Benchmark of WIMS-IST Against MCNP for CANDU Pressure Tube Fast Fluxes
Main Article Content
Abstract
Pressure tube fast-flux data in CANDU are currently calculated using the multi-group neutron transport code WIMS-IST. In this study, the WIMS-IST fast flux calculations are benchmarked against MCNP calculations (a Monte Carlo particle transport code), over the range of fuel burnup and coolant density in CANDU. The comparison shows good agreement between WIMS and MCNP, with WIMS fast fluxes being 1.5% to 4% lower than the MCNP values. The difference is smallest for fresh fuel, and increases with burnup. The fast flux gradient across the pressure tube (factor of 1.23 from inner edge to outer edge) is accurately calculated by WIMS.
When reporting fast fluxes in pressure tubes, these are generally given as >1.000 MeV fluxes. For WIMS, this requires an extra conversion step, since the WIMS ENDF/B libraries do not have a group boundary at 1 MeV. The conversion step is based on a fictitious isotope ONEMEV in the WIMS nuclear data library. The conversion factor in WIMS was found to be about one percent too high. When providing >1 MeV fluxes from WIMS, this partially compensates for the slight under prediction of the fast flux. Pressure tube >1 MeV fluxes from WIMS are therefore 0.5% to 3% lower than MCNP values.
To obtain accurate fast flux data, neutron transport calculations must be performed on a critical cell. For this study, all calculations were performed with radial albedo boundary conditions giving a critical cell. This required the use of an albedo version of MCNP, developed at AECL.
When reporting fast fluxes in pressure tubes, these are generally given as >1.000 MeV fluxes. For WIMS, this requires an extra conversion step, since the WIMS ENDF/B libraries do not have a group boundary at 1 MeV. The conversion step is based on a fictitious isotope ONEMEV in the WIMS nuclear data library. The conversion factor in WIMS was found to be about one percent too high. When providing >1 MeV fluxes from WIMS, this partially compensates for the slight under prediction of the fast flux. Pressure tube >1 MeV fluxes from WIMS are therefore 0.5% to 3% lower than MCNP values.
To obtain accurate fast flux data, neutron transport calculations must be performed on a critical cell. For this study, all calculations were performed with radial albedo boundary conditions giving a critical cell. This required the use of an albedo version of MCNP, developed at AECL.
Article Details
Section
Articles