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.

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.

Numerical evaluation of three dimensional effects in planar flow birefringence

The influence of three dimensional effects on isochromatic birefringence is evaluated for planar flows by means of numerical simulation. Two fluid models are investigated in channel and abrupt contraction geometries. In practice, the flows are confined by viewing windows, which alter the stresses along the optical path. The observed optical properties differ therefore from their counterpart in an ideal two-dimensional flow. To investigate the influence of these effects, the stress optical rule and the differential propagation Mueller matrix are used. The material parameters are selected so that a retardation of multiple orders is achieved, as is typical for highly birefringent melts. Errors due to three dimensional effects are mainly found on the symmetry plane, and increase significantly with the flow rate. Increasing the geometric aspect ratio improve the accuracy provided that the error on the retardation is less than one order. (C) 2004 Elsevier B.V. All rights reserved.

Effect of early marine diagenesis on coral reconstructions of surface-ocean 13C/12C and carbonate saturation state

Recent research suggests that future decreases in the carbonate saturation state of surface seawater associated with the projected build-up of atmospheric CO2 could cause a global decline in coral reef-building capacity. Whether significant reductions in coral calcification are underway is a matter of considerable debate. Multicentury records of skeletal calcification extracted from massive corals have the potential to reconstruct the progressive effect of anthropogenic changes in carbonate saturation on coral reefs. However, early marine aragonite cements are commonly precipitated from pore waters in the basal portions of massive coral skeletons and, if undetected, could result in apparent nonlinear reductions in coral calcification toward the present. To address this issue, we present records of coral skeletal density, extension rate, calcification rate, δ13C, and δ18O for well preserved and diagenetically altered coral cores spanning ∼1830-1994 A.D. at Ningaloo Reef Marine Park, Western Australia. The record for the pristine coral shows no significant decrease in skeletal density or δ13C indicative of anthropogenic changes in carbonate saturation state or δ13C of surface seawater (oceanic Suess effect). In contrast, progressive addition of early marine inorganic aragonite toward the base of the altered coral produces an apparent ∼25% decrease in skeletal density toward the present, which misleadingly matches the nonlinear twentieth century decrease in coral calcification predicted by recent modeling and experimental studies. In addition, the diagenetic aragonite is enriched in 13C, relative to coral aragonite, resulting in a nonlinear decrease in δ13C toward the present that mimics the decrease in δ13C expected from the oceanic Suess effect. Taken together, these diagenetic changes in skeletal density and δ13C could be misinterpreted to reflect changes in surface-ocean carbonate saturation state driven by the twentieth century build-up of atmospheric CO2. Copyright 2004 by the American Geophysical Union.

Modeling the development of acquired clinical immunity to Plasmodium falciparum malaria

Individuals living in regions where malaria is endemic develop an acquired immunity to malaria which enables them to remain asymptomatic while still carrying parasites. Field studies indicate that cumulative exposure to a variety of diverse Plasmodium parasites is required for the transition from symptomatic to asymptomatic malaria. This study used a simulation model of the within-host dynamics of P. falciparum to investigate the development of acquired clinical immunity under different transmission conditions and levels of parasite diversity. Antibodies developed to P. falciparum erythrocyte membrane protein 1 (PfEMP1), a clonally variant molecule, were assumed to be a key human immunological response to P. falciparum infection, along with responses to clonally conserved but polymorphic antigens. The time to the development of clinical immunity was found to be proportional to parasite diversity and inversely proportional to transmission intensity. The effect of early termination of symptomatic infections by chemotherapy was investigated and found not to inhibit the host's ability to develop acquired immunity. However, the time required to achieve this state was approximately double that compared to when no treatment was administered. This study demonstrates that an immune response primarily targeted against PfEMP1 has the ability to reduce clinical symptoms of infections irrespective of whether treatment is administered, supporting its role in the development of acquired clinical immunity. The results also illustrate a novel use for simulation models of P. falciparum infections, investigation of the influence of intervention strategies on the development of naturally acquired clinical immunity.

