High Reynolds number calculations using the dynamic subgrid‐scale stress model

1993 ◽  
Vol 5 (6) ◽  
pp. 1484-1490 ◽  
Author(s):  
Ugo Piomelli
1997 ◽  
Vol 336 ◽  
pp. 151-182 ◽  
Author(s):  
BRANKO KOSOVIĆ

It has been recognized that the subgrid-scale (SGS) parameterization represents a critical component of a successful large-eddy simulation (LES). Commonly used linear SGS models produce erroneous mean velocity profiles in LES of high-Reynolds-number boundary layer flows. Although recently proposed approaches to solving this problem have resulted in significant improvements, questions about the true nature of the SGS problem in shear-driven high-Reynolds-number flows remain open.We argue that the SGS models must capture inertial transfer effects including backscatter of energy as well as its redistribution among the normal SGS stress components. These effects are the consequence of nonlinear interactions and anisotropy. In our modelling procedure we adopt a phenomenological approach whereby the SGS stresses are related to the resolved velocity gradients. We show that since the SGS stress tensor is not frame indifferent a more general nonlinear model can be applied to the SGS parameterization. We develop a nonlinear SGS model capable of reproducing the effects of SGS anisotropy characteristic for shear-driven boundary layers. The results obtained using the nonlinear model for the LES of a neutral shear-driven atmospheric boundary layer show a significant improvement in prediction of the non-dimensional shear and low-order statistics compared to the linear Smagorinsky-type models. These results also demonstrate a profound effect of the SGS model on the flow structures.


2011 ◽  
Vol 110-116 ◽  
pp. 4719-4722
Author(s):  
V. Parthiban ◽  
Ashwin Russelle

In order to predict a turbulent flow around a triangular cylinder a high Reynolds number of 45000 is done in the numerical simulation. In this simulation both steady and unsteady vortex shedding is predicted and various time steps. The numerical method used in this simulation is Reynolds Stress model. For steady and unsteady solution velocity contours and velocity vector plots is to be predicted for the vortex shedding behind the triangular cylinder.


Author(s):  
Peter A. Chang ◽  
Meng Wang ◽  
Jonathan Gershfeld

ATTACHED, wall-bounded flows impose computational requirements on LES that increase drastically with Reynolds number. For that reason, even simple geometries, such as airfoils at small angles of attack, with spanwise uniform section shape, challenge the bounds of LES as chord-based Reynolds numbers increase much above 1 million. Of particular concern is the ability of LES to predict the occurrence, and strength of, weak vortex shedding from the airfoil trailing edge (by weak vortex shedding we mean that the acoustic vortex shedding signature may rise only a few decibels above that for the broadband turbulent boundary layer acoustic sources). Correct prediction of weak vortex shedding may depend on accurately predicting the flow over the entire airfoil that includes the attached, turbulent upstream flow, adverse pressure gradient and separated flow regions and finally, the turbulent wake. This paper compares results of two full-LES and two LES with wall-stress model for the flow about a modified NACA 0016 airfoil with a 41° trailing edge apex angle and a slightly convex pressure side. Comparisons of vortex shedding, as measured by the power spectral density (PSD) of wall pressure fluctuations (WPF) on the pressure side of the TE and the PSD of the vertical velocity fluctuations in the wake are made. The results indicate that vortex shedding predictions are dependent upon the stream-wise and spanwise grid resolution. In order to reduce the large computational times required for simulating the high-Reynolds number flows with fully-resolved LES, a wall-stress model that solves the turbulent boundary layer equations in the near-wall region is applied. Compared with the fully-resolved LES, the LES with wall-stress simulations require about 20 percent the number of grid points and require about 10 percent of the computational time. However, the LES with wall stress model results under-predict the vortex shedding peak in the wake and are not able to predict the vortex shedding signature in TE wall pressure spectra. These results indicate that near-wall turbulence structures need to be resolved in order to correctly predict the occurence and strength of vortex shedding.


Author(s):  
Faustin Ondore

A square duct with a 90-degree streamwise curvature is representative of complex flow domains. Such flow domains are encountered in the designs of fluids engineering systems, especially in the aerospace turbo-machinery components. Examples include the gas turbine engine axial compressor interstage spaces, where the rise in air pressure (and hence compressor efficiency) is dependent on suppression of turbulence. In the case of the centrifugal compressor, pressure rise in the U-shaped diffuser assembly where the suppression of turbulence is critical to the attainable pressure ratio. The results obtained from numerical calculations are analysed and discussed along with the corresponding hot-wire measurements and flow visualization result from a wind-tunnel of identical configuration. Calculations are implemented in four turbulent models, i.e. Standard k-e Module, Algebraic Stress Model (ASM), Non-linear Renormalization Group (RNG) - k-e Model and Differential Stress Model (DSM). The discretization up-winding scheme is the Quadratic Up-winding with Interpolation Kinematics (QUICK). Two high Reynolds number turbulent flows are investigated, with mainstream velocities of 12.3 m/s and 20.4 m/s, representing Re=3.56×105 and Re=6.43×105 respectively. Generally strong correlation between theory and experimental data are recorded. Further, as reported in similar studies, the turbulence modules that are formulated to account for turbulence anisotropy return results that more closely match experimental measurements. Uniquely for this configuration, a massive flow detachment is predicted along the convex wall at about the 90° position. Also, the core of the fluid flow is observed to shift from the outer to the inner areas of the bend in proportion to the secondary (recirculating) flow generated by the bend.


Sign in / Sign up

Export Citation Format

Share Document