基于BEM泊松方程求解的空间电荷效应数值模拟.pdfVIP

第 20卷第 1O期 强 激 光 与 粒 子 束 Vo1．20，No．10 2008年 1O月 HIGH POW ER LASER AND PARTICLE BEAMS Oct．，2008 文章编号 ： 1001—4322(2008)10—1719—05 基于BEM 泊松方程求解的空间电荷效应数值模拟 冯义章， 秦 斌， 樊明武 (华 中科技大学 电气与电子工程学 院，武汉 430074) 摘 要 ： 为了模拟强束流在加速器及其传输线 中的行为，用 c++语言开发了一种包含空间电荷效应的 多粒子跟踪程序(PTP—SC)，它在经典的 PIC方法基础上。基于边界元法 (BEM)和非等距 的网格求解泊松方 程。束流在 自由空间分布的仿真结果与解析结果保持较好的一致性 。给出了一条注人线的模拟计算结果，并 与ORBIT，TRACE3-D的计算结果进行比对。结果表明：该程序与采用数值方法的ORBIT程序的计算结果 有较好 的一致性 。该程序可用于直线加速器及 回旋加速器中的空间电荷效应模拟 。 关键词 ： 空间电荷效应 ； 边界元法 ； 网格 ； 注入线 ； 泊松方程 中圈分类号 ： TL501．5 文献标志码 ： A 质点模拟技术是通过求解一系列质点的运动方程来模拟多体系统 ，其主要应用领域包括高能物理学、流体 动力学和天体物理学 。质点跟踪方法使用了没有太多近似的基本方程 ，这使得模拟能够保持系统的物理特性。 PIC方法是该领域应用较多的方法之一 ，其主要思想是用质点与网格之间的相互作用来代替质点与质点之 间 的相互作用。通过添加数字网格使我们能够更加有效地计算作用在模型质点上的空间电荷力。算法的基本步 骤描述如下：(1)产生初始柬流分布 ，(2)宏粒子轨迹运动方程积分 ，(3)电荷在网格上重新分配，(4)求解 内 部空间电荷场 ，在每一步的计算中重复步骤 (2)至 (4)l1J31。 传统的求解积分是采用简单的二阶蛙跳 ，它能够在保持一定精度的前提下获得较高的效率 ，本文的方法中 采用的是经典的龙格一库塔法 ，并且粒子电荷不再分配在节点之上 ，而是分配在每个网格之上，对粒子所 占据 的 区域进行直接积分来求解空间电荷场[4]，进而计算 出作用在这些宏粒子上的电荷力，运用精确的数值积分算 法，计算结果的精度得到显著提高。 1 自由空间的BEM 计算表达式 首先直接给出势问题泊松方程 “+P一0的基本解 ] 1／27r)ln(1／r．~。) (2一D ) U(X，x。)一I 付 、 (3一D) r(n／2) 1 ( D ) f 。“+ P 一 0 (inn) {U一一U一 ((OilSS1)) ((22)) lau／an— q (onS2) (x)一(，)一一J．fn“pd+-JfS“0U口*ds—JJlSuqds (3) *收稿 日期 ：2008—03—25； 修订 日期：2008—08—23 基金硬 目：国家 自然科学基金项 目 作者简介：冯义章(1982一)，男，硕士，从事加速器结构及束流动力学研究；yzfeng@yahoo．cn。 ● 强 激 光 与 粒 子