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.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Finite element modelling of rock alteration and metamorphic process in hydrothermal systems

We use the finite element method to simulate the rock alteration and metamorphic process in hydrothermal systems. In particular, we consider the fluid-rock interaction problems in pore-fluid saturated porous rocks. Since the fluid rock interaction takes place at the contact interface between the pore-fluid and solid minerals, it is governed by the chemical reaction which usually takes place very slowly at this contact interface, from the geochemical point of view. Due to the relative slowness of the rate of the chemical reaction to the velocity of the pore-fluid flow in the hydrothermal system to be considered, there exists a retardation zone, in which the conventional static theory in geochemistry does not hold true. Since this issue is often overlooked by some purely numerical modellers, it is emphasized in this paper. The related results from a typical rock alteration and metamorphic problem in a hydrothermal system have shown not only the detailed rock alteration and metamorphic process, but also the size of the retardation zone in the hydrothermal system. Copyright (C) 2001 John Wiley & Sons, Ltd.

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.

Scope for further analytical solutions for constant flux infiltration into a semi-infinite soil profile or redistribution in a finite soil profile

[1] We attempt to generate new solutions for the moisture content form of the one-dimensional Richards' [1931] equation using the Lisle [1992] equivalence mapping. This mapping is used as no more general set of transformations exists for mapping the one-dimensional Richards' equation into itself. Starting from a given solution, the mapping has the potential to generate an infinite number of new solutions for a series of nonlinear diffusivity and hydraulic conductivity functions. We first seek new analytical solutions satisfying Richards' equation subject to a constant flux surface boundary condition for a semi-infinite dry soil, starting with the Burgers model. The first iteration produces an existing solution, while subsequent iterations are shown to endlessly reproduce this same solution. Next, we briefly consider the problem of redistribution in a finite-length soil. In this case, Lisle's equivalence mapping is generalized to account for arbitrary initial conditions. As was the case for infiltration, however, it is found that new analytical solutions are not generated using the equivalence mapping, although existing solutions are recovered.

Finite-difference time-domain-based studies of MRI pulsed field gradient-induced eddy currents inside the human body

In modern magnetic resonance imaging (MRI), patients are exposed to strong, rapidly switching magnetic gradient fields that, in extreme cases, may be able to elicit nerve stimulation. This paper presents theoretical investigations into the spatial distribution of induced current inside human tissues caused by pulsed z-gradient fields. A variety of gradient waveforms have been studied. The simulations are based on a new, high-definition, finite-difference time-domain method and a realistic inhomogeneous 10-mm resolution human body model with appropriate tissue parameters. it was found that the eddy current densities are affected not only by the pulse sequences but by many parameters such as the position of the body inside the gradient set, the local biological material properties and the geometry of the body. The discussion contains a comparison of these results with previous results found in the literature. This study and the new methods presented herein will help to further investigate the biological effects caused by the switched gradient fields in a MRI scan. (C) 2002 Wiley Periodicals, Inc.

Potential triggers for akinete differentiation in an Australian strain of the cyanobacterium Cylindrospermopsis raciborskii (AWT 205/1)

Understanding the triggers for some cyanobacteria of the Nostocales and Stigonematales orders to produce specialised reproductive cells termed akinetes, is very important to gain further insights into their ecology. By improving our understanding of their life cycle, appropriate management options may be devised to control the formation of these cells, and therefore the potential bloom inoculum which they are thought to provide, may be reduced. This study investigated the effect of chemical (phosphorus limitation), and environmental variables (temperature shock) on akinete differentiation in the freshwater cyanobacterium Cylindrospermopsis raciborskii (AWT 205/1). From the preliminary results, it is suggested that the availability of phosphorus and changes in temperature were a necessary requirement for the formation of akinetes in this particular strain of C. raciborskii. In the four phosphorus treatments investigated (0, 3, 38 and 75 mug l(-1) P), only the two higher treatments produced akinetes (approximately 220 ml(-1)). When the first akinetes were observed in the 38 and 75 mug l(-1) P treatments, filterable reactive phosphorus (FRP) concentrations in the medium were approximately 22 and 52 mug l(-1) P, respectively, indicating that there was no phosphorus limitation. In the temperature shock experiment, akinetes were observed in the 15 and 20degreesC treatments. However, akinetes were degraded (pale yellow colour, limited swelling and shrivelled edges) and in much lower concentrations, which was thought to be a result of the daily temperature shock. We suggest that the formation of akinetes in C. raciborskii (AWT 205/1) can be triggered by an initial temperature shock and that phosphorus is a necessary requirement to allow further growth and full development of akinetes.

