An Improved Quasi-static Finite-difference Scheme for Induced Field Evaluation in MRI based on the Biconjugate Gradient Method

In modern magnetic resonance imaging (MRI), there are concerns for the health and safety of patients and workers repeatedly exposed to magnetic fields, and therefore accurate and efficient evaluation of in situ electromagnetic field (EMF) distributions has gained a lot of significance. This paper presents a Biconjugate Gradient Method (BiCG) to efficiently implement the quasi-static finite-difference scheme (QSFD), which has been widely utilized to model and analyze magnetically induced electric fields and currents within the human body during the operation of the MRI systems and in other settings. The proposed BiCG method shows computational advantages over the iterative, successive over-relaxation (SOR) algorithm. The scheme has been validated against other known solutions on a lossy, multilayered ellipsoid phantom excited by an ideal loop coil. Numerical results on a 3D human body model demonstrate that the convergence time and memory consumption is significantly reduced using the BiCG method.

[1]  Adnan Trakic,et al.  Numerical study of currents in workers induced by body‐motion around high‐ultrahigh field MRI magnets , 2007, Journal of magnetic resonance imaging : JMRI.

[2]  Stuart Crozier,et al.  Numerical evaluation of the fields induced by body motion in or near high-field MRI scanners. , 2005, Progress in biophysics and molecular biology.

[3]  S. Crozier,et al.  Numerical field evaluation of healthcare workers when bending towards high‐field MRI magnets , 2008, Magnetic resonance in medicine.

[4]  Stuart Crozier,et al.  Transient temperature rise in a mouse due to low-frequency regional hyperthermia , 2006, Physics in medicine and biology.

[5]  P. Dimbylow.  Development of the female voxel phantom, NAOMI, and its application to calculations of induced current densities and electric fields from applied low frequency magnetic and electric fields , 2005, Physics in medicine and biology.

[6]  Hiroshi Kanai,et al.  Impedence Method for Calculation of Power Deposition Patterns in Magnetically Induced Hyperthermia , 1984, IEEE Transactions on Biomedical Engineering.

[7]  F. X. Canning.  Preconditioned iterative solution of IML/moment method problems , 1993 .

[8]  C Gabriel,et al.  The dielectric properties of biological tissues: I. Literature survey. , 1996, Physics in medicine and biology.

[9]  Stuart Crozier,et al.  Calculation of electric fields induced by body and head motion in high-field MRI. , 2003, Journal of magnetic resonance.

[10]  H. V. D. Vorst,et al.  An overview of approaches for the stable computation of hybrid BiCG methods , 1995 .

[11]  S Crozier,et al.  Longitudinal gradient coil optimization in the presence of transient eddy currents , 2007, Magnetic resonance in medicine.

[12]  F Schmitt,et al.  Physiological effects of fast oscillating magnetic field gradients. , 1991, Journal of computer assisted tomography.

[13]  Feng Liu,et al.  On the induced electric field gradients in the human body for magnetic stimulation by gradient coils in MRI , 2003, IEEE Transactions on Biomedical Engineering.

[14]  D. Sorensen.  Numerical methods for large eigenvalue problems , 2002, Acta Numerica.

[15]  Stuart Crozier,et al.  Exposure of workers to pulsed gradients in MRI , 2007, Journal of magnetic resonance imaging : JMRI.

[16]  Seyedeh Mahboubeh Molavi-Arabshahi,et al.  Preconditioned techniques for solving large sparse linear systems arising from the discretization of the elliptic partial differential equations , 2007, Appl. Math. Comput..

[17]  Stuart Crozier,et al.  Finite-difference time-domain-based studies of MRI pulsed field gradient-induced eddy currents inside the human body , 2002 .

[18]  F. Liu,et al.  Electromagnetic fields inside a lossy, multilayered spherical head phantom excited by MRI coils: models and methods. , 2004, Physics in medicine and biology.

[19]  J. Shewchuk.  An Introduction to the Conjugate Gradient Method Without the Agonizing Pain , 1994 .

[20]  I. Duff.  A survey of sparse matrix research , 1977, Proceedings of the IEEE.

[21]  Stuart Crozier,et al.  The Design of Transverse Gradient Coils of Restricted Length by Simulated Annealing , 1994 .

[22]  P. Dimbylow.  Current densities in a 2 mm resolution anatomically realistic model of the body induced by low frequency electric fields. , 2000, Physics in medicine and biology.

[23]  Judith A. Vogel,et al.  Flexible BiCG and flexible Bi-CGSTAB for nonsymmetric linear systems , 2007, Appl. Math. Comput..

[24]  H. Ling,et al.  An efficient wavelet preconditioner for iterative solution of three-dimensional electromagnetic integral equations , 2003 .