电弧加热器试验条件下铁质小行星材料烧蚀机理分析.pdf
《电弧加热器试验条件下铁质小行星材料烧蚀机理分析.pdf》由会员分享,可在线阅读,更多相关《电弧加热器试验条件下铁质小行星材料烧蚀机理分析.pdf(7页珍藏版)》请在咨信网上搜索。
1、电弧加热器试验条件下铁质小行星材料烧蚀机理分析石卫波,党雷宁,罗 跃,孙海浩,黄 洁(中国空气动力研究与发展中心 超高速空气动力研究所,绵阳 621000)摘 要:研究铁质小行星进入地球大气层的烧蚀机理,对评估小行星撞击地球危害具有重要意义。基于试验现象建立了铁质小行星材料烧蚀的熔融烧蚀和熔融层剪切烧蚀模型,采用具有移动边界的气动热、烧蚀与内部热传导耦合求解技术对铁陨石模型烧蚀试验状态进行了计算。计算分析表明建立的剪切烧蚀模型能够得到与试验在定性上一致的驻点烧蚀速率变化规律。若不考虑液态层流失,计算得到的规律与试验相反。计算分析结果显示表面蒸发速率与熔融层质量流失速率相比是小量,表明铁质小行星
2、烧蚀以熔融层剪切流失为主导。可为小行星撞击地球防御提供参考。关键词:铁质小行星;烧蚀;熔融;剪切流失中图分类号: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 o
3、f 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、,引起材料大规模熔融、汽化和烧蚀,导致严重的质量损失、外形变化甚至多次解体,产生的冲击波和热辐射向地面传播,破坏地面设施和生物4-5。小行星极高速进入的烧蚀问题,是小行星进入地球大气的科学问题之一,是研究进入轨迹、光辐射、解体乃至空中爆炸现象的重要基础4。小行星材料大致分为3类6:铁陨石、石铁陨石、石陨石。目前,学术界主要通过小行星进入事件观测、地面试验和理论分析3种手段研究小行星的烧蚀问题。在小行星进入事件观测中,可通过观测到的轨迹分析烧蚀系数7。然而由于铁质小行星占比少 8或缺乏观测数据,铁质小行星的烧蚀系数数据较少。在地面试验方面,美国国家航空航天局(National Aeronauti
5、csand 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。理论分析方面的工作大多以石质小行星为研究对象,针对铁质小行
6、星烧蚀的工作较少。Girin11基于熔融层的失稳提出熔融喷溅的烧蚀模型,计算了6 mm直径铁质微流星的烧蚀,但缺乏与观测或地面试验的对比。apek等12提出另一种熔融喷溅的烧蚀理论,apek等 13对该理论进行了改进。该理论认为小行星表面熔融层在温度达到熔点后立即脱落,分散成液滴在大气中继续减速。该理论对质量0.0061.1 g微流星烧蚀的预测结果与光曲线观测结果符合较好。从研究现状看,国内外对铁质小行星烧蚀问题研究较少,且主要针对微流星,不仅在行星防御范畴之外,建立的烧蚀理论也缺乏地面试验详细测量支撑。本文基于CARDC超高速研究所电弧加热器铁陨石烧蚀试验中观察到的现象,建立了铁质小行星材料
7、烧蚀的 收稿日期: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。电弧加热器产生的高温气体经喷管加速形
8、成高焓高速气流,对试验模型前端面冲击加热,通过调整加热器状态参数模拟小行星极高速进入大气层时的驻点热流、压力、焓值,研究材料的烧蚀特性,试验和测试方法示意如图2所示14。支座安装在试验舱中的快速送进系统上,前端安装试验模型。试验过程中的参数测量:模型背面设置热电偶测量背温;高清摄像机记录流场结构和局部材料瞬时变化过程;比色高温计测量模型驻点附近温度;红外摄像仪获取模型表面温度分布;光栅光谱仪测量流场中的发射光谱;双目视觉系统在线实时测量模型纵截面轮廓变化。试验后的参数测量主要包括烧蚀质量、线烧蚀率。试验视频截图(气流自左向右)如图3所示,可看到熔融层受剪切向后流动并覆盖试验模型全表面。图 3
9、铁石烧蚀过程的视频截图Fig.3 Video screenshot of iron stone ablation process 上述试验现象得到的启发:把熔融、蒸发和剪切作为铁质小行星材料烧蚀的主要因素,熔融是材料在气动加热下由固态到液态的相变,蒸发是由液态到气态的相变,剪切是高速气流对熔融层的机械剥蚀作用。为评估熔融、蒸发和剪切在铁质小行星材料烧蚀中的作用,本文分别建立熔融烧蚀模型和熔融剪切烧蚀模型,模拟电弧射流中的球锥外形铁陨石烧蚀情况。熔融烧蚀模型假设表面达到熔点后熔融层立即去除。熔融剪切烧蚀模型则考虑了表面熔融层的蒸发以及在气流剪切下的流失。烧蚀计算需要和气动热、材料内部热传导耦合起
10、来,是一个非定常问题(图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 试验和测试方法
11、示意图Fig.2 Diagram of experimental set-up 读入材料热物性数据(密度、比热、导热射系数、化学组成)读入喷管出口参数(速度、压温度)读入几何初始外形热环境计算(压力、热流、焓)内部热响应计算、材料烧蚀计算(温度、烧蚀速率)烧蚀外形计算计算数据输出(烧蚀外形、表面热流分布、温度场等参数随试验时间变化数据)t试验结束时间结束是否系数、黏度系数、辐 图 4 铁陨石模型气动热/烧蚀耦合计算流程图Fig.4 Flowchart of aerodynamic heat/ablation coupling calculation of ironmeteorite model
12、第 4 期石卫波,等:电弧加热器试验条件下铁质小行星材料烧蚀机理分析437的边界条件,内部热传导结果是烧蚀计算的输入条件,烧蚀引起的外形变化对气动热产生影响,烧蚀引起的物面后退影响热传导物面边界条件。其中内部热响应和材料烧蚀计算采用紧耦合,其它模块采用松耦合。1.1 气动热计算方法电弧加热器试验条件下模型气动热计算根据模型表面不同区域和边界层不同流动状态选取,驻点热流采用Fay-Riddell公式15;层流热流采用Zoby等16给出的参考焓方法计算简捷,具有较高的精度;湍流热流采用布拉休斯平板表面摩阻关系式、Eckert的参考焓压缩性修正以及Colburn的雷诺比拟关系相结合的方法17。1.2
13、 熔融烧蚀模型与计算方法对于极高速进入地球大气层的铁质小行星,气动加热作用时间短,表面热流只影响局部厚度的温度响应,尚未深入到物体内部,因此,可将铁质小行星的瞬态传热看作半无限大物体的瞬态传热问题。一维烧蚀过程的导热微分方程为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+
14、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其中:为热扩散率;
15、为定压比热。t将式(7)简化后,得到在 时刻小行星表面烧蚀厚度为S(t)=qwtL+cp(TmT0)(8)由式(8)可得金属材料表面熔融烧蚀速率V=qwL+cp(TmT0)(9)1.3 熔融剪切烧蚀模型与计算方法铁质小行星材料在严酷热环境下,一方面与氧发生燃烧,另一方面在气动力/热作用下,达到熔点熔化向后流动。在熔融剪切模型中,假设烧蚀表面存在一层薄的液体层,所有化学反应都发生在液体层表面,且达到平衡。下面建立液体层表面能量守恒、质量守恒和动量守恒关系式。通过求解这些关系式,即可获得烧蚀特性。在固液交界面会发生如下熔化反应Fe(s)Fe(l)(10)在液体层外表面,存在液态层的蒸发和燃烧反应F
16、e(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其中:为静焓;
17、为质量分数。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数情况下,由物面的质量守恒关系的边界层扩散方程可得到 m
18、Fe,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(dwds
19、2d2peds2),驻点(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电弧加热器设备开展。为便于铁陨
20、石烧蚀理论研究与试验测量结果对比分析,根据模拟的参数设计成球锥模型,头部半径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
21、 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
22、 驻点烧蚀速率随时间变化(熔融烧蚀模型)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)计算(熔融剪切模型)计算(熔融模型)试
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电弧 加热器 试验 条件下 铁质 小行星 材料 机理 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。