DenovoA New ThreeDimensional Parallel Discrete Ordinates Code in SCALE
Abstract
Denovo is a new, threedimensional, discrete ordinates (SN) transport code that uses stateoftheart solution methods to obtain accurate solutions to the Boltzmann transport equation. Denovo uses the KochBakerAlcouffe parallel sweep algorithm to obtain high parallel efficiency on O(100) processors on XYZ orthogonal meshes. As opposed to traditional SN codes that use source iteration, Denovo uses nonstationary Krylov methods to solve the withingroup equations. Krylov methods are far more efficient than stationary schemes. Additionally, classic acceleration schemes (diffusion synthetic acceleration) do not suffer stability problems when used as a preconditioner to a Krylov solver. Denovo's generic programming framework allows multiple spatial discretization schemes and solution methodologies. Denovo currently provides diamonddifference, thetaweighted diamonddifference, lineardiscontinuous finite element, trilineardiscontinuous finite element, and step characteristics spatial differencing schemes. Also, users have the option of running traditional source iteration instead of Krylov iteration. Multigroup upscatter problems can be solved using GaussSeidel iteration with transport, twogrid acceleration. A parallel firstcollision source is also available. Denovo solutions to the Kobayashi benchmarks are in excellent agreement with published results. Parallel performance shows excellent weak scaling up to 20000 cores and good scaling up to 40000 cores.
 Authors:

 ORNL
 Oak Ridge Institute for Science and Education (ORISE)
 Publication Date:
 Research Org.:
 Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
 Sponsoring Org.:
 USDOE
 OSTI Identifier:
 984374
 DOE Contract Number:
 DEAC0500OR22725
 Resource Type:
 Journal Article
 Journal Name:
 Nuclear Technology
 Additional Journal Information:
 Journal Volume: 171; Journal Issue: 8; Journal ID: ISSN 00295450
 Country of Publication:
 United States
 Language:
 English
 Subject:
 73 NUCLEAR PHYSICS AND RADIATION PHYSICS; 97 MATHEMATICAL METHODS AND COMPUTING; D CODES; THREEDIMENSIONAL CALCULATIONS; RADIATION TRANSPORT; BOLTZMANN EQUATION; DISCRETE ORDINATE METHOD; PARALLEL PROCESSING; ALGORITHMS; PERFORMANCE
Citation Formats
Evans, Thomas M, Stafford, Alissa, and Clarno, Kevin T. DenovoA New ThreeDimensional Parallel Discrete Ordinates Code in SCALE. United States: N. p., 2010.
Web.
Evans, Thomas M, Stafford, Alissa, & Clarno, Kevin T. DenovoA New ThreeDimensional Parallel Discrete Ordinates Code in SCALE. United States.
Evans, Thomas M, Stafford, Alissa, and Clarno, Kevin T. 2010.
"DenovoA New ThreeDimensional Parallel Discrete Ordinates Code in SCALE". United States.
@article{osti_984374,
title = {DenovoA New ThreeDimensional Parallel Discrete Ordinates Code in SCALE},
author = {Evans, Thomas M and Stafford, Alissa and Clarno, Kevin T},
abstractNote = {Denovo is a new, threedimensional, discrete ordinates (SN) transport code that uses stateoftheart solution methods to obtain accurate solutions to the Boltzmann transport equation. Denovo uses the KochBakerAlcouffe parallel sweep algorithm to obtain high parallel efficiency on O(100) processors on XYZ orthogonal meshes. As opposed to traditional SN codes that use source iteration, Denovo uses nonstationary Krylov methods to solve the withingroup equations. Krylov methods are far more efficient than stationary schemes. Additionally, classic acceleration schemes (diffusion synthetic acceleration) do not suffer stability problems when used as a preconditioner to a Krylov solver. Denovo's generic programming framework allows multiple spatial discretization schemes and solution methodologies. Denovo currently provides diamonddifference, thetaweighted diamonddifference, lineardiscontinuous finite element, trilineardiscontinuous finite element, and step characteristics spatial differencing schemes. Also, users have the option of running traditional source iteration instead of Krylov iteration. Multigroup upscatter problems can be solved using GaussSeidel iteration with transport, twogrid acceleration. A parallel firstcollision source is also available. Denovo solutions to the Kobayashi benchmarks are in excellent agreement with published results. Parallel performance shows excellent weak scaling up to 20000 cores and good scaling up to 40000 cores.},
doi = {},
url = {https://www.osti.gov/biblio/984374},
journal = {Nuclear Technology},
issn = {00295450},
number = 8,
volume = 171,
place = {United States},
year = {2010},
month = {1}
}