Stochastic Analysis and Modeling of Velocity Observations in Turbulent Flows

Highly turbulent water flows,often encountered near human constructions like bridge piers,spillways,and weirs,display intricate dynamics characterized by the formation of eddies and vortices.These formations,varying in sizes and lifespans,significantly influence the distribution of fluid velocities within the flow.Subsequently,the rapid velocity fluctuations in highly turbulent flows lead to elevated shear and normal stress levels.For this reason,to meticulously study these dynamics,more often than not,physical modeling is employed for studying the impact of turbulent flows on the stability and longevity of nearby structures.Despite the effectiveness of physical modeling,various monitoring challenges arise,including flow disruption,the necessity for concurrent gauging at multiple locations,and the duration of measurements.Addressing these challenges,image velocimetry emerges as an ideal method in fluid mechanics,particularly for studying turbulent flows.To account for measurement duration,a probabilistic approach utilizing a probability density function(PDF)is suggested to mitigate uncertainty in estimated average and maximum values.However,it becomes evident that deriving the PDF is not straightforward for all turbulence-induced stresses.In response,this study proposes a novel approach by combining image velocimetry with a stochastic model to provide a generic yet accurate description of flow dynamics in such applications.This integration enables an approach based on the probability of failure,facilitating a more comprehensive analysis of turbulent flows.Such an approach is essential for estimating both short-and long-term stresses on hydraulic constructions under assessment.

Non-equilibrium turbulent phenomena in transitional flat plate boundary-layer flows

Many recent laboratory experiments and numerical simulations support a non-equilibrium dissipation scaling in decaying turbulence before it reaches an equilibrium state.By analyzing a direct numerical simulation(DNS)database of a transitional boundary-layer flow,we show that the transition region and the non-equilibrium turbulence region,which are located in different streamwise zones,present different non-equilibrium scalings.Moreover,in the wall-normal direction,the viscous sublayer,log layer,and outer layer show different non-equilibrium phenomena which differ from those in grid-generated turbulence and transitional channel flows.These findings are expected to shed light on the modelling of various types of non-equilibrium turbulent flows.

Direct numerical simulation of compressible turbulent flows

This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.

Effects of the Reynolds number on the mean skin friction decomposition in turbulent channel flows

As the Reynolds number increases, the skin friction has been identified as the dominant drag in many practical applications. In the present paper, the effects of the Reynolds number on the mean skin friction decomposition in turbulent channel flows up to Reτ= 5 200 are investigated based on two different methods, i.e., the FukagataIwamoto-Kasagi(FIK) identity(FUKAGATA, K., IWAMOTO, K., and KASAGI, N.Contribution of Reynolds stress distribution to the skin friction in wall-bounded flows.Physics of Fluids, 14(11), L73–L76(2002)) and the Renard-Deck(RD) identity(DECK,S., RENARD, N., LARAUFIE, R., and WEISS, P.′E. Large-scale contribution to mean wall shear stress in high-Reynolds-number flat-plate boundary layers up to Reθ= 13 650.Journal of Fluid Mechanics, 743, 202–248(2014)). The direct numerical simulation(DNS) data provided by Lee and Moser(LEE, M. and MOSER, R. D. Direct numerical simulation of turbulent channel flow up to Reτ≈ 5 200. Journal of Fluid Mechanics,774, 395–415(2015)) are used. For these two skin friction decomposition methods, their decomposed constituents are discussed and compared for different Reynolds numbers.The integrands of the decomposed constituents are locally analyzed across the boundary layer to assess the actions associated with the inhomogeneity and multi-scale nature of turbulent motion. The scaling of the decomposed constituents and their integrands are presented. In addition, the boundary layer is divided into three sub-regions to evaluate the contributive proportion of each sub-region with an increase in the Reynolds number.

Drag reduction of turbulent channel flows over an anisotropic porous wall with reduced spanwise permeability