Comparing the concentration of calcium hydroxide in the diffuse double layer by modeling and experiment

The leaching of elements from the surface of charged fly ash particles is known to be an unsteady process. The mass transfer resistance provided by the diffuse double layer has been quantified as one of the reasons for this delayed leaching. In this work, a model based on mass transfer principles for predicting the concentration of calcium hydroxide in the diffuse double layer is presented. The significant difference between predicted calcium hydroxide concentration and the experimentally measured is explained.

Constriction flows of monodisperse linear entangled polymers: Multiscale modeling and flow visualization

We explore both the rheology and complex flow behavior of monodisperse polymer melts. Adequate quantities of monodisperse polymer were synthesized in order that both the materials rheology and microprocessing behavior could be established. In parallel, we employ a molecular theory for the polymer rheology that is suitable for comparison with experimental rheometric data and numerical simulation for microprocessing flows. The model is capable of matching both shear and extensional data with minimal parameter fitting. Experimental data for the processing behavior of monodisperse polymers are presented for the first time as flow birefringence and pressure difference data obtained using a Multipass Rheometer with an 11:1 constriction entry and exit flow. Matching of experimental processing data was obtained using the constitutive equation with the Lagrangian numerical solver, FLOWSOLVE. The results show the direct coupling between molecular constitutive response and macroscopic processing behavior, and differentiate flow effects that arise separately from orientation and stretch. (c) 2005 The Society of Rheology.

Modeling of adsorption in finite cylindrical pores by means of density functional theory

Adsorption of argon at its boiling point infinite cylindrical pores is considered by means of the non-local density functional theory (NLDFT) with a reference to MCM-41 silica. The NLDFT was adjusted to amorphous solids, which allowed us to quantitatively describe argon adsorption isotherm on nonporous reference silica in the entire bulk pressure range. In contrast to the conventional NLDFT technique, application of the model to cylindrical pores does not show any layering before the phase transition in conformity with experimental data. The finite pore is modeled as a cylindrical cavity bounded from its mouth by an infinite flat surface perpendicular to the pore axis. The adsorption of argon in pores of 4 and 5 nm diameters is analyzed in canonical and grand canonical ensembles using a two-dimensional version of NLDFT, which accounts for the radial and longitudinal fluid density distributions. The simulation results did not show any unusual features associated with accounting for the outer surface and support the conclusions obtained from the classical analysis of capillary condensation and evaporation. That is, the spontaneous condensation occurs at the vapor-like spinodal point, which is the upper limit of mechanical stability of the liquid-like film wetting the pore wall, while the evaporation occurs via a mechanism of receding of the semispherical meniscus from the pore mouth and the complete evaporation of the core occurs at the equilibrium transition pressure. Visualization of the pore filling and empting in the form of contour lines is presented.

Theoretical and numerical analysis of large-scale heat transfer problems with temperature-dependent pore-fluid densities

Purpose - In many scientific and engineering fields, large-scale heat transfer problems with temperature-dependent pore-fluid densities are commonly encountered. For example, heat transfer from the mantle into the upper crust of the Earth is a typical problem of them. The main purpose of this paper is to develop and present a new combined methodology to solve large-scale heat transfer problems with temperature-dependent pore-fluid densities in the lithosphere and crust scales. Design/methodology/approach - The theoretical approach is used to determine the thickness and the related thermal boundary conditions of the continental crust on the lithospheric scale, so that some important information can be provided accurately for establishing a numerical model of the crustal scale. The numerical approach is then used to simulate the detailed structures and complicated geometries of the continental crust on the crustal scale. The main advantage in using the proposed combination method of the theoretical and numerical approaches is that if the thermal distribution in the crust is of the primary interest, the use of a reasonable numerical model on the crustal scale can result in a significant reduction in computer efforts. Findings - From the ore body formation and mineralization points of view, the present analytical and numerical solutions have demonstrated that the conductive-and-advective lithosphere with variable pore-fluid density is the most favorite lithosphere because it may result in the thinnest lithosphere so that the temperature at the near surface of the crust can be hot enough to generate the shallow ore deposits there. The upward throughflow (i.e. mantle mass flux) can have a significant effect on the thermal structure within the lithosphere. In addition, the emplacement of hot materials from the mantle may further reduce the thickness of the lithosphere. Originality/value - The present analytical solutions can be used to: validate numerical methods for solving large-scale heat transfer problems; provide correct thermal boundary conditions for numerically solving ore body formation and mineralization problems on the crustal scale; and investigate the fundamental issues related to thermal distributions within the lithosphere. The proposed finite element analysis can be effectively used to consider the geometrical and material complexities of large-scale heat transfer problems with temperature-dependent fluid densities.

