Prediction of Airfoil Stall Based on a Modified k−v2¯−ω Turbulence Model

The accuracy of an airfoil stall prediction heavily depends on the computation of the separated shear layer. Capturing the strong non-equilibrium turbulence in the shear layer is crucial for the accuracy of a stall prediction. In this paper, different Reynolds-averaged Navier–Stokes turbulence models are adopted and compared for airfoil stall prediction. The results show that the separated shear layer fixed k−v2¯−ω (abbreviated as SPF k−v2¯−ω) turbulence model captures the non-equilibrium turbulence in the separated shear layer well and gives satisfactory predictions of both thin-airfoil stall and trailing-edge stall. At small Reynolds numbers (Re~105), the relative error between the predicted CL,max of NACA64A010 by the SPF k−v2¯−ω model and the experimental data is less than 3.5%. At high Reynolds numbers (Re~106), the CL,max of NACA64A010 and NACA64A006 predicted by the SPF k−v2¯−ω model also has an error of less than 5.5% relative to the experimental data. The stall of the NACA0012 airfoil, which features trailing-edge stall, is also computed by the SPF k−v2¯−ω model. The SPF k−v2¯−ω model is also applied to a NACA0012 airfoil, which features trailing-edge stall and an error of CL relative to the experiment at CL>1.0 is smaller than 3.5%. The SPF k−v2¯−ω model shows higher accuracy than other turbulence models.

Large Eddy simulation of tracer gas dispersion in a cavity

This paper assesses the prediction of inert tracer gas dispersion within a cavity of height (H) 1.0 m, and unity aspect ratio, using large Eddy simulation (LES). The flow Reynolds number was 67 000, based on the freestream velocity and cavity height. The flow upstream of the cavity was laminar, producing a cavity shear layer which underwent a transition to turbulence over the cavity. Three distinct meshes are used, with grid spacings of H / 100 (coarse), H / 200 (intermediate), and H / 400 (fine) respectively. The Smagorinsky, WALE, and Germano-Lilly subgrid-scale models are used on each grid to quantify the effects of subgrid-scale modelling on the simulated flow. Coarsening the grid led to small changes in the predicted velocity field, and to substantial over-prediction of the tracer gas concentration statistics. Quantitative metric analysis of the tracer gas statistics showed that the coarse grid simulations yielded results outside of acceptable tolerances, while the intermediate and fine grids produced acceptable output. Interrogation of the fluid dynamics present in each simulation showed that the evolution of the cavity shear layer is heavily influenced by the grid and subgrid scale model. On the coarse and intermediate grids the development of the shear layer is delayed, inhibiting the entrainment and mixing of the tracer gas into the shear layer, reducing the removal of the tracer gas from the cavity. On the fine grid, the shear layer developed more rapidly, resulting in enhanced removal of the tracer gas from the cavity. Concentration probability density functions showed that the fine grid simulations accurately predicted the range, and the most probable value, of the tracer gas concentration towards both walls of the cavity. The results presented in this paper show that the WALE and Germano-Lilly models may be advantageous over the standard Smagorinsky model for simulations of pollutant dispersion in the urban environment.

Measurements of Turbulent Transport in a Square Channel with One Ribbed Wall

The velocity field of the fully developed turbulent flow in a one-sided ribbed square channel (rib-height-to-channel-height ratio of k/h = 0.0667, rib-pitch-to-rib-height ratio of p/k = 9) were measured at Reynolds numbers (based on the channel height h and the mean bulk velocity uB) of Reh = 50 000 and 100 000 by means of Laser-Doppler-Anemometry (LDA). Triple velocity correlations differed slightly between both Reynolds numbers when normalized by the bulk velocity and the channel height, similarly to the first- and second-order statistical moments of the velocity. Their near-wall behavior reflected the crucial role of turbulent transport near the rib crest and within the separated shear layer. Sweep events occurred with the elongated flow structures of the flapping shear layer and gained in importance towards the channel bottom wall, while strong ejection events near the rib leading and trailing edges coincided with flow structures bursting away from the wall. Despite the predominant occurrence of sweep events close to the ribbed wall within the inter-rib spacing, ejection events contributed with higher intensity to the Reynolds shear stress. Ejection and sweep events and their underlying transport phenomena contributing to the Reynolds shear stress were almost Reynolds number-insensitive in the resolved flow range. The invariance to the Reynolds number can be of benefit for the use of scale-resolving simulation methods in the design process of rib structures for heat exchange applications.

Second Moment Closure Modeling and DNS of Stratified Shear Layers

Buoyant shear layers encountered in many engineering and environmental applications have been studied by researchers for decades. Often, these flows have high Reynolds and Richardson numbers, which leads to significant/intractable space-time resolution requirements for DNS or LES. On the other hand, many of the important physical mechanisms, such as stress anisotropy, wake stabilization, and regime transition, inherently render eddy viscosity-based RANS modeling inappropriate. Accordingly, we pursue second-moment closure (SMC), i.e., full Reynolds stress/flux/variance modeling, for moderate Reynolds number non-stratified, and stratified shear layers for which DNS is possible. A range of sub-model complexity is pursued for the diffusion of stresses, density fluxes and variance, pressure strain and scrambling, and dissipation. These sub-models are evaluated in terms of how well they are represented by DNS in comparison to the exact Reynolds averaged terms, and how well they impact the accuracy of full RANS closure. For the non-stratified case, SMC model predicts the shear layer growth rate and Reynolds shear stress profiles accurately. Stress anisotropy and budgets are captured only qualitatively. Comparing DNS of exact and modeled terms, inconsistencies in model performance and assumptions are observed, including inaccurate prediction of individual statistics, non-negligible pressure diffusion, and dissipation anisotropy. For the stratified case, shear layer and gradient Richardson number growth rates, and stress, flux and variance decay rates, are captured with less accuracy than corresponding flow parameters in the non-stratified case. These studies lead to several recommendations for model improvement.

Changes in the Near-surface Shear Layer of the Sun

We use helioseismic data obtained over two solar cycles to determine whether there are changes in the near-surface shear layer (NSSL). We examine this by determining the radial gradient of the solar rotation rate. The radial gradient itself shows a solar-cycle dependence, and the changes are more pronounced in the active latitudes than at adjoining higher latitudes; results at the highest latitudes (≳70°) are unreliable. The pattern changes with depth, even within the NSSL. We find that the near-surface shear layer is deeper at lower latitudes than at high latitudes and that the extent of the layer also shows a small solar-cycle-related change.

Vortex evolution in a rotating tank with an off-axis drain

Fluid entering the periphery of a steadily rotating cylindrical tank exits through an off-axis drain hole, located in the tank's base at the half-radius. Experiments show that, though a concentrated vortex forms over the drain, it soon advects around the tank in what is at first a circular path. Though inviscid vortex dynamics predicts continued motion, our experiments show that the vortex moves inwards from the predicted circular path, finally coming to rest at approximately $50^{\circ }$ from the drain. In this final state, the vorticity is concentrated in a thin shear layer bounding an irrotational core, which passes over the drain. The broadening of the vortex structure and eventual steady-state formation are believed to be due to the growing boundary layer on the outer wall.