A two-dimensional (2D) jet flow and temperature field are simulated by using k- T. C. model and compared with other three nontemperature corrected models, which are standard k-, RNG k-, and SST k- model. Then based on the calculated results, the spectral infrared radiation characteristics within 4∼5 m of the 2D jet flow were calculated. By comparing the computed results of the velocity, temperature field, and infrared radiation with the experimental measurements, it shows that the k- T. C. model predicts mean flow mixing more rapidly and the turbulent kinetic energy dissipates earlier than with no temperature correction; the k- T. C. model could give a good prediction for the velocity and temperature distributions on the centerline of the 2D hot gas jet, but not on the locations off the centerline. The maximum computation error of the 2D hot jet infrared radiation is decreased from 86% to 26%, and the accuracy of the computation is greatly improved.

1. Introduction

The effect of infrared radiation of hot jets is of major importance in the target detection and recognition. Two-dimensional (2D) nozzles which can produce a strong secondary flow are usually applied to reduce the infrared signature of the hot jets [13]. Therefore the improvement of the computation accuracy of 2D hot jet’s infrared radiation is of great significance.

The computation of infrared radiation of hot jets relies on accurate simulation of temperature field. Currently, the most commonly used technique to compute the temperature field is to solve the Reynolds Averaged Navier-Stokes (RANS) equations [4] and scalar conservation equations. These equations are solved with turbulence models, which establish the relation between the Reynolds stress in RANS equations and the turbulent scalar flux in Reynolds averaged scalar equations with the velocity gradient and scalar (such as temperature and concentration) gradient of the flow. Therefore, the turbulence model adopted in the calculation will affect the accuracy of the hot jet temperature filed and in turn affect the calculation accuracy of the infrared radiation.

There are dozens of turbulence models, of them two-equation models, such as standard k-ε model [5], RNG k-ε model [6], and SST k-ω model [7], are currently widely used in engineering [8]. Because the parameters in the equations of these models are calibrated with simple shear flow at the condition of room temperature and low Mach number without considering the variation of temperature, they are often difficult to accurately simulate high temperature jet flow. For example, Dembowski and Georgiadis [9] numericaly simulated the velocity and temperature field of a nozzle with SST k-ω model and Chien k-ε model and compared with Seiner’s et al. [10] experiment data, the results show that accurate results can only be obtained when the temperature of jet is not high, for example, between 313 K and 508 K, and then the deviation will significantly increase while the temperature increases from 508 K to 1366 K.

In order to make the two-equation turbulence models still be applicable to the high temperature jets, Abdol-Hamid et al. [11] proposed a k-ε model based on temperature correction, namely, k-ε T. C. model, in which the eddy viscosity in the model is modified by the local total temperature gradient normalized by the local turbulence length scale. Results show that the model has a good adaptability to both subsonic and supersonic jets with a wide temperature range. Tam and Ganesan [12] proposed a modified k-ε turbulence model to include the effect of the presence of a large-density gradient. In addition, Massey et al. [13] also used k-ε T. C. to study the separation flow of chevron nozzle and got the similar results to Abdol-Hamid. Georgiadis et al. [14] evaluated the k-ε T. C. model (or PAB T. C.) for round jet flow; the results showed that the modified model improved predictions of mean axial velocities for the heated jet but did not improve prediction of the calculated turbulent kinetic energy fields. The study of Li et al. [15, 16] showed that those modified k-ε model provided improved infrared radiation predictions for the same round jet.

While there are considerable experimental and computational studies [1719] on cold jet flow of two-dimensional (rectangular) nozzle, there are comparatively few studies on hot jet flow of this type of nozzle, especially the infrared radiation. The main purpose of this paper is to study how the two-equation turbulence models influence the temperature and infrared radiation of a 2D jet flow. The turbulence models studied includes standard k-ε, RNG k-ε, SST k-ω, and k-ε T. C. models. By comparing the computed results of the velocity, temperature field, and infrared radiation with the experimental measurements, it shows that the calculation with k-ε T. C. model greatly improved the accuracy of the computation.

2. Computation Method of Hot Jet’s Infrared

2.1. Procedure of Hot Jet Infrared Radiation Computation

