公開対象 / 安全性

本記事は MRI / 4D Flow / 統計 の教育目的で作成された一般化された解説です。特定患者、施設、製品仕様、非公開 WIP の代替ではありません。

まず一言で

Vorticity (渦度, ω = ∇ × v) は「流体要素が局所的にどれだけ回転しているか」を表すベクトル場で、単位は 1/s。直感的には「流体の小さなコマを置いたとき、そのコマが回る角速度の 2 倍」である。Helicity (ヘリシティ) はその vorticity と速度ベクトルの内積 h = v · ω で、流れが「直進」しているか「らせん状にねじれて進む」かを定量化する。

4D Flow MRI が画像化する量は voxel 単位の 3 方向速度 v(x,t) だけだが、その空間微分から vorticity / helicity を後処理で導出できる。これは「流速の大きさ」「圧較差」「壁面せん断応力」とは独立に、流れの三次元構造そのもの (回転・ねじれ) を読み取る情報軸を提供する。健常大動脈の生理的 right-handed helical flow二尖弁大動脈 の異常 helical pattern、左室拡張期 vortex の収縮機能との関連、Fontan 循環 での energy loss 増大 — これらはすべて vorticity / helicity を読まなければ捉えられない流れの「形」である (Kilner 1993, Markl 2004, Hope 2007, Bissell 2013, Sengupta 2012)。

ただし — そして本記事の核心はここにある — vorticity は velocity の空間 1 階微分 (curl) であり、WSS と同じく 4D Flow の voxel size・SNR・smoothing に極度に敏感である。|ω| 単独で渦を同定すると shear flow (層流のせん断) を渦として誤検出する。これを避けるため Q-criterion (Hunt et al. 1988) や λ2 criterion (Jeong & Hussain 1995) という古典的 vortex identification 法が必要になる。Helicity は v · ω なので vorticity の誤差を相続し、さらに「左右どちらに巻くか」の符号情報を持つため軽微な vorticity ノイズで符号反転する voxel が頻発する。「数式は単純、解釈は要注意」が本領域でも繰り返し確認されているコンセンサスである (Markl 2012 review; Dyverfeldt 2015 consensus; Soulat 2020 update)。

---

1. Vorticity (渦度) — 速度場の curl

Vorticity ω = ∇ × v の 3D ベクトル幾何と流体要素の局所角速度との対応
図1. Vorticity 定義 — ∇ × v の 3D ベクトル幾何 — vorticity ベクトルの curl 定義と「流体要素の角速度の 2 倍」直感

1.1 ベクトル定義

連続体力学における vorticity は速度場 v(x,t) の curl (回転) として定義される:

ω = ∇ × v,    [ω] = 1/s

成分表示すると 3 成分のベクトル場である:

ω_x = ∂v_z/∂y − ∂v_y/∂z
ω_y = ∂v_x/∂z − ∂v_z/∂x
ω_z = ∂v_y/∂x − ∂v_x/∂y

各成分は「ある軸まわりの局所的回転」を表し、ω ベクトル全体の方向は 回転軸の向き (右手系)、大きさ |ω|回転の強さ に対応する。

1.2 物理的意味 — 「流体要素の角速度の 2 倍」

連続体力学の古典的結果として、流体中の微小な剛体球を浮かべたときの球の自転角速度 Ω は局所 vorticity の半分になる:

Ω_local = (1/2) · ω

この事実から「vorticity は流体要素が回転する角速度の 2 倍」と直感的に呼ばれる。|ω| = 2 1/s ならば、その voxel に置いた仮想コマは 1 rad/s = 約 9.5 rpm で自転している。

1.3 スカラー指標 — vorticity magnitude

しばしば実用上はベクトル ω 全体ではなく大きさ |ω| (vorticity magnitude) を 3D 可視化に使う:

|ω| = √(ω_x² + ω_y² + ω_z²)

