半封闭圆管冲击射流湍流换热数值模拟

采用标准κ-ε模型、RNGκ—ε模型结合壁面函数以及低雷诺数κ—ε模型，对半封闭圆管冲击射流流场的平均速度、湍动能分布和Nu数分布进行了数值计算，并将此3种模型的计算结果与文献中的测量结果进行了比较．结果表明：3种湍流模型都未能完全准确地预测冲击射流场的流动特性与传热特性，其中标准κ—ε模型和低雷诺数κ—ε模型的结果很差，而RNGκ—ε模型的结果与其它两种模型相比更接近实验值。     

基于RNG k-ε湍流模型钝体绕流的数值模拟

采用RNG k-ε湍流模型对钝体绕流流场进行了数值模拟．计算采用非交错网格的有限体积法求解二维不可压N-S方程．计算结果表明，RNG k-ε湍流模型可以成功地模拟绕钝体的不稳定、非定常和剧烈分离流动。     

采用不同湍流模型对半封闭狭缝湍流冲击射流的数值模拟

针对海洋石油工程中废弃平台切割技术的需要,采用充分发展的入口条件结合标准壁面函数,将标准k-ε模型、RNG k-ε模型、Realizable k-ε模型和雷诺应力模型应用于半封闭狭缝湍流冲击射流的数值模拟,采用减小网格尺寸和提高差分格式精度等措施来减小数值伪扩散对冲击射流流动的影响.结果表明:雷诺应力模型和Realizable k-ε湍流模型在湍流冲击射流的数值模拟中具有一定的优越性,Realizable k-ε模型的计算量较小;减小网格尺寸并不能降低数值黏性对计算结果的影响;提高差分格式的精度可以在一定程度上降低数值黏性对计算结果的影响,MUSCL格式并没有显著提高计算结果的精度.     

倾角对湍流狭缝冲击射流流动与传热性能的影响

为了研究倾斜冲击射流的流场结构和传热性能,采用SIMPLE算法和RNGk-ε湍流模型,对三种入射角度下的狭缝湍流冲击射流进行了数值模拟。结果表明,由于射流的卷吸以及上封闭板的作用,在流场两侧的流出区域分别形成了一个环形的回流区。当射流倾角减小时,左侧回流区的范围变小,而右侧的回流区增大。流场中的低温区以及冲击表面压力系数最大值和局部Nu最大值的位置也都随着倾角的减小而向流场的左侧偏移。当倾角减小时,冲击表面压力系数最大值降低,局部Nu最大值没有明显的变化,射流与冲击面之间的平均传热率降低。     

半封闭圆管湍流射流冲击平板的数值研究

采用标准κ-ε模型和壁面函数法,对冲击射流场的平均速度和湍动能分布进行了数值计算,将充分发展和均匀速度分布2种入口条件下的计算结果与文献中激光多普勒测速仪的测量结果进行了比较,结果表明：数值计算能够预测到湍动能分布的趋势,但对滞止点附近湍动能的预测值过高,高达实验值的5倍;在各向同性占优势的壁射流中,平均速度的分布与实验结果具有较好的一致性;射流入口条件对湍动能分布具有较大的影响;要提高冲击射流场总体的预测效果,需要针对滞止区和近壁区内的流动特性,对κ-ε模型进行修正,同时还应采用与实验尽可能一致的入口条件。     

机舱通风口冲击射流湍流特征分析及验证

以机舱个性化通风口为背景,对自由射流加壁面冲击射流,采用几种湍流模型数值分析了10倍冲击高度、雷诺数为9,000的气流遇到冲击板后湍流脉动、流动转折等流动特征.探讨了湍流模型对冲击射流关键位置平均速度、脉动速度以及雷诺应力等流动参数特征的模拟能力,通过与实验测量结果的对比,确定了湍流模型在冲击射流中的适用性和局限性.综合来看,Relizable k-ε模型可以得到与实验值较好的符合度,但在近壁面区域会略有偏差.分析结果认为将RANS和LES组合将是开展实际机舱内通风口射流研究的最佳策略.     

基于CFD的单股冲击射流的数值模拟

运用PHOENICS计算软件,选取标准双方程模型及优选的紊流模型系数,采用有限体积法对轴对称紊流射流流场进行数值计算。计算结果表明,射流轴心速度是逐渐衰减的。选取平面紊动冲击射流问题为研究对象,通过对冲击射流内各区域流动特性的分析,给出了圆形冲击射流冲击区的范围,从喷嘴出口条件入手导出了壁面射流区中壁面切应力和径向速度的沿程分布关系式,该式反映了壁面射流场对初始条件的依赖。     

用非线性κ—ε模型对狭缝冲击射流进行数值计算

采用Ｓｐｅｚｉａｌｅ的非线性κ－ε模型和标准的κ－ε的模型计算了狭缝冲击射流,给出了径向平均速度和脉动速度沿壁面法向的工同实验结果作了对比。结果表明,当其他条件都相同时,线性和非线性模型在冲击射流的计算中相差不大,今后的改进应当在近壁区域对紊流模型本身进行。     

可压缩轴对称冲击射流流场的数值模拟

