二维柱几何非定常中子输运方程基于格式的界面预估校正并行算法

计算物理 ›› 2012, Vol. 29 ›› Issue (2): 198-204.

二维柱几何非定常中子输运方程基于格式的界面预估校正并行算法

魏军侠, 袁光伟, 阳述林, 申卫东

北京应用物理与计算数学研究所, 北京 100088

收稿日期:2011-05-10 修回日期:2011-08-05 出版日期:2012-03-25 发布日期:2012-03-25
作者简介:魏军侠(1975-),女,副研究员,从事中子输运方程计算方法与并行算法研究,E-mail:wei_junxia@iapcm.ac.cn
基金资助:
基础科研计划(B1520110011);国家重点基础研究专项经费(2011CB309705);中国工程物理研究院科学技术发展基金(2008A0202008)资助项目

A Parallel Algorithm with Interface Prediction and Correction for Time-dependent Transport Equation in 2D Cylindrical Geometry

WEI Junxia, YUAN Guangwei, YANG Shulin, SHEN Weidong

Institute of Applied Physics and Computational Mathematics, Beijing 100088, China

Received:2011-05-10 Revised:2011-08-05 Online:2012-03-25 Published:2012-03-25

摘要/Abstract

摘要： 针对二维柱几何非定常中子输运方程的Sn-间断有限元方法,提出基于格式的界面预估校正并行算法.数值算例表明,该并行算法在精度与并行度等诸方面均具有良好的性质,与已有的基于隐式格式的并行扫描算法相比,对于二维中子输运大规模计算问题,并行计算效率较高,并行加速比可增加-倍以上,且可保持原隐式格式的计算精度.

关键词: 中子输运方程, 间断有限元方法, 区域分解, Sn扫描算法, 界面预估校正

Abstract: For discontinuous finite element method of time-dependent neutron transport equations in 2D cylindrical geometry,a parallel algorithm with interface prediction and correction is designed.Numerical experiments show that the parallel algorithm with explicit prediction and implicit correction has good precision,parallelism and simplicity.Compared with parallel algorithm based on implicit scheme,the new parallel algorithm has higher parallel efficiency with accuracy-preserving.Especially,it achieves super-linear speedup on hundreds of processors for large-scale problems.

Key words: transport equation, discontinuous finite element, domain decomposition, sweep algorithm, interface prediction and correction

中图分类号:

O246

魏军侠, 袁光伟, 阳述林, 申卫东. 二维柱几何非定常中子输运方程基于格式的界面预估校正并行算法[J]. 计算物理, 2012, 29(2): 198-204.

WEI Junxia, YUAN Guangwei, YANG Shulin, SHEN Weidong. A Parallel Algorithm with Interface Prediction and Correction for Time-dependent Transport Equation in 2D Cylindrical Geometry[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2012, 29(2): 198-204.

[1]	蔡颖, 张存波, 刘旭, 范征锋, 刘元元, 徐小文, 张爱清. 二维球坐标系中子输运方程的一种并行SN算法[J]. 计算物理, 2022, 39(2): 143-152.
[2]	尹亮, 杨超, 马石庄. 地球外核热对流运动的区域分解多重网格并行数值模拟[J]. 计算物理, 2019, 36(1): 1-14.
[3]	田炜, 任新成. 雪层覆盖地面电磁散射的矩量法[J]. 计算物理, 2018, 35(2): 205-211.
[4]	赵国忠, 蔚喜军, 郭怀民. 二维Lagrangian坐标系下可压气动方程组的间断Petrov-Galerkin方法[J]. 计算物理, 2017, 34(3): 294-308.
[5]	李刚, 张宝印, 邓力, 上官丹骅, 李瑞, 马彦, 付元光, 胡小利. 蒙特卡罗区域分解并行计算中确保串并行结果一致的伪随机数应用[J]. 计算物理, 2017, 34(1): 67-72.
[6]	叶珍宝, 周海京. 高阶间断伽辽金时域有限元方法分析三维谐振腔[J]. 计算物理, 2015, 32(4): 449-454.
[7]	赵国忠, 蔚喜军, 李珍珍. 多介质流模拟的Runge-Kutta控制体积间断有限元方法[J]. 计算物理, 2014, 31(3): 271-284.
[8]	冯涛, 蔚喜军, 安恒斌, 崔霞, 吴迪, 李珍珍. Lax-Wendroff时间离散的自适应间断有限元方法求解三维可压缩欧拉方程[J]. 计算物理, 2013, 30(6): 791-798.
[9]	陈学, 孙创, 夏新林. 基于区域分解和逐级光线分裂的杂散光分析[J]. 计算物理, 2012, 29(6): 876-880.
[10]	梅立泉, 方叶, 原长琦. 中子输运方程源项反演的瑞利商迭代算法[J]. 计算物理, 2012, 29(3): 369-374.
[11]	陈大伟, 秦承森, 王裴, 孙海权, 蔚喜军. 凝聚介质中斜激波的反射[J]. 计算物理, 2011, 28(6): 791-796.
[12]	盛志明, 崔霞, 刘兴平. 扩散方程区域分解的多步算法[J]. 计算物理, 2011, 28(6): 825-830.
[13]	尚月强, 何银年. 非定常Navier-Stokes方程基于完全重叠型区域分解的有限元并行算法[J]. 计算物理, 2011, 28(2): 181-187.
[14]	吴迪, 蔚喜军. 自适应间断有限元方法求解三维欧拉方程[J]. 计算物理, 2010, 27(4): 492-500.
[15]	张磊, 袁礼. 龙格库塔间断有限元方法在计算爆轰问题中的应用[J]. 计算物理, 2010, 27(4): 509-517.