数値計算メモ


  1. semi-implicit スキーム

semi-implicit スキーム

数値モデルを陽解法で計算しようとすれば, 系が含む最も速い伝達速度を持つ波動 によって時間ステップが制約される. プリミティブ方程式系が含む最も速い波は 外部重力波であり, 位相速度は 300 m/s 程度である. このような波動に対し, 陽解法の Courant-Fredrichs-Lewy (CFL) 条件を満たすためには, 短いステップ Δt をとる必要がある. 格子間隔を Δx とすると, CFL 条件は一般に,

begin LATEX

\begin{eqnarray}

\frac{c\delta t}{\delta x} \leq 1

\end{eqnarray}