Continuous-energy Monte Carlo Neutron Research Articles

Neutronic and thermal-hydraulic calculations under steady state and transient conditions for a generic pebble bed fluoride salt cooled high-temperature reactor

Here, we develop multiphysics coupled neutronics/thermal-hydraulics models for a recently available generic Fluoride-salt-cooled High-temperature Reactor (FHR) at equilibrium core burnup conditions with reactivity control mechanisms. We utilized the pebble bed reactor transient tool KP-AGREE with multigroup homogenized cross-sections generated by the continuous energy Monte Carlo neutron transport code Serpent to simulate transient events for an FHR. Comparing global reactor parameters and safety characteristics between Serpent and published results to standalone neutronic calculations in KP-AGREE verifies our initial conditions for our transient models. This work explores three separate transients: a loss of flow accident, an overcooling accident, and a total control rod withdrawal accident. The boundary conditions of these scenarios come from either transient benchmark definitions of a gas-cooled pebble bed reactor or transient analysis of a similar FHR design. For the loss of flow accident, our predicted peak fuel and coolant outlet temperatures show a similar response to published values, albeit slightly higher, with differences between values attributed to our modeling assumptions. The overcooling transient simulation predicts a 52 K temperature increase in the maximum fuel kernel temperature and an approximately 32.6% increase in total core power due to the negative coolant temperature feedback of the FHR design. Our control rod removal transients predicted a 235 K increase in maximum fuel kernel temperature and a 107.5% increase in total core power. For all transient cases, the fuel temperature remained significantly lower than the temperature limit, but coolant outlet temperatures for the reactivity-initiated accident exceeded external core material temperature limits.

Domain decomposition in the GPU-accelerated Shift Monte Carlo code

The GPU solver within the Shift continuous-energy Monte Carlo neutron transport code has been extended to provide domain decomposition in addition to domain replication to enable the solution of problems with memory requirements exceeding the capacity of a single GPU. The strategy follows the Multiple Set, Overlapping Domain (MSOD) approach that is used in Shift’s CPU solver and integrates into the event-based algorithm used for Shift’s GPU solver. Furthermore, the ability to assign processors to spatial domains non-uniformly has been maintained. Two different approaches for communicating particle data between domains are considered, and multiple criteria for load balancing problems have been investigated. Numerical results are presented for both fresh and depleted small modular nuclear reactor (SMR) cores. A parallel efficiency of approximately 80% was achieved with up to 16 spatial domains measured relative to full domain replication. A scaling study on the Summit supercomputer demonstrates a weak scaling parallel efficiency of over 90% on over 24000 GPUs.

Extension of Monte Carlo code MCS to spent fuel cask shielding analysis

The shielding analysis of the spent fuel dry storage cask TN-32 is carried out using the continuous-energy Monte Carlo neutron- and photon-transport code MCS developed by the Computational Reactor Physics and Experiment laboratory of the Ulsan National Institute of Science and Technology. The shielding analysis involves the transport simulation of the decay neutrons and gammas emitted by the spent fuel to obtain the neutron and gamma dose rate values outside of the cask. This deep-penetration problem is solved by means of weight window population control as variance reduction. The TN-32 cask calculations are repeated with the MCNP6 Monte Carlo code for verification purpose and a good agreement within three standard deviations is observed between MCS and MCNP6. The neutron dose rate on the cask surface and 2 m away from the cask can be calculated with MCS, whereas the photon dose rate calculation for positions outside of the cask will be further improved by employing more sophisticated variance reduction techniques.

Validation of Unresolved Neutron Resonance Parameters Using a Thick-Sample Transmission Measurement

