高超声速飞行器相关的摩擦阻力直接测量技术

摩擦阻力是高超声速飞行器气动力的重要组成部分,也是制约高超声速飞行器发展的重要因素,因而对摩擦阻力的准确测量就显得尤为重要。简要介绍了近年国内外与高超声速相关的摩擦阻力直接测量技术的发展状况,同时对中国航天空气动力技术研究院白行研制的两套摩擦阻力测量装置作了介绍,并指出应变式摩阻天平技术是测量摩擦阻力的有效途径之一。文中给出了这两种结构形式摩阻天平的静校结果及其在高超声速风洞中的试验结果,并对结果进行了讨论。     

高超声速二维前体进气道一体化优化设计研究

在飞行器前体进气道的一体化优化设计中,最大总压恢复系数是一个必须考虑的参数.本文从二维高超声速进气道的最大总压恢复系数入手,通过理论分析给出了高超声速飞行器从2波系到6波系的二维高超声速飞行器前体进气道的一体化优化设计计算模型.采用拉格朗日乘子法和序列二次规划法(SQP)分别计算了进气道内一道内激波和两道内激波时的情况,给出了进气道的最大总压恢复系数、进气道内马赫数、激波偏转角和激波强度随来流马赫数的变化关系.比较两种方法的计算结果可知采用的计算方法是合理的.     

二维带动力吸气式高超声速飞行器绕流的PNS-NS混合求解

为了提高吸气式高超声速飞行器绕流的求解效率,采用空间推进求解抛物化PNS(parabolizedNavier-Stokes)方程和时间迭代求解Navier-Stokes(N-S)方程的混合计算流体动力学(CFD)方法来求解高超声速飞行器整机绕流.在超声速占主导的流动区域采用空间推进求解抛物化N-S方程的方法,在亚声速和分离区采用时间迭代求解N-S方程的方法.对于求解二维带化学反应的吸气式高超声速飞行器绕流,混合CFD方法和完全时间迭代方法相比,可得到同等准确的数值模拟结果,并且求解效率提高了数倍.      

复杂外形飞行器热流的NND有限元数值计算方法研究与应用

与微分法相比,有限元方法作为一种加权积分的方法,其特点和权函数的特性使得有限元方法对于复杂外形的边界处理和计算诸如热流(与温度的一阶导数相关)等物理量时有其内在的优越性.本文采用张涵信的NND思想发展了复杂流场的NND有限元计算方法,采用积分的方法提高了物体表面热流的计算精度,对高超声速复杂外形飞行器的绕流获得了良好的数值模拟结果.     

基于欧拉方程的高超声速飞行器的壁面流线生成计算

提出了一种新的利用表面流函数法计算高超声速飞行器表面流线分布的方法.首先提出了表面流函数的概念,并通过理论推导,得到了表面流函数与表面流线的关系;然后运用结构化网格求解三维Eu-ler方程,计算得到高超声速钝头体的边界层外缘外部无粘流场气流参数;最后利用无粘流场气流参数和表面流函数的方法计算了高超声速飞行器的精确表面流线分布.结果表明,在无攻角和攻角小于20°的情况下均可以得到较好的壁面流线分布.高精度的表面流线的得到,为利用边界层内粘性主导区域的积分方程等方法进一步精确预测高超声速飞行器表面的气动加热奠定了基础.      

高超声速复杂外形热流的NND有限元数值计算方法研究

有限元法是一种加权积分的方法,这个特点和权函数的特性使得有限元法在计算诸如热流(与温度的一阶导数相关)等物理量时有其内在的优越性.本文研究了高超声速复杂外形热流的NND有限元数值计算方法,并对高超声速球钝锥外形和带翼飞行器外形的绕流进行了数值模拟,算例结果令人满意.     

再入飞行器高超声速绕流的一步和二步有限元计算

本文针对再入飞行器的空气动力问题，采用一步和二步有限元格式求解可缩守恒型Ｅｕｌｅｒ方程，并结合不同的人工粘性模型来捕捉激波，给出了高超声速二维钝头绕流的结果。     

复杂外形飞行器热流的NND有限元数值计算方法

有限元方法是一种加权积分的方法,其特点和权函数的特性使得有限元方法对于复杂外形的边界处理和计算诸如热流(与温度的一阶导数相关)等物理量时有其内在的优越性.本文采用张涵信的NND思想发展了复杂流场的NND有限元计算方法,算例表明,在相同的网格雷诺数下,采用积分的方法提高了物体表面热流的计算精度,对高超声速复杂外形飞行器的绕流获得了良好的数值模拟结果.     

高超声速全机外形气动加热与结构传热快速计算方法

