On the Influence of Polynomial De-aliasing on Subgrid Scale Models

被引:30
作者
Beck, Andrea D. [1 ]
Flad, David G. [1 ]
Tonhaeuser, Claudia [1 ]
Gassner, Gregor [2 ]
Munz, Claus-Dieter [1 ]
机构
[1] Univ Stuttgart, Inst Aerodynam & Gasdynam, Pfaffenwaldring 21, D-70550 Stuttgart, Germany
[2] Univ Cologne, Math Inst, Weyertal 86-90, D-50923 Cologne, Germany
关键词
Discontinuous Galerkin; Large eddy simulation; Stability; Smagorinsky model; Variational multiscale method; De-aliasing; Collocation; LARGE-EDDY SIMULATION; FINITE-ELEMENT-METHOD; DISCONTINUOUS GALERKIN METHODS; CIRCULAR-CYLINDER; CONSERVATION-LAWS; TURBULENCE; FLOW; ERRORS;
D O I
10.1007/s10494-016-9704-y
中图分类号
O414.1 [热力学];
学科分类号
摘要
In this work we investigate the interplay of polynomial de-aliasing and sub-grid scale models for large eddy simulations based on discontinuous Galerkin discretizations. It is known that stability is a major concern when simulating underresolved turbulent flows with high order nodal collocation type discretizations. By changing the interpolatory character of the nodal collocation type discretization to a projection based discretization by increasing the number of quadrature points (polynomial de-aliasing), one is able to remove the aliasing induced stability problems. We focus on this effect and on the consequence for large eddy simulations with explicit subgrid scale models. Often, subgrid scale models have to achieve two possibly conflicting tasks in a single simulation: firstly stabilizing the numerics and secondly modeling the physical effect of the missing scales. Within a discontinuous Galerkin approach, it is possible to use either a fast (but potentially aliasing afflicted) nodal collocation discretization or a projection-based (but computationally costly) variant in combination with an explicit subgrid scale model. We use this framework to investigate the effect on the appropriate model parameter of a standard Smagorinsky subgrid scale model and of a Variational Multiscale Smagorinsky formulation. For this we first consider the 3-D viscous Taylor-Green vortex example to investigate the impact on the stability of the method and second the turbulent flow past a circular cylinder to investigate and compare the accuracy of the results. We show that the aliasing instabilities of collocative discretizations severely limit the choice of the model constant, in particular for high order schemes, while for de-aliased DG schemes, the closure model parameters can be chosen independently from the numerical scheme. For the cylinder flow, we also find that for the same model settings, the projection-based results are in better agreement with the reference DNS than those of the collocative scheme.
引用
收藏
页码:475 / 511
页数:37
相关论文
共 76 条
[1]  
[Anonymous], 32 AIAA FLUID DYN C
[2]  
[Anonymous], 2013, LARGE EDDY SIMULATIO
[3]  
Arnold DN, 1979, THESIS
[4]   Numerical evaluation of two discontinuous Galerkin methods for the compressible Navier-Stokes equations [J].
Bassi, F ;
Rebay, S .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, 2002, 40 (1-2) :197-207
[5]   A high-order accurate discontinuous finite element method for the numerical solution of the compressible Navier-Stokes equations [J].
Bassi, F ;
Rebay, S .
JOURNAL OF COMPUTATIONAL PHYSICS, 1997, 131 (02) :267-279
[6]  
Bassi F., 1997, 2 EUROPEAN C TURBOMA, P99
[7]  
Baumann CE, 1999, INT J NUMER METH FL, V31, P79, DOI 10.1002/(SICI)1097-0363(19990915)31:1<79::AID-FLD956>3.0.CO
[8]  
2-C
[9]   Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows [J].
Bazilevs, Y. ;
Calo, V. M. ;
Cottrell, J. A. ;
Hughes, T. J. R. ;
Reali, A. ;
Scovazzi, G. .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2007, 197 (1-4) :173-201
[10]  
Beck A., 2013, INT J NUMERICA UNPUB