The Effect of Reynolds Number on Jet in Asymmetric Co-Flows : A CFD Study

In rotary kilns in grate-kiln systems for iron ore pelletizing, a long and stable jet flame is needed to ensure a high quality of the pellets. &e primary jet issuing from the nozzle interacts with two asymmetric co-flows creating a very complex flow. In order to better understand and eventually model this flow with quality and trust, simplified cases need to be studied. In this work, a simplified and virtual model is built based on a down-scaled kiln model established in a previous experimental work.&e aim is to numerically study the jet development as a function of position and Reynolds number (Re).&e numerical simulations are carried out with the standard k-ε model, and quite accurate velocity profiles are obtained while the centerline decays and spreading of the passive scalars are over predicted. &e model is capable of predicting a Re dependency of the jet development. With increasing Re, the jet is longer while it generally decays and spreads faster resulting from the stronger shear between the jet and co-flows and the stronger entrainment from the recirculation zone. &is recirculation found in the simulations restrain the momentum spreading in the spanwise direction, leading to a slower velocity spreading with higher Re. For further validation and understanding, more measurements in the shear layer and simulations with more advanced turbulence models are necessary.


Introduction
A grate kiln can be a part of an iron ore pellet sintering process, and it consists of a travelling grate, a rotary kiln, and a cooler.e rotary kiln can be about 35 m long and 6 m in diameter, and it carries iron ore pellets from the travelling grate to the cooler as the pellets are sintered at temperatures over 1000 °C.At the inlet of the rotary kiln, a burner is placed providing heat to the pellets.e flow from this burner can be modeled as a primary jet flow.Along with the primary jet, there are also two asymmetric co-flows from secondary air channels emanating from the cooler [1].
e two air channels provide excessive air for both the combustion and other important and earlier stages of the pelletizing process.As a result of the high flow rate and irregular shape of the air channels, the kiln performance is highly connected to the kiln aerodynamics, especially the development of the primary jet.e pelletizing process is complex since it is large scale and involves high temperatures and large amounts of air flowing in intricate geometries.Hence, it was realized early that relevant simplified cases should be studied to better understand the process and finally be able to model the whole process [2][3][4][5][6][7][8][9][10].One of these cases have been studied both experimentally [8] and numerically [9] with a cold flow as often done for combustors, especially for diffusion flames where the combustion reaction is diffusion and mixing limited.One of the well-recognized reasons is that the cold flow structure is not significantly different from that with combustion reactions [11].e early measurement by Glass and Bilger [12] indicated that a turbulent jet flame (or as they called "a heated jet-wake") and an isothermal turbulent round jet (isothermal jet-wake) have very similar turbulence quantities, such as the turbulence intensities and shear stress.Cold models provide an easy way to investigate the influence of different variables and facilitate measurements for validation [13]. is and the fact that the results can be compared to experiments, and previous simulations validate further studies on the simplified case in [2][3][4][5][6] now with focus on the jet development as a function of position and Reynolds number (Re).
Jet flow has been widely investigated from the beginning of the 1900s for its numerous applications.Great efforts have been put on the measurements of turbulent axisymmetric jets issuing into a quiescent environment in order to reveal the flow pattern (e.g.[14][15][16][17][18]).Due to the extensive applications in engineering problems, jets ejecting into moving co-flows are of great interest to researchers as well.e interactions between jets and the surrounding flows, which were named as "jet-wake" by Antonia and Bilger [19], make the flow field more complicated.e well-established selfpreserving property of turbulent axisymmetric jets moving into a fluid without any initial motion is no longer valid for jets in co-flows [17,20].Research has focused on jets in axisymmetric co-flows, and the most studied cases are coaxial jets.Less is done on jets with asymmetric co-flows, despite the fact that they are also widely applied in industrial combustors.
It is nowadays well known that the initial conditions of a jet such as nozzle geometry and Re influence the development of the jet downstream.Since about four decades ago, the conclusion of a universal self-similarity behavior of all jets started to be challenged [21].Boersma et al. [22] performed DNS simulations of jets with low Re and different initial velocity profiles.
eir conclusion agrees with the earlier research by George [21] that there is no universal selfsimilarity among jets with different initial conditions.Mi et al. [23] reviewed early experimental research on turbulent round jets issued from different nozzles.ey found that jets issuing from a smooth contraction nozzle have a larger spreading rate than jets issuing from a long straight round pipe.Moreover, the influence of the initial Re on the jet development in the self-similar far field depends on the type of the nozzle as well.e centerline decay of a passive scalar will increase with higher Re if the jet issues from a smooth contraction pipe, while a jet issuing from a long straight pipe is more independent of the initial Re.
By reviewing some of the experimental work regarding turbulent round jets with smoothly contracting jet nozzles with various Re ranges [18,[24][25][26][27][28][29], it can be concluded that there is no universal Re dependency of the jet mean centerline decay either.For jets with low Re, as stated in Todde et al. [24], the mean decay rate 1/B and the virtual origin have different Re dependencies.e decay rate 1/B decreases with Re until about Re � 1600, then remains constant until Re � 4000, and eventually increases again above 4000.In other high Re cases, generally the decay term B increases with Re; that is, the decay rate 1/B decreases, while above the threshold value (Re > 20000 as summarized in Zarruk and Cowen [29]), it is not significantly affected by Re.
e flames in the rotary kilns are confined diffusion jet flames with co-flowing air [5].Co-flows, as well as the jet initial conditions, can significantly influence the jet development.As stated by Yule and Damou [30], the extraturbulence level brought by the co-flows will influence the spreading rate of the confined jet.Moreover, the co-flows in this work are neither annular nor coaxial to the primary jet.us, the high turbulence brought by the irregular co-flows in the rotary kiln may cause an unstable flow field, for example, extra-recirculation of reactants and products.e unstable flow field could lead to problems like combustion instability and thermal fatigue of kiln walls.A long and stable jet flame is needed for an even temperature profile within the pellet bed at the bottom of the kiln.To increase the combustion efficiency and lower the emission of pollutants such as NO x , better mixing between the fuel (coal particles in the real kiln) and air is demanded.To reach these targets, it is necessary to know how to obtain a longer jet and better mixing in the kiln.e first option, based on the literature review of the previous research, is to change the jet initial condition by changing the geometry and see how the jet reacts to the changes.e purpose of this work is to study the effect of initial Re, by altering the nozzle diameter (D j ), on the far field development of the jet in a simplified cold model.e existence of the irregular-shaped secondary air channels brings an extra challenge compared to axisymmetric cases.Moreover, in the rotary kiln in this work, the air channels are separated by a surface called the "back plate."It covers 1/3 of the kiln diameter, and the jet exit is located at the center.e co-flowing air streams flow over the back plate, which acts as a bluff body, and merge together downstream.Flow over bluff bodies are connected to vortex shedding.Previous experimental work by Larsson et al. [8] however showed that the vortex shedding behind the back plate was hard to detect when the jet was introduced.e results clearly show that the shedding process is present and that it breaks up the jet on a large scale and negatively affects the small scale mixing which is promoted by a stable shear layer.In the previous numerical study presented in ASME 2016 International Mechanical Engineering Congress and Exposition (IMECE 2016) by Teng et al. [9], different turbulence models were compared in simulating the interaction between the jet and the asymmetric coflows.It was found that the standard k-ε model performs the best in predicting the mean centerline decay and radial spreading of velocity and scalar concentration.Similar conclusions can be found in early work by Smith et al. [31] where the results of mean and RMS centerline decay and spreading of the jet from k-ε models have good agreement with experimental data.In addition, the research work presented here on rotary kilns is highly linked to the industrial design, for which the robustness and computational economy are significant.Based on these facts, the standard k-ε model was employed in this work to further investigate the effect of initial Re on the mean flow field in the down-scaled kiln model.

Geometry and Mesh.
In previous research regarding kiln aerodynamics, a down-scaled, simplified rotary kiln model was established by Larsson et al. [1,6,8].e rotary kiln was down scaled based on the Craya-Curtet parameter C t � 0.86 as in the study by Larsson et al. [8], which represented a long 2 International Journal of Chemical Engineering and lazy flame.In this work, the geometry as presented in Figure 1 was built based on this down-scaled kiln model within a Cartesian coordinate system.e x-axis represents the streamwise direction.e y-and z-axis represent the spanwise directions.e kiln is simplified as a stationary cylinder with a diameter D kiln � 0.1 m, and the jet inlet is located in the center of the back plate with a diameter D j � 1.5 mm. e height of the back plate is 1/3 of the kiln diameter, and it separates the air channel into two parts.e upper secondary air channel has a 15 °inclining angle, and the lower one is parallel to the kiln.e kiln length is 1.5 m.As stated by Mi et al. [23] and Grandchamp et al. [32], the two commonly used nozzles, contraction nozzle, and straight round pipe, corresponding to a top-hat velocity profile and a parabolic velocity profile similar to the fully developed pipe flow, respectively.In this work, a top-hat velocity profile was approximated by a uniform mass flow rate at the jet inlet, since measurements of the initial velocity profile is lacking.To study the effect of initial Re, four cases were simulated with D j � 0.0012, 0.0015, 0.0018, and 0.0021 m corresponding to Re � 37240, 29792, 24827, and 21280.
With the commercial meshing software ANSYS 16.0 ICEMCFD, a hexagonal mesh was built and multiple o-grids were applied to refine the part near the jet inlet (Figure 1).Extra attention was paid to the grid expansion perpendicular to the jet centerline to avoid too large aspect ratios and to ensure sufficient mesh quality.To enable a grid independence study, three grids with different refinement levels were prepared.e details of this study are presented in Section 2.5.

CFD Setup.
e simulations were performed with the commercial code ANSYS CFX 16.0.Unsteady Reynoldsaveraged Navier-Stokes (uRANS) equations were used since steady state simulations oscillated regarding both the residuals and the monitored values.e time step was set to 1 ms based on the dimensionless time step ∆tU b /D kiln � 0.005 used in Sofia Larsson et al. [6].e simulation time was in its turn set to 5.5 s to ensure that the fluid would cover the entire kiln length, based on the bulk velocity.Transient statistics were taken from 5.5 to 7.5 s to obtain the time-averaged data.
e boundary conditions were the same as in Larsson et al. [8] and Teng et al. [9] with an even mass flow rate of 1.08 kg/s at the two secondary channel inlets.e mass flow rate at the jet nozzle inlet was about 0.031 kg/s.At the kiln outlet, an average static pressure with a relative pressure of 0.0 Pa was applied.At the jet nozzle inlet, a passive scalar (Φ) was initialized as 1.0 while at the inlets of the two co-flows, it was 0.0.e passive scalar helps to visualize the jet development and represents the mixing of the chemical species.
e advection scheme used was high resolution, which includes a blend factor that blends between first-and second-order accurate upwind schemes, ensuring a robust solution also in regions where the flow gradients change rapidly.A second-order backward Euler scheme was used to discretize the transient terms.e convergence criterion for the root-mean-square (RMS) residuals was 10 −6 for tight convergence, and double precision was used.Each transient run was initialized with the corresponding steady state result.

Governing Equations.
e continuity equation and the Navier-Stokes equations describing the mass and momentum transfer are where x i is the spatial coordinate, p is the pressure, t is time, and ] is the kinematic viscosity.e left-hand side in ( 2) is the material derivative of the velocity vector, u i , of the flow field can be expanded as e Reynolds decomposition of the turbulent velocity field gives the mean-momentum or Reynolds equations based on the Navier-Stokes equations as where the brackets denote time-averaged values.Equation (4) can only be solved for very simple cases due to the 6 components of the Reynolds stress tensor u i ′ u j ′ .ese stresses consist of an isotropic and an anisotropic part according to where k(x, t) is the turbulent kinetic energy, δ ij is the Kronecker delta, and a ij is the anisotropic part which, according to the specific assumption in the turbulentviscosity hypothesis, can be determined from where ] T is the turbulent viscosity.
For visualization of the jet and to describe the concentration field, a passive scalar is introduced in the numerical model.e Reynolds-averaged scalar transport equation is where Φ is the quantity of the passive scalar.It could for example be the concentration of a species.Φ ′ is the fluctuating scalar field, and D Φ is the molecular diffusivity, which is not considered in this nonreactive case since the molecular diffusion has quite little effect on the scalar transport compared to the convection.] t /Sc t is the turbulent diffusivity, in which ] t is the turbulent eddy viscosity, and Sc t is the turbulent Schmidt number.Sc t describes the ratio between the turbulent diffusivity and the molecular diffusivity of mass.In this work Sc t is set as a constant 0.7 based on research work regarding axisymmetric free jet flows by Tominaga and Stathopoulos [33] and Mejía et al. [34].S Φ represents the Reynolds-averaged source term, and in this nonreactive case, the source term is zero.

