A one-dimensional local discontinuous Galerkin Richards' equation solution with dual-time stepping

被引:0
|
作者
Xiao, Yilong [1 ]
Kubatko, Ethan J. [1 ]
Conroy, Colton J. [2 ]
机构
[1] Ohio State Univ, Columbus, OH 43210 USA
[2] Columbia Univ, New York, NY USA
基金
美国国家科学基金会;
关键词
Richards' equation; Dual-time stepping; Local discontinuous Galerkin; Hydrostatic pressure; INFILTRATION; MODEL; TRANSIENT;
D O I
10.1007/s10596-021-10098-3
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
We present a compact, high-order Richards' equation solver using a local discontinuous Galerkin finite element method in space and a dual-time stepping method in time. Dual-time stepping methods convert a transient problem to a steady state problem, enabling direct evaluation of residual terms and resolve implicit equations in a step-wise manner keeping the method compact and amenable to parallel computing. Verification of our solver against an analytical solution shows high-order error convergence and demonstrates the solvers ability to maintain high accuracy using low spatial resolution; the method is robust and accurately resolves numerical solutions with time steps that are much larger than what is normally required for lower-order implicit schemes. Resilience of our solver (in terms of nonlinear convergence) is demonstrated in ponded infiltration into homogeneous and layered soils, for which HYDRUS-1D solutions are used as qualitative references to gauge performance of two slope limiting schemes.
引用
收藏
页码:171 / 194
页数:24
相关论文
共 50 条