• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.7
• /
• pp.995-1006
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes for large eddy simulation is evaluated by a dynamic analysis. Large eddy simulation of isotropic turbulence is performed with various dissipative and non-dissipative schemes to investigate the effect of numerical dissipation on the resolved solutions. It is shown by the present dynamic analysis that upwind schemes reduce the aliasing error and increase the finite differencing error. The existence of optimal upwind scheme that minimizes total numerical error is verified. It is also shown that the finite differencing error from numerical dissipation is the leading source of numerical errors by upwind schemes. Simulations of a turbulent channel flow are conducted to show the existence of the optimal upwind scheme.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.7
• /
• pp.973-983
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes is evaluated for large eddy simulation of turbulent flow. Two turbulent flows are considered: turbulent channel flow at Re = 23000 and flow over a circular cylinder at Re = 3900. The effects of numerical dissipation on the finite differencing and aliasing errors and the subgrid-scale stress are investigated. It is shown through the simulations that compact upwind schemes are not suitable for LES, whereas the fourth order-compact centered scheme is a good candidate for LES provided that proper dealiasing of nonlinear terms is performed. The classical issue on the aliasing error and the treatment of nonlinear terms is revisited with compact difference schemes.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.7
• /
• pp.984-994
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes for large eddy simulation is evaluated by a spectral, static error analysis. To investigate the effect of numerical dissipation on LES solutions, power spectra of discretization errors are evaluated for isotropic turbulence models in both continuous and discrete wavevector spaces. Contrary to the common belief, the aliasing errors from upwind-biased schemes are larger than those from comparable non-dissipative schemes. However, this result is the direct consequence of the definition of the power spectral density of the aliasing error, which poses the limitation of the static error analysis for upwind schemes.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.392-397
• /
• 2001

A numerical study is performed to investigate the interaction between subsonic axial turbine blade boundary layer and periodically oncoming rotor induced wakes. An implicit scheme for solving the compressible Navier-Stokes equation is developed, which adopts a 4th-order compact difference for spatial discretiztion, a 2nd order Crank-Nicolson scheme for temporal discretization and the dynamic eddy viscosity model as the subgrid scale model. The efficiency and the accuracy of the proposed method are verified by applying to some benchmark problems such as laminar cylinder flow, laminar airfoil cascade flow and a transitional flat plate boundary layer flow. Computational results show good agreements with previous experimental and numerical results. Finally, flow through a stator cascade is simulated at $Re = 7.5{\times}10^5$ without free-stream turbulence intensity. The velocity fields and skin friction coefficients in the transitional region show similar trends with previous boundary layer natural transition.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.807-810
• /
• 2002

The suitability of high-order accurate, central and upwind-biased compact difference schemes is evaluated for the large-eddy simulations of flows in complex geometry. Two flow geometries are considered: channel and circular cylinder. The effects of numerical dissipation and aliasing error on the evaluation of subgrid scale stress are investigated by extending the analysis by Ghosal (1) to centered and upwind compact schemes. It is shown that the failure of upwind schemes mainly comes from the aliasing error.