Direct numerical simulation of transonic shock/boundary layer interaction under conditions of incipient separation

2010 ◽  
Vol 657 ◽  
pp. 361-393 ◽  
Author(s):  
SERGIO PIROZZOLI ◽  
MATTEO BERNARDINI ◽  
FRANCESCO GRASSO
Keyword(s):  
Numerical Simulation ◽  
Boundary Layer ◽  
Direct Numerical Simulation ◽  
Normal Shock ◽  
Local Regularity ◽  
Numerical Dissipation ◽  
Similar Distribution ◽  
Interaction Zone ◽  
Mean Flow ◽  
The Mean

The interaction of a normal shock wave with a turbulent boundary layer developing over a flat plate at free-stream Mach number M∞ = 1.3 and Reynolds number Reθ ≈ 1200 (based on the momentum thickness of the upstream boundary layer) is analysed by means of direct numerical simulation of the compressible Navier–Stokes equations. The computational methodology is based on a hybrid linear/weighted essentially non-oscillatory conservative finite-difference approach, whereby the switch is controlled by the local regularity of the solution, so as to minimize numerical dissipation. As found in experiments, the mean flow pattern consists of an upstream fan of compression waves associated with the thickening of the boundary layer, and the supersonic region is terminated by a nearly normal shock, with substantial bending of the interacting shock. At the selected conditions the flow does not exhibit separation in the mean. However, the interaction region is characterized by ‘intermittent transitory detachment’ with scattered spots of instantaneous flow reversal throughout the interaction zone, and by the formation of a turbulent mixing layer, with associated unsteady release of vortical structures. As found in supersonic impinging shock interactions, we observe a different amplification of the longitudinal Reynolds stress component with respect to the others. Indeed, the effect of the adverse pressure gradient is to reduce the mean shear, with subsequent suppression of the near-wall streaks, and isotropization of turbulence. The recovery of the boundary layer past the interaction zone follows a quasi-equilibrium process, characterized by a self-similar distribution of the mean flow properties.

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.

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.

Direct Numerical Simulation of By-Pass and Separation-Induced Transition in a Linear Compressor Cascade

2006 ◽  
Author(s):  
Tamer Zaki ◽  
Paul Durbin ◽  
Jan Wissink ◽  
Wolfgang Rodi
Keyword(s):  
Numerical Simulation ◽  
Boundary Layer ◽  
Direct Numerical Simulation ◽  
Pressure Gradient ◽  
Surface Boundary ◽  
Mean Flow ◽  
Compressor Cascade ◽  
Linear Compressor ◽  
Free Stream Turbulence ◽  
Linear Compressor Cascade

Direct Numerical Simulation (DNS) of flow through a linear compressor cascade with incoming free-stream turbulence was performed. On the pressure side, the boundary layer flow is found to undergo by-pass transition: The incident vortical disturbances trigger the formation of elongated boundary layer perturbation jets (or streaks) with amplitudes on the order of 10% of the mean flow. The inception of turbulent spots, which leads to breakdown, is triggered on the backward perturbation jets (negative u-fluctuations). The turbulent patches spread and finally merge into the downstream, fully turbulent region. The suction surface boundary layer is initially subject to a Favorable Pressure Gradient (FPG), followed by a strong Adverse Pressure Gradient (APG). The FPG suppresses the formation of boundary layer streaks. The result is a stabilized boundary layer that does not undergo transition. Farther downstream, the strong APG causes the laminar boundary layer to separate, which is followed by turbulent reattachment.

Direct Numerical Simulation on Transition of an Incompressible Boundary Layer on a Flat Plate

2012 ◽  
Vol 268-270 ◽  
pp. 1143-1147
Author(s):  
Ning Li ◽  
Qi Hong Zeng
Keyword(s):  
Numerical Simulation ◽  
Boundary Layer ◽  
Direct Numerical Simulation ◽  
Flat Plate ◽  
Flow Profile ◽  
Mean Flow ◽  
Turbulent Transition ◽  
Incompressible Boundary Layer ◽  
Incompressible Boundary ◽  
Laminar Turbulent Transition

