Estimation of water retention curve from mercury intrusion porosimetry and van Genuchten model

The water retention curve (WRC) is a hydraulic characteristic of concrete required for advanced modeling of water (and thus solute) transport in variably saturated, heterogeneous concrete. Unfortunately, determination by a direct experimental method (for example, measuring equilibrium moisture levels of large samples stored in constant humidity cells) is a lengthy process, taking over 2 years for large samples. A surrogate approach is presented in which the WRC is conveniently estimated from mercury intrusion porosimetry (MIP) and validated by water sorption isotherms: The well-known Barrett, Joyner and Halenda (BJH) method of estimating the pore size distribution (PSD) from the water sorption isotherm is shown to complement the PSD derived from conventional MIP. This provides a basis for predicting the complete WRC from MIP data alone. The van Genuchten equation is used to model the combined water sorption and MIP results. It is a convenient tool for describing water retention characteristics over the full moisture content range. The van Genuchten parameter estimation based solely on MIP is shown to give a satisfactory approximation to the WRC, with a simple restriction on one. of the parameters.

Effects of temperature-dependent viscosity variation on entropy generation, heat and fluid flow through a porous-saturated duct of rectangular cross-section

Effect of temperature-dependent viscosity on fully developed forced convection in a duct of rectangular cross-section occupied by a fluid-saturated porous medium is investigated analytically. The Darcy flow model is applied and the viscosity-temperature relation is assumed to be an inverse-linear one. The case of uniform heat flux on the walls, i.e. the H boundary condition in the terminology of Kays and Crawford, is treated. For the case of a fluid whose viscosity decreases with temperature, it is found that the effect of the variation is to increase the Nusselt number for heated walls. Having found the velocity and the temperature distribution, the second law of thermodynamics is invoked to find the local and average entropy generation rate. Expressions for the entropy generation rate, the Bejan number, the heat transfer irreversibility, and the fluid flow irreversibility are presented in terms of the Brinkman number, the Péclet number, the viscosity variation number, the dimensionless wall heat flux, and the aspect ratio (width to height ratio). These expressions let a parametric study of the problem based on which it is observed that the entropy generated due to flow in a duct of square cross-section is more than those of rectangular counterparts while increasing the aspect ratio decreases the entropy generation rate similar to what previously reported for the clear flow case.

A numerical study of pore-fluid, thermal and mass flow iu fluid-saturated porous rock basins

We present a numerical methodology for the study of convective pore-fluid, thermal and mass flow in fluid-saturated porous rock basins. lit particular, we investigate the occurrence and distribution pattern of temperature gradient driven convective pore-fluid flow and hydrocarbon transport in the Australian North West Shelf basin. The related numerical results have demonstrated that: (1) The finite element method combined with the progressive asymptotic approach procedure is a useful tool for dealing with temperature gradient driven pore-fluid flow and mass transport in fluid-saturated hydrothermal basins; (2) Convective pore-fluid flow generally becomes focused in more permeable layers, especially when the layers are thick enough to accommodate the appropriate convective cells; (3) Large dislocation of strata has a significant influence off the distribution patterns of convective pore;fluid flow, thermal flow and hydrocarbon transport in the North West Shelf basin; (4) As a direct consequence of the formation of convective pore-fluid cells, the hydrocarbon concentration is highly localized in the range bounded by two major faults in the basin.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Effect of material anisotropy on the onset of convective flow in three-dimensional fluid-saturated faults

In order to investigate the effect of material anisotropy on convective instability of three-dimensional fluid-saturated faults, an exact analytical solution for the critical Rayleigh number of three-dimensional convective flow has been obtained. Using this critical Rayleigh number, effects of different permeability ratios and thermal conductivity ratios on convective instability of a vertically oriented three-dimensional fault have been examined in detail. It has been recognized that (1) if the fault material is isotropic in the horizontal direction, the horizontal to vertical permeability ratio has a significant effect on the critical Rayleigh number of the three-dimensional fault system, but the horizontal to vertical thermal conductivity ratio has little influence on the convective instability of the system, and (2) if the fault material is isotropic in the fault plane, the thermal conductivity ratio of the fault normal to plane has a considerable effect on the critical Rayleigh number of the three-dimensional fault system, but the effect of the permeability ratio of the fault normal to plane on the critical Rayleigh number of three-dimensional convective flow is negligible.

Effects of hot intrusions on pore-fluid flow and heat transfer in fluid-saturated rocks

We use the finite element method to solve coupled problems between pore-fluid flow and heat transfer in fluid-saturated porous rocks. In particular, we investigate the effects of both the hot pluton intrusion and topographically driven horizontal flow on the distributions of the pore-flow velocity and temperature in large-scale hydrothermal systems. Since general mineralization patterns are strongly dependent on distributions of both the pore-fluid velocity and temperature fields, the modern mineralization theory has been used to predict the general mineralization patterns in several realistic hydrothermal systems. The related numerical results have demonstrated that: (1) The existence of a hot intrusion can cause an increase in the maximum value of the pore-fluid velocity in the hydrothermal system. (2) The permeability of an intruded pluton is one of the sensitive parameters to control the pore-fluid flow, heat transfer and ore body formation in hydrothermal systems. (3) The maximum value of the pore-fluid velocity increases when the bottom temperature of the hydrothermal system is increased. (4) The topographically driven flow has significant effects on the pore-fluid flow, temperature distribution and precipitation pattern of minerals in hydrothermal systems. (5) The size of the computational domain may have some effects on the pore-fluid flow and heat transfer, indicating that the size of a hydrothermal system may affect the pore-fluid flow and heat transfer within the system. (C) 2003 Elsevier Science B.V. All rights reserved.

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.

