航天器高超声速气动力的数值计算

本文把航天器的外形用数以百计的平面四边形元素来近似。利用高超声速绕流时面积元素之间的干扰可忽略的假设,分别计算出各面积元素的气动力,积分后得出航天器的气动力,算例的结果是令人满意的。本文也描述了航天器表面网格的生成技术和计算机绘制航天器在不同视角下外形图的原理。本文方法可供航天器初步气动设计时用。     

高超声速飞行器气动力/热参数辨识研究综述

飞行器气动力和气动热参数辨识是高超声速飞行器设计的关键技术之一。笔者对高超声速飞行器的气动力和气动热参数辨识技术进行了综述。介绍了飞行器气动力、气动热参数辨识的基本原理与主要方法,气动力、气动热参数辨识技术在高超声速飞行器研发中的应用情况与发展趋势。同时也简要介绍了中国空气动力研究与发展中心（CARDC）在气动力和气动热参数辨识研究方面的研究概况。     

二元高超声速进气道激波振荡特性实验

通过在进气道/隔离段模型出口附近设置固定的堵塞楔块提高反压,并采用高速纹影拍摄同步壁面压强测量的手段,在马赫数5.9的激波风洞中研究了二元高超声速进气道在不同堵塞度下的流动特征。研究结果表明,在较低的堵塞度下,进气道仍然能够保持起动状态,而在较高的堵塞度下,进气道出现激波振荡。上游产生的压缩波/激波在节流段的反射是出现激波振荡的重要原因之一。随着堵塞度的增加,激波振荡的频率有所升高。     

激波风洞高超声速摩阻直接测量技术研究

介绍了在中国空气动力研究与发展中心（CARDC）激波风洞中进行的摩阻测量技术研究情况。在测量研究中,设计了压电型摩阻天平,为了提高摩阻天平的校准和风洞试验测量结果精度,便于风洞试验和校准之间安装的变换,本项研究的摩阻天平采用一种新结构,也就是测量表面和摩阻天平本体可以分离的分体式结构,由此确保在不同使用场合下,摩阻天平的测量表面或者校准加载块可拆卸和更换。验证性试验是在CARDC0．6m激波风洞中进行的,流场名义马赫数分别为8和10,单位雷诺数分别为2．85×10^7／m和1．58×10^7／m,试验中测量了带压缩拐角的进气道模型表面三个测点的摩擦阻力,也测量了摩阻测点及其附近热流,测量结果表明：模型表面的摩阻和热流与雷诺比拟准则符合得较好。     

高超声速粘性激波层高雷诺数流场的数值求解方法

本文将拉伸变换引入Hosny的高超声速粘性激波层的数值求解方法中,数值计算结果与Little所给出的实验结果相符合。从而不论在雷诺数范围上,或在后身流场的推进计算上,本方法均比Hosny原方法有明显改进。     

高超声速气流中平板上粘滞激波层的数值研究

本文描述绕锐前缘平板高超声速流动的数值模拟结果。     

弯曲激波压缩曲面的二元高超声速进气道研究

以弯曲激波压缩面设计的二元进气道具有较短的长度和较高的性能,研究了以等压力梯度弯曲激波压缩曲面设计的二元进气道一般性能,并以部分等熵压缩前楔板设计的二元进气道的一般性能与之作对比,以探讨前者在Ma 4～7高超声速范围进气道的外压缩面设计上的应用。对以等压力梯度设计的型面开展分析,进行了2D粘性计算并获得了稳定流场,显示进气道有较高的压缩效率。在相同总折转角条件下,通过对比不同起始角方案的进气道总压恢复和长度比,前楔板起始角不宜超过7°,进气道长度比可缩短约10%。在对比计算等压力梯度设计的进气道与部分等熵压缩进气道后,认为前者总压恢复性能要下降3%～5%,但激波封口设计点可以适当提前。最后尝试以流线追踪法扩展到三维乘波前体,在设计点附近获得平直的激波面,显示具较好的乘波特性和二元流动特征。      

