4-6 速度場を得る三つの計測・解析手法

数値流体解析(CFD)の分解能と境界条件

CFD は分解能を自分で決められる唯一の道具です。ただし細かく解けることと正しいことは別で、代わりに手間を払います。前処理は手作業で、1心周期の計算には日単位の時間がかかります。

三つ目の道具は、測るのではなく解く方向です。血管の形をコンピュータ上に作り、そこに流体の方程式を解いて速度場を得ます。計算流体力学(CFD)です。

前の二つと決定的に違う点が一つあります。分解能を自分で決められます。 MRI のボクセルサイズやエコーのサンプリング間隔は装置が決めますが、CFD の格子は自分で切ります。

血管壁の近傍を、三手法はどう刻んでいるか

分解能の話は、壁際に来ると具体的になります。三つの道具は、空間の刻み方そのものが違います。

1-2第1章で定義したとおり、WSS は壁のすぐ内側の速度勾配です。空間微分なので、壁際の刻みの細かさが値をそのまま決めます。

CFD なら壁際に薄い層を何枚も重ねられます。実際、大動脈弓の検証では最外層の厚みを 30 µm に取った例があります。MRI の 1〜2 mm ボクセルとは二桁違います。だから同じ症例でも WSS の絶対値が数倍ずれます。第6章で読む比較検証1は、CFD 側の分解能を段階的に落とすと差が縮むことを示して、この原因を消去法で確かめたものです。

下の図は、血管壁の同じ場所を三通りに刻んだものです。拡大していくと、壁から 1 mm の中に刻みがいくつ入るかが変わります。

壁際の同じ領域を三つの道具で刻んだ比較。心エコーは扇形の極座標、MRI は直方体のボクセル、CFD は壁際に薄い層を重ねられる。破線は壁から 1 mm で、速度勾配が立っている層のおおよその厚み。この中にいくつ刻みが入るかが、そのまま WSS の値を決める。分解能の値は合意声明と当社の比較表に基づく。

MRI の 2 mm ボクセルでは、勾配が立っている層をボクセルが跨いでしまいます。跨いだ分だけ、傾きは緩く出る側に偏ります。 CFD で最外層を 30 µm に取れば、同じ 1 mm の中に十数層が入ります。

CFD の解像度と境界条件の区別

ここを混同すると危険です。CFD は計測ノイズを持たず、いくらでも細かく解けます。しかしそれは「入力が正しければ細かく解ける」という意味にとどまります。

CFD の値は真値ではありません。 粘度の仮定、境界条件、壁のモデル、形状の抽出に依存します。三手法のうち壁近傍の分解能を最も高く設定できるのは CFD ですが、それは「絶対値の基準になる」という意味ではありません。

同じ画像を26チームに渡した国際チャレンジ2では、瘤内平均 WSS のばらつきが変動係数48%に達しました。ソルバーの精度が足りなかったからではありません。形状の切り出し方と流入条件の与え方が割れたからです。

分解能は自分で決められます。境界条件は自分で決めなければなりません。 後者が、CFD の結果を左右します。

CFD の所要時間と前処理の工数

CFD が他の二つと違うもう一つの点は、時間です。

心エコー(VFM) 4D Flow MRI CFD
所要時間 10分程度 数時間 1日〜数日
次元 2次元 2次元・3次元 3次元
壁の動き 扱える 扱える 難しい
空間分解能の目安 0.1〜1 mm 1〜3 mm 0.03〜3 mm(自分で決める)
向く対象 心室・心房 心室・心房・大血管 血管

所要時間の内訳が実務では効きます。Fontan の手術計画に関する総説3では、形状の切り出しに作業時間の約40%がかかり、1つの手術案あたり人手で約10時間、計算に64プロセッサで約48時間という数字が挙げられています。

前処理——形状の抽出と平滑化——は、いまも大半が手作業です。ここが自動化されていないことが、CFD を日常診療に乗せられていない理由の一つです。

流体構造連成(FSI)と剛体壁仮定

表の「壁の動き」の行に注目してください。エコーと MRI は実際の壁の動きを撮っていますが、CFD で壁を動かすのは難しくなります。

