Nuclear Engineering and Technology. 2013. Nov, 45(6): 753-758
Copyright ©2013, Korean Nuclear Society
  • Received : September 09, 2013
  • Accepted : October 10, 2013
  • Published : November 25, 2013
Export by style
Cited by
About the Authors
Research Institute of Nuclear Engineering, University of Fukui
Japan, Egypt Nuclear and Radiological Regulatory Authority 1-2-4 Kanawa-cho, Tsuruga-shi, Fukui, 914-0055, Japan
A new method of calculating sensitivity coefficients of core characteristics relative to infinite-dilution cross sections has been developed. Conventional sensitivity coefficients are evaluated for the changes of effective cross sections which are dependent on individual models of core and cell. Therefore a correction has been derived to the conventional sensitivity coefficients based on the perturbation theory. The accuracy of the present method has been verified by comparing numerical results of sensitivity coefficients with a reference Monte-Carlo method.
Sensitivity coefficients are utilized in the three fields; sensitivity analysis, uncertainty estimation and cross section adjustment. In the sensitivity analysis, the effect of cross section changes, for example from JENDL 3.3 to 4.0, on core characteristics is estimated by calculating the product of sensitivity coefficients and cross section changes, which is important in releasing a new cross section library. In the uncertainty estimation, the prediction uncertainty is evaluated by using the cross section covariance and the sensitivity coefficients 1 . In the cross section adjustment, the sensitivity coefficients play an important role in adjusting the cross sections by using the measured neutronics characteristics based on the Bayes theory. Thus in these three applications, the sensitivity coefficients have important roles, and accurate evaluation is necessary.
Sensitivity coefficients are usually calculated by using the generalized perturbation theory because the calculation requires a lot of computing time. Sensitivity coefficients of individual core characteristics such as k eff , power distribution and control rod worth are calculated relative to cross section changes of each reaction type (fission, capture, scattering) for many energy groups. The generalized perturbation theory shortens the computing time by the use of generalized flux and generalized adjoint flux. TSUNAMI module 2 in SCALE code is a representative sensitivity calculation code. In Japan we have developed SAGEP 3 , SAGEP-T 4 and SAINT 5 codes based on the generalized perturbation theory with use of the diffusion theory, transport S n method, and collision probability method, respectively.
However, the sensitivity coefficients are usually calculated by changing the effective cross sections. The effective cross sections are the product of the infinite- dilution cross sections and the self-shielding factors. The self-shielding factor varies with other cross section changes due to change of the background cross section, and the interaction effect has to be considered in estimating the sensitivities.
The purpose of this paper is to derive a new method to calculate the sensitivity which shows the core characteristic changes due to the change of the infinite dilution cross sections. The present method introduces a correction TERM to the conventional sensitivity coefficient. Therefore, the conventional approximation utilized hitherto with this correction TERM can be used to calculate the improved sensitivities.
To check the accuracy of the present method, numerical calculations are performed for a representative PWR cell model. Approximated (conventional) sensitivity coefficients are calculated by using the SAINT code, which can treat cell geometry using the collision probability method, and then the improved sensitivity coefficients are calculated by using the present method. Finally the improved and approximated sensitivity coefficients are compared with a reference one calculated by MCNP code.
The paper is organized as follow: The relation between improved and approximate sensitivity coefficients is derived in Section II, and numerical results are shown in Sec. III. Section IV summarizes the conclusion.
Core characteristics such as keff, power distribution and control rod worth are calculated by using an effective cross section in deterministic methods. Usually the effective cross sections are calculated by using the Bondarenko selfshielding factor method 1 , the subgroup method1 et al. Here the Bondarenko self-shielding factor method has been used. With that method the effective cross sections are expressed by the infinite dilution cross section and the self-shielding factors ƒ 6 :
Lager Image
The self-shielding factors depend on the background cross section and temperature. The background cross section for nuclide i` in a homogeneous medium is calculated with the formula
Lager Image
where N k is the atomic number density of light nuclide k and σ k t is the microscopic total cross section. In a heterogeneous medium the background cross section for the infinite regular lattice is given by
Lager Image
where C is the Dancoff factor, Σ e is the reciprocal of the fuel diameter, a is the Bell factor, and σ bi` is the background cross defined by Eq. (2).
The sensitivity coefficient is defined by the relative change of the core characteristics due to the relative change of the cross section. Here we consider the following two sensitivities, the sensitivity due to the relative change of the infinite dilution cross sections and the approximate sensitivity due to the relative change of the effective cross sections
Lager Image
Lager Image
The relative change of the infinite dilution cross sections causes the change of the core characteristics as follows:
Lager Image
From Eq. (1), one can obtain
Lager Image
Therefore, the improved sensitivity is expressed by using the approximate sensitivity as follows:
Lager Image
Let us consider the above equation in more detail. The sensitivities and the cross sections are dependent on nuclides, reaction types (such as fission, capture and scattering) and energy groups. Let us consider the case where there is a perturbation in σ of nuclide i, reaction type j in energy group g. This perturbation causes a change in the self-shielding factor ƒ of nuclide i`, reaction type j` in energy group g`. The second term on the right hand side of Eq. (8) has to cover the contributions for all nuclides i`, reaction types j` in energy groups g`, therefore we have to take the summation over i`, j` and g`.
The sensitivity for the nuclide i, reaction type j in energy group g is given by
Lager Image
The first term is the direct contribution to S, which can be calculated using the conventional tools evaluating sensitivity coefficients such as SAGEP, SAGEP-T and SAINT. The second term represents the indirect contribution through the change of self-shielding factor. These coefficients can be calculated as follow:
Here we assume that the wide resonance approximation is applied for heavy nuclides, of-course the narrow resonant approximation is applied for light nuclides. The self-shielding effect depends on the neutron spectrum, where the neutron spectrum for the heavy nuclide i` is written as
Lager Image
Equation (10) indicates that when σ a i` (E) and σ b i` change by the same factor, the neutron spectrum remains the same, This shows that the ratio h i` = σ b i` a i` (E) has the effect on the neutron spectrum, and also on the self-shielding factor.
Now let us consider a case when there are cross section perturbations in σ a i` (E) and σ b i` , this ratio changes to h` i` =σ` b i` /σ` a i` (E), where the superscript ` indicates the perturbed cross sections
When the perturbations of cross sections Δσ = σ` - σ are small compared to σ, so h` can be written as
Lager Image
Equation (11) is interesting, for when the relative changes of Δσ` b i` /σ` b i` and Δσ` a i` (E)/σ` a i` (E) are the same, h`= h, and therefore the self-shielding effect remains the same before perturbation.
The self-shielding effect is dependent on the ratio of σ b to σ a , so the perturbed self-shielding factor can be written as:
Lager Image
Using Eq. (11), Eq. (12) reduces to
Lager Image
When there are some heavy resonant nuclides mixed in fuel (indicated by i and i`), the neutron spectrum in Eq. (10) is replaced by
Lager Image
From Eq. (14) the total background cross section is
Lager Image
Therefore the change in the background cross section due to the cross section perturbation of heavy nuclides i is
Lager Image
because Δσ i` b0 = 0 for these heavy nuclides.
Substituting Eq. (16) to Eq. (13), the change of the self-shielding factor of nuclide i`, reaction j` due to the cross section perturbation of nuclide i, reaction j will be written as
Lager Image
where δ ii` =1 for i=i`, and =0 for i≠i`.
Therefore we obtain:
Lager Image
where dσ i a /dσ i j is usually unity except for scattering reaction;, in which case it is zero.
Usually the ƒ table is described as a function of the background cross section only or a relevant nuclide, and the total cross section is assumed to be fixed in these tables, therefore ∂f i`j` (h i` )/∂ hi` in Eq. (18) can be written as:
Lager Image
∂f/∂ σb is calculated by the finite difference approximation Finally, the coefficient in the second term of the RHS
Lager Image
of Eq. (9) is called TERM and can be written as
Lager Image
When one considers the cross section perturbation of light nuclides k, the narrow resonance approximation has to be used. In this case;, the neutron spectrum is expressed by
Lager Image
and therefore σ a is replaced by σ t in Eq. (20).
The SRAC 7 , SAINT and MCNP 8 codes are used to perform these analyses, where the SRAC code is used for the self-shielding factor calculations, the SAINT code for conventional sensitivity coefficient calculations and the MCNP code for verification. In all calculations, the ENDF/B-VI library was used.
Sensitivities are calculated for the standard PWR –UOX fueled pin-cell. The fuel, clad and water properties are shown in Table 1 .
Isotopic Data for UOX Cell
Lager Image
Isotopic Data for UOX Cell
First, the self-shielding factors are calculated using the SRAC2006 code for different isotopes and reactions. Then the self-shielding gradients with respect to the background cross section are calculated, and from Eq. (20), the TERMs are determined. TERMs are shown in Fig. 1 ~ 4 as a function of energy groups (10 MeV ~ 4.1399E-01 eV).
The TERMs in Figures 1 (a), (b), (c) and (d) represent the changes of U-235 capture, fission, elastic and inelastic self-shielding factors respectively due to perturbations in U-235, U238 and O-16 capture cross sections.
Figures 2 ~ 4 represent the changes of U-238, U-235 and O-16 self-shielding factor due to perturbations in capture, fission and elastic cross sections respectively.
The SAINT code is then used to calculate the conventional approximated k eff sensitivity coefficients for different isotopes and reactions. After calculating the approximated sensitivities and the TERMs, the improved
Lager Image
U-235 / Capture TERM.
Lager Image
U-238 / Capture TERM.
k eff sensitivities can be calculated using Eq. (9). Figures 5 ~ 7 represent the approximated and improved k eff sensitivities for U-235, U-238 and O-16 capture cross sections respectively, while Figures 8 , 9 represent k eff sensitivities for U-235 fission and O-16 elastic cross sections respectively. Figures (a)’s stand for the sensitivity and (b)’s for the differences between approximated and improved sensitivity, where the difference is defined by:
  • Difference= (SImproved–SApprox.)/SImprovedx 100 (%)
It is shown in Figures 5 ~ 9 that the differences between improved and approximated sensitivities are large, at 57 (6.4759 - 8.3153 eV), 53 (1.7604E+01 - 2.2603E+01 eV) and 51 (2.9023E+01 - 3.7266E+01 eV) energy groups, which are the resonant energy groups, and the main contributor to these differences is the U-238 capture cross section.
Lager Image
U-235 / Fission TERM.
Lager Image
O-16 / Elastic TERM.
Lager Image
keff Sensitivity for U-235 Capture Cross Section
Lager Image
keff Sensitivity for U-238 Capture Cross Section
Lager Image
keff Sensitivity for O-16 Capture Cross Section
Lager Image
keff Sensitivity for U-235 Fission Cross Section
Lager Image
keff Sensitivity for O-16 Elastic Cross Section
MCNP5 calculations 9 are performed at these energy groups in order to verify the present formula. The results of MCNP k eff sensitivity coefficients compared to SAINT (Approximated) and improved sensitivity coefficients are shown in Tables 2 , 3 for U-235 Capture and U238-Capture reactions respectively.
Eg is the energy group
Lager Image
keffSensitivity for U-235 Capture Cross Section
Lager Image
keff Sensitivity for U-235 Capture Cross Section
keffSensitivity for U-238 Capture Cross Section
Lager Image
keff Sensitivity for U-238 Capture Cross Section
It can be shown from these tables that Δ2 is smaller than Δ1, which means that the sensitivity calculated by the present model is more accurate and can be used to calculate the sensitivity coefficients.
In order to calculate sensitivity coefficients of neutronics characteristics relative to infinite dilution cross sections, a new formula was introduced. This formula can calculate the sensitivity coefficients using the conventional sensitivity coefficients and the correction term called TERM.
Using the present formula, one can easily revise the conventional sensitivity calculation codes such as SAGEP, SAGEP-T and SAINT.
The present formula has been verified by performing numerical calculations in cell geometry and comparing the results with the reference MCNP code.
The results calculated by the present formula were in good agreement with the reference results, particularly in the resonance energy groups. This agreement shows the usefulness of the present formula, and show the importance of self-shielding factors when calculating sensitivity coefficients in light water reactors.
The authors are very thankful for the kind suggestionsand help to Prof. T. KITADA of Osaka University, andfor kind help to Prof. Van Rooijen of the University ofFukui.
CACUCI D. G. (2010) Handbook of Nuclear Engineering, Volume 3 Springer Science + Business Media LLC Germany Chapter 17
REARDEN B. (2009) TSUNAMI-1D: Control Module for One-Dimensional Cross-Section Sensitivity and UncertaintyAnalysis for Criticality, Version 6 Oak Ridge National Laboratory ORNL0TM-2005039
HARA M. , TAKEDA T. (1984) SAGEP: Two-DimensionalSensitivity Analysis Code Based on General PerturbationTheory Japan JAERI-M 84-027 1 - 30
TAKEDA T. , ASANO K. (2006) Sensitivity Analysis Basedon Transport Theory Journal of Nuclear Science and Technology Japan 43 (7) 743 - 749
TAKEDA T. , NAKANO M. (1986) New Calculation Methodof Sensitivity Coefficients of Cell Parameters and ItsApplication Journal of Nuclear Science and Technology Japan 23 (8) 681 - 694
REUSS P. (2008) Neutron Physics EDP Sciences France Chapter 8
OKUMURA K. , KUGO T. (2007) SRAC2006: A Comprehensive Neutronics Calculation Code System Japan JAEA-Data/Code
X-5 Monte Carlo Team 2004 MCNP — A General Monte CarloN-Particle Transport Code, Version 5 Los Alamos, USA LA-CP-03-0245
ROCHMAN D. , KONING A.J. (2011) Nuclear data uncertaintypropagation: Perturbation vs. Monte Carlo Annals of Nuclear Energy 38 942 - 952