Turbulent boundary layers subject to severe acceleration or strong favorable pressure gradient (FPG) are of fundamental and technological importance. Scientifically, they elicit great interest from the points of view of scaling laws, the complex interaction between the outer and inner regions, and the quasi-laminarization phenomenon. Many flows of industrial and technological applications are subject to strong acceleration such as convergent ducts, turbines blades and nozzles. Our recent numerical predictions (J. Fluid Mech., vol. 775, pp. 189-200, 2015) of turbulent boundary layers subject to very strong FPG with high spatial/temporal resolution, i.e. Direct Numerical Simulation (DNS), have shown a meaningful weakening of the Reynolds shear stresses with an evident logarithmic behavior. In the present study, assessment of three different turbulence models (Shear Stress Transport, k-w and Spalart-Allmaras, henceforth SST, k-w and SA, respectively) in Reynolds-averaged Navier-Stokes (RANS) simulations is performed. The main objective is to evaluate the ability of popular turbulence models in capturing the characteristic features present during the quasi-laminarization phenomenon in highly accelerating turbulent boundary layers. Favorable pressure gradient is prescribed by a top converging surface (sink flow) with an approximately constant acceleration parameter of K = 4.0 x 10^(-6). Furthermore, the quasi-laminarization effect on the temperature field is also examined by solving the energy equation and assuming the temperature as a passive scalar. Validation of RANS results is carried out by means of a large DNS dataset.
more »
« less
Assessment of Turbulence Models over a Curved Hill Flow with Passive Scalar Transport
An incoming canonical spatially developing turbulent boundary layer (SDTBL) over a 2-D curved hill is numerically investigated via the Reynolds-averaged Navier–Stokes (RANS) equations plus two eddy-viscosity models: the K−ω SST (henceforth SST) and the Spalart–Allmaras (henceforth SA) turbulence models. A spatially evolving thermal boundary layer has also been included, assuming temperature as a passive scalar (Pr = 0.71) and a turbulent Prandtl number, Prt, of 0.90 for wall-normal turbulent heat flux modeling. The complex flow with a combined strong adverse/favorable streamline curvature-driven pressure gradient caused by concave/convex surface curvatures has been replicated from wind-tunnel experiments from the literature, and the measured velocity and pressure fields have been used for validation purposes (the thermal field was not experimentally measured). Furthermore, direct numerical simulation (DNS) databases from the literature were also employed for the incoming turbulent flow assessment. Concerning first-order statistics, the SA model demonstrated a better agreement with experiments where the turbulent boundary layer remained attached, for instance, in Cp, Cf, and Us predictions. Conversely, the SST model has shown a slightly better match with experiments over the flow separation zone (in terms of Cp and Cf) and in Us profiles just upstream of the bubble. The Reynolds analogy, based on the St/(Cf/2) ratio, holds in zero-pressure gradient (ZPG) zones; however, it is significantly deteriorated by the presence of streamline curvature-driven pressure gradient, particularly due to concave wall curvature or adverse-pressure gradient (APG). In terms of second-order statistics, the SST model has better captured the positively correlated characteristics of u′ and v′ or positive Reynolds shear stresses ( > 0) inside the recirculating zone. Very strong APG induced outer secondary peaks in and turbulence production as well as an evident negative slope on the constant shear layer.
more »
« less
- PAR ID:
- 10407177
- Date Published:
- Journal Name:
- Energies
- Volume:
- 15
- Issue:
- 16
- ISSN:
- 1996-1073
- Page Range / eLocation ID:
- 6013
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
In this video, we show high-fidelity numerical results of supersonic spatially-developing turbulent boundary layers (SDTBL) under strong concave and concave curvatures and Mach = 2.86. The selected numerical tool is Direct Numerical Simulation (DNS) with high spatial/temporal resolution. The prescribed concave geometry is based on the experimental study by Donovan et al. (J. Fluid Mech., 259, 1-24, 1994). Turbulent inflow conditions are based on extracted data from a previous DNS over a flat plate (i.e., turbulence precursors). The comprehensive DNS information sheds important light on the transport phenomena inside turbulent boundary layers subject to strong deceleration or Adverse Pressure Gradient (APG) caused by concave walls as well as to strong acceleration or Favorable Pressure Gradient (FPG) caused by convex walls at different wall thermal conditions (i.e., cold, adiabatic and hot walls). In this opportunity, the selected scientific visualization tool is Virtual Reality (VR) by extracting vortex core iso-surfaces via the Q-criterion to convert them to a file format readable by the HTC Vive VR toolkit. The reader is invited to visit our Virtual Wind Tunnel (VWT) under a fully immersive environment for further details. The video is available at: https://gfm.aps.org/meetings/dfd-2022/6313a60c199e4c2da9a946bcmore » « less
-
Recent numerical predictions of turbulent boundary layers subject to very strong Favorable Pressure Gradient (FPG) with high spatial/temporal resolution, i.e. Direct Numerical Simulation (DNS), have shown a meaningful weakening of the Reynolds shear stresses with a lengthy logarithmic behavior [1,2]. In the present study, assessment of the Shear Stress Transport and Spalart-Allmaras turbulence models (hence- forth SST and SA, respectively) in Reynolds-averaged Navier-Stokes (RANS) simulations is performed. The main objective is to evaluate the ability of popular turbulence models in capturing the characteristic features present during the quasi-laminarization phenomenon in highly accelerating turbulent boundary layers. A favorable pressure gradient is prescribed by a top converging surface (sink flow) with an approximately constant acceleration parameter of K = 4 . 0 ×10 −6 . Validation of RANS results is carried out by means of a large DNS dataset [1]. Generally speaking, the SA turbulence model has demonstrated the best compromise between accuracy and quick adaptation to the turbulent inflow conditions. Turbulence models properly captured the increasing trend of the freestream and friction velocity in highly accelerated flows; however, they fail to reproduce the decreasing behavior of the skin friction coefficient, which is typical in early stages of the quasi-laminarization process. Both models have shown deficient predictions of the decreasing and logarithmic behavior of Reynolds shear stresses as well as significantly overpredicted the production of Turbulent Kinetic Energy (TKE) in turbulent boundary layers subject to very strong FPG. https://doi.org/10.1016/j.compfluid.2020.104494more » « less
-
Direct Numerical Simulation (DNS) of compressible spatially-developing turbulent boundary layers (SDTBL) is performed at a Mach number of 2.5 and low/high Reynolds numbers over isothermal Zero-Pressure Gradient (ZPG) flat plates. Turbulent inflow information is generated via a dynamic rescaling-recycling approach (J. Fluid Mech., 670, pp. 581-605, 2011), which avoids the use of empirical correlations in the computation of inlet turbulent scales. The range of the low Reynolds number case is approximately 400-800, based on the momentum thickness, freestream velocity and wall viscosity. DNS at higher Reynolds numbers (~3,000, about four-fold larger) is also carried out with the purpose of analyzing the effect of Reynolds number on the transport phenomena in the supersonic regime. Additionally, low/high order flow statistics are compared with DNS of an incompressible isothermal ZPG boundary layer at similar low Reynolds numbers and the temperature regarded as a passive scalar. Peaks of turbulence intensities move closer to the wall as the Reynolds number increases in the supersonic flat plate. Furthermore, Reynolds shear stresses depict a much larger "plateau" (constant shear layer) at the highest Reynolds number considered in present study.more » « less
-
Direct numerical simulations (DNS) are performed to investigate the spatial evolution of flat-plate zero-pressure-gradient turbulent boundary layers over long streamwise domains ( $${>}300\delta _i$$ , with $$\delta _i$$ the inflow boundary-layer thickness) at three different Mach numbers, $2.5$ , $4.9$ and $10.9$ , with the surface temperatures ranging from quasiadiabatic to highly cooled conditions. The settlement of turbulence statistics into a fully developed equilibrium state of the turbulent boundary layer has been carefully monitored, either based on the satisfaction of the von Kármán integral equation or by comparing runs with different inflow turbulence generation techniques. The generated DNS database is used to characterize the streamwise evolution of multiple important variables in the high-Mach-number, cold-wall regime, including the skin friction, the Reynolds analogy factor, the shape factor, the Reynolds stresses, and the fluctuating wall quantities. The data confirm the validity of many classic and newer compressibility transformations at moderately high Reynolds numbers (up to friction Reynolds number $$Re_\tau \approx 1200$$ ) and show that, with proper scaling, the sizes of the near-wall streaks and superstructures are insensitive to the Mach number and wall cooling conditions. The strong wall cooling in the hypersonic cold-wall case is found to cause a significant increase in the size of the near-wall turbulence eddies (relative to the boundary-layer thickness), which leads to a reduced-scale separation between the large and small turbulence scales, and in turn to a lack of an outer peak in the spanwise spectra of the streamwise velocity in the logarithmic region.more » « less