流れと構造の変形を同時に解く流体構造連結解析(FSI)という手法があります。ただし弁の開閉のように変形が大きい問題では、格子の変形が破綻しやすくなります。加えて、血管の弾性特性を一様と仮定せざるを得ないことが多く、動脈硬化のように性質が場所で違う血管では検証そのものが難しくなります。

多くの解析が剛体壁を仮定しているのは、この事情です。剛体壁の仮定で WSS が最大30%変わるという報告があり、これは仮定として明記すべき事項になります。もっとも 頸動脈分岐で両者を比べた研究4は、低ずり域の広さのような指標であれば剛体壁で足りることも示しています。扱いは 6-4第6章で整理します。

測りたい量から逆算する手法の選択

第4章で三つの道具を見ました。整理すると、選択の基準は一つに絞れます。

何を測りたいのか。

流量や逆流分画5なら、実測できる MRI かエコーで足ります。壁近傍の分解能を最も高く設定できるのは CFD です。ベッドサイドで繰り返したいならエコーが向きます。まだ存在しない形状を評価したい6なら CFD だけができます。

そして、どの道具を使っても——ここまで繰り返してきたとおり——値は条件を含んだ量です。 装置が決めた条件か、自分が決めた条件かの違いはありますが、条件を書かずに値だけを比べることはできません。

この記事で出てきた言葉

計算流体力学(CFD)
血管の形状をコンピュータ上に作り、流体の方程式を解いて速度場を求める手法。分解能を自分で決められる代わりに、境界条件も自分で決めなければならない。
空間の離散化
連続した空間を有限の要素に区切ること。エコーは極座標、MRI はボクセル、CFD は多面体格子で区切る。壁際の刻みの細かさが WSS の値を決める。

引用文献(6 本)

  1. 1

    Perinajová R, Juffermans JF, Mercado JL ほか. Assessment of turbulent blood flow and wall shear stress in aortic coarctation using image-based simulations. BioMedical Engineering OnLine. 2021.↩

    条件を完全に揃えても CFD の WSS は MRI の3.7倍 原因は乱流モデルでなく壁際の分解能

  2. 2

    Valen-Sendstad K, Bergersen AW, Shimogonya Y ほか. Real-World Variability in the Prediction of Intracranial Aneurysm Wall Shear Stress: The 2015 International Aneurysm CFD Challenge. Cardiovascular Engineering and Technology. 2018.↩

    同じ画像、26チーム WSS のばらつき48%、母血管で正規化すれば18%

  3. 3

    Trusty PM, Slesnick TC, Wei ZA ほか. Fontan Surgical Planning: Previous Accomplishments, Current Challenges, and Future Directions. Journal of Cardiovascular Translational Research. 2018.↩

    Fontan の CFD、約束できるのは肝静脈血流の分布まで 転帰が良くなるかは未確認

  4. 4

    Zambon S, Mazzi V, Calò K ほか. Influence of wall distensibility on local hemodynamics at the normal carotid bifurcation: a fluid-structure interaction study. Biomechanics and Modeling in Mechanobiology. 2026.↩

    頸動脈分岐は剛体壁 CFD で十分 FSI との面積差は中央値 4.1%

  5. 5

    Dyverfeldt P, Bissell M, Barker AJ ほか. 4D flow cardiovascular magnetic resonance consensus statement. Journal of Cardiovascular Magnetic Resonance. 2015.↩

    WSS・TKE・相対圧を「研究段階」に置いた最初の国際合意 線引きは2026年も現役

  6. 6

    de Leval MR, Dubini G, Migliavacca F ほか. Use of computational fluid dynamics in the design of surgical procedures: application to the study of competitive flows in cavopulmonary connections. The Journal of Thoracic and Cardiovascular Surgery. 1996.↩

    下大静脈吻合を右肺動脈側へ広げると損失最小 CFD 術式設計の起点となった1996年

ここで答えていない問い

  • 計測を最小にする道、あるいは閾値から逃れる道はあるのか。

第4章「速度場を得る三つの計測・解析手法」より 最終更新 2026-09-08