Analytical and Numerical Methods for Nonlinear Fluid Flow Problems in Porous Media

Wenchao Liu,Jun Yao,Weiyao Zhu

DOI: https://doi.org/10.1007/978-981-97-1635-7

2024-01-01

Abstract:This book investigates in detail the mathematical methods and computation methods in efficient solution of some open nonlinear seepage flow problems

What problem does this paper attempt to address?

DEVELOPMENT OF NUMERICAL SIMULATION METHODS FOR FLUID FLOW IN POROUS MEDIA

LI Heng,ZHANG Xi-wen,HE Feng

DOI: https://doi.org/10.3321/j.issn:1000-5870.2000.05.037

2000-01-01

Abstract:The numerical simulation methods for fluid flow in porous media are introduced. Three methods are compared and the advantages and disadvantages of them are also discussed. The numerical simulation methods for porousfluid coupled flow are also reviewed. The porous media model method can be used to simulate the fluid flow in porous materials, especially for the flow in porous-fluid coupled areas.
Analysis and Simulation of Fluid-Particles Interaction with Particle Flow Code

LUO Yong,GONG Xiao-nan,WU Rui-qian

DOI: https://doi.org/10.3785/j.issn.1008-973x.2007.11.032

2007-01-01

Abstract:Based on Darcy law and Navier-Stokes,Ergun and Wen Yu equations,a few of formulas were derived and applied to evaluate the pressure gradient under different porosity and Reynolds numbers.It was proved that viscous drag force and seepage force are the active force and the passive force.By using the fixed coarse-grid fluid scheme for fluid coupling simulation and the Fish language implemented in Particle Flow Code(PFC),codes of a free settling simulation and seepage simulation under different pressures were written,and the numerical results acted in accord with the Darcy law approximately.
A semi-analytical solution for the transient response of one-dimensional unsaturated single-layer poroviscoelastic media

Zhendong Shan,Daosheng Ling,Haojiang Ding,Yungang Wang

DOI: https://doi.org/10.1016/j.compgeo.2016.01.005

IF: 5.218

2016-01-01

Computers and Geotechnics

Abstract:This paper develops a semi-analytical solution for the transient response of an unsaturated single-layer poroviscoelastic medium with two immiscible fluids by using the Laplace transformation and the state-space method. Using the elastic–viscoelastic correspondence principle, we first introduce the Kelvin–Voigt model into Zienkiewicz’s unsaturated poroelastic model. The vibrational response for unsaturated porous material can be obtained by combining these two models and assuming that the wetting and non-wetting fluids are compressible, the solid skeleton and solid particles are viscoelastic, and the inertial and mechanical couplings are taken into account. The Laplace transformation and state-space method are used to solve the basic equations with the associated initial and boundary conditions, and the analytical solution in the Laplace domain is developed. To evaluate the responses in the time domain, Durbin’s numerical inverse Laplace transform method is used to obtain the semi-analytical solution. There are three compressional waves in porous media with two immiscible fluids. Moreover, to observe the three compressional waves clearly, we assume the two immiscible fluids are water and oil. Finally, several examples are provided to show the validity of the semi-analytical solution and to assess the influences of the viscosity coefficients and dynamic permeability coefficients on the behavior of the three compressional waves.
A finite element analysis of the influence of effective stress principle on coupled problems of flow and deformation in porous media

ZX Ding,XN Gong,Z Wang,JL Yu,W Wei

2004-01-01

Abstract:Many different versions of effective stress principle for porous media have been proposed, but none has been acknowledged as being valid in any case. It's reasonable to introduce an abstract comprehensive coefficient tensor of effective stress into Terzaghi's effective stress principle for a variety of porous media. Based on the general effective stress principle, governing equations are given for the coupled mechanism of flow and deformation in fluid saturated porous media in the case of incompressibility, and the corresponding FEM equations are formulated in this paper. It depends on the specific value of the effective stress coefficient tensor whether the system coefficient matrix of the FEM equations is symmetric. A parametric study, of the influence of the effective stress coefficient within the range of possibility, is presented. The results show that the effective stress coefficient tensor has a significant effect on the coupled behavior of porous media. The presented FEM formulations have good prospects in the application to coupled problems of flow and deformation for various porous media.
A Semi-Analytical Solution For One-Dimensional Transient Response Of Single-Layer Unsaturated Porous Media

Daosheng Ling,Yun Wang,Zhendong Shan,Haojiang Ding

DOI: https://doi.org/10.16285/j.rsm.2014.01.009

2014-01-01

Rock and soil mechanics

