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.

Simulating Self-gravitating Hydrodynamic Flows

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

2d Well-Balanced Positivity Preserving Shallow Flow Model With Wetting And Drying On Triangular Grids

A finite element implementation of the incompressible Schrödinger flow method

Three‐dimensional Numerical Modelling of Free Surface Flows with Non‐hydrostatic Pressure

Finite difference method in prolate spheroidal coordinates for freely suspended spheroidal particles in linear flows of viscous and viscoelastic fluids

A 3D Simulation of a Moving Solid in Viscous Free-Surface Flows by Coupling SPH and DEM

A fully Eulerian solver for the simulation of multiphase flows with solid bodies: application to surface gravity waves

A hybrid Eulerian-Lagrangian flow solver

3-D hydrodynamic model for free-surface flows and sediment transport on unstructured σ-coordinate grid. Part I. Numerical scheme

An Implicit Numerical Algorithm for Solving Non-Hydrostatic Free-Surface Flow Problems

A 2D Navier-Stokes Equations Solver Based on Generalized Harmonic Polynomial Cell Method With Non-Uniform Grid

Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers

A numerical strategy to discretize and solve the Poisson equation on dynamically adapted multiresolution grids for time-dependent streamer discharge simulations

Asymptotic Behavior of Incompressible Schrödinger Flow for Small Data in Three Dimensions

Simulation of Fluid–particles Flows: Heavy Particles, Flowing Regime, and Asymptotic-Preserving Schemes

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

A hybrid Cartesian-meshless method for the simulation of thermal flows with complex immersed objects

Aerodynamics Simulations of Three-Dimensional Inviscid Flow Using Curvilinear Discontinuous Galerkin Method on Unstructured Meshes

Simulating Non-Newtonian Flows with the Moving Particle Semi-Implicit Method with an Sph Kernel

Two-dimensional fluids via matrix hydrodynamics