一类用于光滑系数的轴对称静电问题的半粗化多重网格方法.PDF

J ou rna l of A nshan N orm a l Un ivers ity 201012, 12( 6) : 33- 37 一类用于光滑系数的轴对称静电问题的 半粗化多重网格方法 王虹宇 (鞍山师范学院物理系, 辽宁鞍山 114007) : 提出了一种用于求解带光滑介电函数的 对称 Poisson问题的半粗化多重网格方法.利用z 方向的 网格半粗化和 r 方向的线松弛磨光技术, 此方法可以克服多重网格方法用于柱坐标系时的慢收敛问题. 算 例表明此方法具有很高的效率, 并可用于等离子体的隐式 PIC模拟; 通过和并行三对角求解器联用, 此方 法还可以改造为并行形式. :介电 Poisson方程; 多重网格方法; PIC模拟;柱坐标系 : O441. 1 : A : 2010 在进行等离子体的模拟时, 对于大尺寸或者高密度等离子体, 经常需要使用隐式 PIC方法[ 1~ 4]. 种隐式模拟方法的核心问题之一, 是求解带有介电函数的 Poisson方程 Ñ [ ( x ) Ñ ] = - ( x ) 个方程的求解也出现在定常的热传导、流体力学中的不可压流动等问题中. 利用有限差分法离散上述方程就成为标准的稀疏矩阵求解问题. 由于存在介电函数 ( x ) , 通常在 [ 5] [ 6] 规则域上常用的快速 Poisson求解器, 如奇偶循环消去法 、快速傅里叶变换 等方法都不能直接应 用, 通常只能用迭代法求解. 一般情况下, 个方程的差分矩阵是对称正定的, 可以使用各种经典迭代法, 如共轭梯度( CG), 广 义最小残差(GMR S)等方法.但在等离子体仿真或者类似的研究中, 上述方程需要在每一个时间推进 7 步中求解, 而总的推进步可以达到 10 以上. 为了降低计算资源开销, 求解器必须高效. 此外, Poisson方 程的条件数随问题尺寸的增加而增加, 因此经典迭代法在应用到大尺寸问题的时候常常呈现慢收敛. [ 7] 处理 Poisson问题的最有效迭代方法是多重网格方法 . 利用多重网格映射和逐级磨光, 对于标准 的 Poisson问题, 多重网格方法可以达到常数收敛率, 因此常常被作为最重要的 Poisson求解器. 在隐式 等离子体 PIC模拟中, 介电函数 ( x ) 来自等离子体密度的累加, 基本上是一个光滑函数, 因此可以用 多重网格方法有效地处理. 在直角坐标下, 已经有很多关于用多重网格方法处理 Poisson 问题的研 究[ 8, 9] . 尽管如此, 在工程上, 很多等离子体设备具有轴对称性. 例如, 在微电子行业中, 经常使用的容性耦 合放电腔室, 其基本构型为一个圆柱状空腔. 种体系的 2维模拟需要使用柱坐标系, 而从文献和实践 [ 10] 中都发现, 正规的多重网格方法在处理柱坐标系的时候存在明显的慢收敛现象 .例如, 在标准的 2 维直角坐标系中求解常系数 Poisson方程, 使用RBGS方法进