CIP-CSL4法*1,では,CIP法の4拘束条件に関数eqa_phi.gifの積分値eqa_rhoi-12.gifをさらに加える.

eqa_cipcsl4_1.gif

つまり,

eqa_cipcsl4_2.gif

の5拘束条件を満足させるために以下の4次多項式をプロファイルに用いるのが CIP-CSL4法である.

eqa_cipcsl4_3.gif

ここで,eqa_ai.gif,eqa_bi.gif,eqa_ci.gif,eqa_sdi.gif,eqa_ei.gifはグリッドi内での多項式補間パラメータである.

5拘束条件eqa_phii.gif,eqa_gi.gif,eqa_phii-1.gif,eqa_gi-1.gif,eqa_rhoi-12.gifから,

eqa_cipcsl4_4.gif

CIP法では算出した係数eqa_ai.gif,eqa_bi.gif,eqa_ci.gif,eqa_sdi.gif,eqa_ei.gifより, eqa_phi_x.gifeqa_g_x.gifを更新していたが, CIP-CSL4法ではさらに積分値eqa_rho.gifの更新が必要となる.

eqa_rho.gifの時間発展は,

eqa_cipcsl4_5.gif

ここで,

eqa_cipcsl4_6.gif

である.


*1 R. Tanaka, T. Nakamura, T. Yabe, "Constructing exactly conservative scheme in a non-conservative form", Comput. Phys. Commun. 126(3), pp.232-243, 2000.

添付ファイル: fileeqa_gi-1.gif 733件 [詳細] fileeqa_gi.gif 756件 [詳細] fileeqa_phi_x.gif 706件 [詳細] fileeqa_g_x.gif 618件 [詳細] fileeqa_sdi.gif 1009件 [詳細] fileeqa_phii.gif 633件 [詳細] fileeqa_phi.gif 649件 [詳細] fileeqa_phii-1.gif 569件 [詳細] fileeqa_rhoi-12.gif 560件 [詳細] fileeqa_ei.gif 654件 [詳細] fileeqa_rho.gif 635件 [詳細] fileeqa_cipcsl4_6.gif 593件 [詳細] fileeqa_bi.gif 672件 [詳細] fileeqa_cipcsl4_1.gif 618件 [詳細] fileeqa_ci.gif 652件 [詳細] fileeqa_cipcsl4_5.gif 609件 [詳細] fileeqa_cipcsl4_4.gif 550件 [詳細] fileeqa_cipcsl4_2.gif 673件 [詳細] fileeqa_ai.gif 624件 [詳細] fileeqa_cipcsl4_3.gif 623件 [詳細]

トップ   編集 凍結 差分 履歴 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2011-10-28 (金) 13:06:23