基于几何区域分解的三维输运问题并行迭代算法-北京应用物理与计算.pdfVIP

第 21 卷 第 1 期 计 　　算 　　物 　　理 Vol. 21 ,No. 1 　 2004 年 1 月 　CHINESE JOURNAL OF COMPUTATIONAL PHYSICS 　 Jan. , 2004 ( ) [文章编号] 　1001246X 2004 基于几何区域分解的三维输运问题并行迭代算法 阳述林 , 　莫则尧 , 　沈隆钧 (北京应用物理与计算数学研究所 　计算物理实验室 , 北京 　100088) [摘 　要] 　对三维直角坐标下的输运隐式差分方程 , 研究了基于几何区域分解的并行迭代算法 , 给出了串、并 行迭代误差估计. 并对相关数值结果进行了分析 、比较. [ 关键词] 　输运方程 ; 区域分解 ; 并行迭代 [ 中图分类号] 　O246 [文献标识码] 　A 0 　引言 高维粒子输运数值计算是一项复杂的大型科学计算 ,需要高性能计算机和高效率计算方法才能实现 ,而 并行算法就是重要的高效方法之一. 随着大型并行计算机的不断发展 ,特别是国产大型多 CPU 并行机的研 制成功 ,使非定态中子输运计算可以逐步从二维走向三维. 对三维中子输运方程的高效并行算法研究随之提 到了日程上来. ( ) 目前国际上应用较广的并行扫描 Sweep 算法 ,最早是 1992 年由 Koch 、Baker 和 Alcouffe 在正交六面体上 作列分解计算输运方程而提出的[1 ,2 ] ,称为 KBA 算法. 近年来 ,扫描算法得到进一步发展 ,在算法设计中采用 不同的区域分解方法 ,并应用到高维几何空间多种网格上的输运方程数值求解[3 ,4 ] . 由于扫描算法在实际计算中将数值求解过程与数据传送过程耦合在一起 ,需要针对具体并行机处理负 载平衡和极小化通信问题 ,但各处理器数据的非完全独立使并行度受到限制. 另外 ,该算法适用于 MPI 环境 [4 ] 下进行并行计算 ,其并行效率和可扩展性受到通信延迟的极大影响 . 本文试图通过几何区域分解 ,实现按区域分块进行中子输运方程的同步并行迭代计算和块与块之间的 并行数据传送. 这样构造的并行算法及其程序实现 ,具有不依赖于具体并行计算机的内在本性 ,是提高算法 并行度及可扩展性的一条新路. 同时 ,通过在 OpenMP 环境下的数值检测 ,证明该并行迭代计算是方便和高 效的. 1 　三维中子输运方程及其离散格式 三维直角坐标下的非定态多群中子输运方程为 φ φ φ φ 1 g g g g t + ξ + μ + η + σφ = 　　　 v t x y z g g g g 1 2π