二维粒子输运问题GS及SOR迭代方法.docxVIP

二维粒子输运问题G—S及SOR迭代方法 摘要：解决粒子输运问题最基本方法是源迭代。对于简 单系统，源迭代收敛速度很快。然而对于复杂系统用源迭代 求解，收敛速度变得很慢。文章考虑采用高斯-塞德尔(G-S) 迭代方法和超松弛(SOR)迭代方法求解粒子输运问题，加 快迭代收敛速度。数值实验表明，G-S和SOR迭代方法收敛 所需的次数和CPU时间比源迭代要少。 Abstract: The most basic transport iteration method is source iteration. For the simple physical system, source iteration converges fast? But when the physical system is complicated, it is inefficient. In the paper, we shall use Gauss-Seidel (G-S) and Successive Over Relaxation (SOR) iterative methods to solve the particle transport problem to accelerate the convergence speed. Numerical experiments show iterates and CPU time needed for convergence are less than those of the source iteration. 关键词：粒子输运；源迭代；超松弛迭代 Key words： particle transpor t； source it era tion； SOR iterative method 中图分类号：0151.2文献标识码：A文章编号： 1006-4311 （2013） 03-0299-02 0引言 最基本的输运迭代方案为源迭代，对于一个粒子经历很 少碰撞的问题，源迭代收敛得很快，此类问题要么是泄露系 统（粒子可能经历几次碰撞就很快逃离外边界），或者处于 一个系统里到处存在强吸收（粒子可能经过几次碰撞就被俘 获）。然而对于扩散特性的空间区域，其中包含光学厚（泄 露几率很小）及散射居主导地位的复杂系统（俘获几率很 小），大量粒子被俘获或泄露出去之前需要经历多次碰撞， 仍应用源迭代求解，收敛速度会很慢，成本昂贵。而很多实 际问题都属于后一类型，因此，人们一直渴望加快SI的收 敛速度。 1源迭代 柱坐标（r, z； 3, n ）下的输运方程为： ■ ■ （r ?准（r, z； 3, □）） + T] ■ （?准（r, z； 3, n ）） -■ ■ （ 4 ?准（r, z； 3, n ）） +Z ?准（[,z； 3, n ） =S （r, z） （1） 其中，S （r, z） =q （r, z） +■ ?蓦？准（r, z； o , T] ） d T] do U =■ ? U=?cos3, r| =■ ? -1 ?燮 r| ? 燮 1,0?燮 s?燮 ji , r （0, R）, 4 -■ ? 2二■sin 3。 上述方程中各物理量的含义为：工是总截面，o■是散 射截面，q是内部的各向同性源，S是源，？准是待求的角 通量。 为了方便，在简单连续性方程条件下引进一些基本概 念，我们定义如下简化抽象的算子形式： L?准？縮■■?准）+n ■?准-■■（ ￡ ?准）+工？ 准 S?准+q?励 S （r, z） 则方程（1）可以写为L?准二S?准+q。 最基本的输运迭代方案为源迭代，其定义如下： L?准■二S?准?+q 1?叟0 其中？准■由使用者给定，每次迭代的开始，都需把一 个前一时刻的？准■代入上式的右端，得到一个新的估计 值，迭代过程重复执行，直到连续两次估计值小于预先指定 的数值。 2 G-S迭代与SOR迭代 源迭代是解决粒子输运问题的最基本方案，它具有简单 的物理解释。对于大量粒子在被俘获或泄露出去之前只需经 历很少碰撞的简单系统，源迭代很有效。而对于包含有光学 厚（泄露几率很小）及散射居主导地位的复杂系统（俘获几 率很小），大量粒子被俘获或泄露出去之前需要经历多次碰 撞。在这些系统中，应用源迭代收敛很慢，成本昂贵。实际 问题中很多问题属于后一类型，因此长久以来，人们一直渴 望加快源迭代的收敛速度。 从数值分析知道，简单迭代法收敛速度是很慢的，为了 提高收敛速度，自然想到把已算出的值代替旧值。因此，我 们采用Gauss-Seidel迭代法和SOR迭代法对离散纵标方程 求解。从求解格式可以看出，G-S和SOR方法与源迭代相比, 所需的计算量和存储量都是一样的，没有任何明显的增加， 但数值实验表明，迭代次数