Next Article in Journal / Special Issue
Structure-Dependent Doppler Broadening Using a Generalized Thermal Scattering Law
Previous Article in Journal / Special Issue
Generation of the TSL for Zirconium Hydrides from Ab Initio Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Illustrating Important Effects of Second-Order Sensitivities on Response Uncertainties in Reactor Physics

Department of Mechanical Engineering, University of South Carolina, 300 Main Street, Columbia, SC 29208, USA
J. Nucl. Eng. 2021, 2(2), 114-123; https://0-doi-org.brum.beds.ac.uk/10.3390/jne2020012
Submission received: 16 November 2020 / Revised: 11 March 2021 / Accepted: 25 March 2021 / Published: 8 April 2021
(This article belongs to the Special Issue Selected Papers from PHYSOR 2020)

Abstract

:
This paper illustrates the relative importance of the largest first- and second-order sensitivities of the leakage response of an OECD/NEA reactor physics benchmark (a polyethylene-reflected plutonium sphere) to the benchmark’s underlying total cross sections. It will be shown that numerous 2nd-order sensitivities of the leakage response with respect to the total cross sections are significantly larger than the largest corresponding 1st-order sensitivities. In particular, the contributions of the 2nd-order sensitivities cause the mean (expected) value of the response to differ appreciably from its computed value and also cause the response distribution to be skewed towards positive values relative to the mean. Neglecting these large 2nd-order sensitivities would cause very large non-conservative errors by under-reporting the response’s variance and expected value. The results presented in this paper also underscore the need for obtaining reliable cross section covariance data, which are currently unavailable. Finally, comparing the CPU-times needed for computations, this paper demonstrates that the Second-Order Adjoint Sensitivity Analysis Methodology is the only practical method for computing 2nd-order sensitivities exactly, without introducing methodological errors, for large-scale systems characterized by many uncertain parameters.

1. Introduction

The work reported in this work is based on the PHYSOR-2020 plenary invited conference paper entitled “On the importance of second-order response sensitivities to nuclear data in reactor physics uncertainty analysis.” This paper was selected by the PHYSOR-2020 Technical Committee to be published in the Special Issue of the Journal of Nuclear Engineering. To avoid duplicate publication, the contents and the title of the conference paper have been re-written to make the present work suitable for publication in this journal’s Special Issue.
The OECD Nuclear Energy Agency (OECD/NEA) International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook [1] describes several fundamental subcritical reactor physics benchmarks that use a 4.5 kg alpha-phase plutonium sphere constructed at Los Alamos National Laboratory in 1980 [2]. This sphere was initially used as a neutron source for conducting subcritical experiments aimed at estimating the reactivity worth of beryllium reflectors (and was therefore colloquially called the “BeRP ball”) but was subsequently also used for other subcritical experiments using tungsten, nickel and/or polyethylene reflectors. Miller et al. [3] evaluated computationally neutron multiplicity measurements for the “polyethylene-reflected BeRP ball”, noticing significant disagreement between their computational results and the corresponding measurements.
In addition to the nuclear data investigated in [3], the “polyethylene-reflected plutonium sphere” benchmark contains other uncertain data, the most prominent being the total neutron cross sections. The 1st- and 2nd-order sensitivities of this benchmark’s leakage response to the benchmark’s cross sections have been computed by applying Cacuci’s 2nd-ASAM [4]; the most significant results thus obtained are summarized in this work. Many of the 2nd-order sensitivities were found to be much larger than the corresponding 1st-order sensitivities, shifting the leakage response’s mean value significantly away from its computed value and causing asymmetries in the response distribution. All in all, the effects on the benchmark’s leakage response’s expected value, variance, and skewness of the 2nd-order sensitivities with respect to the total cross sections are much larger than the corresponding effects stemming from the corresponding 1st-order sensitivities.

2. Largest 1st- and 2nd-Order Sensitivities of the Perp Benchmark’s Leakage with Respect to Perp’s Nuclear Data

