収束とは何か？数値計算で結果が安定する条件

数値計算における収束とは、刻み幅や格子幅、反復回数などの計算条件を変えたときに、数値解がある値へ安定して近づくことを意味します。

例えば時間刻み $\Delta t$ を小さくしたり、格子幅 $\Delta x$ を細かくしたり、反復回数を増やしたときに、解がほぼ同じ値へ落ち着くなら、その計算は収束していると考えられます。

逆に、計算条件を変えるたびに結果が大きく変わるなら、その計算結果はまだ信用できません。

つまり収束とは、単に「計算が終わること」ではなく、計算結果が安定して近づくことです。

数値計算では、プログラムが最後まで動けば答えが出ます。しかしその値が本当に正しいかどうかは、別の問題です。

刻み幅を変えた瞬間に結果が変わることは珍しくありません。

1) 収束の基本概念

数学では、数列 $\{x_n\}$ がある値 $a$ に近づくとき

$x_n \to a$

と書き、これを収束と呼びます。

数値計算でも同じ考え方が使われます。ただし、近づく対象は反復回数だけではありません。

• 反復回数
• 時間刻み $\Delta t$
• 格子幅 $\Delta x$

これらを変えたときに解が落ち着くかどうかを確認します。

2) 反復法の収束

多くの数値計算では反復法が使われます。例えば解を更新する形で

$x^{(k+1)} = g(x^{(k)})$

のように計算を進めます。

このとき更新量が十分小さくなれば収束と判断します。

典型的な条件は

$|x^{(k+1)}-x^{(k)}| < \varepsilon$

です。

3) 時間刻みと収束

時間発展問題では時間刻みが重要です。

例えば陽的オイラー法では

$u^{n+1} = u^n + \Delta t f(u^n,t^n)$

のように解を更新します。

もし $\Delta t$ が大きすぎると、解が振動したり発散することがあります。

そこで $\Delta t$ を小さくして計算をやり直し、解が落ち着くかを確認します。

4) 格子収束

偏微分方程式を解くときは空間も離散化します。

格子幅 $\Delta x$ が粗いと、流れや温度分布を正しく表現できません。

そこで格子を細かくして再計算し、結果が変わらなくなるかを確認します。

これを格子収束またはメッシュ収束と呼びます。

5) 安定と収束の違い

数値計算では安定性と収束はよく混同されますが、意味は異なります。

数値解析では、整合性と安定性が満たされると収束するという重要な結果があります。

6) 計算例

時間刻みを変えて計算した結果が次のようになったとします。

$\Delta t=0.10 \quad u=1.42$

$\Delta t=0.05 \quad u=1.36$

$\Delta t=0.025 \quad u=1.33$

$\Delta t=0.0125 \quad u=1.32$

刻みを細かくするほど解が一定値に近づいています。

この場合、解はおよそ

$u \approx 1.32$

へ収束していると考えられます。

7) よくあるつまずき

数値計算でよくある誤解は、プログラムが最後まで動けば計算が成功したと思ってしまうことです。

しかし収束していない計算では、結果は刻み幅に強く依存します。

そのため、刻み幅や格子幅を変えて結果が安定するかを必ず確認する必要があります。