具有角点的高超音速轴对称钝体头部绕流问题的一种简易解法

钝体头部与后体相接时,往往要出现角点,如果头部的钝度足够大,则在轴对称绕流情况下,音速点即在此角点处。目前求解此种音速点在角点处的钝体头部绕流的直接问题,只有积分关系式法才较精确,但该法须对沿激波层积分得的近似常微分方程组进行数值解,而且所给的二点边界条件中的一点,须按另一点的边界件条用逐步逼近法确定得,整个求解过程很繁琐且费时,而且在角点隣域进行数值积分时困难很大,见文献[7,8]。 鉴于此种情况,本文提出一种简易而又能具有数值解法精确度的实用计算方法,它系以积分关系式法的一级近似常微分方程组为基础,将头部激波用多项式来逼近,利用通过激波的气体貭量守恒关系式,激波的几何关系式以及物面驻点与音速点处的己知条件,直接确定得这种有角点钝体在高超音速轴对称绕流情况下,头部脱体激波的位置,形状以及头部上的压力、速度与速度梯度的分布。在求解过程中全部采用简单的数学分析与代数运算方法,方便且省时。 本文方法适用于任意形状的头部。文中对平头钝体进行了详细计算,从该算例中可见按本文方法所计算得的激波位置系位于一,二级近似的数值解之间,而所得的压力分布则要较一、二级近似的数值解为佳而更接近于实验值。     

轴对称高超音速绕流平头物体的一个解析解

根据方法,使用级数展开和附面层方法,得到了轴对称高超音速绕流平头物体的一个解析解。方法是直接的,即由给定物体形状可求出表面压力分布和激波形状。对平头物体在M_∞=5.8时所得结果和实验及其他理论结果进行了比较,发现符合得很好。本文所提供的方法也说明了对于高超音速钝体绕流级数求介是可能的。     

基于三维弯曲激波的宽域变马赫数乘波体设计分析

为改善常规乘波体布局的宽域气动性能,提出了一种基于复杂三维弯曲激波面的宽域变马赫数乘波体设计方法,其本质是将局部偏转吻切方法与“并联变马赫数”概念相结合,各流面内的设计马赫数不再保持不变,而可根据需求进行调整,以此提高乘波体的宽速域性能。该方法采用Bezier曲面直接指定所需三维激波形状,使得激波形状的选取更为灵活。结果表明,基于三维弯曲激波的变马赫数乘波体拥有更为均衡的气动与几何特性,更适于宽域飞行,且各流面内设计马赫数的离散单调性对此类乘波体性能有较大影响。相同条件下,马赫数7至12范围内设计的基于三维弯曲激波的变马赫数乘波体相较于吻切锥变马赫数乘波体具有更大的容积与更高的升阻比,但容积率有所下降。     

CO2中速度5~7 km/s自由飞圆球流场参数的实验测量和数值模拟

为研究火星进入条件下的非平衡流动特性，在中国空气动力研究与发展中心超高速空气动力研究所弹道靶上测量了CO₂中针对火星探测器进入速度范围5~7 km/s条件下的自由飞圆球的激波脱体距离。实验数据基于阴影法测量，并将其与数值计算结果进行对比，进一步计算了实验流场温度和组分分布等流场参数。一般认为激波脱体距离随来流速度升高而呈单调减小趋势，但研究结果表明:实验状态下，圆球飞行速度约5.5~7.0 km/s的范围内，圆球激波脱体距离随飞行速度升高而增大；采用Park的双温度非平衡模型和5组分6反应的CO₂化学反应动力模型可基本再现本文自由飞圆球激波脱体距离的实验测量数据；根据计算结果推测，本实验状态下自由飞圆球波后靠近激波一侧区域的流场主要处于热化学非平衡状态；当来流速度在约5.5~7.0 km/s的范围内时，流场组分CO开始发生显著离解，是引起圆球激波脱体距离在该速度范围内随速度升高反而增大的可能原因。     

具有气体动力和物理化学非均匀性的超声速物体绕流(英文)

本文研究的问题是超声速气流流过物体时产生的气体动力和物理化学非均匀性及其对物体空气动力特性的影响。就有关流动非均匀性对温度剖面和密度的效应作了详细探讨,还仔细描述了激波与热层相互干扰时以及在激光辐射脉冲的作用下先行激波生成的情况。值得指出的是超声速气流中加入激光辐射能可能改变流动结构,在局部能量释放区域后面呈现的较高和较低压力,有助于控制飞行中物体的运动。可采取某些措施以改进飞行品质。     

环形激波聚焦流场特性的间断有限元方法(英文)