The “polyethylene-reflected plutonium sphere benchmark” which will henceforth be called using the acronym “PERP” (benchmark) contains two materials, as follows: “Material 1” (core) contains four isotopes: 239Pu, 240Pu, 69Ga, 71Ga; “Material 2” (reflector) contains two isotopes: C and 1H. The neutron flux distribution within the PERP benchmark was computed by using the PARTISN [5] multigroup discrete ordinates particle transport code, which solves the following multi-group approximation of the neutron transport equation with a spontaneous fission source:
B g ( α ) φ g ( r , Ω ) = Q g ( r ) , g = 1 , , G ,
φ g ( r d , Ω ) = 0 , r d S b , Ω n < 0 , g = 1 , , G ,
where r d denotes the radius of the spherical benchmark, and
B g ( α ) φ g ( r , Ω ) Ω · φ g ( r , Ω ) + Σ t g ( r ) φ g ( r , Ω ) g = 1 G 4 π Σ s g g ( r , Ω Ω ) φ g ( r , Ω ) d Ω χ g ( r ) g = 1 G 4 π ( ν Σ ) f g ( r ) φ g ( r , Ω ) d Ω ,
Q g ( r ) k = 1 N f λ k N k , 1 F k S F ν k S F e E g / a k sinh b k E g , g = 1 , , G .
The PARTISN [5] computations of the neutron flux used the MENDF71X [6] 618-group cross section data collapsed to G = 30 energy groups, as well as a P3 Legendre expansion of the scattering cross section, an angular quadrature of S256, and a fine-mesh spacing of 0.005 cm (comprising 759 meshes for the plutonium sphere of radius of 3.794 cm, and 762 meshes for the polyethylene shell of thickness of 3.81 cm). The quantities appearing in Equations (1)–(4) are defined as follows:
  • The quantity φ g ( r , Ω ) is the customary “group-flux” for group g , and is the unknown state-function obtained by solving Equations (1) and (2), where r d is the radius of the PERP sphere while the vector n denotes the outward unit normal vector at each point on the sphere’s outer boundary, denoted as S b .
  • The total cross section Σ t g ( r ) for energy group g , g = 1 , , G , and material m , is computed using the following expression:
    Σ t g ( r ) = m = 1 M = 2 Σ t , m g ( r ) ; Σ t , m g ( r ) = i I N i , m σ t , i g ( r ) = i I N i , m [ σ f , i g ( r ) + σ c , i g ( r ) + g = 1 G σ s , l = 0 , i g g ( r ) ] , m = 1 , 2 ,
    where σ f , i g ( r ) and σ c , i g ( r ) denote, respectively, the tabulated group microscopic fission and neutron capture cross sections for group g , g = 1 , , G . Other nuclear reactions are negligible in the PERP benchmark.
  • The quantity N i , m denotes the atom density of isotope i in material m; i = 1 , , I , m = 1 , , M , where I denotes the total number of isotopes, and M denotes the total number of materials. The computation of N i , m uses the well-known expression:
    N i , m ρ m w i , m N A A i ,
    where ρ m denotes the mass density of material m, m = 1 , , M ; w i , m denotes the weight fraction of isotope i in material m; A i denotes the atomic weight of isotope i , i = 1 , , I ; N A denotes Avogadro’s number. For the PERP benchmark, I = 6 and M = 2 .
  • The quantity Σ s g g ( r , Ω Ω ) represents the scattering transfer cross section from energy group g , g = 1 , , G into energy group g , g = 1 , , G , and is computed in terms of the l -th order Legendre coefficient σ s , l , i g g (of the Legendre-expanded microscopic scattering cross section from energy group g into energy group g for isotope i ), which are tabulated parameters, in the following finite-order expansion:
    Σ s g g ( r , Ω Ω ) = m = 1 M = 2 Σ s , m g g ( r , Ω Ω ) , Σ s , m g g ( r , Ω Ω ) i = 1 I = 6 N i , m l = 0 I S C T = 3 ( 2 l + 1 ) σ s , l , i g g ( r ) P l ( Ω Ω ) , m = 1 , 2 ,
    where I S C T = 3 denotes the order of the respective finite expansion in Legendre polynomial. The expressions in Equations (5) and (7) indicate that the zeroth order (i.e., l = 0 ) scattering cross sections must be considered separately from the higher order (i.e., l 1 ) scattering cross sections, since the former contribute to the total cross sections, while the latter do not.
  • The quantity N f denotes the total number of spontaneous-fission isotopes. The spontaneous-fission isotopes in the PERP benchmark are “isotope 1” (239Pu) and “isotope 2” (240Pu), so N f = 2 , and the spontaneous fission neutron spectrum of 239Pu and 240Pu, respectively, is approximated by a Watts fission spectrum using the evaluated parameters a k and b k . The decay constant for actinide nuclide k is denoted as λ k , and F k S F denotes the fraction of decays that are spontaneous fission (the “spontaneous-fission branching fraction”).
  • PARTISN [5] computes the quantity ( ν Σ f ) g by directly using the quantities ( ν σ ) f , i g , which are provided in nuclear data files for each isotope i , and energy group g , as follows
    ( ν Σ f ) g = m = 1 M = 2 ( ν Σ f ) m g ; ( ν Σ f ) m g = i = 1 I = 6 N i , m ( ν σ f ) i g , m = 1 , 2 .
  • The quantity χ g ( r ) quantifies the fission spectrum in energy group g .
  • The numerical model of the PERP benchmark contains 7477 nonzero parameters which are subject to uncertainties, as follows: (i) 180 group-averaged microscopic total cross sections; (ii) 7101 group-averaged microscopic scattering cross sections; (iii) 60 group-averaged microscopic fission cross sections; (iv) 60 “averaged” number of neutron per fission; (v) 60 group-averaged fission spectrum constants; (vi) 10 external neutron source parameters; (vii) 6 isotopic number densities. The vector α , which appears in the expression of the Boltzmann-operator B g ( α ) , represents the “vector of uncertain model parameters”.
The fundamental quantity underlying the neutron measurements is the total leakage of neutrons leaving the PERP sphere, represented mathematically as follows:
L ( α ) S b d S g = 1 G Ω n > 0 d Ω Ω n φ g ( r , Ω ) .
The expressions of the 1st-order sensitivities of the leakage response to the model parameters underlying the total cross section are as follows:
L ( α ) σ t , i g = g = 1 G d V 4 π d Ω ψ ( 1 ) , g ( r , Ω ) φ g ( r , Ω ) Σ t g ( t ; r ) σ t , i g , i = 1 , , I ; g = 1 , , G ,
where the multigroup adjoint fluxes ψ ( 1 ) , g ( r , Ω ) are the solutions of the following 1st-Level Adjoint Sensitivity System (1st-LASS):
A ( 1 ) , g ( α ) ψ ( 1 ) , g ( r , Ω ) = Ω n δ ( r r d ) , g = 1 , , G ,
ψ ( 1 ) , g ( r d , Ω ) = 0 , Ω n > 0 , g = 1 , , G ,
where
A ( 1 ) , g ( α ) ψ ( 1 ) , g ( r , Ω ) Ω · ψ ( 1 ) , g ( r , Ω ) + Σ t g ( t ; r ) ψ ( 1 ) , g ( r , Ω ) g = 1 G 4 π d Ω Σ s g g ( s ; r , Ω Ω ) ψ ( 1 ) , g ( r , Ω ) ν Σ f g ( f ; r ) g = 1 G 4 π d Ω χ g ( p ; r ) ψ ( 1 ) , g ( r , Ω ) , g = 1 , , G .
The 2nd-order sensitivities of the leakage response to the parameters involved in the definitions of the total cross sections are given by the following expression:
2 L t j t m 2 = V d V 4 π d Ω [ ψ 1 , j ( 2 ) , g m 2 ( r , Ω ) ψ ( 1 ) , g m 2 ( r , Ω ) + ψ 2 , j ( 2 ) , g m 2 ( r , Ω ) φ g m 2 ( r , Ω ) ] N i m 2 , m m 2 , f o r j = 1 , , J σ t ; m 2 = 1 , , J σ t ,
where the 2nd-level adjoint functions ψ 1 , j ( 2 ) , g , j = 1 , , J σ t ; g = 1 , , G , and ψ 2 , j ( 2 ) , g , j = 1 , , J σ t ; g = 1 , , G , are the solutions of the following Second-Level adjoint Sensitivity system (2nd-LASS):
B g ( α 0 ) ψ 1 , j ( 2 ) , g ( r , Ω ) = δ g j g N i j , m j φ g ( r , Ω ) , j = 1 , , J σ t ; g = 1 , , G ,
ψ 1 , j ( 2 ) , g ( r d , Ω ) = 0 , Ω n < 0 ; j = 1 , , J σ t ; g = 1 , , G ,
A ( 1 ) , g ( α 0 ) ψ 2 , j ( 2 ) , g ( r , Ω ) = δ g j g N i j , m j ψ ( 1 ) , g ( r , Ω ) , j = 1 , , J σ t ; g = 1 , , G ,
ψ 2 , j ( 2 ) , g ( r d , Ω ) = 0 , Ω n > 0 ; j = 1 , , J σ t ; g = 1 , , G .
All of the quantities appearing in Equations (14)–(17) are evaluated at the nominal values of all model parameters. The details of these computations are presented in [7]. All of the computations were performed using a DELL computer AMD FX-8350 having an 8-core processor.
The leakage response displayed the largest 1st- and 2nd-order sensitivities with respect to the (group-averaged) total microscopic cross sections for 1H; these sensitivities are presented in Table 1 and Table 2; complete results are presented in [7,8,9,10,11,12].
For the first-order sensitivities, the 2nd-ASAM required ca. 5 min CPU time was needed for computing: (a) the forward group-fluxes φ g ( r , Ω ) by solving Equations (1) and (2); (b) the adjoint group-fluxes ψ ( 1 ) , g ( r , Ω ) by solving Equations (11) and (12); and (c) the 7477 integrals expressed by Equation (10) to obtain all of the first-order sensitivities.
By comparison, it would require 694 h-CPU to compute approximately the first-order responses sensitivities using the two-point finite difference scheme shown below:
R α i R i + 1 R i 1 2 ( δ α i ) = O ( δ α i ) 2 , i = 1 , , N α ,
where R i + 1 R ( α i 0 + δ α i ) and R i 1 R ( α i 0 δ α i ) . Note that Equation (19) introduces its own intrinsic methodological error when approximating the respective derivative; this error would be present even if the numerical computation of the response R ( α i ) at the respective “sampling point” α i were perfect. Additionally, one needs to “play around” to find the “correct” value to use for δ α i , because if δ α i is either too large or too small, Equation (19) would produce erroneous results. The 694 h-CPU required to compute approximately the first-order responses sensitivities using Equation (19) does not include the time needed to find the “appropriate” δ α i which gives a “satisfactory” approximation to the respective derivative.
As indicated by Equation (14), the computation of each of the 27,956,503 distinct non-zero second-order sensitivities of the response (to two non-zero parameters) requires one forward PARTISN computation to obtain the function ψ 1 , j ( 2 ) , g by solving Equations (15) and (16), as well as one adjoint PARTISN computation to obtain the function ψ 2 , j ( 2 ) , g by solving Equations (17) and (18). The total CPU time needed to compute the 27,956,503 distinct non-zero second-order sensitivities is ca. 929 h-CPU, as follows (i) 735 h-CPU needed to perform the 7427 forward PARTISN runs + 7415 PARTISN adjoint runs to compute the 2nd-level adjoint functions ψ 1 , j ( 2 ) , g and ψ 2 , j ( 2 ) , g ; (ii) 194 CPU-hours needed to compute the integrals in Equation (14) to obtain the unmixed and mixed second-order sensitivities.
Consider, for comparison, the simplest finite-difference scheme for computing the second-order responses sensitivities, namely:
2 R α i α i R i + 1 2 R i + R i 1 ( δ α i ) 2 = O ( δ α i ) 2 , i = 1 , , N α ,
2 R α i α j R i + 1 , j + 1 R i 1 , j + 1 R i + 1 , j 1 + R i 1 , j 1 4 δ α i δ α j = O [ ( δ α i ) 2 , ( δ α j ) 2 ] , i = 1 , , N α , j = 1 , , i ,
where R i , j R ( α i 0 + δ α i , α j 0 + δ α j ) , etc., are computed response values at the indicated “sampling points.” Note again that the finite-difference schemes introduce their own intrinsic methodological error when approximating the respective derivative; this error would be present even if the numerical computation of the response R ( α i ) at the respective “sampling point” α i were perfect. Furthermore, one needs to find by “trial and error” the “correct” value to use for δ α i , because if δ α i is either too large or too small, the finite-difference schemes would produce wrong numbers. Not counting the “trial and error” the “correct” value to use for δ α i , the CPU time (using a DELL computer with an 8-core processor, AMD FX-8350) needed to compute the 27,956,503 distinct non-zero 2nd-order sensitivities by using Equations (20) and (21) at the respective 111,811,058 “sampling points” would be 592 YEARS-CPU time! Furthermore, these sensitivities would not be exact (as produced by the 2nd-ASAM), but would contain second-order errors. Evidently, it is not feasible to compute the 2nd-order sensitivities using “sampling approaches”.

