scholarly journals The Impact of Fueling Operations on Full Core Uncertainty Analysis in CANDU Reactors

2020 ◽  
Vol 6 (3) ◽  
Author(s):  
M. D. Tucker ◽  
D. R. Novog

Abstract Within emerging best-estimate-plus-uncertainty (BEPU) approaches, code output uncertainties can be inferred from the propagation of fundamental or microscopic uncertainties. This paper examines the propagation of fundamental nuclear data uncertainties though the entire analysis framework to predict macroscopic reactor physics phenomena, which can be measured in Canada Deuterium Uranium (CANDU) reactors. In this work, 151 perturbed multigroup cross sections libraries, each based on a set of perturbed microscopic nuclear data, were generated. Subsequently, these data were processed into few-group cross sections and used to generate full-core diffusion models in PARCS. The impact of these nuclear data perturbations leads to changes in core reactivity for a fixed set of fuel compositions of 4.5 mk. The impact of online fueling operations was simulated using a series of fueling rules, which attempted to mimic operator actions during CANDU operations such as studying the assembly powers and selecting fueling sites, which would minimize the deviation in power from some desirable reference condition or increasing or decreasing fueling frequency to manage reactivity. An important feature of this analysis was to perform long-transients (1–3 years) starting with each one of the 151 perturbed full core models. It was found that the operational feedback reduced the standard deviation in core reactivity by 99% from 0.0045 to 2.8 × 10−5. Overall, the conclusions demonstrate that while microscopic nuclear data uncertainties may give rise to large macroscopic variability during simple propagation, when important macrolevel feedback are considered the variability is significantly reduced.

2014 ◽  
Vol 2014 ◽  
pp. 1-14
Author(s):  
M. R. Ball ◽  
C. McEwan ◽  
D. R. Novog ◽  
J. C. Luxat

The propagation of nuclear data uncertainties through reactor physics calculation has received attention through the Organization for Economic Cooperation and Development—Nuclear Energy Agency’s Uncertainty Analysis in Modelling (UAM) benchmark. A common strategy for performing lattice physics uncertainty analysis involves starting with nuclear data and covariance matrix which is typically available at infinite dilution. To describe the uncertainty of all multigroup physics parameters—including those at finite dilution—additional calculations must be performed that relate uncertainties in an infinite dilution cross-section to those at the problem dilution. Two potential methods for propagating dilution-related uncertainties were studied in this work. The first assumed a correlation between continuous-energy and multigroup cross-sectional data and uncertainties, which is convenient for direct implementation in lattice physics codes. The second is based on a more rigorous approach involving the Monte Carlo sampling of resonance parameters in evaluated nuclear data using the TALYS software. When applied to a light water fuel cell, the two approaches show significant differences, indicating that the assumption of the first method did not capture the complexity of physics parameter data uncertainties. It was found that the covariance of problem-dilution multigroup parameters for selected neutron cross-sections can vary significantly from their infinite-dilution counterparts.


2010 ◽  
Vol 2 ◽  
pp. 12001 ◽  
Author(s):  
J.N. Wilson ◽  
S. Siem ◽  
S.J. Rose ◽  
A. Georgen ◽  
F. Gunsing ◽  
...  

2021 ◽  
Vol 247 ◽  
pp. 15007
Author(s):  
Liangzhi Cao ◽  
Zhuojie Sui ◽  
Bo Wang ◽  
Chenghui Wan ◽  
Zhouyu Liu

A method of Covariance-Oriented Sample Transformation (COST) has been proposed in our previous work to provide the converged uncertainty analysis results with a minimal sample size. The transient calculation of nuclear reactor is a key part of the reactor-physics simulation, so the accuracy and confidence of the neutron kinetics results have attracted much attention. In this paper, the Uncertainty Quantification (UQ) function of the high fidelity neutronics code NECP-X has been developed based on our home-developed uncertainty analysis code UNICORN, building a platform for the UQ of the transient calculation. Furthermore, the well-known space-time heterogeneous neutron kinetics benchmark C5G7 and its uncertainty propagation from the nuclear data to the interested key parameters of the core have been investigated. To address the problem of “the curse of dimensionality” caused by the large number of input parameters, the COST method has been applied to generate multivariate normal-distribution samples in uncertainty analysis. As a result, the law of the assembly/pin normalized power and their uncertainty with respect to time after introducing an instantaneous perturbation has been obtained. From the numerical results, it can be observed that the maximum relative uncertainties for the assembly normalized power can up to be about 1.65% and the value for the pin-wise power distributions can be about 2.71%.


2021 ◽  
Vol 247 ◽  
pp. 09007
Author(s):  
Isabelle Duhamel ◽  
Nicolas Leclaire ◽  
Luiz Leal ◽  
Atsushi Kimura ◽  
Shoji Nakamura

