Structural ensemble dynamics based closure model for wall-bounded turbulent flow

Wall-bounded turbulent flow involves the development of multi-scale turbulent eddies, as well as a sharply varying boundary layer. Its theoretical descriptions are yet phenomenological. We present here a new framework called structural ensemble dynamics （SED）, which aims at using systematically all relevant statistical properties of turbulent structures for a quantitative description of ensemble means. A new set of closure equations based on the SED approach for a turbulent channel flow is presented. SED order functions are defined, and numerically determined from data of direct numerical simulations （DNS）. Computational results show that the new closure model reproduces accurately the solution of the original Navier-Stokes simulation, including the mean velocity profile, the kinetic energy of the streamwise velocity component, and every term in the energy budget equation. It is suggested that the SED-based studies of turbulent structure builds a bridge between the studies of physical mechanisms of turbulence and the development of accurate model equations for engineering predictions.

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.

Large-eddy simulation of wall-bounded turbulent flow with high-order discrete unified gas-kinetic scheme

In this paper,we introduce the discrete Maxwellian equilibrium distribution function for incompressible flow and force term into the two-stage third-order Discrete Unified Gas-Kinetic Scheme(DUGKS)for simulating low-speed turbulent flows.The Wall-Adapting Local Eddy-viscosity(WALE)and Vreman sub-grid models for Large-Eddy Simulations(LES)of turbulent flows are coupled within the present framework.Meanwhile,the implicit LES are also presented to verify the effect of LES models.A parallel implementation strategy for the present framework is developed,and three canonical wall-bounded turbulent flow cases are investigated,including the fully developed turbulent channel flow at a friction Reynolds number(Re)about 180,the turbulent plane Couette flow at a friction Re number about 93 and lid-driven cubical cavity flow at a Re number of 12000.The turbulence statistics,including mean velocity,the r.m.s.fluctuations velocity,Reynolds stress,etc.are computed by the present approach.Their predictions match precisely with each other,and they are both in reasonable agreement with the benchmark data of DNS.Especially,the predicted flow physics of three-dimensional lid-driven cavity flow are consistent with the description from abundant literature.The present numerical results verify that the present two-stage third-order DUGKS-based LES method is capable for simulating inhomogeneous wall-bounded turbulent flows and getting reliable results with relatively coarse grids.

AModeling Study on ParticleDispersion inWall-Bounded Turbulent Flows

Three physicalmechanismswhichmay affect dispersion of particle’smotion in wall-bounded turbulent flows,including the effects of turbulence,wall roughness in particle-wall collisions,and inter-particle collisions,are numerically investigated in this study.Parametric studies with differentwall roughness extents and with different mass loading ratios of particles are performed in fully developed channel flows with the Eulerian-Lagrangian approach.A low-Reynolds-number k−ǫturbulence model is applied for the solution of the carrier-flow field,while the deterministic Lagrangian method together with binary-collision hard-sphere model is applied for the solution of particle motion.It is shown that the mechanism of inter-particle collisions should be taken into account in the modeling except for the flows laden with sufficiently low mass loading ratios of particles.Influences of wall roughness on particle dispersion due to particle-wall collisions are found to be considerable in the bounded particleladen flow.Since the investigated particles are associated with large Stokes numbers,i.e.,larger than O(1),in the test problem,the effects of turbulence on particle dispersion aremuch less considerable,as expected,in comparison with another two physical mechanisms investigated in the study.

Turbine Passage Secondary Flow Dynamics and Endwall Heat Transfer Under Different Inflow Turbulence

This study presents endwall hydrodynamics and heat transfer in a linear turbine cascade at Re 5×105 at low and high intensities of turbulence.Results are numerically predicted using the standard SST model and Reθ-γtransition model as well as using the high-resolution LES separately.The major secondary flow components,comprising the horseshoe,corner,and passage vortices are recognized and the impact on heat or mass transfer is investigated.The complicated behavior of turbine passage secondary flow generation and establishment are impacted by the perspective of boundary layer attributes and inflow turbulence.The passage vortex concerning the latest big leading-edge vane is generated by the enlargement of the circulation developed at the first instance adjacent to the pressure side becomes powerful and mixes with other vortex systems during its migration towards the suction side.The study conclusions reveal that substantial enhancements are attained on the endwall surface,for the entire spanwise blade extension on the pressure surface,and in the highly 3-D region close to the endwall on the suction surface.The forecasted suction surface thermal exchange depicts great conformity with the measurement values and precisely reproduces the enhanced thermal exchange owing to the development and lateral distribution of the secondary flows along the midspan of the blade passage downstream.The impacts of the different secondary flow structures on the endwall thermal exchange are described in depth.