3. Uncertainty Quantification

The numerical results reported in this paper comprise just the effects of the group-averaged total microscopic cross sections, which will be indicated by using the subscript “t”. Since correlations among the group total cross sections are not available for the PERP benchmark, two extreme situations can be considered: (i) all cross sections are uncorrelated, which will be considered in Section 3.1 below; (ii) all cross sections are fully correlated, which will be considered in Section 3.2 below. The formulas used in this section are as originally used in [12]; for convenient referencing, they are reproduced in Section 3.1 and Section 3.2.

3.1. Uncorrelated Total Microscopic Cross Sections

The expected value of the leakage response has the expression [ E ( L ) ] t ( U ) = L ( α 0 ) + [ E ( L ) ] t ( 2 , U ) , where the superscript “(2,U)” indicates “2nd-order uncorrelated” cross sections and where: (i) L ( α 0 ) denotes the nominal value of the computed leakage response; and (ii) [ E ( L ) ] t ( 2 , U ) = 1 2 g = 1 30 i = 1 6 2 L ( α ) σ t , i g σ t , i g ( s t , i g ) 2 , where s t , i g denotes the standard deviation associated with the group-averaged total microscopic cross sections σ t , i g , for each isotope i. If the total microscopic cross sections are both uncorrelated and normally distributed, which will be indicated using the superscript “(U,N)”, then the following expressions hold: (i) the variance, [ var ( L ) ] t ( U , N ) , of the leakage response of the PERP benchmark becomes [ var ( L ) ] t ( U , N ) = [ var ( L ) ] t ( 1 , U , N ) + [ var ( L ) ] t ( 2 , U , N ) , where the 1st-order term is defined as [ var ( L ) ] t ( 1 , U , N ) g = 1 30 i = 1 6 [ L ( α ) σ t , i g ] 2 ( s t , i g ) 2 and the 2nd-order term is defined as: [ var ( L ) ] t ( 2 , U , N ) 1 2 g = 1 30 i = 1 6 [ 2 L ( α ) σ t , i g σ t , i g ( s t , i g ) 2 ] 2 . Furthermore, the 3rd-order moment, [ μ 3 ( L ) ] t ( U , N ) , of the leakage response has the form [ μ 3 ( L ) ] t ( U , N ) = 3 g = 1 30 i = 1 6 [ L ( α ) σ t , i g ] 2 2 L ( α ) σ t , i g σ t , i g ( s t , i g ) 4 . The skewness, [ γ 1 ( L ) ] t ( U , N ) , of the leakage response, L , due to the variances of uncorrelated and normally distributed microscopic total cross sections is defined as [ γ 1 ( L ) ] t ( U , N ) = [ μ 3 ( L ) ] t ( U , N ) / { [ var ( L ) ] t ( U , N ) } 3 / 2 . If the 2nd-order sensitivities were unavailable, the third moment [ μ 3 ( L ) ] t ( U , N ) and hence the skewness γ 1 of the leakage response would vanish and the response distribution would by default be assumed to be Gaussian.

