Parallel simulations of pore-scale flow through porous media

Smoothed particle hydrodynamics (SPH) is a versatile technique which can be applied to single and multiphase flow through porous media. The versatility of SPH is offset by its computational expense which limits the practicability of SPH for large problems involving low Reynolds number flow. A parallel pore-scale numerical model based on SPH is described for modeling flow phenomena in porous media. Aspects of SPH which complicate parallelization are emphasized. The speed of the method is demonstrated to be proportional to the number of processors for test cases where load balance was achieved. The parallel algorithm permits the application of SPH to more complicated porous media problems than previously considered. For such problems, best performance is achieved when several soil grains are simulated by each processor. Finally, future applications of the method and possible extensions are discussed.     

Multiphase SPH modeling of free surface flow in porous media with variable porosity

This paper presents a numerical model for simulating free surface flow in porous media with spatially varying porosity. The governing equations are based on the mixture theory. The resistance forces between solid and fluid is assumed to be nonlinear. A multiphase SPH approach is presented to solve the governing equations. In the multiphase SPH, water is modeled as a weakly compressible fluid, and solid phase is discretized by fixed solid particles carrying information of porosity. The model is validated by several numerical examples including seepage through specimen, fast flow through rockfill dam and wave interaction with porous structure. Good agreements between numerical results and experimental data are obtained in terms of flow rate and evolution of free surface. Parameter study shows that (1) the nonlinear resistance law provides more accurate results; (2) particle size and porosity have significant influence on the porous flow.     

Smooth particle hydrodynamics simulations of low Reynolds number flows through porous media

In this paper, a three‐dimensional smooth particle hydrodynamics (SPH) simulator for modeling grain scale fluid flow in porous media is presented. The versatility of the SPH method has driven its use in increasingly complex areas of flow analysis, including the characterization of flow through permeable rock for both groundwater and petroleum reservoir research. SPH provides the means to model complex multi‐phase flows through such media; however, acceptance of the methodology has been hampered by the apparent lack of actual verification within the literature, particulary in the three‐dimensional case. In this paper, the accuracy of SPH is addressed via a comparison to the previously recognized benchmarks of authors such as Sangani and Acrivos (Int. J. Multiphase Flow 1982; 8 (4): 343–360), Zick and Homsy (J. Fluid Mech. 1982; 115 :13–26) and Larson and Higdon (Phys. Fluids A 1989; 1 (1):38–46) for the well‐defined classical problems of flow through idealized two‐ and three‐dimensional porous media. The accuracy of results for such low Reynolds number flows is highly dependent on the implementation of no‐slip boundary conditions. A new, robust and numerically efficient, method for implementing such boundaries in SPH is presented. Simulation results for friction coefficient and permeability are shown to agree well with the available benchmarks. Copyright © 2010 John Wiley & Sons, Ltd.     

The numerical modelling of advective transport in the presence of fluid pressure transients

Conventional modelling of transport problems for porous media usually assumes that the Darcy flow velocities are steady. In certain practical situations, the flow velocity can exhibit time‐dependency, either due to the transient character of the flow process or time dependency in the boundary conditions associated with potential flow. In this paper, we consider certain one‐ and three‐dimensional problems of the advective transport of a chemical species in a fluid‐saturated porous region. In particular, the advective flow velocity is governed by the piezo‐conduction equation that takes into account the compressibilities of the pore fluid and the porous skeleton. Time‐ and/or mesh‐refining adaptive schemes used in the computational modelling are developed on the basis of a Fourier analysis, which can lead to accurate and optimal solutions for the advective transport problem with time‐ and space‐dependent advective flow velocity distributions. Copyright © 2006 John Wiley & Sons, Ltd.     

A 2D mixed fracture–pore seepage model and hydromechanical coupling for fractured porous media

