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 Methods and Comparisons for 1D and Quasi 2D Fluid Streamer Propagation Models

Numerical Methods and Comparisons for 1D and Quasi 2D Streamer Propagation Models

Numerical Simulations for the Quasi-3D Fluid Streamer Propagation Model: Methods and Applications

Inflow and Outflow Numerical Simulation Using Least-Square Moving Particle Semi-Implicit Method on GPU

Material Point Method and Smoothed Particle Hydrodynamics Simulations of Fluid Flow Problems: a Comparative Study

An Efficient and Accurate MPI Based Parallel Simulator for Streamer Discharges in Three Dimensions

A Local Discontinuous Galerkin Method for 1.5-Dimensional Streamer Discharge Simulations

Stream Function-Velocity-magnetic Induction Compact Difference Method for the 2D Steady Incompressible Full Magnetohydrodynamic Equations

2-D Discontinuous Galerkin Method for Streamer Discharge Simulations in Nitrogen

Adaptive Mesh Refinement and Embedded Boundary Method for Streamer Discharge Simulations

An efficient compact difference scheme for solving the streamfunction formulation of the incompressible Navier-Stokes equations

An Upwind Compact Difference Scheme for Solving the Streamfunction-Velocity Formulation of the Unsteady Incompressible Navier-Stokes Equation

A Comparison Study of Numerical Methods for Compressible Two-Phase Flows

Comparison Study on the Performances of Finite Volume Method and Finite Difference Method.

A physics-preserving pure streamfunction formulation and high-order compact solver with high-resolution for three-dimensional steady incompressible flows

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

Comparison Investigations Of Numerical Simulations Of Incompressible Viscous Flows By Sph And Mps

A Compact Streamfunction-Velocity Scheme for the 2-D Unsteady Incompressible Navier-Stokes Equations in Arbitrary Curvilinear Coordinates

A high-order upwind compact difference scheme for solving the streamfunction-velocity formulation of the unsteady incompressible Navier–Stokes equations

A Comparison of Two-Fluid Model, Dense Discrete Particle Model and CFD-DEM Method for Modeling Impinging Gas–solid Flows

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