Analysis of particle dispersion in a turbulent flow considering particle rotation

Non-spherical particles exist widely in natural and industrial fluid systems and the motions of nonspherical particles are significantly different from that of spherical particles.In this paper,a simplified model of non-spherical particles considering particle drag correction,lift,and rotation was established.Based on the Eulerian-Lagrangian simulation,the dispersion characteristics of spherical and nonspherical particles with different Stokes numbers in a high-speed turbulent jet were analyzed and compared considering the effect of particle rotation.The results show that,the differences in particle dispersion and radial velocity fluctuation between non-spherical particles and spherical particles in the jet are significant,especially when Stokes number is large.Moreover,the effects of different type of forces on the dispersion of non-spherical particles and spherical particles were compared in detail,which revealed that the change of the Magnus force caused by the increase in the angular velocity of non-spherical particles plays a dominant role in the differences of particle dispersion.

Observed characteristics of flow,water mass,and turbulent mixing in the Preparis Channel

Preparis Channel is the very important exchange path of energy and materials between the northern Bay of Bengal and Andaman Sea(AS).A set of hydrographic measurements,a microstructure profiler,and a deep mooring were used to determine the characteristics of water masses,turbulent mixing,and flows in the Preparis Channel.The unprecedented short-term mooring data reveal that a deep current in the deep narrow passage(below 400 m)of the Preparis Channel flows toward the Bay of Bengal(BoB)with a mean along-stream velocity of 25.26 cm/s at depth of 540 m;above the deep current,there are a relatively weak current flows toward the AS with a mean along-stream velocity of 15.46 cm/s between 500 m and 520 m,and another weak current flows toward the BoB between 430 m and 500 m.Thus,a sandwiched vertical structure of deep currents(below 400 m)is present in the Preparis Channel.The volume transport below 400 m is 0.06 Sv(1 Sv=106 m^(3)/s)from the AS to the BoB.In the upper layer(shallower than 300 m),the sea water of the AS is relatively warmer and fresher than that in the BoB,indicating a strong exchange through the channel.Microstructure profiler observations reveal that the turbulent diffusivity in the upper layer of the Preparis Channel reaches O(10−4 m^(2)/s),one order larger than that in the interior of the BoB and over the continental slope of the northern AS.We speculate that energetic high-mode internal tides in the Preparis Channel contribute to elevated turbulent mixing.In addition,a local“hotspot”of turbidity is identified at the deep mooring site,at depth of about 100 m,which corresponds to the location of elevated turbulent mixing in the Preparis Channel.

Artificial neural network-based subgrid-scale models for LES of compressible turbulent channel flow

Fully connected neural networks(FCNNs)have been developed for the closure of subgrid-scale(SGS)stress and SGS heat flux in large-eddy simulations of compressible turbulent channel flow.The FCNNbased SGS model trained using data with Mach number Ma=3.0 and Reynolds number Re=3000 was applied to situations with different Mach numbers and Reynolds numbers.The input variables of the neural network model were the filtered velocity gradients and temperature gradients at a single spatial grid point.The a priori test showed that the FCNN model had a correlation coefficient larger than 0.91 and a relative error smaller than 0.43,with much better reconstructions of SGS unclosed terms than the dynamic Smagorinsky model(DSM).In a posteriori test,the behavior of the FCNN model was marginally better than that of the DSM in predicting the mean velocity profiles,mean temperature profiles,turbulent intensities,total Reynolds stress,total Reynolds heat flux,and mean SGS flux of kinetic energy,and outperformed the Smagorinsky model.

Simulation Analysis of New Energy Vehicle Engine Cooling System Based on K-E Turbulent Flow Mathematical Model

New energy vehicles have better clean and environmental protection characteristics than traditional fuel vehicles.The new energy engine cooling technology is critical in the design of new energy vehicles.This paper used oneand three-way joint simulation methods to simulate the refrigeration system of new energy vehicles.Firstly,a k-εturbulent flow model for the cooling pump flow field is established based on the principle of computational fluid dynamics.Then,the CFD commercial fluid analysis software FLUENT is used to simulate the flow field of the cooling pump under different inlet flow conditions.This paper proposes an optimization scheme for new energy vehicle engines’“boiling”phenomenon under high temperatures and long-time climbing conditions.The simulation results show that changing the radiator’s structure and adjusting the thermostat’s parameters can solve the problem of a“boiling pot.”The optimized new energy vehicle engine can maintain a better operating temperature range.The algorithm model can reference each cryogenic system component hardware selection and control strategy in the new energy vehicle’s engine.