これは「方向を捨てて回転強度のみを 1 スカラーで読む」操作で、大動脈・心房・心室の vortex 可視化で標準的に用いられる。ただし §4 で詳述するように、|ω| 単独は shear flow を渦として誤検出するので、vortex 同定には Q-criterion / λ2 criterion を併用する。

1.4 2D vs 3D vorticity の違い

教科書では渦と言えば「平面内の回転」(2D vorticity, scalar) で導入されることが多いが、生体血管は本質的に 3D で、2D 投影では捉えられない以下の構造を持つ:

  • Helical flow: 進行方向と回転軸が並行・反並行する、らせん状の流れ。2D 断面では「中心を通る速度」と「周辺の渦」が直交して見えるが、3D で見れば 1 つのコルク栓抜き状の流れ
  • Vortex ring: トーラス状の渦輪 (心室の拡張期 inflow vortex が典型)。2D 平面切断では片側で +ω、対側で −ω に見えるが 3D ではつながった輪
  • Twisted vortex / vortex breakdown: 大動脈 jet が下流で複数の小渦に崩壊する非定常構造

4D Flow MRI は 3D vorticity を直接得られる稀な modality であり、CFD や echocardiography (2D vorticity 主体) に対する強みである (Markl 2012 review)。

---

2. Helicity — 速度と vorticity の内積

Right-handed と left-handed helical flow の比較と LNH 符号の対応
図2. Helicity — right-handed vs left-handed helix — LNH の符号と helical 流れの方向を視覚化

2.1 局所定義

ヘリシティ密度 は速度ベクトルと vorticity ベクトルの内積として定義される:

h = v · ω = v_x · ω_x + v_y · ω_y + v_z · ω_z
[h] = (m/s) × (1/s) = m/s²

物理的には「速度ベクトルと回転軸ベクトルがどれだけ並行しているか」を回転強度で重み付けした量。v ∥ ω のとき最大、v ⊥ ω のとき 0、v ∥ −ω のとき最小 (負の最大)。

2.2 体積積分形 — Helicity total

流れ場全体の helicity は密度の体積積分:

H = ∫∫∫_V (v · ω) dV
[H] = m⁴/s²

これは Moffatt 1969 が示した topological invariant (流体力学的トポロジー不変量) で、理想流体 (粘性ゼロ + 非圧縮) では「絡まった渦の絡まり度」を表し時間保存する量。実際の血流は粘性があるため厳密には保存しないが、「全体としてどれだけねじれているか」の summary 指標として臨床研究で用いられる。

2.3 Normalized helicity — LNH

helicity 密度は単位が m/s² で、流速と vorticity の絶対値に依存するため、施設間 / 患者間の比較が難しい。これを 無次元化 したものが Local Normalized Helicity (LNH, 局所正規化ヘリシティ):

LNH = (v · ω) / (|v| · |ω|),    LNH ∈ [-1, +1]

幾何学的には vω のなす角 θ の cosine。

LNH流れの形物理
+1right-handed helix (右巻きらせん)v と ω が同方向 (進行方向に右手系で回りながら進む)
0純粋並進 or 純粋せん断v と ω が直交 (回転軸が進行方向と垂直)
−1left-handed helix (左巻きらせん)v と ω が反対方向 (進行方向に左手系で回りながら進む)

LNH = ±1 は「コルク栓抜きのように」 進むらせん流。|LNH| が 0.4–0.7 で「中程度の helical 性」、0.1 未満で「ほぼ並進」と臨床的に解釈される (Frydrychowicz 2009, Lorenz 2014)。

2.4 直感 — コルク栓抜きアナロジー

helicity を直感するには「コルク栓抜きを上から押し下げる動作」を考える。手は栓抜きを下方向に押し (進行速度 v) かつ右回りに回している (右手系の角速度 → 上向きの vorticity ω)。v が下向き、ω が上向きなので内積は負。これが left-handed helix。

