skip to main content
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Coupling surface flow and subsurface flow in complex soil structures using mimetic finite differences

Abstract

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 patchymore » 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.« less

Authors:
; ; ; ; ; ; ;
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1644243
Alternate Identifier(s):
OSTI ID: 1649323
Grant/Contract Number:  
LDRD-8872; LDRD-20150397DR; LDRD201200068DR; AC05-00OR22725
Resource Type:
Published Article
Journal Name:
Advances in Water Resources
Additional Journal Information:
Journal Name: Advances in Water Resources Journal Volume: 144 Journal Issue: C; Journal ID: ISSN 0309-1708
Publisher:
Elsevier
Country of Publication:
United Kingdom
Language:
English
Subject:
54 ENVIRONMENTAL SCIENCES; Mimetic finite differences; Integrated hydrology; Complex stratigraphy

Citation Formats

Coon, Ethan T., Moulton, J. David, Kikinzon, Evgeny, Berndt, Markus, Manzini, Gianmarco, Garimella, Rao, Lipnikov, Konstantin, and Painter, Scott L. Coupling surface flow and subsurface flow in complex soil structures using mimetic finite differences. United Kingdom: N. p., 2020. Web. https://doi.org/10.1016/j.advwatres.2020.103701.
Coon, Ethan T., Moulton, J. David, Kikinzon, Evgeny, Berndt, Markus, Manzini, Gianmarco, Garimella, Rao, Lipnikov, Konstantin, & Painter, Scott L. Coupling surface flow and subsurface flow in complex soil structures using mimetic finite differences. United Kingdom. https://doi.org/10.1016/j.advwatres.2020.103701
Coon, Ethan T., Moulton, J. David, Kikinzon, Evgeny, Berndt, Markus, Manzini, Gianmarco, Garimella, Rao, Lipnikov, Konstantin, and Painter, Scott L. Thu . "Coupling surface flow and subsurface flow in complex soil structures using mimetic finite differences". United Kingdom. https://doi.org/10.1016/j.advwatres.2020.103701.
@article{osti_1644243,
title = {Coupling surface flow and subsurface flow in complex soil structures using mimetic finite differences},
author = {Coon, Ethan T. and Moulton, J. David and Kikinzon, Evgeny and Berndt, Markus and Manzini, Gianmarco and Garimella, Rao and Lipnikov, Konstantin and Painter, Scott L.},
abstractNote = {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.},
doi = {10.1016/j.advwatres.2020.103701},
journal = {Advances in Water Resources},
number = C,
volume = 144,
place = {United Kingdom},
year = {2020},
month = {10}
}

Journal Article:
Free Publicly Available Full Text

Save / Share: