一种高超声速稀薄流激波干扰气动热测量技术

2021-01-10 03:27王宏宇石义雷龙正义毛春满

宇航学报 2020年12期

王宏宇，王辉，石义雷，龙正义，毛春满，李杰

(中国空气动力研究与发展中心超高速空气动力研究所，绵阳 621000)

0 引言

高超声速飞行器飞行过程中面临严重的气动加热，可导致飞行器结构在高温下发生烧蚀破坏[1]。气动加热在流场的激波干扰区尤为严重，第IV类干扰模式中，干扰区形成的欠膨胀超声速射流导致局部流动参数幅值发生剧烈变化，对飞行器表面传热造成巨大影响[2-3]。

一些学者在预测飞行器激波干扰区表面压力和热流的CFD和DSMC计算校准方面做了大量工作[4-6]，大量研究结果表明，计算结果大多过高地预测了热流值，与实验结果之间存在差异。迄今为止，国际上尚没有形成一个完备的精确预测激波干扰流场的计算方法，也没有找到影响计算精度的确切原因，因此，发展高精度的实验测量技术对于激波干扰气动热特性问题研究十分重要。

国内在高超声速气动热测量研究方面，杨彦广等[7]和彭治雨等[8]对气动热测量技术的发展情况进行了综合分析探讨，并对其发展趋势进行了讨论。针对激波风洞实验的气动热测量，发展了薄膜热电阻和同轴热电偶热流传感器，通过获取传感器感应面温度-时间数据，根据一维半无限体热传导理论计算出模型表面热流[9-10]。徐大军等[11]在激波风洞实验中，获得了吸气式高超声速飞行器气动热环境；张扣立等[12]针对高超声速边界层转捩研究需求发展了温敏漆技术，验证了温敏漆技术对边界层转捩参数定量测量的适用性。

稀薄流域高超声速飞行器的气动加热问题越来越受到重视，精确预测飞行器在稀薄流域的气动加热及其结构在气动加热作用下的温度变化特性能够为飞行器的结构强度计算和热应力计算提供初始计算依据[13-14]。高马赫数稀薄来流条件下飞行器面临的黏性干扰、激波干扰等现象更加显著，目前对激波干扰气动热特性的预测问题还没有得到有效解决，现有的瞬态热流测量技术无法满足对低密度风洞稀薄来流条件下中低热流量值、长时间的精确测量需求，需要发展新的测量技术。为实现低密度风洞来流激波干扰引起的气动热测量，预测热流随不同流态的变化规律，本文发展了适用于长时间、中低量值热流测量的点式测量技术。

1 适用于长时间中/低热流测量的量热计设计

基于半无限体假设的一维热传导模型，设计了带有封装外套的量热计，通过测量量热计表面的温度-时间曲线，计算表面热流。量热计结构示意图如图1所示，为了提高传感器温度测试信噪比，传热体采用导热系数低的康铜材料，外径为0.5 mm的K型铠装热电偶作为前端面测温元件；为降低热电偶测温点接触热阻和提高热响应速度，采用纯银融接工艺将测温接点与传热体前表面粘结在一起；为减小传热体的侧向传热，实现近似一维传热，在康铜传热体外围固定带有锥口的不锈钢封装外套，使外套与康铜传热体之间存有间隙。这种带有空气间隙的隔热封装外套设计明显增大了传热体侧向热阻系数，使侧向传热显著降低，提高了基于一维半无限体假设的热流测试有效时间和精确度。量热计的公称直径为3.2 mm，长度为8 mm。

获得温度-时间数据后，采用基于一维半无限大体热传导理论的Cook-Felderman公式[15]计算出热流：

(1)

为了便于数据处理，将式(1)转化为：

(2)

式中：c为传热介质的比热，单位为J/(kg·k)；k为传热介质的导热系数，单位为W/(m·K)；Tw为测点的壁温，单位为K；ρ为传热介质密度，单位为kg/m3；t为测量时间，单位为s。在本文的研究中，将新型量热计的测量数据与相同测试原理的同轴热电偶测量数据作对比分析，同轴热电偶的公称直径约为2 mm，使用前需要进行标定试验。

