心臓弁のFSI — 数値解法と実装
離散化手法の選択
心臓弁のFSI解析で使われる数値手法にはどんな種類がありますか？
心臓弁のように開閉を繰り返す場合、ALE法だと毎サイクルリメッシュが大変ですよね。
そのため近年はIB法やoverset mesh法が主流になりつつある。Griffithら（IBAMR）のオープンソースIBコードは心臓弁の研究で広く使われている。
時間積分
心臓の拍動サイクルをどのくらいの時間刻みで解くんですか？
心拍サイクルは約0.8秒。弁の開閉は数十ミリ秒で起こるから、$\Delta t = 0.1$〜$0.5$ ms程度が必要だ。1サイクルで1,600〜8,000ステップになる。
初期過渡を除くために最低3〜5サイクル回し、統計量は安定したサイクルで取得する。流入境界条件にはMRI計測による流速波形や圧力波形を設定するんだ。
壁面せん断応力の評価
WSSの評価が重要とのことですが、どういう指標を使うんですか？
時間平均WSS（TAWSS）と振動せん断指標（OSI）が代表的だ。
低TAWSS（< 0.4 Pa）かつ高OSI（> 0.3）の領域は血栓リスクが高いとされる。FDA（米国食品医薬品局）のガイダンスでもこれらの評価が推奨されている。
心臓シミュレーション——究極のFSI問題
人間の心臓は1日に約10万回拍動し、血液を全身に送り出します。この過程は流体（血液）-構造（心筋・弁）-電気（刺激伝導系）の3場連成問題。心臓のデジタルツインの構築は連成解析の「聖杯」と呼ばれ、世界中の研究者が挑戦しています。実現すれば、手術のシミュレーションや薬の効果予測が患者ごとにカスタマイズできるようになります。
離散化手法の詳細解説
空間離散化における手法選択が数値精度・安定性・計算コストに与える影響を詳述する。
モノリシック法
全物理場を1つの連立方程式系として同時に解く。強い連成に対して安定だが、実装が複雑でメモリ消費が大きい。
パーティション法（分離反復法）
各物理場を独立に解き、界面でデータ交換。実装が容易で既存ソルバーを活用可能。弱い連成に適する。
界面データ転写
最近傍法（最も簡単だが精度低い）、射影法（保存的）、RBF補間（メッシュ非一致に強い）。保存性と精度のバランスが重要。
マトリクスソルバーの選定指針
問題規模と特性に応じた最適なソルバー選択のガイドライン。
| ソルバー種別 | 詳細・推奨条件 |
|---|---|
| ブロック対角前処理 | 各物理場の前処理を対角ブロックとして使用。実装が容易でスケーラビリティが良い。 |
| フィールド分割法 | 速度-圧力分割（流体）、変位-温度分割（熱-構造）等。物理的な意味に基づくブロック分割が安定性に寄与。 |
時間積分法と収束判定
ソルバー内部の制御パラメータと収束判定基準について記述する。
サブイタレーション
各連成ステップ内で十分な反復を行い、界面条件の整合性を確保。残差基準は各物理場の典型値に基づいてスケーリング。
Aitken緩和
連成反復の緩和係数を自動調整。過緩和による発散を防止し、収束を加速する適応的手法。
安定性条件
added mass効果（流体-構造連成で構造密度≈流体密度の場合）に注意。不安定な場合はロビン型界面条件やIQN-ILS法を適用。
数値解法の直感的理解
連成ソルバーのイメージ
モノリシック手法は「バンドの全員が同じ楽譜を見て同時に演奏する」——完璧な同期が取れるが、楽譜（連立方程式系）が巨大で複雑。パーティション手法は「各パートが別々に練習して合わせる」——個々の練習（ソルバー）は簡単だが、合わせ（データ交換）のタイミングと精度が課題。
Aitken緩和のたとえ
Aitken緩和は「シーソーのバランス取り」に似ている。一方が強く押しすぎると反対側が跳ね上がり、その反動でまた強く押しすぎる——この振動を抑えるために、押す力を自動的に調整するのがAitken緩和。連成反復が振動して収束しないとき、前回の修正量を見て次の修正量を自動調整する適応的手法。
連成解析の安定性やデータ転写の精度は、マルチフィジックスの永続的な課題です。 — Project NovaSolverはこの課題に正面から取り組んでいます。
心臓弁のFSIの実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
実務課題アンケートに回答する →