Finite element modelling of reactive fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins

We present the finite element simulations of reactive mineral carrying fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins. In particular we explore the mixing of reactive sulfide and sulfate fluids and the relevant patterns of mineralization for Load, zinc and iron minerals in the regime of temperature-gradient-driven convective flow. Since the mineralization and ore body formation may last quite a long period of time in a hydrothermal basin, it is commonly assumed that, in the geochemistry, the solutions of minerals are in an equilibrium state or near an equilibrium state. Therefore, the mineralization rate of a particular kind of mineral can be expressed as the product of the pore-fluid velocity and the equilibrium concentration of this particular kind of mineral Using the present mineralization rate of a mineral, the potential of the modern mineralization theory is illustrated by means of finite element studies related to reactive mineral-carrying fluids mixing problems in materially homogeneous and inhomogeneous porous rock basins.

Prediction of carbon monoxide in fires by conditional moment closure

Carbon monoxide is the chief killer in fires. Dangerous levels of CO can occur when reacting combustion gases are quenched by heat transfer, or by mixing of the fire plume in a cooled under- or overventilated upper layer. In this paper, carbon monoxide predictions for enclosure fires are modeled by the conditional moment closure (CMC) method and are compared with laboratory data. The modeled fire situation is a buoyant, turbulent, diffusion flame burning under a hood. The fire plume entrains fresh air, and the postflame gases are cooled considerably under the hood by conduction and radiation, emulating conditions which occur in enclosure fires and lead to the freezing of CO burnout. Predictions of CO in the cooled layer are presented in the context of a complete computational fluid dynamics solution of velocity, temperature, and major species concentrations. A range of underhood equivalence ratios, from rich to lean, are investigated. The CMC method predicts CO in very good agreement with data. In particular, CMC is able to correctly predict CO concentrations in lean cooled gases, showing its capability in conditions where reaction rates change considerably.

Theoretical and numerical analyses of convective instability in porous media with temperature-dependent viscosity

Exact analytical solutions of the critical Rayleigh numbers have been obtained for a hydrothermal system consisting of a horizontal porous layer with temperature-dependent viscosity. The boundary conditions considered are constant temperature and zero vertical Darcy velocity at both the top and bottom of the layer. Not only can the derived analytical solutions be readily used to examine the effect of the temperature-dependent viscosity on the temperature-gradient driven convective flow, but also they can be used to validate the numerical methods such as the finite-element method and finite-difference method for dealing with the same kind of problem. The related analytical and numerical results demonstrated that the temperature-dependent viscosity destabilizes the temperature-gradient driven convective flow and therefore, may affect the ore body formation and mineralization in the upper crust of the Earth. Copyright (C) 2003 John Wiley Sons, Ltd.

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.

Numerical modelling of chemical effects of magma solidification problems in porous rocks

The solidification of intruded magma in porous rocks can result in the following two consequences: (1) the heat release due to the solidification of the interface between the rock and intruded magma and (2) the mass release of the volatile fluids in the region where the intruded magma is solidified into the rock. Traditionally, the intruded magma solidification problem is treated as a moving interface (i.e. the solidification interface between the rock and intruded magma) problem to consider these consequences in conventional numerical methods. This paper presents an alternative new approach to simulate thermal and chemical consequences/effects of magma intrusion in geological systems, which are composed of porous rocks. In the proposed new approach and algorithm, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with the proposed mass source and physically equivalent heat source. The major advantage in using the proposed equivalent algorithm is that a fixed mesh of finite elements with a variable integration time-step can be employed to simulate the consequences and effects of the intruded magma solidification using the conventional finite element method. The correctness and usefulness of the proposed equivalent algorithm have been demonstrated by a benchmark magma solidification problem. Copyright (c) 2005 John Wiley & Sons, Ltd.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

Simulation of mixing in unsteady flow through a periodically obstructed channel

Fluid mixing in steady and unsteady Bow through a channel containing periodic square obstructions has been studied using a finite-difference simulation to determine fluid velocities, followed by the use of passive marker particle advection to look at fluid transport out of the cavities formed between each of the obstructions. The geometry and Bow conditions were chosen from the work by Perkins (1989, M.S. Thesis, Lehigh University; 1992, Ph.D. Thesis, Lehigh University); who investigated heat transfer enhancement due to unsteady flow through such an obstructed channel. Particle advection shows that Bow regimes which are predicted to give good mixing based on snapshots of instantaneous streamline contour plots were not necessarily able to efficiently mix fluid which started in the cavity regions throughout the channel. The use of Poincare sections shows regular regions existing under these conditions which inhibit efficient fluid transport. These regular regions are found to disappear when the unsteady Bow velocity is increased. (C) 1997 Elsevier Science Ltd.

Encephalitogenicity of murine, but not bovine, DM20 in SJL mice is due to a single amino acid difference in the immunodominant encephalitogenic epitope

Myelin proteolipid protein (PLP) contains 2 immunodominant encephalitogenic epitopes in SJL mice, namely PLP residues 139-151 and 178-191. DM20, a minor isoform of PLP, lacks residues 116-150 and consequently contains only the single major encephalitogenic epitope 178-191. However, it has been found previously that bovine DM20 is not encephalitogenic in SJL mice. Since residue 188 within peptide 178-191 is phenylalanine (F) in murine DM20 and alanine (A) in bovine DM20, we tested the effect of this difference on the immune responses and induction of EAE. SJL mice were immunized with either highly purified murine or bovine DM20. Residues 178-191 were found to be immunodominant for each, but only murine and not bovine DM20 was encephalitogenic. A synthetic peptide corresponding to the murine 178-191 sequence (F188) was also encephalitogenic, whereas the peptide corresponding to the bovine sequence (A188) was not. Both F188 and A188 bind with high affinity to I-A(s) and both are recognized by the SJL T cell repertoire. A188-specific T cell lines reacted to both A188 and F188, but F188-specific T cell lines were not stimulated by A188. F188-specific T cell lines produced mRNA for the Th1 cytokines IL2 and IFN gamma and, in passive transfer experiments, were encephalitogenic upon stimulation with F188, but not A188. In contrast, A188-specific T cell lines produced mRNA for IL4, IL5 and IL10, in addition to IL2 and IFN gamma, and were not encephalitogenic after stimulation with either F188 or A188. Cotransfer of A188-specific T cell lines with F188-specific T cell lines resulted in protection from EAE. Thus, A188 induces a functionally different phenotype of T cells from that induced by F188. Taken together these data suggest that the failure of bovine DM20 to induce EAE may be attributable to induction of protective rather than pathogenic T cells by the immunodominant epitope.

A compact shock-assisted free-piston driver for impulse facilities

A free-piston driver that employs entropy-raising shock processes with diaphragm rupture has been constructed, which promises significant theoretical advantages over isentropic compression. Results from a range of conditions with helium and argon driver gases are reported. Significant performance gains were achieved in some test cases. Heat losses are shown to have a strong effect on driver processes. Measurements compare well with predictions from a quasi-one-dimensional numerical code.

The application of electron microscopy and microanalysis in conjunction with organic petrology to further the understanding of organic-mineral association: examples from Mount Isa and McArthur Basins, Australia

Organic petrology supported by electron microscopical and micro-analytical techniques was applied to organic matter in Proterozoic sediments to better understand hydrothermal processes responsible for ore-grade mineralisation. It was shown that organic maturation was not only closely linked to the geological history of the sediments, but also highlighted heat transfer by convection as differentiated from conduction solely through sediment burial and step-wise subsidence. Water-rock ratios effect organic maturation in hydrothermal systems, and erratic reflectance profiles are indicators of convective heat transfer. Identification and characterisation of organic materials in terms of source rock and migrated hydrocarbons was shown to be a powerful tool in reconstructing the thermal history of sediments, identifying hydrothermal episodes, fluid pathways and heat source in the northern Australian Proterozoic basins. Higher reflectance of organic matter towards the central parts of the Mount Isa Basin and some of the most northerly parts point to proximity to higher heat flow at times, in contrast to relatively low temperatures prevailing in the western parts of the basin, next to the Murphy Inlier. A close correlation shown between peak organic reflectance values and super-sequence boundaries farther highlighted the valuable information to be gained from organic petrology, by allowing the separation of processes responsible for metal dissolution and transportation from those inducing precipitation. (C) 2001 Elsevier Science B.V All rights reserved.

A generalised dynamic model for char particle gasification with structure evolution and peripheral fragmentation

A generalised model for the prediction of single char particle gasification dynamics, accounting for multi-component mass transfer with chemical reaction, heat transfer, as well as structure evolution and peripheral fragmentation is developed in this paper. Maxwell-Stefan analysis is uniquely applied to both micro and macropores within the framework of the dusty-gas model to account for the bidisperse nature of the char, which differs significantly from the conventional models that are based on a single pore type. The peripheral fragmentation and random-pore correlation incorporated into the model enable prediction of structure/reactivity relationships. The occurrence of chemical reaction within the boundary layer reported by Biggs and Agarwal (Chem. Eng. Sci. 52 (1997) 941) has been confirmed through an analysis of CO/CO2 product ratio obtained from model simulations. However, it is also quantitatively observed that the significance of boundary layer reaction reduces notably with the reduction of oxygen concentration in the flue gas, operational pressure and film thickness. Computations have also shown that in the presence of diffusional gradients peripheral fragmentation occurs in the early stages on the surface, after which conversion quickens significantly due to small particle size. Results of the early commencement of peripheral fragmentation at relatively low overall conversion obtained from a large number of simulations agree well with experimental observations reported by Feng and Bhatia (Energy & Fuels 14 (2000) 297). Comprehensive analysis of simulation results is carried out based on well accepted physical principles to rationalise model prediction. (C) 2001 Elsevier Science Ltd. AH rights reserved.

