An FDTD model for calculation of gradient-induced eddy currents in MRI system

In most magnetic resonance imaging (MRI) systems, pulsed magnetic gradient fields induce eddy currents in the conducting structures of the superconducting magnet. The eddy currents induced in structures within the cryostat are particularly problematic as they are characterized by long time constants by virtue of the low resistivity of the conductors. This paper presents a three-dimensional (3-D) finite-difference time-domain (FDTD) scheme in cylindrical coordinates for eddy-current calculation in conductors. This model is intended to be part of a complete FDTD model of an MRI system including all RF and low-frequency field generating units and electrical models of the patient. The singularity apparent in the governing equations is removed by using a series expansion method and the conductor-air boundary condition is handled using a variant of the surface impedance concept. The numerical difficulty due to the asymmetry of Maxwell equations for low-frequency eddy-current problems is circumvented by taking advantage of the known penetration behavior of the eddy-current fields. A perfectly matched layer absorbing boundary condition in 3-D cylindrical coordinates is also incorporated. The numerical method has been verified against analytical solutions for simple cases. Finally, the algorithm is illustrated by modeling a pulsed field gradient coil system within an MRI magnet system. The results demonstrate that the proposed FDTD scheme can be used to calculate large-scale eddy-current problems in materials with high conductivity at low frequencies.

Analysis of transient eddy currents in MRI using a cylindrical FDTD method

Most magnetic resonance imaging (MRI) spatial encoding techniques employ low-frequency pulsed magnetic field gradients that undesirably induce multiexponentially decaying eddy currents in nearby conducting structures of the MRI system. The eddy currents degrade the switching performance of the gradient system, distort the MRI image, and introduce thermal loads in the cryostat vessel and superconducting MRI components. Heating of superconducting magnets due to induced eddy currents is particularly problematic as it offsets the superconducting operating point, which can cause a system quench. A numerical characterization of transient eddy current effects is vital for their compensation/control and further advancement of the MRI technology as a whole. However, transient eddy current calculations are particularly computationally intensive. In large-scale problems, such as gradient switching in MRI, conventional finite-element method (FEM)-based routines impose very large computational loads during generation/solving of the system equations. Therefore, other computational alternatives need to be explored. This paper outlines a three-dimensional finite-difference time-domain (FDTD) method in cylindrical coordinates for the modeling of low-frequency transient eddy currents in MRI, as an extension to the recently proposed time-harmonic scheme. The weakly coupled Maxwell's equations are adapted to the low-frequency regime by downscaling the speed of light constant, which permits the use of larger FDTD time steps while maintaining the validity of the Courant-Friedrich-Levy stability condition. The principal hypothesis of this work is that the modified FDTD routine can be employed to analyze pulsed-gradient-induced, transient eddy currents in superconducting MRI system models. The hypothesis is supported through a verification of the numerical scheme on a canonical problem and by analyzing undesired temporal eddy current effects such as the B-0-shift caused by actively shielded symmetric/asymmetric transverse x-gradient head and unshielded z-gradient whole-body coils operating in proximity to a superconducting MRI magnet.

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.

Tailoring magnetic field gradient design to magnet cryostat geometry

Eddy currents induced within a magnetic resonance imaging (MRI) cryostat bore during pulsing of gradient coils can be applied constructively together with the gradient currents that generate them, to obtain good quality gradient uniformities within a specified imaging volume over time. This can be achieved by simultaneously optimizing the spatial distribution and temporal pre-emphasis of the gradient coil current, to account for the spatial and temporal variation of the secondary magnetic fields due to the induced eddy currents. This method allows the tailored design of gradient coil/magnet configurations and consequent engineering trade-offs. To compute the transient eddy currents within a realistic cryostat vessel, a low-frequency finite-difference time-domain (FDTD) method using total-field scattered-field (TFSF) scheme has been performed and validated

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.

Crystal structure and magnetic behaviour of a five-coordinate iron(III) complex of pyridoxal-4-methylthiosemicarbazone