SFD方法模拟高超声速进气道中激波与边界层干扰

采用ＳＦＤ方法,在ＥＸＣＥＬ软件上利用ＮＮＤ格式数值模拟了激波与平板边界层之间的相互作用以及激波与管道壁面边界层的相互作用,即高超声速进气道隔离段内激波串流动现象。该方法的特点是简单、直观而且不用常规编程,计算结果基本反映了流场的基本结构,表明ＳＦＤ方法是一种具有潜力的计算方法。     

高超声速升力体气动力气动热数值计算

升力体气动布局，由于其良好的高超声速气动特性、有效的内部空间和有益的防热特性，越来越受到气动工作者的重视，本文应用张涵信院士的NND格式思想，用数值模拟的方法对类似于美国X－33的升力体外形进行了气动力、气动热研究。特别是在粘性项的离散和表面热流率的计算中应用了散度和梯度的积分定义，避免了数值计算的奇异性，保证了通量守恒，提高了计算效率。结果表明气动热计算结果与美国同类飞行器结果接近，驻点热流率与经典理论结果一致。     

熵层对高超声速二维钝楔气动参数的影响

采用薄激波层理论和流管质量守恒相结合的方法,分析了高超声速二维钝楔边界层外缘熵的分布规律,研究了熵层对边界层外缘密度、马赫数以及壁面气动加热等气动参数的影响.结果表明:边界层外缘熵在倒圆-肩部区下降最为剧烈,熵层对气动参数的影响在高超声速下不可忽略,特别是使转捩区和湍流区的气动加热增加约53.6%.因此,将表面流态控制在层流模式对高超声速飞行器热防护具有重要意义.      

Influence of entropy layer for hypersonic on aerodynamics parameters 2-D blunt wedge

采用薄激波层理论和流管质量守恒相结合的方法,分析了高超声速二维钝楔边界层外缘熵的分布规律,研究了熵层对边界层外缘密度、马赫数以及壁面气动加热等气动参数的影响.结果表明：边界层外缘熵在倒圆-肩部区下降最为剧烈,熵层对气动参数的影响在高超声速下不可忽略,特别是使转捩区和湍流区的气动加热增加约53.6%.因此,将表面流态控制在层流模式对高超声速飞行器热防护具有重要意义.     

类乘波体飞行器的气动力工程计算

为估算高超声速类乘波体飞行器的纵向气动力,以三角形面元逼近飞行器外形,根据飞行器迎/背风面情况分别采用达黑姆-巴克法、普朗特-迈耶耳法、切锥法和膨胀波方法计算无粘气动力;通过经验公式估算粘性阻力,并考虑飞行器主要部件之间的气动干扰情况,计算了翼片之间、翼身之间的气动干扰因子,得到整个飞行器的气动力。为验证该方法,以某飞行器为例进行了计算,计算结果与CFD吻合。     

高超声速多体分离过程气动加热计算技术

发展了无粘外流解与气动热工程方法相结合的用于全机外形复杂流动气动加热计算技术,并将其与动态多体分离过程流场解相结合,提出了一种高超声速飞行器多体分离过程气动加热特性计算技术,然后耦合结构传热,实现了飞行器热防护系统中结构材料温度分布特性的数值模拟。采用本文方法,针对假定高超声速多体分离布局及飞行条件开展了分离过程中的结构温度分布特性与热流密度分布特性计算,分析并给出了组合体分离过程中部件干扰对气动热的影响特性。结果表明本文发展的计算技术可为高超声速飞行器的防热方案设计、气动热特性定性分析及热防护系统设计等提供技术支持。     

带减阻杆高超声速飞行器外形气动特性研究

采用高超声速风洞测力试验方法测量钝头飞行器头部减阻杆的高超声速气动特性,研究减阻杆的气动减阻原理,分析了多组不同构型减阻杆的减阻效果。结果表明,减阻杆显著减少了钝头飞行器高超声速的阻力,最大的减阻率达到60％之多;减阻效果与减阻杆构型和迎角状态密切相关;减阻杆会诱发稳定性、“热斑”以及非定常脉动等不利问题。     

