CN106202613B - 一种适用于耦合输运计算的大规模并行方法 - Google Patents

一种适用于耦合输运计算的大规模并行方法 Download PDF

Info

Publication number
CN106202613B
CN106202613B CN201610472928.9A CN201610472928A CN106202613B CN 106202613 B CN106202613 B CN 106202613B CN 201610472928 A CN201610472928 A CN 201610472928A CN 106202613 B CN106202613 B CN 106202613B
Authority
CN
China
Prior art keywords
parallel
dimensional
region
carried out
characteristic curve
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610472928.9A
Other languages
English (en)
Other versions
CN106202613A (zh
Inventor
刘宙宇
赵晨
郑友琦
李云召
吴宏春
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201610472928.9A priority Critical patent/CN106202613B/zh
Publication of CN106202613A publication Critical patent/CN106202613A/zh
Application granted granted Critical
Publication of CN106202613B publication Critical patent/CN106202613B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Complex Calculations (AREA)

Abstract

一种适用于耦合输运计算的大规模并行方法，将径向特征线方法计算、一维离散纵标方法计算和粗网有限差分加速计算均拓展为与径向二维特征线方法相对应的三重并行，彻底消除不匹配现象；本发明方法将在二维一维耦合输运方法的基础上，实现多种方法的并行，并解决多种耦合算法并行不匹配带来的资源浪费问题。

Description

一种适用于耦合输运计算的大规模并行方法
技术领域
本发明涉及核反应堆堆芯设计和反应堆物理计算领域，具体涉及一种适用于耦合输运计算的大规模并行方法。
背景技术
随着计算机计算能力的日益提高，反应堆设计领域希望采用更精确的设计模型，最大程度的减少之前反应堆计算过程中做的人为近似和经验修正，实现高保真度的反应堆堆芯物理计算，从而对反应堆进行精确分析模拟。这种新一代反应堆计算方法称为数值反应堆技术，被认为是未来进行反应堆分析的最重要技术，对于反应堆的安全分析和新型反应堆的设计改进具有十分重要的意义，成为我国核电技术能力达到世界领先水平提供必不可少的基础。
所谓“高保真计算”，即抛弃传统的栅元均匀化、组件均匀化等一系列均匀化过程，对于反应堆进行直接一步计算。这种做法可以避免原有方法无法考虑真实堆芯中组件之间的相互影响，以及组件的燃耗历史效应等等一系列问题。另外，通过对反应堆进行精细建模，充分考虑各种精细效应，可以从模型上提高反应堆计算的精度。
然而现有的计算机水平还不足以直接进行高保真的三维全堆芯直接计算，目前采取的折中办法是采用二维一维或二维三维耦合方法，进行耦合输运计算。即使这种方法已经大幅度减小了计算时间，然而在实际反应堆计算中仍然是极其耗时的，特别是在全堆多群计算的情况下，仍然难以实用。
目前发展中的全堆芯三维输运方法大部分都基于单一变量实现并行，如能群、角度、空间或特征线，但其并行度有限，无法充分利用当前最先进的超级计算机资源。另外，由于二维/一维耦合算法中径向二维特征线方法计算包括空间、角度和特征线可并行，一维计算仅含空间、角度可并行，而粗网加速方法只有空间可并行，三者并行自由度的不匹配将导致一维和粗网加速计算时存在严重的资源浪费现象。
因此需要发明一种适用于耦合输运计算的大规模并行方法。
发明内容
为了实现耦合输运计算的大规模并行，本发明的目的在于提供一种适用于耦合输运计算的大规模并行方法，该方法将在二维一维耦合输运方法的基础上，实现多种方法的并行，并解决多种耦合算法并行不匹配带来的资源浪费问题。
为了实现上述目的，本发明采用如下技术方案：
一种适用于耦合输运计算的大规模并行方法，将径向二维特征线方法计算、轴向一维离散纵标方法计算和粗网有限差分加速计算均拓展为与径向二维特征线方法相对应的三重并行，彻底消除不匹配现象；具体步骤如下：
步骤1：径向二维特征线方法计算包括空间、角度和特征线三个维度，每个维度均能够进行并行；首先将整个计算区域分解为各个子区域，将整个特征线扫描区域按照特征线计算量均分给各个处理器；其次，在每个区域内部再按照角度将特征线计算分给各个处理器，从而实现第二重并行；最后在每个处理器内部通过共享式内存的方式，进行多线程的特征线并行，即第三重并行；
步骤2：一维轴向离散纵标法计算包括栅元和角度两个维度，对一维区域进行区域分解，实现三重并行度；首先进行栅元并行，在一维计算方面将各个栅元计算区域均分给各个处理器；其次进行栅元内部的一维区域分解并行，离散纵标扫描分层进行；最后在各个计算子区域内进行方向并行，在各个子区域内部进行各个方向的方向并行；
步骤3：粗网有限差分CMFD仅存在空间变量的并行，将在利用区域分解的广义最小残差算法GMRES求解粗网有限差分线性方程的基础上，利用径向特征线方法计算时角度和特征线的并行计算资源加速GMRES算法中的矩阵向量积计算，实现二维特征线方法、一维轴向离散纵标法和三维粗网有限差分三者的并行匹配；首先进行区域分解并行，将计算区域均分给各个处理器，在每个区域内部进行粗网有限差分计算，利用计算资源进行矩阵向量积的并行计算。
和现有技术相比较，本发明具备如下优点：
优点1：耦合算法中，径向二维特征线方法计算通常只进行最简单的层并行，并行规模十分有限，无法用于大规模计算；本发明利用空间、角度和特征线三个维度的并行，结合现有大型并行计算机和特征线方法的特点，合理选择分布式内存和共享式内存结构，满足内存限制的情况下尽量提高并行效率，从而实现大规模计算。
优点2：在耦合计算中进行多种耦合算法的匹配，从而避免了计算的浪费。常用的一维轴向离散纵标法计算包括栅元和角度两个维度，本方法在此基础上对一维区域进行区域分解，实现三重并行度，从而与二维特征线方法计算实现匹配。
优点3：在粗网有限差分算法中，本发明加入了广义最小残差算法GMRES中的矩阵向量积并行计算，从而可以更充分利用计算资源。
附图说明
图1为二维三重并行示意图。
图2为一维三重并行示意图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细说明：
考虑到现有并行机的特点，采用分布式内存结构结合共享式内存结构进行计算，充分利用两种并行方案的特点。分布式内存结构最常用的是信息传递接口(MPI)，共享式内存结构最常用的是共享存储并行编程(OpenMP)。MPI具有并行规模大的优点，OpenMP并行规模因为共享内存受限，但具有在节点内部避免通信的优点。在选择通信时需要合理考虑两种并行模式的优缺点，选择合适的并行策略。
下面对本发明一种适用于耦合输运计算的大规模并行方法各步骤进行详细说明：
径向二维特征线方法计算包括空间、角度和特征线三个维度，如图1所示。径向二维特征线方法计算首先采用区域分解算法进行空间并行，每个计算区域作为一个通信域，通信域间通过MPI传递边界角通量；其次在通信域内部采用角度并行，扫描边界上传递反射角通量，通信域内通过MPI进行通信；最后在每个处理器内部通过OpenMP的方式实现多线程特征线并行计算，特征线计算占用内存较少，便于进行共享式内存的并行。
常用的一维轴向离散纵标法计算包括栅元和角度两个维度，本方法在此基础上对一维区域进行区域分解，实现三重并行度；如图2所示，首先进行栅元并行，在一维计算方面将各个栅元计算区域均分给各个处理器；其次进行栅元内部的一维区域分解并行，离散纵标扫描分层进行；最后在各个计算子区域内进行方向并行，在各个子区域内部进行各个方向的方向并行；一维轴向离散纵标法在大规模问题上，每个栅元作为通信域，因为一维计算栅元与栅元之间相互独立，因此通信域之间不需要进行通信；每个栅元内部进行区域分解计算，通过MPI通信边界角通量，在通信域内部进行通信。每个处理器内部开多线程，每个线程进行某一方向的扫描计算，通过OpenMP方式进行共享式内存的并行。常用的粗网有限差分CMFD仅存在空间变量的并行，本发明将在利用区域分解的广义最小残差算法GMRES求解粗网有限差分线性方程的基础上，利用径向特征线方法计算时角度和特征线的并行计算资源加速GMRES算法中的矩阵向量积计算，实现二维特征线方法、一维轴向离散纵标法和三维粗网有限差分三者的并行匹配。在空间上首先采用区域分解算法，每个计算区域作为一个通信域，通信域间通过MPI进行通信；在每个区域内部，通过GMRES算法进行矩阵向量积计算，利用径向MOC计算时角度和特征线的并行计算资源，实现大规模矩阵的并行计算。

Claims (1)

1.一种适用于耦合输运计算的大规模并行方法，其特征在于：将二维特征线方法计算、一维离散纵标方法计算和三维粗网有限差分加速计算均拓展为与二维特征线方法相对应的三重并行，彻底消除不匹配现象；具体步骤如下：
步骤1：二维特征线方法计算包括空间、角度和特征线三个维度，每个维度均能够进行并行；首先将整个计算区域分解为各个子区域，将整个特征线扫描区域按照特征线计算量均分给各个处理器；其次，在每个区域内部再按照角度将特征线计算分给各个处理器，从而实现第二重并行；最后在每个处理器内部通过共享式内存的方式，进行多线程的特征线并行，即第三重并行；
步骤2：一维离散纵标方法计算包括栅元和角度两个维度，对一维区域进行区域分解，实现三重并行度；首先进行栅元并行，在一维计算方面将各个栅元计算区域均分给各个处理器；其次进行栅元内部的一维区域分解并行，一维离散纵标扫描分层进行；最后在各个栅元计算区域内进行方向并行；
步骤3：三维粗网有限差分仅存在空间变量的并行，在利用区域分解的广义最小残差算法GMRES求解粗网有限差分线性方程的基础上，利用二维特征线方法计算时角度和特征线的并行计算资源加速广义最小残差算法GMRES中的矩阵向量积计算，实现二维特征线方法、一维离散纵标方法和三维粗网有限差分三者的并行匹配；首先进行区域分解并行，将计算区域均分给各个处理器，在每个区域内部进行三维粗网有限差分计算，利用计算资源进行矩阵向量积的并行计算。
CN201610472928.9A 2016-06-24 2016-06-24 一种适用于耦合输运计算的大规模并行方法 Active CN106202613B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610472928.9A CN106202613B (zh) 2016-06-24 2016-06-24 一种适用于耦合输运计算的大规模并行方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610472928.9A CN106202613B (zh) 2016-06-24 2016-06-24 一种适用于耦合输运计算的大规模并行方法

Publications (2)

Publication Number Publication Date
CN106202613A CN106202613A (zh) 2016-12-07
CN106202613B true CN106202613B (zh) 2019-05-21

Family

ID=57460796

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610472928.9A Active CN106202613B (zh) 2016-06-24 2016-06-24 一种适用于耦合输运计算的大规模并行方法

Country Status (1)

Country Link
CN (1) CN106202613B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108734181B (zh) * 2018-05-23 2019-04-16 西安交通大学 一种加速在线生成核反应堆堆芯特征线的方法
CN110472187B (zh) * 2019-08-06 2021-05-04 中国原子能科学研究院 一种三维中子输运特征线法的负载均衡并行方法
CN112800641B (zh) * 2020-12-21 2022-02-18 中国核动力研究设计院 一种基于区域分解并行的广义粗网有限差分加速方法
CN112783511B (zh) * 2021-02-05 2023-04-11 成都信息工程大学 一种栅元少群参数计算模块程序的优化方法、系统、终端

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103150424A (zh) * 2013-02-05 2013-06-12 西安交通大学 一种获取反应堆堆芯三维中子通量密度精细分布的方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103150424A (zh) * 2013-02-05 2013-06-12 西安交通大学 一种获取反应堆堆芯三维中子通量密度精细分布的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Solution of the C5G7MOX benchmark three-dimensional extension problems by the DeCART direct whole core caculation code;Jin Young Cho 等;《Progress in Nuclear Energy》;20060731;第48卷(第5期);第456-466页
基于大规模并行计算的2D/1D耦合三维全堆输运程序Tiger-3D;吴文斌 等;《核动力工程》;20141215;第35卷(第S2期);第135-139页
基于模块化特征线方法的二维/一维耦合输运程序开发;梁亮 等;《核动力工程》;20141215;第35卷(第S2期);第131-134页

Also Published As

Publication number Publication date
CN106202613A (zh) 2016-12-07

Similar Documents

Publication Publication Date Title
CN106202613B (zh) 一种适用于耦合输运计算的大规模并行方法
CN110110413B (zh) 一种基于材料场缩减级数展开的结构拓扑优化方法
CN104036095B (zh) 基于区域分解的耦合高精度复杂外形流场快速算法
Gunow et al. Full core 3D simulation of the BEAVRS benchmark with OpenMOC
CN113158527B (zh) 一种基于隐式fvfd计算频域电磁场的方法
CN113158492B (zh) 一种时变电磁场的全隐式双时间步计算方法
Fang et al. An algorithm–hardware co-optimized framework for accelerating n: M sparse transformers
CN107038294A (zh) 针对轻水堆的基于等效一维棒模型的共振自屏计算方法
Sgattoni et al. Optimising piccante-an open source particle-in-cell code for advanced simulations on tier-0 systems
Yamazaki et al. Performance comparison of parallel eigensolvers based on a contour integral method and a Lanczos method
CN106096183B (zh) 一种基于特征线法的多重并行方法
CN107038293A (zh) 针对板状燃料的基于等效一维板模型的共振自屏计算方法
Hao et al. 3D whole-core neutron transport simulation using 2D/1D method via multi-level generalized equivalence theory based CMFD acceleration
Dziekonski et al. Communication and load balancing optimization for finite element electromagnetic simulations using multi-GPU workstation
Bao et al. Parameter estimation method based on parameter function surface
Hou et al. Fast density inversion solution for full tensor gravity gradiometry data
Liang et al. Overlapping communication and computation of GPU/CPU heterogeneous parallel spatial domain decomposition MOC method
CN111339688B (zh) 基于大数据并行算法求解火箭仿真模型时域方程的方法
Galvez et al. Accelerating the solution of families of shifted linear systems with CUDA
Liu et al. A new three-level CMFD method based on the loosely coupled parallel strategy
CN108879691A (zh) 一种大规模连续潮流计算的方法及装置
Agullo et al. Recent advances in sparse direct solvers
Wang et al. Three-dimensional variational nodal method parallelization for pin resolved neutron transport calculations
CN106202862A (zh) 一种针对压水堆栅元非均匀共振积分表的制作方法
Xu et al. Parallelizing a high-order CFD software for 3D, multi-block, structural grids on the TianHe-1A supercomputer

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant