Efficient use of Monte Carlo: the fast correlation coefficient
Department of Physics and Astronomy, Uppsala University,
2 Nuclear Research and Consultancy Group NRG, Petten, The Netherlands
3 Reactor Physics and Thermal Hydraulic Laboratory, Paul Scherrer Institut, Villigen, Switzerland
* e-mail: email@example.com
Received in final form: 16 February 2018
Accepted: 4 May 2018
Published online: 29 June 2018
Random sampling methods are used for nuclear data (ND) uncertainty propagation, often in combination with the use of Monte Carlo codes (e.g., MCNP). One example is the Total Monte Carlo (TMC) method. The standard way to visualize and interpret ND covariances is by the use of the Pearson correlation coefficient,
where x or y can be any parameter dependent on ND. The spread in the output, σ, has both an ND component, σND, and a statistical component, σstat. The contribution from σstat decreases the value of ρ, and hence it underestimates the impact of the correlation. One way to address this is to minimize σstat by using longer simulation run-times. Alternatively, as proposed here, a so-called fast correlation coefficient is used,
In many cases, cov(xstat; ystat) can be assumed to be zero. The paper explores three examples, a synthetic data study, correlations in the NRG High Flux Reactor spectrum, and the correlations between integral criticality experiments. It is concluded that the use of ρ underestimates the correlation. The impact of the use of ρfast is quantified, and the implication of the results is discussed.
© H. Sjöstrand 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.