The essence of hot jet’s infrared computation is solving the radiation transfer equation, which depends on the temperature and species concentration distributions of the jet flow. The temperature and species concentration distributions were obtained by solving the momentum, energy, and species transfer equations by the method of CFD.

The procedure of hot jet infrared radiation computation is a complex process, which involves nozzle geometry modeling, CFD computation domain meshing, CFD computation, IR (infrared radiation) computation, correction of atmosphere transmittance, and results display. The relationship between CFD computation and infrared radiation computation is shown in Figure 1. As shown in the figure, the computation was divided into two modules, which are commercial CFD module and IR module, respectively. The temperature, pressure, and species concentration of hot jet outputted from CFD module were interpolated into the IR module. The results of CFD module were greatly affected by the turbulence model, and four different turbulence models will be discussed in the next section.

2.2. Four Two-Equation Turbulence Models

In the scheme of eddy viscosity models (EVM), the Reynolds stress and turbulent scalar flux are determined with where and are mean and fluctuation velocity in direction, respectively, and are mean and fluctuation scalar, respectively, which could be temperature or species concentration, k is turbulent kinetic energy, and is turbulent Prandtl number or turbulent Schmidt number. is eddy viscosity, which is determined with different turbulent models. The Prandtl number takes 0.85, and the Schmidt number takes 0.7 in this paper.

2.2.1. Standard k-ε Model

In standard k-ε model [5], the expression for is defined by where turbulent kinetic energy and its dissipation rate  ε  are determined by solving their transport equations as follows: where , , , , and .

2.2.2. RNG k-ε Model

The RNG k-ε model [6] was derived using a statistical technique called renormalization group theory. It is similar in form to the standard k-ε model, but it includes the following refinements: the RNG model has an additional term in its ε   equation that improves the accuracy for rapidly strained flows; the effect of swirl on turbulence is included in the RNG k-ε model in order to improve the accuracy for swirling flows; the renormalization group theory provides an analytical formula for turbulent Prandtl number, while the standard k-ε model only uses constant values. With these features, the RNG k-ε model becomes more accurate and reliable for wider class of flows than the standard k-ε model. In RNG model, , , and , .

2.2.3. SST k-ω Model

The Shear-Stress Transport (SST) k-ω model was developed by Menter [7] to effectively combine the robust and accurate formulation of the k-ω model in the near-wall region with the free stream independence of the k-ε model in the distanced field. The definition of eddy viscosity is where . The function of  ω is defined as with . S is strain rate magnitude; . The functions of are

The transport equations of the SST k-ω model are In these equations, represents the generation of turbulence kinetic energy due to mean velocity gradients. The term is defined as The more detailed description of this model could be found in [20].

2.2.4. k-ε T. C. Model

The k-ε T. C. model [12] has the same transport equations (as shown in (3)) as standard k-ε model, but the eddy viscosity is corrected by the local total temperature gradient as follows:

The correction coefficient, , is calculated with (8) as follows: where is stagnation temperature. The function is a correction factor to improve the applicability of the model to high-speed flow; is defined as where , , a is local acoustic speed, and is Heaviside step function.

2.3. The Method to Solve Radiation Transport Equation

Infrared radiation transmission in any direction of is governed by the radiation transport equation [20] (14). Consider where is spectral radiance, s is the local location, is spectral absorption coefficient, is blackbody emissive power, and is local temperature.

By integrating (14) along the ray, the total radiance of the hot jets along the direction of can be obtained. The ray tracing method is used to calculate the infrared radiation transmission of the hot jet. The illustration of the ray tracing method is shown in Figure 2.

The computational expression of is In this work, the narrow band model is used to calculate the absorption coefficient and the spectral step is 10 cm−1. The original data of can be found in NASA SP-3080 or HITEMP database.

The transmittance of the air between the target and the receiver at the band of 4.26 μm is necessary to consider, although the concentration of carbon dioxide in the air is very less. So the total radiation intensity in the field of view (see Figure 2) between wavelength and can be calculated by using where is the distance between the hot jets and the receiver. is the transmissivity of atmosphere between the hot jets and the receiver, and it was calculated by using software Modtran. In the calculation, the atmosphere is treated as a medium without scattering. The composition of the midlatitude summer atmosphere is used to compute atmosphere attenuation.

From the above equations, we can see that the infrared radiation or radiant intensity is mainly dependent on both the temperature and the absorption coefficient of the hot jets.

3. Experimental Measurement of Velocity, Temperature, and Infrared Radiation

3.1. Test Model and Measurement Arrangement

A 2D nozzle with aspect ratio of 6 is used in present study, as shown in Figure 3, the corss section of which is gradually transformed from circular to rectangular. The equivalent diameter of the nozzle outlet is defined as , where is area and is perimeter. The mixer divides the axisymmetric part of the nozzle into two channels, the main channel and bypass channel. The main gas stream has a total temperature of 845 K and mass flow rate of 1 kg/s. The bypass airflow has a total temperature of 342 K and mass flow rate of 1.2 kg/s. The Mach number of nozzle outlet is about 0.15, and the Reynolds number is about 7.1 × 105.

The measurement scenario of pressure, temperature, and infrared radiation of the gas jet is shown in Figures 4(a) and 4(b). The temperature distribution was measured by thermocouples rake, and the total pressure and the static pressure are measured by the pitot tubes rake. The measurement was made within the range from the nozzle exit to 6De downstream. Based on the experiments, the velocity of the gas flow can be calculated. The infrared radiation of the gas jets on the direction of 90° to the axis of the nozzle is measured using a FT-spectroradiometer which has a wavelength band 45 μm and located 34 meters away from the nozzle. The measurements were made in both wide and narrow sides of the 2D nozzle, as shown in Figure 4(b).

3.2. Measurement Results

The measurement results include centerline velocity, temperature, radial velocity, temperature, and spectral infrared radiation, which are given in Figures 10, 14, and 15. The detailed description is omitted here.

3.3. Measurement Uncertainty Analysis

The temperature measurement accuracy with thermocouples is influenced by shunting error, response time error, and radiation error. The type of the thermocouple was commercial K-type sheathed thermocouple. A multichannel temperature inspection instrument was used to acquire the temperature values. Considering the uncertainty of temperature inspection instrument, the total uncertainty of temperature measurement is about ±(1%T + 1.0°C).

The uncertainties of velocity measurement include the uncertainty of both temperature and pressure. The angles between the pitots and the jet axial component were less than 10°; the total uncertainty is about ±2%.

The infrared radiation measurement accuracy with FT-spectroradiometer is influenced by three types of errors. The first type of errors is calibration source errors, which include the error of blackbody temperature, about 2.5%, and the error of blackbody emissivity over the operational spectral range, about 2%. The second type of errors is calibration drift, which is about 2%. Finally, the third type of errors is the spectrometer intrinsic linearity, which is about 1%. These errors are independent of each other, and the total error is calculated by adding up three errors’ contribution in root-sum-square fashion, which is about 4%.

4. Verification of the CFD and IR Computation

The distributions of flow field, temperature field, and species concentration field are computed using the CFD software Fluent [21]. For CFD code, the most important numerical errors and uncertainties are due to use of iterative solution methods and specification of various input parameters such as spatial and time step sizes and other parameters [22]. Since the steady-state simulations are carried out, iteration number and grids size are the major contributions to the uncertainties [23].

4.1. Calculation Domain, Boundary Condition, and Grids Generation

Since the nozzle is symmetrical, in order to save computing time, only 1/4 of the whole region was simulated. The simulation domain and boundary types are shown in Figure 5. The boundary condition of main channel and bypass channel is set to mass flow inlet, and the boundary condition of simulation domain is set to pressure outlet. The boundary condition of symmetric plane is set to symmetry.

The nozzle geometry was complex, and in order to get structured grid, the O-grid method was applied to the meshing domain. The schematic diagram of meshing method is shown in Figure 6. Fan-shaped region is divided into several subblocks. The block has three directions, which are i, j, and .

