Particle clustering in Monte Carlo criticality simulations

Abstract Strong space and cycle correlations affecting the neutron flux have been recently reported to occur in Monte Carlo criticality simulations of pressurized water reactors pin-cells with reflective boundary conditions: when the system size is large, neutrons are observed to gather together and form ‘clusters’ that wander around in space over many cycles. In most cases, the outcome of such clustering is that simulations display wild fluctuations, largely exceeding those expected around the equilibrium distribution. In this paper, we analyze the reasons behind the emergence of these phenomena: we show that the key mechanism is due to the asymmetry between neutron disappearance being uniformly distributed along the pin-cell, and neutron generation being localized at fission sites (i.e., only next to a parent particle). An explanation is provided by resorting to a simplified Brownian transport model coupled with a Galton–Watson birth and death process, which is sufficient to retain the essential features observed in realistic Monte Carlo simulations. An empirical space correlation function is proposed as a diagnostic tool for detecting clustering in criticality simulations.

[1]  Andrea Zoia,et al.  Branching exponential flights: travelled lengths and collision statistics , 2012, 1207.1877.

[2]  G. I. Bell,et al.  Nuclear Reactor Theory , 1952 .

[3]  J. T. Cox,et al.  Occupation Times for Critical Branching Brownian Motions , 1985 .

[4]  Thomas J. Higgins,et al.  Random Processes in Nuclear Reactors , 1975, IEEE Transactions on Nuclear Science.

[5]  Forrest B. Brown,et al.  Time Series Analysis of Monte Carlo Fission Sources—I: Dominance Ratio Computation , 2004 .

[6]  S. Yip,et al.  The Foundations of Neutron Transport Theory , 1963 .

[7]  W. R. Young,et al.  Reproductive pair correlations and the clustering of organisms , 2001, Nature.

[8]  Imre Pázsit,et al.  Neutron Fluctuations: A Treatise on the Physics of Branching Processes , 2007 .

[9]  Mikko Alava,et al.  Branching Processes , 2009, Encyclopedia of Complexity and Systems Science.

[10]  T. Courau,et al.  Dominance Ratio Assessment and Monte Carlo Criticality Simulations: Dealing with High Dominance Ratio Systems , 2010 .

[11]  D. Dawson.  The critical measure diffusion process , 1977 .

[12]  B. Houchmandzadeh,et al.  Neutral clustering in a simple experimental ecological community. , 2008, Physical review letters.

[13]  J. Duderstadt,et al.  Nuclear reactor analysis , 1976 .

[14]  B. Houchmandzadeh.  Clustering of diffusing organisms. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[15]  B. Houchmandzadeh.  Theory of neutral clustering for growing populations. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[16]  Ely M. Gelbard,et al.  Computation of standard deviations in Eigenvalue calculations , 1990 .

[17]  T. Ueki.  Monte Carlo criticality calculation under extreme condition , 2012 .

[18]  A. Mazzolo,et al.  Properties of branching exponential flights in bounded domains , 2012, 1207.6557.

[19]  B. L. Sjenitzer,et al.  Variance reduction for fixed-source Monte Carlo calculations in multiplying systems by improving chain-length statistics , 2011 .

[20]  Yi-Kang Lee,et al.  Source convergence diagnostics using Boltzmann entropy criterion application to different OECD/NEA criticality benchmarks with the 3-D Monte Carlo code Tripoli-4 , 2006 .

[21]  Forrest B. Brown,et al.  Autocorrelation and Dominance Ratio in Monte Carlo Criticality Calculations , 2003 .