稀薄过渡流区横向喷流干扰效应数值模拟研究

为了准确预测稀薄过渡流区横向喷流与稀薄大气的干扰流动特征,本文在建立直角与表面非结构网格混合结构的DSMC数值算法以及碰撞网格自适应算法基础上,提出基于MPI的静态随机负载平衡技术,构建了DSMC并行计算代码。计算的不同压力比条件下三维平板模型横向喷流与稀薄大气干扰的分离长度与低密度风洞试验有较好的一致性,验证了本文算法的可靠性。开展了细长钝双锥外形高超声速稀薄来流与超声速/高超声速横向喷流干扰效应的研究,计算分析了不同飞行高度、不同飞行速度、不同飞行攻角、不同喷流推力下复杂流场结构和对气动力特性的影响规律。考察了RCS喷管出口参数不同(均匀/非均匀)对喷口附近分离涡和分离长度的影响。     

天宫飞行器过渡流区高超声速绕流N-S/DSMC耦合计算

针对类天宫飞行器服役期满再入大气层多次解体过程残骸碎片绕流对气动特性干扰影响的问题,在数值求解N-S方程的CFD和DSMC方法程序研究基础上,采用MPC耦合处理原理,发展了基于Chapman-Enskog非平衡速度分布函数的耦合区域双向信息交换亚松弛计算技术,建立了适于残骸碎片两体不同间隔绕流干扰流场的N-S/DSMC耦合算法。对天宫飞行器两舱结构体在过渡流区低密度风洞试验状态进行了耦合计算,计算结果与风洞试验结果吻合很好,验证了所建立的耦合算法对再入解体残骸碎片气动特性计算的有效性和可靠性。通过对天宫飞行器解体碎片的简化外形球和球柱体在过渡流区高超声速两体干扰气动特性进行不同间隔绕流场耦合计算分析,结果表明,在一定的距离范围内,两体干扰会对残骸碎片的气动特性产生较强的影响,不同的残骸碎片外形在不同距离的干扰规律不一致。计算结果为天宫飞行器再入解体气动融合轨道数值预报提供设计依据。     

高超声速气动热预测技术及发展趋势

高超声速气动热预测技术是高超声速飞行器发展的关键技术之一。对高超声速气动热预测技术的发展情况进行了分析探讨。首先,简要回顾了国内外高超声速气动热理论预测及地面实验技术的发展历程;在此基础上,结合典型外形的计算与风洞试验结果的比较,重点介绍分析了气动热工程计算方法、数值模拟方法、气动热风洞试验设备的模拟能力及目前实验测试技术的研究水平;最后,对气动热预测技术的发展趋势进行了讨论,提出了气动热预测技术应研究解决的问题。     

过渡流区N-S/DSMC耦合计算研究

在已有CFD和DSMC方法程序研究基础上,采用MPC耦合处理原理,发展适于流场分区信息交换的亚松弛技术,建立了N-S/DSMC耦合算法.通过对近连续过渡流区不同Kn数条件下钝锥绕流计算与低密度风洞试验结果比较分析,一方面证明了本文所建立的N-S/DSMC耦合计算模型和方法在低Kn数过渡流区的有效性,另一方面剖析了N-S/DSMC耦合计算误差,可望发展一个工程适用的近连续过渡流区高超声速飞行器气动特性预测分析工具.     

LENS风洞试验返回器模型气动热特性模拟

国家数值风洞高超声速流动模拟软件HyFLOW的研制对打破国外同类软件的技术壁垒具有重要意义。与国外DPLR软件进行了对比研究，同时系统介绍了HyFLOW软件求解器的数值方法、物理化学模型以及壁面催化特性计算模型等主要方法，采用典型算例对有限催化模型进行了数值验证，最后基于LENS风洞试验146 mm返回器标模外形开展了高超声速气动热特性数值模拟。研究结果表明，HyFLOW软件在高超声速热化学非平衡流动模拟与评估方面的气动力计算精度高，与国外同类软件DPLR相当，同时其壁面催化条件下的气动热计算精度可靠，可信度高。     

类AHW助推滑翔飞行器气动布局优化设计分析

