Abstract:This study is concerned with developing a two-dimensional multiphase model that simulates the movement of NAPL in heterogeneous aquifers. Heterogeneity is dealt with in a probabilistic sense by modeling the intrinsic permeability of the porous medium as a stochastic process. The deterministic finite element method is used to spatially discretize the multiphase flow equations. The intrinsic permeability is represented in the model via its Karhunen–Loeve expansion. This is a computationally expedient representation of stochastic processes by means of a discrete set of random variables. Further, the nodal unknowns, water phase saturations and water phase pressures, are represented by their stochastic spectral expansions. This representation involves an orthogonal basis in the space of random variables. The basis consists of orthogonal polynomial chaoses of consecutive orders. The relative permeabilities of water and oil phases, and the capillary pressure are expanded in the same manner, as well. For these variables, the set of deterministic coefficients multiplying the basis in their expansions is evaluated based on constitutive relationships expressing the relative permeabilities and the capillary pressure as functions of the water phase saturations. The implementation of the various expansions into the multiphase flow equations results in the formulation of discretized stochastic differential equations that can be solved for the deterministic coefficients appearing in the expansions representing the unknowns. This method allows the computation of the probability distribution functions of the unknowns for any point in the spatial domain of the problem at any instant in time. The spectral formulation of the stochastic finite element method used herein has received wide acceptance as a comprehensive framework for problems involving random media. This paper provides the application of this formalism to the problem of two-phase flow in a random porous medium.

Modeling groundwater flow with random hydraulic conductivity using radial basis function partition of unity method

A Surrogate Modelling Approach Based on Nonlinear Dimension Reduction for Uncertainty Quantification in Groundwater Flow Models

A Stable Radial Basis Function Partition of Unity Method with d-Rectangular Patches for Modelling Water Flow in Porous Media

Modeling of unsaturated flow through porous media using meshless methods

Localized RBF methods for modeling infiltration using the Kirchhoff-transformed Richards equation

Modeling water relative permeability in uniform grain packs

Stochastic flow simulation in 3D porous media

Solution of the 3D density-driven groundwater flow problem with uncertain porosity and permeability

Stochastic Finite Element Analysis for Multiphase Flow in Heterogeneous Porous Media

A Novel Radial Basis Function Approach for Infiltration-Induced Landslides in Unsaturated Soils

Implicit EXP-RBF techniques for modeling unsaturated flow through soils with water uptake by plant roots

Evaluation of a practical approach for field scale moisture flow modeling in heterogeneous media at a semiarid site

A hyperbolic-elliptic PDE model and conservative numerical method for gravity-dominated variably-saturated groundwater flow

Modeling non‐Fickian solute transport due to mass transfer and physical heterogeneity on arbitrary groundwater velocity fields

Global random walk solvers for fully coupled flow and transport in saturated/unsaturated porous media (extended version)

Modelling of shallow aquifers in interaction with overland water

Stochastic Study on Groundwater Flow and Solute Transport in a Porous Medium with Multi-Scale Heterogeneity

Stochastic finite element modelling of water flow in variably saturated heterogeneous soils

Estimation of 3-D hydraulic conductivity fields from fictive grain-size distributions derived from 3-D geological modeling

An efficient distribution method for nonlinear two-phase flow in highly heterogeneous multidimensional stochastic porous media