Theory manual to OctVCE – a cartesian cell CFD code with special application to blast wave problems, Department of Mechanical Engineering Report 2007/12

OctVCE is a cartesian cell CFD code produced especially for numerical simulations of shock and blast wave interactions with complex geometries, in particular, from explosions. Virtual Cell Embedding (VCE) was chosen as its cartesian cell kernel for its simplicity and sufficiency for practical engineering design problems. The code uses a finite-volume formulation of the unsteady Euler equations with a second order explicit Runge-Kutta Godonov (MUSCL) scheme. Gradients are calculated using a least-squares method with a minmod limiter. Flux solvers used are AUSM, AUSMDV and EFM. No fluid-structure coupling or chemical reactions are allowed, but gas models can be perfect gas and JWL or JWLB for the explosive products. This report also describes the code’s ‘octree’ mesh adaptive capability and point-inclusion query procedures for the VCE geometry engine. Finally, some space will also be devoted to describing code parallelization using the shared-memory OpenMP paradigm. The user manual to the code is to be found in the companion report 2007/13.

User guide for shock and blast simulation with the OctVCE code (version 3.5+), Department of Mechanical Engineering Report 2007/13

OctVCE is a cartesian cell CFD code produced especially for numerical simulations of shock and blast wave interactions with complex geometries. Virtual Cell Embedding (VCE) was chosen as its cartesian cell kernel as it is simple to code and sufficient for practical engineering design problems. This also makes the code much more ‘user-friendly’ than structured grid approaches as the gridding process is done automatically. The CFD methodology relies on a finite-volume formulation of the unsteady Euler equations and is solved using a standard explicit Godonov (MUSCL) scheme. Both octree-based adaptive mesh refinement and shared-memory parallel processing capability have also been incorporated. For further details on the theory behind the code, see the companion report 2007/12.

Electrical impedance particle size method (Coulter Counter) detects the large fat globules in poorly homogenised UHT processed milk

The large fat globules that can be present in UHT milk due to inadequate homogenisation cause a cream layer to form that limits the shelf life of UHT milk. Four different particle size measurement techniques were used to measure the size of fat globules in poorly homogenised UHT milk processed in a UHT pilot plant. The thickness of the cream layer that formed during storage was negatively correlated with homogenisation pressure. It was positively correlated with the mass mean diameter and the percentage volume of particles between 1.5 and 2 mu m diameter, as determined by laser light scattering using the Malvern Mastersizer. Also, the thickness of the cream layer was positively correlated with the volume mode diameter and the percentage volume of particles between 1.5 and 2 mu m diameter, as determined by electrical impedance using the Coulter Counter. The cream layer thickness did not correlate significantly with the Coulter Counter measurements of volume mean diameter, or volume percentages of particles between 2 and 5 mu m or 5 and 10 mu m diameter. Spectroturbidimetry (Emulsion Quality Analyser) and light microscopy analyses were found to be unsuitable for assessing the size of the fat particles. This study suggests that the fat globule size distribution as determined by the electrical impedance method (Coulter Counter) is the most useful for determining the efficiency of homogenisation and therefore for predicting the stability of the fat emulsion in UHT milk during storage.

Data analysis in multiple-frequency bioelectrical impedance analysis

The performance of three analytical methods for multiple-frequency bioelectrical impedance analysis (MFBIA) data was assessed. The methods were the established method of Cole and Cole, the newly proposed method of Siconolfi and co-workers and a modification of this procedure. Method performance was assessed from the adequacy of the curve fitting techniques, as judged by the correlation coefficient and standard error of the estimate, and the accuracy of the different methods in determining the theoretical values of impedance parameters describing a set of model electrical circuits. The experimental data were well fitted by all curve-fitting procedures (r = 0.9 with SEE 0.3 to 3.5% or better for most circuit-procedure combinations). Cole-Cole modelling provided the most accurate estimates of circuit impedance values, generally within 1-2% of the theoretical values, followed by the Siconolfi procedure using a sixth-order polynomial regression (1-6% variation). None of the methods, however, accurately estimated circuit parameters when the measured impedances were low (<20 Omega) reflecting the electronic limits of the impedance meter used. These data suggest that Cole-Cole modelling remains the preferred method for the analysis of MFBIA data.

Modeling ex vivo hematopoiesis using chemical engineering metaphors

Ex vivo hematopoiesis is increasingly used for clinical applications. Models of ex vivo hematopoiesis are required to better understand the complex dynamics and to optimize hematopoietic culture processes. A general mathematical modeling framework is developed which uses traditional chemical engineering metaphors to describe the complex hematopoietic dynamics. Tanks and tubular reactors are used to describe the (pseudo-) stochastic and deterministic elements of hematopoiesis, respectively. Cells at any point in the differentiation process can belong to either an immobilized, inert phase (quiescent cells) or a mobile, active phase (cycling cells). The model describes five processes: (1) flow (differentiation), (2) autocatalytic formation (growth),(3) degradation (death), (4) phase transition from immobilized to mobile phase (quiescent to cycling transition), and (5) phase transition from mobile to immobilized phase (cycling to quiescent transition). The modeling framework is illustrated with an example concerning the effect of TGF-beta 1 on erythropoiesis. (C) 1998 Published by Elsevier Science Ltd. All rights reserved.