A comparison of thermal condition between pilot- and full-scale furnaces for studying slagging and fouling propensity in PF boilers

The Australian Coal Industry Research Laboratory (ACIRL) furnace is scaled to simulate slagging and fouling in operating boilers. This requires that the gas and target temperatures, the heat flux, and the flow pattern be the same as those in real boilers. The gas and target temperatures are maintained by insulating the wall and cooling the target respectively. The flow pattern of a small burner cannot be the same as a large furnace. However, this flow pattern is partially compensated for by placing the slagging panels in three vertical locations. The paper develops the models of radiant heat transfer from the flame to the deposits both in pilot-scale and full-scale furnaces. They are used to compare the effective radiant heat transfer of the pilot- and full-scale furnaces. The experimental data both from the pilot- and full-scale furnaces are used to verify the incident heat flux and temperature profiles in the pilot- and full-scale furnaces. The results showed that the thermal condition in the pilot-scale furnace meets the requirements for studying the slagging regarding the gas temperature and the incident heat flux, particularly for the panel #1. The gas temperature in the convective section also meets the requirement for studying the fouling.

Limitations to photosynthesis of lettuce grown under tropical conditions: alleviation by root-zone cooling

Aerial parts of lettuce plants were grown under natural tropical fluctuating ambient temperatures, but with their roots exposed to two different root-rone temperatures (RZTs): a constant 20 degreesC-RZT and a fluctuating ambient (A-) RZT from 23-40 degreesC, Plants grown at A-RZT showed lower photosynthetic CO2 assimilation (A), stomatal conductance (g(s)), midday leaf relative water content (RWC), and chlorophyll fluorescence ratio F-v/F-m than 20 degreesC-RZT plants on both sunny and cloudy days. Substantial midday depression of A and g(s) occurred on both sunny and cloudy days in both RZT treatments, although F-v/F-m did not vary diurnally on cloudy days. Reciprocal temperature transfer experiments investigated the occurrence and possible causes of stomatal and non-stomatal limitations of photosynthesis. For both temperature transfers, light-saturated stomatal conductance (g(s) (sat)) and photosynthetic CO2 assimilation (A(sat)) were highly correlated with each other and with midday RWC, suggesting that A was limited by water stress-mediated stomatal closure, However, prolonged growth at A-RZT reduced light- and CO2-saturated photosynthetic O-2 evolution (P-max), indicating non-stomatal limitation of photosynthesis. Tight temporal coupling of leaf nitrogen content and P-max during both temperature transfers suggested that decreased nutrient status caused this non-stomatal limitation of photosynthesis.

Transducer for direct measurement of skin friction in hypervelocity impulse facilities

An acceleration compensated transducer was developed to enable the direct measurement of skin friction in hypervelocity impulse facilities. The gauge incorporated a measurement and acceleration element that employed direct shear of a piezoelectric ceramic. The design integrated techniques to maximize rise time and shear response while minimizing the affects of acceleration, pressure, heat transfer, and electrical interference. The arrangement resulted in a transducer natural frequency near 40 kHz. The transducer was calibrated for shear and acceleration in separate bench tests and was calibrated for pressure within an impulse facility. Uncertainty analyses identified only small experimental errors in the shear and acceleration calibration techniques. Although significant errors were revealed in the method of pressure calibration, total skin-friction measurement errors as low as +/-7-12% were established. The transducer was successfully utilized in a shock tunnel, and sample measurements are presented for flow conditions that simulate a flight Mach number near 8.

Numerical investigations of flow and energy fields near a thermoacoustic couple

The flow field and the energy transport near thermoacoustic couples are simulated using a 2D full Navier-Stokes solver. The thermoacoustic couple plate is maintained at a constant temperature; plate lengths, which are short and long compared with the particle displacement lengths of the acoustic standing waves, are tested. Also investigated are the effects of plate spacing and the amplitude of the standing wave. Results are examined in the form of energy vectors, particle paths, and overall entropy generation rates. These show that a net heat-pumping effect appears only near the edges of thermoacoustic couple plates, within about a particle displacement distance from the ends. A heat-pumping effect can be seen even on the shortest plates tested when the plate spacing exceeds the thermal penetration depth. It is observed that energy dissipation near the plate increases quadratically as the plate spacing is reduced. The results also indicate that there may be a larger scale vortical motion outside the plates which disappears as the plate spacing is reduced. (C) 2002 Acoustical Society of America.