An Improved Cylindrical FDTD Algorithm and Its Application to Field-Tissue Interaction Study in MRI

This paper presents a three-dimensional finite-difference time-domain (FDTD) scheme in cylindrical coordinates with an improved algorithm for accommodating the numerical singularity associated with the polar axis. The regularization of this singularity problem is realized based on Faraday's law (for updating the Hr-component) and Ampere's law (for updating the Eφ -, Ez -components). The proposed algorithm has been detailed and verified against a problem with a known solution obtained from a commercial electromagnetic simulation package. The numerical scheme is also illustrated by modeling high-frequency RF field-human body interactions in magnetic resonance imaging (MRI). The results demonstrate the accuracy and capability of the proposed algorithm.

[1]  Jeff H Duyn,et al.  Time-domain finite-difference/finite-element hybrid simulations of radio frequency coils in magnetic resonance imaging. , 2008, Physics in medicine and biology.

[2]  Z. Kancleris.  Handling of Singularity in Finite-Difference Time-Domain Procedure for Solving Maxwell's Equations in Cylindrical Coordinate System , 2008, IEEE Transactions on Antennas and Propagation.

[3]  Jieru Chi,et al.  GPU-Accelerated FDTD Modeling of Radio-Frequency Field–Tissue Interactions in High-Field MRI , 2011, IEEE Transactions on Biomedical Engineering.

[4]  Raj Mittra,et al.  Finite-difference time-domain algorithm for solving Maxwell's equations in rotationally symmetric geometries , 1996 .

[5]  R. Janaswamy,et al.  An unstaggered colocated finite-difference scheme for solving time-domain Maxwell's equations in curvilinear coordinates , 1997 .

[6]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

[7]  Stuart Crozier,et al.  Hybrid numerical techniques for the modelling of radiofrequency coils in MRI , 2009, NMR in biomedicine.

[8]  Feng Liu,et al.  Analysis of Transient Eddy Currents in MRI Using a Cylindrical FDTD Method , 2006, IEEE Transactions on Applied Superconductivity.

[9]  A. Cangellaris,et al.  Analysis of the numerical error caused by the stair-stepped approximation of a conducting boundary in FDTD simulations of electromagnetic phenomena , 1991 .

[10]  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.

[11]  Serkan Aksoy,et al.  An FDTD Model for Calculation of Gradient-Induced Eddy Currents in MRI System , 2011 .

[12]  Dennis W. Prather,et al.  FORMULATION AND APPLICATION OF THE FINITE-DIFFERENCE TIME-DOMAIN METHOD FOR THE ANALYSIS OF AXIALLY SYMMETRIC DIFFRACTIVE OPTICAL ELEMENTS , 1999 .

[13]  S Crozier,et al.  An object-oriented designed finite-difference time-domain simulator for electromagnetic analysis and design in MRI--applications to high field analyses. , 2005, Journal of magnetic resonance.

[14]  Stuart Crozier,et al.  A theoretical comparison of two optimization methods for radiofrequency drive schemes in high frequency MRI resonators , 2005, Physics in medicine and biology.