逆に、右回り (上向き ω) のまま上向きに引き上げる動作 (v 上向き) なら、v · ω > 0 で right-handed helix。

血流の場合: 大動脈で「左心室 outflow tract から ascending aorta に至る間に流れが右回りに捻れる」のが正常 (Kilner 1993)。これは進行方向の血流が右巻きに回る = LNH > 0 (right-handed) を意味する。§7 で詳述する通り、BAV ではこの helical pattern が左右の偏りや弱化・反転を示す (Bissell 2013)。

---

3. 大動脈の生理的 Helical Flow — なぜ正常でも血流はねじれるか

健常大動脈の収縮期 right-handed helical flow pattern (Kilner 1993 系)
図4. 大動脈の正常 right-handed helical flow (Kilner pattern) — 健常者の収縮期 ascending aorta + arch の生理的らせん flow

3.1 Kilner (1993) の発見

Kilner 1993 は当時画期的だった 3D phase contrast MRI で正常ヒト大動脈の血流軌跡を可視化し、収縮期中期から後期にかけて ascending aorta から arch にかけて right-handed helical flow が一過性に出現する ことを示した (*Nature* 1993)。これは病的所見ではなく 健常者の生理的パターン で、その後 100% の正常ボランティアで確認されている (Markl 2004 Circulation; Hope 2007 Radiology)。

3.2 なぜ生理的に helical になるか — 3 つの寄与

1. LV outflow tract の geometry: 左室は球状から outflow tract で捻れた円筒に移行する。LV 内血流の rotation 慣性 (収縮期 vortex の残り) が outflow を通って ascending aorta に伝達される

2. Aortic root の twisted geometry: 大動脈基部は心室から見て約 30–60° 捻れた幾何で、流れに helicity を与える

3. Aortic arch の curvature: arch の S 字状湾曲 + Dean vortex (湾曲管特有の二次流れ) が helical pattern を維持

これら 3 つの寄与により、健常者の ascending aorta では収縮期 peak 前後で LNH = 0.3–0.6 程度の中程度 right-handed helicity が観察される (Frydrychowicz 2009 Eur Radiol; Lorenz 2014 MRM)。

3.3 生理学的意味

helical flow には複数の仮説的役割が提案されている (確立順):

  • Wall washout: らせん状の流れは壁近傍に淀みが生じにくく、内皮環境を均質に保つ (Caro 1996 系の仮説)
  • Energy efficiency: 高 Reynolds 数下で 乱流遷移を抑制する (層流の wider stable range) 効果
  • Atheroprotection: helical pattern と低 OSI / 高 TAWSS 領域が地理的に一致する報告あり (Liu 2009 J Biomech 系)

これらは「helical = 健康な流れ」という単純な等式ではなく、「健常では適切な helicity が維持されているが、病的では失われたり過剰になったり方向反転する」という枠組みで捉えるのが現代的解釈である (Soulat 2020 update review)。

---

4. Vortex Identification — |ω| だけでは渦は同定できない

Vorticity magnitude、Q-criterion、λ2 criterion による vortex 同定の比較
図3. Vortex identification — Q-criterion vs λ2 vs |ω| 比較 — shear flow と vortex の区別、3 手法の同定結果の違い

4.1 問題提起 — Shear flow と vortex の区別

流体力学で歴史的に厄介な問題は 「渦 (vortex) と純粋せん断 (shear) をどう区別するか」 である。例えば壁近傍の Poiseuille 流は中心軸近傍を除いて常に有限の |ω| を持つが、これは「壁付近で流速が急変する shear」であって渦巻きではない。|ω| の閾値で渦を同定すると、壁近傍の薄い shear layer も含めて渦判定されてしまう。

この問題に対し、流体力学コミュニティは 1980 年代後半から複数の vortex identification criterion を提案してきた。代表的な 2 つを以下に示す。

4.2 Q-criterion (Hunt, Wray, Moin 1988)

速度勾配テンソル ∇v を対称部分 (歪み速度テンソル) と反対称部分 (回転テンソル) に分解する:

∇v = S + Ω
S   = (1/2) · (∇v + ∇vᵀ)   [strain rate tensor]
Ω   = (1/2) · (∇v − ∇vᵀ)   [rotation tensor]

Q は両者の Frobenius ノルムの 2 乗の差を半分にした量:

Q = (1/2) · (|Ω|² − |S|²),    [Q] = 1/s²

Q > 0 の領域を渦とみなす。すなわち「局所的に回転が歪みより支配的な領域」を vortex として同定する。

直感: |Ω|² は vorticity の 2 乗に比例し「回転の強さ」、|S|² はせん断歪みの 2 乗に比例し「形を変える強さ」を表す。Q > 0 は「回転 > せん断」、つまり流体要素が変形より回転で動いている領域。これは shear flow を除外し、純粋な渦のみを抽出する。

実用: 4D Flow MRI では Q-criterion による iso-surface 描画 (e.g., Q = 100 1/s² など閾値) が ascending aorta や LV inflow vortex の可視化に標準的に用いられる (Markl 2012 review)。

4.3 λ2 criterion (Jeong & Hussain 1995)

別アプローチ。圧力 Hessian ∂²p/∂x_i∂x_j の固有値解析を行う。Navier-Stokes 方程式の対称部分から:

S² + Ω² = − (1/ρ) · ∇∇p + (粘性項)

左辺 S² + Ω² の 3 固有値を λ_1 ≥ λ_2 ≥ λ_3 と並べたとき、λ_2 < 0 の領域を渦とみなす

直感: vortex core の中心では「圧力極小」が形成され (遠心力との釣り合い)、圧力 Hessian の中間固有値が負になることを利用する。

Q-criterion との違い:

  • Q-criterion は局所的なテンソル不変量 (より計算が単純)
  • λ2 criterion は圧力 (大域情報) を間接的に参照、複雑な flow geometry で robust と報告される (Jeong 1995 J Fluid Mech)

両者は多くの場合似た vortex 同定結果を与えるが、極端な非定常流や強い shear layer 近傍では分岐する。4D Flow MRI 文献では Q-criterion がやや優勢、CFD 文献では λ2 criterion 併用が多い印象。

4.4 |ω| 単独 vs Q vs λ2 — 比較表

手法数式渦判定shear flow を除外実装難度
Vorticity magnitude `ω`√(ω_x² + ω_y² + ω_z²)閾値以上不可 (壁近傍 shear も拾う)最低
Q-criterionQ = (1/2)(Ω² −S²)Q > 0
λ2 criterionλ2(S² + Ω²)λ2 < 0

|ω| 単独表示は誤検出を生む。臨床研究で「vorticity field」を扱う場合は、Q-criterion か λ2 criterion で「渦領域」をマスクしてから可視化するのが standard practice (Markl 2012, Dyverfeldt 2015 consensus)。

4.5 Q-criterion の cancellation 利点

|ω| はベクトル長なので正値しか取らない。一方 Q-criterion は |Ω|² − |S|²正負の符号を持つ scalar field。これにより:

  • 隣接 voxel で正の Q と負の Q が出ても (ノイズで) 平均化されてゼロに近づく → ノイズに stable
  • shear-dominated region では Q < 0 となり、自動的に渦領域から外れる

|ω| は単純に大きさを足すだけなのでノイズが累積するが、Q は cancellation で stable。これも Q-criterion が 4D Flow 解析で好まれる理由の 1 つ。

---

5. 4D Flow MRI による Vorticity / Helicity 推定の落とし穴

4D Flow MRI による vorticity 推定および vortex identification の 5 ステップパイプライン
図6. 4D Flow MRI から vorticity field 推定パイプライン — velocity → smoothing → 空間微分 → vorticity → vortex identification の流れ
LNH map の voxel-wise 符号不安定性と時間平均 / ROI 積分による安定化
図7. Helicity 符号の voxel-wise 不安定性 — ノイズによる LNH 反転 — LNH の voxel 単位 map がノイズで ±1 にランダム反転する artifact の可視化

vorticity は velocity の curl (1 階空間微分) であり、WSS と同じ「微分量の SNR 問題」を抱える。helicity は v · ω なので vorticity の誤差を相続する。

5.1 微分操作によるノイズ増幅

velocity ノイズを σ_v、voxel size を Δx とすると、有限差分で計算される vorticity のノイズは概ね:

σ_ω ~ σ_v / Δx

例: σ_v = 5 cm/s = 0.05 m/s, Δx = 2 mm = 0.002 m のとき σ_ω ~ 25 1/s。生理的 vorticity が 10–50 1/s オーダーなので、SNR は 0.4–2.0 と低い。これは「壁近傍 WSS」と同じ困難で、smoothing なしの直接微分はノイズに埋もれる。

実装上は velocity field を 3D B-spline / 2 次多項式で平滑近似してから解析微分する手法 (Stalder 2008 protocol を vorticity に転用) が多い。Smoothing parameter が結果を支配的に変える点に注意。

5.2 Voxel size 依存性 — 小渦の under-detection

vorticity を取る voxel が大きすぎると、その voxel 内の 正逆 vorticity が打ち消しあい、見かけ上 |ω| が低く出る。特に boundary layer 内の薄い shear vortex や、流速変化の急な領域では voxel 1 個に複数の vortex 構造が混入し、平均化されて vorticity が消える。

Markl 2012 review および Cibis 2015 J Biomech (Fontan 領域) は、voxel size 2–3 mm の典型 4D Flow で 直径 4 mm 以下の vortex は再現性をもって検出できない ことを示している。心室 inflow vortex (直径 20–40 mm) は問題ないが、aortic arch の secondary vortex (直径 5–10 mm) は voxel size に応じて見えたり見えなかったりする。

5.3 時間分解能 — Transient vortex を見逃す

vortex は心拍内で出現・消失する非定常構造である。特に LV inflow vortex は拡張早期 (E wave) と心房収縮期 (A wave) で異なる強度・位置を取る。時相数 (cardiac phase) が 20 未満だと vortex 寿命の半周期以下しかサンプルできず、ピーク強度を見逃す

Dyverfeldt 2015 consensus は 4D Flow MRI で vortex / helicity 解析を行う場合、心拍 1 周期で 25 phase 以上 を推奨している (典型 60 bpm で 1 phase ≈ 40 ms)。それでも fast transient (10 ms オーダー) は捉えられない。

5.4 Helicity 固有の符号不安定性

LNH = (v · ω) / (|v| · |ω|) は分子の符号が vω の相対方向で決まる。vorticity がノイズで方向反転する voxel が頻発するため、LNH map が ±1 の間でランダムに振動する artifact が出やすい。

実用上の対策:

  • LNH を時間平均 (cardiac cycle 全体) で評価し、瞬時値の符号反転を平滑化
  • |ω| < threshold の voxel を mask out (vorticity が弱い領域での LNH は意味なし)
  • 体積積分 helicity H = ∫ v · ω dV を ROI 単位で評価し、voxel-wise の LNH map ではなく integrated metric で議論

これらの対策は OSI の denominator artifact 対策 (参考: OSI denominator artifact) と本質的に同じ「弱い signal voxel での比は信用しない」原則に従う。

5.5 VENC / VNR の重要性

vorticity は velocity の差分なので、低流速領域 (心房 / Fontan / LV apex / 瘤内) で VNR が低いと vorticity がノイズに支配されるMulti-VENCMulti-VENC 完全解説 で議論される VNR 改善は vorticity / helicity の精度を直接律速する。

