2021
DOI: 10.1007/s10596-021-10114-6
|View full text |Cite
|
Sign up to set email alerts
|

Application of an improved P(m)-SOR iteration method for flow in partially saturated soils

Abstract: This paper studies the potential of using the successive over-relaxation iteration method with polynomial preconditioner (P(m)-SOR) to solve variably saturated flow problems described by the linearized Richards' equation. The finite difference method is employed to numerically discretize and produce a system of linear equations. Generally, the traditional Picard method needs to re-evaluate the iterative matrix in each iteration, so it is time-consuming. And under unfavorable conditions such as infiltration int… Show more

Help me understand this report

Search citation statements

Order By: Relevance

Paper Sections

Select...
4
1

Citation Types

0
9
0

Year Published

2023
2023
2023
2023

Publication Types

Select...
5
3

Relationship

0
8

Authors

Journals

citations
Cited by 29 publications
(9 citation statements)
references
References 52 publications
0
9
0
Order By: Relevance
“…Since the granite residual soil slope is an integral structure, the seepage of rainwater on the slope is continuous. Ten, according to the principle of energy conservation, the water fow continuity equation of rainwater on the granite residual soil slope can be expressed as follows [15][16][17]:…”
Section: Teoretical Derivation Of the Water Flow Continuitymentioning
confidence: 99%
“…Since the granite residual soil slope is an integral structure, the seepage of rainwater on the slope is continuous. Ten, according to the principle of energy conservation, the water fow continuity equation of rainwater on the granite residual soil slope can be expressed as follows [15][16][17]:…”
Section: Teoretical Derivation Of the Water Flow Continuitymentioning
confidence: 99%
“…In comparison with (Freni et al, 2009), in which the Green‐Ampt model is used, our approach makes use of the even simpler Darcy's law: where in (Freni et al, 2009) the effective infiltration area is assumed to be “the horizontal area below the structure bottom where the infiltration paths start to be linear and vertical parallel”, here, from observation evidence, a significant role is played by lateral infiltration, both amongst adjacent basins and through lateral external walls of basins system. On the other hand, unlike urban infiltration trenches, designed for handling stormwater runoff, such infiltration basins generally work in saturation conditions, therefore the Darcy approach is even more so justified; furthermore, for such management practices, addressing the solution of Richards equation, with the well‐known numerical difficulties (see, for instance, Abreu et al, 2019; Berardi & Difonzo, 2022; Berardi et al, 2018; Pop et al, 2004; Zhu et al, 2022) is cumbersome: we point out that multidimensional unsaturated flow models are “rarely applied in the design of stormwater infiltration systems because of their data and expertise requirements and complexity” (Browne et al, 2013).…”
Section: Introductionmentioning
confidence: 99%
“…Numerical analysis of unsaturated flow behavior in porous media is fundamental in geotechnical engineering problems involving precipitation-induced landslides [1][2][3], dam failure of a tailing [4,5], seepage control in concrete construction [6][7][8] and CO 2 sequestration [9].…”
Section: Introductionmentioning
confidence: 99%