电弧加热器试验条件下铁质小行星材料烧蚀机理分析.pdf

电弧加热器试验条件下铁质小行星材料烧蚀机理分析石卫波，党雷宁，罗 跃，孙海浩，黄 洁（中国空气动力研究与发展中心 超高速空气动力研究所，绵阳 621000）摘 要：研究铁质小行星进入地球大气层的烧蚀机理，对评估小行星撞击地球危害具有重要意义。基于试验现象建立了铁质小行星材料烧蚀的熔融烧蚀和熔融层剪切烧蚀模型，采用具有移动边界的气动热、烧蚀与内部热传导耦合求解技术对铁陨石模型烧蚀试验状态进行了计算。计算分析表明建立的剪切烧蚀模型能够得到与试验在定性上一致的驻点烧蚀速率变化规律。若不考虑液态层流失，计算得到的规律与试验相反。计算分析结果显示表面蒸发速率与熔融层质量流失速率相比是小量，表明铁质小行星烧蚀以熔融层剪切流失为主导。可为小行星撞击地球防御提供参考。关键词：铁质小行星；烧蚀；熔融；剪切流失中图分类号：P185.7；V11 文献标识码：A 文章编号：2096-9287（2023）04-0436-07DOI:10.15982/j.issn.2096-9287.2023.20230029引用格式：石卫波，党雷宁，罗跃，等.电弧加热器试验条件下铁质小行星材料烧蚀机理分析J.深空探测学报（中英文），2023，10（4）：436-442.Reference format:SHI W B，DANG L N，LUO Y，et al.Analysis of ablation mechanism of iron asteroidmaterials under the condition of arc heater testJ.Journal of Deep Space Exploration，2023，10（4）：436-442.引言近地小行星撞击地球是人类生存和发展的潜在威胁之一。近代以来著名的近地天体撞击事件有1908年通古斯大爆炸1、2003年车里雅宾斯克流星事件等2。学术界已经基本认同：6 500万年前1颗直径10 km的小行星撞击现墨西哥尤卡坦半岛位置，导致包括恐龙在内的全球物种大灭绝，史称“K-T事件”3。小行星以极高速（1270 km/s）进入地球大气层，产生巨大的气动加热，引起材料大规模熔融、汽化和烧蚀，导致严重的质量损失、外形变化甚至多次解体，产生的冲击波和热辐射向地面传播，破坏地面设施和生物4-5。小行星极高速进入的烧蚀问题，是小行星进入地球大气的科学问题之一，是研究进入轨迹、光辐射、解体乃至空中爆炸现象的重要基础4。小行星材料大致分为3类6：铁陨石、石铁陨石、石陨石。目前，学术界主要通过小行星进入事件观测、地面试验和理论分析3种手段研究小行星的烧蚀问题。在小行星进入事件观测中，可通过观测到的轨迹分析烧蚀系数7。然而由于铁质小行星占比少 8或缺乏观测数据，铁质小行星的烧蚀系数数据较少。在地面试验方面，美国国家航空航天局（National Aeronauticsand Space Administration，NASA）艾姆斯研究中心（Ames Research Center）2018年在60 MW电弧风洞中开展了铁质小行星材料进入地球大气层的烧蚀试验9，试验模拟热流32.9 MW/m2、滞止压力126.5 kPa，相当于30 m直径的小行星在65 km海拔高度以20 km/s的速度飞行，试件材料为IAB-MG型Campo Del Cielo铁陨石。中国空气动力研究与发展中心超高速研究所（简称CARDC超高速研究所）2021年在20 MW电弧加热器上开展了球锥外形铁陨石模型的烧蚀试验10。理论分析方面的工作大多以石质小行星为研究对象，针对铁质小行星烧蚀的工作较少。Girin11基于熔融层的失稳提出熔融喷溅的烧蚀模型，计算了6 mm直径铁质微流星的烧蚀，但缺乏与观测或地面试验的对比。apek等12提出另一种熔融喷溅的烧蚀理论，apek等 13对该理论进行了改进。该理论认为小行星表面熔融层在温度达到熔点后立即脱落，分散成液滴在大气中继续减速。该理论对质量0.0061.1 g微流星烧蚀的预测结果与光曲线观测结果符合较好。从研究现状看，国内外对铁质小行星烧蚀问题研究较少，且主要针对微流星，不仅在行星防御范畴之外，建立的烧蚀理论也缺乏地面试验详细测量支撑。本文基于CARDC超高速研究所电弧加热器铁陨石烧蚀试验中观察到的现象，建立了铁质小行星材料烧蚀的 收稿日期：2023-03-20 修回日期：2023-04-20基金项目：国家自然科学基金（11972358）第 10 卷 第 4 期深 空 探 测 学 报（中英文）Vol.10 No.42023 年 8 月Journal of Deep Space ExplorationAugust 2023熔融层剪切烧蚀模型，并采用具有移动边界的气动热、烧蚀与内部热传导耦合求解技术，对烧蚀试验进行了分析。1 铁质小行星材料烧蚀模型与计算方法CARDC超高速研究所于2021年在20 MW电弧加热器上，针对如图1所示的球锥外形，开展了3个状态（见表1）的烧蚀试验10。电弧加热器产生的高温气体经喷管加速形成高焓高速气流，对试验模型前端面冲击加热，通过调整加热器状态参数模拟小行星极高速进入大气层时的驻点热流、压力、焓值，研究材料的烧蚀特性，试验和测试方法示意如图2所示14。支座安装在试验舱中的快速送进系统上，前端安装试验模型。试验过程中的参数测量：模型背面设置热电偶测量背温；高清摄像机记录流场结构和局部材料瞬时变化过程；比色高温计测量模型驻点附近温度；红外摄像仪获取模型表面温度分布；光栅光谱仪测量流场中的发射光谱；双目视觉系统在线实时测量模型纵截面轮廓变化。试验后的参数测量主要包括烧蚀质量、线烧蚀率。试验视频截图（气流自左向右）如图3所示，可看到熔融层受剪切向后流动并覆盖试验模型全表面。图 3 铁石烧蚀过程的视频截图Fig.3 Video screenshot of iron stone ablation process 上述试验现象得到的启发：把熔融、蒸发和剪切作为铁质小行星材料烧蚀的主要因素，熔融是材料在气动加热下由固态到液态的相变，蒸发是由液态到气态的相变，剪切是高速气流对熔融层的机械剥蚀作用。为评估熔融、蒸发和剪切在铁质小行星材料烧蚀中的作用，本文分别建立熔融烧蚀模型和熔融剪切烧蚀模型，模拟电弧射流中的球锥外形铁陨石烧蚀情况。熔融烧蚀模型假设表面达到熔点后熔融层立即去除。熔融剪切烧蚀模型则考虑了表面熔融层的蒸发以及在气流剪切下的流失。烧蚀计算需要和气动热、材料内部热传导耦合起来，是一个非定常问题（图4）。气动热是热传导计算 50R2015812550（参考）单位:mm 图 1 铁陨石烧蚀试验模型Fig.1 Iron meteorite ablation test model 表 1 烧蚀试验模拟状态Table 1 Simulation state of ablation test状态弧室压力/MPa驻点压力/MPa驻点热流/（MWm2）滞止焓（MJkg1）1.020.5113.97.90.660.3415.110.20.560.2819.514.3 高速摄像机比色高温计红外热像仪高清像机光谱仪观察窗试验舱反射镜支座热电偶高温气流排气系统模型 图 2 试验和测试方法示意图Fig.2 Diagram of experimental set-up 读入材料热物性数据（密度、比热、导热射系数、化学组成）读入喷管出口参数（速度、压温度）读入几何初始外形热环境计算（压力、热流、焓）内部热响应计算、材料烧蚀计算（温度、烧蚀速率）烧蚀外形计算计算数据输出（烧蚀外形、表面热流分布、温度场等参数随试验时间变化数据）t试验结束时间结束是否系数、黏度系数、辐 图 4 铁陨石模型气动热/烧蚀耦合计算流程图Fig.4 Flowchart of aerodynamic heat/ablation coupling calculation of ironmeteorite model第 4 期石卫波，等：电弧加热器试验条件下铁质小行星材料烧蚀机理分析437的边界条件，内部热传导结果是烧蚀计算的输入条件，烧蚀引起的外形变化对气动热产生影响，烧蚀引起的物面后退影响热传导物面边界条件。其中内部热响应和材料烧蚀计算采用紧耦合，其它模块采用松耦合。1.1 气动热计算方法电弧加热器试验条件下模型气动热计算根据模型表面不同区域和边界层不同流动状态选取，驻点热流采用Fay-Riddell公式15；层流热流采用Zoby等16给出的参考焓方法计算简捷，具有较高的精度；湍流热流采用布拉休斯平板表面摩阻关系式、Eckert的参考焓压缩性修正以及Colburn的雷诺比拟关系相结合的方法17。1.2 熔融烧蚀模型与计算方法对于极高速进入地球大气层的铁质小行星，气动加热作用时间短，表面热流只影响局部厚度的温度响应，尚未深入到物体内部，因此，可将铁质小行星的瞬态传热看作半无限大物体的瞬态传热问题。一维烧蚀过程的导热微分方程为182Tr2=1aTt(1)Traa=/(c)t其中：为 处的温升；为介质的热扩散系数，；为导热时间。初始条件为t=0，T=T0(2)边界条件为r=S(t)，T=TmTr+qw=LdS(t)dt(3)S(t)Lqw=rS(t)其中：为移动界面；为熔解热（Jkg1）；传入热量等于金属熔化吸收和向小行星内部导热。烧蚀时表面在移动，采用相对动坐标，导热方程经坐标变换为d2T2+1adS(t)tdTd=0(4)Tqw=constV=const其中：是 的函数，偏微分方程转化为常微分方程。由于假定在小的时间步长内小行星表面热流密度为常数，由此可认为小的时间步长内烧蚀表面的后退速度也是常数，即材料在某网格点是逐层烧蚀。边界条件也作相应变换后得到d2Td2+VadTd=0=0,T=TmdTd+qw=VL ,T=T0(5)由式（5）可解得小行星内温度分布为T T0TmT0=exp(Va)=expV(rS(t)/a(6)dTd?=0对式（6）求导，并以代入，由第2个边界条件可得qw=VL+Va(TmT0)=S(t)tL+cp(TmT0)(7)a=/(cp)cp其中：为热扩散率；为定压比热。t将式（7）简化后，得到在 时刻小行星表面烧蚀厚度为S(t)=qwtL+cp(TmT0)(8)由式（8）可得金属材料表面熔融烧蚀速率V=qwL+cp(TmT0)(9)1.3 熔融剪切烧蚀模型与计算方法铁质小行星材料在严酷热环境下，一方面与氧发生燃烧，另一方面在气动力/热作用下，达到熔点熔化向后流动。在熔融剪切模型中，假设烧蚀表面存在一层薄的液体层，所有化学反应都发生在液体层表面，且达到平衡。下面建立液体层表面能量守恒、质量守恒和动量守恒关系式。通过求解这些关系式，即可获得烧蚀特性。在固液交界面会发生如下熔化反应Fe(s)Fe(l)(10)在液体层外表面，存在液态层的蒸发和燃烧反应Fe(l)Fe(g)(11)3Fe(s)+2O2(g)Fe3O4(s)(12)考虑对流换热、辐射传热、反应热与烧蚀带走的热量，则物面的能量平衡方程为Krqc(1hawhr+Qhr)qr=TVCT(TwT0)+H,有液态层ksTy?solid+TVH,无液态层(13)438深空探测学报（中英文）2023年hawhrqrQqcCTTVHksTw其中：为空气在壁温时的焓；为恢复焓；为辐射热；为反应热；为冷壁热流；为熔化前材料的比热；为材料密度；为烧蚀速率；为单位质量熔化潜热；为固体热传导系数；为物面温度。haw空气在壁温时的焓表达式为haw=CO2,ehO2+CN2,ehN2(14)hC其中：为静焓；为质量分数。Q反应热用下式计算Q=BFeQFe=FefFeBQFe(15)BFe=FefFeB(16)QFe=hFe+MOMFehO2MFe3O4MSiO2hFe3O4(17)QFeH其中：在计算中取正值。熔化潜热的表达式为H=fFe2hm(18)hm其中，=269.55 kJ/kg。式（15）中其它各量表达式为B=mKrqc/hr=TVKrqc/hr(19)=11+AfB(20)对式（13）左侧化简得到q=Krqc(1hawhr)+FefFeBQFehr(21)Le=1假定边界层内化学冻结，化学反应在壁面进行，且假定普朗特数和Lewis数情况下，由物面的质量守恒关系的边界层扩散方程可得到 mFe,w=mwKFe,weueCH(KFe,eKFe,w)(22)KFe,e=0eueCH设式（22）中，并除以有BFefFe=(1+fB)KFe,w(23)BB=meueCHff=FefFefFeFe其中：为无量纲烧蚀速率，；为总汽化分数，；为元素的汽化分数。Pi=KiMPeMi从道尔顿定律得1KFe,w=MPePVMFe(24)由式（23）和（24）可得Fe=1fFe1BMMFePePV1(25)MPePVFe其中：为气体平均分子量；为边界层外缘压力；为的蒸气压。PV=101 32510(6.34719 547Tw)(26)在液体层表面有如下动量守恒关系TL(1 f)V=22wL(dwds2d2peds2)，驻点(27)TLws0(1 f)VwLrbds=rb2(w2dpeds)，非驻点(28)TLwrb其中：为固体材料密度；为液体层密度；为液体层厚度；为剪切应力；为物体径向坐标。BTWFeMBB式（13）、（25）、（27）（28）是联系4个未知量、的超越方程组，可用迭代法求解。其中可以从式（13）和的定义得到B=hsCT(TwT0)+H(1hawhr)+FefFeQFeBhrT4wqrKrqw0(29)V联立方程组求解得，然后给出烧蚀速率V=BKr qw0Ths(30)2 铁陨石烧蚀计算结果分析铁陨石烧蚀试验是在CARDC超高速研究所20 MW电弧加热器设备开展。为便于铁陨石烧蚀理论研究与试验测量结果对比分析，根据模拟的参数设计成球锥模型，头部半径R=20 mm，半锥角9，高度70 mm（图1）。铁陨石热物性参数为19：密度7 827.234 kg/m3，熔点1 846 K，比热364 19 J/(kgK)，42.9 15.5 W/(mK)，辐射系数0.671，熔化潜热269.55 kJ/kg。铁陨石模型烧蚀试验状态参数有3个（表1），试验时间4 s。根据1.2节建立的铁陨石熔融烧蚀计算模型，本文分别对表1的3个试验状态进行了计算。第状态驻点平均线烧蚀计算速率1.55 mm/s，第状态驻点平均线烧蚀计算速率2.14 mm/s，第状态驻点平均线烧蚀计算速率2.50 mm/s。烧蚀开始的时间随热流的增加提第 4 期石卫波，等：电弧加热器试验条件下铁质小行星材料烧蚀机理分析439前，得到的驻点烧蚀速率也随驻点热流的增加而增大，而地面试验测得的驻点烧蚀速率随驻点压力的增加而增大（图5，表2）。计算与试验在烧蚀速率差异的原因是铁陨石烧蚀计算模型为熔融模型，假定表面材料温度达到熔点就发生质量损失，没有考虑铁陨石试验高温熔融流动剪切作用。0.003 00.002 5Vw/（ms1）0.002 00.001 50.001 00.000 500 0.5 1.0 1.5 2.0 2.5 3.0 3.5t/s4.0 4.5 5.0 5.5 6.0第状态第状态第状态 图 5 驻点烧蚀速率随时间变化（熔融烧蚀模型）Fig.5 Blation rate of stagnant point changes with time（meltingablation model）表 2 铁陨石球锥试件烧蚀试验烧蚀特征计算与试验对比Table 2 Calculation and experimental comparison of ablationcharacteristics of iron meteorite cone specimen状态驻点烧蚀速率/（mms1）液态层厚度/mm表面蒸发速率/（mms1）液态层质量流失速率/（mms1）计算（熔融剪切模型）计算（熔融模型）试验计算（物面6 mm处）计算计算1.911.559#：1.972.92.0710-191.911.862.1411#：1.552.32.7610-191.861.802.5010#：1.562.23.3410-191.80 根据1.3节建立的铁陨石熔融剪切烧蚀计算模型，本文对铁陨石球锥试件（图1）试验状态开展计算分析，计算与试验测量结果比较见表2。铁陨石球锥试件烧蚀试验的计算结果如图69所示。图67横坐标为时间，图8、图9横坐标为沿物面弧长（从驻点算起），图6纵坐标为驻点温度，图7、图9纵坐标为烧蚀速率，图8纵坐标为液态层厚度。试验状态的试验条件是热流递增、驻点压力递减。热流越高，导致温度越高（图6）、越早烧蚀，状态开始烧蚀的时间分别为3、1.5和1 s（图7）。稳定烧蚀后，驻点烧蚀速率随试验驻点压力增加而增大（图7）。液体层的黏性系数随壁温的增加而降低，因此热流和温度越高（从状态），液体层越薄（图8）。从液体层厚度沿物面的分布看，在头部区域附近较厚，在身部较薄，与试验观察的现象一致（图3）。头部高热流区域，铁陨石材料在高温下熔融并不断被气流吹走，部分覆盖在身部低热流区，而模型身部大部分区域热流低，陨石材料还未发生熔融烧蚀（图9）。2 4002 2002 0001 8001 6001 4001 2001 00080060000.51.01.52.02.53.03.54.0t/sTw/K第状态第状态第状态 图 6 驻点温度随时间变化Fig.6 Temperature of stagnant point changes with time Vw/（ms1）0.002 00.001 50.001 00.000 5002468t/s第状态第状态第状态 图 7 驻点烧蚀速率随时间变化Fig.7 Ablation rate of stagnant point changes with time surface 0.003 00.003 50.002 50.002 00.001 50.001 00.000 5000.010.020.03SF/mDTAL/m0.040.05第状态第状态第状态 图 8 液体层厚度随物面变化Fig.8 Thickness of liquid layer changes with the object surface 440深空探测学报（中英文）2023年00.010.020.03SF/m0.040.050.002 00.001 50.001 00.000 50Vw/（ms1）第状态第状态第状态 图 9 烧蚀速率随物面变化Fig.9 Ablation rate changes with object surface 计算与试验对比如表2所示，表2中包含了熔融剪切烧蚀模型与不考虑剪切的熔融烧蚀模型计算结果。在熔融剪切烧蚀模型下，状态的驻点线烧蚀速率分别为1.91、1.86和 1.80 mm/s，随试验状态的变化规律在定性上与试验相符。而熔融模型得到的驻点烧蚀速率与试验相反，且热流越高，计算结果与试验测量结果相差越大。这说明熔融层的剪切是烧蚀模型中必须考虑的因素。由表2还可看到，表面蒸发速率与液态层流失速率相比是小量，说明熔融层的剪切流失是烧蚀的主要机理。3 结论本文针对铁质小行星以极高速进入地球大气层，在严酷气动加热作用下发生高温熔融烧蚀问题，分别建立了熔融烧蚀计算模型和熔融剪切烧蚀计算模型。通过对电弧加热器试验条件下铁陨石球锥模型试验状态的计算分析，获得了铁陨石高温熔融层剪切流失的烧蚀机理。下一步，将根据建立的铁质小行星进入大气层高温熔融剪切烧蚀模型，评估以极高速进入地球大气层铁质小行星熔融烧蚀情况，为小行星撞击地球防御提供参考。参考文献参考文献 CHYBA C F，THOMAS P J，ZAHNLE K J.The 1908 tunguskaexplosion：atmospheric disruption of a stony asteroidJ.Nature，1993361：40-44.1 POPOVA O P，JENNISKENS P，EMELYANENKO V，et al.Chelyabinsk airburst，damage assessment，meteorite recovery andcharacterizationJ.Science，2013，342（6162）：1069-1073.2 SCHULTE P，ALEGRET L，ARENILLAS I，et al.The chicxulubasteroid impact and mass extinction at the cretaceous-paleogeneboundaryJ.Science，2010，327：1214-1218.3 柳森，党雷宁，赵君尧，等.小行星撞击地球的超高速问题J.力学学报，2018，50（6）：1311-1327.4LIU S，DANG L N，ZHAO J Y，et al.Hypervelocity issues of earthimpact by asteroidsJ.Chinese Journal of Theoretical and AppliedMechanics，2018，50（6）：1311-1327.党雷宁，柳森，白智勇，等.小行星进入与撞击效应评估模型敏感性研究J.力学学报，2021，53（1）：278-292.DANG L N，LIU S，BAI Z Y，et al.Sensitivity research on models ofearth entry and impact effects by asteroidsJ.Chinese Journal ofTheoretical and Applied Mechanics，2021，53（1）：278-292.5 MASON B H.MeteoritesJ.American Scientist，1967，55（4）：429-455.6 EPLECHA Z，BOROVIKA J，ELFORD GW，et al.Meteorphenomena and bodiesJ.Space Science Reviews，1998，84：327-471.7 李世杰.小行星的分类及飞行器探测历史J.大众科学，2021，（1）：56-59.8 AGRAWAL P，JENNISKENS P.Arcjet ablation of stony and ironmeteoritesC/Proceedings of 2018 Aerodynamic MeasurementTechnology and Ground Testing Conference.Georgia，USA：s.n.，2018.9 王磊，罗跃，党雷宁，等.模拟陨石进入的电弧风洞烧蚀试验及辐射特性分析C/第四届全国行星防御研讨会论文集.桂林：第四届全国行星防御研讨会，2021.10 GIRIN O G.A hydrodynamic mechanism of meteor ablation：the melt-spraying modeJ.Astronomy Astrophysics，2017，606：A63.11 APEK D，BOROVIKA J.Ablation of small iron meteoroidsfristresultsJ.Planetary and Space Scicnce，2017，143：159-163.12 APEK D，KOTEN P，BOROVIKA J，et al.Small iron meteoroids-Observation and modeling of meteor light curvesJ.Astronomy&Astrophysics，2019，625：A106.13 罗跃，王磊，党雷宁，等.模拟Chelyabinsk小行星进入的烧蚀实验J.力学学报，2020，52（5）：1362-1370.LUO Y，WANG L，DANG L N，et al.Arcjet ablation experiment tosimulate the chelyabinsk asteroid entryJ.Chinese Journal ofTheoretical and Applied Mechanics，2020，52（5）：1362-1370.14 FAY J A，RIDDELL F R.Theory of stagnation point heat transfer indissociated airJ.Journal of Aerosol Science，1958，25（2）：73-85.15 ZOBY E V，MOSS J N.SUTTON K.Approximate convective heatingequations for hypersonic flowsJ.Journal of Service Research，1981，18：64-70.16 张志成，潘梅林，刘初平，等.高超声速气动热和热防护M.北京：国防工业出版社，2003：222-284.17 吴斌，夏伟，汤勇，等.身管熔化烧蚀的预测数学模型J.火炮发射与控制学报，2002，24（1）：5-10.WU B，XIA W，TANG Y，et al.Mathematical model for predictingmelt-erosion in gun barrelsJ.Journal of Gun Launch&Control，2002，24（1）：5-10.18 OSTROWSKI D，BRYSON K.The physical properties of meteorites J.Planetary and Space Science，2019，165：148-178.19作者简介：石卫波石卫波（1972），男，高级工程师，主要研究方向：气动热与热防护。通信地址：四川省绵阳市二环路南段6号中国空气动力研究与发展中心（621000）电话：[phone]黄洁黄洁（1968），女，研究员，主要研究方向：高超声速空气动力学，超高速碰撞动力学。本文通信作者本文通信作者。通信地址：四川省绵阳市二环路南段6号中国空气动力研究与发展中心（621000）第 4 期石卫波，等：电弧加热器试验条件下铁质小行星材料烧蚀机理分析441Analysis of Ablation Mechanism of Iron Asteroid Materials under theCondition of Arc Heater TestSHI Weibo，DANG Leining，LUO Yue，SUN Haihao，HUANG Jie（Hypervelocity Aerodynamics Institute，China Aerodynamics Research and Development Center，Mianyang 621000，China）Abstract：When the asteroid enters earths atmosphere at a very high speed，the surface temperature rises and the surfacematerial melts and loses under the severe aerodynamic heat.Studying the ablation mechanism of iron asteroids entering the Earthsatmosphere is of great significance to study the ablation mechanism of iron asteroids entering the Earths atmosphere to evaluate theimpact of such asteroids on the earth.In 2021，the Hypervelocity Aerodynamics Institute of China Aerodynamics Research andDevelopment Center carried out the ablation test of the spherical cone-shaped iron meteorite model（head radius 20 mm，half coneangle 9 degrees）on the arc heater.The tested simulated state is that the stagnation heat flux is 13.919.5 MW/m2 and the stagnationpressure is 0.510.28 MPa.In this paper，based on the experimental phenomenon，the melting ablation model and the melting layershear ablation model of iron asteroid material ablation are established.and the ablation test state of meteorite model is calculated byusing the coupled solution methodology of aerodynamic heat，ablation and internal heat conduction with moving boundary.Thecalculation and analysis show that the shear ablation model established in this paper can obtain qualitatively consistent stagnationpoint ablation rate.If the loss of liquid layer is not considered，the calculated law is contrary to the experiment.The calculationshows that the surface evaporation rate is small compared with the mass loss rate of molten layer，which indicates that the ablationof iron asteroids is dominated by the shear loss of molten layer.It can provide reference for asteroid impact Earth defense.Keywords：iron asteroid；ablation；melting；shear lossHighlights：Based on the phenomenon of iron meteorite ablation test,the melting ablation model and the melting layer shear ablation model ofiron asteroid material ablation are established.The coupled solution methodology of aerodynamic heat,ablation and internal heat conduction with moving boundary is used tocalculate the ablation test state of iron meteorite model.The established melting shear ablation model can obtain the change law of stagnation point ablation rate which is qualitativelyconsistent with the experiment.The calculation and analysis show that the ablation of iron asteroids is dominated by the shear loss of molten layer.责任编辑：宋宏，英文审校：宋利辉442深空探测学报（中英文）2023年