Numerical modelling of thermal effects in rats due to high-field magnetic resonance imaging (0.5-1 GHz)

A finite-difference time-domain (FDTD) thermal model has been developed to compute the temperature elevation in the Sprague Dawley rat due to electromagnetic energy deposition in high-field magnetic resonance imaging (MRI). The field strengths examined ranged from 11.75-23.5 T (corresponding to H-1 resonances of 0.5-1 GHz) and an N-stub birdcage resonator was used to both transmit radio-frequency energy and receive the MRI signals. With an in-plane resolution of 1.95 mm, the inhomogeneous rat phantom forms a segmented model of 12 different tissue types, each having its electrical and thermal parameters assigned. The steady-state temperature distribution was calculated using a Pennes 'bioheat' approach. The numerical algorithm used to calculate the induced temperature distribution has been successfully validated against analytical solutions in the form of simplified spherical models with electrical and thermal properties of rat muscle. As well as assisting with the design of MRI experiments and apparatus, the numerical procedures developed in this study could help in future research and design of tumour-treating hyperthermia applicators to be used on rats in vivo.

A distributed equivalent magnetic current based FDTD method for the calculation of E-fields induced by gradient coils

This paper evaluates a new, low-frequency finite-difference time-domain method applied to the problem of induced E-fields/eddy currents in the human body resulting from the pulsed magnetic field gradients in MRI. In this algorithm, a distributed equivalent magnetic current is proposed as the electromagnetic source and is obtained by quasistatic calculation of the empty coil's vector potential or measurements therein. This technique circumvents the discretization of complicated gradient coil geometries into a mesh of Yee cells, and thereby enables any type of gradient coil modelling or other complex low frequency sources. The proposed method has been verified against an example with an analytical solution. Results are presented showing the spatial distribution of gradient-induced electric fields in a multi-layered spherical phantom model and a complete body model. (C) 2004 Elsevier Inc. All rights reserved.

Electromagnetic fields inside a lossy, multilayered spherical head phantom excited by MRI coils: models and methods

The precise evaluation of electromagnetic field (EMF) distributions inside biological samples is becoming an increasingly important design requirement for high field MRI systems. In evaluating the induced fields caused by magnetic field gradients and RF transmitter coils, a multilayered dielectric spherical head model is proposed to provide a better understanding of electromagnetic interactions when compared to a traditional homogeneous head phantom. This paper presents Debye potential (DP) and Dyadic Green's function (DGF)-based solutions of the EMFs inside a head-sized, stratified sphere with similar radial conductivity and permittivity profiles as a human head. The DP approach is formulated for the symmetric case in which the source is a circular loop carrying a harmonic-formed current over a wide frequency range. The DGF method is developed for generic cases in which the source may be any kind of RF coil whose current distribution can be evaluated using the method of moments. The calculated EMFs can then be used to deduce MRI imaging parameters. The proposed methods, while not representing the full complexity of a head model, offer advantages in rapid prototyping as the computation times are much lower than a full finite difference time domain calculation using a complex head model. Test examples demonstrate the capability of the proposed models/methods. It is anticipated that this model will be of particular value for high field MRI applications, especially the rapid evaluation of RF resonator (surface and volume coils) and high performance gradient set designs.

Winding patterns for biplanar MRI shim coils with rectangular and circular target-field regions

A method is presented for calculating the winding patterns required to design independent zonal and tesseral biplanar shim coils for magnetic resonance imaging. Streamline, target-field, Fourier integral and Fourier series methods are utilized. For both Fourier-based methods, the desired target field is specified on the surface of the conducting plates. For the Fourier series method it is possible to specify the target field at additional depths interior to the two conducting plates. The conducting plates are confined symmetrically in the xy plane with dimensions 2a x 2b, and are separated by 2d in the z direction. The specification of the target field is symmetric for the Fourier integral method, but can be over some asymmetric portion pa < x < qa and sb < y < tb of the coil dimensions (-1 < p < q < 1 and -1 < s < t < 1) for the Fourier series method. Arbitrary functions are used in the outer sections to ensure continuity of the magnetic field across the entire coil face. For the Fourier series case, the entire field is periodically extended as double half-range sine or cosine series. The resultant Fourier coefficients are substituted into the Fourier series and integral expressions for the internal and external magnetic fields, and stream functions on both the conducting surfaces. A contour plot of the stream function directly gives the required coil winding patterns. Spherical harmonic analysis of field calculations from a ZX shim coil indicates that example designs and theory are well matched.

Novel target-field method for designing shielded biplanar shim and gradient coils

A new method is presented here for the systematic design of biplanar shielded shim and gradient coils, for use in magnetic resonance imaging (MRI) and other applications. The desired target field interior to the coil is specified in advance, and a winding pattern is then designed to produce a field that matches the target as closely as possible. Both gradient and shim coils can be designed by this approach, and the target region can be located asymmetrically within the coil. The interior target field may be matched at two or more interior locations, to improve accuracy. When shields are present, the winding patterns are designed so that the fields exterior to the biplanar coil are made as small as possible. The method is illustrated here by the design of some transverse gradient and shim coils.

The design of planar gradient coils. Part ll: A weighted superposition method

In this work a superposition technique for designing gradient coils for the purpose of magnetic resonance imaging is outlined, which uses an optimized weight function superimposed upon an initial winding similar to that obtained from the target field method to generate the final wire winding. This work builds on the preliminary work performed in Part I on designing planar insertable gradient coils for high resolution imaging. The proposed superposition method for designing gradient coils results in coil patterns with relatively low inductances and the gradient coils can be used as inserts into existing magnetic resonance imaging hardware. The new scheme has the capacity to obtain images faster with more detail due to the deliver of greater magnetic held gradients. The proposed method for designing gradient coils is compared with a variant of the state-of-the-art target field method for planar gradient coils designs, and it is shown that the weighted superposition approach outperforms the well-known the classical method.

An inverse methodology for high-frequency RF coil design for MRI with de-emphasized B-1 fields

An inverse methodology for the design of biologically loaded radio-frequency (RF) coils for magnetic resonance imaging applications is described. Free space time-harmonic electromagnetic Green's functions and de-emphasized B-1 target fields are used to calculate the current density on the coil cylinder. In theory, with the B-1 field de-emphasized in the middle of the RF transverse plane, the calculated current distribution can generate an internal magnetic field that can reduce the central overemphasis effect caused by field/tissue interactions at high frequencies. The current distribution of a head coil operating at 4 T (170 MHz) is calculated using an inverse methodology with de-emphasized B-1. target fields. An in-house finite-difference time-domain routine is employed to evaluate B-1 field and signal intensity inside a homogenous cylindrical phantom and then a complete human head model. A comparison with a conventional RF birdcage coil is carried out and demonstrates that this method can help in decreasing the normal bright region caused by field/tissue interactions in head images at 170 MHz and higher field strengths.

An object-oriented designed finite-difference time-domain simulator for electromagnetic analysis and design in MRI - application to high field analyses

This paper presents a finite-difference time-domain (FDTD) simulator for electromagnetic analysis and design applications in MRI. It is intended to be a complete FDTD model of an MRI system including all RF and low-frequency field generating units and electrical models of the patient. The pro-ram has been constructed in an object-oriented framework. The design procedure is detailed and the numerical solver has been verified against analytical solutions for simple cases and also applied to various field calculation problems. In particular, the simulator is demonstrated for inverse RF coil design, optimized source profile generation, and parallel imaging in high-frequency situations. The examples show new developments enabled by the simulator and demonstrate that the proposed FDTD framework can be used to analyze large-scale computational electromagnetic problems in modern MRI engineering. (C) 2004 Elsevier Inc. All rights reserved.

A target-field method to design circular biplanar coils for asymmetric shim and gradient fields

The paper presents a method for designing circular, shielded biplanar coils that can generate any desired field. A particular feature of these coils is that the target field may be located asymmetrically within the coil. A transverse component of the magnetic field produced by the coil is made to match a prescribed target field over the surfaces of two concentric spheres (the diameter of spherical volume) that define the target field location. The paper shows winding patterns and fields for several gradient and shim coils. It examines the effect that the finite coil size has on the winding patterns, using a Fourier-transform calculation for comparison.

A time-harmonic target-field method for designing unshielded RF coils in MRI

Time-harmonic methods are required in the accurate design of RF coils as operating frequency increases. This paper presents such a method to find a current density solution on the coil that will induce some desired magnetic field upon an asymmetrically located target region within. This inverse method appropriately considers the geometry of the coil via a Fourier series expansion, and incorporates some new regularization penalty functions in the solution process. A new technique is introduced by which the complex, time-dependent current density solution is approximated by a static coil winding pattern. Several winding pattern solutions are given, with more complex winding patterns corresponding to more desirable induced magnetic fields.

A wave equation technique for designing compact gradient coils

A primary purpose of this research is to design a gradient coil that is planar in construction and can be inserted within existing infrastructure. The proposed wave equation method for the design of gradient coils is novel within the field. it is comprehensively shown how this method can be used to design the planar x-, y-, and z-gradient wire windings to produce the required magnetic fields within a certain domain. The solution for the cylindrical gradient coil set is also elucidated. The wave equation technique is compared with the well-known target held method to gauge the quality of resultant design. In the case of the planar gradient coil design, it is shown that using the new method, a set of compact gradient coils with large field of view can be produced. The final design is considerably smaller in dimension when compared with the design obtained using the target field method, and therefore the manufacturing costs and materials required are somewhat reduced.

A new decoupling method for quadrature coils in magnetic resonance imaging

A powerful decoupling method is introduced to obtain decoupled signal voltages from quadrature coils in magnetic resonance imaging (MRI). The new method uses the knowledge of the position of the signal source in MRI, the active slice, to define a new mutual impedance which accurately quantifies the coupling voltages and enables them to be removed almost completely. Results show that by using the new decoupling method, the percentage errors in the decoupled voltages are of the order of 10(-7)% and isolations between two coils are more than 170 dB.

An inverse method for designing loaded RF coils in MRI

Radio-frequency ( RF) coils are designed such that they induce homogeneous magnetic fields within some region of interest within a magnetic resonance imaging ( MRI) scanner. Loading the scanner with a patient disrupts the homogeneity of these fields and can lead to a considerable degradation of the quality of the acquired image. In this paper, an inverse method is presented for designing RF coils, in which the presence of a load ( patient) within the MRI scanner is accounted for in the model. To approximate the finite length of the coil, a Fourier series expansion is considered for the coil current density and for the induced fields. Regularization is used to solve this ill-conditioned inverse problem for the unknown Fourier coefficients. That is, the error between the induced and homogeneous target fields is minimized along with an additional constraint, chosen in this paper to represent the curvature of the coil windings. Smooth winding patterns are obtained for both unloaded and loaded coils. RF fields with a high level of homogeneity are obtained in the unloaded case and a limit to the level of homogeneity attainable is observed in the loaded case.