A block preconditioner for non-isothermal flow in porous media

Author: 

Roy, T
Jonsthovel, T
Lemon, C
Wathen, A

Publication Date: 

15 October 2019

Journal: 

JOURNAL OF COMPUTATIONAL PHYSICS

Last Updated: 

2019-10-07T15:59:16.767+01:00

Volume: 

395

DOI: 

10.1016/j.jcp.2019.06.038

page: 

636-652

abstract: 

© 2019 Elsevier Inc. In petroleum reservoir simulation, the industry standard preconditioner, the Constrained Pressure Residual method (CPR), is a two-stage process which involves solving a restricted pressure system with Algebraic Multigrid (AMG). Initially designed for isothermal models, this approach is often used in the thermal case. However, it does not have a specific treatment of the additional energy conservation equation and temperature variable. We seek to develop preconditioners which better capture thermal effects such as heat diffusion. In order to study the effects of both pressure and temperature on fluid and heat flow, we consider a model of non-isothermal single phase flow through porous media. For this model, we develop a block preconditioner with an efficient Schur complement approximation. Both the pressure block and the approximate Schur complement are approximately inverted using an AMG V-cycle. The resulting solver is scalable with respect to problem size and parallelization.

Symplectic id: 

1015844

Download URL: 

Submitted to ORA: 

Submitted

Publication Type: 

Journal Article