Theoretical and numerical analyses of convective instability in porous media with temperature-dependent viscosity

Exact analytical solutions of the critical Rayleigh numbers have been obtained for a hydrothermal system consisting of a horizontal porous layer with temperature-dependent viscosity. The boundary conditions considered are constant temperature and zero vertical Darcy velocity at both the top and bottom of the layer. Not only can the derived analytical solutions be readily used to examine the effect of the temperature-dependent viscosity on the temperature-gradient driven convective flow, but also they can be used to validate the numerical methods such as the finite-element method and finite-difference method for dealing with the same kind of problem. The related analytical and numerical results demonstrated that the temperature-dependent viscosity destabilizes the temperature-gradient driven convective flow and therefore, may affect the ore body formation and mineralization in the upper crust of the Earth. Copyright (C) 2003 John Wiley Sons, Ltd.

Analysis of soil wetting and solute transport in subsurface trickle irrigation

The increased use of trickle or drip irrigation is seen as one way of helping to improve the sustainability of irrigation systems around the world. However, soil water and solute transport properties and soil profile characteristics are often not adequately incorporated in the design and management of trickle systems. In this paper, we describe results of a simulation study designed to highlight the impacts of soil properties on water and solute transport from buried trickle emitters. The analysis addresses the influence of soil hydraulic properties, soil layering, trickle discharge rate, irrigation frequency, and timing of nutrient application on wetting patterns and solute distribution. We show that (1) trickle irrigation can improve plant water availability in medium and low permeability fine-textured soils, providing that design and management are adapted to account for their soil hydraulic properties, (2) in highly permeable coarse-textured soils, water and nutrients move quickly downwards from the emitter, making it difficult to wet the near surface zone if emitters are buried too deep, and (3) changing the fertigation strategy for highly permeable coarse-textured soils to apply nutrients at the beginning of an irrigation cycle can maintain larger amounts of nutrient near to and above the emitter, thereby making them less susceptible to leaching losses. The results demonstrate the need to account for differences in soil hydraulic properties and solute transport when designing irrigation and fertigation management strategies. Failure to do this will result in inefficient systems and lost opportunities for reducing the negative environmental impacts of irrigation.

Australia-wide predictions of soil properties using decision trees

This paper describes the construction of Australia-wide soil property predictions from a compiled national soils point database. Those properties considered include pH, organic carbon, total phosphorus, total nitrogen, thickness. texture, and clay content. Many of these soil properties are used directly in environmental process modelling including global climate change models. Models are constructed at the 250-m resolution using decision trees. These relate the soil property to the environment through a suite of environmental predictors at the locations where measurements are observed. These models are then used to extend predictions to the continental extent by applying the rules derived to the exhaustively available environmental predictors. The methodology and performance is described in detail for pH and summarized for other properties. Environmental variables are found to be important predictors, even at the 250-m resolution at which they are available here as they can describe the broad changes in soil property.

Degradation of phenanthrene and pyrene in soil slurry reactors with immobilized bacteria Zoogloea sp.

The environmental fate of polycyclic aromatic hydrocarbons (PAHs) in soils is motivated by their wide distribution, high persistence, and potentially deleterious effect on human health. Polycyclic aromatic hydrocarbons constitute the largest group of environmental contaminants released in the environment. Therefore, the potential biodegradation of these compounds is of vital importance. A biocarrier suitable for the colonization by micro-organisms for the purpose of purifying soil contaminated by polycyclic aromatic hydrocarbons was developed. The optimized composition of the biocarrier was polyvinyl alcohol (PVA) 10%, sodium alginate (SA) 0.5%, and powdered activated carbon (PAC) 5%. There was no observable cytotoxicity of biocarriers on immobilized cells and a viable cell population of 1.86 x 10(10) g(-1) was maintained for immobilized bacterium. Biocarriers made from chemical methods had a higher biodegradation but lower mechanical strengths. Immobilized bacterium Zoogloea sp. had an ideal capability of biodegradation for phenanthrene and pyrene over a relative wide concentration range. The study results showed that the biodegradation of phenanthrene and pyrene reached 87.0 and 75.4%, respectively, by using the optimal immobilized method of Zoogloea sp. cultivated in a sterilized soil. Immobilized Zoogloea sp. was found to be effective for biodegrading the soil contaminated with phenanthrene and pyrene. Even in natural (unsterilized) soil, the biodegradation of phenanthrene and pyrene using immobilized Zoogloea sp. reached 85.0 and 67.1%, respectively, after 168 h of cultivation, more than twice that achieved if the cells were not immobilized on the biocarrier. Therefore, the immobilization technology enhanced the competitive ability of introduced micro-organisms and represents an effective method for the biotreatment of soil contaminated with phenanthrene and pyrene.

