On the efficiency of a matrix-free linearly implicit time integration strategy for high-order Discontinuous Galerkin solutions of incompressible turbulent flows. (15th December 2017)

• Record Type: Journal Article
• Authors: Franciolini, Matteo; Crivellini, Andrea; Nigro, Alessandra
• Is Part Of: Computers & fluids. Volume 159(2017)
• Page Start: 276
• Page End: 294
• Publication Date: 2017-12-15
• Subjects: Discontinuous Galerkin -- Implicit time integration -- Linearly-implicit Rosenbrock schemes -- Matrix-free solvers -- DNS/ILES; Fluid dynamics -- Data processing -- Periodicals; 532.050285
• DOI: 10.1016/j.compfluid.2017.10.008
• Languages: English
• ISSNs: 0045-7930

Abstract

Highlights: Matrix-free linearly implicit time integration for an incompressible DG solver. Adaptive GMRES solution tolerance estimation. Application of the matrix-free solver to 2D laminar and 3D turbulent flows. Memory footprint reduction up 93% for high order polynomials and 3D cases. Matrix-free up to five time faster than matrix-based for DNS/ILES. Abstract: The paper aims at investigating and comparing the efficiency of matrix-based and matrix-free implementations of implicit time integration strategies, here linearly implicit Rosenbrock-type Runge–Kutta schemes, in the context of a high order Discontinous Galerkin space discretization applied to DNS/ILES of turbulent flows. In particular, the effects of time step size, GMRES solution tolerance and preconditioner, mesh size, order of polynomial approximation of the solution and parallel efficiency are here considered. The performance of the solvers have been compared though the solution of the two-dimensional laminar traveling waves test case, the ILES of the Rayleigh–Bénard natural convection up to R a = 10 8, and the channel flow up to R e τ = 950 . The results highlight that while matrix-based strategies are the most adequate for low-order polynomial approximations, considerable improvements in computational efficiency can be achieved, both in terms of CPU time and memory footprint, by a proper use of the matrix-free solver for three-dimensional test cases and very high-order accurate DG space discretizations.