高温化学非平衡湍流边界层脉动量象限分析

高超声速飞行器在较低空域以极高马赫数飞行时，表面会同时存在湍流与化学非平衡流动，但目前针对此类高温化学非平衡湍流边界层流动特性的研究工作还比较有限，对不同湍流特征的主导流动机制的认识还有待于进一步深入。选取高超声速楔形体头部斜激波后的流动状态，设置3种不同的壁面温度，通过直接数值模拟对比了不同壁温条件下的边界层参数分布特性，并采用象限分解技术分析了边界层不同象限流动事件对雷诺剪切应力、湍流热流、湍流质量扩散的贡献。结果显示：在整个边界层中上抛和下扫运动对雷诺剪切应力的贡献占优；冷壁效应会使得流向和法向湍流热流通量的主导流动事件在温度峰值两侧发生改变。O原子组分流向湍流组分扩散主要受到高质量分数流体慢速运动事件和低质量分数流体快速运动事件的影响，而法向湍流组分扩散则主要受到高质量分数流体向上运动事件和低质量分数流体向下运动事件的影响。     

高温气体效应对高超声速磁流体控制的影响

高温气体效应会严重影响高温气体流场的流动特性，进而影响高超声速磁流体控制效率。基于低磁雷诺数假设，通过耦合求解带电磁源项的三维Navier-Stokes流场控制方程和电场泊松方程，开展完全气体模型、平衡气体模型、化学非平衡气体模型、热化学非平衡气体模型等条件下的高超声速磁流体控制数值模拟，分析气体模型对磁流体控制的影响，研究高温气体各种非平衡效应及焦耳热振动能量配比等对高超声速磁流体控制的影响规律。研究表明：化学非平衡效应对高超声速磁流体控制影响显著，采用化学非平衡气体模型模拟得到的磁控增阻特性介于完全气体模型和平衡气体模型之间，平衡气体和完全气体模型磁控热流变化的定性规律，与非平衡气体模型模拟结果差异很大；热力学非平衡效应对高超声速磁流体控制的影响，与焦耳热振动能量作用比率紧密相关，随该配比增大，磁场增阻效果由67%降到约12%；高温气体效应会极大地降低磁控增阻效果，会明显地增强部分表面区域的磁控热流减缓效果，要准确数值模拟高超声速磁流体控制，必须有效地考虑化学和热力学非平衡效应，同时选用接近实际情况的焦耳热振动能量配比。     

热化学非平衡高超声速平板边界层线性稳定性分析

高超声速边界层流动转捩是近期空气动力学研究的热点问题。对于环境扰动较小的自然转捩过程,稳定性分析已被证明是研究扰动演化的重要手段。另一方面,高超声速边界层内的温度会随着马赫数的升高而快速上升,极高的温度会引起所谓的高温真实气体效应,使得量热完全气体假设失效,从而对边界层稳定性和转捩产生影响。本文针对高温热化学非平衡气体,利用空气5组分模型开展了平板边界层的线性稳定性分析,重点研究了热化学过程对模态稳定性的影响,并探究了边界层离散谱模态的演化和同步过程。研究表明,对于由第二模态主导的高超声速二维边界层:(1)扰动相比基本流更趋向于热化学冻结态;(2)扰动方程中新出现的非平衡源项的扰动项对稳定性影响很小,非平衡过程主要是通过改变基本流剖面来间接影响稳定性;(3)声速是影响第二及更高模态的重要参数,热化学平衡态假设引起的声速计算式的变化能够解释边界层温度和厚度降低时第二模态频率反而降低的非常规趋势。     

高超声速尾迹雷达散射湍流场研究

为给高超声速再入尾迹亚密湍流雷达散射分析提供所需的脉动背景场参数 ,本文提出计算非平衡再入湍流尾迹脉动等离子体场的理论方法。在研究高超声速尾迹流动特征的基础上 ,推导、使用包括化学组份浓度脉动强度的k ε g湍流模型 ,用以封闭高超声速尾迹雷诺平均控制方程 ,并用全隐式有限差分求解。以M∞ =2 1 .2 6、Re∞D =1 .33× 1 0 6的小钝锥体流动为例 ,得到的结果说明 :流场的流向和径向参数分布合理 ;在转捩点后较近距离内湍流脉动影响较大 ,随着向下游流动脉动影响迅速减弱 ;本文计算尾迹湍流脉动等离子体场的方法是可行的     

真实气体效应对双楔绕流流动特性的影响

为研究高马赫数下真实气体效应对双楔绕流流动特性的影响，采用量热完全气体模型和化学非平衡气体模型两种气体模型对层流流态、转捩流态、全湍流流态的双楔绕流流场进行了数值模拟，其中转捩流态计算采用基于间歇因子的γ-Re_θt转捩模型，全湍流计算中采用k-ωSST湍流模型。结果表明：双楔拐角附近发生流动分离，与量热完全气体模型相比，相同流态下，化学非平衡气体模型计算得到的分离区尺寸更小。在边界层转捩位置的分析中发现，双楔拐角附近的分离区对边界层转捩有重要影响，分离激波产生的位置与边界层开始转捩的位置高度吻合，与量热完全气体模型相比，化学非平衡气体模型时的边界层转捩起始位置更靠后。针对真实气体效应下单位雷诺数、拐角角度对边界层转捩的影响规律研究发现，真实气体效应在一定程度上抑制了分离区内的转捩流动，随着单位雷诺数减小、拐角角度增大，真实气体效应对边界层转捩的影响随之增大。     

