表面热流可辨识性初步分析

表面热流的可辨识性分析可用于飞行器防热层内温度测量精度和测温点位置的确定,在工程上有较强的实用意义。从无量纲分析和仿真辨识出发,根据防热层材料热物性系数、测点位置、表面热流的频域特性等参数对表面热流辨识结果的影响规律,总结出了表面热流辨识问题的相似参数：基于表面热流频率参数的傅立叶数。此后,以这一傅立叶数为判据,针对不同测量误差值的情况初步建立起了表面热流可辨识性的准则和分析方法。     

高超声速气动热辨识技术实验研究

为估计高超声速飞行中的气动热流,采用热流辨识技术研究了平板表面的气动热流。首先采用仿真数据对辨识方法进行验证,结果表明,在测量噪声较低时热流辨识结果误差较小。此后,开展了平板的高超声速气动热风洞实验,通过内埋热电偶测量平板内部温度,采用热流辨识法获得平板表面的气动热流,采用辅助测点对比验证了内埋测点辨识结果的可信度。最后将辨识结果与平板边界层热流理论估算结果进行对比,分析了两者产生误差的原因。本文的研究为热流辨识技术的实际应用奠定了基础。     

MF-1飞行试验表面热流辨识

利用物体内部测点温度辨识表面热流,是一典型热传导逆问题。利用一维表面热流辨识方法处理分析MF-1航天模型飞行试验的温度测量数据,获得了表面热流数据。辨识结果结合数值模拟预测结果可以有效地判断边界层流态变化(再层流化、转捩)发生的区域,达到了预期目标。一维辨识方法模型过于简单有其局限性,比如一维方法没有考虑到测温结构与飞行器壳体之间的横向传热,辨识结果中出现了不合理的负热流,精度不足以支持深入定量对比分析。为了克服这一缺陷,提高辨识精度,采用三维传热模型来分析处理温度测量数据、辨识表面热流,由此得到的辨识结果比一维辨识结果更加合理。在三维传热模型的表面热流辨识计算中,还考虑了接触热阻对测温部件表面温度辨识结果的影响,并通过模拟计算给出了他们的影响规律。通过这一系列分析可以看到,飞行试验中测温部件传热较为复杂,而基于精细传热模型的表面热流辨识方法可以有效获得精度较高的辨识结果。同时这也预示了这一方法在航天飞行试验中应用的潜力。     

压缩性修正对γ-Re_θ转捩模型的影响研究北大核心CSCD

将基于流场当地变量的γ-Reθ转捩模型加入可压缩RANS方程计算程序,对平均流场控制方程、湍流和转捩控制方程进行基于LU-SGS的隐式紧耦合求解。结合三种压缩性修正方法对高超声速平板、双楔和圆锥流动进行了数值模拟,给出了壁面斯坦顿数、热流值和实验值的比较,以及湍动能和间歇因子分布云图。数值计算结果表明,相同来流湍流度下不同压缩性修正方法得到的转捩位置差别较大。相比于未修正的模型,基于当地马赫数压缩性修正后的转捩位置最靠后,其对带有分离的双楔流动预测较为准确;基于来流马赫数的压缩性修正仅使得转捩位置稍有延迟;对湍动能源项的压缩性修正可提高模型在转捩后湍流段的热流预测精度。     

压缩性修正对γ-Re_θ转捩模型的影响研究

将基于流场当地变量的γ-Reθ转捩模型加入可压缩RANS方程计算程序,对平均流场控制方程、湍流和转捩控制方程进行基于LU-SGS的隐式紧耦合求解。结合三种压缩性修正方法对高超声速平板、双楔和圆锥流动进行了数值模拟,给出了壁面斯坦顿数、热流值和实验值的比较,以及湍动能和间歇因子分布云图。数值计算结果表明,相同来流湍流度下不同压缩性修正方法得到的转捩位置差别较大。相比于未修正的模型,基于当地马赫数压缩性修正后的转捩位置最靠后,其对带有分离的双楔流动预测较为准确;基于来流马赫数的压缩性修正仅使得转捩位置稍有延迟;对湍动能源项的压缩性修正可提高模型在转捩后湍流段的热流预测精度。     

飞行试验热流辨识和边界层转捩滞后现象

