超声速欠膨胀冲击射流数值模拟
|
摘要
超声速欠膨胀冲击射流有着重要的实际应用价值,如S/TOVL飞行器、火箭发射、除尘等。其流场结构复杂,包含间断激波、反射激波、马赫盘、滞止泡以及冲击平面传热不同于亚声速冲击射流的特点。为了分析超声速欠膨胀冲击射流流场和传热,采用有限体积法,结合k-l湍流模型以及二阶精度的TVD格式进行数值模拟:对比实验和k-ε湍流模型的努赛尔数,得出k-l湍流模型在传热问题中更具优势;对比阴影图和计算密度云图以及对比冲击平面压力系数的实验值和计算值,验证了k-l湍流模型模拟超声速欠膨胀冲击射流流场的合理性;采用k-l湍流模型研究3种冲击高度(3D,6D,10D),3种压比(2.0,3.4,4.4),3种喷管总温(493 K,591 K,580 K)下,冲击平面温度分布。数值研究结果对分析超声速欠膨胀冲击平面的烧蚀有一定的指导意义。
Under-expanded supersonic impinging jets are widely used in many engineering fields,such as S/TOVL aircraft,the rocket launch,and the removal of dust,et al. Supersonic impinging jet also has complicated structures,reflected shock,mash disk,stagnation bubble,and the local heat transfer on an impinged surface has a different characteristic than a subsonic impinging jet. In order to analyze the flow-field and heat transfer of under-expanded supersonic impinging jet,it is that used a finite volume approach combined with k-l turbulence model and second-order TVD scheme. The main investigation:Compared to the nusselt number of experimental value and k-ε turbulence model,which showed that k-l turbulence model had an advantage in heat transfer problems. Compared to the shadowgraph and density contour,the pressure coefficient of calculated value and experimental value,which verified and validated the k-l turbulence model. The k-l turbulence model was then used to investigate the distribution of temperature on three impingement plane heights(3 D,5 D,7 D,l0 D),three nozzle pressure ratios(2.0,3.4,4.4) and three nozzle total temperatures(493 K,591 K,580 K). An understanding of these flow characteristics is essential since they contribute directly to the study of ground erosion.
Under-expanded supersonic impinging jets are widely used in many engineering fields,such as S/TOVL aircraft,the rocket launch,and the removal of dust,et al. Supersonic impinging jet also has complicated structures,reflected shock,mash disk,stagnation bubble,and the local heat transfer on an impinged surface has a different characteristic than a subsonic impinging jet. In order to analyze the flow-field and heat transfer of under-expanded supersonic impinging jet,it is that used a finite volume approach combined with k-l turbulence model and second-order TVD scheme. The main investigation:Compared to the nusselt number of experimental value and k-ε turbulence model,which showed that k-l turbulence model had an advantage in heat transfer problems. Compared to the shadowgraph and density contour,the pressure coefficient of calculated value and experimental value,which verified and validated the k-l turbulence model. The k-l turbulence model was then used to investigate the distribution of temperature on three impingement plane heights(3 D,5 D,7 D,l0 D),three nozzle pressure ratios(2.0,3.4,4.4) and three nozzle total temperatures(493 K,591 K,580 K). An understanding of these flow characteristics is essential since they contribute directly to the study of ground erosion.
引文
[1]徐惊雷,徐忠,肖敏,等.冲击射流的研究概述[J].力学与实践,1999,21(6):8-17.
[2]DONALDSON C D,SNEDEKER R S,MARGOLIS D P. Astudy of free jet impingement. part 2. free jet turbulent struc-ture and impingement heat transfer[J]. Journal of Fluid Me-chanics,1971,45(03):477-512.
[3]MACKIE S,TAGHAVI R. Supersonic impinging jets:a com-putational investigation[C]//40th AIAA Aerospace SciencesMeeting and Exhibit,2001:2002-0671.
[4]LOU H,SHIH C,ALVI F S. A PIV study of supersonic im-pinging jet[J]. AIAA,2003,3263(9):1003-1007.
[5]LIMAYE M D,VEDULA R P,PRABHU S V. Local heattransfer distribution on a flat plate impinged by a compress-ible round air jet[J]. International Journal of Thermal Sci-ences,2010,49(11):2157-2168.
[6]何枫,谢峻石.超声速尔膨胀冲击射流的数值模拟[J].推进技术,2002,23(2):96-99.
[7]刘海,姚朝晖.超声速高度欠膨胀冲击射流的大涡模拟[J].空气动力学学报,2012,30(1):80-85.
[8]姚朝晖,侯修洲,郝鹏飞.超声速冲击射流的PIV实验研[J].实验流体力学,2007,21(4):32-35.
[9]刘小军,傅德彬,牛青林,等.燃气射流冲击传热特性的数值模拟[J].爆炸与冲击,2015,35(2):229-235.
[10]ZUCKERMAN N,LIOR N. Jet impingement heat transfer:physics,correlations,and numerical modeling[J]. Advancesin Heat Transfer,2006,39:565-631.
[11] ANGIOLETTI M,NINO E,RUOCCO G. CFD turbulentmodelling of jet impingement and its validation by particleimage velocimetry and mass transfer measurements[J]. In-ternational Journal of Thermal Sciences,2005,44(4):349-356.
[12]CHUNG Y M,LUO K H. Unsteady heat transfer analysis ofan impinging jet[J]. Journal of Heat Transfer,2002,124(6):1039-1048.
[13]GOLDBERG U,CHAKRAVARTHY S. A k-l turbulenceclosure for wall-bounded flows[J]. AIAA,2005,4638:2005.
[14]刘超群.多重网格法及其在计算流体力学中的应用[M].北京:清华大学出版社,1995.
[15]BAUGHN J W,HECHANOVA A E,YAN X. An experi-mental study of entrainment effects on the heat transfer froma flat surface to a heated circular impinging jet[J]. Journalof Heat Transfer,1991,113(4):1023-1025.
[16] YAN X,BAUGHN J W,MESBAH M. The effect ofReynolds number on the heat transfer distribution from aflat plate to an impinging jet[J]. Fundamental and AppliedHeat Transfer Research for Gas Turbine Engines,1992:1-7.
[17]BEHNIA M,PARNEIX S,DURBIN P A. Prediction of heattransfer in an axisymmetric turbulent jet impinging on a flatplate[J]. International Journal of Heat and Mass Transfer,1998,41(12):1845-1855.
[18]ALVI F S,LADD J A,BOWER W W. Experimental andcomputational investigation of supersonic impinging jets[J].AIAA Journal,2002,40(4):599-609.
[19]MESSERSMITH N L,MURTHY S N B. Gas dynamics andheat transfer from impinging underexpanded jets[C]//AIAAand DGLR,International Aerospace Planes and HypersonicsTechnologies Conference,5 th,Munich,Germany,1993.
[2]DONALDSON C D,SNEDEKER R S,MARGOLIS D P. Astudy of free jet impingement. part 2. free jet turbulent struc-ture and impingement heat transfer[J]. Journal of Fluid Me-chanics,1971,45(03):477-512.
[3]MACKIE S,TAGHAVI R. Supersonic impinging jets:a com-putational investigation[C]//40th AIAA Aerospace SciencesMeeting and Exhibit,2001:2002-0671.
[4]LOU H,SHIH C,ALVI F S. A PIV study of supersonic im-pinging jet[J]. AIAA,2003,3263(9):1003-1007.
[5]LIMAYE M D,VEDULA R P,PRABHU S V. Local heattransfer distribution on a flat plate impinged by a compress-ible round air jet[J]. International Journal of Thermal Sci-ences,2010,49(11):2157-2168.
[6]何枫,谢峻石.超声速尔膨胀冲击射流的数值模拟[J].推进技术,2002,23(2):96-99.
[7]刘海,姚朝晖.超声速高度欠膨胀冲击射流的大涡模拟[J].空气动力学学报,2012,30(1):80-85.
[8]姚朝晖,侯修洲,郝鹏飞.超声速冲击射流的PIV实验研[J].实验流体力学,2007,21(4):32-35.
[9]刘小军,傅德彬,牛青林,等.燃气射流冲击传热特性的数值模拟[J].爆炸与冲击,2015,35(2):229-235.
[10]ZUCKERMAN N,LIOR N. Jet impingement heat transfer:physics,correlations,and numerical modeling[J]. Advancesin Heat Transfer,2006,39:565-631.
[11] ANGIOLETTI M,NINO E,RUOCCO G. CFD turbulentmodelling of jet impingement and its validation by particleimage velocimetry and mass transfer measurements[J]. In-ternational Journal of Thermal Sciences,2005,44(4):349-356.
[12]CHUNG Y M,LUO K H. Unsteady heat transfer analysis ofan impinging jet[J]. Journal of Heat Transfer,2002,124(6):1039-1048.
[13]GOLDBERG U,CHAKRAVARTHY S. A k-l turbulenceclosure for wall-bounded flows[J]. AIAA,2005,4638:2005.
[14]刘超群.多重网格法及其在计算流体力学中的应用[M].北京:清华大学出版社,1995.
[15]BAUGHN J W,HECHANOVA A E,YAN X. An experi-mental study of entrainment effects on the heat transfer froma flat surface to a heated circular impinging jet[J]. Journalof Heat Transfer,1991,113(4):1023-1025.
[16] YAN X,BAUGHN J W,MESBAH M. The effect ofReynolds number on the heat transfer distribution from aflat plate to an impinging jet[J]. Fundamental and AppliedHeat Transfer Research for Gas Turbine Engines,1992:1-7.
[17]BEHNIA M,PARNEIX S,DURBIN P A. Prediction of heattransfer in an axisymmetric turbulent jet impinging on a flatplate[J]. International Journal of Heat and Mass Transfer,1998,41(12):1845-1855.
[18]ALVI F S,LADD J A,BOWER W W. Experimental andcomputational investigation of supersonic impinging jets[J].AIAA Journal,2002,40(4):599-609.
[19]MESSERSMITH N L,MURTHY S N B. Gas dynamics andheat transfer from impinging underexpanded jets[C]//AIAAand DGLR,International Aerospace Planes and HypersonicsTechnologies Conference,5 th,Munich,Germany,1993.