罗跃^,2), 王磊, 党雷宁, 刘进博, 张军, 柳森^,3)中国空气动力研究与发展中心超高速空气动力研究所,四川绵阳 621000

ARCJET ABLATION EXPERIMENT TO SIMULATE THE CHELYABINSK ASTEROID ENTRY ¹⁾

Luo Yue^,2), Wang Lei, Dang Leining, Liu Jinbo, Zhang Jun, Liu Sen^,3)Hypervelocity Aerodynamics Institute of China Aerodynamics Research and Development Center,Mianyang 621000,Sichuan,China

通讯作者: 2)罗跃,高级工程师,主要研究方向：高超声速飞行器热防护试验技术. E-mail:myheater@sohu.com;3)柳森,研究员,主要研究方向：超高速空气动力学和超高速碰撞动力学. E-mail:liusen@cardc.cn

收稿日期:2020-03-11接受日期:2020-03-11网络出版日期:2020-09-18

基金资助:

1)国家自然科学基金资助项目.11972358 国家自然科学基金资助项目.11802321

Received:2020-03-11Accepted:2020-03-11Online:2020-09-18
作者简介 About authors


摘要
烧蚀是小行星极高速进入地球大气层后最重要的现象之一,在很大程度上决定了小行星的质量/尺寸变化、飞行轨迹、甚至光辐射特性. 为观测小行星材料在超高速高温流场中的烧蚀现象,在电弧加热器上开展了模拟Chelyabinsk小行星事件典型弹道状态(速度约5.6 km/s,高度17 km,流星体直径1 m)的烧蚀实验. 试件为钝头外形,头部半径20 mm,半锥角18$^\circ$. 作为对比,试件分别采用玄武岩和碳钢材料. 成功记录了清晰的烧蚀动态过程,观察到两种材料试件表面的熔融损失流动、以及玄武岩试件的蒸发喷射和崩裂剥落等现象,全程测得烧蚀气体发射光谱、试件实时外形变化、表面热图变化等数据. 分析结果显示了两种材料明显不同的烧蚀现象和质量损失机制：碳钢在高温气流冲击作用下溅射成大量微小液滴,跟随气流高速流失;玄武岩质量损失以熔融物剪切流动为主,并伴随少量块状剥落及蒸发喷射. 烧蚀时间为4 s,玄武岩和碳钢的质量损失及驻点后退量分别为37.9 g,72.7 g以及7.3 mm,13.1 mm,估算玄武岩材料的有效烧蚀焓约为2.6 MJ/kg,两种材料的烧蚀光谱测量组分与电镜能谱扫描结果吻合.
关键词： 小行星;进入;烧蚀;电弧加热器;实验

Abstract
Ablation is one of the most important phenomenon when an asteroid enters the earth atmosphere at hypervelocity, which largely determines the mass loss, flight trajectory, and even radiation characteristics of the asteroid. To research the typical ablation process of asteroids when entering the earth atmosphere, experiments were conducted in an arc heater to simulate the typical conditions (velocity: 6 km/s, height: 17 km, diameter: 1 m) of Chelyabinsk asteroid event. The blunt-shaped test samples with the head radius of 20mm were made by carbon steel and basalt. In this work, the ablation process of test samples were clearly recorded, in which the melt flow of two different materials and the spallation of fragments as well as vaporization of basalt were observed. The evolutions of emission spectroscopy, recession profile and surface temperature profile during the whole process were acquired. The results indicate that the ablation phenomenon and the mechanism of mass loss of two materials are obviously different: Under the impact of the high-temperature flow, the carbon steel was sputtering into mass of tiny droplets which were washed away by the flow rapidly, while the mass loss of basalt were the shear flow of molten matter with small amount of massive spalling and evaporation spraying. All samples were exposed for four seconds in the plasma stream, mass loss and recession of the basalt and carbon steel were [37.9 g, 7.3 mm] and [72.7 g, 13.1 mm] respectively. The estimated effective enthalpy of ablation was 2.6 MJ/kg, the component measured by emission spectroscopy conforms to the electron microscopy (EDS) scanning.
Keywords：asteroid;entry;ablation;arc-heater;experiment


