Journal
ADVANCES IN WATER RESOURCES
Volume 144, Issue -, Pages -Publisher
ELSEVIER SCI LTD
DOI: 10.1016/j.advwatres.2020.103701
Keywords
Mimetic finite differences; Integrated hydrology; Complex stratigraphy
Categories
Funding
- LANL [LDRD-20150397DR, LDRD201200068DR]
- Interoperable Design of Extreme-scale Application Software (IDEAS) project
- Next Generation Ecosystem Experiment (NGEE-Arctic) project
- Office of Biological and Environmental Research in the DOE Office of Science
- [ORNLLDRD-8872]
Ask authors/readers for more resources
We explore the coupling of surface and subsurface flows on fully unstructured meshes that conform to complex soil structures. To accommodate the distorted meshes that inevitably result from explicit representation of complex soil structures, we leverage the structure of the Mimetic Finite Difference (MFD) spatial discretization scheme to couple surface and subsurface flows. The MFD method achieves second-order accuracy and maintains local mass conservation on distorted meshes. We couple the diffusion wave approximation for surface flows to the Richards equation for subsurface flow, ensuring continuity of both pressure and flux between the surface and subsurface. The MFD method is particularly convenient for this coupling because it uses face-based constraints in the subsurface system that can be expressed as face-pressure unknowns. Those unknowns are coincident with surface cell-based unknowns, thus allowing the discrete surface system to be directly substituted into the subsurface system and solved implicitly as a global system. Robust representation of the transition between wet and dry surface conditions requires upwinding of the relative permeability and is facilitated by globalization in the nonlinear solver. The approach and its implementation in the Advanced Terrestrial Simulator (ATS) are evaluated by comparison to previously published benchmarks. Using runoff from soils with patchy groundcover (duff) as an example, we show that the new method converges significantly faster in mesh convergence tests than the commonly used two-point flux approximation.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available