Massively parallel, three-dimensional transport solutions for the k-eigenvalue problem

Gregory G. Davidson, Thomas M. Evans, Joshua J. Jarrell, Steven P. Hamilton, Tara M. Pandya, Rachel N. Slaybaugh

Research output: Contribution to journalArticlepeer-review

25 Scopus citations

Abstract

We have implemented a new multilevel parallel decomposition in the Denovo discrete ordinates radiation transport code. In concert with Krylov subspace iterative solvers, the multilevel decomposition allows concurrency over energy in addition to space-angle, enabling scalability beyond the limits imposed by the traditional Koch-Baker-Alcouffe (KBA) space-angle partitioning. Furthermore, a new Arnoldi-based k-eigenvalue solver has been implemented. The added phase-space concurrency combined with the high-performance Krylov and Arnoldi solvers has enabled weak scaling to 0(105) cores on the Titan XK7 supercomputer. The multilevel decomposition provides a mechanism for scaling to exascale computing and beyond.

Original languageEnglish
Pages (from-to)111-125
Number of pages15
JournalNuclear Science and Engineering
Volume177
Issue number2
DOIs
StatePublished - Jun 2014

Fingerprint

Dive into the research topics of 'Massively parallel, three-dimensional transport solutions for the k-eigenvalue problem'. Together they form a unique fingerprint.

Cite this