Turbulence Models.
e standard k-ε model is widely used in industrial problems for its robustness.In previous work by Granström et al. [10] and Teng et al. [9], the standard k-ε model was proved to predict the mean flow field in the kiln in a proper way.Since the detailed structures of the jet, especially the small scale flow motions, is not taken into account in this work, the standard k-ε model, again, was used to close the governing equations.e k-ε model solves two equations for turbulent kinetic energy k and the energy dissipation rate ε: where E ij represents the component of rate of deformation.e turbulent eddy viscosity μ t is specified as e standard k-ε model is typically referred to that developed by Launder and Sharma [35].e model constants used are

Grid Study.
To perform the grid independence study, three grids were prepared.e grid refinement factor calculated based on the procedure in Celik [36] was about 1.3.As mentioned in Section 2.1, the mesh was built by a structured mesh with multiple o-grids in order to refine the grid near the jet inlet.e refinement was done in a structured way to ensure a systematic discretization error estimation.
Two variables were chosen for the error estimation, the time-averaged scalar concentration at 60D j on the geometrical centerline, and the radial spreading rate S calculated based on the jet half width r 1/2 (x) from 0 to 100 D j .
e jet half width r 1/2 (x) is defined as the spanwise location where the mean velocity or scalar concentration equals to half of the jet centerline value.Ideally, a round jet develops symmetrically, and r 1/2 (x) should be isotropic at each YOZ plane.For a self-preserving jet, r 1/2 (x) has a linear relationship with the streamwise location x as follows [18]: e term S is the spreading rate for momentum (velocity) or a passive scalar in the spanwise direction.To distinguish between them, S denotes the spreading rate of momentum, and S Φ denotes the spreading rate of a passive scalar.At the right-hand side, x is the streamwise location of the observation point, and x 0 is the location of the virtual origin.For ideal cases, S and S Φ are treated as constants, and they can be calculated from (11) based on r 1/2 .e spreading rates represent how fast the jet spreads in the spanwise direction.e larger spreading rates a jet has, the faster it spreads, and as a result, the better it mixes with the surrounding flow.erefore, spreading rates are often used to describe the development of a turbulent round jet in the spanwise direction.
It is worth to note that, from previous simulations by Teng et al. [9], it was found that the jet develops asymmetrically because of the geometry.As a result, the jet half width r 1/2 (x) is no longer isotropic.e visualization of the velocity fields where the streamwise velocity U ≥ U c /2 from the steady state simulation of the case with D j � 0.0015 m clearly shows an asymmetric jet development downstream (Figure 2).U c is the jet centerline velocity at each streamwise location.As the jet develops, the shapes of the velocity contours are nearly elliptical, and they are defined as the jet 4 International Journal of Chemical Engineering half width regions in this work.Moreover, they are nonaxisymmetric as depicted in Figure 2. erefore, in order to quantify the jet spreading in the spanwise direction, the jet half widths in this work are estimated by the areas of these jet half regions taken every 2 D j downstream in the X-direction.
Richardson extrapolation was performed based on the procedure described by Celik [36] using the results from steady state simulations with three grids of different densities.e grid information and the Richardson extrapolation results show an oscillatory convergence with these three girds (Table 1).e errors between the extrapolated value and the simulated values, e ext,1 and e ext,2 , respectively, are quite small.It is not necessary to always use the finest mesh as long as the errors between the results with the second grid and the extrapolated value are small enough.In this work, the second grid was therefore chosen for further investigation to save computational time.

