高超声速飞行器前缘缝隙流动数值模拟研究

通过分析高超声速飞行器前缘防热瓦结构,建立了一种开缝前缘的简化模型。针对这一模型的流场通过求解三维可压缩Navier Stokes方程进行了数值模拟。研究了缝隙诱导形成的三维旋涡的空间分布特征和旋涡运动对物面气动加热的影响规律。模型圆弧段缝隙肩部倒圆区因存在较强的三维效应形成“常规”高热流区,而缝隙内主旋涡再附致使侧壁上存在一个“非常规”高热流区;模型平直段展向流动诱导缝隙上方出现较强的旋涡运动,同时流动在缝隙倒圆区形成分离涡并于缝隙侧壁面再附,受这些旋涡运动的影响,缝隙肩部倒圆区转变为局部热流低值区,缝隙侧壁上存在局部热流高值区。     

MSL火星探测器高超声速流场预测及气动性分析

为获取火星探测器升力-弹道式进入火星大气层的气动特性,以火星科学实验室探测器（MSL）为研究对象,开展基于有效比热比方法的高超声速流动的数值模拟研究,并建立探测器进入过程的气动性能分析方法。计算分析表明,探测器绕流流场因有效比热比较小而具有很薄的激波层,物面气动力与文献值吻合良好;气动分析表明,对特定飞行状态,质心位置决定配平升阻比和俯仰静稳定度,合理的质心法向偏移和纵向移动可满足探测器的升阻特性和静稳定性,同时给定的配平特性和静稳定性要求也可指导质心位置调整,进而为探测器的气动布局设计和舱内设备布置提供建议。     

电热除冰传热特性的结冰风洞实验研究

利用结冰风洞设备和电加热除冰装置,采用实验的方法研究了不同加热模式、冷却时间、加热功率和冰脱落对传热特性的影响.研究表明:设置合理的冷却时间和加热功率,采用高功率的周期性加热模式比采用低功率的连续性加热模式更优越,不仅除冰时间更少,而且能量消耗也更少,从而为电热除冰系统加热模式的选取和传热特性的优化提供了实验依据.     

结冰风洞水滴直径标定方法研究

过冷水滴粒径大小是重要的结冰云雾参数,获知结冰风洞中的水滴直径,是得到定量结冰风洞实验结果的基础.对于结冰风洞内水滴直径单一或者分布比较集中的情况,提出了一种采用数值计算和结冰风洞实验相结合的手段标定水滴直径的方法.该方法首先采用拉格朗日法数值计算水滴运动轨迹,得到撞击极限随水滴直径变化的关系曲线,在此基础上,进行结冰风洞实验,测量实验得到的水滴撞击极限,通过在撞击极限与水滴直径关系曲线上进行插值,进而得到实验水滴直径大小.采用该方法对0.3m×0.2m结冰风洞内的水滴直径进行了标定,分别计算和测量了25m/s和35m/s两种速度条件下的水滴撞击极限,得到的水滴直径值相差不超过1μm,初步说明该方法的合理性.同时,对于结冰风洞内水滴粒径多尺寸分布的情况,还提出了相应的标定其容积平均直径MVD的方法,该方法在计算水滴收集率的基础上,通过测量驻点处的结冰厚度,实现对MVD的测量.采用本文提出的两种方法进行结冰风洞水滴粒径标定,只需要一般的长度测量工具即可进行,操作方便,成本低廉,克服了常规的水滴直径测量或标定需要专门设备的不足.     

热气动弹性变形对飞行器结构温度场的影响研究

气动加热造成的结构温升可能造成飞行器结构失效从而带来安全隐患。准确预测结构温度场在高超声速飞行器热防护系统与结构设计中显得尤为重要。气动热与传热耦合是提高结构温度场预测精度的有效手段,经长期研究与发展,不管是耦合方法研究还是实际工程应用都已开展了大量工作。但这些研究工作均未考虑结构变形对气动加热和结构温度场的影响。而在实际飞行过程中,特别是长时间飞行后,结构变形对结构温度场的影响往往是不能忽略的,对气动力/热环境也都有直接的影响。本文以飞行器静热气动弹性计算方法为基础,对高超声速飞行器机翼模型进行了考虑热气动弹性变形影响的气动热与传热耦合计算,并与不考虑变形对热环境影响情况的计算结果进行了对比分析。结果表明,虽然对于大面积区域变形对气动热/结构温度场的影响较小,但对于热防护结构重点关注且精度要求较高的前缘驻点附近区域计算结果变化明显。由此,也说明了考虑弹性变形对结构温度场预测的重要性。该研究工作为进一步提高飞行器结构温度场预测精度和结构热安全性能评估能力奠定了基础。     

