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.

Characterisation and environmental application of an Australian natural zeolite for basic dye removal from aqueous solution

An Australian natural zeolite was collected, characterised and employed for basic dye adsorption in aqueous solution. The natural zeolite is mainly composed of clinoptiloite, quartz and mordenite and has cation-exchange capacity of 120 meq/100 g. The natural zeolite presents higher adsorption capacity for methylene blue than rhodamine B with the maximal adsorption capacity of 2.8 x 10(-5) and 7.9 x 10(-5) Mot/g at 50 degrees C for rhodamine B and methylene blue, respectively. Kinetic studies indicated that the adsorption followed the pseudo second-order kinetics and could be described as two-stage diffusion process. The adsorption isotherm could be fitted by the Langmuir and Freundlich models. Thermodynamic calculations showed that the adsorption is endothermic process with Delta H degrees at 2.0 and 8.7 kJ/mol for rhodamine B and methylene blue. It has also found that the regenerated zeolites by high-temperature calcination and Fenton oxidation showed similar adsorption capacity but lower than the fresh sample. Only 60% capacity could be recovered by the two regeneration techniques. (c) 2006 Elsevier B.V. All rights reserved.

Comparison of blast furnace raceway size with theory

A previously developed one-dimensional mathematical model, to explain raceway hysteresis, is used to predict the raceway diameter in operating blast furnaces and hot models. Raceway size obtained from the open literature under various conditions for various blast furnaces are compared with computed predictions. In addition the predictions are also compared with published outcomes from other hot models. Simulated results on raceway diameter are in very good agreement with published operating blast furnace and hot model data. The effect of various parameters such as tuyere and hearth diameter, coke size and density, void fraction and bed height on raceway diameter has been studied.

Comparison of cellulose solubilisation rates in rumen and landfill leachate inoculated reactors

The aim of this study was to conduct a number of controlled digestions to obtain easily comparable cellulose solubilisation rates and to compare these rates to those found in the literature to see which operational differences were significant in affecting cellulose degradation during anaerobic digestion. The results suggested that differences in volumetric cellulose solubilisation rates were not indicative of the true performance of cellulose digestion systems. When cellulose solubilisation rates were normalised by the mass of cellulose in the reactor at each time step, the comparison of the rates became more meaningful. Cellulose solubilisation was surface area limited. Therefore, changes in the loading rate of cellulose to the reactor altered the volumetric solubilisation rate without changing the mass normalised rate. Comparison of mass normalised solubilisation rates from this study and the literature demonstrated that differences in reactor configuration and operational conditions did not significantly impact on the solubilisation rate whereas the difference in composition of the microbial communities showed a marked effect. This work highlights the importance of using appropriately normalised data when making comparisons between systems with differing operational conditions. (c) 2005 Elsevier Ltd. All rights reserved.

Control of culture environment for improved polyethylenimine-mediated transient production of recombinant monoclonal antibodies by CHO cells

In this study we describe optimization of polyethylenimine (PEI)-mediated transient production of recombinant protein by CHO cells by facile manipulation of a chemically defined culture environment to limit accumulation of nonproductive cell biomass, increase the duration of recombinant protein production from transfected plasmid DNA, and increase cell-specific production. The optimal conditions for transient transfection of suspension-adapted CHO cells using branched, 25 kDa PEI as a gene delivery vehicle were experimentally determined by production of secreted alkaline phosphatase reporter in static cultures and recombinant IgG(4) monoclonal antibody (Mab) production in agitated shake flask cultures to be a DNA concentration of 1.25 mu g 10(6) cells(-1) mL(-1) at a PEI nitrogen: DNA phosphate ratio of 20:1. These conditions represented the optimal compromise between PEI cytotoxicity and product yield with most efficient recombinant DNA utilization. Separately, both addition of recombinant insulin-like growth factor (LR3-IGF) and a reduction in culture temperature to 32 degrees C were found to increase product titer 2- and 3-fold, respectively. However, mild hypothermia and LR3-IGF acted synergistically to increase product titer 11-fold. Although increased product titer in the presence of LR3-IGF alone was solely a consequence of increased culture duration, a reduction in culture temperature post-transfection increased both the integral of viable cell concentration (IVC) and cell-specific Mab production rate. For cultures maintained at 32 degrees C in the presence of LR3-IGF, IVC and qMab were increased 4- and 2.5-fold, respectively. To further increase product yield from transfected DNA, the duration of transgene expression in cell populations maintained at 32 C in the presence of LR3-IGF was doubled by periodic resuspension of transfected cells in fresh media, leading to a 3-fold increase in accumulated Mab titer from similar to 13 to similar to 39 mg L-1. Under these conditions, Mab glycosylation at Asn297 remained essentially constant and similar to that of the same Mab produced by stably transfected GS-CHO cells. From these data we suggest that the efficiency of transient production processes (protein output per rDNA input) can be significantly improved using a combination of mild hypothermia and growth factor(s) to yield an extended activated hypothermic synthesis.

