Numerical solution of the Richards equation based catchment runoff model with dd-adaptivity algorithm

2016 ◽  
Author(s):  
Michal Kuraz
2008 ◽  
Vol 10 (3) ◽  
pp. 227-244 ◽  
Author(s):  
Olaf Kolditz ◽  
Jens-Olaf Delfs ◽  
Claudius Bürger ◽  
Martin Beinhorn ◽  
Chan-Hee Park

In this paper we present an object-oriented concept for numerical simulation of multi-field problems for coupled hydrosystem analysis. Individual (flow) processes modelled by a particular partial differential equation, i.e. overland flow by the shallow water equation, variably saturated flow by the Richards equation and saturated flow by the groundwater flow equation, are identified with their corresponding hydrologic compartments such as land surface, vadose zone and aquifers, respectively. The object-oriented framework of the compartment approach allows an uncomplicated coupling of these existing flow models. After a brief outline of the underlying mathematical models we focus on the numerical modelling and coupling of overland flow, variably saturated and groundwater flows via exchange flux terms. As each process object is associated with its own spatial discretisation mesh, temporal time-stepping scheme and appropriate numerical solution procedure. Flow processes in hydrosystems are coupled via their compartment (or process domain) boundaries without giving up the computational necessities and optimisations for the numerical solution of each individual process. However, the coupling requires a bridging of different temporal and spatial scales, which is solved here by the integration of fluxes (spatially and temporally). In closing we present three application examples: a benchmark test for overland flow on an infiltrating surface and two case studies – at the Borden site in Canada and the Beerze–Reusel drainage basin in the Netherlands.


1994 ◽  
Vol 29 (1-2) ◽  
pp. 267-275 ◽  
Author(s):  
S. Herath ◽  
K. Musiake

A modelling approach is presented to simulate infiltration systems in urban areas. The model consists of a hydrological sub-model and an infiltration system sub-model. Infiltration characteristics of individual facilities are first established using steady state numerical simulation of Richards' equation. These are represented as linear relations between the facility water head and infiltration rate for given facility widths. The infiltration system model is obtained by applying continuity equation to infiltration facilities lumped over a sub-catchment. This model is then coupled to a catchment runoff model to simulate storm runoff with infiltration systems. The model is applied to an infiltration system installation in a residential area, where stormwater runoff is monitored in a pilot area and a comparative area. The observed results suggest the method is adequate to evaluate the performance of infiltration systems. Except for the catchment storage routing parameter, all model parameters are determined from physical catchment characteristics.


2011 ◽  
Vol 49 (6) ◽  
pp. 2576-2597 ◽  
Author(s):  
Heiko Berninger ◽  
Ralf Kornhuber ◽  
Oliver Sander

2021 ◽  
Vol 233 ◽  
pp. 03042
Author(s):  
Yan SU ◽  
Yan SU ◽  
Zhi-ming ZHENG ◽  
Cheng-yu GU ◽  
Long-teng ZHANG

In order to solve the characteristics of low accuracy and slow efficiency in traditional numerical solution the free surface problem, the multiquardatic radial base function collocation method(MQ RBF) is used to analyze the constant seepage and unsteady seepage of the homogeneous earth dam. Computation of transient problem of free surface of earth dam by the linear derivation of Richards equation. The results show that the calculation accuracy of the MQRBF is higher than that of the traditional numerical method. The solution process does not involve numerical integral calculation and grid reorganization, which greatly reduces the calculation amount. Compared with the Trefftz method, it has the advantage of solving boundary values and internal values at the same time. It is not limited by the solution of the Laplace equation, and its application is wider and simpler.


2004 ◽  
Vol 9 (2) ◽  
pp. 21-24 ◽  
Author(s):  
Yong-Qiang Xu ◽  
Koichi Unami ◽  
Toshihiko Kawachi ◽  
Shuhei Yoshimoto

2018 ◽  
Vol 101 ◽  
pp. 168-175 ◽  
Author(s):  
C. Chávez-Negrete ◽  
F.J. Domínguez-Mota ◽  
D. Santana-Quinteros

Author(s):  
Liyuan Qiu ◽  
Yu Zhang ◽  
Sheng Zhang ◽  
Jingwei Zhao ◽  
Tengfei Wang ◽  
...  

Abstract In urban areas, the buildings and pavements make it hard for rainwater to infiltrate into the ground. The hardened underlaying sub-crust has increased the total rainfall runoff, pushing up the peak flood flow. Drawing on the construction concept of sponge city, this paper probes deep into the materials in each layer of permeable pavement for sidewalks. Specifically, a runoff model was constructed for sidewalk pavements under rainfall conditions through numerical simulation and model testing. Using the precipitation pattern of Qingdao, China, several combinations of materials were subject to rainfall simulations, revealing how each permeable pavement controls and affects the surface runoff. The results show that the permeability of surface course and sub-crust directly bear on the starting time, peak flow, total runoff and runoff time of sub-catchment runoff; and the latter has a greater impact than the former on sub-catchment runoff.


2019 ◽  
Vol 6 (5) ◽  
Author(s):  
Yuanyuan Zha ◽  
Jinzhong Yang ◽  
Jicai Zeng ◽  
Chak‐Hau M. Tso ◽  
Wenzhi Zeng ◽  
...  

Water ◽  
2020 ◽  
Vol 12 (6) ◽  
pp. 1780
Author(s):  
Dariusz Gąsiorowski ◽  
Tomasz Kolerski

Research on seepage flow in the vadose zone has largely been driven by engineering and environmental problems affecting many fields of geotechnics, hydrology, and agricultural science. Mathematical modeling of the subsurface flow under unsaturated conditions is an essential part of water resource management and planning. In order to determine such subsurface flow, the two-dimensional (2D) Richards equation can be used. However, the computation process is often hampered by a high spatial resolution and long simulation period as well as the non-linearity of the equation. A new highly efficient and accurate method for solving the 2D Richards equation has been proposed in the paper. The developed algorithm is based on dimensional splitting, the result of which means that 1D equations can be solved more efficiently than as is the case with unsplit 2D algorithms. Moreover, such a splitting approach allows any algorithm to be used for space as well as time approximation, which in turn increases the accuracy of the numerical solution. The robustness and advantages of the proposed algorithms have been proven by two numerical tests representing typical engineering problems and performed for typical properties of soil.


Sign in / Sign up

Export Citation Format

Share Document