超声速/高超声速风洞试验段结构形式对流场性能的影响研究

通过求解轴对称 N-S 方程,对Φ1 m 高超声速风洞马赫数3和6状态下的流场进行了模拟,计算结果与试验数据基本一致,验证了所用数值方法的可信性。在此基础上,对比研究了马赫数3和6状态下采用闭口等直圆截面和开口自由射流两种试验段结构形式的超声速/高超声速风洞在起动条件下的稳态流场性能。结果表明:采用闭口等直圆截面试验段和开口自由射流试验段的流场均匀区内速度场性能指标均满足相关标准要求;马赫数3喷管采用闭口试验段时,沿风洞轴向-300mm~900mm 截面范围内的流场均匀区直径均保持在Φ882mm 以上,均匀区面积较开口试验段增加了约31.57%;马赫数6喷管采用闭口试验段时,均匀区面积比开口试验段仅增加了约8.24%,流场品质略为提高。超声速条件下,闭口试验段的流场均匀区增加明显;但在高超声速条件下,闭口试验段的流场均匀区增加比较有限。     

超声速开式空腔流动的数值模拟

研究了不同长深比、不同雷诺数条件下,超声速开式空腔流动特性的变化。算法采用空间上5阶的WENO格式,时间上具有TVD性质的3阶Runge-Kutta方法,湍流模型选用B－L模型。数值模拟获得了不同的时刻的注场图像,压力波、涡以及空腔结构之间的相互耦合作用,使得空腔内存在典型的周期性压力脉动,但不同长深比空腔的流动性存在较大差别。     

四喷管射流流场的数值模拟及分析

本文分析了ENO格式的特点，并应用于全NS方程的迁移项和压力项，模拟了三维射流，首先计算了四喷管射流干扰流场，给出了相应的计算结果，指出了其四股射流外部有形成类似一股射流波系的趋势，而内部四股射流始终存在并相互作用，且伴随流对内部干扰流场有重要影响；然后计算了四股射流流入圆管内的流场，给出了壁面和截面物理量的分布，分析了其拓扑结构。     

吸气-振荡射流激励器振荡特性

高效的振荡射流系统需要振荡特性优良的激励器,这要求减少激励器内流损失以提高出流流速、减少射流在0°偏角处的停滞时间以利于在流场中产生非定常旋涡、增大射流偏角以增大控制范围、有效调节振荡频率以接近最优控制频率。现以能大幅减少引入高压气源气体流量的吸气-振荡射流激励器为研究对象,通过数值模拟研究了不同几何外形激励器的起振、出流和频率特性。结果表明：只有当扣除射流宽度后的喉道高度大于反馈通道宽度的1.2倍,且反馈段长度足够、扩张角大小合适时,射流才会稳定振荡,并能与扩张段壁面相切;截短扩张段可使出口处射流中心速度提高67.3%;减小扩张段内分离涡的长宽比最大可使射流扫掠角达到±110°;改变反馈通道的宽度和长度会通过改变通道内通流面积和沿程损失以改变流量,从而影响频率。     

各向异性网格自适应计算在超燃模拟中的应用

采用各向异性网格自适应求解技术,将其应用于DLR超燃冲压发动机燃烧室中的超燃模拟.开展了3个算例,包括采用滑移壁面条件的燃烧室冷流场模拟、采用无滑移壁面条件的冷流场模拟及采用无滑移壁面条件的反应流模拟.模拟中,各向异性网格自适应计算捕捉到了如激波、射流、边界层、火焰面等具有各向异性特征的大梯度区域,并利用各向异性网格进行了很好的加密.对比利用各向同性网格的初始流场计算,各向异性网格自适应计算使基于滑移条件、无滑移条件的冷流计算及反应流计算的网格单元数量分别下降了36.2%,36.4%和36.8%,有效降低了计算规模,而且流场大梯度区域的计算结果更准确,辨析度更好.结果表明:对于像超燃这类具有各向异性特征的问题,各向异性网格系统比各向同性网格系统有更好的计算效率及准确性,同时也表明基于Mach数场构造的各向异性网格系统可以有效应用于超燃计算.      

脉冲射流强化圆柱喷流混合的三维数值模拟

对利用脉冲射流强化圆柱喷流混合的流场进行了数值模拟,这项工作将高频高强度信号射流射入亚音剪切层以强化喷流混合.利用有限体积法和重整化群(RNG)k-ε湍流模型求解N-S方程.分别进行了低亚音马赫数及高亚音马赫数喷流直接排入空气的数值模拟.主要计算3个状态:未受激励的状态,斯德鲁哈尔数St分别为0.2和0.4时的状态.数值模拟结果显示,脉冲激励引起喷流在与激励源平行的平面迅速扩展,在与之垂直的平面收缩,当脉冲射流激励的St在0.2附近时,相关的混合较好,而且少量的射流流量可以极大地增加速度的衰减率.      

燃油射流横流穿透深度试验和数值模拟

采用激光粒子成像测速仪(PIV),试验研究了不同动量比下燃油射流在横向空气流中的穿透深度特性。同时,采用数值计算方法,对射流穿透特性进行了模拟,并将试验和计算结果,分别与已有的经验关系式,及考虑气动韦伯数影响修改后的关系式进行对比。结果表明:燃油射流上边界深度与幂指数关系式较吻合;采用VOF两相流模型能较准确地模拟出燃油射流的喷雾核心深度;修改后的关系式与不同燃油-空气动量比范围下的喷雾核心深度较吻合;燃油-空气动量比和气动韦伯数,是影响燃油射流横流穿透深度的主要参数。     

超音通流风扇叶栅流场数值模拟

用矢通量分裂法求解Ｎ－Ｓ方程和Ｅｕｌｅｒ方程，对超音通流风扇叶栅流场进行了数值模拟，由计算结果给出的物理图象表明，叶栅通道内没有强激波，亦未出现分离回流，全叶栅流场均为超音速，说明计算的两种叶型符合设计要求，其性能有待实验检验。     

横向气流中非牛顿液体射流直接数值模拟

采用直接数值模拟研究了动量比为6的非牛顿液体射流在横向气流中的破碎特征,重点分析了表面波发展、弯曲特性和展向扩散等射流结构及其非牛顿特征。在非牛顿液体射流的近喷嘴位置,其表面迹线随时间摆动,并有逆横向气流方向运动的趋势。射流的展向扩散仅在喷入初期较为明显,此后基本维持在接近35°的扩散半角。射流的一次破碎,尤其是液丝与液滴、甚而卫星液滴的生成过程都被精细地捕捉和描述。在近喷嘴的射流柱附近,横向气流流动具有一定的圆柱扰流相似性,在其他区域则表现为极其复杂的紊流特征。射流的非牛顿特性主要体现在黏性系数,不同位置的液体黏性系数相差超过20%。相较牛顿液体射流,剪切稀化的非牛顿射流具有更易破碎的特征。     

二元双喉道射流推力矢量喷管流动参数影响的数值研究

采用数值模拟方法研究了不同流动参数对二元双喉道射流推力矢量喷管(Dual-throat fluidic Thrust-vectoring Nozzle,DTN)内流特性和推力矢量控制效果的影响。结果表明,DTN在非推力矢量情况下,NPR在3~4范围时,推力系数较大,达到0.968,而流量系数较小,仅为0.93;NPR再增大,推力系数迅速下降。在推力矢量情况下,落压比一定时,随着次流流量比的增加,推力矢量角增加,而流量系数、推力系数、推力矢量效率减小;次流流量比一定时,随着落压比的增加,推力矢量角减小,系统推力系数先增加后减小,流量系数略微增加。