Numerical simulation of a combined oxidation ditch flow using 3D k-ε turbulence model

The standard three dimensional（3D） k-ε turbulence model was applied to simulate the flow field of a small scale combined oxidation ditch. The moving mesh approach was used to model the rotor of the ditch. Comparison of the computed and the measured data is acceptable. A vertical reverse flow zone in the ditch was found, and it played a very important role in the ditch flow behavior. The flow pattern in the ditch is discussed in detail, and approaches are suggested to improve the hydrodynamic performance in the ditch.

Turbulence,aeration and bubble features of air-water flows in macro-and intermediate roughness conditions

Free surface flows aeration potential causing the in macro- and intermediate flow characteristics to vary roughness conditions have a high with slopes and discharges. The underlying mechanism of two-phase flow characteristics in macro- and intermediate roughness conditions were analyzed in an experimental setup assembled at the Laboratory of Hydraulic Protection of the Territory （PITLAB） of the University of Pisa, Italy. Crushed angular rocks and hemispherical boulders were used to intensify the roughness of the bed. Flow rates per unit width ranging between 0.03 m^2/s and 0.09 m^2/s and slopes between 0.26 and 0.46 were tested over different arrangements of a rough bed. Analyses were mainly carried out in the inner flow region, which consists of both bubbly and intermediate flow regions. The findings revealed that the two-phase flow properties over the rough bed were much affected by rough bed arrangements. Turbulence features of two-phase flows over the rough bed were compared with those of the stepped chute data under similar flow conditions. Overall, the results highlight the flow features in the inner layers of the two-phase flow, showing that the maximum turbulence intensity decreases with the relative submergence, while the bubble frequency distribution is affected by the rough bed elements.

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.

A dual-scale turbulence model for gas–liquid bubbly flows

A dual-scale turbulence model is applied to simulate cocurrent upward gas-liquid bubbly flows and validated with available experimental data. In the model, liquid phase turbulence is split into shear-induced and bubble- induced turbulence. Single-phase standard k-e model is used to compute shear-induced turbulence and another transport equation is added to model bubble-induced turbulence. In the latter transport equation, energy loss due to interface drag is the production term, and the characteristic length of bubble-induced turbulence, simply the bubble diameter in this work, is introduced to model the dissipation term. The simulated results agree well with experimental data of the test cases and it is demonstrated that the proposed dual-scale turbulence model outperforms other models. Analysis of the predicted turbulence shows that the main part of turbulent kinetic en- ergy is the bubble-induced one while the shear-induced turbulent viscosity predominates within turbulent vis- cosity, especially at the pipe center. The underlying reason is the apparently different scales for the two kinds of turbulence production mechanisms： the shear-induced turbulence is on the scale of the whole pipe while the bubble-induced turbulence is on the scale of bubble diameter. Therefore, the model reflects the multi-scale phe- nomenon involved in gas-liquid bubbly flows.

Effects of sawtooth heat pulses on edge flows and turbulence in a tokamak plasma

Enhancements of edge zonal flows,radial electric fields,and turbulence are observed in electron cyclotron resonance heating-heated plasmas(Zhao et al 2013 Nucl.Fusion 53083011).In this paper,the effects of sawtooth heat pulses on flows and turbulence are presented.These experiments are performed using multiple Langmuir probe arrays in the edge plasmas of the HL-2A tokamak.The edge zonal flows,radial electric fields,and turbulence are all enhanced by sawteeth.Propagation of the zonal flow and turbulence intensities is also observed.The delay time of the maximal intensity of the electric fields,zonal flows,and turbulence with respect to the sawtooth crashes is estimated as~1 ms and comparable to that of the sawtooth-triggered intermediate phases.Not only the zonal flows but also the radial electric fields lag behind the turbulence.Furthermore,the intensities of both the zonal flows and electric fields nearly linearly increase/decrease with the increase/decrease of the turbulence intensity.A double-source predator-prey model analysis suggests that a relatively strong turbulence source may contribute to the dominant zonal flow formation during sawtooth cycles.

Turbulence modulation model for gas–particle flow based on probability density function approach

he paper focuses on the turbulence modulation problem in gas–particle flow with the use of probability density function（PDF） approach. By means of the PDF method, a general statistical moment turbulence modulation model without considering the trajectory difference between two phases is derived from the Navier–Stokes equations. A new turbulence production term induced by the dispersed-phase is analyzed and considered. Furthermore, the trajectory difference between two media is taken into account. Subsequently, a new k–ε turbulence modulation model in dilute particle-laden flow is successfully set up. Then, the changes to several terms, including the turbulence production, dissipation, and diffusion terms, are well described consequently. The promoted model provides a more probable explanation for the modification of particles on the turbulence. Finally, we applied the model to simulate a gas–particle turbulence flow case in a wall jet, and found that the simulation results agree well with the experimental data.

NUMERICAL SIMULATION OF WING-BODY JUNCTION TURBULENCE FLOW

Wing-body junction turbulence flow is simulated by using RANS equation and boundary fitted coordinate technique. Three order differential scheme is used in the computation of convection term and two layers turbulence model are employed in the calculation.

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.

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.

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.

An improved large eddy simulation of two-phase flows in a pump impeller

