电磁波在等离子体中的传播.ppt

主题 电磁波在等离子体中的传播 时域有限差分法（FDTD）的基本思想 电磁波在等离子体传播的FDTD模拟计算方法 大气压下二维电磁波在等离子体中的传播特性 电磁波在等离子体中的传播 一维波动方程以及耦合电子运动方程求解 二维情形下波动方程求解的困难 应用FDTD进行电磁波在等离子体中传播的数值模拟计算 一维波动方程以及耦合电子运动方程求解 差分方程的边界条件 金属边界条件 二维波动方程求解的困难 FDTD算法进行等离子体数值模拟 时域有限差分法（FDTD）基本思想 FDTD示意图 二维TM波的麦克斯韦方程组 二维TM波的差分格式 吸收边界条件 由于计算机容量的限制，FDTD只能在有限区域进行，为了能够模拟开域的电磁散射过程，在计算区域的外边要求给出吸收边界条件，我们采用被广为应用的二阶Mur吸收边界条件，TM波其吸收边界条件在边界上为： 二维TM角点的吸收条件 二维TM波的连接边界条件 平面波的加入 电磁波在等离子体传播的FDTD模拟计算方法 麦克斯韦方程以及电子流体方程组 二维TM波的差分格式 数值实验 应用FDTD方法进行数值计算，要验证计算结果的有效性，需要进行数值实验，数值实验主要有以下几个方面： 连接边界上加入射波，不加入等离子体以及金属板散射体，计算散射场区以及总场区的电场强度。 连接边界上加入射波，有金属板散射体，不加入等离子体计算散射场区以及总场区的电场强度。 FDTD算例参数 FDTD数值实验结果 FDTD示意图 数值实验计算结果－900入射角 数值实验计算结果－600入射角 数值实验结果－450 入射角 电磁波在等离子体中传播特性 * 华中科技大学电气与电子工程学院 * 江中和 胡希伟 刘明海 Huazhong University of Science and Technology 波动方程： 电子流体方程： 入射边界条件 当入射边界取距离入射点足够远时，使得反射回来的电磁波在等离子体边界面附近经过自由空间一个或者几个周期还没有到达起始点，可以不考虑入射边界对电磁波传播的影响。 差分格式要采用二维的差分格式，这决定差分得到的方程较为复杂。 为了避免边界截断造成的非物理型反射影响，使得整个计算域的空间要取得足够大，这对计算机软硬件都存在很大的依赖，最明显的是计算机内存要足够大。 FDTD算法，其空间节点采用Yee元胞的方法，电场和磁场节点空间与时间上都采用交错抽样，因而使得麦克斯韦旋度方程离散后构成显式差分方程，相比较与前面的波动方程求解，计算得到大大简化。 由于FDTD采用吸收边界条件的方法，使得计算可以在有限的空间范围内进行，这样就可以降低程序对计算机硬件的要求。 FDTD算法是K.S.Yee于1966年提出的、直接对麦克斯韦方程作差分处理、来解决电磁脉冲在电磁介质中传播和反射问题的算法。基本思想是：FDTD计算域空间节点采用Yee元胞的方法，同时电场和磁场节点空间与时间上都采用交错抽样；把整个计算域划分成包括散射体的总场区以及只有反射波的散射场区，这两个区域是以连接边界相连接，最外边是采用特殊的吸收边界，同时在这两个边界之间有个输出边界，用于近、远场转换；在连接边界上采用连接边界条件加入入射波，从而使得入射波限制在总场区域；在吸收边界上采用吸收边界条件，尽量消除反射波在吸收边界上的非物理性反射波。 在二维矩形FDTD计算域内的角点，上述的吸收边 界条件无法应用，根据一阶Mur吸收条件可以得到 在角点时吸收边界条件的差分格式： 式中Δt、Δδ分别为时间步长与空间步长。 连接边界条件的目的就是要使得在散射场区域内只存在反射波，而在总场区域则入射波与反射波同时存在，因此连接边界条件就是对于总场区域的节点，如果其差分格式中有散射区域的节点，需要在散射场节点加上入射波，而对于散射场区域的节点在差分格式中有总场节点，则总场节点要扣除入射波值。因此对于前面提到的TM波的麦克斯韦方程的差分格式在连接边界附近作上述修正。此即为连接边界条件。 入射波是加在连接边界附近上，由于是平面波，故而采用一维FDTD随时间逐步推进地在总场区引进入射波，这样可以减小散射场入射波的泄漏。对于与底边成θ角平面波，可以通过连接边界上向一维平面波作垂线，根据平面波特性以及插值方法得到边界上的入射波。 虽然目前的FDTD算法可以解决给定介电系数和磁导率的介质对入射电磁波的响应（包括传播、反射、吸收和透射），但无法正确地反映大气环境中等离子体对入射电磁波的响应。在这种等离子体中，电子与中性原子（分子）的弹性碰撞频率将等于或大于入射电磁波的频率，因此仅用复介电常数已经不能真实地反映等离子体的介电性质。使用将等离子体视为复介电常数介质