Spanwise reflection symmetry breaking and turbulence control: plane Couette flow

2014 ◽  
Vol 745 ◽  
pp. 300-320
Author(s):  
G. Chagelishvili ◽  
G. Khujadze ◽  
H. Foysi ◽  
M. Oberlack
Keyword(s):  
Kinetic Energy ◽  
Shear Flow ◽  
Turbulent Kinetic Energy ◽  
Couette Flow ◽  
Promising Result ◽  
Spanwise Direction ◽  
Reflection Symmetry ◽  
Plane Couette Flow ◽  
Mean Flow ◽  
Turbulence Control

AbstractWe propose and analyse a new strategy of shear flow turbulence control that can be realized by the following steps: (i) imposing specially designed seed velocity perturbations, which are non-symmetric in the spanwise direction, at the walls of a flow; (ii) the configuration of the latter ensures a gain of shear flow energy and the breaking of turbulence spanwise reflection symmetry: this leads to the generation of spanwise mean flow; (iii) that changes the self-sustained dynamics of turbulence and results in a considerable reduction of the turbulence level and the production of turbulent kinetic energy. In fact, by this strategy the shear flow transient growth mechanism is activated and the formed spanwise mean flow is an intrinsic, nonlinear composition of the controlled turbulence and not directly introduced in the system. In the present paper, a weak near-wall volume forcing is designed to impose the velocity perturbations with required characteristics in the flow. The efficiency of the proposed scheme has been demonstrated by direct numerical simulation using plane Couette flow as a representative example. A promising result was obtained: after a careful parameter selection, the forcing reduces the turbulence kinetic energy and its production by up to one-third. The strategy can be naturally applied to other wall-bounded flows, e.g. channel and boundary-layer flows. Of course, the considered volume force is theoretical and hypothetical. Nevertheless, it helps to gain knowledge concerning the design of the seed velocity field that is necessary to be imposed in the flow to achieve a significant reduction of the turbulent kinetic energy. This is convincing with regard to a new control strategy, which could be based on specially constructed blowing/suction or riblets, by employing the insight gained by the comprehension of the results obtained using the investigated methodology in this paper.

Secondary flow in weakly rotating turbulent plane Couette flow

1996 ◽  
Vol 317 ◽  
pp. 195-214 ◽  
Author(s):  
Knut H. Bech ◽  
Helge I. Andersson
Keyword(s):  
Kinetic Energy ◽  
Secondary Flow ◽  
Couette Flow ◽  
Three Dimensional ◽  
Wall Turbulence ◽  
Plane Couette Flow ◽  
Mean Flow ◽  
Comparable Amount ◽  
Flow Vorticity ◽  
Turbulent Plane

As in the laminar case, the turbulent plane Couette flow is unstable (stable) with respect to roll cell instabilities when the weak background angular velocity Ωk is antiparallel (parallel) to the spanwise mean flow vorticity (-dU/dy)k. The critical value of the rotation number Ro, based on 2Ω and dU/dy of the corresponding laminar flow, was estimated as 0.0002 at a low Reynolds number with fully developed turbulence. Direct numerical simulations were performed for Ro = ±0.01 and compared with earlier results for non-rotating Couette flow. At the low rotation rates considered, both senses of rotation damped the turbulence and the number of near-wall turbulence-generating events was reduced. The destabilized flow was more energetic, but less three-dimensional, than the non-rotating flow. In the destabilized case, the two-dimensional roll cells extracted a comparable amount of kinetic energy from the mean flow as did the turbulence, thereby decreasing the turbulent kinetic energy. The turbulence anisotropy was practically unaffected by weak spanwise rotation, while the secondary flow was highly anisotropic due to its inability to contract and expand in the streamwise direction.

Experimental investigation of coherent structures of a three-dimensional separated turbulent boundary layer

2018 ◽  
Vol 859 ◽  
pp. 1-32 ◽  
Author(s):  
Mohammad Elyasi ◽  
Sina Ghaemi
Keyword(s):  
Boundary Layer ◽  
Kinetic Energy ◽  
Turbulent Boundary Layer ◽  
Saddle Point ◽  
Turbulent Kinetic Energy ◽  
Coherent Structures ◽  
Three Dimensional ◽  
Spanwise Direction ◽  
Mean Flow ◽  
Elliptical Cross

