A New Direct Calculation Method of Response Matrices Using a Monte Carlo Calculation

A novel direct calculation method of response matrices on heterogeneous lattices by using the Monte Carlo method is proposed. These direct response matrices (DRMs) can be used in core calculations in place of the conventional homogenized lattice constants. The DRMs are formalized by four sub response matrices (sub-RMs) in order to respond to a core eigenvalue, k; thus the DRMs can be re-evaluated on each outer iteration in the core calculations. The sub-RMs can be evaluated by analyzing each neutron's trajectory from ordinary lattice calculations with the Monte Carlo code. Since these sub-RMs are calculated directly under an actual complex assembly geometry, i.e., without a homogenization process, intra-assembly heterogeneous effects can be reflected on global partial current balance calculations. With using two of the sub-RMs, which deal with neutron production probabilities for each fuel pin, and the obtained partial current balance, pin-wise neutron production distributions can also be reconstructed. The DRM method is applied to infinite lattice calculations and multi-assembly calculations. As a result, it is shown that the derivation of the method is valid and intra-assembly heterogeneity is reflected in the results correctly.