Robust and accurate schemes are proposed to couple subsurface and overland
flows by enforcing the continuity of the normal flux and the pressure. Richards’ equation
governing the subsurface flow is discretized using a Backward Differentiation Formula in
time and a symmetric interior penalty Discontinuous Galerkin method in space. The kinematic
wave equation governing the overland flow is discretized using a Godunov scheme.
Both schemes are individually mass conservative and can be used within coupling algorithms
that ensure overall mass conservation owing to a specific design of the interface
fluxes in the multi-step case. For field drainage problems, we also propose a method for
representing drain tubes using Signorini type conditions. Numerical results are presented
to illustrate the performances of the proposed algorithms.
Keywords: Backward Differentiation Formula, Beavers-Joseph-Saffman condition, field drainage problems,
Godunov scheme, high-order initialization, hydraulic conductivity, kinematic wave equation, mass
conservative schemes, matrix renumerotation, Richards’ equation, saturated porous media, Signorini type
conditions, subsurface flow.