AHW作为美国首飞即成功完成高超声速助推滑翔飞行试验的飞行器受到越来越多的关注。本文针对该类钝双锥十字形小尺寸弹翼外形气动布局阶段的共性问题进行了研究。基于参数化方法建立的基本外形和工程气动力估算模型,采用正交设计方法进行了参数敏感性分析,并对正交优化结果进行了数值模拟验证分析,在此基础上利用多目标优化方法完成了弹体布局的进一步优化。同时在优化外形基础上考虑气动特性以及总体和防热需求,对操纵面的类型、质心与操纵面尺寸配置以及操纵效率进行了探讨分析,给出了气动布局建议。研究表明,该类布局方式可以获得较高的配平升力、配平升阻比及容积率,并且合理的质心位置/舵面尺寸的组合可以实现操纵性需求,是高超声速滑翔飞行器的一种潜在布局方案。     

风洞非均匀流场中的一体化高超声速飞行器缩比模型气动性能研究

准确预测气动推进性能是吸气式高超声速飞行器研究的重要挑战之一。针对CARDC吸气式高超声速实验室（AHL）自主设计的一体化高超声速飞行器风洞试验模型,通过数值模拟计算,研究了CARDC600mm脉冲燃烧风洞的流场,并与试验结果做了对比,确定了试验模型在风洞中的合理安装位置,分析了带舵面飞行器在进气道打开、发动机不工作情况下的气动性能,对比研究了试验模型部分处于风洞流场非均匀区时,风洞结果对模型气动性能产生的影响,对比了数值计算结果和风洞试验结果。结果为利用风洞试验结果准确分析飞行器气动性能提供了重要依据。     

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

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

锥导乘波构型设计、优化与分析

乘波构型是高超声速飞行器高升阻比气动布局设计的参考外形之一,设计中需要综合考虑升阻比、容积率和容积等要求。在锥导乘波构型参数化设计的基础上,采用工程估算和计算流体力学相结合的方法,通过正交试验设计分析了不同参数对目标影响的敏感性,合理选择设计参数优化区间,应用改进的多目标遗传算法对乘波构型进行了优化设计,针对优化外形开展了气动性能的数值模拟研究,并在高超声速炮风洞中完成了缩比模型的验证性实验。结果表明:优化设计外形具有良好的升阻比,且在一定攻角范围内升阻比较高,数值模拟和实验分析基本吻合。研究结果可为高超声速滑翔式飞行器的设计提供参考。     

稀薄过渡流中简单外形体气动特性的实验与工程计算方法

本文扼要介绍简单外形体在低密度高超声速风洞中的实验及稀薄过渡流态下工程计算方法的特点和类型。 在实验数据和有关文献的基础上,本文对零攻角尖、钝锥阻力系数、横流圆柱体的驻点Nusselt数、零攻角半球柱体的驻点Stanton数等,给出了较为简便的经验公式或近似公式,计算结果与实验数据作了比较和讨论。     

高超声速再入飞行器气动特性的快速预示—局部方法的推广应用

本文将计算高超声速稀薄气流过渡领域中气动特性的局部方法,推广应用到连续介质中弹头型高超声速再入飞行器气动力特性的快速估算。由激波风洞中M_∞=9.9时,一个8°钝锥的气动力测量结果,导出这一实验条件下的领域系数,并以此来估算不同锥角、不同钝度比及不同外形弹头型再入飞行器的气动力和力矩系数,其结果与无粘数值解及实验结果作了比较,在攻角2°～14°范围内吻合得很好。局部方法可用于弹头型高超声速再入飞行器气动特性的快速预示。     

火星进入器高空稀薄气动特性

针对火星稀薄大气环境进入器气动特性问题,以类火星科学实验室外形为例,计算分析火星稀薄大气真实气体效应对气动特性的影响,给出火星高空稀薄环境下的气动特性规律。研究发现,随着飞行高度的增加,稀薄度增加,激波脱体距离、激波厚度增大,激波强度减弱,明显的激波结构逐渐消失,流场等值线更趋于圆弧状分布;真实气体效应使得迎风面压缩及背风面膨胀增强,轴向力、法向力及顶点力矩系数等预测结果与完全气体模型预测结果相比绝对值偏大;随着稀薄度增大,轴向力、法向力及顶点力矩系数等绝对值增大,在同样的迎角下,随稀薄度的增加,纵向压心前移,进入器的静稳定性变差。     

电离对高超声速稀薄流飞行器气动热影响

