A low-dissipation DG method for the under-resolved simulation of low Mach number turbulent flows

Abstract In recent years the use of high-order Discontinuous Galerkin (DG) methods for the under-resolved direct numerical simulations (uDNS) of turbulent flows has received special attention. The suitability of the approach for this kind of applications is related to the dissipation and dispersion proprieties of the scheme: while the dispersion errors are small over a broad range of frequencies, a relevant dissipation error mainly acts at the smallest under-resolved scales, resembling a high frequency filter. Nevertheless, it was recognized (Flad and Gassner, 2017; Mengaldo et al., 2018) that the choice of the interface convective numerical flux strongly affects this dissipation behaviour and ultimately the success of the uDNS approach. In this regard, the excess of numerical dissipation caused by some upwind numerical convective fluxes must be avoided, in particular when dealing with low-speed flows, since this behaviour is exacerbated approaching the incompressibility limit. Fixes for the excess of numerical dissipation of these schemes have been proposed by several authors in the context of different numerical methods, see for example Weiss and Smith (1995). In this work a simple modification of the dissipation term of the low Mach preconditioned Roe scheme proposed by Weiss and Smith is considered. The aim is to reduce further the amount of numerical dissipation with the intent of improving the results of uDNS. A spatial DG discretization coupled with a linearly-implicit Rosenbrock-type time integrator is here considered as a numerical framework perfectly suited for the assessment and comparison of different numerical flux functions. Results on canonical turbulent flow problems as the Taylor–Green vortex and the flow in a straight sided channel are presented. The improved accuracy of the proposed flux function is demonstrated. The new low-dissipation flux can be useful also in the context of standard, lower order, finite volume methods.

[1]  Alessandro Colombo,et al.  Linearly implicit Rosenbrock-type Runge–Kutta schemes applied to the Discontinuous Galerkin solution of compressible and incompressible unsteady flows , 2015 .

[2]  Leland Jameson,et al.  Numerical Convergence Study of Nearly Incompressible, Inviscid Taylor–Green Vortex Flow , 2005, J. Sci. Comput..

[3]  Alessandro Colombo,et al.  On the development of an implicit high-order Discontinuous Galerkin method for DNS and implicit LES of turbulent flows , 2016 .

[4]  Hester Bijl,et al.  A comparison of Rosenbrock and ESDIRK methods combined with iterative solvers for unsteady compressible flows , 2016, Advances in Computational Mathematics.

[5]  H. Guillard,et al.  On the behavior of upwind schemes in the low Mach number limit: II. Godunov type schemes , 2004 .

[6]  A. Nigro,et al.  On the efficiency of a matrix-free linearly implicit time integration strategy for high-order Discontinuous Galerkin solutions of incompressible turbulent flows , 2017 .

[7]  C. Munz,et al.  On the Effect of Flux Functions in Discontinuous Galerkin Simulations of Underresolved Turbulence , 2014 .

[8]  Vít Dolejší,et al.  Adaptive backward difference formula–Discontinuous Galerkin finite element method for the solution of conservation laws , 2008 .

[9]  Alexander Jaust,et al.  Implicit Multiderivative Collocation Solvers for Linear Partial Differential Equations with Discontinuous Galerkin Spatial Discretizations , 2017, Journal of Scientific Computing.

[10]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..

[11]  Ralf Hartmann,et al.  Discontinuous Galerkin solution of preconditioned Euler equations for very low Mach number flows , 2009 .

[12]  Ralf Hartmann.  Error estimation and adjoint based refinement for an adjoint consistent DG discretisation of the compressible Euler equations , 2007, Int. J. Comput. Sci. Math..

[13]  Francesco Bassi,et al.  Modified extended BDF scheme for the discontinuous Galerkin solution of unsteady compressible flows , 2014 .

[14]  Gregor Gassner,et al.  On the use of kinetic energy preserving DG-schemes for large eddy simulation , 2017, J. Comput. Phys..

