FEniCS/FEniCSx — Python×FEMの研究プラットフォーム徹底解説

FEniCSとは何か

例えばポアソン方程式 \(-\nabla^2 u = f\) の弱形式は：

$$\int_\Omega \nabla u \cdot \nabla v \, dx = \int_\Omega f v \, dx$$

これをFEniCSでは a = inner(grad(u), grad(v)) * dx と書く。数学の教科書とコードがほぼ一対一対応するのが最大の特徴だ。

UFL — 数学をそのままコードにする

• grad(u) — 勾配 \(\nabla u\)
• div(u) — 発散 \(\nabla \cdot u\)
• curl(u) — 回転 \(\nabla \times u\)
• inner(a, b) — 内積 \(a \cdot b\)
• dx — 体積積分の測度
• ds — 境界積分の測度

例えば線形弾性体の変分形式なら、応力テンソル \(\sigma = \lambda \mathrm{tr}(\varepsilon) I + 2\mu \varepsilon\) をそのままUFLで書ける。商用ソルバーでは「線形弾性モデルを選択」するだけだが、FEniCSでは構成則そのものをカスタマイズできるのが強みだ。

自動コード生成の仕組み

ユーザーから見ると「Pythonで3行書いただけ」なのに、裏では何千行ものC++コードが自動生成されて高速実行されている。この「高レベル記述 → 低レベル最適化」の自動変換がFEniCSのアーキテクチャの核心だよ。

FEniCSxへの進化

• 並列性能：MPIベースの分散メモリ並列がより洗練され、数千コアへのスケーリングが改善
• メッシュ対応：高次要素メッシュのネイティブサポート、混合要素の取り扱い改善
• 外部メッシャー連携：GmshやMeshio経由でのメッシュインポートが容易に
• APIの見直し：より明示的でPythonic な設計。古い暗黙の挙動が排除された

新規プロジェクトではFEniCSxを使うべきだ。旧版のチュートリアルが多いので混乱しやすいが、公式ドキュメントはFEniCSxに統一されつつある。

マルチフィジックス研究での活用

• 生体力学：心臓の電気伝導シミュレーション、脳組織の力学モデル。複雑な異方性材料の構成則をUFLで直接実装
• 相場モデル（Phase-field）：Cahn-Hilliard方程式やAllen-Cahn方程式による破壊伝播・相分離のモデリング
• 流体-構造連成（FSI）：ALE定式化やイマーズド法による流体と構造の相互作用
• トポロジー最適化：密度法やレベルセット法による形状最適化の研究実装

Google Scholarで "FEniCS" を検索すると、年間1,000件以上の論文がヒットする。研究用FEMプラットフォームとしては世界最大級の利用実績だ。

商用ソルバー・他OSSとの比較

• Abaqus/ANSYS：「すぐに使える」完成品。GUIで設定して計算する。カスタマイズ性は限定的
• FEniCS：「研究者が弱形式から自由に組める」フレームワーク。GUIなし。Pythonプログラミング必須
• deal.II：C++ベースのFEMライブラリ。適応メッシュ細分化（AMR）が最強。FEniCSより低レベルだが、性能チューニングの自由度が高い
• CalculiX：Abaqus互換の.inpファイルで動く汎用ソルバー。FEniCSとは設計思想が全く異なる

FEniCSは「アイデアをすぐコードにしたい研究者」に最適化されたツールだ。産業用途で定型的な解析を回すなら商用ソルバー、数値手法の研究開発ならFEniCS、という棲み分けになる。