辐射加热对返回舱气动热环境影响的数值研究

基于高超声速再入飞行器气动热环境预测分析的需要,建立了高温非平衡气体辐射加热对飞行器热环境影响的计算分析手段。采用数值求解化学非平衡N-S方程的方法,对返回舱绕流流场进行模拟,获得高温空气组分质量分数和温度等流场参数分布。基于辐射传输方程,考虑高温气体组分的主要辐射机制,计算分析高温流场气体辐射加热对返回舱热环境的影响。分析表明,在同一飞行弹道条件下,返回舱大底半径尺寸对气动加热的影响较大,在再入热环境严酷区,辐射加热对物面总热流的贡献达30!;产生辐射加热效应的主要机制是高温流场中O和N原子产生连续谱和线状谱以及N2的第一正带系;物面催化效应对辐射加热影响不大。     

基于轴对称比拟的高超声速飞行器表面热环境数值与工程耦合算法研究

基于轴对称比拟概念,发展了一套笛卡尔坐标系下的高超声速飞行器气动热环境计算方法。首先利用有限体积法数值求解Euler方程获得较为准确的边界层外缘无粘流场参数,然后基于有限元的四节点单元变换方法,直接利用笛卡尔坐标系下的三维速度分量计算无粘表面流线和尺度因子;在获取无粘表面流线和尺度因子的基础上,利用Zoby、Moss和Sutton提出的热流公式计算表面热流,从而实现数值算法和工程算法的耦合。将上述方法用于求解球锥在攻角分别为0°、8°和16°时的表面热流,并将计算结果同经典流线法结果及实验值进行比较,从而对方法进行考核验证。结果显示:本文计算结果与实验值吻合的较好,计算精度较经典流线法有较大提高。     

利用二维平面和轴对称逆向喷流减阻和降低热流的计算研究

对近年来国内外学者在减阻和降低热流方面的工作做了简要回顾,并对带逆向喷流钝体的高超声速绕流流场进行了数值计算,对二维圆柱和球头物形的差异给喷流减阻和降低热流效能带来的影响进行了研究。结果表明逆向喷流对二维平面物形具有良好的减阻和降低热流作用,但对于轴对称的球头由于肩部流动再附处的强烈激波干扰,会形成热斑,其局部热流反而高于无喷流时的热流值,适当加大喷流压比可以避免激波在肩部再附区的干扰,达到降低表面热流的目的。本文认为在逆向喷流的研究中,采用真实维数的物形对于得到正确的研究结论是非常重要的。     

结冰风洞水滴直径标定方法研究

过冷水滴粒径大小是重要的结冰云雾参数,获知结冰风洞中的水滴直径,是得到定量结冰风洞实验结果的基础。对于结冰风洞内水滴直径单一或者分布比较集中的情况,提出了一种采用数值计算和结冰风洞实验相结合的手段标定水滴直径的方法。该方法首先采用拉格朗日法数值计算水滴运动轨迹,得到撞击极限随水滴直径变化的关系曲线,在此基础上,进行结冰风洞实验,测量实验得到的水滴撞击极限,通过在撞击极限与水滴直径关系曲线上进行插值,进而得到实验水滴直径大小。采用该方法对0.3m×0.2m结冰风洞内的水滴直径进行了标定,分别计算和测量了25m/s和35m/s两种速度条件下的水滴撞击极限,得到的水滴直径值相差不超过1μm,初步说明该方法的合理性。同时,对于结冰风洞内水滴粒径多尺寸分布的情况,还提出了相应的标定其容积平均直径MVD的方法,该方法在计算水滴收集率的基础上,通过测量驻点处的结冰厚度,实现对MVD的测量。采用本文提出的两种方法进行结冰风洞水滴粒径标定,只需要一般的长度测量工具即可进行,操作方便,成本低廉,克服了常规的水滴直径测量或标定需要专门设备的不足。     

导弹投放分离的数值模拟

采用非结构网格、有限体积方法,利用准定常假设,对两种飞机在不同状态下与导弹的分离进行了数值模拟,并将计算结果和试验进行了比较。数值模拟的结果和试验结果表明,对于升力系数和俯仰力矩系数随攻角的变化、导弹轨迹随时间的变化两者吻合较好。