Calculating current densities and fields due to shielded bi-planar radio-frequency coils

A method for the accurate computation of the current densities produced in a wide-runged bi-planar radio-frequency coil is presented. The device has applications in magnetic resonance imaging. There is a set of opposing primary rungs, symmetrically placed on parallel planes and a similar arrangement of rungs on two parallel planes surrounding the primary serves as a shield. Current densities induced in these primary and shielding rungs are calculated to a high degree of accuracy using an integral-equation approach, combined with the inverse finite Hilbert transform. Once these densities are known, accurate electrical and magnetic fields are then computed without difficulty. Some test results are shown. The method is so rapid that it can be incorporated into optimization software. Some preliminary fields produced from optimized coils are presented.

Percolative fragmentation of char particles during gasification

Percolative fragmentation was confirmed to occur during gasification of three microporous coal chars. Indirect evidence obtained by the variation of electrical resistivity (ER) with conversion was supported by direct observation of numerous fragments during gasification. The resistivity increases slowly at low conversions and then sharply after a certain conversion value, which is a typical percolation phenomenon suggesting the occurrence of internal fragmentation at high conversion. Two percolation models are applied to interpret the experimental data and determine the percolation threshold. A percolation threshold of 0.02-0.07 was found, corresponding to a critical conversion of 92-96% for fragmentation. The electrical resistivity variation at high conversions is found to be very sensitive to diffusional effects during gasification. Partially burnt samples with a narrow initial particle size range were also observed microscopically, and found to yield a large number of small fragments even when the particles showed no disintegration and chemical control prevailed. It is proposed that this is due to the separation of isolated clusters from the particle surface. The particle size distribution of the fragments was essentially independent of the reaction conditions and the char type, and supported the prediction by percolation theory that the number fraction distribution varies linearly with mass in a log-log plot. The results imply that perimeter fragmentation would occur in practical combustion systems in which the reactions are strongly diffusion affected.

In vivo measurement of the effect of intra-abdominal pressure on the human spine

In humans, intra-abdominal pressure (IAP) is elevated during many everyday activities. This experiment aimed to investigate the extent to which increased IAP-without concurrent activity of the abdominal or back extensor muscles-produces an extensor torque. With subjects positioned in side lying on a swivel table with its axis at L3, moments about this vertebral level were measured when IAP was transiently increased by electrical stimulation of the diaphragm via the phrenic nerve. There was no electromyographic activity in abdominal and back extensor muscles. When IAP was increased artificially to similar to 15% of the maximum IAP amplitude that could be generated voluntarily with the trunk positioned in flexion, a trunk extensor moment (similar to6 Nm) was recorded. The size of the effect was proportional to the increase in pressure. The extensor moment was consistent with that predicted from a model based on measurements of abdominal cross-sectional area and IAP moment arm. When IAP was momentarily increased while the trunk was flexed passively at a constant velocity, the external torque required to maintain the velocity was increased. These results provide the first in vivo data of the amplitude of extensor moment that is produced by increased IAP. Although the net effect of this extensor torque in functional tasks would be dependent on the muscles used to increase the IAP and their associated flexion torque, the data do provide evidence that IAP contributes, at least in part, to spinal stability. (C) 2001 Elsevier Science Ltd. All rights reserved.

A time-harmonic inverse methodology for the design of RF coils in MRI

An inverse methodology is described to assist in the design of radio-frequency (RF) coils for magnetic resonance imaging (MRI) applications. The time-harmonic electromagnetic Green's functions are used to calculate current on the coil and shield cylinders that will generate a specified internal magnetic field. Stream function techniques and the method of moments are then used to implement this theoretical current density into an RF coil. A novel asymmetric coil operating for a 4.5 T MRI machine was designed and constructed using this methodology and the results are presented.

Using multidimensional projection techniques for reaching a high distinguishing ability in biosensing