Measurements of seismic attenuation and transient fluid pressure in partially saturated Berea sandstone: Evidence of fluid flow on the mesoscopic scale

A novel laboratory technique is proposed to investigate wave-induced fluid flow on the mesoscopic scale as a mechanism for seismic attenuation in partially saturated rocks. This technique combines measurements of seismic attenuation in the frequency range from 1 to 100?Hz with measurements of transient fluid pressure as a response of a step stress applied on top of the sample. We used a Berea sandstone sample partially saturated with water. The laboratory results suggest that wave-induced fluid flow on the mesoscopic scale is dominant in partially saturated samples. A 3-D numerical model representing the sample was used to verify the experimental results. Biot's equations of consolidation were solved with the finite-element method. Wave-induced fluid flow on the mesoscopic scale was the only attenuation mechanism accounted for in the numerical solution. The numerically calculated transient fluid pressure reproduced the laboratory data. Moreover, the numerically calculated attenuation, superposed to the frequency-independent matrix anelasticity, reproduced the attenuation measured in the laboratory in the partially saturated sample. This experimental?numerical fit demonstrates that wave-induced fluid flow on the mesoscopic scale and matrix anelasticity are the dominant mechanisms for seismic attenuation in partially saturated Berea sandstone.

Condensational growth of atmospheric aerosol particles in an expanding water saturated air flow

Aerosol particles and water vapour are two important constituents of the atmosphere. Their interaction, i.e. thecondensation of water vapour on particles, brings about the formation of cloud, fog, and raindrops, causing the water cycle on the earth, and being responsible for climate changes. Understanding the roles of water vapour and aerosol particles in this interaction has become an essential part of understanding the atmosphere. In this work, the heterogeneous nucleation on pre-existing aerosol particles by the condensation of water vapour in theflow of a capillary nozzle was investigated. Theoretical and numerical modelling as well as experiments on thiscondensation process were included. Based on reasonable results from the theoretical and numerical modelling, an idea of designing a new nozzle condensation nucleus counter (Nozzle-CNC), that is to utilise the capillary nozzle to create an expanding water saturated air flow, was then put forward and various experiments were carried out with this Nozzle-CNC under different experimental conditions. Firstly, the air stream in the long capillary nozzle with inner diameter of 1.0~mm was modelled as a steady, compressible and heat-conducting turbulence flow by CFX-FLOW3D computational program. An adiabatic and isentropic cooling in the nozzle was found. A supersaturation in the nozzle can be created if the inlet flow is water saturated, and its value depends principally on flow velocity or flow rate through the nozzle. Secondly, a particle condensational growth model in air stream was developed. An extended Mason's diffusion growthequation with size correction for particles beyond the continuum regime and with the correction for a certain particle Reynolds number in an accelerating state was given. The modelling results show the rapid condensational growth of aerosol particles, especially for fine size particles, in the nozzle stream, which, on the one hand, may induce evident `over-sizing' and `over-numbering' effects in aerosol measurements as nozzle designs are widely employed for producing accelerating and focused aerosol beams in aerosol instruments like optical particle counter (OPC) and aerodynamical particle sizer (APS). It can, on the other hand, be applied in constructing the Nozzle-CNC. Thirdly, based on the optimisation of theoretical and numerical results, the new Nozzle-CNC was built. Under various experimental conditions such as flow rate, ambient temperature, and the fraction of aerosol in the total flow, experiments with this instrument were carried out. An interesting exponential relation between the saturation in the nozzle and the number concentration of atmospheric nuclei, including hygroscopic nuclei (HN), cloud condensation nuclei (CCN), and traditionally measured atmospheric condensation nuclei (CN), was found. This relation differs from the relation for the number concentration of CCN obtained by other researchers. The minimum detectable size of this Nozzle-CNC is 0.04?m. Although further improvements are still needed, this Nozzle-CNC, in comparison with other CNCs, has severaladvantages such as no condensation delay as particles larger than the critical size grow simultaneously, low diffusion losses of particles, little water condensation at the inner wall of the instrument, and adjustable saturation --- therefore the wide counting region, as well as no calibration compared to non-water condensation substances.

Viscous-Flow Approach to In Situ Infiltration and In Vitro Saturated Hydraulic Conductivity Determination

Infiltration is dominantly gravity driven, and a viscous-flow approach was developed. Laminar film flow equilibrates gravity with the viscous force and a constant flow velocity evolves during a period lasting 3/2 times the duration of a constant input rate, qS. Film thickness F and the specific contact area L of the film per unit soil volume are the key parameters. Sprinkler irrigation produced in situ time series of volumetric water contents, θ(z,t), as determined with TDR probes. The wetting front velocity v and the time series of the mobile water content, w(z,t) were deduced from θ(z,t). In vitro steady flow in a core of saturated soil provided volume flux density, q(z,t), and flow velocity, v, as determined from a heat front velocity. The F and L parameters of the in situ and the in vitro experiments were compared. The macropore-flow restriction states that, for a particular permeable medium, the specific contact area L must be independent from qS i.e., dL/dqS = 0. If true, then the relationship of qS ∝ v3/2 could scale a wide range of input rates 0 ≤ qS ≤ saturated hydraulic conductivity, Ksat, into a permeable medium, and kinematic-wave theory would become a versatile tool to deal with non-equilibrium flow. The viscous-flow approach is based on hydromechanical principles similar to Darcy’s law, but currently it is not suited to deduce flow properties from specified individual spatial structures of permeable media.