埋め込み境界法（IBM） — トラブルシューティング
より充実した内容を immersed-boundary.html でご覧いただけます。
界面のスムージング（ぼやけ）
Peskin型IBMで界面がぼやけるんですが、改善方法はありますか？
Peskin型の界面のぼやけ（diffuse interface）はデルタ関数のサポート幅に起因する。改善策は以下の通り。
| 対策 | 効果 | 副作用 |
|---|---|---|
| 格子を細かくする | 界面が相対的にシャープに | 計算コスト増大 |
| 狭いデルタ関数（2点） | ぼやけ減少 | 不安定性増加 |
| Sharp interface法に切替 | 根本的に解決 | 実装が複雑 |
| AMR（適合細分化） | 界面近傍のみ高解像 | AMR実装が必要 |
Sharp interface法とDiffuse interface法のどちらを選ぶべきですか？
壁面の精度が重要（Re数が高い、壁面熱伝達が必要）ならSharp interface法（Ghost cell、Cut-cell）を選ぶ。構造が薄くて柔軟（弁、膜）で、壁面境界層の精密な解像がそこまで重要でないならPeskin型でも十分だ。
体積保存の問題
IBMで閉じた物体内部の体積が保存されないんですが...
Direct Forcing型IBMでは、物体内部のセルを「Dead cell」として速度をゼロにするが、速度場の連続性が完全には保証されないため、微小な質量流入/流出が生じることがある。対策として、
1. 質量ソース補正: 物体表面を通過する質量フラックスを毎ステップ計算し、補正ソース項で打ち消す
2. Pressure correction: 物体内部の圧力方程式を修正して質量保存を強制
3. Cut-cell法への移行: 幾何学的に正確なセル体積を使うため保存性が高い
この問題は長時間計算で蓄積しますか？
蓄積する。特に周期運動する物体（振動する翼、回転するプロペラ）では、サイクルごとに微小な体積誤差が蓄積して非物理的な圧力場の変動が生じる。質量保存の補正は長時間計算では必須と考えた方がいい。
IBMのメッシュ解像度指針
IBMで必要なメッシュ解像度の目安を教えてください。
物体のサイズ $D$ に対する格子幅 $h$ の比が重要だ。
| 用途 | $D/h$ の目安 | 備考 |
|---|---|---|
| 定性的な流れパターン | 10-20 | ストロハル数は±10% |
| 定量的な力の評価 | 30-50 | $C_D$ は±5%程度 |
| 壁面せん断応力 | 50-100 | 境界層の解像が必要 |
| DNS品質 | 100+ | 壁面摩擦の精密計算 |
$D/h = 50$ として3D球だと、球の直径方向に50セル分ですね。計算領域全体だとかなりの規模になりそうです。
そう。だからこそAMR（Adaptive Mesh Refinement）との組み合わせが重要なんだ。Basiliskは8レベルの八分木AMRを標準搭載していて、界面近傍だけ高解像度にできる。等間隔格子比で実効的に10-100倍の効率化が可能だよ。
IBMの「スプリアス圧力振動」——誰もが一度は悩む壁付近の乱れ
IBMを使い始めたエンジニアが最初に悩むのが、壁面付近の圧力場に生じる非物理的な振動（スプリアス振動）だ。固体と流体の境界で格子を「切り替える」際に生じる数値誤差で、流れ場の可視化ビジュアルにモワレ模様のような圧力縞が現れる。これを見た人は「計算がおかしい」と焦るが、多くは空間解像度の問題か強制項の実装方法の問題だ。美しいビジュアルが出るまでがIBMの「洗礼」と呼ばれています。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——埋め込み境界法（IBM）の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告