Prediction and investigation of the turbulent flow over a rotating disk

2000 ◽  
Vol 418 ◽  
pp. 231-264 ◽  
Author(s):  
XIAOHUA WU ◽  
KYLE D. SQUIRES

Large-eddy simulation (LES) has been used to predict the statistically three-dimensional turbulent boundary layer (3DTBL) over a rotating disk. LES predictions for six parameter cases were compared to the experimental measurements of Littell & Eaton (1994), obtained at a momentum thickness Reynolds number of 2660. A signal-decomposition scheme was developed by modifying the method of Spalart (1988) to prescribe time-dependent boundary conditions along the radial direction, entrainment towards the disk surface was prescribed by satisfying global mass conservation. Predictions of the mean velocities and r.m.s. fluctuations are in good agreement with data, with the largest discrepancy occurring in the prediction of the wall-normal intensities. The primary and two secondary shear stresses are also in good agreement with the measurements and one-dimensional energy spectra of the velocity fluctuations agree well with established laws, i.e. a −1 slope in the buffer region and −5/3 slope near the edge of the boundary layer.Conditionally averaged velocities provide new evidence in support of the structural model of Littell & Eaton (1994) concerning the interaction of mean-flow three-dimensionality and shear-stress producing structures. Inside the buffer region under strong ejections, the conditionally averaged crossflow (radial) velocity is larger than the unconditioned mean, and the profile conditioned on strong sweeps is smaller than the mean. This is consistent with the notion that streamwise vortices having the same sign as the mean streamwise vorticity, and beneath the peak crossflow location, are mostly responsible for strong sweep events; streamwise vortices with opposite sign as the mean streamwise vorticity promote strong ejections. Comparison of two-point spatial correlations with previous measurements in two-dimensional turbulent boundary layers (2DTBLs) indicates interesting structural similarities, e.g. the correlation of wall pressure and surface-normal velocity fluctuations is an odd function of streamwise separation, being positive downstream and negative upstream. These similarities offer quantitative indirect support to the hypothesis advanced by Littell & Eaton (1994) and Johnston & Flack (1996) that structural models describing 2DTBLs may be employed as a baseline in (equilibrium) 3DTBL structural studies.

1975 ◽  
Vol 42 (3) ◽  
pp. 591-597 ◽  
Author(s):  
D. H. Wood ◽  
R. A. Antonia

Mean velocity and turbulence intensity measurements have been made in a fully developed turbulent boundary layer over a d-type surface roughness. This roughness is characterised by regular two-dimensional elements of square cross section placed one element width apart, with the cavity flow between elements being essentially isolated from the outer flow. The measurements show that this boundary layer closely satisfies the requirement of exact self-preservation. Distribution across the layer of Reynolds normal and shear stresses are closely similar to those found over a smooth surface except for the region immediately above the grooves. This similarity extends to distributions of third and fourth-order moments of longitudinal and normal velocity fluctuations and also to the distribution of turbulent energy dissipation. The present results are compared with those obtained for a k-type or sand grained roughness.


2018 ◽  
Vol 849 ◽  
pp. 192-230 ◽  
Author(s):  
Dominic A. van der A ◽  
Pietro Scandura ◽  
Tom O’Donoghue

Turbulence characteristics of an asymmetric oscillatory boundary layer flow are analysed through two-component laser-Doppler measurements carried out in a large oscillatory flow tunnel and direct numerical simulation (DNS). Five different Reynolds numbers, $R_{\unicode[STIX]{x1D6FF}}$, in the range 846–2057 have been investigated experimentally, where $R_{\unicode[STIX]{x1D6FF}}=\tilde{u} _{0max}\unicode[STIX]{x1D6FF}/\unicode[STIX]{x1D708}$ with $\tilde{u} _{0max}$ the maximum oscillatory velocity in the irrotational region, $\unicode[STIX]{x1D6FF}$ the Stokes length and $\unicode[STIX]{x1D708}$ the fluid kinematic viscosity. DNS has been carried out for the lowest three $R_{\unicode[STIX]{x1D6FF}}$ equal to 846, 1155 and 1475. Both experimental and numerical results show that the flow statistics increase during accelerating phases of the flow and especially at times of transition to turbulent flow. Once turbulence is fully developed, the near-wall statistics remain almost constant until the late half-cycle, with values close to those reported for steady wall-bounded flows. The higher-order statistics reach large values within a normalized wall distance of approximately $y/\unicode[STIX]{x1D6FF}=0.2$ at phases corresponding to the onset of low-speed streak breaking, because of the intermittency of the velocity fluctuations at these times. In particular, the flatness of the streamwise velocity fluctuations reaches values of the order of ten, while the flatness of the wall-normal velocity fluctuations reaches values of several hundreds. Far from the wall, at locations where the vertical gradient of the streamwise velocity is zero, the skewness is approximately zero and the flatness is approximately equal to 3, representative of a normal distribution. At lower elevations the distribution of the fluctuations deviate substantially from a normal distribution, but are found to be well described by other standard theoretical probability distributions.