Abstract:Based on the theory of unsaturated porous media proposed by Zienkiewicz, considering the inertia, viscous and mechanical couplings and the compressibility of fluid and solid particles, semi-analytical solutions for one-dimensional transient response of single-layer unsaturated porous media are obtained with the example of a typical boundary condition problem. During the solution procedure, the dimensionless displacement governing equations in matrix form are derived firstly. Then the the nonhomogeneous boundary conditions are transformed into homogeneous boundary conditions and the eigenvalue problems of the governing equations are solved. After that, the variation coefficient method and the orthogonality of the characteristics functions are utilized to obtain a series of ordinary differential equations with their initial conditions. Applying the precise time-integration method, the semi-analytical solutions of the ordinary differential equations are developed. Finally, several numerical examples are provided to prove the correctness of the present solution and investigate the features of one-dimensional transient response of unsaturated single-layer porous media. This method can be extended to any other boundary conditions.
Numerical methods for multiscale transport equations and application to two-phase porous media flow

Xingye Yue,Weinan E

DOI: https://doi.org/10.1016/j.jcp.2005.05.009

IF: 4.645

2005-01-01

Journal of Computational Physics

Abstract:We discuss numerical methods for linear and nonlinear transport equations with multiscale velocity fields. These methods are themselves multiscaled in nature in the sense that they use macro and micro grids, multiscale test functions. We demonstrate the efficiency of these methods and apply them to two-phase flow in heterogeneous porous media.
Dynamics of Fluid-saturated Porous Media and Its Numerical Solution

Li Jin-jun

Abstract:The purpose of this paper is to give a review of dynamics of fluid-saturated porous media and its numerical solution.Emphasis is given to the discussions on the basic formulation of dynamics of fluid-saturated porous elastic solid,elastic wave velocity and attenuation in saturated soils and time-domain numerical solution procedures in saturated soil dynamics with application to earthquake response.Topics for further researches needed are suggested.

Geology,Engineering,Materials Science,Environmental Science
ANALYSIS OF ELASTIC WAVE PROPAGATION IN SANDSTONE SATURATED BY TWO IMMISCIBLE FLUIDS

CAI Yuanqiang,LI Baozhong,XU Changjie

DOI: https://doi.org/10.3321/j.issn:1000-6915.2006.10.010

2006-01-01

Abstract:Based on the mechanics of fluids in porous media and the continuum mechanics,a continuum theory of multiphase porous media is developed to analyze the propagation of elastic waves in a porous medium saturated by a mixture of two immiscible,viscousc,ompressible fluids.The behavior of porous media can be described by three equations of motion,two continuity equations of flow,and a supplementary equation that describes the constitutive relationships.A closed-form solution is derived for a fully coupled model considering the capillary pressure and irreducible/residual saturation,which demonstrates the existence of three types of compressional waves and one type of shear wave.The theoretical results are plotted and compared with the experimental data available in the previous studies.This solution is highly versatile in that it considers compression of the solid grain,compression of both pore fluids,deformation of the porous skeleton,and spatial damping.In addition,all the material parameters in this model are measurable.
Exact Solutions for One‐dimensional Transient Response of Fluid‐saturated Porous Media

Zhendong Shan,Daosheng Ling,Haojiang Ding

DOI: https://doi.org/10.1002/nag.904

IF: 4

2011-01-01

International Journal for Numerical and Analytical Methods in Geomechanics

Abstract:Based on the Biot theory, the exact solutions for one‐dimensional transient response of single layer of fluid‐saturated porous media and semi‐infinite media are developed, in which the fluid and solid particles are assumed to be compressible and the inertial, viscous and mechanical couplings are taken into account. First, the control equations in terms of the solid displacement u and a relative displacement w are expressed in matrix form. For problems of single layer under homogeneous boundary conditions, the eigen‐values and the eigen‐functions are obtained by means of the variable separation method, and the displacement vector u is put forward using the searching method. In the case of nonhomogeneous boundary conditions, the boundary conditions are first homogenized, and the displacement field is constructed basing upon the eigen‐functions. Making use of the orthogonality of eigen‐functions, a series of ordinary differential equations with respect to dimensionless time and their corresponding initial conditions are obtained. Those differential equations are solved by the state‐space method, and the series solutions for three typical nonhomogeneous boundary conditions are developed. For semi‐infinite media, the exact solutions in integral form for two kinds of nonhomogeneous boundary conditions are presented by applying the cosine and sine transforms to the basic equations. Finally, three examples are studied to illustrate the validity of the solutions, and to assess the influence of the dynamic permeability coefficient and the fluid inertia to the transient response of porous media. Copyright © 2010 John Wiley & Sons, Ltd.
Exact Solutions of A Third Grade Fluid Flow on A Porous Plate

