情報科学 › 偏微分方程式の離散化 › 差分法の整合性・安定性・収束

4 つのスキームの安定性と精度

定理 7.1(「差分法の整合性・安定性・収束」§7) ・ ✓公理まで証明済み

a≠0a \ne 0、格子比 λ=Δt/h\lambda = \Delta t/h を固定し、ν=aλ\nu = a\lambda とする。u0∈C4(R)u_0 \in C^4(\mathbb{R}) は 4 階までの導関数が有界とし、Um=sup⁡∣u0(m)∣U_m = \sup|u_0^{(m)}|、u(x,t)=u0(x−at)u(x,t) = u_0(x - at) とする。s=sin⁡(θ/2)s = \sin(\theta/2) と書く。

  1. 風上差分(a>0a > 0)。∣g(θ)∣2=1−4ν(1−ν)s2|g(\theta)|^2 = 1 - 4\nu(1-\nu)s^2 で、ℓ2\ell^2 で安定であることと 0<ν≤10 < \nu \le 1 は同値。このとき ∥un+1∥∞≤∥un∥∞\|u^{n+1}\|_\infty \le \|u^n\|_\infty でもある。∣τ∣≤ah2(1+ν)U2|\tau| \le \frac{ah}{2}(1 + \nu)U_2(1 次精度)。
  2. Lax–Friedrichs。∣g(θ)∣2=1−(1−ν2)sin⁡2θ|g(\theta)|^2 = 1 - (1 - \nu^2)\sin^2\theta で、ℓ2\ell^2 で安定であることと ∣ν∣≤1|\nu| \le 1 は同値。このとき ∥un+1∥∞≤∥un∥∞\|u^{n+1}\|_\infty \le \|u^n\|_\infty でもある。∣τ∣≤h2(a2λ+1λ)U2+∣a∣h26U3|\tau| \le \frac{h}{2}\bigl(a^2\lambda + \frac{1}{\lambda}\bigr)U_2 + \frac{|a|h^2}{6}U_3(1 次精度)。
  3. Lax–Wendroff。∣g(θ)∣2=1−4ν2(1−ν2)s4|g(\theta)|^2 = 1 - 4\nu^2(1 - \nu^2)s^4 で、ℓ2\ell^2 で安定であることと ∣ν∣≤1|\nu| \le 1 は同値。∣τ∣≤∣a∣h26(1+ν2)U3+a2Δt h224U4|\tau| \le \frac{|a|h^2}{6}(1 + \nu^2)U_3 + \frac{a^2\Delta t\,h^2}{24}U_4(2 次精度)。
  4. FTCS。∣g(θ)∣2=1+ν2sin⁡2θ|g(\theta)|^2 = 1 + \nu^2\sin^2\theta で、すべての λ>0\lambda > 0 について ℓ2\ell^2 で不安定。∣τ∣≤a2Δt2U2+∣a∣h26U3|\tau| \le \frac{a^2\Delta t}{2}U_2 + \frac{|a|h^2}{6}U_3 なので整合的ではある。

本文で読む(証明と前後の説明)→

閉包の状態

✓ 公理まで証明済み

証明が使う定理を 前提 としてたどると、Rikai の中の証明だけで定義・公理(物理では基本仮定)までつながっています(前提の定理類 90 件、定義なども含めて 105 件)。

印の意味: ✓公理まで証明済み ・ △証明あり(閉包未完)・ ※主張のみ。テーマごとの集計は 閉包の状況。

前提をたどる 証明が使う定理・定義(直接 4 件)

▸ を開くと、その項目の証明が使う前提が出ます。深い先は各項目のページでたどれます。

これを使っている定理 証明で参照している 2 件

関連するノック 4 問