针对环形激波聚焦时产生的高温、高压特性,采用间断有限元方法对环形激波在同轴圆柱形激波管内的聚焦流场进行了数值模拟。为了验证采用方法的有效性,首先对经典激波管问题进行了求解,计算结果与解析解吻合很好。其次,对入射激波马赫数为2.0到5.0的环形激波聚焦流场进行了计算,获得了环形激波在圆柱形激波管内聚焦过程的气动特性。模拟结果表明,采用间断有限元方法能够有效地捕捉激波聚焦过程形成的绕射、反射和聚焦等主要流动特征,而且在聚焦点附近,数值解具有较大的梯度变化,表明该方法对间断解具有较强的捕捉能力。最后比较了网格尺度对聚焦点峰值压力的影响。     

激波在通道口部绕射反射研究

在地下工程研究和设计中,地面运动的高速空气冲击波进入通道后的波动问题是非常重要和急待解决的问题。为了弄清激波进入通道后的传播规律,我们利用新近研制的光学激波管和其他激波管设备,对侧向激波进入通道后的各种模型进行了系统试验研究,观察了激波进入通道后的波动图像,获得了一些很有价值的成果。本文介绍了试验方法和设备,描述了地面激波进入通道后在口部的绕射、反射波动图像和波后多旋涡运动规律,定量给出了通道口部绕射区压力变化曲线。     

高速列车绕流的数值模拟

给出了用低速N－S方程求解高速列车绕流的数值计算方法，采用显隐式结合的方法来加快收敛速度，基于该方法编制了计算软件系统Train3D，用该软件系统对高速列车外形的外流场进行了数值模拟，获得了该外形的空气动力特性和流场结构。     

高超音速钝体小攻角绕流的直接问题

本文建立了高超音速钝体小攻角的绕流积分关系式法的一级近似常微分方程组。确定了钝体对称轴上(s=0)物面与激波处的所有流动参数,使该方程组的边界条件成为一点边值问题。因此对它的求解,要较相应另攻角轴对称绕流的二点边值问题简便。文中还给出了该方程组的一种近似简便解法。     

根据风洞试验结果建立有尾翼导弹数学模型

通过分析有尾翼导弹气动力的对称特性,建立了有尾翼导弹空气动力系数的三角级数模型,给出了模型结构确定方法,并对某导弹风洞试验数据进行了建模研究,建模结果验证了所建模型的有效性和可行性。该套方法可以减少风洞试验量、帮助总结气动力规律、校正风洞试验的误差,并能为控制规律设计和性能分析提供方便。     

亚、超音速细长翼身组合体大迎角气动特性的计算研究

本文提出了一种适于初步设计使用、具有良好精度的亚、超音速细长翼身组合体大迎角气动特性的综合性计算方法。对大迎角情况下的涡升力,采用吸力比拟原理计算;位流升力的计算,采用基本解的数值计算方法。关于机翼翼剖面头部圆度和涡破碎对涡升力的影响,进行经验性修正。翼身干扰的贡献,通过翼身干扰系数进行计算。并按文[4]原理,将亚音速计算方法推广到亚音速前缘的超音速情况。对几种机翼与翼身组合体的计算结果表明,本文方法具有方法简便、计算快速和符合设计精度要求的优点。     

航天飞机在再入段的亚音速空气动力特性的一种有效算法

本文采用吸力比拟原理,结合基本解的数值计算方法,用来计算航天飞机机翼从小迎角到大迎角(a=0°～30°)的亚音速纵向气动特性;而对零升阻力和机身气动特性,则用工程估算方法计算。由于目前的航天飞机,一般为下单翼的复杂外形翼-身组合体,根据文[9]的原理,可忽略翼-身干扰对纵向气动特性的影响。 本文导得可以计及涡效应的任意平面形状边条机翼的亚音速气动特性的计算公式,亦可计算尖梢机翼的展向升力分布。公式中所需的位流系数可采用涡格面元法进行数值计算来获得,压缩性效应则通过位流系数来计及。 本文计算了多种机翼和航天飞机的气动特性。与实验数据比较表明,本方法具有方法简便、计算快速和计算结果具有设计精度的优点,是计算航天飞机亚音速气动特性的一种有效方法。可供航天飞机初步设计使用,亦可作为航天飞机气动优化设计系统中的子系统。经过适当推导,本方法可推广应用于亚音速前缘的超音速情况。     

三维旋转效应对叶片非定常气动特性的影响

