期刊
ADVANCES IN WATER RESOURCES
卷 94, 期 -, 页码 11-22出版社
ELSEVIER SCI LTD
DOI: 10.1016/j.advwatres.2016.04.016
关键词
Richards' equation; Variable saturated flows; Nonlinear solver; Picards' method
资金
- National Nuclear Security Administration of the U.S. Department of Energy at Los Alamos National Laboratory [DE-AC52-06NA25396]
- US Department of Energy Office of Science Advanced Scientific Computing Research (ASCR) Program in Applied Mathematics Research
- DOE Office of Environmental Management Advanced Simulation Capability for Environmental Management (ASCEM) Program
We develop a new approach for solving the nonlinear Richards' equation arising in variably saturated flow modeling. The growing complexity of geometric models for simulation of subsurface flows leads to the necessity of using unstructured meshes and advanced discretization methods. Typically, a numerical solution is obtained by first discretizing PDEs and then solving the resulting system of nonlinear discrete equations with a Newton-Raphson-type method. Efficiency and robustness of the existing solvers rely on many factors, including an empiric quality control of intermediate iterates, complexity of the employed discretization method and a customized preconditioner. We propose and analyze a new preconditioning strategy that is based on a stable discretization of the continuum Jacobian. We will show with numerical experiments for challenging problems in subsurface hydrology that this new preconditioner improves convergence of the existing Jacobian-free solvers 3-20 times. We also show that the Picard method with this preconditioner becomes a more efficient nonlinear solver than a few widely used Jacobian-free solvers. (C) 2016 Published by Elsevier Ltd.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据