基于高效特征线追踪和实时生成的直接三维特征线计算方法

俞荣君, 杨睿, 胡赟, 赵守智

俞荣君, 杨睿, 胡赟, 赵守智. 基于高效特征线追踪和实时生成的直接三维特征线计算方法[J]. 原子能科学技术, 2025, 59(1): 118-126. DOI: 10.7538/yzk.2024.youxian.0115
引用本文: 俞荣君, 杨睿, 胡赟, 赵守智. 基于高效特征线追踪和实时生成的直接三维特征线计算方法[J]. 原子能科学技术, 2025, 59(1): 118-126. DOI: 10.7538/yzk.2024.youxian.0115
YU Rongjun, YANG Rui, HU Yun, ZHAO Shouzhi. Direct 3D Method of Characteristic Based on Efficient Characteristic Tracing and On-the-fly Generation[J]. Atomic Energy Science and Technology, 2025, 59(1): 118-126. DOI: 10.7538/yzk.2024.youxian.0115
Citation: YU Rongjun, YANG Rui, HU Yun, ZHAO Shouzhi. Direct 3D Method of Characteristic Based on Efficient Characteristic Tracing and On-the-fly Generation[J]. Atomic Energy Science and Technology, 2025, 59(1): 118-126. DOI: 10.7538/yzk.2024.youxian.0115

基于高效特征线追踪和实时生成的直接三维特征线计算方法

基金项目: 中核集团“青年英才”启明星项目
详细信息
    通讯作者:

    赵守智

  • 中图分类号: TL329