The direct numerical simulation (DNS) is carried out for the incompressible viscous turbulent flows over an anisotropic porous wall. Effects of the anisotropic porous wall on turbulence modifications as well as on the turbulent drag reduction are investigated. The simulation is carried out at a friction Reynolds number of 180, which is based on the averaged friction velocity at the interface between the porous medium and the clear fluid domain. The depth of the porous layer ranges from 0.9 to 54 viscous units. The permeability in the spanwise direction is set to be lower than the other directions in the present simulation. The maximum drag reduction obtained is about 15.3% which occurs for a depth of 9 viscous units. The increasing of drag is addressed when the depth of the porous layer is more than 25 wall units. The thinner porous layer restricts the spanwise extension of the streamwise vortices which suppresses the bursting events near the wall. However, for the thicker porous layer, the wall-normal fluctuations are enhanced due to the weakening of the wall-blocking effect which can trigger strong turbulent structures near the wall.

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.

Advances in Studies on Turbulent Dispersed Multiphase Flows

Dispersed multiphase flows,including gas-particle(gas-solid),gas-spray,liquid-particle(liquid-solid) ,liquid-bubble,and bubble-liquid-particle flows,are widely encountered in power,chemical and metallurgical,aeronautical and astronautical,transportation,hydraulic and nuclear engineering. In this paper,advances and re-search needs in fundamental studies of dispersed multiphase flows,including the particle/droplet/bubble dynamics,particle-particle,droplet-droplet and bubble-bubble interactions,gas-particle and bubble-liquid turbulence interac-tions,particle-wall interaction,numerical simulation of dispersed multiphase flows,including Reynolds-averaged modeling(RANS modeling),large-eddy simulation(LES) and direct numerical simulation(DNS) are reviewed. The research results obtained by the present author are also included in this review.

LES prediction of space-time correlations in turbulent shear flows

We compare the space-time correlations calculated from direct numerical simulation （DNS） and large-eddy simulation （LES） of turbulent channel flows. It is found from the comparisons that the LES with an eddy-viscosity subgrid scale （SGS） model over-predicts the space-time corre- lations than the DNS. The overpredictions are further quantified by the integral scales of directional correlations and convection velocities. A physical argument for the overpre- diction is provided that the eddy-viscosity SGS model alone does not includes the backscatter effects although it correctly represents the energy dissipations of SGS motions. This argument is confirmed by the recently developed elliptic model for space-time correlations in turbulent shear flows. It suggests that enstrophy is crucial to the LES prediction of spacetime correlations. The random forcing models and stochastic SGS models are proposed to overcome the overpredictions on space-time correlations.

IMPROVED SUBGRID SCALE MODEL FOR DENSE TURBULENT SOLID-LIQUID TWO-PHASE FLOWS

The dense solid-phase governing equations for two-phase flows are obtained by using the kinetic theory of gas molecules.Assuming that the solid-phase velocity distributions obey the Maxwell equations,the collision term for particles under dense two-phase flow conditions is also derived. In comparison with the governing equations of a dilute two-phase flow,the solid-particle's governing equations are developed for a dense turbulent solid-liquid flow by adopting some relevant terms from the dilute two-phase governing equations.Based on Cauchy-Helmholtz theorem and Smagorinsky model, a second-order dynamic sub-grid-scale(SGS)model,in which the sub-grid-scale stress is a function of both the strain-rate tensor and the rotation-rate tensor,is proposed to model the two-phase governing equations by applying dimension analyses.Applying the SIMPLEC algorithm and staggering grid system to the two-phase discretized governing equations and employing the slip boundary conditions on the walls,the velocity and pressure fields,and the volumetric concentration are calculated.The simulation results are in a fairly good agreement with experimental data in two operating cases in a conduit with a rectangular cross-section and these comparisons imply that these models are practical.

MASS TRANSFER IN TURBULENT PULSATING FLOWS

The effect of flow oscillation to the mass transfer between turbulent fluid and solid wall was investigatedby measuring the mass transfer rate between fluid and pipe wall with imposed oscillating flow usingelectrochemical method.The velocity and concentration field in the viscous sublayer which controls the mass trans-fer in such a process was simulated by a simple wave model of single harmonics.Experimental results confirmthat the flow oscillation has no influene on time averaged mass transfer rate,but the phase difference betweenphase averaged velocity field and concentration field shifts with the frequency of imposed oscillating flow.Numeri-cal analysis reveals that the concentration boundarylayer which is responsible for the mass transfer is muchthinner than the viscous sublayer which greatly weakens the influence of imposed oscillating flow on mass transfer.

Basic Laws in Mechanics of Turbulent Flows

Turbulent flow is a basic form of fluid motion widely observed in nature. In hydraulic engineering, especially in the study of sediment movement, turbulence is a key problem. In this paper, based on the stochastic theory of wall turbulence developed by the author and the results by other investigators, fluc-tuation and mean structures and drag coefficient for Newtonian and drag reduction flows in all states (laminar, transitional, turbulent) and in all regions (smooth, transitional, rough) are theoretically discussed in detail. General laws for laminar and turbulent flows obtained by the author are verified by the experimental results obtained by others, and there is good agreement between them.

The dispersion of particles in turbulent semi-circular duct flows

The flow field in a semi-circular duct is simulated by Large Eddy Simulation(LES)and its particle field is simulated by Lagrange particle tracking method.Reynolds number Reb(based on bulk velocity and hydraulic diameter)is 80,000 and Ret(based on friction velocity and hydraulic diameter)is 3528.Particle diameter dpis chosen as 10,50,100,500 mm corresponding to St as 0.10,2.43,9.72,243.05.The results show that the intensity of the secondary flow near the ceiling is less than that near the floor because the ceiling is curved and able to inhibit the secondary flow.It is found that the difference between the semicircular duct and the square duct is that the secondary flow in a corner of the semi-circular duct is not symmetrical along the diagonal although they have the same generation mechanism.Regarding the particles,small particles(dp≤10 mm)are found to uniformly distribute in the duct,while large particles(dp≥50 mm)preferentially distribute in the corner and floor center.The maximum particles(dp=500mm)fall on the floor quickly and their dispersion mainly depends on the secondary flow near the floor.Particle deposition in the corner depends on particle size due to the effect of secondary flow and gravity.The effect of lift force on particles becomes more significant for 50 and 100 mm particles in comparison with other smaller particles.In the end,the effect of secondary flow is found to be more significant to dominate particle behavior than that of flow fluctuation.

NUMERICAL STUDY OF TURBULENT FLOWS OVER WAVY BOUNDARIES

The fully developed túrbulent flows over wavy boundaries are investigated by means of the k-ε model.Predicted flow characteristics over rigid wavy walls are in good agreement with the availa- ble experimental data.Moreover drag reduction has been found in a 2-dimensional channel with periodi- cal wavy walls.The energy input from turbulent wind to regular waves is also studied in the paper by the same turbulence model with carefully posed boundary conditions at wind-wave interface.Better agreement has been obtained in the prediction of the growth rates of wind waves as compared with the previous theoretical and numerical results.

Multi-scale Equations for Incompressible Turbulent Flows

The short-range property of interactions between scales in incompressible turbulent flow was examined. Some formulae for the short-range eddy stress were given. A concept of resonant-range interactions between extremely contiguous scales was introduced and some formulae for the resonant-range eddy stress were also derived. Multi-scale equations for the incompressible turbulent flows were proposed. Key words turbulence - incompressible flow - interactions between scales - multi-scale equations MSC 2000 76F70

Experimental investigation of turbulent flows around high-rise structure foundations and implications on scour

Many studies have been undertaken to predict local scour around offshore high-rise structure foundations(HRSFs),which have been used in constructing the Donghai Wind Farm in China.However,there have been few works on the turbulent flow that drives the scour process.In this study,the characteristics of the turbulent flow fields around an HRSF were investigated using the particle image velocimetry technique.The mean flow,vorticity,and turbulence intensity were analyzed in detail.The relationship between the flow feature and scour development around an HRSF was elaborated.The results showed that the flow velocity increased to its maximum value near the third row of the pile group.The shear layer and wake vortices could not be fully developed downstream of the last row of the piles at small Reynolds numbers.The strong flow and turbulent fluctuation near the third piles explained the existence of a longtail scour pattern starting from the HRSF shoulders and a trapezoidal deposition region directly downstream of HRSF.This laboratory experiment gains insight into the mechanism of the turbulent flow around HRSFs and provides a rare dataset for numerical model verifications.

