ダムブレイク流体-構造連成 — トラブルシューティングガイド
よくある問題と対策
ダムブレイクFSIの典型的なトラブルを教えてください。
1. VOFの界面拡散
症状: 自由表面がぼやけて、水と空気の境界が不明瞭になる。
原因: 数値拡散、メッシュが粗い、時間刻みが大きい。
対策:
- OpenFOAMではcAlpha（界面圧縮係数）を1〜2に設定
- HRIC（High Resolution Interface Capturing）スキームを使用
- AMR（Adaptive Mesh Refinement）で界面付近を自動細分化
2. SPHの引張不安定性
症状: 粒子間に空隙が発生し、非物理的な粒子分布になる。
原因: SPHのテンシル不安定性（負圧下で粒子が離散する）。
対策:
- 人工粘性の追加（Monaghan型）
- XSPH速度補正の適用
- δ-SPH（密度拡散項追加）の使用
3. LS-DYNAのALE-FSIでリーク
流体が構造を通り抜けてしまうのはなぜですか？
症状: *CONSTRAINED_LAGRANGE_IN_SOLIDで流体が構造壁を貫通する。
原因: ペナルティ係数（PFAC）が小さい、またはALEメッシュが粗い。
対策:
- PFACを0.1から1.0に段階的に増加させてテスト
- 構造壁の厚さ方向にALE要素を最低2層配置
- NQUAD（流体-構造間の積分点数）を増加
4. 衝撃圧のメッシュ非収束
症状: メッシュを細かくするほど衝撃圧ピークが上昇し続ける。
対策:
- インパルス（力の時間積分）の収束性を代わりに確認
- 空気の圧縮性を考慮したモデルに切り替え
- 実験データのばらつき範囲と比較して妥当性を判断
| チェック項目 | 基準 |
|---|---|
| 洪水波先端の到達時間 | 理論解（Ritter解）と比較 |
| 水位時刻歴 | 実験データと10%以内 |
| 衝撃圧インパルス | メッシュ3水準で収束確認 |
| 構造変位 | 実験データと20%以内 |
衝撃圧の絶対値よりもインパルスの収束を確認するのが実務的なコツですね。
リバティ船の溶接割れ——連成問題の教訓
第二次世界大戦中、アメリカは「リバティ船」を溶接で大量生産し、戦争の物流を支えました。しかし約1,500隻のうち約400隻に船体の亀裂が発生。原因は溶接残留応力と低温脆性の連成——溶接時の急激な温度変化が残留応力を生み、北大西洋の冷たい海水で鋼材が脆くなり、亀裂が伝播したのです。現代の溶接シミュレーションは、この「温度→残留応力→破壊」の連鎖を予測できます。
トラブル解決の考え方
デバッグのイメージ
連成解析のトラブルシューティングは「チームプレーの問題解決」に似ている。まず「どのチーム（物理場）に問題があるか」を切り分け、次に「チーム間の連携（データ転写）に問題がないか」を確認する。各物理場を単独で動かして問題がなければ、連成の設定が原因。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ダムブレイク流体-構造連成の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
連成解析の安定性やデータ転写の精度は、マルチフィジックスの永続的な課題です。 — Project NovaSolverはこの課題に正面から取り組んでいます。
CAEの未来を、実務者と共に考える
Project NovaSolverは、ダムブレイク流体-構造連成における実務課題の本質に向き合い、エンジニアリングの現場を支える道具づくりを目指す研究開発プロジェクトです。
プロジェクトの最新情報を見る →