A novel two-dimensional mixed fracture–pore seepage model for fluid flow in fractured porous media is presented based on the computational framework of finite-discrete element method (FDEM). The model consists of a porous seepage model in triangular elements bonded by unbroken joint elements, as well as a fracture seepage model in broken joint elements. The principle for determining the fluid exchange coefficient of the unbroken joint element is provided to ensure numerical accuracy and efficiency. The mixed fracture–pore seepage model provides a simple but effective tool for solving fluid flow in fractured porous media. In this paper, examples of 1D and 2D seepage flow in porous media and porous media with a single fracture or multiple fractures are studied. The simulation results of the model match well with theoretical solutions or results obtained by commercial software, which verifies the correctness of the mixed fracture–pore seepage model. Furthermore, combining FDEM mechanical calculation and the mixed fracture–pore seepage model, a coupled hydromechanical model is built to simulate fluid-driven dynamic propagation of cracks in the porous media, as well as its influence on pore seepage and fracture seepage.

     

Numerical simulation of flow processes in liquefied soils using a soil–water-coupled smoothed particle hydrodynamics method

Liquefaction can result in the damage or collapse of structures during an earthquake and can therefore be a great threat to life and property. Many site investigations of liquefaction disasters are needed to study the large-scale deformation and flow mechanisms of liquefied soils that can be used for performance assessments and infrastructure improvement. To overcome the disadvantages of traditional flow analysis methods for liquefied soils, a soil–water-coupled smoothed particle hydrodynamics (SPH) modeling method was developed to analyze flow in liquefied soils. In the proposed SPH method, water and soil were simulated as different layers, while permeability, porosity, and interaction forces could be combined to model water-saturated porous media. A simple shear test was simulated using the SPH method with an elastic model to verify its application to solid phase materials. Subsequently, the applicability of the proposed SPH modeling method to the simulation of interaction forces between water and soil was verified by a falling-head permeability test. The coupled SPH method produced good simulations for both the simple shear and falling-head permeability tests. Using a fit-for-purpose experimental apparatus, a physical flow model test of liquefied sand has been designed and conducted. To complement the physical test, a numerical simulation has been undertaken based on the soil–water-coupled SPH method. The numerical results correspond well with the physical model test results in observed configurations and velocity vectors. An embankment failure in northern Sweden was selected so that the application of the soil–water-coupled SPH method could be extended to an actual example of liquefaction. The coupled SPH method simulated the embankment failure with the site investigation well. They have also estimated horizontal displacements and velocities, which can be used to greatly improve the seismic safety of structures.     

Flow nets from finite element data

Finite element solutions for steady-state flow of fluids through porous media do not include an explicit develpoment of the flow net, and many engineers find this lack frustrating. To solve for the stream function, the analyst can use essentially the same finite element routines as he used for the potential, but the matrix corresponding to the hydraulic conductivity is not the inverse of the conductivity. Contouring the potential and the stream function is done element by element; certain anomalous situations are described in an appendix.     

多孔介质渗透系数的空间尺度效应研究进展

多孔介质渗透系数的空间尺度问题是一个与地下流体运动和溶质运移的数值模拟密切相关的应用性课题,广泛的应用需求和新的计算方法使其成为近年的热门课题之一。它涉及到相互联系的两个方面：①非均质介质场渗透系数空间尺度行为的分析与模拟;②将局部测量尺度下的试验参数转化为数值模拟网格尺度下的参数输入值的升尺度（ｕｐｓｃａｌｉｎｇ）方法和计算模型。首先介绍了该课题在概念上的拓展及其物理含义,进而以方法为主线,对这一领域具有代表性的研究成果进行了分类和评述,讨论了该课题的研究对地下水流和溶质运移的模拟分析乃至整个多孔介质流体运动研究的意义。     

Modelling Two-Phase Incompressible Flow in Porous Media Using Mixed Hybrid and Discontinuous Finite Elements

In this paper, we present a numerical model for simulating two-phase (oil–water and air–water) incompressible and immiscible flow in porous media. The mathematical model which is based on a fractional flow formulation is formed of two nonlinear partial differential equations: a mean pressure equation and a water saturation equation. These two equations can be solved in a sequential manner. Two numerical methods are used to discretize the equations of the two-phase flow model: mixed hybrid finite elements are used to treat the pressure equation, h-based Richards' equation and the diffusion term in the saturation equation, the advection term in the saturation equation is treated with the discontinuous finite elements. We propose a better way to calculate the nonlinear coefficients contained in our equations on each element of the discretized domain. In heterogeneous porous media, the saturation becomes discontinuous at the interface between two porous media. We show in this paper how to use the capillary pressure–saturation relationship in order to handle the saturation jump in the mixed hybrid finite element method. The two-phase flow simulator is verified against analytical solutions for some flow problems treated by other authors.     

