[論文レビュー] A distribution-free lattice Boltzmann method for compartmental reaction-diffusion systems with application to epidemic modelling

Alessandro De Rosis|arXiv (Cornell University)|Mar 20, 2026

Lattice Boltzmann Simulation Studies被引用数 0

ひとこと要約

要約: 論文は SEIRD 反応拡散システムの直接 macroscopic compartment densities を進化させる単一ステップの簡略化された格子Boltzmann法（SSLBM）を導入し、第四階の有限差分ソルバーと検証し、BGK LBM と比較してあらゆる領域で精度が向上することを示す。

ABSTRACT

We introduce a distribution-free lattice Boltzmann formulation for general compartmental reaction--diffusion systems arising in mathematical epidemiology. The proposed scheme, termed a single-step simplified lattice Boltzmann method (SSLBM), evolves directly macroscopic compartment densities, eliminating the need for particle distribution functions and explicit streaming operations. This yields a compact and computationally efficient framework while retaining the kinetic consistency of lattice Boltzmann methodologies. The approach is applied to a SEIRD (Susceptible-Exposed-Infected-Recovered-Deceased) reaction-diffusion model as a representative case. The resulting discrete evolution equations are derived and shown to recover the target macroscopic dynamics. The method is systematically validated against a fourth-order finite difference reference solution and compared with a standard BGK lattice Boltzmann formulation. Numerical results demonstrate that the SSLBM consistently improves accuracy across all compartments and norms. The error reduction is robust with respect to both the basic reproduction number and diffusion strength, typically ranging between factors of approximately two and five depending on the regime. In particular, the method shows enhanced control of localised errors in regimes characterised by strong nonlinear coupling and steep spatial gradients. Our findings indicate that the proposed formulation provides an accurate and efficient alternative to classical lattice Boltzmann approaches for reaction-diffusion systems, with particular advantages in stiff and nonlinear epidemic dynamics.

研究の動機と目的

• 難病モデルを ODE のみでなく compartmental PDE を用いて空間的に拡張する動機付け。
• PDFなしで macroscopic densities を進化させ、メモリ効率の良い SSLBM を開発する。
• 高次有限差分法と比較して一貫性、安定性、質量保存、精度を示す。
• SEIRD への適用と拡散レジームおよび非線形結合全体での性能分析。

提案手法

• 拡散係数を持つ S, E, I, R 区分と D の ODE を含む一般的な compartmental 反応–拡散 PDE を定式化。
• 各移動可能区分に対する BGK 格子ボルツマン表現を導出し、拡散係数を d=c_s^2(τ−1/2) に結びつけた緩和関係を導出。
• バックワード/フォワードのスタencil を用いて密度更新を計算する単一ステップ SSLBM 更新と、予測子-修正子を用いた単一パスへの還元を導入。
• Explicit SSLBM 更新: ρ^C(x,t+Δt)=ρ_b^C+(τ^C−1)(ρ_f^C−2ρ_c^C+ρ_b^C)+Ψ^C(x,t)Δt (Eq. 18).
• SSLBM の拡張演算子は構造上第2次中心差分ラプラシアンと一致するが、運動論的導出と質量保存を保持。
• 整合性（二次空間）、安定性（DFE 付近の拡散と線形化反応安定性）、各方面性（D2Q5 ステンシル）、および離散レベルでの厳密な質量保存を分析。

(a) 0.01

実験結果

リサーチクエスチョン

• RQ1 SSLBM は高次有限差分参照と比較して SEIRD 反応–拡散ダイナミクスを正確に再現するか。
• RQ2 diffusion レジームと初期撹乱の異なる状況で、BGK LBM と比べた精度とメモリフットプリントの性能はどうか。
• RQ3 区分モデルに対する SSLBM の安定性、異方性、質量保存特性は。
• RQ4 SEIRD を超える一般的な compartmental PDE へ SSLBM フレームワークを拡張可能か。
• RQ5 拡散の強さと非線形結合は SSLBM と BGK LBM の局所誤差制御にどう影響するか。

主な発見

χ	Compartment	L2-norm SSLBM	L2-norm BGK	Ratio (BGK/SSLBM)	L∞-norm SSLBM	L∞-norm BGK	Ratio (BGK/SSLBM)
0.01	S	0.0364	0.0624	1.71	0.0482	0.0799	1.66
0.01	E	0.0719	0.1208	1.68	0.0755	0.1154	1.53
0.01	I	0.0711	0.1158	1.63	0.0750	0.1097	1.46
0.01	R,D	0.0284	0.0475	1.67	0.0431	0.0696	1.61
0.05	S	0.0241	0.0603	2.50	0.0252	0.0603	2.39
0.05	E	0.0384	0.0856	2.23	0.0468	0.0856	1.83
0.05	I	0.0413	0.0836	2.02	0.0511	0.0836	1.64
0.05	R,D	0.0137	0.0508	3.71	0.0233	0.0508	2.18
0.10	S	0.0204	0.0591	2.90	0.0162	0.0549	3.39
0.10	E	0.0303	0.0692	2.28	0.0369	0.0740	2.01
0.10	I	0.0342	0.0717	2.10	0.0436	0.0994	2.28
0.10	R,D	0.0104	0.0239	2.30	0.0158	0.0441	2.79
0.20	S	0.0234	0.0696	2.97	0.0121	0.0609	5.03
0.20	E	0.0291	0.0601	2.07	0.0345	0.1168	3.39
0.20	I	0.0332	0.0725	2.18	0.0384	0.1811	4.72
0.20	R,D	0.0094	0.0212	2.26	0.0135	0.0424	3.14

• SSLBM は SEIRD ダイナミクスを第四階の有限差分参照と優れた一致度で再現し、ピークの時刻と振幅を含む。
• 初期曝露割合 χ ∈ {0.01,0.05,0.10,0.20} に対して、SSLBM は全区分と各ノルムで一貫して BGK LBM を上回り、精度向上は χ が大きくなるにつれて拡大。
• χ=0.01 での誤差削減はおおよそ 1.5–1.7 倍、多くのレジームで 2–5 倍、χ=0.20 では L∞ で >4 に達することもある。
• 方法は離散レベルで総人口を厳密に保存し、空間精度は第2次で、拡張は τ≥1/2 かつ Δt ≤ 4/((α+γ)−√((α−γ)^2+4αβ)) の標準 BGK LBM に整合する安定条件を示す。
• 異方性解析では D2Q5 拡散演算子は一階近似で等方性だが四次の異方性項を持つ。SSLBM は運動学的フレームワークを介してより複雑な輸送へ拡張可能。

(b) 0.05