K. Fakhar,Zhenli Xu,Cheng Yi

DOI: https://doi.org/10.1016/j.amc.2008.02.029

IF: 4.397

2008-01-01

Applied Mathematics and Computation

Abstract:In this paper, the exact solutions for an unsteady flow of an incompressible fluid of third grade occupying space over an infinite porous plate are obtained for a given velocity component. Translational type symmetries are employed to reduce the non-linear equations governing the flow under discussion to two ordinary differential equations (with and without dispersion terms). Graphs representing the solutions are discussed, and appropriate conclusions drawn.
Boundary Element Method for Linear Problems of Fluid-Saturated Porous Media

HUANG Quan-zhang,ZHENG Xiao-ping,WANG Bin,WANG Yun,YAO Zhen-han

2011-01-01

Abstract:In this paper,a new method to solve the linear problem of fluid-saturated porous media is presented by using the boundary element method.Firstly,the whole problem can be decomposed into a series of similar sub-problems,each of which contains only one fluid inclusion.Furthermore,the relation between the fluid volume variation and pressure variation for each of the sub-problems can be obtained,and based on these relations,a set of algebraic equations dealing with the fluid pressure are constructed consequently.After solving the algebraic equations,all the fluid pressure in different pores as well as the displacement fields in the whole domain can be solved.In order to illustrate the applications of the presented method,one example of a plane strain problem is presented;and the effective elastic modulus is computed under different porosities,which shows good agreement with the results in the references.
Two Phase Compressible Flow in Porous Media

Steven H. Collicott,Daniel T. Valentine,E.L. Houghton,P.W. Carpenter

DOI: https://doi.org/10.1016/s0252-9602(11)60391-3

2011-01-01

Abstract:We study the mathematical model of two phase compressible flows through porous media. Under the condition that the compressibility of rock, oil, and water is small, we prove that the initial-boundary value problem of the nonlinear system of equations admits a weak solution.
Numerical Simulation of Thermoelastic Nonlinear Waves in Fluid Saturated Porous Media with Non-local Darcy Law

Miglena N. Koleva,Yuri Poveschenko,Lubin G. Vulkov

DOI: https://doi.org/10.1007/978-3-030-55347-0_24

2020-08-08

Abstract:In this work we present a numerical study of a non-linear space-fractional parabolic system, describing thermoelastic waves propagation in fluid saturated porous media with non-local Darcy law. We use implicit cell-centered finite difference method for numerical discretizations, combined with $$L1-2$$ formula for the fractional derivative approximation. We apply Newton’s method to compute the finite difference solution. We solve a linear system of algebraic equation in block-matrix form on each iteration, using Schur complement. Numerical experiments attest good attributes of the proposed numerical method.
Numerical analysis for generalized Forchheimer flows of slightly compressible fluids in porous media

Thinh Kieu

DOI: https://doi.org/10.48550/arXiv.1512.03951

2015-12-12

Numerical Analysis

Abstract:In this paper, we will consider the generalized Forchheimer flows for slightly compressible fluids. Using Muskat's and Ward's general form of Forchheimer equations, we describe the fluid dynamics by a nonlinear degenerate parabolic equation for density. The long time numerical approximation of the nonlinear degenerate parabolic equation with time dependent boundary conditions is studied. The stability for all positive time is established in both a continuous time scheme and a discrete backward Euler scheme. A Gronwall's inequality-type is used to study the asymptotic behavior of the solution. Error estimates for the solution are derived for both continuous and discrete time procedures. Numerical experiments confirm the theoretical analysis regarding convergence rates.
NUMERICAL AND EXPERIMENTAL STUDY OF COMPRESSIBLE GAS FLOW THROUGH A POROUS/FLUID–COUPLED AREA

Xiwen Zhang,Zhaohui Yao,Pengfei Hao

DOI: https://doi.org/10.1615/jpormedia.v21.i4.20

2018-01-01

Journal of Porous Media

Abstract:Numerical simulations and theoretical investigations of fluid flow and heat transfer in porous/fluid-coupled areas have increased greatly in recent years. Scant research has concentrated on compressible gas flow in porous/fluid-coupled areas. First, an effective experiment method is proposed that can be used to accurately measure the permeability and inertial coefficient of porous media when gas flows at a high speed and high differential pressure. Using the method, the parameters of the non-Darcy flow through polyvinyl formal porous materials were tested in an actual experiment. Additionally, modified governing equations were used to solve the non-Darcy flow in porous media and turbulent flow in compressible air-coupled regimes. A method to determine the source terms for flow in porous media is presented. A robust numerical scheme was used to discretize the equations, and time-dependent boundary conditions were used to treat the boundary conditions. A detailed numerical and experimental investigation of compressible gas flow in a straight, round pipe with porous/fluid-coupled areas and backward- and forward-facing steps is given. The computational results were in strong agreement with the experimental data.
Pore-Scale Numerical Investigations of Fluid Flow in Porous Media Using Lattice Boltzmann Method

