高超声速绕钝体熵层基本流特性研究

基于N-S方程,针对超声速气体绕0°迎角钝锥流动采用三阶WENO(Weighed Essentially Non-Oscillatory)格式进行数值模拟,研究了熵层的基本流特性。研究给出了边界层及熵层外缘的确定方法,并分析了其适用性与局限性。数值计算结果显示基本流剖面存在"双广义拐点"现象,N-S方程与Euler方程数值结果对比发现,熵层区出现的广义拐点是由熵层区有旋运动引起的,而边界层的广义拐点由粘性作用导致。这是进行流动稳定性分析的基础。文章还分析了球头半径与马赫数对熵层特性的影响。结果表明,马赫数越大,熵层区的涡量越大,熵层影响的范围越大;雷诺数对无粘流的特征影响很小,熵层大小随球头半径几乎呈线性关系。     

有局部稀薄气体效应的高超声速流动数值模拟

近空间高超声速飞行器当飞行高度和速度足够高时,其流场计算可能要考虑稀薄气体效应,传统的计算流体力学(CFD)方法预测的阻力和升阻比将不够准确。而现有的模拟稀薄气体流动的计算方法由于其计算量巨大,难以在工程实际中应用。因此需要发展能用于近空间高超声速飞行器流场的可行、可靠的计算方法。陈杰和赵磊在文献[1]中针对边界层中既有强剪切而气体分子自由程又相对较大的情况进行分析,提出了刻画此类局部稀薄效应的无量纲参数Zh,并提出了在传统CFD中通过采用依赖于Zh参数的等效黏性系数考虑局部稀薄效应对阻力计算影响的研究思路。因此,本文尝试将此等效黏性系数纳入CFD模型中,以在70km高空,以马赫数15飞行的小迎角钝平板为例,来检验计算方法是否合理可行。结果表明:和传统的CFD方法所得结果相比,新模型计算的阻力减小,升阻比增加,其改进的方向与现有飞行试验结果定性相符,且所增加的计算时间非常有限,可方便地应用于现有的计算空气动力学中。     

气体稀薄效应对热流计算的影响

临近空间高超声速飞行器的研发越来越受到关注。周恒和张涵信在《中国科学》上的《空气动力学的新问题》一文中分析了在此情况下现有的空气动力学的不足,其中之一就是需要考虑流场中可能出现的局部稀薄气体效应,并且针对临近空间高超声速飞行器边界层进行了详细的分析。陈杰和赵磊在此基础上研究了强剪切下的气体稀薄效应,给出了判别气体稀薄效应的无量纲参数Zh,并提出对传统连续介质模型中的黏性系数基于Zh参数进行修正。本文将通过类似的研究方法,采用DSMC (Direct Simulation Monte Carlo)研究纯导热问题中的气体稀薄效应,通过对粒子速度分布函数的分析提出了相应的刻画气体稀薄效应的参数ZhT,获得了依赖参数ZhT的导热系数修正规律,并进一步将该修正规律纳入CFD (Computational Fluid Dynamics)计算,验证了该算法对圆柱绕流问题表面热流预测的精度。