An implicit high-order spectral difference approach for large eddy simulation

The filtered fluid dynamic equations are discretized in space by a high-order spectral difference (SD) method coupled with large eddy simulation (LES) approach. The subgrid-scale stress tensor is modelled by the wall-adapting local eddy-viscosity model (WALE). We solve the unsteady equations by advancing in time using a second-order backward difference formulae (BDF2) scheme. The nonlinear algebraic system arising from the time discretization is solved with the nonlinear lower-upper symmetric Gauss-Seidel (LU-SGS) algorithm. In order to study the sensitivity of the method, first, the implicit solver is used to compute the two-dimensional (2D) laminar flow around a NACA0012 airfoil at Re=5x10^5 with zero angle of attack. Afterwards, the accuracy and the reliability of the solver are tested by solving the 2D ''turbulent'' flow around a square cylinder at Re=10^4 and Re= 2.2x10^4. The results show a good agreement with the experimental data and the reference solutions.

[1]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[2]  A. Mochida,et al.  On turbulent vortex shedding flow past 2D square cylinder predicted by CFD , 1995 .

[3]  Chi-Wang Shu,et al.  The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems , 1998 .

[4]  Chris Lacor,et al.  Dispersion and dissipation properties of the 1D spectral volume method and application to a p-multigrid algorithm , 2007, J. Comput. Phys..

[5]  Zhi J. Wang,et al.  Efficient Implicit LU-SGS Algorithm for High-Order Spectral Difference Method on Unstructured Hexahedral Grids , 2007 .

[6]  C. Fureby,et al.  Mathematical and Physical Constraints on Large-Eddy Simulations , 1997 .

[7]  Chris Lacor,et al.  An accuracy and stability study of the 2D spectral volume method , 2007, J. Comput. Phys..

[8]  Jan G. Wissink,et al.  DNS of 2D turbulent flow around a square cylinder , 1997 .

[9]  T. Iliescu,et al.  A Numerical Study of a Class of LES Models , 2003 .

[10]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[11]  Marcel Vinokur,et al.  Spectral difference method for unstructured grids I: Basic formulation , 2006, J. Comput. Phys..

[12]  Marcel Vinokur,et al.  Spectral (finite) volume method for conservation laws on unstructured grids V: Extension to three-dimensional systems , 2006, J. Comput. Phys..

[13]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[14]  T. A. Zang,et al.  Direct and large-eddy simulations of three-dimensional compressible Navier-Stokes turbulence , 1992 .

[15]  Yuzhi Sun,et al.  Spectral (finite) volume method for conservation laws on unstructured grids VI: Extension to viscous flow , 2006, J. Comput. Phys..

[16]  Chi-Wang Shu,et al.  The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case , 1990 .

[17]  Zhi Jian Wang,et al.  Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids III: One Dimensional Systems and Partition Optimization , 2004, J. Sci. Comput..

[18]  A. Jameson,et al.  Implicit schemes and LU decompositions , 1981 .

[19]  Shmuel Einav,et al.  A laser-Doppler velocimetry study of ensemble-averaged characteristics of the turbulent near wake of a square cylinder , 1995, Journal of Fluid Mechanics.

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

[21]  Antony Jameson,et al.  Spectral Difference Method for Unstructured Grids II: Extension to the Euler Equations , 2007, J. Sci. Comput..

[22]  Douglas N. Arnold,et al.  Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems , 2001, SIAM J. Numer. Anal..

[23]  D. Bouris,et al.  2D LES of vortex shedding from a square cylinder , 1999 .

[24]  Haiyang Gao,et al.  A High-Order Lifting Collocation Penalty Formulation for the Navier-Stokes Equations on 2-D Mixed Grids , 2009 .

[25]  Z. Wang,et al.  Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids , 2002 .

[26]  Wolfgang Rodi,et al.  On the Simulation of Turbulent Flow Past Bluff Bodies , 1993 .

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

[28]  P. Lax.  Weak solutions of nonlinear hyperbolic equations and their numerical computation , 1954 .

[29]  Richard D. Sandberg,et al.  Direct numerical simulations of tonal noise generated by laminar flow past airfoils , 2009 .

[30]  W. Rodi,et al.  Calculation of Vortex Shedding Past a Square Cylinder with Various Turbulence Models , 1993 .

[31]  H. T. Huynh,et al.  A Flux Reconstruction Approach to High-Order Schemes Including Discontinuous Galerkin Methods , 2007 .

[32]  Zhi J. Wang,et al.  Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids. Basic Formulation , 2002 .

[33]  Zhi J. Wang,et al.  Spectral (finite) volume method for conservation laws on unstructured grids IV: extension to two-dimensional systems , 2004 .

[34]  Manuel V. Heitor,et al.  Measurements of turbulent and periodic flows around a square cross-section cylinder , 1988 .

[35]  Eli Turkel,et al.  Implicit LU-SGS algorithm for high-order methods on unstructured grid with p-multigrid strategy for solving the steady Navier-Stokes equations , 2010, J. Comput. Phys..

[36]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems , 1989 .

[37]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[38]  Antony Jameson,et al.  Lower-upper implicit schemes with multiple grids for the Euler equations , 1987 .

[39]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

[40]  Michael Breuer,et al.  First Experiences with LES of Flows past Bluff Bodies , 1996 .

[41]  Zhi J. Wang,et al.  Extension of the spectral volume method to high-order boundary representation , 2006 .

[42]  Eli Turkel,et al.  Simulation of Compressible Turbulent Flows with an Implicit LU-SGS Algorithm for High-Order Spectral Difference Method on Unstructured Grids , 2009 .

[43]  James V. Lambers.  Practical Implementation of Krylov Subspace Spectral Methods , 2007, J. Sci. Comput..

[44]  Chris Lacor,et al.  On the Stability and Accuracy of the Spectral Difference Method , 2008, J. Sci. Comput..

[45]  H. C. Yee,et al.  Entropy Splitting for High Order Numerical Simulation of Compressible Turbulence , 2002 .

[46]  F. Nicoud,et al.  Subgrid-Scale Stress Modelling Based on the Square of the Velocity Gradient Tensor , 1999 .