电磁波在等离子体中传播的二维FDTD数值计算课件.ppt

1、江中和 胡希伟 刘明海Huazhong University ofScience and Technology主题p电磁波在等离子体中的传播p时域有限差分法（FDTD）的基本思想 p电磁波在等离子体传播的FDTD模拟计算方法p大气压下二维电磁波在等离子体中的传播特性 电磁波在等离子体中的传播p一维波动方程以及耦合电子运动方程求解p二维情形下波动方程求解的困难p应用FDTD进行电磁波在等离子体中传播的数值模拟计算一维波动方程以及耦合电子运动方程求解eeueEutm 222222140,eeEEJJen uxctct 波动方程：电子流体方程：差分方程的边界条件p金属边界条件p入射边界条件(0,)0

2、Et 当入射边界取距离入射点足够远时，使得反射回来的电磁波在等离子体边界面附近经过自由空间一个或者几个周期还没有到达起始点，可以不考虑入射边界对电磁波传播的影响。二维波动方程求解的困难22222222140zzzzeEEEunezyctctp差分格式要采用二维的差分格式，这决定差分得到的方程较为复杂。p为了避免边界截断造成的非物理型反射影响，使得整个计算域的空间要取得足够大，这对计算机软硬件都存在很大的依赖，最明显的是计算机内存要足够大。22222222140yyyyeEEEunezyctctFDTD算法进行等离子体数值模拟pFDTD算法，其空间节点采用Yee元胞的方法，电场和磁场节点空间与时

3、间上都采用交错抽样，因而使得麦克斯韦旋度方程离散后构成显式差分方显式差分方程程，相比较与前面的波动方程求解，计算得到大大简化。p由于FDTD采用吸收边界条件的方法，使得计算可以在有限的空间范围内有限的空间范围内进行，这样就可以降低程序对计算机硬件的要求。时域有限差分法（FDTD）基本思想 FDTD算法是K.S.Yee于1966年提出的、直接对麦克斯韦方程作差分处理、来解决电磁脉冲在电磁介质中传播和反射问题的算法。基本思想是：FDTD计算域空间节点采计算域空间节点采用用Yee元胞的方法，同时电场和磁场节点空间与时间上都元胞的方法，同时电场和磁场节点空间与时间上都采用交错抽样；把整个计算域划分成包

4、括散射体的总场采用交错抽样；把整个计算域划分成包括散射体的总场区以及只有反射波的散射场区，这两个区域是以连接边区以及只有反射波的散射场区，这两个区域是以连接边界相连接，最外边是采用特殊的吸收边界，同时在这两界相连接，最外边是采用特殊的吸收边界，同时在这两个边界之间有个输出边界，用于近、远场转换；在连接个边界之间有个输出边界，用于近、远场转换；在连接边界上采用连接边界条件加入入射波，从而使得入射波边界上采用连接边界条件加入入射波，从而使得入射波限制在总场区域；在吸收边界上采用吸收边界条件，尽限制在总场区域；在吸收边界上采用吸收边界条件，尽量消除反射波在吸收边界上的非物理性反射波量消除反射波在吸收

5、边界上的非物理性反射波。FDTD示意图.Z ZY YX X吸收边界吸收边界输出边界输出边界连接边界连接边界等离子体层等离子体层金属板层金属板层二维TM波的麦克斯韦方程组yzHExtzxEHyt yxzzHHEJxyt二维TM波的差分格式11,1,2211,22nnnnz j kz j kxxj kj kEEtHHy111,2211,22nnnnz jkz j kyyjkjkEEtHHx1111222211111,122222,()nnnnyyxxjkjkj kj knnnz j kz j kz j kHHHHtEEJxy吸收边界条件 由于计算机容量的限制，FDTD只能在有限区域进行，为了能够模

6、拟开域的电磁散射过程，在计算区域的外边要求给出吸收边界条件，我们采用被广为应用的二阶Mur吸收边界条件，TM波其吸收边界条件在边界上为：y102zzzbEEcHxcty 下 边 界：yb102zzzEEcHxcty上边界：a102zzzxEEcHxcty右 边 界：102zzzxaEEcHxcty 左 边 界：二维TM角点的吸收条件在二维矩形FDTD计算域内的角点，上述的吸收边界条件无法应用，根据一阶Mur吸收条件可以得到在角点时吸收边界条件的差分格式：n+1nn+1n000000002E(,)E(1,1)E(1,1)E(,)2zzzzc ti jijiji jc t 式中t、分别为时间步长与

