for timestep and grid spacing場合、Crank-Nicolsonの最大原理が成り立ち
ます。一般に、の形式の -scheme
ここで、は標準のラプラシアン行列で、です。もし、次にスキームが安定しています。(これはフーリエ法で簡単に示すことができます。)ただし、一般に最大原理が成立するためには、というより強い基準が必要です。のkHθUN+1=UN+μ
μ ≐ Kh2≤ 1
khθA0≤θ≤1μ(1-2θ)≤1あなたn + 1= uん+ μ2((1 - θ )A Uん+ θ A Un + 1)
あ0 ≤ θ ≤ 1 μ(1-θ)≤1μ (1 - 2 θ )≤ 12μ (1 - θ )≤ 12
証明については、KW Mortonによる偏微分方程式の数値解を参照してください。特に、セクション2.10および2.11と定理2.2をご覧ください。
制約がない場合、最大原理がCrank-Nicolsonに対して一般的に成り立たないことを確認する良い方法もあります。μ
境界を含む3点を含む離散化を使用して、の熱方程式を考えます。タイムステップとグリッドポイントでの離散化を表すとしましょう。すべてのとなるようにディリクレ境界を仮定します。次に、クランクニコルソンは減少し
これはさらに削減できます
Uはkはiは kをI uはkは0 = uはkは2 = 0 K ( 1 - μを[ 0 、1 ]あなたk私k私あなたk0= uk2= 0kU N + 1 1 =(1-μ
( 1 - μ2(− 2 ))un + 11= ( 1 + μ2(− 2 ))uん1、
あなたn + 11= (1 - μ1 + μ)uん1。
我々はの初期条件を検討している場合、我々は持っている
それが常になりますけれどもケースという、我々はそれにもかかわらず持っているだろうという奇数のためない限り、。したがって、ない限り、最大/最小原理に違反します。これは、クランクニコルソンが任意のに対して安定しているという事実に照らして特に注目に値します。あなた01= 1
あなたん1= (1 - μ1 + μ)ん、
あなたん1≤ 1あなたん1< 0んμ ≤ 1μ ≤ 1μ
foobarbazの要求に応えて、証明のスケッチを追加しました。
重要なのは、スキームをの形式で記述することです。
(1 + 2 θ μ )のun + 1j= θ μ (Un + 1j − 1+ un + 1j + 1)+ (1 - θ )μ (Uんj − 1+ uんj + 1)+ [ 1 - 2 (1 - θ )μ ] Uんj
という仮説は、上記のすべての係数が非負であるという事実とまったく同じです。μ (1 - θ )≤ 12
あなたn + 1jあなたn + 1j − 1あなたn +1j + 1あなたんj − 1あなたんj + 1あなたんjあなたn + 1jあなたn + 1jの場合、上記の等式と係数の非負性は、
(1 + 2 θ μ )のun + 1j> θ μ (Un + 1j − 1+ un + 1j + 1)+ (1 - θ )μ (Uんj − 1+ uんj + 1)+ [ 1 - 2 (1 - θ )μ ] Uんj= (1 + 2 θ μ )Un + 1j
あなたn + 1jあなた