[论文解读] 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.

研究动机与目标

• 为对称（圆形）空间区域内单种生物种群动态构建非局部时滞反应-扩散模型。
• 解决尽管已有越来越多关于此类栖息地的实证研究，但在对称区域内建模仍显不足的问题。
• 通过使用贝塞尔函数表示出生函数，将完整的年龄结构模型简化为径向对称形式。
• 对不同初始条件和参数组合下的简化模型进行数值模拟与长期行为分析。
• 展示空间区域几何形状与边界条件对模型推导及解动态的影响。

提出的方法

• 通过在径向对称条件下简化年龄结构模型，推导出二维圆形区域内非局部时滞反应-扩散模型。
• 假设无方位依赖性（因各向同性扩散），将偏微分方程简化为径向形式。
• 在简化模型中使用第一类贝塞尔函数表示出生函数。
• 施加零通量边界条件，并使用有限差分法或谱方法对所得偏微分方程进行数值求解。
• 对同质（无出生函数）与非同质（含密度依赖性出生函数）两种情况分别进行数值模拟。
• 通过随时间比较解的变化，评估其是否收敛至平凡解（灭绝）或正平衡态。

实验结果

研究问题

• RQ1 对称（圆形）区域的几何形状如何影响非局部时滞反应-扩散模型在单种生物扩散中的推导与结构？
• RQ2 当出生函数缺失与存在且具有密度依赖性时，圆形区域内的种群动态将如何变化？
• RQ3 边界条件与区域大小如何影响种群模型的长期行为？
• RQ4 能否通过贝塞尔函数将模型简化为径向对称形式？该简化的意义是什么？
• RQ5 在不同参数组合与初始条件下，种群密度随时间演变将呈现出何种数值模式？

主要发现

• 当无出生函数时，种群密度均匀收敛至零，表明发生灭绝，且在 t > 400 时最大密度降至 5.3×10⁻⁶。
• 当采用密度依赖性出生函数 b(w) = 0.25w²exp(−0.1w) 时，种群稳定在正的常数平衡态。
• 模型行为对出生函数的存在极为敏感：若无出生函数，由于高死亡率且无繁殖，种群将灭绝。
• 数值模拟证实，即使存在扩散与时滞，当出生函数缺失时系统仍收敛至平凡解。
• 引入非线性出生函数可导致稳定且正的种群平衡态，表明种群具有持续性。
• 当区域半径 R → ∞ 时，所推导模型正确地退化为整个平面上的标准非局部时滞反应-扩散模型。