心臓弁のFSI — トラブルシューティングガイド
より充実した内容を heart-valve.html でご覧いただけます。
弁葉接触時の計算破綻
弁が閉じるときにメッシュが潰れて計算が止まるんですが。
ALE法の最大の弱点だ。対策は3つある。
1. 接触ギャップの設定: 弁葉間に最小ギャップ（0.05〜0.1 mm）を維持し、完全閉鎖を避ける
2. overset mesh: STAR-CCM+やAnsys Fluentのoverlapping grid機能で弁葉ごとに独立メッシュを使用
3. IB法への切り替え: 格子を固定することで接触問題を回避
最小ギャップを設けると逆流量が過大評価されませんか？
良い指摘だ。ギャップ設定は結果に影響するから、ギャップ感度解析を実施すべきだ。複数のギャップ値（0.02, 0.05, 0.1 mm）で解析して、regurgitant volumeへの影響を定量評価する。
非ニュートンモデルの収束問題
Carreau-Yasudaモデルを使うと収束が悪くなるんですが。
せん断率がゼロに近い領域で粘度が急激に上昇するのが原因だ。対策として、
- 粘度の上限値を設定する（例：$\mu_{max} = 0.1$ Pa·s）
- 初期はニュートン流体で計算し、途中から非ニュートンに切り替える
- Under-relaxation factorを0.5〜0.7程度に下げる
生理学的妥当性の確認
計算結果が正しいかどうかの判断基準はありますか？
以下の生理学的パラメータと比較するのが標準だ。
| パラメータ | 正常範囲 | 確認方法 |
|---|---|---|
| 大動脈弁口面積（EOA） | 3.0〜4.0 cm² | 弁口でのjet面積から算出 |
| 最大ジェット速度 | 1.0〜1.5 m/s | 弁口中央の最大速度 |
| 逆流率 | < 5% | 閉鎖時逆流量/拍出量 |
| 圧力損失 | < 15 mmHg | 弁前後の圧力差 |
| TAWSS（大動脈洞） | 0.5〜2.0 Pa | 壁面せん断応力の時間平均 |
これらの値から大きく外れていたら、モデルを見直す必要があるんですね。
その通り。特にEOAとpressure gradientは臨床のエコー検査データと直接比較できるから、V&Vの最初のステップとして有用だ。
「弁が完全に閉じないまま計算が止まった」——接触収束失敗あるある
心臓弁FSI解析でよく遭遇するのが「弁尖の接触判定が収束しない」という問題です。弁が閉じる瞬間、圧力と変位が急変するため、連成反復が振動して収束しません。まず確認すべきは「時間刻み幅（dt）が閉鎖時刻をまたいでいないか」です。実際の弁閉鎖は1〜2ms以内に完了するため、dt=1ms以上だと閉鎖の瞬間を計算が飛び越えてしまいます。次に「接触剛性（ペナルティ剛性）が高すぎないか」——高すぎると数値振動が起き、低すぎると弁尖が貫通します。経験則では、接触剛性は弁尖の曲げ剛性の100〜1000倍程度に設定するのが良いとされています。それでも収束しない場合は、弁尖に数値減衰（Rayleigh減衰）を少し加えるのも有効な手です。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——心臓弁のFSIの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告