Estimating the potential refolding yield of recombinant proteins expressed as inclusion bodies

Recombinant protein production in bacteria is efficient except that insoluble inclusion bodies form when some gene sequences are expressed. Such proteins must undergo renaturation, which is an inefficient process due to protein aggregation on dilution from concentrated denaturant. In this study, the protein-protein interactions of eight distinct inclusion-body proteins are quantified, in different solution conditions, by measurement of protein second virial coefficients (SVCs). Protein solubility is shown to decrease as the SVC is reduced (i.e., as protein interactions become more attractive). Plots of SVC versus denaturant concentration demonstrate two clear groupings of proteins: a more aggregative group and a group having higher SVC and better solubility. A correlation of the measured SVC with protein molecular weight and hydropathicity, that is able to predict which group each of the eight proteins falls into, is presented. The inclusion of additives known to inhibit aggregation during renaturation improves solubility and increases the SVC of both protein groups. Furthermore, an estimate of maximum refolding yield (or solubility) using high-performance liquid chromatography was obtained for each protein tested, under different environmental conditions, enabling a relationship between yield and SVC to be demonstrated. Combined, the results enable an approximate estimation of the maximum refolding yield that is attainable for each of the eight proteins examined, under a selected chemical environment. Although the correlations must be tested with a far larger set of protein sequences, this work represents a significant move beyond empirical approaches for optimizing renaturation conditions. The approach moves toward the ideal of predicting maximum refolding yield using simple bioinformatic metrics that can be estimated from the gene sequence. Such a capability could potentially screen, in silico, those sequences suitable for expression in bacteria from those that must be expressed in more complex hosts. (C) 2004 Wiley Periodicals, Inc.

Production of fibrolytic enzymes by Aspergillus japonicus C03 using agro-industrial residues with potential application as additives in animal feed

Solid-state fermentation obtained from different and low-cost carbon sources was evaluated to endocellulases and endoxylanases production by Aspergillus japonicus C03. Regarding the enzymatic production the highest levels were observed at 30 degrees C, using soy bran added to crushed corncob or wheat bran added to sugarcane bagasse, humidified with salt solutions, and incubated for 3 days (xylanase) or 6 days (cellulase) with 70% relative humidity. Peptone improved the xylanase and cellulase activities in 12 and 29%, respectively. The optimum temperature corresponded to 60 degrees C and 50-55 degrees C for xylanase and cellulase, respectively, both having 4.0 as optimum pH. Xylanase was fully stable up to 40 degrees C, which is close to the rumen temperature. The enzymes were stable in pH 4.0-7.0. Cu(++) and Mn(++) increased xylanase and cellulase activities by 10 and 64%, respectively. A. japonicus C03 xylanase was greatly stable in goat rumen fluid for 4 h during in vivo and in vitro experiments.

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.

The therapeutic potential of recombinant BCG expressing the antigen S1PT in the intravesical treatment of bladder cancer

