有限体積法

有限体積法（FVM ）は、偏微分方程式を代数方程式の形で表現し、評価する方法である。 ^{[1]有限体積法では、}発散項 を含む偏微分方程式の体積積分は、発散定理を用いて面積分に変換される。これらの項は、各有限体積の表面におけるフラックスとして評価される。与えられた体積に入るフラックスは隣接する体積から出るフラックスと同一であるため、これらの方法は保存的である。有限体積法のもう一つの利点は、非構造メッシュを考慮に入れて容易に定式化できることである。この方法は、多くの数値流体力学パッケージで使用されている。「有限体積」とは、メッシュ上の各節点を取り囲む小さな体積を指す。^[2]

有限体積法は、節点値を用いて導関数を近似する有限差分法や、局所データを用いて解の局所近似値を作成し、それらをつなぎ合わせて全体近似値を構築する有限要素法と比較対照することができます。一方、有限体積法は、ある体積における解の平均値の正確な式を評価し、このデータを用いてセル内の解の近似値を構築します。^{[3] [4]}

例

単純な 1D移流問題を考えてみましょう。

${\displaystyle {\frac {\partial \rho }{\partial t}}+{\frac {\partial f}{\partial x}}=0,\quad t\geq 0.}$

1

ここで、は状態変数、は のフラックスまたは流れを表す。慣例的に、正は右方向の流れ、負は左方向の流れを表す。式( 1 )が一定面積の流れ媒体を表すと仮定すると、空間領域 を有限体積、つまりセル中心を でインデックス付けしたセルに分割することができる。特定のセル について、時刻および におけるの体積平均値を次のように 定義できる。${\displaystyle \rho =\rho \left(x,t\right)}$ ${\displaystyle f=f\left(\rho \left(x,t\right)\right)}$ ${\displaystyle \rho }$ ${\displaystyle f}$ ${\displaystyle f}$ ${\displaystyle x}$ ${\displaystyle i}$ ${\displaystyle i}$ ${\displaystyle {\rho }_{i}\left(t\right)=\rho \left(x,t\right)}$ ${\displaystyle {t=t_{1}}}$ ${\displaystyle {x\in \left[x_{i-1/2},x_{i+1/2}\right]}}$

${\displaystyle {\bar {\rho }}_{i}\left(t_{1}\right)={\frac {1}{x_{i+1/2}-x_{i-1/2}}}\int _{x_{i-1/2}}^{x_{i+1/2}}\rho \left(x,t_{1}\right)\,dx,}$

2

そして、時には、 ${\displaystyle t=t_{2}}$

${\displaystyle {\bar {\rho }}_{i}\left(t_{2}\right)={\frac {1}{x_{i+1/2}-x_{i-1/2}}}\int _{x_{i-1/2}}^{x_{i+1/2}}\rho \left(x,t_{2}\right)\,dx,}$

3

ここで、およびは、それぞれセルの上流​​および下流の面またはエッジの位置を表します。 ${\displaystyle x_{i-1/2}}$ ${\displaystyle x_{i+1/2}}$ ${\displaystyle i^{\text{th}}}$

式（１）を時間に関して積分すると次のようになる。

${\displaystyle \rho \left(x,t_{2}\right)=\rho \left(x,t_{1}\right)-\int _{t_{1}}^{t_{2}}f_{x}\left(x,t\right)\,dt,}$

4

どこ。 ${\displaystyle f_{x}={\frac {\partial f}{\partial x}}}$

時刻における の体積平均を得るには、セルの体積について積分し、その結果を で割ります。つまり、 ${\displaystyle \rho \left(x,t\right)}$ ${\displaystyle t=t_{2}}$ ${\displaystyle \rho \left(x,t_{2}\right)}$ ${\displaystyle \left[x_{i-1/2},x_{i+1/2}\right]}$ ${\displaystyle \Delta x_{i}=x_{i+1/2}-x_{i-1/2}}$

