[논문 리뷰] 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명|arXiv (Cornell University)|2019. 02. 06.

Microwave Imaging and Scattering Analysis참고 문헌 68인용 수 9

한 줄 요약

이 논문은 고대비, 비균질적인 MRI 환경에서 정확한 전자기(EM)장 계산을 위한 불연속(piecewise) 선형(PWL) 기저 함수를 사용한 빠르고 안정적인 체적 적분방정식(VIE) 해법을 제시한다. FFT 가속화된 행렬-벡터 곱과 새로운 조절자(preconditioner)를 활용한 갈레르킨 방법의 모멘트(MoM)를 통해, 특히 굵은 해상도에서 뛰어난 정확도와 수렴성 확보를 달성하였으며, 더 적은 자유도로도 기존의(piecewise constant, 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장 계산 문제를 해결한다.
• 특히 굵은 메esh 해상도에서 정확도와 수렴 안정성 측면에서 저차수(piecewise constant) 기저 함수의 한계를 극복한다.
• 메쉬 세분화가 종종 불가능한 초고장력 MRI에서 환자 맞춤형 SAR 계산 및 RF 코ils 모델링에 적합한 안정적이고 빠른 VIE 해법을 개발한다.
• 계산 비용이나 자유도 증가 없이도 정확도를 확보하기 위해 h-세분화 대신 p-세분화(고차수 기저 함수)를 적용할 수 있도록 한다.
• 기존의 SIE 기반 프레임워크에 통합하여 실제 인간 신체 모델(RHBMs)에서의 EM 상호작용을 하이브리드 모델링할 수 있도록 한다.

제안 방법

• 안정적인 체적 적분방정식(VIE) 형식화를 위해 등가 극화 및 자화 전류를 사용하여 EM 산란 문제를 수립한다.
• 균일한 볼록체(voxelized) 격자에 정의된 불연속(piecewise) 선형(PWL) 기저 및 시험 함수를 사용하여 VIE를 이산화함으로써,(piecewise constant) 함수보다 더 나은 장 표현을 확보한다.
• 갈레르킨 방법의 모멘트를 적용하여 갈레르킨 수직성과 수치적 안정성 및 정확도를 향상시킨다.
• 반복 해법(GMRES 등)에서 행렬-벡터 곱을 빠른 푸리에 변환(FFT)을 활용해 가속화함으로써 대규모 문제에 대한 효율적 계산을 가능하게 한다.
• 특수 설계된 조절자를 구현하여 반복 해법의 안정성을 제고하고, 특히 해상도가 높은 경우에 발산하는 문제를 해결한다.
• 균일한 볼록체 격자를 활용하여 FFT 기반의 효율적 행렬-벡터 곱산술을 실현하고, 기존의 MRI 모델링 파이프라인과의 호환성을 유지한다.

실험 결과

연구 질문

• RQ1 표준(piecewise constant, PWC) 함수에 비해 불연속(piecewise) 선형(PWL) 기저 함수가 고대비, 비균질 MRI 문제에서 VIE 해법의 정확도와 수렴 안정성에 뚜렷한 향상을 이끌 수 있는가?
• RQ2 PWL 기저 함수를 사용한 p-세분화가 PWC 기저 함수를 사용한 h-세분화와 비교해 자유도를 줄이면서도 동일한 정확도를 달성할 수 있는가?
• RQ3 PWL 기저 함수를 사용한 현재 기반 VIE 형식이 매우 높은 해상도(예: 1 mm)에서도 발산 없이 수치적 안정성과 수렴성을 유지하는가?
• RQ4 고유성율이 높은 셰이밍 패드와 같은 도전적인 상황에서, 제안된 조절자가 반복 해법의 수렴 속도를 얼마나 효과적으로 가속화하는가?
• 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장 계산을 가능하게 하여 고대비 상황에서의 강건성(robustness)을 입증한다.