The crystal structure and magnetic properties of a penta-coordinate iron(III) complex of pyridoxal-4-methylthiosemicarbazone, [Fe(Hmthpy)Cl](CHCHSO), are reported. The synthesised ligand and the metal complex were characterised by spectroscopic methods (H NMR, IR, and mass spectroscopy), elemental analysis, and single crystal X-ray diffraction. The complex crystallises as dark brown microcrystals. The crystal data determined at 100(1) K revealed a triclinic system, space group P over(1, ¯) (Z = 2). The ONSCl geometry around the iron(III) atom is intermediate between trigonal bipyramidal and square pyramidal (t = 0.40). The temperature dependence of the magnetic susceptibility (5-300 K) is consistent with a high spin Fe(III) ion (S = 5/2) exhibiting zero-field splitting. Interpretation of these data yielded: D = 0.34(1) cm and g = 2.078(3). © 2007 Elsevier B.V. All rights reserved.

State detection of bond wires in IGBT modules using eddy current pulsed thermography

Insulated gate bipolar transistor (IGBT) modules are important safety critical components in electrical power systems. Bond wire lift-off, a plastic deformation between wire bond and adjacent layers of a device caused by repeated power/thermal cycles, is the most common failure mechanism in IGBT modules. For the early detection and characterization of such failures, it is important to constantly detect or monitor the health state of IGBT modules, and the state of bond wires in particular. This paper introduces eddy current pulsed thermography (ECPT), a nondestructive evaluation technique, for the state detection and characterization of bond wire lift-off in IGBT modules. After the introduction of the experimental ECPT system, numerical simulation work is reported. The presented simulations are based on the 3-D electromagnetic-thermal coupling finite-element method and analyze transient temperature distribution within the bond wires. This paper illustrates the thermal patterns of bond wires using inductive heating with different wire statuses (lifted-off or well bonded) under two excitation conditions: nonuniform and uniform magnetic field excitations. Experimental results show that uniform excitation of healthy bonding wires, using a Helmholtz coil, provides the same eddy currents on each, while different eddy currents are seen on faulty wires. Both experimental and numerical results show that ECPT can be used for the detection and characterization of bond wires in power semiconductors through the analysis of the transient heating patterns of the wires. The main impact of this paper is that it is the first time electromagnetic induction thermography, so-called ECPT, has been employed on power/electronic devices. Because of its capability of contactless inspection of multiple wires in a single pass, and as such it opens a wide field of investigation in power/electronic devices for failure detection, performance characterization, and health monitoring.

Studies of turbulent flows in continuous casting of steel with and without magnetic field

This thesis develops and tests various transient and steady-state computational models such as direct numerical simulation (DNS), large eddy simulation (LES), filtered unsteady Reynolds-averaged Navier-Stokes (URANS) and steady Reynolds-averaged Navier-Stokes (RANS) with and without magnetic field to investigate turbulent flows in canonical as well as in the nozzle and mold geometries of the continuous casting process. The direct numerical simulations are first performed in channel, square and 2:1 aspect rectangular ducts to investigate the effect of magnetic field on turbulent flows. The rectangular duct is a more practical geometry for continuous casting nozzle and mold and has the option of applying magnetic field either perpendicular to broader side or shorter side. This work forms the part of a graphic processing unit (GPU) based CFD code (CU-FLOW) development for magnetohydrodynamic (MHD) turbulent flows. The DNS results revealed interesting effects of the magnetic field and its orientation on primary, secondary flows (instantaneous and mean), Reynolds stresses, turbulent kinetic energy (TKE) budgets, momentum budgets and frictional losses, besides providing DNS database for two-wall bounded square and rectangular duct MHD turbulent flows. Further, the low- and high-Reynolds number RANS models (k-ε and Reynolds stress models) are developed and tested with DNS databases for channel and square duct flows with and without magnetic field. The MHD sink terms in k- and ε-equations are implemented as proposed by Kenjereš and Hanjalić using a user defined function (UDF) in FLUENT. This work revealed varying accuracies of different RANS models at different levels. This work is useful for industry to understand the accuracies of these models, including continuous casting. After realizing the accuracy and computational cost of RANS models, the steady-state k-ε model is then combined with the particle image velocimetry (PIV) and impeller probe velocity measurements in a 1/3rd scale water model to study the flow quality coming out of the well- and mountain-bottom nozzles and the effect of stopper-rod misalignment on fluid flow. The mountain-bottom nozzle was found more prone to the longtime asymmetries and higher surface velocities. The left misalignment of stopper gave higher surface velocity on the right leading to significantly large number of vortices forming behind the nozzle on the left. Later, the transient and steady-state models such as LES, filtered URANS and steady RANS models are combined with ultrasonic Doppler velocimetry (UDV) measurements in a GaInSn model of typical continuous casting process. LES-CU-LOW is the fastest and the most accurate model owing to much finer mesh and a smaller timestep. This work provided a good understanding on the performance of these models. The behavior of instantaneous flows, Reynolds stresses and proper orthogonal decomposition (POD) analysis quantified the nozzle bottom swirl and its importance on the turbulent flow in the mold. Afterwards, the aforementioned work in GaInSn model is extended with electromagnetic braking (EMBr) to help optimize a ruler-type brake and its location for the continuous casting process. The magnetic field suppressed turbulence and promoted vortical structures with their axis aligned with the magnetic field suggesting tendency towards 2-d turbulence. The stronger magnetic field at the nozzle well and around the jet region created large scale and lower frequency flow behavior by suppressing nozzle bottom swirl and its front-back alternation. Based on this work, it is advised to avoid stronger magnetic field around jet and nozzle bottom to get more stable and less defect prone flow.

