限制器对高超声速体表面热流数值模拟的影响

采用Roe的FDS（Flux Difference Splitting）差分格式结合min mod,van Leer和Osher-C限制器,对高超声速体外部流场进行数值模拟.研究了限制器对高超声速体气动力/热数值模拟的区别,对比分析了3种限制器对三维模型表面热流数值模拟结果的影响.研究表明,限制器对数值模拟精度有重要影响,且对气动热数值模拟的影响大于对气动力的影响;对简单或较复杂的气动外形,采用Roe的FDS格式结合min mod限制器获得的表面热特性比它结合van Leer和Osher-C限制器时要好,可为热防护设计提供可靠结果.      

高超声速热化学非平衡对气动热环境影响

高超声速气动加热严重,考虑热化学非平衡对气动热环境影响,可以为热防护系统设计提供有效保障。采用Park和Gupta热化学非平衡模型,数值计算研究5组元(N_2,O_2,N,O,NO),17组化学反应的热化学非平衡效应对高超声速飞行器气动热环境影响,并与完全气体和热化学平衡模型进行对比分析。热化学非平衡模型流场温度及激波距离均比完全气体模型小。激波后气体密度因离解、化学反应而增大,且气体密度越大,激波距离越小,热化学平衡模型激波距离最小。完全气体和热化学平衡模型热流载荷计算值均比实验值偏大。Park和Gupta热化学非平衡模型数值计算激波距离及气动力载荷差别小。Park模型热流载荷计算值偏大,Gupta模型与实验结果相符,它可对气动热环境可靠预测。     

高超声速气动热数值计算壁面网格准则

首先对高超声速气动热数值模拟中的3种壁面法向网格准则(平均自由程(MFP)、自由来流参数网格雷诺数和壁面参数网格雷诺数准则)进行了分析.随后,提出了壁面参数预估方法,使得壁面参数网格雷诺数和MFP两种准则无需进行试算即可直接预估壁面法向网格尺度.其次,运用提出的壁面参数预估方法结合不同网格准则,在同一实验条件下确定多个壁面法向网格尺度,通过与实验数据对比研究各网格准则确定的壁面法向网格尺度是否满足热环境模拟精度要求,并对提出的壁面参数预估方法进行验证.最后,开展了高超声速二维完全气体及真实气体效应算例、三维钝双锥算例的数值模拟.研究表明,壁面参数预估方法是可靠的;自由来流参数网格雷诺数准则随来流静温变化所确定的网格尺度变化趋势与物理推理相悖;基于壁面参数网格雷诺数和MFP两种准则主要受壁面温度影响,其在同一条件下所确定的网格尺度基本一致,并满足热环境模拟的精度要求.      

高超声速热流计算湍流模型性能评估

采用计算流体力学方法,通过对高超声速来流下双椭球热流计算,综合分析了影响壁面热流预测精度的因素.针对工程上4种常用的湍流模型:BL(Baldwin-Lomax)模型、SA(Spalart-Allmaras)模型、k-ω模型及SST(Shear-Stress Transport)模型,评估了其在高超声速热流预测中的性能,获得了湍流模型下壁面法向网格雷诺数对热流计算的影响规律.此外,对高马赫数下不同后掠角钝舵模型进行热流计算,获得了后掠角对钝舵前缘热流峰值的影响规律.研究表明,两方程湍流模型相比一方程和零方程模型更适合高超声速气动热计算,其中SST模型对热流预测精度更高,总体性能更为优异;钝舵驻点线最大热流随后掠角增大呈非单调变化,在后掠角约为22°时热流达到峰值.      

高超声速二方程湍流模型的数值模拟对比

对基于Reynolds和Favré(密度加权)混合平均的二方程湍流模型进行了修正,同时根据数值模拟高超声速流动时必须具有高分辨率捕捉间断面与在边界层内抑制数值粘性的能力的要求,提出了新的总变差减小(TVD)格式熵修正函数.在此基础上,通过对压缩拐角的高超声速湍流的数值模拟,对基于Reynolds和Favré混合平均的二方程湍流模型,以及其它不可压缩模型及可压缩性修正模型进行了对比,显示了不同湍流模型及可压缩性修正在计算壁面压力分布和热流分布上的特点,说明了对高超声速压缩拐角型流动,湍流模型可压缩修正的必要性,得到了基于Reynolds和Favré混合平均的二方程湍流模型的计算结果最接近实验结果.      