Purpose: Bacillus Calmette-Guerin (BCG) continues to be employed as the most effective immunotherapy against superficial bladder cancer. We have developed an rBCG-S1PT strain that induces a stronger cellular immune response than BCG. This preclinical study was designed to test the potential of rBCG-S1PT as an immunotherapeutic agent for intravesical bladder cancer therapy. Materials and methods: A tumor was induced in C57BL/6 mice after chemical cauterization of the bladder and inoculation of the tumor cell line MB49. Next, mice were treated by intravesical instillation with BCG, rBCG-S1PT, or PBS once a week for 4 weeks. After 35 days, the bladders were removed and weighed, Th1 (IL-2, IL-12, INOS, INF-gamma, TNF-alpha), and Th2 (IL-5, IL-6, IL-10, TGF-beta) cytokine mRNA responses in individual mice bladders were measured by quantitative real time PCR, and the viability of MB49 cells in 18-hour coculture with splenocytes from treated mice was assessed. In an equivalent experiment, animals were observed for 60 days to quantify their survival. Results: Both BCG and rBCG-S1PT immunotherapy resulted in bladder weight reduction, and rBCG-S1PT increased survival time compared with the control group. There were increases in TNF-alpha in the BCG treated group, as well as increases in TNF-alpha and IL-10 mRNA in the rBCG-S1PT group. The viability of MB49 cells cocultured with splenocytes from rBCG-S1PT-treated mice was lower than in both the BCG and control groups. Conclusions: rBCG-S1PT therapy improved outcomes and lengthened survival times. These results indicate that rBCG could serve as a useful substitute for wild-type BCG. (C) 2010 Elsevier Inc. All rights reserved.

Quantum dynamics of an atomic Bose-Einstein condensate in a double-well potential

We consider the quantum dynamics of a neutral atom Bose-Einstein condensate in a double-well potential, including many-body hard-sphere interactions. Using a mean-field factorization we show that the coherent oscillations due to tunneling are suppressed when the number of atoms exceeds a critical value. An exact quantum solution, in a two-mode approximation, shows that the mean-field solution is modulated by a quantum collapse and revival sequence.

Potential GABAB receptor antagonists .9. The synthesis of 3-amino-3-(4-chlorophenyl)propanoic acid, 2-amino-2-(4-chlorophenyl)ethylphosphonic acid and 2-amino-2-(4-chlorophenyl)ethanesulfonic acid

This paper describes the synthesis of 3-amino-3-(4-chlorophenyl)propanoic acid and the corresponding phosphonic and sulfonic acids, lower homologues of baclofen, phaclofen and saclofen respectively. The chlorinated acids were all weak specific antagonists of GABA at the GABAB receptor, with the sulfonic acid (pA(2) 4.0) being stronger than the phosphonic acid (pA(2) 3.8) and carboxylic acid (pA(2) 3.5).

Phase-dependent fluorescence linewidth narrowing in a three-level atom damped by a finite-bandwidth squeezed vacuum

We examine subnatural phase-dependent linewidths in the fluorescence spectrum of a three-level atom damped by a narrow-bandwidth squeezed vacuum in a cavity. Using the dressed-atom model approach of a strongly driven three-level cascade system, we derive the master equation of the system from which we obtain simple analytical expressions for the fluorescence spectrum. We show that the phase effects depend on the bandwidths of the squeezed vacuum and the cavity relative to the Rabi frequency of the driving fields. When the squeezing bandwidth is much larger than the Rabi frequency, the spectrum consists of five lines with only the central and outer sidebands dependent on the phase. For a squeezing bandwidth much smaller than the Rabi frequency the number of lines in the spectrum and their phase properties depend on the frequency at which the squeezing and cavity modes are centered. When the squeezing and cavity modes are centered on the inner Rabi sidebands, the spectrum exhibits five lines that are completely independent of the squeezing phase with only the inner Rabi sidebands dependent on the squeezing correlations. Matching the squeezing and cavity modes to the outer Rabi sidebands leads to the disappearance of the inner Rabi sidebands and a strong phase dependence of the central line and the outer Rabi sidebands. We find that in this case the system behaves as an individual two-level system that reveals exactly the noise distribution in the input squeezed vacuum. [S1050-2947(97)00111-X].