RETRACTED ARTICLE: Numerical analysis of multiphase flow in porous material

Recent developments in the application of x-ray micro-tomography in laboratory geomechanics have allowed all the individual grains of sand in a test sample to be seen and identified uniquely in 3D. Combining such imaging capabilities with experiments carried out “in situ” within an imaging set-up has led to the possibility of directly observing the mechanisms of deformation as they happen. The challenge has thus become extracting pertinent, quantified information from these rich time-lapse 3D images to elucidate the mechanics at play. This paper presents a new approach (ID-Track) for the quantification of individual grain kinematics (displacements and rotations) of large quantities of sand grains (tens of thousands) in a test sample undergoing loading. With ID-Track, grains are tracked between images based on some geometrical feature(s) that allow their unique identification and matching between images. This differs from Digital Image Correlation (DIC), which makes measurements by recognising patterns between images. Since ID-Track does not use the image of a grain for tracking, it is significantly faster than DIC. The technique is detailed in the paper, and is shown to be fast and simple, giving good measurements of displacements, but suffering in the measurement of rotations when compared with Discrete DIC. Subsequently, results are presented from successful applications of ID-track to triaxial tests on two quite different sands: the angular Hostun sand and the rounded Caicos Ooids. This reveals details on the performance of the technique for different grain shapes and insight into the differences in the grain-scale mechanisms occurring in these two sands as they exhibit strain localisation under triaxial loading.     

Numerical modeling of multiphase fluid flow in deforming porous media: A comparison between two- and three-phase models for seismic analysis of earth and rockfill dams

In this paper, a fully coupled numerical model is presented for the finite element analysis of the deforming porous medium interacting with the flow of two immiscible compressible wetting and non-wetting pore fluids. The governing equations involving coupled fluid flow and deformation processes in unsaturated soils are derived within the framework of the generalized Biot theory. The displacements of the solid phase, the pressure of the wetting phase and the capillary pressure are taken as the primary unknowns of the present formulation. The other variables are incorporated into the model using the experimentally determined functions that define the relationship between the hydraulic properties of the porous medium, i.e. saturation, relative permeability and capillary pressure. It is worth mentioning that the imposition of various boundary conditions is feasible notwithstanding the choice of the primary variables. The modified Pastor–Zienkiewicz generalized constitutive model is introduced into the mathematical formulation to simulate the mechanical behavior of the unsaturated soil. The accuracy of the proposed mathematical model for analyzing coupled fluid flows in porous media is verified by the resolution of several numerical examples for which previous solutions are known. Finally, the performance of the computational algorithm in modeling of large-scale porous media problems including the large elasto-plastic deformations is demonstrated through the fully coupled analysis of the failure of two earth and rockfill dams. Furthermore, the three-phase model is compared to its simplified one which simulates the unsaturated porous medium as a two-phase one with static air phase. The paper illustrates the shortcomings of the commonly used simplified approach in the context of seismic analysis of two earth and rockfill dams. It is shown that accounting the pore air as an independent phase significantly influences the unsaturated soil behavior.     

复杂水力路径下非饱和流动的数值分析

多孔介质中的流动问题,与孔隙介质的特性,含水量状态以及含水量的变化历史密切相关。基于毛细循环滞回理论模型,考虑含水量变化历史对土水特征关系的影响,在开发的U-DYSAC2有限元程序中进行了相应的数值实施。在试验给定的初边值条件下进行了非饱和渗流模拟分析,并将模拟结果与实测数据比较,表明在压力边界条件反复变化下,考虑滞回效应能获得更接近实测的结果,证实该模型在模拟各种循环变化条件下非饱和土渗流初边值问题的适用性与必要性。对入渗重分布反复变化条件下非饱和土柱流动的数值模拟表明,考虑滞回与不考虑滞回条件下,含水量、孔隙水压力和湿峰的迁移的预测在入渗后的重分布过程差异较大。考虑滞回效应时,土柱上部的脱湿速率、下部的吸湿速率比不考虑滞回时要低。从而证实了非饱和多孔介质中的土水状态依赖于含水量变化,而且强烈依赖于土体的水力路径变化。因此,循环边界条件变化下,毛细滞回效应在非饱和渗流模拟中的影响显著,必须加以考虑。     

Image analysis algorithms for estimating porous media multiphase flow variables from computed microtomography data: a validation study

