InterPore2018 New Orleans

A posteriori error estimates, stopping criteria, and adaptivity for a two phase flow with exchange between phases as a nonlinear complementarity problem in porous media

17 May 2018, 10:47

15m

New Orleans

Oral 20 Minutes GS 2: Computational challenges in porous media simulation Parallel 10-G

Speaker

Mr Jad Dabaghi (Inria Paris)

Description

In this work we develop an a posteriori-steered algorithm for a two phase compositional flow with exchange of components between the phases in porous media. The discretization of our model is based on a backward Euler scheme in time and a finite volume scheme in space. The phase transition is treated introducing a formulation based on Henry's law. The resulting nonlinear system is solved via an inexact semi-smooth Newton method. The key ingredient for the a posteriori analysis are the discretization, linearization, and algebraic flux reconstructions allowing to devise estimators for each error component. These enable to formulate criteria for stopping the iterative algebraic solver and the iterative linearization solver whenever the corresponding error components do not affect significantly the overall error. Numerical experiments are performed using the semi-smooth Newton-min algorithm as well as the Fischer--Burmeister algorithm and the GMRES iterative linear solver to show the efficiency of the method.