Uncertainties due to grid size are estimated using systematically refined grids with constant refinement ration, . Grid studies were conducted using three grids. Grid dimensions for grid refinement studies are given in Table 1, and a comparison of the three grids at symmetry plane is shown in Figure 7. The total number of cells of grid 2 and 3 is about 2.7 and 2.9 times to that of gird 1 and 2, respectively. The product of does not equal the total number of cells and this is due to the way of meshing. In all cases, the simulation domain extended roughly 22De downstream of the nozzle exit and 4.2De from the centerline in the radial direction.

4.2. Numerical Solution of RANS CFD Code

The governing equations described above were solved numerically using commercial software, Fluent. Standard k-ε model, RNG k-ε model, and SST k-ω model can be directly selected in the software. The k-ε T. C. model was achieved by embedding a UDF function in the software. The density is calculated using the ideal gas model, and the molecule viscosity is calculated using the Sutherland model. An implicit solution based on the finite volume method (FVM) along with a quadrilateral nonuniform structured grid in the Cartesian coordinate system was used. The flow, turbulent kinetic energy, and specific dissipation rate equations were discretized using the second-order upwind scheme at the cell faces.

4.3. Numerical Uncertainties Analysis of CFD Computation

The convergence criterion of numerical computation is that the solution residuals and nozzle outlet temperature do not change significantly with the number of iterations. Iteration history for 2D nozzle on gird 3 is shown in Figure 8. As shown in the figure, the solution changes drop three to five orders of magnitude. The variation in dimensional temperature at nozzle outlet is about 0.15% (where is the solution on the finest-grid 3) over the last period of oscillation. Iterative uncertainty is estimated as half of the range of the maximum and minimum values over the last two periods of oscillation, so the iterative uncertainty is about 0.075%.

Grid convergence is assessed through core lengths solutions on three systematically refined grids with constant refinement ratio (see Table 1). The core length is defined for the in which (where is the temperature at the nozzle exit) along the axis of the jet. Simulated core lengths using k-ε T. C. and RNG k-ε model are shown in Table 2. The core lengths using k-ε T. C. model are much shorter than RNG k-ε model and display monotonic convergence with the increase in the number of gird cells. According to Stern and Wilson theory [23], the estimated grid uncertainty is less than 0.82% and 1.4% for k-ε T. C. model and RNG k-ε model, respectively. The core lengths using standard k-ε and SST k-ω modelare between k-ε T. C. and RNG k-ε model, and the uncertainties of these two models may be between 0.82% and 1.4%.

The difference of core length with different turbulence models exceeds 100%, which is obviously greater than the numerical uncertainty using the finest-grid 3. The following computation and analysis were all carried out on grid 3.

4.4. Numerical Uncertainties Analysis of IR Computation

An important factor that affects the IR computation is the size of ray discrete step (see Figure 2). The influence of ray discrete step on radiation intensity and consumed time is shown in Figure 9. The shorter the , the more stable the radiation intensity, and the more time spent. The numerical uncertainties of IR computation are about 0.5% when . In the following calculations, were taken as 0.05.

5. Results and Analysis

5.1. Velocity and Turbulent Kinetic Energy Distribution of the Jet

A comparison of centerline axial velocities profiles obtained from the four turbulence models is made with experimental data for the 2D jet in Figure 10. It shows that the RNG k-ε model gives the longest length of potential core region which is approximately 1.8 times longer than the measurement values. Beyond the potential core, however, the rate of jet decay is slower than that produced by the other models and the experimental data. The standard k-ε and SST k-ω   models get similar predicted length of the core region that is about 1.2 times longer and then downstream mixing rates that are also too slow relative to the experimental data. The k-ε T. C. model dramatically improves the accuracy of the predicted length of the velocity core region because the correction of eddy viscosity is only 0.2 times longer than the measured data. The k-ε T. C. model results have a slightly slow rate of jet decay compared to experimental data.

The computed centerline turbulent kinetic energy profiles and the turbulent kinetic energy contours of the two symmetric planes are compared in Figure 11 to Figure 13. The decline change which appears between is caused by reflow phenomenon after central cone. The RNG k-ε model results show the slowest propagation of turbulence to the jet centerline. The standard k-ε and SST k-ω models approaches enable a faster transport of turbulent kinetic energy to RNG k-ε model. The k-ε T. C. model gives the fastest transport of turbulent kinetic energy. This could be more obviously observed in the turbulent kinetic energy contours shown in Figure 12 and Figure 13. The peak turbulent kinetic energy region produced in the jet shear layer by k-ε T. C. model mostly appears before , while for the other models the peak region mostly appears after . It can also be seen that the diffusion of turbulent kinetic energy region in direction in narrow side predicted by k-ε T. C. model is slightly quicker than the models with no temperature correction.

