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.

A novel method for tailoring the pore-opening size of MCM-41 materials

A novel pore tailoring method is proposed by which the pore-opening sizes of MCM-41 materials can be finely tuned without significant loss in pore volume and surface area.

Numerical and experimental study of seepage in unconfined aquifers with a periodic boundary condition

The assessment of groundwater conditions within an unconfined aquifer with a periodic boundary condition is of interest in many hydrological and environmental problems. A two-dimensional numerical model for density dependent variably saturated groundwater flow, SUTRA (Voss, C.I., 1984. SUTRA: a finite element simulation model for saturated-unsaturated, fluid-density dependent ground-water flow with energy transport or chemically reactive single species solute transport. US Geological Survey, National Center, Reston, VA) is modified in order to be able to simulate the groundwater flow in unconfined aquifers affected by a periodic boundary condition. The basic flow equation is changed from pressure-form to mixed-form. The model is also adjusted to handle a seepage-face boundary condition. Experiments are conducted to provide data for the groundwater response to the periodic boundary condition for aquifers with both vertical and sloping faces. The performance of the numerical model is assessed using those data. The results of pressure- and mixed-form approximations are compared and the improvement achieved through the mixed-form of the equation is demonstrated. The ability of the numerical model to simulate the water table and seepage-face is tested by modelling some published experimental data. Finally the numerical model is successfully verified against present experimental results to confirm its ability to simulate complex boundary conditions like the periodic head and the seepage-face boundary condition on the sloping face. (C) 1999 Elsevier Science B.V. All rights reserved.

A method of analysing stream channel response to environmental change: gauge data for the Tully River

Gauging data are available from numerous streams throughout Australia, and these data provide a basis for historical analysis of geomorphic change in stream channels in response to both natural phenomena and human activities. We present a simple method for analysis of these data, and a briefcase study of an application to channel change in the Tully River, in the humid tropics of north Queensland. The analysis suggests that this channel has narrowed and deepened, rather than aggraded: channel aggradation was expected, given the intensification of land use in the catchment, upstream of the gauging station. Limitations of the method relate to the time periods over which stream gauging occurred; the spatial patterns of stream gauging sites; the quality and consistency of data collection; and the availability of concurrent land-use histories on which to base the interpretation of the channel changes.

Encapsulation of lemon oil by paste method using beta-cyclodextrin: Encapsulation efficiency and profile of oil volatiles

Microencapsulation of lemon oil was undertaken by kneading with beta-cyclodextrin, at a beta-cyclodextrin to lemon oil ratio of 88:12 (w/w). The resulting paste samples of the complex were vacuum- or spray-dried. Ten selected lemon oil flavor volatiles (alpha-pinene, sabinene, beta-pinene, beta-myrcene, limonene, gamma-terpinene, terpinolene, linalool, neral, and geranial) in the complex were analyzed periodically after 1, 2, 5, 10, 15, 20, and 30 min of kneading time. The results indicated that the levels of these volatiles were not significantly different (P > 0.05) irrespective of mixing time or type of the drying (vacuum- or spray-drying) used. An optimum mixing time was found to be 15 min, at which time the maximum encapsulation of lemon oil (97.7 mg/g of beta-cyclodextrin) was obtained in the complex powder.

Comparison of trends in method-specific suicide rates in Australia and England & Wales, 1968-97

Objective: To compare secular trends in method-specific suicide rates among young people in Australia and England & Wales between 1968 and 1997. Methods: Australian data were obtained from the Australian Bureau of Statistics, and for England & Wales from the Office for National Statistics. Overall and method-specific suicide rates for 15-34 year old males and females were calculated using ICD codes E950-9 and E980-9 except E988.8. Results: In both settings, suicide rates have almost doubled in young males over the past 30 years (from 16.8 to 32.9 per 100,000 in Australia and from 10.1 to 19.0 in England & Wales). Overall rates have changed little in young females. In both sexes and in both settings there have been substantial increases in suicide by hanging (5-7 fold increase in Australia and four-fold increase in England & Wales). There have also been smaller increases in gassing in the 1980s and '90s. In females, the impact of these increases on overall rates has been offset by a decline in drug overdose, the most common method in females. Conclusions: Rates of male suicide have increased substantially in both settings in recent years, and hanging has become an increasingly common method of suicide. The similarity in observed trends in both settings supports the view that such changes may have common causes. Research should focus on understanding why hanging has increased in popularity and what measures may be taken to diminish it.

