New Developments in Fluid Mechanics and Its Engineering Applications
View this Special IssueResearch Article  Open Access
Numerical Simulation on the Effect of Turbulence Models on Impingement Cooling of Double Chamber Model
Abstract
Investigation of the effects of impingement cooling for the different turbulence models and study of the aerodynamic behavior of a simplified transition piece model (TP) are the two themes of this paper. A model (double chamber model) of a onefourth cylinder is designed which could simulate the transition piece structure and performance. The relative strengths and drawbacks of renormalization group theory (RNG), the realizable (RKE), the , the shear stress transport (SST), and largeeddy simulation (LES) models are used to solve the closure problem. The prediction of the inner wall temperature, cooling effectiveness, and velocity magnitude contours in various conditions are compared in five different turbulence models. Surprisingly, the and SST models can produce even better predictions of fluid properties in impinging jet flows. It is recommended as the best compromise between solution speed and accuracy.
1. Introduction
Given the large number of sustained operational hours required for industrial turbines, two important demands placed on such engines are component life and overall engine performance. These demands are somewhat conflicting because high temperatures are required at turbine combustor in order to achieve high performance; however, increasing combustor outlet temperature in turn reduced component life, high repair costs, and downtime costs. Impingement cooling is an enhanced heat transfer method capable of cooling a transition piece (TP) without injecting cool air directly into the gas chamber. Cooling the transition piece from the gas inlet enables engineers to dissipate the heat load and maintains more uniform temperatures in the turbine region needed for efficient turbine [1].
In the example of turbine cooling applications [2], impinging jet flows may be used to cool several different sections of the engine such as the combustor case (combustor can walls), turbine case/liner, and the critical high temperature turbine blades. General applications and performance of impinging jets had been discussed in a number of reviews [3–6]. The jet impingement angle has an effect on heat transfer and was studied frequently [3, 4]. Göppert et al. [5] investigated the effects of an unstable precessing jet over a fixed target plate. Hwang et al. [6] altered the flow pattern in the initial shearing layer by using coaxial jets.
Numerical modeling of impinging jet flows and heat transfer is employed widely for prediction, sensitivity analysis, and device design. Finite element, finite difference, and finite volume computational fluid dynamics (CFD) models of impinging jets have succeeded in making rough predictions of heat transfer coefficients and velocity fields. Turbulent impinging jet CFD employs practically all available numerical methods that will be critically reviewed in the following sections.
An earlier critical review of this topic was conducted by Polat et al. [7] in 1989. Since that date, the variety of numerical models have been established and computational research in predicting the physical behavior of impinging jets. Tzeng et al. [8] compared seven low Re modifications of the model using a confined turbulent slot jet array problem with three adjacent jets at . Heck et al. [9] showed the RNG model that provided a close match of Nu in the walljet region but an error up to 10% in the stagnation region. Modeling by Cziesla et al. [10] demonstrated the ability of LES to predict local Nu under a slot jet within 10% of experimental measurements. Silieti et al. [11] investigated the numerical prediction of cooling effectiveness in 2D and 3D gas turbine end walls/shrouds for the cases of conjugate and adiabatic heat transfer models. They considered different cooling hole geometries, that is, cooling slots, cylindrical, and fanshaped cooling holes at different blowing ratios. They incorporated the effects of different turbulence models in predicting the surface temperature and hence the cooling effectiveness.
There is a great interest in the application of impingement cooling to protect the transition piece from high temperature gas streams. The model created in the paper looks like a quarter torus with the curved double chambers simulating the structure of the transition piece. The relative strengths and drawbacks of renormalization group theory model (RNG), the realizable model (RKE), the model, the shear stress transport model (SST), and largeeddy simulation model (LES) for impinging jet flow and heat transfer are compared. As well, the velocity and temperature fields, in addition to centerline and twodimensional impingement cooling effectiveness, will be presented.
2. Turbulence Models
2.1. The Renormalization Group Theory Model (RNG)
The RNG turbulence model is derived from the instantaneous NavierStokes equations, using a mathematical technique called renormalization group (RNG) methods, borrowed from quantum mechanics. The analytical derivation results in a model with constants different from those in the standard model, and it results in additional terms and functions in the transport equations for the turbulent kinetic energy () and for dissipation rate (). A more comprehensive description of RNG theory and its application to turbulence can be found in [12, 13]. The governing equations for this model are as follows:
equation equation where is the convective (or substantial time) derivative, and is related to the mean strain tensor, , as , , , , , and .
2.2. The Realizable Model (RKE)
The term realizable means that the model satisfies certain mathematical constraints on the normal stresses consistent with the physics of turbulent flows. In this model, the equation is the same as in RNG model; however is not a constant and varies as a function of mean velocity field and turbulence (0.09 in loglayer , 0.05 in shear layer of ). The equation is based on a transport equation for the meansquare vorticity fluctuation [14] as where and . This model is designed to avoid unphysical solutions in the flow field.
2.3. The Model
Durbin's model, also known as the “normal velocity relaxation model,” has shown some of the best predictions to date, with calculated Nu values falling within the spread of experimental data [15]. The model uses an eddy viscosity to increase stability with two additional differential equations beyond those of the model, forming a fourequation model. The additional equations are defined as where , , , , , , , and . Similar equations exist for predicting the transport of a scalar (e.g., thermal energy) with a as a funtion of and .
2.4. The Shear Stress Transport Model (SST)
There are two major ways in which the SST model differs from the standard model. The first is the gradual change from the standard model in the inner region of the boundary layer to a high Reynoldsnumber version of the model in the outer part of the boundary layer. The second is the modified turbulent viscosity formulation to account for the transport effects of the principal turbulent shear stress. The SST (shear stress transport) model consists of the zonal (blended) equations and clips of turbulent viscosity so that turbulent stress stays within what is dictated by structural similarity constant. The equation is the same as the standard model whereas the resulting blended equation for is [16] where in the inner layer and in the outer layer, and .
2.5. LargeEddy Simulation (LES)
Largeeddy simulation (LES) is a technique intermediate between the direct simulation of turbulent flows and the solution of the Reynoldsaveraged equations. In LES the contribution of the large, energycarrying structures to momentum and energy transfer is computed exactly; well only the effect of the smallest scales of turbulence is modeled. Since the small scales tend to be more homogeneous and universal and less affected by the boundary conditions than the large ones, there is hope that their models can be simpler and require fewer adjustments when applied to different flows than similar models for the RANS equations [17, 18].
3. Computational Model, Boundary Conditions, and Grid
3.1. Model Description
Transition piece develops heat transformation on both internal and external walls to eliminate resonant frequency concerns. As well the transition piece conducts gas flow directly from the corresponding combustion liners toward the first stage of the gas turbine (stator). The discrete coolant jets, forming a protective film chamber on the side of transition piece, are drawn from the upstream compressor in an operational gas turbine engine. From the supply plenum, the coolant is ejected through the several rows of discrete holes over the external boundary layer against the local high thermal conduction on the other side of the transition pieces [19]. The model, as a onefourth cylinder, could simulate the transition piece’s structure and performance, to elaborate the turbulence effect on the impingementcooling performance over a concave surface (see Figure 1) [4].
A schematic diagram of the flow domain along with boundary conditions and dimensions is given in Figure 2(a). As shown in the figure, the onefourth cylinder model has two layers of chambers with a length of 1050 mm and an outer radius and an inner radius of 200 mm and 162 mm, respectively. In the diagram, one side of the outer chamber called the coolant chamber is closed; contrarily, both sides of the mainstream chamber as the inner chamber are unfolding in which gas could flow through from one side to the other. There are 18 holes distributed uniformly in three rows on the surface of the outer wall. The distance between the two rows is 68 mm and the diameter of all the holes is about 10.26 mm. In this scenario, the velocity of the coolant flow is set at 6 m/s; the temperatures of the coolant flow and the mainstream flow are set as 300 K and 1300 K, respectively.
(a)
(b)
3.2. Boundary Conditions
Boundary conditions are applied to specific faces within the domain to specify the flow and thermal variables that dictate conditions within the model. Figure 2(a) discloses the boundary conditions used in the model in which the cooling air and gas are moving along, respectively, in the two layers of chambers in opposite directions. In the cooling chamber, the simulation is performed by using air as the cooling flow, while velocity and temperature are set on the jet holes, with pressure on the exit. In another chamber, it is assumed that the mainstream is a mixture of O_{2}, H_{2}O, CO_{2}, N_{2}, and some rare gases. The model created is considered as boundary condition (Table 1) [4]. Assumption of the solid wall of the quarter torus is formed with a hypothesis of negligible thermal resistance by conduction; the thermal properties of the material were considered by Nimonic 263.

3.3. Meshing and Simulation Procedures
The computational domain incorporates the model, the HEXA mesh in the software, ICEM/CFD, used to generate the structured multiblock and the bodyfitted grid system. In this study, the grid system associated with the parts of the mainstream and the coolant supply plenum is Htype. Figure 3 shows the grids of the computational domain, and the total number of the cells for the 3D domain is 776,828. Local grid refinement is used near the hole regions.
The wall plus (7) is a nondimensional number similar to local Reynolds number, determining whether the influences in the walladjacent cells are laminar or turbulent, hence indicating the part of the turbulent boundary layer that they resolve:
The subdivisions of the nearwall region in a turbulent boundary layer can be summarized as follows [20]:(a): in the viscous sublayer region (velocity profiles are assumed to be laminar and viscous stress dominates the wall shear),(b): buffer region (both viscous and turbulent shear dominate),(c): fully turbulent portion or loglaw region (corresponds to the region where turbulent shear predominates).
As it can be seen in Figure 4, for all cases, all nodes on the inner wall surface have the plus value smaller than 300.
This study is using a commercial CFD code based on the controlvolume method, ANSYSFLUENT 12.0.16, which in order to predict temperature, impingementcooling effectiveness, and velocity fields. All runs were made on a PC cluster with sixteen Pentium4 3.0 GHz personal computers. The convergence criteria of the steadystate solution are judged by the reduction in the mass residual by a factor of 6, typically, in 2000 iterations.
4. Results and Discussion
4.1. Inner Wall Temperature
Figure 5 shows the temperature contours predicted by five turbulence models. The figure illustrates that the temperature at the starting point of the wall is high and then starts to go down. The starting point is cooled by the coolant holes at the curve mm on the outer wall while the same temperature is maintained throughout the coolant holes. After the coolant strikes the inner wall, there are vortexes formed. The jet impingement and the vortex formed out of the coolant flow cool the surface of the inner wall. The delicate color region (orange, yellow, and green region) is approximately the region where the coolant strikes the inner wall after being reflected by the coolant holes, in which the inner wall is cooled purely by impingement cooling. These figures confirm that the model has well simulation of heat transfer in the double chamber model. The calculated temperature at the center line of the inner wall is presented in Figure 5(b); there is a difference among the predictions of SST, LES, and RKE models. Both the temperature contours and the color lines have the same trend in the five turbulence models whereas the calculated temperatures between and 600 mm are diverse for all turbulence models. All turbulence models continue to predict similar wall temperatures, while the RNG model underpredicts the temperature by up to 5%–8% in the same location.
(a)
(b)
4.2. Cooling Effectiveness
To define cooling effectiveness, the surface temperature downstream of the cooling hole has to be measured. The adiabatic cooling effectiveness () is used to examine the performance of cooling. The definition of is where is the mainstream hot gas inlet temperature, which is a fixed value for calculation of the adiabatic cooling effectiveness of any location, and is the temperature of the coolant, which is assigned as a constant of 300 K in this issue. is the adiabatic wall temperature [21].
Figures 6(a) and 6(b) exhibit the comparison of cooling effectiveness on distribution and averaged cooling effectiveness along the axis, respectively. It could be concluded that the effectiveness is high at the beginning of the film while decreasing gradually in downstream. The distribution of the cooling effectiveness of the and SST model cases is significantly different. The centerline effectiveness by using the five turbulence model cases is shown in Figure 6(b). Overall, particularly the model gave better results compared to other models. Values of cooling effectiveness compare with the RNG model case, there is an approximately 30% difference in model case.
(a)
(b)
4.3. Characteristics of the Flow Field
Since the thermal field of a jetincrossflow interaction is dictated by the hydrodynamics, the flow field results were predicted by five turbulence models. Figure 7(a) shows the velocity contours at the section plane of coolant chamber. As it can be seen, all models predicted the low momentum region along the downstream edge and the corresponding high momentum or jetting region along the upstream edge within the cooling holes. The predicted reattachment region by SST and model cases is approximately at the inlet of the cooling chamber, whereas LES, RKE, and RNG predict a larger separation bubble and the flow separates behind the cooling holes. The computed near inner wall velocity contours (m/s) along the centerline plane are shown in Figure 7(b), where the turbulence closure was simulated using the five different turbulence models. All predictions are extremely close to each other with a skewed upstream velocity profile.
(a)
(b)
The actual computational cost will of course vary with model complexity and computing power. With the parallel computing resources of a desktop computer available at the time of writing, six Pentium4 typically 3 GHz processors, for a highresolution twodimensional problem, the steady timeaveraged eddy viscosity models (RNG, RKE) will have computation times of a few hours (0.5–1.5). In comparison, the more complex SST and could take 2–4 hours depending on how smoothly the model converges. Based on recent work, unsteady LES models have computation times at least two orders of magnitude higher; a wellresolved threedimensional LES impinging jet model could take a day to provide a solution.
5. Conclusion
A numerical simulation has been performed to study the flow and heat transfer of impinging cooling on the double chamber model, and a comparative study, indicating the ability of five turbulence models, is presented. The research of turbulence model tasks is important to improve the design and resulting performance of impinging jets.
During this investigation, numerical simulation is impacted with five turbulence models, which has some practical value for real processing and guiding significance for theory. To date, the SST and models offer the best results for the least amount of computation time. It is very important to consider the effect of heat conduction within the metal on the predictions of an accurate surface temperature and hence impingementcooling effectiveness. The validation of the present study has confirmed angle cases and will be employed in future studies of impingementcooling parameters optimization.
Nomenclature
:  Diameter of coolant chamber 
:  Diameter of mainstream chamber 
:  Length of the model 
:  Absolute static temperature 
:  Nondimensional coordinates in diameter, spanwise, and mainstream directions. 
:  Cooling effectiveness. 
:  Mainstream flow 
:  Coolant flow 
:  Adiabatic wall. 
Conflict of Interests
The authors declare that there is no conflict of interests.
Acknowledgment
This research is supported by the Technology Development of Jilin Province (no. 20126001).
References
 J. C. Han, S. Dutta, and S. V. Ekkad, Gas Turbine Heat Transfer and Cooling Technology, Taylor and Francis, New York, NY, USA, 2000.
 N. Zuckerman and N. Lior, “Impingement heat transfer: correlations and numerical modeling,” Journal of Heat Transfer, vol. 127, no. 5, pp. 544–552, 2005. View at: Publisher Site  Google Scholar
 A. A. Tawfek, “Heat transfer studies of the oblique impingement of round jets upon a curved surface,” Heat and Mass Transfer, vol. 38, no. 6, pp. 467–475, 2002. View at: Publisher Site  Google Scholar
 Z. L. Yu, T. Xu, J. L. Li, L. Ma, and T. S. Xu, “Comparison of a series of double chamber model with various hole angles for enhancing cooling effectiveness,” International Communications in Heat and Mass Transfer, vol. 44, pp. 38–44, 2013. View at: Publisher Site  Google Scholar
 S. Göppert, T. Gürtler, H. Mocikat, and H. Herwig, “Heat transfer under a precessing jet: effects of unsteady jet impingement,” International Journal of Heat and Mass Transfer, vol. 47, no. 1213, pp. 2795–2806, 2004. View at: Publisher Site  Google Scholar
 S. D. Hwang, C. H. Lee, and H. H. Cho, “Heat transfer and flow structures in axisymmetric impinging jet controlled by vortex pairing,” International Journal of Heat and Fluid Flow, vol. 22, no. 3, pp. 293–300, 2001. View at: Publisher Site  Google Scholar
 S. Polat, B. Huang, A. S. Mujumdar, and W. J. M. Douglas, “Numerical flow and heat transfer under impinging jets: a review,” Annual Review of Numerical Fluid Mechanics and Heat Transfer, vol. 2, pp. 157–197, 1989. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. Y. Tzeng, C. Y. Soong, and C. D. Hsieh, “Numerical investigation of heat transfer under confined impinging turbulent slot jets,” Numerical Heat Transfer A, vol. 35, no. 8, pp. 903–924, 1999. View at: Publisher Site  Google Scholar
 U. Heck, U. Fritsching, and K. Bauckhage, “Fluid flow and heat transfer in gas jet quenching of a cylinder,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 11, no. 1, pp. 36–49, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 T. Cziesla, G. Biswas, H. Chattopadhyay, and N. K. Mitra, “Largeeddy simulation of flow and heat transfer in an impinging slot jet,” International Journal of Heat and Fluid Flow, vol. 22, no. 5, pp. 500–508, 2001. View at: Publisher Site  Google Scholar
 M. Silieti, E. Divo, and A. J. Kassab, “Film cooling effectiveness from a single scaledup fanshaped hole a CFD simulation of adiabatic and conjugate heat transfer models,” in Proceedings of the ASME Turbo Expo: Power for Land, Sea, and Air, pp. 431–441, June 2005, paper no. GT200568431. View at: Google Scholar
 S. H. Park, S. W. Park, S. H. Rhee, S. B. Lee, J.E. Choi, and S. H. Kang, “Investigation on the wall function implementation for the prediction of ship resistance,” International Journal of Naval Architecture and Ocean Engineering, vol. 5, no. 1, pp. 33–46, 2013. View at: Publisher Site  Google Scholar
 S.E. Kim and S. H. Rhee, “Efficient engineering prediction of turbulentwing tip vortex flows,” Computer Modeling in Engineering and Sciences, vol. 62, no. 3, pp. 291–309, 2010. View at: Google Scholar
 Z. L. Yu, T. Xu, J. L. Li, T. S. Xu, and Y. Tatsuo, “Computational analysis of droplet mass and size effect on mist/air impingement cooling performance,” Advances in Mechanical Engineering, vol. 2013, Article ID 181856, 8 pages, 2013. View at: Publisher Site  Google Scholar
 M. Behnia, S. Parneix, Y. Shabany, and P. A. Durbin, “Numerical study of turbulent heat transfer in confined and unconfined impinging jets,” International Journal of Heat and Fluid Flow, vol. 20, no. 1, pp. 1–9, 1999. View at: Publisher Site  Google Scholar
 T. H. Park, H. G. Choi, J. Y. Yoo, and S. J. Kim, “Streamline upwind numerical simulation of twodimensional confined impinging slot jets,” International Journal of Heat and Mass Transfer, vol. 46, no. 2, pp. 251–262, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Y. Huang, S. Wang, and V. Yang, “Systematic analysis of leanpremixed swirlstabilized combustion,” AIAA Journal, vol. 44, no. 4, pp. 724–740, 2006. View at: Publisher Site  Google Scholar
 V. Moureau, P. Domingo, and L. Vervisch, “From largeeddy simulation to direct numerical simulation of a lean premixed swirl flame: filtered laminar flamePDF modeling,” Combustion and Flame, vol. 158, no. 7, pp. 1340–1357, 2011. View at: Publisher Site  Google Scholar
 J. Benoit, C. Johnston, and M. Zingg, “Enhancing gas turbine power plant profitability: chronic transition piece and turbine part failures in some 501F gas turbines led to a replacement part redesign,” Power Engineering, vol. 111, no. 11, pp. 140–144, 2007. View at: Google Scholar
 M. Moshfeghi, Y. J. Song, and Y. H. Xie, “Effects of nearwall grid spacing on SSTKω model using NREL Phase VI horizontal axis wind turbine,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 107, pp. 94–105, 2012. View at: Publisher Site  Google Scholar
 H. Wu, J. Wang, and Z. Tao, “Evaluation of predicted heat transfer on a transonic vane using ${v}^{2}f$ turbulence models,” Journal of Thermal Science and Technology, vol. 6, no. 3, pp. 424–435, 2011. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Zhenglei Yu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.