Direct 3D Method of Characteristic Based on Efficient Characteristic Tracing and On-the-fly Generation

  • 摘要:

    直接三维特征线方法(3D MOC)具有强大的几何处理能力,但是其对计算机的性能要求苛刻,尤其是极大的存储规模会限制机器适用性。为此,本文提出了一种基于层次建模树和改进R函数方法的特征线追踪方法,并在此基础上采用特征线信息实时生成(on-the-fly generation)技术,可有效降低计算时间和存储空间。对Takeda和C5G7算例的数值结果表明,在特征值等问题计算上与参考解相比具有良好的计算精度,提出的高效特征线追踪和实时生成方法能够显著降低存储空间,在分别付出增加约45%、23%的时间成本后,两个算例的存储空间分别最大可降低原来的99.51%、99.92%,证明了该方法可大大提高直接3D MOC的机器适用性。

     

    Abstract:

    High-precision full reactor simulation is extremely important for reactor design and safety analysis. In recent years, three-dimensional (3D) full reactor transport calculation becomes one of the key development directions of deterministic calculation in reactor physics. 3D full reactor transport calculation requires the calculation program to have strong anisotropy, fine geometric processing and large-scale parallel computing capabilities. The method of characteristics (MOC) is one of the excellent implementations of the next generation of neutronic deterministic calculations due to its powerful geometric processing capabilities and natural parallelism. However, the huge requirements of the MOC for computing time and space storage are also important factors restricting development. The main technical routes of MOC include 2D/1D method and direct 3D method. Among them, the direct 3D MOC has powerful geometric processing capabilities, but its performance requirements for computers are harsh, especially the large storage scale will limit the applicability of the machine. To this end, a characteristic tracking method based on hierarchical modeling tree and improved R-function method was proposed, and on this basis, on-the-fly generation of characteristic information technology was adopted, which could effectively reduce computation time and spatial storage. The numerical results of the Takeda and C5G7 benchmarks show that the method have good computational accuracy compared to the reference solution in the calculation of eigenvalues and other benchmarks. The proposed efficient characteristic tracking and on-the-fly generation technology can significantly reduce spatial storage. After increasing the time cost by about 45% and 23% respectively, the storage of the two benchmarks can be reduced by up to 99.51% and 99.92%, respectively. It is proved that this method can greatly improve the machine applicability of the direct 3D MOC.

     

  • 图  1   rU函数图像

    Figure  1.   rU function image

    图  2   实时生成计算流程图

    Figure  2.   Flowchart of on-the-fly generation method

    图  3   Takeda算例

    Figure  3.   Takeda benchmark

    图  4   C5G7算例的堆芯径向排布和边界条件

    UO2—— UO2组件;MOX——MOX组件

    Figure  4.   Radial core arrangement and boundary condition of C5G7 benchmark

    图  5   C5G7算例栅元归一化功率的分布

    Figure  5.   Pin normalized power distribution of C5G7 benchmark

    图  6   C5G7算例的相对误差

    Figure  6.   Relative error distribution of C5G7 benchmark

    表  1   Takeda算例的keff计算结果

    Table  1   keff result of Takeda benchmark

    对比数据 case1 case2
    keff 误差/pcm keff 误差/pcm
    参考解 0.977 80(±60 pcm) 0.962 40(±60 pcm)
    计算值 0.977 60 −20 0.962 57 17
    下载: 导出CSV

    表  2   Takeda case1算例各区域的平均通量

    Table  2   Average flux of each region in Takeda case1 benchmark

    对比数据 能群 空区 燃料区 反射层
    平均通量 相对误差/% 平均通量 相对误差/% 平均通量 相对误差/%
    参考解 1 1.450 0×10−3 4.750 9×10−3 5.9251×10−4
    2 9.740 6×10−4 8.6998×10−4 9.1404×10−4
    计算值 1 1.464 3×10−3 0.98 4.7241×10−3 −0.56 5.9489×10−4 0.40
    2 9.537 2×10−4 −2.09 8.7009×10−4 0.01 8.8284×10−4 −3.41
    下载: 导出CSV

    表  3   Takeda case2算例各区域的平均通量

    Table  3   Average flux of each region in Takeda case2 benchmark

    对比数据 能群 控制棒 燃料区 反射层
    平均通量 相对误差/% 平均通量 相对误差/% 平均通量 相对误差/%
    参考解 1 1.224 7×10−3 4.9125×10−3 5.910 9×10−4
    2 2.460 4×10−4 8.6921×10−4 8.7897×10−4
    计算值 1 1.224 9×10−3 0.01 4.8794×10−3 −0.67 5.9305×10−4 0.33
    2 2.514 9×10−4 2.21 8.6934×10−4 0.02 8.5151×10−4 −3.12
    下载: 导出CSV

    表  4   改变方位角的测试结果

    Table  4   Test result of changed azimuth angle

    极角
    数量
    方位角
    数量
    保存分段信息的
    存储空间/MB
    实时生成的
    存储空间/MB
    减少的空间
    存储比/%
    保存分段信息的
    计算时间/s
    实时生成的
    计算时间/s
    时间增加
    比例/%
    keff计算偏差/
    pcm
    4 4 698.5 35.1 94.97 1 361.52 1 940.254 42.51 −166
    4 8 1 218.2 35.1 97.12 3 145.536 4 590.991 45.95 −206
    4 16 2 336.8 35.1 98.50 5 414.808 7 800.432 44.06 −199
    4 32 4 589.6 35.1 99.24 11 571.95 16 628.77 43.70 −205
    下载: 导出CSV

    表  5   改变密度测试结果

    Table  5   Test result of changed density

    (边长/特征线密度)/cm−1 保存分段信息的
    存储空间/MB
    实时生成的
    存储空间/MB
    减少的空间
    存储比/%
    保存分段信息的
    计算时间/s
    实时生成的
    计算时间/s
    时间增加
    比例/%
    keff计算偏差/
    pcm
    5 339 35.1 89.65 745.362 1 066.985 7 43.15 −204
    10 1 218.2 35.1 97.12 3 145.536 4 584.304 2 45.74 −206
    15 2 661.1 35.1 98.68 5 816.61 8 399.184 8 44.40 −207
    20 4 665.3 35.1 99.25 11 759.769 16 876.444 43.51 −207
    25 7 223.9 35.1 99.51 18 146.043 25841.78 42.41 −207
    下载: 导出CSV

    表  6   C5G7算例计算结果汇总

    Table  6   Summary of calculation result for C5G7 benchmark

    参数 参数值
    算例参考解(MCNP) MOC
    keff 1.186 55 1.185 19
    keff误差,pcm −136
    AVG,% 0.32 0.435
    RMS,% 0.34 0.596
    MRE,% 0.27 0.376
    最大栅元归一化功率 2.498 2.505
    最大栅元功率相对误差,% ±0.16 0.28
    最小栅元归一化功率 0.232 0.236
    最小栅元功率相对误差,% ±0.58 1.72
    栅元最大相对误差,% 3.90
    下载: 导出CSV

    表  7   2个极角和4个方位角下的计算时间和存储空间

    Table  7   Calculation time and storage requirement of 2 polar and 4 azimuth angles

    求积组 保存分段信息的
    存储空间/MB
    实时生成的
    存储空间/MB
    减少的空间
    存储比/%
    保存分段信息的
    计算时间/s
    实时生成的
    计算时间/s
    时间增加
    比例/%
    keff计算偏差/
    pcm
    (2,4)
    (8,16)
    6 940
    111 040
    84
    84
    98.79
    99.92
    82 576
    1 215 955
    101 386
    1 492 950
    22.78
    22.78
    −434
    −136
    下载: 导出CSV
  • [1]

    CHO J Y, JOO H G, KIM K S, et al. Cell based CMFD formulation for acceleration of whole-core method of characteristics calculations[J]. Journal of the Korean Nuclear Society, 2002, 34(3): 250-258.

    [2]

    YUK S, CHO N Z. Whole-core transport solutions with 2-D/1-D fusion kernel via p-CMFD acceleration and p-CMFD embedding of nonoverlapping local/global iterations[J]. Nuclear Science and Engineering, 2015, 181(1): 1-16. doi: 10.13182/NSE14-88

    [3]

    WU W, YU Y, LI Q, et al. Verification of the 2D/1D coupling 3D transport code TIGER with C5G7 benchmarks[J]. International Journal of Advanced Nuclear Reactor Design and Technology, 2019, 1: 43-50. doi: 10.1016/j.jandt.2019.10.001

    [4]

    KELLEY B W, LARSEN E W. A consistent 2D/1D approximation to the 3D neutron transport equation[J]. Nuclear Engineering and Design, 2015, 295: 598-614. doi: 10.1016/j.nucengdes.2015.07.026

    [5]

    CHO J Y, KIM K S, LEE C C, et al. Axial SPN and radial MOC coupled whole core transport calculation[J]. Journal of Nuclear Science and Technology, 2007, 44(9): 1156-1171. doi: 10.1080/18811248.2007.9711359

    [6]

    JOO H G, CHO J Y, KIM K S, et al. Methods and performance of a three-dimensional whole-core transport code DeCART[C]. Chicago: [s. n.], 2004.

    [7]

    BAN Y S, MASIELLO E, LENAIN R, et al. Code-to-code comparisons on spatial solution capabilities and performances between nTRACER and the standalone IDT solver of APOLLO3®[J]. Annals of Nuclear Energy, 2018, 115: 573-594. doi: 10.1016/j.anucene.2018.02.011

    [8]

    CHEN J, LIU Z, ZHAO C, et al. A new high-fidelity neutronics code NECP-X[J]. Annals of Nuclear Energy, 2018, 116: 417-428. doi: 10.1016/j.anucene.2018.02.049

    [9]

    LIU Z, WU H, CAO L, et al. A new three-dimensional method of characteristics for the neutron transport calculation[J]. Annals of Nuclear Energy, 2011, 38(2/3): 447-454.

    [10]

    KOCHUNAS B M. A hybrid parallel algorithm for the 3-D method of characteristics solution of the Boltzmann transport equation on high performance compute clusters[R]. [S. l.]: [s. n.], 2013.

    [11]

    SANTANDREA S, GRAZIANO L, SCIANNANDRONE D. Accelerated polynomial axial expansions for full 3D neutron transport MOC in the APOLLO3® code system as applied to the ASTRID fast breeder reactor[J]. Annals of Nuclear Energy, 2018, 113: 194-236. doi: 10.1016/j.anucene.2017.11.010

    [12]

    SCIANNANDRONE D, SANTANDREA S, SANCHEZ R. Optimized tracking strategies for step MOC calculations in extruded 3D axial geometries[J]. Annals of Nuclear Energy, 2016, 87: 49-60. doi: 10.1016/j.anucene.2015.05.014

    [13]

    GRAZIANO L, SANTANDREA S, SCIANNANDRONE D. Polynomial axial expansion in the method of characteristics for neutron transport in 3D extruded geometries[J]. EPJ Web of Conferences, 2017, 153: 06027. doi: 10.1051/epjconf/201715306027

    [14]

    WU W, YU Y, LUO Q, et al. Verification of the 3D capability of OpenMOC with the C5G7 3D extension benchmark[J]. Annals of Nuclear Energy, 2020, 140: 107293. doi: 10.1016/j.anucene.2019.107293

    [15]

    OpenMOC Development Team. OpenMOC: A method of characteristic code for nuclear reactor physics calculations[EB/OL]. https://github.com/mit-crpg/OpenMOC.

    [16] 杨睿, 胡赟, 单浩栋, 等. 适用于任意几何的特征线边界条件处理方法[J]. 原子能科学技术, 2020, 54(11): 2145-2152. doi: 10.7538/yzk.2019.youxian.0849

    YANG Rui, HU Yun, SHAN Haodong, et al. Boundary condition processing method for MOC calculation in arbitrary geometry[J]. Atomic Energy Science and Technology, 2020, 54(11): 2145-2152 (in Chinese). doi: 10.7538/yzk.2019.youxian.0849

    [17]

    TAKEDA T, IKEDA H. 3-D neutron transport benchmarks[J]. Journal of Nuclear Science and Technology, 1991, 28(7): 656-669. doi: 10.1080/18811248.1991.9731408

    [18]

    SMITH M, LEWIS E, NA B. Benchmark on deterministic 3-D MOX fuel assembly transport calculations without spatial homogenization[J]. Progress in Nuclear Energy, 2004, 45: 383-393.

图(6)  /  表(7)
计量
  • 文章访问数:  47
  • HTML全文浏览量:  8
  • PDF下载量:  13
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-03-11
  • 修回日期:  2024-06-09
  • 网络出版日期:  2024-10-15
  • 刊出日期:  2025-01-19

目录

    /

    返回文章
    返回