A scalable fully implicit framework for reservoir simulation on parallel computers

Handle URI:
http://hdl.handle.net/10754/626151
Title:
A scalable fully implicit framework for reservoir simulation on parallel computers
Authors:
Yang, Haijian; Sun, Shuyu ( 0000-0002-3078-864X ) ; Li, Yiteng; Yang, Chao
Abstract:
The modeling of multiphase fluid flow in porous medium is of interest in the field of reservoir simulation. The promising numerical methods in the literature are mostly based on the explicit or semi-implicit approach, which both have certain stability restrictions on the time step size. In this work, we introduce and study a scalable fully implicit solver for the simulation of two-phase flow in a porous medium with capillarity, gravity and compressibility, which is free from the limitations of the conventional methods. In the fully implicit framework, a mixed finite element method is applied to discretize the model equations for the spatial terms, and the implicit Backward Euler scheme with adaptive time stepping is used for the temporal integration. The resultant nonlinear system arising at each time step is solved in a monolithic way by using a Newton–Krylov type method. The corresponding linear system from the Newton iteration is large sparse, nonsymmetric and ill-conditioned, consequently posing a significant challenge to the fully implicit solver. To address this issue, the family of additive Schwarz preconditioners is taken into account to accelerate the convergence of the linear system, and thereby improves the robustness of the outer Newton method. Several test cases in one, two and three dimensions are used to validate the correctness of the scheme and examine the performance of the newly developed algorithm on parallel computers.
KAUST Department:
Physical Sciences and Engineering (PSE) Division
Citation:
Yang H, Sun S, Li Y, Yang C (2017) A scalable fully implicit framework for reservoir simulation on parallel computers. Computer Methods in Applied Mechanics and Engineering. Available: http://dx.doi.org/10.1016/j.cma.2017.10.016.
Publisher:
Elsevier BV
Journal:
Computer Methods in Applied Mechanics and Engineering
KAUST Grant Number:
BAS/1/1351-01-01
Issue Date:
10-Nov-2017
DOI:
10.1016/j.cma.2017.10.016
Type:
Article
ISSN:
0045-7825
Sponsors:
The authors would like to express their appreciations to 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 the state key laboratory program of LASG (20170062). S. Sun was also supported by KAUST through the grant BAS/1/1351-01-01.
Additional Links:
http://www.sciencedirect.com/science/article/pii/S0045782517306862
Appears in Collections:
Articles; Physical Sciences and Engineering (PSE) Division

Full metadata record

DC FieldValue Language
dc.contributor.authorYang, Haijianen
dc.contributor.authorSun, Shuyuen
dc.contributor.authorLi, Yitengen
dc.contributor.authorYang, Chaoen
dc.date.accessioned2017-11-14T12:46:05Z-
dc.date.available2017-11-14T12:46:05Z-
dc.date.issued2017-11-10en
dc.identifier.citationYang H, Sun S, Li Y, Yang C (2017) A scalable fully implicit framework for reservoir simulation on parallel computers. Computer Methods in Applied Mechanics and Engineering. Available: http://dx.doi.org/10.1016/j.cma.2017.10.016.en
dc.identifier.issn0045-7825en
dc.identifier.doi10.1016/j.cma.2017.10.016en
dc.identifier.urihttp://hdl.handle.net/10754/626151-
dc.description.abstractThe modeling of multiphase fluid flow in porous medium is of interest in the field of reservoir simulation. The promising numerical methods in the literature are mostly based on the explicit or semi-implicit approach, which both have certain stability restrictions on the time step size. In this work, we introduce and study a scalable fully implicit solver for the simulation of two-phase flow in a porous medium with capillarity, gravity and compressibility, which is free from the limitations of the conventional methods. In the fully implicit framework, a mixed finite element method is applied to discretize the model equations for the spatial terms, and the implicit Backward Euler scheme with adaptive time stepping is used for the temporal integration. The resultant nonlinear system arising at each time step is solved in a monolithic way by using a Newton–Krylov type method. The corresponding linear system from the Newton iteration is large sparse, nonsymmetric and ill-conditioned, consequently posing a significant challenge to the fully implicit solver. To address this issue, the family of additive Schwarz preconditioners is taken into account to accelerate the convergence of the linear system, and thereby improves the robustness of the outer Newton method. Several test cases in one, two and three dimensions are used to validate the correctness of the scheme and examine the performance of the newly developed algorithm on parallel computers.en
dc.description.sponsorshipThe authors would like to express their appreciations to 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 the state key laboratory program of LASG (20170062). S. Sun was also supported by KAUST through the grant BAS/1/1351-01-01.en
dc.publisherElsevier BVen
dc.relation.urlhttp://www.sciencedirect.com/science/article/pii/S0045782517306862en
dc.rightsNOTICE: this is the author’s version of a work that was accepted for publication in Computer Methods in Applied Mechanics and Engineering. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Computer Methods in Applied Mechanics and Engineering, [, , (2017-11-10)] DOI: 10.1016/j.cma.2017.10.016 . © 2017. This manuscript version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/en
dc.subjectReservoir simulationen
dc.subjectFully implicit methoden
dc.subjectSchwarz preconditioneren
dc.subjectParallel computingen
dc.titleA scalable fully implicit framework for reservoir simulation on parallel computersen
dc.typeArticleen
dc.contributor.departmentPhysical Sciences and Engineering (PSE) Divisionen
dc.identifier.journalComputer Methods in Applied Mechanics and Engineeringen
dc.eprint.versionPost-printen
dc.contributor.institutionCollege of Mathematics and Econometrics, Hunan University, Changsha, Hunan 410082, PR Chinaen
dc.contributor.institutionState Key Laboratory of Computer Science, Chinese Academy of Sciences, Beijing 100190, PR Chinaen
dc.contributor.institutionInstitute of Software, Chinese Academy of Sciences, Beijing 100190, PR Chinaen
kaust.authorSun, Shuyuen
kaust.authorLi, Yitengen
kaust.grant.numberBAS/1/1351-01-01en
All Items in KAUST are protected by copyright, with all rights reserved, unless otherwise indicated.