Three-dimensional h-adaptivity for the multigroup neutron diffusion equations

Yaqi Wang, Wolfgang Bangerth, Jean Ragusa

Research output: Contribution to journalArticlepeer-review

53 Scopus citations


Adaptive mesh refinement (AMR) has been shown to allow solving partial differential equations to significantly higher accuracy at reduced numerical cost. This paper presents a state-of-the-art AMR algorithm applied to the multigroup neutron diffusion equation for reactor applications. In order to follow the physics closely, energy group-dependent meshes are employed. We present a novel algorithm for assembling the terms coupling shape functions from different meshes and show how it can be made efficient by deriving all meshes from a common coarse mesh by hierarchic refinement. Our methods are formulated using conforming finite elements of any order, for any number of energy groups. The spatial error distribution is assessed with a generalization of an error estimator originally derived for the Poisson equation. Our implementation of this algorithm is based on the widely used Open Source adaptive finite element library deal.II and is made available as part of this library's extensively documented tutorial. We illustrate our methods with results for 2-D and 3-D reactor simulations using 2 and 7 energy groups, and using conforming finite elements of polynomial degree up to 6. © 2008 Elsevier Ltd. All rights reserved.
Original languageEnglish (US)
Pages (from-to)543-555
Number of pages13
JournalProgress in Nuclear Energy
Issue number3
StatePublished - Apr 2009
Externally publishedYes

Bibliographical note

KAUST Repository Item: Exported on 2020-10-01
Acknowledged KAUST grant number(s): KUS-C1-016-04
Acknowledgements: Part of this research was funded through DOE grants DE-FG07-071D14767 and DE-FG07-051D74692. This publication is also partly based on work supported by Award No. KUS-C1-016-04, made by the King Abdullah University of Science and Technology (KAUST).
This publication acknowledges KAUST support, but has no KAUST affiliated authors.


Dive into the research topics of 'Three-dimensional h-adaptivity for the multigroup neutron diffusion equations'. Together they form a unique fingerprint.

Cite this