TY - JOUR
T1 - Scalable semismooth Newton methods with multilevel domain decomposition for subsurface flow and reactive transport in porous media
AU - Cheng, Tianpei
AU - Yang, Haijian
AU - Yang, Chao
AU - Sun, Shuyu
N1 - KAUST Repository Item: Exported on 2022-09-14
Acknowledged KAUST grant number(s): BAS/1/1351-01, URF/1/3769-01, URF/1/4074-01
Acknowledgements: The authors would like to express their appreciations to the anonymous reviewer for the invaluable comments that have greatly improved the quality of the manuscript. This work is partially supported by the National Natural Science Foundation of China (No. 12131002 and No. 11971006), the Hunan Province Natural Science Foundation of China (No. 2020JJ2002). The fourth author also greatly thanks for the support from King Abdullah University of Science and Technology (KAUST) through the grants BAS/1/1351-01, URF/1/4074-01, and URF/1/3769-01.
PY - 2022/7/14
Y1 - 2022/7/14
N2 - Large-scale modeling and predictive simulations of the subsurface flow and reactive transport system in porous media is significantly challenging, due to the high nonlinearity of the governing equations and the strong heterogeneity of material coefficients. The design of novel mathematical models and state-of-the-art methods for the flow simulation through porous media typically needs to satisfy the so-called bound-preserving property, i.e., the computed solution should stay inside a physically meaningful range. This paper presents a robust, scalable numerical framework based on the variational inequality formulation and the semismooth Newton method in a fully implicit manner, to model and simulate highly nonlinear flows without violating the boundedness requirement of the solution. Rigorous theoretical analysis for the variational inequality formulation of the problem is provided for facilitating the design of algorithms. Specifically, our approach further enhances the numerical formulation by utilizing a family of multilevel monolithic overlapping Schwarz methods for efficiently preconditioning, and the parallel implementation of the bound-preserving solver is based on the fast and robust domain decomposition technique. Numerical experiments are presented to demonstrate the efficiency and parallel scalability of the solution strategy for both standard benchmarks as well as realistic flow problems involving strong heterogeneity and high nonlinearity. We also show that the proposed framework is more robust and efficient than the commonly used inexact Newton algorithm in terms of the bound-preserving property.
AB - Large-scale modeling and predictive simulations of the subsurface flow and reactive transport system in porous media is significantly challenging, due to the high nonlinearity of the governing equations and the strong heterogeneity of material coefficients. The design of novel mathematical models and state-of-the-art methods for the flow simulation through porous media typically needs to satisfy the so-called bound-preserving property, i.e., the computed solution should stay inside a physically meaningful range. This paper presents a robust, scalable numerical framework based on the variational inequality formulation and the semismooth Newton method in a fully implicit manner, to model and simulate highly nonlinear flows without violating the boundedness requirement of the solution. Rigorous theoretical analysis for the variational inequality formulation of the problem is provided for facilitating the design of algorithms. Specifically, our approach further enhances the numerical formulation by utilizing a family of multilevel monolithic overlapping Schwarz methods for efficiently preconditioning, and the parallel implementation of the bound-preserving solver is based on the fast and robust domain decomposition technique. Numerical experiments are presented to demonstrate the efficiency and parallel scalability of the solution strategy for both standard benchmarks as well as realistic flow problems involving strong heterogeneity and high nonlinearity. We also show that the proposed framework is more robust and efficient than the commonly used inexact Newton algorithm in terms of the bound-preserving property.
UR - http://hdl.handle.net/10754/679872
UR - https://linkinghub.elsevier.com/retrieve/pii/S0021999122005022
UR - http://www.scopus.com/inward/record.url?scp=85134401679&partnerID=8YFLogxK
U2 - 10.1016/j.jcp.2022.111440
DO - 10.1016/j.jcp.2022.111440
M3 - Article
SN - 1090-2716
VL - 467
SP - 111440
JO - Journal of Computational Physics
JF - Journal of Computational Physics
ER -