Often discrepancies can be found in the corresponding cross sections of different evaluated nuclear data libraries. Traditional integral benchmarks that are used to validate such libraries are sensitive to cross-section values across many different energies. This means an erroneously low cross section at one energy may compensate for an erroneously high cross section at another energy, and the integral benchmark value may still be met. While the evaluated cross section may agree with that single benchmark, it could affect other systems differently. To reduce the potential for this error, an energy differential validation method is proposed herein for continuous energy Monte Carlo neutron transport models in the resolved resonance region and the unresolved resonance region (URR). The proposed method exposes the underlying physics of the URR and validates both the average cross section and resonance self-shielding effect driven by the fluctuations in that cross section. This is done by measuring the neutron transmission of a thick sample that, by its nature, exaggerates the resonance self-shielding effect. This validation method is shown to be very sensitive to the cross-section model used (resolved versus unresolved) and the fluctuation correction employed, allowing it to probe the validity of the previously mentioned cross-section evaluations. Tantalum-181 is used as an example to demonstrate the impact of different resonance evaluations. It was found that the JEFF-3.3 and JENDL-4.0u evaluations made reasonable choices for cross-section models of 181Ta; none of the current evaluations, however, can be used to properly model the validation transmission over all energies. It was also found that updating resonance parameters in the URR provided better agreement with the validation transmission.

Continuous-energy Monte Carlo neutron transport on GPUs in the Shift code

A continuous-energy Monte Carlo neutron transport solver executing on GPUs has been developed within the Shift code. Several algorithmic approaches are considered, including both history-based and event-based implementations. Unlike in previous work involving multigroup Monte Carlo transport, it is demonstrated that event-based algorithms significantly outperform a history-based approach for continuous-energy transport as a result of increased device occupancy and reduced thread divergence. Numerical results are presented for detailed full-core models of a small modular reactor (SMR), including a model containing depleted fuel materials. These results demonstrate the substantial gains in performance that are possible with the latest-generation of GPUs. On the depleted SMR core configuration, an NVIDIA P100 GPU with 56 streaming multiprocessors provides performance equivalent to 90 CPU cores, and the latest V100 GPU with 80 multiprocessors offers the performance of more than 150 CPU cores.

Physics-oriented optimization strategy for the energy lookup algorithm in continuous energy Monte Carlo neutron transport simulation

The continuous energy Monte Carlo method is a most high-fidelity and high-resolution method for neutron transport simulations in reactor physics with minimal approximations. However, one of the major disadvantages is that it is very time-consuming and computational-intensive for large scale whole-core simulations, especially for coupling with depletion analysis for realistic reactors. Some recent researches indicate that one of the principle performance bottlenecks for the problem lies in the energy lookup algorithm during the calculation of energy-dependent material cross sections. Therefore, two physics-oriented optimization strategies, based on making use of the physical characteristics of neutron transport behaviors, are developed to optimize the run-time performance of the algorithm for accelerating the energy lookup without any loss in precision and accuracy. The first optimization strategy is called Neighbored Material Cascade Grid (NMCG) which is a hybrid approach utilizing the key features of the cascade grid and double indexing method. The second optimization strategy is called Adaptive Optimal Logarithmic Grid (AOLG) which is a variation of the conventional logarithmic energy grid method utilizing the advantages of energy hash tables. The strategies are incorporated into a continuous energy Monte Carlo neutron transport code and tested on realistic whole-core reactor systems. The computational performance as measured by memory usage, elapsed runtime and overall speedup, associated with each of the optimization strategies are demonstrated in the whole-core Monte Carlo simulations. Depending on the complexity of the models, the number of nuclides in the material compositions and the utilization of different optimization strategies, overall speedup ratios of 1.2-1.7, relative to the conventional binary lookup algorithm, are routinely observed. Furthermore, the numerical results indicate that the run-time performance of the new physics-oriented optimization strategies performs a bit better than that of conventional optimization methods with existing approaches.

MC21/CTF and VERA multiphysics solutions to VERA core physics benchmark progression problems 6 and 7

Abstract The continuous energy Monte Carlo neutron transport code, MC21, was coupled to the CTF subchannel thermal-hydraulics code using a combination of Consortium for Advanced Simulation of Light Water Reactors (CASL) tools and in-house Python scripts. An MC21/CTF solution for VERA Core Physics Benchmark Progression Problem 6 demonstrated good agreement with MC21/COBRA-IE and VERA solutions. The MC21/CTF solution for VERA Core Physics Benchmark Progression Problem 7, Watts Bar Unit 1 at beginning of cycle hot full power equilibrium xenon conditions, is the first published coupled Monte Carlo neutronics/subchannel T-H solution for this problem. MC21/CTF predicted a critical boron concentration of 854.5 ppm, yielding a critical eigenvalue of 0.99994 ± 6.8E-6 (95% confidence interval). Excellent agreement with a VERA solution of Problem 7 was also demonstrated for integral and local power and temperature parameters.