Magnetic multilayer fabrication technology with selective activation of SU-8 films

Developing magnetic multilayers are essential for reducing the core eddy current losses in the integrated power magnetic components (inductors/transformers). PVD based processes are typically used to achieve the multilayers with thin dielectric spacers. However, those processes are costly, and can be difficult to integrate. It is evident that cost effective alternative is needed. In recent years, electrochemical processes have been investigated to address these issues. One such method would be to successive metallization of insulating photoresists acting as spacer layer (such as SU-8) with soft magnetic films (such as Ni-Fe-Co alloys). This paper describes an experimental procedure to fabricate magnetic multilayers with a thin variant of SU-8 2 (< 1.5 µm) as inter-layers for integrated micro-inductors/transformers for power conversion applications.

Structural Connectivity Of The Default Mode Network And Cognition In Alzheimer׳s Disease.

Disconnectivity between the Default Mode Network (DMN) nodes can cause clinical symptoms and cognitive deficits in Alzheimer׳s disease (AD). We aimed to examine the structural connectivity between DMN nodes, to verify the extent in which white matter disconnection affects cognitive performance. MRI data of 76 subjects (25 mild AD, 21 amnestic Mild Cognitive Impairment subjects and 30 controls) were acquired on a 3.0T scanner. ExploreDTI software (fractional Anisotropy threshold=0.25 and the angular threshold=60°) calculated axial, radial, and mean diffusivities, fractional anisotropy and streamline count. AD patients showed lower fractional anisotropy (P=0.01) and streamline count (P=0.029), and higher radial diffusivity (P=0.014) than controls in the cingulum. After correction for white matter atrophy, only fractional anisotropy and radial diffusivity remained significantly lower in AD compared to controls (P=0.003 and P=0.05). In the parahippocampal bundle, AD patients had lower mean and radial diffusivities (P=0.048 and P=0.013) compared to controls, from which only radial diffusivity survived for white matter adjustment (P=0.05). Regression models revealed that cognitive performance is also accounted for by white matter microstructural values. Structural connectivity within the DMN is important to the execution of high-complexity tasks, probably due to its relevant role in the integration of the network.

Cardiac Magnetic Resonance Assessment Of Interstitial Myocardial Fibrosis And Cardiomyocyte Hypertrophy In Hypertensive Mice Treated With Spironolactone.

Nearly 50% of patients with heart failure (HF) have preserved LV ejection fraction, with interstitial fibrosis and cardiomyocyte hypertrophy as early manifestations of pressure overload. However, methods to assess both tissue characteristics dynamically and noninvasively with therapy are lacking. We measured the effects of mineralocorticoid receptor blockade on tissue phenotypes in LV pressure overload using cardiac magnetic resonance (CMR). Mice were randomized to l-nitro-ω-methyl ester (l-NAME, 3 mg/mL in water; n=22), or l-NAME with spironolactone (50 mg/kg/day in subcutaneous pellets; n=21). Myocardial extracellular volume (ECV; marker of diffuse interstitial fibrosis) and the intracellular lifetime of water (τic; marker of cardiomyocyte hypertrophy) were determined by CMR T1 imaging at baseline and after 7 weeks of therapy alongside histological assessments. Administration of l-NAME induced hypertensive heart disease in mice, with increases in mean arterial pressure, LV mass, ECV, and τic compared with placebo-treated controls, while LV ejection fraction was preserved (>50%). In comparison, animals receiving both spironolactone and l-NAME (l-NAME+S) showed less concentric remodeling, and a lower myocardial ECV and τic, indicating decreased interstitial fibrosis and cardiomyocyte hypertrophy (ECV: 0.43 ± 0.09 for l-NAME versus 0.25 ± 0.03 for l-NAME+S, P<0.001; τic: 0.42 ± 0.11 for l-NAME groups versus 0.12 ± 0.05 for l-NAME+S group). Mice treated with a combination of l-NAME and spironolactone were similar to placebo-treated controls at 7 weeks. Spironolactone attenuates interstitial fibrosis and cardiomyocyte hypertrophy in hypertensive heart disease. CMR can phenotype myocardial tissue remodeling in pressure-overload, furthering our understanding of HF progression.

Emergent Antiferromagnetism Out Of The Hidden-order" State In Uru2si2: High Magnetic Field Nuclear Magnetic Resonance To 40 T."

Very high field (29)Si-NMR measurements using a fully (29)Si-enriched URu(2)Si(2) single crystal were carried out in order to microscopically investigate the hidden order (HO) state and adjacent magnetic phases in the high field limit. At the lowest measured temperature of 0.4 K, a clear anomaly reflecting a Fermi surface instability near 22 T inside the HO state is detected by the (29)Si shift, (29)K(c). Moreover, a strong enhancement of (29)K(c) develops near a critical field H(c) ≃ 35.6 T, and the ^{29}Si-NMR signal disappears suddenly at H(c), indicating the total suppression of the HO state. Nevertheless, a weak and shifted (29)Si-NMR signal reappears for fields higher than H(c) at 4.2 K, providing evidence for a magnetic structure within the magnetic phase caused by the Ising-type anisotropy of the uranium ordered moments.

Feasibility And Reproducibility Of Diffusion-weighted Magnetic Resonance Imaging Of The Fetal Brain In Twin-twin Transfusion Syndrome.

The aim of this study is to test the feasibility and reproducibility of diffusion-weighted magnetic resonance imaging (DW-MRI) evaluations of the fetal brains in cases of twin-twin transfusion syndrome (TTTS). From May 2011 to June 2012, 24 patients with severe TTTS underwent MRI scans for evaluation of the fetal brains. Datasets were analyzed offline on axial DW images and apparent diffusion coefficient (ADC) maps by two radiologists. The subjective evaluation was described as the absence or presence of water diffusion restriction. The objective evaluation was performed by the placement of 20-mm(2) circular regions of interest on the DW image and ADC maps. Subjective interobserver agreement was assessed by the kappa correlation coefficient. Objective intraobserver and interobserver agreements were assessed by proportionate Bland-Altman tests. Seventy-four DW-MRI scans were performed. Sixty of them (81.1%) were considered to be of good quality. Agreement between the radiologists was 100% for the absence or presence of diffusion restriction of water. For both intraobserver and interobserver agreement of ADC measurements, proportionate Bland-Altman tests showed average percentage differences of less than 1.5% and 95% CI of less than 18% for all sites evaluated. Our data demonstrate that DW-MRI evaluation of the fetal brain in TTTS is feasible and reproducible.

Beam-energy Dependence Of Charge Separation Along The Magnetic Field In Au+au Collisions At Rhic.

Local parity-odd domains are theorized to form inside a quark-gluon plasma which has been produced in high-energy heavy-ion collisions. The local parity-odd domains manifest themselves as charge separation along the magnetic field axis via the chiral magnetic effect. The experimental observation of charge separation has previously been reported for heavy-ion collisions at the top RHIC energies. In this Letter, we present the results of the beam-energy dependence of the charge correlations in Au+Au collisions at midrapidity for center-of-mass energies of 7.7, 11.5, 19.6, 27, 39, and 62.4 GeV from the STAR experiment. After background subtraction, the signal gradually reduces with decreased beam energy and tends to vanish by 7.7 GeV. This implies the dominance of hadronic interactions over partonic ones at lower collision energies.