CST自带案例:粒子枪仿真和Track Solver追踪求解
这一期我们一起看一下CST自带案例之一,电子枪。ComponentLibrary 中搜Gun,然后选ElectronGun模型,其他粒子枪模型也是不错的学习案例。直接看模型不容易理解,这里我们从模板开始,重新建模仿真学习粒子工作室的Trk-solver追踪求解器。
Step 1. 粒子追踪模板
粒子枪可以属于加速器件类,也可以属于真空电子装置类,总之Trk-solver是推荐,辅助为静电场求解器Es-solver和静磁场求解器Ms-solver。
Step 2. 阴极阳极建模
先画第一个圆柱作为离子源和阴极,圆柱取名为Emission,以坐标原点为中心,半径1mm,Z+方向高为0.1mm,材料为默认PEC,在粒子工作室中的PEC默认不会让粒子通过,也不会发生电子撞击造成的电离,可在属性中粒子界面根据需要更改,这里我们用默认就可以。
然后利用本地坐标,在Z=20mm的平面上画第二个空心圆柱,内半径5mm,外半径6mm,高1mm,材料同为PEC,取名Anode,作为阳极。两个圆柱中心为同轴Z轴。
然后定义两个电势差,圆柱Emission为potential1 = -5000V,蓝色阴极;圆环Anode为potential2=0V,红色阳极。
然后定义粒子源,选Particle Circular Source,再选阴极圆柱的上沿,这里我们假设电子以同心圆的形式产生。默认Emissiondensity 发散密度中的Lines为5,就是5个同心圆区域,points为81,就是图中蓝点个数81。
粒子源定义追踪算法,选Space Charge, 编辑中发射电势选potential1,参考电势自动是potential2,就是-5000V,下面的粒子属性选electron,也就是我们的模型是个电子枪。如下图:
CST还有很多追踪算法,这里不过多解释。简单地说,Space Charge是空间粒子模型,粒子形成的电流密度满足Child-Langmuir模型;空间中的电子会将理想线性电动势拉低成并变成非线性,甚至将发射电极的位置等效地向参考电极位置靠近,如图中的蓝色点面和红网格就不在一个平面上。
这里还可以定义二维平面的粒子监视器,可同时定义多个平面,这里我们定义一个在Z=30的平面上。当然在电子束路径上每隔一段距离加一个平面对后处理有好处,比如可以用来提取包络形状,或者像医院CT扫描一样切片看粒子分布。
边界条件为open (add space if), 并且我们手动在Z+方向上多加一些距离,因为电子路径会朝着Z+方向延伸。
Step 3. 仿真和结果
接下来我们直接开始仿真,看看这样的阴极阳极能有什么效果。Trk-solver中选Low Frequency和E-static静电场,因为我们的阴阳极都是静态,这里会自动调用Es-solver。
仿真结束后,三维结果查看Trajectories,粒子路径。
可以看到粒子颜色由蓝变红,能量增加,被圆环阳极吸引发散。
再看静电场E-field,也就是电势对距离的偏导,大部分集中在阴极附近,库仑力加速基本也就集中在靠近阴极的区域。
还有粒子密度,由于是电子,颜色表示相反,蓝色的库伦密度是更密集电子。
因为目前只有电场加速没有磁场改变方向来聚集波束,目前还不算是个电子枪,最多算个简单的电子加速器,还是发散的。
Step 4. 加线圈静磁场仿真
接下来我们在阳极圆环前方加个同轴线圈,电流为0.4A,1000匝,type为stranded,也就是电流均匀分布于截面。建模可以用sweep curve(画个圆圈曲线,方型截面curve,sweep curve功能)。
有了线圈,Trk-solver 中就可以加上静磁场了,直接仿真。
仿真结束我们先看B-field磁感应强度,不出所料的回路分布。这里可以用洛伦兹力预测一下电子的受力和路径。
比如拿阴极圆柱上半部为例,也就是图中Y-方向,红色箭头指尖位置。生成电子右移(Z+),电流向左(Z-),磁场有少量向下分量(主要分量在Z+),根据左手定则,受力由纸面向外(读者应该是手心朝上的“无辜”手势,拇指向自己),所以从右看是顺时针旋转地向线圈和阳极移动。
好了不为难个别左右手定则搞不清的朋友,直接看路径结果:
从阴极电子源看进去顺时针旋转。确实电子束有收敛,但是聚焦点很远,下一步我们要再加强磁场。
Step 5. 再加线圈静磁场仿真
将线圈coil1 复制一个在阳极圆环后面。开始仿真,仿真结束后可看到电子束焦点。
我们的二维监视器平面也是接收到了81个电子。
其他后处理包括提取不同截面的粒子分布,焦点的坐标,能量直方图等等,这里我们就留给读者自己探索。
最后再提个小技巧,可以在求解器列表选转换Es-solver 或Ms-solver做一些设置,再回来Trk-solver继续仿真。比如去静电场Es-solver的specials 设置不同的边界的电势处理。