Direct Numerical Simulation(DNS) was carried out for laminar-turbulent transition of an incompressible boundary layer on a flat plate based on disturbance Navier-Stokes(N-S) equation in spatial mode with Massage Passing Interface(MPI) technology. Study on breakdown mechanism of laminar-turbulent transition was carried on. The effect of mean flow distortion on the process of breakdown in laminar-turbulent transition was investigated. Results indicate that change of instability characteristic of mean flow profile plays a key role during process of breakdown.

Direct numerical simulation of turbulent duct flow with inclined or V-shaped ribs mounted on one wall

2021 ◽  
Vol 932 ◽  
Author(s):  
S.V. Mahmoodi-Jezeh ◽  
Bing-Chen Wang
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Secondary Flows ◽  
Wall Turbulence ◽  
Spanwise Direction ◽  
Duct Flow ◽  
Mean Flow ◽  
Square Duct ◽  
Turbulence Structures ◽  
The Mean

In this research, highly disturbed turbulent flow of distinct three-dimensional characteristics in a square duct with inclined or V-shaped ribs mounted on one wall is investigated using direct numerical simulation. The turbulence field is highly sensitive to not only the rib geometry but also the boundary layers developed over the side and top walls. In a cross-stream plane secondary flows appear as large longitudinal vortices in both inclined and V-shaped rib cases due to the confinement of four sidewalls of the square duct. However, owing to the difference in the pattern of cross-stream secondary flow motions, the flow physics is significantly different in these two ribbed duct cases. It is observed that the mean flow structures in the cross-stream directions are asymmetrical in the inclined rib case but symmetrical in the V-shaped rib case, causing substantial differences in the momentum transfer across the spanwise direction. The impacts of rib geometry on near-wall turbulence structures are investigated using vortex identifiers, joint probability density functions between the streamwise and vertical velocity fluctuations, statistical moments of different orders, spatial two-point autocorrelations and velocity spectra. It is found that near the leeward and windward rib faces, the mean inclination angle of turbulence structures in the V-shaped rib case is greater than that of the inclined rib case, which subsequently enhances momentum transport between the ribbed bottom wall and the smooth top wall.

scholarly journals Direct numerical simulation of a self-similar adverse pressure gradient turbulent boundary layer at the verge of separation

2017 ◽  
Vol 829 ◽  
pp. 392-419 ◽  
Author(s):  
V. Kitsios ◽  
A. Sekimoto ◽  
C. Atkinson ◽  
J. A. Sillero ◽  
G. Borrell ◽  
...  
Keyword(s):  
Numerical Simulation ◽  
Boundary Layer ◽  
Reynolds Number ◽  
Direct Numerical Simulation ◽  
Pressure Gradient ◽  
Adverse Pressure Gradient ◽  
Streamwise Velocity ◽  
Far Field ◽  
Self Similar ◽  
The Mean

