Abstract:Understanding fluid flows in heterogeneous porous media is fundamental to applied geosciences. The wide connectivity variations in the natural aquifer or oil reservoirs make the equivalent permeability have strong spatial variations. When performing the simulations for subsurface flows, the permeabilities may have strong discontinuities across the interfaces between different grid cells. Utilizing the traditional numerical schemes to simulate flows in strong heterogeneous media, the refinement ratio for the grid cell needs to increase dramatically to get an accurate result. Recently, we proposed a finite analytic numerical scheme to solve the two-dimensional fluid flows in heterogeneous porous media. With only 2×2 or 3×3 subdivisions, this scheme can provide rather accurate solutions. In this paper, we develop the finite analytic numerical method for solving the three-dimensional fluid flows in heterogeneous porous media. For the rectangular grid system, it is generally proposed that the pressure gradient in a plane normal to the edge joining different permeability regions will tend to infinite as approaching the edge according to a typical power-law solution and the tangential derivate of the pressure along the edge must be of limited value due to the pressure continuity. Consequently, the three-dimensional flow will reduce to the two-dimensional one in the neighborhood around each edge. Such quasi-two-dimensional behavior is then applied to construct a finite analytic numerical scheme. Numerical examples show that the proposed scheme can provide rather accurate solutions with only 2×2×2 or 3×3×3 subdivisions and the convergent speed is independent of the permeability heterogeneity. Due to its high calculation efficiency, the proposed scheme is utilized to test the well known LLM (Landau, Lifshitz and Matheron) conjecture, which provides keq/kG=exp⁡(16σln⁡k2) for the isotropic log-normal porous medium. The numerical results do not support this conjecture for large σln⁡k, but strongly suggest the linear relation keq/kG=1+16σln⁡k2.

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

Properties of Generalized Forchheimer Flows in Porous Media

Anisotropic flows of Forchheimer-type in porous media and their steady states

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

Finite Analytic Numerical Method for Two-Dimensional Fluid Flow in Heterogeneous Porous Media

Numerical Analysis for a Non-isothermal Incompressible Navier–Stokes–Allen–Cahn System

Numerical simulations of unsteady viscous incompressible flows using general pressure equation

Novel numerical analysis of multi-term time fractional viscoelastic non-Newtonian fluid models for simulating unsteady MHD Couette flow of a generalized Oldroyd-B fluid

Homogenization of Non-Local Navier-Stokes-Korteweg Equations for Compressible Liquid-Vapour Flow in Porous Media

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

Numerical analysis of a porous natural convection system with vorticity and viscous dissipation

Natural convection flow and heat transfer of generalized Maxwell fluid with distributed order time fractional derivatives embedded in the porous medium

Mathematical Analysis of the Poiseuille Flow of a Fluid with Temperature-Dependent Properties

Steady, quasi-one-dimensional, internal compressible flow with area change, heat addition and friction

Muskat-Leverett two-phase flow in thin cylindric porous media: Asymptotic approach

Numerical analysis of two-dimensional MHD flow and heat transfer of generalized Maxwell fluid through a rectangular pipe

Numerical approximation for hybrid‐dimensional flow and transport in fractured porous media

A Numerical Method for a Model of Two-Phase Flow in a Coupled Free Flow and Porous Media System.

Mathematical Analysis of a Diffuse Interface Model for Multi-phase Flows of Incompressible Viscous Fluids with Different Densities

Numerical Analysis of a Periodically Varying Flow in a Circular Tube Containing a Slightly Compressible Fluid

A robust conservative mixed finite element method for compressible flow on pipe networks