On the generation of drift flows in wall-bounded flows transiting to turbulence

Despite recent progress, laminar-turbulent coexistence in transitional planar wall-bounded shear flows is still not well understood. Contrasting with the processes by which chaotic flow inside turbulent patches is sustained at the local(minimal flow unit) scale, the mechanisms controlling the obliqueness of laminar-turbulent interfaces typically observed all along the coexistence range are still mysterious. An extension of Waleffe's approach [Waleffe, 1997] is used to show that,already at the local scale, drift flows breaking the problem's spanwise symmetry are generated just by slightly detuning the modes involved in the self-sustainment process. This opens perspectives for theorizing the formation of laminar-turbulent patterns.

STUDY ON TURBULENT SPOTS IN PLANE COUETTE FLOW

Turbulent spots play a key role in the formation of the turbulence and the transition. The generation and evolution of turbulent spots using the wall impulse model in the plane Couette flow are studied by direct numerical simulation of Navier-Stokes equations. A group of three-dimensional coupling compact difference schemes with high accuracy and high resolution is used in the numerical calculation. The important characteristics of turbulent spots based on the results of examples are analyzed, including the formation of random pulse, the generation of Reynolds stress, the growth of disturbance amplitude, and the continuous change of spot shape, especially the complex evolution process of the streamwise vortices. Computational results confirm that basic properties of turbulent spots in the laminar flow are similar to those in the turbulent flow.

Numerical simulation of gas–liquid flow in the bubble column using Wray–Agarwal turbulence model coupled with population balance model

In this paper,an improved computational fluid dynamic(CFD)model for gas-liquid flow in bubble column was developed using the one-equation Wary-Agarwal(WA)turbulence model coupled with the population balance model(PBM).Through 18 orthogonal test cases,the optimal combination of interfacial force models,including drag force,lift force,turbulent dispersion force.The modified wall lubrication force model was proposed to improve the predictive ability for hydrodynamic behavior near the wall of the bubble column.The values simulated by optimized CFD model were in agreement with experimental data,and the errors were within±20%.In addition,the axial velocity,turbulent kinetic energy,bubble size distribution,and the dynamic characteristic of bubble plume were analyzed at different superficial gas velocities.This research work could provide a theoretical basis for the extension of the CFD-PBM coupled model to other multiphase reactors..

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.

NUMERICAL SIMULATION OF 3-D DENSE SOLID-LIQUID TWO-PHASE TURBULENT FLOW IN A NON-CLOGGING MUD PUMP

A mathematical model is set to evaluate the 3-D dense solid-liquid two-phaseturbulent flow in a non-clogging mud pump, the flow feature in the impeller channel is simulatedwith the tool of IPSA. Meanwhile, resort to TECPLOT as the post-processor, the simulation results isvisualized. The results show the main flow characteristics: There exists backflow and aberrantvelocities at inlet area and a relative velocity slip between two phases; A jet-wake flow pattern isdiscerned around the shroud-suction side area; The relative velocity vector of solid phase iscloser to the pressure surface than that of liquid phase and the trend is more obvious with theincrease of diameter; The kinetic energy of turbulence k and the dissipation rate e reach theirpeaks at the corner of pressure and suction surface. The simulation results show a good agreementwith the experimental flow features in the impeller channel, which prove the turbulent model used isvalid and provide a theoretical design basis to non-clogging pumps.

Re-understanding the law-of-the-wall for wall-bounded turbulence based on in-depth investigation of DNS data

Based on direct numerical simulation (DNS) data of the straight ducts,namely square and rectangular annular ducts,detailed analyses were conducted for the mean streamwise velocity,relevant velocity scales,and turbulence statistics.It is concluded that turbulent boundary layers (TBL) should be broadly classified into three types (Type-A,-B,and-C) in terms of their distribution patterns of the time-averaged local wall-shear stress (τw) or the mean local frictional velocity (uτ).With reference to the Type-A TBL analysis by von Karman in developing the law-of-the-wall using the time-averaged local frictional velocity (uτ) as scale,the current study extended the approach to the Type-B TBL and obtained the analytical expressions for streamwise velocity in the inner-layer using ensemble-averaged frictional velocity (ūτ) as scale.These analytical formulae were formed by introducing the general damping and enhancing functions.Further,the research applied a near-wall DNS-guided integration to the governing equations of Type-B TBL and quantitatively proved the correctness and accuracy of the inner-layer analytical expressions for this type.

