Simulation of underresolved turbulent flows by adaptive filtering using the high order discontinuous Galerkin spectral element method

Scale-resolving simulations of turbulent flows in complex domains demand accurate and efficient numerical schemes, as well as geometrical flexibility. For underresolved situations, the avoidance of aliasing errors is a strong demand for stability. For continuous and discontinuous Galerkin schemes, an effective way to prevent aliasing errors is to increase the quadrature precision of the projection operator to account for the non-linearity of the operands (polynomial dealiasing, overintegration). But this increases the computational costs extensively. In this work, we present a novel spatially and temporally adaptive dealiasing strategy by projection filtering. We show this to be more efficient for underresolved turbulence than the classical overintegration strategy. For this novel approach, we discuss the implementation strategy and the indicator details, show its accuracy and efficiency for a decaying homogeneous isotropic turbulence and the transitional Taylor-Green vortex and compare it to the original overintegration approach and a state of the art variational multi-scale eddy viscosity formulation.

[1]  S. Collis,et al.  Monitoring unresolved scales in multiscale turbulence modeling , 2001 .

[2]  E. Toro Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

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

[4]  Analysis of Relaxation Regularization as Subgrid-Scale Model for Large-Eddy Simulation , 2004 .

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

[6]  Leonhard Kleiser,et al.  High-pass filtered eddy-viscosity models for large-eddy simulations of transitional and turbulent flow , 2005 .

[7]  Wei Liu,et al.  An analysis of subgrid-scale interactions in numerically simulated isotropic turbulence , 1993 .

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

[9]  Steven A. Orszag,et al.  On the Elimination of Aliasing in Finite-Difference Schemes by Filtering High-Wavenumber Components , 1971 .

[10]  David A. Kopriva,et al.  Implementing Spectral Methods for Partial Differential Equations: Algorithms for Scientists and Engineers , 2009 .

[11]  S. A. Orsag,et al.  Small-scale structure of the Taylor-Green vortex , 1984 .

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

[13]  P. Schlatter Large-eddy simulation of transition and turbulence in wall-bounded shear flow , 2005 .

[14]  M. Carpenter,et al.  Fourth-order 2N-storage Runge-Kutta schemes , 1994 .

[15]  J. A. Domaradzki,et al.  Truncated Navier–Stokes equations with the automatic filtering criterion: Reynolds stress and energy budgets , 2008 .

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

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

[18]  J. A. Domaradzki,et al.  Large eddy simulations using truncated Navier–Stokes equations with the automatic filtering criterion , 2008 .

[19]  P. Swarztrauber The Vector Harmonic Transform Method for Solving Partial Differential Equations in Spherical Geometry , 1993 .

[20]  Jan S. Hesthaven,et al.  A spectral multidomain penalty method model for the simulation of high Reynolds number localized incompressible stratified turbulence , 2005 .

[21]  M. Brachet Direct simulation of three-dimensional turbulence in the Taylor–Green vortex , 1991 .

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

[23]  J. Domaradzki,et al.  Large eddy simulations of decaying rotating turbulence , 2004 .

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

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

[26]  George Em Karniadakis,et al.  A Spectral Vanishing Viscosity Method for Large-Eddy Simulations , 2000 .

[27]  David A. Kopriva,et al.  Metric Identities and the Discontinuous Spectral Element Method on Curvilinear Meshes , 2006, J. Sci. Comput..

[28]  Hyeong-Bin Cheong,et al.  Fourier-series representation and projection of spherical harmonic functions , 2012, Journal of Geodesy.

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