MAGIC第一讲分解.pptVIP

* * Magic 软件是一种粒子模拟(PIC 全称 particle-in-cell)软件，主要用于模拟电真空器件中电磁波与空间电荷的互作用过程；计算和分析有空间电荷存在的复杂电磁问题。实际上，该软件提供了一种在计算机上对电真空器件进行模拟实验的方法。 Magic软件的概述 基本理论基础 Magic软件中，求解一个电磁学问题的基本方法表示如下： Maxwell方程组 Lorenz力学方程 电流连续性方程 具体边界条件 具体问题中的 电场和磁场 有源情况 注意：Magic软件中求解Maxwell方程组是采用时域有限差分的方法，简称FDTD 英文表示：finite-difference-time-domain Maxwell 方程 Lorenz力学方程 电流连续性方程 Maxwell 方程的分量形式 差分法原理 前向差分 后向差分 中心差分 二维电磁场泊松方程的差分格式 设 在一个由边界C限定的二维场域D内满足泊松方程 首先将场域D离散化，如图a；场域内部某节点0附近的各节点如图b所示（步长相等） D x y (i,j) (i,j-1) (i+1,j) (i-1,j) (i,j+1) (i+1,j+1) (i-1,j+1) (i-1,j-1) (i+1,j-1) 0 1 2 3 4 h h h h 图(a) 图(b) X (i) y (j) 以 分别代表在节点0、1、2、3、4处的 值 采用中心差分公式： 得差分方程(用节点表示)： 从而，连续的泊松方程就变成了离散的差分方程。只要求出各网格点上的 也就求出了泊松方程 值 有限差分法的原理： 把电磁场连续域内的问题变为离散系统的问题，即用各离散点上的数值 解来逼近连续场域内的真实值。 在求解Maxwell方程组时，将对时间的偏微分也进行差分，所以又叫时域有限 差分法。同时，要对三维空间差分，因此要采用三维网格。 Yee网格的划分 1966年，K.S.Yee提出了一种如右图 的差分格式，网格节点与一组相应 的整数标号一一对应，同时将空间任 一矩形网格上的E和H的6个分量 如右图放置。 z x y B1 B3 (i,j,k+1) (i+1,j,k+1) (i+1,j,k) (i+1,j+1,k) (i,j+1,k) (i,j+1,k+1) 电场和磁场各分量在空间的取值点被交叉地放置，使得在每个坐标平面上,每个电场分量的四周由磁场分量环绕，同时每个磁场分量四周由电场分量环绕。这样的电磁场空间配置符合电磁场的基本规律——Faraday电磁感应定律和Ampere环流定律。因而也符合电磁波在空间传播的规律。 Yee网格的特点 “蛙跳”格式 为了得到精度较高的解，Yee将E和H在时 间上相差半个步长交替计算。即用E的n时 刻值和B的n＋1/2时刻的值，得到E的n＋1 时刻的值。同理，用E的n＋1时刻的值和B 的n＋1/2时刻的值，可以得到B的n＋3/2时 刻的值。 我们把这种计算方法叫“蛙跳”格式。 离散的电磁场方程 离散方程的特点 由上面的6个方程，可以看出在时间和空间上都采用了中心差分，任一网络点上的电场分量只与它上一个时间步的值及四周环绕它的磁场分量有关。同样，任一网格点上的磁场分量也只与它上一时间步时的值及四周环绕它的电场分量有关。 计算电磁学 王秉中 著 （科学出版社） 电磁波时域有限差分方法 葛德彪 著（西安电子科技大学出版社） 参考书： 场 粒子受力 粒子动量，位置 电流密度 1.求解推动粒子: 已知n时刻的粒子位置, 电磁场和n－1/2时刻的粒子速度，求解洛仑兹方程，得到n＋1/2时刻的粒子动量。再由n＋1/2时刻的粒子动量和n时刻的粒子位置求解运动方程,得到n＋1时刻的粒子位置。 自洽场的求解过程 2.由粒子的坐标和速度求得新的电荷密度和电流密度; 3.由电流密度求解Maxwell方程得到电磁场。 1）已知n时刻的电场强度，n＋1/2时刻的磁场强度和电流密度J，求出n+1时刻的电场强度。 2）由n+1时刻的电场强度和n＋1/2时刻的磁场强度，求出n＋3/2时刻的磁场强度。 在Magic软件中我们要做的事情： a) 给出一个具体模型 （在此举例说明） c) 给具体物体赋材料特性 其它具体计算工作完全由Magic软件完成！！ e) 给出具体算法 f) 设定观察命令 b) 在模拟区域内划分网格 g) 分析结果 d) 定义边界条件 编写Magic程序的基本流程： 打开编辑器 定义要模拟的相关信息 (header命令完成) 定义坐标系 (System 完成) 定义变量 定义物体 (point, line, area, volume命令完成） 划分网格 (mark, a