广义Burgers方程声爆传播模型高阶格式离散

声爆预测技术是研制新一代环保型超声速民机的关键技术之一。针对当前基于广义Burgers方程的声爆远场传播模型的数值求解格式精度较低的情况,开展了高阶格式离散方法研究。通过分析该模型中各效应项的物理性质,应用合适的高阶精度格式对各项分别进行离散求解,并采用经典几何声学射线法计算声爆传播路径,实现了对地面声爆波形的精准预测。通过对美国F-5E机型的声爆飞行试验和第二届国际声爆预测研讨会的典型算例进行数值模拟,验证了本文高阶离散方法的可靠性。进一步的结果分析表明,采用高阶离散格式的预测方法相比传统二阶精度预测方法具有明显的低耗散特性,在同等网格数目上可获得更高分辨率的计算结果,且其网格收敛性远优于传统二阶精度格式,同时在计算效率上也有一定优势;在Burgers方程的各效应项中,非线性项的影响更为明显,因而采用高阶离散格式时其优势更为突出;热黏性吸收项对数值结果影响较小,采用高阶离散格式对计算精度的提升效果也不明显,实际计算中仍可以采用传统离散格式,甚至忽略该项的贡献。     

火星大气高温光谱建模与非平衡辐射热流预测

火星科学实验室成功着陆后的热环境重构数据表明,气动辐射加热在火星进入防热设计中具有不同于以往认识的重要影响,未知机制和模型不确定性等问题有待进一步研究。探测器高速进入火星大气产生极高温非平衡气动环境,造成火星气动辐射与常规CO₂红外辐射研究显著不同。针对火星大气高温光谱和辐射热流预测,首先,建立适用于火星大气的高温非平衡光谱辐射模型,获得典型高温条件CO₂光谱结构和辐射强度,与NASA和JAXA试验结果对比,结果显示符合较好。其次,依靠激波管和发射光谱测量技术,开展典型进入条件的辐射强度测量试验,数值结果、试验值和NASA试验结果相互符合较好,验证了光谱模型。最后,对探路者号进行气动辐射加热分析,完成典型进入条件下的非平衡流动和辐射特性计算,基于光线法得到光谱辐射强度沿驻点线的变化,表明高速与低速条件的气动辐射机制存在显著差异;基于有限体积法获得进入器表面辐射热流分布,结果显示辐射热流的分布及变化规律与地球再入显著不同,进入速度6 km/s以下时辐射热流随进入速度增加而减小,同时进入器锥身及肩部的辐射热流高于驻点区域。     

SABRE预冷器计算模型及其在整机模型中的应用

预冷器的性能对SABRE等预冷组合循环发动机具有重要影响,为实现发动机方案设计阶段预冷器的快速设计与评估,建立了预冷器准二维快速评估模型。将SABRE预冷器的几何结构简化为一个扇环形区域,沿径向和周向将该区域划分为二维节点。应用守恒方程及传热关联式完成单个节点计算,再求解节点矩阵的平衡方程组,计算内外流体特定节点上参数的二维分布,得到预冷器出口参数。将模型嵌入发动机总体性能程序中,实现了发动机设计及非设计点的预冷器性能计算的功能。与文献数据对比结果表明,预冷器模型传热计算误差小于5%,摩擦阻力误差小于10%。整机计算结果显示,Ma0~5范围内,预冷器空气侧温降范围为143K至932K,温降随飞行马赫数升高单调上升。预冷器传热有效度范围为0.896-0.945,空气侧总压恢复系数范围为0.852-0.904。     

飞行器太阳能源计算方法CSCD

介绍了空间飞行器典型电源系统组成,通过对影响太阳能源输出的各项参数因素分析,提出了空间飞行器太阳能源计算方法,该方法包括由功率计算所需太阳阵面积的公式和由太阳阵面积计算输出功率的公式,达到了准确计算太阳能源的效果,解决了以往空间飞行器太阳能源估算误差大,在总体回路设计中电源参数难以闭环的问题。     

前体尾流对降落伞工作性能的影响

为研究前体尾流对降落伞工作性能的非定常影响,基于Realizablek-ε湍流模型采用PISO算法开展了物伞系统的非定常绕流数值计算,获得了精细的流场旋涡结构。在此基础上,研究了不同拖曳比下物伞系统的尾涡演变规律、流场分布规律以及伞衣气动特性变化。结果表明:前体尾涡导致伞衣入口处的涡量大小和方向时刻变化,随拖曳比增加,涡量黏性耗散增强,进入伞衣的旋涡强度逐渐减弱,伞衣入口形成稳定的负涡量区,伞衣尾涡脱离周期随之延长;拖曳比对尾涡区后端（伞衣入口处）流场压力的影响远大于前端,随拖曳比增加,流动形式逐渐由闭式转变为开式,流场的速度分布和压力分布更为对称,伞衣入口形成稳定的正压区,内外压差增加;当拖曳比大于9时,前体尾流对降落伞阻力系数和表面压强系数的影响减小。      

矩形断面非线性驰振自激力测量及间接验证中若干重要问题的讨论

采用专门研制的小型动态测力天平,通过弹簧悬挂节段模型内置天平同步测力测振风洞试验,对3:2矩形断面的非线性驰振自激力进行了测量。比较了基于实测自激力重构的节段模型位移响应时程与试验结果,从而对自激力测量精度进行了间接的验证;讨论了在进行这种验证时考虑节段模型系统等效阻尼和刚度参数非线性特性的重要性、非风致附加自激力和风致自激力在动态力中的占比、忽略非风致附加气动阻尼力和惯性力的非线性特性对自激力测量精度的影响等若干重要问题。结果显示:对于3:2矩形断面,非风致附加自激力在测得的总动态力中的占比超过了风致自激力的占比,因此从测得的总动态力中提取自激力时必须扣除非风致附加自激力;非风致附加气动阻尼力和惯性力的非线性对驰振自激力测量精度有一定影响,值得考虑;节段模型系统等效阻尼和刚度参数的非线性对节段模型驰振位移响应的重构精度有明显影响,在验证自激力测量精度时必须加以考虑。     

探测器安装结构对航天器壁温测量的影响分析

针对飞行试验中飞行器薄壁壳体测量温度与预测温度存在较大差异这一问题,采用气动热工程算法结合热传导计算方法,分析了测温探测器安装结构对测点温度的影响,并提出了改进措施。结果表明:对于薄壁结构飞行器在上升段有气动加热、其表面处于升温过程或热量由壳体表面向内部传导时,测温探测器安装结构对测点温度基本无影响。但当飞行器处于飞行中段,在辐射散热、表面温度低于壳体内部温度造成热量由壳体内部向外表面传导时,测点温度受原探测器安装结构影响明显,测量温度明显低于不装探测器时的预测温度;而采用本文提出的探测器安装方案,可明显降低对测点温度的影响,在飞行器的测点位置最大影响小于0.5K。     

附属构件及桥面粗糙度对近流线型宽体箱梁气动静力系数的影响

针对高宽比为1/12的近流线型宽体箱梁进行了一系列参数化的三分力测力试验,研究了栏杆透风率、检修车轨道位置、桥面粗糙度对三分力系数的影响.运用CFD工程软件Fluent 15.0模拟二维风场,对三分力系数变化的可能机理进行讨论.结果表明,对于近流线型宽体箱梁:(1)随迎角由负到正加大,CD先降后升,CL与Cm总体趋于升高,但均不关于零迎角对称或反对称;(2)栏杆透风率的减小一方面加大了断面挡风面积,使尾流变宽,CD升高;另一方面扩大了断面上方负压区范围(正迎角工况),提高了上表面压力系数(负迎角工况),使正、负迎角下的升力均增加,即CL绝对值升高;(3)正迎角来流时,检修车轨道的内移使断面下方的负压区缩小、压力系数减小、尾流变窄,则CD、CL均降低;但负迎角来流时,其位置的移动未对分离区造成明显影响,则CD与CL均未出现显著波动;(4)正迎角来流时,桥面粗糙度的提高减小了上表面的分离角,使负压区缩小、尾流变窄,则CD降低;且负压区的缩小使得朝上的升力减小,即CL降低.负迎角来流时,粗糙度的提高使上表面的阻力与压力系数均增加,则CD升高、CL的绝对值升高;(5)栏杆透风率、检修车轨道位置及桥面粗糙度对Cm的影响可以忽略不计.     

尾身组合体不同类型栅格翼气动特性试验研究

对于不同类型的栅格翼,采取了一种模型设计简洁、克服了在较小试验模型上安装铰链力矩天平的困难,同时又能测量栅格翼全部气动力和力矩分量的试验方法,获得了简单框架式、正置蜂窝式以及不同网格密度的斜置蜂窝式4种栅格翼在弹身影响下的亚声速和跨声速气动特性,并与平板尾翼的气动特性进行了对比,同时研究了栅格翼在尾身组合体迎风侧和背风侧不同位置时的气动特性。结果表明:栅格翼与常规平板翼各有优缺点,栅格翼结构形式、网格密度和尺寸对其气动特性有重要影响,位置对不同类型栅格翼气动特性的影响也不尽相同。     

卫星介质深层充电电场时域差分计算方法研究

根据卫星介质深层充电的物理机制,分别对应介质板正面接地、背面接地和正面背面同时接地三种情况建立相应的介质板深层充电模型,研究了利用时域差分和数值积分方法求解卫星介质深层充电电场的算法,计算了在单能电子束流垂直入射作用下卫星内部介质板(Teflon)深层充电电场,计算结果表明,文章建立的数值计算方法能求解出介质深层充电电场,且计算结果能反映充电电场在介质板中的时空分布特征.