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.

On the steady-state assumption and its application to the rotating disk voltammetry of adsorbed enzymes

Rotating disk voltammetry is routinely used to study electrochemically driven enzyme catalysis because of the assumption that the method produces a steady-state system. This assumption is based on the sigmoidal shape of the voltammograms. We have introduced an electrochemical adaptation of the King-Altman method to simulate voltammograms in which the enzyme catalysis, within an immobilized enzyme layer, is steadystate. This method is readily adaptable to any mechanism and provides a readily programmable means of obtaining closed form analytical equations for a steady-state system. The steady-state simulations are compared to fully implicit finite difference (FIFD) simulations carried out without any steady-state assumptions. On the basis of our simulations, we conclude that, under typical experimental conditions, steady-state enzyme catalysis is unlikely to occur within electrode-immobilized enzyme layers and that typically sigmoidal rotating disk voltammograms merely reflect a mass transfer steady state as opposed to a true steady state of enzyme intermediates at each potential.

Effective cell-centred time-domain Maxwell's equations numerical solvers

This research work analyses techniques for implementing a cell-centred finite-volume time-domain (ccFV-TD) computational methodology for the purpose of studying microwave heating. Various state-of-the-art spatial and temporal discretisation methods employed to solve Maxwell's equations on multidimensional structured grid networks are investigated, and the dispersive and dissipative errors inherent in those techniques examined. Both staggered and unstaggered grid approaches are considered. Upwind schemes using a Riemann solver and intensity vector splitting are studied and evaluated. Staggered and unstaggered Leapfrog and Runge-Kutta time integration methods are analysed in terms of phase and amplitude error to identify which method is the most accurate and efficient for simulating microwave heating processes. The implementation and migration of typical electromagnetic boundary conditions. from staggered in space to cell-centred approaches also is deliberated. In particular, an existing perfectly matched layer absorbing boundary methodology is adapted to formulate a new cell-centred boundary implementation for the ccFV-TD solvers. Finally for microwave heating purposes, a comparison of analytical and numerical results for standard case studies in rectangular waveguides allows the accuracy of the developed methods to be assessed. © 2004 Elsevier Inc. All rights reserved.

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.

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.

Stochastic modelling and simulations for solute transport in porous media

A stochastic model for solute transport in aquifers is studied based on the concepts of stochastic velocity and stochastic diffusivity. By applying finite difference techniques to the spatial variables of the stochastic governing equation, a system of stiff stochastic ordinary differential equations is obtained. Both the semi-implicit Euler method and the balanced implicit method are used for solving this stochastic system. Based on the Karhunen-Loeve expansion, stochastic processes in time and space are calculated by means of a spatial correlation matrix. Four types of spatial correlation matrices are presented based on the hydraulic properties of physical parameters. Simulations with two types of correlation matrices are presented.

An object-orient designed FDTD simulator: Applications to high field systems

An object-oriented finite-difference time-domain (FDTD) simulator has been developed for electromagnetic study and design applications in Magnetic Resonance Imaging. It is aimed to be a complete FDTD model of an MRI system including all high and low-frequency field generating units and electrical models of the patient. The design method is described and MRI-based numerical examples are presented to illustrate the function of the numerical solver, particular emphasis is placed on high field studies.

Computational modeling of optical trapping of vaterite microspheres

The ability to grow microscopic spherical birefringent crystals of vaterite, a calcium carbonate mineral, has allowed the development of an optical microrheometer based on optical tweezers. However, since these crystals are birefringent, and worse, are expected to have non-uniform birefringence, computational modeling of the microrheometer is a highly challenging task. Modeling the microrheometer - and optical tweezers in general - typically requires large numbers of repeated calculations for the same trapped particle. This places strong demands on the efficiency of computational methods used. While our usual method of choice for computational modelling of optical tweezers - the T-matrix method - meets this requirement of efficiency, it is restricted to homogeneous isotropic particles. General methods that can model complex structures such as the vaterite particles, such as finite-difference time-domain (FDTD) or finite-difference frequency-domain (FDFD) methods, are inefficient. Therefore, we have developed a hybrid FDFD/T-matrix method that combines the generality of volume-discretisation methods such as FDFD with the efficiency of the T-matrix method. We have used this hybrid method to calculate optical forces and torques on model vaterite spheres in optical traps. We present and compare the results of computational modelling and experimental measurements.

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

FDTD study of an UWB radar technique for breast tumor detection

