Maximum kinetic energy dissipation and the stability of turbulent Poiseuille flow

2015 ◽  
Vol 767 ◽  
pp. 342-363 ◽  
Author(s):  
J. Bertram

AbstractFollowing Malkus’s (J. Fluid Mech., vol. 1, 1956, pp. 521–539) proposal that turbulent Poiseuille channel flow maximises total viscous dissipation $D$, a variety of variational procedures have been explored involving the maximisation of different flow quantities under different constraints. However, the physical justification for these variational procedures has remained unclear. Here we address more recent claims that mean flow viscous dissipation $D_{m}$ should be maximised on the basis of a statistical stability argument, and that maximising $D_{m}$ yields realistic mean velocity profiles (Malkus, J. Fluid Mech., vol. 489, 2003, pp. 185–198). We clarify the connection between maximising $D_{m}$ and other flow quantities, verify Malkus & Smith’s, (J. Fluid Mech., vol. 208, 1989, pp. 479–507) claim that maximising the ‘efficiency’ yields realistic profiles and show that, in contrast, maximising $D_{m}$ does not yield realistic mean velocity profiles as recently claimed. This leads us to revisit Malkus’s statistical stability argument for maximising $D_{m}$ and to address some of its limitations. We propose an alternative statistical stability argument leading to a principle of minimum kinetic energy for fixed pressure gradient, which suggests a principle of maximum $D$ for fixed Reynolds number under certain conditions. We discuss possible ways to reconcile these conflicting results, focusing on the choice of constraints.

2018 ◽  
Vol 857 ◽  
pp. 345-373 ◽  
Author(s):  
Davide Gatti ◽  
Andrea Cimarelli ◽  
Yosuke Hasegawa ◽  
Bettina Frohnapfel ◽  
Maurizio Quadrio

This paper addresses the integral energy fluxes in natural and controlled turbulent channel flows, where active skin-friction drag reduction techniques allow a more efficient use of the available power. We study whether the increased efficiency shows any general trend in how energy is dissipated by the mean velocity field (mean dissipation) and by the fluctuating velocity field (turbulent dissipation). Direct numerical simulations (DNS) of different control strategies are performed at constant power input (CPI), so that at statistical equilibrium, each flow (either uncontrolled or controlled by different means) has the same power input, hence the same global energy flux and, by definition, the same total energy dissipation rate. The simulations reveal that changes in mean and turbulent energy dissipation rates can be of either sign in a successfully controlled flow. A quantitative description of these changes is made possible by a new decomposition of the total dissipation, stemming from an extended Reynolds decomposition, where the mean velocity is split into a laminar component and a deviation from it. Thanks to the analytical expressions of the laminar quantities, exact relationships are derived that link the achieved flow rate increase and all energy fluxes in the flow system with two wall-normal integrals of the Reynolds shear stress and the Reynolds number. The dependence of the energy fluxes on the Reynolds number is elucidated with a simple model in which the control-dependent changes of the Reynolds shear stress are accounted for via a modification of the mean velocity profile. The physical meaning of the energy fluxes stemming from the new decomposition unveils their inter-relations and connection to flow control, so that a clear target for flow control can be identified.


Author(s):  
G. K. Batchelor