3.2. Fully Correlated Total Microscopic Cross Sections

The effects of correlations among the group-averaged microscopic total cross sections are transmitted to the response moments (mean value, variance, skewness) through the 2nd-order mixed sensitivities of the leakage response to these cross sections. The exact effect of such correlations cannot be assessed exactly since they are unavailable. When the cross sections are fully correlated (an extreme case denoted by using the superscript “FC”), the effects of the 2nd-order sensitivities can be quantified as follows:
(i)
The mean value of the leakage response becomes [ E ( L ) ] t ( F C ) = L ( α 0 ) + [ E ( L ) ] t ( 2 , F C ) , where [ E ( L ) ] t ( 2 , F C ) 1 2 g = 1 30 g = 1 30 i = 1 6 k = 1 6 2 L ( α ) σ t , i g σ t , k g ( s t , i g s t , k g ) . The quantity [ E ( L ) ] t ( 2 , M S C ) = [ E ( L ) ] t ( 2 , F C ) [ E ( L ) ] t ( 2 , U ) , where the superscript “MSC” stands for “mixed second-order correlated,” quantifies the 2nd-order contributions to the mean value of the stemming from the mixed 2nd-order sensitivities when the total cross sections are fully correlated.
(ii)
When the cross sections are normally distributed and fully correlated, the variance of the response is [ var ( L ) ] t ( F C , N ) = [ var ( L ) ] t ( 1 , F C , N ) + [ var ( L ) ] t ( 2 , F C , N ) , where the contributions from the 1st-order sensitivities are contained in the term [ var ( L ) ] t ( 1 , F C , N ) g = 1 30 g = 1 30 i = 1 6 k = 1 6 [ L ( α ) σ t , i g L ( α ) σ t , k g ] ( s t , i g s t , k g ) while the contributions from the 2nd-order sensitivities are contained in the term [ var ( L ) ] t ( 2 , F C , N ) 1 2 [ g = 1 30 g = 1 30 i = 1 6 k = 1 6 2 L ( α ) σ t , i g σ t , k g s t , i g s t , k g ] [ g = 1 G g = 1 G μ = 1 I ν = 1 I 2 L ( α ) σ t , μ g σ t , ν g s t , μ g s t , ν g ] . The superscript “(FC,N)” indicates “fully correlated, normally distributed.” The contributions to [ var ( L ) ] t ( F C , N ) involving the first-order sensitivities will be denoted as [ var ( L ) ] t ( 1 , M S C ) and are obtained by subtracting the uncorrelated terms from the fully correlated ones, i.e., [ var ( L ) ] t ( 1 , M S C ) = [ var ( L ) ] t ( 1 , F C , N ) [ var ( L ) ] t ( 1 , U , N ) . The superscript “(1, MSC)” denotes “first-order, mixed sensitivities, correlated.” Similarly, the quantity [ var ( L ) ] t ( 2 , M S C ) = [ var ( L ) ] t ( 2 , F C ) [ var ( L ) ] t ( 2 , U , N ) represents the contributions to [ var ( L ) ] t ( F C , N ) involving the 2nd-order mixed and correlated sensitivities.
When the total cross sections are fully correlated, the skewness, [ γ 1 ( L ) ] t ( F C , N ) = [ μ 3 ( L ) ] t ( F C , N ) / { [ var ( L ) ] t ( F C , N ) } 3 / 2 of the response distribution response involve a third order cumulant that has the following expression:
[ μ 3 ( L ) ] t ( F C , N ) = [ g = 1 30 g = 1 30 i = 1 6 k = 1 6 L ( α ) σ t , i g s t , i g L ( α ) σ t , k g s t , k g ] [ g = 1 30 g = 1 30 μ = 1 6 ν = 1 6 2 L ( α ) σ t , μ g σ t , ν g s t , μ g s t , ν g ] + [ g = 1 30 g = 1 30 i = 1 6 ν = 1 6 L ( α ) σ t , i g s t , i g L ( α ) σ t , ν g s t , ν g ] [ g = 1 30 g = 1 30 k = 1 6 μ = 1 6 2 L ( α ) σ t , k g σ t , μ g s t , k g s t , μ g ] + [ g = 1 30 g = 1 30 i = 1 6 k = 1 6 2 L ( α ) σ t , i g σ t , k g s t , i g s t , k g ] [ g = 1 30 g = 1 30 μ = 1 6 ν = 1 6 L ( α ) σ t , μ g s t , μ g L ( α ) σ t , ν g s t , ν g ] .
When the total cross sections are fully correlated and normally distributed, the quantity [ μ 3 ( L ) ] t ( M S C , N ) = [ μ 3 ( L ) ] t ( F C , N ) [ μ 3 ( L ) ] t ( U , N ) provides the contributions from the “mixed second-order sensitivities and fully correlated normally distributed parameters,” which is indicated using the superscript “(MSC,N)”.
Table 3 presents illustrative results when considering uniform relative standard deviations s t , i g of 5% and 10%, respectively. These numerical results highlight the importance of both the unmixed and mixed 2nd-order sensitivities for: (a) causing the mean value of the leakage response to differ from its computed value; (b) contributing significantly to increase the total response variance [ var ( L ) ] t ( U ) ; (c) causing the response distribution to become non-Gaussian.