Shear Flows in the Near-Turbulent Wake Region of High Speed Trains

Two flow cases for scaled high speed train models with different length are numerically analyzed in the framework of the improved delayed detachededdy simulation model.Specific attention is paid to the shear flows and related mechanisms in the near turbulent wake created by these moving models.In particular,a comparative analysis is made on the distributions of turbulent kinetic energy(TKE)and turbulence production(TP)in planes perpendicular to the streamwise direction.The numerical results suggest that,in the wake region very close to the tail,significant TKE and TP can be ascribed to the dynamic interaction between powerful eddies and strong shear,which explain why these quantities are sensitive to the shear strength.The shear flows are essentially governed by the boundary layers developing along the streamwise direction on the train surfaces,especially from the under-body region and the side walls.For other positions located in the downstream direction away from the tail,the interaction of vortices with the non-slip ground serves as a mechanism to promote transfer of energy from weak eddies to turbulence through the shear present in planes parallel to the ground.

On the capability of the curvilinear immersed boundary method in predicting near-wall turbulence of turbulent channel flows

The immersed boundary method has been widely used for simulating flows over complex geometries.However,its accuracy in predicting the statistics of near-wall turbulence has not been fully tested.In this work,we evaluate the capability of the curvilinear immersed boundary(CURVIB)method in predicting near-wall velocity and pressure fluctuations in turbulent channel flows.Simulation results show that quantities including the time-averaged streamwise velocity,the rms(root-mean-square)of velocity fluctuations,the rms of vorticity fluctuations,the shear stresses,and the correlation coefficients of u'and v"computed from the CURVIB simulations are in good agreement with those from the body-fitted simulations.More importantly,it is found that the time-averaged pressure,the rms and wavenumber-frequency spectra of pressure fluctuations computed using the CURVIB method agree well with the body-fitted results.

REYNOLDS-STRESS MODELLING OF TURBULENT ROTATING FLOWS

The present analysis shows that the EVM can not reflect the turbulence physics in non-inertial frame. The effects of Coriolis force on turbulence is embodied naturally in the Reynolds-stress transport equation. It is observed that the existing second-moment closure model with appropriate near-wall treatment can adequately predict flows in rotating channel and in axially rotating pipe for moderate rotation rate.

RESEARCH ON THE EFFECT OF CYLINDER PARTICLES ON THE TURBULENT PROPERTIES IN PARTICULATE FLOWS

The fluid fluctuating velocity equations which include the term of cylinder particles were established. The turbulent intensity and Reynolds stress of fluid were obtained by averaging fluctuating velocity based on the solution of the fluctuating velocity equations. Above approach was used to solve the channel turbulent flows, and computational results were compared with the experimental ones for the case of single phase flow. The effects of volume fraction of particles, the ratio of particle length to diameter and the particle relaxation time on turbulent properties were illustrated by changing cylinder particle parameters. It is shown that particles play a restraining role to turbulent properties in the flows. The degree of restraint is directly proportional to the volume fraction of particle, the ratio of particle length to diameter and inversely proportional to particle relaxation time.

THE NUMERICAL SIMULATION OF TURBULENT FLOWS OF NEWTONIAN AND A SORT OF NON-NEWTONIAN FLUID IN 180-DEGREE SQUARE SECTIONED BEND

Using k- model of turbulence and measured wall functions, turbulent flows of Newtonian (pure water) andasort of non-Newtonian fluid (dilute, drag-reduction solution of polymer) in a 180-degree curved bend were simulated numerically. The calculated results agreed well with the measured velocity profiles. On the basis of calculation and measurement, appropriateness of turbulence model to complicated flow in which the large-scale vortex exists was analyzed and discussed.