Second-Order MaxEnt Predictive Modelling Methodology. III: Illustrative Application to a Reactor Physics Benchmark
- 1 Center for Nuclear Science and Energy, Department of Mechanical Engineering, University of South Carolina, Columbia, SC, USA
- 2 Center for Nuclear Science and Energy, Department of Mechanical Engineering, University of South Carolina, Columbia, SC, USA
Abstract
This work illustrates the innovative results obtained by applying the recen tly developed the 2 nd -order predictive modeling methodology called “2 nd - BERRU-PM”, where the acronym BERRU denotes “best-estimate results with reduced uncertainties” and “PM” denotes “predictive modeling.” The physical system selected for this illustrative application is a polyethylene-reflected plutonium (acronym: PERP) OECD/NEA reactor physics benchmark. This benchmark is modeled using the neutron transport Boltzmann equation (involving 21,976 uncertain parameters), the solution of which is representative of “large-scale computations.” The results obtained in this work confirm the fact that the 2 nd -BERRU-PM methodology predicts best-estimate results that fall in between the corresponding computed and measured values, while reducing the predicted standard deviations of the predicted results to values smaller than either the experimentally measured or the computed values of the respective standard deviations. The obtained results also indicate that 2 nd -order response sensitivities must always be included to quantify the need for including (or not) the 3 rd - and/or 4 th -order sensitivities. When the parameters are known with high precision, the contributions of the higher-order sensitivities diminish with increasing order, so that the inclusion of the 1 st - and 2 nd -order sensitivities may suffice for obtaining accurate predicted best- estimate response values and best-estimate standard deviations. On the other hand, when the parameters’ standard deviations are sufficiently large to approach (or be outside of) the radius of convergence of the multivariate Taylor-series which represents the response in the phase-space of model parameters, the contributions stemming from the 3 rd - and even 4 th -order sensitivities are necessary to ensure consistency between the computed and measured response. In such cases, the use of only the 1 st -order sensitivities erroneously indicates that the computed results are inconsistent with the respective measured response. Ongoing research aims at extending the 2 nd -BERRU-PM methodology to fourth-order, thus enabling the computation of third-order response correlations (skewness) and fourth-order response correlations (kurtosis).
- Cacuci, D.G. (2023) Second-Order MaxEnt Predictive Modeling Methodology. I: Deterministically Incorporated Computational Model Responses (2nd-BERRU-PMD). American Journal of Computational Mathematics.
- Cacuci, D.G. (2023) Second-Order MaxEnt Predictive Modelling Methodology. II: Probabilistically Incorporated Computational Model (2nd-BERRU-PMP). American Journal of Computational Mathematics.
- Jaynes, E.T. (1957) Information Theory and Statistical Mechanics. Physical Review Journals Archive, 106, 620-630. https://doi.org/10.1103/PhysRev.106.620
- Shannon, C.E. (1948) A Mathematical Theory of Communication. The Bell System Technical Journal, 27, 379-423. https://doi.org/10.1002/j.1538-7305.1948.tb01338.x
- Lewis, J.M., Lakshmivarahan, S. and Dhall, S.K. (2006) Dynamic Data Assimilation: A Least Square Approach. Cambridge University Press, Cambridge. https://doi.org/10.1017/CBO9780511526480
- Valentine, T.E. (2006) Polyethylene-Reflected Plutonium Metal Sphere Subcritical Noise Measurements, SUB-PU-METMIXED-001. International Handbook of Evaluated Criticality Safety Benchmark Experiments, NEA/NSC/DOC(95)03/I-IX, Organization for Economic Co-Operation and Development, Nuclear Energy Agency, Paris.
- Alcouffe, R.E., Baker, R.S., Dahl, J.A., Turner, S.A. and Ward, R. (2008) PARTISN: A Time-Dependent, Parallel Neutral Particle Transport Code System. Los Alamos National Laboratory, Los Alamos.
- Wilson, W.B., Perry, R.T., Shores, E.F., Charlton, W.S., Parish, T.A., Estes, G.P., Brown, T.H., Arthur, E.D., Bozoian, M., England, T.R., et al. (2002) SOURCES 4C: A Code for Calculating (α, n), Spontaneous Fission, and Delayed Neutron Sources and Spectra. Proceedings of the American Nuclear Society/Radiation Protection and Shielding Division 12th Biennial Topical Meeting, Santa Fe, 14-18 April 2002.
- Conlin, J.L., Parsons, D.K., Gardiner, S.J., Beth Lee, M.G.M. and White, M.C. (2013) MENDF71X: Multigroup Neutron Cross-Section Data Tables Based upon ENDF/B-VII.1X. Los Alamos National Laboratory Report LA-UR-5-29571, Los Alamos National Laboratory, Los Alamos.
- Cacuci, D.G. and Fang, R. (2023) The nth-Order Comprehensive Adjoint Sensitivity Analysis Methodology, Volume II: Overcoming the Curse of Dimensionality: Large-Scale Application. Springer, Cham. https://doi.org/10.1007/978-3-031-19635-5
- Fang, R. and Cacuci, D.G. (2022) Fourth-Order Adjoint Sensitivity and Uncertainty Analysis of an OECD/NEA Reactor Physics Benchmark: II. Computed Response Uncertainties. Journal of Nuclear Engineering, 3, 1-16. https://doi.org/10.3390/jne3010001