Parallel multilevel restricted Schwarz preconditioners for implicit simulation of subsurface flows with Peng-Robinson equation of state

by   Rui Li, et al.

Parallel algorithms and simulators with good scalabilities are particularly important for large-scale reservoir simulations on modern supercomputers with a large number of processors. In this paper, we introduce and study a family of highly scalable multilevel restricted additive Schwarz (RAS) methods for the fully implicit solution of subsurface flows with Peng-Robinson equation of state in two and three dimensions. With the use of a second-order fully implicit scheme, the proposed simulator is unconditionally stable with the relaxation of the time step size by the stability condition. The investigation then focuses on the development of several types of multilevel overlapping additive Schwarz methods for the preconditioning of the resultant linear system arising from the inexact Newton iteration, and some fast solver technologies are presented for the assurance of the multilevel approach efficiency and scalability. We numerically show that the proposed fully implicit framework is highly efficient for solving both standard benchmarks as well as realistic problems with several hundreds of millions of unknowns and scalable to 8192 processors on the Tianhe-2 supercomputer.


page 14

page 15

page 16

page 17

page 19


Parallel energy-stable phase field crystal simulations based on domain decomposition methods

In this paper, we present a parallel numerical algorithm for solving the...

Energy stable schemes for gradient flows based on the DVD method

The existing discrete variational derivative method is only second-order...

A shortcut for IMEX methods: integrate the residual explicitly

In numerical time-integration with implicit-explicit (IMEX) methods, a w...

An Improved Multi-Stage Preconditioner on GPUs for Compositional Reservoir Simulation

The compositional model is often used to describe multicomponent multiph...

A staggered-grid multilevel incomplete LU for steady incompressible flows

Algorithms for studying transitions and instabilities in incompressible ...

Large scale simulation of pressure induced phase-field fracture propagation using Utopia

Non-linear phase field models are increasingly used for the simulation o...

Please sign up or login with your details

Forgot password? Click here to reset