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 difference time domain (FDTD) method for modeling the effect of switched gradients on the human body in MRI

Numerical modeling of the eddy currents induced in the human body by the pulsed field gradients in MRI presents a difficult computational problem. It requires an efficient and accurate computational method for high spatial resolution analyses with a relatively low input frequency. In this article, a new technique is described which allows the finite difference time domain (FDTD) method to be efficiently applied over a very large frequency range, including low frequencies. This is not the case in conventional FDTD-based methods. A method of implementing streamline gradients in FDTD is presented, as well as comparative analyses which show that the correct source injection in the FDTD simulation plays a crucial rule in obtaining accurate solutions. In particular, making use of the derivative of the input source waveform is shown to provide distinct benefits in accuracy over direct source injection. In the method, no alterations to the properties of either the source or the transmission media are required. The method is essentially frequency independent and the source injection method has been verified against examples with analytical solutions. Results are presented showing the spatial distribution of gradient-induced electric fields and eddy currents in a complete body model.

The source of the Great Dyke, Zimbabwe, and its tectonic significance: Evidence from Re-Os isotopes

Re-Os data for chromite separates from 10 massive chromitite seams sampled along the 550-km length of the 2.58-Ga Great Dyke layered igneous complex, Zimbabwe, record initial 187Os/188Os ratios in the relatively narrow range between 0.1106 and 0.1126. This range of initial 187Os/188Os values is only slightly higher than the value for the coeval primitive upper mantle (0.1107) as modeled from the Re-Os evolution of chondrites and data of modern mantle melts and mantle derived xenoliths. Analyses of Archean granitoid and gneiss samples from the Zimbabwe Craton show extremely low Os concentrations (3-9 ppt) with surprisingly unradiogenic present-day 187Os/188Os signatures between 0.167 and 0.297. Only one sample yields an elevated 187Os/188Os ratio of 1.008. Using these data, the range of crustal contamination of the Great Dyke magma would be minimally 0%-33% if the magma source was the primitive upper mantle, whereas the range estimated from Nd and Pb isotope systematics is 5%-25%. If it is assumed that the primary Great Dyke magma derived from an enriched deep mantle reservoir (via a plume), a better agreement can be obtained. A significant contribution from a long-lived subcontinental lithospheric mantle (SCLM) reservoir with subchondritic Re/Os to the Great Dyke melts cannot be reconciled with the Os isotope results at all. However, Os isotope data on pre-Great Dyke ultramafic complexes of the Zimbabwe Craton and thermal modeling show that such an SCLM existed below the Zimbabwe Craton at the time of the Great Dyke intrusion. It is therefore concluded that large melt volumes such as that giving rise to the Great Dyke were able to pass lithospheric mantle keels without significant contamination in the late Archean. Because the ultramafic-mafic melts forming the Great Dyke must have originated below the SCLM (which extends to at least a 200-km depth ), the absence of an SCLM signature precludes a subduction-related magma-generation process.

Enhanced biological phosphorus removal in a sequencing batch reactor using propionate as the sole carbon source

An enhanced biological phosphorus removal (EBPR) system was developed in a sequencing batch reactor (SBR) using propionate as the sole carbon source. The microbial community was followed using fluorescence in situ hybridization (FISH) techniques and Candidatus 'Accumulibacter phosphatis' were quantified from the start up of the reactor until steady state. A series of SBR cycle studies was performed when 55% of the SBR biomass was Accumulibacter, a confirmed polyphosphate accumulating organism (PAO) and when Candidatus 'Competibacter phosphatis,' a confirmed glycogen-accumulating organism (GAO), was essentially undetectable. These experiments evaluated two different carbon sources (propionate and acetate), and in every case, two different P-release rates were detected. The highest rate took place while there was volatile fatty acid (VFA) in the mixed liquor, and after the VFA was depleted a second P-release rate was observed. This second rate was very similar to the one detected in experiments performed without added VFA. A kinetic and stoichiometric model developed as a modification of Activated Sludge Model 2 (ASM2) including glycogen economy, was fitted to the experimental profiles. The validation and calibration of this model was carried out with the cycle study experiments performed using both VFAs. The effect of pH from 6.5 to 8.0 on anaerobic P-release and VFA-uptake and aerobic P-uptake was also studied using propionate. The optimal overall working pH was around 7.5. This is the first study of the microbial community involved in EBPR developed with propionate as a sole carbon source along with detailed process performance investigations of the propionate-utilizing PAOs. (C) 2004 Wiley Periodicals, Inc.

Forecasting Helicoverpa populations in Australia: A comparison of regression based models and a bio-climatic based modeling approach

Long-term forecasts of pest pressure are central to the effective management of many agricultural insect pests. In the eastern cropping regions of Australia, serious infestations of Helicoverpa punctigera (Wallengren) and H. armigera (Hübner)(Lepidoptera: Noctuidae) are experienced annually. Regression analyses of a long series of light-trap catches of adult moths were used to describe the seasonal dynamics of both species. The size of the spring generation in eastern cropping zones could be related to rainfall in putative source areas in inland Australia. Subsequent generations could be related to the abundance of various crops in agricultural areas, rainfall and the magnitude of the spring population peak. As rainfall figured prominently as a predictor variable, and can itself be predicted using the Southern Oscillation Index (SOI), trap catches were also related to this variable. The geographic distribution of each species was modelled in relation to climate and CLIMEX was used to predict temporal variation in abundance at given putative source sites in inland Australia using historical meteorological data. These predictions were then correlated with subsequent pest abundance data in a major cropping region. The regression-based and bioclimatic-based approaches to predicting pest abundance are compared and their utility in predicting and interpreting pest dynamics are discussed.