数值模拟一种可压缩轴对称冲击射流。所构造的数值模拟方法是:直接求解柱坐标系下的二维可压缩Navier-Stokes方程的差分离散方程,其中对流项采用基于非等距网格上的五阶精度迎风紧致型差分格式,黏性项采用基于非等距网格上的六阶精度对称紧致型差分格式,时间项采用3步三阶精度Runge-Kutta方法。模拟不同雷诺数、马赫数条件下冲击射流大尺度涡结构的演化过程。结果表明:流体从喷嘴射出后卷起形成一个独立的大尺度负涡,即初生漩涡,它会在壁面处逐渐激发出一个具有正涡量的壁面二次生成涡;初生漩涡和二次生成涡互相旋转挤压,壁面二次生成涡的力量很快占优势,带动初生漩涡向流场内部发展;随马赫数的增大,初生漩涡具有更强的力量,抑制了壁面二次生成涡和其他小尺度负涡的发展;随雷诺数的增大,初生漩涡的力量有所减弱,促进了壁面二次生成涡和其他小尺度负涡的发展。     

一个湍流重正化群二阶矩封闭模式的数值模拟

通过对充分发展旋转槽道和后台阶流动两个二维定常湍流问题的计算,比较了标准k-ε模型、Gibson-Launder二阶矩和一个基于湍流重正化群理论推导的二阶矩模型模拟复杂湍流运动的能力,尤其是检验这种新型的湍流重正化群二阶矩模型的性能．计算结果表明,与标准k-ε模型相比,Reynolds应力模型能够捕捉到旋转流动中由旋转带来的湍流流场结构的改变,能更准确地预测后台阶流动中的回流区长度,对于各物理量的计算,湍流重正化群二阶矩模型和Gibson-Launder模型的精度相当．     

湍流射流的数值模拟

用Ｋ－ε和Ｋ－Ｗ湍流模型及亚格子涡模型进行湍流射流的数值模拟。方法采用４步Ｒｕｎｇｅ－Ｋｕｔｔａ方法离散时间导数项,３阶ＥＮＯ格式离散对流通量项,中心差分格式离散粘性通量项,数值求解Ｒｅｙｎｏｌｄｓ平均可压缩Ｎ－Ｓ方程。     

微小孔径旋转冲击射流换热特性的数值模拟

采用数值模拟方法模拟了旋转冲击射流的换热过程,分析了换热过程中喷射孔径、喷射间距、旋转角速度以及流场分布特性对冲击射流换热的表面传热系数与平均换热效果的影响.结果表明:相同雷诺数下,3 mm孔径射流的平均换热效果比6 mm孔径的强;大孔径射流时的平均传热系数受角速度的影响要比小孔径时大;角速度的增加使换热板上的最大换热系数减小且由驻点向外偏移;旋转使板上的换热更加均匀,表现为角速度越高,平均表面传热系数曲线越平坦.     

三维引射流动数值模拟及紊流模型选择

针对引射流动三维湍流流动的特点,应用κ-ε、S-A以及k-ω湍流模型分别对三维引射流动进行了数值模拟。比较引射流量的数值结果与实验结果表明,湍流模型的选取对数值结果具有较大影响,采用k-ω模型进行模拟可以较为准确地预测引射流量。     

湍流浸没射流反应区结构的数值计算

介绍HC1气体以湍流浸没喷入NH_3水溶液容器内而发生化学反应的反应区结构的数值计算。通过化学平衡的计算,求得以混合分数为基础的状态关系;采用局部均相流动模型(LHF)、κ-ε-g湍流模型、截尾正态分布模型和法夫雷(Favre)均化的榆运方程组,得到反应区内空隙率、温度、气体穿透距离和各组分的浓度等参数值。     

聚能射流问题的数值模拟

基于多物质流体的Euler算法，用面向对象的C 语言自行编制了M-MMIC通用多物质二维流体弹塑性程序，对锥形聚能装药形成过程进行了数值模拟，并用VISC2D可视化软件对射流的形成过程进行动画演示，计算结果符合聚能射流形成的物理现象和规律，说明该物理模型和数值算法比较合理，可用于指导聚能破甲战斗部的工程设计。     

横流中三维圆形垂直浮力射流特性数值分析

采用κ-ε双方程模型，结合混合有限分析方法对横流中三维圆形垂直热水浮力射流进行了数值模拟．计算所得的温升分布与试验值较为吻合．对射流的流动结构、分叉现象、射流轨迹线等进行了预报，得出了基于速度场和温度场的射流轨迹线不同，横断面上的两个高温区的夹角范围与空气中所得结论一致，影响射流的主要因素是射流对来流的流速比等结论，从而证明了运用该数学模型对横流中三维圆形垂直热水浮力射流进行数值研究的有效性．     

用新紊流模式研究圆孔自由紊流射流

采用笔者建议的紊汉模式研究了圆孔自由紊流射流，指出脉动粘性系数Ｖｒ、基本脉动长度ｌｍ的均方差Ｌｒ和脉动能Ｋ是描述紊乱统计特征的三个基本量，首次证明了Ｖｒ与主流速Ｕｘ在断面上的变化趋势同步，提出计算Ｖｒ的新公式，用来求解流速分布和脉动切应力分布，均符合实际，文中还研究了Ｌｒ的取值问题，建议了脉动分布的公式和它在三向的分配公式，也符合实测资料，说明了模式的适用性。