Finite element modelling of reactive mass transport problems in fluid-saturated porous media

We use the finite element method to solve reactive mass transport problems in fluid-saturated porous media. In particular, we discuss the mathematical expression of the chemical reaction terms involved in the mass transport equations for an isothermal, non-equilibrium chemical reaction. It has turned out that the Arrhenius law in chemistry is a good mathematical expression for such non-equilibrium chemical reactions especially from the computational point of view. Using the finite element method and the Arrhenius law, we investigate the distributions of PH (i.e. the concentration of H+) and the relevant reactive species in a groundwater system. Although the main focus of this study is on the contaminant transport problems in groundwater systems, the related numerical techniques and principles are equally applicable to the orebody formation problems in the geosciences. Copyright (C) 1999 John Wiley & Sons, Ltd.

Approach to designing rotating drum bioreactors for solid-state fermentation on the basis of dimensionless design factors

The development of large-scale solid-stale fermentation (SSF) processes is hampered by the lack of simple tools for the design of SSF bioreactors. The use of semifundamental mathematical models to design and operate SSF bioreactors can be complex. In this work, dimensionless design factors are used to predict the effects of scale and of operational variables on the performance of rotating drum bioreactors. The dimensionless design factor (DDF) is a ratio of the rate of heat generation to the rate of heat removal at the time of peak heat production. It can be used to predict maximum temperatures reached within the substrate bed for given operational variables. Alternatively, given the maximum temperature that can be tolerated during the fermentation, it can be used to explore the combinations of operating variables that prevent that temperature from being exceeded. Comparison of the predictions of the DDF approach with literature data for operation of rotating drums suggests that the DDF is a useful tool. The DDF approach was used to explore the consequences of three scale-up strategies on the required air flow rates and maximum temperatures achieved in the substrate bed as the bioreactor size was increased on the basis of geometric similarity. The first of these strategies was to maintain the superficial flow rate of the process air through the drum constant. The second was to maintain the ratio of volumes of air per volume of bioreactor constant. The third strategy was to adjust the air flow rate with increase in scale in such a manner as to maintain constant the maximum temperature attained in the substrate bed during the fermentation. (C) 2000 John Wiley & Sons, Inc.

Finite element modelling of dissipative structures for nonequilibrium chemical reactions in fluid-saturated porous media

We use the finite element method to model and predict the dissipative structures of chemical species for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. In particular, we explore the conditions under which dissipative structures of the species may exist in the Brusselator type of nonequilibrium chemical reaction. Since this is the first time the finite element method and related strategies have been used to study the chemical instability problems in a fluid-saturated porous medium, it is essential to validate the method and strategies before they are put into application. For this purpose, we have rigorously derived the analytical solutions for dissipative structures of chemical species in a benchmark problem, which geometrically is a square. Comparison of the numerical solutions with the analytical ones demonstrates that the proposed numerical method and strategy are robust enough to solve chemical instability problems in a fluid-saturated porous medium. Finally, the related numerical results from two application examples indicate that both the regime and the magnitude of pore-fluid flow have significant effects on the nature of the dissipative structures that developed for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. The motivation for this study is that self-organization under conditions of pore-fluid flow in a porous medium is a potential mechanism of the orebody formation and mineralization in the upper crust of the Earth. (C) 2000 Elsevier Science S.A. All rights reserved.

Finite element analysis of heat transfer and mineralization in layered hydrothermal systems with upward throughflow

We use the finite element method to solve the coupled problem between convective pore-fluid flow, heat transfer and mineralization in layered hydrothermal systems with upward throughflow. In particular, we present the improved rock alteration index (IRAI) concept for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in the systems. To validate the numerical method used in the computation, analytical solutions to a benchmark problem have been derived. After the numerical method is validated, it is used to investigate the pattern of pore-fluid Aom, the distribution of temperature and the mineralization pattern of gold minerals in a layered hydrothermal system with upward throughflow. The related numerical results have demonstrated that the present concept of IRAI is useful and applicable for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in hydrothermal systems. (C) 2000 Elsevier Science S.A. All rights reserved.

Influence of Mg content on the microstructure and solid solution chemistry of Al-7%Si-Mg casting alloys during solution treatment

The solution treatment stage of the T6 heat-treatment of Al-7%Si-Mg foundry alloys influences microstructural features such as Mg2Si dissolution, and eutectic silicon spheroidisation and coarsening. Microstructural and microanalytical studies have been conducted across a range of Sr-modified Al-7%Si alloys, with an Fe content of 0.12% and Mg contents ranging from 0.3-0.7wt%. Qualitative and quantitative metallography have shown that, in addition to the above changes, solution treatment also results in changes to the relative proportions of iron-containing intermetallic particles and that these changes are composition-dependent. While solution treatment causes a substantial transformation of pi phase to beta phase in low Mg alloys (0.3-0.4%), this change is not readily apparent at higher Mg levels (0.6-0.7%). The pi to beta transformation is accompanied by a release of Mg into the aluminum matrix over and above that which arises from the rapid dissolution of Mg2Si. Since the level of matrix Mg retained after quenching controls an alloy's subsequent precipitation hardening response, a proper understanding of this phase transformation is crucial if tensile properties are to be maximised.

Mesoionic pyridopyrimidinylium and pyridooxazinylium olates and non-mesoionic pyridopyrimidinones. Structures in the solid state, solution, and matrices

X-Ray crystal structures, C-13 NMR spectra and theoretical calculations (B3LYP/6-31G*) are reported for the mesoionic (zwitterionic) pyridopyrimidinylium- and pyridooxazinyliumolates 2a, 3a and 5a,b as well as the enol ether 11b and the enamine 11c. The 1-NH compounds like 1a, 2a and 3a exist in the mesoionic form in the crystal and in solution, but the OH tautomers such as 1b and 2b dominate in the gas phase as revealed by the Ar matrix IR spectra in conjunction with DFT calculations. All data indicate that the mesoionic compounds can be regarded as intramolecular pyridine-ketene zwitterions (cf. 16 --> 17) with a high degree of positive charge on the pyridinium nitrogen, a long pyridinium N-CO bond (ca. 1.44-1.49 Angstrom), and normal C=O double bonds (ca. 1.22 Angstrom). All mesoionic compounds exhibit a pronounced tilting of the olate C=O groups (the C=O groups formally derived from a ketene) towards the pyridinium nitrogen, giving NCO angles of 110-118 degrees. Calculations reveal a hydrogen bond with 6-CH, analogous to what is found in ketene-pyridine zwitterions and the C3O2-pyridine complex. The 2-OH tautomers of type 1b, 2b, and 11 also show a high degree of zwitterionic character as indicated by the canonical structures 11 12.

Finite element modelling of three-dimensional convection problems in fluid-saturated porous media heated from below

We use the finite element method to model three-dimensional convective pore-fluid flow in fluid-saturated porous media when they are heated from below. In particular, we employ the particle-tracking technique to mimic the trajectories of particles in three-dimensional fluid flow problems. The related numerical results demonstrated that: (1) The progressive asymptotic approach procedure, which was previously developed for the finite element modelling of two-dimensional convective pore-fluid flow problems, is equally applicable to the finite element modelling of three-dimensional convective pore-fluid flow in fluid-saturated porous media heated from below. (2) The perturbation of gravity at different planes has a significant effect on the pattern of three-dimensional convective pore-fluid flow and therefore, may influence the pattern of orebody formation and mineralization in three-dimensional hydrothermal systems. Copyright (C) 2001 John Wiley & Sons, Ltd.

Change in the rate of shell deposition and shell microstructure in response to shell borers in the abalone Haliotis rubra

