scholarly journals Transition to turbulence in the rotating disk boundary layer of a rotor–stator cavity

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.

2013 ◽  
Vol 718 ◽  
pp. 596-621 ◽  
Author(s):  
J. C. Klewicki

AbstractThis study investigates how and why dynamical self-similarities emerge with increasing Reynolds number within the canonical wall flows beyond the transitional regime. An overarching aim is to advance a mechanistically coherent description of turbulent wall-flow dynamics that is mathematically tractable and grounded in the mean dynamical equations. As revealed by the analysis of Fife, Klewicki & Wei (J. Discrete Continuous Dyn. Syst.A, vol. 24, 2009, pp. 781–807), the equations that respectively describe the mean dynamics of turbulent channel, pipe and boundary layer flows formally admit invariant forms. These expose an underlying self-similar structure. In all cases, two kinds of dynamical self-similarity are shown to exist on an internal domain that, for all Reynolds numbers, extends from$O(\nu / {u}_{\tau } )$to$O(\delta )$, where$\nu $is the kinematic viscosity,${u}_{\tau } $is the friction velocity and$\delta $is the half-channel height, pipe radius, or boundary layer thickness. The simpler of the two self-similarities is operative on a large outer portion of the relevant domain. This self-similarity leads to an explicit analytical closure of the mean momentum equation. This self-similarity also underlies the emergence of a logarithmic mean velocity profile. A more complicated kind a self-similarity emerges asymptotically over a smaller domain closer to the wall. The simpler self-similarity allows the mean dynamical equation to be written as a closed system of nonlinear ordinary differential equations that, like the similarity solution for the laminar flat-plate boundary layer, can be numerically integrated. The resulting similarity solutions are demonstrated to exhibit nearly exact agreement with direct numerical simulations over the solution domain specified by the theory. At the Reynolds numbers investigated, the outer similarity solution is shown to be operative over a domain that encompasses${\sim }40\hspace{0.167em} \% $of the overall width of the flow. Other properties predicted by the theory are also shown to be well supported by existing data.


2019 ◽  
Vol XVI (2) ◽  
pp. 13-22
Author(s):  
Muhammad Ehtisham Siddiqui

Three-dimensional boundary-layer flow is well known for its abrupt and sharp transition from laminar to turbulent regime. The presented study is a first attempt to achieve the target of delaying the natural transition to turbulence. The behaviour of two different shaped and sized stationary disturbances (in the laboratory frame) on the rotating-disk boundary layer flow is investigated. These disturbances are placed at dimensionless radial location (Rf = 340) which lies within the convectively unstable zone over a rotating-disk. Mean velocity profiles were measured using constant-temperature hot-wire anemometry. By careful analysis of experimental data, the instability of these disturbance wakes and its estimated orientation within the boundary-layer were investigated.


2021 ◽  
Vol 929 ◽  
Author(s):  
N. Agastya Balantrapu ◽  
Christopher Hickling ◽  
W. Nathan Alexander ◽  
William Devenport

Experiments were performed over a body of revolution at a length-based Reynolds number of 1.9 million. While the lateral curvature parameters are moderate ( $\delta /r_s < 2, r_s^+>500$ , where $\delta$ is the boundary layer thickness and r s is the radius of curvature), the pressure gradient is increasingly adverse ( $\beta _{C} \in [5 \text {--} 18]$ where $\beta_{C}$ is Clauser’s pressure gradient parameter), representative of vehicle-relevant conditions. The mean flow in the outer regions of this fully attached boundary layer displays some properties of a free-shear layer, with the mean-velocity and turbulence intensity profiles attaining self-similarity with the ‘embedded shear layer’ scaling (Schatzman & Thomas, J. Fluid Mech., vol. 815, 2017, pp. 592–642). Spectral analysis of the streamwise turbulence revealed that, as the mean flow decelerates, the large-scale motions energize across the boundary layer, growing proportionally with the boundary layer thickness. When scaled with the shear layer parameters, the distribution of the energy in the low-frequency region is approximately self-similar, emphasizing the role of the embedded shear layer in the large-scale motions. The correlation structure of the boundary layer is discussed at length to supply information towards the development of turbulence and aeroacoustic models. One major finding is that the estimation of integral turbulence length scales from single-point measurements, via Taylor's hypothesis, requires significant corrections to the convection velocity in the inner 50 % of the boundary layer. The apparent convection velocity (estimated from the ratio of integral length scale to the time scale), is approximately 40 % greater than the local mean velocity, suggesting the turbulence is convected much faster than previously thought. Closer to the wall even higher corrections are required.


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.


