Numerical prediction of inner turbulent flow in conical diffuser by using a new five-point scheme and DLR k-ε turbulence model

The internal turbulent flow in conical diffuser is a very complicated adverse pressure gradient flow.DLR k-ε turbulence model was adopted to study it.The every terms of the Laplace operator in DLR k-ε turbulence model and pressure Poisson equation were discretized by upwind difference scheme.A new full implicit difference scheme of 5-point was constructed by using finite volume method and finite difference method.A large sparse matrix with five diagonals was formed and was stored by three arrays of one dimension in a compressed mode.General iterative methods do not work wel1 with large sparse matrix.With algebraic multigrid method(AMG),linear algebraic system of equations was solved and the precision was set at 10-6.The computation results were compared with the experimental results.The results show that the computation results have a good agreement with the experiment data.The precision of computational results and numerical simulation efficiency are greatly improved.

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.

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.

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.

Effects of Noncircular Inlet on the Flow Structures in Turbulent Jets

Turbulent jet flows with noncircular nozzle inlet are investigated by using a Reynolds Stress Model. In order to analyze the effects of noncircular inlet, the cross section of inlet are selected as circular, square, and equilateral triangular shape. The jet half-width, vorticity thickness, and developments of the secondary flow are presented. From the result, it is confirmed that the secondary flows of square and equilateral triangular nozzle are more vigorous than that of the circular jet. This development of secondary flows is closely related to the variations of vortical motions in axial and azimuthal directions.

Diffusion Process of High Concentration Spikes in a Quasi-Homogeneous Turbulent Flow

When a mass spreads in a turbulent flow, areas with obviously high concentration of the mass compared with surrounding areas are formed by organized structures of turbulence. In this study, we extract the high concentration areas and investigate their diffusion process. For this purpose, a combination of Planar Laser Induced Fluorescence (PLIF) and Particle Image Velocimetry (PIV) techniques was employed to obtain simultaneously the two fields of the concentration of injected dye and of the velocity in a water turbulent channel flow. With focusing on a quasi-homogeneous turbulence in the channel central region, a series of PLIF and PIV images were acquired at several different downstream positions. We applied a conditional sampling technique to the PLIF images to extract the high concentration areas, or spikes, and calculated the conditional-averaged statistics of the extracted areas such as length scale, mean concentration, and turbulent diffusion coefficient. We found that the averaged length scale was constant with downstream distance from the diffusion source and was smaller than integral scale of the turbulent eddies. The spanwise distribution of the mean concentration was basically Gaussian, and the spanwise width of the spikes increased linearly with downstream distance from the diffusion source. Moreover, the turbulent diffusion coefficient was found to increase in proportion to the spanwise distance from the source. These results reveal aspects different from those of regular mass diffusion and let us conclude that the diffusion process of the spikes differs from that of regular mass diffusion.

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.

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.

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.

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.

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

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.

Numerical Simulation of Turbulent Flow of Hydraulic Oil through 90° Circular-sectional Bend

Oil flow through pipe bends is found in many engineering applications. However, up to now, the studies of oil flow field in the pipe bend appear to be relatively sparse, although the oil flow field and the associated losses of pipe bend are very important in practice. In this paper, the relationships between the turbulent flow of hydraulic oil in a bend and the Reynolds number Re and the curvature ratio δare studied by using computational fluid dynamics （CFD）. A particular emphasis is put on hydraulic oil, which differs from air or water, flowing through 90° circular-sectional bend, with the purpose of determining the turbulent flow characteristics as well as losses. Three turbulence models, namely, RNG κ-ε model, realizable k-ε model, and Reynolds stress model （RSM）, are used respectively. The simulation results in the form of contour and vector plots for all the three turbulence models for pipe bends having curvature ratio of δ=0.5, and the detailed pressure fields and total pressure losses for different Re and δ for RSM are presented. The RSM can predict the stronger secondary flow in the bend better than other models. As Re increases, the pressure gradient changes rapidly, and the pressure magnitude increases at inner and outer wall of the bend. When δ decreases, two transition points or transition zones of pressure gradient arise at inner wall, meanwhile, the transition point moves towards the inlet at outer wall of the bend. Owing to secondary flow, the total pressure loss factor k increases as the bend tightens, on the contrary, as Re increases, factor k decreases due to higher velocity heads, and the rapid change of pressure gradient on the surface of the bend leads to increasing of friction and separation effects, and magnified swirl intensity of secondary flow. A new mathematical model is proposed for predicting pressure loss in terms of Re and δ in order to provide support to the one-dimensional simulation software. The proposed research provides reference for the analysis of oil flow with higher Re in the large bends.

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.

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.

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.

Numerical simulation on turbulent flow field in convergentdiveroent nozzle

Because of the complication of turbulence＇s mechanism and law as well as the jet pressure in nozzle is difficult to test by experiment, five turbulent models were applied to numerically simulate the turbulent flow field in convergent-divergent nozzle. Theory analysis and experiment results of mass flow rates conclude that the RNG k-ε model is the most suitable model. The pressure distribution in the convergent-divergent nozzle was revealed by computational fluid dynamic （CFD） simulating on the turbulent flow field under different pressure conditions. The growing conditions of cavitation bubbles were shown; meanwhile, the phenomena in the experiment could be explained. The differential pres- sure between the upstream and downstream in nozzle throat section can improve the cavitating effect of cavitation water jet.

Simultaneous determination of 15 pesticide residues in Chinese cabbage and cucumber by liquid chromatography-tandem mass spectrometry utilizing online turbulent flow chromatography

In this experiment,a liquid chromatography tandem mass spectrometry method was built to determine 15 pesticide residues in Chinese cabbage and cucumber samples based on online turbulent flow chromatography purification.After modified quick,easy,cheap,effective,rugged,and safe(QuEChERS)extraction,extracts were directly injected to the TLX(TurboFlow Liquid Xcalibur)system and brought to TurboFlow™columns for on-line purification and then transferred to analytical column for further separation and analysis.TurboFlow™columns types,transfer flow rate,and transfer time were optimized.Limits of detection and limits of quantification of the method obtained for 15 pesticide residues were ranged between 0.2–1.0μg/kg and 0.5–2.0μg/kg in Chinese cabbage and cucumber samples.Recoveries of pesticide residues were in range of 75.3%–103.7%.Matrix effects for 15 pesticides were in range of 5.6%–106.6%.The developed method has been successfully used for the determination of 15 pesticide residues in real samples.

Investigation of scale efect for the computation of turbulent flow around a circular cylinder

In order to investigate the scale effect of turbulent flow around a circular cylinder, two similarity numbers （criteria） based on turbulent kinetic and dissipation rates associ- ated with the fluctuation characteristics of turbulence wake are deduced by analyzing the Reynolds averaged NavierStokes equations （RANS）. The RNG k-s models and finite volume method are used to solve the governing equations and the second-order implicit time and upwind space discretization algorithms are used to discrete the governing equations. A numerical computation of flow parameters around a two-dimensional circular cylinder with Reynolds numbers ranging from 102 to l07 is accomplished and the result indicates that the fluctuation of turbulence flow along the center line in the wake of circular cylinder can never be changed with increasing Reynolds numbers when Re ≥ 3 × 10^6. This conclusion is useful for controlling the scale of numerical calculations and for applying model test data to engineering practice.

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.