A finite difference time domain (FDTD) method is applied to investigate capabilities of an ultra-wide band (UWB) radar system to detect a breast tumor. The first part of the investigations concerns FDTD simulations of a phantom formed by a plastic container with liquid and a small reflecting target. The second part focuses on a three-dimensional numerical breast model with a small tumor. FDTD simulations are carried out assuming a planar incident wave. Various time snap shots of the electromagnetic field are recorded to learn about the physical phenomenon of reflection and scattering in different layers of the phantom.

The steady performance of a porous and compliant aerostatic thrust bearing

Recent developments in aerostatic thrust bearings have included: (a) the porous aerostatic thrust bearing containing a porous pad and (b) the inherently compensated compliant surface aerostatic thrust bearing containing a thin elastomer layer. Both these developments have been reported to improve the bearing load capacity compared to conventional aerostatic thrust bearings with rigid surfaces. This development is carried one stage further in a porous and compliant aerostatic thrust bearing incorporating both a porous pad and an opposing compliant surface. The thin elastomer layer forming the compliant surface is bonded to a rigid backing and is of a soft rubber like material. Such a bearing is studied experimentally and theoretically under steady state operating conditions. A mathematical model is presented to predict the bearing performance. In this model is a simplified solution to the elasticity equations for deflections of the compliant surface. Account is also taken of deflections in the porous pad due to the pressure difference across its thickness. The lubrication equations for flow in the porous pad and bearing clearance are solved by numerical finite difference methods. An iteration procedure is used to couple deflections of the compliant surface and porous pad with solutions to the lubrication equations. Comparisons between experimental results and theoretically predicted bearing performance are in good agreement. However these results show that the porous and compliant aerostatic thrust bearing performance is lower than that of a porous aerostatic thrust bearing with a rigid surface in place of the compliant surface. This discovery is accounted to the recess formed in the bearing clearance by deflections of the compliant surface and its effect on flow through the porous pad.

The strength of column-to-foundation joints in reinforced concrete

This thesis is concerned with the experimental and theoretical investigation into the compression bond of column longitudinal reinforcement in the transference of axial load from a reinforced concrete column to a base. Experimental work includes twelve tests with square twisted bars and twenty four tests with ribbed bars. The effects of bar size, anchorage length in the base, plan area of the base, provision of bae tensile reinforcement, links around the column bars in the base, plan area of column and concrete compressive strength were investigated in the tests. The tests indicated that the strength of the compression anchorage of deformed reinforcing steel in the concrete was primarily dependent on the concrete strength and the resistance to bursting, which may be available within the anchorage . It was shown in the tests without concreted columns that due to a large containment over the bars in the foundation, failure occurred due to the breakdown of bond followed by the slip of the column bars along the anchorage length. The experimental work showed that the bar size , the stress in the bar, the anchorage length, provision of the transverse steel and the concrete compressive strength significantly affect the bond stress at failure. The ultimate bond stress decreases as the anchorage length is increased, while the ultimate bond stress increases with increasing each of the remainder parameters. Tests with concreted columns also indicated that a section of the column contributed to the bond length in the foundation by acting as an extra anchorage length. The theoretical work is based on the Mindlin equation( 3), an analytical method used in conjunction with finite difference calculus. The theory is used to plot the distribution of bond stress in the elastic and the elastic-plastic stage of behaviour. The theory is also used to plot the load-vertical displacement relationship of the column bars in the anchorage length, and also to determine the theoretical failure load of foundation. The theoretical solutions are in good agreement with the experimental results and the distribution of bond stress is shown to be significantly influenced by the bar stiffness factor K. A comparison of the experimental results with the current codes shows that the bond stresses currently used are low and in particular, CPIlO(56) specifies very conservative design bond stresses .

A Note on the “Constructing” of Nonstationary Methods for Solving Nonlinear Equations with Raised Speed of Convergence

This paper is partially supported by project ISM-4 of Department for Scientific Research, “Paisii Hilendarski” University of Plovdiv.

Optimal Control of a Second Order Parabolic Heat Equation

In this paper, we are concerned with the optimal control boundary control of a second order parabolic heat equation. Using the results in [Evtushenko, 1997] and spatial central finite difference with diagonally implicit Runge-Kutta method (DIRK) is applied to solve the parabolic heat equation. The conjugate gradient method (CGM) is applied to solve the distributed control problem. Numerical results are reported.

A Gradient-Type Optimization Technique for the Optimal Control for Schrodinger Equations

In this paper, we are considered with the optimal control of a schrodinger equation. Based on the formulation for the variation of the cost functional, a gradient-type optimization technique utilizing the finite difference method is then developed to solve the constrained optimization problem. Finally, a numerical example is given and the results show that the method of solution is robust.