The Monte Carlo code MCNPX has been utilized to calculate the effective delayed neutron fraction and reactivity by using the area-ratio method. The effective delayed neutron fractioneff has been calculated with the fission probability method proposed by Meulekamp and van der Marck. MCNPX was used to calculate separately the fission probability of the delayed and the prompt neutrons by using the TALLYX user subroutine of MCNPX. In this way, �eff was obtained from the one criticality (k-code) calculation without performing an adjoint calculation. The traditional k-ratio method requires two criticality calculations to calculateeff, while this approach utilizes only one MCNPX criticality calculation. Therefore, the approach described here is referred to as a one-run method. In subcritical systems driven by a pulsed neutron source, the area-ratio method is used to calculate reactivity (in dollar units) as the ratio between the prompt and delayed areas. These areas represent the integral of the reaction rates induced from the prompt and delayed neutrons during the pulse period. Traditionally, application of the area-ratio method requires two separate fixed- source MCNPX simulations: one with delayed neutrons and the other without. The number of source particles in these two simulations must be extremely high in order to obtain accurate results with low statistical errors because the values of the total and prompt areas are very close. Consequently, this approach is time consuming and suffers from the statistical errors of the two simulations. The present paper introduces a more efficient method for estimating the reactivity calculated with the area method by taking advantage of the TALLYX user subroutine of MCNPX. This subroutine has been developed for separately scoring the reaction rates caused by the delayed and the prompt neutrons during a single simulation. Therefore the method is referred to as a one- run calculation. These methodologies have been utilized for several published benchmark cases, and a comparison of the results shows good agreement.
[1]
N. G. Sjostrand.
MEASUREMENTS ON A SUBCRITICAL REACTOR USING A PULSED NEUTRON SOURCE
,
1956
.
[2]
G. R. Keepin,et al.
Physics of Nuclear Kinetics
,
1967
.
[3]
J. F. Briesmeister.
MCNP-A General Monte Carlo N-Particle Transport Code
,
1993
.
[4]
M. M. Bretscher.
Evaluation of reactor kinetic parameters without the need for perturbation codes.
,
1998
.
[5]
Viktoriya V. Kulik,et al.
Space-Time Correction for Reactivity Determination in Source-Driven Subcritical Systems
,
2006
.
[6]
R. K. Meulekamp,et al.
Calculating the Effective Delayed Neutron Fraction with Monte Carlo
,
2006
.
[7]
M. Iqbal,et al.
Flow of kinetic parameters in a typical swimming pool type research reactor
,
2008
.
[8]
Ted Taylor.
Los Alamos National Laboratory
,
2005
.
[9]
G. Chiba,et al.
Calculation of Effective Delayed Neutron Fraction Using a Modified k-Ratio Method
,
2009
.
[10]
Alberto Talamo,et al.
Deterministic and Monte Carlo modeling and analyses of YALINA-THERMAL subcritical assembly.
,
2010
.
[11]
J. B. Briggs,et al.
Identification of Integral Benchmarks for Nuclear Data Testing Using DICE (Database for the International Handbook of Evaluated Criticality Safety Benchmark Experiments)
,
2011
.