4. Concluding Remarks

The following conclusions can be drawn based on the results presented in Table 3:
  • The 2nd-order sensitivities become more important than the 1st-order sensitivities when the parameters’ standard deviations and correlation become larger.
  • The 2nd-order sensitivities cause the mean value of the leakage response to differ from its computed value.
  • The 2nd-order sensitivities cause the leakage distribution in parameter space to be skewed towards positive values relative to the response’s mean value.
  • In reality, the total cross sections are partially correlated, but these correlations are currently unavailable. The results (presented in Table 3) involving the mixed 2nd-order sensitivities highlight the need for future experimental research for quantifying these currently unavailable correlations.
  • It has also been shown that the Second Order Adjoint Sensitivity Analysis Methodology (2nd-ASAM) conceived by Cacuci [13,14,15] is the only practical method for computing the second-order sensitivities exactly (i.e., without introducing methodological errors) for realistic and practical models of large-scale systems, which invariably involve many uncertain parameters. In particular, it has been shown that the simplest “2-point sampling method” requires for computing the 1st-order sensitivities (while introducing second-order errors!) almost as much time as the 2nd-ASAM requires for computing (without methodological errors!) all of the 2nd-order sensitivities. Additionally, it has also been shown that even the simplest “sampling approach” requires unthinkably large amount of computational time, while still producing approximate results, subject to “sampling approach errors”.
The complete 2nd-order sensitivity analysis of the PERP benchmark can be found in References [7,8,9,10,11,12]. The complete neutron transport computational model of this benchmark comprises 21,886 imprecisely known nuclear data parameters, which give rise to 21,886 first-order and 478,996,996 second-order sensitivities. The analysis presented in [7,8,9,10,11,12] indicates that many of the 2nd-order sensitivities of the leakage response to the nuclear data underlying this benchmark are more important than the 1st-order ones. While the effects of 2nd-order response sensitivities may be less markedly important for other reactor physics systems than they are for the leakage response of the PERP benchmark, the results presented in this work and in [7,8,9,10,11,12] clearly underscore the need for computing them, rather than considering them a priori to be negligible, as has been the practice thus far in the published literature.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. International Handbook of Evaluated Criticality Safety Benchmark Experiments; Nuclear Energy Agency: Paris, France, 2016.
  2. Brandon, E. Assembly of 239Pu Ball for Criticality Experiment; CMB-11-FAB-80-65; Los Alamos National Laboratory: Los Alamos, NM, USA, 1980. [Google Scholar]
  3. Miller, E.C.; Mattingly, J.K.; Clarke, S.D.; Solomon, C.J.; Dennis, B.D.; Meldrum, A.; A Pozzi, S. Computational Evaluation of Neutron Multiplicity Measurements of Polyethylene-Reflected Plutonium Metal. Nucl. Sci. Eng. 2014, 176, 167–185. [Google Scholar] [CrossRef]
  4. Cacuci, D.G. Application of the Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology to Compute 1st- and 2nd-Order Sensitivities of Flux Functionals in a Multiplying System with Source. Nucl. Sci. Eng. 2019, 193, 555–600. [Google Scholar] [CrossRef] [Green Version]
  5. Alcouffe, R.E.; Baker, R.S.; Dahl, J.A.; Turner, S.A.; Ward, R. PARTISN: A Time-Dependent, Parallel Neutral Particle Transport Code System; LA-UR-05-3925; Los Alamos National Laboratory: Los Alamos, NM, USA, 2008. [Google Scholar]
  6. Conlin, J.L.; Parsons, D.K.; Gardiner, S.J.; Gray, M.G.; Lee, M.B.; White, M.C. MENDF71X: Multigroup Neutron Cross-Section Data Tables Based upon ENDF/B-VII.1X; LA-UR-15-29571; Los Alamos National Laboratory: Los Alamos, NM, USA, 2013. [Google Scholar]
  7. Cacuci, D.G.; Fang, R.; Favorite, J.A. Comprehensive Second-Order Adjoint Sensitivity Analysis Methodology (2nd-ASAM) Applied to a Subcritical Experimental Reactor Physics Benchmark: I. Effects of Imprecisely Known Microscopic Total and Capture Cross Sections. Energies 2019, 12, 4219. [Google Scholar] [CrossRef] [Green Version]
  8. Fang, R.; Cacuci, D.G. Comprehensive Second-Order Adjoint Sensitivity Analysis Methodology (2nd-ASAM) Applied to a Subcritical Experimental Reactor Physics Benchmark: II. Effects of Imprecisely Known Microscopic Scattering Cross Sections. Energies 2019, 12, 4114. [Google Scholar] [CrossRef] [Green Version]
  9. Cacuci, D.G.; Fang, R.; Favorite, J.A.; Badea, M.C.; Di Rocco, F. Rocco Comprehensive Second-Order Adjoint Sensitivity Analysis Methodology (2nd-ASAM) Applied to a Subcritical Experimental Reactor Physics Benchmark: III. Effects of Imprecisely Known Microscopic Fission Cross Sections and Average Number of Neutrons per Fission. Energies 2019, 12, 4100. [Google Scholar] [CrossRef] [Green Version]
  10. Fang, R.; Cacuci, D.G. Comprehensive Second-Order Adjoint Sensitivity Analysis Methodology (2nd-ASAM) Applied to a Subcritical Experimental Reactor Physics Benchmark: IV. Effects of Imprecisely Known Source Parameters. Energies 2020, 13, 1431. [Google Scholar] [CrossRef] [Green Version]
  11. Fang, R.; Cacuci, D.G. Comprehensive Second-Order Adjoint Sensitivity Analysis Methodology (2nd-ASAM) Applied to a Subcritical Experimental Reactor Physics Benchmark: V. Computation of Mixed 2nd-Order Sensitivities Involving Isotopic Number Densities. Energies 2020, 13, 2580. [Google Scholar] [CrossRef]
  12. Cacuci, D.G.; Fang, R.; Favorite, J.A. Comprehensive Second-Order Adjoint Sensitivity Analysis Methodology (2nd-ASAM) Applied to a Subcritical Experimental Reactor Physics Benchmark. VI: Overall Impact of 1st- and 2nd-Order Sensitivities on Response Uncertainties. Energies 2020, 13, 1674. [Google Scholar] [CrossRef] [Green Version]
  13. Cacuci, D.G. Second-order adjoint sensitivity analysis methodology (2nd-ASAM) for computing exactly and efficiently first- and second-order sensitivities in large-scale linear systems: I. Computational methodology. J. Comput. Phys. 2015, 284, 687–699. [Google Scholar] [CrossRef] [Green Version]
  14. Cacuci, D.G. Second-Order Adjoint Sensitivity Analysis Methodology (2nd-CASAM) for Large-Scale Nonlinear Systems: I. Theory. Nucl. Sci. Eng. 2016, 184, 16–30. [Google Scholar] [CrossRef]
  15. Cacuci, D.G. The Second-Order Adjoint Sensitivity Analysis Methodology; CRC Press: Boca Raton, FL, USA; Taylor & Francis Group: Oxfordshire, UK, 2018. [Google Scholar]
