Fatemeh MOHAMMADHASANİ, Ahmad PİROUZMAND

Multi-Group Neutron Cross Sections and Scattering Matrix Calculations Via Monte Carlo Method

This study describes a new methodology to estimate multi-group neutron cross sections and scattering matrix elements from a Monte Carlo simulation, particularly from MCNPX 2.7 code. The geometric flexibility associated with the Monte Carlo method makes it very suitable for the generation of highly accurate multi-group constants. While the deterministic lattice codes are not capable of dealing with energy as a continuous variable the Monte Carlo codes such as MCNPX make use of a continuous energy cross sections data for neutron transport calculations. To determine the scattering matrix, an output file of a MCNPX run (socalled PTRAC file) with all relevant source, collision and terminal events of the simulation is used. First, by a separately special program in MATLAB, the PTRAC file is read and tracks are selected in the geometrical region for which one wants to calculate the multi-group constants. Then, information such as coordinates of the particle position, particle direction with coordinates axes, energy and weight of the particle are extracted. The multi-group scattering matrix elements are generated via the weight-to-flux ratio method using the above data available in the PTRAC file. The proposed method is validated using three benchmark problems involve a slab, a pin cell, and a fuel assembly of Tehran research reactor (TRR). The generated multi-group constants via presented method and multiplication factor calculation presents a good comparison to the reference values.

PDF

___

Duderstadt, J.J. and Hamilton, L.J., "Nuclear Reactor Analysis", 1st Edition, John Wiley & Sons. (1976)

Hoogenboom, J.H., "Methodology of continuous-energy adjoint Monte Carlo for neutron, photon, and coupled neutron-photon transport", Nuclear Science and Engineering. 143 (2), 99-120, (2003)

Hoogenboom J.E, Khotylev V.A, and Tholammakkil J. M, "Generation of Multi-Group Cross Sections and Scattering Matrices with the Monte Carlo Code MCNP5", Joint International Topical Meeting on Mathematics & Computation and Supercomputing in Nuclear Applications, 15-19 April 2007, California, American Nuclear Society.

Pounders, J. M., "Stochastically Generated Multigroup Diffusion Coefficients", MSc. Thesis, Georgia Institute of Technology (2006).

Thiyagarajan T. K., Ganesan S, Jagannathan V, and Karthikeyan R, "A Computer Program with Graphical User Interface to Plot the Multigroup Cross Sections of WIMS-D Library", Annals of Nuclear Energy. 29 (14),1735-1745, (2002)

Thiyagarajan T. K., Ganesan S, Jagannathan V, and Karthikeyan R, "Generation of Handbook of Multi-group Cross Sections of WIMS-D Libraries by Using the XnWlup2.0", Annals of Nuclear Energy, 30 (7), 793-809, (2003)

Rantamäki K. M., Häkkinen S, and Anttila M, "Impact of Cross Section Libraries on CASMO-4E and CASMO-4 Results", Annals of Nuclear Energy. 38 (12), 2830-2835,(2011)

Redmond, E.L., "Multigroup Cross Section Generation via Monte Carlo Methods", Ph.D. Thesis, MIT University.(1997)

Yoshioka K, and Ando Y, "Multigroup scattering matrix generation method using weight-to-flux ratio based on a continuous energy Monte Carlo technique" J. nuclear science and technology. 47 (10), 908-916, (2010)

Duderstadt J.J, and Hamilton L.J, "Nuclear Reactor Analysis", 1st Edition, Wiley (1978)

Yamamoto T.Y, and Nomura Y, "OECD/NEA Source Convergence Benchmark 3: Three thick one-dimensional slabs", Japan Atomic Energy Research Institute, Japan.

Atomic Energy Organization of Iran, "Safety Analysis Report for Tehran Research Reactor (TRR)", Tehran (2006)