Workshop on coupled hydrological modeling, Padova, September 23-24, 2015 LSCE Darcy multi-domain approach for coupling surface-subsurface flows: Application to benchmark problems Claude MUGLER, Emmanuel MOUCHE Laboratoire des Sciences du Climat et de l’Environnement UMR 8212 CEA/CNRS/UVSQ, Orme des Merisiers, 91191 Gif-sur-Yvette, France 1/17
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Workshop on coupled hydrological modeling, Padova, September 23-24, 2015 LSCE
Darcy multi-domain approach for coupling surface-subsurface flows:
Application to benchmark problems
Claude MUGLER, Emmanuel MOUCHE
Laboratoire des Sciences du Climat et de l’Environnement UMR 8212 CEA/CNRS/UVSQ, Orme des Merisiers,
91191 Gif-sur-Yvette, France
1/17
Workshop on coupled hydrological modeling, Padova, September 23-24, 2015 LSCE
Summary
• The integrated model: Description and validation
• Integrated Hydrologic Model Intercomparison • Conclusion
2/17
Workshop on coupled hydrological modeling, Padova, September 23-24, 2015 LSCE
– Unsaturated Zone (UZ): Richards equation
– Saturated Zone (SZ): Darcy equation
Pressure head h as the main variable è unified description of flow in the UZ and SZ
)))((.()( zhhKthhC ∇+∇∇=
!!!
∂∂
))(.( zhKthS sat ∇+∇∇=
!!!
∂∂
⎩⎨⎧
=∂
∂=
SZinSUZinhC
hhC sub
sub
)()( θ
h 0
Ksub(h)
Csub(h)
⎩⎨⎧
=SZinKUZinhK
hKsat
sub
)()(
Ksat
K(h)
C(h)
S
Subsurface Model
Le Potier, CMWR XII (1998)
3/17
Workshop on coupled hydrological modeling, Padova, September 23-24, 2015 LSCE
diffusive wave + Manning formula
))((3/5
sss
ss zhxSn
hxt
h+
∂∂
∂∂
=∂∂ hs = runoff water depth
zs = soil surface elevation n = Manning’s coefficient Ss = soil slope
Surface-subsurface coupling: Introduction of runoff
4/17
Workshop on coupled hydrological modeling, Padova, September 23-24, 2015 LSCE
diffusive wave + Manning law
))((3/5
sss
ss zhxSn
hxt
h+
∂∂
∂∂
=∂∂ hs = runoff water depth
zs = soil surface elevation n = Manning’s coefficient Ss = soil slope
Surface-subsurface coupling: Introduction of runoff
5/17
è same type of equation as Richards and Darcy equations
è Runoff modeled as Darcean flow in a porous layer
Weill, PhD thesis (2008) Weill et al., J. Hydrol. (2009)
Workshop on coupled hydrological modeling, Padova, September 23-24, 2015 LSCE
• Unified equation
)())(.()( zhHqHHKtHHC uniuniuni +==∇−∇−
!!
∂∂
Integrated model: Darcy multidomain
• Physical laws for the whole domain
⎪⎩
⎪⎨⎧
=surfacehK
subsurfacehKHK
ss
subuni )(
)()(
⎩⎨⎧
=surfacehsubsurfaceh
Hss
subuni )(
)()(
θ
θθ
hHC uni
uni ∂θ∂
=)( with
A single equation describes the whole set of surface & subsurface processes and their interactions
6/17
Workshop on coupled hydrological modeling, Padova, September 23-24, 2015 LSCE
• Resolution of a single nonlinear system with domain dependent parameters (Darcean continuum)
• Natural continuity of pressure and flux at the soil surface
• Runoff / infiltration partitioning naturally controlled by pressure at the soil surface
• Same formalism to describe runoff and streams
• Can take into account any friction law
Integrated model: Advantages of the approach
7/17
Workshop on coupled hydrological modeling, Padova, September 23-24, 2015 LSCE
● Cast3M simulation platform (www-cast3m.cea.fr)
● Spatial scheme:
- Mixed Hybrid Finite Elements - Finite Volumes
● Time scheme:
- Iterative Picard algorithm for nonlinear terms (n: time index, i: iteration index)
- Underrelaxation for nonlinear laws
))(.()( 1,11,1
1,1,1 zhK
tHHhC in
in
ninin ∇+∇∇=
Δ− ++
++
+++
!!!
)10()()1()( 1,1,11,1 <<−+= −++++ ααα ininin hKhKK
Integrated model: Numerics
8/17
Workshop on coupled hydrological modeling, Padova, September 23-24, 2015 LSCE
Abdul and Gillham (WRR, 1984)
Ogden and Watts (WRR, 2000) Govindaraju and Kavvas (WRR, 1991)
Di Giammarco et al (J Hydrol 1996)
Mugler et al ( sub. J Hydrol)
Vauclin et al (WRR, 1978)
Subsurface flow and transport
Overland flow model Integrated surface/subsurface model
outlet
saturated zone
unsaturated zone
Rainfall
prescribed head boundary
no flow boundaries
saturated length
3D configuration
Validation & Application
9/17
Weill, PhD thesis (2008); Weill et al., J. Hydrol. (2009)
Workshop on coupled hydrological modeling, Padova, September 23-24, 2015 LSCE
2nd phase of the « Integrated Hydrologic Model Intercomparison Project » Maxwell et al., WRR 2014; Kollet et al., EGU 2015; www.hpsc-terrsys.de/intercomparison-project
- Organizers: S. Kollet (Forschungszentrum Jülich GmbH), R. Maxwell (Colorado School of Mines), M. Putti (Univ. of Padova), C. Paniconi (Univ. of Québec)
- Models: CATHY, Cast3M, HydroGeoSphere, OpenGeoSys, MIKE SHE,
ParFlow, PAWS, PIHM
- Focus: - 3D surface-subsurface flow interactions - more complex heterogeneity - a field experiment
Bonn meeting, 2013
Application to benchmark problems (1/2)
10/17
Workshop on coupled hydrological modeling, Padova, September 23-24, 2015 LSCE
1 Tilted v-catchment: 3D, homogenous subsurface, recession and rain/recession
3 Borden field experiment: 3D, real topography, rain/recession experiment