特に LNH の符号情報 は VNR < 10 では信頼性が低く、Fontan や stagnation 領域では Multi-VENC が事実上必須となる (Bock 2010, Soulat 2020 review)。

5.6 Wall motion の無視

大動脈・心室は心拍で 5–15% 径変化する。多くの vorticity 解析は壁を時間固定とみなし、結果として wall-region の vorticity が見かけ上 fluctuating に出ることがある。動的 segmentation を組み合わせた解析 (Bock 2010 系) は研究レベルで存在するが計算コストが高い。

---

6. 派生指標 — Helicity / Vorticity と他の hemodynamic indices

vorticity / helicity は流れ場全体の指標で、壁ローカル指標 (WSS, OSI, RRT) と相補的位置づけ。

6.1 Helicity と Energy Loss

理論的には helicity が高い (helical flow) と viscous dissipation (粘性散逸) が増えるという主張と、減るという主張の両方がある。

  • 増える説: 3D の helical motion は単純並進より速度勾配が複雑で、Φ = μ · (∇v + ∇vᵀ) : ∇v が増える
  • 減る説: 適度な helicity は乱流遷移を抑制し、結果として乱流性散逸が減る

実測 (4D Flow MRI) では Fontan circulation や BAV 大動脈で H 増 ↔ Energy Loss 増 の正相関が観察される (Hsiao 2009, Stankovic 2014)。これは「過剰な helicity = 効率の悪い流れ」の解釈と整合するが、健常範囲内の生理的 helicity については議論継続中。

6.2 Vorticity と WSS

壁直近の vorticity 成分 ω_t (壁面接線成分) は、wall shear rate (壁面せん断率) と直接比例する:

wall_shear_rate ≈ |ω_t|_wall
WSS_magnitude  ≈ μ · |ω_t|_wall

つまり「壁近傍 vorticity の接線成分」は WSS そのものに近い情報を持つ。逆に言えば、4D Flow で WSS を推定するときの誤差源 (boundary layer 過小サンプル) はそのまま wall vorticity 推定の誤差源にもなる。

6.3 Helicity と OSI / RRT の同時上昇

頚動脈分岐部・大動脈瘤・BAV jet 領域などの disturbed flow region では、低 TAWSS + 高 OSI + 高 RRT の「動脈硬化好発指紋」とともに異常 helicity pattern (反転 / 弱化) が観察される。これらを組み合わせて「disturbed flow signature」として multi-axis に評価する研究が増えている (Soulat 2020 update)。

6.4 4 指標 + helicity の使い分け表

指標単位何を読むか主たる臨床応用
WSS (TAWSS)Pa壁を擦る平均強度BAV aortopathy、低 WSS atherogenesis
OSI無次元WSS 方向振動性動脈硬化好発部位同定
RRT1/Pa滞留性古典的好発部位 summary
Vorticity `ω`1/s局所回転強度vortex 可視化 (要 Q-criterion 併用)
Helicity H / LNHm⁴/s² / 無次元流れのねじれBAV pattern, helical wash, Fontan
Energy LossmW (積分)流れ場全体の粘性散逸弁逆流効率、Fontan

「壁を読む」 (WSS/OSI/RRT) と「流れ場全体を読む」 (vorticity/helicity/EL) は補完関係。両者を組み合わせて初めて hemodynamic 環境が記述できる。

---

7. 臨床応用 — 現状エビデンス

BAV cusp fusion 形態別 (R-L vs R-N) の異常 helical pattern と jet 偏倚の模式図
図5. BAV における異常 helical pattern (RL vs RN fusion) — 弁形態別 jet 偏倚と helical pattern 異常
Helicity / Vorticity の 5 主要臨床応用領域 (BAV / LV vortex / PAH / Fontan / LAA) のマップ
図8. Helicity / Vorticity の臨床応用マップ — BAV / LV vortex / PAH / Fontan / Atria の応用領域を 1 枚に整理

7.1 BAV (二尖弁大動脈) と ascending aortopathy