有翼再入飞行器的超/高超声速气动力工程方法

有翼再入飞行器气动力的工程计算方法尚有待进一步完善,计算精度有待进一步提高。文章根据有翼再入飞行器的气动布局特点,改进和发展了一套适于其气动布局的部件划分策略和压强计算选取准则。以此对航天飞机轨道器和类X-34飞行器的纵向气动力特性进行工程计算,并与其风洞试验和数值模拟结果进行对比分析。结果表明,文中的工程计算方法在马赫数3以上的超/高超声速范围内可准确预测有翼再入飞行器的升阻特性,其中升力系数和阻力系数误差小于10%;并可准确预测俯仰力矩特性随迎角改变的变化规律。与现有方法相比,文中的部件划分策略和压强计算选取准则在升阻特性精度相当情况下,可明显提高俯仰力矩的预测精度,并反映飞行器表面的压强分布特征。     

高超声速再入飞行器气动布局多目标优化

气动布局的多目标优化是飞行器设计中的关键技术。提出一种新的高超声速再入飞行器气动外形参数的多目标优化方法,证明外形优化对高超声速流下飞行器性能的影响。通过实例仿真对飞行器所受阻力和升力对制导性能影响进行详细验证分析,将飞行器落点圆概率偏差、末速大于500 m/s的占比、最大飞行过载小于60g的占比这3个性能指标作为优化目标,将升力特性作为中间参数,将气动布局优化问题分解为2个子问题,通过基于搜索算法的升力特性优化和基于改进的模拟退火算法的外形参数优化,减少优化计算时间、提升计算效率、实现对飞行器主体和襟翼的气动布局优化、获得高超声速流下的最佳飞行器外形。仿真结果表明：在确定的约束条件下,优化算法增加了飞行器在超音速流下的气动升力,有效提高了升阻比。在不影响最大飞行过载的前提下,优化后的飞行器表现出更高的气动性能,显著提升了命中精度,同时末速也满足指标要求,制导系统性能得到有效改善。     

一种高超声速飞行器的非线性再入姿态控制方法

针对高超声速飞行器的再入非线性动力学模型,利用SDRE（state dependent Riccati equation）设计姿态控制器。基于奇异摄动理论,把姿态动力学分解成姿态角和姿态角速度跟踪内、外环回路,同时把非线性动力学伪线性化。每个跟踪回路用SDRE获得控制律,考虑到SDRE局部渐近稳定的特点,可以保证系统闭环稳定。最后设计高超声速飞行器飞行控制系统,并在高超声速条件下进行仿真,验证了该方案的有效性。     

高超声速伸缩式变形飞行器再入制导方法

针对高超声速变形飞行器再入制导问题,提出了一种采用伸缩式机翼的高超声速变形飞行器外形方案,建立了含有展长变形量的气动模型和动力学模型。将该变形飞行器的展长变形量扩展为控制变量,分析了倾侧角、展长变形量和终端航程、高度之间的关系。在此基础上,利用倾侧角和展长变形量在线预测剩余航程和终端高度,通过数值方法校正2个控制量以满足航程约束和高度约束,通过航向角走廊确定倾侧角符号。仿真结果表明:该变形飞行器再入制导方法制导精度高,相比于传统固定外形飞行器终端约束能力更强、轨迹更加平滑,且在扰动条件下具有一定鲁棒性。      

高超声速流动湍流模式评估

选取了超声速二维膨胀-压缩拐角和高超声速双椭球绕流作为基准流动,考察了几种当前CFD(Computational Fluid Dynamics)工程应用较流行的湍流模式:BL(Baldwin-Lomax)模式、SA(Spalart-Allmaras)模式、k-ω模式以及SST(Shear-Stress Transport)模式,通过数值计算结果和实验结果的对比及分析,对有关的湍流模式进行了评估,得到一些有意义的结论,并指出受流动可压缩效应的影响,计算高超声速湍流壁面热流时要对现有的湍流模式进行相应的可压缩修正.      

超声速场中的反向喷流数值模拟

采用了高分辨率的N-S(Navier-Stockes) 方程数值模拟方法,对钝头体头部反向喷流进行研究.计算并分析了头部反向喷流现象中喷流马赫数、来流马赫数、攻角等因素对流场细致结构的影响,并对反向喷流减少阻力和减少气动加热的原理进行了深入分析和探讨.研究分析表明这些因素通过决定弓形激波、马赫盘的强度和位置,来影响回流区的大小和位置,从而使有反向喷流时阻力系数所受到的影响比无喷流时更大.      