Computation of bifurcation boundaries for power systems: A new Delta-plane method

This paper is devoted to the problems of finding the load flow feasibility, saddle node, and Hopf bifurcation boundaries in the space of power system parameters. The first part contains a review of the existing relevant approaches including not-so-well-known contributions from Russia. The second part presents a new robust method for finding the power system load flow feasibility boundary on the plane defined by any three vectors of dependent variables (nodal voltages), called the Delta plane. The method exploits some quadratic and linear properties of the load now equations and state matrices written in rectangular coordinates. An advantage of the method is that it does not require an iterative solution of nonlinear equations (except the eigenvalue problem). In addition to benefits for visualization, the method is a useful tool for topological studies of power system multiple solution structures and stability domains. Although the power system application is developed, the method can be equally efficient for any quadratic algebraic problem.

A two-dimensional analytical solution of groundwater responses to tidal loading in an estuary and ocean

Previous studies on tidal dynamics of coastal aquifers have focussed on the inland propagation of oceanic tides in the cross-shore direction, a configuration that is essentially one-dimensional. Aquifers at natural coasts can also be influenced by tidal waves in nearby estuaries, resulting in a more complex behaviour of head fluctuations in the aquifers. We present an analytical solution to the two-dimensional depth-averaged groundwater flow equation for a semi-infinite aquifer subject to oscillating head conditions at the boundaries. The solution describes the tidal dynamics of a coastal aquifer that is adjacent to a cross-shore estuary. Both the effects of oceanic and estuarine tides on the aquifer are included in the solution. The analytical prediction of the head fluctuations is verified by comparison with numerical solutions computed using a standard finite-difference method. An essential feature of the present analytical solution is the interaction between the cross- and along-shore tidal waves in the aquifer area near the estuary's entry. As the distance from the estuary or coastline increases, the wave interaction is weakened and the aquifer response is reduced, respectively, to the one-dimensional solution for oceanic tides or the solution of Sun (Sun H. A two-dimensional analytical solution of groundwater response to tidal loading in an estuary, Water Resour Res 1997;33:1429-35) for two-dimensional non-interacting tidal waves. (C) 2000 Elsevier Science Ltd. All rights reserved.

Neurons in the hilus region of the rat hippocampus are depleted in number by exposure to alcohol during early postnatal life

We have previously shown that exposing rats to a relatively high dose of ethanol during early postnatal life resulted in a deficit in spatial learning ability. This ability is controlled, at least in part, by the hippocampal formation. The purpose of the present study was to determine whether exposure of rats to ethanol during early postnatal life affected the number of specific neurons in the hippocampus. Wistar rats were exposed to a relatively high daily dose of ethanol between postnatal days 10 and 15 by placing them for 3 h each day in a chamber containing ethanol vapor. The blood ethanol concentration was about 430 mg/dl at the end of the exposure period. Groups of ethanol-treated (ET) rats, separation controls (SC), and mother-reared controls (MRC) were anesthetized and killed at 16 days of age by perfusion with phosphate-buffered glutaraldehyde (2.5%). The Cavalieri principle was used to determine the volume of various subdivisions of the hippocampal formation (CA1, CA2+CA3, hilus, and granule cell layer), and the physical disector method was used to estimate the numerical densities of neurons within each subdivision. The total number of neurons was calculated by multiplying estimates of the numerical density with the volume. There were, on average, about 441,000 granule cells in the granule cell layer and 153,000 to 177,000 pyramidal cells in both the CA1 and CA2+CA3 regions in all three treatment groups. In the hilus region, ET rats had about 27,000 neuronal cells. This was significantly fewer than the average of 38,000 such neurons estimated to be present in both MRC and SC animals. Thus, neurons in the hilus region may be particularly vulnerable to the effects of a high dose of ethanol exposure during early postnatal life. (C) 2000 Wiley-Liss, Inc.

