有限要素法(FEM)と有限体積法(FVM)はどちらも偏微分方程式を数値的に解くための代表的な離散化手法ですが、数学的定式化・適用分野・計算特性に明確な違いがあります。
構造解析では有限要素法が主流、流体解析(CFD)では有限体積法が主流というのが現在の工業界の大まかな使い分けですが、両手法の理論的背景と特性を理解することで、問題に応じた適切な選択が可能になります。
この記事では、有限要素法と有限体積法の定式化の違い・保存則の扱い・適用分野・代表的ソフトウェアについて詳しく解説していきます。
有限要素法と有限体積法の基本的な違いをまず理解しよう
それではまず、有限要素法と有限体積法それぞれの基本的な定式化と根本的な違いについて解説していきます。
最も本質的な違いは有限要素法が弱形式(変分形式)への変換を基礎とする手法であり、有限体積法が保存則の積分形(検査体積に対する収支バランス)を基礎とする手法であるという点です。
有限要素法では解析領域を要素に分割し、各要素内で形状関数による補間を行って全体の剛性方程式を組み立てます。
有限体積法では解析領域を制御体積(コントロールボリューム)に分割し、各制御体積の境界を通過するフラックス(流束)の収支バランスを離散化方程式として定式化します。
有限体積法の基本原理
有限体積法(FVM:Finite Volume Method)は保存型の偏微分方程式(質量保存・運動量保存・エネルギー保存)を制御体積に対する積分形として表現し、フラックスの離散化によって数値解を求める手法です。
保存量(質量・運動量・エネルギー)の制御体積内での変化が、境界を通過するフラックスの正味収支に等しいという積分保存則が離散化の出発点です。
有限体積法の基本方程式(1次元定常移流拡散の例):d/dx(ρuφ) = d/dx(Γ dφ/dx) + S。これを制御体積[w,e]で積分すると:(ρuφA)_e – (ρuφA)_w = (ΓA dφ/dx)_e – (ΓA dφ/dx)_w + SΔV。各境界面での値を補間スキームで離散化します。
有限体積法は各制御体積で物理量の保存則が厳密に満足されるため、保存性が重要な流体解析・熱移動解析に特に適した手法です。
有限差分法との三者比較
有限差分法(FDM)・有限要素法(FEM)・有限体積法(FVM)の三手法の主な特性を比較します。
| 特性 | 有限差分法(FDM) | 有限要素法(FEM) | 有限体積法(FVM) |
|---|---|---|---|
| 離散化の基礎 | テイラー展開による微分近似 | 弱形式・変分法 | 積分保存則 |
| メッシュの柔軟性 | 低い(構造格子) | 高い(非構造格子可) | 高い(非構造格子可) |
| 複雑形状への対応 | 難しい | 優れている | 良好 |
| 保存性 | 局所的に非保存の場合あり | 大域的保存性 | 局所的・大域的保存性 |
| 主な適用分野 | 流体・熱・地震波伝播 | 構造・電磁気・熱 | 流体(CFD)・熱・多相流 |
| 代表ソフトウェア | FDTD法ツール | ANSYS・ABAQUS | OpenFOAM・FLUENT |
適用分野による使い分けの実際
工業界での使い分けとして、構造強度解析・固体力学・線形・非線形変形解析では有限要素法が圧倒的に主流です。
流体力学(CFD:Computational Fluid Dynamics)・燃焼解析・多相流解析では有限体積法が主流で、OpenFOAM・FLUENT・StarCCM+などが代表的なFVMソルバーです。
電磁気解析では有限要素法(ANSYS Maxwell・COMSOL)と有限差分時間領域法(FDTD)が分野によって使い分けられています。
熱流体連成解析(FSI:流体構造連成)では、流体側にFVM・構造側にFEMを組み合わせたマルチフィジクスアプローチが採用されます。
有限要素法と有限体積法の数学的定式化の違い
続いては、両手法の数学的定式化の詳細な違いを確認していきます。
定式化の違いを理解することで、各手法の強みと限界がより明確になります。
弱形式と積分保存則の違い
有限要素法の出発点は支配方程式の強形式(微分形式)を弱形式(変分形式)に変換することです。
弱形式では解の連続性要件が緩和され(C⁰連続性で十分)、複雑な境界条件・材料不均質性の取り扱いが容易になります。
有限体積法は積分保存則を出発点とするため、物理量の局所的な保存が各制御体積で厳密に保証されます。
この局所保存性はシミュレーション全体での質量・運動量・エネルギーの過不足が生じないことを保証する重要な特性です。
有限体積法の局所保存性は、非定常流体解析・衝撃波・接触不連続面を含む問題での数値的信頼性において特に重要な特性です。
補間スキームと数値拡散の問題
有限体積法では制御体積境界でのフラックスを計算する際に補間スキームの選択が精度と安定性に大きく影響します。
一次精度の迎え風差分スキーム(Upwind Scheme)は安定性が高いですが数値拡散(過剰な数値的散逸)が大きい問題があります。
二次精度の中心差分スキーム(Central Differencing Scheme)は精度が高いですが高ペクレ数(対流支配)で振動(数値不安定)が生じる場合があります。
TVD(Total Variation Diminishing)スキームやHRスキーム(High Resolution Scheme)は安定性と精度を両立させた高性能スキームです。
非構造メッシュへの対応
両手法とも非構造メッシュ(三角形・四面体など任意形状の要素)に対応していますが、アプローチが異なります。
有限要素法は元々任意形状要素への対応が基本設計に組み込まれており、アイソパラメトリック要素によって複雑形状を高精度に表現できます。
有限体積法も非構造メッシュへの拡張が進んでいますが、非直交メッシュでの拡散項の離散化に非直交補正が必要で実装が複雑になります。
どちらの手法も現代のソフトウェアではGUI上での自動メッシュ生成・品質評価ツールが整備されており、実用上の障壁は大幅に低下しています。
代表的なFEM・FVMソフトウェアと選び方
続いては、代表的な有限要素法・有限体積法ソフトウェアの特徴と選び方を確認していきます。
適切なソフトウェア選択が解析効率・品質に大きく影響します。
主要FEMソフトウェア
ANSYS Mechanicalは世界最大シェアの汎用構造解析ソフトウェアで、線形・非線形・動的・疲労解析に幅広く対応しています。
ABAQUSは特に材料非線形・接触・大変形解析に強く、自動車・航空宇宙・製造業での高度解析に多用されています。
NX Nastranは航空宇宙業界で標準的に使用される高信頼性FEMソルバーです。
オープンソースのCode_Aster(フランス電力EDF開発)・CalculiXは無料で使用可能な高機能FEMソルバーとして研究・中小企業に普及しています。
FEMソフトウェアの選定では解析の種類・必要精度・ライセンスコスト・サポート体制・社内の習熟度を総合的に考慮することが重要です。
主要FVMソフトウェア(CFDツール)
OpenFOAMはオープンソースの汎用CFDソルバーで、研究・教育・中小企業に広く普及しており、ユーザーコミュニティが活発です。
ANSYS FLUENTは世界最大シェアの商用CFDソフトウェアで、豊富な物理モデル・使いやすいGUIが特徴です。
StarCCM+は自動メッシュ生成・連成解析機能に優れた商用CFDソフトウェアで、自動車・船舶・エネルギー分野で多用されています。
マルチフィジクス解析ツール
COMSOLはFEM・FVM・FDMを統合したマルチフィジクス解析プラットフォームで、流体・構造・電磁気・化学反応の連成解析を一つの環境で実施できます。
ANSYS WorkbenchはFEM(Mechanical)とFVM(Fluent)を同一プラットフォームで統合管理でき、流体構造連成(FSI)解析に強みを持ちます。
マルチフィジクス解析への需要は今後も増大することが見込まれており、FEM・FVMを統合したプラットフォームの重要性が高まっています。
まとめ
この記事では、有限要素法と有限体積法の基本的な定義・定式化の違い・保存性・適用分野・代表的ソフトウェアについて解説しました。
有限要素法は弱形式への変換を基礎とし複雑形状・構造解析に強み、有限体積法は保存則の積分形を基礎とし流体解析での局所保存性に優れるという根本的な特性の違いがあります。
実際の工業解析では問題の種類(構造・流体・熱・電磁気)に応じた使い分けと、必要に応じた連成解析アプローチが最適解を導きます。
両手法を理解した上で適切なツールを選択・活用することが、高品質な数値解析を行う上での重要な技術リテラシーといえるでしょう。