DOI QR코드

DOI QR Code

Performing linear regression with responses calculated using Monte Carlo transport codes

  • Price, Dean (Department of Nuclear Engineering and Radiological Science, University of Michigan) ;
  • Kochunas, Brendan (Department of Nuclear Engineering and Radiological Science, University of Michigan)
  • Received : 2021.08.18
  • Accepted : 2021.11.03
  • Published : 2022.05.25

Abstract

In many of the complex systems modeled in the field of nuclear engineering, it is often useful to use linear regression-based analyses to analyze relationships between model parameters and responses of interests. In cases where the response of interest is calculated by a simulation which uses Monte Carlo methods, there will be some uncertainty in the responses. Further, the reduction of this uncertainty increases the time necessary to run each calculation. This paper presents some discussion on how the Monte Carlo error in the response of interest influences the error in computed linear regression coefficients. A mathematical justification is given that shows that when performing linear regression in these scenarios, the error in regression coefficients can be largely independent of the Monte Carlo error in each individual calculation. This condition is only true if the total number of calculations are scaled to have a constant total time, or amount of work, for all calculations. An application with a simple pin cell model is used to demonstrate these observations in a practical problem.

Keywords

Acknowledgement

We would like to thank Professor Brian Kiedrowski and David Griesheimer for their helpful thoughts during our discussions concerning the subject of this work. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship under Grant No. DGE 1256260.

References

  1. K. Dalbey, M.S. Eldred, G. Geraci, J.D. Jakeman, K.A. Maupin, J.A. Monschke, D.T. Seidl, L.P. Swiler, A. Tran, F. Menhorn, et al., Dakota, a Multilevel Parallel Object-Oriented Framework for Design Optimization, Parameter Estimation, Uncertainty Quantification, and Sensitivity Analysis: Version 6.12theory Manual, SAND2020-4987, Sandia National Laboratory, 2020.
  2. A. Alfonsi, C. Wang, P. Talbot, M.M. Abdo, M. Gamal, D. Mandelli, C. Rabiti, J.J. Cogliati, R.A. Kinoshita, Raven User Manual, Tech. rep., Idaho National Lab.(INL), 2021. Idaho Falls, ID (United States).
  3. M. Reyes-Fuentes, E. del Valle-Gallegos, J. Duran-Gonzalez, J. Or t iz-Villafuerte, R. Castillo-Dur an, A.G. omez-Torres, C. Queral, Aztusia: a new application software for un certainty and sensitivity analysis for nuclear reactors, Reliab. Eng. Syst. Saf. 209 (2021) 107441. https://doi.org/10.1016/j.ress.2021.107441
  4. C. Yu, W. Yao, Robust linear regression: a review and comparison, Commun. Stat. Simulat. Comput. 46 (8) (2017) 6261-6282. https://doi.org/10.1080/03610918.2016.1202271
  5. W. Zwermann, N. Berner, A. Aures, K. Velkov, Sensitivity and uncertainty analysis for the uam-sfrsub-exercises with linear regression from random sampling, Ann. Nucl. Energy 149 (2020) 107832. https://doi.org/10.1016/j.anucene.2020.107832
  6. N.R. Draper, R.C. Van Nostrand, Ridge regression and james-stein estimation: review and comments, Technometrics 21 (4) (1979) 451-466. https://doi.org/10.1080/00401706.1979.10489815
  7. R. Tibshirani, Regression shrinkage and selection via the lasso: a retrospective, J. Roy. Stat. Soc. B 73 (3) (2011) 273-282. https://doi.org/10.1111/j.1467-9868.2011.00771.x
  8. D. Price, A. Maile, J. Peterson-Droogh, D. Blight, A Methodology for Uncertainty Quantification and Sensitivity Analysis for Responses Subject to Monte Carlo Uncertainty with Application to Fuel Plate Char-Acteristics in the Atrc, Nuclear Engineering and Technology.
  9. G.A. Seber, A.J. Lee, Linear Regression Analysis, vol. 329, John Wiley & Sons, 2012.
  10. Y. Dodge, The Concise Encyclopedia of Statistics, Springer New York, New York, NY, 2008.
  11. F.B. Brown, A Review of Best Practices for Monte Carlo Criticality Calculations, Los Alamos National Lab.(LANL), Los Alamos, NM (United States), 2009.
  12. L. Wasserman, All of Statistics: a Concise Course in Statistical Inference, Springer Science & Business Media, 2013.
  13. J. Leppanen, M. Pusa, T. Viitanen, V. Valtavirta, T. Kaltiaisenaho, The serpent Monte Carlo code: status, development and applications in 2013, Ann. Nucl. Energy 82 (2015) 142-150, joint International Conference on Supercomputing in Nuclear Applications and Monte Carlo 2013, SNA + MC 2013. Pluri-and Trans-disciplinarity, Towards New Modeling and Numerical Simulation Paradigms. https://doi.org/10.1016/j.anucene.2014.08.024