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 »
- Authors:
- Publication Date:
- Research Org.:
- Oak Ridge National Laboratory (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. doi: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 = {Thu Oct 01 00:00:00 EDT 2020},
month = {Thu Oct 01 00:00:00 EDT 2020}
}
https://doi.org/10.1016/j.advwatres.2020.103701
Works referenced in this record:
Mathematical Modeling of Watershed Hydrology
journal, July 2002
- Singh, Vijay P.; Woolhiser, David A.
- Journal of Hydrologic Engineering, Vol. 7, Issue 4
Weighted Implicit Finite-Volume Model for Overland Flow
journal, September 1998
- Lal, A. M. Wasantha
- Journal of Hydraulic Engineering, Vol. 124, Issue 9
A comparison of two physics-based numerical models for simulating surface water–groundwater interactions
journal, April 2010
- Sulis, Mauro; Meyerhoff, Steven B.; Paniconi, Claudio
- Advances in Water Resources, Vol. 33, Issue 4
Hyperresolution global land surface modeling: Meeting a grand challenge for monitoring Earth's terrestrial water: OPINION
journal, May 2011
- Wood, Eric F.; Roundy, Joshua K.; Troy, Tara J.
- Water Resources Research, Vol. 47, Issue 5
Support-Operator Finite-Difference Algorithms for General Elliptic Problems
journal, April 1995
- Shashkov, Mikhail; Steinberg, Stanly
- Journal of Computational Physics, Vol. 118, Issue 1
The integrated hydrologic model intercomparison project, IH-MIP2: A second set of benchmark results to diagnose integrated hydrology and feedbacks: INTEGRATED HYDROLOGIC MODEL INTERCOMPARISON, IH-MIP2
journal, January 2017
- Kollet, Stefan; Sulis, Mauro; Maxwell, Reed M.
- Water Resources Research, Vol. 53, Issue 1
An introduction to the European Hydrological System — Systeme Hydrologique Europeen, “SHE”, 1: History and philosophy of a physically-based, distributed modelling system
journal, October 1986
- Abbott, M. B.; Bathurst, J. C.; Cunge, J. A.
- Journal of Hydrology, Vol. 87, Issue 1-2
Application of the mixed hybrid finite element approximation in a groundwater flow model: Luxury or necessity?
journal, November 1994
- Mosé, R.; Siegel, P.; Ackerer, P.
- Water Resources Research, Vol. 30, Issue 11
Improved numerical solvers for implicit coupling of subsurface and overland flow
journal, December 2014
- Osei-Kuffuor, D.; Maxwell, R. M.; Woodward, C. S.
- Advances in Water Resources, Vol. 74
A general mass-conservative numerical solution for the unsaturated flow equation
journal, July 1990
- Celia, Michael A.; Bouloutas, Efthimios T.; Zarba, Rebecca L.
- Water Resources Research, Vol. 26, Issue 7
Hydrologic conditions controlling runoff generation immediately after wildfire: HYDROLOGIC CONDITIONS CONTROLLING RUNOFF AFTER WILDFIRE
journal, March 2012
- Ebel, Brian A.; Moody, John A.; Martin, Deborah A.
- Water Resources Research, Vol. 48, Issue 3
Modeling Coupled Surface-Subsurface Flow Processes: A Review
journal, May 2008
- Furman, Alex
- Vadose Zone Journal, Vol. 7, Issue 2
A Mathematical Model for the Coupled Simulation of Surface Runoff and Infiltration
journal, January 2000
- Gandolfi, C.; Savi, F.
- Journal of Agricultural Engineering Research, Vol. 75, Issue 1
Physically based modeling in catchment hydrology at 50: Survey and outlook: SURVEY OF PHYSICALLY BASED MODELING IN CATCHMENT HYDROLOGY
journal, September 2015
- Paniconi, Claudio; Putti, Mario
- Water Resources Research, Vol. 51, Issue 9
HydroGeoSphere: A Fully Integrated, Physically Based Hydrological Model
journal, December 2011
- Brunner, Philip; Simmons, Craig T.
- Ground Water, Vol. 50, Issue 2
Modeling of Conjunctive Two-Dimensional Surface-Three-Dimensional Subsurface Flows
journal, February 2002
- Morita, Masaru; Yen, Ben Chie
- Journal of Hydraulic Engineering, Vol. 128, Issue 2
Unsteady one‐dimensional flow over a plane: Partial equilibrium and recession hydrographs
journal, April 1980
- Morris, E. M.; Woolhiser, D. A.
- Water Resources Research, Vol. 16, Issue 2
Simulating Surface and Subsurface Water Balance Changes Due to Burn Severity
journal, January 2018
- Atchley, Adam L.; Kinoshita, Alicia M.; Lopez, Sonya R.
- Vadose Zone Journal, Vol. 17, Issue 1
A continuous/discontinuous Galerkin framework for modeling coupled subsurface and surface water flow
journal, May 2008
- Dawson, Clint
- Computational Geosciences, Vol. 12, Issue 4
Global threats to human water security and river biodiversity
journal, September 2010
- Vörösmarty, C. J.; McIntyre, P. B.; Gessner, M. O.
- Nature, Vol. 467, Issue 7315
Nonlinear Krylov acceleration applied to a discrete ordinates formulation of the k-eigenvalue problem
journal, April 2013
- Calef, Matthew T.; Fichtl, Erin D.; Warsa, James S.
- Journal of Computational Physics, Vol. 238
Managing complexity in simulations of land surface and near-surface processes
journal, April 2016
- Coon, Ethan T.; David Moulton, J.; Painter, Scott L.
- Environmental Modelling & Software, Vol. 78
First‐order exchange coefficient coupling for simulating surface water–groundwater interactions: parameter sensitivity and consistency with a physics‐based approach
journal, March 2009
- Ebel, Brian A.; Mirus, Benjamin B.; Heppner, Christopher S.
- Hydrological Processes, Vol. 23, Issue 13
A framework for classifying and comparing distributed hillslope and catchment hydrologic models: DISTRIBUTED MODEL REVIEW
journal, May 2007
- Kampf, Stephanie K.; Burges, Stephen J.
- Water Resources Research, Vol. 43, Issue 5
New preconditioning strategy for Jacobian-free solvers for variably saturated flows with Richards’ equation
journal, August 2016
- Lipnikov, Konstantin; Moulton, David; Svyatskiy, Daniil
- Advances in Water Resources, Vol. 94
Hydrologic-Response simulations for the R-5 catchment with a comprehensive physics-based model
journal, April 2001
- VanderKwaak, Joel E.; Loague, Keith
- Water Resources Research, Vol. 37, Issue 4
Integrated surface–groundwater flow modeling: A free-surface overland flow boundary condition in a parallel groundwater flow model
journal, July 2006
- Kollet, Stefan J.; Maxwell, Reed M.
- Advances in Water Resources, Vol. 29, Issue 7
A Second-Order Accurate, Finite Volume-Based, Integrated Hydrologic Modeling (FIHM) Framework for Simulation of Surface and Subsurface Flow
journal, November 2009
- Kumar, Mukesh; Duffy, Christopher J.; Salvage, Karen M.
- Vadose Zone Journal, Vol. 8, Issue 4
The role of bedrock topography on subsurface storm flow: THE ROLE OF BEDROCK TOPOGRAPHY ON SUBSURFACE STORM FLOW
journal, December 2002
- Freer, Jim; McDonnell, J. J.; Beven, K. J.
- Water Resources Research, Vol. 38, Issue 12
Surface-subsurface flow modeling with path-based runoff routing, boundary condition-based coupling, and assimilation of multisource observation data: SURFACE-SUBSURFACE FLOW MODELING
journal, February 2010
- Camporese, M.; Paniconi, C.; Putti, M.
- Water Resources Research, Vol. 46, Issue 2
Surface-subsurface model intercomparison: A first set of benchmark results to diagnose integrated hydrology and feedbacks
journal, February 2014
- Maxwell, Reed M.; Putti, Mario; Meyerhoff, Steven
- Water Resources Research, Vol. 50, Issue 2
Design and Application of a Gradient-Weighted Moving Finite Element Code II: in Two Dimensions
journal, May 1998
- Carlson, Neil N.; Miller, Keith
- SIAM Journal on Scientific Computing, Vol. 19, Issue 3
A fully coupled physically-based spatially-distributed model for evaluating surface/subsurface flow
journal, April 2004
- Panday, Sorab; Huyakorn, Peter S.
- Advances in Water Resources, Vol. 27, Issue 4
An Investigation of Numerical Grid Effects in Parameter Estimation
journal, November 2006
- Zyvoloski, George A.; Vesselinov, Velimir V.
- Ground Water, Vol. 44, Issue 6 Understanding
Design and Application of a Gradient-Weighted Moving Finite Element Code I: in One Dimension
journal, May 1998
- Carlson, Neil N.; Miller, Keith
- SIAM Journal on Scientific Computing, Vol. 19, Issue 3
A scalable coupled surface–subsurface flow model
journal, August 2015
- De Maet, T.; Cornaton, F.; Hanert, E.
- Computers & Fluids, Vol. 116
A Wetland Simulation Module for the MODFLOW Ground Water Model
journal, September 1998
- Restrepo, Jorge I.; Montoya, Angela M.; Obeysekera, Jayantha
- Groundwater, Vol. 36, Issue 5
Mixed-hybrid finite element approximations of second-order elliptic boundary-value problems
journal, May 1977
- Babuska, I.; Oden, J. T.; Lee, J. K.
- Computer Methods in Applied Mechanics and Engineering, Vol. 11, Issue 2
Mimetic finite difference method
journal, January 2014
- Lipnikov, Konstantin; Manzini, Gianmarco; Shashkov, Mikhail
- Journal of Computational Physics, Vol. 257
Connections between groundwater flow and transpiration partitioning
journal, July 2016
- Maxwell, R. M.; Condon, L. E.
- Science, Vol. 353, Issue 6297
Hyper-resolution global hydrological modelling: what is next?: “Everywhere and locally relevant”
journal, December 2014
- Bierkens, Marc F. P.; Bell, Victoria A.; Burek, Peter
- Hydrological Processes, Vol. 29, Issue 2
Analysis of coupling errors in a physically-based integrated surface water–groundwater model
journal, December 2012
- Dagès, Cécile; Paniconi, Claudio; Sulis, Mauro
- Advances in Water Resources, Vol. 49
Overland Flow on an Infiltrating Surface
journal, August 1971
- Smith, Roger E.; Woolhiser, David A.
- Water Resources Research, Vol. 7, Issue 4
A generalized Richards equation for surface/subsurface flow modelling
journal, March 2009
- Weill, S.; Mouche, E.; Patin, J.
- Journal of Hydrology, Vol. 366, Issue 1-4
Newton–Krylov-multigrid solvers for large-scale, highly heterogeneous, variably saturated flow problems
journal, July 2001
- Jones, Jim E.; Woodward, Carol S.
- Advances in Water Resources, Vol. 24, Issue 7
The mimetic finite difference method for elliptic and parabolic problems with a staggered discretization of diffusion coefficient
journal, January 2016
- Lipnikov, Konstantin; Manzini, Gianmarco; Moulton, J. David
- Journal of Computational Physics, Vol. 305
Comparative Study of Coupling Approaches for Surface Water and Subsurface Interactions
journal, May 2009
- Huang, Guobiao; Yeh, Gour-Tsyh
- Journal of Hydrologic Engineering, Vol. 14, Issue 5
A Closed-form Equation for Predicting the Hydraulic Conductivity of Unsaturated Soils1
journal, January 1980
- van Genuchten, M. Th.
- Soil Science Society of America Journal, Vol. 44, Issue 5
Fill and Spill Hillslope Runoff Representation With a Richards Equation‐Based Model
journal, November 2019
- Camporese, M.; Paniconi, C.; Putti, M.
- Water Resources Research, Vol. 55, Issue 11
Relationships among some locally conservative discretization methods which handle discontinuous coefficients
journal, December 2004
- Klausen, R. A.; Russell, T. F.
- Computational Geosciences, Vol. 8, Issue 4