锥形管路系统内紊流场的数值仿真与预测研究
|
摘要
用DLR型k-ε紊流模型·BFC(边界拟合曲线坐标变换)法,对扩散角为4°,扩散度为4的锥形渐扩管路系统内具有逆压梯度的充分发展的不可压粘性紊流进行数值仿真与诊断。应用实例的入口雷诺数分别为1.16×10~5和2.93×10~5。在近壁密集型径向适体非均匀网格系统下,通过36种算例的数值实验,研究模型函数f_μ和f_2的不同组合,改变网格布局,简化ε方程和压力泊松方程及改变每时间步的迭代次数I和时间步长Δt等条件对计算结果的影响,将计算结果与实验结果相比较,绘制出流场各流动参数的分布图。其中包括时均流速(?),(?),压力(?),紊流动能k,紊流耗散率ε,涡动粘性系数v_t和壁面压力系数C_(pw)。对数值仿真的结果进行诊断分析,获得难以用实物实验或模型实验测得的相对精确的数据,为该种紊流理论的研究和工程应用提供比较可靠的参考依据。
在文章的结尾,作者对本文所做的工作进行了总结和概括,并且指出了今后的研究工作。
Fully developed incompressible viscous turbulent flow with adverse pressure gradient in a conical diffuser tube system having a divergence of 4°and an area ratio of 4 has been simulated and diagnosed by a DLR k-εturbulent model and it's BFC(Boundary-Fitted Coordinate Transformation) method.The research has been done for pipe entry Reynolds number of 1.16×10~5 and 2.93×10~5.Based on the near wall denseness model radial fittedbody nonuniform grid system,The calculated results of 36 cases have been given under the different calculated conditions,such as different model function combinations,different grid dispositions,predigest theεfunction with pressure function and different I with different△t.They were compared with the experimental results respectively.The parameter distributions of the area have been given respectively,including the distributions of the mean flow velocity(?)(?),mean pressure(?),turbulent energy k,turbulent dissipationε, eddy kinematic viscosity v_t and the coefficient of wall pressure C_(pw).The results of numerical simulation have been diagnosed in order to obtain the more accurately numerical prediction results,which haven't been gotten by using material experiment and model experiment,the results have offered relatively credible reference for turbulent theoretical research and engineering technical actual application.
At last,the researcher has been summarized this article and the work in the future have been put forward.
在文章的结尾,作者对本文所做的工作进行了总结和概括,并且指出了今后的研究工作。
Fully developed incompressible viscous turbulent flow with adverse pressure gradient in a conical diffuser tube system having a divergence of 4°and an area ratio of 4 has been simulated and diagnosed by a DLR k-εturbulent model and it's BFC(Boundary-Fitted Coordinate Transformation) method.The research has been done for pipe entry Reynolds number of 1.16×10~5 and 2.93×10~5.Based on the near wall denseness model radial fittedbody nonuniform grid system,The calculated results of 36 cases have been given under the different calculated conditions,such as different model function combinations,different grid dispositions,predigest theεfunction with pressure function and different I with different△t.They were compared with the experimental results respectively.The parameter distributions of the area have been given respectively,including the distributions of the mean flow velocity(?)(?),mean pressure(?),turbulent energy k,turbulent dissipationε, eddy kinematic viscosity v_t and the coefficient of wall pressure C_(pw).The results of numerical simulation have been diagnosed in order to obtain the more accurately numerical prediction results,which haven't been gotten by using material experiment and model experiment,the results have offered relatively credible reference for turbulent theoretical research and engineering technical actual application.
At last,the researcher has been summarized this article and the work in the future have been put forward.
引文
[1]孔珑(主编).流体力学[M].北京:高等教育出版社,2003:1-5.
[2]陈懋章(编著).粘性流体动力学理论及紊流工程计算[M].北京:北京航空航天大学出版社,1986:1-3.
[3]陈懋章(编著).粘性流体动力学基础[M].北京:高等教育出版社,2002:1-4,46-50,438-446.
[4]付德薰,马延文.计算流体力学[M].北京:高等教育出版社,2002:1-4,191-197.
[5]张兆顺,崔桂香,许春晓.湍流理论与模拟[M].北京:清华大学出版社,2005:1,209-213.
[6]陶文铨.数值传热学[M].西安:西安交通大学出版社,2001:332-335,362-370.
[7]赵新雅.紊流数学模型的发展[J].太原:山西建筑,2005,31(10):133-134.
[8]阎超,钱翼稷,连祺祥(编著).粘性流体力学[M].北京:北京航空航天大学出版社,2005:130-132,164-170.
[9]Lander B E,Spalding D B.The Numerical Computation of Turbulent Flows[M].New York:Northholland Publishing Company,1974:269.
[10]华绍增,杨学宁等(编译).实用流体阻力手册[M].北京:国防工业出版社,1985:178-198.
[11]何永森,刘邵英.机械管内流体数值预测[M].北京:国防工业出版社,1999:19-21,201-230,251.
[12]何永森,小林敏雄,森西洋平.用近壁低Re数型k-ε模型对锥形扩散器内紊流的数值预测[C].第四届亚洲流体机械国际会议论文集,1993:417-424.
[13]何永森,蒋光彪.锥形渐扩管内紊流数值预测诊断系统的研究-模型函数、网格配置和雷诺数的影响[J].湘潭:湘潭大学自然科学学报,2006,(2):28.
[14]A.N.Rousseau,L.D.Albright,and K.E.Torrance.A Short Comparison of Damping Functions of Standard Low-Reynolds-Number Models[J].JUNE 1997,Vol.119,460.
[15]张涵信,沈孟育.计算流体力学—差分方法的原理和应用[M].北京:国防工业出版社,2003:347-380.
[16]任安禄.不可压缩粘性流场计算方法[M].北京:国防工业出版社,2003:23-33.
[17]J.F.Thompson,F.C.Thames,C.W.Mastin.TOMCAT-A Code for Numerical Generation of Boundary-Fitted Curvilinear Coordinate Systems on Fields Containing Any Number of Arbitrary Two-Dimensional Bodies[J].Computational Physics,1977,24:274-302.
[18]王德意,魏文礼,刘玉玲,等.正交曲线网格生成技术研究[J].西安理工大学学报,2000,2(16):188-191.
[19]Wei Wenli.A new method to determine the source terms of poisson equations for grid generation[J].Hydrodynamics,2000,12(1):16-21.
[20]刘哲,魏文礼,郭永涛,等.边界拟合坐标网格生成方法研究[J].西安理工大学学报,2004,4(20):413-415.
[21]Thompson J F.Numerical grid generation.Amsterdam,Elsevier Science Publishing Co.1982:1-30.
[22]J.L.Steger,R.L.Sorenson.Automatic Mesh-Point Clustering Near a Boundary in Grid Generation with Elliptic partial Differential Equations[J].Computational Physics,1979,33:405-410.
[23]刘玉玲,魏文礼,沈永明,等.边界处正交的曲线网格生成技术合理调节因子选取的研究[J].计算力学学报,2000,2(19):195-197.
[24]张涵信.网格和高精度差分计算问题[C].全国第九届计算流体力学会议,1989:1-9.
[25]张国强,吴家鸣.流体力学[M].北京:机械工业出版社,2006:235-242.
[26]P.J.罗奇(著),钟锡昌,刘学宗(译).计算流体动力学[M].北京:科学出版社,1983:257-275.
[27]林建中,阮晓东,陈邦国等(编著).流体力学[M].北京:清华大学出版社,2005:508-516.
[28]孙勇,何永森.锥形渐扩管内紊流数值预测诊断系统的研究-网格数和差分方法的影响[J].湘潭:湘潭大学自然科学学报,2005,(1):47.
[29]J.Laufer.The Structure of Turbulence in Fully Developed Pipe Flow[J].NACA Rep (1174),Washington:U.S.G.P.O,1955:6-15.
[30]徐士良FORTRAN常用算法程序集[M].北京:清华大学出版社,1992:26-39.
[31]何光渝,高水利.Visual Fortran常用数值算法集[M].北京:科学出版社,2003:62-65.
[32]Virendra C.Patel,Wolfgang Rodi,Georg Scheuerer.Turbulence Models for Near-Wall and Low Reynolds Number Flows:A Review[J].AIAA Journal,1984,23(9):1308-1319.
[33]S.J.Wang,A.S.Mujumdar.A comparative study of five low Reynolds number k-ε models for impingement heat transfer[J].Applied Thermal Engineering,2005,25:31-44.
[34]C.M.Hrenya,E.J.Bolio,D.Chakrabarti,etal.Comparison of low Reynolds number k-εturbulence models in predicting fully developed pipe flow[J].Chemical Engineering Science,1995,50(12):1923-1941.
[35]P.A.C.Okwuobi,R.S.Azad.Turbulence in a Conical Diffuser with Fully Developed Flow at Entry[J].Fluid Mech,1973,57(3):603-622.
[36]D.Singh,R.S.Azad.Turbulent Kinetic Energy Balance in a Conical Diffuser[J].Proc.of Turbulence,1981:21-33.
[37]He Yongsen,Kobayashi T,Morinishi Y.Numerical Prediction of Turbulent Flow in a Conical Diffuser Using k-ε Model[J].Acta Meshanica Sinica,1992,8(2):117-126.
[38]He Yongsen,Kobayashi T,Morinishi Y.Numerical Prediction of Turbulent Flow in a Conical Diffuser Using Model for Near wall and Low Re Number[J].Suzhou,CMES,1993:593-599.
[39]Myong,H.K.and Kasagi,N.A New Approach to the Improvement of k-εTurbulence Model for wall-Bunded shear Flows[J].ASME,Int.J.SerⅡ,1990(33):63.
[40]Georgi Kalitzin,Gorezd Medic,Gianluca Laccarino,etal.Near-Wall behavior of RANS turbulence models and implications for wall functions[J].Journal of Computational Physics,2004:1-27.
[41]Nagano Y,Hishida M.Improved From of the k-ε Model for Wall Turbulent Shear Flows[J].TRrans ASME,Fluids Eng,1987:109-156.
[42]卢小勇,骆坚,何永森,等.锥形渐扩管内紊流计算机仿真诊断系统的研究[J].矿业研究与开发,2001,21(1):25-28.
[2]陈懋章(编著).粘性流体动力学理论及紊流工程计算[M].北京:北京航空航天大学出版社,1986:1-3.
[3]陈懋章(编著).粘性流体动力学基础[M].北京:高等教育出版社,2002:1-4,46-50,438-446.
[4]付德薰,马延文.计算流体力学[M].北京:高等教育出版社,2002:1-4,191-197.
[5]张兆顺,崔桂香,许春晓.湍流理论与模拟[M].北京:清华大学出版社,2005:1,209-213.
[6]陶文铨.数值传热学[M].西安:西安交通大学出版社,2001:332-335,362-370.
[7]赵新雅.紊流数学模型的发展[J].太原:山西建筑,2005,31(10):133-134.
[8]阎超,钱翼稷,连祺祥(编著).粘性流体力学[M].北京:北京航空航天大学出版社,2005:130-132,164-170.
[9]Lander B E,Spalding D B.The Numerical Computation of Turbulent Flows[M].New York:Northholland Publishing Company,1974:269.
[10]华绍增,杨学宁等(编译).实用流体阻力手册[M].北京:国防工业出版社,1985:178-198.
[11]何永森,刘邵英.机械管内流体数值预测[M].北京:国防工业出版社,1999:19-21,201-230,251.
[12]何永森,小林敏雄,森西洋平.用近壁低Re数型k-ε模型对锥形扩散器内紊流的数值预测[C].第四届亚洲流体机械国际会议论文集,1993:417-424.
[13]何永森,蒋光彪.锥形渐扩管内紊流数值预测诊断系统的研究-模型函数、网格配置和雷诺数的影响[J].湘潭:湘潭大学自然科学学报,2006,(2):28.
[14]A.N.Rousseau,L.D.Albright,and K.E.Torrance.A Short Comparison of Damping Functions of Standard Low-Reynolds-Number Models[J].JUNE 1997,Vol.119,460.
[15]张涵信,沈孟育.计算流体力学—差分方法的原理和应用[M].北京:国防工业出版社,2003:347-380.
[16]任安禄.不可压缩粘性流场计算方法[M].北京:国防工业出版社,2003:23-33.
[17]J.F.Thompson,F.C.Thames,C.W.Mastin.TOMCAT-A Code for Numerical Generation of Boundary-Fitted Curvilinear Coordinate Systems on Fields Containing Any Number of Arbitrary Two-Dimensional Bodies[J].Computational Physics,1977,24:274-302.
[18]王德意,魏文礼,刘玉玲,等.正交曲线网格生成技术研究[J].西安理工大学学报,2000,2(16):188-191.
[19]Wei Wenli.A new method to determine the source terms of poisson equations for grid generation[J].Hydrodynamics,2000,12(1):16-21.
[20]刘哲,魏文礼,郭永涛,等.边界拟合坐标网格生成方法研究[J].西安理工大学学报,2004,4(20):413-415.
[21]Thompson J F.Numerical grid generation.Amsterdam,Elsevier Science Publishing Co.1982:1-30.
[22]J.L.Steger,R.L.Sorenson.Automatic Mesh-Point Clustering Near a Boundary in Grid Generation with Elliptic partial Differential Equations[J].Computational Physics,1979,33:405-410.
[23]刘玉玲,魏文礼,沈永明,等.边界处正交的曲线网格生成技术合理调节因子选取的研究[J].计算力学学报,2000,2(19):195-197.
[24]张涵信.网格和高精度差分计算问题[C].全国第九届计算流体力学会议,1989:1-9.
[25]张国强,吴家鸣.流体力学[M].北京:机械工业出版社,2006:235-242.
[26]P.J.罗奇(著),钟锡昌,刘学宗(译).计算流体动力学[M].北京:科学出版社,1983:257-275.
[27]林建中,阮晓东,陈邦国等(编著).流体力学[M].北京:清华大学出版社,2005:508-516.
[28]孙勇,何永森.锥形渐扩管内紊流数值预测诊断系统的研究-网格数和差分方法的影响[J].湘潭:湘潭大学自然科学学报,2005,(1):47.
[29]J.Laufer.The Structure of Turbulence in Fully Developed Pipe Flow[J].NACA Rep (1174),Washington:U.S.G.P.O,1955:6-15.
[30]徐士良FORTRAN常用算法程序集[M].北京:清华大学出版社,1992:26-39.
[31]何光渝,高水利.Visual Fortran常用数值算法集[M].北京:科学出版社,2003:62-65.
[32]Virendra C.Patel,Wolfgang Rodi,Georg Scheuerer.Turbulence Models for Near-Wall and Low Reynolds Number Flows:A Review[J].AIAA Journal,1984,23(9):1308-1319.
[33]S.J.Wang,A.S.Mujumdar.A comparative study of five low Reynolds number k-ε models for impingement heat transfer[J].Applied Thermal Engineering,2005,25:31-44.
[34]C.M.Hrenya,E.J.Bolio,D.Chakrabarti,etal.Comparison of low Reynolds number k-εturbulence models in predicting fully developed pipe flow[J].Chemical Engineering Science,1995,50(12):1923-1941.
[35]P.A.C.Okwuobi,R.S.Azad.Turbulence in a Conical Diffuser with Fully Developed Flow at Entry[J].Fluid Mech,1973,57(3):603-622.
[36]D.Singh,R.S.Azad.Turbulent Kinetic Energy Balance in a Conical Diffuser[J].Proc.of Turbulence,1981:21-33.
[37]He Yongsen,Kobayashi T,Morinishi Y.Numerical Prediction of Turbulent Flow in a Conical Diffuser Using k-ε Model[J].Acta Meshanica Sinica,1992,8(2):117-126.
[38]He Yongsen,Kobayashi T,Morinishi Y.Numerical Prediction of Turbulent Flow in a Conical Diffuser Using Model for Near wall and Low Re Number[J].Suzhou,CMES,1993:593-599.
[39]Myong,H.K.and Kasagi,N.A New Approach to the Improvement of k-εTurbulence Model for wall-Bunded shear Flows[J].ASME,Int.J.SerⅡ,1990(33):63.
[40]Georgi Kalitzin,Gorezd Medic,Gianluca Laccarino,etal.Near-Wall behavior of RANS turbulence models and implications for wall functions[J].Journal of Computational Physics,2004:1-27.
[41]Nagano Y,Hishida M.Improved From of the k-ε Model for Wall Turbulent Shear Flows[J].TRrans ASME,Fluids Eng,1987:109-156.
[42]卢小勇,骆坚,何永森,等.锥形渐扩管内紊流计算机仿真诊断系统的研究[J].矿业研究与开发,2001,21(1):25-28.