当前位置：文档之家 › 中子输运方程SP2方法的并行计算1

中子输运方程SP2方法的并行计算1

  • 格式：pdf
  • 大小：132.35 KB
  • 文档页数：4

下载文档原格式

  / 4

合集下载

中子输运方程特征线解法及嵌入式组件均匀化方法的研究

中子输运方程特征线解法及嵌入式组件均匀化方法的研究

中子输运方程特征线解法及嵌入式组件均匀化方法的研究目前一代轻水堆堆芯物理分析理论和数值计算方法形成于上个世纪七十年代末。

受限于当时计算机的硬件条件,理论上存在近似和不足的现行堆芯物理分析方法已难以满足当今日益复杂的堆芯设计需求。

因此,近些年国际上广泛地开展了下一代(第三代)轻水堆堆芯数值计算理论和方法(Next Generation Methods, NGM)的研究。

归纳起来,现有NGM研究的主要思路是完全抛弃组件均匀化和堆芯扩散计算而直接进行Pin-by-Pin的三维全堆输运计算,这样势必会造成对多年积累的堆芯物理分析方法、经验、程序等资源的浪费。

为此,本文提出了一套全新的NGM方案,目标是将组件输运计算产生等效均匀化参数的过程嵌入三维扩散计算中。

本文通过改进输运计算方法、组件均匀化方法和等效均匀化参数产生方法,满足未来复杂堆芯物理计算精度与速度的要求。

本文的主要工作集中于中子输运方程特征线解法的研究、程序PEACH的研制以及嵌入式组件均匀化方法的可行性分析。

首先,本文基于特征线法研制了NGM的输运计算工具PEACH程序。

程序PEACH的几何处理,采用了以组件为模块产生特征线的方法(Assembly Modular Ray Tracing, AMRT)。

这种特征线的产生方法既具备计算速度快和节省存储量的优点,又可以比较充分地考虑组件内子区的任意性。

程序PEACH的特征线跟踪计算方法,分别采用了平源近似的步特征线法、菱形差分的特征线法和线性源近似的特征线法。

本文在探讨后两种数值计算模型时,提出了相关负中子角通量和负中子源的处理方法,并通过OECD/NEA C5G7-MOX 2D 等基准题的数值检验,证明线性源近似的特征线法在相同精度的前提下占用更少的系统内存和计算时间。

程序PEACH中运用了由多群粗网有限差分(Coarse Mesh Finite Difference, CMFD)与少群粗网有限差分组合而成的两重粗网有限差分加速方法,并建立了指数函数插值表,呈现出良好的加速效果。

以栅元为模块进行特征线跟踪的中子输运方程解法

以栅元为模块进行特征线跟踪的中子输运方程解法

以栅元为模块进行特征线跟踪的中子输运方程解法
汤春桃;张少泓
【期刊名称】《核动力工程》
【年(卷),期】2009(0)4
【摘要】为解决复杂几何条件下中子输运方程的求解问题,分析了特征线法理论模型,探讨了以栅元为模块的高效特征线产生方法,以及与之相关的空间角度离散和边界条件处理问题。

采用自行研制的特征线法数值计算软件——PEACH,对经济合作组织核能机构(OECD/NEA)UO2和MOX燃料混合装载的7群(C5G7MOX)基准问题的数值进行了检验。

结果表明,无论是计算keff还是棒功率分布该方法都具有很高的精度。

【总页数】5页(P32-36)
【关键词】特征线法;跟踪;中子输运方程
【作者】汤春桃;张少泓
【作者单位】上海交通大学核科学与工程学院
【正文语种】中文
【中图分类】TL329.2
【相关文献】
1.采用栅元重组CMFD加速技术的六边形几何模块化特征线方法及程序验证 [J], 韩宇;蒋校丰
2.基于非均匀射线追踪的栅元模块化特征线方法研究 [J], 彭良辉;陈笑松;陈耀东;
刚直
3.基于 OpenMP 的中子输运方程特征线法并行计算研究 [J], 于锐;赵强
4.线性源近似的中子输运方程特征线解法 [J], 汤春桃
5.基于AutoCAD二次开发实现中子输运方程特征线法求解 [J], 陈其昌;吴宏春;曹良志
因版权原因，仅展示原文概要，查看原文内容请购买。

第三章 中子扩散理论

