ALE法によるFSI

Eulerian記述では格子が固定されていて流体が格子を通過する。大変形に強いけど界面追跡が苦手。Lagrangian記述では格子が物質と一緒に動くから界面は鮮明だけど、大変形でメッシュが潰れる。ALE法はこの中間で、格子速度を任意（Arbitrary）に設定できる。界面では格子を物質と一緒に動かし（Lagrangian的）、内部では格子を適度に動かしてメッシュ品質を維持する（Eulerian的）。

ALE記述でのNavier-Stokes方程式の運動量保存則はこうなる。

その通り。FSI問題では流体-構造界面上で $\hat{\mathbf{u}} = \dot{\mathbf{d}}_s$（構造の速度）とし、内部のメッシュ速度はスムージングアルゴリズムで決定する。

最も重要なのがGCL（幾何学的保存則）の満足だ。メッシュが動く場合、保存則の離散化が一貫していないと、一様流に対して偽のソース項が発生する。

一様流を入力しても圧力やエネルギーに非物理的な変動が生じる。特に時間精度が1次だとGCL誤差が顕著になる。2次精度時間積分を使い、メッシュ速度の計算に前ステップのメッシュ位置を正しく使うことが重要だよ。

代表的な手法を比較しよう。

連成の強さによって安定性と精度が大きく変わる。

弱連成は1ステップの遅れ（time lag）がある。密度比 $\rho_s/\rho_f$ が大きい場合（金属構造+空気など）は弱連成でも安定だが、$\rho_s/\rho_f \approx 1$（血管+血液、ゴム+水など）では付加質量不安定性が発生して弱連成は使えない。

空気中の構造（$\rho_s/\rho_f \gg 1$）なら弱連成でも多くの場合安定だ。ただし精度を重視するなら強連成を使う。水中構造（$\rho_s/\rho_f \sim 1-10$）は必ず強連成が必要だよ。

Aitken $\Delta^2$ 加速法が最もシンプルかつ効果的だ。連成反復の緩和係数を動的に最適化する。

preCICEには標準搭載。Ansys System Couplingにも最近追加された。OpenFOAMではsolidDisplacementFoamとpimpleFoamの連成でカスタム実装が必要だけど、preCICEアダプターを使えば簡単だよ。

構造の変形量が元のメッシュサイズと同程度になると、スムージングだけでは対応できない。そのときは自動リメッシング（re-meshing）が必要だ。

Turek & Hron (2006)のFSIベンチマークが最も広く使われている。円柱の後ろに弾性フラグを取り付けた2次元チャネル流れの問題だ。

その通り。密度比1は付加質量効果が最大になるケースで、弱連成では確実に発散する。強連成 + Aitken加速（またはIQN-ILS）が必須だ。このベンチマークでの参照解はフラグ先端のx方向変位・y方向変位と抗力・揚力の時刻歴で、複数の研究グループが高精度の参照値を公開している。

Ansys Workbenchでのワークフローはこうだ。

構造の固有振動周期 $T_n$ を基準にして、$\Delta t < T_n / 20$ とするのが安全だ。流体側のCFL条件も同時に確認する。FSI問題では構造の時間スケールが流体より短いことが多いから、構造側の要求で時間刻みが決まることが多い。

OpenFOAM + preCICE + CalculiX（またはdeal.II）の組み合わせが定番だ。preCICEのチュートリアルにTurekベンチマークが含まれているので、まずはそれを動かしてみるのがいいよ。

preCICEの設定ファイル（precice-config.xml）で重要なパラメータは以下だ。

以下にまとめるよ。

COMSOLのFSIモジュールは設定が非常に簡単で、教育用途や小規模問題には最適だ。ただしメッシュ規模が大きくなると計算効率でFluent/STAR-CCM+に劣る。また、乱流モデルの種類が限られるのが産業用途での弱点だ。

以下の判断基準で選択するといいよ。

変形量が翼弦長の10%以下なら問題なくALE法で対応できる。それを超える場合（例えば大振幅のLCO）や、制御舵面の大角度回転を含む場合はOverset meshの方が安全だ。Fluent 2022以降はOverset meshの精度が大幅に改善されていて、FSIとの組み合わせも安定しているよ。

いくつかの代表的な事例だ。

血管壁は超弾性体（Mooney-Rivlin等）でモデル化し、血液は非ニュートン流体（Carreau-Yasudaモデル等）で扱う。密度比が1に近いから強連成が必須で、計算は難しいけど臨床応用への期待が大きい分野だよ。

Tezduyarらが提唱したDSD/SST（Deforming-Spatial-Domain/Stabilized Space-Time）法が有名だ。時間と空間を統合した4次元時空有限要素法で、メッシュの変形を自然に取り込める。GCLも自動的に満足される美しい定式化だよ。

確かに実装の複雑さと計算コストが障壁で、商用ソフトには搭載されていない。研究コードレベルではTeam for Advanced Flow Simulation and Modeling（TAFSM、Tezduyarグループ）が精力的に開発している。タイヤと地面の接触、パラシュートの展開など、大変形FSIの高精度解析に使われているよ。

CutFEM（Cut Finite Element Method）やXFEM（eXtended FEM）を流体側に使う手法が研究されている。固定のEulerianメッシュ上で流体を解き、構造界面でメッシュを「切断」して境界条件を適用する。

利点は以下の通りだ。

もちろん。FSI + 熱連成（Conjugate Heat Transfer + FSI）は、以下のような応用で重要だ。

Ansys System Couplingでは3ソルバー同時連成（Fluent + Mechanical + ICEPACK等）が可能だ。STAR-CCM+のCo-SimulationもAbaqusとの熱-構造-流体3者連成に対応している。ただし3者連成では収束性がさらに厳しくなるから、連成反復回数の増加とunder-relaxationの調整が重要だよ。

デジタルツインの実現に向けて、流体+構造+熱+電磁場+化学反応のフルマルチフィジックス連成が研究の方向性だ。計算コストの問題はROMやAI代替モデルで対処しつつ、高忠実度シミュレーションを検証データとして使うハイブリッドアプローチが現実的だね。

ALE法で最も頻繁に遭遇するエラーだね。メッシュの変形が局所的に大きすぎてセルが裏返ってしまう状態だ。

デフォルトのUniform diffusivityだと、壁面変位が均等に伝播する。Boundary Distanceだと壁面から離れるほど拡散率が大きくなり、壁面近傍のセル（プリズム層など）の変形が抑制される。FSI問題ではこの設定がほぼ必須だよ。Fluent GUIでは Dynamic Mesh → Smoothing → Diffusion → Diffusivity Based On: Boundary Distance だ。

以下を順番にチェックしよう。

System CouplingのGUIで「Mapping Visualization」を有効にすると、マッピング元と先のメッシュの対応関係が可視化される。未マッピング領域（Orphan point）が多い場合はメッシュの再配置が必要だ。

簡単なテストがある。一様流をFSI界面に処方的な振動変位を与えて（構造ソルバーなしで）計算する。メッシュが動いているだけで流れ場に変化がなければGCLが正しく満足されている。圧力やエネルギーに非物理的な変動が出たらGCL問題だ。

Fluent/STAR-CCM+では標準的にGCLが満足されるから問題になることは稀だけど、OpenFOAMのカスタムソルバーや独自コードでは要注意。特に時間積分が1次精度の場合にGCL誤差が顕在化しやすい。2次精度時間積分（Crank-Nicolson, BDF2等）を使うことで改善される。