Numerical Simulation of the Effect of Air Distribution on Turbulent Flow and Combustion in a Tubular Heating Furnace

A three-dimension full-size numerical simulation of the effect of air distribution on turbulent flow and combustion in a tubular heating furnace was carried out. A standard k –ε turbulent model, a simplified PDF combustion model and a discrete ordinate transfer radiation model were used. The hybrid grid combining a structured and a non-structured grid was generated without any simplification of the complicated geometric configuration around the burner. It was found that the multistage combustion could reduce and control the peak value of temperature. At the same time, it was concluded that the amount of primary air had little effect on the global distribution of velocity and temperature in the furnace, but a great effect on that around the burner. It is recommended that 45% - 65% of the total amount of air be taken in in primary air inlets in the furnace. All the results are important to optimize the combustion progress.

Comparison Between LES and RANS Modeling of Turbulent Swirling Flows and Swirling Diffusion Combustion

Turbulent swirling flows and methane-air swirling diffusion combustion are simulated by both large-eddy simulation (LES) using a Smagorinsky-Lilly subgrid-scale (SGS) turbulence model, a second-order moment (SOM) subgrid-scale combustion model and an eddy break up (EBU) combustion model and Reynolds-averaged NavierStokes (RANS) modeling using the Reynolds stress equation model and a second-order moment (SOM) combustion model. For swirling flows, the LES statistical results give better agreement with the experimental results than the RANS modeling, indicating that the adopted subgrid-scale turbulence model is suitable for swirling flows. For swirling combustion, both the proposed SOM SGS combustion model and the RANS-SOM model give the results in good agreement with the experimental results, but the LES-EBU modeling results are not in agreement with the experimental results.

Large-eddy structures of turbulent swirling flows and methane-air swirling diffusion combustion

Turbulent swirling flows and methane-air swirling diffusion combustion are studied by large-eddy simulation （LES） using a Smagorinsky-Lilly subgrid scale turbulence model and a second-order moment （SOM） SGS combustion model, and also by RANS modeling using the Reynolds Stress equation model with the IPCM＋wall and IPCM pressure-strain models and SOM combustion model. The LES statistical results for swirling flows give good agreement with the experimental results, indicating that the adopted subgrid-scale turbulence model is suitable for swirling flows. The LES instantaneous results show the complex vortex shedding pattern in swirling flows. The initially formed large vortex structures soon break up in swirling flows. The LES statistical results of combustion modeling are near the experimental results and are as good as the RANS-SOM modeling results. The LES results show that the size and range of large vortex structures in swirling combustion are different from those of isothermal swirling flows, and the chemical reaction is intensified by the large-eddy vortex structures.

Validation of the RANS-SOM Combustion Model Using Direct Numerical Simulation of Incompressible Turbulent Reacting Flows

The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation （DNS） of incompressible turbulent reacting channel flows. The instantaneous DNS results show the near-wall strip structures of concentration and temperature fluctuations. The DNS statistical results give the budget of the terms in the correlation equations, showing that the production and dissipation terms are most important. The DNS statistical data are used to validate the closure model in RANS second-order moment （SOM） combustion model. It is found that the simulated diffusion and production terms are in agreement with the DNS data in most flow regions, except in the near-wall region, where the near-wall modification should be made, and the closure model for the dissipation term needs further improvement. The algebraic second-order moment （ASOM） combustion model is well validated by DNS.

Simulation of gas-solid flow characteristics of the circulating fluidized bed boiler under pure-oxygen combustion conditions

Under the pressure of carbon neutrality,many carbon capture,utilization and storage technologies have witnessed rapid development in the recent years,including oxy-fuel combustion(OFC)technology.However,the conventional OFC technology usually depends on the flue gas recirculation system,which faces significant investment,high energy consumption,and potential low-temperature corrosion problem.Considering these deficiencies,the direct utilization of pure oxygen to achieve particle fluidization and fuel combustion may reduce the overall energy consumption and CO_(2)-capture costs.In this paper,the fundamental structure of a self-designed 130 t·h^(-1) pure-oxygen combustion circulating fluidized bed(CFB)boiler was provided,and the computational particle fluid dynamics method was used to analyze the gas-solid flow characteristics of this new-concept boiler under different working conditions.The results indicate that through the careful selection of design or operational parameters,such as average bed-material size and fluidization velocity,the pure-oxygen combustion CFB system can maintain the ideal fluidization state,namely significant internal and external particle circulation.Besides,the contraction section of the boiler leads to the particle backflow in the lower furnace,resulting in the particle suspension concentration near the wall region being higher than that in the center region.Conversely,the upper furnace still retains the classic core-annulus flow structure.In addition to increasing solid circulation rate by reducing the average bed-material size,altering primary gas ratio and bed inventory can also exert varying degrees of influence on the gas-solid flow characteristics of the pure-oxygen combustion CFB boiler.

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

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.

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.

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.

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.

Turbulent boundary layers and hydrodynamic flow analysis of nanofluids over a plate