第三章 中子扩散理论
J z (r ) J z (r ) J z (r )
s (r )
3 z
J z (r ) 叫做z方向的中子流密度或净中子流密度，若dA的取向
与x轴垂直，沿着x方向穿过dA平面单位面积净中子数Jx为
J x (r )
s ( r )
3 x
同样，沿着y方向穿过dA平面单位面积净中子数Jy为
推导中并没有考虑中子源的贡献，中子流密度的贡献只 是来自中子与介质核的散射碰撞 在强中子源两三个平均自由程的区域内，斐克定律不 适用。
3.2 非增殖介质内中子扩散方程的解
稳态单能的中子扩散方程
D2 (r ) a (r ) S (r ) 0
无源情况下，即除中子源所在的位置以外的无源区域,扩散 方程有以下形式: 1 D (r ) 2 2 2 2 L ( r ) 0 (r ) (r ) 0 或 k2 L2
第三章 中子扩散理论
中子在介质中的输运过程中 的运动状态由位置矢量r(x,y,z), 能量 E, 和运动方向Ω表示。 Ω通过极角θ 和方位角φ 来 表示
dS r 2 sin dd d 2 sin dd 2 r r
中子角密度函数n(r,E, Ω)定义： 方向 Ω的表示 在r处单位体积内和能量为E的 单位能量间隔内，运动方向为 Ω的单位立体角内的中子数目。 中子角通量密度定义为: ( r, E, ) n( r , E, )v( E ) 对中子角密度和中子角通量密度积分便可得到与运动方向无 关的标量中子密度和标量中子通量密度
4

tr d
6 dx
0
d dx
x 0
30 2tr
应用输运理论和扩散理论的 外推距离求得的扩散方程的解

基于特征线方法的三维中子输运程序(Ⅰ)——边界条件的插值处理

基于特征线方法的三维中子输运程序(Ⅰ)——边界条件的插值处理
的方 法来 解 决 空 间 问题 和角 度 问题 ，采 用 这 种 方
收 稿 日期 ：２ ０ — １ ９ ０ －７
而且对方程 （ ） ２ 积分可以得到射线线段 ｋ 上的平 均角 通量 ：
ｌ ２
核 动 力 工 程
Ｖ Ｉ３ ． ． ． ０ ０ ｂ． １Ｎｏ２ ２ ｌ
十 （） （） （） ：Ｑ 式中 ， （） １
特征线方法在处理三维 中子输运问题 的反射 边界 条件 时存 在 空 间 问题 和角 度 问题 ，然 而 过 去 的研究主要集中于二维问题 的边界条件处理 ， ： 如 ＲｃａｄＳｎｈｚ等人 系 统地研 究 了精 确处 理 二 维 ｉ ｒ ａ ｃｅ ｈ 问题边界条件的方法 【，但是精确处理边界条件 ６ 】 的方 法对 求 积 组选 取 及 射线 布 置均 有 限制 ，这 种
摘要 ：特征线方法 在处理三维 中子输运 问题 的反射边界条 件时存在空间问题和角度 问题 ，本文提出 了采 用 平面插值和球面插值 方法来分别处理这两个 问题 ，插值方 法不 仅可以有效解决空 间问题 和角度问题 ，而且 三维特 征线 方法在采用插 值方法处理边界条件后对 离散求积组和射线 布置 没有 要求 。
ｏ
法 后 ，三 维特 征 线方 法 对 离散 求 积组 和射 线 布置 没 有要 求 。
２ 特 征 线 方 法简 介
特 征线方法是在 空间离散角 的直线上求解 中子输运方程 ， 这些直线可以认为是中子 的飞行 路线 ，在这些直线上 ，中子输运方程变为全微分 的形式 。假设几何 区域被分成小网格 ，而且每一 个 网格 只含 一 种材 料 ，网格 内部通 量 和源 是不 随 空 间变 化 的 ， 也就 是说 使 用平 源 近似 。 在方 向 上 ，输 运 方程 为 ：

蒙特卡洛方法在中子输运中的应用

蒙特卡洛方法在中子输运中的应用

《中子输运理论与数值方法》课程作业——蒙特卡洛方法目录1. 前言 (3)2. 蒙特卡洛方法概述 (3)2.1 蒙特卡洛方法的基本思想 (4)2.2 蒙特卡洛方法的收敛性、误差 (4)2.2.1 蒙特卡洛方法的收敛性 (4)2.2.2 蒙特卡洛方法的误差 (5)2.3 蒙特卡洛方法的特点 (6)2.4 蒙特卡洛方法的主要应用范围 (7)3. 随机数 (7)3.1 线性乘同余方法 (9)3.2 伪随机数序列的均匀性和独立性 (9)3.2.1 伪随机数的均匀性 (9)3.2.2 伪随机数的独立性 (10)4. 蒙特卡洛方法在粒子输运上的应用 (10)4.1 屏蔽问题模型 (10)4.2 直接模拟方法 (11)4.2.1 状态参数与状态序列 (11)4.2.2 模拟运动过程 (12)4.2.3 记录结果 (15)4.3 蒙特卡洛方法的效率 (16)5. 蒙特卡洛方法应用程序—MCNP (17)5.1 MCNP简述 (17)5.2 MCNP误差的估计 (18)5.3 MCNP效率因素 (19)6. 结论 (19)参考文献 (20)1.前言半个多世纪以来，由于科学技术的发展和电子计算机的发明，蒙特卡洛(Monte Carlo)方法作为一种独立的方法被提出来，并首先在核武器的试验与研制中得到了应用。

