Navier-Stokes方程式での移流項 (u・▽)u について

----
#contents
----

**移流方程式 [#v073943b]
移流方程式
#ref(eqa_advec_equation1.gif,,50%,移流方程式)

*空間微分の離散化 [#h165c435]
**風上差分(Upwind differencing) [#tac8e468]
#include(風上差分,title)

-例:矩形波の移流
#include(矩形波 - 設定,notitle)
風上差分では下の中心差分のように振動が全くない代わりに,形状がなまっている.
これは数値拡散と呼ばれる現象によるためである.
#include(矩形波 - 風上差分,notitle)


**中心差分(Central differencing) [#pa3383ab]
#include(中心差分,title)

-例:矩形波の移流
#include(矩形波 - 設定2,notitle)
風上差分ではt=5まで移流させた結果のみ示したが,中心差分では非常に振動が大きくなるため,
t=0.5, 1.0, 5.0の3つのグラフを示した.流れの後方に多くの振動が現れ成長していく様子が分かる.
#include(矩形波 - 中心差分,notitle)


**Lax-Wendroff法 [#g3a9f016]
#include(Lax-Wendroff法,title)

-例:矩形波の移流
#include(矩形波 - 設定2,notitle)
#include(矩形波 - Lax-Wendroff法,notitle)


**QUICK(Quadratic Upstream Interpolation for Convective Kinematics) [#la926299]
#include(QUICK,title)

-例:矩形波の移流
#include(矩形波 - 設定2,notitle)
#include(矩形波 - QUICK,notitle)


**QUICKEST(QUICK with estimated streaming terms) [#u6eb7363]
#include(QUICKEST,title)

-例:矩形波の移流
#include(矩形波 - 設定2,notitle)
#include(矩形波 - QUICKEST,notitle)


**河村・桑原スキーム(Kawamura-Kuwabara scheme) [#y2dbc484]
#include(河村・桑原スキーム,title)

-例:矩形波の移流
#include(矩形波 - 設定2,notitle)
#include(矩形波 - 河村・桑原スキーム,notitle)



**ENO(Essentially Non-Oscillatory polynomial interpolation) [#ufe6b856]
#include(ENO,title)

-例:矩形波の移流
#include(矩形波 - 設定,notitle)
#include(矩形波 - ENO,notitle)



**WENO(Weighted Essentially Non-Oscillatory polynomial interpolation) [#ofec5d0e]
#include(WENO,title)

-例:矩形波の移流
#include(矩形波 - 設定,notitle)
#include(矩形波 - WENO,notitle)



**セミラグランジュ(semi-Lagrangian)法 [#d5baf792]
#include(セミラグランジュ法,title)

-例:矩形波の移流
#include(矩形波 - 設定,notitle)
#include(矩形波 - セミラグランジュ法,notitle)



**BFECC(Back and Forth Error Compensation and Correction) [#q7fba469]
#include(BFECC,title)

-例:矩形波の移流
#include(矩形波 - 設定,notitle)
#include(矩形波 - BFECC,notitle)



**MacCormack [#w4d9a820]
参考文献1((A. Selle, R. Fedkiw, B. Kim, Y. Liu, J. Rossignac, "An unconditionally stable MacCormack method", J. Sci. Comput. 35(2), pp.350-371, 2008.)),
参考文献2((R. MacCormack, "The effect of viscosity in hypervelocity impact cratering", AIAA Hypervelocity Impact Conference, 1969.))




**CIP法(Constrained Interpolation Profile scheme) [#h2dafb81]
#include(CIP法,title)

-例
#include(矩形波 - 設定,notitle)
#include(矩形波 - CIP,notitle)



**RCIP法(Rational Constrained Interpolation Profile scheme) [#f1688623]
#include(RCIP法,title)

-例
#include(矩形波 - 設定,notitle)
#include(矩形波 - TCIP,notitle)
#include(矩形波 - RCIP,notitle)



**CIP-CSL法 [#df441b69]
CIP法では格子iのプロファイル&ref(eqa_phii.gif,,70%);の決定のために,
4拘束条件から3次多項式の4係数を求めたが,
この条件では値とその勾配値は保存されるが体積は保存されない.
保存系CIP移流法を用いることでこれを解決できる.
保存系でのCIP法の実装として,CIP-CSL2法やCIP-CSL4法, CIP-CSL3法などが提案されている.

***CIP-CSL4 [#x72994d6]
#include(CIP-CSL4法,title)

***CIP-CSL2 [#t2bee67e]
#include(CIP-CSL2法,title)

***CIP-CSL3 [#sfe8a80b]
#include(CIP-CSL3法,title)

-例
#include(矩形波 - 設定,notitle)
#include(矩形波 - CIP-CSL4,notitle)
#include(矩形波 - CIP-CSL2,notitle)

体積保存を確かめるためにSin速度場で移流させた結果を示す.
#include(Sin速度場 - 設定,notitle)
#include(Sin速度場 - CIP,notitle)
#include(Sin速度場 - CIP-CSL4,notitle)
#include(Sin速度場 - CIP-CSL2,notitle)


**保存系セミラグランジュ法(Conservative Semi-Lagrangian scheme) [#d30ce6ad]
#include(保存系セミラグランジュ法,title)

-例
#include(矩形波 - 設定,notitle)
#include(矩形波 - 保存系セミラグランジュ法,notitle)

体積保存を確かめるためにSin速度場で移流させた結果を示す.
#include(Sin速度場 - 設定,notitle)
#include(Sin速度場 - セミラグランジュ法,notitle)
#include(Sin速度場 - 保存系セミラグランジュ法,notitle)



*時間微分の離散化 [#ze36fe32]
**修正オイラー法 [#u074bb60]
#include(修正オイラー法,title)
**前進オイラー法 [#u074bb60]
#include(前進オイラー法,title)

**ルンゲクッタ法(Runge-Kutta scheme) [#h0540b9d]
#include(RK,title)

**TVDルンゲクッタ(Total-Variation Diminishing Runge-Kutta) [#w5918b01]
#include(TVD RK,title)

**アダムス・バシュフォース法(Adams-Bashforth scheme) [#n6b9dee6]
#include(アダムス・バシュフォース法,title)

**完全陰解法(full implicit scheme) [#o1d773cf]
#include(完全陰解法,title)

**クランク・ニコルソン法(Crank-Nicolson scheme) [#b7d95388]
#include(クランク・ニコルソン法,title)

**予測子・修正子法(predictor-corrector scheme) [#o951b385]
#include(予測子・修正子法,title)

トップ   編集 差分 履歴 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS