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.

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.

Synthesis of gadolinium-labeled shell-crosslinked nanoparticles for magnetic resonance imaging applications

Shell-crosslinked knedel-like nanoparticles (SCKs; knedel is a Polish term for dumplings) were derivatized with gadolinium Shell chelates and studied as robust magnetic-resonance-imaging-active structures with hydrodynamic diameters of 40 +/- 3 nm. SCKs possessing an amphiphilic core-shell morphology were produced from the aqueous assembly of diblock copolymers of poly(acrylic acid) (PAA) and poly(methyl acrylate) (PMA), PAA(52)-b-PMA(128), and subsequent covalent crosslinking by amidation upon reaction with 2,2'-(ethylenedioxy)bis(ethylamine) throughout the shell layer. The properties of these materials, including non-toxicity towards mammalian cells, non-immunogenicity within mice, and capability for polyvalent targeting, make them ideal candidates for utilization within biological systems. The synthesis of SCKs derivatized with Gd-III and designed for potential use as a unique nanometer-scale contrast agent for MRI applications is described herein. Utilization of an amino-functionalized diethylenetriaminepentaacetic acid-Gd analogue allowed for direct covalent conjugation throughout the hydrophilic shell layer of the SCKs and served to increase the rotational correlation lifetime of the Gd. In addition, the highly hydrated nature of the shell layer in which the Gd was located allowed for rapid water exchange; thus, the resulting material demonstrated large ionic relaxivities (39 s(-1) mM(-1)) in an applied magnetic field of 0.47 T at 40 degrees C and, as a result of the large loading capacity of the material, also demonstrated high molecular relaxivities (20 000 s(-1) mM(-1)).

Numerical evaluation of the fields induced by body motion in or near high-field MRI scanners

In modern magnetic resonance imaging, both patients and health care workers are exposed to strong. non-uniform static magnetic fields inside and outside of the scanner. In which body movement may be able to induce electric currents in tissues which could be potentially harmful. This paper presents theoretical investigations into the spatial distribution of induced E-fields in a tissue-equivalent human model when moving at various positions around the magnet. The numerical calculations are based on an efficient. quasi-static, finite-difference scheme. Three-dimensional field profiles from an actively shielded 4 T magnet system are used and the body model projected through the field profile with normalized velocity. The simulation shows that it is possible to induce E-fields/currents near the level of physiological significance under some circumstances and provides insight into the spatial characteristics of the induced fields. The methodology presented herein can be extrapolated to very high field strengths for the evaluation of the effects of motion at a variety of field strengths and velocities. (C) 2004 Elsevier Ltd. 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.

Magnetic susceptibility measurement applied to the minerals industry

The efficiency of numerous mineral processing operations can be determined by measuring the magnetic properties of the process streams. This measurement, if done at all, is currently performed by laboratory testing of spot samples. This is an inherently slow process and, where feed grades are variable, optimum performance is generally not achieved. This paper describes the rapid measurement of the magnetic properties of minerals. AC measurement techniques, including the analysis of the phase component of the magnetic vector; frequency dependent magnetic susceptibility and the effect of applied magnetic field strength will be discussed. Industrial applications in mineral sands, copper smelting, ferrosilicon testing and drill core scanning will be reported. (C) 2006 Elsevier Ltd. All rights reserved.

Modulation of channel activity and gadolinium block of MscL by static magnetic fields

The magnetic field of the Earth has for long been known to influence the behaviour and orientation of a variety of living organisms. Experimental studies of the magnetic sense have, however, been impaired by the lack of a plausible cellular and/or molecular mechanism providing meaningful explanation for detection of magnetic fields by these organisms. Recently, mechanosensitive (MS) ion channels have been implied to play a role in magnetoreception. In this study we have investigated the effect of static magnetic fields (SMFs) of moderate intensity on the activity and gadolinium block of MscL, the bacterial MS channel of large conductance, which has served as a model channel to study the basic physical principles of mechanosensory transduction in living cells. In addition to showing that direct application of the magnetic field decreased the activity of the MscL channel, our study demonstrates for the first time that SMFs can reverse the effect of gadolinium, a well-known blocker of MS channels. The results of our study are consistent with a notion that (1) the effects of SMFs on the MscL channels may result from changes in physical properties of the lipid bilayer due to diamagnetic anisotropy of phospholipid molecules and consequently (2) cooperative superdiamagnetism of phospholipid molecules under influence of SMFs could cause displacement of Gd3+ stop ions from the membrane bilayer and thus remove the MscL channel block.