Gastropod shells consist of two crystal types of calcium carbonate, an outer, prismatic calcite layer and an inner nacreous layer made of aragonite. In cross-section, the nacre of the nacreous layer appears to have a regular brick-like microstructure composed of thin laminae of aragonite crystals, separated by very thin sheets of protein (Lutz and Rhoads, 1980; Nakahara, 1983). In abalone (Genus, Haliotis) and other gastropods, thin layers of non-lamellar pigmented material occur within the nacre and have been termed alternatively, fine lines, growth rings or growth lines (Shepherd et al., 1995). It has been suggested that these pigmented layers are small, prismatic, calcite layers (Shepherd and Avalos-Borja, 1997; Zaremba et al., 1996) but investigations using a Raman laser in Haliotis rubra show that they contain aragonite rather than calcite (Hawkes et al, 1996). Day and Fleming (1992) suggest that the occurrence of pigmented layers is correlated with regular exogenous cues such as reproduction or temperature changes and indeed in some species, pigmented layers in the shell can be used to age abalone (review: Shepherd and Triantafillos, 1997). However, McShane and Smith (1992) suggest that pigmented layers can occur irregularly and therefore may be unreliable indicators of age.

Finite element modelling of heat transfer through permeable cracks in hydrothermal systems with upward throughflow

We use the finite element method to model the heat transfer phenomenon through permeable cracks in hydrothermal systems with upward throughflow. Since the finite element method is an approximate numerical method, the method must be validated before it is used to soh,e any new, kind of problem. However, the analytical solution, which can be used to validate the finite element method and other numerical methods, is rather limited in the literature, especially, for the problem considered here. Keeping this in mind, we have derived analytical solutions for the temperature distribution along the vertical axis of a crack in a fluid-saturated porous layer. After the finite element method is validated by comparing the numerical solution with the analytical solution for the same benchmark problem, it is used to investigate the pore-fluid flow and heat transfer in layered hydrothermal systems with vertical permeable cracks. The related analytical and numerical results have demonstrated that vertical cracks are effective and efficient members to transfer heat energy from the bottom section to the top section in hydrothermal systems with upward throughflow.

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.

Thermal analysis, nuclear magnetic resonance spectroscopy, and impedance spectroscopy of N,N-dimethyl-pyrrolidinium iodide: An ionic solid exhibiting rotator phases

N,N-dimethyl-pyrrolidinium iodide has been investigated using differential scanning calorimetry, nuclear magnetic resonance (NMR) spectroscopy, second moment calculations, and impedance spectroscopy. This pyrrolidinium salt exhibits two solid-solid phase transitions, one at 373 K having an entropy change, Delta S, of 38 J mol(-1) K-1 and one at 478 K having Delta S of 5.7 J mol(-1) K-1. The second moment calculations relate the lower temperature transition to a homogenization of the sample in terms of the mobility of the cations, while the high temperature phase transition is within the temperature region of isotropic tumbling of the cations. At higher temperatures a further decrease in the H-1 NMR linewidth is observed which is suggested to be due to diffusion of the cations. (C) 2005 American Institute of Physics.

Solution structure of χ-conopeptide MrIA, a modulator of the human norepinephrine transporter

The chi-conopeptides MrIA and MrIB are 13-residue peptides with two disulfide bonds that inhibit human and rat norepinephrine transporter systems and are of significant interest for the design of novel drugs involved in pain treatment. In the current study we have determined the solution structure of MrIA using NMR spectroscopy. The major element of secondary structure is a hairpin with the two strands connected by an inverse gamma-turn. The residues primarily involved in activity have previously been shown to be located in the turn region (Sharpe, I. A.; Palant, E.: Schroder, C. L; Kaye, D. M.; Adams, D. I.; Alewood, P. F.; Lewis, R. J. J Biol Client 2003, 278, 40317-40323), which appears to be more flexible than the beta-strands based on disorder in the ensemble of calculated structures. Analogues of MrIA with N-terminal truncations indicate that the N-terminal residues play a role in defining a stable conformation and the native disulfide connectivity. In particular, noncovalent interactions between Val3 and Hypl2 are likely to be involved in maintaining a stable conformation. The N-terminus also affects activity, as a single N-terminal deletion introduced additional pharmacology at rat vas deferens, while deleting the first two amino acids reduced chi-conopeptide potency. This article was originally published online as an accepted preprint. The Published Online date corresponds to the preprint version. You can request a copy of the preprint by entailing the Biopolymers editorial office at biopolymers@wiley.com (c) 2005 Wiley Periodicals, Inc.

