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

江中和胡希伟刘明海HuazhongUniversityofScienceandTechnology华中科技大学电气与电子工程学院主题电磁波在等离子体中的传播时域有限差分法（FDTD）的基本思想电磁波在等离子体传播的FDTD模拟计算方法大气压下二维电磁波在等离子体中的传播特性华中科技大学电气与电子工程学院电磁波在等离子体中的传播一维波动方程以及耦合电子运动方程求解二维情形下波动方程求解的困难应用FDTD进行电磁波在等离子体中传播的数值模拟计算华中科技大学电气与电子工程学院一维波动方程以及耦合电子运动方程求解eeueEutm222222140,eeEEJJenuxctct波动方程：电子流体方程：华中科技大学电气与电子工程学院差分方程的边界条件金属边界条件入射边界条件(0,)0Et当入射边界取距离入射点足够远时，使得反射回来的电磁波在等离子体边界面附近经过自由空间一个或者几个周期还没有到达起始点，可以不考虑入射边界对电磁波传播的影响。华中科技大学电气与电子工程学院二维波动方程求解的困难22222222140zzzzeEEEunezyctct差分格式要采用二维的差分格式，这决定差分得到的方程较为复杂。为了避免边界截断造成的非物理型反射影响，使得整个计算域的空间要取得足够大，这对计算机软硬件都存在很大的依赖，最明显的是计算机内存要足够大。22222222140yyyyeEEEunezyctct华中科技大学电气与电子工程学院FDTD算法进行等离子体数值模拟FDTD算法，其空间节点采用Yee元胞的方法，电场和磁场节点空间与时间上都采用交错抽样，因而使得麦克斯韦旋度方程离散后构成显式差分方程，相比较与前面的波动方程求解，计算得到大大简化。由于FDTD采用吸收边界条件的方法，使得计算可以在有限的空间范围内进行，这样就可以降低程序对计算机硬件的要求。华中科技大学电气与电子工程学院时域有限差分法（FDTD）基本思想FDTD算法是K.S.Yee于1966年提出的、直接对麦克斯韦方程作差分处理、来解决电磁脉冲在电磁介质中传播和反射问题的算法。基本思想是：FDTD计算域空间节点采用Yee元胞的方法，同时电场和磁场节点空间与时间上都采用交错抽样；把整个计算域划分成包括散射体的总场区以及只有反射波的散射场区，这两个区域是以连接边界相连接，最外边是采用特殊的吸收边界，同时在这两个边界之间有个输出边界，用于近、远场转换；在连接边界上采用连接边界条件加入入射波，从而使得入射波限制在总场区域；在吸收边界上采用吸收边界条件，尽量消除反射波在吸收边界上的非物理性反射波。华中科技大学电气与电子工程学院FDTD示意图.ZYX吸收边界输出边界连接边界等离子体层金属板层华中科技大学电气与电子工程学院二维TM波的麦克斯韦方程组yzHExtzxEHytyxzzHHEJxyt华中科技大学电气与电子工程学院吸收边界条件由于计算机容量的限制，FDTD只能在有限区域进行，为了能够模拟开域的电磁散射过程，在计算区域的外边要求给出吸收边界条件，我们采用被广为应用的二阶Mur吸收边界条件，TM波其吸收边界条件在边界上为：y102zzzbEEcHxcty下边界：＋yb102zzzEEcHxcty上边界：＋a102zzzxEEcHxcty右边界：＋＋102zzzxaEEcHxcty左边界：华中科技大学电气与电子工程学院二维TM角点的吸收条件在二维矩形FDTD计算域内的角点，上述的吸收边界条件无法应用，根据一阶Mur吸收条件可以得到在角点时吸收边界条件的差分格式：n+1nn+1n000000002E(,)E(1,1)E(1,1)E(,)2zzzzctijijijijct式中Δt、Δδ分别为时间步长与空间步长。华中科技大学电气与电子工程学院二维TM波的连接边界条件连接边界条件的目的就是要使得在散射场区域内只存在反射波，而在总场区域则入射波与反射波同时存在，因此连接边界条件就是对于总场区域的节点，如果其差分格式中有散射区域的节点，需要在散射场节点加上入射波，而对于散射场区域的节点在差分格式中有总场节点，则总场节点要扣除入射波值。因此对于前面提到的TM波的麦克斯韦方程的差分格式在连接边界附近作上述修正。此即为连接边界条件。华中科技大学电气与电子工程学院平面波的加入.ZYX吸收边界输出边界连接边界金属板层..p+1p入射波是加在连接边界附近上，由于是平面波，故而采用一维FDTD随时间逐步推进地在总场区引进入射波，这样可以减小散射场入射波的泄漏。对于与底边成θ角平面波，可以通过连接边界上向一维平面波作垂线，根据平面波特性以及插值方法得到边界上的入射波。华中科技大学电气与电子工程学院电磁波在等离子体传播的FDTD模拟计算方法虽然目前的FDTD算法可以解决给定介电系数和磁导率的介质对入射电磁波的响应（包括传播、反射、吸收和透射），但无法正确地反映大气环境中等离子体对入射电磁波的响应。在这种等离子体中，电子与中性原子（分子）的弹性碰撞频率将等于或大于入射电磁波的频率，因此仅用复介电常数已经不能真实地反映等离子体的介电性质。使用将等离子体视为复介电常数介质的传统FDTD算法、以及其它时域算法已经很难保证计算结果的可靠性与合理性。正确的作法是将等离子体的电子流体运动方程直接和麦克斯韦方程组耦合求解，由此得出电磁场在这种特殊等离子体中的传播特性。华中科技大学电气与电子工程学院麦克斯韦方程以及电子流体方程组()zezzzJenuJEyzHExtzxEHytyxzzHHEJxytzzzeueEutm华中科技大学电气与电子工程学院二维TM波的差分格式11,1,2211,,22nnnnzjkzjkxxjkjkEEtHHy111,,2211,,22nnnnzjkzjkyyjkjkEEtHHx1111222211111,,,,122222,,,()nnnnyyxxjkjkjkjknnnzjkzjkzjkHHHHtEEJxy111222,,,nnnzjkejkzjkJenu1122,,,12[(2)](2)nnnzjkzjkzjkeetutuEtm华中科技大学电气与电子工程学院数值实验应用FDTD方法进行数值计算，要验证计算结果的有效性，需要进行数值实验，数值实验主要有以下几个方面：连接边界上加入射波，不加入等离子体以及金属板散射体，计算散射场区以及总场区的电场强度。连接边界上加入射波，有金属板散射体，不加入等离子体计算散射场区以及总场区的电场强度。华中科技大学电气与电子工程学院FDTD算例参数根据FDTD思想进行算例计算，其中电磁波频率在2GHz至20GHz范围，空间步长取Δh=0.25mm，时间步长根据计算稳定性要(ΔtΔh/c)Δt=0.5Δh/c(c为真空中光速)，吸收边界为830×160，连接边界为810×140，金属板层尺寸为800×10，等离子体电子密度随等离子体厚度为双曲正切上升函数，等离子体厚度变化范围为10mm至30mm。华中科技大学电气与电子工程学院FDTD数值实验结果当不加入等离子体以及金属板散射体，散射场区电场场强很小可看成0，总场区的电场强度随时间变化与入射波随时间变化一致。说明连接边界条件起到很好的防止入射波泄漏到散射场作用。有金属板散射体，不加入等离子体计算散射场区以及总场区的电场强度，数值实验结果表明，在垂直入射时沿x轴方向，在靠近金属板两端面附近的散射波振幅较小，而在Y轴附近电场的散射波振幅基本保持不变，随着入射角变小，靠近入射端的散射波振幅变小区域变大，而出射端的散射波振幅变小区域变小，这是由于金属板全反射引起的。华中科技大学电气与电子工程学院FDTD示意图.ZYX吸收边界输出边界连接边界等离子体层金属板层华中科技大学电气与电子工程学院数值实验计算结果－900入射角020040060080010001200-1.2-1.0-0.8-0.6-0.4-0.20.00.20.40.60.81.01.2ElectricFieldTime020040060080010001200-1.2-1.0-0.8-0.6-0.4-0.20.00.20.40.60.81.01.2ElectricFieldTime(-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.2ElectricFieldTime020040060080010001200-1.2-1.0-0.8-0.6-0.4-0.20.00.20.40.60.81.01.2ElectricFieldTime(-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.2ElectricFieldTime(-400,-95)(-350,-95)(-300,-95)(-250,-95)(-200,-95)(-100,-95)(0,-95)(250,-95)(400,-95)华中科技大学电气与电子工程学院电磁波在等离子体中传播特性当等离子体厚度为30mm时，在不同入射角情况下散射场中电场强度随入射波频率衰减图02468101214161820225101520253035attenuation(dB)f(GHz)30mm90degree30mm60degree30mm45degree华中科技大学电气与电子工程学院电磁波在等离子体中传播特性0246810121416182022051015202530354045attenuation(dB)f(GHz)20mm90degree20mm60degree20mm45degree当等离子体厚度为20mm时，在不同入射角情况下散射场中电场强度随入射波频率衰减图华中科技大学电气与电子工程学院024681012141618202201234567attenuation(dB)f(GHz)10mm90degree10mm60degree10mm45degree电磁波在等离子体中传播特性当等离子体厚度为10mm时，在不同入射角情况下散射场中电场强度随入射波频率衰减图华中科技大学电气与电子工程学院电磁波在等离子体中传播特性0246810121416182022051015202530attenuation(dB)f(GHz)30mm90degree20mm90degree10mm90degree当入射角为90时，在不同等离子体厚度情况下散射场中电场强度随入射波频率衰减图华中科技大学电气与电子工程学院0246810121416182022051015202530354045attenuation(dB)f(GHz)30mm60degree20mm60degree10mm60degree当入射角为60时，在不同等离子体厚度情况下散射场中电场强度随入射波频率衰减图电磁波在等离子体中传播特性华中科技大学电气与电子工程学院电磁波在等离子体中传播特性024681012141618202205101520253035atte