## On the efficiency of a matrix-free linearly implicit time integration strategy for high-order discontinuous Galerkin solutions of incompressible turbulent flows.(English)Zbl 1390.76312

Summary: 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 $$\mathrm{Ra} = 10^8$$, and the channel flow up to $$\mathrm{Re}_{\tau}= 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.

### MSC:

 76M10 Finite element methods applied to problems in fluid mechanics 76F65 Direct numerical and large eddy simulation of turbulence 65M60 Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs 76D05 Navier-Stokes equations for incompressible viscous fluids

### Software:

PETSc; Spalart-Allmaras; ROS3P; NITSOL
Full Text: