“Virtual density” and traditional boundary perturbation theories: analytic equivalence and numeric comparison
Name
Forget_Virtual density and traditional.pdf
Description
Accepted version
Size
517.68 KB
Format
Adobe PDF
Checksum (MD5)
52c27b3e51aaf3240ac413c2ee0a271b
Author(s)
Reed, Mark
Smith, Kord
Forget, Benoit
Alternative Title
Virtual density and traditional boundary perturbation theories: analytic equivalence and numeric comparison
Date Issued
February 2018
Journal
Annals of nuclear energy
Publisher
Elsevier BV
Citation
Reed, Mark, Kord Smith, and Benoit Forget, "'Virtual density' and traditional boundary perturbation theories: analytic equivalence and numeric comparison." Annals of nuclear energy 112 (2018): p. 531-48 doi 10.1016/J.ANUCENE.2017.09.028 ©2018 Author(s)
Version
Author's final manuscript
Abstract
We compare and contrast “virtual density” perturbation theory with the traditional boundary perturbation theory developed by Pomraning, Larsen, and Rahnema in the context of diffusion theory. First, after reviewing that literature, we mathematically prove that virtual density perturbations and traditional boundary perturbations are precisely equivalent for arbitrary 1-D problems, which constitute non-uniform isotropic expansions. We also mathematically prove that these two perturbation theories are equivalent for 2-D boundary shift problems, which constitute non-uniform anisotropic expansions. Extension of this proof to swellings or 3-D problems is straightforward. We compare the two theories numerically for a series of alternating uranium and sodium 1-D slabs in finite difference diffusion, and we show that virtual density theory predicts reactivities much more accurately and efficiently than traditional boundary perturbation theory. Boundary perturbation theory is often very inaccurate on a coarse mesh but converges to the virtual density solution as the mesh becomes finer. We also compare the two theories for axial assembly swelling in an abbreviated FFTF benchmark with a coarse mesh. Here we find that reactivity coefficients obtained via virtual density perturbation theory agree with reference solutions to within 0.1%, while those obtained via boundary perturbation theory exhibit sporadic accuracy – sometimes in the range of 1–5% error, more frequently in the range 5–20% error, and occasionally well over 100% error in control rod assemblies. We conclude that although virtual density perturbation theory and boundary perturbation theory are analytically equivalent, boundary perturbations in diffusion theory are often thwarted in coarse mesh finite difference solutions due to inaccurate flux gradients along mesh cell surfaces in heterogeneous cores. ©2018
MIT Department
Massachusetts Institute of Technology. Department of Nuclear Science and Engineering
Terms of Use
Creative Commons Attribution-NonCommercial-NoDerivs License
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.1016/J.ANUCENE.2017.09.028