蒙特卡洛方法是一种计算方法，但与一般数值计算方法有很大区别。

它是以概率统计理论为基础的一种方法。

由于蒙特卡洛方法能够比较逼真地描述事物的特点及物理实验过程，解决一些数值方法难以解决的问题，因而该方法的应用领域日趋广泛。

蒙特卡洛模拟计算是解决中子在介质中输运较为成熟、有效的方法，对于原子能、辐射防护、剂量学和辐射生物物理学等研究领域实际问题的计算，都可以利用蒙特卡洛方法予以实现。

粒子输运过程可以用玻耳兹曼方程加以描述，然而，以此基础上发展起来的近似数值方法如扩散近似法、离散坐标方法在处理截面与能量相关以及散射各向异性介质、复杂几何条件问题时碰到了较大困难。

基于特征线方法的三维中子输运程序(Ⅱ)——数值验证

基于特征线方法的三维中子输运程序(Ⅱ)——数值验证

精度略高 ；相对 于角度一次插值 ，角度二次插值 得到的计算结果 的精度有大幅提高。角度二次插 值能够大幅提高计算精度 的原因是 ：Ｓ或者 ６ ｓ ｌ 离 散求积组中的角度数 目 仍然较少 ；由于二次插值 的代数精度高于一次插值 的代数精度 ，因此角度 二次插值后得到的计算结果 ， 其精度有较大提高。 当选用 ８ 离散求积组时 ， 求积组中的角度数 目较 多 ，采用一次插值即可达到很好 的计算精度 。空 间二次插值计算精度提高较小的原 因是 ：射线之 间的间距较小 ， 一次插值已达到很好 的计算精度 ， 采 用 二次插 值后 精度 提高有 限。
Ｆ ｇ ４ Ｓ ｈ ｍａ ｉ ａ ｒ ｍ ｆ ｇ ｌｒ ｅｒ ｈ ｄ ｏ ｉ ． ｃ ｅ ｔ Ｄｉｇ ａ ｏ Ｒｅ ｕ ａ ｔ ｅ ｒ ｎ ｃ Ｔ ａ
图５ 正 四面体基准题的ｅｒ ｈ ｄ ｏ ｉ ． ｓ ｏ ｆ ｕ ａｉ ｎ ｏ Ｒｅ ｕ ａ Ｔ ｔａ ｅ ｒ ｎ ｉ ｒ
第 ３ 卷 第 ２ 期 １
２Ｏ １ ０
核 动 力 工 程
Ｎｕｃｅ ｒＰｏ ｒＥｎｇ ｎ ｅ ｉ ｇ ｌａ ｗｅ ｉ ｅ ｒｎ
Ｖ． ．３ ．Ｎ０ ２ ０ １ １ ． Ａｐ ．２ ０ ｒ １ ０
年 ４月
文章 编号 ：０ ５ ．９ ６２ １）２０ １．５ ２ ８０２ （０ ００ ．０ ６０
摘要 ：介绍 了开发 的三维 中子输运特征线方法程序 Ｔ Ｍ。使用数值实例对 Ｔ Ｍ 的准确性进行 了验证 ， Ｃ Ｃ
同时比较 了不 同插值 阶数 、 离散角度数 目和 网格大小对计算结果的影响 。 结果表 明 ， 三维特征线方法程序 Ｔ Ｍ Ｃ 能够求解任意几何形状及反射边界条件的中子输运 问题 ，计算结果具有较高的精度 。 关键词 ：输运方程 ；特征线方法 ；反射边 界条 件 ；插值 ；基准题

基于%._d二次开发实现中子输运方程特征线法求解

