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.

Effect of removing user fees on attendance for curative and preventive primary health care services in rural South Africa

User fees are used to recover costs and discourage unnecessary attendance at primary care clinics in many developing countries. In South Africa, user fees for children aged under 6 yea rs and pregnant women were removed in 1994, and in 1997 all user fees at all primary health care clinics were abolished. The intention of these policy changes was to improve access to health services for previously disadvantaged communities. We investigated the impact of these changes on clinic attendance patterns in Hlabisa health district. Average quarterly new registrations and total attendances for preventive services (antenatal care, immunization, growth monitoring) and curative services (treatment of ailments) at a mobile primary health care unit were studied from 1992 to 1998. Regression analysis was undertaken to assess whether trends were statistically significant. There was a sustained increase in new registrations (P = 0.0001) and total attendances (P = 0.0001)for curative services, and a fall in new registrations (P = 0.01) and total attendances for immunization and growth monitoring (P = 0.0002) over the study period. The upturn in demand for curative services started at the time of the first policy change. The decreases in antenatal registrations (P = 0.07) and attendances (P = 0.09) were not statistically significant The number of new registrations for immunization and growth monitoring increased following the first policy change but declined thereafter. We found no evidence that the second policy change influenced underlying trends. The removal of user fees improved access to curative services but this may have happened at the expense of some preventive services. Governments should remain vigilant about the effects of new health policies in order to ensure that objectives are being met.

An equivalent algorithm for simulating thermal effects of magma intrusion problems in porous rocks

An equivalent algorithm is proposed to simulate thermal effects of the magma intrusion in geological systems, which are composed of porous rocks. Based on the physical and mathematical equivalence, 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 a physically equivalent heat source. From the analysis of an ideal solidification model, the physically equivalent heat source has been determined in this paper. The major advantage in using the proposed equivalent algorithm is that the fixed finite element mesh with a variable integration time step can be employed to simulate the thermal effect of the intruded magma solidification using the conventional finite element method. The related numerical results have demonstrated the correctness and usefulness of the proposed equivalent algorithm for simulating the thermal effect of the intruded magma solidification in geological systems. (C) 2003 Elsevier B.V. All rights reserved.

Machining of micro aspherical mould inserts

This paper reports a parametric investigation and development of grinding technologies for micro aspherical mould inserts using parallel grinding method. The parametric investigation revealed that at nanometric scale the undeformed chip thickness has little influence on the surface finish of ground inserts. The grinding trace spacing has a slightly larger influence on the surface finish. A new technique was developed to true and dress the resin bonded micro wheels with mesh size of #3000, which produced a satisfactory wheel form accuracy and relatively high grain packing density. A form error compensation technique was also developed, with which mould inserts of submicron form accuracy were consistently produced. Using the developed technologies, micro aspherical inserts of diameters ranging from 200 mu m to 1000 mu m with surface finish of around 10 nm and form error of similar to 0.2-0.4 mu m were successfully fabricated. (c) 2005 Elsevier Inc. All rights reserved.

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.

Application of Petrov-Galerkin methods to transient boundary value problems in chemical engineering: adsorption with steep gradients in bidisperse solids

Petrov-Galerkin methods are known to be versatile techniques for the solution of a wide variety of convection-dispersion transport problems, including those involving steep gradients. but have hitherto received little attention by chemical engineers. We illustrate the technique by means of the well-known problem of simultaneous diffusion and adsorption in a spherical sorbent pellet comprised of spherical, non-overlapping microparticles of uniform size and investigate the uptake dynamics. Solutions to adsorption problems exhibit steep gradients when macropore diffusion controls or micropore diffusion controls, and the application of classical numerical methods to such problems can present difficulties. In this paper, a semi-discrete Petrov-Galerkin finite element method for numerically solving adsorption problems with steep gradients in bidisperse solids is presented. The numerical solution was found to match the analytical solution when the adsorption isotherm is linear and the diffusivities are constant. Computed results for the Langmuir isotherm and non-constant diffusivity in microparticle are numerically evaluated for comparison with results of a fitted-mesh collocation method, which was proposed by Liu and Bhatia (Comput. Chem. Engng. 23 (1999) 933-943). The new method is simple, highly efficient, and well-suited to a variety of adsorption and desorption problems involving steep gradients. (C) 2001 Elsevier Science Ltd. All rights reserved.

