1. Introduction
Magnetic resonance imaging (MRI) is a non-ionizing and non-invasive diagnostic technique based on the nuclear magnetic resonance (NMR) phenomenon.
In MRI systems, the radiofrequency (RF) field is generated and picked up by using coils, which should be able to support the wide field of view (FOV) with high magnetic field homogeneity in transmission and to achieve a high signal-to-noise ratio (SNR) in reception [
1]. The SNR, which depends on the hardware, on the acquisition of sequence parameters and on the tissue relaxation properties, is an accepted standard for quality measurements in MR [
2]. In order to design RF coils with optimized performance for a given application, an accurate simulation process is a necessary constraint. Such simulation has to permit the selection of coil parameters (sizes and conductor geometry) for ensuring the optimal SNR, which depends on the total losses (coil resistance and sample-induced resistance) [
3]. The sample-induced losses are caused by RF currents in the sample, induced by the fluctuating magnetic field, and by electric fields in the sample, mainly generated by the coil capacitors, while the coil losses comprise conductor resistance and radiation resistance, which take into account cross-sectional conductor geometry and the so-called “antenna effect”, respectively. Losses in the conductor can become the dominant losses contribution, especially for low-field MR systems, while radiative losses increase with the radiation frequency and coil size [
4]. As described in the literature [
5], sample-induced resistance can be estimated with numerical tools employing the finite-difference time-domain (FDTD) algorithm, while for coil losses, estimation, analytical formulations [
6] and finite element method (FEM) [
7] were employed. In particular, the FEM provides the great advantage of appropriately modeling structures with a small radius of curvature [
8].
FDTD modeling has been used in many designs. Plasmonic non-uniform nano-gratings for surface plasmon resonance (SPR) sensing and imaging were studied with FDTD to determine their optical characteristics [
9], as well as the optical response of a mechanically tunable ultra-narrowband optical filter based on a one-dimensional quasi-periodic photonic crystal (QPPC), and an on-chip integrated MEMS actuator was estimated with the same method [
10].
FDTD was even employed for modeling the optical properties of metal nanoparticles (MNPs) and recently a 3D-optimized FDTD (OFDTD) version was proposed for reducing computational requirements and simulation time by introducing new FDTD approximation terms based on the physical events occurring during the plasmonic oscillations in MNP [
11].
Even more recently [
12], a least-squares finite-difference time-domain method (LS-FDTD) was formulated in order to attenuate the high-frequency non-physical modes, superimposed to the physical solution, produced by Yee’s space discretization every time that time step is larger than the Courant–Friedrichs–Lewy (CFL) limit. Such a method, which provides simple central approximations for spatial derivatives in Maxwell’s equations, permitted us to obtain computer implementation simplicity and considerable processing time gains.
Further developments and applications in electromagnetic simulations will be possible thanks to recent upgrades, such as the leapfrog scheme for the unconditionally stable complying divergence implicit (CDI) finite-difference time-domain (FDTD) method, characterized by unconditional stability, complying divergence and an efficient leapfrog scheme with reduced right-hand side (RHS) flops [
13].
In this paper, we propose the application of an FDTD algorithm for separately estimating conductor and radiative losses in a circular loop constituted by a wire (cylindrical rod shape) conductor from 21 to 128 MHz, corresponding to a static field from 0.5 to 3.0 T. To our knowledge, this simulation approach has not been fully theoretically and experimentally validated yet for MR coil loss contribution estimation, since copper parts are usually modeled as perfect electric conductors in FDTD simulations. The final objective of this work was to demonstrate that an FDTD-based coil SNR model is able to supply all loss contributions (coil- and sample-induced resistance values) and magnetic field pattern calculation, without approximations in sample and coil geometries.
Results provided by such FDTD simulations for conductor and radiation resistance estimations were compared with analytical, FEM and workbench results obtained with a home-built coil prototype.
3. Results
Table 1 and
Table 2 show, respectively, conductor resistance and radiation resistance results for the non-segmented loop (one-feed coil) obtained by FDTD simulations at different frequencies compared with the FEM (CST—computer simulation technology, AG, Darmstadt, Germany) values [
7].
By naming
Rtot the sum of conductor resistance
Rcond and radiation resistance
Rrad, the experimental measurements, published previously, provided values of
Rtot equal to 260 mΩ at 63.9 MHz and 950 mΩ at 127.8 MHz [
18]. These experimental values included further resistive losses attributable to the solder joints between the coil and the cable for the connection with the analyzer (
Rsol = 25 mΩ and 60 mΩ estimated at 64 MHz and 128 MHz [
16], respectively). Total resistances
Rtot estimated by FDTD and FEM along with experimental values reduced by soldering losses are listed in
Table 3.
Table 4 and
Table 5 summarize, respectively, conductor resistance and radiation resistance results for the segmented loop obtained by FDTD simulations at different frequencies and compared with the analytical calculated values.
Table 6 and
Table 7 show, respectively, conductor resistance and radiation resistance results for the
n = 2-segment loop obtained by FDTD simulations at different frequencies and compared with the analytical calculated values.
Table 8 and
Table 9 show, respectively, conductor resistance and radiation resistance results for the
n = 4 segmented loop obtained by FDTD simulations at different frequencies and compared with the analytical calculated values.
Figure 2 depicts the plots of electric (E) and magnetic field (H) variation across the coil for the
n = 2, 4 and 8 breaks at 127.8 MHz.
4. Discussion
In general, FDTD-based simulation tools permit the inclusion of the computational space complex structures, such as a part of the human body, to simulate various geometry systems, without approximations in sample and coil geometries [
21]. The literature described coil simulations where an SNR coil model was developed by using FDTD for sample-induced resistance and magnetic-field pattern estimation, but the coil conductor resistance was calculated by Ohm’s law, and radiation resistance was neglected [
22]. However, Ohm’s law can be easily applied for conductor resistance calculation only for simple coil geometries, while in some cases the radiation losses cannot be neglected, and they can be significantly different for the loaded coils and when it is placed in the scanner bore.
Although FDTD was not considered good at estimating coil losses and much inferior to FEM, we demonstrated that with an optimal simulation setup, it can provide results with a good agreement with FEM, analytical and measurement results.
In particular, for a non-segmented loop, which essentially acts as a folded dipole [
23], when the loop circumference of the loop approaches a significant fraction of the wavelength, a non-uniform current flows on it. In these conditions, FDTD was able to provide results similar to the ones obtained with FEM (relative differences of <4.35% and <6.62% for, respectively, conductor and radiation resistances). Moreover, the sum of conductor resistance and radiation resistance
Rtot calculated with FDTD provides values closer to experimental results with respect to FEM analysis, as listed in
Table 3. Conversely, by segmenting this loop with
n = 8 ports and feeding them concurrently, both conductor and radiation resistances are extremely close to the analytical method results, which contemplates a uniform current in the coil path (relative differences of <1.59% and <0.50% for, respectively, conductor and radiation resistances).
Simulation results performed with
n = 2 and
n = 4 breaks demonstrated that the relative difference of the FDTD-estimated coil conductor and radiation resistances with respect to analytical calculation increase when the break number diminishes, because the current distribution in the loop becomes less uniform, as can be guessed from the E field and H field patterns shown in
Figure 2. In particular, for
n = 2 the relative differences were 25.21% and 10.13% for, respectively, conductor and radiation resistances at 127.8 MHz, while the differences were reduced to 16.31% (conductor resistance) and 2.11% (radiation resistance) for
n = 4 at the same frequency.
As predicted by theory, the radiation losses were higher in the non-segmented loop case due to the presence of larger field heterogeneities induced by the RF current phase shift on longer conductors with respect to the segmented loop. The segmentation also has the advantage of distributing the electric field all around the coil and reducing the coil antenna mode contribution [
4]. We believe that the optimal number of breaks can be chosen by considering that each coil segment has to be a very small fraction (<1/20) of the wavelength associated with the highest frequency (about
λ/40 in our
n = 8 case).