ダムブレイク流体-構造連成 — トラブルシューティングガイド
この記事は統合版に移行しました
より充実した内容を dam-break-fsi.html でご覧いただけます。
より充実した内容を dam-break-fsi.html でご覧いただけます。
よくある問題と対策
ダムブレイク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%以内 |
衝撃圧の絶対値よりもインパルスの収束を確認するのが実務的なコツですね。
Coffee Break よもやま話
「衝撃圧が実験の10倍出る」——数値的な圧力スパイクの正体
ダムブレイクFSI解析でよくある困り事が「衝突直後に非物理的な圧力スパイクが出る」問題です。実験では数kPaの衝撃圧なのにシミュレーションでは数十kPaが出て構造が壊れてしまう、というケース。原因の多くは自由表面の数値的な不安定性で、VOF法でのinterfaceが「数値的に飛ぶ」現象です。対策はタイムステップをCFL数0.1以下に絞ること、界面の人工的な平滑化を控えること、そして圧力ポアソン方程式のソルバー収束基準を厳しくすることです。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ダムブレイク流体-構造連成の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
この記事の評価
ご回答ありがとうございます！
参考に
なった
なった
もっと
詳しく
詳しく
誤りを
報告
報告