Recent advances in the control of molecular engineering architectures have allowed unprecedented ability of molecular recognition in biosensing, with a promising impact for clinical diagnosis and environment control. The availability of large amounts of data from electrical, optical, or electrochemical measurements requires, however, sophisticated data treatment in order to optimize sensing performance. In this study, we show how an information visualization system based on projections, referred to as Projection Explorer (PEx), can be used to achieve high performance for biosensors made with nanostructured films containing immobilized antigens. As a proof of concept, various visualizations were obtained with impedance spectroscopy data from an array of sensors whose electrical response could be specific toward a given antibody (analyte) owing to molecular recognition processes. In addition to discussing the distinct methods for projection and normalization of the data, we demonstrate that an excellent distinction can be made between real samples tested positive for Chagas disease and Leishmaniasis, which could not be achieved with conventional statistical methods. Such high performance probably arose from the possibility of treating the data in the whole frequency range. Through a systematic analysis, it was inferred that Sammon`s mapping with standardization to normalize the data gives the best results, where distinction could be made of blood serum samples containing 10(-7) mg/mL of the antibody. The method inherent in PEx and the procedures for analyzing the impedance data are entirely generic and can be extended to optimize any type of sensor or biosensor.

Impedance characteristics of the diamond/carbon fiber electrodes for electrical double-layer capacitor

The electrochemical performance of carbon fibers (CF) and boron-doped diamond electrodes grown on carbon fiber substrate (BDD/CF) was studied. CF substrates were obtained from polyacrylonitrile precursor heat treated at two different temperatures of 1000 and 2000 degrees C to produce the desirable CF carbon graphitization index. This graphitization process influenced the CF conductivity and its chemical surface, also analyzed from X-ray photoelectron spectroscopy measurements. These three-dimensional CF structures allowed a high incorporation of diamond films compared to other carbon substrates such as glass carbon or HOPG. The electrochemical responses, from these four classes of electrodes, were evaluated focusing their application as electrical double-layer capacitors using cyclic voltammetry and impedance measurements. Cyclic voltammetry results revealed that the electrode formed from BDD grown on CF-2000 presented a typical capacitor behavior with the best rectangular shape, compared to those electrodes of CF or BDD/CF-1000. Furthermore, the BDD/CF-2000 electrode presented the lowest impedance, associated to its significant capacitance value of 1940 mu F/cm(2) taking into account the BDD films. This behavior was attributed to the strong dependence between diamond coating texture and the CF graphitization temperature. The largest surface area of BDD/CF-2000 was promoted by its singular film growth mechanism associated to the substrate chemical surface. (c) 2008 Elsevier B.V. All rights reserved.

Investigation of the electrical and electrochemical properties of nanocomposites from V2O5, polypyrrole, and polyaniline

In this work, an investigation of the electrical and electrochemical properties responsible for the energy storage capability of nanocomposites has been carried out. We demonstrate that, in the case of the V2O5 xerogel and the nanocomposites polypyrrole (Ppy)/V2O5 and polyaniline (PANI)/V2O5, the quadratic logistic equation (QLE) can be used to fit the inverse of the resistance values as a function of the injected charge in non-steady-state conditions. This contributes to a phenomenological understanding of the lithium ion and electron transport. The departure of the experimental curve from the fitting observed for the V2O5 xerogel can be attributed to the trapping sites formed during the lithium electroinsertion, which was observed by electrochemical impedance spectroscopy. The amount of trapping sites was obtained on the basis of the QLE. Similar values used to fit the inverse of the resistance were also used to fit the absorbance changes, which is also associated with the small polaron hopping from the V(IV) to the V(V) sites. On the other hand, there was good agreement between the experimental and the theoretical data when the profile of the inverse of the resistance as a function of the amount of inserted lithium ions of the nanocomposites Ppy/V2O5 and PANI/ V2O5 was concerned. We suggest that the presence of the conducting polymers is responsible for the different electrical profile of the V2O5 xerogel compared with those of the nanocomposites. In the latter case, interactions between the lithium ions and oxygen atoms from V2O5 are shielded, thus decreasing the trapping effect of lithium ions in the V2O5 sites. The different values of the lithium ion diffusion coefficient into these intercalation materials are in agreement with this hypothesis.

Modeling of columnar-to-equiaxed transition in solidified Al-Si alloys

A Cellular-Automaton Finite-Volume-Method (CAFVM) algorithm has been developed, coupling with macroscopic model for heat transfer calculation and microscopic models for nucleation and growth. The solution equations have been solved to determine the time-dependent constitutional undercooling and interface retardation during solidification. The constitutional undercooling is then coupled into the CAFVM algorithm to investigate both the effects of thermal and constitutional undercooling on columnar growth and crystal selection in the columnar zone, and formation of equiaxed crystals in the bulk liquid. The model cannot only simulate microstructures of alloys but also investigates nucleation mechanisms and growth kinetics of alloys solidified with various solute concentrations and solidification morphologies.

Signature of mott-insulator transition with ultracold fermions in a one-dimensional optical lattice

Using the exact Bethe ansatz solution of the Hubbard model and Luttinger liquid theory, we investigate the density profiles and collective modes of one-dimensional ultracold fermions confined in an optical lattice with a harmonic trapping potential. We determine a generic phase diagram in terms of a characteristic filling factor and a dimensionless coupling constant. The collective oscillations of the atomic mass density, a technique that is commonly used in experiments, provide a signature of the quantum phase transition from the metallic phase to the Mott-insulator phase. A detailed experimental implementation is proposed.