Available nuclear data for molybdenum included in the nuclear data libraries are not of sufficient quality for reactor physics or criticality safety issues and indeed information about uncertainties and covariance is either missing or leaves much to be desired. Therefore, IRSN and JAEA performed experimental measurements on molybdenum at the J-PARC (Japan Proton Accelerator Research Complex) facility in Japan. The aim was to measure capture cross section and transmission of natural molybdenum at the ANNRI (Accurate Neutron-Nucleus Reaction measurement Instrument) in the MLF (Material Life and science Facility) of J-PARC. The measurements were performed on metallic natural molybdenum samples with various thicknesses. A NaI detector, placed at a flight-path length of about 28 m, was used for capture measurements and a Li-glass detector (flight-path length of about 28.7 m) for transmission measurements. Following the data reduction process, the measured data are being analyzed and evaluated to produce more accurate cross sections and associated uncertainties.


2008 ◽  
Vol 2008 ◽  
pp. 1-5
Author(s):  
I. Kodeli

An experiment on a mockup of the test blanket module based on helium-cooled lithium lead (HCLL) concept will be performed in 2008 in the Frascati Neutron Generator (FNG) in order to study neutronics characteristics of the module and the accuracy of the computational tools. With the objective to prepare and optimise the design of the mockup in the sense to provide maximum information on the state-of-the-art of the cross-section data the mockup was pre-analysed using the deterministic codes for the sensitivity/uncertainty analysis. The neutron fluxes and tritium production rate (TPR), their sensitivity to the underlying basic cross-sections, as well as the corresponding uncertainties were calculated using the deterministic transport codes (DOORS package), the sensitivity/uncertainty code package SUSD3D, and the VITAMINJ/ COVA covariance matrix libraries. The cross-section reactions with largest contribution to the uncertainty of the calculated TPR were identified to be (n,2n) and (n,3n) reactions on lead. The conclusions of this work support the main benchmark design and suggest some modifications and improvements. In particular this study recommends the use, as far as possible, of both natural and enriched lithium pellets for the TRP measurements. The combined use is expected to provide additional and complementary information on the sensitive cross-sections.


2021 ◽  
Vol 247 ◽  
pp. 07003
Author(s):  
A. Sargeni ◽  
E. Ivanov

The paper presents our first results of the exercise III-I-2c from the OECD-NEA UAM-LWR benchmark intended to an elaboration of the methodology of uncertainty propagation. The considered case studied a full PWR core behavior in fast (~0.1 sec) rod ejection transient. According to the benchmark, the core represented a Hot Zero Power state. Authors used brute-force sampling propagating nuclear data and thermo-fluid uncertainties using 3D computational IRSN chain HEMERA. It couples the reactor physics code CRONOS and thermal-hydraulic core code FLICA4. The nuclear data uncertainties were represented in a form of cross sections standard deviations (in percentage of the mean cross sections values) supplied by the UAM team. In addition to the original benchmark, the study includes a case with an increased power peak by supplementary rod ejection, i.e. with higher reactivity. Both the results are similar to what we obtained in the mini-core rod ejection: the power standard deviation follows, in percentage of the mean power, the mean power curve. We split the variance with a direct calculation: once the cross sections are modified and the thermal-hydraulics inputs are kept constant, another time the contrary. The results show that uncertainties dues to nuclear data dominate over ones due to the thermal-flow area. Furthermore, the major contributors in peak-of-power variance lie in a fast group of cross sections.


2021 ◽  
pp. 113-131
Author(s):  
Wei Shen ◽  
Benjamin Rouben

Reactor physics aims to understand accurately the reactivity and the distribution of all the reaction rates (most importantly of the power), and their rate of change in time, for any reactor configuration. To do this, the multiplication factor (or, equivalently, reactivity) and the neutron-flux distribution under various operating conditions and at different times need to be calculated repeatedly. Most of the other parameters of interest (such as neutron reaction rates, power, heat deposition, etc.) are derived from them. They are governed by the geometry, the material composition and the nuclear data (i.e., the neutron cross sections, their energy dependence, the energy spectra and the angular distributions of secondary particles, etc.). For radiation-shielding calculations, additional photon interactions and coupled neutron-photon interaction data are required.


2021 ◽  
Vol 2 (3) ◽  
pp. 281-308 ◽  
Author(s):  
Ruixian Fang ◽  
Dan Gabriel Cacuci

