On the use of kinetic energy preserving DG-schemes for large eddy simulation

Abstract Recently, element based high order methods such as Discontinuous Galerkin (DG) methods and the closely related flux reconstruction (FR) schemes have become popular for compressible large eddy simulation (LES). Element based high order methods with Riemann solver based interface numerical flux functions offer an interesting dispersion dissipation behavior for multi-scale problems: dispersion errors are very low for a broad range of scales, while dissipation errors are very low for well resolved scales and are very high for scales close to the Nyquist cutoff. In some sense, the inherent numerical dissipation caused by the interface Riemann solver acts as a filter of high frequency solution components. This observation motivates the trend that element based high order methods with Riemann solvers are used without an explicit LES model added. Only the high frequency type inherent dissipation caused by the Riemann solver at the element interfaces is used to account for the missing sub-grid scale dissipation. Due to under-resolution of vortical dominated structures typical for LES type setups, element based high order methods suffer from stability issues caused by aliasing errors of the non-linear flux terms. A very common strategy to fight these aliasing issues (and instabilities) is so-called polynomial de-aliasing, where interpolation is exchanged with projection based on an increased number of quadrature points. In this paper, we start with this common no-model or implicit LES (iLES) DG approach with polynomial de-aliasing and Riemann solver dissipation and review its capabilities and limitations. We find that the strategy gives excellent results, but only when the resolution is such, that about 40% of the dissipation is resolved. For more realistic, coarser resolutions used in classical LES e.g. of industrial applications, the iLES DG strategy becomes quite inaccurate. We show that there is no obvious fix to this strategy, as adding for instance a sub-grid-scale models on top doesn't change much or in worst case decreases the fidelity even more. Finally, the core of this work is a novel LES strategy based on split form DG methods that are kinetic energy preserving. The scheme offers excellent stability with full control over the amount and shape of the added artificial dissipation. This premise is the main idea of the work and we will assess the LES capabilities of the novel split form DG approach when applied to shock-free, moderate Mach number turbulence. We will demonstrate that the novel DG LES strategy offers similar accuracy as the iLES methodology for well resolved cases, but strongly increases fidelity in case of more realistic coarse resolutions.

[1]  Gregor Gassner,et al.  A Comparison of the Dispersion and Dissipation Errors of Gauss and Gauss-Lobatto Discontinuous Galerkin Spectral Element Methods , 2011, SIAM J. Sci. Comput..

[2]  T. Hughes,et al.  Large Eddy Simulation and the variational multiscale method , 2000 .

[3]  Praveen Chandrashekar,et al.  Kinetic energy preserving and entropy stable finite volume schemes for compressible Euler and Navier-Stokes equations , 2012, ArXiv.

[4]  Heinz Pitsch,et al.  High order conservative finite difference scheme for variable density low Mach number turbulent flows , 2007, J. Comput. Phys..

[5]  Philipp Birken,et al.  L2Roe: a low dissipation version of Roe's approximate Riemann solver for low Mach numbers , 2016 .

[6]  George Em Karniadakis,et al.  De-aliasing on non-uniform grids: algorithms and applications , 2003 .

[7]  R. Kraichnan.  Eddy Viscosity in Two and Three Dimensions , 1976 .

[8]  Hester Bijl,et al.  Fourth-Order Runge–Kutta Schemes for Fluid Mechanics Applications , 2005, J. Sci. Comput..

[9]  Victor Lee,et al.  New Trends in Large-Eddy Simulations of Turbulence , 2011 .

[10]  Nikolaus A. Adams,et al.  An adaptive local deconvolution method for implicit LES , 2005, J. Comput. Phys..

[11]  Laslo T. Diosady,et al.  A Space-Time Discontinuous-Galerkin Approach for Separated Flows , 2016 .

[12]  Jung Yul Yoo,et al.  Discretization errors in large eddy simulation: on the suitability of centered and upwind-biased compact difference schemes , 2004 .

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

[14]  Gregor Gassner,et al.  A Skew-Symmetric Discontinuous Galerkin Spectral Element Discretization and Its Relation to SBP-SAT Finite Difference Methods , 2013, SIAM J. Sci. Comput..

[15]  Steven H. Frankel,et al.  Entropy Stable Spectral Collocation Schemes for the Navier-Stokes Equations: Discontinuous Interfaces , 2014, SIAM J. Sci. Comput..

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

[17]  Claus-Dieter Munz,et al.  Simulation of underresolved turbulent flows by adaptive filtering using the high order discontinuous Galerkin spectral element method , 2016, J. Comput. Phys..

[18]  P. Moin,et al.  A dynamic subgrid‐scale eddy viscosity model , 1990 .

[19]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .

[20]  J. Chollet.  Two-Point Closure Used for a Sub-Grid Scale Model in Large Eddy Simulations , 1985 .

[21]  Per-Olof Persson,et al.  Implicit Large Eddy Simulation of transition to turbulence at low Reynolds numbers using a Discontinuous Galerkin method , 2011 .

[22]  Claus-Dieter Munz,et al.  High‐order discontinuous Galerkin spectral element methods for transitional and turbulent flow simulations , 2014 .

[23]  Stefan Hickel,et al.  Subgrid-scale modeling for implicit large eddy simulation of compressible flows and shock-turbulence interaction , 2014 .

[24]  P. Moin,et al.  On the Effect of Numerical Errors in Large Eddy Simulations of Turbulent Flows , 1997 .

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

[26]  Spencer J. Sherwin,et al.  Linear dispersion-diffusion analysis and its application to under-resolved turbulence simulations using discontinuous Galerkin spectral/hp methods , 2015, J. Comput. Phys..

[27]  Laslo T. Diosady,et al.  Higher-Order Methods for Compressible Turbulent Flows Using Entropy Variables , 2015 .

[28]  Gregor Gassner,et al.  Split form nodal discontinuous Galerkin schemes with summation-by-parts property for the compressible Euler equations , 2016, J. Comput. Phys..

[29]  Joel H. Ferziger,et al.  A robust high-order compact method for large eddy simulation , 2003 .

[30]  Gregor Gassner,et al.  The BR1 Scheme is Stable for the Compressible Navier–Stokes Equations , 2017, J. Sci. Comput..

[31]  Claus-Dieter Munz,et al.  On the Influence of Polynomial De-aliasing on Subgrid Scale Models , 2016 .

[32]  Spencer J. Sherwin,et al.  Dealiasing techniques for high-order spectral element methods on regular and irregular grids , 2015, J. Comput. Phys..

[33]  Travis C. Fisher,et al.  High-order entropy stable finite difference schemes for nonlinear conservation laws: Finite domains , 2013, J. Comput. Phys..

[34]  D. Lilly,et al.  A proposed modification of the Germano subgrid‐scale closure method , 1992 .

[35]  Sergio Pirozzoli,et al.  Numerical Methods for High-Speed Flows , 2011 .

[36]  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..

[37]  L. Bricteux,et al.  Implicit LES of free and wall‐bounded turbulent flows based on the discontinuous Galerkin/symmetric interior penalty method , 2015 .

[38]  D. Drikakis,et al.  Large eddy simulation using high-resolution and high-order methods , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[39]  P. Moin,et al.  Suitability of upwind-biased finite difference schemes for large-eddy simulation of turbulent flows , 1997 .

[40]  Y. Kaneda,et al.  Effects of Wavenumber Truncation on High-Resolution Direct Numerical Simulation of Turbulence , 2002 .

[41]  Christopher A. Kennedy,et al.  Reduced aliasing formulations of the convective terms within the Navier-Stokes equations for a compressible fluid , 2008, J. Comput. Phys..