有限体積法（FVM）｜制御体積で保存則を満たす離散化手法

有限体積法(FVM)

有限体積法（FVM）とは、支配方程式を体積要素（コントロールボリューム）にわたって積分し、面を通過する流束（フラックス）の収支として離散化する数値計算法である。保存則の形を保ったまま空間離散を行うため、質量・運動量・エネルギーなどの局所保存性に優れ、圧縮性・非圧縮性流れ、熱伝導、対流拡散現象などの計算流体力学（CFD）で広く用いられる。構造格子・非構造格子を問わず適用でき、複雑形状の離散化に強いこと、界面や衝撃波のような不連続の取り扱いに適していることが特長である。

Table Of Contents

有限体積法(FVM)

基本概念と保存則

有限体積法（FVM）の出発点は、支配方程式を保存形で表し、各セルにガウスの発散定理を適用して体積積分を面積分に変換する点にある。これにより、セル境界を通る数値流束の総和がセル内の物理量の変化率に一致する。セル中心配置（collocated）やスタガード配置を選択でき、圧力–速度の結合には補間（例: Rhie–Chow）が用いられる。

離散化の手順

ジオメトリと格子の準備（構造・非構造、直交・非直交、ボディフィット）
体積積分の面積分化とセル平均量の定義
対流・拡散・源項の評価と境界条件の付与
時間積分（陽解法・陰解法）と収束判定

数値流束の評価

対流項は上流差分的な考え方を取り入れ、Godunov系スキームやRoe、AUSM、HLLCなどで面を跨ぐ一次元リーマン問題を近似的に解く。二次精度化にはMUSCL再構築やTVD/TVB制限子（minmod、van Leerなど）を用い、傾き制限により非物理的振動を抑える。拡散項は面法線方向の勾配を中央差分で評価するのが一般的である。

時間積分と安定性

陽解法は各タイムステップが軽量で実装が容易である一方、CFL条件による制約を受ける。陰解法は大きなΔtが選べるが線形（または非線形）方程式系の解法を要し、前処理付き反復法などのソルバが鍵となる。定常問題では疑似時間進行や残差低減を指標に収束を図る。

圧力–速度連成（非圧縮流）

非圧縮性流れでは連続の式を満たす圧力場が必要になる。代表的な枠組みとしてSIMPLE、PISO、投影法があり、運動量式の予測解から圧力補正方程式（ポアソン方程式）を解いて発散のない速度場を得る。離散ポアソン問題は多重格子法やIC/ILU前処理CG/GMRESで効率化する。

境界条件の扱い

有限体積法（FVM）では、面に課す量が直接流束に影響するため、壁面の無滑り・等温、入口の速度・乱流量、出口の圧力、対称・周期などの境界条件の整合性が重要である。壁面近傍は壁関数または十分に細かい格子（y+管理）でせん断応力と熱流束を評価する。

格子とデータ配置

非構造格子（三角形・四面体・多面体）は形状追従性に優れ、局所細分で境界層や渦構造の解像度を高められる。構造格子（直交・曲線座標）は離散誤差が制御しやすく、メモリアクセス効率が高い。セル中心配置は汎用的で、スタガード配置は圧力–速度のデカップリングを避けやすい。

利点と欠点

利点：厳密な局所保存、格子の柔軟性、衝撃波や界面の表現に強い、離散の物理的解釈が明瞭
欠点：高次精度には再構築・制限子が必須で実装が複雑、非直交格子での勾配評価が難しい、陰解法では大規模連立解法がボトルネック

適用分野

圧縮性・非圧縮性CFD、熱流体の対流拡散、燃焼、気液二相流、微粒子輸送、放射伝熱などに広く用いられる。固体内熱伝導や電磁気の保存形方程式にも適用可能である。乱流モデル（RANS、LES、DES）や反応流モデルとも親和性が高い。

精度検証と格子収束

解の信頼性は格子収束試験で担保する。代表量の格子依存性を解析し、リチャードソン外挿やGCIで格子誤差を定量化する。スキームの形式精度と実効精度は異なることが多く、非直交やセル歪みが大きい場合は二次精度の維持に注意を要する。

実務での実装上の注意

ソルバ選択（CG/GMRES/BiCGStab）、前処理（ILU/AMG）、並列化（MPI/スレッド）の設計が解の時間を支配する。入出力のバイナリ化、疎行列フォーマット、キャッシュ効率を高めるデータ構造など、数値アルゴリズム外の最適化も重要である。ベンチマークはCFL、残差履歴、計算資源の実測で比較する。

代表的ソフトウェアと利用

オープンソースではOpenFOAMやSU2が著名で、教育・研究・産業で広く使われる。商用ではANSYS FluentやSTAR-CCM+などがあり、メッシュ生成から前後処理まで統合される。モデル選択・メッシュ独立性・検証問題のセットアップを標準作業とし、ブラックボックス化を避けることが肝要である。