This work illustrates the innovative
results obtained by applying the recently developed the 2nd-order
predictive modeling methodology called “2nd-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 2nd-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 2nd-order
response sensitivities must always be included to quantify the need for
including (or not) the 3rd- and/or 4th-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 1st- and 2nd-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 3rd- and even 4th-order
sensitivities are necessary to ensure consistency between the computed and
measured response. In such cases, the use of only the 1st-order
sensitivities erroneously indicates that the computed results are inconsistent
with the respective measured response. Ongoing research
aims at extending the 2nd-BERRU-PM methodology to fourth-order, thus
enabling the computation of third-order response
References
[1]
Cacuci, D.G. (2023) Second-Order MaxEnt Predictive Modeling Methodology. I: Deterministically Incorporated Computational Model Responses (2nd-BERRU-PMD). American Journal of Computational Mathematics.
[2]
Cacuci, D.G. (2023) Second-Order MaxEnt Predictive Modelling Methodology. II: Probabilistically Incorporated Computational Model (2nd-BERRU-PMP). American Journal of Computational Mathematics.
[3]
Jaynes, E.T. (1957) Information Theory and Statistical Mechanics. Physical Review Journals Archive, 106, 620-630. https://doi.org/10.1103/PhysRev.106.620
[4]
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
[5]
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
[6]
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.
[7]
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.
[8]
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.
[9]
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.
[10]
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
[11]
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
[12]
Cacuci, D.G. and Fang, R. (2021) Fourth-Order Adjoint Sensitivity Analysis of an OECD/NEA Reactor Physics Benchmark: I. Mathematical Expressions and CPU-Time Comparisons for Computing 1st-, 2nd-and 3rd-Order Sensitivities. American Journal of Computational Mathematics, 11, 94-132. https://doi.org/10.4236/ajcm.2021.112009
[13]
Cacuci, D.G. and Fang, R. (2021) Fourth-Order Adjoint Sensitivity Analysis of an OECD/NEA Reactor Physics Benchmark: II. Mathematical Expressions and CPU-Time Comparisons for Computing 4th-Order Sensitivities. American Journal of Computational Mathematics, 11, 133-156. https://doi.org/10.4236/ajcm.2021.112010