[15]  P. Tesini,et al.  On the flexibility of agglomeration based physical space discontinuous Galerkin discretizations , 2012, J. Comput. Phys..

[16]  Alessandro Colombo,et al.  Time Integration in the Discontinuous Galerkin Code MIGALE - Unsteady Problems , 2015 .

[17]  Petros Koumoutsakos,et al.  A comparison of vortex and pseudo-spectral methods for the simulation of periodic vortical flows at high Reynolds numbers , 2011, J. Comput. Phys..

[18]  Spencer J. Sherwin,et al.  On the eddy-resolving capability of high-order discontinuous Galerkin approaches to implicit LES / under-resolved DNS of Euler turbulence , 2017, J. Comput. Phys..

[19]  P. Roe.  CHARACTERISTIC-BASED SCHEMES FOR THE EULER EQUATIONS , 1986 .

[20]  A. Ern,et al.  Mathematical Aspects of Discontinuous Galerkin Methods , 2011 .

[21]  J. Williamson.  Low-storage Runge-Kutta schemes , 1980 .

[22]  A. Beck,et al.  On the accuracy of high-order discretizations for underresolved turbulence simulations , 2013 .

[23]  S. Rebay,et al.  A High-Order Accurate Discontinuous Finite Element Method for the Numerical Solution of the Compressible Navier-Stokes Equations , 1997 .

[24]  Andrea Crivellini,et al.  Galerkin solutions of unsteady compressible flows 2 , 2017 .

[25]  Alessandro Colombo,et al.  h-multigrid agglomeration based solution strategies for discontinuous Galerkin discretizations of incompressible flow problems , 2017, J. Comput. Phys..

[26]  Francesco Bassi,et al.  Up to sixth-order accurate A-stable implicit schemes applied to the Discontinuous Galerkin discretized Navier-Stokes equations , 2014, J. Comput. Phys..

[27]  Clinton P. T. Groth,et al.  Assessment of Riemann solvers for unsteady one-dimensional inviscid flows for perfect gases , 1988 .

[28]  Steven J. Ruuth,et al.  A New Class of Optimal High-Order Strong-Stability-Preserving Time Discretization Methods , 2002, SIAM J. Numer. Anal..

[29]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

[30]  Rémi Abgrall,et al.  High‐order CFD methods: current status and perspective , 2013 .

[31]  Ralf Hartmann,et al.  A high-order accurate discontinuous Galerkin finite element method for laminar low Mach number flows , 2013 .

[32]  Vít Dolejší,et al.  Discontinuous Galerkin Method: Analysis and Applications to Compressible Flow , 2015 .

[33]  Spencer J. Sherwin,et al.  Spatial eigensolution analysis of discontinuous Galerkin schemes with practical insights for under-resolved computations and implicit LES , 2017, Computers & Fluids.

[34]  John Kim,et al.  DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOWS UP TO RE=590 , 1999 .

[35]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[36]  Ralf Hartmann,et al.  A discontinuous Galerkin method for inviscid low Mach number flows , 2009, J. Comput. Phys..

[37]  Koen Hillewaert,et al.  Assessment of a discontinuous Galerkin method for the simulation of vortical flows at high Reynolds number , 2014 .

[38]  Per-Olof Persson,et al.  Stage-parallel fully implicit Runge-Kutta solvers for discontinuous Galerkin fluid simulations , 2017, J. Comput. Phys..

[39]  Thomas J. R. Hughes,et al.  A comparative study of different sets of variables for solving compressible and incompressible flows , 1998 .

[40]  S. Orszag,et al.  Small-scale structure of the Taylor–Green vortex , 1983, Journal of Fluid Mechanics.

[41]  Per-Olof Persson,et al.  Newton-GMRES Preconditioning for Discontinuous Galerkin Discretizations of the Navier--Stokes Equations , 2008, SIAM J. Sci. Comput..

[42]  Wayne A. Smith,et al.  Preconditioning Applied to Variable and Constant Density Flows , 1995 .

[43]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .