Various aspects of fluid vortices

A general review of the current research in vortex dynamics is presented, based on contributions given during a workshop held in May 2003 at Porquerolles, France. This article aims at providing a picture of the work performed on this subject in the French community. Various cases are covered, from 2D vortex patches to 3D vortex tubes; from isolated vortices to shear flows. Different contexts are considered: pure Euler and Navier–Stokes flows as well as stratified, rotating and magnetic flows. To cite this article: I. Delbende et al., C. R. Mecanique 332 (2004).     

Dynamical Equations and Turbulent Closures in Geophysics

Based on scaling arguments the governing equations for turbulent flows are classified. The similarity for stratified and rotating flows is characterized and the conditions for a hydrostatic assumption are shown for several flow regimes. For stratified rotating flow a scale analysis of the turbulent stresses exhibits different classes of second order closure. The complete sets of the governing equations for second and third order turbulent closures are presented. The evolution of the equations is embedded into a historical chronology. Received September 15, 1997     

A Direct Numerical Study of Transport and Anisotropy in a Stably Stratified Turbulent Flow with Uniform Horizontal Shear

Direct numerical simulations of homogeneous turbulence in stably stratified shear flow have been performed to aid the understanding of turbulence and turbulent mixing in geophysical flow. Two cases are compared. In the first case, which has been studied in the past, the mean velocity has vertical shear and the mean density is vertically stably stratified. In the second case, which has not been studied systematically before, the mean velocity has horizontal shear and the mean density is again vertically stably stratified. The critical value of the gradient Richardson number, for which a constant turbulence level is obtained, is found to be an order of magnitude larger in the horizontal shear case. The turbulent transport coefficients of momentum and vertical mass transfer are also an order of magnitude larger in the horizontal shear case. The anisotropy of the turbulence intensities are found to be in the range expected of flows with mean shear with no major qualitative change in the range of Richardson numbers studied here. However, the anisotropy of the turbulent dissipation rate is strongly affected by stratification with the vertical component dominating the others. This revised version was published online in July 2006 with corrections to the Cover Date.     

稳定分层流动中湍流特性的研究综述

本文总结了近60 年来分层流动中湍流特性研究的成果. 主要从两个方面进行了综述:(1) 分层流动中湍流场的演变和混合. 在这方面主要分析稳定分层对湍流混合和湍流结构的影响, 以及混合层内湍流结构的特性和混合层的演化规律. (2) 分层流动中湍流的扩散和输运. 动量和标量的逆梯度输运特性是分层湍流研究的一个重要方向;分析分层对湍流扩散的影响. 并指出了一些值得今后进一步研究的方向.     

A Two-Time-Scale Model for Turbulent Mixing Flows Induced by Rayleigh–Taylor and Richtmyer–Meshkov Instabilities

A two-time-scale closure model for compressible flows previously developed is extended to turbulent Rayleigh-Taylor and Richtmyer-Meshkov driven flows where mixing coexists with mean pressure gradients. Two model coefficients are calibrated with the help of Canuto-Goldman's model. For several Rayleigh-Taylor configurations, it is shown that the characteristic lengths scale as t ² while the kinetic energies and spectral transfers behave as t ² and t, respectively. The computed phenomenological coefficients of Youngs' scaling law are compared with experimental data ones. Comparisons with Youngs' three-dimensional numerical simulation (The Physics of Fluids A 3 (1991) 1312) are also performed. Finally three shock tube experiments, where the Richtmyer-Meshkov instability initiates the mixing, are simulated. The mixing thickness evolution is well reproduced while the turbulence levels seem to be overestimated with such first order models. The capability of the two-time-sale model to recover available data for different turbulent flows allows us to conclude to a more universal behavior in comparison with single-time-scale models. This revised version was published online in July 2006 with corrections to the Cover Date.     

Toward a Turbulence Constitutive Relation for Geophysical Flows

Rapidly rotating turbulent flows are frequently in approximate geostrophic balance. Single-point turbulence closures, in general, are not consistent with a geostrophic balance. This article addresses and resolves the possibility of a constitutive relation for single-point second-order closures for classes of rotating and stratified flows relevant to geophysics. Physical situations in which a geostrophic balance is attained are described. Closely related issues of frame-indifference, horizontal divergence, and the Taylor–Proudman theorem are discussed. It is shown that, in the absence of vortex stretching along the axis of rotation, turbulence is frame-indifferent. Unfortunately, no turbulence closures are consistent with this frame-indifference that is frequently an important feature of rotating or quasi-geostrophic flows. A derivation and discussion of the geostrophic constraint which ensures that the modeled second-moment equations are frame-invariant, in the appropriate limit, is given. It is shown that rotating, stratified, and shallow water flows are situations in which such a constitutive relation procedure is useful. A nonlinear nonconstant coefficient representation for the rapid-pressure strain covariance appearing in the Reynolds stress and heat flux equations, consistent with the geostrophic balance, is described. The rapid-pressure strain closure features coefficients that are not constants determined by numerical optimization but are functions of the state of turbulence as parametrized by the Reynolds stresses and the turbulent heat fluxes as is required by tensor representation theory. These issues are relevant to baroclinic and barotropic atmospheric and oceanic flows. The planetary boundary layers in which there is a transition, with height or depth, from a thermally or shear driven turbulence to a geostrophic turbulence is a classic geophysical example to which the considerations in this article are relevant. Received 14 October 1996 and accepted 9 June 1997     

工程湍流模式理论综述及展望

本文讨论了国内外湍流模拟的现状和发展趋势．指出湍流模式的建立除了应遵循理性力学原则外,还必须密切结合工程流体的复杂流动现象,如对具有浮力的回流、分离流及强旋流的模拟,对逆梯度的输运模拟,对单相流、多相流,单流体、多流体的湍流牛顿流体及湍流的非牛顿流体的模拟．应加深对湍流机理的认识,改进湍流模拟手段,结合工程实际,提出较为通用的工程湍流模式．      

Recent progress in the development of the Elliptic Blending Reynolds-stress model

The Elliptic Blending Reynolds Stress Model (EB-RSM), originally proposed by Manceau and Hanjalić (2002) to extend standard, weakly inhomogeneous Reynolds stress models to the near-wall region, has been subject to various modifications by several authors during the last decade, mainly for numerical robustness reasons. The present work revisits all these modifications from the theoretical standpoint and investigates in detail their influence on the reproduction of the physical mechanisms at the origin of the influence of the wall on turbulence. The analysis exploits recent DNS databases for high-Reynolds number channel flows, spanwise rotating channel flows with strong rotation rates, up to complete laminarization, and the separated flow after a sudden expansion without and with system rotation. Theoretical arguments and comparison with DNS results lead to the selection of a recommended formulation for the EB-RSM model. This formulation shows satisfactory predictions for the configurations described above, in particular as regards the modification of the mean flow and turbulent anisotropy on the anticyclonic or pressure side.     

Mixing in a stably stratified shear layer: two- and three-dimensional numerical experiments

A direct method for analyzing diapycnal mixing in a stably stratified fluid (Winters et al., 1995) has been applied to the stably stratified shear layer. The diapycnal flux and mixing efficiency are computed as functions of time, whatever the turbulent activity in the fluid. The mixing properties of two- and three-dimensional numerical simulations of the Boussinesq equations are analyzed and compared. The interest of the former simulations is to emphasize the fundamental role of three-dimensional effects in fluid mixing and to quantify it. We focus on the influence of stratification (measured by the minimum Richardson number J) and changes in Prandtl number on the overall mixing that occurs as the computed flows evolve from unstable initial conditions.

In three dimensions, the flow dynamics exhibit three successive stages, each with different mixing properties. During the first stage, a primarily two-dimensional Kelvin–Helmholtz instability develops and the mixing efficiency is high (the flux Richardson number Rf_b ranges between 0.37 and 0.68, decreasing as J increases). The second stage is characterized by the development of small-scale three-dimensional instabilities. These motions result in significantly higher diapycnal flux than during the first stage but in only moderate mixing efficiency (Rf_b0.32), as the rate of kinetic energy dissipation is also high during this stage. Finally, the turbulent activity is progressively expulsed toward the outer regions of the shear layer and decays in time while the central region relaminarizes. During this final stage, Rf_b approaches an asymptotic value close to 0.25 and the diapycnal diffusivity displays a clear functional dependence on a gradient Richardson number Ri_b of the form Ri_b⁻².

As expected, the two-dimensional flows are unable to reproduce the mixing properties of the flow, except during the first stage. During the subsequent turbulent regime, both the diapycnal flux and the dissipation rate of kinetic energy are too small (because, for the latter quantity, of the nonlinear enstrophy conservation constraint). The final stage consists in a quasi-stationary weakly turbulent regime, for which the diapycnal diffusivity behaves as Ri_b⁻¹. It should be noted that, despite these differences, Rf_b relaxes toward the 0.25 value found in three dimensions.     

Water entrainment due to spillway surface jets

Strong flow entrainment has been observed downstream of spillways constructed with flow deflectors. This water entrainment has important environmental and ecological impacts because it improves the mixing of powerhouse and spillway flows, but may negatively impact fish migration or create adverse flow conditions.

Most studies found in the literature attempt to explain this entrainment with turbulent mixing. Both reduced-scale hydraulic models and single-phase, isotropic RANS models grossly under-predict the degree of entrainment observed in prototypes. In this paper, an anisotropic model that accounts for the bubble volume fraction and attenuation of the normal velocity fluctuations at the free surface is presented. The model adequately predicts the main mechanisms causing water entrainment and compares well against experimental data for a round surface jet and for Brownlee Dam at model scale. It is shown that appropriate entrainment can only be captured if the turbulence anisotropy and the two-phase nature of the flow are modelled.     

Thermal-field structure of a non-premixed turbulent flame formed in a strong pressure-gradient flow

Statistical characteristics of a non-premixed turbulent flame formed in a curved-rectangular duct and spatio-temporal structures of the thermal field were investigated experimentally. The flame was much affected by a strong pressure gradient in the radial direction of the duct curvature, which caused strong gradient diffusion in turbulent heat transfer on the inner-wall side of the flame and, in contrast, counter-gradient heat transfer on the outer-wall side. Two-point correlation measurement of temperature fields revealed that, in the strong gradient diffusion region, a spatial thermal pattern generated by turbulent mixing of high- and low-temperature fluid parcels was advected downstream with little diffusion. In contrast, the pattern was attenuated and diffused rapidly in the counter-gradient diffusion region. These results accurately correspond to the generation mechanism of the counter-gradient heat transport so far observed in stably stratified turbulent flows.     

Advances in studies on two-phase turbulence in dispersed multiphase flows

Particle/droplet/bubble fluctuation and dispersion are important to mixing, heat and mass transfer, combustion and pollutant formation in dispersed multiphase flows, but are insufficiently studied before the 90 years of the last century. In this paper, the present author reports his systematic studies within nearly 20 years on two-phase turbulence in dispersed multiphase flows, including particle fluctuation in dilute gas-particle and bubble-liquid flows, particle-wall collision effect, coexistence of particle turbulence and inter-particle collisions, fluid turbulence modulation due to the particle wake effect and validation of the two-fluid RANS modeling using large-eddy simulation.     

On basic equations of turbulent swirling gas-solid flows and their application in cyclones

The basic equations of turbulent gas-solid flows are derived by using the pseudo-fluid model of particle phase with a refined two-phase turbulence model. These equations are then applied to swirling gas-particle flows for analyzing the collection efficiency in cyclone separators.     

Parallel computation of turbulent flows using moment base lattice Boltzmann method

This paper describes parallel computing approach for simulating turbulent flows using a moment base lattice Boltzmann method. The distribution functions of the lattice Boltzmann method are expressed by corresponding moments. Choosing proper relaxation times for higher order moments, a minimum numerical dissipation is implicitly added to stabilise the method at high Reynolds numbers. Validation of the method is made by computing free decaying periodic turbulent flows and fully developed turbulent channel flows on a GPU platform. Though the present method requires additional work to calculate the higher order moments, it is shown that additional computational cost is negligible in the GPU computing. The numerical results stably obtained for the turbulent flows are in good agreement with those of a pseudo-spectral method and corresponding DNS database.     

Direct numerical simulation of low Reynolds number flows in an open‐channel with sidewalls

A direct numerical simulation of low Reynolds number turbulent flows in an open‐channel with sidewalls is presented. Mean flow and turbulence structures are described and compared with both simulated and measured data available from the literature. The simulation results show that secondary flows are generated near the walls and free surface. In particular, at the upper corner of the channel, a small vortex called inner secondary flows is simulated. The results show that the inner secondary flows, counter‐rotating to outer secondary flows away from the sidewall, increase the shear velocity near the free surface. The secondary flows observed in turbulent open‐channel flows are related to the production of Reynolds shear stress. A quadrant analysis shows that sweeps and ejections are dominant in the regions where secondary flows rush in toward the wall and eject from the wall, respectively. A conditional quadrant analysis also reveals that the production of Reynolds shear stress and the secondary flow patterns are determined by the directional tendency of the dominant coherent structures. Copyright © 2009 John Wiley & Sons, Ltd.     

A new evaluation method for Nusselt numbers in naphthalene sublimation experiments in rotating-disk systems

A comparative analysis of theoretical and empirical relations, as well as experimental data for mass transfer of a rotating disk in laminar, transitional and turbulent flows for naphthalene sublimation in air was done. New correlations between local and average Sherwood numbers for the entire disk were offered. A new evaluation approach for Nusselt numbers based on the experimental data for naphthalene sublimation in laminar, transitional and turbulent flows was developed.     

Interactions of particles and microbubbles with turbulence

Dilute, dispersed two-phase flows arise in many contexts ranging from solid particles or droplets in gas flows to bubbles in liquids. Many of the flows of interest are turbulent, which presents a complex problem to analyze or to determine the dominant physical processes contributing to the observed phenomena. Advances in experimental techniques have made it possible to measure directly turbulent and particle velocity fluctuations in dilute systems. This has provided a counterpart to advances in computational and analytical models and a basis on which to test these models. Three specific areas are considered: the fluctuating forces on an individual particle in an unsteady flow, the response of a solid particle to a turbulent air flow, and the corresponding response of a small bubble in turbulent liquid flows. Results from direct numerical simulations are presented for each of these, including the nonuniform distribution of particles generated by local instantaneous features of the flow. The issue of turbulence modulation at low to moderate void fractions is discussed.     

Experiments to study interactions between baroclinic lower flows and a stably stratified upper layer

Experiments were conducted on a rotating fluid annulus to study the basic interactions between baroclinic lower flows and a stably stratified upper layer. Sufficiently stable stratification is necessary for steady flows to emerge in the lower layer. Upward fluid motions make the baroclinic flows permeate into the upper layer. The stable stratification, however, suppresses upward motions so that zonal fluid velocities decrease with height. In fact, their maximum appears at the top level of the baroclinic lower layer and the sign of the radial temperature gradient changes there; namely, it is warmer on the inner side of the annulus in the upper layer. This temperature profile is reflected in a meridional fluid circulation mixing both layers. In the upper layer of the wave flow, there exists a critical level below and above which the zonal fluid velocities have opposite directions for the wave to have a phase shift of half a wavelength in appearance. The experimental results correspond to real atmospheric phenomena.     

Precessing Vortex Structures in Turbulent Flow within Rotor-Stator Disc Cavities

Several studies of enclosed turbulent flows within rotating discs or cylinders (e.g. [6, 7]) have revealed that, while the geometry may be strictly axisymmetric, it is possible for non-axisymmetric flow patterns to be created within the space. Here we report a visualization study off low induced in the cavity formed between two discs, one rotating, the other stationary. This is an idealization of the flow configuration that occurs between successive stages in the `hot section' of a gas turbine. Such rotor-stator cavities have hitherto been regarded as creating asymmetric flow pattern but Owen [8] has conjectured that the failure to predict heat transfer coefficients accurately for certain radius-to-height ratios may indicate that here, too, organized rotating vortex structures were playing a crucial role. The present study has made an experimental visualization of this flow over a range of conditions in order to test this conjecture and to help guide future numerical explorations. The apparatus comprised a rotating disc over which is fitted a Perspex stationary disc and shroud. The lower disc was rotated for a number of distinct speeds between 30 and 120 rpm and for two ratios of gap-height to radius (H/R). The spin Reynolds number based on gap height and maximum rotational speed, ρΩRh/μ, ranged from 3.7 × 10E4 to 2.24 × 10E5. The flow structures were visualized by injecting ink through a small hypodermic tube at various radii and depths within the cavity and recording the ensuing dye streaks with a video camera mounted above the discs. The results show that, for a wide range of conditions,structured flow with large-scale vortices does indeed arise, the number of vortices diminishing as the spin Reynolds number is increased. The two-vortex S-shaped pattern is stable over a wide range of conditions but three, five and seven vortices have also been observed. These results suggest that an accurate numerical simulation of the flow within rotor-stator disc cavities may require unsteady,three-dimensional CFD modelling over at least certain ranges of flow parameters. This revised version was published online in July 2006 with corrections to the Cover Date.     

Experimental study of stratified jet by simultaneous measurements of velocity and density fields

Stratified flows with small density difference commonly exist in geophysical and engineering applications, which often involve interaction of turbulence and buoyancy effect. A combined particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF) system is developed to measure the velocity and density fields in a dense jet discharged horizontally into a tank filled with light fluid. The illumination of PIV particles and excitation of PLIF dye are achieved by a dual-head pulsed Nd:YAG laser and two CCD cameras with a set of optical filters. The procedure for matching refractive indexes of two fluids and calibration of the combined system are presented, as well as a quantitative analysis of the measurement uncertainties. The flow structures and mixing dynamics within the central vertical plane are studied by examining the averaged parameters, turbulent kinetic energy budget, and modeling of momentum flux and buoyancy flux. At downstream, profiles of velocity and density display strong asymmetry with respect to its center. This is attributed to the fact that stable stratification reduces mixing and unstable stratification enhances mixing. In stable stratification region, most of turbulence production is consumed by mean-flow convection, whereas in unstable stratification region, turbulence production is nearly balanced by viscous dissipation. Experimental data also indicate that at downstream locations, mixing length model performs better in mixing zone of stable stratification regions, whereas in other regions, eddy viscosity/diffusivity models with static model coefficients represent effectively momentum and buoyancy flux terms. The measured turbulent Prandtl number displays strong spatial variation in the stratified jet.