CIP-CSL2法*1,*2,*3は, 値eqa_phi_x.gifとその積分値eqa_cipcsl2_1.gifを用いる. 積分値eqa_D.gifに対する移流方程式は,

eqa_cipcsl2_2.gif

$D$の微分eqa_cipcsl2_3.gifと定義すると

eqa_cipcsl2_4.gif

これは1次元の移流方程式と同じであり, CIP法でのeqa_cipcsl2_5.gifと置き換えることで CIP法と同様に計算ができる.

格子点i内でのeqa_Di.gifを以下のように定義する.

eqa_cipcsl2_6.gif

eqa_Di.gifeqa_xi.gifから補間点までの積分値を表し, 3次関数で補間する.

eqa_cipcsl2_7.gif

CIP法における空間微分値gはeqa_phi.gifで置き換えられるので, [eqa_xi.gif, eqa_xi+1.gif]のプロファイルは,

eqa_cipcsl2_8.gif

eqa_D.gifの定義より,

eqa_cipcsl2_9.gif

また,

eqa_cipcsl2_10.gif

これらの条件により係数eqa_ai.gif,eqa_bi.gifは以下のように求められる.

eqa_cipcsl2_11.gif

eqa_rho.gifの時間発展は,eqa_xi.gif,eqa_xi+1.gifをバックトレースした位置eqa_cipcsl2_12.gifに挟まれた部分の積分値となる(下図).

eqa_cipcsl2_13.gif
csl2_int.gif

*1 T. Yabe, R. Tanaka, T. Nakamura, F. Xiao, "An exactly conservative semi-Lagrangian scheme (CIP-CSL) in one dimension", Mon. Weather Rev. 129, pp.332-344, 2001.
*2 T. Nakamura, R. Tanaka, T. Yabe and K. Takizawa, "Exactly Conservative Semi-Lagrangian Scheme for Multi-dimensional Hyperbolic Equations with Directional Splitting Technique", J. Comput. Phys. 174, pp.171-207, 2001.
*3 K. Takizawa, T. Yabe, T. Nakamura, "Multi-dimensional semi-Lagrangian scheme that guarantees exact conservation", Comput. Phys. Commun. 148(2), pp.137-159, 2002.

添付ファイル: fileeqa_cipcsl2_4.gif 540件 [詳細] fileeqa_phi.gif 581件 [詳細] fileeqa_cipcsl2_3.gif 519件 [詳細] fileeqa_rho.gif 539件 [詳細] fileeqa_xi.gif 576件 [詳細] fileeqa_cipcsl2_8.gif 503件 [詳細] fileeqa_cipcsl2_9.gif 463件 [詳細] fileeqa_cipcsl2_5.gif 484件 [詳細] fileeqa_cipcsl2_6.gif 464件 [詳細] fileeqa_phi_x.gif 489件 [詳細] fileeqa_cipcsl2_7.gif 450件 [詳細] fileeqa_xi+1.gif 469件 [詳細] fileeqa_Di.gif 589件 [詳細] fileeqa_cipcsl2_2.gif 467件 [詳細] fileeqa_ai.gif 543件 [詳細] fileeqa_cipcsl2_10.gif 467件 [詳細] filecsl2_int.gif 1362件 [詳細] fileeqa_cipcsl2_13.gif 441件 [詳細] fileeqa_cipcsl2_11.gif 454件 [詳細] fileeqa_bi.gif 556件 [詳細] fileeqa_D.gif 560件 [詳細] fileeqa_cipcsl2_12.gif 457件 [詳細] fileeqa_cipcsl2_1.gif 483件 [詳細]

トップ   編集 凍結 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2011-10-28 (金) 13:07:05 (3257d)