The Development and Implementation of a One-Dimensional Sn Method in the 2D-1D Integral Transport Solution

被引:11
|
作者
Hursin, M. [1 ]
Collins, B. [1 ]
Xu, Y. [1 ]
Downar, T. [1 ]
机构
[1] Univ Michigan, Dept Nucl Engn & Radiol Sci, Ann Arbor, MI 48109 USA
关键词
BENCHMARK;
D O I
10.13182/NSE12-4
中图分类号
TL [原子能技术]; O571 [原子核物理学];
学科分类号
0827 ; 082701 ;
摘要
During the last several years, a class of algorithms has been developed based on two-dimensional-one-dimensional (2D-1D) decomposition of the reactor transport problem. The current 2D-1D algorithm implemented in the DeCART (Deterministic Core Analysis based on Ray Tracing) code solves a set of coupled 2D planar transport and 1D axial diffusion equations. This method has been successfully applied to several light water reactor analysis problems. However, applications with strong axial heterogeneities have exposed the limitations of the current diffusion solvers used for the axial solution. The work reported in this paper is the implementation of a discrete ordinates (S-n)-based axial solver in DeCART. An S-n solver is chosen to preserve the consistency of the angular discretization between the radial method of characteristics and axial solvers. This paper presents the derivation of the nodal expansion method (NEM)-S S-n equations and its implementation in DeCART. The subplane spatial refinement method is introduced to reduce the computational cost and improve the accuracy of the calculations. The NEM-S-n axial solver is tested using the C5G7 benchmark. The DeCART results with the axial diffusion solver shows k(eff) errors of approximately -95, -74, and -110 pcm for the unrodded configuration, rodded configuration A, and rodded configuration B, respectively. These errors decrease to approximately -40, -11, and -12 pcm by using the NEM-S-n solver. In terms of pin power distribution, the use of the NEM-S-n solver has a small effect, except for the heavily rodded configuration. The implementation of the subplane scheme makes it possible to maintain a coarse axial mesh and therefore to reduce the computational cost of the three-dimensional calculations without reducing the accuracy of the solution.
引用
收藏
页码:186 / 200
页数:15
相关论文
共 50 条
  • [1] Strongly-interacting bosons at 2D-1D dimensional crossover
    Yao, Hepeng
    Pizzino, Lorenzo
    Giamarchi, Thierry
    SCIPOST PHYSICS, 2023, 15 (02):
  • [2] Elastodynamic 2D-1D coupling using the DtN method
    Rabinovich, Daniel
    Givoli, Dan
    JOURNAL OF COMPUTATIONAL PHYSICS, 2022, 448
  • [3] 2D-1D Dimensional reduction in a toy model for magnetoelastic interactions
    Mouhcine Tilioua
    Applications of Mathematics, 2011, 56
  • [5] 2D-1D mixed-dimensional heterostructures: progress, device applications and perspectives
    Huang, Pei-Yu
    Qin, Jing-Kai
    Zhu, Cheng-Yi
    Zhen, Liang
    Xu, Cheng-Yan
    JOURNAL OF PHYSICS-CONDENSED MATTER, 2021, 33 (49)
  • [6] One-Dimensional Integral Imaging 3D Display Systems
    Hirayama, Yuzo
    IDW'11: PROCEEDINGS OF THE 18TH INTERNATIONAL DISPLAY WORKSHOPS, VOLS 1-3, 2011, : 531 - 534
  • [7] A modified boundary integral solution of coupled linear and nonlinear one-dimensional transport equations
    Onyejekwe, OO
    BOUNDARY ELEMENTS XXII, 2000, 8 : 471 - 480
  • [8] One-dimensional (1-D) nanoscale heterostructures
    Shen, Guozhen
    Chen, Di
    Bando, Yoshio
    Golberg, Dmitri
    JOURNAL OF MATERIALS SCIENCE & TECHNOLOGY, 2008, 24 (04) : 541 - 549
  • [9] One-Dimensional(1-D) Nanoscale Heterostructures
    Yoshio BANDO
    Dmitri GOLBERG
    JournalofMaterialsScience&Technology, 2008, (04) : 541 - 549
  • [10] Vectorization of a 2D-1D Iterative Algorithm for the 3D Neutron Transport Problem in Prismatic Geometries
    Moustafa, Salli
    Fevotte, Francois
    Lathuiliere, Bruno
    Plagne, Laurent
    SNA + MC 2013 - JOINT INTERNATIONAL CONFERENCE ON SUPERCOMPUTING IN NUCLEAR APPLICATIONS + MONTE CARLO, 2014,