高超声速风洞变雷诺数试验技术研究

为满足高超声速飞行器气动力雷诺数效应研究需求,在CARDC的Φ1米高超声速风洞中开展了变雷诺数试验技术研究.该项试验技术是利用Φ1米高超声速风洞采用高压下吹-真空抽吸驱动运行方式、风洞运行参数范围宽的特点,通过宽范围内调节风洞运行总压而大幅改变模拟雷诺数.研究采用了单点变雷诺数试验技术和连续变雷诺数试验技术两种手段来开展高超声速飞行器气动力雷诺数效应模拟.单点变雷诺数试验是通过一系列不同雷诺数条件、不同试验车次的试验结果,获得气动特性随雷诺数的变化规律;连续变雷诺数试验时,控制风洞总压从高到低连续变化,测量获取模型处于某一姿态角条件时气动力随雷诺数的变化规律.本文介绍了变雷诺数试验的风洞开车方式、试验及数据处理方法等,并开展了某升力体飞行器和某弹头模型雷诺数效应试验研究.研究结果表明:采用单点和连续变雷诺数试验技术相结合的方式,能较为完整、准确地获得飞行器模型气动力随雷诺数的变化规律.     

高超声速飞行器研制系统工程风险概率分析

针对高超声速飞行器研制工程的高风险特点,对工程决策方和研制方面临的不同类型风险进行了建模分析。基于系统方案或关键技术的固有风险概率和抽象化的研发与验证过程,分别计算分析了工程决策方误判验收通过不合格产品和研制方过度研发或重复验证较低失败概率产品的风险概率,进而提出了高风险研发项目中研制方过度研发风险的概念,明确了工程决策误判风险与研制方过度研发风险的相互影响规律,并基于概率方法建立了一种可以综合权衡决策方风险和研制方风险、合理确定研制周期的系统工程优化方法。     

高超声速边界层的转捩及预测

首先描述了边界层转捩的基本过程及研究内容。在此基础上,指出了高超声速边界层不同于不可压缩边界层的流动不稳定性特性,并介绍了边界层的转捩机理与感受性特征;给出了高超声速三维边界层中预测转捩的常用方法,并着重介绍了多用于工程实际的e N方法以及对e N方法的理性改进,同时列举了在高超声速三维边界层中应用e N方法实现转捩预测的多个实例。最后,分析并总结了高超声速边界层转捩预测所存在的困难及需要解决的问题。     

高超声速数值激波失稳的网格依赖性

高超声速强激波的稳定捕捉仍极具挑战性。目前工程计算中普遍应用的数值格式在模拟多维强激波时通常会遭遇明显的激波不稳定现象,且数值格式的激波稳定性对计算网格表现出严重的依赖性。基于矩阵稳定性分析法,对比了具有不同耗散性质的数值格式稳定捕捉激波的能力,分析了空间二阶精度格式的激波稳定性及限制器对激波稳定性的影响。在此基础上,重点探究了计算网格对激波稳定性的影响规律,研究了网格长宽比和畸变角度对激波稳定性的影响。结果显示,在激波附近通过增大网格长宽比或改变网格畸变角度可有效改善激波捕捉的稳定性;相比于增大网格长宽比,改变网格畸变角度提升激波捕捉稳定性的效果更加明显。在此基础上,结合数值耗散分析对高马赫数下数值激波失稳现象的网格依赖性机制进行了探讨。     

高超弯曲激波二维进气道初步研究

研究了一种压升规律的曲面压缩面,设计了高超弯曲激波二维进气道,并用数值模拟手段对该进气道和同等条件下的常规高超二维三楔进气道、楔+等熵进气道的性能进行了比较.数值模拟表明:通过给定合理压缩面压升规律来设计压缩面并改善压缩面附面层稳定性是可行的,弯曲激波二维进气道的长度比同等条件下的常规二维三楔、楔+等熵进气道分别缩短12%和10%,并且对来流 Ma 变化不敏感,综合性能优势明显、应用前景大.