In this work, a time-accurate integration of the discontinuous Galerkin space-discretized Navier-Stokes equations is performed exploiting the matrix-free (MF) approach to speed up the solution process of the modified extended backward differentiation formulae (MEBDF) schemes. MEBDF are high-order accurate implicit multistep schemes composed by three nonlinear stages. The proposed algorithm consists in solving the resulting nonlinear sys- tem of each stage with a preconditioned MF Newton/Krylov method using a frozen preconditioner strategy to improve its efficiency. Numerical results for compressible inviscid and viscous test cases, both with a known analytical solu- tion, aim at assessing the performance of the proposed MF-MEBDF algorithm, by comparing it with the one obtained by using its matrix-explicit counterpart or the explicit strong stability preserving Runge-Kutta scheme. In particular, the influence of some relevant physical (low-speed flows) and discretization (aspect ratio, polynomial degree) aspects on the performance of the different time integration schemes are investigated, highlighting the pros and cons of the proposed algorithm and its effectiveness in solving nonstiff and stiff systems. Furthermore, the scalability in parallel computations of the proposed algorithm is investigated and, finally, its potential for efficient long-time simulations is demonstrated computing a laminar vortex shedding behind a circular cylinder at different Reynolds numbers and by comparing its effectiveness with that of the strong stability preserving Runge-Kutta scheme.

Matrix‐free modified extended BDF applied to the discontinuous Galerkin solution of unsteady compressible viscous flows

Nigro, A.
;
De Bartolo, C.;
2018

Abstract

In this work, a time-accurate integration of the discontinuous Galerkin space-discretized Navier-Stokes equations is performed exploiting the matrix-free (MF) approach to speed up the solution process of the modified extended backward differentiation formulae (MEBDF) schemes. MEBDF are high-order accurate implicit multistep schemes composed by three nonlinear stages. The proposed algorithm consists in solving the resulting nonlinear sys- tem of each stage with a preconditioned MF Newton/Krylov method using a frozen preconditioner strategy to improve its efficiency. Numerical results for compressible inviscid and viscous test cases, both with a known analytical solu- tion, aim at assessing the performance of the proposed MF-MEBDF algorithm, by comparing it with the one obtained by using its matrix-explicit counterpart or the explicit strong stability preserving Runge-Kutta scheme. In particular, the influence of some relevant physical (low-speed flows) and discretization (aspect ratio, polynomial degree) aspects on the performance of the different time integration schemes are investigated, highlighting the pros and cons of the proposed algorithm and its effectiveness in solving nonstiff and stiff systems. Furthermore, the scalability in parallel computations of the proposed algorithm is investigated and, finally, its potential for efficient long-time simulations is demonstrated computing a laminar vortex shedding behind a circular cylinder at different Reynolds numbers and by comparing its effectiveness with that of the strong stability preserving Runge-Kutta scheme.
discontinuous Galerkin, frozen preconditioner strategy, matrix-free approximation, modified extended backward differentiation formulae, stiff and nonstiff system of equations, unsteady flows
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/20.500.11770/291293
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus ND
  • ???jsp.display-item.citation.isi??? 1
social impact