${\displaystyle {\bar {\rho }}_{i}\left(t_{2}\right)={\frac {1}{\Delta x_{i}}}\int _{x_{i-1/2}}^{x_{i+1/2}}\left\{\rho \left(x,t_{1}\right)-\int _{t_{1}}^{t_{2}}f_{x}\left(x,t\right)dt\right\}dx.}$

5

は適切に振る舞い、積分の順序を逆にできると仮定します。また、流れはセルの単位面積に垂直であることを思い出してください。ここで、1次元 では発散定理を適用できます。つまり、となり、 を有限体積のセル表面（端と）で評価したの値を発散の体積積分に代入すると、次のようになります。 ${\displaystyle f\ }$ ${\displaystyle f_{x}\triangleq \nabla \cdot f}$ ${\displaystyle \oint _{v}\nabla \cdot fdv=\oint _{S}f\,dS}$ ${\displaystyle f(x)}$ ${\displaystyle x_{i-1/2}}$ ${\displaystyle x_{i+1/2}}$

${\displaystyle {\bar {\rho }}_{i}\left(t_{2}\right)={\bar {\rho }}_{i}\left(t_{1}\right)-{\frac {1}{\Delta x_{i}}}\left(\int _{t_{1}}^{t_{2}}f_{i+1/2}dt-\int _{t_{1}}^{t_{2}}f_{i-1/2}dt\right).}$

6

どこ。 ${\displaystyle f_{i\pm 1/2}=f\left(x_{i\pm 1/2},t\right)}$

したがって、セル中心を、セル端のフラックスを として、( 6 )を時間について微分して以下の式を得ることにより、 上記の問題に対する半離散数値スキームを導くことができる。${\displaystyle i}$ ${\displaystyle i\pm 1/2}$

${\displaystyle {\frac {d{\bar {\rho }}_{i}}{dt}}+{\frac {1}{\Delta x_{i}}}\left[f_{i+1/2}-f_{i-1/2}\right]=0,}$

7

ここで、エッジフラックスの値は、セル平均の内挿または外挿によって再構成できます。式( 7 )は体積平均に対して正確です。つまり、導出の過程で近似は行われていません。 ${\displaystyle f_{i\pm 1/2}}$

この方法は、ノードの周囲の東西面に加えて、北面と南面を考慮することで 2D状況にも適用できます。

一般保存則

次のような偏微分方程式で表される一般保存則問題を考えることもできる。

${\displaystyle {\frac {\partial \mathbf {u} }{\partial t}}+\nabla \cdot {\mathbf {f} }\left({\mathbf {u} }\right)={\mathbf {0} }.}$

8

ここで、 は状態ベクトルを表し、は対応するフラックステンソルを表します。ここでも、空間領域を有限体積、つまりセルに分割することができます。特定のセル について、セル全体の体積 に対する体積積分 をとると、次の式が得られます。 ${\displaystyle \mathbf {u} }$ ${\displaystyle \mathbf {f} }$ ${\displaystyle i}$ ${\displaystyle v_{i}}$

${\displaystyle \int _{v_{i}}{\frac {\partial \mathbf {u} }{\partial t}}\,dv+\int _{v_{i}}\nabla \cdot {\mathbf {f} }\left({\mathbf {u} }\right)\,dv={\mathbf {0} }.}$

9

最初の項を積分して体積平均を求め、発散定理を2番目の項に適用すると、次の式が得られる。

${\displaystyle v_{i}{{d{\mathbf {\bar {u}} }_{i}} \over dt}+\oint _{S_{i}}{\mathbf {f} }\left({\mathbf {u} }\right)\cdot {\mathbf {n} }\ dS={\mathbf {0} },}$

10

ここで、は細胞の全表面積を表し、は表面に垂直で外側を向く単位ベクトルである。したがって、最終的に（ 8） と同等の一般的な結果を示すことができる。${\displaystyle S_{i}}$ ${\displaystyle {\mathbf {n} }}$

