A multithreaded parallel upwind sweep algorithm for the S_(N) transport equations discretized with discontinuous finite elements  

作　　者：Zhi‑Wei Zong Mao‑Song Cheng Ying‑Chi Yu Zhi‑Min Dai 

机构地区：[1]Shanghai Institute of Applied Physics,Chinese Academy of Sciences,Shanghai 201800,China [2]University of Chinese Academy of Sciences,Beijing 100049,China

出　　处：《Nuclear Science and Techniques》2023年第12期229-241,共13页核技术（英文）

摘　　要：The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.

关 键 词：Shielding calculation Discrete ordinates method Discontinuous Galerkin finite element method Unstructured meshes 

分 类 号：TL329[核科学技术—核技术及应用]