场协同原理在高超声速化学非平衡流动中的推广

应用理论分析与数值模拟方法,将对流传热的场协同原理从不可压缩流动推广至高超声速化学非平衡流动中。结果表明,高超声速化学非平衡层流与湍流的热流密度取决于流动的当地单位体积的动量与单位质量总焓梯度的协同。用当地单位体积的动量与单位质量总焓梯度的协同研究高超声速化学非平衡流动的壁面传热问题,对层流流动下的对流传热,不但计及了高超声速化学非平衡流的密度变化对热流密度的影响,而且包括了静焓梯度、压力梯度、边界层内的分子黏性剪切效应对热流密度的作用;对湍流问题,除了上述层流流动各项对热流密度的影响外,还计及了雷诺剪切应力对热流密度的作用。考虑到高超声速化学非平衡流静焓的定义,高超声速化学非平衡层流及湍流的场协同同时计及所有组分的平动能、转动能、振动能及电子能等梯度的贡献。     

壁面温度分布的高超声速平板边界层稳定性分析及转捩预测

针对来流马赫数为4.5、6.0和7.0的高超声速平板边界层,取30km高空处的气体参数,壁面为等温、绝热和温度分布等3种不同条件,采用eN方法进行转捩预测。其中,壁面温度分布条件下,在等温壁（冷壁）和绝热壁之间,给出4种流向分布进行分析。取扰动的初始幅值为0.3‰,以幅值达到1.5%作为转捩判断依据。结果表明:当温度为来流温度时,等温壁面条件的转捩位置最靠前,并随马赫数增大更加靠前;绝热壁面条件的转捩位置随马赫数增大而推后;壁面温度分布条件下,在相同时刻,马赫数越大,转捩位置越靠前。相同马赫数下,壁面温度较高者转捩位置较靠后（马赫数为7.0时,不完全满足此规律）。在马赫数为4.5和6.0时,绝热壁面条件转捩由第一模态主导,其余情况主导转捩的都是第二模态。     

计算气动光学研究进展

计算气动光学是采用数值计算的方法研究高速空气动力流场对光波传输和光学成像影响的一门交叉学科,在红外成像制导设计等工程应用中发挥着重要的作用。总结了计算气动光学研究的三种基本计算方法,包括基于CFD简化方法的光学统计估算、基于RANS流场的计算和基于LES/DNS瞬态流场的计算,对这三种计算方法的基本理论、计算过程进行了介绍,综述了新近的一些计算气动光学研究进展,最后指出了未来的研究重点和方向。     

超声速混合层中PIV粒子的湍流变动作用研究

对二维超声速气固两相混合层进行双向耦合,研究了粒子图像测速技术（PIV）中示踪粒子对超声速混合层的湍流变动作用。超声速气固两相混合层的气相采用大涡模拟,离散相采用拉格朗日颗粒轨道模型求解。结果表明:与无负载示踪粒子时的超声速混合层相比,小Stokes数示踪粒子在超声速混合层中的布撒减弱了流向湍流,而强化了法向湍流,使雷诺应力峰值增大了9.68%;大Stokes数示踪粒子对混合层的湍流脉动起到了一定的削弱作用,最大雷诺应力值只有无负载时的41.74%。大质量载荷时,大量示踪粒子的运动尾迹抹平了部分法向速度脉动,使最大法向速度脉动只有无负载粒子时的38.63%;中等质量载荷时,超声速混合层的法向速度脉动和雷诺应力峰值与无负载粒子时相近;而小质量载荷时,超声速混合层中心线及其附近的法向速度脉动得到较小的增强,而最大流向速度脉动却被削弱了19.29%。小Stokes数和中等质量载荷示踪粒子对原始无负载粒子时的流场影响相对较小,研究结论对高速流动PIV测试有着重要的参考价值。     

球形粒子在流体中的跟随性

本文从Maxey和Riley的粒子运动基本方程出发,得到了球形固体粒子在相对流动雷诺数很小情形下的分析解。讨论了粒子跟随性对外力、初始条件和流场性质的依赖关系;对均匀湍流场,分析了不同密度比和扰动频率对跟随性的影响     

几种常用格式和湍流模型在高速流动中的适用性研究

针对高超声速飞行器中存在的强激波、激波/边界层干扰、分离、湍流等复杂流动现象,对比分析了当前计算流体力学中的主要空间离散格式及湍流模型,发现不同格式对强激波的分辨率基本相同,Roe和LDE(low diffusion E-CUSP(convective upwind split pressure))格式对摩擦因数和传热系数的模拟优于其他格式;S-A(Spalart-Allmaras)一方程湍流模型计算的摩擦因数比k-ωSST(shear stress transport)两方程湍流模型高10%左右,而后者预测的分离区约为前者的2倍,且分离点靠前.     

近壁湍流相干结构的数值模拟

提出了近壁湍流对称单个相干结构的理论模型,采用紧致有限差分和Fourier谱展开相结合的数值方法研究它的演化问题,得到雷诺应力分布及高剪切层的形成过程与实验结果一致。     

轴对称射流矢量喷管的试验和数值模拟

对基于激波来实现推力矢量的轴对称射流矢量喷管的缩比模型进行了测力和测压试验，用推广到可计算可压缩流的SIMPLE方法其内外流场进行了数值模拟，根据试验和数值模拟结果分析了喷管主流与次流相互作用产生的复杂流场结构、二次流流量和落压比对气动矢量角的影响，在落压比3～6范围内，二次流流量和喷管主流比值增大，气动矢量角增大，两者比值相同时，落压比增大，气动矢量角减小。     

混合解析/数值方法及其在湍流数值模拟上的应用

针对带小时间尺度的源项的方程描述的流动问题,提出了混合解析/数值方法。混合解析/数值方法的基本思想是:分裂原始方程组为对流-扩散部分的偏微分方程和源项的常微分方程。偏微分方程采用合适的数值方法求解,而常微分方程采用解析方式积分。模型方程的理论误差研究表明,混合方法提高了源项处理的精度,降低了混合方法的整体数值误差。分析同时表明,基于时间分裂的算法在求解含源项双曲系统的定常类型问题,会存在数值振荡。为此发展了非分裂方式的混合解析/数值方法,在湍流模型数值计算中提高了数值稳定性,而且加快了计算的收敛速度。     

可压缩流动激波装配在格心型有限体积法中的应用

发展了一种基于格心型有限体积方法(FVM)的激波装配算法。通过定义网格节点属性可以灵活调用激波装配和激波捕捉计算方法。在使用激波装配方法时,激波节点运动速度和下游运动速度通过Rankine-Hugoniot(R-H)关系式获得,同时采用非结构动网格技术描述激波的运动以及调整其他网格节点的位置。流过激波面元的通量为上游单元的基本通量,物理概念更加清晰,通量计算也更为准确。在计算过程中,网格节点属性可以发生变化,以此实现对带有拓扑变化流场的描述。数值试验表明:本文提出的计算方法不但具有较高的计算精度,同时能有效地避免由于捕捉激波而出现的数值问题。