埋め込み境界法（IBM） — 先端手法と応用
弾性構造のIBM-FSI
弾性体の大変形FSIにIBMを使う場合、構造側はどう解くんですか？
構造のモデリングは問題によって異なるが、代表的な3パターンがある。
| 構造モデル | 対象 | 解法 | 代表例 |
|---|---|---|---|
| 1D弾性繊維 | フィラメント、旗 | 有限差分/FEM | 旗のはためき |
| 2Dシェル | 膜、弁 | MITC shell要素 | 心臓弁 |
| 3Dソリッド | 厚肉構造 | FEM（Abaqus等） | 血管壁 |
薄い弾性構造（旗、膜、弁）のIBM-FSIでは、構造をLagrangianメッシュで表現し、Peskin型のデルタ関数で流体とカップリングするのが伝統的なアプローチだ。構造の支配方程式は
ここで $\mathbf{T}$ は膜張力テンソル、$\mathbf{F}_{fluid}$ は流体から受ける力だ。
旗のはためき問題は教科書でよく見ますね。
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パラメータで振動特性が決まる。
接触を含むFSI
IBMで接触問題が解けるというのは具体的にどういうことですか？
心臓弁の開閉を考えてみよう。弁尖（leaflet）が閉じるとき、2つの弁尖が接触する。ALE法ではメッシュが潰れてしまうけど、IBMでは弁尖がEulerian格子上を自由に動くから接触自体は問題にならない。
ただし、接触時の流体の扱いには工夫が必要だ。2つの弁尖の間に流体が残っていると非物理的な圧力が発生する。反発力モデルや最小ギャップモデルで弁の完全閉鎖を模擬する手法が使われている。
粒子浮遊系のIBM
流体中の粒子群の運動にもIBMが使えますか？
代表的な研究コードを紹介しよう。
| コード | 開発元 | 手法 | 粒子数 | 応用 |
|---|---|---|---|---|
| PUReIBM | University of Iowa | Direct Forcing | ~10,000 | 流動層、沈降 |
| Basilisk | CNRS | VOF + IBM | ~1,000 | 液滴、気泡 |
| CaNS | TU Delft | IBM + FFT | ~100,000 | チャネル流中粒子 |
10万粒子のPR-DNSってすごい計算量ですね。
直交等間隔格子 + FFTベースの高速Poisson解法のおかげで、CaNSは非常に効率が良い。GPUにも対応していて、単一GPUで数億格子点の計算が可能だ。粒子-流体相互作用の基礎研究に威力を発揮しているよ。
ライト兄弟は最初の「CFDエンジニア」だった？
ライト兄弟は1901年に自作の風洞で200以上の翼型を試験しました。当時のコンピュータは？ もちろん存在しません。彼らは手作業で揚力と抗力を測定し、最適な翼型を見つけ出した。現代のCFDエンジニアがFluent1発で計算する揚力係数を、ライト兄弟は何百回もの風洞実験で手に入れたのです。
先端技術を直感的に理解する
この分野の進化のイメージ
CFDの最先端は「天気予報の進化」に似ている。かつての天気予報（RANS）は平均的な傾向しか分からなかったが、最新の数値天気予報（LES/DNS）は個々の雲の動きまでシミュレーションできる。AIとの融合により「数秒で近似予測」も可能になりつつある。
なぜ先端技術が必要なのか — 埋め込み境界法（IBM）の場合
従来手法で埋め込み境界法（IBM）を解析すると、計算時間・精度・適用範囲に限界がある。例えば、設計パラメータを100通り試したい場合、従来手法では100回の解析が必要だが、サロゲートモデルを使えば数回の解析結果から100通りの予測が可能になる。「全部試す」から「賢く推測する」への転換が先端技術の本質。
CFDメッシュの品質管理や乱流モデルの選定に悩む時間を、もっと創造的な設計作業に使えたら。 — Project NovaSolverはそんな実務者の声から生まれました。
Project NovaSolver — CAE実務の課題に向き合う研究開発
「埋め込み境界法（IBM）をもっと効率的に解析できないか？」——私たちは実務者の声に耳を傾け、既存ワークフローの改善を目指す次世代CAEプロジェクトに取り組んでいます。具体的な機能はまだ公開前ですが、開発の進捗をお届けします。
進捗通知を受け取る →