" / / 1 2 . ) ) 0 ) ; 2 0 # ' *1 2 3 # 4 2 #' 4 !# 2 . 4 ) 0 ) K 3 ' 4L 3 ' ) + ) 4 4 3 ? # * : 3 + K 3 ' 4F < ) ) ! *! B . 3 4 '# ,J 8 %! 8> %!
% : 9 5 . ( + * . : 22 , : $ 8$ Q/ : 7 3 7 / , 2 3 . P , . / P" ;? ># P $ % O 2 P, : 21 2 J , 3 $ 1, 3 7 1 P $ 3 ,2 J 7 , . $ 1 9 R f ! 7 % $ 1 7 / :/ : 7 3 7 / , 2 3 . P , . / P % . 1 2 & 9 : 2 $ 3 2 , . / 7 % % . , . P . 1 8 2 2 1 8 2 1 , $ Q , : 24 2 $ 2 , 3 : 7 2$ Q R 5P R 42 5 ! !7 S $ J 1 8 7 3 1 83 2 . $ 1&_ $ T 2 O 2 3 , : 24 2 $ 2 , 3 : $ J % 8S 28 2 P / 3 . S 2 8 1 83 7 3 7 / . 1 57 4 5P 5, 4 !, P $ J % 8S 2/ 7 3 3 . 2 8$ J ,Q . 3 P , % &H 1 $ P ,$ Q, : 22 N . P , . 1 >/ $ 8 2 P : 24 2 $ 2 , 3 5 4 ;? 5 8 2 P / 3 . , . $ 17 1 83 7 , 3 7 / . 1 , . % % : 7 O 2$ 3 2$ 3 % 2 P P % . . , 7 , . $ 1& 9 $T . 2$ Q Q , : 2 P 24 2 $ 2 , 3 R 5 4P R 5 % . . , 7 , . $ 1 P Q $ 3J P . 1 >. 1$ 3 2/ $ % 2 NQ J 2 % 7 P P 2 S % 7 % / J % 7 , . $ 1! 71 2 T ;? >/ $ 8 2 4 ;? R 5/ + J , $ ;? >S 7 P 2 8$ 1, : 2/ J P , $ . U 7 , . $ 1$ Q+ J , $ > + \ T 7 P8 2 O 2 % $ 2 8 &9 : 21 J 2 3 . / 7 % R 3 2 P J % , PP : $ T, : 7 ,+ J , $ ;? >1 $ ,$ 1 % 7 P: . :Q % 2 N . S . % . , 14 2 $ 2 , 3 J ,7 % P $4 . O 2 P 5: 4 5. 5S 7 / / J 3 7 , 2P $ % J , . $ 1 P . 1/ $ 7 3 . P $ 1T . , :P , 7 , 2 D $ Q D , : 2 D 7 3 ,2 , : $ 8 $ % $ . 2 P & R 4 % & & & & ; ' / ( > 5 2 , : $ 8$ Q/ : 7 3 7 / , 2 3 . P , . / P + J , $ > + \ / J P , $ . U 7 , . $ 1 / $ % 2 N4 2 $ 2 , 3 R 5 <= 1 2 J , 3 $ 1, 3 7 1 P $ 3 ,2 J 7 , . $ 1 R f 在 未 来 的 先 进 反 应 堆 的 组 件 设 计 计 算 中! 可能需精确求解复杂几何情况下的中子输运方

蒙特卡罗中子输运程序的MPI／OpenMP混合并行研究


Ａｂｓ ｔ ｒ ａ ｃ ｔ ： Ｐａ ｒ ａ ｌ ｌ ｅ ｌ ｐｒ ｏ ｇ ｒ ａ ｍｍ ｉ ｎ ｇ ｗｉ ｔ ｈ ｍｉ ｘ ｅ ｄ ｍ ｏ ｄ ｅ ｏ ｆ ｍ ｅ ｓ ｓ ａ ｇ ｅ ｓ — ｐ ａ ｓ ｓ ｉ ｎ ｇ ａ ｎ ｄ ｓ ｈａ ｒ ｅ ｄ —
ｍｅ ｍ ｏｒ ｙ ｈ ａ ｓ ｓ ｅ ｖｅ ｒ ａ ｌ ａ ｄｖ ａ ｎ ｔ ａ ｇ ｅ ｓ ｗ ｈｅ ｎ ｕｓ ｅ ｄ ｉ ｎ Ｍｏ ｎｔ ｅ Ｃａ ｒ ｌ ｏ ｎ ｅ ｕ ｔ ｒ ｏｎ ｔ ｒ ａ ｎｓ ｐｏ ｒ ｔ ｃ ｏ ｄｅ，ｓ ｕｃ ｈ ａ ｓ ｆ ｉ ｔ ｔ ｉ ｎｇ ｈａ ｒ ｄ ｗａ ｒ ｅ ｏ ｆ ｄ ｉ ｓ ｔ ｒ ｉ ｂ ｕｔ ｅ ｄ — ｓ ｈａ ｒ ｅ ｄ ｃ ｌ ｕｓ ｔ ｅ ｒ ｓ， ｅ ｃ ｏ ｎｏ ｍｉ ｚ ｉ ｎｇ ｍ ｅ ｍ ｏｒ ｙ ｄｅ ｍａ ｎｄ ｏ ｆ
ｔ ｈ ａ ｔ ｔ ｈ ｅ ｈ ｙ ｂ ｒ ｉ ｄ ｐ ａ ｒ ａ ｌ ｌ ｅ ｌ ｃ ｏ ｄ ｅ ｃ ａ ｎ ｒ ｅ ａ ｃ ｈ ｇ ｏ ｏ ｄ ｐ ｅ ｒ ｆ ｏ ｒ ｍａ ｎ ｃ ｅ ｊ ｕ ｓ ｔ ａ ｓ ｐ ｕ ｒ ｅ ＭＰ Ｉ ｐ ａ ｒ ａ ｌ ｌ ｅ ｌ
Ｍｏ ｎ ｔ ｅ Ｃａ ｒ ｌ ｏ ｔ ｒ ａ ｎ ｓ ｐ ｏ ｒ ｔ ，ｉ ｍｐ ｒ ｏ ｖ ｉ ｎ ｇ ｐ ａ ｒ ａ ｌ ｌ ｅ ｌ ｐ ｅ ｒ ｆ ｏ ｒ ｍａ ｎ ｃ ｅ ，ａ ｎ ｄ Ｓ Ｏ ｏ ｎ ． ＭＰ Ｉ ／ Ｏｐ ｅ ｎ ＭＰ
ｈ ｙｂ ｒ ｉ ｄ ｐ ａ ｒ ａ ｌ ｌ ｅ ｌ ｉ ｓ ｍ ｗａ ｓ ｉ ｍ ｐｌ ｅ ｍｅ ｎ ｔ ｅ ｄ ｂ ａ ｓ ｅ ｄ ｏ ｎ ａ ｏ ｎｅ ｄ ｉ ｍｅ ｎｓ ｉ ｏ ｎ Ｍ ｏｎ ｔ ｅ Ｃａ ｒ ｌ ｏ ｎｅ ｕｔ ｒ ｏ ｎ

一种求解反应堆中子扩散方程和输运方程的方法

! 第!"卷 第#期 !")*#年#月
原子能科学技术 +,%-./0123456/.21/271892/:1%&%45
$%&'!"(%'# +;4'")*#
一种求解反应堆中子扩散方程 和输运方程的方法
胡文军曹永刚
中国原子能科学研究院 反应堆工程技术研究部北京!*)"=*<
摘要结合反应堆中子物理方程的特点和非结构网格技术提出了一种在任意网格下求 解 反 应 堆 中 子 物 理方程的方法以反应堆中子扩散方程和一阶离散纵标的输运方程为例从空间离散方 程 离 散 和 边 界 条件实施等方面介绍了该 方 法 的 实 施 过 程 利 用 该 方 法 编 制 的 程 序 计 算 了 N(K?)) 基 准 题AG@*))) 概念堆和 NU@ 栅 元 < 种 情 况 的 有 效 增 殖 因 数并 与 其 他 程 序 的 计 算 结 果 进 行 了 对 比初 步 验 证 了 该 方 法的合理性 关 键 词 反 应 堆 物 理 数 值 方 法 区 域 离 散 方 程 离 散 边 界 条 件 实 施 中图分类号9Q<?='!!!! 文献标志码+!!!文章编号*)))K?><*")*#)#K*<#*K)L &'(*)'L!<#5RI'")*L'5%;S.71')L=!
本文在研究反应堆物理扩散和输运方程的 特征后!提出 一 种 空 间 离 散 的 方 法% 不 同 于 有 限 差 分 法 和 有 限 元 方 法!该 方 法 从 方 程 积 分 形 式 出 发 !结 合 中 子 守 恒 特 性 !可 对 三 维 空 间 区 域 进行任意形状网格的离散%本文将以扩散方程 和 一 阶 离 散 纵 标 的 输 运 方 程 稳 态 形 式 为 例!从 区 域 离 散'方 程 离 散 和 边 界 条 件 实 施 等 方 面 介 绍该方法的实施过程%