This work extends the investigation of higher-order sensitivity and uncertainty analysis from 3rd-order to 4th-order for a polyethylene-reflected plutonium (PERP) OECD/NEA reactor physics benchmark. Specifically, by applying the 4th-order comprehensive adjoint sensitivity analysis methodology (4th-CASAM) to the PERP benchmark, this work presents the numerical results of the most important 4th-order sensitivities of the benchmark’s total leakage response with respect to the benchmark’s 180 microscopic total cross sections, which includes 180 4th-order unmixed sensitivities and 360 4th-order mixed sensitivities corresponding to the largest 3rd-order ones. The numerical results obtained in this work reveal that the number of 4th-order relative sensitivities that have large values (e.g., greater than 1.0) is far greater than the number of important 1st-, 2nd- and 3rd-order sensitivities. The majority of those large sensitivities involve isotopes 1H and 239Pu contained in the PERP benchmark. Furthermore, it is found that for most groups of isotopes 1H and 239Pu of the PERP benchmark, the values of the 4th-order relative sensitivities are significantly larger than the corresponding 1st-, 2nd- and 3rd-order sensitivities. The overall largest 4th-order relative sensitivity S(4)σt,6g=30,σt,6g=30,σt,6g=30,σt,6g=30=2.720×106 is around 291,000 times, 6350 times and 90 times larger than the corresponding largest 1st-order, 2nd-order and 3rd-order sensitivities, respectively, and the overall largest mixed 4th-order relative sensitivity S(4)σt,630,σt,630,σt,630,σt,530=2.279×105 is also much larger than the largest 2nd-order and 3rd-order mixed sensitivities. The results of the 4th-order sensitivities presented in this work have been independently verified with the results obtained using the well-known finite difference method, as well as with the values of the corresponding symmetric 4th-order sensitivities. The 4th-order sensitivity results obtained in this work will be subsequently used on the 4th-order uncertainty analysis to evaluate their impact on the uncertainties they induce in the PERP leakage response.


2020 ◽  
Vol 2020 ◽  
pp. 1-14
Author(s):  
Ishita Trivedi ◽  
Jason Hou ◽  
Giacomo Grasso ◽  
Kostadin Ivanov ◽  
Fausto Franceschini

In this study, the Best Estimate Plus Uncertainty (BEPU) approach is developed for the systematic quantification and propagation of uncertainties in the modelling and simulation of lead-cooled fast reactors (LFRs) and applied to the demonstration LFR (DLFR) initially investigated by Westinghouse. The impact of nuclear data uncertainties based on ENDF/B-VII.0 covariances is quantified on lattice level using the generalized perturbation theory implemented with the Monte Carlo code Serpent and the deterministic code PERSENT of the Argonne Reactor Computational (ARC) suite. The quantities of interest are the main eigenvalue and selected reactivity coefficients such as Doppler, radial expansion, and fuel/clad/coolant density coefficients. These uncertainties are then propagated through safety analysis, carried out using the MiniSAS code, following the stochastic sampling approach in DAKOTA. An unprotected transient overpower (UTOP) scenario is considered to assess the effect of input uncertainties on safety parameters such as peak fuel and clad temperatures. It is found that in steady state, the multiplication factor shows the most sensitivity to perturbations in 235U fission, 235U ν, and 238U capture cross sections. The uncertainties of 239Pu and 238U capture cross sections become more significant as the fuel is irradiated. The covariance of various reactivity feedback coefficients is constructed by tracing back to common uncertainty contributors (i.e., nuclide-reaction pairs), including 238U inelastic, 238U capture, and 239Pu capture cross sections. It is also observed that nuclear data uncertainty propagates to uncertainty on peak clad and fuel temperatures of 28.5 K and 70.0 K, respectively. Such uncertainties do not impose per se threat to the integrity of the fuel rod; however, they sum to other sources of uncertainties in verifying the compliance of the assumed safety margins, suggesting the developed BEPU method necessary to provide one of the required insights on the impact of uncertainties on core safety characteristics.


2021 ◽  
Vol 247 ◽  
pp. 15021
Author(s):  
Dan G. Cacuci

This invited keynote presentation compares the relative importance of 1st-order versus 2nd-order sensitivities of the leakage response of an OECD/NEA benchmark (polyethylene-reflected plutonium sphere) to the nuclear data characterizing this benchmark. The imprecisely known parameters underlying the neutron transport computational model for this benchmark include 180 group-averaged total microscopic cross sections, 21600 group-averaged scattering microscopic cross sections, 60 parameters describing the fission process, 30 parameters describing the fission spectrum, 10 parameters describing the system’s sources, and 6 isotopic number densities. Thus, this benchmark comprises 21886 1st-order sensitivities of the leakage response with respect to the model parameters, and 478,996,996 2nd-order sensitivities, of which 239,509,441 are distinct. The exact deterministic computation of all of these 1st- and 2nd-order sensitivities was made possible by the application of the Second-Order Adjoint Sensitivity Analysis Methodology (2nd-ASAM) developed by Cacuci. Thousands (out of the 32 400 elements) of the 2nd-order sensitivities of the leakage response with respect to the total cross sections turned out to be significantly larger than the largest corresponding 1st-order sensitivities, contrary to some previously held beliefs in the reactor physics community. Hence, it will be shown that neglecting the 2nd-order sensitivities to total cross sections would cause very large non-conservative errors by under-reporting the response’s variance and expected value. The 2nd-order sensitivities also cause the response distribution to be skewed towards positive values relative to the expected value, which, in turn, is significantly larger than the computed value of the leakage response. The result presented in this paper also underscore the need for obtaining reliable cross section covariance data, which are not available at this time.


Sign in / Sign up

Export Citation Format

Share Document