5.2. Temperature Distribution of the Jet

Figure 14(a) shows the computed temperature distribution along the centerline of the jet from the 2D nozzle. Similar improvement on temperature core region length prediction is observed as that of velocity core region prediction. The measured temperature core length is about 0.5, while the simulated core lengths are 0.593, 0.987, 1.221, and 1.125 for k-ε T. C. model, Standard k-ε, RNG k-ε, and SST k-ω model, respectively. Although all of the four turbulence models give a longer core region prediction, the k-ε T. C. model predicted a core region only 0.2 times longer than the measured value; it produces a much better result compared to the other three models which produce a core region about 11.4 times longer than the measured value. The differences caused by different turbulence models are about 8.5%–106%, and the errors between predictions and experiments are about 19%–144%.

Figures 14(b) and 14(c) show the static temperature distribution of middle section from wide side and narrow side, respectively. In the core region, all the four turbulence models can give a reasonable result. But in the mixing region, the deviations are increased obviously, and the computational high temperature region is significantly larger than the experimental values. The computed temperature profiles by k-ε T. C. model are more flat than other nontemperature turbulence models for the rapid rate of jet decay. And from Figure 14(b), it can be seen that the results of small are not very good when the k-ε T. C. model is used. From the viewpoint of the In general, the computed temperature of the jet from 2D nozzle using k-ε T. C. model is comparatively better than the other models.

5.3. Infrared Radiation of the Jet in 90° to the Nozzle Axis

Figure 15 shows the comparison of computed and measured spectral radiation intensity in the direction 90° to the nozzle axis, and all the data are normalized. Only the maximum error bars are shown, since the small error bar will coincide with the symbols representing the experiment values. It can be observed that the spectral radiation intensity of the hot jet is mainly concentrated within a wavelength band of 4.15 μm4.6 μm. The appearance of the double peaks is caused by the emission of the high temperature combustion production of CO2 in the jet and the absorption of the low temperature CO2 in atmosphere. With the use of the standard k-ε, RNG k-ε, and SST k-ω model, the predicted peak values of the spectral radiation are significantly higher than the measured data, and the predicted wavelength band widths of the spectrum of the radiation near the wavelength of 4.5 μm are also obviously wider than the measured data. However, with the use of the k-ε T. C. model, a significantly improved prediction is observed.

Table 3 gives the comparison of computed and measured integral intensity integrated from wavelength 4 μm5 μm, where Ic represents a calculated value and IE represents a measured value. It can be seen that, for the standard k-ε model, RNG k-ε model, and SST k-ω model, the error of prediction in wide and narrow side is larger than 70.7% and 27.9%, respectively, while the k-ε T. C. model reduces the error to 26.0% and −5.8%, respectively.

6. Conclusion

The predicted lengths of the velocity and temperature core region by standard k-ε model, RNG k-ε model, and SST k-ω model are obviously longer than the measured data. For this reason, the computed infrared radiation intensities are significantly higher than the measured data. The maximum error is up to 80% or more. Because the numerical uncertainty and the experiment uncertainty are small, the errors between predictions and experiments are mainly caused by different turbulence models.

The temperature corrected k-ε T. C. model directly increases the turbulent viscosity, causing more rapid mixing. This model dramatically improves the computed velocity distribution and temperature distribution. In turn, the accuracy of the predicted infrared radiation is improved evidently.

Although the k-ε T. C. model gives a good prediction for the centerline velocity and temperature distribution of the 2D hot gas jet, it cannot make satisfactory prediction on the locations off the centerline due to the flow complexity. The computation error of the 2D hot jet infrared radiation using k-ε T. C. model is the smallest but still up to 26%. To further improve the accuracy, a more accurate turbulence model is needed.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.