一阶中子输运方程的离散纵标有限元方法研究

一阶中子输运方程的离散纵标有限元方法研究肖锋【摘要】利用离散纵标(SN)方法离散SN一阶多群中子输运方程,建立了基于剖分单元的有限元离散与非结构网格扫描方法相结合的求解模型.针对给定的三角形单元应用Galerkin变分,获得线性方程组,通过引入定解条件,求解线性代数方程组,获得该三角形单元所有节点上的角通量,然后对其他三角形单元进行扫描,从而解出所有节点处中子通量密度.根据上述理论模型,编制了相应的计算程序FEGT,对一系列例题进行校验的数值结果表明,该程序的计算精度满足要求.%With the first order neutron multi-group transport equation discretized using the discrete ordinates solution (SN ) ,a calculation model was established by using finite element method (FEM ) based on single element and sweep method based on unstructured meshes . In a preset triangular element , the flux was discretized using Galerkin variation ,and a linear system of equations was obtained .The angular flux of the nodes in the triangle element can be got by solving the linear system of equations , with offering some conditions .With the sweep method ,all triangular elements will be swept and the differential equation can be solved .Based on the above theoretical model , the program FEGT was encoded .In order to validate this program ,a series of examples were calculated .The results show that FEGT code meets the precision requirement .【期刊名称】《原子能科学技术》【年(卷),期】2017(051)003【总页数】6页(P444-449)【关键词】一阶中子输运方程;有限元方法;离散纵标【作者】肖锋【作者单位】中国核动力研究设计院核反应堆系统设计技术重点实验室,四川成都610041【正文语种】中文【中图分类】TL329非结构网格的有限元方法[1]能适应几何结构不规则问题的数值计算，国内外已开展了大量使用有限元方法计算中子输运方程的研究。

  1. 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

