[論文レビュー] A Fast Volume Integral Equation Solver with Linear Basis Functions for the Accurate Computation of Electromagnetic Fields in MRI

Ioannis P. Georgakis, Ilias I. Giannakopoulos+2 more|arXiv (Cornell University)|Feb 6, 2019

Microwave Imaging and Scattering Analysis参考文献 68被引用数 9

ひとこと要約

本稿では、高対比で不均一なMRI環境における電磁界（EM）計算の高精度を実現するため、不連続な区分線形（PWL）基底関数を用いた高速で安定した体積積分方程式（VIE）ソルバを提案する。ガラーキン法のモーメント法とFFTを用いた行列・ベクトル積、および新規のプリコンディショナを活用することで、特に粗いメッシュ解像度下でも優れた精度と収束性を達成し、自由度を減らした状態で従来の区分定数（PWC）スキームを上回る性能を発揮する。

ABSTRACT

A stable volume integral equation (VIE) solver based on polarization/magnetization currents is presented, for the accurate and efficient computation of the electromagnetic scattering from highly inhomogeneous and high contrast objects.We employ the Galerkin Method of Moments to discretize the formulation with discontinuous piecewise linear basis functions on uniform voxelized grids, allowing for the acceleration of the associated matrix-vector products in an iterative solver, with the help of FFT. Numerical results illustrate the superior accuracy and more stable convergence properties of the proposed framework, when compared against standard low order (piecewise constant) discretization schemes and a more conventional VIE formulation based on electric flux densities. Finally, the developed solver is applied to analyze complex geometries, including realistic human body models, typically used in modeling the interactions between electromagnetic waves and biological tissue.

研究の動機と目的

• 高磁場MRIに特徴的な高不均一性・高対比な生物学的組織における、正確で効率的なEM界計算の課題に対処すること。
• 特に粗いメッシュ解像度下で精度と収束安定性に欠ける低次元（区分定数）基底関数の限界を克服すること。
• メッシュの細分化がしばしば不実用的である超高磁場MRIにおける、患者個別化SAR計算およびRFコイルモデル化に適した安定的で高速なVIEソルバを開発すること。
• 計算コストや自由度の増加なしに精度を向上させるためのh-Refinementの代替手段として、p-Refinement（高次元基底関数）を可能とすること。
• 既存のSIEベースのフレームワークに統合し、現実的な人体モデル（RHBMs）におけるEM相互作用のハイブリッドモデリングを可能とすること。

提案手法

• 等価な極化電流および磁化電流を用いてEM散乱問題を定式化し、安定な体積積分方程式（VIE）の定式化を可能にする。
• 均一なボクセルグリッド上に定義された不連続な区分線形（PWL）基底関数およびテスト関数を用いてVIEを離散化し、区分定数関数よりも優れた場の表現を実現する。
• ガラーキン法のモーメント法を適用することで、ガラーキン直交性を確保し、数値的安定性および精度を向上させる。
• 繰り返し解法（例：GMRES）における行列・ベクトル積を高速フーリエ変換（FFT）を用いて加速し、大規模問題における効率的計算を実現する。
• 特に細かい解像度下で未処理系が発散する状況を考慮し、繰り返し解法の安定化を図る専用プリコンディショナを実装する。
• FFTを用いた行列・ベクトル乗算の効率化と、既存のMRIモデリングパイプラインとの互換性を維持するため、均一なボクセルグリッドを活用する。

実験結果

リサーチクエスチョン

• RQ1 不連続な区分線形（PWL）基底関数は、標準的な区分定数（PWC）関数と比較して、高対比・不均一なMRI問題におけるVIEソルバの精度と収束安定性を顕著に向上させるか？
• RQ2 PWL基底関数によるp-Refinementは、PWC関数を用いたh-Refinementと同等の精度を達成しつつ、自由度の増加を抑えることができるか？
• RQ3 PWL関数を用いた現在基準のVIE定式化は、極めて細かい解像度（例：1 mm）でも数値的安定性と収束性を維持できるか、発散を回避できるか？
• RQ4 提案されたプリコンディショナは、高誘電率のシェーマンパッドのような困難な状況下でも、PWLベースのVIE系の繰り返し解法の収束をどれほど効果的に加速できるか？
• RQ5 PWLベースのVIEソルバは、粗い解像度下でも現実的人体モデル（RHBMs）におけるEM界および磁界エネルギーを信頼性高く計算でき、実用的な患者個別化SAR解析を可能にするか？

主な発見

• PWLベースのVIEソルバは、粗い解像度下でもPWCベースのソルバを上回る高い精度を達成する。2 mmのPWLにおける磁界エネルギーは1 mmのPWCを下回り、誤差が小さいことを示している。
• 1 mm PWCの自由度の半分しか持たない2 mmのPWLでも、磁界エネルギー（error_h1）、磁界エネルギー密度（error_h2）、磁界ノルム（error_h3）の誤差が1 mm PWCを下回る。
• 5 mmのPWL解は、自由度が約3.75倍多い2 mm PWCと同等またはそれ以上の精度を達成する。
• 未処理のPWLソルバは1 mm解像度で発散するが、提案されたプリコンディショナにより収束が安定化され、反復回数が3900回から700回に減少し、PWCソルバと同等の性能を達成する。
• PWLソルバはプリコンディショナを用いることで安定的かつ効率的に収束し、メッシュの細分化が不実用な高負荷なMRIシミュレーションにおいても実用的である。
• 提案されたフレームワークは、高誘電率のシェーマンパッドを含む複雑なRHBMsにおいても、高対比な状況下で頑健に動作し、正確なEM界計算を可能にする。