Ferrolysis induced soil transformation by natural drainage in Vertisols of sub-humid South India

In sub-humid South India, recent studies have shown that black soil areas (Vertisols and vertic Intergrades), located on flat valley bottoms, have been rejuvenated through the incision of streambeds, inducing changes in the pedoclimate and soil transformation. Joint pedological, geochemical and geophysical investigations were performed in order to better understand the ongoing processes and their contribution to the chemistry of local rivers. The seasonal rainfall causes cycles of oxidation and reduction in a perched watertable at the base of the black soil, while the reduced solutions are exported through a loamy sand network. This framework favours a ferrolysis process, which causes low base saturation and protonation of clay, leading to the weathering of 2:1 then 1:1 clay minerals. Maximum weathering conditions occur at the very end of the wet season, just before disappearance of the perched watertable. Therefore, the by-products of soil transformation are partially drained off and calcareous nodules, then further downslope, amorphous silica precipitate upon soil dehydration. The ferrolysed area is fringing the drainage system indicating that its development has been induced by the streambed incision. The distribution of (14)C ages of CaCO(3) nodules suggests that the ferrolysis process started during the late Holocene, only about 2 kyr B.P. at the studied site and about 5 kyr B.P. at the watershed outlet. The results of this study are applied to an assessment of the physical erosion rate (4.8x10(-3) m/kyr) since the recent reactivation of the erosion process. (C) 2010 Elsevier B.V. All rights reserved.

Soil mineral genesis and distribution in a saline lake landscape of the Pantanal Wetland, Brazil

The origin of the saline lakes in the Pantanal wetland has been classically attributed to processes occurring in past periods. However, recent studies have suggested that saline water is currently forming from evaporative concentration of fresh water, which is provided annually by seasonal floods. Major elements (Ca, Mg, K) and alkalinity appear to be geochemically controlled during the concentration of waters and may be involved in the formation of carbonates and clay minerals around the saline lakes. The mineralogy of soils associated with a representative saline lake was investigated using XRD, TEM-EDS, and ICP-MS in order to identify the composition and genesis of the secondary minerals suspected to be involved in the control of major elements. The results showed that Ca, Mg, and K effectively undergo oversaturation and precipitation as the waters become more saline. These elements are incorporated in the authigenically formed carbonates, smectites, and micas surrounding the saline lake. The control of Ca occurs by precipitation of calcite and dolomite in nodules while Mg and K are mainly involved in the neoformation of Mg-smectites (stevensitic and saponitic minerals) and, probably, iron-enriched micas (ferric-illite) in surface and subsurface horizons. Therefore, our study confirms that the salinity of Pantanal, historically attributed to inheritance from former regimes, has a contribution of current processes. (C) 2009 Elsevier B.V. All rights reserved.

Soil morphological control on saline and freshwater lake hydrogeochemistry in the Pantanal of Nhecolandia, Brazil