第七届反应堆数值计算与粒子拾运学术交流会中子输运方程ＳＰ：方法的并行计算１ 沈智军袁光伟沈隆钧目匕京应用物理与计体致学研究计算物理实脸室，周毓麟

北京１０００８８）

本文将研究中子输运方程ＳＰｚ通近方程的并行差分格式，证明了一般并行差分格式的稳定性，用数值例子检脸了一维和二维的并行格式，比较了格式的加速比，所使用的并行算法具有较好的可扩展性，并且与ｓ．方法相比，计算精度相近，但这

里给出的格式的计算ｔ少得多． 羊滋１５１中子输运方程并行计算差分方法稳定性

１引言 由于中子枪运方程的数值求解运算量惊人，加上当今大型并行计算机发展迅速，因

此，研究输运方程的并行计算方法是十分迫切的重要间题，已经得到广泛关注。在［２１中研究了标准的多群方法的并行化。在｛（３１中对求解离散纵标方程的源迭代和扩散模拟加速算法提出了空间区域分解法，其结果表明，直到子区域变得光学薄为止，空间分解对收故率影响不大。最近，文［［ｌ１中对非定常中子枪运方程研究了简化的Ｐ２（即ＳＰ２）方法，发现ＳＰ２方程对输运方程是强壮的高阶渐近逼近。在｛’］中对ＳＰ２方法提出了一类修改的Ｍａｃｓｈａｋ边界条件，并证明了所得到的ＳＰ２方程的定解间题是适定的，井构造了纯显式和纯隐式格式，证明了格式的收敛性。虽然纯隐式格式是无条件稳定的，但每一时间层需迭代求解，运算量大，尤其是在多维情形。而纯显式格式对时间步长大小有限制，如果时间步长大的话，解将出现振荡（见［’］）。因此需要研究并行差分方法． 本文将对ＳＰ２方法提出一般的具有并行本性差分格式，证明其稳定性。这里构造并行差分格式以及稳定性证明的方法可参考文献［５１和｛６１。所讨论的ＳＰ２方程以及初边值条件见［（４１。 数值试验结果表明所构造的并行差分方法具有较好的可扩展性，并且与５８方法相比，计算精度相近，但这里给出的格式的计算量少得多。由于数值例子的计算粒度较小，在四个处理器的倩形下，所得加速比并不十分理想，但在两个处理器的情形下，表明这里

的井行算法还是比较令人满意的。

２一般并行差分格式用点二ｉ二７△二（７＝０，１，－，Ｊ）剖分区间风Ｌ１，　Ｊ４ｘ＝Ｌ，记ｘ－ｚ＝ｘｏ，ｘＪ＋；二ｘ７，

｀，￥苏蔺而不￥（４１０９５２１｝．　．｝藕科学基金和大规模科学与工程计算的方法和理论资助项目３一５５数值计算方法月！１，

马砖二合〔几＋令＋：），△云为时间步长，Ｎ山二Ｔ，Ｊ和Ｎ是正整数。构造如下具有并行本性的差分格式 些俨雌豁兰十、科

些了丛十菩噬愁烂十会竺遭影笠翌Ｊ二０工，…，Ｊ（１）

十。。，裂＝。，ｊ二０，１，…Ｊ一１

（２）

十、、罗尹＝０，ｊ二。，：，…，了

󰀀󰀀。。

、

󰀀｜！。󰀀｜｜ｌｄ．．．．．．

．．．．

其中才气Ｑ（。，。勺，，拼二咧拱＋（，一ａ７）屹十、，其它记号类似定义。 （３）假定．

十卜一 叻一扫汗一 举一

，‘

－‘

口

＋

争

格式中出现的参数。，口，７浦，刀，（均一致有界。当某些。，口等于。，即可得到并行格式。 离散初值条件为

崎二功从＊＝劝易二０，夕＝Ｄ，１，…，Ｊ（’）离散边界条件为 叻戮一妇”几势豁十告二儿（劝品＋２劝升）（５）其中几是给定的非负常数。 类似地可构造两维５几方程初边值间题的具有并行本性的差分格式。