Shear deformation at 29% solid during solidification of magnesium alloy AZ91 and aluminium alloy A356

Partially solid commercial Al-Si and Mg-Al alloys have been deformed in shear during solidification using vane rheometry. The dendritic mush was deformed for a short period at 29% solid and allowed to cool naturally after deformation. Both alloys exhibited yield point behaviour and deformation was highly localised at the surface of maximum shear stress. The short period of deformation was found to have a distinct impact on the as-cast microstructure leading to fragmented dendrites in the deformation region of both alloys. In the case of the Mg-Al alloy, a concentrated region of interdendritic porosity was also observed in the deformation region. Concentrated porosity was not observed in the Al-Si alloy. (c) 2005 Elsevier B.V. All rights reserved.

Molecular orbital calculations of two-electron states for P-donor solid-state spin qubits

We theoretically study the Hilbert space structure of two neighboring P-donor electrons in silicon-based quantum computer architectures. To use electron spins as qubits, a crucial condition is the isolation of the electron spins from their environment, including the electronic orbital degrees of freedom. We provide detailed electronic structure calculations of both the single donor electron wave function and the two-electron pair wave function. We adopted a molecular orbital method for the two-electron problem, forming a basis with the calculated single donor electron orbitals. Our two-electron basis contains many singlet and triplet orbital excited states, in addition to the two simple ground state singlet and triplet orbitals usually used in the Heitler-London approximation to describe the two-electron donor pair wave function. We determined the excitation spectrum of the two-donor system, and study its dependence on strain, lattice position, and interdonor separation. This allows us to determine how isolated the ground state singlet and triplet orbitals are from the rest of the excited state Hilbert space. In addition to calculating the energy spectrum, we are also able to evaluate the exchange coupling between the two donor electrons, and the double occupancy probability that both electrons will reside on the same P donor. These two quantities are very important for logical operations in solid-state quantum computing devices, as a large exchange coupling achieves faster gating times, while the magnitude of the double occupancy probability can affect the error rate.

Pectinase production by solid fermentation from Aspergillus niger by a new prescription experiment

The Ordos Plateau in China is covered with up to 300,000 ha of peashrub (Caragana) which is the dominant natural vegetation and ideal for fodder production. To exploit peashrub fodder, it is crucially important to optimize the culture conditions, especially culture substrate to produce pectinase complex. In this study, a new prescription process was developed. The process, based on a uniform experimental design, first optimizes the solid substrate and second, after incubation, applies two different temperature treatments (30 degrees C for the first 30 h and 23 degrees C for the second 42 h) in the fermentation process. A multivariate regression analysis is applied to a number of independent variables (water, wheat bran, rice dextrose, ammonium sulfate, and Tween 80) to develop a predictive model of pectinase activity. A second-degree polynomial model is developed which accounts for an excellent proportion of the explained variation (R-2 = 97.7%). Using unconstrained mathematical programming, an optimized substrate prescription for pectinase production is subsequently developed. The mathematical analysis revealed that the optimal formula for pectinase production from Aspergillus niger by solid fermentation under the conditions of natural aeration, natural substrate pH (about 6.5), and environmental humidity of 60% is rice dextrose 8%, wheat bran 24%, ammonium sulfate ((NH4)(2)SO4) 6%, and water 61%. Tween 80 was found to have a negative effect on the production of pectinase in solid substrate. With this substrate prescription, pectinase produced by solid fermentation of A. niger reached 36.3IU/(gDM). Goats fed on the pectinase complex obtain an incremental increase of 0.47 kg day(-1) during the initial 25 days of feeding, which is a very promising new feeding prospect for the local peashrub. It is concluded that the new formula may be very useful for the sustainable development of and and semiarid pastures such as those of the Ordos Plateau. (c) 2005 Elsevier Inc. All rights reserved.