Novel silica gel supported TiO2 photocatalyst synthesized by CVD method

TiO2 in anatase crystal phase is a very effective catalyst in the photocatalytic oxidation of organic compounds in water. To improve the recovery rate of TiO2 photocatalysts, which in most cases are in fine powder form, the chemical vapor deposition (CVD) method was used to load TiO2 onto a bigger particle support, silica gel. The amount of titania coating was found to depend strongly on the synthesis parameters of carrier gas flow rate and coating time. XPS and nitrogen ads/desorption results showed that most of the TiO2 particles generated from CVD were distributed on the external surface of the support and the coating was stable. The photocatalytic activities of TiO2/silica gel with different amounts of titania were evaluated for the oxidation of phenol aqueous solution and compared with that of Degussa P25. The optimum titania loading rate was found around 6 wt % of the TiO2 bulk concentration. Although the activity of the best TiO2/silica gel sample was still lower than that of P25, the synthesized TiO2/silica gel catalyst can be easily separated from the treated water and was found to maintain its TiO2 content and catalytic activity.

A simplified, sequential, phosphorus fractionation method

Hedley er al. (1982) developed what has become the most widely used land modified), phosphorus (P) fractionation technique. It consists of sequential extraction of increasingly less phytoavailable P pools. Extracts are centrifuged at up to 25000 g (RCF) and filtered to 0.45 mu m to ensure that soil is not lost between extractions. In attempting to transfer this method to laboratories with limited facilities, it was considered that access to high-speed centrifuges, and the cost of frequent filtration may prevent adoption of this P fractionation technique. The modified method presented here was developed to simplify methodology, reduce cost, and therefore increase accessibility of P fractionation technology. It provides quantitative recovery of soil between extractions, using low speed centrifugation without filtration. This is achieved by increasing the ionic strength of dilute extracts, through the addition of NaCl, to flocculate clay particles. Addition of NaCl does not change the amount of P extracted. Flocculation with low speed centrifugation produced extracts comparable with those having undergone filtration (0.025 mu m). A malachite green colorimetric method was adopted for inorganic P determination, as this simple manual method provides high sensitivity with negligible interference from other anions. This approach can also be used for total P following digestion, alternatively non-discriminatory methods, such as inductively coupled plasma atomic emission spectroscopy, may be employed.

Numerical modelling of single-layer folding: clarification of an issue regarding the possible effect of computer codes and the influence of initial irregularities

The influence of initial perturbation geometry and material propel-ties on final fold geometry has been investigated using finite-difference (FLAC) and finite-element (MARC) numerical models. Previous studies using these two different codes reported very different folding behaviour although the material properties, boundary conditions and initial perturbation geometries were similar. The current results establish that the discrepancy was not due to the different computer codes but due to the different strain rates employed in the two previous studies (i.e. 10(-6) s(-1) in the FLAC models and 10(-14) s(-1) in the MARC models). As a result, different parts of the elasto-viscous rheological field were bring investigated. For the same material properties, strain rate and boundary conditions, the present results using the two different codes are consistent. A transition in Folding behaviour, from a situation where the geometry of initial perturbation determines final fold shape to a situation where material properties control the final geometry, is produced using both models. This transition takes place with increasing strain rate, decreasing elastic moduli or increasing viscosity (reflecting in each case the increasing influence of the elastic component in the Maxwell elastoviscous rheology). The transition described here is mechanically feasible but is associated with very high stresses in the competent layer (on the order of GPa), which is improbable under natural conditions. (C) 2000 Elsevier Science Ltd. All rights reserved.