Coherent structures of a three-dimensional (3D) separation due to an adverse pressure gradient are investigated experimentally. The flow set-up consists of a flat plate to develop a turbulent boundary layer upstream of an asymmetric two-dimensional diffuser with one diverging surface. The diffuser surface has an initial mild curvature followed by a flat section where flow separation occurs. The top and the two sidewalls of the diffuser are not equipped with any flow control mechanism to form a 3D separation. Planar particle image velocimetry (PIV) using four side-by-side cameras is applied to characterize the flow with high spatial resolution over a large streamwise-wall-normal field of view (FOV). Tomographic PIV (tomo-PIV) is also applied for volumetric measurement in a domain flush with the flat surface of the diffuser. The mean flow obtained from averaging instantaneous velocity fields of this intermittent unsteady flow appears as a vortex with an elliptical cross-section. The major axis of the ellipse is tilted with respect to the streamwise direction. As a result, the average velocity in the mid-span of the diffuser has an upstream forward flow and a downstream backward flow, separated by a point of zero wall shear stress. Sweep motions mainly carry out transport of turbulent kinetic energy upstream of this point, while ejections dominate at the downstream region. In the instantaneous flow fields, forward and backward flows have equivalent strength, and the separation front is extended in the spanwise direction. The conditional average of the separation instants forms a saddle-point structure with streamlines converging in the spanwise direction. Proper orthogonal decomposition (POD) of the tomo-PIV data demonstrates that about 42 % of the turbulent kinetic energy is present in the first pair of modes, with a strong spanwise component. The spatial modes of POD also show focus, node and saddle-point structures. The average of the coefficients of the dominant POD modes during the separation events is used to develop a reduced-order model (ROM). Based on the ROM, the instantaneous 3D separation over the diffuser is a saddle-point structure interacting with focus-type structures.

A Comparison of Cavitation Inception Index Measurements to the Spatial Pressure Distribution Within a 2D Cavity Shear Flow

2007 ◽  
Author(s):  
Xiaofeng Liu ◽  
Joseph Katz
Keyword(s):  
Kinetic Energy ◽  
Shear Flow ◽  
Pressure Distribution ◽  
Turbulent Kinetic Energy ◽  
Mean Flow ◽  
Local Pressure ◽  
Trailing Edge ◽  
Integration Algorithm ◽  
Cavitation Inception ◽  
Pressure Diffusion

Direct cavitation inception index measurements and observation on occurrence of cavitation are compared to results of novel spatial pressure distribution measurements in a 2D cavity shear flow. This non-intrusive technique utilizes four-exposure PIV to measure the distribution of material acceleration, and integrating it by means of omni-directional virtual boundary integration algorithm to obtain the pressure distribution (Liu and Katz, 2006). Consequently, it provides the instantaneous spatial distributions of velocity, material acceleration and pressure over a sample area along with their statistics. The present Reynolds numbers based on the cavity length vary from 1.7×105 to 3.4×105. High-speed imaging of cavitation inception, recorded at 30,000 fps, indicates that for this 2D cavity flow, the onset of cavitation always occurs on the top of the cavity trailing edge, regardless of the free stream speed. With decreasing pressure cavitation intermittently expand to the region located just in front of the cavity. The time-averaged spatial pressure distribution has a minimum just above the trailing edge due to the interaction of the impinging shear layer with the trailing wall. Around the cavity trailing edge, the mean flow first decelerates due to the impingement, but then accelerates right above the trailing edge, creating a local pressure minimum there. RMS values and PDFs of pressure fluctuations show that the highest fluctuations occur around the cavity trailing edge, and that the pressure peaks are consistent with the measured cavitation inception indices. There is also agreement between pressure statistics and conditions of appearance of cavitation in front of the trailing edge. The paper also provides the first directly measured experimental data on pressure-velocity correlation and pressure diffusion terms that appear in the evolution equation for turbulent kinetic energy. Results compared to other terms that act as sources and sinks in the turbulent kinetic energy balance. It is evident that near the trailing edge of the cavity, the contribution of pressure diffusion is comparable to that of turbulent kinetic energy production rate, and is much larger than the turbulent diffusion rate. Trends and spatial distribution of pressure diffusion also differs from those of turbulence diffusion.