Results and Discussion
To investigate the effect of initial Re on the jet development, four simulations with different nozzle diameters D j � 0.0012, 0.0015, 0.0018, and 0.0021 m were run with grid 2 based on the grid study results.With the same initial mass flow rate 0.031 kg/s as in the PIV measurement by Larsson et al. [8], the jet initial velocities are U j � 27.70, 17.73, 12.31, and 9.04 m/s, respectively, for the four cases.e corresponding initial Re are 37240, 29792, 24827, and 21280.Among these four cases, the one with D j � 0.0015 m (Re � 29792) is seen as the standard case since it is directly built based on the work in [8].
e instantaneous streamwise velocity fields (Figure 3) at the last time step show the major difference brought by altering the nozzle diameter D j .In order to visualize the jet in a better way, the streamwise velocity is normalized as U/U b .U b is the bulk velocity in the whole kiln calculated by the total flow rate and the kiln cross-sectional area.e high-velocity region where U/U b > 5 can represent the jet length, and it increases with Re, which is beneficial to the sintering of pellets.On the contrary, the length of the region where 1 < U/U b < 5 significantly decreases with Re, indicating that the velocity decays faster in higher Re cases.
is region can be called the fast-decay region.It can be seen that a low-velocity region where U/U b < 1 forms at the upper kiln wall after the upper co-flow starts to interact with the primary jet in all cases.is low-velocity region orients the jet by pushing it downwards, restraining the jet from spreading upwards.Inside this region, recirculation may happen.e recirculation is marked by black contour lines in Figure 3 to visualize the location and area of the back flow, which is the region where the flow moves in the negative X-direction.
From the turbulent kinetic energy profile taken at x � 0.09 m (corresponding to 60 D j in the standard case), the existence and location of the recirculation zone disclosed in Figure 3 can be identified as well.In Figure 4, the turbulence kinetic energy is first scaled by the square of centerline velocity U c at x � 0.09 m (Figure 4(a)) and then expressed in the logarithm scale (Figure 4(b)) to highlight the region where the weak co-flows affect the flow field.e x-axis in Figure 4 is the normalized turbulence kinetic energy.e yaxis is the spanwise distance y from the jet centerline in Ydirection.e profiles of the kinetic energy depict that there is a large scale flow motion outside the jet, though it is weaker quantitatively compared to the jet itself.It corresponds to the upper recirculation zone seen in Figure 3.With increasing Re (decreasing D j ), the boundary of the recirculation zone moves inwards to the jet centerline in the Y-direction and the turbulence kinetic energy increases, indicating that the recirculation zone is enlarged and enhanced.
In combustion reactors, there are two types of recirculation zones: the internal recirculation zone (IRZ) and the external recirculation zone (ERZ).ey are of great significance for diffusion flames since they directly influence the mixing between the fuel and oxidizer.e IRZ is located inside the jet and has been proven to be able to improve the internal mixing and flame stability according to Driscoll et al. [37] and Vanoverberghe et al. [38] for instance.e ERZ is located outside the jet, and it can be generated by shear and the entrainment of ambient flow by the jet in a confined environment.Rabacal et al. [39]   International Journal of Chemical Engineering concluded that, based on large eddy simulation (LES) results of a pulverized coal jet flame, the reactant particles recirculate inside the ERZ until being fully burned.As a consequence, the ERZ may increase the residence time and the burnout rate of the fuel particles and then further affects the flame stability.In this work, the four cases are nonreactive and singlephase, and hence, the passive scalar field is used to show the effect of the ERZ formed near the upper kiln wall on mixing.A passive scalar is a diffusive contaminant without dynamic effect on the flow motion.e transport of a passive scalar is governed by Equation ( 7), in which the diffusion term consists of turbulent and molecular diffusion.
e distribution of the passive scalar can represent the mixing and entrainment of species and can be useful for the understanding of the mixing status in a fluid flow.
It can be seen from the scalar concentration field (Figure 5) that, at the same location as the ERZ shown in Figure 3, the passive scalar is entrained up towards the upper kiln wall.It indicates that the appearance of the ERZ promotes the spreading of the passive scalar in the spanwise direction in the upper part of the kiln.In the reactive case, the ERZ can improve the mixing between the reactants, assist in the flame-anchoring, and further benefit  Besides the spreading rate mentioned earlier, another important characterization of the mean flow field of a turbulent round jet is the centerline decay.e decay of the mean velocity of such a jet at the centerline U c , with the streamwise distance from the nozzle x, can be expressed with a linear relationship as in [18]: Similarly, it has been suggested that the decay of the centerline concentration Φ c of a turbulent round jet can be written as in [23]: e terms 1/B and 1/B Φ , the decay rates of momentum (velocity), and the passive scalar, for ideal cases, are treated as constants.U j is the jet initial velocity, and Φ j is the passive scalar concentration at the jet inlet, which was defined as 1 in the mass fraction form in this work.
e decay rates of the mean streamwise velocity 1/B and the scalar concentration 1/B Φ are calculated by the curve fitting the time-averaged data from simulation results based on Equations ( 12) and (13).For validation, data from 30 to 50 D j were used for calculation in order to be consistent with Larsson et al. [8].From the profiles of U j /U c and Φ j /Φ c shown in Figure 6, it can be seen that, for x < 6 D j , there is no decay for any Re. is is the so-called potential core of the jet with a virtually uniform mean velocity.e decay rate of the mean streamwise velocity 1/B increases with Re (Figure 6(a)), while 1/B Φ of the mean scalar is not significantly affected (Figure 6(b)). is independency of Re of the mean scalar decay agrees with the findings in the coupled PIV-LIF measurements by Zarruk and Cowen [29].It is also clear that 1/B u is dependent on x for all Re investigated, and hence, 1/B u � f 1 (Re, x).
e spreading rates of the streamwise velocity S and the passive scalar concentration S Φ in the spanwise direction are calculated based on Equation (11).As described in Section 2.5, the jet half widths r 1/2 are estimated from the areas of jet half width regions taken from the time-averaged data.From Figure 7, it can be seen that the spreading rate of velocity S slightly decreases with Re (Figure 7(a)), but it almost remains constant.e spreading rate of the scalar S Φ has an obvious different trend in that it increases with Re (Figure 7(b)), meaning that S Φ � f 2 (Re, x).In addition, a slower increase of S Φ is seen after about 25 D j in all cases, indicating that the scalar spreads slower.is agrees with Larsson et al. [23], who concluded that there is less mixing after the wake region extending to about x/D j � 25.With higher Re (smaller D j ), this effect tends to be eliminated.
All the terms calculated from the centerline decay and spreading are summarized in Table 2 together with the experimental data from Larsson et al. [8].A comparison of the simulation results to the measurement shows that the standard k-ε model overpredicts the centerline decay term B by about 45% and the spreading rate of scalar concentration S Φ by about 20%.e prediction of the spreading rate of the velocity S is relatively accurate with an error of about 3%.
e k-ε models can be improved for better prediction of the centerline decay and spreading of turbulent round jets by modifying the model constants based on individual decay and spreading data of each jet [38].International Journal of Chemical Engineering Zarruk and Cowen [29] also concluded that the decay term B of the mean velocity increases with initial Re and then tends to be stable at very high Re (>20,000) by reviewing the studies of turbulent round jets without co-flow.An opposite trend is seen in this work where B decreases with Re. ere are several possible explanations to the discrepancy in Re dependency of mean centerline decay and spreading between this work and those studies without co-flow.Firstly, in this   When Re increases, the velocity difference between the jet and the coflows increase as well, resulting in stronger shear between them.Entrainment of the surrounding flows into the primary jet is promoted by the shear layers at the interfaces between the jet and co-flows formed by the large velocity gradient there.Strong shear can be identified from the radial profile of the shear stress uv taken at x � 0.015, 0.045, and 0.09 m (Figure 8).e x-axis, stress uv, is scaled by the square of the centerline velocity U c at each streamwise location.e y-axis, the spanwise distance from the centerline y, is scaled by the jet half widths calculated previously.It can be seen that, with increasing Re, the shear layers are stronger and located more outwards in the spanwise direction.In words, the jet is wider with higher Re.In addition, it is found before that the ERZ formed near the upper kiln wall suppresses the jet downwards and this zone is enlarged and reinforced at higher Re.At the same time, it entrains the passive scalar upwards.In other words, the ERZ restrains the momentum spreading while it prompts the scalar spreading in the spanwise direction.
In earlier work by Larsson et al. [5], it was found that the back plate between the two co-flows causes vortex shedding.Furthermore, in the measurements by Larsson et al. [8], the vortex shedding behind the back plate was hard to detect when the primary jet was introduced.e results clearly shows that the shedding process is present and that it breaks up the jet on a large scale and negatively affects the small scale mixing which is promoted by a stable shear layer.However, the effect from such underlying flow structures cannot be modeled by k-ε models.is might be another possible reason for the opposite trend of the mean flow field dependence on Re.Two-equation models are not capable of predicting anisotropy and different turbulence International Journal of Chemical Engineering length scales.As stated in Smith et al. [40], even though k-ε models perform better than other two-equation models in predicting the spreading rate of round jets, the underlying flow structures are not taken into account.It results in failure of predicting the correct trends of the centerline decay rate of jets with different nozzle geometries. is problem might be solved by finding suitable model constants for jets with different initial conditions where the constants are functions of the centerline decay rates and the jet half widths [41].Despite the fact that the spreading issue can be taken care of, k-ε models still cannot capture the effect from vortex shedding.
is implies that more advanced turbulence models are necessary, as well as more measurements, for a more thorough understanding of how the jet initial conditions affect the flow field.
is work including the published prestudy [9] are well summarized in thesis [42].