Image analysis of three-dimensional microtomographic image data has become an integral component of pore scale investigations of multiphase flow through porous media. This study focuses on the validation of image analysis algorithms for identifying phases and estimating porosity, saturation, solid surface area, and interfacial area between fluid phases from gray-scale X-ray microtomographic image data. The data used in this study consisted of (1) a two-phase high precision bead pack from which porosity and solid surface area estimates were obtained and (2) three-phase cylindrical capillary tubes of three different radii, each containing an air–water interface, from which interfacial area was estimated. The image analysis algorithm employed here combines an anisotropic diffusion filter to remove noise from the original gray-scale image data, a k-means cluster analysis to obtain segmented data, and the construction of isosurfaces to estimate solid surface area and interfacial area. Our method was compared with laboratory measurements, as well as estimates obtained from a number of other image analysis algorithms presented in the literature. Porosity estimates for the two-phase bead pack were within 1.5% error of laboratory measurements and agreed well with estimates obtained using an indicator kriging segmentation algorithm. Additionally, our method estimated the solid surface area of the high precision beads within 10% of the laboratory measurements, whereas solid surface area estimates obtained from voxel counting and two-point correlation functions overestimated the surface area by 20–40%. Interfacial area estimates for the air–water menisci contained within the capillary tubes were obtained using our image analysis algorithm, and using other image analysis algorithms, including voxel counting, two-point correlation functions, and the porous media marching cubes. Our image analysis algorithm, and other algorithms based on marching cubes, resulted in errors ranging from 1% to 20% of the analytical interfacial area estimates, whereas voxel counting and two-point correlation functions overestimated the analytical interfacial area by 20–40%. In addition, the sensitivity of the image analysis algorithms on the resolution of the microtomographic image data was investigated, and the results indicated that there was little or no improvement in the comparison with laboratory estimates for the resolutions and conditions tested.     

RETRACTED ARTICLE: Enhancement of bio-sludge gravitational thickening with weak ultrasound

Gravitational thickening is the prevailing method to reduce bio-sludge volume though the process is slow and usually requires addition of polyelectrolyte(s). This paper investigated the potential benefits of sonication on enhancing the bio-sludge gravitational thickening with very low energy dose, thereby called ??weak ultrasound??. Results showed that weak sonication significantly changed the bio-sludge settleability and the main mechanism was release of the loosely bound extracellular polymeric substances. The changes were strongly influenced by power density and sonication duration. Low frequency was better but the difference was insignificant. Weak sonication (<680?kJ/kg?DS) improved the sludge gravitational thickening while high ultrasonic energy deteriorated the process. Considering both the sludge thickening efficiency and energy consumption, the optimum conditions were 0.15?W/ml, 7?s and 25?kHz. Under such conditions, the energy dose was only 155?kJ/kg?DS, much lower than literature reports and the sludge settling time shortened from the original 24 to 12?h. Weak sonication could substitute expensive polyelectrolyte coagulant for bio-sludge thickening. Combination of weak sonication and polyelectrolyte could further reduce the settling time to 6?h. The final water content of the thickened bio-sludge was not changed after sonication or polyelectrolyte addition.     

Numerical solutions for flow in porous media

A numerical approach is proposed to model the flow in porous media using homogenization theory. The proposed concept involves the analyses of micro‐true flow at pore‐level and macro‐seepage flow at macro‐level. Macro‐seepage and microscopic characteristic flow equations are first derived from the Navier–Stokes equation at low Reynolds number through a two‐scale homogenization method. This homogenization method adopts an asymptotic expansion of velocity and pressure through the micro‐structures of porous media. A slightly compressible condition is introduced to express the characteristic flow through only characteristic velocity. This characteristic flow is then numerically solved using a penalty FEM scheme. Reduced integration technique is introduced for the volumetric term to avoid mesh locking. Finally, the numerical model is examined using two sets of permeability test data on clay and one set of permeability test data on sand. The numerical predictions agree well with the experimental data if constraint water film is considered for clay and two‐dimensional cross‐connection effect is included for sand. Copyright © 2003 John Wiley & Sons, Ltd.     

RETRACTED ARTICLE: The hydrochemical framework of surface water basins in southern Ghana

Surface water resources play a crucial role in the domestic water delivery system in Ghana. In addition, sustainable food production is based on the quality and quantity of water resources available for irrigation purposes to supplement rain-fed agricultural activities in the country. The objective of this research was to determine the main controls on the hydrochemistry of surface water resources in the southern part of Ghana and assess the quality of water from these basins for irrigation activities in the area. R-mode factor and cluster analyses were applied to 625 data points from 6 river basins in southern Ghana after the data had been log transformed and standardized for homogeneity. This study finds that surface water chemistry in the south is controlled by the chemistry of silicate mineral weathering, chemistry of rainfall, fertilizers from agricultural activities in the area, as well as the weathering of carbonate minerals. A Gibb’s diagram plotted with total dissolved solids (TDS) on the vertical axis against (Na⁺ + K⁺)/(Ca²⁺ + K⁺ + Na⁺) on the horizontal axis indicates that rock weathering plays a significant role in the hydrochemistry. Activity diagrams for the CaO–Na₂O–Al₂O–SiO₂–H₂O and CaO–MgO–Al₂O₃–SiO₂–H₂O systems suggest that kaolinite is the most stable clay mineral phase in the system. In addition, an assessment of the irrigation quality of water from these basins suggests that the basins are largely low sodium—low to medium salinity basins, delivering water of acceptable quality for irrigation purposes.     

