高马赫数层流摩阻数值计算精度

飞行器在高空高速飞行时黏性效应显著,摩阻预示精度对飞行器的关键气动性能意义重大。目前,摩阻预示主要依赖数值计算,但高马赫数层流摩擦阻力计算与试验测量结果仍存在差距。以具有高速飞行器典型部件特征的球锥、三角翼为对象,结合风洞试验摩阻测量结果,使用国家数值风洞（NNW）数值计算软件和自研CFD程序研究了数值计算中影响摩阻计算精度的格式数值耗散及壁面温度边界条件等重要因素。研究结果表明：数值耗散越小,表面摩阻的计算精度越高;在速度较低的边界层近壁区内关闭熵修正,将有助于提高表面摩阻的预示精度。此外,在高马赫数流动问题的数值模拟中,壁面温度条件对表面摩阻计算同样具有重要影响。最后,基于分析结果和工程需要提出了对高精度摩阻数值预示的研究需求。     

超声速尾迹-剪切流的混合增强

超声速条件下燃料和空气之间的高效混合是超然冲压发动机技术上的主要挑战。基于大涡模拟和流动稳定性分析,针对超声速尾迹-剪切流动开展了混合增强方法研究。尾迹的存在改变了混合层流动的速度剖面,对流动稳定性产生了重要影响,使混合层由三维最不稳定变为二维最不稳定,最不稳定扰动波频率和增长率增大。基于流动稳定性结果引入扰动的混合增强方式依然有效,根据稳定性结果设计了波纹隔板。数值结果表明：二维波纹壁引入的扰动未能增长,不具备混合强化效果,而三维波纹壁引入的扰动能够快速增长,具有混合强化效果,且波纹壁参数越接近最不稳定扰动波参数,混合强化效果越明显。     

高超声速条件下星光传输模拟方法与仿真

以高超声速飞行器为研究对象,建立了高速流场下的星光传输模拟仿真方法,定性、定量地研究平台气动光学效应对星光导航的影响。采用基于k-ω/SST两方程湍流模型的雷诺平均法计算星敏感器安装在飞行器不同位置处（头部、中部和尾部）的外流场密度分布,进而得到流场的空间折射率分布。在此基础上,利用几何光学法计算星光由高速层流流场引起的光程差、点扩散函数、像偏移等光学传输效应,得出经过层流流场的星光降晰图。同时运用统计光学理论计算高速湍流流场引起的星光出瞳面的波像差、密度均方差、相位均方差、Strehl比等光学传输效应,得出经过湍流流场的星光降晰图。仿真结果与理论分析吻合,从而为研究星光在高超声速流场中传输的机理提供了可靠的技术手段。     

激波引燃冲压发动机一体化流场数值模拟

激波引燃冲压发动机是一种采用爆轰形式组织燃烧的吸气式高超声速飞行器动力系统。采用AUSMPW+迎风格式、氢氧7组分8步基于反应模型,在非结构网格离散域上求解二维多组分化学非平衡流Euler方程。采用发展的数值方法求解了激波引燃冲压发动机内外一体化流场,研究了台阶长度、斜劈尖角度对发动机流场结构和发动机性能参数的影响规律。计算结果表明,所发展的数值方法能用于激波引燃冲压发动机的一体化流场和性能预示计算;台阶长度和斜劈尖角度影响发动机流场结构和预混气体能量释放程度。推力和燃料比冲均随台阶长度的增大而增大,随斜劈尖角度的增大而减少。     

带抑流楔升力前体设计初探

升力前体由于前缘激波不能附着于前体前缘,其上下表面流场存在相互耦合,进而导致其下表面流场存在较强的横向流动,这会对腹部进气系统的设计造成不利影响.为了使前体出口流场能较好地保持二维流动特征,通过在升力前体两侧增加小角度抑流斜楔,以抑制横向流动.这种升力前体及一体化进气系统的数值模拟表明,该升力前体出口流场二维特征良好,采用该前体的进气系统能较好地维持理想的流场结构,并具有较优的性能.     

基于POD方法的弯曲扩压通道分离流控制的时空特性分析

为了分析一种运用于压气机内分离流控制的无源脉冲射流控制技术的特点,基于弯曲扩压通道试验模型进行了脉冲射流控制的试验和数值模拟研究,结果均表明当射流频率接近通道内分离涡主频时控制效果最为明显;引入了本征正交分解（POD）技术对无控状态下通道内流场结构进行分析,得到了POD各阶模态的流动结构特征。在此基础上对比分析了定常及非定常控制特点,结果表明：非定常控制方式主要是重分配各阶模态之间的能量,有选择性地强化或削弱某阶模态;定常射流控制则是整体削弱高阶模态,压制通道内复杂流动现象;合理地构建脉冲射流可使能量从高阶模态向平均流模态进行转移,能量的转移通过空间流场结构的重构和模态时间演化特性的序化实现。最后针对POD分析结果进行了验证性试验研究,试验结果部分反映了时空特性的变化规律,提升了POD分析结果的可信度。     

翼型前缘变形对动态失速效应影响的数值计算

翼型或机翼的动态失速效应所引起的低头力矩和正气动阻尼限制了飞行器气动性能的提高,甚至可能诱导发生不稳定运动。应用于小尺寸机翼的前缘动态变形（DDLE）技术,通过实时改变前缘形状,能够改善翼型前缘区域的速度梯度,进而抑制动态失速效应。采用转捩剪切应力输运（SST）黏性模型结合分区混合动态网格技术,研究了这种前缘变形对机翼俯仰运动所引起的非定常流动的影响,得到通过小幅度前缘变形抑制和延迟动态失速的方法,从而提高翼型的气动性能。翼型NAC A0012的数值模拟结果与动态失速风洞试验结果比较表明：所使用的数值计算方法能够较为准确地模拟翼型在动态失速过程中升力系数与俯仰力矩系数的变化情况,可用于研究前缘变形对翼型俯仰运动所引起的非定常流动的影响。前缘动态变形翼型俯仰运动过程的非定常流场的数值模拟表明：在大迎角下不同幅度的前缘下垂运动能够抑制流动分离的发生,从而抑制动态失速,但在大迎角下小幅度高频率的前缘下垂变形能更高效地抑制动态失速;前缘变形幅度以及变形沿中弧线的分布对升力系数和俯仰力矩系数的影响并不明显。     

火箭发动机气体-颗粒两相流双流体模型研究

探讨了两相流双流体模型中颗粒相压强、两相界面压强以及两相界面速度等若干关键问题,从理论上推导了双流体模型的控制方程,并讨论了颗粒相控制方程的双曲性.把LU隐式时间格式和Van-Leer矢通量分裂格式应用到两相流双流体模型中,并对一维等截面传热摩擦管道和轴对称JPL喷管中气固两相流进行了计算,计算结果与文献结果吻合良好,所建数学模型正确,所用数值计算方法可靠.     

二维微槽道电渗流的增强混合

微尺度流动中的混合一直是微流系统研制中的一个大问题,对壁面zeta势非均匀分布情况下的二维微槽道电渗流进行了数值分析,重点分析了这种情况下的微流动混合增强.基于电渗流的Helmholtz-Smoluchowski滑移速度模型,应用SIMPLE方法计算了4种滑移速度在整个壁面上的分布不均匀一致时的微槽道电渗流动的流场及被动标量在其中的混合.计算结果表明,滑移速度的非均匀一致导致微槽道电渗流中的流线发生剧烈的转折,从而使得微槽道流动的混合效率有很大的提高.流场中流线的转折越剧烈、转折的次数越多,微槽道流动混合效率的提高就越大.