主权项 |
一种预测腔体滤波器微放电阈值的方法,其特征在于,包括下列步骤:a.设置输入腔体滤波器的初始电磁信号及其频率f<sub>0</sub>、仿真时间T、时间步长t<sub>0</sub>、总时间步N<sub>t</sub>、输入功率P<sub>0</sub>及其功率扫描范围P<sub>1</sub>~P<sub>2</sub>,P<sub>1</sub>≤P<sub>0</sub>≤P<sub>2</sub>;b.对腔体滤波器进行三维几何建模,设置腔体滤波器的金属材料与金属材料的二次电子发射理论模型,在腔体滤波器中空部分建立粒子模拟区域,将三维几何模型与粒子模拟区域剖分成多个六面体网格,粒子随机分布在粒子模拟区域;c.按时间步迭代进行微放电数值模拟,达到仿真时间T后结束数值模拟,然后转入步骤d;在每一个时间步n,1≤n≤N<sub>t</sub>,进行的微放电数值模拟包括:(i)求解时域有限差分离散形式的麦克斯韦方程得到每一个网格节点和网格线处的电场值和磁场值,求解时域差分离散形式的牛顿‑洛伦兹方程得到每一个粒子的位移与动量;(ii)根据粒子的位移判断粒子是否与腔体滤波器碰撞,若发生碰撞,根据碰撞时粒子的动量与金属材料的二次电子发射理论模型得到出射的二次粒子产额和二次粒子的能量与相位分布,累计所有发生碰撞的粒子数M<sub>i</sub>(n),所有出射的二次粒子的总数M<sub>s</sub>(n);d.求解平均二次电子发射M<sub>a</sub>(n),其中<img file="FDA0000560974210000011.GIF" wi="359" he="140" />1≤n≤N<sub>t</sub>,对平均二次电子发射进行内插均值处理,具体方法如下:设置内插时间步N<sub>1</sub>,2≤N<sub>1</sub><N<sub>t</sub>,从M<sub>a</sub>(n)开始,依次取其及以后N<sub>1</sub>‑1个时间步的平均值作为内插均值二次电子发射M<sub>a1</sub>(n),其中<img file="FDA0000560974210000012.GIF" wi="538" he="216" />1≤n≤(N<sub>t</sub>‑N<sub>1</sub>+1);e.根据初始电磁信号的频率f<sub>0</sub>设置粒子雪崩效应判定起始点N<sub>m</sub>,设置微放电阈值判断基数M<sub>2</sub>,M<sub>2</sub>≥1,从粒子雪崩效应判定起始点N<sub>m</sub>开始对M<sub>a1</sub>(n)与M<sub>2</sub>进行比较,若在N<sub>m</sub>≤n≤(N<sub>t</sub>‑N<sub>1</sub>+1)的范围内,M<sub>a1</sub>(n)小于M<sub>2</sub>,输入功率为P<sub>0</sub>时不发生微放电,若在N<sub>m</sub>≤n≤(N<sub>t</sub>‑N<sub>1</sub>+1)的范围内,M<sub>a1</sub>(n)大于M<sub>2</sub>,输入功率为P<sub>0</sub>时发生微放电;f.令输入功率P<sub>0</sub>在设置的功率扫描范围内由P<sub>1</sub>开始以第一间隔步进,分别进行步骤c至步骤e,若扫描至输入功率为P<sub>3</sub>时不发生微放电,P<sub>1</sub>≤P<sub>3</sub>≤P<sub>2</sub>,扫描至P<sub>4</sub>时发生微放电,P<sub>3</sub><P<sub>4</sub>≤P<sub>2</sub>,从P<sub>3</sub>开始,以第二间隔步进,分别进行步骤c至步骤e,若输入功率为P<sub>5</sub>时不发生微放电,为P<sub>6</sub>时发生微放电,则预测P<sub>5</sub>为腔体滤波器的微放电阈值,若输入功率步进至P<sub>2</sub>时仍未发生微放电,则预测在功率扫描范围P<sub>1</sub>~P<sub>2</sub>内不发生微放电。 |