TY - JOUR

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

AU - Parsani, M.

AU - Ghorbaniasl, G.

AU - Lacor, C.

AU - Turkel, E.

N1 - Funding Information:
This research was funded by IWT under Project SBO 050163 . This funding is gratefully acknowledged.

PY - 2010/7

Y1 - 2010/7

N2 - 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=5×105 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=104 and Re= 2.2×104. The results show a good agreement with the experimental data and the reference solutions.

AB - 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=5×105 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=104 and Re= 2.2×104. The results show a good agreement with the experimental data and the reference solutions.

KW - High-order spectral difference method

KW - Implicit LU-SGS algorithm

KW - Large eddy simulation

KW - Wall-adapting local eddy-viscosity model

UR - http://www.scopus.com/inward/record.url?scp=77952882930&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2010.03.038

DO - 10.1016/j.jcp.2010.03.038

M3 - מאמר

AN - SCOPUS:77952882930

VL - 229

SP - 5373

EP - 5393

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 14

ER -