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_sdi.gif 986件 [詳細] fileeqa_phii.gif 623件 [詳細] fileeqa_g_x.gif 588件 [詳細] fileeqa_phi_x.gif 687件 [詳細] fileeqa_gi-1.gif 697件 [詳細] fileeqa_phi.gif 635件 [詳細] fileeqa_gi.gif 722件 [詳細] fileeqa_rhoi-12.gif 553件 [詳細] fileeqa_rho.gif 614件 [詳細] fileeqa_phii-1.gif 543件 [詳細] fileeqa_cipcsl4_5.gif 585件 [詳細] fileeqa_cipcsl4_6.gif 574件 [詳細] fileeqa_cipcsl4_2.gif 657件 [詳細] fileeqa_ai.gif 589件 [詳細] fileeqa_ci.gif 622件 [詳細] fileeqa_cipcsl4_3.gif 599件 [詳細] fileeqa_ei.gif 621件 [詳細] fileeqa_bi.gif 643件 [詳細] fileeqa_cipcsl4_1.gif 589件 [詳細] fileeqa_cipcsl4_4.gif 538件 [詳細]

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