Joint pedological, geochemical, hydrological and geophysical investigations were performed to study the coexistence or saline and freshwater lakes in close proximity and similar climatic conditions in the Nhecolandia region, Pantanal wetlands in Brazil. The saline lakes are concentrically surrounded by green sandy loam horizons, which cause differential hydrological regimes. Mg-calcite, K-silicates, and amorphous silica precipitate in the soil cover, whereas Mg-silicates and more soluble Na-carbonates are concentrated in the topsoil along the shore of the saline lake. In saline solutions, some minor elements (As, Se) reach values above the water quality recommendations, whereas others are controlled and incorporated in solid phases (Ba, Sr). Locally, the destruction of the sandy loam horizons generates very acidic soil solution (pH similar to 3.5) through a process not yet understood. The soil distributions indicate that some freshwater lakes are former saline lakes. They are invaded by freshwater after destruction of the sandy loam green horizons, then the freshwater becomes enriched in K(+), SO(4)(2-), Fe, Al, and a stream of minor and trace elements. The formation of these green sandy loam horizons in the saline environment and their destruction in the non-saline one emphasizes the dynamic nature of this environment (C) 2008 Elsevier B.V. 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.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

Numerical simulation of dendrite growth during solidification

A comprehensive probabilistic model for simulating dendrite morphology and investigating dendritic growth kinetics during solidification has been developed, based on a modified Cellular Automaton (mCA) for microscopic modeling of nucleation, growth of crystals and solute diffusion. The mCA model numerically calculated solute redistribution both in the solid and liquid phases, the curvature of dendrite tips and the growth anisotropy. This modeling takes account of thermal, curvature and solute diffusion effects. Therefore, it can simulate microstructure formation both on the scale of the dendrite tip length. This model was then applied for simulating dendritic solidification of an Al-7%Si alloy. Both directional and equiaxed dendritic growth has been performed to investigate the growth anisotropy and cooling rate on dendrite morphology. Furthermore, the competitive growth and selection of dendritic crystals have also investigated.

Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. (c) 2004 Elsevier Ltd. All rights reserved.

A numerical study of effect of grain boundaries on elastic and plastic properties in nanocomposite materials

Nanocomposite materials have received considerable attention in recent years due to their novel properties. Grain boundaries are considered to play an important role in nanostructured materials. This work focuses on the finite element analysis of the effect of grain boundaries on the overall mechanical properties of aluminium/alumina composites. A grain boundary is incorporated into the commonly used unit cell model to investigate its effect on material properties. By combining the unit cell model with an indentation model, coupled with experimental indentation measurements, the ''effective'' plastic property of the grain boundary is estimated. In addition, the strengthening mechanism is also discussed based on the Estrin-Mecking model.

Numerical modelling of tide-induced beach water table fluctuations

Field studies have shown that the elevation of the beach groundwater table varies with the tide and such variations affect significantly beach erosion or accretion. In this paper, we present a BEM (Boundary Element Method) model for simulating the tidal fluctuation of the beach groundwater table. The model solves the two-dimensional flow equation subject to free and moving boundary conditions, including the seepage dynamics at the beach face. The simulated seepage faces were found to agree with the predictions of a simple model (Turner, 1993). The advantage of the present model is, however, that it can be used with little modification to simulate more complicated cases, e.g., surface recharge from rainfall and drainage in the aquifer may be included (the latter is related to beach dewatering technique). The model also simulated well the field data of Nielsen (1990). In particular, the model replicated three distinct features of local water table fluctuations: steep rising phase versus flat falling phase, amplitude attenuation and phase lagging.

Numerical and experimental study of a homogenizer impinging jet

High-pressure homogenization is a key unit operation used to disrupt cells containing intracellular bioproducts. Modeling and optimization of this unit are restrained by a lack of information on the flow conditions within a homogenizer value. A numerical investigation of the impinging radial jet within a homogenizer value is presented. Results for a laminar and turbulent (k-epsilon turbulent model) jet are obtained using the PHOENICS finite-volume code. Experimental measurement of the stagnation region width and correlation of the cell disruption efficiency with jet stagnation pressure both indicate that the impinging jet in the homogenizer system examined is likely to be laminar under normal operating conditions. Correlation of disruption data with laminar stagnation pressure provides a better description of experimental variability than existing correlations using total pressure drop or the grouping 1/Y(2)h(2).