边界元素法及在NSRL束测调试中的应用.pdfVIP

边界元素法及其在NSRL束测系统调试中的应用 刘建宏． 王筠华． 刘祖平， 孙葆根． 裴元吉 中国科学技术大学国家同步辐射实验室 合肥230029 [摘要] 我们用它计算了BPM探头的灵敏度，并以此检验模拟定标测试数据的可信程度。作为该方法 的进一步探讨．我们尝试．由模拟程序直接计算出有机械加工误差时的灵敏度，给出了机械 加工要求的一个参考。 关键字：边界元素法 柬流位置监测器模拟定标 一．前言 储存环束流位置监测器(BPM)，是用以确定环中柬流沿环形轨道，诒：横截面方向上位簧 的束流探测器。对于NSRL电子存储环图型真空室，倾斜45”放置的圆型钮扪电极(见圈3)． 计算定点的束流位置与四个感应电压的关系： 燃×糍y咄×黼㈩ ^，，k．与探头的几何形状有关．称之为探头标定系数u ，： JJ (2) 1+I，+ y一一 K一屹 型“黼 令。=寺=岳舞1}、 I’：+I = (3) 则灵敏度被定义为：S，：掣 ， 坐秒 二．边界元素法 2．1理论简介 在光滑束流管道中，相对论带电粒子束流在横向j因为“洛仑兹”效应而受到压缩，橙 截面很小。在稳定的-隋况下，它的电场相当于无限长线电荷的静电场。当然，束流也有磁场 产生．但钮扣电极对其不敏感，所以只需考虑电场部分【2J。 荷来等效地替代它，如罔Figl(b)所示。 一， t+ ．—i／p：．7-9 ，p 十 ／，1…c譬≯ ，／ ++一≥譬 + (二二，一≥’擘。 十：+十 (b)空间自由电荷描述的等价问题 Figl(a)场满足的方程wV二巾=0和边肄条件 用解析法直接解这样的微分方程通常是不可能的，利用BEM法将其转化成直接解积分 一177一 方程的近似计算方法使问题迎刃而解。引用流行的加权余量法【3l的结论得到，在边界内任何 一点场分布如下式表达： (4) (p(r1=f。二P【，·。)I【(r．r。)dr．，+』。d(s)∥(，，j)c缸 情况可简单取为格林函数 G： ! 】。上 (5) 2Ⅱ￡c。 r 由于要利用数值积分手段，边界被划分成 2所示。 若干直线段(边界元素)．1埘图Fig 此处面和(3-分别为第i个元素的场势和感 应面电菏密度， 2边界划分示意图 现在，可以由方程(4)与(5j推出第i个元素 Fig 凼}㈣ 2Ⅱsn 的电势腻龇。，=j鼍m去+釉掣k丽1r2∞。 {)篇‘o l，(s) E¨电荷束流的线电茼密度：r¨从束流到第i个元素中心的距离：o，(s)：第j个元素上的 感应电荷密度；s：第j个元素上的长度积分参数；r．，(s)：从第i个元素的中心到第j个元素 上的各积分点的距离 式(6)中．第一项代表柬流‘。在面元之处产生的势，第二项是所有边界元素上的感应电 荷对面元之处的贡献，当边界元素划分足够小时．O．町近似看作在元素j上是常量。将边界 分为N个元素．由方程(6)得到