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

Coon, Ethan T.; Moulton, J. David; Kikinzon, Evgeny; Berndt, Markus; Manzini, Gianmarco; Garimella, Rao V.; Lipnikov, Konstantin; Painter, S. L.

doi:10.1016/j.advwatres.2020.103701

Cited by 31 publications

(23 citation statements)

References 52 publications

Supporting

Mentioning

Contrasting

Unclassified

Order By: Relevance

“…ATS is an integrated, distributed hydrologic code that solves the diffusion wave approximation of the St-Venant equations for surface flow coupled to the Richards equation for flow in variably saturated porous media in the subsurface (Coon et al, 2019(Coon et al, , 2020. The Richards equation is described as…”

Section: Ats Model Setupmentioning

confidence: 99%

The effects of spatial and temporal resolution of gridded meteorological forcing on watershed hydrological responses

Shuai

Mital

Coon

et al. 2022

Hydrol. Earth Syst. Sci.

Self Cite

View full text Add to dashboard Cite

Abstract. Meteorological forcing plays a critical role in accurately simulating the watershed hydrological cycle. With the advancement of high-performance computing and the development of integrated watershed models, simulating the watershed hydrological cycle at high temporal (hourly to daily) and spatial resolution (tens of meters) has become efficient and computationally affordable. These hyperresolution watershed models require high resolution of meteorological forcing as model input to ensure the fidelity and accuracy of simulated responses. In this study, we utilized the Advanced Terrestrial Simulator (ATS), an integrated watershed model, to simulate surface and subsurface flow and land surface processes using unstructured meshes at the Coal Creek Watershed near Crested Butte (Colorado). We compared simulated watershed hydrologic responses including streamflow and distributed variables such as evapotranspiration, snow water equivalent (SWE), and groundwater table driven by three publicly available, gridded meteorological forcings (GMFs) – Daily Surface Weather and Climatological Summaries (Daymet), the Parameter-elevation Regressions on Independent Slopes Model (PRISM), and the North American Land Data Assimilation System (NLDAS). By comparing various spatial resolutions (ranging from 400 m to 4 km) of PRISM, the simulated streamflow only becomes marginally worse when spatial resolution of meteorological forcing is coarsened to 4 km (or 30 % of the watershed area). However, the 4 km-resolution has much worse performance than finer resolution in spatially distributed variables such as SWE. Using the temporally disaggregated PRISM, we compared models forced by different temporal resolutions (hourly to daily), and sub-daily resolution preserves the dynamic watershed responses (e.g., diurnal fluctuation of streamflow) that are absent in results forced by daily resolution. Conversely, the simulated streamflow shows better performance using daily resolution compared to that using sub-daily resolution. Our findings suggest that the choice of GMF and its spatiotemporal resolution depends on the quantity of interest and its spatial and temporal scale, which may have important implications for model calibration and watershed management decisions.

show abstract

Section: Ats Model Setupmentioning

confidence: 99%

The effects of spatial and temporal resolution of gridded meteorological forcing on watershed hydrological responses

Shuai

Mital

Coon

et al. 2022

Hydrol. Earth Syst. Sci.

Self Cite

View full text Add to dashboard Cite

show abstract

“…The (near-elliptic) PDE written in mixed form is discretized fully implicitly using a hybrid MFD scheme [13,14] involving cell-centered degrees of freedom-the cell-centered pressure, p-as well as face-centered degrees of freedom-the one-sided face flux, w, and the Lagrange multiplier, π-as described in detail in Appendix A. The key feature of the MFD scheme is that it provides a consistent and accurate discretization for complex meshes made of highly distorted cells often encountered in flow simulation in complex geological formations [14,15]. The scheme yields large saddle-point block linear systems [16,17] that cannot be easily handled by black-box preconditioners.…”

Section: Hybrid Mfd Methods For Single-phase Flowmentioning

confidence: 99%

Multigrid reduction preconditioning framework for coupled processes in porous and fractured media

Bui,

Hamon,

Castelletto

et al. 2021

Preprint

View full text Add to dashboard Cite

Many subsurface engineering applications involve tight-coupling between fluid flow, solid deformation, fracturing, and similar processes. To better understand the complex interplay of different governing equations, and therefore design efficient and safe operations, numerical simulations are widely used. Given the relatively long time-scales of interest, fully-implicit time-stepping schemes are often necessary to avoid time-step stability restrictions. A major computational bottleneck for these methods, however, is the linear solver. These systems are extremely large and ill-conditioned. Because of the wide range of processes and couplings that may be involved-e.g. formation and propagation of fractures, deformation of the solid porous medium, viscous flow of one or more fluids in the pores and fractures, complicated well sources and sinks, etc.-it is difficult to develop general-purpose but scalable linear solver frameworks. This challenge is further aggravated by the range of different discretization schemes that may be adopted, which have a direct impact on the linear system structure. To address this obstacle, we describe a flexible framework based on multigrid reduction that can produce purely algebraic preconditioners for a wide spectrum of relevant physics and discretizations. We demonstrate its broad applicability by constructing scalable preconditioners for several problems, notably: a hybrid discretization of single-phase flow, compositional multiphase flow with complex wells, and hydraulic fracturing simulations. Extension to other systems can be handled quite naturally. We demonstrate the efficiency and scalability of the resulting solvers through numerical examples of difficult, field-scale problems.

show abstract

“…An extensive overview of MFD history may be found in [29], while theoretical foundations, derivation of mimetic schemes and their applications are presented in [30]. Recently, mimetic schemes have been applied in subsurface flows simulation, including groundwater flow described by Richards equation [31,32] and multiphase flows [33].…”

Section: Mimetic Finite Difference Schemementioning

confidence: 99%

First-order continuation method for steady-state variably saturated groundwater flow modeling

Anuprienko¹

2022

Preprint

View full text Add to dashboard Cite

Recently, the nonlinearity continuation method has been used to numerically solve boundary value problems for steady-state Richards equation. The method can be considered as a predictor-corrector procedure with the simplest form which has been applied to date having a trivial, zeroth-order predictor. In this article, effect of a more sophisticated predictor technique is examined. Numerical experiments are performed with finite volume and mimetic finite difference discretizations on various problems, including reallife examples.

show abstract

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

Cited by 31 publications

References 52 publications

The effects of spatial and temporal resolution of gridded meteorological forcing on watershed hydrological responses

The effects of spatial and temporal resolution of gridded meteorological forcing on watershed hydrological responses

Multigrid reduction preconditioning framework for coupled processes in porous and fractured media

First-order continuation method for steady-state variably saturated groundwater flow modeling

Contact Info

Product

Resources

About