On-the-fly Doppler broadening of unresolved resonance region cross sections

Two methods for computing temperature-dependent unresolved resonance region cross sections on-the-fly within continuous-energy Monte Carlo neutron transport simulations are presented. The first method calculates Doppler broadened cross sections directly from zero-temperature average resonance parameters. In a simulation, at each event that requires cross section values, a realization of unresolved resonance parameters is generated about the desired energy and temperature-dependent single-level Breit-Wigner resonance cross sections are computed directly via the analytical ψ−χ Doppler integrals. The second method relies on the generation of equiprobable cross section magnitude bands on an energy-temperature mesh. Within a simulation, the bands are sampled and interpolated in energy and temperature to obtain cross section values on-the-fly. Both of the methods, as well as their underlying calculation procedures, are verified numerically in extensive code-to-code comparisons. Energy-dependent pointwise cross sections calculated with the newly-implemented procedures are shown to be in excellent agreement with those calculated by a widely-used nuclear data processing code. Relative differences at or below 0.1% are observed. Integral criticality benchmark results computed with the proposed methods are shown to reproduce those computed with a state-of-the-art processed nuclear data library very well. In simulations of fast spectrum systems which are highly-sensitive to the representation of cross section data in the unresolved region, k-eigenvalue and neutron flux spectra differences of <10 pcm and <1.0% are observed, respectively. The direct method is demonstrated to be well-suited to the calculation of reference solutions — against which results obtained with a discretized representation may be assessed — as a result of its treatment of the energy, temperature, and cross section magnitude variables as continuous. Also, because there is no pre-processed data to store (only temperature-independent average resonance parameters) the direct method is very memory-efficient. Typically, only a few kB of memory are needed to store all required unresolved region data for a single nuclide. However, depending on the details of a particular simulation, performing URR cross section calculations on-the-fly can significantly increase simulation times. Alternatively, the method of interpolating equiprobable probability bands is demonstrated to produce results that are as accurate as the direct reference solutions, to within arbitrary precision, with high computational efficiency in terms of memory requirements and simulation time. Analyses of a fast spectrum system show that interpolation on a coarse energy-temperature mesh can be used to reproduce reference k-eigenvalue results obtained with cross sections calculated continuously in energy and directly at an exact temperature to within <10 pcm. Probability band data on a mesh encompassing the range of temperatures relevant to reactor analysis usually require around 100 kB of memory per nuclide. Relative to the case in which probability table data generated at a single, desired temperature are used, minor increases in simulation times are observed when probability band interpolation is employed.

Monte Carlo Sensitivity and Uncertainty Analysis with Continuous-Energy Covariance Data

The continuous-energy Monte Carlo (MC) sensitivity and uncertainty (S/U) analysis conducted using the multigroup covariance matrices has a theoretical pitfall in that it is inconsistent with the principle of continuous-energy MC neutronics calculations because the use of the multigroup covariance matrices means treating covariance data as multigroup variables rather than continuous-energy variables. As a way to get around this deficiency and perform the MC S/U analysis on the theoretically consistent principle, this paper presents a new continuous-energy MC S/U formulation which directly utilizes the continuous-energy covariance data in the evaluated nuclear data libraries instead of the multigroup covariance matrices produced by nuclear data processing codes. The validity of the new MC S/U formulation is examined in terms of the input-nuclear-data-induced k uncertainty of the Godiva critical assembly and the TMI-1 pin cell problem by inputting the continuous-energy covariance data of nuclides involved directly into the continuous-energy MC transport calculations by a Seoul National University MC code, McCARD.

Performance and accuracy of criticality calculations performed using WARP – A framework for continuous energy Monte Carlo neutron transport in general 3D geometries on GPUs