对中国空气动力研究与发展中心马赫数为5左右的球锥模型在首次航天模型飞行试验中的温度测量数据进行了分析,通过辨识获得热流分布,发现飞行试验的测热数据后处理方法与地面风洞试验有很大差别,必须考虑温度变化历史,并考虑测温单元与周围飞行器壳体的三维传热才能得到正确的热流结果。采用工程计算方法对模型表面热流分布进行了计算,通过与飞行试验测量结果对比分析,发现测温点在发射上升段由湍流完全变为层流和在再入下降段由层流向湍流转捩具有不同的转捩准则数,边界层转捩存在滞后现象;根据地面风洞试验拟合出的转捩准则受到风洞噪声等因素的影响,预测的转捩位置比实际情况靠前;对于球钝锥外形,当x/R50时,流场和热流趋于锥形流结果。本次模型飞行试验还首次验证了气动热工程方法对于马赫数小于5情况的适应性。     

考虑烧蚀情况下的表面热流辨识

针对烧蚀传热问题,在热解面模型的基础上通过伴随方程推导建立了基于测点温度辨识表面热流的方法,并进行了算例考核。结果表明:热流辨识结果与真值符合较好,辨识结果与真值之间的偏差随测量误差的增加而增加;烧蚀后退量测量结果的误差对辨识结果有较为显著的影响。然后,将该辨识方法用于钝头型碳酚醛材料Narmco4028试件在陶瓷加热风洞中的烧蚀试验结果分析,结果表明辨识出的表面热流与加热功率基本符合,辨识方法是有效的,在工程实际中有较好的应用前景。     

MF-1模型飞行试验表面压力与温度测量技术研究

MF-1是我国首次以高超声速空气动力学基础问题研究为目的的航天模型飞行试验,试验模型为锥-柱-裙体,主要研究0°迎角圆锥边界层转捩和压缩拐角激波/边界层干扰现象。针对飞行试验转捩区测量需求,引入和改进了风洞试验中常用薄壁测热技术,设计了一种新型变厚度薄壁测温结构,有效抑制了侧向导热损失,可基于一维热流辨识方法获取可靠的表面热流数据;与现有风洞试验薄壁测热技术相比,该方法可提高有效测量时间,降低时间延迟效应,适于长时间飞行试验测量。针对柱-裙压缩拐角激波/边界层干扰区压力测量需求,采用了风洞试验中常用的基于引压管和电子压力扫描阀的测量方案,通过改进装配工艺,提高了系统耐压能力,实现了模型飞行试验全弹道表面压力测量。模型飞行试验结果表明:MF-1模型飞行试验测量系统可靠,获得了可供边界层转捩和激波/边界层干扰研究分析及CFD验证的可信数据;在热流急剧下降时一维热流辨识存在较大误差,以及压力测量中的时间延迟和低压测量准确度存在不足,是需要进一步改进的问题。     

固定转捩和自由转捩流动预测方法验证

采用固定转捩位置和基于SST k-ω湍流模型的γ-Reθt转捩模型自由转捩模拟两种方法研究了RAE2822翼型流动转捩问题。先采用γ-Reθt转捩模型确定转捩位置,再采用所获得的转捩位置实施固定转捩位置计算,并与全湍流计算结果和其他软件计算结果进行了对比。表明两种方法均能有效捕捉超临界翼型流动转捩现象,固定转捩位置模拟结果基本上重现了自由转捩预测结果。     

后掠机翼人工转捩最佳粗糙带高度数值预测

介绍了基于当地变量的γ-Re_θ转捩模型,并将该模型应用到后掠机翼的转捩预测和人工转捩最佳粗糙带高度以及人工转捩技术能够模拟的大气飞行雷诺数的确定中。为检验γ-Re_θ转捩模型对后掠机翼转捩的预测能力,对ONERA M6机翼和DLR-F4标模机翼进行了边界层转捩预测,采用结构化网格和有限体积法求解雷诺平均Navier-Stokes(RANS)方程,得到了机翼表面的摩擦阻力系数分布,从而可以得到相应的转捩位置,预测得到的转捩位置与试验结果比较吻合,说明该模型对后掠机翼转捩预测是可信的。最后在DLR-F4标模机翼上表面固定了粗糙带,通过相同的方法得到了转捩位置,从而确定了马赫数为0.785、雷诺数为3.0×10⁶时最佳粗糙带高度为0.11 mm;通过不断增大雷诺数使自由转捩位置不断向前缘移动,验证了人工转捩对大气飞行雷诺数的模拟能力。结果表明,在最佳粗糙带高度为0.11 mm下,可以实现对大气飞行高雷诺数的模拟。     

