A block preconditioner for non-isothermal flow in porous media

JOURNAL OF COMPUTATIONAL PHYSICS

卷:395

Article

Roy, Thomas¹  Jonsthovel, Tom B.²  Lemon, Christopher³  Wathen, Andrew J.¹ 

[1] Univ Oxford, Math Inst, Oxford, England

[2] Schlumberger Houston Prod Technol Ctr, Houston, TX USA

[3] Schlumberger Abingdon Technol Ctr, Abingdon, Oxon, England

关键词: Preconditioning;    Iterative solvers;    Porous media;    Thermal reservoir simulation;   

DOI  :  10.1016/j.jcp.2019.06.038

来源: Elsevier

【 摘 要 】

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. (C) 2019 Elsevier Inc. All rights reserved.

【 授权许可】

Free