An improved large eddy simulation using a dynamic second-order sub-grid-scale （SGS） stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating coordinate system, and continuity is conserved by a mass-weighted method to solve the filtered governing equations. In the cur- rent second-order SGS model, the SGS stress is a function of both the resolved strain-rate and rotation-rate tensors, and the model parameters are obtained from the dimensional consistency and the invariants of the strain-rate and the rotation-rate tensors. In the numerical calculation, the finite volume method is used to discretize the governing equations with a staggered grid system. The SIMPLEC algorithm is applied for the solution of the discretized governing equations. Body- fitted coordinates are used to simulate the two-phase flows in complex geometries. Finally the second-order dynamic SGS model is successfully applied to simulate the dense turbu-lent particle-liquid two-phase flows in a centrifugal impeller. The predicted pressure and velocity distributions are in good agreement with experimental results.

Effect of Agitator’s Types on the Hydrodynamic Flow in an Agitated Tank

The aim of this paper is to study the effect of agitator’s types on the turbulent flows in stirred tanks without and with baffles.The hydrodynamics behavior induced by four different agitator’s types:a Rushton turbine(RT),a circular blade turbine(CBT),a diverging triangular blade turbine(DTBT)and converging triangular blade turbine(CTBT)are numerically predicted by solving the Navier-Stokes equations and RNGκ–εturbulent model.The simulations are carried out using the Multi Reference Frame(MRF)approach.The numerical results showed good agreement with experiment.We find that the agitator CTBT gives an important profit on the power consumption per report/ratio the others and DTBT give a good reduction of the vortex size of the impeller angles.

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.

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.

Characteristics of turbulent flow distribution in branch piping system

Flow distribution in branch piping system is affected by flow characteristics and different geometric variations. Most of the flow distribution studies are performed with one-dimensional analysis to yield overall information only. However, detailed analysis is required to find effects of design parameters on the flow distribution. For this aspect, three-dimensional turbulent flow analysis was performed to assess turbulence model performance and effects of upstream pressure and branch pipe geometry. Three different turbulence models of standard k-e model, realizable k-e model and standard k-co yield similar results, indicating small effects of turbulence models on flow characteristics analysis. Geometric variations include area ratio of main and branch pipes, branch pipe diameter, and connection shape of main and branch pipes. Among these parameters, area ratio and branch diameter and shape show strong effect on flow distribution due to high friction and minor loss. Uniform flow distribution is one of common requirements in the branch piping system and this can be achieved with rather high total loss design.

EXPERIMENTAL STUDY ON TURBULENT FEATURES IN THE NEGATIVE TRANSPORT REGION OF ASYMMETRIC PLANE CHANNEL FLOW

Turbulent features of streamwise and vertical components of velocity in the negative transport region of asymmetric plane channel flow have been studied experimentally in details. Experiments show that turbulent fluctuations in negative transport region are suppressed, and their probability distributions are far from Gaussian. Besides, the skewness factors attain their negative maxima at the position of the maximum mean velocity, whereas the flatness factors attain their positive maxima at the same position.

NUMERICAL SIMULATION FOR INCOMPRESSIBLE TURBULENT FLOW IN THE IMPELLER OF SEWAGE PUMP

The 3-D turbulent flow in the impeller of sewage pump is simulated. Thetime-averaged N-S equations and the kappa-epsilon turbulent model is modified. The calculation iscarried out in body-fitted coordinated grid by applying SIMPLE-C algorithm. The calculated velocity,pressure distributions of the turbulent flow in the sewage pump are obtained for the first time,which will be helpful for the optimal design and performance prediction of sewage pumps on the basisof flow field simulation.

Assessment of k–ε models using tetrahedral grids to describe the turbulent flow field of a PBT impeller and validation through the PIV technique

In turbulence modeling, the RNG and Realizable models have important improvements in the turbulent production and dissipation terms in comparison to the Standard. The selection of the appropriate turbulence model has an impact on the convergence and solution in STRs, and they are used in mixing, multiphase modeling or as starting solution of transient models as DES and LES. Although there are several studies with the pitched blade turbine(PBT) impeller, most of them used the Standard model as representative of all k–ε models, using structured hexahedral grids composed of low number of cells, and in some cases under axial symmetry assumptions.Accordingly, in this work the assessment of the Standard, RNG and Realizable models to describe the turbulent flow field of this impeller, using the Multiple Reference Frame(MRF) and Sliding Mesh(SM) approaches with tetrahedral domains in dense grids, is presented. This kind of cell elements is especially suitable to reproduce complex geometries. Flow velocities and turbulent parameters were verified experimentally by PIV and torque measurements. The three models were capable of predicting fairly the pumping number, the power number based on torque, and velocities. Although the RNG improved the predictions of the turbulent kinetic energy and dissipation rate, the Realizable model presented better performance for both approaches. All models failed in the prediction of the total dissipation rate, and a dependence of its value on the number of cells for the MRF was found.

Effect of Downward Seepage on Turbulent Flow Characteristics and Bed Morphology around Bridge Piers

In this work, experimental investigations have been pursued to analyse the influence of downward seepage on the turbulent characteristics of flow and corresponding changes in vortex structure around circular bridge pier in alluvial channel. Experiments were conducted in sand bed channel with circular piers of different sizes for no seepage, 10% seepage and 20% seepage cases. The measurement of turbulent flow statistics such as velocity and Reynolds stresses is found to be negative within the scour hole at upstream of the pier whereas application of downward seepage retards the reversal of the flow causing a decrement in the velocity and Reynolds stresses. Higher Reynolds shear stress prevails at the downstream side because of the production of wake vortices. Contribution of all bursting events to the total Reynolds shear stress production has been observed to increase with downward seepage. The analysis of integral scale suggest that size of eddies increases with seepage, which is responsible for increase in particle mobility. Initially rate of scouring is more which abatements gradually with expanding time as well as with the increased of downward seepage. Presence of downward seepage reduces the depth and length of vortex and shifts towards downstream side of the pier.

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