Fully Implicit Multidimensional Hybrid Upwind Scheme for Coupled Flow and Transport

Abstract

Robust and accurate fully implicit finite-volume schemes applied to Darcy-scale multiphase flow and transport in porous media are highly desirable. Recently, a smooth approximation of the saturation-dependent flux coefficients based on Implicit Hybrid Upwinding (IHU) has been proposed to improve the nonlinear convergence in fully implicit simulations with buoyancy. Here, we design a truly multidimensional extension of this approach that retains the simplicity and robustness of IHU while reducing the sensitivity of the results to the orientation of the computational (Cartesian) grid. This is achieved with the introduction of an adaptive, local coupling between the fluxes that takes the flow pattern into account. We analyze the mathematical properties of the proposed methodology to show that the scheme is monotone in the presence of competing viscous and buoyancy forces and yields saturations remaining between physical bounds. Finally, we demonstrate the efficiency and accuracy of the scheme on challenging two-dimensional two-phase examples with buoyancy, with an emphasis on the reduction of the grid orientation effect.