Cheng Gao,Rui-Na Xu,Pei-Xue Jiang

DOI: https://doi.org/10.1108/hff-07-2014-0202

2015-01-01

International Journal of Numerical Methods for Heat & Fluid Flow

Abstract:Purpose – Lattice Boltzmann method (LBM) is employed to explore friction factor of single-phase fluid flow through porous media and the effects of local porous structure including geometry of grains in porous media and specific surface of porous media on two-phase flow dynamic behavior, phase distribution and relative permeability. The paper aims to discuss this issue. Design/methodology/approach – The 3D single-phase LBM model and the 2D multi-component multi-phase Shan-Chen LBM model (S-C model) are developed for fluid flow through porous media. For the solid site, the bounce back scheme is used with non-slip boundary condition. Findings – The predicted friction factor for single-phase fluid flow agrees well with experimental data and the well-known correlation. Compared with porous media with square grains, the two-phase fluids in porous media with circle grains are more connected and continuous, and consequently the relative permeability is higher. As for the factor of specific porous media surface, the relative permeability of wetting fluids varies a little in two systems with different specific surface areas. In addition, the relative permeability of non-wetting fluid decreases with the increasing of specific surface of porous media due to the large flow resistance. Originality/value – Fluid-fluid interaction and fluid-solid interaction in the SC LBM model are presented, and schemes to obtain immiscible two-phase flow and different contact angles are discussed. Two-off mechanisms acting on the wetting fluids is proposed to illustrate the relative permeability of wetting fluids varies a little in two systems with different specific surface.
Numerical analysis for two-phase flow with non-equilibrium capillary pressure in anisotropic porous media

Khaled Bouadjila,Ali Samir Saad,Mazen Saad,Wissal Mesfar

DOI: https://doi.org/10.1007/s10444-022-09972-0

2022-08-21

Advances in Computational Mathematics

Abstract:In this paper, we are concerned with the convergence analysis of a combined finite volume-non-conforming finite element scheme, to approximate the two incompressible phase flow with dynamic capillary pressure in anisotropic porous media. All diffusion terms are anisotropic and heterogeneous and they are discretized by piecewise linear non-conforming triangular finite elements whereas the mobilities are discretized on dual diamond mesh. The mobilities are approximated by a non-standard way considering the inverse of the mean value of the inverse of mobilities across the interface of the dual mesh. This approximation ensures a priori estimates on the global pressure and on the water saturation. Under the assumption of non-degeneracry of mobilities, we prove the convergence of the combined scheme for the complete two-phase flow including the non-equilibrium capillary pressure.

mathematics, applied
Nonlinearly Preconditioned Semismooth Newton Methods for Variational Inequality Solution of Two-Phase Flow in Porous Media.

Haijian Yang,Shuyu Sun,Chao Yang

DOI: https://doi.org/10.1016/j.jcp.2016.11.036

IF: 4.645

2017-01-01

Journal of Computational Physics

Abstract:Most existing methods for solving two-phase flow problems in porous media do not take the physically feasible saturation fractions between 0 and 1 into account, which often destroys the numerical accuracy and physical interpretability of the simulation. To calculate the solution without the loss of this basic requirement, we introduce a variational inequality formulation of the saturation equilibrium with a box inequality constraint, and use a conservative finite element method for the spatial discretization and a backward differentiation formula with adaptive time stepping for the temporal integration. The resulting variational inequality system at each time step is solved by using a semismooth Newton algorithm. To accelerate the Newton convergence and improve the robustness, we employ a family of adaptive nonlinear elimination methods as a nonlinear preconditioner. Some numerical results are presented to demonstrate the robustness and efficiency of the proposed algorithm. A comparison is also included to show the superiority of the proposed fully implicit approach over the classical IMplicit Pressure-Explicit Saturation (IMPES) method in terms of the time step size and the total execution time measured on a parallel computer.
An advanced algorithm for solving incompressible fluid dynamics: from Navier–Stokes to Poisson equations

Mutaz Mohammad,Alexander Trounev

DOI: https://doi.org/10.1140/epjs/s11734-024-01209-z

2024-06-26

The European Physical Journal Special Topics

