相对运动等离子体平板与电磁波相互作用研究_满良.pdf
《相对运动等离子体平板与电磁波相互作用研究_满良.pdf》由会员分享,可在线阅读,更多相关《相对运动等离子体平板与电磁波相互作用研究_满良.pdf(9页珍藏版)》请在咨信网上搜索。
1、满良,邓浩川,薄勇,等.相对运动等离子体平板与电磁波相互作用研究J.电波科学学报,2023,38(1):164-172.DOI:10.12265/j.cjors.2022003MAN L,DENG H C,BO Y,et al.Interaction between relative moving plasma plate and electromagnetic waveJ.Chinese journal of radio science,2023,38(1):164-172.(in Chinese).DOI:10.12265/j.cjors.2022003相对运动等离子体平板与电磁波相互作用
2、研究满良1,2邓浩川1,2薄勇3肖志河1,2杨利霞3*(1.中国航天科工集团有限公司第二研究院,北京 100854;2.电磁散射重点实验室,北京 100854;3.安徽大学 信息材料与智能感知安徽省实验室,合肥 230601)摘要 高超声速飞行器在飞行过程中,由于高温高压的作用,会在飞行器表面形成一层等离子体鞘层,飞行器及其等离子体鞘层相对于地面测控中心做高速相对运动,对目标物电磁回波会带来严重影响.此外,等离子体鞘层的时变特性,也会对目标物雷达回波进行调制,使测控中心难以识别、跟踪目标.本文利用 Lorentz-FDTD方法研究了相对运动等离子体鞘层与电磁波之间的相互作用,并分析了等离子体鞘
3、层的相对运动特性和时变特性对电磁波造成的影响,发现运动的时变等离子体除了对电磁波造成多普勒频移外,还会对入射波频谱进行调制.关键词时域有限差分;Lorentz 变换;相对运动;时变等离子体中图分类号O441/V11文献标志码A文章编号1005-0388(2023)01-0164-09DOI 10.12265/j.cjors.2022003Interaction between relative moving plasma plate and electromagnetic waveMAN Liang1,2DENG Haochuan1,2BO Yong3XIAO Zhihe1,2YANG Lix
4、ia3*(1.The Second Academy of China Aerospace Science and Industry Corporation,Beijing 100854,China;2.Science and Technologyon Electromagnetic Scattering Laboratory,Beijing 100854,China;3.Information Materials and IntelligentSensing Laboratory of Anhui Province,Anhui University,Hefei 230601,China)Abs
5、tractDuring the flight of hypersonic aircraft,due to the impact of high temperature and high pressure,alayer of plasma sheath is formed on the surface of the aircraft.With the aircraft and its plasma sheath motion at highspeed relative to the ground measurement and control center,a serious impact on
6、 the electromagnetic echo of the targetwill be caused.In addition,the time-varying characteristics of plasma sheath also modulate the radar echo of the target,making it difficult for the center to identify and track the target.In this paper,Lorentz-FDTD method is used to studythe interaction between
7、 the relatively moving plasma sheath and the electromagnetic wave,and the influence of therelative motion and time-varying characteristics of the plasma sheath on the electromagnetic wave is analyzed.It isfound that the moving time-varying plasma not only causes Doppler shift to the electromagnetic
8、wave but alsomodulates the spectrum of incident wave.KeywordsFDTD;Lorentz transformation;relative motion;time varying plasma 引言随着航空航天及深空探测技术的不断发展,飞行器飞行速度也在不断地刷新记录,当飞行速度到达一定程度后,高速运动对电磁波传播、散射产生的相对论效应影响不能忽略.近年来,高速飞行目标物 收稿日期:2022-01-04资助项目:国家自然科学基金(62071003,41874174,61901004);国家国防科工局基础研究重点项目(稳定支持项目);安徽高
9、校自然科学研究项目(KJ2020A0024);安徽省重点科研平台协同创新项目(GXXT-2020-050);安徽省高校协同创新计划项目(GXXT-2021-028);脉冲功率激光技术国家重点实验室开放研究基金(AHL 2020 KF04);计算智能与信号处理教育部重点实验室开放课题(2020A008)通信作者:杨利霞 E-mail: 第 38 卷第 1 期电波科学学报Vol.38,No.12023 年 2 月CHINESE JOURNAL OF RADIO SCIENCEFebruary,2023 的电磁散射特性研究引起了人们越来越多的关注.高超声速飞行器飞行过程中,飞行器的表面会覆盖一层相对
10、运动的等离子体流场,高速的等离子体流场会对飞行目标的探测和识别带来困难.高速运动目标电磁散射特性仿真计算研究主要有两种途径:Lorentz 坐标系变换方法1-6和相对论边界条件方法7-9.Lorentz 坐标系变换法通过 Lorentz 变换将所要求解的电磁波传输、散射问题(麦克斯韦方程、边界条件)变换到与目标物拥有相同速度的运动坐标系中进行计算,此时由于运动目标和运动坐标系之间相对静止,电磁场值可以利用传统的时域有限差分(finite-difference time-domain,FDTD)方法进行求解;计算得到的场值结果再进一步经过相对论 Lorentz逆变换,回到与地面测控系统相对静止的
11、实验坐标系中.基于 Lorentz 变换建立的 Lorentz-FDTD 算法1-6能够和经典的 FDTD 方法保持相同的精度和稳定性,是分析高速运动目标电磁散射问题的有效途径之一.相对论边界条件方法在处理运动目标电磁散射计算时,需要不断地更新边界处的场值,易造成计算结果的不收敛,但在处理非线性运动的目标物电磁散射时有一定的优势7-9.文献 1 通过在静止坐标系和运动坐标系之间做 Lorentz 变换计算了运动的金属锥目标的雷达散射截面,并研究了微运动对雷达散射的影响.文献 2 将 Lorentz-FDTD 方法应用到三维复杂运动目标雷达回波信号的模拟中,并与理论计算结果进行对比,发现结果吻合
12、良好.文献 3利用 Lorentz-FDTD 方法计算了入射波为高斯脉冲,相对运动的多层介质板的反射波频谱.文献 4 在仿真高速运动一维导体目标散射时,发现反射波的幅度和频率由于相对论效应都受到调制作用.在此基础上,文献 5-6又研究了高速运动的介质平板对反射波幅度和频率的调制效果,并与解析方法的计算结果进行了对比.随着国内外专家对高超声速飞行器外等离子体绕流场研究的不断深入,发现流场中等离子体电子密度与流场的状态紧密相关,当流场发生湍动时,电子密度也会发生抖动.20 世纪 50 年代末,Kistler 就发现了等离子体湍流中温度抖动现象10,Demetriade利用温度抖动计算了电子密度的抖
13、动11.文献 12-13从仿真模拟和风洞实验的角度验证了飞行器再入过程中等离子体鞘层的时变特性,发现等离子体湍流会造成电子密度的振荡.文献 14 又指出等离子体流场中电子密度的时变频率近似等于湍流的猝发频率,约为 20 60 kHz,但具体服从的抖动过程和规律尚不明确.文献 15-16 分别通过仿真和实验研究发现快速生成的等离子体会对入射的电磁波造成频率漂移现象.文献 17 利用火箭推进器产生动态等离子体喷焰,研究了时变的等离子体对电磁波的寄生调制效应,发现被寄生调制的载波信号频点出现了偏移,偏移量约为等离子体的时变频率.综上可知,国内外研究者们已经对运动目标物与电磁波之间的相互作用做了一些研
14、究,对飞行器等离子体鞘层的时变特性也有了一定的认识.本文在前人工作的基础上进一步研究了运动的等离子体与电磁波之间的相互作用,并考虑等离子体鞘层运动和时变两种特殊情况共同作用下对电磁波传输特性带来的影响.1 理论推导 1.1 相对运动等离子体介质建模kv当电磁波在运动各向同性介质(等离子体)中传播时,电磁波频率会发生多普勒频移,假设电磁波传播矢量 与目标运动速度 之间的夹角为,求得频移后的频率为=(1cos).(1):=1/12=v/c vvcvkk=式中;,为 的模值,为真空中的光速1-3.当地面测控中心发射电磁波方向与目标运动方向一致时,与 位于同一方向,会发生频率下移或红移;当目标物向地面
15、测控中心靠近时,会发生频率上移或蓝移;当接收机垂直于 运动时,有横向多普勒位移,此时属于纯相对论效应.常用 Appleton 公式来描述磁化冷等离子体的色散关系,根据文献 18-19 可得r()=12p2+2enjen2p2+2en.(2)rpen式中:为等离子的相对介电常数;为等离子体截止频率;为电子与中性粒子之间的碰撞频率.由于相对论效应,当频率为 的入射电磁波进入运动的等离子体中时,入射频率变为,对应的运动等离子体相对介电常数变为r()=12p2+2enjen2p2+2en.(3)r()为使用色散介质移位算子 FDTD 方法对运动等离子体进行计算,此处将运动等离子体的相对介电常数写成有理
16、分式的形式:r()=0r()=0Nn=0pn(j)nNn=0qn(j)n,(4)则沿 x 方向本构方程的时域形式可写为Dx(t)=0r(t)Ex(t).(5)0式中,为真空中的介电常数.C=(1cos)令,则第 1 期满良,等:相对运动等离子体平板与电磁波相互作用研究165 j=Cj.(6)运动等离子体的相对介电常数在频域和时域可以分别表示为:r()=Nn=0pnCn(j)nNn=0qnCn(j)n,(7)r(t)=Nn=0pnCn(t)nNn=0qnCn(t)n.(8)将式(8)带入本构方程(5)中,可得电场本构关系的时域表示:|Nn=0qn(t)n|Dx(t)=0|Nn=0pn(t)n|E
17、x(t).(9)qn=qnCnpn=pnCn式中:;.下面推导离散时域的时间导数19,设函数y(t)=f(t)t.(10)(n+1/2)t对该函数在 处进行中心差分近似,可得yn+1+yn2=fn+1 fnt.(11)zt引入离散时域的移位算子:ztfn=fn+1,(12)则yn=2tzt1zt+1fn.(13)(/t)h(zt1)/(zt+1)由式(13)可得从时域微分算子过渡到离散时域移 位 算 子 的 表 达 式,对 于m 阶微分算子有(t)mh(zt1)(zt+1)m.(14)h=2/t式中,.将式(14)带入本构关系式(9),可得|Nl=0ql(2tzt1zt+1)l|Dnx=0|N
18、l=0pl(2tzt1zt+1)l|Enx.(15)(zt+1)N式(15)两边同乘以得|Nl=0ql(2t)l(zt+1)Nl(zt1)l|Dnx=0|Nl=0pl(2t)l(zt+1)Nl(zt1)l|Enx.(16)将式(16)求和形式展开可得|q0+q12t+q2(2t)2|z2t+|2q02q2(2t)2|zt+|q0q12t+q2(2t)2|Dnx=|p0+p12t+p2(2t)2|z2t+|2p02p2(2t)2|zt+|p0 p12t+p2(2t)2|0Enx,(17)从而得出运动色散介质中电场的更新方程:En+1x=1b0a0(Dn+1x0)+a1(Dnx0)+a2(Dn1x
19、0)b1Enxb2En1x.(18)式中:|a0=q0+q12t+q2(2t)2a1=2q02q2(2t)2a2=q0q12t+q2(2t)2;(19)|b0=p0+p12t+p2(2t)2b1=2p02p2(2t)2b2=2p02p2(2t)2.(20)通过式(1)(20)即可完成 Dx到 Ex的迭代过程,获得考虑相对论效应的电场的更新方程.电场在y 方向和 z 方向的更新方程可以通过类似的方法推导出来.1.2 运动坐标系反推静止坐标系中的电磁场麦克斯韦方程在时空变换下是 Lorentz 协变的,即无论是在相对运动还是相对静止的坐标系中,麦克斯韦方程组的形式都是不变的,而场量会发生变化,两种
20、坐标系中场量之间满足 Lorentz 变换式20.根据相对论中的协变性原理,即在给定坐标系下的分量在坐标变换中满足协同变换规律,可知运动坐标系中等离子体的旋度方程和静止坐标系中旋度方程形式是一样的,因此在运动坐标系中 FDTD迭代式在形式上基本一样.图 1 为 Lorentz-FDTD 方法设计流程图,步骤如下:静止直角坐标系运动直角坐标系相对论协变性运动等离子体中的FDTD迭代公式Lorentz变换运动坐标系中的场量静止坐标系中的场量H=E=DtBtD=0rEE=0HH,D,E,B,rH,D,E,B,rH=E=DtBtD=0rEE=0H图 1 Lorentz-FDTD 方法计算运动目标电磁问
21、题流程图Fig.1 Flowchart of Lorentz-FDTD method forelectromagnetic problems of moving targets 166电波科学学报第 38 卷1)首先,根据相对论协变性(麦克斯韦方程形式不变性),将要求解的位于实验坐标系中的电磁问题转化到与目标物有着相同运动速度的相对静止的运动坐标系中;2)其次,在运动坐标系中,目标物处于相对静止的状态,而等离子体则处于相对运动的状态,通过公式(8)对运动坐标系中等离子体介电常数进行处理,并利用Lorentz-FDTD 方法计算运动坐标系中的电磁场;3)最后,对于上一步得到的运动坐标系中的电磁场
22、值,进一步利用 Lorentz 逆变换,得到静止坐标系(相对地面测控中心静止)中电磁场值,并进行后续处理计算.运动坐标系下麦克斯韦方程组的两个旋度方程可以表示为:Dt=H,(21)Ht=10E.(22)等离子体中本构关系如式(23)所示:D()=0r()E().(23)EDHBEDHB对式(21)(23)进行差分离散,得到在运动坐标系中 FDTD 的迭代公式.通过 FDTD 计算可以得到运动坐标系中的电场、电位移矢量以及磁场、磁通量,然后利用 Lorentz 变换得到相对静止坐标系中电场、电位移矢量和磁场、磁通量.下面推导运动坐标系到静止坐标系中场量变换关系的 FDTD 的迭代式.电场的 Lo
23、rentz 变换为|ExEyEz|=|ExEyEz|+(1)Exvx+Eyvy+Ezvzv2v.(24)式中,vx、vy、vz分别表示沿 x、y、z 方向的速度分量.Ex对式(24)中进行整理后可得Ex=Ex+(vyBzvzBy)+(1)Exvx+Eyvy+Ezvzv2vx,(25)ExEx进而得到离散时域与之间的关系式:En+1x=v2(v2y+v2z)+v2xEn+1x+(+1)vx(v2y+v2z)+v2x(Enyvy+Enzvz)+v2(v2y+v2z)+v2x(vzBn+12yvyBn+12z).(26)ExEyEzByBz式(26)中,和、不在同一个空间节点位置,可以通过空间插值处
24、理.同理我们可以再利用BxBxDxDxHxHxLorentz 变换公式求解与、与、与之间的 FDTD 迭代式,如式(27)、(30)、(33)所示.|BxByBz|=|BxByBz|c2|xyzvxvyvzBxByBz|+(1)Exvx+Eyvy+Ezvzv2v.(27)对式(27)进行整理可得Bx=Bxc2(vyEzvzEy)+(1)Bxvx+Byvy+Bzvzv2vx.(28)Bx在离散时域,可以得到的 FDTD 迭代式:Bn+12x=v2(v2y+v2z)+v2xBn+12x+(1)vx(v2y+v2z)+v2x(Bn+12yvy+Bn+12zvz)+v2(v2y+v2z)c2+v2xc
25、2(vyEnzvzEny).(29)|DxDyDz|=DxDyDz+c2|xvxHxyvyHyzvzHz|+(1)Dxvx+Dyvy+Dzvzv2v.(30)对式(30)进行整理后可得Dx=Dx+c2(vyHzvzHy)+(1)Dxvx+Dyvy+Dzvzv2vx.(31)在离散时域,可以得到 Dx的 FDTD 迭代式:Dnx=v2(v2y+v2z)+v2xDnx+(1)(v2y+v2z)+v2x(Dn+12yvy+Dn1zvz)+v2(v2y+v2z)c2+v2xc2(vzHn12yvyHn12z).(32)|HxHyHz|=|HxHyHz|xvxDxyvyDyzvzDz|+(1)Hxvx+
- 配套讲稿:
如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。