Development of Monte Carlo Machine for Particle Transport Problem

Monte Carlo machine, Monte-4 has been developed to realize high performance computing of Monte Carlo codes for particle transport. The calculation for particle tracking in a complex geometry requires (1) classification of particles by the region types using multi-way conditional branches, and (2) determination whether intersections of particle paths with surfaces of the regions are on the boundaries of the regions or not, using nests of conditional branches. How-ever, these procedures require scalar operations or unusual vector operations. Thus the speedup ratios have been low, i.e. nearly two times, in vector processing of Monte Carlo codes for particle transport on conventional vector processors. The Monte Carlo machine Monte-4 has been equipped with the special hardware called Monte Carlo pipelines to process these procedures with high performance. Additionally Monte-4 has been equipped with enhanced load/store pipelines to realize fast transfer of indirectly addressed data for the purpose of resolving imbalances between the performance of data transfers and arithmetic operations in vector processing of Monte Carlo codes on conventional vector processors. Finally, Monte-4 has a parallel processing capability with four processors to multiply the performance of vector processing. We have evaluated the effective performance of Monte-4 using production-level Monte Carlo codes such as vectorized KENO-IV and MCNP. In the performance evaluation, nearly ten times speedup ratios have been obtained, compared with scalar processing of the original codes.     

Neutron Transport Calculations by Using Double-Differential Cross Sections

Some test calculations were carried out to demonstrate the usefulness of double-differential cross sections for neutron transport calculations including anisotropic scattering. A transport code system NITRAN was applied for the purpose. In NITRAN, the anisotropy of elastic and inelastic scattering can be treated in a general form by double-differential total neutron-emission cross sections, which are generated from single-differential and/or original double-differential cross section data base.

The test calculations were performed for neutron flux spectra in aluminum and lead slabs, and also for tritium production rates in a natural lithium sphere. Since the treatment free from collision kinematics is possible by using the double-differential cross sections in the Sⁿcalculations, the discretization of secondary neutron energy distribution becomes independent of the segmentation of angular distribution. A significant improvement due to this independence can be seen in calculating the anisotropy of general inelastic scattering and the extreme anisotropy of elastic scattering by heavy nuclei. For precise anisotropic transport calculations, it is therefore concluded that the nuclear data of double-differential type are more suitable than those of single-differential type.     

Cross Section Assessment by Measurement of Time Dependent Neutron Spectra from Lithium Assemblies

Time dependent neutron spectra from lithium assemblies were measured to assess the neutron cross sections of ⁷Li in ENDF/B-IV, which is important nuclide for the D-T fusion reactor blanket material. Pulsed neutrons produced by D-D or D-T reaction were used to measure leakage neutron spectra from cubical lithium assemblies as a function of time by the use of NE213 liquid scintillator. Calculations of time dependent neutron spectra were carried out by the Monte Carlo code SIMON, which was prepared for this study. The group constants used in these calculations were processed from ENDF/B-IV data. The calculated and the measured neutron spectra were compared for the following three; a stationary spectrum, spectra at each time interval and decay curves for specified energy groups. Discrepancies between the measured and the calculated neutron spectra were found in these comparisons. In order to assure the cause of these discrepancies, some calculations were carried out with recently measured cross sections of inelastic scattering which excite 0.478 and 4.63 MeV level of ⁷Li. It was concluded that some of the neutron cross section data of ⁷Li in ENDF/B-IV should be ameliorated.     

Analysis of Critical Experiments Using Medium-Enriched-Uranium Fuel in Kyoto University Critical Assembly (KUCA)

The critical experiments using medium-enriched-uranium (MEU) fuel in the Kyoto University Critical Assembly (KUCA), a light-water-moderated and heavy-water-reflected cylindrical core, were started in May 1981, as a part of the international Reduced Enrichment for Research and Test Reactors (RERTR) program.

The following KUCA critical experiments were analyzed: (1) the criticality measurements for high-enriched-uranium (HEU) and MEU cores and (2) the reactivity effect measurements of boron burnable-poison (BP) for MEU cores. Five-group constants were generated using the EPRI-CELL code, and two-dimensional diffusion calculations were performed using a conventional finite-difference code DIF3D(2D), and a finite-element code 2D-FEM-KUR. Some of the results from the two diffusion codes were compared with each other. Advantage was taken of the finite-element method for the application of the 2D-FEM-KUR code to a detailed analysis of the BP effect measurements.

Differences between the results of calculations and experiments were less than 1.8% in C/E ratios for eigenvalues. The agreement between the results obtained using the DIF3D(2D) code and the 2D-FEM-KUR code was excellent. The calculated results of the BP effects with use of the 2D-FEM-KUR code approximately agreed with the experiments.     

Evaluation of Critical Bucklings of Light-Water Moderated Low Enriched UO2 Cores by the Variable Loading Method

Abstract

Horizontal extrapolated distances Xh and critical bucklings B² _C for light-water moderated low enriched UO² cores were evaluated using the variable loading method. A set of critical height (water levels) of the cores having various horizontal dimensions was used in this method. Errors of the present method were smaller than the previous ones by the flux shape method in which the horizontal extrapolated distance was determined from the horizontal power distribution. The presently obtained Xh and B² _C reproduced well the change in the critical water level for the whole range of the horizontal core size, while the previous ones could be used only for a limited range. Moreover, the effective multiplication factors obtained through the cell calculation of the SRAC code by using the presently evaluated critical bucklings agreed much better with the rigorous ones by the continuous energy Monte Carlo calculation using a full core model, than those by using the previous values. Therefore, it was concluded that the variable loading method can give more accurate values for Xh and B² _C than the flux shape method for the cores investigated.