Numerical Simulation on Gas-Solid Two-Phase Turbulent Flow in FCC Riser Reactors(Ⅰ) Turbulent Gas-Solid Flow-Reaction Model

Gas-solid two-phase turbulent flows,mass transfer,heat transfer and catalytic cracking reactions areknown to exert interrelated influences in commercial fluid catalytic cracking(FCC)riser reactors.In the presentpaper,a three-dimensional turbulent gas-solid two-phase flow-reaction model for FCC riser reactors was devel-oped.The model took into account the gas-solid two-phase turbulent flows,inter-phase heat transfer,masstransfer,catalytic cracking reactions and their interrelated influence.The k-V-k_P two-phase turbulence modelwas employed and modified for the two-phase turbulent flow patterns with relatively high particle concentration.Boundary conditions for the flow-reaction model were given.Related numerical algorithm was formed and a nu-merical code was drawn up.Numerical modeling for commercial FCC riser reactors could be carried out with thepresented model.

The investigation and optimization of drag reduction in turbulent flow of Newtonian fluid passing through horizontal pipelines using functionalized magnetic nanophotocatalysts and lecithin

The nanophotocatalysts were synthesized in four stages and evaluated by FTIR, FESEM and VSM analysis. The influence of nanofluids containing functionalized magnetic Ti O2 nanophotocatalyst and dipalmitoylphosphatidylcholine lecithin in drag reduction of turbulent flow in four horizontal pipelines was studied. The effective parameters on drag reduction(nanoparticle concentration, surfactant concentration, p H and Re number) were investigated and optimized in each pipeline using response surface method. The drag reduction in 1/2 " galvanized, 3/4" galvanized, 1/2 "five-layer and 1/2" cuprous pipelines was found 99.1%, 92.5%, 87.6% and 85.2%, respectively. The model adequacy was measured using ANOVA. Based on the high determination coefficient, more than 95% of variance of experimental data in all pipelines was described by quadratic model.

Numerical investigation of velocity distribution of turbulent flow through vertically double-layered vegetation

The velocity structures of flow through vertically double-layered vegetation(VDLV)as well as single-layered rigid vegetation(SLV)were investigated computationally with a three-dimensional(3D)Reynolds stress turbulence model,using the computational fluid dynamics(CFD)code FLUENT.The detailed velocity distribution was explored with a varying initial Froude number(Fr),with consideration of the steady subcritical flow conditions of an inland tsunami.In VDLV flows,the numerical model successfully captured the inflection point in the profiles of mean streamwise velocities in the mixing-layer region around the top of short submerged vegetation.An upward and downward movement of flow occurred at the positions located just behind the tall and short vegetation,respectively.Overall,higher streamwise velocities were observed in the upper vegetation layer due to high porosity,with Pr=98%(sparse vegetation,where Pr is the porosity),as compared to those in the lower vegetation layer,which had comparatively low porosity,with Pr=91%(dense vegetation).A rising trend of velocities was found as the flow passed through the vegetation region,followed by a clear sawtooth distribution,as compared to the regions just upstream and downstream of vegetation where the flow was almost uniform.In VDLV flows,a rising trend in the flow resistance was observed with the increase in the initial Froude number,i.e.,Fr?0.67,0.70,and 0.73.However,the flow resistance in the case of SLV was relatively very low.The numerical results also show the flow structures within the vicinity of short and tall vegetation,which are difficult to attain through experimental measurements.

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.

Study of Turbulent Flow of Film Cooling Holes with Lateral Expanded Exits

Hot wire measurements and flow visualization are presented for studying the turbulent flow field over a flat gas turbine film cooling blade with lateral expanded holes. Three mass flux ratios of jet to free stream, M = 0.5, 0.89, 1.5, are tested. The streamwise velocity, the turbulent intensities and the Reynolds shear stress are measured. The effect of the lateral expanded holes on the improvement of the turbulent flow field for film cooling of gas turbines can be analyzed from the measured spatial di...