Browse > Article
http://dx.doi.org/10.1016/j.net.2017.02.008

Stabilization effect of fission source in coupled Monte Carlo simulations  

Olsen, Borge (Division of Nuclear Reactor Technology, KTH Royal Institute of Technology, AlbaNova University Center)
Dufek, Jan (Division of Nuclear Reactor Technology, KTH Royal Institute of Technology, AlbaNova University Center)
Publication Information
Nuclear Engineering and Technology / v.49, no.5, 2017 , pp. 1095-1099 More about this Journal
Abstract
A fission source can act as a stabilization element in coupled Monte Carlo simulations. We have observed this while studying numerical instabilities in nonlinear steady-state simulations performed by a Monte Carlo criticality solver that is coupled to a xenon feedback solver via fixed-point iteration. While fixed-point iteration is known to be numerically unstable for some problems, resulting in large spatial oscillations of the neutron flux distribution, we show that it is possible to stabilize it by reducing the number of Monte Carlo criticality cycles simulated within each iteration step. While global convergence is ensured, development of any possible numerical instability is prevented by not allowing the fission source to converge fully within a single iteration step, which is achieved by setting a small number of criticality cycles per iteration step. Moreover, under these conditions, the fission source may converge even faster than in criticality calculations with no feedback, as we demonstrate in our numerical test simulations.
Keywords
Monte Carlo; Coupled simulation; Fixed-point iteration; Numerical stability; Feedback;
Citations & Related Records
연도 인용수 순위
  • Reference
1 M. Mori, W. Maschek, E. Laurien, K. Morita, Monte-Carlo/Simmer-III Reactivity Coefficients Calculations for the Supercritical Water Fast Reactor, Proceedings of the ANS/ENS Topical Meeting GLOBAL, New Orleans, 2003.
2 H.G. Joo, J.Y. Cho, K.Y. Kim, M.H. Chang, B.S. Han, C.H. Kim, Consistent Comparison of Monte Carlo and Whole-core Transport Solutions for Cores with Thermal Feedback. PHYSOR 2004-the Physics of Fuel Cycles and Advanced Nuclear Systems, Global Developments, Chicago, IL, April 2004.
3 R.L. Moore, B.G. Schnitzler, C.A. Wemple, MOCUP: MCNP-ORIGEN2 coupled utility program, Technical Report INEL-95/0523, 130667, September 1995.
4 D.P. Griesheimer, D.F. Gill, B.R. Nease, T.M. Sutton, M.H. Stedry, P.S. Dobreff, D.C. Carpenter, T.H. Trumbull, E. Caro, H. Joo, D.L. Millman, MC21 v.6.0-a continuous-energy Monte Carlo particle transport code with integrated reactor feedback capabilities, Ann. Nucl. Energy 82 (2015) 29-40.   DOI
5 J. Leppanen, Serpent-a Continuous-energy Monte Carlo Reactor Physics Burnup Calculation Code, Technical Report, VTT Technical Research Centre of Finland, August 16, 2012.
6 K. Wang, Z. Li, D. She, J. Liang, Q. Xu, Y. Qiu, J. Yu, J. Sun, X. Fan, G. Yu, RMC-a Monte Carlo code for reactor core analysis, Ann. Nucl. Energy 82 (2015) 121-129.   DOI
7 B.L. Sjenitzer, J.E. Hoogenboom, J. Jimnez Escalante, V. Sanchez Espinoza, Coupling of dynamic Monte Carlo with thermal-hydraulic feedback, Ann. Nucl. Energy 76 (2015) 27-39.   DOI
8 A.E. Isotalo, J. Leppanen, J. Dufek, Preventing xenon oscillations in Monte Carlo burnup calculations by enforcing equilibrium xenon distribution, Ann. Nucl. Energy 60 (2013) 78-85.   DOI
9 J. Dufek, D. Kotlyar, E. Shwageraus, The stochastic implicit Euler method-a stable coupling scheme for Monte Carlo burnup calculations,, Ann. Nucl. Energy 60 (2013) 295-300.   DOI
10 J. Dufek, H. Anglart, Derivation of a stable coupling scheme for Monte Carlo burnup calculations with the thermal-hydraulic feedback, Ann. Nucl. Energy 62 (2013) 260-263.   DOI
11 Roots of nonlinear equations, in: J.R. Hauser (Ed.), Numerical Methods for Nonlinear Engineering Models, Springer, Netherlands, 2009, pp. 43-76.
12 C. Kelley, Iterative Methods for Linear and Nonlinear Equations, Frontiers in Applied Mathematics, Society for Industrial and Applied Mathematics, Philadelphia, January 1995.
13 D.F. Griffiths, D.J. Higham, Linear multistep methods-V: solving implicit methods, in: Numerical Methods for Ordinary Differential Equations, Springer London, London, 2010, pp. 109-121.
14 G. Alefeld, J. Herzberger, Fixed point iteration for nonlinear systems of equations, in: W. Rheinboldt (Ed.), Introduction to Interval Computation, Academic Press, San Diego, 1983, pp. 131-142.
15 J. Dufek, W. Gudowski, Stochastic approximation for Monte Carlo calculation of steady-state conditions in thermal reactors, Nucl. Sci. Eng. 152 (2006) 274-283.   DOI
16 E.M. Gelbard, A.G. Gu, Biases in Monte Carlo eigenvalue calculations, Nucl. Sci. Eng. 117 (1994) 1-9.   DOI
17 J. Dufek, J.E. Hoogenboom, Description of a stable scheme for steady-state coupled Monte Carlo-thermal-hydraulic calculations, Ann. Nucl. Energy 68 (2014) 1-3.   DOI
18 H. Etherington, Nuclear Engineering Handbook, first ed., McGraw-Hill, New York, 1958.
19 U.S. Department of Energy, Technical Report DOE-HDBK-1019/2-93, DOE Fundamentals Handbook-Nuclear Physics and Reactor Theory, Vol. 2 of 2, U.S. Department of Energy, Washington, D.C., January 1993.
20 H.J. Shim, C.H. Kim, Elimination of variance bias in Monte Carlo eigenvalue calculations, Trans. Am. Nucl. Soc. 97 (2007) 653-655.
21 H.J. Shim, C.H. Kim, Real variance estimation using an inter-cycle fission source correlation for Monte Carlo eigenvalue calculations, Nucl. Sci. Eng. 162 (2009) 98-108.   DOI
22 D.F. Gill, D.P. Griesheimer, D.L. Aumiller, Numerical methods in coupled Monte Carlo and thermal-hydraulic calculations, Nucl. Sci. Eng. 185 (2017) 194-205.   DOI