静电模型PIC方法的Matlab仿真设计

任务要求：

采用PIC模拟方法的静电模型来模拟多环形电子注在圆柱系统中的运动轨迹。模拟电子枪阴极表面发射电子注，电子在静态场的作用下运动直至稳定的运动过程。其中
系统长：0.01m
系统半径：0.005m
入射的每个宏电子电流：-1/310-3A
初始5圈电子以同样的速度从左平面入射进系统，一圈36个，且横向分布满足轴对称性，轴向速度为107m/s。
最内圈宏电子注入半径：2.010-3m，以后均匀递增2.010-4m
每隔相同的时间间隔从阴极注入一批电子。令Nr=101，Nz=201，Nθ=36,
径向网格数 ，角向网格数 ，轴向网格数 。
总共6000个步长，每个步长为1.6710-13s。
每十个步长注入一次电子。
左平面固定电位为100V，右平面接地电位为0，两平面之间电位均匀变化，且横向等电位。
结果：
根据任务要求的得出仿真结果：

由上面图片不容易看出仿真轨迹，因此对参数进行修改，主要是两个参数，一是宏电子电流大小，二是左右平面所加的电压。其中宏电子电流能够使电子轨迹更加明显，左右平面电压能够使电子在z方向上加速行进。
更改参数为宏电子电流为-1/3*10-3A，左端电压为800V。
得出仿真图：

结果讨论：

分析最初的结果可以看出，由于受电子之间相互作用力的影响，会出现比较小的弧度，但是由于每个宏电子电流大小比较小，所以受到的力不是很明显，当增加宏电子电流大小时，适当增加导体两端电压可以得出非常明显的电子轨迹图。首先可以看出由于电子的初速度及施加的初始电场的作用下，电子沿着z方向移动；然后由于电子之间的相互之间的斥力，我们进行电荷平均分配，经过一定的计算之后得出开始电子会向着导体外方向移动，随着移动距离的增加，电子逐渐分散，外方向的电子密度更大，就会导致电子再次汇聚，最终形成以上的轨迹图。
需要补充的是在程序中没有进行边界的收敛，即边界对电子的吸收处理，因此在两端电压较小并且宏电子电流过大时会导致少量电子在运行一段时间后朝着左平面行进，出现错误。如果两端电压过大，会使电子接触到右平面，出现错误。这个也是可以通过对电压进行合理设置进行避免的，当然也可以进行边界的收敛程序设置，由于任务没有要求，所以仅仅得出大致轨迹。

收获：

在本次设计中，由于对MATLAB软件学习不够，一些能够比较简单的处理在写的时候做的过于复杂，比如构造一个三维矩阵，可以直接用zeros（x,y,z）的方式来建立，但是我却用了一个三重循环的方式，这大大增加了编程的难度。因此我需要在这些的基础上来深入挖掘常见语句的运用，为以后的使用打下基础。
另一方面，我在各个模块化处理方面可以完成任务，但在整个系统的整体性方面做得不够好，一个是没有正确的认识到dt和1/2*dt的关系仅仅是一个时间段或者一个循环内的先后关系，而不是分割成两个循环，再一个是在模块化处理时没有合理的考虑到系统性处理的要求，这使我在最后整合工作时无从下手，感到困难。因此在以后设计中我会牢记模块化与系统性的内在联系性。
最后，在通过将程序与课程的结合，我对粒子模拟这门课程有了初步的了解，我也会在这方面继续前进。
源码下载
https://download.csdn.net/download/Artisan_w/88468540
课程设计文章下载
https://download.csdn.net/download/artisan_w/11156092