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.gifとeqa_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_sdi.gif 1043件 [詳細] fileeqa_phii.gif 679件 [詳細] fileeqa_g_x.gif 648件 [詳細] fileeqa_phi_x.gif 739件 [詳細] fileeqa_gi-1.gif 747件 [詳細] fileeqa_phi.gif 689件 [詳細] fileeqa_gi.gif 780件 [詳細] fileeqa_rhoi-12.gif 609件 [詳細] fileeqa_rho.gif 677件 [詳細] fileeqa_phii-1.gif 600件 [詳細] fileeqa_cipcsl4_5.gif 633件 [詳細] fileeqa_cipcsl4_6.gif 630件 [詳細] fileeqa_cipcsl4_2.gif 707件 [詳細] fileeqa_ai.gif 640件 [詳細] fileeqa_ci.gif 676件 [詳細] fileeqa_cipcsl4_3.gif 657件 [詳細] fileeqa_ei.gif 674件 [詳細] fileeqa_bi.gif 694件 [詳細] fileeqa_cipcsl4_1.gif 648件 [詳細] fileeqa_cipcsl4_4.gif 596件 [詳細]
トップ   編集 凍結 差分 履歴 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2022-11-30 (水) 13:48:06