有限体積法のソースを表示

{{Differential equations}}
{{計算物理学}}
'''有限体積法'''（ゆうげんたいせきほう、{{lang-en|finite volume method}}、FVM）とは、[[数値解析]]手法の一つである。領域を有限個のコントロールボリューム（{{en|control volume}}）に分割し、各ボリュームに対して積分形の物理量の[[保存則|保存方程式]]を適用するものである<ref name=Ferziger>{{cite|和書 |title=コンピュータによる流体力学 |author=Joel H. Ferziger |author2=Milovan Peri&#x107; |translator=小林敏雄、谷口伸行、坪倉誠 |publisher=シュプリンガー・フェアラーク東京 |year=2003 |isbn=4-431-70842-1 |pages=36-37}}</ref><ref name="intro"/><ref name="cfd"/>。

[[1960年代]]に[[ロスアラモス国立研究所]]において{{仮リンク|非構造格子|en|unstructured grid}}に基づく[[数値流体力学|流体解析]]手法として開発され<ref name="intro"/><ref name="cfd"/>、現在では、多くの商用の流体解析コードに標準的な離散化解析手法として採用されている<ref name="intro">Versteeg, H. K., & Malalasekera, W. (2007). An introduction to computational fluid dynamics: the finite volume method. Pearson education.</ref><ref name="cfd">Moukalled, F., Mangani, L., & Darwish, M. (2016). The finite volume method in computational fluid dynamics (Vol. 113). Berlin, Germany:: Springer.</ref><ref name=JSCE>{{cite|和書 |editor=[[土木学会]] 応用力学委員会 計算力学小委員会 |title=いまさら聞けない計算力学の常識 |publisher=丸善 |year=2008 |isbn=978-4-621-08022-1 |page=9-10}}</ref>。

== 概要 ==
[[有限差分法]]と[[有限要素法]]の両方の特徴を合わせ持つ手法と言える<ref name=JSCE/><ref>Idelsohn, S. R., & Onate, E. (1994). Finite volumes and finite elements: two ‘good friends’. International journal for numerical methods in engineering, 37(19), 3323-3341.</ref>。

解析領域をセル（{{en|cell}}）と呼ばれる小領域に分割し、セルの格子点を中心とする領域であるコントロールボリュームあるいは検査領域''D<sub>e</sub>'' を定義する。そして、有限要素法と同様にその離散化には[[重み付き残差法]]<ref>Finlayson, B. A., & Scriven, L. E. (1966). The method of weighted residuals—a review. Appl. Mech. Rev, 19(9), 735-748.</ref>を適用する。ただし有限体積法では、コントロールボリューム''D<sub>e</sub>'' ごとに、重み関数を 1 として重み付き残差式を離散化する。
:<math>\int_{D_{e}} u^{*} f(u)d\Omega =\int_{D_{e}} f(u)d\Omega = 0</math>

;長所
*保存方程式をコントロールボリュームで積分するので、この積分領域内の物理量の[[保存則|保存]]が満足される。コントロールボリュームが重ならないかぎり、領域全体での保存性も満足される。
*非構造格子など、どのようなタイプの[[計算格子]]にも適用できるため、任意形状への適合性が良い<ref name=Ferziger/>。
*上式における積分中の微分の近似には中点公式（または中点差分近似<ref name="yama"> 山本哲朗『数値解析入門』サイエンス社〈サイエンスライブラリ 現代数学への入門 14〉、2003年6月、増訂版</ref>）
:<math>\frac{df}{dx}\simeq\frac{f(x+h)-f(x-h)}{2h}\quad (h<<1)</math>
:が一般に用いられるため、離散化が簡便で有限要素法に比べて計算時間の面で有利である<ref name="yama"/>。なお、微分の近似表現に中点公式などを用いているため、構造格子を用いた場合には、離散化された代数方程式は有限差分法を適用して導かれたそれと一致することがある<ref name=JSCE/>。

;短所
*高次精度化が煩雑あるいは困難である<ref name=Ferziger/><ref name="intro"/><ref name="cfd"/>。有限体積法は3段階の近似（補間、微分、積分）を必要とするために、3次元で2次精度よりも高い精度の方法を実現することが難しい<ref name=Ferziger/><ref name="intro"/><ref name="cfd"/>。
==参考文献==
* Eymard, R. Gallouët, T. R. Herbin, R. (2000) The finite volume method, in Handbook of Numerical Analysis, Vol. VII, 2000, p. 713–1020. Editors: P.G. Ciarlet and J.L. Lions.
* LeVeque, Randall (2002), Finite Volume Methods for Hyperbolic Problems, [[:en:Cambridge University Press]].
== 脚注 ==
{{reflist|2}}
{{偏微分方程式の数値解法}}

{{Applied-math-stub}}
{{Physics-stub}}

{{DEFAULTSORT:ゆうけんたいせきほう}}
[[Category:物理数学]]
[[Category:力学]]
[[Category:数値解析]]
[[Category:数学に関する記事]]