Multiscale finite-volume formulation for multiphase flow in porous media: black oil formulation of compressible,three-phase flow with gravity

Most practical reservoir simulation studies are performed using the so-called black oil model, in which the phase behavior is represented using solubilities and formation volume factors. We extend the multiscale finite-volume (MSFV) method to deal with nonlinear immiscible three-phase compressible flow in the presence of gravity and capillary forces (i.e., black oil model). Consistent with the MSFV framework, flow and transport are treated separately and differently using a sequential implicit algorithm. A multiscale operator splitting strategy is used to solve the overall mass balance (i.e., the pressure equation). The black-oil pressure equation, which is nonlinear and parabolic, is decomposed into three parts. The first is a homo geneous elliptic equation, for which the original MSFV method is used to compute the dual basis functions and the coarse-scale transmissibilities. The second equation accounts for gravity and capillary effects; the third equation accounts for mass accumulation and sources/ sinks (wells). With the basis functions of the elliptic part, the coarse-scale operator can be assembled. The gravity/capillary pressure part is made up of an elliptic part and a correction term, which is computed using solutions of gravity-driven local problems. A particular solution represents accumulation and wells. The reconstructed fine-scale pressure is used to compute the fine-scale phase fluxes, which are then used to solve the nonlinear saturation equations. For this purpose, a Schwarz iterative scheme is used on the primal coarse grid. The framework is demonstrated using challenging black-oil examples of nonlinear compressible multiphase flow in strongly heterogeneous formations.     

Numerical modeling of strongly coupled microscale multiphase flow and solid deformation

When fluid flows in porous media under subsurface conditions, significant deformation can occur. Such deformation is dependent on structural and phase characteristics. In this paper, we investigate the effect of multiphase flow on the deformation of porous media at the pore scale by implementing a strongly coupled partitioned solver discretized with finite volume (FV) technique. Specifically, the role of capillary forces on grain deformation in porous media is investigated. The fluid and solid subdomains are meshed using unstructured independent grids. The model is applied for solving multiphase coupled equations and is capable of capturing pore scale physics during primary drainage by solving the Navier-Stokes equation and advecting fluid indicator function using volume of fluid (VOF) while the fluid is interacting with a nonlinear elastic solid matrix. The convergence of the coupled solver is accelerated by Aitken underrelaxation. We also reproduce geomechanical stress conditions, at the pore scale, by applying uniaxial stress on the solid while simultaneously solving the multiphase fluid-solid interaction problem to investigate the effect of external stress on fluid occupancy, velocity-field distribution, and relative permeability. We observe that the solid matrix exhibits elasto-capillary behavior during the drainage sequence. Relative permeability endpoints are shifted on the basis of the external stress exerted.     

Solving stiff mass transfers in compositional multiphase flow models: Numerical stability and spurious solutions

This paper points out two numerical problems linked to the resolution of compositional multiphase flow models for porous media with the finite‐volume technique. In particular, we consider fluid mixtures featuring fast mass transfers between the phases, hence stiff. In this context, we show how the computation of mass exchange kinetics can be expensive and that erroneous saturation front locations arise. A numerical splitting method is developed which is proven to be stable with advection‐type time steps, whatever the stiffness of the mass transfer. The erroneous front location problem is illustrated and shown to be intrinsically linked to the numerical diffusion. If we assume that the fluids are in thermodynamical equilibrium, we find that spurious solutions can be avoided by deriving and solving a new uncoupled hyperbolic equation for the saturation. This revised version was published online in July 2006 with corrections to the Cover Date.     

Node-centered finite volume discretizations for the numerical simulation of multiphase flow in heterogeneous porous media

Three node-centered finite volume discretizations for multiphase porous media flow are presented and compared. By combination of these methods two additional discretization methods are generated. The ability of these schemes to describe flows at textural interfaces of different geologic formations is investigated. It was found that models with nonzero-entry pressures for the capillary pressure-saturation relationship in conjunction with the Box discretization may give rise to spurious oscillations for flows around low permeable lenses. Furthermore, the applicability and sensitivity of the discretization methods with regard to the used computational grids is discussed. The schemes are used for the numerical study of two-phase flow in porous media with zones of different material properties. This revised version was published online in July 2006 with corrections to the Cover Date.     

