ABSTRACT
The Richards equation is a nonlinear degenerate advection diffusion equation that models flow in saturated/unsaturated porous media, it's crucially important for prediction of disasters when heavy rain attacks. Efficient and precise linearized numerical schemes are necessary, but there is few study related it, and the numerical theory is incomplete because of the degeneracy and strong nonlinearity. In this paper, we establish a linearized Crank–Nicolson finite difference scheme which is a three-level scheme with almost second-order accuracy. In stability analysis, we develop a creative technique to overcome the degeneracy by adding a small positive perturbation ϵ. We also propose the error estimates by applying Young's inequality and prove the convergence order is approximate to second-order. Numerical examples are also provided to verify our main results and show the relationship between the computational error and ϵ is linear.
Acknowledgments
The draft of this paper was finished when FZ 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.
Disclosure statement
No potential conflict of interest was reported by the author(s).