[논문 리뷰] Modeling and Numerical Simulations of Single Species Dispersal in Symmetrical Domains

Majid Bani-Yaghoub, Guangming Yao외 1명|arXiv (Cornell University)|2014. 12. 15.

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차원 원형 영역에서 비국소 지연 반응-확산 모델을 유도한다.
• 등방성 확산에 기인한 각도 의존성이 없음을 가정하여 PDE를 반경 방향 형태로 단순화한다.
• 감소된 모델에서 생식 함수를 제1종 베르누이 함수로 표현한다.
• 0-유속 경계 조건을 적용하고, 유한 차분 또는 스펙트럼 방법을 사용하여 유도된 PDE를 수치적으로 해결한다.
• 균일한(생식 함수 없음) 및 비균일한(밀도에 따라 변하는 생식 함수 있음) 경우에 대해 수치 시뮬레이션을 수행한다.
• 시간에 따른 해를 비교하여 자취(멸종) 또는 정상 상태로 수렴하는지 평가한다.

실험 결과

연구 질문

• RQ1 대칭(원형) 영역의 기하학적 특성이 단일 종 확산을 위한 비국소 지연 반응-확산 모델의 유도 및 구조에 어떻게 영향을 미치는가?
• RQ2 생식 함수가 없을 경우와 밀도에 따라 변하는 생식 함수가 있을 경우 원형 영역에서 인구 역학은 어떻게 변화하는가?
• RQ3 경계 조건과 영역 크기가 인구 모델의 장기적 행동에 어떻게 영향을 미치는가?
• RQ4 베르누이 함수를 사용하여 모델을 반경 방향 대칭 형태로 단순화할 수 있으며, 이러한 단순화의 의미는 무엇인가?
• RQ5 다양한 파라미터 조합과 초기 조건 하에서 시간에 따른 인구 밀도에서 나타나는 수치적 패턴은 무엇인가?

주요 결과

• 생식 함수가 없을 경우, 인구 밀도는 균일하게 0으로 수렴하며 멸종을 나타내며, t > 400일 때 최대 밀도가 5.3×10⁻⁶로 감소한다.
• 밀도에 따라 변하는 생식 함수 b(w) = 0.25w²exp(−0.1w)가 존재할 경우, 인구는 정상 상태의 양의 일정한 값으로 정착한다.
• 모델의 행동은 생식 함수의 유무에 매우 민감하다: 생식 함수가 없을 경우 고사율과 번식 없음으로 인해 멸종이 발생한다.
• 수치 시뮬레이션은 생식 함수가 없을 경우에도 확산과 지연가 있음에도 불구하고 시스템이 자취 해로 수렴함을 확인한다.
• 비선형 생식 함수의 포함은 안정적인 양의 인구 평형 상태를 이끌어내며, 인구 지속성을 입증한다.
• 영역 반경 R → ∞일 때 유도된 모델은 전체 평면에서 표준 비국소 지연 반응-확산 모델로 정확히 단순화된다.