发展了一种无黏流场解与工程计算方法相结合的高超声速全机外形气动加热与结构传热快速计算方法。该计算方法结合了三维块结构网格无黏流场数值计算技术可处理复杂外形流动的优点与工程计算方法效率高的特点,将气动热的计算简化为绕飞行器的无黏外流(边界层以外)数值解和边界层内热流求解两个部分,同时耦合了防热结构传热计算模型、高温化学非平衡热效应估算方法以及弹道状态动态插值方法,可用于快速计算与分析三维复杂外形高超声速飞行器在弹道飞行状态下全机热环境参数、防热结构内温度场等随飞行时间的变化特性。以RAM-CⅡ、类Ⅹ-37B等典型高超声速飞行器为研究对象,在设定的飞行条件及热防护方案下,进行了气动加热与结构传热问题的求解,给出了全机表面热流密度与防热结构材料温度的时变特性。结果对比表明,所发展的方法具有快速、高效的特点,且计算精度可满足工程设计初期选型需求,可为高超声速飞行器的热防护系统初期设计及热环境特性快速计算分析提供技术支持。     

机动飞行器中等攻角高超声速无粘绕流数值模拟

本文采用ＮＮＤ差分格式，利用推进－迭代方法，数值模拟了机动飞行器的中等攻角高超声速无粘绕流流场，发展了一种适用于处理中等攻下背风面流场的技术，大量计算表明，结果准确、可靠。气动力力系数和壁面压力分布与实验数据一致，流场波系结构模拟正确。     

侧向喷流与飞行器三维绕流干扰流场特性的DSMC/EPSM数值模拟

采用DSMC/EPSM混合算法,研究了带有侧向喷流的复杂外形高超声速飞行器稀薄区和连续区三维混合流场的干扰特性。采用一种飞行器物面网格与DSMC计算域网格分别标识的方法,通过判断模拟分子与表面碰撞来完成飞行器物面网格与DSMC计算域网格间的信息传递和信息存贮,对于复杂外形飞行器的精确描述的物面网格不需做进一步处理,直接应用于不依赖于飞行器外形的DSMC计算的通用子程序中。引入碰撞数做为混合参数,对流动区域进行划分,并通过网格中碰撞数的统计和处理,有效地将DSMC和EPSM方法结合在一起。仿真计算了三维复杂外形飞行器带喷流流场压力、热流分布量,飞行器表面气动力、气动力矩和气动热参数,对喷流与高超声速空气绕流相互作用以及它们与物面之间的相互扰动进行了分析,证明了采用的方案和技术的有效性。     

高超声速边界层转捩机理及应用的若干进展回顾

高超声速边界层转捩对飞行器的热传递、表面摩阻和流动分离等有重要影响,尤其是再入飞行器和吸气式巡航飞行器。然而,人们对边界层转捩机理中的很多问题认识还不清楚,或存在争议。本文从扰动波演化的角度回顾了高超声速边界层感受性、线性稳定性和非线性作用的国内若干研究进展,并以基于谐波共振的人工转捩技术为例示范了这些机理认识在转捩控制上的应用。扰动的产生和发展是认识边界层转捩机理的核心。通过研究扰动波来认识边界层转捩机理,开展应用创新研究对提升飞行器性能具有重要意义。     

跨流域高超声速绕流Boltzmann模型方程并行算法

通过对Boltzmann方程碰撞积分进行模型化处理,提出了统一描述各流域复杂高超声速流动输运现象的气体分子速度分布函数控制方程,使用离散速度坐标法对分布函数方程所依赖的速度空间离散降维,构造出直接求解分子速度分布函数的气体动理论耦合迭代数值格式,研制了复杂飞行器高超声速绕流气动热力学计算模型。基于对气体动理论数值计算方法内在并行性、变量依赖关系、数据通信与并行可扩展性的分析研究,使用区域分解并行化方法提出了新型的气体动理论数值算法并行方案;研究了数据的并行分布与并行执行特征,开展了大规模的并行化程序设计,构造了可稳定运行于成千上万CPU的高性能并行算法,用以模拟各流域复杂飞行器的高超声速绕流问题。以稀薄流到连续流环境下不同Knudsen数、不同马赫数的可重复使用类球锥卫星体及翼身组合复杂飞行器等气动力、热绕流问题为研究对象展开大规模并行计算,并进行算法验证,所得计算结果与理论分析、直接模拟蒙特卡罗方法(DSMC)的模拟值及有关实验数据吻合较好,揭示了飞行器跨流域高超声速下的复杂流动机理与变化规律,提供了一条能够可靠模拟高超声速飞行器跨流域气动力及热问题的统一的算法应用研究途径。     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

高超声速壁面黏性力快速计算方法

