Third Order Adjoint Sensitivity and Uncertainty Analysis of an OECD/NEA Reactor Physics Benchmark: III. Response Moments

The (180)3 third-order mixed sensitivities of the leakage response of a polyethylene-reflected plutonium (PERP) experimental benchmark with respect to the benchmark’s 180 microscopic total cross sections have been computed in accompanying works [1] [2]. This work quantifies the contributions of these (180)3 third-order mixed sensitivities to the PERP benchmark’s leakage response distribution moments (expected value, variance and skewness) and compares these contributions to those stemming from the corresponding first- and second-order sensitivities of the PERP benchmark’s leakage response with respect to the total cross sections. The numerical results obtained in this work reveal that the importance of the 3rd-order sensitivities can surpass the importance of the 1st- and 2nd-order sensitivities when the parameters’ uncertainties increase. In particular, for a uniform standard deviation of 10% of the microscopic total cross sections, the 3rd-order sensitivities contribute 80% to the response variance, whereas the contribution stemming from the 1st- and 2nd-order sensitivities amount only to 2% and 18%, respectively. Consequently, neglecting the 3rd-order sensitivities could cause a very large non-conservative error by under-reporting the response variance by a factor of 506%. The results obtained in this work also indicate that the effects of the 3rd-order sensitivities are to reduce the response’s skewness in parameter space, rendering the distribution of the leakage response more symmetric about its expected value. The results obtained in this work are the first such results ever published in reactor physics. Since correlations among the group-averaged microscopic total cross sections are not available, only the effects of typical standard deviations for these cross sections could be considered. Due to this lack of correlations among the cross sections, the effects of the mixed 3rd-order sensitivities could not be quantified exactly at this time. These effects could be quantified only when correlations among the group-averaged microscopic total cross sections would be obtained experimentally by the nuclear physics community.

Third-Order Adjoint Sensitivity Analysis of an OECD/NEA Reactor Physics Benchmark: I. Mathematical Framework

This work extends to third-order previously published work on developing the adjoint sensitivity and uncertainty analysis of the numerical model of a polyethylene-reflected plutonium (acronym: PERP) OECD/NEA reactor physics benchmark. The PERP benchmark comprises 21,976 imprecisely known (uncertain) model parameters. Previous works have used the adjoint sensitivity analysis methodology to compute exactly and efficiently all of the 21,976 first-order and (21,976)2 second-order sensitivities of the PERP benchmark’s leakage response to all of the benchmark’s uncertain parameters, showing that the largest and most consequential 1st- and 2nd-order response sensitivities are with respect to the total microscopic cross sections. These results have motivated extending the previous adjoint-based derivations to third-order, leading to the derivation, in this work, of the exact mathematical expressions of the (180)3 third-order sensitivities of the PERP leakage response with respect to these total microscopic cross sections. The formulas derived in this work are valid not only for the PERP benchmark but can also be used for computing the 3rd-order sensitivities of the leakage response of any nuclear system involving fissionable material and internal or external neutron sources. Subsequent works will use the adjoint-based mathematical expressions obtained in this work to compute exactly and efficiently the numerical values of these (180)3 third-order sensitivities (which turned out to be very large and consequential) and use them for a third-order uncertainty analysis of the PERP benchmark’s leakage response.

Third Order Adjoint Sensitivity and Uncertainty Analysis of an OECD/NEA Reactor Physics Benchmark: II. Computed Sensitivities

This work presents the results of the exact computation of (180)3 = 5,832,000 third-order mixed sensitivities of the leakage response of a polyethylene-reflected plutonium (PERP) experimental benchmark with respect to the benchmark’s 180 microscopic total cross sections. This computation was made possible by applying the Third-Order Adjoint Sensitivity Analysis Methodology developed by Cacuci. The numerical results obtained in this work revealed that many of the 3rd-order sensitivities are significantly larger than their corresponding 1st- and 2nd-order ones, which is contrary to the widely held belief that higher-order sensitivities are all much smaller and hence less important than the first-order ones, for reactor physics systems. In particular, the largest 3rd-order relative sensitivity is the mixed sensitivity of the PERP leakage response with respect to the lowest energy-group (30) total cross sections of 1H (“isotope 6”) and 239Pu (“isotope 1”). These two isotopes are shown in this work to be the two most important parameters affecting the PERP benchmark’s leakage response. By comparison, the largest 1st-order sensitivity is that of the PERP leakage response with respect to the lowest energy-group total cross section of isotope 1H, having the value , while the largest 2nd-order sensitivity is . The 3rd-order sensitivity analysis presented in this work is the first ever such analysis in the field of reactor physics. The consequences of the results presented in this work on the uncertainty analysis of the PERP benchmark’s leakage response will be presented in a subsequent work.

Fourth-Order Adjoint Sensitivity Analysis of an OECD/NEA Reactor Physics Benchmark: II. Mathematical Expressions and CPU-Time Comparisons for Computing 4^th-Order Sensitivities

This work extends to fourth-order previously published work on developing the adjoint sensitivity and uncertainty analysis of the numerical model of a polyethylene-reflected plutonium (acronym: PERP) OECD/NEA reactor physics benchmark. Previous works showed that the third-order sensitivities of the PERP leakage response with respect to these total microscopic cross sections are far larger than the corresponding 1st-order and 2nd-order ones, thereby having the largest impact on the uncertainties induced in the PERP benchmark’s response. This finding has motivated the development of the original 4th-order formulas presented in this work, which are valid not only for the PERP benchmark but can also be used for computing the 4th-order sensitivities of response of any nuclear system involving fissionable material and internal or external neutron sources. Subsequent works will use the adjoint-based mathematical expressions obtained in this work to compute exactly and efficiently the numerical values of the largest fourth-order sensitivities of the PERP benchmark’s response to the total microscopic cross sections, and use them for a pioneering fourth-order uncertainty analysis of the PERP benchmark’s response.