Abstract:In this study, an extensive exploration of numerical methods for tackling the complexities of unsteady incompressible flow dynamics is undertaken. The investigation encompasses a range of cutting-edge techniques, including the Euler wavelets collocation method, nonlinear finite element method (FEM), linear FEM, and FEM with a projection step, all meticulously applied to the Navier–Stokes equation. The primary objective is to achieve an unprecedented level of accuracy in capturing fluid behavior, validated against existing numerical algorithms. The results of simulations reveal remarkable fidelity, particularly with the Euler wavelets collocation method, showcasing its exceptional capability to depict intricate velocity components and pressure fields. The nonlinear FEM technique is adept at unraveling phenomena such as bathtub vortex formation, while linear FEM and FEM with a projection step excel in capturing flow patterns around obstacles. Intriguingly, the study delves into blowup solutions, where the Euler wavelets collocation method shines brightest, predicting velocity components amidst transient phenomena with unparalleled precision. Ultimately, the study underscores the indispensable role of numerical methods in understanding the nuances of fluid dynamics, each method offering unique insights into the intricate dance of incompressible flows.

physics, multidisciplinary
A Novel Numerical Model for Fluid Flow in 3D Fractured Porous Media Based on an Equivalent Matrix-Fracture Network

Chi Yao,Chen He,Jianhua Yang,Qinghui Jiang,Jinsong Huang,Chuangbing Zhou

DOI: https://doi.org/10.1155/2019/9736729

IF: 2.006

2019-01-03

Geofluids

Abstract:An original 3D numerical approach for fluid flow in fractured porous media is proposed. The whole research domain is discretized by the Delaunay tetrahedron based on the concept of node saturation. Tetrahedral blocks are impermeable, and fluid only flows through the interconnected interfaces between blocks. Fractures and the porous matrix are replaced by the triangular interface network, which is the so-called equivalent matrix-fracture network (EMFN). In this way, the three-dimensional seepage problem becomes a two-dimensional problem. The finite element method is used to solve the steady-state flow problem. The big finding is that the ratio of the macroconductivity of the whole interface network to the local conductivity of an interface is linearly related to the cubic root of the number of nodes used for mesh generation. A formula is presented to describe this relationship. With this formula, we can make sure that the EMFN produces the same macroscopic hydraulic conductivity as the intact rock. The approach is applied in a series of numerical tests to demonstrate its efficiency. Effects of the hydraulic aperture of fracture and connectivity of the fracture network on the effective hydraulic conductivity of fractured rock masses are systematically investigated.

geochemistry & geophysics,geology

Analytical and Numerical Methods for Nonlinear Fluid Flow Problems in Porous Media

DEVELOPMENT OF NUMERICAL SIMULATION METHODS FOR FLUID FLOW IN POROUS MEDIA

Analysis and Simulation of Fluid-Particles Interaction with Particle Flow Code

A semi-analytical solution for the transient response of one-dimensional unsaturated single-layer poroviscoelastic media

A finite element analysis of the influence of effective stress principle on coupled problems of flow and deformation in porous media

A Semi-Analytical Solution For One-Dimensional Transient Response Of Single-Layer Unsaturated Porous Media

Numerical methods for multiscale transport equations and application to two-phase porous media flow

Dynamics of Fluid-saturated Porous Media and Its Numerical Solution

ANALYSIS OF ELASTIC WAVE PROPAGATION IN SANDSTONE SATURATED BY TWO IMMISCIBLE FLUIDS

Exact Solutions for One‐dimensional Transient Response of Fluid‐saturated Porous Media

Exact Solutions of A Third Grade Fluid Flow on A Porous Plate

Boundary Element Method for Linear Problems of Fluid-Saturated Porous Media

Two Phase Compressible Flow in Porous Media

Numerical Simulation of Thermoelastic Nonlinear Waves in Fluid Saturated Porous Media with Non-local Darcy Law

Numerical analysis for generalized Forchheimer flows of slightly compressible fluids in porous media

NUMERICAL AND EXPERIMENTAL STUDY OF COMPRESSIBLE GAS FLOW THROUGH A POROUS/FLUID–COUPLED AREA

Pore-Scale Numerical Investigations of Fluid Flow in Porous Media Using Lattice Boltzmann Method

Numerical analysis for two-phase flow with non-equilibrium capillary pressure in anisotropic porous media

Nonlinearly Preconditioned Semismooth Newton Methods for Variational Inequality Solution of Two-Phase Flow in Porous Media.

An advanced algorithm for solving incompressible fluid dynamics: from Navier–Stokes to Poisson equations

A Novel Numerical Model for Fluid Flow in 3D Fractured Porous Media Based on an Equivalent Matrix-Fracture Network