变状态测量中薄壁量热计后壁导热损失修正

针对超燃冲压发动机脉冲燃烧试验需求,基于能量平衡原理给出了对变状态条件下修正薄壁量热计导热损失的一种方法.对薄壁量热计进行了数值分析和试验验证,对比分析了修正前后的薄壁量热计测量结果.有限元分析结果表明,该方法可以有效补偿后壁导热损失,算例中修正前的热流测量误差不小于20%,修正后不大于4%.在超燃冲压发动机脉冲燃烧试验中,利用薄壁量热计与同轴热电偶热流测量了轴对称模型壁面上相同轴向位置的热流,冷态试验中修正前后的薄壁热流与同轴热电偶热流测量结果的相对偏差分别为-11.2%和3.3%.热态试验结果表明,试验中存在振荡燃烧现象,薄壁量热计的响应时间约为0.02s,与振荡周期(约为0.03s)接近,不能正确反映燃烧状态的变化.      

“Ⅳ型”激波干扰中流-热-固耦合问题一体化计算分析

针对高超声速进气道前缘"Ⅳ型"激波干扰产生的气动加热与结构传热多物理场耦合计算问题，发展了一种基于有限体积法的流-热-固一体化计算方法。该方法采用一体化控制方程组统一离散求解外部高速流场与内部结构温度场，规避了传统分区耦合算法在时间域内交替迭代的繁琐数据交换策略。另外，提出一种新的双温阻模型计算流-固交界面的物性参数以保证计算准确性，采用LU-SGS隐式时间迭代和自适应时间步长以提高计算效率。采用经典高超声速二维圆管流-热-固耦合算例对该一体化方法进行验证，计算结果与试验值和参考文献数据吻合较好，证明了该方法的可靠性和正确性。利用一体化方法对高超声速前缘"Ⅳ型"激波干扰流-热-固耦合问题进行定常/非定常计算与分析，给出了温度与热流的时变特性，计算结果表明，激波干扰作用产生的超声速"喷流"不断冲击壁面，使得壁面最大压力系数增大约9倍，壁面最大热流增大约4.7倍，给高速飞行器的热防护设计与选材带来严峻挑战。同时，也表明了一体化计算方法可以较好地用于长航时飞行条件下与复杂飞行环境下的高超声速热防护系统的热环境特性分析与综合性能评估。     

高超声速飞行器表面测热技术综述

高超声速飞行中飞行器表面气动加热量是飞行器热防护系统最为关键的设计输入,在理论计算与地面模拟的有限近似条件下,通过飞行试验实时获取真实环境下飞行器表面的气动加热量,并在此基础上完成对计算模型与地面试验的验证与改进具有重要意义。详细列举了20世纪50年代以来,国外具有代表性的飞行试验与测热方案。以“内置式”与“嵌入式”作为测热技术的分类特征,介绍了各类测试设备及相应的飞行试验结果。着重分析了“热匹配性”与“结构匹配性”作为关键因素对飞行测热技术的影响,通过飞行试验实例介绍了该问题的解决方法及工程经验。归纳了飞行测热技术发展的共性、特点与未来趋势,并结合当前我国发展现状,对该领域未来研究提出建议。     

典型尖楔结构燃气流热试验工况确定方法

在传统冷壁热流模拟方法的基础上,进一步提出以热壁温度及热流密度的时序变化曲线为控制目标的燃气流热试验工况确定方法,即利用壁温控制目标与实测值的偏差对热壁热流控制目标做一定修正,以尽可能消除和弥补前期试验误差,同时利用300K冷壁边界热流密度数据库插值迭代方法,快速确定一定气动热模拟所需燃气流温度,解决了沿飞行轨迹瞬态热试验技术难题之一。利用CFD数值模拟方法,建立了典型尖楔结构高/中温双路燃气流组合热试验300K冷壁边界热流密度数据库,并针对典型尖楔结构沿某飞行轨迹9个典型状态气动热模拟需求,确定相应双路燃气流热模拟参数。相关数值计算结果显示,驻点区域热流密度平均模拟偏差为4.5%,平板区热流密度平均模拟偏差为4.6%,两者最大模拟偏差均不大于8%,满足工程试验精度要求。同时,瞬态热分析结果显示第45s时,距驻点1mm处最大温度梯度达到21K/mm,距驻点10.1mm处最大温度梯度达到18K/mm,满足气动热大温度梯度效应需求。      