A set of model equations is given to describe the gross features of a statistically steady or 'slowly varying’ inhomogeneous field of turbulence and the mean velocity distribution. The equations are based on the idea that turbulence can be characterized by ‘densities’ which obey nonlinear diffusion equations. The diffusion equations contain terms to describe the convection by the mean flow, the amplification due to interaction with a mean velocity gradient, the dissipation due to the interaction of the turbulence with itself, and the dif­fusion also due to the self interaction. The equations are similar to a set proposed by Kolmo­gorov (1942). It is assumed that both an ‘energy density’ and a ‘vorticity density’ satisfy diffusion equations, and that the self diffusion is described by an eddy viscosity which is a function of the energy and vorticity densities; the eddy viscosity is also assumed to describe the diffu­sion of mean momentum by the turbulent fluctuations. It is shown that with simple and plausible assumptions about the nature of the interaction terms, the equations form a closed set. The appropriate boundary conditions at a solid wall and a turbulent interface, with and without entrainment, are discussed. It is shown that the dimensionless constants which appear in the equations can all be estimated by general arguments. The equations are then found to predict the von Kármán constant in the law of the wall with reasonable accuracy. An analytical solution is given for Couette flow, and the result of a numerical study of plane Poiseuille flow is described. The equations are also applied to free turbulent flows. It is shown that the model equations completely determine the structure of the similarity solutions, with the rate of spread, for instance, determined by the solution of a nonlinear eigenvalue problem. Numerical solutions have been obtained for the two-dimensional wake and jet. The agreement with experiment is good. The solutions have a sharp interface between turbulent and non-turbulent regions and the mean velocity in the turbulent part varies linearly with distance from the interface. The equations are applied qualitatively to the accelerating boundary layer in flow towards a line sink, and the decelerating boundary layer with zero skin friction. In the latter case, the equations predict that the mean velocity should vary near the wall like the 5/3 power of the distance. It is shown that viscosity can be incorporated formally into the model equations and that a structure can be given to the interface between turbulent and non-turbulent parts of the flow.


An experimental investigation has been made of turbulent boundary layer response to harmonic oscillations associated with a travelling wave imposed on an otherwise constant freestream velocity and convected in the freestream direction. The tests covered oscillation frequencies of 4-12 Hz for freestream amplitudes of up to 11% of the mean velocity. Additional steady flow measurements were used to infer the quasi-steady response to freestream oscillations. The results show a welcome insensitivity of the mean flow and turbulent intensity distributions to the freestream oscillations tested. An approximate analysis based on these results has been developed. It is probably of limited validity but it does provide a useful guide to the physical processes involved. The effects on boundary layer response of varying the travelling wave convection velocity and frequency of oscillation are illustrated by the analysis and show a behaviour broadly similar to that of laminar boundary layers. The travelling wave convection velocity exhibits a dominant influence on the turbulent boundary layer response to freestream oscillations.


1982 ◽  
Vol 119 ◽  
pp. 121-153 ◽  
Author(s):  
Udo R. Müller

An experimental study of a steady, incompressible, three-dimensional turbulent boundary layer approaching separation is reported. The flow field external to the boundary layer was deflected laterally by turning vanes so that streamwise flow deceleration occurred simultaneous with cross-flow acceleration. At 21 stations profiles of the mean-velocity components and of the six Reynolds stresses were measured with single- and X-hot-wire probes, which were rotatable around their longitudinal axes. The calibration of the hot wires with respect to magnitude and direction of the velocity vector as well as the method of evaluating the Reynolds stresses from the measured data are described in a separate paper (Müller 1982, hereinafter referred to as II). At each measuring station the wall shear stress was inferred from a Preston-tube measurement as well as from a Clauser chart. With the measured profiles of the mean velocities and of the Reynolds stresses several assumptions used for turbulence modelling were checked for their validity in this flow. For example, eddy viscosities for both tangential directions and the corresponding mixing lengths as well as the ratio of resultant turbulent shear stress to turbulent kinetic energy were derived from the data.


2000 ◽  
Vol 407 ◽  
pp. 145-166 ◽  
Author(s):  
LIAN SHEN ◽  
GEORGE S. TRIANTAFYLLOU ◽  
DICK K. P. YUE