在旋转的作用下，流动分离向下游延迟，结果大大改变了风力机叶片的气动特性，这种由旋转所带来的影响被称为三维旋转效应。本文以叶片表面边界层理论分析的结果为基础，给出了三维旋转效应对气流分离影响的解析关系，使得这一复杂问题可以方便地解决，为准确计算旋转叶片的气动栽荷提供基础。通过算例表明，当这一三维旋转效应应用于叶片非定常气动特性计算时，能大大改善计算结果。     

Re数对细长旋成体非对称涡及气动力特性影响的实验研究

通过对细长拱柱旋成体大迎角绕流不同截面测压结果分析,探讨了绕流Re数对非对称涡结构和气动特性的影响,得出Re数不仅影响分离线位置和绕流流态结构,而且影响边界层的绕流特征及其分离涡的强度,非对称性的出现与细长体两侧边界层的绕流特征和分离涡的强度不等存在密切的关系。特别是在同种流态下,两侧边界层的绕流特征和分离涡强度不等是造成侧向力的主要原因;在两侧不同的流态下,转捩不对称是产生大侧向力的主要原因。     

边条机翼气动特性的一种解析计算法

本文以文[1]为基础,研究了平板锐缘边条机翼亚音速气动特性的解析估算,导出了曲线前缘边条机翼气动特性的通用计算式,并可计算尖梢机翼的展向升力分布;所需的位流常数采用涡格面元法来确定。本文就多种机翼进行了计算,与实验结果的比较表明,本方法具有计算简单快速、计算结果具有实用精度的优点,可供初步设计与性能分析时使用。     

尖前缘任意翼型的超音速和高超音速俯仰安定性导数

本文推广了文[1]的欧拉方程摄动法,讨论了尖前缘任意翼型俯仰安定性导数的计算方法,并应用该方法计算了双圆弧翼型的俯仰安定性导数,其数值结果与文[2]的超音速线化位流方程有限差分法的数值结果比较,更接近于实验值。在激波附体的条件下,本文的方法可用于计算大迎角尖前缘任意翼型的超音速和高超音速俯仰安定性导数。     

一种折叠弹翼悬挂物的分离轨迹试验技术

提出了一种新的折叠弹翼悬挂物机弹分离轨迹试验技术,重点解决了在机弹分离过程中折叠弹翼动态展开时悬挂物的气动力获取问题。研究表明,提出的试验技术通过将悬挂物气动力修正方法引入到悬挂物分离安全性研究当中,准确地得到悬挂物的分离特性,解决了折叠弹翼悬挂物分离轨迹风洞试验技术瓶颈,为折叠弹翼悬挂物的投放分离安全性提供一套工程实用的解决方案。     

非烧蚀防热材料地面模拟试验有效性分析

流场化学非平衡度与材料表面催化反应的耦合控制着服役于化学非平衡流场中非烧蚀防热材料表面的气动热载荷,若在该类防热材料性能模拟研究中忽略这一耦合效应,地面模拟试验将无法获得材料的有效使用性能。为此,本文依据钝头体高超声速飞行器边界层驻点热流关系式,分析了影响驻点热流的主要流场参数、地面高焓模拟设备所提供的高焓超声速流场特点以及与飞行热环境之间的主要差异,采用CFD分析了"三参数"模拟方法的有效性。针对化学非平衡边界层驻点传热分析,提出"四参数"模拟方法并分析了"四参数"模拟方法中离解焓无法模拟时的防热材料性能并提出初步解决方法。     

三维高超声速底部喷流干扰流场数值模拟与试验研究

通过数值模拟和风洞试验两种手段对来流马赫数M∞=4、喷流压比Pj/P∞=156.8、不同迎角下的三维高超声速底部喷流干扰流场进行了研究.研究结果表明超声速底部喷流干扰流场结构复杂,有、无喷流时底部流场有很大不同,对气动力系数影响显著;在大喷流压比情况下,喷流干扰使导弹纵向气动力系数下降、压心前移.最后,对数值模拟与风洞试验在结果上的差异进行了分析.     

地面效应下的不同翼型亚声速气动特性分析

采用数值方法研究了亚声速地面效应条件下不同翼型的气动特性,进一步以Ma=0.5来流工况为例,研究了翼型参数和飞行高度对气动特性的影响。计算结果表明在Ma为0.5、迎角为6°的地效情况下,翼型弯度减小,更容易在翼型前缘产生激波阻力;翼型下翼面后缘弯度增大使得后缘压力更高,升力系数和低头力矩相应增大;随着飞行高度的减小,地效作用加强,翼型下翼面压力增大,下翼面的升力增量大于上翼面吸力损失,机翼升力系数和升阻比增加越来越显著。