A new and fruitful theory of turbulent motion was published in 1941 by A. N. Kolmogoroff. It does not seem to be as widely known outside the U.S.S.R. as its importance warrants, and the present paper therefore describes the theory in some detail before presenting a number of extensions and making a comparison of experimental results with some of the theoretical predictions.Kolmogoroff's basic notion is that at high Reynolds number all kinds of turbulent motion, of arbitrary mean-flow characteristics, show a similar structure if attention is confined to the smallest eddies. The motion due to these eddies of limited size is conceived to be isotropic and statistically steady. Within this range of eddies we recognize two limiting processes. The influence of viscosity on the larger eddies of the range is negligible if the Reynolds number is large enough, so that their motion is determined entirely by the amount of energy which they are continually passing on to smaller eddies. This quantity of energy is the local mean energy dissipation due to turbulence. On the other hand, the smaller eddies of the range dissipate through the action of viscosity a considerable proportion of the energy which they receive, and the motion of the very smallest eddies is entirely laminar. The analytical expression of this physical picture is that the motion due to eddies less than a certain limiting size in an arbitrary field of turbulence is determined uniquely by two quantities, the viscosity and the local mean energy dissipation per unit mass of the fluid.The mathematical method of describing the motion due to eddies of a particular size is to construct correlations between the differences of parallel-velocity components at two points at an appropriate distance apart. Kinematical results analogous to those for turbulence which is isotropic in the ordinary sense are obtained, and then the scalar functions occurring in the expressions for the correlations are determined by dimensional analysis. The consequences of the theory in the case of turbulence which possesses ordinary isotropy are analysed and various predictions are made. One of these, namely that dimensionless ratios of moments of the probability distribution of the rate of extension of the fluid in any direction are universal constants, is confirmed by recent experiments, so far as the second and third moments are concerned. In several other cases it can be said that relations predicted by the theory have the correct form, but further experiments at Reynolds numbers higher than those hitherto used will be required before the theory can be regarded as fully confirmed. If valid, Kolmogoroff's theory of locally isotropic turbulence will provide a powerful tool for the analysis of problems of non-uniform turbulent flow, and for the determination of statistical characteristics of space and time derivatives of quantities influenced by the turbulence.


2001 ◽  
Vol 439 ◽  
pp. 403-412 ◽  
Author(s):  
RAMA GOVINDARAJAN ◽  
R. NARASIMHA

We first demonstrate that, if the contributions of higher-order mean flow are ignored, the parabolized stability equations (Bertolotti et al. 1992) and the ‘full’ non-parallel equation of Govindarajan & Narasimha (1995, hereafter GN95) are both equivalent to order R−1 in the local Reynolds number R to Gaster's (1974) equation for the stability of spatially developing boundary layers. It is therefore of some concern that a detailed comparison between Gaster (1974) and GN95 reveals a small difference in the computed amplitude ratios. Although this difference is not significant in practical terms in Blasius flow, it is traced here to the approximation, in Gaster's method, of neglecting the change in eigenfunction shape due to flow non-parallelism. This approximation is not justified in the critical and the wall layers, where the neglected term is respectively O(R−2/3) and O(R−1) compared to the largest term. The excellent agreement of GN95 with exact numerical simulations, on the other hand, suggests that the effect of change in eigenfunction is accurately taken into account in that paper.


Water ◽  
2020 ◽  
Vol 12 (3) ◽  
pp. 750
Author(s):  
Ender Demirel ◽  
Mustafa M. Aral

Interaction of recirculating and mean flow downstream of a submerged gate may form significant vortex structures, which may affect the stability of the gate. Although these flow structures that appear in submerged hydraulic jumps received considerable attention in the literature, relatively less work was devoted to the analysis and suppression of the vortex structures downstream of a submerged gate. In this work, internal flow structure and vortex dynamics around a submerged gate were investigated through laboratory tests and large-eddy simulation (LES) using computational fluid dynamics (CFD). It is shown that numerical results obtained for mean velocity field are in good agreement with the experimental measurements. A helical vortex pair connected with a horseshoe vortex system was identified within the roller region using high-resolution numerical simulations. Damping performance of different types of anti-vortex elements placed on the downstream face of the gate are evaluated based on numerical studies. It is shown that the horizontal porous baffle mounted at an elevation below the free surface reduced the vortex magnitudes in the roller region by 26.8%. With the implementation of the proposed vortex breaker, lift forces acting on the gate lip were reduced by 9.4% and drag forces acting on the downstream face of the gate were reduced by 8.6%. Finally, in this study, we assess the performance of the vortex breaker under different flow conditions.


1967 ◽  
Vol 27 (4) ◽  
pp. 657-689 ◽  
Author(s):  
R. E. Kelly

