HE Qin,YANG Hao,ZHOU Shijie,et al.Parallel Optimization of MOC/SN Coupled Three-Dimensional Neutron Transport Program[J].Journal of Chengdu University of Information Technology,2024,39(01):37-42.[doi:10.16836/j.cnki.jcuit.2024.01.007]
MOC/SN耦合三维中子输运程序的并行优化
- Title:
- Parallel Optimization of MOC/SN Coupled Three-Dimensional Neutron Transport Program
- 文章编号:
- 2096-1618(2024)01-0037-06
- 关键词:
- KYCORE程序; MPI二维拓扑通信结构; 特征线并行算法
- Keywords:
- KYCORE program; MPI 2D topology communication structure; parallel algorithm of method of characteristics
- 分类号:
- TP301.6
- 文献标志码:
- A
- 摘要:
- KYCORE程序是中国核动力研究设计院开发的径向特征线法与轴向离散纵标法耦合三维中子输运程序。KYCORE程序采用几何区域并行方法,可并行计算不同区域的数据。但由于其庞大的计算量以及内存占用量,KYCORE基于几何区域并行取得的时间增益不高,且程序无法再对区域中输运模块二级并行。MPI二维拓扑通信结构将不同维度的计算划分至不同通信域,适用于二维计算问题。将其应用至KYCORE程序,使几何区域和输运模块独立开,实现了在不破坏原有一维几何区域并行结构的情况下,可以进一步实现区域内部输运模块即二维维度上的并行优化。针对输运模块中的特征线法求解模块,特征线两级并行优化算法自适应地为各个进程划分能群计算任务,且各个进程动态负载均衡。针对C5G7基准例题测试,对于不同大小的算例,输运特征线法模块的并行效率达70%及以上,测试表明该结构与方法取得了明显的加速效果,具有精确性、高可扩展性。
- Abstract:
- KYCORE program is a three-dimensional neutron transport program developed by Nuclear Power Research Institute of China, which is coupled with radial characteristic line method and axial discrete longitudinal scale method. KYCORE program adopts geometric domain parallel method, which can compute data of different domains in parallel. However, due to its huge computational load and memory occupation, the time gain obtained from KYCORE based on geometric domain parallelism is insufficient, and the program can no longer parallel the transport module level in the domain. To address this issue, MPI 2D topological communication structure divides the computation of different dimensions into different communication domains, which is suitable for 2D computation problems. By applying it to KYCORE program, the geometric domain and the transport module become independent, and the parallel optimization of the transport module within the domain, namely the two-dimensional dimension, can be further realized without destroying the original parallel structure of the one-dimensional geometric domain. For the method of characteristics module in the transport module, the two-level parallel optimization algorithm adaptively divides the energy group computing tasks for each process, each process dynamically balances the load., Finally, for the C5G7 benchmark example test, the parallel efficiency of the transport of method of characteristics module reaches 70% or more for the cases of varing sizes.
参考文献/References:
[1] Mokhtar-Kharroubi M.Mathematical topics in neutron transport theory:new aspects[M].World Scientific,1997.
[2] 唐霄,李庆,柴晓明,等.MOC/SN耦合三维中子输运程序KYCORE开发与初步验证[J].强激光与粒子束,2017,29(3):5.
[3] Zhang G,Li Z,Liu J.A linearized coarse mesh finite difference preconditioner for the within-group Krylov subspace iteration based on two-dimensional method of characteristics[J].Annals of Nuclear Energy,2020,151.
[4] Masiello E.“On the fly”stabilization of the Coarse-Mesh Finite Difference acceleration for multidimensional discrete-ordinates transport calculations[J].Journal of Computational Physics,2018,373:1-27.
[5] Shen W.On the better performance of the coarse-mesh finite-difference method for CANDU-type reactors[J].Annals of Nuclear Energy,2012,46:169-178.
[6] Shapira Y.Domain decomposition multigrid for unstructured grids[R].Los Alamos National Lab(LANL),Los Alamos,NM(United States),1997.
[7] Toselli A,Widlund O B.Domain Decomposition Methods——Algorithms and Theory[M].北京:科学出版社,2006.
[8] Zhao J,Lianghui P,Zhaoyuan L,et al.Calculation and Analysis of C5G7 Benchmark by SN2D Code[J].Atomic Energy ence & Technology,2013,47(增刊).
[9] Forum T M.MPI:a message passing interface[M].Oregon Graduate Institute School of Science & Engineering,1994.
[10] Gropp,Lusk W,Skjellum E,et al.Using MPI-Portable Parallel Programming with the Message Passing Interface 2e[J].scientific programming,2014,5(3):275-276.
[11] Chandra R,Dagum L,Kohr D,et al.Parallel Programming in OpenMP[M].Morgan Kaufmann Publishers,2001.
[12] Hoskin N E,Abbott M B.An Introduction to the Method of Characteristics[J].Mathematical Gazette,1968,52(380):207.
[13] Ganapol B D.A more efficient implementation of the discrete-ordinates method for an approximate model of particle transport in a duct[J].Annals of Nuclear Energy,2015,86:13-22.
[14] Sayyadikord Abadi R,Shojaei A F,Tatafei F E,et al.Theoretical Study of Octreotide Derivatives as Anti-Cancer Drugs using QSAR,Monte Carlo Method and formation of Complexes[J].Russian Journal of Physical Chemistry B,2022,16(1):127-137.
[15] Karpeev A.Dynamic load balancing algorithm for continuum mechanics problems with essential redistribution of workloads among the processes[J].Journal of Physics:Conference Series,2022,2154(1):12007.
[16] Pacheco P.An introduction to parallel programming[M].Elsevier,2011.
备注/Memo
收稿日期:2023-03-28
基金项目:四川省科技厅重点研发资助项目(2022YFS0542)