We study numerically and analytically the turbulent diffusion characteristics in a low-Froude-number turbulent shear flow beneath a free surface. In the numerical study, the Navier–Stokes equations are solved directly subject to viscous boundary conditions at the free surface. From an ensemble of such simulations, we find that a boundary layer develops at the free surface characterized by a fast reduction in the value of the eddy viscosity. As the free surface is approached, the magnitude of the mean shear initially increases over the boundary (outer) layer, reaches a maximum and then drops to zero inside a much thinner inner layer. To understand and model this behaviour, we derive an analytical similarity solution for the mean flow. This solution predicts well the shape and the time-scaling behaviour of the mean flow obtained in the direct simulations. The theoretical solution is then used to derive scaling relations for the thickness of the inner and outer layers. Based on this similarity solution, we propose a free-surface function model for large-eddy simulations of free-surface turbulence. This new model correctly accounts for the variations of the Smagorinsky coefficient over the free-surface boundary layer and is validated in both a priori and a posteriori tests.


2016 ◽  
Vol 789 ◽  
pp. 127-165 ◽  
Author(s):  
Xiang I. A. Yang ◽  
Jasim Sadique ◽  
Rajat Mittal ◽  
Charles Meneveau

We conduct a series of large-eddy simulations (LES) to examine the mean flow behaviour within the roughness layer of turbulent boundary layer flow over rough surfaces. We consider several configurations consisting of arrays of rectangular-prism roughness elements with various spacings, aspect ratios and height distributions. The results provide clear evidence for exponential behaviour of the mean flow with respect to the wall normal distance. Such behaviour has been proposed before (see, e.g., Cionco, 1966 Tech. Rep. DTIC document), and is represented as $U(z)/U_{h}=\exp [a(z/h-1)]$, where $U(z)$ is the spatially/temporally averaged fluid velocity, $z$ is the wall normal distance, $h$ represents the height of the roughness elements and $U_{h}$ is the velocity at $z=h$. The attenuation factor $a$ depends on the density of the roughness element distribution and details of the roughness distribution on the wall. Once established, the generic velocity profile shape is used to formulate a fully analytical model for the effective drag exerted by turbulent flow on a surface covered with arrays of rectangular-prism roughness elements. The approach is based on the von Karman–Pohlhausen integral method, in which a shape function is assumed for the mean velocity profile and its parameters are determined based on momentum conservation and other fundamental constraints. In order to determine the attenuation parameter $a$, wake interactions among surface roughness elements are accounted for by using the concept of flow sheltering. The model transitions smoothly between ‘$k$’ and ‘$d$’ type roughness conditions depending on the surface coverage density and the detailed geometry of roughness elements. Comparisons between model predictions and experimental/numerical data from the existing literature as well as LES data from this study are presented. It is shown that the analytical model provides good predictions of mean velocity and drag forces for the cases considered, thus raising the hope that analytical roughness modelling based on surface geometry is possible, at least for cases when the location of flow separation over surface elements can be easily predicted, as in the case of wall-attached rectangular-prism roughness elements.


2008 ◽  
Vol 130 (6) ◽  
Author(s):  
Martin Agelinchaab ◽  
Mark F. Tachie

A particle image velocimetry is used to study the mean and turbulent fields of separated and redeveloping flow over square, rectangular, and semicircular blocks fixed to the bottom wall of an open channel. The open channel flow is characterized by high background turbulence level, and the ratio of the upstream boundary layer thickness to block height is considerably higher than in prior experiments. The variation of the Reynolds stresses along the dividing streamlines is discussed within the context of vortex stretching, longitudinal strain rate, and wall damping. It appears that wall damping is a more dominant mechanism in the vicinity of reattachment. In the recirculation and reattachment regions, profiles of the mean velocity, turbulent quantities, and transport terms are used to document the salient features of block geometry on the flow. The flow characteristics in these regions strongly depend on block geometry. Downstream of reattachment, a new shear layer is formed, and the redevelopment of the shear layer toward the upstream open channel boundary layer is studied using the boundary layer parameters and Reynolds stresses. The results show that the mean flow rapidly redeveloped so that the Clauser parameter recovered to its upstream value at 90 step heights downstream of reattachment. However, the rate of development close to reattachment strongly depends on block geometry.


Sign in / Sign up

Export Citation Format

Share Document