[論文レビュー] Modeling and Numerical Simulations of Single Species Dispersal in Symmetrical Domains

Majid Bani-Yaghoub, Guangming Yao+1 more|arXiv (Cornell University)|Dec 15, 2014

Mathematical and Theoretical Epidemiology and Ecology Models参考文献 18被引用数 7

ひとこと要約

本稿では、円形（対称的）領域における単一種の拡散を対象とした非局所遅延反応拡散モデルを構築し、数値的にシミュレートした。問題は径方向対称性に還元され、出生関数にはベッセル関数が用いられた。主な発見は、出生関数が存在しない場合には種の絶滅が生じるが、密度に依存する出生関数が存在する場合には正の平衡状態に収束することである。

ABSTRACT

We develop a class of nonlocal delay Reaction-Diffusion (RD) models in a circular domain. Previous modeling efforts include RD population models with respect to one-dimensional unbounded domain, unbounded strip and rectangular spatial domain. However, the importance of an RD model in a symmetrical domain lies in the increasing number of empirical studies conducted with respect to symmetrical natural habitats of single species. Assuming that the single species has no directional preference to spread in the symmetrical domain, the RD model is reduced to an equation with no angular dependance. The model can be further reduced by considering the birth function in the form of the Bessel function of the first kind. We numerically simulate the reduced forms of the nonlocal delay RD model to study the dispersal and growth of behaviors of the single species in a circular domain. Although spatial patterns of population densities are gradually developed, it is numerically shown that the single species population goes extinct in the absence of the birth function or it may converge to a positive equilibrium in the presence of the birth function.

研究の動機と目的

• 対称的（円形）空間領域における単一種の集団動態の非局所遅延反応拡散モデルの構築を目的とする。
• このような生息環境に関する実証的研究が増加しているにもかかわらず、そのような領域におけるモデル化の不足に取り組むことを目的とする。
• 出生関数にベッセル関数を用いることで、年齢構造モデルを径方向対称形に簡略化することを目的とする。
• 異なる初期条件およびパラメータ設定下での、簡略化モデルの長期的挙動の数値的シミュレーションと解析を目的とする。
• 空間領域の幾何構造および境界条件がモデルの導出と解のダイナミクスに与える影響を示すことを目的とする。

提案手法

• 径方向対称性の下で年齢構造モデルを簡略化することにより、2次元の円形領域における非局所遅延反応拡散モデルを導出する。
• 等方的拡散のための角度依存性がないと仮定し、偏微分方程式を径方向形式に簡略化する。
• 簡略化モデルにおける出生関数として、第一種ベッセル関数を用いる。
• ゼロフラックス境界条件を適用し、有限差分法またはスペクトル法を用いて得られた偏微分方程式を数値的に解く。
• 均一（出生関数なし）および非均一（密度に依存する出生関数あり）の両ケースについて数値的シミュレーションを実施する。
• 時間経過に伴う解の比較を通じて、自明解（絶滅）または正の平衡状態への収束を評価する。

実験結果

リサーチクエスチョン

• RQ1 対称的（円形）領域の幾何構造が、単一種の拡散のための非局所遅延反応拡散モデルの導出および構造にどのように影響を与えるか？
• RQ2 出生関数が存在しない場合と、密度に依存する出生関数が存在する場合とで、円形領域における集団動態はどのように変化するか？
• RQ3 境界条件および領域サイズが、集団モデルの長期的挙動にどのように影響を与えるか？
• RQ4 ベッセル関数を用いることで、モデルを径方向対称形に簡略化できるか？また、その簡略化の意義は何か？
• RQ5 異なるパラメータ設定および初期条件下で、時間経過に伴う集団密度にどのような数値的パターンが現れるか？

主な発見

• 出生関数が存在しない場合、集団密度は均一にゼロに収束し、絶滅を示す。t > 400 の時点で最大密度は 5.3×10⁻⁶ にまで低下した。
• 密度に依存する出生関数 b(w) = 0.25w²exp(−0.1w) が存在する場合、集団は正の定常状態に到達する。
• モデルの挙動は出生関数の有無に極めて敏感である。出生関数がなければ、高い死亡率と再生不能性により絶滅が生じる。
• 数値的シミュレーションにより、出生関数が存在しない場合でも、拡散と遅延が存在するにもかかわらず、システムが自明解に収束することが確認された。
• 非線形の出生関数を組み込むことで、安定した正の集団平衡状態が得られ、集団の持続可能性が示された。
• 領域半径 R → ∞ の極限において、導出されたモデルは、平面全体における標準的な非局所遅延RDモデルに正しく還元される。