图1 量热计示意图

2 地面标定实验

为提高风洞实验的气动热测量精度，量热计在使用前需进行标定，本文采用空气动力研究与发展中心超高速所自主研制的热流传感器标定系统对量热计进行标定，其原理图如图2所示。标定系统由弧光灯光源、位移机构和数据采集系统组成。弧光灯光源用于产生均匀缝补的稳定辐射热流，通过改变弧光灯的输入电流可以产生不同的辐射热流值。入射的辐射热流经过已知吸收率的感应面涂层进入热流传感器中。位移机构调节传感器的位置，使弧光灯标定光源装置的中心对准量热计和戈登计中心。数据采集系统采集温度信号随时间变化，采样频率为1 kHz。将已校准的标准戈登计作为基准传感器测量入射热流，并根据量热计感应面涂层吸收率，计算出量热计净标定入射热流，然后根据式(1)，由量热计所测温度信号计算出热流示值。最后，根据不同弧光灯电流值下的净标定入射热流和热流示值，采用最小二乘优化算法，拟合一条截距为零的直线，其斜率即为量热计的标定系数。

图2 热流传感器标定实验原理图

图3～图6对比了不同辐射热流条件下(光源恒定电流为50 A，70 A，90 A，110 A)量热计和同轴热电偶所测得的温升和热流变化。由于外界气流干扰带来的噪声较小，量热计和同轴热电偶所获得的温升曲线均比较平滑，热流信号的噪声水平较低。其中，量热计的响应时间较慢，大约为1 s。两种测量方式获得的初始热流值较为接近，表明两种传感器均能够对中低热流进行测量。同轴热电偶的响应时间较快，但由于受到侧向传热的影响，其获得的热流信号很快呈现出下降的趋势，且热流值越大，下降越快。增加了封装外套的量热计，由于减小了侧向传热，所测得的热流长时间保持恒定，有效降低了测量误差。标定实验验证了量热计具备对中低量值的热流长时间测量的能力。

图3 量热计温升曲线

图4 同轴热电偶温升曲线

图5 通过量热计温升计算的热流

图6 通过同轴热电偶温升计算的热流

图7给出了量热计热流示值q和戈登计获得的净标定入射热流q*比对数据，由二者线性回归的直线斜率获得量热计的测量数据标定系数k，则量热计的校准值q*为热流示值q与k的乘积。即q*=

图7 量热计热流标定曲线

kq。标定了16个量热计用于风洞实验热流测量。

3 高超声速低密度风洞激波干扰实验

3.1 实验设置

为验证所设计的量热计用于测量高超声速低密度来流条件下模型表面热流的可行性，在中国空气动力研究与发展中心超高速所Φ1 m低密度风洞上开展了双锥模型激波干扰气动热测量实验。模型尺寸采用国际上的标准双锥模型[16]，如图8所示：模型一级锥角和二级锥角分别为25°和55°，总直径为261.85 mm；测量时模型的攻角为0，在模型的两条特征线上设置直径为3.3 mm和2.2 mm的测量孔，分别用于安装量热计和同轴热电偶。模型底部安装法兰盘将模型固定于风洞实验段的快送机构上。当风洞运行流场稳定后，快送机构将模型快速送进至均匀来流的中心位置，送进时间约为0.2 s。模型被送进后，传感器与纹影系统同步采集数据。为对热流分布的规律性进行定性研究，实验采用了高灵敏度的双光程纹影成像系统捕捉流场激波图像，纹影系统配备LED光源，纹影镜，刀口和相机等部件，相机的采样速率为30帧/秒，分辨率为1920×1080，与采用单次平行光通过实验流场的单光程纹影系统相比，双光程纹影系统采用单球面反射镜，利用发散光通过实验流场，由于光束受到两次扰动，光线的偏折角加倍，流场成像的灵敏度更高[17]。

风洞前室总压p0采用测量范围为0～10 MPa绝对压力传感器测量，精度0.2%；前室总温T0采用铂铑-铂B型热电偶测量，测量精度为0.5%。实验采用M10和M12两种喷管，通过典型状态的流场马赫数校测，喷管核心区马赫数分别为9.82(M10)和11.54(M12)。实验过程中，每种工作状态开展三次实验，以考核传感器测量的重复性。

