| 研究生: |
謝銘峻 Hsieh, Ming-June |
|---|---|
| 論文名稱: |
含吸力面薄膜冷卻渦輪靜子葉片的熱傳模擬 Numerical Thermal Simulation of a Turbine Nozzle with Suction Side Film Cooling |
| 指導教授: |
陳世雄
Chen, Shih-Hsiung |
| 學位類別: |
博士 Doctor |
| 系所名稱: |
工學院 - 航空太空工程學系 Department of Aeronautics & Astronautics |
| 論文出版年: | 2008 |
| 畢業學年度: | 96 |
| 語文別: | 英文 |
| 論文頁數: | 183 |
| 中文關鍵詞: | 熱傳模擬 、CFD 、渦輪葉片 、薄膜冷卻 |
| 外文關鍵詞: | CFD, heat transfer simulation, film cooling, turbine blades |
| 相關次數: | 點閱:116 下載:1 |
| 分享至: |
| 查詢本校圖書館目錄 查詢臺灣博碩士論文知識加值系統 勘誤回報 |
為了提高渦輪引擎燃油使用效率與輸出功率,工程上普遍透過提高渦輪進口溫度來達到這個目標,然而也因燃氣溫度遠高於渦輪葉片所能承受的範圍,因此需要導引壓縮機的低溫氣流來冷卻以保護葉片。本研究的目的是利用數值方法來模擬分析含薄膜冷卻渦輪靜子葉片的熱傳,並探討冷卻氣流量與壁面熱傳的關係。使用的數值計算工具為CFX-TASCflow軟體,為了驗證數值計算工具的準確度,研究中首先針對二維含肋條管道、旋轉中的方型管與無薄膜冷卻渦輪靜子葉片的熱傳現象進行計算分析,並與相關實驗數據比較後得到不錯的結果。在旋轉方形管熱傳模擬中,為了減少因原模型方型轉角造成流體急速轉彎分離而形成的壓力損失,將方管轉角改成圓弧型來模擬其流場。雖然圓弧型轉角成功地減少方形管內的壓力損失,但轉角處與第二個通道中的壁面熱傳卻也因此有了些許的下降。
在含薄膜冷卻渦輪葉片熱傳模擬研究中,為了準確模擬壁面熱傳分布,紊流的計算採用SST紊流模型,計算區域則包含葉片流道、冷卻噴孔與氣室。計算模型為一渦輪靜子葉片,並在吸力面上鑽有兩排相鄰交錯的圓錐形噴孔。模擬結果發現吹出比m=0.45時的壁面平均熱傳係數分布趨勢與實驗相符,然而在吹出比m=0.6與m=1.0時噴孔附近的熱傳分布則有較大的差異,即使在遠離噴孔處有較接近的結果,在吹出比m=1.0時噴孔附近的熱傳結果甚至比實驗值高。在高吹出比下數值結果與實驗值的差異應來自於熱傳模擬中冷卻氣流已抬升遠離壁面,但在實驗中並無此現象的結果。實驗中此冷卻氣流抬升的現象可能在吹出比大於1.0的情況下發生,但在此熱傳模擬相對較低的吹出比範圍內卻提早發生。
為了探討冷卻氣流量與壁面熱傳的關係,本研究熱傳模擬的吹出比範圍為m=0.15-1.0。模擬結果發現吹出比越高,遠離噴孔的平均壁面熱傳係數會越低,這是由於保護葉片不受外界熱氣影響的冷卻氣流量增加的緣故。但在靠近噴孔處,平均壁面熱傳係數則呈現先降低後增加的趨勢。除此之外,在吹出比大於0.45之後平均壁面熱傳係數的分佈會有一峰值與谷值出現。在分析過靠近噴孔處的溫度場分佈後,發現冷卻氣流在噴出之後的下游處形成一對轉渦旋,且渦旋的尺度隨著吹出比的增加而增大。而在m=1.0時帶有高動量的冷卻氣流被噴離了壁面,使得外界熱氣鑽入冷卻氣流的下方來接觸壁面。這代表過高的吹出比會導致氣流抬升現象,使得靠近噴孔處的熱傳增加、冷卻效率下降。這也表示了增加吹出比並不完全是提昇冷卻效率的好方法,選擇一個適當的吹出比來操作可更有效率地冷卻渦輪葉片。另外,因冷卻氣流噴出所形成的對轉渦旋而被捲入的外界熱氣也會導致冷卻效率的下降。由於冷卻氣流的抬升與外界熱氣的捲入會造成在高吹出比時葉片壁面熱傳的增加與冷卻效率的下降,如何找出較佳的冷卻方式來避免這些現象的發生將是未來的研究重點。
Raising the turbine entry temperature is a common method to increase the performance and power output of a gas turbine engine. However, the level of turbine entry temperature is much higher than the blade material can endure. Thus, it is necessary to guard the blade against the entering hot gas by guiding cooling air from the compressor into the turbine. The purpose of this research is to numerically simulate the flowfield and temperature distribution of gas turbine blades with film cooling, and then to investigate the influence of coolant quantity on the blade heat transfer. The finite element based finite volume method CFD software, CFX-TASCflow, was adopted as the numerical tool. To validate the numerical tool, thermal simulations in a 2D channel with symmetric ribs on the walls, a rotating square channel, and a turbine nozzle without film cooling were conducted, and the numerical results were in good agreement with the experimental measurements. In the thermal simulations of a rotating square channel, the models with round corners were also tested to reduce the pressure loss resulting from the separations at the sharp corners in the original model. Though it was successful in decreasing the pressure loss by replacing the sharp corners with the round corners, the wall heat transfer was also reduced at the turn and in the second pass with the modifications.
In order to improve the numerical prediction of heat transfer on a turbine blade with film cooling, the SST turbulence model was adopted to simulate the turbulent flow. The computational domain included the blade passage, the coolant holes, and the plenum. The test model was a turbine nozzle blade and two staggered rows of conical holes drilled on the suction side as the cooling holes. The results indicated that the averaged heat transfer coefficient distribution agreed with the experimental measurements satisfactorily at a blowing ratio of m=0.45. Noticeable discrepancies on are discovered near the rear hole at m=0.6 and m=1.0, even though the results agree with the measured data downstream away from the holes. The average heat transfer coefficient in the region near the holes is even overpredicted at m=1.0. The discrepancies resulted from the lift-off phenomenon which occurred in the heat transfer prediction at high blowing ratios, but not in the heat transfer measurements. It would be discovered at a higher blowing ratio than 1.0, but occurs too early at this relatively low blowing ratio in this computation.
In order to investigate the influence of the blowing ratio on the turbine blade heat transfer, the numerical predictions were performed in the range of a blowing ratio m=0.15-1.0. It was found that with increasing the blowing ratio, the average heat transfer coefficient decreased on the wall downstream away from the holes because the cooling air injected for wall protection from the ambient hot gas was increased. But the average heat transfer coefficient decreased first and then increased in the vicinity of the holes. In addition, a peak and a valley value also appeared in the average heat transfer coefficient distribution as m>0.45. After analyzing the thermal field in the vicinity of the holes, it was discovered that a pair of counter-rotating vortices formed while the cooing air was ejected out of the holes. The scales of the induced vortices became greater as the blowing ratio was increased. At m=1.0, the cooling air with high momentum lifts off the surface, and the ambient hot gas completely flows below the cooling air. This implies that an excess blowing ratio causes wall heat transfer to increase while poor cooling effectiveness near the film cooling holes is due to the lift-off phenomenon. This also means that increasing the blowing ratio is not always a good way to improve cooling at any given position on the suction side. Choosing an adequate blowing ratio can make the cooling air cool the turbine blades more effectively. In addition, the entrainment by the induced counter-rotating vortices also degrades the cooling effectiveness near the film cooling holes. Cooling air lift-off and hot gas entrainment by induced vortices result in high wall heat transfer and poor cooling effectiveness at high blowing ratios. Finding the cooling configuration to prevent cooling air lift-off and hot gas entrainment will be the focus for the future research.
[1] Han, J. C., “ Heat Transfer and Friction in Channels with Two Opposite
Rid-Roughened Walls,” ASME Journal of Heat Transfer, Vol. 106, pp. 774-782, 1984.
[2] Han, J. C.,“Heat Transfer and Friction Characteristics in Rectangular Channels with
Rid-Roughened Turbulators,” ASME Journal of Heat Transfer, Vol. 110, pp. 321-328,
1988.
[3] Liou, T. M., and Lin, J.,“Measurements of Turbulent Flow in a Duct with Repeated
Ribs Applied to Two Opposite Walls,” Journal of Chinese Institute Engineers, Vol. 11,
No. 4, pp. 319-326, 1988.
[4] Liou, T. M., and Hwang, J. J.,“Effect of Ridge Shapes on Turbulent Heat Transfer
and Friction in Rectangular Channel,” International Journal of Heat and Mass
Transfer, Vol. 36, No. 8, pp. 931-940, 1993.
[5] Patankar, S. V., Liu, C. H., and Sparrow, E. M.,“Fully Developed Flow and Heat
Transfer in Ducts Having Streamwise-Periodic Variations of Cross-Sectional Area,”
ASME Journal of Heat Transfer, Vol. 99, pp. 180-186, 1977.
[6] Acharya, S., Dutta, S., Myrum, T. A., and Baker, R. S.,“Periodically Developed
Flow and Heat Transfer in a Ribbed Duct,”International Journal of Heat Mass
Transfer, Vol. 36, No. 8, pp. 2069-2082, 1993.
[7] Liou, T. M., Hwang, J. J., and Chen, S. H.,“Simulation and Measurement of
Enhanced Turbulent Heat Transfer in a Channel with Periodic Ribs on One Principal
Wall,”International Journal of Heat and Mass Transfer, Vol. 36, No. 4, pp. 507-517,
1993.
[8] Wagner, J. H., Johnson, B. V., and Kopper, F. C.,“Heat Transfer in Rotating
Serpentine Passages with Smooth Walls,” ASME Journal of Turbomachinery, Vol. 113,
pp. 321-330, 1991.
[9] Dutta, S., and Han, J. C.,“Local Heat Transfer in Rotating Smooth and Ribbed
Two-Pass Square Channels with Three Channel Orientations,” ASME Journal of
Heat Transfer, Vol. 118, pp. 578-584, 1996.
[10] Liou, T. M., and Chen, C. C.,“LDV Study of Developing Flows Through a Smooth
86
Duct with a 180 Deg Straight-Corner Turn,”ASME Journal of Turbomachinery, Vol.
121, pp. 167-174, 1999.
[11] Kirillov, A. I., Ris, V. V., Smirnov, E. M. and Zaitsev, D. K.,“Numerical Simulation
of Local Heat Transfer in Rotating Two-Pass Square Channels,” Heat Transfer in Gas
Turbine Systems, pp. 456-463, 2000.
[12] Chen, H. C., Jang, Y. J., and Han, J. C.,“Computation of Heat Transfer in Rotating
Two-Pass Square Channels by a Second-Moment Closure Model,”International
Journal of Heat and Mass Transfer, Vol. 43, pp. 1603-1616, 2000.
[13] Camci, C. and Arts, T.,“An experiment Convective Heat Transfer Investigation
Around a Film-Cooled Gas Turbine Blade,”Journal of Turbomachinery, Vol. 112, pp.
497-503, 1990.
[14] Consigny, H., and Richards, B. E.,“Short Duration Measurements of Heat Transfer
Rate to a Gas Turbine Rotor Blade,”ASME Journal of Engineering for Power, Vol.
104, No.3, pp. 542-551, 1982.
[15] Medic, G., and Durbin, P. A.,“Toward Improved Prediction of Heat Transfer on
Turbine Blades,”Journal of Turbomachinery, Vol. 124, pp. 187-192, 2002.
[16] Calzada, P. D. L., and Alonso, A.,“Numerical Investigation of Heat Transfer in
Turbine Cascades with Separated Flows,” ASME Journal of Turbomachinery, Vol. 125,
pp. 260-266, 2003.
[17] Hermanson, K., Kern S., Picker, G., and Parneix, S.,“Prediction of External Heat
Transfer for Turbine Vanes and Blades with Secondary Flowfields,”ASME Journal
of Turbomachinery, Vol. 125, pp. 107-113, 2003.
[18] Goldstein, R. J., Eckert, E. R. G., and Ramsey, J., W.,“Film Cooling With Injection
Through a Circular Hole,”ASME Journal of Engineering for Power, Vol. 90, pp.
384-395, 1968.
[19] Eriksen, V. L., and G.,Goldstein, R. J.,“Heat Transfer and Film Cooling Injection
Through Inclined Circular Tubes,”ASME Journal of Heat Transfer, Vol. 96, pp.
239-245, 1974.
[20] Eriksen, V. L., and G.,Goldstein, R. J.,“Heat Transfer and Film Cooling Injection
Through a Round Hole,”ASME Paper No. 74-GT-6, 1974.
[21] Bergeles, G., Gosman, A. D., and Launder, B. E.,“The Near-Field Character of a Jet
87
Discharged Normal to a Mainstream,” ASME Journal of Heat Transfer, Vol. 90, pp.
384-395, 1968.
[22] Bergeles, G., Gosman, A. D., and Launder, B. E., “ The Prediction of
Three-Dimensional Discrete-Hole Cooling Processes, Part I: Laminar Flow,” ASME
Journal of Heat Transfer, Vol. 98, pp. 379-386, 1976.
[23] Launder, B. E., and York, J.,“Discrete Hole Cooling in the Presence of Freestream
Turbulence and Strong Favourable Pressure Gradient,”Int. J. Heat Mass Transfer,
Vol. 17, pp. 1403-1409, 1974.
[24] Kadotani, K., and Goldstein, R. J.,“Effect of Mainstream Variables on Jet Issuing
From a Row of Inclined Round Holes,”ASME Journal of Engineering for Power, Vol.
101, pp. 298-304, 1979.
[25] Andreopoulos, J., and Rodi, W., “ Experimental Investigation of Jets in a
Crossflow,”Journal of Fluid Mechanics, Vol. 138, pp. 92-127, 1984.
[26] Petersen, D. R., Eckert, E. R. G., and Goldstein, R. J.,“Film Cooling With Large
Density Difference Between the Mainstream and Secondary Fluid Measured by
Heat-Mass Transfer Analogy,” ASME Journal of Heat Transfer, Vol. 99, pp. 620-627,
1977.
[27] Foster, N. W., and Lampard, D., “The Flow and Film Cooling Effectiveness
Following Injection Through a Row of Holes,” ASME Journal of Engineering for
Power, Vol. 102, pp 584-588, 1980.
[28] Sinha, A. K., Bogard, D. G., and Crawford, M. E.,“Film–Cooling Effectiveness
Downstream of a Single Row of Holes With Variable Density Ratio,” ASME Journal
of Turbomachinery, Vol. 113, pp. 442-449, 1991.
[29] Ammari, H. D., Hay, N., and Lampard, D.,“The Effect of Density Ratio on the Heat
Transfer Coefficient from a Film Cooled Flat Plate,” ASME Journal of
Turbomachiney, Vol. 113, pp. 464-470, 1990.
[30] Pietrzyk, J. R., Bogard, D. G., and Crawford, M. E.,“Hydrodynamic Measurements
of Jets in Crossflow for Gas Turbine Film Cooling Applications,”ASME Journal of
Turbomachinery, Vol. 111, pp. 139-145, 1989.
[31] Goldstein, R. J., and Taylor, J. R.,“The Influence of a Laminar Boundary Layer and
Laminar Injection on Film Cooling Performance,” ASME Journal of Heat Transfer,
Vol. 104, pp. 715-721, 1982.
88
[32] Baldauf, S., Scheurlen, M., Schlz, A., and Wittig, S.,“Correlation of Film-Cooling
Effectiveness from Thermographic Measurements at Enginelike Conditions,”ASME
J. Turbomachinery, Vol. 124, pp. 686-698, 2002.
[33] Hay, N., Lampard, D., and Saluja, C. L.,“Effect of Cooling Films on the Heat
Transfer Coefficient on a Flat Plate with Zero Mainstream Pressure Gradient,” ASME
Journal of Heat Transfer, Vol. 104, pp. 715-721, 1985.
[34] Gritsch, M., Schulz, A., and Wittig, S.,“Adiabatic Wall Effectiveness Measurements
of Film Cooling Holes with Expanded Exits,” ASME Journal of Turbomachinery, Vol.
120, pp. 549-556, 1998b.
[35] Saumweber, C., Schulz, A., and Wittig, A.,“Free-Stream Turbulence Effects on Film
Cooling With Shaped Holes,” ASME Journal of Turbomachinery, Vol. 125, pp. 65-73,
2003.
[36] Ligrani, P. M., Wigle, J. M., Ciriello, S., and Jackson, S., M.,“Film-Cooling from
Holes with Compound Angle Orientations. Part I: Results Downstream of Two
Staggered Rows of Holes with 3d Spanwise Spacing,” ASME Journal of Heat
Transfer, Vol. 116, pp. 341-352, 1994a.
[37] Sen, B., Schmidt, D. L., and Bogard, D. G.,“Film Cooling with Compound Angle
Holes: Heat Transfer,” ASME Journal of Turbomachinery, Vol. 118, pp. 800-806,
1996.
[38] Ekkad, S. V., Zapata, D., and Han, J. C.,“Heat Transfer Coefficient Over a Flat
Surface with Air and CO2 Film Injection Through Compound Angle Holes Using a
Transient Liquid Crystal Image Method,” ASME Journal of Turbomachinery, Vol. 119,
pp. 580-586, 1997a.
[39] Ekkad, S. V., Zapata, D., and Han, J. C.,“Film Effectiveness Over a Flat Surface
with Air and CO2 Film Injection Through Compound Angle Holes Using a Transient
Liquid Crystal Image Method,” ASME Journal of Turbomachinery, Vol. 119, pp.
594-600, 1997b.
[40] Maska, J. F., Fish, R. W., and Suo, M.,“The Additive Nature of Film Cooling from
Rows of Holes,” ASME Journal of Engineering for Power, Vol. 98, pp. 457-464,
1976.
[41] Sasaki, M., Takahara, K., Kumagai, T., and Hamano, M., “ Film Cooling
Effectiveness for Injection from Multirow Holes,” ASME Journal of Engineering for
Power, Vol. 101, pp. 101-108, 1979.
89
[42] Ito, S., Goldstein, R. J., and Eckert, E. R. G.,“Film Cooling of a Gas Turbine Blade,”
ASME Journal of Engineering for Power, Vol. 100, pp. 476-481, 1978.
[43] Schwarz, S. G., and Goldstein, R. J.,“The Two-Dimensional Behavior of Film
Cooling Jets on Concave Surfaces,” ASME Journal of Turbomachinery, Vol. 111, pp.
124-130, 1989.
[44] Schwarz, S. G., Goldstein, R. J., and Eckert, E. R. G.,“The Influence of Curvature on
Film Cooling Performance,” ASME Journal of Turbomachinery, Vol. 113, pp. 472-478,
1991.
[45] Goldstein, R. J., and Stone, L. D.,“Row-of-Holes Film Cooling of Curved Walls at
Low Injection Angles,” ASME Journal of Turbomachinery, Vol. 119, pp. 574-579,
1997.
[46] Chen, P. H., Hung, M. S., and Ding, P. P.,“A Transient Method Using Liquid Crystal
for Film Cooling Over a Convex Surface,” International Journal of Rotating
Machinery, Vol. 7, No. 3, pp. 153-164, 2001.
[47] Behle, M., K., and Patankar, S. V.,“Curvature Effects on Discrete-Hole Film
Cooling,” ASME Journal of Turbomachinery, Vol. 121, pp. 792-803, 1999a.
[48] Camci, C., and Arts, T.,“An Experimental Convective Heat Transfer Investigation
Around a Film-Cooled Gas Turbine Blade,” ASME Journal of Turbomachinery, Vol.
112, pp. 497-503, 1990.
[49] Arts, T., and Bourguignon, A. E.,“Behavior of Coolant Film With Two Rows of
Holes Along the Pressure Side of a High-Pressure Nozzle Guide Vane,” ASME
Journal of Turbomachinery, Vol. 112, pp. 497-503, 1990.
[50] Drost, U., and Bölcs, A.,“Investigation of Detailed Film Cooling Effectiveness and
Heat Transfer Distributions on a Gas Turbine Airfoil,” ASME Journal of
Turbomachinery, Vol. 121, pp. 233-242, 1999.
[51] Crawford, M. E., Kays, W. M., and Moffat, R. J.,“Full-Coverage Film Cooling, Part
I: Comparison of Heat Transfer Data for Three Injection Angles,”ASME Journal of
Engineering for Power, Vol. 102, pp. 1000-1005, 1980.
[52] Tafti, D. K., and Yavuzkurt, S.,“Prediction of Heat transfer Characteristics for
Discrete Hole Film Cooling for Turbine Blade Applications,” ASME Journal of
Turbomachinery, Vol. 112, pp. 504-511, 1990.
90
[53] Bergeles, G., Gosman, A. D., and Launder, B. E.,“Double-Row Discrete-Hole
Cooling: an Experimental and Numerical Study,” ASME Journal of Engineering for
Power, Vol. 102, pp. 498-503, 1980.
[54] Ameri, A. A., and Arnone, A.,“Transition of Modeling Effects on Turbine Rotor
Blade Heat Transfer Predictions,” ASME Journal of Turbomachinery, Vol. 118, pp.
307-314, 1996.
[55] Leylek, J. H., and Zerkle, R. D.,“Discrete-Jet Film Cooling: A Comparison of
Computational Results with Experiments,” ASME Journal of Turbomachinery, Vol.
116, pp. 358-368, 1994.
[56] Lin, Y. L., and Shih, I. P.,“Film Cooling over Flat, Convex, and Concave Surfaces,”
AIAA Paper 99-0344, 1999.
[57] Dorney, D. J., and Davis, R. L.,“Numerical Simulation of Turbine Hot Spot
Alleviation Using Film Cooling,” Journal of Propulsion and Power, Vol. 9, pp.
329-336, 1993.
[58] Abhari, R. S.,“Impact of Rotor-Stator Interaction on Turbine Blade Film Cooling,”
ASME Journal of Turbomachinery, Vol. 118, pp. 123-133, 1996.
[59] Garg, V. K., and Gaugler, R. E.,“Effect of Coolant Temperature and Mass Flow on
Film Cooling of Turbine Blades,” Int. J. Heat Mass Transfer, Vol. 40, No. 2, pp.
435-445, 1997.
[60] Garg, V. K., and Gaugler, R. E.,“Effect of Velocity and Temperature Distribution at
the Hole Exit on Film Cooling of Turbine Blades,” ASME Journal of Turbomachinery,
Vol. 119, pp. 343-351, 1997.
[61] Medic, G., and Durbin, P. A.,“Toward Improved Film Cooling Prediction,” ASME
Journal of Turbomachinery, Vol. 124, pp. 193-199, 2002.
[62] Heidmann, J. D., Rigby, D. L., and Ameri, A. A.,“A Three-Dimensional Coupled
Internal/External Simulation of a Film-Cooled Turbine Vane,” ASME Journal of
Turbomachinery, Vol. 122, pp. 348-359, 2000.
[63] Shina, A. K., Bogard, D. G., and Crawford, M. E.,“Film-Cooling Effectiveness
Downstream of a Single Row of Holes with Variable Density Ratio,” ASME Journal
of Turbomachinery, Vol. 113, pp. 442-449, 1991.
[64] Goldstein, R. J., and Yoshida, T.,“The Influence of a Laminar Boundary Layer and
91
Laminar Injection on Film Cooling Performance,” ASME Journal of Heat Transfer,
Vol. 104, pp. 715-721, 1982.
[65] Fric, T. F., and Roshko, A.,“Vortical Structure in the Wake of a Transverse Jet,”
Journal of Fluid Mechanics, Vol. 279, pp. 1-47, 1994.
[66] Crabb, D., Durao, D. F. G., and Whitelaw, J. H.,“A round Jet Normal to a
Crossflow,” ASME Journal of Fluids Engineering, Vol. 103, pp. 142-152, 1981.
[67] Chima, R. V., Giel, P. W., and Boyle, R. J.,“An Algebraic Turbulence Mod el for
Three-Dimensional Viscous Flows,” AIAA paper 93-0083 (NASA TM-105931),
1993.
[68] Radomsky, R., and Thole, K. A.,“Highly Turbulent Flowfield Measurement a Stator
Vane,” ASME J. of Turbomachinery, Vol. 122, pp. 255-262, 2000.
[69] Cebeci, T., and Smith. A. M. O., “ Analysis of Turbulent Boundary Layers,”
Academic Press, New York, U.S.A, 1974.
[70] Rodi, W.,“Experience with Two-Layer Models Combining the k-ε model with a
One-Equation Model near the Wall,” AIAA paper 91-0216, 1991.
[71] Launder, B. E., and Spalding, D. B.,“The Numerical Computation of Turbulent
Flows,” Comp. Meth. Appl. Mech. Eng., pp. 269-289, 1974.
[72] Yakhot, V., Orszag, S. A., Tangham, S., Gatski, T. B., and Speciale C. G.
“Development of Turbulence Models for Shear Flows by a Double Expansion
Technique,” Physics of Fluids A, vol. 4, pp. 1510-1520, 1992.
[73] Forster, D., Nelson, D. R., and Stephen M. J.“Large Distance and Long Time
Properties of a Randomly Stirred Field,” Phys. Rev. A, vol. 16, pp. 732-749, 1977.
[74] Wilcox, D. C.“Multiscale model for turbulent flows,” AIAA 24th Aerospace
Sciences Meeting American Institute of Aeronautics and Astronautics, 1986.
[75] Menter, F. R., “Influence of free stream values on k-ω turbulence model predictions,”
AIAA Journal, Vol. 30, No. 6, pp 1657-1659, 1992.
[76] Johnson, D. A., and King, L. S., “Mathematically Simple Turbulence Closure Model
for Attached and Separated Turbulent Boundary Layers,” AIAA Journal, Vol. 23, No.
11, pp 1684-1692, 1985.
[77] Bardina, J. E., Huang, P. G., and Coakley, T.,“Turbulence Modeling Validation,”
AIAA. Paper 97-2121, 1997.
92
[78] Rodi, W.,“A New Algebraic Relation for Calculation for the Reynolds Stresses,” Z.
Ange. Math. Mech., Vol.56, pp. 219-221, 1976.
[79] Speciale, C. G., and Gatski, T. B.,“Explicit Algebraic Stress Models for Complex
Turbulent Flows,” Journal of Fluid Mechanics, Vol.254, pp. 59-78, 1993.
[80] Kato, M., and Launder, B. E.,“The Modeling of Turbulent Flow Around Stationary
and Vibrating Square Cylinders,” Proc. 9th Symposium on Turbulent Shear Flows,
Kyoto, 1993.
[81] Baliga, B. R., and Patankar, S. V.,“A Control Volume Finite-Element Method for
Two-Dimensional Fluid Flow and Heat Transfer,” Numerical Heat Transfer, Vol. 6, pp.
245-261, 1983.
[82] Raw, M. R., Galpin, P. E., and Hutchinson, B. R.,“A Co-located Finite Volume
Method for Solving the Naiver-Stokes Equations for Incompressible and
Compressible Flows in Turbomachinery: Results and Application,” Canadian
Aeronautics and Space Journal, Vol. 35, No. 4, December, 1989.
[83] Van Doormaal, J. P., Turan, A. and Raithby, G. D.,“Evaluation of New Techniques
for the Calculation of Internal Recirculating Flows,” AIAA paper 87-0059, 1987.
[84] Raithby, G. D., “Skewed Upstream Differencing Schemes for Problems Involving
Fluid Flow,” Comp. Meth. Appl. Mech. Eng., Vol.9, No. 2, pp. 153-164, 1976.
[85] Patankar, S. V., “Numerical Heat Transfer,” Hemisphere, Washington, D.C., 1980.
[86] Schneider, G. E. and Raw, M. J., “Control Volume Finite Element Method for Heat
Transfer and Fluid Flow Using Colocated Variables Computational Procedure,”
Numerical Heat Transfer, Vol.11, pp. 363-390, 1987.
[87] Rhie , C. M. and Chow , W. L., “A Numerical Study of the Turbulent Flow Past an
Isolated Airfoil with Trailing Edge Separation,” AIAA paper 82-0998, 1982.
[88] Prakash , C. and Patankar , S. V., “A Control Volume Based Finite Element Method
for Solving the N-S Equations Using Equal Order Velocity-Pressure Interpolation,”
Numerical Heat Transfer, Vol.8, pp. 259-280, 1985.
[89] Ferziger, J. H. and Peric, M., “Computational Methods for Fluid Dynamics,” Springer,
Germany, 2002.
[90] Sorensen, R. L., “A Computer Program to Generate Two-Dimensional Grids About
93
Airfoils and Other Arbitrary Shapes by the Use of Poisson’s Equation,” NASA
TM-81198, 1980.
[91] Beach, T. A., “An Interactive Grid Generation Procedure for Axial and Radial Flow
Turbomachinery,” AIAA Paper 90-0344. Reno. NV. Jan. 1990.
[92] Lee, L. C., Chen, S. H., Su, S. P., and Hsu, W. Y., “An Interactive Grid Generation
Method for Turbomachinery Applications,” The 3rd Computational Fluid Dynamics
Conference, Taiwan, August 1995.
[93] Grotjans, H., and Menter, F., “Wall Functions for Industrial Applications,”
Computational Fluid Dynamics '98, Proceedings of the Fourth ECCOMAS
Conference, Papailiou K. D. et al., eds., John Wiley and Sons Ltd., Chichester, Vol. 1,
Part 2, pp.1112-1117, 1998.