[논문 리뷰] High order discretization methods for spatial-dependent SIR models

Bálint Takács, Yiannis Hadjimichael|arXiv (Cornell University)|2021. 01. 01.

Mathematical and Theoretical Epidemiology and Ecology Models인용 수 5

한 줄 요약

이 논문은 적분미분방정식으로 지배되는 공간에 따라 변화하는 SIR 모델을 위한 고차수 공간 이산화 방법을 제안하며, 수치 해에서 보존성, 양성, 단조성을 보장한다. 고차수 스킴이 주요 모델 성질을 유지하기 위한 충분조건을 도출하고, 계산 실험을 통해 정확도와 수렴성을 검증한다.

ABSTRACT

In this paper, an SIR model with spatial dependence is studied and results regarding its stability and numerical approximation are presented. We consider a generalization of the original Kermack and McKendrick model in which the size of the populations differs in space. The use of local spatial dependence yields a system of integro-differential equations. The uniqueness and qualitative properties of the continuous model are analyzed. Furthermore, different choices of spatial and temporal discretizations are employed, and step-size restrictions for population conservation, positivity, and monotonicity preservation of the discrete model are investigated. We provide sufficient conditions under which high order numerical schemes preserve the discrete properties of the model. Computational experiments verify the convergence and accuracy of the numerical methods.

연구 동기 및 목표

• 공간에 따라 다름을 보이는 인구의 이질성과 함께 공간적으로 의존하는 SIR 모델의 유일성과 정성적 행동을 분석하기 위해.
• 유도된 적분미분방정식계를 위한 고차수 시간 및 공간 이산화 방법을 개발하기 위해.
• 이산 해가 인구 보존성, 양성, 단조성을 유지하도록 보장하기 위한 스텝 크기 제약 조건을 유도하기 위해.
• 고차수 스킴이 연속 모델의 구조적 성질을 유지할 수 있는 충분조건을 설정하기 위해.

제안 방법

• 공간적으로 변화하는 인구 규모를 고려한 일반화된 Kermack-McKendrick SIR 프레임워크를 사용하여 인구 역학을 모델링하기 위해.
• 지역적 공간 의존성을 반영하기 위해 모델을 적분미분방정식계로 공식화하기 위해.
• 고차수 시간 통합 방법(예: 룬게-쿠타 방법)을 사용하고, 유한차분법 또는 스펙트럼 방법을 통해 공간 이산화를 수행하기 위해.
• 이산 보존성, 양성, 단조성을 유지하기 위한 스텝 크기 제약 조건에 기반한 안정성 조건을 도출하기 위해.
• 제안된 스킴의 수렴성 및 정확도를 검증하기 위해 계산 실험을 수행하기 위해.
• 다양한 공간 및 시간 이산화 선택에 따른 이산 모델의 행동 분석하기 위해.

실험 결과

연구 질문

• RQ1 고차수 수치 스킴이 공간에 따라 변화하는 SIR 모델에서 인구 보존성을 유지하기 위한 조건은 무엇인가?
• RQ2 시간 단계를 거쳐 이산 해의 양성과 단조성을 어떻게 유지할 수 있는가?
• RQ3 이산 모델에서 보존성 및 양성을 유지하기 위해 필요한 스텝 크기 제약 조건은 무엇인가?
• RQ4 다양한 공간 및 시간 이산화 선택이 수치 해의 정확도와 수렴성에 미치는 영향은 어떠한가?
• RQ5 고차수 스킴이 연속 공간 SIR 모델의 정성적 행동을 어느 정도 유지하는가?

주요 결과

• 고차수 수치 스킴이 이산 인구 총량의 보존성을 유지할 수 있는 충분조건가 도출되었다.
• 이산 해의 양성과 단조성을 보장하기 위한 스텝 크기 제약 조건이 설정되었다.
• 제안된 고차수 스킴은 공간과 시간 모두에서 최적의 수렴 속도를 보이며, 이론적 기대를 확인하였다.
• 계산 실험을 통해 다양한 공간 구성에서 방법의 정확성과 안정성이 확인되었다.
• 적절한 고차수 이산화를 적용할 경우 이산 모델이 연속 시스템의 핵심 정성적 특성을 유지한다.
• 분석 결과 고차수 방법이 생물학적 제약 조건을 위반하지 않으면서도 공간 이질성 SIR 모델의 역학을 효과적으로 포착할 수 있음을 보여주었다.