Discontinuous Galerkin method for two-component liquid–gas porous media flows

We consider two-component (typically, water and hydrogen) compressible liquid–gas porous media flows including mass exchange between phases possibly leading to gas-phase (dis)appearance, as motivated by hydrogen production in underground repositories of radioactive waste. Following recent work by Bourgeat, Jurak, and Smaï, we formulate the governing equations in terms of liquid pressure and dissolved hydrogen density as main unknowns, leading mathematically to a nonlinear elliptic–parabolic system of partial differential equations, in which the equations degenerate when the gas phase disappears. We develop a discontinuous Galerkin method for space discretization, combined with a backward Euler scheme for time discretization and an incomplete Newton method for linearization. Numerical examples deal with gas-phase (dis)appearance, ill-prepared initial conditions, and heterogeneous problem with different rock types.