The Application of Eddy-Viscosity Stress Limiters for Modeling Cross-Flow Separation

2010 ◽  
Vol 132 (9) ◽  
Author(s):  
P. A. Gregory ◽  
P. N. Joubert ◽  
M. S. Chong ◽  
A. Ooi
Keyword(s):  
Kinetic Energy ◽  
Skin Friction ◽  
Turbulent Kinetic Energy ◽  
Flow Separation ◽  
Eddy Viscosity ◽  
Cross Flow ◽  
The Body ◽  
Mean Flow ◽  
Viscosity Models ◽  
Experimental Apparatus

The ability of eddy-viscosity models to simulate the turbulent wake produced by cross-flow separation over a curved body of revolution is assessed. The results obtained using the standard k−ω model show excessive levels of turbulent kinetic energy k in the vicinity of the stagnation point at the nose of the body. Additionally, high levels of k are observed throughout the wake. Enforcing laminar flow upstream of the nose (which replicates the experimental apparatus more accurately) gives more accurate estimates of k throughout the flowfield. A stress limiter in the form of Durbin’s T-limit modification for eddy-viscosity models is implemented for the k−ω model, and its effect on the computed surface pressures, skin friction, and surface flow features is assessed. Additionally, the effect of the T-limit modification on both the mean flow and the turbulent flow quantities within the wake is also examined. The use of the T-limit modification gives significant improvements in predicted levels of turbulent kinetic energy and Reynolds stresses within the wake. However, predicted values of skin friction in regions of attached flow become up to 50% greater than the experimental values when the T-limit is used. This is due to higher values of near-wall turbulence being created with the T-limit.

On the structure of shear stress and turbulent kinetic energy flux across the roughness layer of a gravel-bed channel flow

2009 ◽  
Vol 638 ◽  
pp. 423-452 ◽  
Author(s):  
EMMANUEL MIGNOT ◽  
D. HURTHER ◽  
E. BARTHELEMY
Keyword(s):  
Shear Stress ◽  
Kinetic Energy ◽  
Turbulent Kinetic Energy ◽  
Channel Flow ◽  
Mean Flow ◽  
Local Flow ◽  
Gravel Bed ◽  
Wall Flows ◽  
The Mean ◽  
Conditional Statistics

This study examines the structure of shear stress and turbulent kinetic energy (TKE) flux across the roughness layer of a uniform, fully rough gravel-bed channel flow (ks+ ≫ 100, δ/k = 20) using high-resolution acoustic Doppler velocity profiler measurements. The studied gravel-bed roughness layer exhibits a complex random multi-scale roughness structure in strong contrast with conceptualized k- or d-type roughness in standard rough-wall flows. Within the roughness layer, strong spatial variability of all time-averaged flow quantities are observed affecting up to 40% of the boundary layer height. This variability is attributed to the presence of bed zones with emanating bed protuberances (or gravel clusters) acting as local flow obstacles and bed zones of more homogenous roughness of densely packed gravel elements. Considering the strong spatial mean flow variability across the roughness layer, a spatio-temporal averaging procedure, called double averaging (DA), has been applied to the analysed flow quantities. Three aspects have been addressed: (a) the DA shear stress and DA TKE flux in specific bed zones associated with three classes of velocity profiles as previously proposed in Mignot, Barthélemy & Hurther (J. Fluid Mech., vol. 618, 2009, p. 279), (b) the global and per class DA conditional statistics of shear stress and associated TKE flux and (c) the contribution of large-scale coherent shear stress structures (LC3S) to the TKE flux across the roughness layer. The mean Reynolds and dispersive shear structure show good agreement between the protuberance bed zones associated with the S-shape/accelerated classes and recent results obtained in standard k-type rough-wall flows (Djenidi et al., Exp. Fluids, vol. 44, 2008, p. 37; Pokrajac, McEwan & Nikora, Exp. Fluids, vol. 45, 2008, p. 73). These gravel-bed protuberances act as local flow obstacles inducing a strong turbulent activity in their wake regions. The conditional statistics show that the Reynolds stress contribution is fairly well distributed between sweep and ejection events, with threshold values ranging from H = 0 to H = 8. However, the TKE flux across the roughness layer primarily results from the residual shear stress between ejection and sweep of very high magnitude (H = 10–20) and of small turbulent scale. Although LC3S are seen to penetrated the interfacial roughness layer, their TKE flux contribution is found to be negligible compared to the very energetic small-scale sweep events. These sweeps are dominantly produced in the bed zones of local gravel protuberances where the velocity profiles are inflexional of S-shape type and the mean flow properties are of mixing-layer flow type as previously shown in Mignot et al. (2009).