将电离反应模型扩展到（direct simulation Monte Carlo,DSMC）方法中,研究了电离反应效应对高超声速稀薄流飞行器气动热的影响特性.针对稀薄流场中电子出现带来的实际困难,引入“捆绑法”思想处理电子在流场中的运动,并给出了电离反应模型及电离反应处理方法.在以RAM-C Ⅱ飞行器外形为例对增加了电离反应的DSMC代码进行验证的基础上,以“星尘号”探测器外形为研究对象,针对不同飞行高度下5组元混合气体模型（无电离）和11组元混合气体模型（含电离）的化学非平衡流动开展了数值模拟,细致分析和对比了电离反应效应对探测器气动热的影响规律.研究结果表明:采用的电离反应处理方案能够模拟带电离反应的高超声速化学非平衡稀薄流动.在飞行高度为60km时电离反应对探测器气动热的影响最为强烈,使探测器的驻点热流密度降低了5.12%,电离反应对探测器气动热的影响随气体稀薄程度增加而减弱.      

A parameterized geometry design method for inward turning inlet compatible waverider

Intensive studies have been carried out on generations of waverider geometry and hypersonic inlet geometry. However, integration efforts of waverider and related air-intake system are restricted majorly around the X43A-like or conical flow field induced configuration, which adopts mainly the two-dimensional air-breathing technology and limits the judicious visions of developing new aerodynamic profiles for hypersonic designers. A novel design approach for integrating the inward turning inlet with the traditional parameterized waverider is proposed. The proposed method is an alternative means to produce a compatible configuration by linking the off-the-shelf results on both traditional waverider techniques and inward turning inlet techniques. A series of geometry generations and optimization solutions is proposed to enhance the lift-to-drag ratio. A quantitative but efficient aerodynamic performance evaluation approach （the hypersonic flow panel method） with lower computational cost is employed to play the role of objective function for opti- mization purpose. The produced geometry compatibility with a computational fluid dynamics （CFD） solver is also verified for detailed flow field investigation. Optimization results and other numerical validations are obtained for the feasibility demonstration of the proposed method.     

高超声速稀薄流的气粒多相流动DSMC算法建模研究

基于直接模拟Monte Carlo(DSMC)方法,构造适用于DSMC算法的固态和液态颗粒碰撞、聚合和分离模型,发展稀薄条件下双向耦合作用的气粒多相流的DSMC算法,在此基础上初步实现高超声速稀薄流环境中的气粒多相喷流流场数值模拟.算例结果表明该方法能为稀薄过渡区气粒多相流动提供一种新的应用研究手段.     

升力体机动飞行器气动布局概念设计

开展高超声速机动飞行器的气动布局概念设计研究,提出一种新的升力体机动飞行器气动布局方案,进行了机体的优化设计,并进行了控制舵的匹配设计,研究了飞行器的气动特性和操纵效率问题。研究表明,该方案可以获得较高的配平升阻比及配平攻角,有较高的机动控制效率,是高超声速飞行器实现高机动的潜在可行方案。     

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

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

Control-oriented Modeling for Air-breathing Hypersonic Vehicle Using Parameterized Configuration Approach

This article presents a parameterized configuration modeling approach to develop a 6 degrees of freedom (DOF) rigid-body model for air-breathing hypersonic vehicle (AHV). The modeling process involves the parameterized configuration design, inviscous hypersonic aerodynamic force calculation and scramjet engine modeling. The parameters are designed for airframe-propulsion integration configuration, the aerodynamic force calculation is based on engineering experimental methods, and the engine model is acquired from gas dynamics and quasi-one dimensional combustor calculations. Multivariate fitting is used to obtain analytical equations for aerodynamic force and thrust. Furthermore, the fitting accuracy is evaluated by relative error (RE). Trim results show that the model can be applied to the investigation of control method for AHV during the cruise phase. The modeling process integrates several disciplines such as configuration design, aerodynamic calculation, scramjet modeling and control method. Therefore the modeling method makes it possible to conduct AHV aerodynamics/propulsion/control integration design.     

吸气式高超声速飞行器内外流同时测力试验

为研究吸气式高超声速飞行器新型测力试验技术的可行性,在马赫数Ma=6、高度H=27km模拟飞行环境下开展了风洞试验。以轴对称吸气式高超声速飞行器无舵翼简化构型为研究对象,通过合理的试验系统设计解决内外流解耦关键技术,采用两个环式六分量天平在同一车次上测量模型内外流气动力载荷。试验结果表明:内外流气动力天平测量数据和全模气动力数据准确反映了内外流解耦、内外流窜流、进气道起动/不起动以及溢流影响等物理现象;从物理上实现了内外流解耦;内流气动力对全机气动性能贡献大;内外流解耦、进气道起动时,气动力数据准确、重复性好;窜流产生的传力属于内力。试验证明,同时测力试验技术可行,为解决内外流气动数据源于不同试验所致的不确定性问题提供了有效手段。