Parallel iterative linear solvers for oil reservoir models

1986 ◽  
Vol 30 (2) ◽  
pp. 184-192 ◽  
Author(s):  
I. Efrat ◽  
M. Tismenetsky
1997 ◽  
Vol 85 (1) ◽  
pp. 145-167 ◽  
Author(s):  
Eleonora Messina ◽  
Jacques J.B de Swart ◽  
Wolter A van der Veen

2014 ◽  
Author(s):  
Seonkyoo Yoon ◽  
Zeid Alghareeb ◽  
John Williams

1995 ◽  
Vol 21 (4) ◽  
pp. 669-682 ◽  
Author(s):  
A. Asenov ◽  
D. Reid ◽  
J.R. Barker

2020 ◽  
Vol 77 (2) ◽  
pp. 597-616
Author(s):  
Siamak Akbarzadeh ◽  
Jan Hückelheim ◽  
Jens-Dominik Müller

Abstract Algorithmic differentiation (AD) is a widely-used approach to compute derivatives of numerical models. Many numerical models include an iterative process to solve non-linear systems of equations. To improve efficiency and numerical stability, AD is typically not applied to the linear solvers. Instead, the differentiated linear solver call is replaced with hand-produced derivative code that exploits the linearity of the original call. In practice, the iterative linear solvers are often stopped prematurely to recompute the linearisation of the non-linear outer loop. We show that in the reverse-mode of AD, the derivatives obtained with partial convergence become inconsistent with the original and the tangent-linear models, resulting in inaccurate adjoints. We present a correction term that restores consistency between adjoint and tangent-linear gradients if linear systems are only partially converged. We prove the consistency of this correction term and show in numerical experiments that the accuracy of adjoint gradients of an incompressible flow solver applied to an industrial test case is restored when the correction term is used.


Sign in / Sign up

Export Citation Format

Share Document