A stabilized finite element method for a flow problem arising from 4D flow magnetic resonance imaging

技術要約：4D Flow MRIから派生する流体問題のための安定化有限要素法

問題提起
本論文は、ノイズを含む速度測定値から、圧力などの血行動態バイオマーカーを再構成する際の、4D Flow MRIデータの品質評価という課題に取り組んでいる。4D Flow MRIは、心血管系の血流を定量化するための非侵襲的なゴールドスタンダードであるが、取得時間を短縮するために用いられる積極的なアンダーサンプリング技術は、アーティファクトやノイズを導入する。さらに、標準的な再構成アルゴリズムは、流体力学の物理法則（特に質量保存則およびナビエ・ストークス方程式）を無視することが多く、圧力のような導出量に系統的な誤差をもたらす。

著者らは、真の物理的速度場 $\mathbf{u}$ を、MRIで再構成された速度 $\mathbf{u}_m$（データとして扱う）と未知のノイズ成分 $\mathbf{w}$ の和としてモデル化している。この分解を非圧縮性ナビエ・ストークス方程式に代入することで、ノイズ速度 $\mathbf{w}$ と流体圧力 $p$ を解く問題へと定式化を書き換えている。このアプローチにより、圧力再構成問題は、ノイズを含むMRIデータ $\mathbf{u}_m$ に依存するソース項を持つ境界値問題へと変換される。

手法
著者らは、この再定式化された問題の特性を扱うために設計された、安定化有限要素法（FEM）を提案している。

1. モデルの導出： ナビエ・ストークス方程式から出発し、ノイズ速度 $\mathbf{w}$ と圧力 $p$ に関する修正されたOseen型問題（線形化バージョン）を導出している。この定式化は、ノイズにより一般にゼロではない $\mathbf{u}_m$ の発散を考慮に入れている。
2. 等次補間： MRIデータの処理における実用的な効率性の必要性と性質に対応するため、本手法では速度と圧力の両方に等次有限要素空間（例：共に $P_k$）を採用している。この選択は、標準的なinf-sup（Ladyzhenskaya-Babuška-Brezzi）安定条件に抵触する。
3. 安定化戦略： 安定性を回復し、等次ペアの使用を可能にするために、著者らは「Unusual Stabilized Finite Element Methods」（具体的には文献[5]の研究を参照）に着想を得た安定化定式化を導入している。この安定化には以下が含まれる：
   • 運動方程式の残差に対する最小二乗項。
   • 対流微分の項。
   • 圧力安定化項。
   • 要素サイズ $h_T$、粘性 $\mu$、および時間ステップに相当するパラメータ $\sigma$ に依存するパラメータ $\tau_T$。
4. 線形化： 理論的解析は、線形化された問題（修正Oseen方程式）に対して行われる。非線形なケース（摂動を加えたナビエ・ストークス）については、ピカール反復を用いて数値的に対処している。

主要な貢献と理論的結果

• ウェルポーズド性（良定常性）： 著者らは、特定の条件下（特に $\sigma > 4\rho \|\nabla \mathbf{u}_m\|_\infty$ という $\sigma$ に関する制約）において、線形化された修正Oseen問題が一意の解を持つことを証明している。
• 安定性と誤差評価： 提案された安定化有限要素スキームについて、以下の事項を確立している：
  • 強圧性（Coercivity）： 双線形形式が特定のエネルギーノルムにおいて安定性評価を満たすこと。
  • 最適収束性： 速度と圧力の両方について、最適な次数の誤差評価を導出している。誤差境界は、メッシュサイズ $h$、解の正則性、および対流速度場の近似精度に依存する。
  • 整合誤差（Consistency Error）： 本スキームは、非線形な文脈における対流項の近似により、厳密には整合的ではないことを認めているが、この整合誤差が最適に抑えられることを示している。
• 圧力再構成： 本手法は、ノイズと圧力を同時に推定する逆問題を解くことにより、非侵襲的な圧力再構成の手順を提供する。

数値実験
著者らは、FEniCSライブラリを用いて3つの数値実験を行い、手法を検証している。

1. Kovasznay Flow（学術的テスト）： 解析解を用いて、様々な粘性値において、線形（$k=1$）、二次（$k=2$）、および三次（$k=3$）要素で最適な収束率を示すことを実証している。これらの結果は、線形モデルと非線形摂動ナビエ・ストークスの両方で成立する。
2. 屈曲領域上のランダム速度： 四角形メッシュ上で区分定数速度データを用いた4D Flow MRIの設定をシミュレートし、本手法が圧力場を正常に再構成することを示している。一部のぼけ（区分定数データと連続有限要素間の補間誤差に起因）は見られるものの、圧力の回復は堅牢である。
3. 補間されたナビエ・ストークス速度： 解像度の低いナビエ・ストークスシミュレーションから生成され、補間された速度データを、圧力を再構成するために使用する現実的なシナリオ。再構成された圧力および修正された速度場は、高解像度の参照解と高い一致を示しており、特定の断面における最大誤差は約1.6%である。

意義と主張
本論文は、提案された安定化有限要素法が、以下のための堅牢な数学的枠組みを提供することを主張している：

• 品質評価： 「ノイズ」成分 $\mathbf{w}$ を定量化することによる、測定された速度場の品質推定。
• 圧力再構成： 純粋に現象論的な圧力再構成手法を回避し、侵襲的な処置を避けるための、物理に基づいた圧力再構成の代替手段の提供。
• 効率性： 複雑な3次元流体問題において極めて重要となる、安定性を損なうことなく計算効率の高い等次有限要素ペアの使用を可能にすること。

著者らは、線形化されたケースについては完全に解析されているものの、完全な時間依存問題への拡張や、圧力ロバストまたは発散フリーな離散化の探索は、今後の研究における未解決の課題であると述べている。本研究は、臨床的な4D Flow MRIデータからより信頼性の高い血行動態バイオマーカーを算出するための、一歩となるものである。