Application of a coupled ground-surface water flow model to simulate periodic groundwater flow influenced by a sloping boundary, capillarity and vertical flows

Theoretical developments as well as field and laboratory data have shown the influence of the capillary fringe on water table fluctuations to increase with the fluctuation frequency. The numerical solution of a full, partially saturated flow equation can be computationally expensive. In this paper, the influence of the capillary fringe on water table fluctuations is simplified through its parameterisation into the storage coefficient of a fully-saturated groundwater flow model using the complex effective porosity concept [Nielsen, P., Perrochet, P., 2000. Water table dynamics under capillary fringes: experiments and modelling. Advances in Water Resources 23 (1), 503-515; Nielsen, P., Perrochet, P., 2000. ERRATA: water table dynamics under capillary fringes: experiments and modelling (Advances in Water Resources 23 (2000) 503-515). Advances in Water Resources 23, 907-908]. The model is applied to sand flume observations of periodic water table fluctuations induced by simple harmonic forcing across a sloping boundary, analogous to many beach groundwater systems. While not providing information on the moisture distribution within the aquifer, this approach can reasonably predict the water table fluctuations in response to periodic forcing across a sloping boundary. Furthermore, he coupled ground-surface water model accurately predicts the extent of the seepage face formed at the sloping boundary. (C) 2005 Elsevier Ltd. All rights reserved.

A review of flow modeling for dense medium cyclones

A critical assessment is presented for the existing fluid flow models used for dense medium cyclones (DMCs) and hydrocyclones. As the present discussion indicates, the understanding of dense medium cyclone flow is still far from the complete. However, its similarity to the hydrocyclone provides a basis for improved understanding of fluid flow in DMCs. The complexity of fluid flow in DMCs is basically due to the existence of medium as well as the dominance of turbulent particle size and density effects on separation. Both the theoretical and experimental analysis is done with respect to two-phase motions and solid phase flow in hydrocyclones or DMCs. A detailed discussion is presented on the empirical, semiempirical, and the numerical models based upon both the vorticity-stream function approach and Navier-Stokes equations in their primitive variables and in cylindrical coordinates available in literature. The existing equations describing turbulence and multiphase flows in cyclone are also critically reviewed.

Modeling hydrogen separation in high temperature silica membrane systems

In this work, a working model is proposed of molecular sieve silica (MSS) multistage membrane systems for CO cleanup at high temperatures (up to 500 degrees C) in a simulated fuel cell fuel processing system. Gases are described as having little interactions with each other relative to the pore walls due to low isosteric heat of adsorption on silica surfaces and high temperatures. The Arrhenius function for activated transport of pure gases was used to predict mixture concentration in the permeate and retentate streams. Simulation predicted CO could be reduced to levels below the required 50 ppmv for polymer electrolyte membrane fuel cell anodes at a stage H-2/CO selectivity of higher than 40 in 4 series membrane units. Experimental validation showed predicting mixture concentrations required only pure gas permeation data. This model has significant application for setting industrial stretch targets and as a robust basis for complex membrane model configurations. (c) 2006 American Institute of Chemical Engineers.

Simulation of Brugada syndrome using cellular and three-dimensional whole-heart modeling approaches

