强流电子枪数值分析的虚边界元方法研究.pdfVIP

殷海荣宫玉彬王文祥魏彦玉 (电子科技大学物理电子学院，成都610054) [摘要]以虚边界元方法为基础，利用FORTRAN语言开发了一个计算强流电子枪的程序， 对强流电子．睑进行了计算。计算结果表明，选取适当的虚实边界最小距离，虚边界元解法 比同等条件下的实边界元法求解精度要高。适当的虚实边界最小距离是获得最好计算精度 的必要条件，最好精度还与边界单元长度值有关，只有在最佳距离与最佳边界单元长度达 到最佳组合时，才能获得最好的求解精度。虚实界元法的计算速度比同等条件下的边界元 法略快。 [关键词]强流电子枪虚边界元方法边界元方法数值计算 强流电子枪的数值分析中，早期常用的基本方法是有限差分法和有限元法。这两种方 法需要对整个求解区域进行计算，计算速度慢，对于电子枪及其漂移区的电磁场分析来说， 这是不必要的。而且，这两种方法无法直接求解无限域，这个缺陷无可避免地会给强流电 子枪的分析带来误差。在进行强流电子枪的初步设计时，需要进行大量重复的计算，基于 有限元法和有限差分法的程序由于其速度关系无法满足这种大量计算的需要，寻求—种耗 时少且求解精度高的方法将会推动强流电子枪的设计和新结构的发展。 边界元法只需对电极边界进行数值分1析，再利用其结果解析求出域内的电位，它不需 要对整个区域作全面的计算，也不要求区域必须封闭。从边界元法的这些好处来看，它十 分适宜于用在强流电子光学系统的分析之中。但是，边界元法存在着大量的奇异积分，如 何处理这些奇异积分成为它解的精度的关键因素，也是它的运用受到限制的最大原因。 虚边界元法[1]可以完全避免边界分析．亡的奇异积分。它自1984年被提出以来，在力 学领域内一直受到重视和研究，但在电磁场溅尚未有应用的报道。虚边界元法具有边 1爱釜灏目：孝幻争溯丛纭如|c秀人才j￡持项目资助，重点婆习筮童基金，514400∞104D囱0囱砣。 353 界元法的所有优点，然而在实际计算中，虚实边界最小距离的确定却存在着实质性的困难， 这就使得解的稳定性无法得到保证。本文推导了求解电子枪电磁场的虚边界元法的边界离 散积分方程，利用所编写的轨迹模拟程序对强流电子枪进行了计算。 1、虚边界=耐蝴憾钥瓣 虚边界元法的楱心思想是将场延拓到无限域，并在实边界之外虚设—场域边界，假设 虚边界匕分布有场源函数Q(r)，对于位势问题，若待求问题满足第一=类边界条件的泊 松方程，可以将虚边界的域内任何一点的位势值及其法向导数摹沥专为： (1) C n,rw)6(，一)dQ(，一) ￡，(r，)；Jr，￡，‘(r．r，)Q r)dr，(，)+J。￡，·(r 鼍竽：J，．等∽，，)Q㈩∥㈩+J口簪∥／川一dQ∥， ‘2’ 上式中，U‘表示满足控银防．程位势问题的基本解，r7表示场点，，．和，．’分别表示虚 边界匕的虚源点和实场域内的源点，r1表示虚边界，Q表示实边界所围自勺=场域，n表示 沿谋卜．方向的单檀乡遏，b是区域内的源函数。 2、用虚逝孵磁持峥阿虢电荆曲嘱睛嗡匿 在金属电极边界条僻下，由于不存在偶电层的作用，由式(1)，取基本解为1／R，因 此有： 帅，=亡上最静：∥，+亡』器an／ Q) U(r)Ih岫=，(r) 将虚边界单元匕的虚设电觥由两部分眵i：组成，仃肋由实边界匕的电位感应．所产 生，仃p由空间电荷2窗鳓听产生，令： (4) 仃=仃朋+仃p 这样，将电位U也看成由两部分揪： c，．：亡上}：’ u‘h岫2，