We present a compact, high-order Richards’ equation solver using a local discontinuous Galerkin finite element method in space and a dual-time stepping method in time. Dual-time stepping methods convert a transient problem to a steady state problem, enabling direct evaluation of residual terms and resolve implicit equations in a step-wise manner keeping the method compact and amenable to parallel computing. Verification of our solver against an analytical solution shows high-order error convergence and demonstrates the solvers ability to maintain high accuracy using low spatial resolution; the method is robust and accurately resolves numerical solutions with time steps that are much larger than what is normally required for lower-order implicit schemes. Resilience of our solver (in terms of nonlinear convergence) is demonstrated in ponded infiltration into homogeneous and layered soils, for which HYDRUS-1D solutions are used as qualitative references to gauge performance of two slope limiting schemes.

In this work, we measure the performance of the fixed stress split algorithm for the immiscible water-oil flow coupled with linear poromechanics. The two-phase flow equations are solved on general hexahedral elements using the multipoint flux mixed finite element method whereas the poromechanics equations are discretized using the conforming Galerkin method. We introduce a rigorous calculation of the update in poroelastic properties during the iterative solution of the coupled system equations. The effects of the coupling parameter on the performance of the fixed stress algorithm is demonstrated in two field studies: the Frio oil reservoir and the Cranfield injection site.