Accurate velocity reconstruction for Discontinuous Galerkin approximations of two-phase porous media flows

We consider Discontinuous Galerkin approximations of two-phase, immiscible porous media flows in the global pressure/fractional flow formulation with capillary pressure. A sequential approach is used with a backward Euler step for the saturation equation, equal-order interpolation for the pressure and the saturation, and without any limiters. An accurate total velocity field is recovered from the global pressure equation to be used in the saturation equation. Numerical experiments show the advantages of the proposed reconstruction. To cite this article: A. Ern et al., C R. Acad. Sci. Paris, Ser. 1347 (2009). (C) 2009 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.

Consolidation of elastic-plastic saturated porous media by the boundary element method

This paper presents a domain boundary element formulation for inelastic saturated porous media with rate-independent behavior for the solid skeleton. The formulation is then applied to elastic-plastic behavior for the solid. Biot`s consolidation theory, extended to include irreversible phenomena is considered and the direct boundary element technique is used for the numerical solution after time discretization by the implicit Euler backward algorithm. The associated nonlinear algebraic problem is solved by the Newton-Raphson procedure whereby the loading/unloading conditions are fully taken into account and the consistent tangent operator defined. Only domain nodes (nodes defined inside the domain) are used to represent all domain values and the corresponding integrals are computed by using an accurate sub-elementation scheme. The developments are illustrated through the Drucker-Prager elastic-plastic model for the solid skeleton and various examples are analyzed with the proposed algorithms. (c) 2008 Elsevier B.V. All rights reserved.

Discontinuous Galerkin approximation of two-phase flows in heterogeneous porous media with discontinuous capillary pressures

We design and investigate a sequential discontinuous Galerkin method to approximate two-phase immiscible incompressible flows in heterogeneous porous media with discontinuous capillary pressures. The nonlinear interface conditions are enforced weakly through an adequate design of the penalties on interelement jumps of the pressure and the saturation. An accurate reconstruction of the total velocity is considered in the Raviart-Thomas(-Nedelec) finite element spaces, together with diffusivity-dependent weighted averages to cope with degeneracies in the saturation equation and with media heterogeneities. The proposed method is assessed on one-dimensional test cases exhibiting rough solutions, degeneracies, and capillary barriers. Stable and accurate solutions are obtained without limiters. (C) 2010 Elsevier B.V. All rights reserved.

Multiphase flow in the vascular system of wood: From microscopic exploration to 3-D Lattice Boltzmann experiments

This paper provides insights into liquid free water dynamics in wood vessels based on Lattice Boltzmann experiments. The anatomy of real wood samples was reconstructed from systematic 3-D analyses of the vessel contours derived from successive microscopic images. This virtual vascular system was then used to supply fluid-solid boundary conditions to a two-phase Lattice Boltzmann scheme and investigate capillary invasion of this hydrophilic porous medium. Behavior of the liquid phase was strongly dependent on anatomical features, especially vessel bifurcations and reconnections. Various parameters were examined in numerical experiments with ideal vessel bifurcations, to clarify our interpretation of these features. (c) 2010 Elsevier Ltd. All rights reserved.

Sphericity of apatite particles determined by gas permeability through packed beds

Because shape is an assessment of the three-dimensional form of a particle, it may be described in terms of sphericity (Psi), which is a measure of how closely a particle approaches a spherical configuration. In this study, Darcy`s law and the Kozeny-Carman model for fluid flow through porous media were applied to packed beds to determine the sphericity (Psi) of apatite particles. The beds were composed of glass spheres or particles of apatite (igneous from Brazil and sedimentary from the United States) of three classes of size (Class 1: -297 +210 mu m; Class 2: -210 +149 mu m; Class 3: -149 +105 mu m). Glass spheres were used to validate the model because of its known sphericity (Psi = 1.00). Apatite particles, either igneous or sedimentary, showed very close values for particle sphericity (Psi approximate to 0.6). Observations on particle images conducted by scanning electron microscopy illustrated that igneous (Psi = 0.623) and sedimentary (Psi = 0.644) particles of apatite of Class 2 predominantly exhibit elongated shape. The close value of particle sphericity (Psi approximate to 0.6) showed by either igneous or sedimentary apatite may be justified by the similarity in particle shape.