The statistical properties are presented for the direct numerical simulation of a self-similar adverse pressure gradient (APG) turbulent boundary layer (TBL) at the verge of separation. The APG TBL has a momentum thickness-based Reynolds number range from $Re_{\unicode[STIX]{x1D6FF}_{2}}=570$ to 13 800, with a self-similar region from $Re_{\unicode[STIX]{x1D6FF}_{2}}=10\,000$ to 12 300. Within this domain the average non-dimensional pressure gradient parameter $\unicode[STIX]{x1D6FD}=39$, where for a unit density $\unicode[STIX]{x1D6FD}=\unicode[STIX]{x1D6FF}_{1}P_{\!e}^{\prime }/\unicode[STIX]{x1D70F}_{w}$, with $\unicode[STIX]{x1D6FF}_{1}$ the displacement thickness, $\unicode[STIX]{x1D70F}_{w}$ the mean shear stress at the wall and $P_{\!e}^{\prime }$ the far-field pressure gradient. This flow is compared with previous zero pressure gradient and mild APG TBL ($\unicode[STIX]{x1D6FD}=1$) results of similar Reynolds number. All flows are generated via the direct numerical simulation of a TBL on a flat surface with far-field boundary conditions tailored to apply the desired pressure gradient. The conditions for self-similarity, and the appropriate length and velocity scales, are derived. The mean and Reynolds stress profiles are shown to collapse when non-dimensionalised on the basis of these length and velocity scales. As the pressure gradient increases, the extent of the wake region in the mean streamwise velocity profiles increases, whilst the extent of the log-layer and viscous sublayer decreases. The Reynolds stress, production and dissipation profiles of the APG TBL cases exhibit a second outer peak, which becomes more pronounced and more spatially localised with increasing pressure gradient. This outer peak is located at the point of inflection of the mean velocity profiles, and is suggestive of the presence of a shear flow instability. The maximum streamwise velocity variance is located at a wall normal position of $\unicode[STIX]{x1D6FF}_{1}$ of spanwise wavelength of $2\unicode[STIX]{x1D6FF}_{1}$. In summary as the pressure gradient increases the flow has properties less like a zero pressure gradient TBL and more akin to a free shear layer.

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.

Data Analysis Methods for Stereo PIV of a High Aspect Ratio Flexible Cylinder

2007 ◽  
Author(s):  
D. Furey ◽  
P. Atsavapranee ◽  
K. Cipolla
Keyword(s):  
Boundary Layer ◽  
Data Analysis ◽  
Aspect Ratio ◽  
High Aspect Ratio ◽  
Particle Image Velocimetry Data ◽  
Flow Fields ◽  
Mean Flow ◽  
Boundary Flow ◽  
The Mean ◽  
Cylinder Flow

Stereo Particle Image velocimetry data was collected over high aspect ratio flexible cylinders (L/a = 1.5 to 3 × 105) to evaluate the axial development of the turbulent boundary layer where the boundary layer thickness becomes significantly larger than the cylinder diameter (δ/a>>1). The flexible cylinders are approximately neutrally buoyant and have an initial length of 152 m and radii of 0.45 mm and 1.25 mm. The cylinders were towed at speeds ranging from 3.8 to 15.4 m/sec in the David Taylor Model Basin. The analysis of the SPIV data required a several step procedure to evaluate the cylinder boundary flow. First, the characterization of the flow field from the towing strut is required. This evaluation provides the residual mean velocities and turbulence levels caused by the towing hardware at each speed and axial location. These values, called tare values, are necessary for comparing to the cylinder flow results. Second, the cylinder flow fields are averaged together and the averaged tare fields are subtracted out to remove strut-induced ambient flow effects. Prior to averaging, the cylinder flow fields are shifted to collocate the cylinder within the field. Since the boundary layer develops slowly, all planes of data occurring within each 10 meter increment of the cylinder length are averaged together to produce the mean boundary layer flow. Corresponding fields from multiple runs executed using the same experimental parameters are also averaged. This flow is analyzed to evaluate the level of axisymmetry in the data and determine if small changes in cylinder angle affect the mean flow development. With axisymmetry verified, the boundary flow is further averaged azimuthally around the cylinder to produce mean boundary layer profiles. Finally, the fluctuating velocity levels are evaluated for the flow with the cylinder and compared to the fluctuating velocity levels in the tare data. This paper will first discuss the data analysis techniques for the tare data and the averaging methods implemented. Second, the data analysis considerations will be presented for the cylinder data and the averaging and cylinder tracking techniques. These results are used to extract relevant boundary layer parameters including δ, δ* and θ. Combining these results with wall shear and momentum thickness values extracted from averaged cylinder drag data, the boundary layer can be well characterized.

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.

Sign in / Sign up

Export Citation Format

Share Document