再入飞行器鼻锥逆向喷流对流场及气动热的影响

 使用计算流体力学（CFD）方法研究逆向喷流热防护系统对降低再入飞行器鼻锥物面热流的效果,获得了流场参数,回流再附点位置,物面压力分布以及热流分布。分析了逆向喷流对降低物面热流的物理机理,喷流通过与来流相互作用形成马赫盘,将来流导流到四周,不与物面直接作用形成气动加热,同时喷流回流形成低温区,降低物面与接触气体的温差,进而降低了物面热流。随着总压比率增大,这种效果越明显,气动加热越轻。为更合理分析喷流强度对流场及传热量的影响,将总压比率和流量相结合,提出了新的参数R PA。分析该参数的应用效果,结果发现不同的流量与总压比率组合成相同的参数R PA,可以实现相同的激波位置、再附点位置、表面热流峰值位置和总传热量。这说明该参数可用于表征喷流强度,用以分析喷流对流场及传热量的影响。     

大面积气动加热的石英灯阵模拟优化设计

针对结构热试验模拟的精细化需求,对大面积气动加热的石英灯阵模拟优化设计进行了研究。发展了石英灯阵辐射热流模拟程序,分析了采用传统单灯热流分布数据库插值叠加获得的石英灯阵热流分布的适用范围,基于遗传算法,发展了以灯阵中各灯功率为优化参数的石英灯阵热流模拟优化设计方法,并基于所建方法对某飞行器结构部件迎风面气动加热进行了灯阵模拟,获得了灯阵加热和气动加热条件下迎风面温度变化特性。结果表明,基于本文方法对石英灯阵中各灯功率进行优化设计,采用简单石英灯阵即可有效地模拟大面积非均匀气动加热,从而有效提高试验模拟精度,但前缘等位置的高热流区模拟精度有待进一步提高。     

高超声速圆锥边界层失稳条纹结构实验研究

边界层转捩的准确预测是高超声速飞行面临的关键气动问题之一。为研究高超声速边界层失稳和转捩机理,以前缘半径1.6mm、半锥角7°的圆锥模型为研究对象,在FD-07高超声速风洞中采用红外热图技术开展边界层转捩实验测量。通过与工程计算结果对比,确认模型表面边界层流态。实验结果表明:有迎角条件下,模型表面中后段出现条纹结构,条纹结构的起始位置随着周向角的增加而向上游移动;随着迎角的增加,条纹起始位置向上游移动,条纹强度差异和条纹与模型中心线的夹角越来越大。实验获得的条纹结构与不同频率扰动波相互作用直接数值模拟获得的条纹结构现象一致。通过对比分析,认为边界层内不同频率扰动波相互作用是产生条纹结构的一种机制。     

狭窄空间内沸腾的几个特征

将石英或瓷质板用电热膜技术加工成加热板。将两块加热板对称地放入饱合氟利昂－１１３（Ｒ－１１３）中形成竖直窄通道。对通道内沸腾的始沸点、流型和传热规律等基本特性进行实验研究。结果表明：始沸热流随加热板通道间距减小而降低;沸腾流型与加热过程、通道间距等因素有关;存在最佳传热通道间距;最后提出小通道内始沸热流、传热规律的半经验关系式。     

Investigation of thermal protection system by forward-facing cavity and opposing jet combinatorial configuration

 This paper focuses on the usage of the forward-facing cavity and opposing jet combinatorial configuration as the thermal protection system (TPS) for hypersonic vehicles. A hemispherecone nose-tip with the combinatorial configuration is investigated numerically in hypersonic free stream. Some numerical results are validated by experiments. The flow field parameters, aerodynamic force and surface heat flux distribution are obtained. The influence of the opposing jet stagnation pressure on cooling efficiency of the combinatorial TPS is discussed. The detailed numerical results show that the aerodynamic heating is reduced remarkably by the combinatorial system. The recirculation region plays a pivotal role for the reduction of heat flux. The larger the stagnation pressure of opposing jet is, the more the heating reduction is. This kind of combinatorial system is suitable to be the TPS for the high-speed vehicles which need long-range and long time flight.