３差分格式的稳定性 本节将证明并行格式（ｌ）一（５）对任意固定的备都是稳定的。 用知孙此十、和知黔分别乘以方程（ｌ），（ｚ）和（ｓ）１将所得三式分别对ｊ求和，再将它们加起来，得忐力心王一、，、十态到（‘拼一‘冷命（‘群一心可叽＋、）

扩护笋＿话”十盯 “十全‘工Ｊ一古

３△公哈＋艺粤、扩寸屹十云馨（‘狱一‘城’心‘

十忿剧（‘黔一俨、・＊十（‘拼一‘冷叼＋（‘摄一心沁・蕙，一裂，尔、・念睿（，叮二碧）、艺。灯，扩叮吩铭第七属反应堆数值计界与粗子翰运学术报告二；Ｉ．，｝乙ｎ＋（｝ｗ０７４７ｊ＝０

注意到（Ｙ＇０＋１一０；　）Ｖ，Ｏｊ＝蠢（ＩＶ，ｏｉ　１１一｝，ＰＯｊ　１一｝，Ｐｎ０＋１一Ｖ＇Ｏｊ　１），

（ ０　．＋．ｊ１ｊｊ一‘＋）３ｊ＋ｏｋ．＂＋，６１，一－＋ＭＡｂｏｊ　，　）ＩｆｉＩｊ＋ｉ二（，Ｏ－，ｊ　＋、一，＃，－ｊ一尹凡＋　（０ｏｊ十，一吩）屹＋、十喇举ｎ＋１｝｝（｝ｌｉ＋；一０ｉｊ＋Ｉ）１Ｇｏ－ｊ卜－Ｒ；　（，ｋｎ＋１一Ｖ＇，ｊ－　ｉ　）＇ｈｏ　ｊ十ｎ｝　（＇｝ｏｊ＋ｌ｝一＇Ｐｏｉ＋１）＇０１　ｊ十＊一Ｑｉ　（，ｏ０１１一ｐｏｉ　）＇ｉｉｉ＋｝，　， 一。｝：２　（１ｘ１２＋；。一），Ｖａ，　ｂ，并对所有出现的亡方向差商项，不等式

如ｍｏ　＊＇－ｓａ；，再应用Ｐｆ方程组（：）一（３）。于Ａ％出如下

众（１１ｏｎ　１１１１＋３１１＋Ｇｉｎ　１１１２＋５１１Ｇｚｎ　１１１，一（！。，；！。‘＋３１｝ｉｈ１１＋５１１Ｎ２ｎｎ１２：。（备）（１１４＂󰀀＋０：　ＩＩ，十“＇ｈｏｔｉ　１１１　＋３１０　ｉｎ　１１１，十５　１１０　ｓｅ　１１１，十“＋ｈｏｎｌ２十、１Ｙ１ｅ１１２十５１１０２ｈ１１２）．利用离散Ｇｒｏｎｗａｌｌ不等式，推出如下不等式成立０＜ｎ＜Ｎ｝１１　｝＇Ｏｈ　１２　＋　１１＋Ｌ１ｎ　１１２　＋　１Ｖ，：＂　１１２：Ｃ艺１ｐｎ＋０＂　１１２Ａｔ．（６）

因此，巳经证明对任意固定的器，具有并行本性的差分格式（１）－（５）是稳定的． 与一维间题类似，对二维和三维ＳＰ方程可构造具有并行本性的差分格式，并证明其稼定性。

４救值结果 我们已经计算了一维平几何和二维Ｘ一Ｙ几何Ｓ乃方程并行格式的数值解。实际计算的一维并行格式是在差分格式（１）－（５）中取一个显式点，将区域分成两个不重亚的子区域。所计算的二维并行格式是将区域分成四块不重处的子区域，在子区域的交界线上显式计算。下面仅给出二维问题的计算结果，其中。：＝　１．０，　ｖ。二０．１，Ｑ＝　１．０，空间正方形区域的边长为１０。初始通量为０，区域左边和下边的边界条件为反射条件，右边和上边的边界条件是真空条件。３一５７数位计算方法闷

！｜！翎』．ＪＪ．『１，．１｜Ｊ飞１（２一一１

１１１１

时间步数２００２００２００４００２００２００２００时间步长０．０５ ０．１０．０３０．０３空间步数 ４　ＣＰＵ时间（ｓ） ２　ＣＰＵ时间（ｓ）串行时间（ｓ）加速比（４　ＣＰＵ）加速比

