Constrained Pressure-Temperature Residual (CPTR) Preconditioner Performance for Large-Scale Thermal CO2 Injection Simulation

by   Matthias A. Cremon, et al.

This work studies the performance of a novel preconditioner, designed for thermal reservoir simulation cases and recently introduced in Roy et al. (2020) and Cremon et al. (2020), on large-scale thermal CO2 injection cases. For Carbon Capture and Sequestration (CCS) projects, injecting CO2 under supercritical conditions is typically tens of degrees colder than the reservoir temperature. Thermal effects can have a significant impact on the simulation results, but they also add many challenges for the solvers. More specifically, the usual combination of an iterative linear solver (such as GMRES) and the Constrained Pressure Residual (CPR) physics-based block-preconditioner is known to perform rather poorly or fail to converge when thermal effects play a significant role. The Constrained Pressure-Temperature Residual (CPTR) preconditioner retains the 2x2 block structure (elliptic/hyperbolic) of CPR but includes the temperature in the elliptic subsystem. The elliptic subsystem is now formed by two equations, and is dealt with by the system-solver of BoomerAMG (from the HYPRE library). Then a global smoother, ILU(0), is applied to the full system to handle the local, hyperbolic temperature fronts. We implemented CPTR in the multi-physics solver GEOS and present results on various large-scale thermal CCS simulation cases, including both Cartesian and fully unstructured meshes, up to tens of millions of degrees of freedom. The CPTR preconditioner severely reduces the number of GMRES iterations and the runtime, with cases timing out in 24h with CPR now requiring a few hours with CPTR. We present strong scaling results using hundreds of CPU cores for multiple cases, and show close to linear scaling. CPTR is also virtually insensitive to the thermal Peclet number (which compares advection and diffusion effects) and is suitable to any thermal regime.


page 5

page 6

page 7

page 8


A constrained pressure-temperature residual (CPTR) method for non-isothermal multiphase flow in porous media

In petroleum reservoir simulation, the standard preconditioner, the Cons...

Multi-Stage Preconditioners for Thermal-Compositional-Reactive Flow in Porous Media

We present a family of multi-stage preconditioners for coupled thermal-c...

A block preconditioner for non-isothermal flow in porous media

In petroleum reservoir simulation, the industry standard preconditioner,...

A Scalable Thermal Reservoir Simulator for Giant Models on Parallel Computers

This paper introduces the model, numerical methods, algorithms and paral...

Development of A Scalable Thermal Reservoir Simulator for Distributed-memory Parallel Computers

This paper introduces the model, numerical methods, algorithms and paral...

Estimating a pressure dependent thermal conductivity coefficient with applications in food technology

In this paper we introduce a method to estimate a pressure dependent the...

Parallelized Discrete Exterior Calculus for Three-Dimensional Elliptic Problems

A formulation of elliptic boundary value problems is used to develop the...

Please sign up or login with your details

Forgot password? Click here to reset