Brugada syndrome (BS) is a genetic disease identified by an abnormal electrocardiogram ( ECG) ( mainly abnormal ECGs associated with right bundle branch block and ST-elevation in right precordial leads). BS can lead to increased risk of sudden cardiac death. Experimental studies on human ventricular myocardium with BS have been limited due to difficulties in obtaining data. Thus, the use of computer simulation is an important alternative. Most previous BS simulations were based on animal heart cell models. However, due to species differences, the use of human heart cell models, especially a model with three-dimensional whole-heart anatomical structure, is needed. In this study, we developed a model of the human ventricular action potential (AP) based on refining the ten Tusscher et al (2004 Am. J. Physiol. Heart Circ. Physiol. 286 H1573 - 89) model to incorporate newly available experimental data of some major ionic currents of human ventricular myocytes. These modified channels include the L-type calcium current (ICaL), fast sodium current (I-Na), transient outward potassium current (I-to), rapidly and slowly delayed rectifier potassium currents (I-Kr and I-Ks) and inward rectifier potassium current (I-Ki). Transmural heterogeneity of APs for epicardial, endocardial and mid-myocardial (M) cells was simulated by varying the maximum conductance of IKs and Ito. The modified AP models were then used to simulate the effects of BS on cellular AP and body surface potentials using a three-dimensional dynamic heart - torso model. Our main findings are as follows. (1) BS has little effect on the AP of endocardial or mid-myocardial cells, but has a large impact on the AP of epicardial cells. (2) A likely region of BS with abnormal cell AP is near the right ventricular outflow track, and the resulting ST-segment elevation is located in the median precordium area. These simulation results are consistent with experimental findings reported in the literature. The model can reproduce a variety of electrophysiological behaviors and provides a good basis for understanding the genesis of abnormal ECG under the condition of BS disease.

High-field magnetic resonance imaging with reduced field/tissue RF artefacts - A modeling study using hybrid MoM/FEM and FDTD technique

Due to complex field/tissue interactions, high-field magnetic resonance (MR) images suffer significant image distortions that result in compromised diagnostic quality. A new method that attempts to remove these distortions is proposed in this paper and is based on the use of transceiver-phased arrays. The proposed system uses, in the examples presented herein, a shielded four-element transceive-phased array head coil and involves performing two separate scans of the same slice with each scan using different excitations during transmission. By optimizing the amplitudes and phases for each scan, antipodal signal profiles can be obtained, and by combining both the images together, the image distortion can be reduced several fold. A combined hybrid method of moments (MoM)/finite element method (FEM) and finite-difference time-domain (FDTD) technique is proposed and used to elucidate the concept of the new method and to accurately evaluate the electromagnetic field (EMF) in a human head model. In addition, the proposed method is used in conjunction with the generalized auto-calibrating partially parallel acquisitions (GRAPPA) reconstruction technique to enable rapid imaging of the two scans. Simulation results reported herein for 11-T (470-MHz) brain imaging applications show that the new method with GRAPPA reconstruction theoretically results in improved image quality and that the proposed combined hybrid MoM/FEM and FDTD technique is. suitable for high-field magnetic resonance imaging (MRI) numerical analysis.

Modelling trickle irrigation: Comparison of analytical and numerical models for estimation of wetting front position with time

Irrigation practices that are profligate in their use of water have come under closer scrutiny by water managers and the public. Trickle irrigation has the propensity to increase water use efficiency but only if the system is designed to meet the soil and plant conditions. Recently we have provided a software tool, WetUp (http://www.clw.csiro.au/products/wetup/), to calculate the wetting patterns from trickle irrigation emitters. WetUp uses an analytical solution to calculate the wetted perimeter for both buried and surface emitters. This analytical solution has a number of assumptions, two of which are that the wetting front is defined by water content at which the hydraulic conductivity (K) is I mm day(-1) and that the flow occurs from a point source. Here we compare the wetting patterns calculated with a 2-dimensional numerical model, HYDRUS2D, for solving the water flow into typical soils with the analytical solution. The results show that the wetting patterns are similar, except when the soil properties result in the assumption of a point source no longer being a good description of the flow regime. Difficulties were also experienced with getting stable solutions with HYDRUS2D for soils with low hydraulic conductivities. (c) 2005 Elsevier Ltd. All rights reserved.