图8 风洞实验模型示意图

3.2 风洞实验结果及分析

图9为实验纹影图，给出了两种工况下双锥模型诱导的激波及其相互干扰状态：包括分离激波、分离区、再附激波和剪切层等典型流场结构，研究表明两种来流工况下模型诱导的流场结构相似。由风洞实验纹影图可知，两种工况下分离点的位置几乎相同，大致位于第一锥体的中心，且激波干涉点的位置也大致相同。由纹影图可预测热流的分布情况，由激波干扰的流动特性可知，位于第一锥体头部附近的边界层黏性干扰区和第二锥体的激波干扰区存在较大的热流。

图9 实验纹影

图10～图13给出了M10和M12工况下量热计和同轴热电偶所测量的模型特征线上的热流分布。由图10～13可知，一级锥的大部分区域热流较低，只有位于上游黏性干扰区的热流较高，自模型拐点起，表面热流呈上升趋势，在激波干扰区再附点处热流达到峰值。再附点下游，热流逐渐减小，出现了一个极低值，而后又逐渐升高，该结果与激波/边界层相互干扰的流动机制完全相符，说明测量结果能够真实反映气动加热过程。三个车次实验结果的对比表明两种传感器测量重复性较好，最大重复性误差小于6%。值得注意的是，同轴热电偶的测量结果较量热计偏大，分析原因，新型量热计的导热体与其封装外壳通过前端刃口配合密封，最大限度减小侧向传热接触面积，也保证了传热体有效受热面积在标定实验和风洞实验中的恒定；而同轴热电偶在风洞实验中，由于模型与安装孔壁面之间存在较明显的缝隙，产生了局部侧向热输入，导致其实际有效换热面积大于同轴热电偶端面的测热面积，从而导致测量结果偏高。因此，同轴热电偶在使用时应与模型测量孔无缝配合以减小测量误差[18]。

图10 M10风洞实验模型热流分布(量热计测量)

图11 M10风洞实验模型热流分布(同轴热电偶测量)

图12 M12风洞实验模型热流分(量热计测量)

图13 M12风洞实验模型热流分布(同轴热电偶测量)

图14～图17给出了量热计和同轴热电偶所测量热流的结果，分别选取低热流值的测点4和高热流值的测点11和测点12进行分析。由图14和图15可知，两种传感器均可以很好捕捉到低热流值(q<10 kW/m2)。由于低热流时侧向传热不明显，同轴热电偶也可以对低热流较好地捕捉，且热流信号的信噪比较量热计高。值得注意的是，两种传感器捕捉的信号中均有一个初始峰值，且同轴热电偶的信号较为明显，这是因为模型送进过程中受到来流的瞬时冲击作用。待流场稳定后，呈现较为均匀的变化。通过地面标定实验不难发现同轴热电偶测量低热流的阈值，大概为q=20 kW/m2。对于较高热流的测点(测点11和测点12)来说，可以明显看到两种传感器的差异性，量热计的响应时间明显变长，测点12经过两秒热流才基本达到恒定，这可能是由于局部流场加热的非定常性引起的。由于同轴热电偶本身的响应时间较快，可以捕捉到振荡较为明显的热流变化，但两侧点的热流值随测试时间延长均显著下降，与地面标定实验的测量结果相似。因此，本文量热计的设计可以满足低密度来流条件下中低量值热流测量，而对于低量值的热流测量，同轴热电偶也能满足测量精度要求，考虑到实验成本，测量的热流较低时，采用同轴热电偶是合适的。

4 结论

本文针对高超声速低密度风洞复杂激波干扰气动热问题的研究，给出了一种长时间、中低量值热流测量方法。采用空气隔热设计的量热计可有效改善量热计侧向传热的影响，实现更高的热流测量精度。通过开展量热计的地面标定实验和M10、M12高超声速低密度风洞激波边界层气动热测量实验，验证了所研制的量热计具备长时间、中低量值热流的测量能力，并且可用于高超声速稀薄来流条件下激波干扰复杂气动热问题的研究。通过本文研究可得出