TY - JOUR
T1 - Parallel reservoir simulators for fully implicit complementarity formulation of multicomponent compressible flows
AU - Yang, Haijian
AU - Sun, Shuyu
AU - Li, Yiteng
AU - Yang, Chao
N1 - KAUST Repository Item: Exported on 2020-10-01
Acknowledgements: The authors would like to express their appreciationto the anonymous reviewers for the invaluable comments that have greatly improved the quality of the manuscript. The work was supported in part by the National Natural Science Foundation of China (11571100 and 11871069).
PY - 2019/7/30
Y1 - 2019/7/30
N2 - The numerical simulation of multicomponent compressible flow in porous media is an important research topic in reservoir modeling. Traditional semi-implicit methods for such problems are usually conditionally stable, suffer from large splitting errors, and may accompany with violations of the boundedness requirement of the numerical solution. In this study we reformulate the original multicomponent equations into a nonlinear complementarity problem and discretize it using a fully implicit finite element method. We solve the resultant nonsmooth nonlinear system of equations arising at each time step by a parallel, scalable, and nonlinearly preconditioned semismooth Newton algorithm, which is able to preserve the boundedness of the solution and meanwhile treats the possibly imbalanced nonlinearity of the system. Some numerical results are presented to demonstrate the robustness and efficiency of the proposed algorithm on the Tianhe-2 supercomputer for both standard benchmarks as well as realistic problems in highly heterogeneous media.
AB - The numerical simulation of multicomponent compressible flow in porous media is an important research topic in reservoir modeling. Traditional semi-implicit methods for such problems are usually conditionally stable, suffer from large splitting errors, and may accompany with violations of the boundedness requirement of the numerical solution. In this study we reformulate the original multicomponent equations into a nonlinear complementarity problem and discretize it using a fully implicit finite element method. We solve the resultant nonsmooth nonlinear system of equations arising at each time step by a parallel, scalable, and nonlinearly preconditioned semismooth Newton algorithm, which is able to preserve the boundedness of the solution and meanwhile treats the possibly imbalanced nonlinearity of the system. Some numerical results are presented to demonstrate the robustness and efficiency of the proposed algorithm on the Tianhe-2 supercomputer for both standard benchmarks as well as realistic problems in highly heterogeneous media.
UR - http://hdl.handle.net/10754/656717
UR - https://linkinghub.elsevier.com/retrieve/pii/S001046551930222X
UR - http://www.scopus.com/inward/record.url?scp=85071355472&partnerID=8YFLogxK
U2 - 10.1016/j.cpc.2019.07.011
DO - 10.1016/j.cpc.2019.07.011
M3 - Article
SN - 0010-4655
VL - 244
SP - 2
EP - 12
JO - Computer Physics Communications
JF - Computer Physics Communications
ER -