A numerical analysis of the log-law behavior for the turbulent boundary layer of a wall-bounded flow is performed over a flat plate immersed in three nanofluids(Zn O-water,SiO_(2)-water,TiO_(2)-water).Numerical simulations using CFD code are employed to investigate the boundary layer and the hydrodynamic flow.To validate the current numerical model,measurement points from published works were used,and the compared results were in good compliance.Simulations were carried out for the velocity series of 0.04,0.4 and 4 m/s and nanoparticle concentrations0.1% and 5%.The influence of nanoparticles’ concentration on velocity,temperature profiles,wall shear stress,and turbulent intensity was investigated.The obtained results showed that the viscous sub-layer,the buffer layer,and the loglaw layer along the potential-flow layer could be analyzed based on their curving quality in the regions which have just a single wall distance.It was seen that the viscous sub-layer is the biggest area in comparison with other areas.Alternatively,the section where the temperature changes considerably correspond to the thermal boundary layer’s thickness goes a downward trend when the velocity decreases.The thermal boundary layer gets deep away from the leading edge.However,a rise in the volume fraction of nanoparticles indicated a minor impact on the shear stress developed in the wall.In all cases,the thickness of the boundary layer undergoes a downward trend as the velocity increases,whereas increasing the nanoparticle concentrations would enhance the thickness.More precisely,the log layer is closed with log law,and it is minimal between Y^(+)=50 and Y^(+)=95.The temperature for nanoparticle concentration φ=5%is higher than that for φ=0.1%,in boundary layers,for all studied nanofluids.However,it is established that the behavior is inverted from the value of Y^(+)=1 and the temperature for φ =0.1% is more important than the case of φ =5%.For turbulence intensity peak,this peak exists at Y^(+)=100 for v=4 m/s,Y^(+)=10 for v=0.4 m/s and Y^(+)=8 for v=0.04 m/s.

Numerical Investigations on the Impact of Turbulent Prandtl Number and Schmidt Number on Supersonic Combustion

The flow field inside the combustor of a scramjet is highly complicated and the related turbulent Prandtl and Schmidt numbers have a significant impact on the effective numerical prediction of such dynamics.As in many cases researchers set these parameters on the basis of purely empirical laws,assessing their impact(via parametric numerical simulations)is a subject of great importance.In the present work,in particular,two test cases with different characteristics are selected for further evaluation of the role played by these non-dimensional numbers:Burrows-Kurkov case and DLR case.The numerical results indicate that these parameters influence ignition location.Moreover,the temperature distribution is more sensitive to them than to H2O mass fraction and velocity distributions.

NUMERICAL STUDY OF AN OSCILLATORY TURBULENT FLOW OVER A FLAT PLATE

Oscillatory turbulent flow over a flat plate was studied by using large eddy simulation (LES) and Reynolds-average Navier-Stokes (RANS) methods. A dynamic subgrid-scale model was employed in LES and Saffman's turbulence model was used in RANS. The flow behaviors were discussed for the accelerating and decelerating phases during the oscillating cycle. The friction force on the wall and its phase shift from laminar to turbulent regime were also investigated for different Reynolds numbers. (Edited author abstract) 11 Refs.

Three-dimensional turbulent model of heat transfer and fluid flow in GTAW process

A two-equation K-ε turbulent fluid flow model is built to model the heat transfer and fluid flow in gas tungsten arc welding （GTAW） process of stainless steel S US310 and S US316. This model combines the buoyancy force, lorentz force and marangni force as the driving forces of thefluidflow in the weld pool. The material properties are functions of temperature in this model. The simulated results show that the molten metal flowing outward is mainly caused by the marangoni convection, which makes the weld pool become wider and shallower. The comparison of the weld pool shape of SUS310 and SUS316 shows that the slight differences of the value of thermal conductivity mainly attributes to the difference of the weld pool shape and the distinction of heat transport in laminar and turbulent model makes large diversity in the simulated results.

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.

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.

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.

TURBULENT SEPARATED REATTACHED FLOW IN A TWO-DIMENSIONAL CURVED-WALL DIFFUSER

A turbulent separation-rcattachment flow in a two-dimensional asymmetrical curved-wall diffuser is studied by a two-dimensional laser doppler velocimeter.The turbulent boundary layer separates on the lower curved wall under strong pressure gradient and then reattaches on a parallel channel.At the inlet of the diffuser,Reynolds number based on the diffuser height is 1.2×10~5 and the velocity is 25.2m/s.The re- sults of experiments are presented and analyzed in new defined streamline-aligned coordinates.The experiment shows that after Transitory Detachment Reynolds shear stress is negative in the near-wall backflow region. Their characteristics are approximately the same as in simple turbulent shear layers near the maximum Reynolds shear stress.A scale is formed using the maximum Reynolds shear stresses.It is found that a Reynolds shear stress similarity exists from separation to reattachment and the Schofield-Perry velocity law ex- ists in the forward shear flow.Both profiles are used in the experimental work that leads to the design of a new eddy-viscosity model.The length scale is taken from that developed by Schofield and Perry.The composite velocity scale is formed by the maximum Reynolds shear stress and the Schofield Perry velocity scale as well as the edge velocity of the boundary layer.The results of these experiments are presented in this paper