最も研究蓄積が厚い領域。BAV では弁形態 (right-left fusion vs right-non fusion) に応じて、ascending aorta 内の helical pattern が健常 right-handed pattern から偏移する。

Hope 2007 (*Radiology*) は 4D Flow MRI で正常 / BAV / Marfan の helical flow pattern を比較し、BAV で right-handed helical flow の amplification と局所反転 を初めて報告した。Bissell 2013 (*Circ CV Imaging*) はより大規模コホートで BAV cusp fusion 形態と helical pattern の関係 を定量化し、right-left fusion で右側 ascending aorta に jet が偏倚、right-non fusion で前壁に偏倚することを示した。

Mahadevia 2014 (*Circulation*) と Guzzardi 2015 (*JACC*) は同じ BAV cohort で WSS pattern + helical pattern の両方を解析し、異常 helical flow region と局所 WSS 異常 region が地理的に一致 することを示した。これは「helical pattern も大動脈リモデリングの hemodynamic 駆動因子」という仮説の根拠となっている。

7.2 LV (左室) inflow vortex

Sengupta 2012 (*JACC Imaging*) は echocardiography の vortex 解析を 4D Flow MRI に拡張し、拡張期早期に出現する LV vortex の寿命・強度・位置が LV function と相関 することを示した。具体的には:

  • 健常者: 拡張期早期 vortex が安定して形成、収縮期に解離
  • 拡張機能不全 (e.g., HFpEF): vortex 強度低下、位置が apex 寄りに偏倚
  • 拡張型心筋症 (DCM): vortex が disorganized、収縮期にも残存

vortex の Q-criterion ベース可視化、vortex 中心 trajectory、vortex 体積積分 helicity を用いた評価が行われている。臨床的価値は研究段階だが、HFpEF の早期診断 marker として注目されている (Soulat 2020 update)。

7.3 Pulmonary artery (PA) — PAH

肺高血圧症 (PAH) では主肺動脈 (MPA) 拡張に伴い helical / vortex pattern が顕著化する。Geiger 2011 (*Eur Radiol*) は PAH 患者で 4D Flow MRI による MPA helical flow を評価し、helical flow の強度と pulmonary vascular resistance (PVR) が相関 することを示した。PAH 治療効果 (例: prostacyclin) の longitudinal evaluation に応用研究が進んでいる。

7.4 Fontan circulation — energy loss と helicity

Hsiao 2009 (*Radiology*) は Fontan 循環で 4D Flow MRI による helical flow および energy loss を評価し、post-surgical helical pattern が energy loss と強く相関 することを示した。Fontan 外科手技 (extracardiac vs intracardiac conduit, lateral tunnel) の選択で hemodynamic outcome が異なり、helicity / energy loss を術後評価で用いる試みが行われている。

ただし Fontan の流速は 5–20 cm/s と低く、§5.5 で議論した通り Multi-VENC + VNR 改善が前提条件となる。

7.5 Atria — LAA 内 helical flow と stasis

左心耳 (LAA) は心房細動患者で thrombus 形成の主要部位。LAA 内の血流は通常 helical pattern を持ち、AF や LAA morphology (chicken wing vs cactus vs windsock vs cauliflower) で pattern が変化する (Stankovic 2014 JCMR 系)。Helicity 低下と stagnation の関連が研究されているが、LAA は小構造のため 4D Flow voxel size の制約が厳しく、研究段階に留まる。

7.6 Cerebral aneurysm

脳動脈瘤内の vortex 構造 (single vortex vs multiple vortices, stable vs unstable) と破裂リスクの関連が CFD 主体で研究され、4D Flow MRI でも検証が始まっている。Vortex の数、不安定性、helicity の正負交代頻度などが評価指標。臨床利用には程遠いが、術前計画への活用が研究段階で進む。

---

8. 限界・落とし穴 (チェックリスト)