1996 ◽  
Vol 326 ◽  
pp. 151-179 ◽  
Author(s):  
Junhui Liu ◽  
Ugo Piomelli ◽  
Philippe R. Spalart

The interaction between a zero-pressure-gradient turbulent boundary layer and a pair of strong, common-flow-down, streamwise vortices with a sizeable velocity deficit is studied by large-eddy simulation. The subgrid-scale stresses are modelled by a localized dynamic eddy-viscosity model. The results agree well with experimental data. The vortices drastically distort the boundary layer, and produce large spanwise variations of the skin friction. The Reynolds stresses are highly three-dimensional. High levels of kinetic energy are found both in the upwash region and in the vortex core. The two secondary shear stresses are significant in the vortex region, with magnitudes comparable to the primary one. Turbulent transport from the immediate upwash region is partly responsible for the high levels of turbulent kinetic energy in the vortex core; its effect on the primary stress 〈u′v′〉 is less significant. The mean velocity gradients play an important role in the generation of 〈u′v′〉 in all regions, while they are negligible in the generation of turbulent kinetic energy in the vortex core. The pressure-strain correlations are generally of opposite sign to the production terms except in the vortex core, where they have the same sign as the production term in the budget of 〈u′v′〉. The results highlight the limitations of the eddy-viscosity assumption (in a Reynolds-averaged context) for flows of this type, as well as the excessive diffusion predicted by typical turbulence models.


1994 ◽  
Vol 116 (4) ◽  
pp. 586-596 ◽  
Author(s):  
P. L. Andrew ◽  
Wing-fai Ng

The turbulent character of the supersonic wake of a linear cascade of fan airfoils has been studied using a two-component laser-doppler anemometer. The cascade was tested in the Virginia Polytechnic Institute and State University intermittent wind tunnel facility, where the Mach and Reynolds numbers were 2.36 and 4.8 × 106, respectively. In addition to mean flow measurements, Reynolds normal and shear stresses were measured as functions of cascade incidence angle and streamwise locations spanning the near-wake and the far-wake. The extremities of profiles of both the mean and turbulent wake properties´ were found to be strongly influenced by upstream shock-boundary -layer interactions, the strength of which varied with cascade incidence. In contrast, the peak levels of turbulence properties within the shear layer were found to be largely independent of incidence, and could be characterized in terms of the streamwise position only. The velocity defect turbulence level was found to be 23 percent, and the generally accepted value of the turbulence structural coefficient of 0.30 was found to be valid for this flow. The degree of similarity of the mean flow wake profiles was established, and those profiles demonstrating the most similarity were found to approach a state of equilibrium between the mean and turbulent properties. In general, this wake flow may be described as a classical free shear flow, upon which the influence of upstream shock-boundary-layer interactions has been superimposed.


1994 ◽  
Vol 266 ◽  
pp. 175-207 ◽  
Author(s):  
Howard S. Littell ◽  
John K. Eaton