New Analytic Solution Related to the Richards, Philip, and Green-Ampt Equations for Infiltration

Based on physical laws of similarity, an analytic solution of the soil water potential form of the Richards equation was derived for water infiltration into a homogeneous sand. The derivation assumes a similarity between the soil water retention function and that of the soil water content profiles taken at fixed times. The new solution successfully described soil water content profiles experimentally measured for water infiltrating downward, upward, and horizontally into a homogeneous sand and agrees with that presented by Philip in 1957. The utility of this analysis is still to be verified, but it is expected to hold for soils that have a narrow pore-size distribution before wetting and that manifest a sharp increase of water content at the wetting front during infiltration. The effect of van Genuchten`s parameters alpha and n on the application of the solution to other porous media was investigated. The solution also improves and provides a more realistic description of the infiltration process than that pioneered by Green and Ampt in 1911.

Alternative Analytical Expressions for the General van Genuchten-Mualem and van Genuchten-Burdine Hydraulic Conductivity Models

The van Genuchten expressions for the unsaturated soil hydraulic properties, first published in 1980, are used frequently in various vadose zone flow and transport applications assuming a specific relationship between the m and n soil hydraulic parameters. By comparison, probably because of the complexity of the hydraulic conductivity equations, the more general solutions with independent m and n values are rarely used. We expressed the general van Genuchten-Mualem and van Genuchten-Burdine hydraulic conductivity equations in terms of hypergeometric functions, which can be approximated by infinite series that converge rapidly for relatively large values of the van Genuchten-Mualem parameter n but only very slowly when n is close to one. Alternative equations were derived that provide very close approximations of the analytical results. The newly proposed equations allow the use of independent values of the parameters m and n in the soil water retention model of van Genuchten for subsequent prediction of the van Genuchten-Mualem and van Genuchten-Burdine hydraulic conductivity models, thus providing more flexibility in fitting experimental pressure-head-dependent water content, theta(h), and hydraulic conductivity, K(h), or K(theta) data.

Modeling of transpiration reduction in van Genuchten-Mualem type soils

We derive an analytic expression for the matric flux potential (M) for van Genuchten-Mualem (VGM) type soils which can also be written in terms of a converging infinite series. Considering the first four terms of this series, the accuracy of the approximation was verified by comparing it to values of M estimated by numerical finite difference integration. Using values of the parameters for three soils from different texture classes, the proposed four-term approximation showed an almost perfect match with the numerical solution, except for effective saturations higher than 0.9. Including more terms reduced the discrepancy but also increased the complexity of the equation. The four-term equation can be used for most applications. Cases with special interest in nearly saturated soils should include more terms from the infinite series. A transpiration reduction function for use with the VGM equations is derived by combining the derived expression for M with a root water extraction model. The shape of the resulting reduction function and its dependency on the derivative of the soil hydraulic diffusivity D with respect to the soil water content theta is discussed. Positive and negative values of dD/d theta yield concave and convex or S-shaped reduction functions, respectively. On the basis of three data sets, the hydraulic properties of virtually all soils yield concave reduction curves. Such curves based solely on soil hydraulic properties do not account for the complex interactions between shoot growth, root growth, and water availability.

NMR imaging of fluid pathways during drainage of softwood in a pressure membrane chamber

An experiment was implemented to study fluid flow in a pressure media. This procedure successfully combines nuclear magnetic resonance imaging with a pressure membrane chamber in order to visualize the non-wetting and wetting fluid flows with controlled boundary conditions. A specially designed pressure membrane chamber, made of non-magnetic materials and able to withstand 4 MPa, was designed and built for this purpose. These two techniques were applied to the drainage of Douglas fir sapwood. In the study of the longitudinal flow, narrow drainage fingers are formed in the latewood zones. They follow the longitudinal direction of wood and spread throughout the sample length. These fingers then enlarge in the cross-section plane and coalesce until drainage reaches the whole latewood part. At the end of the experiments, when the drainage of liquid water in latewood is completed, just a few sites of percolation appear in earlywood zones. This difference is a result of the wood anatomical structure, where pits, the apertures that allow the sap to flow between wood cells, are more easily aspirated in earlywood than in latewood. (C) 2007 Elsevier Ltd. All rights reserved.

Basic oxygen furnace slag as a treatment material for pathogens: Contribution of inactivation and attachment in virus attenuation

Basic oxygen furnace (BOF) slag media were studied as a potential treatment material in on-site sanitation systems. Batch and column studies were conducted to evaluate attenuation of the bacteriophage PR772 and 0.190 mu m diameter microspheres by BOF media, and to delineate the relative contributions of two principle processes of virus attenuation: inactivation and attachment. In the batch studies, conducted at 4 degrees C, substantial inactivation of PR772 did not occur in the pH 7.6 and 9.5 suspensions. At pH 11.4, bimodal inactivation of PR772 was observed, at an initial rate of 2.1 log C/C(0) day(-1) for the first two days, followed by a much slower rate of 0.124 log C/C(0) day(-1) over the following 10 days. Two column studies were conducted at 4 degrees C at a flow rate of 1 pore volume day(-1) using two slag sources (Stelco, Ontario; Tubarao, Brazil) combined with sand and pea gravel. In both column experiments, the effluent microsphere concentration approached input concentrations over time (reductions of 0.1-0.2 log C/C(0)), suggesting attachment processes for microspheres were negligible. Removal of PR772 virus was more pronounced both during the early stages of the experiments, but also after longer transport times (0.5-1.0 log C/C(0)). PR772 reduction appeared to be primarily as a result of virus inactivation in response to the elevated pH conditions generated by the BOF mixture (10.6-11.4). On-site sanitation systems using BOF media should be designed to maintain sufficient contact time between the BOF media and the wastewater to allow sufficient residence time of pathogens at elevated pH conditions. (C) 2009 Published by Elsevier Ltd.

Sorption aspects for unconsolidated materials of the Sao Carlos region (SP), Brazil

The transport of liquid and gaseous pollutants through porous geological media depends on the physical and chemical characteristics of the unconsolidated material, rocks and water associated with the characteristics of the pollutants. Of these characteristics, the sorption aspect is of fundamental importance and is a function of the mineral proportions, pH, Eh and void aspects encountered in the porous media. In the Sao Carlos region, located in the eastern-central part of the 9 ate of Sao Paulo, Brazil, there are basically two types of unconsolidated materials: the first is a residual from sandstones cemented with fines and the secondarily composed of claystones, siltstones and conglomerates from the Cretaceous Period that constitute the Itaqueri Formation; the second is a sandy sediment of the Tertiary Period. These geological conditions are found in areas where chemical products are disposed of characterized as either diffuse or point pollutions sources. Because of this situation, a study was developed to evaluate the sorption aspects of some inorganic cations that are frequently found in these sources, in varied concentrations. Taken into consideration were their physical/chemical properties, such as: specific weight, grain size, mineralogy, cationic exchange capacity, pH, hydraulic conductivity. Batch tests were run using solutions of KCl, ZnCl(2) and CuCl center dot H(2)O at three different pH values, and then with a combined solution (KCl + ZnCl(2) + CuCl center dot H(2)O), also at three different pH values.

An alternative approach for the determination of soil water mobility

A new laboratory method was proposed to establish an easily performed standard for the determination of mobile soil water close to real conditions during the infiltration and redistribution of water in a soil. It consisted of applying a water volume with a tracer ion on top of an undisturbed ring sample on a pressure plate under a known suction or pressure head. Afterwards, soil water mobility was determined by analyzing the tracer-ion concentration in the soil sample. Soil water mobility showed to be a function of the applied water volume. No relation between soil water mobility and applied pressure head could be established with data from the present experiment. A simple one- or two-parameter equation can be fitted to the experimental data to parameterize soil water mobility as a function of applied solute volume. Sandy soils showed higher mobility than loamy soils at low values of applied solute volumes, and both sandy and loamy soils showed an almost complete mobility at high applied solute volumes.

Culture and exopolysaccharide production from sugarcane molasses by Gordonia polyisoprenivorans CCT 7137, isolated from contaminated groundwater in Brazil

Gordonia polyisoprenivorans CCT 7137 was isolated from groundwater contaminated with leachate in an old controlled landfill (Sauo Paulo, Brazil), and cultured in GYM medium at different concentrations of sugarcane molasses (2%, 6%, and 10%). The strain growth was analyzed by monitoring the viable cell counts (c.f.u. mL(-1)) and optical density and EPS production was evaluated at the end of the exponential phase and 24 h after it. The analysis of the viable cell counts showed that the medium that most favored bacterial growth was not the one that favored EPS production. The control medium (GYM) was the one that most favored the strain growth, at the maximum specific growth rate of 0.232 h(-1). Differences in bacterial growth when cultured at three different concentrations of molasses were not observed. Production of EPS, in all culture media used, began during the exponential phase and continued during the growth stationary phase. The highest total EPS production, after 24 h of stationary phase, was observed in 6% molasses medium (172.86 g L(-1)) and 10% (139.47 g L(-1)) and the specific total EPS production was higher in 10% molasses medium (39.03 x10(-11)g c.f.u.(-1)). After the exponential phase, in 2%, 6%, and 10% molasses media, a higher percentage of free exopolysaccharides (EPS) was observed, representing 88.4%, 62.4%, and 64.2% of the total, respectively. A different result was observed in pattern medium, which presented EPS made up of higher percentage of capsular EPS (66.4% of the total). This work is the first study on EPS production by G. polyisoprenivorans strain in GYM medium and in medium utilizing sugarcane molasses as the sole nutrient source and suggests its potential use for EPS production by G. polyisoprenivorans CCT 7137 aiming at application in biotechnological processes.

Nonlinear diffraction of light beams propagating in photorefractive media with embedded reflecting wire

The theory of nonlinear diffraction of intensive light beams propagating through photorefractive media is developed. Diffraction occurs on a reflecting wire embedded in the nonlinear medium at a relatively small angle with respect to the direction of the beam propagation. It is shown that this process is analogous to the generation of waves by a flow of a superfluid past an obstacle. The ""equation of state"" of such a superfluid is determined by the nonlinear properties of the medium. On the basis of this hydrodynamic analogy, the notion of the ""Mach number"" is introduced where the transverse component of the wave vector plays the role of the fluid velocity. It is found that the Mach cone separates two regions of the diffraction pattern: inside the Mach cone oblique dark solitons are generated and outside the Mach cone the region of ""optical ship waves"" (the wave pattern formed by a two-dimensional packet of linear waves) is situated. Analytical theory of the ""optical ship waves"" is developed and two-dimensional dark soliton solutions of the generalized two-dimensional nonlinear Schrodinger equation describing the light beam propagation are found. Stability of dark solitons with respect to their decay into vortices is studied and it is shown that they are stable for large enough values of the Mach number.

Catalytic oxidation of ethanol on gold electrode in alkaline media

This work presents a study of the catalytic oxidation of ethanol on polycrystalline gold electrode in alkaline media. The investigation was carried out by means of chronoamperometry, cyclic voltammetry, and in situ FTIR spectroscopy. The main goal was to investigate the early stages of ethanol electrooxidation, namely at fairly low potentials (E = 600 mV vs. RHE) and for moderate reaction times (t < 300 s). Chronoamperometric experiments show a current increase accompanying the increasing in the ethanol concentration up to about 2 M and then a slight decrease at 3 M. Adsorbed CO has been observed as early as about 200 mV vs. RHE and indicates that the cleavage of the C-C bond might occur, probably to a small extent, at very low overpotentials during ethanol adsorption on gold surface. The amount of dissolved acetate ions produced during the chronoamperomentry was followed by the asymmetric stretching band at 1558 cm(-1) as a function of time, and found to increase linearly with time up to 300 s. This allowed estimating the reaction order of acetate formation with respect to ethanol concentration.