高超声速滑翔飞行器再入段制导方法综述

针对高超声速滑翔飞行器再入飞行段,回顾了制导技术的发展历程和研究现状。建立了高超声速滑翔飞行器运动模型,并分析了再入段的路径约束、终端约束和地理约束。将再入制导方法分为三类：标准轨迹制导方法、预测-校正制导方法、混合制导方法,分别对研究现状进行了综述。然后,专门针对侧向平面制导方法进行了讨论和分类,根据飞行任务不同分为了常规约束的制导问题与附加地理约束的制导问题两类。最后,对再入制导方法进行了总结,并结合未来高超声速滑翔飞行器的任务需求,展望了再入制导技术的发展方向。     

高超声速三维化学非平衡流动的数值模拟

高温非平衡化学反应对高超声速再入飞行器绕流流场的气动热特性有较大影响.通过采用有限体积NND(Non-oscillatory, containing No free parameters and Dissipative scheme)格式求解三维N-S方程,分别计算了小攻角条件下球锥的完全气体粘性流场和化学反应非平衡气体流场.比较并分析了化学非平衡作用对三维高超声速钝头体绕流流场的影响.所得结果与国外文献结果符合的较好,证明了计算方法与程序的有效性和实用性.      

TVD格式与高超音速气动加热数值模拟

为了在化学非平衡流动中获得准确的流场解以及表面热流分布,将总变差减小TVD(Total Variation Diminishing)格式中的熵修正函数,由各向同性分布改为各向异性分布,同时让熵修正函数中的参数与流场中的压力梯度分布相关.将改进后的熵修正函数运用到高超音速化学非平衡绕流流场的数值模拟中,获得了较使用原有熵修正函数更为准确的流场参数和表面热流分布.采用改进的熵修正函数,可以提高壁面附近的粘性分辨率,降低热流计算结果对壁面附近法向网格尺度的敏感性.      

数值模拟高超音速化学非平衡流的隐式格式

针对高超音速化学非平衡流提出了两种隐式格式.一种是近似因子分解算法对角化形式在隐式TVD(全变差递减)格式上的推广,另一种可视为隐式TVD格式的LU分解形式.数值试验表明,本文给出的算法保持了TVD格式的高分辨率特性,同时计算量大大减少,是高超音速化学非平衡流数值模拟的有效方法.     

战术导弹横向喷流数值模拟

为了研究喷流对导弹不同升力面的干扰影响,通过数值求解N-S方程来模拟超声速外流场中横向喷流的干扰流场,采用分块对接网格和"O"型网格技术,精确模拟喷口截面及弹翼形状,生成高质量的贴体计算网格.通过几种升力面和弹身组合的超声速横向喷流干扰流场的数值模拟,计算的马赫数为2和4.5,计算结果显示具有升力面的喷流干扰比光弹身的喷流干扰大.研究和分析了喷口附近流场的涡系结构和波系结构,并将导弹喷流干扰力放大因子及喷流引起的压心变化的计算结果与目前国外公布的试验值进行比较分析,二者显示了良好的一致性.     

高温燃气风洞加热特性仿真分析

提出了利用亚音速高温燃气流进行近空间高超飞行器热环境地面模拟的试验方案,在试验装置试验段,通过高温高速的燃气流引射低速的冷气流,达到仅使飞行器头锥驻点附近区域产生局部高温而其余头锥蒙皮表面低温的目的.对某型高超音速飞行器的头锥利用高温燃气进行加热并利用CFD(Computational Fluid Dynamics)方法对13种模拟工况进行了数值仿真分析.将数值模拟计算结果与飞行器在高超音速飞行状态下对应机体部位气动热的理论计算值进行了对比,证实了亚声速高温燃气热环境模拟方法的可行性,为高温燃气地面模拟设备技术方案论证提供了依据.      

航天飞机气动加热计算

介绍了国内外预测航天飞机气动加热的工程计算方法.给出了轴对称、非轴对称和机翼前缘等各种外形驻点热流密度的计算,分析和比较了平板、锥体的层流和湍流热流密度计算的各种方法,概述了航天飞机迎风面中心线和离开中心线横向热流密度的计算.计算结果表明,本文方法用在航天飞机气动热环境初步设计中比较简单并且有足够精度.