Three-dimensional temporally dependent dispersion through porous media: analytical solution

A three-dimensional model for non-reactive solute transport in physically homogeneous subsurface porous media is presented. The model involves solution of the advection-dispersion equation, which additionally considered temporally dependent dispersion. The model also account for a uniform flow field, first-order decay which is inversely proportional to the dispersion coefficient and retardation factor. Porous media with semi-infinite domain is considered. Initially, the space domain is not solute free. Analytical solutions are obtained for uniform and varying pulse-type input source conditions. The governing solute transport equation is solved analytically by employing Laplace transformation technique (LTT). The solutions are illustrated and the behavior of solute transport may be observed for different values of retardation factor, for which simpler models that account for solute adsorption through a retardation factor may yield a misleading assessment of solute transport in ‘‘hydrologically sensitive’’ subsurface environments.     

Modeling of dynamic cohesive fracture propagation in porous saturated media

In this paper, a mathematical model is presented for the analysis of dynamic fracture propagation in the saturated porous media. The solid behavior incorporates a discrete cohesive fracture model, coupled with the flow in porous media through the fracture network. The double‐nodded zero‐thickness cohesive interface element is employed for the mixed mode fracture behavior in tension and contact behavior in compression. The crack is automatically detected and propagated perpendicular to the maximum effective stress. The spatial discretization is continuously updated during the crack propagation. Numerical examples from the hydraulic fracturing test and the concrete gravity dam show the capability of the model to simulate dynamic fracture propagation. The comparison is performed between the quasi‐static and fully dynamic solutions, and the performance of two analyses is investigated on the values of crack length and crack mouth opening. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

Hydrodynamic model of cells for designing systems of urban groundwater drainage

An improved mathematical hydrodynamic quasi-two-dimensional model of cells, CELSUB3, is presented for simulating drainage systems that consist of pumping well fields or subsurface drains. The CELSUB3 model is composed of an assemblage of algorithms that have been developed and tested previously and that simulate saturated flow in porous media, closed conduit flow, and flow through pumping stations. A new type of link between aquifer cells and drainage conduits is proposed. This link is verified in simple problems with well known analytical solutions. The correlation between results from analytical and mathematical solutions was considered satisfactory in all cases. To simulate more complex situations, the new proposed version, CELSUB3, was applied in a project designed to control the water-table level within a sewer system in Cha?ar Ladeado Town, Santa Fe Province, Argentina. Alternative drainage designs, which were evaluated under conditions of dynamic recharge caused by rainfall in a critical year (wettest year for the period of record) and a typical year, are briefly described. After analyzing ten alternative designs, the best technical–economic solution is a subsurface drainage system of closed conduits with pumping stations and evacuation channels. Electronic Publication     

Modelling the flow of self‐compacting concrete

A Lagrangian particle‐based method, smooth particle hydrodynamics (SPH), is used in this paper to model the flow of self‐compacting concretes (SCC) with or without short steel fibres. An incompressible SPH method is presented to simulate the flow of such non‐Newtonian fluids whose behaviour is described by a Bingham‐type model, in which the kink in the shear stress vs shear strain rate diagram is first appropriately smoothed out. The viscosity of the SCC is predicted from the measured viscosity of the paste using micromechanical models in which the second phase aggregates are treated as rigid spheres and the short steel fibres as slender rigid bodies. The basic equations solved in the SPH are the incompressible mass conservation and Navier–Stokes equations. The solution procedure uses prediction–correction fractional steps with the temporal velocity field integrated forward in time without enforcing incompressibility in the prediction step. The resulting temporal velocity field is then implicitly projected on to a divergence‐free space to satisfy incompressibility through a pressure Poisson equation derived from an approximate pressure projection. The results of the numerical simulation are benchmarked against actual slump tests carried out in the laboratory. The numerical results are in excellent agreement with test results, thus demonstrating the capability of SPH and a proper rheological model to predict SCC flow and mould‐filling behaviour. Copyright © 2010 John Wiley & Sons, Ltd.     

