1887

Abstract

Summary

We explore and develop POD-based deflation methods to accelerate the solution of large-scale linear systems resulting from two-phase flow simulation.

The techniques here presented collect information from the system in a POD basis, which is later used in a deflation scheme.

The snapshots required to obtain the POD basis are captured in two ways: a moving window approach, where the most recently computed solutions are used, and a training phase approach, where a full pre-simulation is run. We test this methodology in highly heterogeneous porous media: a full SPE 10 model containing O(10^6) cells, and in an academic layered problem presenting a contrast in permeability layers up to 10^6. Among the experiments, we study cases including gravity and capillary pressure terms.

With the POD-based deflated procedure, we accelerate the convergence of a Preconditioned Conjugate Gradient (PCG) method, reducing the required number of iterations to around 10–30 %, i.e., we achieve speed-ups of factors three to ten.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201802122
2018-09-03
2024-03-28
Loading full text...

Full text loading...

References

  1. Astrid, P., Papaioannou, G., Vink, J. and Jansen, J.
    [2011] Pressure Preconditioning Using Proper Orthogonal Decomposition. In: 2011 SPE Reservoir Simulation Symposium, The Woodlands, Texas, USA. 21–23.
    [Google Scholar]
  2. Aubry, R., Mut, F., Lohner, R. and Cebral, J.
    [2008] Deflated preconditioned conjugate gradient solvers for the Pressure-Poisson equation. Journal of Computational Physics, 227(24), 10196–10208.
    [Google Scholar]
  3. Christie, M. and Blunt, M.
    [2001] Tenth SPE Comparative Solution Project: a Comparison of Upscaling Techniques. SPE Reservoir Engineering and Evaluation, 4(4), 308–317.
    [Google Scholar]
  4. Clemens, M., Wilke, M., Schuhmann, R. and Weiland, T.
    [2004] Subspace projection extrapolation scheme for transient field simulations. IEEE Transactions on Magnetics, 40(2), 934–937.
    [Google Scholar]
  5. Diaz-Cortes, G., Vuik, C. and Jansen, J.
    [2016] Physics-based Pre-conditioners for Large-scale Subsurface Flow Simulation. In: Proceedings of the 15th European Conference on the Mathematics of Oil Recovery, ECMOR XV.1031–1051.
    [Google Scholar]
  6. [2017] On POD-based Deflation Vectors for DPCG applied to porous media problems. Report 17-1, Delft University of Technology, Delft Institute of Applied Mathematics, Delft. http://ta.twi.tudelft.nl/nw/users/vuik/papers/Dia17VJ.pdf.
  7. [2018] On POD-based Deflation Vectors for DPCG applied to porous media problems.Journal of Computational and Applied Mathematics, 330(Supplement C), 193–213.
    [Google Scholar]
  8. Golub, G. and Loan, C.V.
    [1996] Matrix Computations.Johns Hopkins University Press, Baltimore, MD, USA, 3rd edn.
    [Google Scholar]
  9. Kahl, K. and Rittich, H.
    [2017] The deflated conjugate gradient method: Convergence, perturbation and accuracy.Linear Algebra and its Applications, 515, 111–129.
    [Google Scholar]
  10. Lie, K.
    [2013] An Introduction to Reservoir Simulation Using MATLAB: User guide for the Matlab Reservoir Simulation Toolbox (MRST). SINTEF ICT.
    [Google Scholar]
  11. Lohner, R. F., J.R.Cebral, J.R., Aubry, R. and Houzeaux, G.
    [2011] Deflated preconditioned conjugate gradient solvers for the pressure-Poisson equation: Extensions and improvements. International Journal for Numerical Methods in Engineering, 87(1-5), 2–14.
    [Google Scholar]
  12. Markovinović, R. and Jansen, J.
    [2006] Accelerating iterative solution methods using reduced-order models as solution predictors. International journal for numerical methods in engineering, 68(5), 525–541.
    [Google Scholar]
  13. Napov, A. and Notay, Y.
    [2012] An algebraic multigrid method with guaranteed convergence rate. SIAM journal on scientific computing, 34(2), A1079–A1109.
    [Google Scholar]
  14. Notay, Y.
    [2010] An aggregation-based algebraic multigrid method. Electronic transactions on numerical analysis, 37(6), 123–146.
    [Google Scholar]
  15. [2012] Aggregation-based algebraic multigrid for convection-diffusion equations. SIAM journal on scientific computing, 34(4), A2288–A2316.
    [Google Scholar]
  16. Pasetto, D., Ferronato, M. and Putti, M.
    [2017] A reduced order model-based preconditioner for the efficient solution of transient diffusion equations. International Journal for Numerical Methods in Engineering, 109(8), 1159–1179.
    [Google Scholar]
  17. Saad, Y., Yeung, M., Erhel, J. and Guyomarc’h, F.
    [2000] A deflated version of the conjugate gradient algorithm. SIAM Journal on Scientific Computing, 21(5), 1909–1926.
    [Google Scholar]
  18. Smith, B., Bjorstad, P. and Gropp, W.
    [1996] Domain decomposition: parallel multilevel methods for elliptic partial differential equations.Cambridge University PressNew York.
    [Google Scholar]
  19. Tang, J.
    [2008] Two-Level Preconditioned Conjugate Gradient Methods with Applications to Bubbly Flow Problems. Ph.D. thesis, Delft University of Technology.
    [Google Scholar]
  20. Tang, J., Nabben, R., Vuik, C. and Erlangga, Y.
    [2009] Comparison of two-level preconditioners derived from deflation, domain decomposition and multigrid methods. Journal of scientific computing, 39(3), 340–370.
    [Google Scholar]
  21. der Vorst, H.V.
    [2003] Iterative Krylov methods for large linear systems, 13. Cambridge University Press.
    [Google Scholar]
  22. der Vorst, H.V. and Dekker, C.
    [1988] Conjugate gradient type methods and preconditioning.Journal of Computational and Applied Mathematics, 24, 73–87.
    [Google Scholar]
  23. Vuik, C., Segal, A. and Meijerink, J.
    [1999] An Efficient Preconditioned CG Method for the Solution of a Class of Layered Problems with Extreme Contrasts in the Coefficients.Journal of Computational Physics, 152, 385.
    [Google Scholar]
  24. Vuik, C., Segal, A., Yaakoubi, L. and Dufour, E.
    [2002] A comparison of various deflation vectors applied to elliptic problems with discontinuous coefficients. Applied Numerical Mathematics, 41(1), 219–233.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201802122
Loading
/content/papers/10.3997/2214-4609.201802122
Loading

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error