EPJ Nuclear Sci. Technol.
Volume 4, 2018
|Number of page(s)||6|
|Published online||17 October 2018|
Extension of Bayesian inference for multi-experimental and coupled problem in neutronics − a revisit of the theoretical approach
French Atomic Energy and Alternative Energies Commission CEA, DEN, CAD, Reactor Studies Department,
Saint Paul-Lez-Durance, France
* e-mail: firstname.lastname@example.org
Received in final form: 22 April 2018
Accepted: 9 July 2018
Published online: 17 October 2018
Bayesian methods are known for treating the so-called data re-assimilation. The Bayesian inference applied to core physics allows to get a new adjustment of nuclear data using the results of integral experiments. This theory leading to reassimliation encompasses a broader approach. In previous papers, new methods have been developed to calculate the impact of nuclear and manufacturing data uncertainties on neutronics parameters. Usually, adjustment is performed step by step with one parameter and one experiment by batch. In this document, we rewrite Orlov theory to extend to multiple experimental values and parameters adjustment. We found that the multidimensional system expression looks like can be written as the monodimensional system in a matrix form. In this extension, correlation terms appears between experimental processes (manufacturing and measurements) and we discuss how to fix them. Then formula are applied to the extension to the Boltzmann/Bateman coupled problem, where each term could be evaluated by computing depletion uncertainties, studied in previous papers.
© T. Frosio et al., published by EDP Sciences, 2018
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The required accuracies on the nuclear data are difficult to reach using only differential experiments, even if innovative experimental techniques are used. The use of integral experiments has been essential in the past to ensure enhanced predictions for power fast reactor cores. In some cases, these integral experiments have been documented in an effective manner and the associated uncertainties are well understood. A combined use of scientifically based covariance data and of integral experiments can be made using advanced statistical adjustment techniques (see, e.g. [1–4]). These techniques can provide in a first step adjusted nuclear data [5,6] for a wide range of applications with new and improved covariance data and bias factors (with reduced uncertainties) for the required design parameters, in order to meet design target accuracies. Moreover, the purpose of cross sections adjustment is more and more perceived as of providing useful feedback to evaluators and differential measurement experimentalists. It then allows to improve the knowledge of neutron cross sections to be used in a wider range of applications.
However, results concerning the nuclear data propagation on GEN-II and GEN-III reactor parameters, and even more for MTR are scarce: if uncertainty quantification methods are well established for Boltzmann problem, as well as separated Bateman problems, an accurate and rigorous treatment of nuclear data uncertainty propagation in coupled problem is still missing.
An important gap between “step 0” uncertainty calculation and depletion uncertainty calculation must then be filled. The major unknown, uncertainties on isotopes concentrations in the reactor core can be estimated by decorrelating sources of uncertainties.
The large amount of integral experiments performed on different nuclear reactors could be taken into account to improve the knowledge of the nuclear data and covariance matrix, in addition to nuclear data differential experiments. These integral experiments could be based on different integral parameters such as reactivity but also radionuclide concentrations, power distribution or reactivity coefficients. They could be simultaneously used for nuclear data adjustment as we will explain in this document.
After the description of the representativity theory, we extend the concept to multiple experiments and then we take into account Boltzmann/Bateman coupled problem.
Sensitivity coefficients and uncertainties can be used to adjust nuclear data and have a feedback on neutronics parameters of interest from integral experiments. An example of adjustment with all the necessary steps (sensitivity analysis, interest parameters uncertainties, experiment analysis) can be found in .
We consider a nuclear data set σ0 associated to a covariance matrix Mσ from a data evaluation. The representativity or data adjustment has been introduced by Usachev . This method allows to reduce the uncertainties specified in Mσ called “a priori” to get a new “a posteriori” covariance matrix taking into account the information contained in the experiments as used in  for super-Phenix reactor. This covariance matrix is associated to σ* the posterior nuclear data re-evaluated during the same process of adjustment.
One of the important hypothesis done for this equation (1) is based on a first order Taylor expansion to linearize the C parameters for small nuclear data perturbations around their expectancies. This can be translated by equation (2): (2)
C(σ*) refers to the value of C computed with σ* and C is linked to σ0.
Different reasoning ways  can give the next equation (3). Considering this set of integral experiments measured in E experiments, we can show using Bayes theorem and maximizing the likelihood, with generalized least-squares, and applying Sherman-Morrison-Woodbury formulation of matricial inverse: (3) with:
σ0 is the evaluation vector of the nuclear data (prior). The dimension of σ0 is q;
C is the integral parameters vector calculated with a simulation code and with the σ0 nuclear data. The dimension of C is p;
E is the experimental vector of the integral parameters. The dimension of E is p;
SC is the sensitivity matrix of the C interest parameters to σ0. The dimensions of SC is q × p;
SE is the sensitivity matrix of the E interest parameters to σ0. The dimensions of SE is q × q;
ME is the covariance matrix associated to the integral parameters observed in the E experiments. The dimension of ME is p × p, and represents the experimental uncertainty. In a general meaning, ME defined as the experimental uncertainty represents:
the knowledge we have about the results given by the measurement (uncertainty of the measurement device or measurement process);
the knowledge we have on the manufacturing parameters such as geometry, position of assemblies, enrichment, mass balance.
σ* is the adjuster vector of the nuclear data (posterior). The dimension of σ* is q.
is the new covariance matrix associated to σ*. The dimension of is q × q.
σ* and are built during the same mathematical process and associated in a system of equations. They are fully linked and cannot be used independently.
The transposition is an application of the representativity. It allows to transmit the result of an experiment E modelled in a calculation code C to another application R evaluated with a calculation code in case that the representativity between the experience and the other application is high (close to 1).
whit R* represents the integral parameters estimated for another design with the adjusted set σ* of nuclear data and is the posterior covariance matrix for the integral parameters of R. We recognize MR = SRMσSRT the a priori covariance matrix of the integral parameters of R.
Considering one integral experiment and one parameter of interest associated to this experiment, we find the Orlov formulation : (5) with:
rE,R the representativity between the reactor and the experiment E for the neutronics parameter R. It is a Pearson correlation coefficient;
ω is a neutronics weight describing the quality of the experiment. ε2(R) and ε*2(R) are respectively the uncertainties a priori and a posteriori of the parameter R. It provides information on the contribution of the experiment in terms of uncertainty knowledge.
ε(X) representing the standard deviation of X, X being any parameter;
δ2 represents the experimental uncertainty. We use this notation instead of ε to not overload the equations with indices to distinguish nuclear data uncertainties and experimental uncertainties.
In the following, we will work with the transposition process starting with equations (4). However, the nuclear data adjustment is following exactly the same assumptions and transformations, starting with equation (3) instead of (4). A work has already been done in  to consider different experiments for a set of nuclear data, but we want here to give details on the experimental covariance matrix ME and express experimental correlations.
In the case where experiments are performed for example on local neutronics parameters such as power factors, or reactivities from different benchmark, the number of experimental results can be important. This is the reason why we now want to extend the data assimilation to cases with multiple experiments and for different neutronics parameters. These parameters can be those studied in  and : isotopic concentrations, power factors and reactivity. Different design studies and different experiments not necessarily correlated can be used for the adjustment. Starting from equation (4), we will express . This matrix is depending on the experimental data base. We need to express ME and will use the methodology proposed by Dos Santos  for this purpose.
ME is a matrix containing the experimental sources of uncertainties that we will write δ2 to not have redundancy with the ε for nuclear data uncertainties. The difficulty is to express these terms because they are evaluated by the experimenter in a first part and by the manufacturing uncertainty propagation in the second part. We consider that these two sources of uncertainties are not correlated and then, the Pearson coefficient between measurements and manufacturing data is null. Manufacturing data consist in all the data associated to the construction of the nuclear core such as mass balance, geometry of the core or position of elements of the core.
represents the variances/covariances coming from Manufacturing Data (MD);
the uncertainty coming from the measurement method.
If we consider different experiments, these coefficients need to take into account correlations between measurements and manufacturing processes.
The evaluation of can be performed using the method described in .
For measurement uncertainties, the formulation is similar. The difficulty associated to these equations is the quantification of the Pearson correlations rMD(Ei, Ej) and rMEAS(Ei, Ej). A method has been proposed in .
To calculate the correlation factors, it is necessary in a first step to identify systematic uncertainties between Ei and Ej and statistical uncertainties. Systematic uncertainties are characterized by a correlation factor of 1 whereas statistical uncertainties have a correlation factor of 0. Then, correlation factor can be estimated as follow: (9)
The formulation is similar for measurement uncertainties.
Examples of systematic and statistical uncertainties are given in .
In the case of impossibilities to describe uncertainties as systematic and statistical, the following method is applied:
The value of rMD(Ei, Ej) can be estimated as follow:
if the experiments Ei and Ej are performed on reactor cores which are manufactured in the same way, by the same manufacturer, we assume ;
otherwise, if the reactor cores are built with different manufacturing processes, .
To express measurement variances/covariances, we consider that the measures are fully correlated if the measurement technique for two experiments is based on the same process/measurement device. At the opposite, if the measurement techniques are different, we assume that the correlations are null. Then:
In the next paragraph, we apply the adjustment system for multidimensional experiments.
Now we consider we know N integral parameters measurements and the associated simulations and we want to calculate the posterior uncertainties of Q integral parameters (more general case). Starting with equation (4), we want to describe the system more precisely using matrix notation to make appear the Orlov coefficients ω and r: (10)
Detailing the matrix with indices, equation (10) can be written with the following notations:
Multiplying equation (10) by and we transform the system as follow: (11)
Finally, in this equation we recognize the Orlov representativity in the matrix C and the neutronic weight in the matrix B. If we come back to the assumptions of equation (5), considering one experiment and one application design, we find the monodimensional expression (Eq. (5)).
In the case where some experiments of the database are performed for example at different burnup to take into account the irradiation process, the method needs to be extended. Nuclide transmutation has to be added in the system.The new uncertainty terms coming from fission yield and transmutation have to be taken into account. These methods are described in  for transmutation terms, and in  for the fission yield uncertainties. Depletion perturbation theory  can also be used.
MR,σ and MR,FY the variances/covariances matrix of coming respectively from the cross-sections and from the fission yields;
SR(D), SR(T) and SR(FY) the sensitivities matrix of to respectively the cross-sections (direct effect), the isotopic concentration (transmutation effect) and the fission yields as defined in ;
Fσ and FFY the Jacobian matrix describing the variations of isotopic concentrations respectively to variations of the cross-sections and the fission yields.
Restarting from equation (4), (13)
With the new notation,, and SE is equivalent to SR but associated to the set of E experiments and MR = MR,σ + MR,FY.
The total uncertainty ε(Ei) of a neutronics parameter Ei now depends on the cross-sections during core evolution and the fission yields. We consider they are independent because the correlations do not exist in the literature between cross sections and fission yields. The new (4) matricial system allows to take into account the isotopic composition of the core during time evolution.
In the case that some experiments represent an isotopic concentration evolution measurement, then the direct term SR(D) = SE(D) = 0 and SR(T) = SR(FY) = I, I being the identity matrix. The matrix system (4) written adding this new information and performing the same transformations than from the previous paragraph becomes: (14)
The different quantities SR, SE, , , have to be calculated taking into account the evolution step of the core.
Accurate knowledge of nuclear data represents a major challenge in numerical simulations for reactor physics, since they are clearly identified as being the major source of propagated uncertainties on reactor integral parameters, such as reactivity, power distributions, and isotopic inventories. These nuclear data, and in particular microscopic cross sections, evaluated through differential experiments sometimes suffer from uncertainties currently not compatible with target uncertainties on some integral parameters, such as reactivity, reactivity effects, or reaction/transmutation rates. The combination of uncertainties propagation results described in [13–15], with inferential methods called “representativity” or “transposition”, extended to the use of several integral parameters and with the results of experiments enables the re-estimation (sometimes still called adjustment) of cross sections of which remains to validate the effectiveness.
The construction of the posterior covariance matrix, i.e. taking into account measurement, is associated with adjusted cross sections that can be called posterior cross sections. On the other hand, the adjustment needed the consideration of the full set of cross sections which are impacting the interest parameters uncertainties . Indeed, experimental data are collected relating to all cross sections since nuclear reactors are macroscopic systems configured by the substance they contain. These sections are all therefore affecting the quality of the result.
The extension of the inference method performed in this paper deals with two new items:
the fuel depletion, which slightly modifies the variables used in the calculation of uncertainties;
the cumulative uncertainty from transmutation effect. These uncertainties affect isotopic concentrations which in return modify the neutronics parameters of interest. If the parameter of interest is the isotopic concentration, then the direct effects are vanishing and some parameters of the equation do not have to be evaluated.
The second aspect to consider is the inclusion of measurements of local parameters as power factors or isotopic concentrations, which can reveal experimental correlations, because the measurement method is the same. It is then possible to extend the inference method to multiple and correlated experiments.
All the authors have contributed equally to the work presented in this manuscript.
- V.V. Orlov, A.A. Van'kov, V.I. Voropaev, Y.A. Kazansky, V.I. Matveev, V.M. Murogov, Problems of fast reactor physics related to breeding, Energy Rev. 4, 991 (1980) [Google Scholar]
- A. Gandini, On transposition of experimental reactor data to reference design, Technical report, Comitato nazionale per la ricerca e per lo sviluppo dell energia nucleare e delle energie alternative, 1983 [Google Scholar]
- P. Blaise, E. Fort, Resonance parameter adjustement methodoloby based on integral experiment analysis, Nucl. Sci. Eng. 133, 235 (1999) [CrossRef] [Google Scholar]
- T. Ivanova, E. Ivanov, F. Ecrabet, Uncertainty assessment for fast reactors based on nuclear data adjustment, Nucl. Data Sheets 118, 592 (2014) [CrossRef] [Google Scholar]
- C. De Saint Jean et al., Evaluation of cross section uncertainties using physical constraints: focus on integral experiments, Nucl. Data Sheets 123, 178 (2015) [CrossRef] [Google Scholar]
- E. Privas et al., Generation of 238U covariance matrices by using the integral data assimilation technique of the CONRAD code, EPJ Web Conf. 106, 04015 (2016) [CrossRef] [Google Scholar]
- G. Palmiotti et al., Combined use of integral experiments and covariance data, Nucl. Data Sheets 118, 596 (2014) [CrossRef] [Google Scholar]
- L.N. Usachev, Yu.G. Bokkov, INDC(CCP)-19/U, International Nuclear Data Committee, Vienna, 1972 [Google Scholar]
- G. Palmiotti, M. Salvatores, Use of integral experiments in the assessment of large liquid-metal fast breeder reactor basic design parameters, Nucl. Sci. Eng. 87, 333 (1984) [CrossRef] [Google Scholar]
- C. De Saint Jean, E. Dupont, M. Ishikawa, G. Palmioti, M. Salvatores, Assessment of existing nuclear data adjustment methodologies, Nuclear Science NEA/WPEC-33, 2010 [Google Scholar]
- P. Archier et al., COMAC, Nuclear data covariance matrices library for reactor applications, in Physor 2014 [Google Scholar]
- N. Terranova et al., Covariance generation and uncertainty propagation for thermal and fast neutron induced fission yields, EPJ Web Conf. 146, 02013 (2017) [CrossRef] [Google Scholar]
- T. Frosio, T. Bonaccorsi, P. Blaise, Nuclear data uncertainties propagation methods in Boltzmann/Bateman coupled problem: application to reactivity in MTR, Ann. Nucl. Energy 90, 303 (2016) [CrossRef] [Google Scholar]
- T. Frosio, T. Bonaccorsi, P. Blaise, Fission yields and cross section uncertainty propagation in Boltzmann/Bateman coupled problems: global and local parameters analysis with a focus on MTR, Ann. Nucl. Energy 98, 43 (2016) [CrossRef] [Google Scholar]
- T. Frosio, T. Bonaccorsi, P. Blaise, Manufacturing data uncertainties propagation method in burn-up problems, Sci. Technol. Nucl. Ins. 2017, 7275346 (2017) [Google Scholar]
- Nuclear Science, Methods and issues for the combined use of integral experiments and covariance data, NEA/NSC/WPEC/DOC(2013)445, 2013 [Google Scholar]
- M. Salvatores et al., Needs and issues of covariance data application, Nucl. Data Sheets 109, 2725 (2008) [CrossRef] [Google Scholar]
- N. Dossantos, Optimisation de l'approche de représentativité et de transposition pour la conception neutronique de programmes expérimentaux dans les maquettes critiques, PHD 2013, GRENI033 [Google Scholar]
- W.S. Yang, T.J. Downar, Depletion perturbation theory for the constrained equilibrium cycle, Nucl. Sci. Eng. 102, 365 (1989) [CrossRef] [Google Scholar]
- G. Palmiotti, M. Salvatores, Developments in sensitivity methodologies and the validation of reactor physics calculations, Sci. Technol. Nucl. Ins. 2012 529623 (2012) [Google Scholar]
Cite this article as: Thomas Frosio, Thomas Bonaccorsi, Patrick Blaise, Extension of Bayesian inference for multi-experimental and coupled problem in neutronics − a revisit of the theoretical approach, EPJ Nuclear Sci. Technol. 4, 19 (2018)
Initial download of the metrics may take a while.