TY - JOUR
T1 - A Linearized Finite Difference Scheme for the Richards Equation Under Variable-Flux Boundary Conditions
AU - Fengnan, Liu
AU - Fukumoto, Yasuhide
AU - Zhao, Xiaopeng
N1 - Funding Information:
LF was supported by Fundamental Research Funds for the Central Universities (No. DUT19RC(4)038); YF was supported in part by a Grant-in-Aid for Scientific Research from the Japan Society for the Promotion of Science (Grant No. 16K05476); XZ was supported by the China Postdoctoral Science Foundation (Grant No. 2015M581689). The draft of this paper was finished when LF and XZ visited the Institute of Mathematics for Industry of Kyushu University in the summer of 2018. They wish to appreciate the hospitality of Ms. Sasaguri.
Funding Information:
LF was supported by Fundamental Research Funds for the Central Universities (No. DUT19RC(4)038); YF was supported in part by a Grant-in-Aid for Scientific Research from the Japan Society for the Promotion of Science (Grant No. 16K05476); XZ was supported by the China Postdoctoral Science Foundation (Grant No. 2015M581689). The draft of this paper was finished when LF and XZ visited the Institute of Mathematics for Industry of Kyushu University in the summer of 2018. They wish to appreciate the hospitality of Ms. Sasaguri.
Publisher Copyright:
© 2020, Springer Science+Business Media, LLC, part of Springer Nature.
PY - 2020/4/1
Y1 - 2020/4/1
N2 - The Richards equation is a degenerate nonlinear PDE that models a flow through saturated/unsaturated porous media. Research on its numerical methods has been conducted in many fields. Implicit schemes based on a backward Euler format are widely used in calculating it. However, it is difficult to obtain stability with a numerical scheme because of the strong nonlinearity and degeneracy. In this paper, we establish a linearized semi-implicit finite difference scheme that is faster than backward Euler implicit schemes. We analyze the stability of this scheme by adding a small positive perturbation ϵ to the coefficient function of the Richards equation. Moreover, we show that there is a linear relationship between the discretization error in the L∞-norm and ϵ. Numerical experiments are carried out to verify our main results.
AB - The Richards equation is a degenerate nonlinear PDE that models a flow through saturated/unsaturated porous media. Research on its numerical methods has been conducted in many fields. Implicit schemes based on a backward Euler format are widely used in calculating it. However, it is difficult to obtain stability with a numerical scheme because of the strong nonlinearity and degeneracy. In this paper, we establish a linearized semi-implicit finite difference scheme that is faster than backward Euler implicit schemes. We analyze the stability of this scheme by adding a small positive perturbation ϵ to the coefficient function of the Richards equation. Moreover, we show that there is a linear relationship between the discretization error in the L∞-norm and ϵ. Numerical experiments are carried out to verify our main results.
UR - http://www.scopus.com/inward/record.url?scp=85082663830&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85082663830&partnerID=8YFLogxK
U2 - 10.1007/s10915-020-01196-y
DO - 10.1007/s10915-020-01196-y
M3 - Article
AN - SCOPUS:85082663830
VL - 83
JO - Journal of Scientific Computing
JF - Journal of Scientific Computing
SN - 0885-7474
IS - 1
M1 - 16
ER -