On the interactions between large-scale structure and fine-grained turbulence in a free shear flow I. The development of temporal interactions in the mean

1976 ◽  
Vol 352 (1669) ◽  
pp. 213-247 ◽  
Keyword(s):  
Energy Transfer ◽  
Kinetic Energy ◽  
Shear Layer ◽  
Turbulent Kinetic Energy ◽  
Viscous Dissipation ◽  
Mean Flow ◽  
Fine Grained ◽  
Large Eddy ◽  
The Mean ◽  
Three Components

In this problem a mean turbulent shear layer originally exists, homogeneous in the streamwise direction, formed perhaps by previous instabilities, but in equilibrium with the fine-grained turbulence. At a given time, a large eddy of a fixed horizontal wavenumber is initiated. We study the subsequent time development of the non-equilibrium interactions between the three components of flow as they adjust towards ultimate simultaneous equilibrium, using the integrated energy-balance conservation equations to derive the amplitude equations. This necessarily involves the usual averaging procedure and a conditional or phase-averaging procedure by which the large structure motion is educed from the total fluctuations. In general, the mean flow growth is due to the energy transfer to both fluctuating components, the large eddy gains energy from the mean motion and exchanges energy with the fine-grained turbulence, while the fine-grained turbulence gains energy from the mean flow and exchanges with the large eddy and converts its energy to heat through viscous dissipation of the smallest scales. The closure problem is obtained via the shape assumptions which enter into the interaction integrals. The situation in which the fine-grained turbulent kinetic energy production and viscous dissipation are in local balance is considered, the displacement from equilibrium being due only to the energy transfer from the large eddy. The large eddy shape is taken to be two-dimensional, instability-wavelike, with its vorticity axis perpendicular to the direction of the mean outer stream. Prior to averaging, detailed but approximate calculations of the wave-induced turbulent Reynolds stresses are obtained; the product of these stresses with the appropriate large-eddy rates of strain give the energy transfer mechanism between the two disparate scales of fluctuations. Coupled, nonlinear amplitude or energy density equations for the three components of motion are obtained, the coefficients of which are the interaction integrals guided by the shape assumptions. It is found that for the special case of parallel flow, the energy of the large eddy first undergoes a hydrodynamic-instability type of amplification but eventually decays due to the energy transfer to the fine-grained turbulence, while the turbulent kinetic energy is displaced from an original level of equilibrium to a new one because of the ability of the large eddy to negotiate an indirect energy transfer from the mean flow. For the growing shear layer, approximate considerations show that if the mechanism of energy transfer from the large to the small scale is eventually weakened by the shear layer growth compared to the large-eddy production mechanism so that the amplification and decay process repeats, ‘bursts’ of the remnant of the same large eddy will occur repeatedly until an ultimate equilibrium is reached among the three interacting components of motion. However, for the large eddy whose wavenumber corresponds to that of the initially most amplified case, the ‘bursting’ phenomenon is much less pronounced and equilibrium is very nearly reached at the end of the very first ‘burst’.

The evolution of a uniformly sheared thermally stratified turbulent flow