7、空间步长。二维TM波的连接边界条件连接边界条件的目的就是要使得在散射场区域内只存在反射波，而在总场区域则入射波与反射波同时存在，因此连接边界条件就是对于总场区域对于总场区域的节点，如果其差分格式中有散射区域的节点，的节点，如果其差分格式中有散射区域的节点，需要在散射场节点加上入射波，而对于散射场区需要在散射场节点加上入射波，而对于散射场区域的节点在差分格式中有总场节点，则总场节点域的节点在差分格式中有总场节点，则总场节点要扣除入射波值要扣除入射波值。因此对于前面提到的TM波的麦克斯韦方程的差分格式在连接边界附近作上述修正。此即为连接边界条件。平面波的加入.Z ZY YX X吸收边界吸收边界输出

8、边界输出边界连接边界连接边界金属板层金属板层.p+1pp 入射波是加在连接边界附近上，由于是平面波，故而采用一维FDTD随时间逐步推进地在总场区引进入射波，这样可以减小散射场入射波的泄漏。对于与底边成角平面波，可以通过连接边界上向一维平面波作垂线，根据平面波特性以及插值方法得到边界上的入射波。电磁波在等离子体传播的FDTD模拟计算方法 虽然目前的FDTD算法可以解决给定介电系数和磁导率的介质对入射电磁波的响应（包括传播、反射、吸收和透射），但无法正确地反映大气环境中等离子体对入射电磁波的响应。在这种等离子体中，电子与中性原子（分子）的弹性碰撞频率将等于或大于入射电磁波的频率，因此仅用复介电常数

9、已经不能真实地反映等离子体的介电性质。使用将等离子体视为复介电常数介质的传统FDTD算法、以及其它时域算法已经很难保证计算结果的可靠性与合理性。正确的作法是将等离子体的电子流体运动方程直接将等离子体的电子流体运动方程直接和麦克斯韦方程组耦合求解，由此得出电磁场在这种特和麦克斯韦方程组耦合求解，由此得出电磁场在这种特殊等离子体中的传播特性殊等离子体中的传播特性。麦克斯韦方程以及电子流体方程组()zezzzJen uJE yzHExtzxEHyt yxzzHHEJxytzzzeueEutm 二维TM波的差分格式11,1,2211,22nnnnz j kz j kxxj kj kEEtHHy111,

10、2211,22nnnnz jkz j kyyjkjkEEtHHx1111222211111,122222,()nnnnyyxxjkjkj kj knnnz j kz j kz j kHHHHtEEJxy111222,nnnz j ke j kz j kJenu 1122,12(2)(2)nnnz j kz j kz j kee tut uEtm数值实验 应用FDTD方法进行数值计算，要验证计算结果的有效性，需要进行数值实验，数值实验主要有以下几个方面：p连接边界上加入射波，不加入等离子体以及金属板散射体，计算散射场区以及总场区的电场强度。p连接边界上加入射波，有金属板散射体，不加入等离子体计算

11、散射场区以及总场区的电场强度。FDTD算例参数p根据FDTD思想进行算例计算，其中电磁波频率在2GHz至20GHz范围，空间步长取h=0.25mm，时间步长根据计算稳定性要(th/c)t=0.5h/c(c为真空中光速)，吸收边界为830160，连接边界为810140，金属板层尺寸为80010，等离子体电子密度随等离子体厚度为双曲正切上升函数，等离子体厚度变化范围为10mm至30mm。FDTD数值实验结果 p当不加入等离子体以及金属板散射体，散射场区电场场强很小可看成0，总场区的电场强度随时间变化与入射波随时间变化一致。说明连接边界条件起到很好的防止入射波泄漏到散射场作用。p有金属板散射体，不加

12、入等离子体计算散射场区以及总场区的电场强度，数值实验结果表明，在垂直入射时沿x轴方向，在靠近金属板两端面附近的散射波振幅较小，而在Y轴附近电场的散射波振幅基本保持不变，随着入射角变小，靠近入射端的散射波振幅变小区域变大，而出射端的散射波振幅变小区域变小，这是由于金属板全反射引起的。FDTD示意图.Z ZY YX X吸收边界吸收边界输出边界输出边界连接边界连接边界等离子体层等离子体层金属板层金属板层数值实验计算结果900入射角020040060080010001200-1.2-1.0-0.8-0.6-0.4-0.20.00.20.40.60.81.01.2 Electric FieldTime0

13、20040060080010001200-1.2-1.0-0.8-0.6-0.4-0.20.00.20.40.60.81.01.2 Electric FieldTime(-400,-95)(-350,-95)(-300,-95)(-250,-95)(-200,-95)数值实验计算结果600入射角020040060080010001200-1.2-1.0-0.8-0.6-0.4-0.20.00.20.40.60.81.01.2 Electric FieldTime020040060080010001200-1.2-1.0-0.8-0.6-0.4-0.20.00.20.40.60.81.01.2

14、Electric FieldTime(-400,-95)(-350,-95)(-300,-95)(-250,-95)(-200,-95)(-150,-95)(-100,-95)数值实验结果450 入射角020040060080010001200-1.2-1.0-0.8-0.6-0.4-0.20.00.20.40.60.81.01.2 Electric FieldTime(-400,-95)(-350,-95)(-300,-95)(-250,-95)(-200,-95)(-100,-95)(0,-95)(250,-95)(400,-95)电磁波在等离子体中传播特性 当等离子体厚度为30mm时，在

15、不同入射角情况 下散射场中电场强度随入射波频率衰减图02468101214161820225101520253035 attenuation(dB)f(GHz)30mm90degree 30mm60degree 30mm45degree电磁波在等离子体中传播特性0246810121416182022051015202530354045 attenuation(dB)f(GHz)20mm90degree 20mm60degree 20mm45degree 当等离子体厚度为20mm时，在不同入射角情况 下散射场中电场强度随入射波频率衰减图024681012141618202201234567 at

16、tenuation(dB)f(GHz)10mm90degree 10mm60degree 10mm45degree电磁波在等离子体中传播特性当等离子体厚度为10mm时，在不同入射角情况下散射场中电场强度随入射波频率衰减图 电磁波在等离子体中传播特性0246810121416182022051015202530 attenuation(dB)f(GHz)30mm90degree 20mm90degree 10mm90degree当入射角为90时，在不同等离子体厚度情况下散射场中电场强度随入射波频率衰减图0246810121416182022051015202530354045 attenuati

17、on(dB)f(GHz)30mm60degree 20mm60degree 10mm60degree当入射角为60时，在不同等离子体厚度情况下散射场中电场强度随入射波频率衰减图电磁波在等离子体中传播特性电磁波在等离子体中传播特性024681012141618202205101520253035 attenuation(dB)f(GHz)30mm45degree20mm45degree10mm45degree当入射角为45时，在不同等离子体厚度情况下散射场中电场强度随入射波频率衰减图电磁波在等离子体中传播特性 当等离子体厚度为20mm时，电子碰撞频率为23.6GHz在不同入射角情况下散射场中电场

18、强度随入射波频率衰减图024681012141618202202468101214161820 Attenuation(dB)f(GHz)20mm90degree 20mm60degree 20mm45degree电磁波在等离子体中传播特性 p在大气条件下，电磁波在一定厚度的等离子体传播，在频率较低时（小于10GHz），随着斜入射角度的变化反射的电磁波衰减相差不大，当频率较高时斜入射角度越小，电磁波衰减的越大，p等离子体的厚度对电磁波反射的衰减存在较大的影响，在等离子体厚度为10mm时，电磁波衰减的较小，而在等离子体厚度为20mm和30mm时存在较大衰减峰值，峰值的位置随等离子体厚度以及斜入射角变化而发生改变。这是由于当等离子体厚度与其波长相比拟时等离子体对电磁波就有比较明显的吸收效果，斜入射角越大，电磁波在等离子体中传播的距离就越大。p电子碰撞频率对电磁波在等离子体中的吸收有较大的影响。Huazhong University ofScience and Technology