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_gi.gif 811件 [詳細] fileeqa_phi.gif 723件 [詳細] fileeqa_phi_x.gif 769件 [詳細] fileeqa_phii-1.gif 633件 [詳細] fileeqa_phii.gif 710件 [詳細] fileeqa_rho.gif 709件 [詳細] fileeqa_rhoi-12.gif 638件 [詳細] fileeqa_sdi.gif 1074件 [詳細] fileeqa_ai.gif 669件 [詳細] fileeqa_bi.gif 724件 [詳細] fileeqa_ci.gif 708件 [詳細] fileeqa_cipcsl4_1.gif 677件 [詳細] fileeqa_cipcsl4_2.gif 740件 [詳細] fileeqa_cipcsl4_3.gif 687件 [詳細] fileeqa_cipcsl4_4.gif 626件 [詳細] fileeqa_cipcsl4_5.gif 665件 [詳細] fileeqa_cipcsl4_6.gif 665件 [詳細] fileeqa_ei.gif 704件 [詳細] fileeqa_g_x.gif 680件 [詳細] fileeqa_gi-1.gif 779件 [詳細]
トップ   編集 凍結 差分 履歴 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2024-03-08 (金) 18:06:02