Effect of adsorption deformation on thermodynamic characteristics of a fluid in slit pores at sub-critical conditions

Solvation. pressure due to adsorption of fluids in porous materials is the cause of elastic deformation of an adsorbent, which is accessible to direct experimental measurements. Such a deformation contributes to the Helmholtz free energy of the whole adsorbent-adsorbate system due to accumulation of compression or tension energy by the solid. It means that in the general case the solid has to be considered as not solely a source of the external potential field for the fluid confined in the pore volume, but also as thermodynamically nonmert component of the solid-fluid system. We present analysis of nitrogen adsorption isotherms and heat of adsorption in slit graphitic pores accounting for the adsorption deformation by means of nonlocal density functional theory. (c) 2006 Elsevier Ltd. All rights reserved.

Effects of rock properties on specific cutting energy in linear cutting of sandstones by picks

Specific cutting energy (SE) has been widely used to assess the rock cuttability for mechanical excavation purposes. Some prediction models were developed for SE through correlating rock properties with SE values. However, some of the textural and compositional rock parameters i.e. texture coefficient and feldspar, mafic, and felsic mineral contents were not considered. The present study is to investigate the effects of previously ignored rock parameters along with engineering rock properties on SE. Mineralogical and petrographic analyses, rock mechanics, and linear rock cutting tests were performed on sandstone samples taken from sites around Ankara, Turkey. Relationships between SE and rock properties were evaluated using bivariate correlation and linear regression analyses. The tests and subsequent analyses revealed that the texture coefficient and feldspar content of sandstones affected rock cuttability, evidenced by significant correlations between these parameters and SE at a 90% confidence level. Felsic and mafic mineral contents of sandstones did not exhibit any statistically significant correlation against SE. Cementation coefficient, effective porosity, and pore volume had good correlations against SE. Poisson's ratio, Brazilian tensile strength, Shore scleroscope hardness, Schmidt hammer hardness, dry density, and point load strength index showed very strong linear correlations against SE at confidence levels of 95% and above, all of which were also found suitable to be used in predicting SE individually, depending on the results of regression analysis, ANOVA, Student's t-tests, and R-2 values. Poisson's ratio exhibited the highest correlation with SE and seemed to be the most reliable SE prediction tool in sandstones.

Efficient and stable nitritation and denitritation of ammonium-rich sludge dewatering liquor using an SBR with continuous loading

Separate treatment of dewatering liquor from anaerobic sludge digestion significantly reduces the nitrogen load of the main stream and improves overall nitrogen elimination. Such ammonium-rich wastewater is particularly suited to be treated by high rate processes which achieve a rapid elimination of nitrogen with a minimal COD requirement. Processes whereby ammonium is oxidised to nitrite only (nitritation) followed by denitritation with carbon addition can achieve this. Nitrogen removal by nitritation/denitritation was optimised using a novel SBR operation with continuous dewatering liquor addition. Efficient and robust nitrogen elimination was obtained at a total hydraulic retention time of 1 day via the nitrite pathway. Around 85-90% nitrogen removal was achieved at an ammonium loading rate of 1.2 g NH4+-N m(-3) d(-1). Ethanol was used as electron donor for denitritation at a ratio of 2.2gCODg(-1) N removed. Conventional nitritation/denitritation with rapid addition of the dewatering liquor at the beginning of the cycle often resulted in considerable nitric oxide (NO) accumulation during the anoxic phase possibly leading to unstable denitritation. Some NO production was still observed in the novel continuous mode, but denitritation was never seriously affected. Thus, process stability can be increased and the high specific reaction rates as well as the continuous feeding result in decreased reactor size for full-scale operation. (c) 2006 Elsevier Ltd. All rights reserved.

Electrodynamic heart model construction and ECG simulation