Making Fedora easier to implement with Fez: A free open source content model and workflow management front-end to Fedora

The University of Queensland, Australia has developed Fez, a world-leading user-interface and management system for Fedora-based institutional repositories, which bridges the gap between a repository and users. Christiaan Kortekaas, Andrew Bennett and Keith Webster will review this open source software that gives institutions the power to create a comprehensive repository solution without the hassle..

Modeling control objectives for business process compliance

Business process design is primarily driven by process improvement objectives. However, the role of control objectives stemming from regulations and standards is becoming increasingly important for businesses in light of recent events that led to some of the largest scandals in corporate history. As organizations strive to meet compliance agendas, there is an evident need to provide systematic approaches that assist in the understanding of the interplay between (often conflicting) business and control objectives during business process design. In this paper, our objective is twofold. We will firstly present a research agenda in the space of business process compliance, identifying major technical and organizational challenges. We then tackle a part of the overall problem space, which deals with the effective modeling of control objectives and subsequently their propagation onto business process models. Control objective modeling is proposed through a specialized modal logic based on normative systems theory, and the visualization of control objectives on business process models is achieved procedurally. The proposed approach is demonstrated in the context of a purchase-to-pay scenario.

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.

Reactivity of chars and carbons: New insights through molecular modeling

A new model proposed for the gasification of chars and carbons incorporates features of the turbostratic nanoscale structure that exists in such materials. The model also considers the effect of initial surface chemistry and different reactivities perpendicular to the edges and to the faces of the underlying crystallite planes comprising the turbostratic structure. It may be more realistic than earlier models based on pore or grain structure idealizations when the carbon contains large amounts of crystallite matter. Shrinkage of the carbon particles in the chemically controlled regime is also possible due to the random complete gasification of crystallitic planes. This mechanism can explain observations in the literature of particle size reduction. Based on the model predictions, both initial surface chemistry and the number of stacked planes in the crystallites strongly influence the reactivity and particle shrinkage. Its test results agree well with literature data on the air-oxidation of Spherocarb and show that it accurately predicts the variation of particle size with conversion. Model parameters are determined entirely from rate measurements.

A state-contingent production approach to principal-agent problems with an application to point-source pollution control

Two basic representations of principal-agent relationships, the 'state-space' and 'parameterized distribution' formulations, have emerged. Although the state-space formulation appears more natural, analytical studies using this formulation have had limited success. This paper develops a state-space formulation of the moral-hazard problem using a general representation of production under uncertainty. A closed-form solution for the agency-cost problem is derived. Comparative-static results are deduced. Next we solve the principal's problem of selecting the optimal output given the agency-cost function. The analysis is applied to the problem of point-source pollution control. (C) 1998 Published by Elsevier Science S.A. All rights reserved.

Modeling the optical constants of GaP, InP, and InAs

An extension of the Adachi model with the adjustable broadening function, instead of the Lorentzian one, is employed to model the optical constants of GaP, InP, and InAs. Adjustable broadening is modeled by replacing the damping constant with the frequency-dependent expression. The improved flexibility of the model enables achieving an excellent agreement with the experimental data. The relative rms errors obtained for the refractive index equal 1.2% for GaP, 1.0% for InP, and 1.6% for InAs. (C) 1999 American Institute of Physics. [S0021-8979(99)05807-7].

Modeling the fluid-flow-induced stress and collapse in a dendritic network

An analytical approach to the stress development in the coherent dendritic network during solidification is proposed. Under the assumption that stresses are developed in the network as a result of the friction resisting shrinkage-induced interdendritic fluid flow, the model predicts the stresses in the solid. The calculations reflect the expected effects of postponed dendrite coherency, slower solidification conditions, and variations of eutectic volume fraction and shrinkage. Comparing the calculated stresses to the measured shear strength of equiaxed mushy zones shows that it is possible for the stresses to exceed the strength, thereby resulting in reorientation or collapse of the dendritic network.

Modeling the optical constants of AlxGa1-xAs alloys

The extension of Adachi's model with a Gaussian-like broadening function, in place of Lorentzian, is used to model the optical dielectric function of the alloy AlxGa1-xAs. Gaussian-like broadening is accomplished by replacing the damping constant in the Lorentzian line shape with a frequency dependent expression. In this way, the comparative simplicity of the analytic formulas of the model is preserved, while the accuracy becomes comparable to that of more intricate models, and/or models with significantly more parameters. The employed model accurately describes the optical dielectric function in the spectral range from 1.5 to 6.0 eV within the entire alloy composition range. The relative rms error obtained for the refractive index is below 2.2% for all compositions. (C) 1999 American Institute of Physics. [S0021-8979(99)00512-5].

Modeling the optical properties of AlSb, GaSb, and InSb

Optical constants of AlSb, GaSb, and InSb are modeled in the 1-6 eV spectral range. We employ an extension of Adachi's model of the optical constants of semiconductors. The model takes into account transitions at E-0, E-0 + Delta(0), E-1, and E-1 + Delta(1) critical points, as well as higher-lying transitions which are modeled with three damped harmonic oscillators. We do not consider indirect transitions contribution, since it represents a second-order perturbation and its strength should be low. Also, we do not take into account excitonic effects at E-1, E-1 + Delta(1) critical points, since we model the room temperature data. In spite of fewer contributions to the dielectric function compared to previous calculations involving Adachi's model, our calculations show significantly improved agreement with the experimental data. This is due to the two main distinguishing features of calculations presented here: use of adjustable line broadening instead of the conventional Lorentzian one, and employment of a global optimization routine for model parameter determination.