Simulation of RF Electromagnetic Fields for MRI: Methods and Applications
Christopher Collins1
1New York University School of Medicine, United States

Synopsis

Keywords: Physics & Engineering: RF Safety, Physics & Engineering: Hardware, Physics & Engineering: Physics

Simulations of electromagnetic fields have many valuable uses for design and evaluation of RF coils with consideration of both performance and safety. What methods are best depend on the intended application, available resources, and preferences of the researcher. Successful use of RF simulations for MRI requires good understanding of the simulation process, the expected behavior of RF coils and the fields they produce, and how these fields relate to signal and/or heating during MRI. In this presentation, we will cover important concepts and methods for simulating RF coils, fields, and their effects in MRI.

Introduction:

RF fields in MRI are applied to manipulate nuclear magnetic moments, but in doing so cause heating of tissues. RF coils are used to detect signal from nuclear magnetic moments, but are also sensitive to noise induced by random motion of ions and molecules throughout the body. While we typically desire a homogeneous field for nuclear excitation, this is increasingly difficult to achieve as we move to higher frequencies in search of better SNR. And while we desire a configuration of receive coils with high sensitivity to signal and field distributions optimal for SNR and accelerated imaging, design of such arrays can be confounded by the variation of field distributions with subject geometry.
Simulation of RF fields in MRI is increasingly used to understand observed effects, ensure patient safety (including in the presence of active and passive metallic implants), design transmit and receive coils and arrays, and design pulses and sequences that maximize desired image characteristics, minimize undesired effects and artifacts. There are a wide variety of approaches to simulation and interpretation of results, and the best methods are often dependent on the particular application.
In this presentation we will review a wide variety of methods for field calculation, representation of coils and the human body, interpretation of results, and compromises between these approaches for a wide array of engineering problems in MRI.

Consideration of Commonly-used Simulation Methods and Settings for MRI:

The following figures summarize some considerations in selecting methods and settings for RF field simulations. Simulation methods included here are Finite Difference Time Domain (FDTD), Finite Integration Technique (FIT), Finite Element Method (FEM), Boundary Element Method (BEM), and Method of Moments. When applicable, commercial tools in each category are listed in alphabetical order by manufacturer name.
Regarding commercially available simulation methods, FDTD and FIT methods (Microwave Studio, Sim4Life, xFDTD) most often use a rectilinear grid and can simulate high-resolution multi-tissue anatomical models with relatively little computational resources (memory and time), whereas FEM (HFSS), BEM (COMSOL), and MoM (Feko), not being constrained to rectlinear representations, can simulate conductive surfaces in coils, shields, and implants with arbitrary orientation, but generally have limitations with the anatomical detail that can be simulated with reasonable computational resources. Between FEM, BEM, and MoM, generally FEM has higher memory requirements for a given problem, but also better ability to adapt to non-sparse problem geometries, including human body models with a surface-based format.
Regarding coil driving methods, if coil current and/or voltage distributions are well understood and the goal of the simulations is to understand field distributions in the sample or patient, directly implementing those current or voltage distributions (especially with sources across gaps where capacitors would be placed in the physical coil) can save significant time. On the other hand, if the purpose of the simulation is more about designing the coil or array to perform in a desired manner or to relate the forward power from the amplifier to the power absorbed in the subject, modeling the coil circuitry very exactly becomes important. Depending on the number of electrical components (especially capacitor values) to be optimized, simulating the coil with multiple ports in the place of all capacitors followed by use of a circuit co-simulator can be very valuable.
Often a major consideration regarding time and memory requirements for a simulation, as well as the storage space for results, is the spatial resolution of the grid or mesh used. It often takes a significant amount of experience with a software package to design efficient grids or meshes. When it comes to evaluating the results of a simulation, while the various available software packages include a differing degree of tools for analysis relevant to MRI (sometimes including temperature simulation tools or even a Bloch-based MRI simulator), often advanced analysis is better accomplished in other environments, such as with Matlab. In these cases, the ability to easily export data for use in another environment is an important feature.
Besides the commercial packages listed above, a number of freely-available options exist (with the understanding that ease of use – especially regarding GUIs for setting up simulations) exist. Limited versions of some of the commercial packages are also available for use at academic institutions free of charge.

Acknowledgements

No acknowledgement found.

References

Recommended reading:

Collins CM and Wang Z. Calculation of radiofrequency electromagnetic fields and their effects in MRI of human subjects. Magn Reson Med 2011;65:1470–1482

Hand JW. Modelling the interaction of electromagnetic fields (10 MHz–10 GHz) with the human body: methods and applications. Phys. Med. Biol. 2008;53:R243–R286.

Hoult DI and Lauterbur PC. The sensitivity of the zeugmatographic experiment involving human samples. J Magn Reson 1979; 34: 425-433

Edelstein WA, Glover GH, Hardy CJ, Redington RW. The intrinsic signal-to-noise ratio in NMR imaging. Magn. Reson. Med. 1986;3:604-618.

Roemer PB, Edelstein WA, Hayes CE, Souza SP, Mueller OM. The NMR phased array. Magn Reson Med 1990; 16: 192-225. Kozlov M, Turner R. Fast MRI coil analysis based on 3-D electromagnetic and RF circuit co-simulation. Journal of Magnetic Resonance 2009;200:147–152.

Kozlov M, Turner R. Fast MRI coil analysis based on 3-D electromagnetic and RF circuit co-simulation. Journal of Magnetic Resonance 2009;200:147–152.

Eichfelder G, Gebhardt M. Local specific absorption rate control for parallel transmission by virtual observation points. Magnetic resonance in medicine. 2011 Nov;66(5):1468-76.

Collins CM, Liu W, Wang JH, Gruetter R, Vaughan JT, Ugurbil K, and Smith MB, “Temperature and SAR calculations for a human head within volume and surface coils at 64 and 300 MHz.” J. Magn. Reson. Imag., 19:650-656, 2004.

Wiesinger F, Boesiger P, Pruessmann KP. Electrodynamics and Ultimate SNR in Parallel MR Imaging. Magnetic Resonance in Medicine 2004;52:376–390.

Wiesinger F, Van de Moortele P-F, Adriany G, De Zanche N, Ugurbil K, Pruessmann KP. Parallel imaging performance as a function of field strength—an experimental investigation using electrodynamic scaling. Magn Reson Med 2004; 52: 953–964.

Christ A, Kainz W, Hahn EG, Honegger K, Zefferer M, Neufeld E, Rascher W, Janka R, Bautz W, Chen J, Kiefer B, Schmitt P, Hollenbach H-P, Shen J, Oberle M, Szczerba D, Kam A, Guag JW, Kuster N. The virtual family – development of surface-based anatomical models of two adults and two children for dosimetric simulations. Phys Med Biol 2010;55:23-38.

Metzger GJ, Snyder C, Akgun C, Vaughan JT, Ugurbil K, Van de Moortele P-F. Local B1+ shimming for prostate imaging with transceiver arrays at7T based on subject-dependent transmit phase measurements. Magnetic Resonance in Medicine 2008;59:396–409.

Cloos MA, Luong M, Ferrand G, Amadon A, LeBihan D, Boulant N. Local SAR reduction in parallel excitation based on channel-dependent tikhonov parameters. J Magn Reson Imaging 2010;32:1209-1216.

Li BK, Liu F, Crozier S. Focused, eight-element transceive phased array coil for parallel magnetic resonance imaging of the chest - theoretical considerations. Magn Reson Med 2005; 53: 1251–1257.

Bassen H, Kainz W, Mendoza G, Kellom T. MRI-induced heating of thin wire metallic implants – laboratory and computational studies – findings and new questions raised. Minimally Invasive Therapy and Allied Technologies, 2006;15:76 – 84.

Bankson JA, Wright SM. Simulation-Based Investigation of Partially Parallel Imaging With a Linear Array at High Accelerations. Magnetic Resonance in Medicine 2002;47:777–786.

Cao Z, Oh S, Sica CT, McGarrity JM, Horan T, Luo W, Collins CM. Bloch-based MRI System Simulator Considering Realistic Electromagnetic Fields for Calculation of Signal, Noise, and SAR. Magn Reson Med 2014;72:237-247.

Proc. Intl. Soc. Mag. Reson. Med. 31 (2023)