An inverse methodology for high frequency RF head coil design with pre-emphasized B1 field in MRI

An inverse methodology to assist in the design of radio-frequency (RF) head coils for high field MRI application is described in this work. Free space time-harmonic electromagnetic Green's functions and preemphasized B1 field are used to calculate the current density on the coil cylinder. With B1 field preemphasized and lowered in the middle of the RF transverse plane, the calculated current distribution can generate an internal magnetic field that can reduce the EM field/tissue interactions at high frequencies. The current distribution of a head coil operating at 4 T is calculated using inverse methodology with preemphasized B1 fields. FDTD is employed to calculate B1 field and signal intensity inside a homogenous cylindrical phantom and human head. A comparison with conventional RF birdcage coil is reported here and demonstrated that inverse-method designed coil with preemphasized B1 field can help in decreasing the notorious bright region caused by EM field/tissue interactions in the human head images at 4 T.

An equivalent distributed magnetic current based FDTD method for the calculation of e-fields induced by gradient coils in MRI

This paper evaluates a low-frequency FDTD 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 (DEMC) 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 discretizing of complicated gradient coil geometries into a mesh of Yee cells, and thereby enables any type of gradient coil modeling 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 multilayered spherical phantom model and a complete body model.

A study of magnetic non-linearity and finite length effects in solid iron subjected to a travelling MMF wave

This thesis describes an experimental and analytic study of the effects of magnetic non-linearity and finite length on the loss and field distribution in solid iron due to a travelling mmf wave. In the first half of the thesis, a two-dimensional solution is developed which accounts for the effects of both magnetic non-linearity and eddy-current reaction; this solution is extended, in the second half, to a three-dimensional model. In the two-dimensional solution, new equations for loss and flux/pole are given; these equations contain the primary excitation, the machine parameters and factors describing the shape of the normal B-H curve. The solution applies to machines of any air-gap length. The conditions for maximum loss are defined, and generalised torque/frequency curves are obtained. A relationship between the peripheral component of magnetic field on the surface of the iron and the primary excitation is given. The effects of magnetic non-linearity and finite length are combined analytically by introducing an equivalent constant permeability into a linear three-dimensional analysis. The equivalent constant permeability is defined from the non-linear solution for the two-dimensional magnetic field at the axial centre of the machine to avoid iterative solutions. In the linear three-dimensional analysis, the primary excitation in the passive end-regions of the machine is set equal to zero and the secondary end faces are developed onto the air-gap surface. The analyses, and the assumptions on which they are based, were verified on an experimental machine which consists of a three-phase rotor and alternative solid iron stators, one with copper end rings, and one without copper end rings j the main dimensions of the two stators are identical. Measurements of torque, flux /pole, surface current density and radial power flow were obtained for both stators over a range of frequencies and excitations. Comparison of the measurements on the two stators enabled the individual effects of finite length and saturation to be identified, and the definition of constant equivalent permeability to be verified. The penetration of the peripheral flux into the stator with copper end rings was measured and compared with theoretical penetration curves. Agreement between measured and theoretical results was generally good.

The influence of age on the pattern and flash visual evoked magnetic response (VEMR)

The visual evoked magnetic response (VEMR) was measured over the occipital cortex to pattern and flash stimuli in 86 normal subjects aged 15-86 years. The latency of the major positive component (outgoing magnetic field) to the pattern reversal stimulus (P100M) increased with age, particularly after 55 years, while the amplitude of the P100M decreased more gradually over the lifespan. By contrast, the latency of the major positive component to the flash stimulus (P2M) increased more slowly with age after about 50 years, while its amplitude may have decreased in only a proportion of the elderly subjects. The changes in the P100M with age may reflect senile changes in the eye and optic nerve, e.g. senile miosis, degenerative changes in the retina or geniculostriate deficits. The P2M may be more susceptible to senile changes in the visual cortex. The data suggest that the contrast channels of visual information processing deteriorate more rapidly with age than the luminance channels.

Developments and extensions of Roth's method of double Fourier series with applications to the solution of electromagnetic field problems

The first part of the thesis compares Roth's method with other methods, in particular the method of separation of variables and the finite cosine transform method, for solving certain elliptic partial differential equations arising in practice. In particular we consider the solution of steady state problems associated with insulated conductors in rectangular slots. Roth's method has two main disadvantages namely the slow rate of con­vergence of the double Fourier series and the restrictive form of the allowable boundary conditions. A combined Roth-separation of variables method is derived to remove the restrictions on the form of the boundary conditions and various Chebyshev approximations are used to try to improve the rate of convergence of the series. All the techniques are then applied to the Neumann problem arising from balanced rectangular windings in a transformer window. Roth's method is then extended to deal with problems other than those resulting from static fields. First we consider a rectangular insulated conductor in a rectangular slot when the current is varying sinusoidally with time. An approximate method is also developed and compared with the exact method.The approximation is then used to consider the problem of an insulated conductor in a slot facing an air gap. We also consider the exact method applied to the determination of the eddy-current loss produced in an isolated rectangular conductor by a transverse magnetic field varying sinusoidally with time. The results obtained using Roth's method are critically compared with those obtained by other authors using different methods. The final part of the thesis investigates further the application of Chebyshdev methods to the solution of elliptic partial differential equations; an area where Chebyshev approximations have rarely been used. A poisson equation with a polynomial term is treated first followed by a slot problem in cylindrical geometry.

A critical assessment of the finite element method for calculating electric and magnetic fields

The present dissertation is concerned with the determination of the magnetic field distribution in ma[.rnetic electron lenses by means of the finite element method. In the differential form of this method a Poisson type equation is solved by numerical methods over a finite boundary. Previous methods of adapting this procedure to the requirements of digital computers have restricted its use to computers of extremely large core size. It is shown that by reformulating the boundary conditions, a considerable reduction in core store can be achieved for a given accuracy of field distribution. The magnetic field distribution of a lens may also be calculated by the integral form of the finite element rnethod. This eliminates boundary problems mentioned but introduces other difficulties. After a careful analysis of both methods it has proved possible to combine the advantages of both in a .new approach to the problem which may be called the 'differential-integral' finite element method. The application of this method to the determination of the magnetic field distribution of some new types of magnetic lenses is described. In the course of the work considerable re-programming of standard programs was necessary in order to reduce the core store requirements to a minimum.

The electron optical limits of performance of single-polepiece magnetic electron lenses

The thesis is concerned with the electron properties of single-polepiece magnetic electron lenses especially under conditions of extreme polepiece saturation. The electron optical properties are first analysed under conditions of high polepiece permeability. From this analysis, a general idea can be obtained of the important parameters that affect ultimate lens performance. In addition, useful information is obtained concerning the design of improved lenses operating under conditions of extreme polepiece saturation, for example at flux densities of the order of 10 Tesla. It is shown that in a single-polepiece lens , the position and shape of the lens exciting coil plays an important role. In particular, the maximum permissible current density in the windings,rather than the properties of the iron, can set a limit to lens performance. This factor was therefore investigated in some detail. The axial field distribution of a single-polepiece lens, unlike that of a conventional lens, is highly asymmetrical. There are therefore two possible physical arrangements of the lens with respect to the incoming electron beam. In general these two orientations will result in different aberration coefficients. This feature has also been investigated in some detail. Single-pole piece lenses are thus considerably more complicated electron- optically than conventional double polepiece lenses. In particular, the absence of the usual second polepiece causes most of the axial magnetic flux density distribution to lie outside the body of the lens. This can have many advantages in electron microscopy but it creates problems in calculating the magnetic field distribution. In particular, presently available computer programs are liable to be considerably in error when applied to such structures. It was therefore necessary to find independent ways of checking the field calculations. Furthermore, if the polepiece is allowed to saturate, much more calculation is involved since the field distribution becomes a non-linear function of the lens excitation. In searching for optimum lens designs, care was therefore taken to ensure that the coil was placed in the optimum position. If this condition is satisfied there seems to be no theoretical limit to the maximum flux density that can be attained at the polepiece tip. However , under iron saturation condition, some broadening of the axial field distribution will take place, thereby changing the lens aberrations . Extensive calculations were therefore made to find the minimum spherical and chromatic aberration coefficients . The focal properties of such lens designs are presented and compared with the best conventional double-polepiece lenses presently available.