ＣＰＵ１１巧石ｎ‘０甘ｎ曰３４ｎＪ介Ｊｌ』Ｊ９９９５６３０２４３蛇￣日甘，Ｊ月０．．占目』ｄ马」‘．１１，ｌｑ白，曰介」０．０５０．０３０．０４３．６６７．４８９．０５１５．１８１３．１２１５．９４２１．５１５．９５７．８４１５．５７１９．６９３７．１１３１．６２３４．７８４９．１９２．１４２．０８２．１７２．４４２．４１２．１８２．２８｝一１｛七口

户勺

‘〕

已ｄ

已Ｊ

工合

民０

４４ｔ‘间了月才Ｑ户

ｎＵ

从表中可以看出，当空间步长固定，随着时间步长增加，加速比随之增加。其原因是，当时间步长增加时，求解的迭代次数增加引起计算规模增大。

参考文献１．　Ｕ．　Ｓｈｉｎ　ａｎｄ　Ｗ．　Ｆ．　Ｍｉｌｌｅｒ，　Ｊｒ．，　Ｔｈｅ　Ｔｉｍｅ－Ｄｅｐｅｎｄｅｎｔ　Ｓｉｍｐｌｉｆｉｅｄ几Ｅｑｕａｔｉｏｎｓ：　Ａｓｙｍｐ｝ｔｏｔｉｃ　Ａｎａｌｙｓｅｓ　ａｎｄ　Ｎｕｍｅｒｉｃａｌ　Ｅｘｐｅｒｉｍｅｎｔｓ，　Ｎｕｃｌ．　Ｓｃｉ．　Ｅｎｇ．，　１２８，　２７－４６，　１９９８．

Ｂ．　Ｒ．　Ｗｉｅｎｋｅ　ａｎｄ　Ｒ．　Ｅ．　Ｈｉｒｏｍｏｔｏ，　Ｐａｒａｌｌｅｌ　Ｓ󰀀　Ｔｒａｎｓｐｏｒｔ　Ａｌｇｏｒｉｔｈｍｓ，　Ｔｒａｎｓｐ．　Ｔｈｅｏｒｙ女Ｓｔａｔ．　Ｐ妙ｓ．，　１５（１＆２），　４９－５９，　１９８６．

Ｍ．　Ｙａｖｕｚ　ａｎｄ　Ｅ．　Ｗ．　Ｌａｒｓｅｎ，　Ｓｐａｔｉａｌ　Ｄｏｍａｉｎ　Ｄｅｃｏｍｐｏｓｉｔｉｏｎ　ｆｏｒ　Ｎｅｕｔｒｏｎ　ＴｒａｎｓｐｏｒｔＰｒｏｂｌｅｍｓ，　Ｔｒａｎｓｐ．　Ｔｈｅｏｒｙ＆Ｓｔａｔ．　Ｐｈｙｓ．，　１８（２）；　２０５－２１９，　１９８９．

４袁光伟， １９９８．沈智军，沈隆钧，周毓麟，中子抬运方程ＳＰ；方法稳定性分析，Ｐｒｅｐｒｉｎｔ，

５．　Ｚｈｏｕ　Ｙｕ－ｌｉｎ，　Ｙｕａｎ　Ｇｕａｎｇ－ｗｅｉ；　Ｇｅｎｅｒａｌ　ｄｉｆｅｒｅｎｃｅ　ｓｃｈｅｍｅｓ　ｗｉｔｈ　ｉｎｔｒｉｎｓｉｃ　ｐａｒａｌｌｅｌｉｓｍ ｆｏｒ　ｎｏｎｌｉｎｅａｒ　ｐａｒａｂｏｌｉｃ　ｓｙｓｔｅｍｓ，　Ｓｃｉｅｎｃｅ　ｉｎ　Ｃｈｉｎａ　（ｓｅｒ．　Ａ），　Ｖｏｌ．　４０，　Ｎｏ．　４，　３５７－３６５， １９９７．

６．　Ｚｈｏｕ　Ｙｕ－＂，　Ｓｈｅｎ　Ｌｏｎｇ　ｊｕｎ，　Ｙｕａｎ　Ｇｕａｎｇ－ｗｅｉ，　Ｓｏｍｅ　ｐｒａｃｔｉｃａｌ　ｄｉｆｅｒｅｎｃｅ，　ｓｃｈｅｍｅｓ　ｗｉｔｈ ｉｎｔｒｉｎｓｉｃ　ｐａｒａｌｌｅｌｉｓｍ　ｆｏｒ　ｎｏｎｌｉｎｅａｒ　ｐａｒａｂｏｌｉｃ　ｓｙｓｔｅｍｓ，　Ｃｈｉｎｅｓｅ　Ｊ．　Ｎｕｍｅｒ．　Ｍａｔｈ．　ａｎｄ Ａｐｐｌ．，　Ｖｏｌ．　１９，　Ｎｏ．　３，　４６－５７，　１９９７．

３一５８