vorticity / helicity を含む 4D Flow 研究を読む / 書くときの最低限チェック:

  • VENC / VNR: 関心領域の velocity に対する VENC マージン、Multi-VENC 使用有無 (helicity 符号情報は VNR > 10 が事実上下限)
  • Voxel size: ≤ 2 mm 等方性が望ましい (大動脈)、心室・Fontan・LAA はさらに小
  • 時間分解能: cardiac phase ≥ 25 (60 bpm で 1 phase ≈ 40 ms)、transient vortex 評価には 30+ phase 推奨
  • Vortex identification 法: |ω| 単独か、Q-criterion か、λ2 criterion か (|ω| 単独は shear flow を誤検出)
  • Helicity 正規化: 絶対 helicity か LNH か (施設間比較には LNH 必須)
  • Smoothing parameter: velocity field の事前平滑化 (B-spline order, kernel size) が結果を支配的に変える
  • Wall motion modeling: 静止壁仮定か、動的 segmentation か
  • 時間積分の有無: 瞬時 vorticity / helicity 評価か、cardiac cycle 平均か
  • Mask 閾値: 弱 vorticity voxel の LNH をどう扱うか (mask out か、表示するか)
  • 比較対象: 絶対値ではなく群間差・部位間差・atlas 偏差で議論しているか

「絶対値は ±30–50% 動く可能性がある + helicity 符号は voxel-wise には信頼度低い」という前提を共有した上で、relative / integrated metric として使うのが現時点のベストプラクティスである (Markl 2012, Dyverfeldt 2015 consensus, Soulat 2020 update)。

---

9. ZK 化と次の記事への橋

この記事は次の atomic note に分解する (詳細は zk_notes/):

  • Vorticity 定義 — curl of velocity と物理直感 (260516_vorticity_curl_velocity.md)
  • Helicity 密度 — v · ω 内積と normalized helicity (260516_helicity_density_local_definition.md)
  • Vortex identification — Q-criterion vs λ2 vs |ω| (260516_q_criterion_lambda2_vortex_identification.md)
  • 大動脈生理的 helical flow — Kilner 1993 起源と right-handed pattern (260516_helical_flow_aorta_physiology.md)
  • 4D Flow vorticity 推定の限界 — 微分量の SNR、voxel/時相依存 (260516_4dflow_vorticity_estimation_limits.md)
  • BAV helicity pattern — Hope 2007 / Bissell 2013 系 (260516_bav_helicity_pattern.md)
  • Helicity の臨床的意義 — LV/PA/Fontan/Atria の応用 summary (260516_helicity_clinical_significance.md)
  • MOC: 上記 7 atomic を orchestrate する Map of Content (260516_helicity_vorticity_moc.md)

連続記事の予定:

  • WSS パック — 兄弟記事 (壁ローカル指標)
  • OSI / RRT パック — 振動性 / 滞留性
  • Multi-VENC パック — vorticity / helicity 精度の前提
  • Energy_Loss_4DFlow — helicity と energy loss の関連
  • Bland-Altman / ICC パック — helicity の test-retest reproducibility

---

10. 公開対象 / 安全性

  • 施設名、患者情報、内部契約情報、メール、個人情報は含まない
  • 製品名 (装置・WIP) は含まない、引用ソフトは公開既知のもののみ
  • 一次文献の DOI / PubMed は references.md に集約
  • 画像生成プロンプトは image_prompts.json に分離 (Codex 側で GPT Image2 実行)
  • HTML 化、Cloudflare deploy、Discord 通知は Codex 担当

---

Disclaimer (推奨掲載)

本記事は教育目的の解説であり、特定の臨床判断・治療方針を推奨するものではない。vorticity / helicity は依然として研究段階の指標が中心で、ガイドライン推奨は限定的である。臨床応用は施設の倫理審査・専門医判断に基づくこと。引用文献の DOI / PubMed は記事末尾の references.md に集約してある。