Solving the Neutron Transport Equation for Microreactor Modeling Using Unstructured Meshes and Exascale Computing Architectures

The Microreactor Exascale eZ CALculation (MEZCAL) tool has been developed to accurately and efficiently solve the neutron transport equation in general, unstructured meshes to support the design and modeling of microreactors. MEZCAL solves the self-adjoint angular flux form of the neutron transport equation using the finite element method. As the neutron transport equation is computationally expensive to solve, MEZCAL is designed to efficiently use exascale computing architectures, with an emphasis on graphics processing unit computing. To leverage existing tools, MEZCAL is built using the MFEM library and uses solvers from HYPRE, PETSc, and SLEPc. Verification of the neutron transport solver in MEZCAL is demonstrated with the solution to a one-dimensional cylindrical problem that has a semi-analytic solution. After verification, a realistic microreactor based on the MARVEL microreactor design is modeled using MEZCAL. Spatial and angular refinement results are presented for a two-dimensional model of the MARVEL microreactor, and the eigenvalue is converged to approximately 60 pcm. This convergence required a very fine mesh and more than 3.76 Billion Degrees Of Freedom (BDOF). Preliminary results are also presented for a three-dimensional model of the MARVEL microreactor. Finally, a weak scaling study is performed to investigate how the methods in MEZCAL will scale for larger problems with the next generation of exascale computing architectures.