当前高超声速飞行器壁面黏性力的预测方法主要有精确的Novier-Stokes方程CFD方法(计算流体力学)和快速工程计算方法。前者计算精度高,但效率很低;后者效率极高,但精度难以满足工程应用要求。如何将两种方法有效结合,发展具有一定精度且效率也较高的壁面黏性力计算方法,一直是计算空气动力学研究的重要方向之一。为此,本文基于无黏Euler方程CFD方法获得的壁面流场参数,发展了一种更加精确的高超声速飞行器壁面黏性力快速工程计算方法。该方法以经典参考温度法为基础,以壁面当地流场参数替代自由来流参数,更加充分地考虑了高超声速复杂流动的特性,从理论上提高了壁面黏性力的预测精度。为验证该方法的准确性,以典型高超声速高升阻比飞行器外形为研究对象,分别采用精确Navier-Stokes方程CFD方法和本文发展的快速工程计算方法,预测了高空高马赫数飞行环境下飞行器壁面的黏性力;通过对比分析表明,本文所建立的快速计算方法预测偏差为10~20%,能够满足工程初步设计的要求。     

高超声速滑翔飞行器表面加热特点研究

高超声速滑翔飞行器具有高机动和远程快速到达能力,是高超声速技术应用的前沿领域。快速准确地预测飞行器表面受热变化特征,对高超声速滑翔飞行器热防护系统设计十分重要。以基于锥导乘波构型的高超声速滑翔飞行器和跳跃飞行弹道为研究对象,对其表面不同特征区域的加热开展了快速预测方法研究;采用选用方法分析驻点辐射平衡温度随飞行弹道的变化规律;在弹道特征位置上针对飞行器前缘及上、下表面的加热情况进行分析,获得高超声速滑翔飞行器受热的整体特征,可为分区域选择热防护措施提供参考。     

基于流线跟踪法的气动热工程计算研究

运用结构化网格求解三维Euler方程,计算得到边界层外缘无粘流场气流参数;利用无粘流场气流参数和表面流函数的方法计算了飞行器无粘表面流线分布;在理论和半经验公式的基础上,计算了定比热比和变比热比情况下驻点热流密度,非驻点区域采用参考焓、局部相似性等方法来确定飞行器表面的气动加热,实现了数值算法与工程算法的耦合.上述方法用于求解高超声速钝双锥的表面热流分布,计算结果与经典的热流公式和实验结果进行对比,平均精度为10%左右,满足高超声速飞行器概念研究和初步设计的需要.      

网格对高超声速钝头体表面热流数值模拟结果的影响

高超声速钝头体热流的精确预估，对高超声速飞行器热防护的设计具有重要意义，而网格划分尺度对预估飞行器表面热流精度具有重要的影响。应用理论分析及数值模拟方法，提出了一种基于壁面网格雷诺数及基于钝体特征长度的来流雷诺数网格划分方法，并给出了壁面网格雷诺数的取值范围。应用所提出的网格划分方法，对不同来流条件的高超声速半圆柱及球头表面热流进行了数值模拟，结果表明，应用所提出的网格划分方法，及给出的可供参考的壁面网格雷诺数的取值范围划分网格，在满足高超声速钝体表面热流精度要求的同时，有利于网格合理分布并提高数值计算效率。     

壁面催化对高超声速飞行器气动特性影响

针对高超声速流动中的高温真实气体效应,采用数值模拟求解三维N-S方程的方法研究了壁面催化对高超声速飞行器气动特性的影响规律。研究发现:对于文中所选两类高超声速飞行器———大钝头CEV再入飞行器和仿HTV2高升阻比升力体飞行器,壁面催化对表面压力影响均较小,对剪切应力的影响在飞行器不同部位表现不同:在头部和前缘等强压缩区域,壁面催化对表面剪切力影响明显,在大面积和背风区位置,壁面催化对表面剪切力影响微弱。这是因为壁面催化使得具有更大惯性的大分子气体在近壁处聚集,从而导致更高的速度梯度。由于大钝头外形波阻在整体气动特性中占优,而高升阻比外形头部强压缩区域面积较小,头部强压缩区域剪切力对整体气动特性贡献较小,最终体现为壁面催化对整体气动力影响微弱。     

基于蚁群算法的高超声速飞行器气动布局优化设计

针对目前飞行器优化设计领域遗传算法计算量大、效率低的情况,提出了基于连续空间蚁群算法的高超声速飞行器气动布局优化设计方法.蚁群算法是近年来发展的一种全新仿生算法,具有全局性和高效性等特点,已经成功地应用于离散空间的优化设计.采用连续空间蚁群算法,对高超声速飞行器进行了多变量、多约束下的气动布局优化设计,并与采用遗传算法和约束可变多面体法的优化结果进行了对比,指出了蚁群算法的优点.本文的研究可为蚁群算法应用于复杂、高维的大规模飞行器设计问题提供参考.