激波不稳定性效应的k-ε可压缩修正湍流模型

针对高速可压缩湍流流动, 在已有的压力膨胀项和可压缩耗散率的可压缩性修正湍流模型基础上, 引入激波不稳定效应修正, 发展了一个新的可压缩性修正k-ε湍流模型. 新模型采用抑制湍流动能和耗散率方程中湍流动能产生项的方法模化激波不稳定性效应, 压力膨胀项和可压缩耗散率的可压缩性修正采用广泛使用的Sarkar修正模型. 新模型物理意义明确, 形式简单, 可适用于超声速复杂湍流流动. 对自由流动中超声速混合层和复杂的超声速横侧射流干扰流场的多个工况进行计算分析以及与实验结果的比较, 表明本文发展的k-ε模型能抑制过大的湍流动能增长, 预测结果显著优于标准k-ε模型. 对超声速混合层流动, 新模型准确预测到了混合层增长速率随对流马赫数增加而减小的趋势, 与实验结果符合地较好. 对复杂横侧射流干扰流场中的分离流动, 激波不稳定性修正抑制激波区域过大的湍流动能增长, 计算出较宽的激波区域, 从而显著改善了对强分离流动的预测结果. 流体分离越强, 修正模型效果越明显, 即使在强分离情况下, 新模型的预测结果也与实验结果较好吻合.     

带挡板圆柱可压缩绕流的大涡模拟

采用大涡模拟方法数值研究了来流马赫数为0.75、雷诺数为2×105的带挡板圆柱可压缩绕流.由于带挡板圆柱绕流的被动控制,与圆柱相比,带挡板圆柱的阻力减小了33%,升力系数的振荡几乎被完全抑制.在此类流动中,由剪切效应引起的减阻大于可压缩效应的影响.通过分析挡板的被动控制机理,我们发现抑制剪切层的失稳可以提高柱体的背压,这与柱体减阻和升力系数振荡抑制密切相关.通过分析沿剪切层的对流马赫数和自激振荡现象,发现在可压缩效应和高频力效应共同作用下,带挡板圆柱的剪切层失稳得到了抑制.     

激波与湍流相互作用的实验研究

在超声速湍流混合层风洞中, 以NPLS技术为基础研究了激波与湍流的相互作用, 观察到了斜激波与湍流边界层的相互作用以及湍流混合层中的大涡对斜激波的影响. 所得到的NPLS实验图像不仅再现了复杂的流场结构, 而且不同图像之间还具有时间相关性. 在NPLS图像基础上探讨了激波与湍流相互作用的规律.     

受损伴流湍流氢气射流火焰的直接数值模拟

采用直接数值模拟DNS的方法对受损伴流湍流氢气射流火焰进行了数值模拟, 采用16步的氢气燃烧详细化学反应机理, 冷的高速H₂/N₂燃料射流喷入热的低速伴随流, 伴随流由贫燃氢气预混火焰燃烧产生, 温度1045 K, 氧量较低. 化学反应源项由主程序在每一时间步长内动态调用CHEMKIN库函数获得. 计算采用消息传递MPI的并行计算方法, 采用12颗CPU在并行计算机上完成. 作为与实验对比的Faver平均结果由DNS瞬态结果做长时间的统计平均后获得. 火焰中涡结构的卷起以及发展过程均能很好地被捕捉, 可以观察到同旋向涡结构之间的相互吸引和反旋向涡结构之间的相互排斥过程, 伴随射流两侧涡结构彼此复杂的吸引、合并、挤压和撕裂过程, 湍流拟序结构由最初的轴对称模式开始向非对称模式演化. 流场中5.67 ms时刻瞬态的H, OH和H₂O分布, 表征了燃料射流自点燃过程中的详细火焰结构. 计算中获取的火焰抬升高度为9d ~ 11d, 与实验结果相吻合. 计算发现由OH和H粒子表征的火焰锋面中, 在火焰锋面转角位置, 燃烧过程得到强化, 可能与火焰面的拉伸以及较长的停留时间有关. 从湍流强度的分布曲线来看, 火焰的传播应该是从两侧向中心发展的. 这里的DNS结果可以作为今后发展更准确通用湍流燃烧模型的参考.     

利用双流体模型研究柴油高压喷嘴内部的空化流动

利用双流体模型对一单孔柴油喷嘴内部复杂的空化流动进行数值模拟, 并分别考察进口压力恒定和进口压力瞬变这两种不同的边界条件所对应的空化过程及流动特征. 模型预测的结果与文献中已有的实验结果吻合较好. 通过数值分析发现, 超空化的出现引起了流场结构及出流条件的显著改变, 超空流状态下的液相湍动能分布特征和喷孔出口速度类型有利于柴油射流的一次雾化过程. 此外, 喷孔上游的压力波动对空化过程有着重大的影响, 无论是部分空化还是超空化都会随着上游压力的急剧升高或降低而表现出不稳定的特性.