In this companion paper to “Algorithmic Choices in WARP – A Framework for Continuous Energy Monte Carlo Neutron Transport in General 3D Geometries on GPUs” (doi:http://dx.doi.org/10.1016/j.anucene.2014.10.039), the WARP Monte Carlo neutron transport framework for graphics processing units (GPUs) is benchmarked against production-level central processing unit (CPU) Monte Carlo neutron transport codes for both performance and accuracy. Neutron flux spectra, multiplication factors, runtimes, speedup factors, and costs of various GPU and CPU platforms running either WARP, Serpent 2.1.24, or MCNP 6.1 are compared. WARP compares well with the results of the production-level codes, and it is shown that on the newest hardware considered, GPU platforms running WARP are between 0.8 and 7.6 times as fast as CPU platforms running production codes. The GPU platforms running WARP are also shown to be between 15% and 50% as expensive to purchase and between 80% and 90% as expensive to operate as equivalent CPU platforms performing at an equal simulation rate.

Neutron Cross Section Processing Methods for Improved Integral Benchmarking of Unresolved Resonance Region Evaluations

In this work we describe the development and application of computational methods for processing neutron cross section data in the unresolved resonance region (URR). These methods are integrated with a continuous-energy Monte Carlo neutron transport code, thereby enabling their use in high-fidelity analyses. Enhanced understanding of the effects of URR evaluation representations on calculated results is then obtained through utilization of the methods in Monte Carlo integral benchmark simulations of fast spectrum critical assemblies. First, we present a so-called on-the-fly (OTF) method for calculating and Doppler broadening URR cross sections. This method proceeds directly from ENDF-6 average unresolved resonance parameters and, thus, eliminates any need for a probability table generation pre-processing step in which tables are constructed at several energies for all desired temperatures. Significant memory reduction may be realized with the OTF method relative to a probability table treatment if many temperatures are needed. Next, we examine the effects of using a multi-level resonance formalism for resonance reconstruction in the URR. A comparison of results obtained by using the same stochastically-generated realization of resonance parameters in both the single-level Breit-Wigner (SLBW) and multi-level Breit-Wigner (MLBW) formalisms allows for the quantification of level-level interference effects on integrated tallies such as k eff and energy group reaction rates. Though, as is well-known, cross section values at any given incident energy may differ significantly between single-level and multi-level formulations, the observed effects on integral results are minimal in this investigation. Finally, we demonstrate the calculation of true expected values, and the statistical spread of those values, through independent Monte Carlo simulations, each using an independent realization of URR cross section structure throughout. It is observed that both probability table and OTF treatments reproduce the true expected values, calculated by averaging the results of many independent simulations, quite well. However, the spread of independent calculation results is shown to be relatively significant. The k eff eigenvalues for fast spectrum systems can differ by more than 250 pcm from one simulation to the next.

Impact of new evaluated nuclear data libraries on core characteristics of innovative reactor designs

The impact of new evaluated nuclear data libraries (ENDF/B, JENDL and JEFF) on core characteristics of thermal and fast spectra innovative reactor designs was investigated. The innovative reactor design with thermal neutron spectrum is represented by small-sized, long-life CANDLE (Constant Axial shape of Neutron flux, nuclide densities and power shape During Life of Energy producing reactors) HTGRs (High Temperature Gas-cooled Reactors) with uranium and thorium fuel cycles, while the one with the fast neutron spectrum is represented by compact, sodium-cooled B&BRs (Breed-and-Burn Reactors). The CANDLE HTGR core characteristics (i.e. kinf, keff, discharge burnup, burning region moving velocity, core life time and axial power peaking factor) were evaluated by a dedicated deterministic neutronics and depletion code based on few-group neutron diffusion theory in 2-D RZ geometry, while the required microscopic cross sections were prepared by using SRAC2006 with JENDL-4.0, ENDF/B-VI and JEFF-3.1 based new SRAC2006 libraries. The B&BR core characteristics were evaluated using a continuous energy Monte Carlo neutron transport code, SERPENT, with JENDL-4.0, ENDF/B-VII.0, ENDF/B-VII.1 and JEFF-3.1.1 nuclear data based new libraries. The impact of new evaluated nuclear data libraries for innovative CANDLE HTGR core characteristics was found significantly larger for thorium fuel cycle than the one for uranium fuel cycle, especially in the discharge burnup, burning region moving velocity which affected the core life time significantly. The findings indicated the needs of more accurate nuclear data libraries for nuclides involved in thorium fuel cycle. The impact of new evaluated nuclear data libraries for innovative B&BR core characteristics was found in a slightly large variation for keff evolution during burnup which in turn affected the estimated core life time. However, a good agreement within the statistical error on the integral kinetic parameters, Doppler reactivity coefficients, coolant density reactivity coefficients, axial and radial expansion reactivity coefficients, power density distributions can be observed among the libraries. In addition, the buildup of fission products and minor actinides were also similar.

Algorithmic choices in WARP – A framework for continuous energy Monte Carlo neutron transport in general 3D geometries on GPUs

In recent supercomputers, general purpose graphics processing units (GPGPUs) are a significant faction of the supercomputer’s total computational power. GPGPUs have different architectures compared to central processing units (CPUs), and for Monte Carlo neutron transport codes used in nuclear engineering to take advantage of these coprocessor cards, transport algorithms must be changed to execute efficiently on them. WARP is a continuous energy Monte Carlo neutron transport code that has been written to do this. The main thrust of WARP is to adapt previous event-based transport algorithms to the new GPU hardware; the algorithmic choices for all parts of which are presented in this paper. It is found that remapping history data references increases the GPU processing rate when histories start to complete. The main reason for this is that completed data are eliminated from the address space, threads are kept busy, and memory bandwidth is not wasted on checking completed data. Remapping also allows the interaction kernels to be launched concurrently, improving efficiency. The OptiX ray tracing framework and CUDPP library are used for geometry representation and parallel dataset-side operations, ensuring high performance and reliability.

Monte Carlo Fuel Temperature Coefficient Estimation by an Adjoint-Weighted Correlated Sampling Method

It is very time-consuming to obtain a high-precision Monte Carlo (MC) estimate of the fuel temperature reactivity coefficient (FTC) through direct subtraction of two reactivity values from MC calculations at two different fuel temperatures. As an alternative to the direct subtraction MC estimate of the FTC, this paper presents a new method based on the adjoint-weighted correlated sampling technique. The new method translates the change in fuel temperature as the corresponding changes in both the microscopic cross sections and the transfer probabilities in scattering kernels described by the free gas model. The effectiveness of the new method is examined through continuous-energy MC neutronics calculations for pressurized water reactor pin cell and CANDU pressurized heavy water reactor lattice problems. The isotope-wise and reaction-type–wise contributions to the FTCs in the two problems are examined for two free gas models: the constant-cross-section and the resonance-cross-section models. It is demonstrated that the new MC method can predict the reactivity change due to fuel temperature variation as accurately as the conventional, more time-consuming direct subtraction MC method.

Benchmarking ENDF/B-VII.1, JENDL-4.0 and JEFF-3.1.1 with MCNP6

Recent releases of three major world nuclear reaction data libraries, ENDF/B-VII.1, JENDL-4.0, and JEFF-3.1.1, have been tested extensively using benchmark calculations. The calculations were performed with the latest release of the continuous energy Monte Carlo neutronics code MCNP, i.e. MCNP6. Three types of benchmarks were used, viz. criticality safety benchmarks, (fusion) shielding benchmarks, and reference systems for which the effective delayed neutron fraction is reported.For criticality safety, more than 2000 benchmarks from the International Handbook of Criticality Safety Benchmark Experiments were used. Benchmarks from all categories were used, ranging from low-enriched uranium, compound fuel, thermal spectrum ones (LEU-COMP-THERM), to mixed uranium-plutonium, metallic fuel, fast spectrum ones (MIX-MET-FAST). For fusion shielding many benchmarks were based on IAEA specifications for the Oktavian experiments (for Al, Co, Cr, Cu, LiF, Mn, Mo, Si, Ti, W, Zr), Fusion Neutronics Source in Japan (for Be, C, N, O, Fe, Pb), and Pulsed Sphere experiments at Lawrence Livermore National Laboratory (for 6Li, 7Li, Be, C, N, O, Mg, Al, Ti, Fe, Pb, D2O, H2O, concrete, polyethylene and teflon). The new functionality in MCNP6 to calculate the effective delayed neutron fraction was tested by comparison with more than thirty measurements in widely varying systems. Among these were measurements in the Tank Critical Assembly (TCA in Japan) and IPEN/MB-01 (Brazil), both with a thermal spectrum, two cores in Masurca (France) and three cores in the Fast Critical Assembly (FCA, Japan), all with fast spectra.The performance of the three libraries, in combination with MCNP6, is shown to be good. The results for the LEU-COMP-THERM category are on average very close to the benchmark value. Also for most other categories the results are satisfactory. Deviations from the benchmark values do occur in certain benchmark series, or in isolated cases within benchmark series. Such instances can often be related to nuclear data for specific non-fissile elements, such as C, Fe, or Gd. Indications are that the intermediate and mixed spectrum cases are less well described.The results for the shielding benchmarks are generally good, with very similar results for the three libraries in the majority of cases. Nevertheless there are, in certain cases, strong deviations between calculated and benchmark values, such as for Co and Mg. Also, the results show discrepancies at certain energies or angles for e.g. C, N, O, Mo, and W.The functionality of MCNP6 to calculate the effective delayed neutron fraction yields very good results for all three libraries.

Two practical methods for unionized energy grid construction in continuous-energy Monte Carlo neutron transport calculation

A considerable speed-up in continuous-energy Monte Carlo neutron transport calculation can be achieved by using the same unionized energy grid for all point-wise reaction cross sections. This speed-up results from the fact that time-consuming grid iteration is reduced to minimum, and if the unionized grid is constructed by combining the grids of all nuclides, there is no loss of data or accuracy in the calculation. The drawback of this approach is that computer memory is wasted for storing a large number of redundant data points. Memory usage may become a problem, especially in burnup calculation, in which the irradiated materials consist of several hundred actinide and fission product nuclides. The grid size easily increases to over 1 million points, requiring tens of gigabytes of memory for storing the cross section data. This paper presents two practical methods for reducing the memory demand, while trying preserve the accuracy of the original data. The calculation routines are included in the PSG2/Serpent Monte Carlo reactor physics burnup calculation code and the methods are tested in a BWR assembly burnup calculation.

Benchmarking ENDF/B-VII.0

The new major release VII.0 of the ENDF/B nuclear data library has been tested extensively using benchmark calculations. These were based upon MCNP-4C3 continuous-energy Monte Carlo neutronics simulations, together with nuclear data processed using the code NJOY. Three types of benchmarks were used, viz., criticality safety benchmarks, (fusion) shielding benchmarks, and reference systems for which the effective delayed neutron fraction is reported. For criticality safety, more than 700 benchmarks from the International Handbook of Criticality Safety Benchmark Experiments were used. Benchmarks from all categories were used, ranging from low-enriched uranium, compound fuel, thermal spectrum ones (LEU-COMP-THERM), to mixed uranium-plutonium, metallic fuel, fast spectrum ones (MIX-MET-FAST). For fusion shielding many benchmarks were based on IAEA specifications for the Oktavian experiments (for Al, Co, Cr, Cu, LiF, Mn, Mo, Si, Ti, W, Zr), Fusion Neutronics Source in Japan (for Be, C, N, O, Fe, Pb), and Pulsed Sphere experiments at Lawrence Livermore National Laboratory (for 6Li, 7Li, Be, C, N, O, Mg, Al, Ti, Fe, Pb, D2O, H2O, concrete, polyethylene and teflon). For testing delayed neutron data more than thirty measurements in widely varying systems were used. Among these were measurements in the Tank Critical Assembly (TCA in Japan) and IPEN/MB-01 (Brazil), both with a thermal spectrum, and two cores in Masurca (France) and three cores in the Fast Critical Assembly (FCA, Japan), all with fast spectra. In criticality safety, many benchmarks were chosen from the category with a thermal spectrum, low-enriched uranium, compound fuel (LEU-COMP-THERM), because this is typical of most current-day reactors, and because these benchmarks were previously underpredicted by as much as 0.5% by most nuclear data libraries (such as ENDF/B-VI.8, JEFF-3.0). The calculated results presented here show that this underprediction is no longer there for ENDF/B-VII.0. The average over 257 benchmarks deviates only 0.017% from the measured benchmark value. Moreover, no clear trends (with e.g. enrichment, lattice pitch, or spectrum) have been observed. Also for fast spectrum benchmarks, both for intermediately or highly enriched uranium and for plutonium, clear improvements are apparent from the calculations. The results for bare assemblies have improved, as well as those with a depleted or natural uranium reflector. On the other hand, the results for plutonium solutions (PU-SOL-THERM) are still high, on average (over 120 benchmarks) roughly 0.6%. Furthermore there still is a bias for a range of benchmarks based on cores in the Zero Power Reactor (ANL) with sizable amounts of tungsten in them. The results for the fusion shielding benchmarks have not changed significantly, compared to ENDF/B-VI.8, for most materials. The delayed neutron testing shows that the values for both thermal and fast spectrum cases are now well predicted, which is an improvement when compared with ENDF/B-VI.8.

Critical and Subcritical Mass Calculations of Curium-243 to -247 Based on JENDL-3.2 for Revision of ANSI/ANS-8.15

Critical and subcritical masses were calculated for a sphere of five curium isotopes from 243Cm to 247Cm in metal and in metal-water mixtures considering three reflector conditions: bare, with a water reflector or a stainless steel reflector. The calculation were made mainly with a combination of a continuous energy Monte Carlo neutron transport calculation code, MCNP, and the Japanese Evaluated Nuclear Data Library, JENDL-3.2. Other evaluated nuclear data files, ENDF/B-VI and JEF-2.2, were also applied to find differences in calculation results of the neutron multiplication factor originated from different nuclear data files. A large dependence on the evaluated nuclear data files was found in the calculation results: more than 10%Δk/k relative differences in the neutron multiplication factor for a homogeneous mixture of 243Cm metal and water when JENDL-3.2 was replaced with ENDF/B-VI and JEF-2.2, respectively; and a 44% reduction in the critical mass by changing from JENDL-3.2 to ENDF/B-VI for 246Cm metal. The present study supplied basic information to the ANSI/ANS-8.15 Working Group for revision of the standard for nuclear criticality control of special actinide elements. The new or revised values of the subcritical mass limits for curium isotopes accepted by the ANSI/ANS-8.15 Working Group were finally summarized.

Critical and Subcritical Mass Calculations of Curium-243 to-247 Based on JENDL-3.2 for Revision of ANSI/ANS-8.15.

Critical and subcritical masses were calculated for a sphere of five curium isotopes from 243Cm to 247Cm in metal and in metal-water mixtures considering three reflector conditions: bare, with a water reflector or a stainless steel reflector. The calculation were made mainly with a combination of a continuous energy Monte Carlo neutron transport calculation code, MCNP, and the Japanese Evaluated Nuclear Data Library, JENDL-3.2. Other evaluated nuclear data files, ENDF/B-VI and JEF-2.2, were also applied to find differences in calculation results of the neutron multiplication factor originated from different nuclear data files. A large dependence on the evaluated nuclear data files was found in the calculation results: more than 10%Δk/k relative differences in the neutron multiplication factor for a homogeneous mixture of 243Cm metal and water when JENDL-3.2 was replaced with ENDF/B-VI and JEF-2.2, respectively; and a 44% reduction in the critical mass by changing from JENDL-3.2 to ENDF/B-VI for 246Cm metal. The present study supplied basic information to the ANSI/ANS-8.15 Working Group for revision of the standard for nuclear criticality control of special actinide elements. The new or revised values of the subcritical mass limits for curium isotopes accepted by the ANSI/ANS-8.15 Working Group were finally summarized.