Experimental study on the velocity-dependent dispersion of the solute transport in different porous media

The hydrodynamic dispersion is an important factor influencing the reactive solute transport in the porous media, and many previous studies assumed that it linearly varied with the average velocity of the groundwater flow. Actually, such linear relationship has been challenged by more and more experimental observations, even in homogeneous media. In this study, we aim to investigate the relationship between hydrodynamics dispersion and the flow velocity in different types of porous media through a laboratory-controlled experiment. The results indicate that (1) the dispersion coefficient should not be a linear function of the flow velocity when the relationship between the flow velocity and the hydraulic gradient can be described by Darcy’s law satisfactorily; (2) Power function works well in describing the dispersion coefficient changing with the flow velocity for different types of porous media, and the power value is between 1.0-2.0 for different particle sizes.     

Experimental study on the velocity-dependent dispersion of the solute transport in different porous media

The hydrodynamic dispersion is an important factor influencing the reactive solute transport in the porous media, and many previous studies assumed that it linearly varied with the average velocity of the groundwater flow. Actually, such linear relationship has been challenged by more and more experimental observations, even in homogeneous media. In this study, we aim to investigate the relationship between hydrodynamics dispersion and the flow velocity in different types of porous media through a laboratory-controlled experiment. The results indicate that(1) the dispersion coefficient should not be a linear function of the flow velocity when the relationship between the flow velocity and the hydraulic gradient can be described by Darcy's law satisfactorily;(2) Power function works well in describing the dispersion coefficient changing with the flow velocity for different types of porous media, and the power value is between 1.0-2.0 for different particle sizes.     

Heterogeneous modelling behaviour at an interface in porous media

The fluid flow induced by an incident wave at a discontinuity separating two porous media is governed by the hydraulic permeabilities of both media and that of the interface. In the context of Biot’s theory, we derive the time-harmonic Green’s function for the two half-space problem allowing incident fast and slow dilatational waves to assess the heterogeneous modelling behaviour for diverse hydraulic conditions. It is found that when at least one of the media is permeated with inviscid fluids, heterogeneous modelling simulates open boundary conditions. On the other hand, when the model is saturated with viscous fluids, the modelling reproduces restrained fluid flow whose values correspond to sealed pore interface conditions, in agreement with the theoretical results. Therefore the numerical technique models correctly the wave diffusion and propagation phenomena attendant at the boundary.     

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 simulation of NAPL flow in the subsurface

A three-dimensional, three-phase numerical model is presented for simulating the movement of immiscible fluids, including nonaqueous-phase liquids (NAPLs), through porous media. The model is designed to simulate soil flume experiments and for practical application to a wide variety of contamination scenarios involving light or dense NAPLs in heterogeneous subsurface systems. The model is derived for the three-phase flow of water, NAPL, and air in porous media. The basic governing equations are based upon the mass conservation of the constitutents within the phases. The descretization chosen to transform the governing equations into the approximating equations, although logically regular, is very general. The approximating equations are a set of simultaneous coupled nonlinear equations which are solved by the Newton-Raphson method. The linear system solutions needed for the Newton-Raphson method are obtained using a matrix of preconditioner/accelerator iterative methods. Because of the special way the governing equations are implemented, the model is capable of simulating many of the phenomena considered necessary for the sucessful simulation of field problems including entry pressure phenomena, entrapment, and preferential flow paths. The model is verified by comparing it with several exact analytic test solutions and three soil flume experiments involving the introduction and movement of light nonaqueous-phase liquid (LNAPL) or dense nonaqueous-phase liquid (DNAPL) in heterogeneous sand containing a watertable. This revised version was published online in August 2006 with corrections to the Cover Date.     

Forward and Inverse Problems in Modeling of Multiphase Flow and Transport Through Porous Media

The paper deals with numerical simulation techniques for forward and inverse modelling in multiphase (multicomponent) flow through porous media. The forward simulation software system MUFTE-UG uses recent discretization techniques and fast solvers. The efficient integration of optimization strategies for the solution of the inverse problems is demonstrated in detail and also applied to practical numerical examples.