scholarly journals Direct numerical simulation of ‘short’ laminar separation bubbles with turbulent reattachment

2000 ◽  
Vol 403 ◽  
pp. 223-250 ◽  
Author(s):  
M. ALAM ◽  
N. D. SANDHAM
Keyword(s):  
Numerical Simulation ◽  
Boundary Layer ◽  
Direct Numerical Simulation ◽  
Laminar Boundary Layer ◽  
Stokes Equations ◽  
Separation Bubble ◽  
Three Dimensional ◽  
Wall Turbulence ◽  
Mean Flow ◽  
Laminar Separation

Direct numerical simulation of the incompressible Navier-Stokes equations is used to study flows where laminar boundary-layer separation is followed by turbulent reattachment forming a closed region known as a laminar separation bubble. In the simulations a laminar boundary layer is forced to separate by the action of a suction profile applied as the upper boundary condition. The separated shear layer undergoes transition via oblique modes and Λ-vortex-induced breakdown and reattaches as turbulent flow, slowly recovering to an equilibrium turbulent boundary layer. Compared with classical experiments the computed bubbles may be classified as ‘short’, as the external potential flow is only affected in the immediate vicinity of the bubble. Near reattachment budgets of turbulence kinetic energy are dominated by turbulence events away from the wall. Characteristics of near-wall turbulence only develop several bubble lengths downstream of reattachment. Comparisons are made with two-dimensional simulations which fail to capture many of the detailed features of the full three-dimensional simulations. Stability characteristics of mean flow profiles are computed in the separated flow region for a family of velocity profiles generated using simulation data. Absolute instability is shown to require reverse flows of the order of 15–20%. The three-dimensional bubbles with turbulent reattachment have maximum reverse flows of less than 8% and it is concluded that for these bubbles the basic instability is convective in nature.

scholarly journals Direct numerical simulation of conical shock wave–turbulent boundary layer interaction

2019 ◽  
Vol 877 ◽  
pp. 167-195 ◽  
Author(s):  
Feng-Yuan Zuo ◽  
Antonio Memmolo ◽  
Guo-ping Huang ◽  
Sergio Pirozzoli
Keyword(s):  
Numerical Simulation ◽  
Boundary Layer ◽  
Shock Wave ◽  
Turbulent Boundary Layer ◽  
Direct Numerical Simulation ◽  
Pressure Gradient ◽  
Stokes Equations ◽  
Shear Stresses ◽  
Mean Flow ◽  
Conical Shock

Direct numerical simulation of the Navier–Stokes equations is carried out to investigate the interaction of a conical shock wave with a turbulent boundary layer developing over a flat plate at free-stream Mach number $M_{\infty }=2.05$ and Reynolds number $Re_{\unicode[STIX]{x1D703}}\approx 630$, based on the upstream boundary layer momentum thickness. The shock is generated by a circular cone with half opening angle $\unicode[STIX]{x1D703}_{c}=25^{\circ }$. As found in experiments, the wall pressure exhibits a distinctive N-wave signature, with a sharp peak right past the precursor shock generated at the cone apex, followed by an extended zone with favourable pressure gradient, and terminated by the trailing shock associated with recompression in the wake of the cone. The boundary layer behaviour is strongly affected by the imposed pressure gradient. Streaks are suppressed in adverse pressure gradient (APG) zones, but re-form rapidly in downstream favourable pressure gradient (FPG) zones. Three-dimensional mean flow separation is only observed in the first APG region associated with the formation of a horseshoe vortex, whereas the second APG region features an incipient detachment state, with scattered spots of instantaneous reversed flow. As found in canonical geometrically two-dimensional wedge-generated shock–boundary layer interactions, different amplification of the turbulent stress components is observed through the interacting shock system, with approach to an isotropic state in APG regions, and to a two-component anisotropic state in FPG. The general adequacy of the Boussinesq hypothesis is found to predict the spatial organization of the turbulent shear stresses, although different eddy viscosities should be used for each component, as in tensor eddy-viscosity models, or in full Reynolds stress closures.

Direct Numerical Simulation and Large Eddy Simulation of Laminar Separation Bubbles at Moderate Reynolds Numbers

2014 ◽  
Vol 136 (6) ◽  
Author(s):  
Francois Cadieux ◽  
Julian A. Domaradzki ◽  
Taraneh Sayadi ◽  
Sanjeeb Bose
Keyword(s):  
Heat Transfer ◽  
Numerical Simulation ◽  
Boundary Layer ◽  
Direct Numerical Simulation ◽  
Separation Bubble ◽  
Leading Edge ◽  
Transition To Turbulence ◽  
Laminar Separation Bubble ◽  
Laminar Separation ◽  
Large Eddy