Objectives: In this paper, we present a unified electrodynamic heart model that permits simulations of the body surface potentials generated by the heart in motion. The inclusion of motion in the heart model significantly improves the accuracy of the simulated body surface potentials and therefore also the 12-lead ECG. Methods: The key step is to construct an electromechanical heart model. The cardiac excitation propagation is simulated by an electrical heart model, and the resulting cardiac active forces are used to calculate the ventricular wall motion based on a mechanical model. The source-field point relative position changes during heart systole and diastole. These can be obtained, and then used to calculate body surface ECG based on the electrical heart-torso model. Results: An electromechanical biventricular heart model is constructed and a standard 12-lead ECG is simulated. Compared with a simulated ECG based on the static electrical heart model, the simulated ECG based on the dynamic heart model is more accordant with a clinically recorded ECG, especially for the ST segment and T wave of a V1-V6 lead ECG. For slight-degree myocardial ischemia ECG simulation, the ST segment and T wave changes can be observed from the simulated ECG based on a dynamic heart model, while the ST segment and T wave of simulated ECG based on a static heart model is almost unchanged when compared with a normal ECG. Conclusions: This study confirms the importance of the mechanical factor in the ECG simulation. The dynamic heart model could provide more accurate ECG simulation, especially for myocardial ischemia or infarction simulation, since the main ECG changes occur at the ST segment and T wave, which correspond with cardiac systole and diastole phases.

Field observations of instantaneous water slopes and horizontal pressure gradients in the swash-zone

Field observations of instantaneous water surface slopes in the swash zone are presented. For free-surface flows with a hydrostatic pressure distribution the surface slope is equivalent to the horizontal pressure gradient. Observations were made using a novel technique which in its simplest form consists of a horizontal stringline extending seaward from the beach face. Visual observation, still photography or video photography is then sufficient to determine the surface slope where the free-surface cuts the line or between reference points in the image. The method resolves the mean surface gradient over a cross-shore distance of 5 m or more to within +/- 0.001, or 1/20th -1/100th of typical beach gradients. In addition, at selected points and at any instant in time during the swash cycle, the water surface slope can be determined exactly to be dipping either seaward or landward. Close to the location of bore collapse landward dipping water surface slopes of order 0.05-0.1 occur over a very small region (order 0.5 m) at the blunt or convex leading edge of the swash. In the middle and upper swash the water surface slope at this leading edge is usually very close to horizontal or slightly seaward. Behind the leading edge, the water surface slope was observed to be very close to horizontal or dipping seaward at all times throughout the swash uprush. During the backwash the water surface slope was observed to be always dipping seaward, approaching the beach slope, and remained seaward until a new uprush edge or incident bore passed any particular cross-shore location of interest. The observations strongly Suggest that the swash boundary layer is subject to an adverse pressure gradient during uprush and a favourable pressure gradient during the backwash. Furthermore, assuming Euler's equations are a good approximation in the swash, the observations also show that the total fluid acceleration is negative (offshore) for almost the whole of the uprush and for the entire backwash. The observations are contrary to recent work suggesting significant shoreward directed accelerations and pressure gradients occur in the swash (i.e., delta u/delta t > 0 similar to delta p/delta x < 0), but consistent with analytical and numerical solutions for swash uprush and backwash. The results have important implications for sediment transport modelling in the swash zone.

Influence of adsorbate interaction on transport in confined spaces

This article provides a review of the recent theory of transport in nanopores developed in the author's laboratory. In particular the influence of fluid-solid interactions on the transport coefficient is examined, showing that such interactions reduce the value of the coefficient by almost an order of magnitude in comparison to the Knudsen theory for non-interacting systems. The activation energy and potential energy barriers for diffusion in smooth pores with a one-dimensional potential energy profile are also discussed, indicating the inadequacy of the commonly used assumption of proportionality between the activation energy and heat of adsorption or the minimum pore potential energy. A further feature affected by fluid-solid interactions is the nature of the reflection of fluid molecules colliding with a pore wall surface, varying from being nearly specular - such as in carbon nanotubes - to nearly diffuse for amorphous solids. Diffuse reflection leads to momentum loss and reduced transport coefficients. However, fluid-solid interactions do not affect the transport coefficient in the single-file diffusion regime when the surface reflection is diffuse, and the transport coefficient in this case is largely independent of the adsorbed density.

Kinetic study of the thermal degradation of high density polyethylene