1997 ◽  
Vol 334 ◽  
pp. 61-86 ◽  
Author(s):  
PAUL PICCIRILLO ◽  
CHARLES W. VAN ATTA
Keyword(s):  
Kinetic Energy ◽  
Turbulent Kinetic Energy ◽  
Richardson Number ◽  
Initial Conditions ◽  
Velocity Shear ◽  
Spectral Energy ◽  
Small Scale ◽  
Mean Flow ◽  
Mean Velocity ◽  
Small Scales

Experiments were carried out in a new type of stratified flow facility to study the evolution of turbulence in a mean flow possessing both uniform stable stratification and uniform mean shear.The new facility is a thermally stratified wind tunnel consisting of ten independent supply layers, each with its own blower and heaters, and is capable of producing arbitrary temperature and velocity profiles in the test section. In the experiments, four different sized turbulence-generating grids were used to study the effect of different initial conditions. All three components of the velocity were measured, along with the temperature. Root-mean-square quantities and correlations were measured, along with their corresponding power and cross-spectra.As the gradient Richardson number Ri = N2/(dU/dz)2 was increased, the downstream spatial evolution of the turbulent kinetic energy changed from increasing, to stationary, to decreasing. The stationary value of the Richardson number, Ricr, was found to be an increasing function of the dimensionless shear parameter Sq2/∈ (where S = dU/dz is the mean velocity shear, q2 is the turbulent kinetic energy, and ∈ is the viscous dissipation).The turbulence was found to be highly anisotropic, both at the small scales and at the large scales, and anisotropy was found to increase with increasing Ri. The evolution of the velocity power spectra for Ri [les ] Ricr, in which the energy of the large scales increases while the energy in the small scales decreases, suggests that the small-scale anisotropy is caused, or at least amplified, by buoyancy forces which reduce the amount of spectral energy transfer from large to small scales. For the largest values of Ri, countergradient buoyancy flux occurred for the small scales of the turbulence, an effect noted earlier in the numerical results of Holt et al. (1992), Gerz et al. (1989), and Gerz & Schumann (1991).

scholarly journals The streamwise turbulence intensity in the intermediate layer of turbulent pipe flow

2015 ◽  
Vol 774 ◽  
pp. 324-341 ◽  
Author(s):  
J. C. Vassilicos ◽  
J.-P. Laval ◽  
J.-M. Foucaut ◽  
M. Stanislas
Keyword(s):  
Reynolds Number ◽  
Kinetic Energy ◽  
Turbulent Kinetic Energy ◽  
Intermediate Layer ◽  
High Reynolds Number ◽  
Logarithmic Derivative ◽  
Turbulent Pipe Flow ◽  
Mean Flow ◽  
The Mean ◽  
High Reynolds

The spectral model of Perryet al. (J. Fluid Mech., vol. 165, 1986, pp. 163–199) predicts that the integral length scale varies very slowly with distance to the wall in the intermediate layer. The only way for the integral length scale’s variation to be more realistic while keeping with the Townsend–Perry attached eddy spectrum is to add a new wavenumber range to the model at wavenumbers smaller than that spectrum. This necessary addition can also account for the high-Reynolds-number outer peak of the turbulent kinetic energy in the intermediate layer. An analytic expression is obtained for this outer peak in agreement with extremely high-Reynolds-number data by Hultmarket al. (Phys. Rev. Lett., vol. 108, 2012, 094501;J. Fluid Mech., vol. 728, 2013, pp. 376–395). Townsend’s (The Structure of Turbulent Shear Flows, 1976, Cambridge University Press) production–dissipation balance and the finding of Dallaset al. (Phys. Rev. E, vol. 80, 2009, 046306) that, in the intermediate layer, the eddy turnover time scales with skin friction velocity and distance to the wall implies that the logarithmic derivative of the mean flow has an outer peak at the same location as the turbulent kinetic energy. This is seen in the data of Hultmarket al. (Phys. Rev. Lett., vol. 108, 2012, 094501;J. Fluid Mech., vol. 728, 2013, pp. 376–395). The same approach also predicts that the logarithmic derivative of the mean flow has a logarithmic decay at distances to the wall larger than the position of the outer peak. This qualitative prediction is also supported by the aforementioned data.

Sign in / Sign up

Export Citation Format

Share Document