近壁湍流边界层壁面局部扰动诱导形成相干结构的动力学模型

在壁面局部脉冲扰动作用下,采用直接数值模拟的方法,数值研究了湍流边界层近壁区两个相干结构形成的理论机制和非线性作用问题.然后,通过数值计算的手段,进一步探讨了近壁湍流边界层流中两个相干结构之间的非线性作用机理,分别确定它们在流向、展向方向上排列的不同结构类型分布的两个相干结构之间相互作用的效果有何区别以及它们之间间距的大小对相干结构的产生和发展有什么影响,从而可能合理地解释平面剪切湍流中低速条纹结构出现的原因.     

颗粒粒径和气载比改变时湍流调制的变化

采用相多普勒颗粒分析仪(PDPA)测试了两相圆湍射流速度场, 研究了颗粒雷诺数较小情况下, 粒径和气载比变化时颗粒对气固两相圆湍射流的湍流调制作用.实验结果表明:颗粒粒径的减小或气载比的增加均使颗粒对气相湍流度的削弱作用增强.通过无量纲化的分析, 表明颗粒粒径的减小对气相湍流度的削弱程度强于气载比增加时对湍流度的削弱程度.      

Numerical simulation of axial liquid film cooling in rocket combustor

Numerical simulation has been done for liquid film cooling in liquid rocket combustor.Multiple species of axial Navier-Stokes equations have been solved for liquid-film/hot-gas flow field,and k-ε equations have been used for compressible turbulent flow.The results of the model agree well with the results of software FLUENT.The results show that:(1) Liquid film can decrease the wall heat flux and temperature effectively,and the cold border area formed by the film covers the whole combustor and nozzle wall.(2) The turbulent viscosity is higher than the physical viscosity, and its biggest value is in the border area of the convergent area in nozzle.The effect of turbulent flow on the whole simulation field can not be ignored.(3) The mass fraction of kerosene at the film inlet is 1,but it decreases along the nozzle wall and achieves its lowest value at the outlet.However,the mass fraction of kerosene near the wall is the biggest at any axial location.      

火箭发动机燃烧室液膜-再生复合冷却数值仿真

对液体火箭发动机燃烧室液膜-再生复合冷却进行了数值计算,针对液膜-燃气流场区多组分、轴对称Navier-Stokes(N-S)方程和再生冷却区单组分N-S方程进行求解,并使用k-ε方程求解湍流流动.对文献中的某液氧/煤油火箭发动机燃烧室进行了数值模拟,该模型的计算结果能够与文献中的计算结果较好地吻合.计算结果表明:①液膜-再生复合冷却能有效地减少壁面热流密度和降低壁面温度,且其形成的冷气边区覆盖了整个燃烧室及喷管壁面;②再生冷却液入口质量流量越大,复合冷却作用越明显,壁面温度越低;③随再生冷却液质量流量的不同其温升在450~600K之间,且质量流量越大,再生冷却液的温升越小.④壁面煤油的质量分数不断下降,在喷管出口壁面处达到最低值,但含有煤油的区域不断变大.      

基于Delaunay背景网格插值和局部网格重构的变形体动态混合网格生成技术

建立了一种基于Delaunay背景网格插值方法和局部网格重构方法相结合的变形体动态混合网格生成方法.首先利用作者发展的定态混合网格生成技术生成变形体初始网格,即变形体附近利用结构化的四边形网格,外场采用自适应Cartesian网格,中间由三角形网格过渡.当物体运动或变形时, 首先利用Liu和Qin提出的基于Delaunay背景网格的插值方法对变形体附近的网格进行变形;当物体运动或变形位移很大导致局部网格质量急剧下降、甚至网格相交时,则在局部重新生成网格.利用该方法生成了多种变形体外形的动态混合网格,如单个鱼体、双鱼串列巡游,昆虫单翼、双翼扑动,多段翼型变形等,并在鱼体巡游的非定常数值模拟中得到了应用.     

利用分区ALE算法数值模拟圆柱湍流涡致振动

用分区算法结合任意拉格朗日-欧拉法(ALE)数值模拟了圆柱湍流涡致振动.求解基于非交错网格系统,压力求解采用压力泊松方程提法,湍流模型采用标准k-ε模型和重整化群(RNG)k-ε模型.计算取中等雷诺数Re=5000、10000、15000、25000、50000等,质量系数M=10.阻尼系数ζ=0.00331,自振频率fn=0.18315、0.1628.计算结果表明:湍流涡致振动下圆柱时均阻力系数大于孤立圆柱绕流,而升力系数(振幅)值都小于孤立圆柱绕流.随着雷诺数增大,湍流粘性系数随之增大,但湍动能和湍流耗散率变化趋势不明显.对孤立圆柱绕流,研究结果与前人实验结果基本相符.     

用电子束荧光技术研究高超声速边界层中的扰动现象

利用电子束荧光技术对高超声速平板边界层中的扰动现象进行实验研究,实验在炮风洞中进行,自由流马赫数为7．8,单位长度雷诺数为3．5×107/m,测量了平均密度分布、脉动密度分布,并获得了湍流密度脉动的互相关分布和频谱分布。实验结果表明在过渡区中有序与随机现象并存。     

圆转方非轴对称喷管无粘流动的多层网格法求解

本文利用多层网格法求解三维无粘亚音和跨音速圆转方非轴对称喷管的流动,基本格式是把Ni[1]的二维有限体积积分格式推广到三维流场。边界条件采用文[2]中的“预测一修正”处理方法。为了进一步提高格式稳定性,在计算流动变量的一阶时间变化量时,采用单元体积平均法而不是文[1]中的算术平均法。引入当地时间步及多层网格方法,以加快格式的收敛速度。文中还对多层网格构造形式及其收敛效果作了讨论。通过对两种喷管模型流动的计算和与试验数据的对比,验证了该算法的精度。     

偏心底吹气体搅拌下冶金熔池水模型内三维湍流的实验研究

用一个圆柱形水模型模拟偏心底吹气体搅拌的冶金熔池,对模型内的三个典型剖面测量了三维速度场和湍流参数分布。方法是用二维激光多普勒测速仪在两个方向上测两次合成。面内的平行速度分布表明小流量偏心底吹熔池内的流谱为正循环流动。垂直速度分布表明,在气水两相流气泡柱附近有两个较强的圆周方向的环流,这种现象在中心底吹轴对称二维熔池流场中是不存在的。在气泡柱附近的区域内,速度模与湍动能值与较远的区域内的值相比要大得多,这就增强了这里的传热、传质和化学反应过程。在这里,气泡柱与壁面间的距离对湍流度的大小有影响,在速度值相近的情况下,较宽的区域内具有较大的湍流度值。     

二维不可压内流湍流附面层分离流的预测

本文采用适合于内流的准同时粘性/无粘性流匹配模型,快速求解二维通道内的不可压湍流附面层分离流动。无粘流采用由复变函数论的Plemeli公式推导的面元法;粘性流动区采用Head积分法求解。文中给出了两个算例。结果表明:所采用的准同时粘性/无粘性流匹配模型,不但可有效地克服分离形成的奇异性,而且收敛速度快,稳定性好。在湍流分离阶段,其预测结果与实验值基本吻合。文中还讨论了形状因子及表面摩擦系数关系式对预测结果的影响,形状因子经验关系对预测有比较明显的影响,而表面摩擦系数经验关系对预测无显著影响。文中还提出了改进模型的方向。