« 拡散方程式を高次ルンゲクッタで計算してみる(その1)。 | トップページ | 移流方程式を高次ルンゲクッタで計算してみる(その1) »

2009年12月 9日 (水)

拡散方程式を高次ルンゲクッタで計算してみる(その2)。

昨日は時間微分のみ5次のルンゲクッタで計算して、空間は2次のままだった。バランスとるために時間も高次にしてみよう。

u"(x) ≒ (-u(x-2Δx) + 16u(x-Δx) -30u(x) + 16u(x+Δx) -u(x+2Δx) )/(12*Δx*Δx)

だと4次まで上げられる。で時間は同様に5次で計算してみた。

Diffusion003

普通に計算できてるけど、実は安定性が落ちていて、

空間2次:DΔt/Δx^2 < 0.9

空間4次:DΔt/Δx^2 < 0.7

くらいになった。あちらを立てればこちらが立たず。まさに人生。

« 拡散方程式を高次ルンゲクッタで計算してみる(その1)。 | トップページ | 移流方程式を高次ルンゲクッタで計算してみる(その1) »

学問・資格」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

トラックバック

この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/512682/46860100

この記事へのトラックバック一覧です: 拡散方程式を高次ルンゲクッタで計算してみる(その2)。:

« 拡散方程式を高次ルンゲクッタで計算してみる(その1)。 | トップページ | 移流方程式を高次ルンゲクッタで計算してみる(その1) »

最近のコメント

2018年11月
        1 2 3
4 5 6 7 8 9 10
11 12 13 14 15 16 17
18 19 20 21 22 23 24
25 26 27 28 29 30  
フォト
無料ブログはココログ