Conclusions
Based on the results from CFD simulations with the standard k-ε model and comparisons to other research work, some conclusions can be drawn as below.
e standard k-ε model gave the closest prediction to the measurements in previous work.
e prediction of the spreading rate S of the streamwise velocity is quite accurate with an error of about 3%.ere is an overprediction of the centerline decay term B of the velocity by about 45% and of the spreading rate S Φ of the passive scalar by about 20% compared to the results in previous measurements.e model is capable of predicting a Re dependency of the jet development as the jet with higher initial Re is longer, which is beneficial to the sintering of pellets, but the centerline velocity decays faster downstream.e spanwise spreading rate S Φ of the passive scalar increases with Re, but the spreading rate S of the velocity almost remains constant when Re increases.One of the reasons behind this Re dependency could be the enhanced shear layer between the jet and co-flows resulting from higher velocity difference when Re increases.is stronger shear also improves the scalar mixing in the flow field.
An external recirculation zone (ERZ) forms at the upper kiln wall after the inclined upper co-flow and starts to interact with the primary jet as a result of the shear between them.It orients the jet slightly downwards, restraining the jet from spreading towards the upper kiln wall, and preventing possible thermal fatigue there.e ERZ entrains the passive scalar upwards, promoting the spreading of scalar in the spanwise direction.Moreover, the ERZ grows when Re increases.In a reactive case, the ERZ can improve the mixing between the reactants, assists in the flame anchoring, and hence further benefit the flame stability.Since a long jet stably located in the center of the kiln is expected in this type of kiln for iron ore pellets sintering, the existence of the ERZ is therefore significantly advantageous.
e k-ε model can be improved for better accuracy but no matter how accurate it can predict the mean properties of the jet development, two equation models are still not capable of capturing the underlying flow structures such as vortex shedding and its effect on small-scale mixing.For better validation and understanding of the flow field in this virtual kiln model, more measurements and advanced turbulence models or methods will be considered in future work.

Data Availability
e data supporting the conclusions in this work are included in this manuscript.And the other datasets generated and/or analysed during the current study are available from the corresponding author on reasonable request.

Figure 1 :
Figure 1: Geometry (part) of the down-scaled, simplified kiln model in Cartesian coordinate system and the hexagonal mesh (fine grid) with refinement near the jet inlet with multiple o-grids.

Figure 8 :
Figure 8: e profile of shear stress uv/U 2 c ∼ y/r 1/2 at (a) x � 0.015 m, (b) x � 0.045 m, and (c) x � 0.09 m (corresponding to 10, 30, and 45 D j in the standard case) in the Y-direction.

Table 1 :
Grid information and Richardson extrapolation resultswith three grids and steady state simulations.

Table 2 :
e calculated centerline decay and radial spreading terms.