Flows over airfoils and blades in rotating machinery for unmanned and microaerial vehicles, wind turbines, and propellers consist of different flow regimes. A laminar boundary layer near the leading edge is often followed by a laminar separation bubble with a shear layer on top of it that experiences transition to turbulence. The separated turbulent flow then reattaches and evolves downstream from a nonequilibrium turbulent boundary layer to an equilibrium one. Typical Reynolds-averaged Navier–Stokes (RANS) turbulence modeling methods were shown to be inadequate for such laminar separation bubble flows (Spalart and Strelets, 2000, “Mechanisms of Transition and Heat Transfer in a Separation Bubble,” J. Fluid Mech., 403, pp. 329–349). Direct numerical simulation (DNS) is the most reliable but is also the most computationally expensive alternative. This work assesses the capability of large eddy simulations (LES) to reduce the resolution requirements for such flows. Flow over a flat plate with suitable velocity boundary conditions away from the plate to produce a separation bubble is considered. Benchmark DNS data for this configuration are generated with the resolution of 59 × 106 mesh points; also used is a different DNS database with 15 × 106 points (Spalart and Strelets, 2000, “Mechanisms of Transition and Heat Transfer in a Separation Bubble,” J. Fluid Mech., 403, pp. 329–349). Results confirm that accurate LES are possible using O(1%) of the DNS resolution.

scholarly journals Compressible Direct Numerical Simulation of Low-Pressure Turbines: Part I — Methodology

2014 ◽  
Author(s):  
Richard D. Sandberg ◽  
Richard Pichler ◽  
Liwei Chen ◽  
Roderick Johnstone ◽  
Vittorio Michelassi
Keyword(s):  
Numerical Simulation ◽  
Boundary Layer ◽  
Direct Numerical Simulation ◽  
Environmental Parameters ◽  
Boundary Layer Separation ◽  
Low Pressure ◽  
Boundary Layer Transition ◽  
Layer Transition ◽  
Laminar Separation ◽  
Background Turbulence

Modern low pressure turbines (LPT) feature high pressure ratios and moderate Mach and Reynolds numbers, increasing the possibility of laminar boundary-layer separation on the blades. Upstream disturbances including background turbulence and incoming wakes have a profound effect on the behavior of separation bubbles and the type/location of laminar-turbulent transition and therefore need to be considered in LPT design. URANS are often found inadequate to resolve the complex wake dynamics and impact of these environmental parameters on the boundary layers and may not drive the design to the best aerodynamic efficiency. LES can partly improve the accuracy, but has difficulties in predicting boundary layer transition and capturing the delay of laminar separation with varying inlet turbulence levels. Direct Numerical Simulation (DNS) is able to overcome these limitations but has to date been considered too computationally expensive. Here a novel compressible DNS code is presented and validated, promising to make DNS practical for LPT studies. Also, the sensitivity of wake loss coefficient with respect to freestream turbulence levels below 1% is discussed.

scholarly journals Self-similar compressible turbulent boundary layers with pressure gradients. Part 1. Direct numerical simulation and assessment of Morkovin’s hypothesis

2019 ◽  
Vol 880 ◽  
pp. 239-283 ◽  
Author(s):  
Christoph Wenzel ◽  
Tobias Gibis ◽  
Markus Kloker ◽  
Ulrich Rist
Keyword(s):  
Numerical Simulation ◽  
Boundary Layer ◽  
Mach Number ◽  
Direct Numerical Simulation ◽  
Boundary Layers ◽  
Turbulent Boundary Layers ◽  
Shear Stresses ◽  
Pressure Gradients ◽  
Mean Flow ◽  
Self Similar

A direct numerical simulation study of self-similar compressible flat-plate turbulent boundary layers (TBLs) with pressure gradients (PGs) has been performed for inflow Mach numbers of 0.5 and 2.0. All cases are computed with smooth PGs for both favourable and adverse PG distributions (FPG, APG) and thus are akin to experiments using a reflected-wave set-up. The equilibrium character allows for a systematic comparison between sub- and supersonic cases, enabling the isolation of pure PG effects from Mach-number effects and thus an investigation of the validity of common compressibility transformations for compressible PG TBLs. It turned out that the kinematic Rotta–Clauser parameter $\unicode[STIX]{x1D6FD}_{K}$ calculated using the incompressible form of the boundary-layer displacement thickness as length scale is the appropriate similarity parameter to compare both sub- and supersonic cases. Whereas the subsonic APG cases show trends known from incompressible flow, the interpretation of the supersonic PG cases is intricate. Both sub- and supersonic regions exist in the boundary layer, which counteract in their spatial evolution. The boundary-layer thickness $\unicode[STIX]{x1D6FF}_{99}$ and the skin-friction coefficient $c_{f}$, for instance, are therefore in a comparable range for all compressible APG cases. The evaluation of local non-dimensionalized total and turbulent shear stresses shows an almost identical behaviour for both sub- and supersonic cases characterized by similar $\unicode[STIX]{x1D6FD}_{K}$, which indicates the (approximate) validity of Morkovin’s scaling/hypothesis also for compressible PG TBLs. Likewise, the local non-dimensionalized distributions of the mean-flow pressure and the pressure fluctuations are virtually invariant to the local Mach number for same $\unicode[STIX]{x1D6FD}_{K}$-cases. In the inner layer, the van Driest transformation collapses compressible mean-flow data of the streamwise velocity component well into their nearly incompressible counterparts with the same $\unicode[STIX]{x1D6FD}_{K}$. However, noticeable differences can be observed in the wake region of the velocity profiles, depending on the strength of the PG. For both sub- and supersonic cases the recovery factor was found to be significantly decreased by APGs and increased by FPGs, but also to remain virtually constant in regions of approximated equilibrium.

Sign in / Sign up

Export Citation Format

Share Document