摘要
为建立一种可以准确高效预测通气超空化流动的湍流模型,结合数值计算和实验结果,对绕锥头回转体通气超空化流动特性进行研究.实验采用高速录像观察了通气空化随时间的流场变化;数值计算中,分别应用标准k-ε湍流模型和密度分域的湍流模型计算了绕锥头回转体通气空化流动.其中,密度分域的湍流模型是在实验观测的基础上建立,即在空泡的前端含气量较大的区域应用DCM模型,以体现附着型空穴的可压缩性;在空泡尾部含气量较大的雾状空泡区域应用FBM湍流模型,以捕捉多尺度的空泡涡团结构.研究结果表明:与标准k-ε湍流模型相比,基于密度分域的湍流模型计算的结果与实验观测的现象基本吻合,有效减小了通气空化空泡区域内的湍流黏性,可以捕捉空泡区域内多尺度旋涡结构的演化过程,进而可以准确地预测通气超空化空泡断裂脱落的非定常流动细节.
To build a turbulence model which can accurately and effectively simulate unsteady ventilated supercavity,the characteristics of ventilated supercavity flows around a axisymmetric body were studied based on both of experimental and numerical methods.A high speed camera was used to visualize the flow structures.A standard k-εturbulence model and a modified density based turbulence model were used to simulate ventilated supercavity.The modified density based turbulence model was built based on the experimental phenomenon.In the experiment,a DCM turbulence model was applied to reflect compressibility of the attaching cavity in the front of ventilated cavity where more gas contained.And a FBM turbulence model was applied to catch multi-vortex in the rear of ventilated cavity where more fogging gas existed.The results show that compared with k-εturbulence model,the modified density based turbulence model can get better results consisting with experimental results.So the modified density based turbulence model can easily catch the evolution process of complex vortex structures in the ventilated cavity because of reducing turbulence viscosity in the domain of ventilated cavity.
引文
[1]Gopalan S,Katz J.Flow structure and modeling issues in the closure region of attached cavitation[J].Physics of Fluids,2000,12(4):3414-3431.
[2]Breuer M,Jovicic N,Mazaev K.Comparison of DES,RANS and LES forth separated flow around a flat plate at high incidence[J].International Journal Numerical Methods Fluids,2003,41:357-388.
[3]Koutmos P,Mavridis C.A computational investigation of unsteady separated flows[J].International Journal Heat and Fluid Flow,1997,18:297-306.
[4]Mavridis C,Bakrozis A,Koutmos P,et al.Isothermal and nonpremixed turbulent reacting wake flows past a two-dimensional square cylinder[J].Experimental Thermal and Fluid Sci,1998,17:90-99.
[5]Wu Jiongyang.Filter-based modeling of unsteady turbulent cavitating flow computations[D].Gainesville,Florida,USA:University of Florida,2005.
[6]余志毅.FBM模型在栅中翼型空化流动计算中的应用[J].工程热物理学报,2010,31(5):777-780.Yu Zhiyi.Application of the FBM model for the computations of cavitation flow around a cascade hydrofoil[J].Journal of Engineering Thermophysics,2010,31(5):777-780.(in Chinese)
[7]张博.亚空化非定常流动机理及动力特性研究[D].北京:北京理工大学,2009.Zhang Bo.The Hydrodynamic investigation of unsteady cloud cavitating flows[D].Beijing:Beijing Institute of Technology,2009.(in Chinese)
[8]黄彪.非定常空化数值计算模型及流动机理研究[D].北京:北京理工大学,2012.Huang Biao.Physical and numerical investigation of unsteady cavitating flows[D].Beijing:Beijing Institute of Technology,2012.(in Chinese)
[9]黄彪,王国玉,王复峰,等.非定常空化流场结构的实验研究[J].实验力学,2011,26(4):417-424.Huang Biao,Wang Guoyu,Wang Fufeng,et al.Experimental investigation on flow filed structure of unsteady cavitation flow[J].Journal of Experimental Mechanics,2011,26(4):417-424.(in Chinese)
[10]Launder B E,Spalding D B.The numerical computation of turbulent flows[J].Comp Meth Appl Mech Eng,1974,3:269-289.
[11]Coutier-Delgosha O.Numerical prediction of cavitation flow on a two-dimensional symmetrical hydrofoil and comparison to experiments[J].Journal of Fluids Engineering,2007,129(3):279-291.