ALE法によるFSI — トラブルシューティング
より充実した内容を ale-fsi.html でご覧いただけます。
負体積セルの発生
「Negative cell volume detected」というエラーが出て計算が止まるんですが...
ALE法で最も頻繁に遭遇するエラーだね。メッシュの変形が局所的に大きすぎてセルが裏返ってしまう状態だ。
| 対策 | 効果 | 実装方法 |
|---|---|---|
| Diffusivity=Boundary Distance | 壁面近傍のメッシュ変形を抑制 | Fluent Dynamic Mesh設定 |
| 時間刻みを小さくする | 1ステップあたりの変形量を減らす | $\Delta t$ を半分に |
| リメッシング有効化 | 品質劣化セルを再生成 | Fluent Remeshing ON |
| Overset meshに切り替え | メッシュ変形自体を回避 | Overset Interface作成 |
| 層数の増加 | プリズム層の余裕を増やす | メッシュ再生成 |
DiffusivityをBoundary Distanceにすると何が変わるんですか？
デフォルトのUniform diffusivityだと、壁面変位が均等に伝播する。Boundary Distanceだと壁面から離れるほど拡散率が大きくなり、壁面近傍のセル（プリズム層など）の変形が抑制される。FSI問題ではこの設定がほぼ必須だよ。Fluent GUIでは Dynamic Mesh → Smoothing → Diffusion → Diffusivity Based On: Boundary Distance だ。
連成反復が収束しない
System Couplingの反復が最大回数に達しても収束しないんですが...
以下を順番にチェックしよう。
1. Under-relaxation factorを下げる: 0.5→0.3→0.1 と段階的に。密度比が小さいほど小さい値が必要
2. マッピングを確認: 構造メッシュと流体メッシュのFSI界面が空間的に一致しているか。ギャップがあると力の転送が正しくない
3. 初期条件を改善: 静的釣り合い状態（静荷重のみの定常解）からスタートする
4. 時間刻みを小さくする: 非定常問題の各ステップでの変化量を減らす
5. Ramping: 初期数ステップで荷重を線形に増加させる
マッピングのギャップはどうやって確認しますか？
System CouplingのGUIで「Mapping Visualization」を有効にすると、マッピング元と先のメッシュの対応関係が可視化される。未マッピング領域（Orphan point）が多い場合はメッシュの再配置が必要だ。
GCL誤差の確認
GCL誤差が疑わしい場合、どうやって確認すればいいですか？
簡単なテストがある。一様流をFSI界面に処方的な振動変位を与えて（構造ソルバーなしで）計算する。メッシュが動いているだけで流れ場に変化がなければGCLが正しく満足されている。圧力やエネルギーに非物理的な変動が出たらGCL問題だ。
Fluent/STAR-CCM+では標準的にGCLが満足されるから問題になることは稀だけど、OpenFOAMのカスタムソルバーや独自コードでは要注意。特に時間積分が1次精度の場合にGCL誤差が顕在化しやすい。2次精度時間積分（Crank-Nicolson, BDF2等）を使うことで改善される。
ALE解析が「謎の発散」で止まる——あるあるな落とし穴
ALE-FSI解析の現場でよく聞く失敗談がある。「設定はちゃんとしたはずなのに、数百ステップで計算が突然発散した」というものだ。大抵の原因はメッシュスムージングの頻度が低すぎて、境界付近の要素が徐々に歪み、ある瞬間に負体積が発生すること。「エラーメッセージを読まずに時間刻みを細かくして再計算」を繰り返す人がいるが、根本はメッシュ品質の問題なので全く解決しない。エラーメッセージは本当のことを言っています。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ALE法によるFSIの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告