A consistent point-searching algorithm for solution interpolation in unstructured meshes consisting of 4-node bilinear quadrilateral elements

To translate and transfer solution data between two totally different meshes (i.e. mesh 1 and mesh 2), a consistent point-searching algorithm for solution interpolation in unstructured meshes consisting of 4-node bilinear quadrilateral elements is presented in this paper. The proposed algorithm has the following significant advantages: (1) The use of a point-searching strategy allows a point in one mesh to be accurately related to an element (containing this point) in another mesh. Thus, to translate/transfer the solution of any particular point from mesh 2 td mesh 1, only one element in mesh 2 needs to be inversely mapped. This certainly minimizes the number of elements, to which the inverse mapping is applied. In this regard, the present algorithm is very effective and efficient. (2) Analytical solutions to the local co ordinates of any point in a four-node quadrilateral element, which are derived in a rigorous mathematical manner in the context of this paper, make it possible to carry out an inverse mapping process very effectively and efficiently. (3) The use of consistent interpolation enables the interpolated solution to be compatible with an original solution and, therefore guarantees the interpolated solution of extremely high accuracy. After the mathematical formulations of the algorithm are presented, the algorithm is tested and validated through a challenging problem. The related results from the test problem have demonstrated the generality, accuracy, effectiveness, efficiency and robustness of the proposed consistent point-searching algorithm. Copyright (C) 1999 John Wiley & Sons, Ltd.

Hybrid modelling of coupled pore fluid-solid deformation problems

A hybrid formulation for coupled pore fluid-solid deformation problems is proposed. The scheme is a hybrid in the sense that we use a vertex centered finite volume formulation for the analysis of the pore fluid and a particle method for the solid in our model. The pore fluid formally occupies the same space as the solid particles. The size of the particles is not necessarily equal to the physical size of materials. A finite volume mesh for the pore fluid flow is generated by Delaunay triangulation. Each triangle possesses an initial porosity. Changes of the porosity are specified by the translations of the mass centers of particles. Net pore pressure gradients are applied to the particle centers and are considered in the particle momentum balance. The potential of our model is illustrated by means of a simulation of coupled fracture and fluid flow developed in porous rock under biaxial compression condition.

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.

Copolymers obtained by the radiation-induced grafting of styrene onto poly(tetrafluoroethylene-co-perfluoropropylvinyl ether) substrates. 1. Preparation and structural investigation

A study has been made to investigate the radiation grafting of styrene onto poly(tetrafluoroethylene-co-perfluoropropylvinyl ether) (PFA) substrates, using the simultaneous irradiation method. Two PFA polymers of different comonomer perfluoropropyl vinyl ether (PPVE) content and degree of crystallinity were used. Effects of grafting conditions such as monomer concentrations, type of solvent, dose rate, and irradiation dose on the grafting yield were investigated. Of the six different solvents used, the most efficient in terms of increasing grafting yield were dichloromethane, benzene, and methanol. The degree of grafting increased with increasing radiation dose up to 500 kGy, stabilizing above this dose. However, the grafting yield decreased with an increase in the dose rate. The grafting of styrene onto the PFA substrates was confirmed by FTIR-ATR and micro-Raman spectroscopy, The increase in the overall grafting yield was accompanied by a proportional increase in the penetration depth of the grafts into the substrate.

Spectroscopic study of the penetration depth of grafted polystyrene onto poly (tetrafluoroethylene-co-perfluoropropylvinylether) substrates. 1. Effect of grafting conditions

This study concerns the radiation grafting of styrene onto poly(tetrafluoroethylene-co-perfluoropropylvinylether) (PFA) substrates and the penetration depth of the graft. Grafting was obtained by the simultaneous irradiation method, and the spectroscopic analysis was made with the micro-Raman technique. Effects of grafting conditions such as the type of solvent, dose rate, and irradiation dose on the grafting yield were investigated. Of the different solvents used, the most efficient in terms of increasing grafting yield were dichloromethane, benzene, and methanol, respectively. A mixture of methanol and dichloromethane used as a solvent for styrene achieved a higher degree of grafting and concentration of grafted polystyrene onto the surface of PFA substrates than solutions of the monomer in the separate solvents. The degree of grafting increased with increasing radiation dose up to 500 kGy, stabilizing above this dose. However, the grafting yield decreased with an increase in the dose rate. The increase in the overall grafting yield was accompanied by a proportional increase in the penetration depth of the grafts into the substrate. (C) 2002 Wiley Periodicals, Inc.

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.

Characterisation of grafted supports used for solid-phase synthesis

A series of 'pellicular' type supports were fabricated by direct gamma-radiation-mediated graft polymerisation of styrene onto polypropylene, followed by aminomethylation. Raman spectroscopy was used for measuring the level of penetration of polystyrene graft into polypropylene, and other structural features such as density of graft and depth of functionalisation. The kinetics of the coupling of fluorenylmethylcarbamate (Fmoc)-labelled amino acids, to the aminomethylated polystyrene grafts have been measured by UV absorption followed cleavage of the Fmoc chromophore. The Raman spectroscopy results showed that for this series of experiments the calculated rate coefficient for coupling of Fmoc-labelled amino acids was primarily dependent on graft thickness, but was also influenced by the proportion of polystyrene graft to polypropylene. In general, it was also shown that with increasing loading capacity of support the calculated rate coefficient for amino-acid coupling decreased correspondingly. In addition, a support that had both a high rate coefficient and a high loading capacity was prepared from polypropylene base material with a co-continuous porous structure (high surface area). (C) 2003 Society of Chemical Industry.