EPJ Web Conf.
Volume 247, 2021PHYSOR2020 – International Conference on Physics of Reactors: Transition to a Scalable Nuclear Future
|Number of page(s)||8|
|Published online||22 February 2021|
PARAMETRIC STUDY OF PARALLEL BLOCK JACOBI / SOURCE ITERATION HYBRID METHODS IN 2-D CARTESIAN GEOMETRY AND CONSTRUCTION OF THE INTEGRAL TRANSPORT MATRIX METHOD MATRICES VIA GREEN’S FUNCTIONS
NC State University, Department of Nuclear Engineering 2500 Stinson Drive, Raleigh, NC, 27607, United States
Published online: 22 February 2021
Parallel Block Jacobi (PBJ)  is an asynchronous spatial domain decomposition with application in solving the neutron transport equation due to its extendibility to massively parallel solution in unstructured spatial meshes (grids) without the use of the computationally complex and expensive sweeps required by the Source Iteration (SI) method in these applications.  However, PBJ iterative methods suffer a lack of iterative robustness in problems with optically thin cells,  which we have previously demonstrated to be a consequence of PBJ’s asynchronicity. To mitigate this effect, we have developed multiple PBJ / SI hybrid methods which employ a PBJ method (Parallel Block Jacobi - Integral Transport Matrix Method (PBJ-ITMM) or Inexact Parallel Block Jacobi (IPBJ)) along with SI. [3,4] In this work, we perform a parametric study to determine performance of numerous PBJ / SI hybrid methods as a function of multiple problem parameters. This parametric study reached 5 main conclusions: 1) our hybrid approach is more effective with PBJ-ITMM than with IPBJ, 2) for PBJ-ITMM, there is a hybrid method that mitigates the aforementioned iterative slowdown in optically thin cells without diminishing the method’s potential parallelism in unstructured grids, 3) this hybrid method is most effective in problems with large, continuous regions of very thin cells, 4) the best performing hybrid method consistently executes within a factor of ten slower than current state-of-the-art acceleration methods that are not efficiently extendable to the massively parallel regime, and 5) both PBJ-ITMM and IPBJ are observed to be viable approaches for our desired applications. In the pursuit of implementing PBJ-ITMM in unstructured grids, we conclude with a description of the Green’s Function ITMM Construction (GFIC) algorithm, which allows for the ITMM matrices to be constructed using the pre-existing SI sweep algorithm already present in unstructured grid SN transport codes.
Key words: Neutron Transport / Parallel / Parallel Block Jacobi / Integral Transport Matrix Method
© The Authors, published by EDP Sciences, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.