熱方程式の最大/最小原理は、クランクニコルソン離散化によって維持されますか?


10

クランクニコルソン有限差分スキームを使用して、1D熱方程式を解いています。熱方程式の最大/最小原理(つまり、最大/最小が初期条件または境界で発生する)も離散化された解に当てはまるかどうか疑問に思っています。

これはおそらく、クランクニコルソンが安定した収束型スキームであることによって暗示されています。しかし、Crank-Nicolsonステンシルから作成された行列を使用して、線形代数引数を介してこれを直接証明できる可能性があるようです。

これに関する文献へのポインタをいただければ幸いです。ありがとう。


こんにちはfoobarbaz、そしてscicompへようこそ!私はあなたが解決している問題にソース用語がないと思いますよね?
ポール

回答:


8

for timestep and grid spacing場合、Crank-Nicolsonの最大原理が成り立ち ます。一般に、の形式の -scheme ここで、は標準のラプラシアン行列で、です。もし、次にスキームが安定しています。(これはフーリエ法で簡単に示すことができます。)ただし、一般に最大原理が成立するためには、というより強い基準が必要です。のkHθUN+1=UN+μ

μkh21
khθA0θ1μ1-2θ1
あなた+1=あなた+μ21θあなた+θあなた+1
0θ1 μ1-θ1μ12θ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 - μを[01]あなたkkあなた0k=あなた2k=0kU N + 1 1 =1-μ

1μ22あなた1+1=1+μ22あなた1
あなた1+1=1μ1+μあなた1

我々はの初期条件を検討している場合、我々は持っている それが常になりますけれどもケースという、我々はそれにもかかわらず持っているだろうという奇数のためない限り、。したがって、ない限り、最大/最小原理に違反します。これは、クランクニコルソンが任意のに対して安定しているという事実に照らして特に注目に値します。あなた10=1

あなた1=1μ1+μ
あなた11あなた1<0μ1μ1μ

foob​​arbazの要求に応えて、証明のスケッチを追加しました。

重要なのは、スキームをの形式で記述することです。

1+2θμあなたj+1=θμあなたj1+1+あなたj+1+1+1θμあなたj1+あなたj+1+[121θμ]あなたj

という仮説は、上記のすべての係数が非負であるという事実とまったく同じです。μ1θ12

あなたj+1あなたj1+1あなたj+1+1あなたj1あなたj+1あなたjあなたj+1あなたj+1の場合、上記の等式と係数の非負性は、

1+2θμあなたj+1>θμあなたj1+1+あなたj+1+1+1θμあなたj1+あなたj+1+[121θμ]あなたj=1+2θμあなたj+1

あなたj+1あなた


ありがとう!たまたまモートン以外の参考文献を知っていますか?Googleブックのプレビューでこれらのセクションまたは定理にアクセスできません。証拠を理解したいのですが。
foob​​arbaz 2013年

@foobarbaz他に便利なリファレンスはありませんが、証明の概要を追加しました。明確にできるかどうか教えてください。
ベン

0

安定性とは、摂動が時間的に制限されたままであることを意味します。最大の原則が離散レベルで満たされるという意味ではありません。それは別の問題です。離散最大原理を満たすことで十分ですが、安定性には必要ありません。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.