An efficient, high-order probabilistic collocation method on sparse grids for three-dimensional flow and solute transport in randomly heterogeneous porous media

In this study, a probabilistic collocation method (PCM) on sparse grids is used to solve stochastic equations describing flow and transport in three-dimensional, saturated, randomly heterogeneous porous media. The Karhunen–Loève decomposition is used to represent log hydraulic conductivity Y=lnKsY=lnKs. The hydraulic head h   and average pore-velocity vv are obtained by solving the continuity equation coupled with Darcy’s law with random hydraulic conductivity field. The concentration is computed by solving a stochastic advection–dispersion equation with stochastic average pore-velocity vv computed from Darcy’s law. The PCM approach is an extension of the generalized polynomial chaos (gPC) that couples gPC with probabilistic collocation. By using sparse grid points in sample space rather than standard grids based on full tensor products, the PCM approach becomes much more efficient when applied to random processes with a large number of random dimensions. Monte Carlo (MC) simulations have also been conducted to verify accuracy of the PCM approach and to demonstrate that the PCM approach is computationally more efficient than MC simulations. The numerical examples demonstrate that the PCM approach on sparse grids can efficiently simulate solute transport in randomly heterogeneous porous media with large variances.