${\displaystyle {{d{\mathbf {\bar {u}} }_{i}} \over {dt}}+{{1} \over {v_{i}}}\oint _{S_{i}}{\mathbf {f} }\left({\mathbf {u} }\right)\cdot {\mathbf {n} }\ dS={\mathbf {0} }.}$

11

繰り返しになりますが、エッジフラックスの値は、セル平均値の内挿または外挿によって再構成できます。実際の数値計算手法は、問題の形状とメッシュ構成によって異なります。MUSCL再構成は、解に衝撃波や不連続性が存在する 高解像度の計算手法でよく使用されます。

有限体積法は、セル平均値がエッジフラックスを通じて変化するため、保守的です。言い換えれば、あるセルの損失は常に別のセルの利得となります。

参照

• 有限要素法
• フラックスリミッター
• ゴドゥノフの計画
• ゴドゥノフの定理
• 高解像度スキーム
• KIVA（ソフトウェア）
• MIT一般循環モデル
• MUSCLスキーム
• セルゲイ・K・ゴドゥノフ
• 総変動は減少している
• 非定常流れに対する有限体積法

参考文献

1. ルヴェック、ランドール（2002年）。双曲型問題のための有限体積法。ISBN 9780511791253。
2. Wanta, D.; Smolik, WT; Kryszyn, J.; Wróblewski, P.; Midura, M. (2021年10月). 「電気容量トモグラフィーのモデリングにおける四分木非均一構造メッシュを用いた有限体積法」. Proceedings of the National Academy of Sciences, India Section A: Physical Sciences . 92 (3): 443– 452. doi : 10.1007/s40010-021-00748-7 .
3. Fallah, NA; Bailey, C.; Cross, M.; Taylor, GA (2000-06-01). 「幾何学的非線形応力解析における有限要素法と有限体積法の適用の比較」.応用数学モデリング. 24 (7): 439– 455. doi : 10.1016/S0307-904X(99)00047-5 . ISSN  0307-904X.
4. Ranganayakulu, C. (Chennu) (2018年2月2日). 「第3章 3.1節」.コンパクト熱交換器：FEMおよびCFDアプローチを用いた解析、設計、最適化. Seetharamu, KN, Hoboken, NJ. ISBN 978-1-119-42435-2. OCLC  1006524487.

さらに読む

• Eymard, R. Gallouët, TR, Herbin, R. (2000)有限体積法 『数値解析ハンドブック』第7巻、2000年、p. 713–1020。編集者：PG CiarletおよびJL Lions。
• Hirsch, C. (1990)、「内部流れと外部流れの数値計算、第 2 巻: 非粘性流れと粘性流れの計算方法」、Wiley。
• Laney, Culbert B. (1998)、「計算ガスダイナミクス」、ケンブリッジ大学出版局。
• LeVeque, Randall (1990)、「保存則の数値的手法」、ETH 数学講義シリーズ、Birkhauser-Verlag。
• LeVeque, Randall (2002)、「双曲型問題のための有限体積法」、ケンブリッジ大学出版局。
• Patankar, Suhas V. (1980)、「数値的熱伝達と流体の流れ」、Hemisphere。
• Tannehill, John C.、他 (1997)、「計算流体力学および熱伝達」、第 2 版、Taylor and Francis。
• Toro, EF (1999)、「Riemann Solvers and Numerical Methods for Fluid Dynamics」、Springer-Verlag。
• Wesseling, Pieter (2001)、「計算流体力学の原理」、Springer-Verlag。

外部リンク

• R. Eymard、T Gallouët、 R. Herbinによる有限体積法、2000年の数値解析ハンドブックに掲載された記事の更新版
• Rübenkönig, Oliver. 「有限体積法（FVM）入門」。2009年10月2日時点のオリジナルよりアーカイブ。GFDLに基づいて利用可能です。
• FiPy: NIST の Python を使用した有限体積 PDE ソルバー。
• CLAWPACK: 波動伝播アプローチを用いて双曲型偏微分方程式の数値解を計算するために設計されたソフトウェアパッケージ