1 概説
計算流体力学は、流体のふるまいを数値計算によって扱う学問分野である。流れ場、圧力、温度、濃度などを方程式として表し、計算機上で近似的に解くことで、観測や実験だけでは捉えにくい現象を解析する。
この方法は、流体そのものの観察に加えて、設計条件の比較や性能予測にも用いられる。航空宇宙、自動車、エネルギー、気象、海洋、医療機器など、対象は多岐にわたる。
1.1 定義
計算流体力学は、流体の支配方程式を離散化し、数値解を求める技術と理論の総称である。対象は気体、液体、混相系を含み、熱や物質の移動を伴う場合も扱う。
解析の中心には、連続体としての流体を仮定し、空間と時間を細かい要素に分けて計算する手続きがある。これにより、複雑な形状や条件下の流れを再現しやすくなる。
1.2 研究対象
研究対象は、速度や圧力の分布だけでなく、温度差による対流、拡散、界面の移動、反応に伴う変化なども含む。現象によっては、固体の変形や電磁場との相互作用を組み合わせて扱うこともある。
1.2.1 流体の運動
流体の運動は、速度場、渦、せん断、圧力勾配などによって特徴づけられる。これらは物体まわりの抗力や揚力、配管内の圧力損失、噴流の拡がりなどに関係する。
1.2.2 熱移動と物質移動
熱移動と物質移動では、温度分布や濃度分布が時間とともに変化する。対流、伝導、拡散の組み合わせによって、冷却、混合、蒸発、吸収などの挙動が決まる。
1.3 歴史的背景
計算流体力学は、流体力学と数値解析、計算機科学の発展とともに形成された。初期には単純な近似や限定的な格子計算が中心だったが、計算機性能の向上により、より複雑な流れの解析が可能になった。
その後、有限体積法や高精度スキーム、乱流モデル、多相流モデルなどが整備され、工学設計の実用手段として広く普及した。
2 基礎理論
計算流体力学の基盤は、流体を連続体として記述する力学である。そこでは、質量、運動量、エネルギーの保存を満たす方程式が中心となる。
2.1 連続体力学
連続体力学では、流体を分子の集合ではなく、密度や速度、圧力を連続的に定義できる媒質として扱う。この仮定により、微分方程式を用いた記述が可能になる。
2.1.1 質量保存則
質量保存則は、流体の質量が生成も消滅もしないという原理である。流入と流出の差、あるいは密度変化によって、ある領域内の質量収支が決まる。
2.1.2 運動量保存則
運動量保存則は、外力と応力が流体の速度変化を生むことを表す。圧力、粘性、重力などが主要な要因であり、流れの加速や減速、方向転換に関与する。
2.1.3 エネルギー保存則
エネルギー保存則は、内部エネルギー、運動エネルギー、熱の移動を統一的に記述する。温度変化や圧縮、粘性散逸を考慮することで、熱流動の解析が可能になる。
2.2 支配方程式
支配方程式は、流体の状態を定める中心的な数式である。問題設定に応じて、非粘性近似、粘性流体、乱流、反応を含む形へ拡張される。
2.2.1 オイラー方程式
オイラー方程式は、粘性を無視した流れを表す。高速流れや、粘性の影響が相対的に小さい場面で用いられ、衝撃波を伴う問題にも現れる。
2.2.2 ナビエ・ストークス方程式
ナビエ・ストークス方程式は、粘性流体の運動を記述する基本式である。実用上の多くの流体現象はこの式を基礎に扱われ、数値計算の中心的対象となっている。
2.2.3 乱流方程式
乱流方程式は、乱れた速度変動の影響を平均的に表現するための式群である。直接解く方法のほか、平均化やモデル化を通じて、複雑な流れの性質を近似する。
2.3 境界条件と初期条件
境界条件は、壁面、入口、出口、対称面などで流体が満たすべき制約を与える。初期条件は、計算開始時点の速度、圧力、温度などの状態を定める。
これらの設定は解の妥当性に大きく影響する。条件が不適切だと、数値解は実際の現象からずれやすくなる。
3 数値手法
数値手法は、連続な方程式を計算可能な形に置き換える技術である。離散化、時間更新、格子設計、安定化の工夫が組み合わされて構成される。
3.1 離散化の考え方
離散化は、空間や時間を有限個の点や領域に分け、微分を差分や積分の近似に変換する操作である。これにより、方程式は代数方程式の集合として扱えるようになる。
3.1.1 有限差分法
有限差分法は、導関数を格子点間の差で近似する方法である。実装が比較的簡潔で、規則的な格子に適している。
3.1.2 有限体積法
有限体積法は、各セルにおける保存量の収支を基礎にする。保存則を自然に満たしやすく、工学分野で広く用いられている。
3.1.3 有限要素法
有限要素法は、解を要素ごとの基底関数で表現する手法である。複雑な形状への適応力が高く、流体と構造の連成にも利用される。
3.2 時間積分法
時間積分法は、離散化された方程式を時間方向に進める手続きである。精度、安定性、計算負荷のバランスが重要になる。
3.2.1 陽解法
陽解法は、現時点の値から次の時刻の値を明示的に計算する方法である。扱いやすい一方、時間刻みに制約が生じやすい。
3.2.2 陰解法
陰解法は、未知の次時刻の値を含む式を解く方法である。安定性に優れる場合が多いが、各ステップで反復計算が必要になることがある。
3.3 格子生成
格子生成は、計算領域を要素分割する作業である。格子の品質は、解の精度と計算効率に直接影響する。
3.3.1 構造格子
構造格子は、規則的な番号付けを持つ格子である。計算手順が整理しやすく、整った領域では高効率に使える。
3.3.2 非構造格子
非構造格子は、三角形や四面体など不規則な要素で構成される。複雑な幾何形状に対応しやすい点が利点である。
3.3.3 適合格子と適応格子
適合格子は境界や領域の形に合わせて整えられた格子であり、適応格子は解の変化が大きい部分を自動的に細分化する。後者は局所的な精度向上に有効である。
3.4 安定性と収束性
安定性は計算誤差が増幅しにくい性質を指し、収束性は格子や刻みを細かくしたときに解が真の値へ近づく性質を意味する。両者は実用計算の信頼性を左右する。
3.4.1 数値拡散
数値拡散は、離散化により実際より滑らかな分布が生じる現象である。急な勾配や界面がぼやける原因となる。
3.4.2 数値振動
数値振動は、解の近傍に不自然な揺れが現れる現象である。高勾配領域や不適切なスキームで生じやすい。
3.4.3 収束判定
収束判定は、計算をいつ終了するかを決める基準である。残差、保存量の変化、観測量の安定化などが目安として使われる。
4 解析対象と応用
計算流体力学は、単純な流れから複雑な相変化や化学反応まで幅広く対応する。対象に応じてモデルや数値手法が選ばれる。
4.1 単相流
単相流は、ひとつの相のみからなる流れである。最も基本的な対象であり、多くの理論や計算法の出発点となる。
4.1.1 層流
層流は、流線が比較的整い、流体層が滑らかに移動する状態である。低速、低レイノルズ数の条件で現れやすい。
4.1.2 乱流
乱流は、渦や速度変動が複雑に入り混じる流れである。工学上の多くの場面で重要だが、再現にはモデル化が必要になることが多い。
4.2 多相流
多相流は、気体、液体、固体など複数の相が同時に存在する流れである。界面の追跡や相互作用の記述が大きな課題となる。
4.2.1 気液流
気液流は、泡、噴霧、気泡流、液滴流などを含む。相の分布や界面の移動が、圧力損失や熱交換に影響する。
4.2.2 固液流
固液流は、粒子が液体中を移動する現象である。土砂移動、スラリー、沈降などの解析に関係する。
4.2.3 気固流
気固流は、気体中の粒子輸送を扱う。集じん、粉体搬送、噴霧乾燥などで重要になる。
4.3 反応流
反応流は、流れと化学変化が同時に進む系である。温度上昇、組成変化、発熱反応などが流動場に影響を与える。
4.3.1 燃焼
燃焼は、燃料と酸化剤の反応によって熱と生成物が生じる過程である。火炎構造、混合、反応速度の評価が中心となる。
4.3.2 化学反応
化学反応は、複数成分の濃度変化とエネルギー収支を伴う。反応機構が複雑な場合、簡略化されたモデルが用いられることがある。
4.4 連成解析
連成解析は、流体の挙動を他の物理現象と結びつけて扱う方法である。相互作用を含めることで、実際の装置や構造の挙動に近づけられる。
4.4.1 流体構造連成
流体構造連成は、流れと固体の変形が相互に影響し合う解析である。柔軟な構造物や振動を伴う機器で重要となる。
4.4.2 流体熱連成
流体熱連成は、流れと温度場を同時に解く考え方である。冷却、加熱、熱交換器の設計に広く使われる。
5 実務上の課題
実務での利用では、理論式をそのまま解くだけでは不十分である。モデル選択、計算時間、精度保証の各要素を総合的に考える必要がある。
5.1 モデル化誤差
モデル化誤差は、現象を単純化したことによって生じるずれである。未解明の物理や省略した要因が、結果の差として現れる。
5.1.1 乱流モデルの選択
乱流モデルの選択は、計算精度と負荷の妥協点を決める作業である。対象の流れに応じて、平均化モデルやより詳細な手法が使い分けられる。
5.1.2 近似の限界
近似の限界は、計算モデルが再現できる範囲の上限を意味する。極端に複雑な流れや微細構造では、誤差が増えやすい。
5.2 計算資源
計算資源は、処理時間、メモリ、計算機台数などを含む。大規模問題では、資源配分が解析の成否を左右する。
5.2.1 並列計算
並列計算は、複数の計算装置に処理を分担させる方法である。大きな領域や高解像度計算で特に有効である。
5.2.2 高性能計算
高性能計算は、専用の大規模計算機や高速演算環境を利用する。詳細な三次元解析や長時間計算を支える基盤となる。
5.3 検証と妥当性確認
検証と妥当性確認は、計算結果が正しく実装されているか、また現実を適切に表しているかを確かめる工程である。信頼できる解析には不可欠である。
5.3.1 実験との比較
実験との比較は、数値結果を観測データと照らし合わせる作業である。流速、圧力、温度などの一致度からモデルの妥当性を評価する。
5.3.2 ベンチマーク問題
ベンチマーク問題は、標準的な試験例として広く参照される。異なる手法の性能比較や、コードの品質確認に用いられる。
6 関連分野
計算流体力学は、他の流体研究や計算科学と密接に結びついている。目的に応じて、実験、理論、数値、設計の各分野と連携する。
6.1 実験流体力学
実験流体力学は、実測によって流れを調べる分野である。可視化や計測を通じて、数値解析の検証材料を提供する。
6.2 理論流体力学
理論流体力学は、流体現象を数式で理解しようとする分野である。支配方程式の性質や解の構造を明らかにする点に重点がある。
6.3 数値解析
数値解析は、連続的な問題を計算可能な形へ変換する理論と技術の総称である。誤差評価、収束、安定性の研究は、この分野の重要な柱である。
6.4 工学シミュレーション
工学シミュレーションは、実機試験の前段階で性能や挙動を見積もる方法である。流体解析は、その中核技術の一つとして位置づけられる。
</INTERNAL_LINK_CANDIDATES> 連続体力学（流体を連続的な媒質として扱う理論） 保存則（質量や運動量などが保たれる原理） オイラー方程式（粘性を無視した流れの支配方程式） ナビエ・ストークス方程式（粘性流体の運動を表す基本方程式） 乱流（速度や渦が不規則に変動する流れ） 有限差分法（導関数を格子点の差で近似する手法） 有限体積法（各領域の収支で保存量を計算する手法） 有限要素法（要素ごとの基底関数で解を表す手法） 陽解法（現時点の値から次時刻を直接求める時間積分法） 陰解法（未知の次時刻を含む方程式を解く時間積分法） 構造格子（規則的な番号付けを持つ格子） 非構造格子（不規則な要素で構成される格子） 適応格子（解の変化に応じて細分化する格子） 数値拡散（離散化で生じる解の平滑化） 数値振動（離散化で生じる不自然な揺れ） ベンチマーク問題（手法比較や品質確認に用いる標準例） 流体構造連成（流れと固体変形が相互作用する解析） 流体熱連成（流れと温度場を同時に解く解析） 実験流体力学（実測で流れを調べる分野） 高性能計算（大規模な計算を支える計算環境）