埋め込み境界法（IBM）

ALE法ではメッシュが構造界面に適合（body-conforming）するのに対し、IBMでは固定の直交格子（Cartesian grid）上に構造を「埋め込む」。構造の境界はメッシュと無関係に定義され、境界条件を体積力（forcing term）や補間で課すんだ。

Peskin（1972）が心臓弁の血流シミュレーションのために考案した手法で、元々の定式化はこうなる。

Peskinの正規化デルタ関数を使う。4点離散デルタ関数が標準的で、

大きく2系統に分けられる。

Continuous forcing（Peskin型）は実装が簡単だけど、界面がデルタ関数のサポート幅分だけぼやける。壁面の $y^+$ 制御が難しいため、高Re数の乱流壁面流れには不向きだ。

厳密にはCut-cell法はIBMの一種で、界面がセルを横切る場所でセルを幾何学的に切断して、切断セル上で保存則を満たす。保存性が高く精度も良いが、切断セルの体積が非常に小さくなる問題（small cell problem）への対策が必要だ。Cell mergingやflux redistributionで対処する。

Direct Forcing法（Mohd-Yusof, 1997; Fadlun et al., 2000）が最も広く使われている。界面近傍のEulerianセルで、速度を所望の境界条件値に強制する。

界面がセルの中間にある場合は1次精度。改良版では界面からの距離に基づく補間を使って2次精度を達成する。Ghost cell法では界面の内側に仮想セル（ghost cell）を置いて、反射補間で境界条件を課す。

IBM-FSI連成の流れはこうだ。

ALE法と比較したIBM-FSIの利点と欠点は以下の通りだ。

その通り。心臓弁の開閉、パラシュートの展開、旗のはためきなど、構造が大きく変形してトポロジーが変わるような問題はIBMの独壇場だ。

限定的だが存在する。

厳密なIBMが必要な場合は、研究コード（Nek5000、CaNS、PeleLM、AFiDなど）を使うか、OpenFOAMにカスタム実装するのが現実的だよ。

以下が代表的な検証問題だ。

Lai & Peskin (2000)、Uhlmann (2005)の論文がIBM文脈での標準的な参照だ。$Re = 100$ でのストロハル数 $St \approx 0.164$、抗力係数 $C_D \approx 1.35$、揚力係数振幅 $C_L' \approx 0.34$ が目安だよ。格子幅 $h$ に対して円柱直径が $D/h = 30-50$ 程度あれば十分な解像度だ。

ESI版OpenFOAM（v2006以降）にはimmersedBoundaryFvPatchが含まれている。Ghost cell型のIBMで、以下の手順で使える。

sixDoFRigidBodyMotionと組み合わせて、物体の並進・回転運動を6DoFソルバーで計算し、IBMマスクを毎ステップ更新する。FSIの場合は構造ソルバーとの連成が必要で、preCICEとの組み合わせが有力だ。

LBMは直交等間隔格子上で動作するから、IBMとの組み合わせが非常に自然だ。Palabos（オープンソース）やPowerFLOW（商用、ダッソー・システムズ）がLBM + IBMを採用している。

PowerFLOWは自動車の空力騒音解析で広く使われていて、内部的にはLattice Boltzmann + IBM（Very Large Eddy Simulation）を基盤としている。複雑な車体形状のメッシュ生成が自動化されるのが大きな利点だよ。

Overset（Chimera）法は複数のbody-conformingメッシュを重ね合わせる手法で、各メッシュは物体形状に適合している。IBMとの根本的な違いは壁面メッシュの有無だ。

その通り。高Re数の乱流壁面流れ（航空機翼など）ではOversetが事実上の標準だ。一方、低〜中Re数のFSI（生体流れ、マイクロ流体）や、接触を含む問題（心臓弁など）ではIBMが有利。

以下の手順だ。

Overset FSIの注意点は、コンポーネントメッシュが回転や並進しても背景メッシュとの重なり領域が十分に確保されること。重なり領域が狭すぎると補間が破綻する。最低3-5セル幅のオーバーラップ領域を確保しよう。

Fluent 2023R2以降で対応している。Oversetメッシュのコンポーネント面をSystem Coupling Regionとして定義し、Ansys Mechanicalと力・変位を交換できる。翼の大振幅フラッターやタービンブレードの変形解析に有効だよ。

以下の3方向が活発だ。

NVIDIA Omniverse上のPhysX Flow（旧Flow Simulator）はGPUネイティブのIBMベースCFDソルバーで、リアルタイムの流体可視化に使われている。産業用途ではまだ精度検証が不十分だけど、デジタルツインのリアルタイム流体シミュレーションへの応用が期待されているよ。

構造のモデリングは問題によって異なるが、代表的な3パターンがある。

薄い弾性構造（旗、膜、弁）のIBM-FSIでは、構造をLagrangianメッシュで表現し、Peskin型のデルタ関数で流体とカップリングするのが伝統的なアプローチだ。構造の支配方程式は

Zhu & Peskin (2002)、Connell & Yue (2007)の研究が代表的だ。旗のはためきはIBM-FSIの標準ベンチマークの一つで、質量比 $M^* = \rho_s h_s / (\rho_f L)$ と無次元曲げ剛性 $KB = EI / (\rho_f U^2 L^3)$ の2パラメータで振動特性が決まる。

心臓弁の開閉を考えてみよう。弁尖（leaflet）が閉じるとき、2つの弁尖が接触する。ALE法ではメッシュが潰れてしまうけど、IBMでは弁尖がEulerian格子上を自由に動くから接触自体は問題にならない。

ただし、接触時の流体の扱いには工夫が必要だ。2つの弁尖の間に流体が残っていると非物理的な圧力が発生する。反発力モデルや最小ギャップモデルで弁の完全閉鎖を模擬する手法が使われている。

Particle-Resolved DNS（PR-DNS）はIBMの重要な応用分野だ。数百から数千の粒子を流体中に配置し、各粒子周りの流れを完全に解像する。

代表的な研究コードを紹介しよう。

直交等間隔格子 + FFTベースの高速Poisson解法のおかげで、CaNSは非常に効率が良い。GPUにも対応していて、単一GPUで数億格子点の計算が可能だ。粒子-流体相互作用の基礎研究に威力を発揮しているよ。

Peskin型の界面のぼやけ（diffuse interface）はデルタ関数のサポート幅に起因する。改善策は以下の通り。

壁面の精度が重要（Re数が高い、壁面熱伝達が必要）ならSharp interface法（Ghost cell、Cut-cell）を選ぶ。構造が薄くて柔軟（弁、膜）で、壁面境界層の精密な解像がそこまで重要でないならPeskin型でも十分だ。

Direct Forcing型IBMでは、物体内部のセルを「Dead cell」として速度をゼロにするが、速度場の連続性が完全には保証されないため、微小な質量流入/流出が生じることがある。対策として、

蓄積する。特に周期運動する物体（振動する翼、回転するプロペラ）では、サイクルごとに微小な体積誤差が蓄積して非物理的な圧力場の変動が生じる。質量保存の補正は長時間計算では必須と考えた方がいい。

物体のサイズ $D$ に対する格子幅 $h$ の比が重要だ。

そう。だからこそAMR（Adaptive Mesh Refinement）との組み合わせが重要なんだ。Basiliskは8レベルの八分木AMRを標準搭載していて、界面近傍だけ高解像度にできる。等間隔格子比で実効的に10-100倍の効率化が可能だよ。