PDF (21285KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
罗跃, 王磊, 党雷宁, 刘进博, 张军, 柳森. 模拟Chelyabinsk小行星进入的烧蚀实验 ¹⁾. 力学学报[J], 2020, 52(5): 1362-1370 DOI:10.6052/0459-1879-20-081
Luo Yue, Wang Lei, Dang Leining, Liu Jinbo, Zhang Jun, Liu Sen. ARCJET ABLATION EXPERIMENT TO SIMULATE THE CHELYABINSK ASTEROID ENTRY ¹⁾. Chinese Journal of Theoretical and Applied Mechanics[J], 2020, 52(5): 1362-1370 DOI:10.6052/0459-1879-20-081

引言

烧蚀是小行星极高速进入地球大气层后最重要的现象之一,在很大程度上决定了小行星的质量/尺寸变化、飞行轨迹、甚至光辐 射特性^[1-2]. 开展小行星进入地球大气烧蚀机理研究,有助于更好地理解小行星进入地球大气的烧蚀现象,可为解体爆炸、进入轨迹和光辐射 研究奠定基础,支撑建立小行星防御地面风险评估系统^[3].

近十年来,随着行星防御研究的开展以及人类空间探测活动的深入,空气动力学方面的研究者们开始注重小行星材料烧蚀问题的研究,这 些工作都是基于半个多世纪以来在航天飞行器热防护材料烧蚀研究中建立的理论和设备而开展^[4-11],存在极高速进入热环境、天然材质烧蚀机理、不稳定热结构等诸多问题亟待探索.目前,国外机构主要依靠激光加热、电弧(高焓等离子体)风洞两类设备在这方面开展初步实验研究^[12-19] (图1). 以NASA的Ames研究中心为代表,基于小行星威胁评估计划项目(ATAP),从2014年开始在LHMEL 设备上开展陨石激光烧蚀实验研究^[13],相比于Milley^[12],其设备能量和测试仪器都有较大提升,表面热流达到150 MW/m$^{2}$. 分析结果显示,激光烧蚀实验中石陨石材料损失以汽化蒸发为主,而玄武岩材料以熔融物溅射为主. 随后的激光烧蚀/低速剪切联合实验得到了不同辐射热流、不同剪切力对烧蚀光谱信号影响的定性结论,其实验数据提供给ATAP以研究数值计算方法. 2016年及2017年,烧蚀实验在60 MW电弧风洞上开展,试件材料包括石陨石(Tamdakht)、铁陨石(Campo Del Ceilo)、玄武岩及石英,调试了驻点热流40 MW/m$^{2}$、驻点压力140 kPa,总焓20 MJ/kg的来流条件,模拟直径30 m的陨石在65 km高空以20 km/s速度飞行时的热环境. 陨石材料均在5 s时间内完全烧尽,没有残留实验后的模型,但记录了重要的动态过程影像. 基于实验结果分析,Ames认为小行星材料的烧蚀主要包含蒸发、熔融和机械剥蚀等机理,估算有效烧蚀焓为2 MJ/kg,初步数值计算与测量的发射谱线有明显偏差,部分原因可能是缺少熔化模型^[14].

图1

新窗口打开|下载原图ZIP|生成PPT
图1国外开展的陨石材料烧蚀实验研究^[18]

Fig. 1Ablation Experiments of meteorite abroad ^[18]



国内之前尚未开展小行星材料烧蚀机理研究. 为观测小行星材料在超高速高温流场中的烧蚀现象,中国空气动力研 究与发展中心超高速空气动力研究所首次开展了模拟小行星进入的烧蚀实验研究,在电弧加热设备平台上构建Chelyabinsk小行星极高速进入地球大气层的典型气动加热环境,针对玄武岩和碳钢实验件进行烧蚀测试,测量了烧蚀动态过程、气体发射光谱、实时外形变化、表面温度变化等数据.

1 实验方法

1.1 实验方案和测量设备

该项实验在中国空气动力研究与发展中心的20 MW片式电弧加热器上采用自由射流驻点烧蚀的方法^[20-21]进行,电弧加热器产生的高温气体经喷管加速形成高焓高速气流,对实验模型前端面冲击加热,通过调整加热器状态参数模拟小行星极高速进入大气层时的驻点热流、压力、焓值,研究材料的烧蚀特性,实验和测试方法示意见图2.

图2

新窗口打开|下载原图ZIP|生成PPT
图2实验和测试方法示意图

Fig. 2Diagram of experimental set-up



高清摄像机用于记录流场结构和局部材料瞬时变化过程,比色高温计通过获取材料双波段辐射信号测量模型前端温度,红外摄像仪用于获取模型表面温度分布,光栅光谱仪通过聚焦透镜和光纤探头测量流场中的发射光谱,基于结构光的视觉系统在线实时测量模型纵 截面轮廓变化.

实验后的参数测量主要包括烧蚀外形、表面形态、烧蚀产物成分：利用三维激光扫描仪测量烧后模型三维外形及尺寸,通过对比原始外形获取总的质量损失;扫描电镜能谱仪(SEM/EDS)通过捕捉材料的二次电子和X射线连续谱,获得烧蚀模型表面材料化学成分,测试仪器配置见表1所示.

Table 1
表1
表1测试仪器配置
Table 1Measurement instrument

新窗口打开|下载CSV

1.2 实验状态和模型

材料驻点烧蚀特性地面考核实验通常模拟试件的驻点热流、焓值以复现烧蚀过程^[22],根据公开发表已观测到的Chelyabinsk小行星进入大气层轨迹 ^[23] (图3),结合现有电弧风洞能力(焓值小于30 MJ/kg,功率小于50 MW),此次实验模拟Chelyabinsk小行星在5.6 km/s的速度下,在17 km高度 解体为1 m直径后的驻点热流和焓值分别为25 MW/m$^2$和16 MJ/kg.

图3

新窗口打开|下载原图ZIP|生成PPT
图3Chelyabinsk小行星进入轨迹 ^[23]

Fig. 3Entry trajectory of Chelyabinsk asteroid ^[23]



实验模型材质为玄武岩和碳钢,分别模拟石质陨石和铁陨石,外形尺寸根据来流条件和需模拟的参数而设计(图4),头部半径为20 mm,半锥角9$^\circ$. 模型底部加工内螺纹孔,与支座连接,由快速送进系统控制试件进、出高温射流区,以确保在稳定的流场中烧蚀4 s时间,图4中黑色圆圈为比色高温计和光谱仪聚焦测量点.

图4

新窗口打开|下载原图ZIP|生成PPT
图4实验模型尺寸及照片

Fig. 4Model dimentions and image



喷管出口尺寸为$\phi $42mm,出口名义马赫数 1.8,前段弧室总压0.56 MPa,气流总焓16 MJ/kg. 焓值由总压流量法计算 ^[20],波后驻点压力和热流分别采用总压探针和零点量热计测量,流场结构如图5所示.

图5

新窗口打开|下载原图ZIP|生成PPT
图5流场结构示意

Fig. 5Flow field structure diagram

2 验结果及分析

2.1 烧蚀过程及形貌

碳钢模型烧蚀时自发光过于饱和,摄像机没有记录到较为清晰的烧蚀过程. 但通过模型刚退出流场一半时刻的视频截图(图6),能够观察到模型烧蚀的两个明显特征：整个视场大量火花飞溅、熔液在球、锥相切线堆积. 模型烧蚀后(图7)锥体氧化为褐黑色,表面无明显烧蚀痕迹和产物堆积;头部端面光亮平滑,熔融物凝固堆积在球、锥相切线位置,与图6特征对应,头部熔化的钢水在高压气流作用下被不断剥离,表层来不及氧化,因而呈光亮银白色,与锥面形成鲜明对比. 碳钢头部熔化为钢水后,一部分向后流动并驻留在球、锥相切线位置,一部分在气流冲击和剪切作用下溅射成微小液滴并被气流带走,从而表现为图6中的大量火花,后者也是碳钢烧蚀质量损失的主要机制. 利用激光三维扫描仪对烧蚀后的模型外形进行了测量,如图8所示,头部曲率半径为45 mm,远大于初始的20 mm,顶点后退13.1 mm,烧蚀质量为72.7 g.

图6

新窗口打开|下载原图ZIP|生成PPT
图6碳钢烧蚀结束时刻视频截图

Fig. 6Image of carbon steel at the end of ablation

图7

新窗口打开|下载原图ZIP|生成PPT
图7碳钢烧蚀前后外貌

Fig. 7Pre- and post-test image of carbon steel

图8

新窗口打开|下载原图ZIP|生成PPT
图8碳钢烧蚀后三维轮廓

Fig. 8Recession profile for carbon steel



图9为玄武岩烧蚀过程视频截图. 可以看出,玄武岩试件刚进入流场时(0 s)前端就被融化,相对于钢水来说,玄武岩熔融物黏性较大,并没有形成大量细小液珠溅射,而是受剪切向后流动. 随着烧蚀进行,熔融层逐渐覆盖整个球头并向锥身流动. 球头驻点持续后退,外形不断变化,在1.5 s及3 s时刻均有明显片状物剥落现象,在烧蚀结束时刻(4 s)熔融层已覆盖整个锥身,且表面明显可见由于不均匀剥落造成的凹坑. 烧蚀后的玄武岩(图10)近似双锥外形,越靠近底部熔融层凝固覆盖现象越明显,表面多处大小不一的块状崩裂痕迹,且在底部已有较深的贯穿性裂纹. 局部照片(图11)可见部分熔融液滴被蒸发喷射气体冲破,较厚的熔融层表面形成密集气孔. 由此可见,在实验条件下,玄武岩烧蚀质量损失以熔融物剪切流失为主,并伴随少量块状剥落及蒸发喷射. 通过测量,玄武岩试件顶点后退7.3 mm,烧蚀质量为37.9 g.

图9

新窗口打开|下载原图ZIP|生成PPT
图9玄武岩烧蚀过程视频截图

Fig. 9Image of ablated basalt at different time

图10

新窗口打开|下载原图ZIP|生成PPT
图10玄武岩烧蚀前后外貌

Fig. 10Pre- and post-test image of basalt

图11

新窗口打开|下载原图ZIP|生成PPT
图11玄武岩烧蚀后局部外貌

Fig. 11Local detals of basalt after ablation

2.2 表面温度

高温计采用德国IMPAC公司的ISR12-LO M33比色高温计,测试波段分别为0.8 $\mu$m和1.05 $\mu$m,量程为1000${^\circ}$C $\sim$ 3300${^\circ}$C,通过获取双波段辐射信号计算图4中球、锥相切线聚焦点(约$\phi $8 mm)的温度,因此不受测量窗口透射率的影响,在其测量波段受高温空气辐射影响也较小(高温空气在0.9 $\sim$ 3.2 $\mu$m范围辐射能量仅占5.9%).热像仪采用德国INFRA公司的ImageIR 8300 hp,测试波段为2 $\sim$ 5.7 $\mu$m,参考NASA的试验^[14]设置碳钢和玄武岩的发射率为0.4和0.8,但其实际发射率受窗口衰减和温度变化并不能确定. 图12和图13显示了碳钢和玄武岩的比色高温计所测聚焦点温度、热像仪所测试件驻点温度以及从热图中提取的聚焦点平均温度的历程变化. 碳钢的驻点位置在开始0.5 s内快速升温,并基本稳定在约1500${^\circ}$C,热像仪所测聚焦点温度经历约1.6 s时间上升后基本稳定在约1600${^\circ}$C,并在结束前0.5 s时刻快速下降,显示了随着材料焼蚀后退,$\phi $8 mm测温点的移动过程：上升阶段是从低温区移向前端高温区,稳定阶段是在高温区范围内,下降阶段是部分移出试件表面. 比色高温计和热像仪所测聚焦点温度均高于驻点温度,与温度分布热图相符,且远高于熔点1300${^\circ}$C. 玄武岩驻点温升速度明显大于碳钢,由于焼蚀后退量较小,聚焦点温度没有出现类似碳钢的后段快速下降现象,在稳定阶段驻点温度和聚焦点温度相近,约为1900${^\circ}$C,同样远高于熔点980${^\circ}$C $\sim$ 1200${^\circ}$C^[24]. 比色高温计的测量数据均只有后半段的约2.5 s历程,即使测点温度高于测量下限(1000${^\circ}$C)也没有信号输出,且在输出信号开始时突然跃升,可能是由于前1.5 s时间段高温计所获取能量没有达到仪器所预设能量阈值,没有及时输出测算信号.

图12

新窗口打开|下载原图ZIP|生成PPT
图12碳钢烧蚀过程表面温度历程

Fig. 12Temperature history of carbon steel from pyrometer

图13

新窗口打开|下载原图ZIP|生成PPT
图13玄武岩烧蚀过程表面温度历程

Fig. 13Temperature history of carbon steel from pyrometer



红外热像仪所测表面温度分布变化历程 (图14和图15),印证了前文分析的模型烧蚀过程. 碳钢最高温度始终位于球、锥相切线附近,由于驻点附近的半球面熔液被不断快速剥离,半球面的温度低于球、锥相切线位置,而锥身位置没有被直接加热,温度远低于前端. 玄武岩的表面高温区域由端头局部逐渐向后扩展至整体,也是玄武岩熔融物不断流动覆盖锥面的过程.

图14

新窗口打开|下载原图ZIP|生成PPT
图14碳钢烧蚀表面温度分布云图

Fig. 14Temperature field history of carbon steel from infrared imaging

图15

新窗口打开|下载原图ZIP|生成PPT
图15玄武岩烧蚀表面温度分布云图

Fig. 15Temperature field history of basalt from infrared imaging

2.3 外形变化及有效烧蚀焓

对光学视觉形貌测量系统所测玄武岩动态烧蚀过程进行图像数字处理,提取了不同时刻模型外形轮廓,见图16所示. 经过烧蚀后的玄武岩试件头部外形有较大变化,驻点等效半径增大,球面变为锥面,整个模型轮廓近似为双锥形. 通过前端点相对位置变化,获取了驻点后退量的变化历程,如图17所示,玄武岩试件后退量近似为线性增长,并没有明显的初始缓变过程,说明玄武岩在极短的时间内完成吸热升温(小于0.3 s),并迅速熔化后退.

图16

新窗口打开|下载原图ZIP|生成PPT
图16玄武岩烧蚀轮廓动态变化

Fig. 16Time evolution of profile for basalt

图17

新窗口打开|下载原图ZIP|生成PPT
图17玄武岩驻点后退历程及线性拟合

Fig. 17Recession history of Stagnation point on basalt



有效烧蚀焓是反映牺牲性防热材料耐烧蚀性能的重要参数,其定义为 ^[25]

(1) $ h_{\rm eff} = \dfrac{q_{\rm hw} }{m} = \dfrac{q_{\rm cw} }{m} \left[ {\left( {h_{\rm e} - h_{\rm hw} } \right) / \left(h_{\rm e} - h_{\rm cw} \right ) } \right] $
其中,$q_{\rm hw} $为热壁热流,$q_{\rm cw} $为冷壁热流,$h_{\rm e} $为边界层外缘焓值,$h_{\rm hw} $为热壁焓,$h_{\rm cw}$为冷壁焓,$m$为单位面积质量损失率. 在此次实验条件下,$h_{\rm cw} $远小于$h_{\rm e} $,且驻点位置的$h_{\rm e}$可近似为总焓 ,因此上式可简化为

(2) $ h_{\rm eff} = \dfrac{q_{\rm cw} }{m}\left( {1 - \dfrac{h_{\rm hw} }{h_0 }} \right) $
式中,$h_0 $为已知的实验测量参数16 MJ/kg,$h_{\rm hw} $由所测驻点温度确定为2.5 MJ/kg,根据前文得到的烧蚀外形变化数据,驻点呈线性后退,半径由20 mm变化至30 mm,由Fey-Riddle公式计算得到时间平均驻点热流为22.5 MW/m$^{2}$. 质量损失在半球面上并非均匀分布,而这种分布主要由热流分布导致,因此可根据球头层流无量纲热流分布情况对时间平均驻点热流作面积加权平均,得到等效冷壁热流$q_{\rm sw}$为14 MW/m$^{2}$. 这样,$m$则可直接利用质量烧蚀率; $m_{\rm t}$计算

(3) $ h_{\rm eff} = \dfrac{q_{\rm sw} }{\left( {m_{\rm t} / s} \right)}\left( {1 - \dfrac{h_{\rm hw} }{h_0 }} \right) $
根据上述方法估算的玄武岩材料有效烧蚀焓约为2.6 MJ/kg,与文献^[14]的2 MJ/kg接近,且远小于文献^[26]的蒸发烧蚀焓,可见地面模拟条件下玄武岩试件的质量损失是以熔融层剪切流失为主导.

2.4 烧蚀发射光谱

碳钢和玄武岩试件烧蚀过程中的发射光谱测量结果如图18所示,(a)为碳钢,(b)为玄武岩. 从图18中可以看出碳钢和玄武岩在烧蚀过程中的发射光谱包含了O,N,Cu,Fe,K,Na,Cr,Si,Ca等元素的原子辐射谱线,这些元素大致可以分为3类：N和O原子来自实验气体中氮气和氧气的解离;Cu原子来自电弧加热器电极的烧蚀;Fe,K,Na,Cr,Si,Ca等为实验件材料元素,其中碳钢中的Fe,Cr为主要材质元素和添加剂,Na,K可能是炼制过程中炉渣带入的少量硅酸盐杂质,玄武岩中的Si,Ca为主要材质元素(未测到Al),Na,K为碱金属杂质,也是人工陨石的典型添加元素 ^[13].

图18

新窗口打开|下载原图ZIP|生成PPT
图18碳钢及玄武岩发射光谱测量结果

Fig. 18Emission spectrum of carbon steel and basalt



图19显示了碳钢和玄武岩试件烧蚀过程中各元素谱线的时间演化,从图中看出在碳钢和玄武岩的烧蚀中Na,K的谱线都是最先出现并迅速升高,反映了碱金属在高温下易挥发不稳定性. 烧蚀过程中模型端头随着烧蚀量的增加而后退,使得测量点更加接近驻点,因此各元素谱线强度随着烧蚀的进行都有一定的上升趋势.

图19

新窗口打开|下载原图ZIP|生成PPT
图19碳钢及玄武岩烧蚀过程光谱强度时间演化

Fig. 19Evolution of spectral line intensity



利用碳钢实验中的多条谱线,通过多谱线相对强度法(玻尔兹曼分布法)获取了实验气体的总温. 根据原子光谱理论,假设模型端头为局部热力学平衡态,各原子能级分布满足玻尔兹曼分布,原子谱线强度和上能级能量满足下式

(4) $ {\rm In} \left( {\dfrac{I\lambda }{gA}} \right) = - \left( {\dfrac{E}{kT}} \right) + {\rm In} \left( {\dfrac{N}{U}} \right) $
式中,$E, g, A, \lambda $和$I$分别为谱线的上能级能量、统计权重、跃迁几率、波长和强度,$k$为玻尔兹曼常数,$N$为总态数,$U$为布居函数.对于选定原子的谱线,式中第3项为常数,因此温度可由斜率 $(-1/kT)$得到,选择多条谱线对${\rm In} (I \lambda/gA) $和$E$进行直线拟合,得到斜率进而求出温度,根据图18, 选定Fe原子的4条谱线和铜原子的3条谱线进 行作图(图20).根据斜率分别算出气体总温为 6634.4 K(铁)和6754.1 K(铜),与气流总焓对应的温度相符.

图20

新窗口打开|下载原图ZIP|生成PPT
图20铁原子($\times $)和铜原子($\square $)玻尔兹曼图

Fig. 20Boltzmann plot of Fe($\times $) and Cu($\square $) lines

2.5 电镜能谱测量

利用扫描电镜能谱仪对烧蚀后玄武岩试件进行了表面材料成分定性测量. 从测量结果(图21)看,玄武岩主要含有C,Ca,O,Na,Mg,Al,Si等元素成分. 对比光谱测量数据,除红外波段的C原子谱线和紫外波段的Al和Mg原子谱线处于仪器测量范围(450 $\sim$ 1100 nm)边缘,没有有效捕捉,其余元素的两种测量结果是对应一致的.

图21

新窗口打开|下载原图ZIP|生成PPT
图21烧蚀后玄武岩电镜能谱扫描

Fig. 21Post-test SEM/EDS result of basalt

3 结语

(1) 通过实验记录了玄武岩和碳钢试件烧蚀动态过程,观察到了清晰的烧蚀现象,分析揭示了两者不同的质量损失机制：碳钢端头熔化为钢水,在气流冲击作用下溅射成大量微小液滴,跟随气流高速流失;玄武岩端头熔化为黏性较大的熔融层,在气流剪切作用下以表面流动损失为主,伴随少量块状剥落及蒸发喷射.

(2) 利用实验获取的玄武岩试件外形轮廓变化过程、冷壁热流、表面温度、气流总焓、质量损失等数据,估算得到实验条件下玄武岩有效烧蚀焓为2.6 MJ/kg,远小于文献^[26]的蒸发烧蚀焓8 MJ/kg,揭示了类陨石材料在低空高压与 高空低压环境下的烧蚀机理差异. 玄武岩和碳钢试件的质量损失及驻点后退量分别为37.9 g,72.7 g以及7.3 mm,13.1 mm.

(3) 在线光谱测量获得多种原子发射光谱,包括材料蒸发的Fe,K,Na,Cr,Si,Ca,空气离解的N和O,以及设备电极烧损的Cu,与电镜能谱扫描结果比对一致. 对测试数据采用多谱线相对强度法计算气流总温约为6700 K,与流场调试校测参数相符.

致谢

本文在研究中得到中国空气动力研究与发展中心超高速空气动力研究所的石卫波、杨远剑、陈卫等同志的帮助,周平、杨靖康等参与实施了实验,在此向他们表示感谢.

参考文献 View Option 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子

[1] 柳森, 党雷宁, 赵君尧 等. 小行星撞击地球的超高速问题
力学学报, 2018,50(6):1311-1327
[本文引用: 1]

( Liu Sen, Dang Leining, Zhao Junyao, et al. Hypervelocity issues of Earth impact by asteroids
Chinese Journal of Theoretical and Applied Mechanics, 2018,50(6):1311-1327 (in Chinese)) )
[本文引用: 1]