In experiments concerning the instability of free shear layers, oscillations have been observed in the downstream flow which have a frequency exactly half that of the dominant oscillation closer to the origin of the layer. The present analysis indicates that the phenomenon is due to a secondary instability associated with the nearly periodic flow which arises from the finite-amplitude growth of the fundamental disturbance.At first, however, the stability of inviscid shear flows, consisting of a non-zero mean component, together with a component periodic in the direction of flow and with time, is investigated fairly generally. It is found that the periodic component can serve as a means by which waves with twice the wavelength of the periodic component can be reinforced. The dependence of the growth rate of the subharmonic wave upon the amplitude of the periodic component is found for the case when the mean flow profile is of the hyperbolic-tangent type. In order that the subharmonic growth rate may exceed that of the most unstable disturbance associated with the mean flow, the amplitude of the streamwise component of the periodic flow is required to be about 12 % of the mean velocity difference across the shear layer. This represents order-of-magnitude agreement with experiment.Other possibilities of interaction between disturbances and the periodic flow are discussed, and the concluding section contains a discussion of the interactions on the basis of the energy equation.


Author(s):  
Sylvain C. Humbert ◽  
Jonas Moeck ◽  
Alessandro Orchini ◽  
Christian Oliver Paschereit

Abstract Thermoacoustic oscillations in axisymmetric annular combustors are generally coupled by degenerate azimuthal modes, which can be of standing or spinning nature. Symmetry breaking due to the presence of a mean azimuthal flow splits the degenerate thermoacoustic eigenvalues, resulting in pairs of counter-spinning modes with close but distinct frequencies and growth rates. In this study, experiments have been performed using an annular system where the thermoacoustic feedback due to the flames is mimicked by twelve identical electroacoustic feedback loops. The mean azimuthal flow is generated by fans. We investigate the standing/spinning nature of the oscillations as a function of the Mach number for two types of initial states, and how the stability of the system is affected by the mean azimuthal flow. It is found that spinning, standing or mixed modes can be encountered at very low Mach number, but increasing the mean velocity promotes one spinning direction. At sufficiently high Mach number, spinning modes are observed in the limit cycle oscillations. In some cases, the initial conditions have a significant impact on the final state of the system. It is found that the presence of a mean azimuthal flow increases the acoustic damping. This has a beneficial effect on stability: it often reduces the amplitude of the self-sustained oscillations, and can even suppress them in some cases. However, we observe that the suppression of a mode due to the mean flow may destabilize another one. We discuss our findings in relation with an existing low-order model.


2009 ◽  
Vol 23 (03) ◽  
pp. 509-512 ◽  
Author(s):  
SUHUA SHEN ◽  
JIANZHONG LIN

To explore the rheological property in turbulent channel flow of fiber suspensions, the equation of probability distribution function for mean fiber orientation and the Reynolds averaged Navier-Stokes equation with the term of additional stress resulted from fibers were solved with numerical methods to get the distributions of the mean velocity and turbulent kinetic energy. The simulation results show that the effect of fibers on turbulent channel flow is equivalent to an additional viscosity. The turbulent velocity profiles of fiber suspension become gradually sharper by increasing the fiber concentration and/or decreasing the Reynolds number. The turbulent kinetic energy will increase with increasing Reynolds number and fiber concentration.


1972 ◽  
Vol 52 (3) ◽  
pp. 401-423 ◽  
Author(s):  
Timothy W. Kao ◽  
Cheol Park

The stability of the laminar co-current flow of two fluids, oil and water, in a rectangular channel was investigated experimentally, with and without artificial excitation. For the ratio of viscosity explored, only the disturbances in water grew in the beginning stages of transition to turbulence. The critical water Reynolds number, based upon the hydraulic diameter of the channel and the superficial velocity defined by the ratio of flow rate of water to total cross-sectional area of the channel, was found to be 2300. The behaviour of damped and growing shear waves in water was examined in detail using artificial excitation and briefly compared with that observed in Part 1. Mean flow profiles, the amplitude distribution of disturbances in water, the amplification rate, wave speed and wavenumbers were obtained. A neutral stability boundary in the wave-number, water Reynolds number plane was also obtained experimentally.It was found that in natural transition the interfacial mode was not excited. The first appearance of interfacial waves was actually a manifestation of the shear waves in water. The role of the interface in the transition range from laminar to turbulent flow in water was to introduce and enhance spanwise oscillation in the water phase and to hasten the process of breakdown for growing disturbances.


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’.


1997 ◽  
Vol 334 ◽  
pp. 61-86 ◽  
Author(s):  
PAUL PICCIRILLO ◽  
CHARLES W. VAN ATTA

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).


Sign in / Sign up

Export Citation Format

Share Document