  The compositional model is often used to describe multicomponent multiphase
porous media flows in the petroleum industry. The fully implicit method with
strong stability and weak constraints on time-step sizes is commonly used in
the mainstream commercial reservoir simulators. In this paper, we develop an
efficient multi-stage preconditioner for the fully implicit compositional flow
simulation. The method employs an adaptive setup phase to improve the parallel
efficiency on GPUs. Furthermore, a multi-color Gauss-Seidel algorithm based on
the adjacency matrix is applied in the algebraic multigrid methods for the
pressure part. Numerical results demonstrate that the proposed algorithm
achieves good parallel speedup while yields the same convergence behavior as
the corresponding sequential version.

Numerical Analysis;Numerical Analysis