RETRACTED ARTICLE: Space and time distributions of major winter storms in the United States

Winter storms are a major weather problem in the United States and their losses have been rapidly increasing. A total of 202 catastrophic winter storms involving ice storms, blizzards, and snowstorms, each causing >$5 million in damages, occurred during 1949–2003, and their losses totaled $35.2 billion (2003 dollars). Catastrophic winter storms occurred in most parts of the contiguous United States, but were concentrated in the eastern half of the nation where 88% of all storm losses occurred. They were most frequent in the Northeast climate district (95 storms), and were least frequent in the West district (14 catastrophic storms). The annual average number of storms is 3.7 with a 1-year high of nine storms, and one year had no storms. Temporal distributions of storms and their losses exhibited considerable spatial variability across the nation. For example, when storms were very frequent in the Northeast, they were infrequent elsewhere, a result of spatial differences in storm-producing weather conditions over time. The time distribution of the nation’s 202 storms during 1949–2003 had a sizable downward trend, whereas the nation’s storm losses had a major upward trend for the 55-year period. This increase over time in losses, given the decrease in storm incidences, was a result of significant temporal increases in storm sizes and storm intensities. Increases in storm intensities were small in the northern sections of the nation, but doubled across the southern two-thirds of the nation, reflecting a climatic shift in conditions producing intense winter storms.     

Smoothed particle hydrodynamics and its applications for multiphase flow and reactive transport in porous media

Smoothed particle hydrodynamics (SPH) is a Lagrangian method based on a meshless discretization of partial differential equations. In this review, we present SPH discretization of the Navier-Stokes and advection-diffusion-reaction equations, implementation of various boundary conditions, and time integration of the SPH equations, and we discuss applications of the SPH method for modeling pore-scale multiphase flows and reactive transport in porous and fractured media.     

Numerical modelling of multiphase immiscible flow in double‐porosity featured groundwater systems

A numerical model describing the flow of multiphase, immiscible fluids in a deformable, double‐porosity featured soil has been developed. The model is focused on the modelling of the secondary porosity features in soil, which is more relevant to groundwater contamination problems. The non‐linear saturation and relative permeabilities were expressed as functions of the capillary pressures. The governing partial differential equations in terms of soil displacement and fluid pressures were solved numerically. Galerkin's weighted‐residual finite element method was employed to obtain the spatial discretization whereas temporal discretization was achieved using a fully implicit scheme. The model was verified against established, peer‐reviewed works, and the assumption that the immiscible fluids (non‐aqueous phase liquids) will flow preferentially through the secondary porosity features in soil was validated. Copyright © 2011 John Wiley & Sons, Ltd.     

A novel finite element double porosity model for multiphase flow through deformable fractured porous media

Based on the theory of double-porosity, a novel mathematical model for multiphase fluid flow in a deforming fractured reservoir is developed. The present formulation, consisting of both the equilibrium and continuity equations, accounts for the significant influence of coupling between fluid flow and solid deformation, usually ignored in the reservoir simulation literature. A Galerkin-based finite element method is applied to discretize the governing equations both in the space and time domain. Throughout the derived set of equations the solid displacements as well as the fluid pressure values are considered as the primary unknowns and may be used to determine other reservoir parameters such as stresses, saturations, etc. The final set of equations represents a highly non-linear system as the elements of the coefficient matrices are updated during each iteration in terms of the independent variables. The model is employed to solve a field scale example where the results are compared to those of ten other uncoupled models. The results illustrate a significantly different behaviour for the case of a reservoir where the impact of coupling is also considered. © 1997 by John Wiley & Sons, Ltd.     

Numerical methods for flow in a porous medium with internal boundaries

Flow of fluids and transport of solutes in porous media are subjects of wide interest in several fields of applications: reservoir engineering, subsurface hydrology, chemical engineering, etc. In this paper we will study two-phase flow in a model consisting of two different types of sediments. Here, the absolute permeability, the relative permeabilities and the capillary pressure are discontinuous functions in space. This leads to interior boundary value problems at the interface between the sediments. The saturation S_w will be discontinuous or experience large gradients at the interface. A new solution procedure for such problems will be presented. The method combines the modified method of characteristics with a weak formulation where the basis functions are discontinuous at the interior boundary. The modified method of characteristics will provide a good first approximation for the jump in the discontinuous basis functions, which leads to a fast converging iterative solution scheme for the complete problem. The method has been implemented in a two-dimensional simulator, and results from numerical experiments will be presented. This revised version was published online in August 2006 with corrections to the Cover Date.