壁の伸展性は頸動脈分岐の局所血行動態をどこまで変えるか — FSI と剛体壁 CFD の対比較
Influence of wall distensibility on local hemodynamics at the normal carotid bifurcation: a fluid-structure interaction study
頸動脈の WSS 解析は剛体壁の CFD で足りるのか。壁の動きを入れた FSI にする必要があるのか。
健常形態の頸動脈分岐10例で二方向連成 FSI と剛体壁 CFD を対で比較したところ、低 TAWSS・高 OSI・高 TSVI にさらされる面積の差は中央値で 4.1%・1.4%・2.3%、分布の重なりを表す類似度指数は 0.83・0.79・0.68 でした。臨床的に意味のある血行動態の特徴は剛体壁 CFD で捉えられます。FSI が要るのは壁応力そのものを見たいときです。
当社の視点
「壁を固定して計算していて本当に良いのか」は、当社が繰り返し受けるご相談です。この文献はそれに正面から答えたもので、同じ10形状に剛体壁 CFD と FSI を対で走らせ、差を数字で示しています。結論は「頸動脈分岐の範囲では剛体壁で足りる」。特筆すべきは、比較対象に TSVI（WSS 位相幾何の時間変動）やヘリシティといった、壁の動きに敏感と考えられていた新しい指標まで含めている点です。同時に著者らは、この結論が変形量に依存することも明記しています。肺動脈自家移植片のように大きく動く構造では剛体壁仮定が破綻した先行報告を引き、頸動脈で成立したことをそのまま他部位に持ち込まないよう注意を促しています。
扱う範囲
-
01
対象と設計
内腔形状が正常とみなせるヒト頸動脈分岐10例。過去の CFD 研究で得た高 RRT 露出面積の十分位ごとに1例ずつ選び、正常から乱れた流れまでの幅をわざと確保しています。同一形状に剛体壁 CFD と FSI を対で適用した比較としては最大規模と著者らは述べています。
-
02
壁のモデル化
膠原線維を埋め込んだ異方性の超弾性材料で壁を表現し、拡張期のプレストレスと周囲組織の粘弾性支持まで入れています。線維方向は最大主応力に沿うという仮定から反復計算で決定。壁厚は黒血 MRI で計測した実測値（0.76〜1.07 mm）を被験者ごとに一様に与えています。
-
03
流れの条件
血液は非圧縮ニュートン流体（密度 1060 kg/m³、粘度 0.0035 Pa·s）。総頸動脈入口に PC-MRI 実測流量、内頸・外頸の出口に3要素 Windkessel。ソルバは SimVascular の svFSI（ALE 定式化）。3心拍を計算して周期収束を確認しています。
-
04
比較した量
壁面側は TAWSS・OSI に加えて WSS 位相幾何（DIVw による安定／不安定マニフォールド）と TSVI。内腔側はヘリシティ指標（h2・h4）と、周期平均の軸方向速度が負になる逆流体積の割合です。
-
05
妥当性の確認
総頸動脈断面積の心周期内最大変化は 9.1〜20.8%（平均 13.2 ± 3.4%）で、MRI 実測の 15.1 ± 5.1% と整合。径変化も超音波実測 6.6 ± 2.2% に対し 6.4 ± 1.5% でした。FSI が生理的に妥当な変形量を再現できていることの確認です。
剛体壁 CFD と FSI の差（10例の中央値）
| 項目 | 値 | 補足 |
|---|---|---|
| 低 TAWSS 露出面積の差 | 4.1% | 四分位範囲 1.63〜4.94%。3指標のなかで最も大きい |
| 高 OSI 露出面積の差 | 1.4% | 四分位範囲 1.09〜3.85% |
| 高 TSVI 露出面積の差 | 2.3% | 四分位範囲 1.69〜3.90% |
| 類似度指数（低 TAWSS） | 0.83 | 領域が空間的にどれだけ重なるか。1 が完全一致 |
| 類似度指数（高 OSI） | 0.79 | |
| 類似度指数（高 TSVI） | 0.68 | 3指標で最も低く、TSVI が壁の動きに最も敏感 |
| ヘリシティ強度 h2 の差 | 0.95 m·s⁻² | FSI のほうが概して低い |
| らせんの左右バランス h4 の差 | 0.01 | 壁の伸展性はほぼ影響しない |
| 逆流体積割合の差 | 0.96% | 四分位範囲 0.71〜1.26% |
| 断面積の心周期内変化 | 9.1〜20.8% | 総頸動脈での最大値。FSI の生理的妥当性の確認に使われた |
いずれもこの文献が報告した、同一形状に対する対比較の実測値です。四分位範囲も原著の記載によります。
実務チェック
- 頸動脈・大血管のように変形が2割程度に収まる部位では、剛体壁 CFD を選んで計算時間を優先する。
- 壁の変形が大きい構造を扱うときは、この結論をそのまま持ち込まず、対象形状での感度確認から始める。
- TSVI を主要指標にする場合は、他の WSS 指標より壁の動きに敏感であることを踏まえ、結果の解釈に幅を持たせる。
- 壁応力と血行動態の関係そのものを論じたい研究では、剛体壁では原理的に足りないため FSI を選ぶ。
- 他施設の FSI 結果と比べるときは、壁の材料モデル・壁厚・プレストレスの扱いが揃っているかを先に確認する。
この文献では分からないこと
- 対象は内腔形状が正常とみなせる頸動脈分岐のみ。プラークのある狭窄病変で剛体壁仮定が成り立つかは扱っていません。
- Windkessel の係数は剛体壁計算で調整したものを FSI にもそのまま使っており、内頸／総頸の流量配分に最大 5.2% の差が残っています。著者らはこれが一部の症例の OSI の差を説明しうると認めています。
- 壁の材料定数は被験者ごとではなく文献の組織学的データに基づく共通値です。壁厚だけが実測値で、しかも部位によらず一様と仮定されています。
- 10例のうち、分岐部の広がりが最大で心拍数と入口レイノルズ数も最高だった1例では、内腔側の指標すべてで最大の差が出ています。形状と流れ条件しだいで差が無視できなくなります。
書誌情報
| 原題 | Influence of wall distensibility on local hemodynamics at the normal carotid bifurcation: a fluid-structure interaction study |
|---|---|
| 著者 | Zambon S, Mazzi V, Calò K, Arminio M, Nocerino S, Chiastra C, Steinman DA, Morbiducci U, Gallo D. |
| 掲載誌 | Biomechanics and Modeling in Mechanobiology（2026） |
| DOI | 10.1007/s10237-026-02103-4 |
| PMID | 42479204 |
| 本文の入手 | オープンアクセス（CC BY）。全文を読んで記載しています。 |