Table 1. First-order relative sensitivities, ( L / σ t , i = 6 g ) ( σ t , i = 6 g / L ) , g = 1 , , 30 of the leakage to the total cross sections of 1H (isotope i = 6).
Table 1. First-order relative sensitivities, ( L / σ t , i = 6 g ) ( σ t , i = 6 g / L ) , g = 1 , , 30 of the leakage to the total cross sections of 1H (isotope i = 6).
g1st-OrderG1st-Orderg1st-Orderg1st-OrderG1st-Order
1–6<−0.0111−0.1916−1.1621−9.6926−0.65
7−0.0712−0.4417−1.1722−0.9227−0.58
8−0.0913−0.5218−1.1423−0.8928−0.55
9−0.1414−0.5719−1.0924−0.7529−0.55
10−0.1715−0.5820−1.0325−0.7130−9.366
Table 2. Selected 2nd-order relative sensitivities ( 2 L / σ t , i = 6 g σ t , k = 6 g ) ( σ t , i = 6 g σ t , k = 6 g / L ) , g = 1 , , 30 , having values greater than 1.0, for 1H.
Table 2. Selected 2nd-order relative sensitivities ( 2 L / σ t , i = 6 g σ t , k = 6 g ) ( σ t , i = 6 g σ t , k = 6 g / L ) , g = 1 , , 30 , having values greater than 1.0, for 1H.
g′ = 12131415161718192021222324252627282930
g = 120.6530.3150.3400.3560.7400.7630.7510.7250.6880.6480.5970.5540.5020.4770.4400.3930.3690.3726.432
130.3150.9740.4710.4710.9761.0050.9880.9530.9040.8510.7840.7280.6590.6270.5770.5160.4840.4878.424
140.3400.4711.2610.5791.1581.1921.1721.1301.0721.0090.9300.8630.7820.7420.6840.6110.5740.5769.968
150.3560.4710.5791.3911.2551.2771.2551.2101.1481.0810.9960.9240.8370.7950.7330.6550.6150.61710.67
160.7400.9761.1581.2554.4612.7002.6472.5532.4212.2802.1001.9491.7671.6771.5461.3811.2961.30022.48
170.7631.0051.1921.2772.7004.8532.7892.6842.5462.3982.2092.0501.8581.7641.6251.4521.3631.36723.62
180.7510.9881.1721.2552.6472.7894.8282.6892.5462.3992.2102.0511.8591.7641.6261.4531.3631.36723.62
190.7250.9531.1301.2102.5532.6842.6894.6192.4982.3492.1652.0101.8221.7291.5941.4241.3361.34023.15
200.6880.9041.0721.1482.4212.5462.5462.4984.2842.2662.0851.9361.7551.6661.5351.3721.2871.29022.29
210.6480.8511.0091.0812.2802.3982.3992.3492.2663.9372.0041.8571.6841.5991.4741.3171.2361.23821.40
220.5970.7840.9300.9962.1002.2092.2102.1652.0852.0043.5151.7601.5931.5121.3941.2461.1691.17120.24
230.5540.7280.8630.9241.9492.0502.0512.0101.9361.8571.7603.1771.5211.4401.3281.1871.1141.11619.28
240.5020.6590.7820.8371.7671.8581.8591.8221.7551.6841.5931.5212.7921.3581.2491.1171.0481.04918.13
250.4770.6270.7420.7951.6771.7641.7641.7291.6661.5991.5121.4401.3582.6041.2141.0821.0161.01717.58
260.4400.5770.6840.7331.5461.6251.6261.5941.5351.4741.3941.3281.2491.2142.3491.0370.9710.97216.79
270.3930.5160.6110.6551.3811.4521.4531.4241.3721.3171.2461.1871.1171.0821.0372.0390.9130.91215.76
280.3690.4840.5740.6151.2961.3631.3631.3361.2871.2361.1691.1141.0481.0160.9710.9131.8850.88815.30
290.3720.4870.5760.6171.3001.3671.3671.3401.2901.2381.1711.1161.0491.0170.9720.9120.8881.89115.39
306.4328.4249.9710.6722.4823.6223.6223.1522.2921.4020.2419.2818.1317.5816.7915.7615.3015.39429.6
Table 3. Response moments for various cases.
Table 3. Response moments for various cases.
Fully Correlated Cross SectionsUncorrelated Cross Sections
Rel. St. Dev.10%5%Rel. St. Dev.10%5%
L ( α 0 ) 1.7648 × 10 6 1.7648 × 10 6 L ( α 0 ) 1.7648 × 10 6 1.7648 × 10 6
[ E ( L ) ] t ( 2 , F C ) 2.9451 × 10 7 7.3627 × 10 6 [ E ( L ) ] t ( 2 , U ) 4.5980 × 10 6 1.1495 × 10 6
[ E ( L ) ] t ( 2 , M S C ) 2.4853 × 10 7 6.2132 × 10 6 [ E ( L ) ] t ( U ) 6.3628 × 10 6 2.9143 × 10 6
[ E ( L ) ] t ( F C ) 3.1216 × 10 7 9.1275 × 10 6 [ var ( L ) ] t ( 1 , U , N ) 3.4196 × 10 12 8.5490 × 10 11
[ var ( L ) ] t ( 1 , F C , N ) 4.7601 × 10 13 1.1900 × 10 13 [ var ( L ) ] t ( 2 , U , N ) 2.8789 × 10 13 1.7993 × 10 12
[ var ( L ) ] t ( 1 , M S C ) 4.4181 × 10 13 1.1045 × 10 13 [ var ( L ) ] t ( U , N ) 3.2208 × 10 13 2.6542 × 10 12
[ var ( L ) ] t ( 2 , F C , N ) 1.7347 × 10 15 1.0842 × 10 14 [ μ 3 ( L ) ] t ( U , N ) 6.2267 × 10 19 3.8917 × 10 18
[ var ( L ) ] t ( 2 , M S C ) 1.7059 × 10 15 1.0662 × 10 14 [ γ 1 ( L ) ] t ( U , N ) 0.3407 0.8999
[ var ( L ) ] t ( F C , N ) 1.7823 × 10 15 1.2302 × 10 14
[ μ 3 ( L ) ] t ( F C , N ) 8.4113 × 10 21 5.2571 × 10 20
[ μ 3 ( L ) ] t ( M S C , N ) 8.3490 × 10 21 5.2181 × 10 20
[ γ 1 ( L ) ] t ( F C , N ) 0.1118 0.3983
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cacuci, D.G. Illustrating Important Effects of Second-Order Sensitivities on Response Uncertainties in Reactor Physics. J. Nucl. Eng. 2021, 2, 114-123. https://0-doi-org.brum.beds.ac.uk/10.3390/jne2020012

AMA Style

Cacuci DG. Illustrating Important Effects of Second-Order Sensitivities on Response Uncertainties in Reactor Physics. Journal of Nuclear Engineering. 2021; 2(2):114-123. https://0-doi-org.brum.beds.ac.uk/10.3390/jne2020012

Chicago/Turabian Style

Cacuci, Dan G. 2021. "Illustrating Important Effects of Second-Order Sensitivities on Response Uncertainties in Reactor Physics" Journal of Nuclear Engineering 2, no. 2: 114-123. https://0-doi-org.brum.beds.ac.uk/10.3390/jne2020012

Article Metrics

Back to TopTop