甲烷/空气超声速燃烧流动数值模拟

用流动项与化学反应生成源项解耦处理的化学非平衡流动计算方法，从薄层近似三维N-S方程出发，采用ENO差分格式数值模拟了超声速冲压加速器简化模型中，高速甲烷气流从后体多个喷口射入高超声速空气流形成流场，研究了高温异质气体效应和19组分65反应模型的非平衡效应对冲压加速器表面压力分布的影响。计算表明新型解耦方法适合反应机理复杂的碳氢燃料超燃冲压发动机内部流动模拟，为开发应用软件系统打下基础。     

简化三方程转捩模型在高超声速流动中的应用

在雷诺平均方法的基础上,通过耦合k_(-ω)SST湍流模型和间歇因子转捩模型,引入湍流模型和转捩模型的可压缩修正方法,对高超声速平板、双楔、尖锥三类模型边界层转捩流动开展了数值模拟研究。与实验结果的对比分析表明基于压力梯度表征参数T_w=R_(TΩ/ω)的简化三方程转捩模型,能够准确捕捉高超声速平板边界层流动的转捩起始位置、转捩区域长度以及湍流区壁面热流。而对于双楔、尖锥模型,改进前的简化三方程转捩模型由于受到流动可压缩效应的影响,边界层转捩后湍流区的壁面热流模拟预测结果明显高于实验值。在添加模型可压缩修正方法后,转捩区域长度和湍流区壁面热流模拟结果得到有效改善,与实验值吻合较好。可见,简化三方程转捩模型在添加可压缩修正方法后具备准确模拟预测高超声速边界层流动转捩的潜力。     

电导率模拟对高超声速MHD控制影响

高温气体电导率是高超声速电磁流动控制数值模拟最重要的参数之一。针对电导率模拟准确性及其对高超声速磁流体控制影响的问题，考虑高超声速飞行器流场中高温气体热化学非平衡效应，采用三维低磁雷诺数磁流体动力学（MHD）数值模拟方法及程序，结合国内外常见电导率处理方法开展典型状态高超声速MHD控制数值模拟，分析电导率模拟对高超声速磁流体流场分布、气动力/热特性的影响。研究表明：磁控热流减缓效果与电导率呈非线性关系，电导率较大时将出现电导率的磁控热饱和现象，其产生的原因可能与化学反应趋向于平衡态存在一定关系；采用定电导率方法，会人为放大磁场洛仑兹力的磁阻力效果，使阻力系数的预测值偏大；不同电导率模型计算得到的电导率分布差异很大，甚至存在数量级的差别，显著影响了磁流体的控制效果，这与电导率模型的适用范围、参数选取原则存在很大关联；对于含多种离解、电离组分的高温气体流动来说，采用基于多电离组分迁移碰撞的电导率模型（本文模型M8），计算与试验一致性最好。     

用PIV实验研究边界层内结构函数标度律

湍流标度律的变化情况显示了湍流边界层内的流动状况。利用高分辨率、高帧率的粒子图像测速技术（PIV）对平板湍流边界层进行了实验研究。研究结果显示,当脉动速度增量结构函数的阶数为正整数且逐渐增大时,由实验所得的脉动速度增量结构函数的标度指数逐渐偏离K41和SL标度指数,出现奇异标度律;而当结构函数的阶数为正分数时,实验所得的标度指数和K41及SL标度指数较为接近;当结构函数的阶数由0趋向于-1时,实验的标度指数明显偏离K41和SL标度指数,出现奇异标度律。这也说明了用低阶结构函数的标度指数也能说明标度律的奇异性及边界层内存在间歇性。对于法向速度增量和涡量分量来说,不管结构函数的高阶或低阶的标度指数都不随壁面位置的不同而改变,具有普适性。     

高超声速边界层转捩的若干问题及工程应用研究进展综述

高超声速飞行器绕流存在着激波、边界层、流动分离、稀薄气体效应和高温气体效应等多种复杂流动现象的空气动力学问题,其中高超声速边界层转捩既是空气动力学的基础问题,也是高超声速流动研究的热点和难点。若能对边界层转捩进行准确预示及有效控制,则可以实现对飞行器气动力热特性的精细设计,改进飞行器性能,提高任务执行能力。文章针对工程中具有复杂外形飞行器存在的典型失稳特征进行了研究进展回顾,提出了工程实际中亟需解决的复杂边界层转捩问题,明确了高超声速边界层转捩研究的工程应用方向。文章最后还对高超声速边界层的流动控制进行了回顾,以期在今后高超声速飞行器设计中实现对边界层的流动控制,提高飞行器的飞行性能。