[論文レビュー] High order discretization methods for spatial-dependent SIR models

Bálint Takács, Yiannis Hadjimichael|arXiv (Cornell University)|Jan 1, 2021

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モデルのダイナミクスを効果的に捉えることができることが示された。