Surface-Subsurface Flow Modeling
Shallow water equations
Surface flow can be modeled using the shallow water equations. They are used to model rivers, floods, flow in channels; in dam break problems; in oceanic modeling for tsunamis, tides etc.

Simulation created by the author of this article using Amanzi and Blender
I have worked on the shallow water equation solver in Amanzi. Amanzi is a framework which allows users to model phenomena involving coupling of flow and reactive transport.

Here are some aspects of the numerical implementation of shallow water equations which I have worked on:
- Well-balancedness: a well-balanced numerical scheme is that which preserves the stationary steady-state i.e. lake at rest solutions
- Extension to polygonal meshes: achieved using piecewise linear interpolation of bathymetry after it is defined on cell nodes.
- Central upwind flux: more accurate than the Rusanov flux.
- Higher order time stepping.
- Depth positivity: achieved using reconstruction of total water depth in partially wet/ dry cells.

Simulation created by the author of this article using Amanzi and Blender
- Kurganov, Finite-volume schemes for shallow-water equations, 2018
- Beljadid, Mohammadian, Kurganov, Well- balanced positivity preserving cell-vertex central-upwind scheme for shallow water flows, 2016
Surface-subsurface coupling
The coupling between the shallow water equations and the Richards equation can be achieved through subcycling.

Through subcycling the different time scales at which the two flows occur in most practical situations is exploited. Essentially, the more physically transient nature of surface flow compared to subsurface flow is taken into account.
Example:
An example showing 1D dam break over a fully saturated subsurface is shown below. Also shown is the pressure and flux distribution in the subsurface.

Simulation created by the author of this article using Amanzi