Measurements of the boundary layer on an effectively infinite rotating disk in a quiescent environment are described for Reynolds numbers up to Reδ2 = 6000. The mean flow properties were found to resemble a ‘typical’ three-dimensional crossflow, while some aspects of the turbulence measurements were significantly different from two-dimensional boundary layers that are turned. Notably, the ratio of the shear stress vector magnitude to the turbulent kinetic energy was found to be at a maximum near the wall, instead of being locally depressed as in a turned two-dimensional boundary layer. Also, the shear stress and the mean strain rate vectors were found to be more closely aligned than would be expected in a flow with this degree of crossflow. Two-point velocity correlation measurements exhibited strong asymmetries which are impossible in a two-dimensional boundary layer. Using conditional sampling, the velocity field surrounding strong Reynolds stress events was partially mapped. These data were studied in the light of the structural model of Robinson (1991), and a hypothesis describing the effect of cross-stream shear on Reynolds stress events is developed.


2000 ◽  
Vol 420 ◽  
pp. 47-83 ◽  
Author(s):  
NIKOLAUS A. ADAMS

The turbulent boundary layer along a compression ramp with a deflection angle of 18° at a free-stream Mach number of M = 3 and a Reynolds number of Reθ = 1685 with respect to free-stream quantities and mean momentum thickness at inflow is studied by direct numerical simulation. The conservation equations for mass, momentum, and energy are solved in generalized coordinates using a 5th-order hybrid compact- finite-difference-ENO scheme for the spatial discretization of the convective fluxes and 6th-order central compact finite differences for the diffusive fluxes. For time advancement a 3rd-order Runge–Kutta scheme is used. The computational domain is discretized with about 15 × 106 grid points. Turbulent inflow data are provided by a separate zero-pressure-gradient boundary-layer simulation. For statistical analysis, the flow is sampled 600 times over about 385 characteristic timescales δ0/U∞, defined by the mean boundary-layer thickness at inflow and the free-stream velocity. Diagnostics show that the numerical representation of the flow field is sufficiently well resolved.Near the corner, a small area of separated flow develops. The shock motion is limited to less than about 10% of the mean boundary-layer thickness. The shock oscillates slightly around its mean location with a frequency of similar magnitude to the bursting frequency of the incoming boundary layer. Turbulent fluctuations are significantly amplified owing to the shock–boundary-layer interaction. Reynolds-stress maxima are amplified by a factor of about 4. Turbulent normal and shear stresses are amplified differently, resulting in a change of the structure parameter. Compressibility affects the turbulence structure in the interaction area around the corner and during the relaxation after reattachment downstream of the corner. Correlations involving pressure fluctuations are significantly enhanced in these regions. The strong Reynolds analogy which suggests a perfect correlation between velocity and temperature fluctuations is found to be invalid in the interaction area.


Author(s):  
Witold Elsner ◽  
Stephane Vilmin ◽  
Stanislaw Drobniak ◽  
Wladyslaw Piotrowski

The paper presents an experimental and numerical analysis of the interaction between wakes and boundary layers on aerodynamic blade profiles. The experiment revealed that incoming wakes interact with boundary layers and cause the significant increase of velocity fluctuations in the boundary layer and in consequence shift the transition zone towards the leading edge. The full time evolution of periodic wake induced transition was reproduced from measurements. The numerical simulation of the flow around the blade profile has been performed with the use of the adaptive grid viscous flow unNEWT PUIM solver with a prescribed unsteady intermittency method (PUIM) developed at Cambridge University, UK. The results obtained give evidence that the turbulence transported within the wake is mainly responsible for the transition process. The applied CFD solver was able to reproduce some essential flow features related to the bypass and wake-induced transitions and the simulations reveal good agreement with the experimental results in terms of localisation and extent of wake-induced transition.


Author(s):  
Barbara L. da Silva ◽  
David Sumner ◽  
Don Bergstrom

Abstract The flow downstream of surface-mounted finite-height square prisms with aspect ratio AR = 1 (cube) and 0.5 (block) was investigated experimentally in a low-speed wind tunnel, to determine the overall structure and dynamic behavior of the wake and the source of the streamwise vorticity. The Reynolds number based on the prisms' width D was Re = 7.5×104 and the boundary layer thickness at the location of the prisms was d/D = 0.73. A vortex shedding frequency was found in the wake of the cube, but no periodicity was found in the wake of the block. The mean wake of the cube showed features of prisms below the critical AR, but the wake of the block had a distinct behavior due to the dominant shear flow from the boundary layer. The shear changed the downwash and, consequently, the streamwise vorticity distribution in the wake, in addition to reducing the magnitude of the Reynolds stresses. The phase-average analysis for the cube revealed the alternate shedding of inclined structures related to the streamwise vorticity in the upper part of the wake. These vorticity regions were caused by the alternate bending and entrainment of the side flow, caused by the downwash. The periodic component of the total Reynolds stresses was, however, significantly smaller than the turbulence-related stresses. The present study showed that the wake had a transitional behavior for the cube, but became fundamentally different for the block, when compared with prisms of higher AR.


2011 ◽  
Vol 688 ◽  
pp. 258-296 ◽  
Author(s):  
I. Jacobi ◽  
B. J. McKeon

AbstractThe zero-pressure-gradient turbulent boundary layer over a flat plate was perturbed by a temporally oscillating, spatial impulse of roughness, and the downstream response of the flow field was interrogated by hot-wire anemometry and particle-image velocimetry. The key features common to impulsively perturbed boundary layers, as identified in Jacobi & McKeon (J. Fluid Mech., 2011), were investigated, and the unique contributions of the dynamic perturbation were isolated by contrast with an appropriately matched static impulse of roughness. In addition, the dynamic perturbation was decomposed into separable large-scale and small-scale structural effects, which in turn were associated with the organized wave and roughness impulse aspects of the perturbation. A phase-locked velocity decomposition of the entire downstream flow field revealed strongly coherent modes of fluctuating velocity, with distinct mode shapes for the streamwise and wall-normal velocity components. Following the analysis of McKeon & Sharma (J. Fluid Mech., vol. 658, 2010, pp. 336–382), the roughness perturbation was treated as a forcing of the Navier–Stokes equation and a linearized analysis employing a modified Orr–Sommerfeld operator was performed. The experimentally ascertained wavespeed of the input disturbance was used to solve for the most amplified singular mode of the Orr–Sommerfeld resolvent. These calculated modes were then compared with the streamwise and wall-normal velocity fluctuations. The discrepancies between the calculated Orr–Sommerfeld resolvent modes and those experimentally observed by phase-locked averaging of the velocity field were postulated to result from the violation of the parallel flow assumption of Orr–Sommerfeld analysis, as well as certain non-equilibrium effects of the roughness. Additionally, some difficulties previously observed using a quasi-laminar eigenmode analysis were also observed under the resolvent approach; however, the resolvent analysis was shown to provide reasonably accurate predictions of velocity fluctuations for the forced Orr–Sommerfeld problem over a portion of the boundary layer, with potential applications to designing efficient flow control strategies. The combined experimental and analytical effort provides a new opportunity to examine the non-equilibrium and forcing effects in a dynamically perturbed flow.


2018 ◽  
Vol 848 ◽  
pp. 631-647 ◽  
Author(s):  
Eunok Yim ◽  
J.-M. Chomaz ◽  
D. Martinand ◽  
E. Serre

The transition to turbulence in the rotating disk boundary layer is investigated in a closed cylindrical rotor–stator cavity via direct numerical simulation (DNS) and linear stability analysis (LSA). The mean flow in the rotor boundary layer is qualitatively similar to the von Kármán self-similarity solution. The mean velocity profiles, however, slightly depart from theory as the rotor edge is approached. Shear and centrifugal effects lead to a locally more unstable mean flow than the self-similarity solution, which acts as a strong source of perturbations. Fluctuations start rising there, as the Reynolds number is increased, eventually leading to an edge-driven global mode, characterized by spiral arms rotating counter-clockwise with respect to the rotor. At larger Reynolds numbers, fluctuations form a steep front, no longer driven by the edge, and followed downstream by a saturated spiral wave, eventually leading to incipient turbulence. Numerical results show that this front results from the superposition of several elephant front-forming global modes, corresponding to unstable azimuthal wavenumbers $m$, in the range $m\in [32,78]$. The spatial growth along the radial direction of the energy of these fluctuations is quantitatively similar to that observed experimentally. This superposition of elephant modes could thus provide an explanation for the discrepancy observed in the single disk configuration, between the corresponding spatial growth rates values measured by experiments on the one hand, and predicted by LSA and DNS performed in an azimuthal sector, on the other hand.


Sign in / Sign up

Export Citation Format

Share Document