Solution techniques for transport problems involving steep concentration gradients: application to noncatalytic fluid solid reactions

Some efficient solution techniques for solving models of noncatalytic gas-solid and fluid-solid reactions are presented. These models include those with non-constant diffusivities for which the formulation reduces to that of a convection-diffusion problem. A singular perturbation problem results for such models in the presence of a large Thiele modulus, for which the classical numerical methods can present difficulties. For the convection-diffusion like case, the time-dependent partial differential equations are transformed by a semi-discrete Petrov-Galerkin finite element method into a system of ordinary differential equations of the initial-value type that can be readily solved. In the presence of a constant diffusivity, in slab geometry the convection-like terms are absent, and the combination of a fitted mesh finite difference method with a predictor-corrector method is used to solve the problem. Both the methods are found to converge, and general reaction rate forms can be treated. These methods are simple and highly efficient for arbitrary particle geometry and parameters, including a large Thiele modulus. (C) 2001 Elsevier Science Ltd. All rights reserved.

The efficiency and condition of oysters and macroalgae used as biological filters of shrimp pond effluent

Current shrimp pond management practices generally result in elevated concentrations of nutrients, suspended solids, bacteria and phytoplankton compared with the influent water. Concerns about adverse environmental impacts caused by discharging pond effluent directly into adjacent waterways have prompted the search for cost-effective methods of effluent treatment. One potential method of effluent treatment is the use of ponds or raceways stocked with plants or animals that act as natural biofilters by removing waste nutrients. In addition to improving effluent water quality prior to discharge, the use of natural biofilters provides a method for capturing otherwise wasted nutrients. This study examined the potential of the native oyster, Saccostrea commercialis (Iredale and Roughley) and macroalgae, Gracilaria edulis (Gmelin) Silva to improve effluent water quality from a commercial Penaeus japonicus (Bate) shrimp farm, A system of raceways was constructed to permit recirculation of the effluent through the oysters to maximize the filtration of bacteria, phytoplankton and total suspended solids. A series of experiments was conducted to test the ability of oysters and macroalgae to improve effluent water quality in a flow-through system compared with a recirculating system. In the flow-through system, oysters reduced the concentration of bacteria to 35% of the initial concentration, chlorophyll a to 39%, total particulates (2.28-35.2 mum) to 29%, total nitrogen to 66% and total phosphorus to 56%. Under the recirculating flow regime, the ability of the oysters to improve water quality was significantly enhanced. After four circuits, total bacterial numbers were reduced to 12%, chlorophyll a to 4%, and total suspended solids to 16%. Efforts to increase biofiltration by adding additional layers of oyster trays and macroalgae-filled mesh bags resulted in fouling of the lower layers causing the death of oysters and senescence of macroalgae. Supplementary laboratory experiments were designed to examine the effects of high effluent concentrations of suspended particulates on the growth and condition of oysters and macroalgae. The results demonstrated that high concentrations of particulates inhibited growth and reduced the condition of oysters and macroalgae. Allowing the effluent to settle before biofiltration improved growth and reduced signs of stress in the oysters and macroalgae. A settling time of 6 h reduced particulates to a level that prevented fouling of the oysters and macroalgae.

Predicting cave initiation and propagation

The most widely used method for predicting the onset of continuous caving is Laubscher's caving chart. A detailed examination of this method was undertaken which concluded that it had limitations which may impact on results, particularly when dealing with stronger rock masses that are outside current experience. These limitations relate to inadequate guidelines for adjustment factors to rock mass rating (RMR), concerns about the position on the chart of critical case history data, undocumented changes to the method and an inadequate number of data points to be confident of stability boundaries. A review was undertaken on the application and reliability of a numerical method of assessing cavability. The review highlighted a number of issues, which at this stage, make numerical continuum methods problematic for predicting cavability. This is in particular reference to sensitivity to input parameters that are difficult to determine accurately and mesh dependency. An extended version of the Mathews method for open stope design was developed as an alternative method of predicting the onset of continuous caving. A number of caving case histories were collected and analyzed and a caving boundary delineated statistically on the Mathews stability graph. The definition of the caving boundary was aided by the existence of a large and wide-ranging stability database from non-caving mines. A caving rate model was extrapolated from the extended Mathews stability graph but could only be partially validated due to a lack of reliable data.