[2] 柳森, 黄洁, 李毅 等. 行星防御中的超高速问题
// 香山科学会议,北京, 2018
[本文引用: 1]

( Liu Sen, Huang Jie, Li Yi, et al. Hyervelocity issues of planetary defense
// Xiangshan Scientific Conference, Beijing, 2018 (in Chinese))
[本文引用: 1]

[3] Dang LN, Liu S, Bai ZY. Hazard estimate of 2019 PDC impact scenario
// 6th IAA Planetary Defense Conference, Washington DC, USA, 2019
[本文引用: 1]

[4] Chen YK. Thermal ablation modeling for silicate materials
// 54th AIAA Aerospace Science Meeting, 2016
[本文引用: 1]

[5] Johnston CO, Samareh J. Aerothermodynamic characteristics of 16-22 km/s earth entry. AIAA 2015-3110


[6] Johnston CO, Stern EC. Impact of coupled radiation and ablation on the aerothermodynamics of meteor entries
// 47th AIAA Thermophysics Conference, Denver, USA, 2017


[7] Johnston CO, Gnoffo PA, Sutton K. The influence of ablation on radiative heating for earth entry
Journal of Spacecraft and Rocket, 2009,46(3):481-491


[8] Dias B, Magin TE. Stagnation-line simulations of meteor ablation
// Annual BRAIN-BE Meeting, Brussels, Belgium, 2015


[9] Girin OG. A hydrodynamic mechanism of meteor ablation: The melt-spraying mode
Astronomy Astrophysics, 2017,606, A63


[10] Shi WB, Dang LN, Xiao Y, et al. Aerodynamic heating/ablation/conduction calculation of iron asteroid entering the earth's atmosphere at hypervelocity
// 6th IAA Planetary Defense Conference, Washington DC, USA, 2019


[11] 白智勇, 柳森, 党雷宁 等. 高温气体热力学特性及其在小行星进入问题中的应用
// 全国第二届行星防御研讨会, 北京, 2019
[本文引用: 1]

( Bai Zhiyong, Liu Sen, Dang Leining, et al. thermodynamic characteristics of high temperature gas and its application at the problem of asteroid entry
// The 2nd Planetary Defense Symposium, Beijing, 2019 (in Chinese))
[本文引用: 1]

[12] Milley EP, Hawkes RL, Ehrman JM. Meteor luminosity simulation through laser ablation of meteorites
Monthly Notices of the Royal Astronomical Society: Letter, 2007,382(1):L67-L71
[本文引用: 2]

[13] White S, Stern EC. Laser ablation experiments on the Tamdakht H5 Chondrite
// 80th Annual Meeting of the Meteoritical Society, 2017
[本文引用: 2]

[14] Agrawal P, Jenniskens P, Stern EC, et al. Arcjet ablation of stony and iron Meteorites
// 2018 Aerodynamic Measurement Technology and Ground Testing Conference, Georgia, USA, 2018
[本文引用: 3]

[15] Fretter E. Ames arc jets in support for planetary defance
// 1st International Workshop on PHA Characterization, Atmospheric Entry and Risk Assessment, 2015,09


[16] Stern EC, Chen Y-K, White S, et al. Numerical and experimental investigation of meteoroid ablation
// 1st International Workshop on PHA Characterization, Atmospheric Entry and Risk Assessment, 2015,09


[17] Loehle S, Zander F, Hermann T, et al. Experimental simulation of meteorite ablation during earth entry using a plasma wind tunnel
Astrophysical Journal, 2017,837(2):112


[18] Magin T. Investigation of the meteor phenomenon plasmatron wind-tunnel experiments and ground observations
// 61st Course of the International School of Quantum Eletroncs, Erice, 2017
[本文引用: 2]

[19] Watanabe M, Hironori S, Shinsuke A, et al. Development of artificial meteor for observation of upper atmosphere
Acta Astronautica, 2016,121:172-178
[本文引用: 1]

[20] 张松贺, 杨远剑, 王茂刚 等. 电弧风洞热/透波联合试验技术研究及应用
空气动力学学报, 2017,35(1):141-145
[本文引用: 2]

( Zhang Songhe, Yang Yuanjian, Wang Maogang, et al. Studies and applications of thermal/wave-transmission test technique in arc-heated wind tunnel
Acta Aerodynamica Sinica, 2017,35(1):141-145 (in Chinese))
[本文引用: 2]

[21] 罗跃, 周玮, 杨鸿 等. 电弧加热器湍流平板试验流场计算分析
实验流体力学, 2017, ( 2):86-92
[本文引用: 1]

( Luo Yue, Zhou Wei, Yang Hong, et al. CFD analysis of the arc heater turbulent flow field of flat plate testing
Journal of Experiments in Fluid Mechanics, 2017, ( 2):86-92 (in Chinese))
[本文引用: 1]

[22] Prabhu D, Saunders D, Oishi T, et al. CFD analysis framework for arc-heated flowfields, I:stagnation testing in arc-jets at NASA ARC
// 41st AIAA Thermophysics Conference, San Antonio, Texas, USA, 2009
[本文引用: 1]

[23] Borovicka J, Spurny P, Brown P, et al. The trajectory, structure and origin of the chelyabinsk asteroidal impactor
Nature, 2013,503(7475):235-237
DOIURL [本文引用: 3]
Earth is continuously colliding with fragments of asteroids and comets of various sizes. The largest encounter in historical times occurred over the Tunguska river in Siberia in 1908, producing(1,2) an airburst of energy equivalent to 5-15 megatons of trinitrotoluene (1 kiloton of trinitrotoluene represents an energy of 4.185 x 10(12) joules). Until recently, the next most energetic airburst events occurred over Indonesia(3) in 2009 and near the Marshall Islands(4) in 1994, both with energies of several tens of kilotons. Here we report an analysis of selected video records of the Chelyabinsk superbolide(5) of 15 February 2013, with energy equivalent to 500 kilotons of trinitrotoluene, and details of its atmospheric passage. We found that its orbit was similar to the orbit of the two-kilometre-diameter asteroid 86039 (1999 NC43), to a degree of statistical significance sufficient to suggest that the two were once part of the same object. The bulk strength-the ability to resist breakage-of the Chelyabinsk asteroid, of about one megapascal, was similar to that of smaller meteoroids(6) and corresponds to a heavily fractured single stone. The asteroid broke into small pieces between the altitudes of 45 and 30 kilometres, preventing more-serious damage on the ground. The total mass of surviving fragments larger than 100 grams was lower than expected(7).

[24] Larsen ES. The temperature of magmas
American Mineralogist, 1929,14:81-94
[本文引用: 1]

[25] ASTM. E 458-72 Standard Test Method for Heat of Ablation.ASTM Standard.
[本文引用: 1]

[26] ?pik EJ. Physics of Meteor Flight in the Atmosphere. Dover Press, 1958
[本文引用: 2]