The thermal degradation of high density polyethylene has been modelled by the random breakage of polymer bonds, using a set of population balance equations. A model was proposed in which the population balances were lumped into representative sizes so that the experimentally determined molecular weight distribution of the original polymer could be used as the initial condition. This model was then compared to two different cases of the unlumped population balance which assumed unimolecular initial distributions of 100 and 500 monomer units, respectively. The model that utilised the experimentally determined molecular weight distribution was found to best describe the experimental data. The model fits suggested a second mechanism in addition to random breakage at slow reaction rates. (c) 2005 Elsevier Ltd. All rights reserved.

London-van der Waals adhesiveness of rough particles

Van der Waals forces often dominate interactions and adhesion between fine particles and, in turn, decisively influence the bulk behaviour of powders. However, so far there is no effective means to characterize the adhesive behaviour of such particles. A complication is that most powder particles have rough surfaces, and it is the asperities on the surfaces that touch, confounding the actual surface that is in contact. Conventional approaches using surface energy provide limited information regarding adhesion, and pull-off forces measured through atomic force microscope (AFM) are highly variable and difficult to interpret. In this paper we develop a model which combines the Rumpf-Rabinovich and the JKR-DMT theories to account simultaneously for the effects of surface roughness and deformation on adhesion. This is applied to a 'characteristic asperity' which may be easily obtained from AFM measurements. The concept of adhesiveness, a material property reflecting the influences of elastic deformability, surface roughness, and interfacial surface energy, is introduced as an efficient and quantitative measure of the adhering tendency of a powder. Furthermore, a novel concept of specific adhesiveness is proposed as a convenient tool for characterizing and benchmarking solid materials. This paper provides an example to illustrate the use of the proposed theories. (c) 2005 Elsevier B.V. All rights reserved.

Long wave generation by the shoaling and breaking of transient wave groups on a beach

This paper presents new laboratory data on the generation of long waves by the shoaling and breaking of transient-focused short-wave groups. Direct offshore radiation of long waves from the breakpoint is shown experimentally for the first time. High spatial resolution enables identification of the relationship between the spatial gradients of the short-wave envelope and the long-wave surface. This relationship is consistent with radiation stress theory even well inside the surf zone and appears as a result of the strong nonlinear forcing associated with the transient group. In shallow water, the change in depth across the group leads to asymmetry in the forcing which generates significant dynamic setup in front of the group during shoaling. Strong amplification of the incident dynamic setup occurs after short-wave breaking. The data show the radiation of a transient long wave dominated by a pulse of positive elevation, preceded and followed by weaker trailing waves with negative elevation. The instantaneous cross-shore structure of the long wave shows the mechanics of the reflection process and the formation of a transient node in the inner surf zone. The wave run-up and relative amplitude of the radiated and incident long waves suggests significant modification of the incident bound wave in the inner surf zone and, the dominance of long waves generated by the breaking process. It is proposed that these conditions occur when the primary short waves and bound wave are not shallow water waves at the breakpoint. A simple criterion is given to determine these conditions, which generally occur for the important case of storm waves.

Multidimensional modelling to investigate interspecies hydrogen transfer in anaerobic biofilms

Anaerobic digestion is a multistep process, mediated by a functionally and phylogenetically diverse microbial population. One of the crucial steps is oxidation of organic acids, with electron transfer via hydrogen or formate from acetogenic bacteria to methanogens. This syntrophic microbiological process is strongly restricted by a thermodynamic limitation on the allowable hydrogen or formate concentration. In order to study this process in more detail, we developed an individual-based biofilm model which enables to describe the processes at a microbial resolution. The biochemical model is the ADM1, implemented in a multidimensional domain. With this model, we evaluated three important issues for the syntrophic relationship: (i) is there a fundamental difference in using hydrogen or formate as electron carrier? (ii) Does a thermodynamic-based inhibition function produced substantially different results from an empirical function? and; (iii) Does the physical colocation of acetogens and methanogens follow directly from a general model. Hydrogen or formate as electron carrier had no substantial impact on model results. Standard inhibition functions or thermodynamic inhibition function gave similar results at larger substrate field grid sizes (> 10 mu m), but at smaller grid sizes, the thermodynamic-based function reduced the number of cells with long interspecies distances (> 2.5 mu m). Therefore, a very fine grid resolution is needed to reflect differences between the thermodynamic function, and a more generic inhibition form. The co-location of syntrophic bacteria was well predicted without a need to assume a microbiological based mechanism (e.g., through chemotaxis) of biofilm formation.