血行動態シミュレーション — 数値解法と実装

カテゴリ: 連成解析 | 2026-01-20
hemodynamics-method
数値解法の舞台裏

安定化有限要素法

🧑‍🎓

血流CFDではSUPG/PSPG安定化が必要だと聞きましたが、なぜですか?


🎓

等次補間(P1-P1)はinf-sup条件を満たさないため、圧力にスプリアス振動が出る。SUPG/PSPG安定化を追加することで安定な解が得られる。


$$ \sum_e \int_{\Omega_e} \tau_{SUPG} (\mathbf{u} \cdot \nabla \mathbf{w}) \cdot \mathbf{R} \, d\Omega = 0 $$

🎓

血管分岐部やステント周辺では局所Re数が高く、安定化パラメータの設計が計算安定性に直結する。


連成アルゴリズム

🧑‍🎓

血管のFSIは弱連成で大丈夫ですか?


🎓

血管壁は薄く付加質量効果が顕著だから、弱連成では発散しやすい。Robin-Neumann分割やMonolithic法が推奨される。Generalized Robin条件を使った半陰的連成法が注目されている。


$$ \alpha_f \mathbf{u}_f + \boldsymbol{\sigma}_f \cdot \mathbf{n} = \alpha_f \dot{\mathbf{d}}_s + \boldsymbol{\sigma}_s \cdot \mathbf{n} $$

画像ベースモデリングのワークフロー

🧑‍🎓

臨床画像からメッシュを作るまでの流れを教えてください。


🎓

1. CT血管造影またはMRAの取得

2. セグメンテーション(Mimics, ITK-SNAP, 3D Slicer)

3. 表面スムージングと欠損修復

4. 体積メッシュ生成(TetGen, VMTK, ICEM CFD

5. 境界層メッシュ挿入(VMTKのboundary layer機能)

6. 入口・出口面の特定と延長管追加


VMTK(Vascular Modeling Toolkit)は血管専用OSSで、中心線抽出から境界層メッシュまでを一貫処理できる。

Coffee Break よもやま話

心臓シミュレーション——究極のFSI問題

人間の心臓は1日に約10万回拍動し、血液を全身に送り出します。この過程は流体(血液)-構造(心筋・弁)-電気(刺激伝導系)の3場連成問題。心臓のデジタルツインの構築は連成解析の「聖杯」と呼ばれ、世界中の研究者が挑戦しています。実現すれば、手術のシミュレーションや薬の効果予測が患者ごとにカスタマイズできるようになります。

離散化手法の詳細解説

空間離散化における手法選択が数値精度・安定性・計算コストに与える影響を詳述する。

モノリシック法

全物理場を1つの連立方程式系として同時に解く。強い連成に対して安定だが、実装が複雑でメモリ消費が大きい。

パーティション法(分離反復法

各物理場を独立に解き、界面でデータ交換。実装が容易で既存ソルバーを活用可能。弱い連成に適する。

界面データ転写

最近傍法(最も簡単だが精度低い)、射影法(保存的)、RBF補間(メッシュ非一致に強い)。保存性と精度のバランスが重要。

マトリクスソルバーの選定指針

問題規模と特性に応じた最適なソルバー選択のガイドライン。

ソルバー種別詳細・推奨条件
ブロック対角前処理各物理場の前処理を対角ブロックとして使用。実装が容易でスケーラビリティが良い。
フィールド分割法速度-圧力分割(流体)、変位-温度分割(熱-構造)等。物理的な意味に基づくブロック分割が安定性に寄与。

時間積分法と収束判定

ソルバー内部の制御パラメータと収束判定基準について記述する。

サブイタレーション

各連成ステップ内で十分な反復を行い、界面条件の整合性を確保。残差基準は各物理場の典型値に基づいてスケーリング。

Aitken緩和

連成反復の緩和係数を自動調整。過緩和による発散を防止し、収束を加速する適応的手法。

安定性条件

added mass効果(流体-構造連成で構造密度≈流体密度の場合)に注意。不安定な場合はロビン型界面条件やIQN-ILS法を適用。

数値解法の直感的理解

連成ソルバーのイメージ

モノリシック手法は「バンドの全員が同じ楽譜を見て同時に演奏する」——完璧な同期が取れるが、楽譜(連立方程式系)が巨大で複雑。パーティション手法は「各パートが別々に練習して合わせる」——個々の練習(ソルバー)は簡単だが、合わせ(データ交換)のタイミングと精度が課題。

Aitken緩和のたとえ

Aitken緩和は「シーソーのバランス取り」に似ている。一方が強く押しすぎると反対側が跳ね上がり、その反動でまた強く押しすぎる——この振動を抑えるために、押す力を自動的に調整するのがAitken緩和。連成反復が振動して収束しないとき、前回の修正量を見て次の修正量を自動調整する適応的手法。

連成解析の安定性やデータ転写の精度は、マルチフィジックスの永続的な課題です。 — Project NovaSolverはこの課題に正面から取り組んでいます。

CAEの未来を、実務者と共に考える

Project NovaSolverは、血行動態シミュレーションにおける実務課題の本質に向き合い、エンジニアリングの現場を支える道具づくりを目指す研究開発プロジェクトです。

プロジェクトの最新情報を見る →