Search NASA⌕ Search

SEARCH · Search NASA

Results for “MCNP Model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

High-level fuel fabrication facility designs from discrete-event simulation

Like other industrial processes, the production of metallic nuclear fuels (MNF) requires that fabrication facilities be able to reliably meet production demands, operate efficiently, and adhere to federal and local safety regulations. In turn, the set of design variables employed by such a facility, such as operations policies, infrastructure and machinery purchased, and the type and number of staff hired, directly impact a facility’s ability to satisfy these goals. Therefore, facility designers must carefully determine which set of design variable values optimally satisfies these constraints. In this paper, we explore how values for these high-level design variables, namely hiring requirements, can be determined in the context of nuclear fuel manufacturing through the coupling of physics-based and discrete-event simulation technologies. Using the Versatile Test Reactor (VTR) program as a case study, we demonstrate how SCALE and MCNP nuclear physics model outputs can be integrated into ExtendSim discrete-event simulation (DES) models of the fuel fabrication process to determine the optimal number of staff hired to ensure fuel production goals are met, operations comply with effective dose limit regulations, and overall project costs are reduced.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Challenges and Opportunities Specific to Microreactors: Compact and Light-Weight Shielding During Operation and Transportation

The problem of designing a shield for a mobile microreactor has been addressed. The principal challenges encountered in this endeavor are the dual goals of increasing radiation protection and decreasing weight. It is recognized that these two goals are contradictory. Increasing shielding naturally increases weight. Conversely, decreasing weight naturally calls for decreasing material inventory, including possibly decreasing shield mass. The operational and regulatory drivers for shielding and transportation constraints are identified and summarized. The main constraints pertain to volume, mass, radiation dose, material feasibility and availability, and decommissioning expectations. The latter are limits on residual contamination and dose rates after the reactor is removed. Progress toward solving the challenge of shielding a transportable microreactor is made by considering a tungsten tetraboride (WB4) shell surrounding the reactor vessel. The shield and a generic microreactor were modeled using the MCNP code. It was determined that the presence of a 16-cm-thick WB4 shield reduces the dose equivalent rate during operations to 3.6 mrem/h just outside of a concrete vault and dome and soil berm surrounding the reactor, which is very near the occupational dose rate limits. Without the WB4 layer, the dose rate at the same location outside the berm would be 120 mrem/h. Most importantly, the presence of the WB4 shield substantially protects the surrounding concrete and soil from activation, thus decreasing the duration of administrative control for the site following removal of the reactor from 16 to 3.4 years. At that point, the site may be decommissioned and returned to unrestricted public access. A limitation on the incorporation of a WB4 shield as an intrinsic component of a microreactor is that it results in a design that is transportable only by using the largest road conveyances and that the load would require special permits for transit through the U.S. highway system. New directions for further improvements and possibly more effective solutions to the shielding problem of a mobile microreactor are discussed in the conclusion section.

Microreactor↗

Calibration of ITER Instant Power Neutron Monitors: Recommended Scenario of Experiments at the Reactor

Instant power is a key parameter of the ITER. Its monitoring with an accuracy of a few percent is an urgent and challenging aspect of neutron diagnostics. In a series of works published in Problems of Atomic Science and Technology, Series: Thermonuclear Fusion under a common title, the step-by-step neutronics analysis was given to substantiate a calibration technique for the DT and DD modes of the ITER. A Gauss quadrature scheme, optimal for processing “expensive” experiments, is used for numerical integration of {sup 235}U and {sup 238}U detector responses to the point sources of 14-MeV neutrons. This approach allows controlling the integration accuracy in relation to the number of coordinate mesh points and thus minimizing the number of irradiations at the given uncertainty of the full monitor response. In the previous works, responses of the divertor and blanket monitors to the isotropic point sources of DT and DD neutrons in the plasma profile and to the models of real sources were calculated within the ITER model using the MCNP code. The neutronics analyses have allowed formulating the basic principles of calibration that are optimal for having the maximum accuracy at the minimum duration of in situ experiments at the reactor. In this work, scenarios of the preliminary and basic experimental ITER runs are suggested on the basis of those principles. It is proposed to calibrate the monitors only with DT neutrons and use correction factors to the DT mode calibration for the DD mode. It is reasonable to perform full calibration only with {sup 235}U chambers and calibrate {sup 238}U chambers by responses of the {sup 235}U chambers during reactor operation (cross-calibration). The divertor monitor can be calibrated using both direct measurement of responses at the Gauss positions of a point source and simplified techniques based on the concepts of equivalent ring sources and inverse response distributions, which will considerably reduce the amount of measurements. It is shown that the monitor based on the average responses of the horizontal and vertical neutron chambers remains spatially stable as the source moves and can be used in addition to the staff monitor at neutron fluxes in the detectors four orders of magnitude lower than on the first wall, where staff detectors are located. Owing to low background, detectors of neutron chambers do not need calibration in the reactor because it is actually determination of the absolute detector efficiency for 14-MeV neutrons, which is a routine out-of-reactor procedure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Upgrade the EMPIRE code as needed for evaluation and production of covariances for actinides

Treatment of actinides in the EMPIRE code has been updated by adding capabilities necessary for effective assimilation of nuclei that undergo fission. The upgrades comprise calculations of Prompt Fission Neutron Spectra (PFNS), and the average cosine of elastically scattered neutrons ($\overline{µ}$) and include sensitivity and covariance computations for both quantities. In addition, EMPIRE outputs the average number of neutrons released per fission event ($\overline{ν}$) that is read from the existing evaluated file and scaled by a factor. Kalman filter is used to fit calculated PFNS, $\overline{µ}$, and $\overline{ν}$ to the experimental data and to produce covariances for the quantities and respective model parameters. These updates expand EMPIRE assimilation possibilities by adding PFNS, $\overline{ν}$, and $\overline{µ}$ to the physical quantities that might be perturbed along with the cross sections during the adjustment. We describe PFNS models used in the EMPIRE code, the modifications required to accommodate these models, and detail updates to the EMPIRE scripts used for fitting, and illustrate the new features applying them to specific nuclei. Finally, we report on actual assimilation of 239 Pu performed by calculating sensitivity of the JEZEBEL $k_{eff}$ to the full set of reaction model parameters (EMPIRE+MCNP) and using Kalman filter to perform adjustment of the parameters. This work was performed according to the goals set forth in milestones M3FT-12BN0210011 and M3FT-12BN0210012.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Integrated Approach to Fluoride High Temperature Reactor Technology and Licensing Challenges (FHR-IRP)

This document constitutes progress made on this project by the Georgia Tech led Integrated Approach to Fluoride High Temperature Reactor Technology and Licensing Challenges (FHR-IRP). This is the first draft of the project final report. Chapter 3 is focused on Tritium management. Tritium management in Fluoride-salt High-temperature Reactors (FHRs) presents a major operational concern and a significant challenge. Tritium is generated via multiple pathways including various interactions in fuel and structural materials, as well as coolant (primarily interaction with lithium and beryllium), leading to a complex evolution and spatial distribution of tritium generation rate. FHR generates orders of magnitude more tritium than in LWRs, and has potential for a significant release, via permeation through structural materials (due to high temperature), and through heat exchanger(s). Therefore, accurate estimate is necessary. For that purpose relevant tritium production pathways were identified. A methodology to obtain space-time dependent tritium generation rate was developed and implemented, employing SCALE and MCNP codes. Detailed models of the reference Advanced High Temperature Reactor (AHTR) design were developed. A multi-parametric metrics was introduced to capture relevant tritium production indicators, including the instantaneous tritium production rate, total tritium production (neglecting subsequent decay), and accumulated tritium production (accounting for the decay). The obtained results will inform tritium management, including the sizing of the components, and to gain insight on the impact of the core design parameters on tritium generation rates. Research effort on Liquid Salt Impurity Removal, Redox and Corrosion Control is discussed in Chapter 4. Section 4.1 is an introduction chapter. Based on the current available knowledge and operation experience on previous MSRs and other types of nuclear reactors, the report identifies the major impurities that are related to either material compatibility, salt properties or radioactive source terms, and the potential contamination sources of these impurities, their transport in the primary coolant loop and available cleanup methodologies for a particular group of impurity. In section 4.2, the available purification methods were discussed, including the liquid/liquid extraction using liquid Bi-Li which has been applied in MSBR and electrochemical purification which is developed based on Pyroprocessing. In this section, the new multi-stage and multi-salt-component was developed and discussed. Section 4.3 discussed the redox potential control for salt chemistry and corrosion mitigation, including the control theory and the soluble/soluble buffer control method. The section also developed method on how to analyze soluble/soluble redox reaction in a molten salt system and designed an electrochemical sensor for redox potential control. Corrosion models were also reported in this section. Section 4.4 focused on how to measure the fundamental properties in a molten salt. As an example, the section focused on La fluoride in the molten salt. The section developed methods on how analyze electrochemical measurement data to obtain the fundamental properties including the redox potential, diffusion coefficient and exchange current density. All these fundamental information is important for salt purification system development. Section 4.5 focused on measurements of corrosion products in the salts. Because of the different behaviors between Cr and La in the salt, new analysis method for Cr fluorides in the salt was developed. Some results for Fe fluoride in the salt and corrosion measurements were also reported. Section 4.6 demonstrated the separation of oxide from the salt. This section also measured the solubility of La 2 O 3 in the salt as function of temperature. Section 4.7 focused on the separation of La and Ce from the molten salt. These two elements were selected because La was found to be separated by using inert electrode, while for separation of Ce, reactive electrode must be selected. Section 4.8 demonstrated the separation iodine ions from the molten salt. Two kinds of working electrode materials (metal W and graphite) were investigated. Section 4.9 focused on the fundamental data of fission products and actinides in liquid Bi. The chapter reviewed the available data, then based a CALPHAD model, the available data were broadened to the concentration-dependent cases.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Neutron Transmission Imaging with a Portable D-T Neutron Generator

Fast-neutron transmission imaging provides complementary information to x-ray transmission imaging. While fast neutron imaging resolution is generally below x-ray imaging, 14-MeV neutrons have an advantage over portable x-ray systems. Neutrons have higher transmission through high-Z materials due to a more uniform attenuation as a function of material atomic number Z compared to X-rays, and can therefore image low-Z materials inside high-Z materials. As a result, portable neutron transmission imaging has many applications, including inspection of concrete and welds for corrosion in vehicles, bridges, and other infrastructure, measurement of material levels in containers, and inspection of suspicious packages. Fast-neutron imaging is also more practical for field use than thermal-neutron imaging due to the size and shielding requirements typical of thermal-imaging systems compared to the availability of small 14.1 MeV D-T neutron generators. However, there are limitations in portable fast-neutron imaging systems, including limited neutron output, limited light produced by neutron scintillators, and lower resolution due to neutron source spot size and 2-3 mm scintillator thickness. In addition, digital-panel dark-noise is roughly 100x higher than neutron scintillator light, and variations in noise across the panel and in time is comparable to the imaging signal. Here we discuss recent efforts in developing a portable fast-neutron radiography system, including an improved neutron scintillator, mitigation of panel noise, and new commercial portable D-T neutron generators. We also present MCNP efforts to model neutron imaging, including scintillator resolution and the effects of neutron scattering from the object and surrounding materials.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Secondary Pion Production using the LINAC Beam on Graphite and Tungsten Targets

The Fermilab 400-MeV Linac primary proton beam can produce precision, single-species secondary beams using a production target in the MeV Test Area experimental hall; specifically, pions from 4 up to 120 MeV kinetic energy. (These low energy pions rapidly decay into a muon beam in a secondary collection and transport beamline.) Although graphite is a standard production target material given its ultra-high temperature tolerance and low impact on the post-target primary beam, secondary pion/muon total production cross sections increase approximately as Z1/3 for positive muons and N2/3 for negative muons. Higher mass targets such as tungsten can potentially increase low-energy m+ and m- rates by factors of 3 and 8, respectively. Initial production studies have confirmed higher muon yield from heavy targets (tungsten vs carbon), but also, surprisingly, show significant differences in pion production between modern hadronic models (GENIEhad) and between GEANT and MCNP, state of the art production modeling codes. This project will further investigate the production discrepancies between these two target materials and between different physics models with high-statistics runs using the GEANT-based code, G4Beamline to better quantify the results and compare with production data which will be obtained using the MTA secondary beamline.

Ahmed, Shiza↗

Demonstrating Safeguards Applications of DRiFT Gas Detector Capabilities [Slides]

DRiFT, a Detector Response Function Toolkit, post-processes MCNP output to model detector response. Three primary components: scintillator (focus of FY 21 TED work), gas detectors (focus of this FY), and semiconductor detectors. Capability to model scintillator detectors has previously been demonstrated.

47 OTHER INSTRUMENTATION↗

Design and testing of a combined neutron and gamma assay system for measuring fissile material in fuel and waste items - 15340

A combined neutron and gamma ray assay system for the measurement of fissile material in both fuel and waste items has been designed, modelled, developed and tested. The purpose of the system is to measure and characterize a wide range of fissile materials including uranium, plutonium, americium and mixed oxide fuel arising from a variety of processes and facilities across the United Kingdom. Dounreay Site Restoration Limited at Dounreay currently stores these items form the United Kingdom Nuclear Decommissioning Authority's un-irradiated fuel inventory. This inventory comprises historic materials from ex-United Kingdom Atomic Energy Authority nuclear fuel projects and experiments together with fast reactor fuel materials and some commercial fuel items as delivered to Dounreay for recovery. The combined assay system has been designed and constructed to support the operation of a fuel processing facility, which will be operated at Dounreay to treat these historic legacy fuel materials, making them suitable for long term storage. Appropriate characterization of the fissile material is required for safety, criticality control and nuclear materials safeguards purposes prior to final disposition. The system consists of both gamma ray and neutron measurement components. The results of a neutron measurement are combined with the gamma ray isotopic ratio measurement results to generate the assay result for the fissile material being characterized. The gamma ray component of the system employs a safeguards quality high purity Germanium detector with electro-mechanical cooling and the measurements of fissile material containers are made in a shielded chamber with a tin and copper graded lining. Containers of fissile material are rotated during the gamma ray measurement and a variable automatic steel collimator aperture mechanism is employed to adjust the measured count rate to control detector dead time. Gamma ray spectra are analysed using the plutonium and uranium isotopic ratio code PC/FRAM. The gamma ray system can also be configured to perform a direct assay of small quantities of uranium bearing material. For this assay measurement a series of cadmium filters can be positioned in front of the detector to reduce the contribution from 59 keV gamma rays from Am-241. In addition to gamma ray isotopic measurements, both active and passive neutron assay is performed by the system. The neutron-measuring component consists of a thick walled cylindrical polyethylene chamber, which provides neutron moderation for 30 one-inch diameter He-3 detector tubes located in re-entrant channels in the polyethylene moderator. Two removable polyethylene-inserts incorporating lower plug assemblies are provided, one with and one without a cadmium lining. Also two top plug units are provided, again one with cadmium and one without a cadmium liner. The Neutron Subsystem can be operated both as a passive and as an active neutron coincidence counter. In both fast and thermal active neutron mode, the neutron component employs two americium-beryllium (Am-Be) alpha-n neutron sources, one positioned in the top and the other in the bottom plug unit. The neutron sources are removed when the counter is operated in passive mode. Passive neutron detection efficiency and active neutron performance has been modelled using the MCNP Monte Carlo computer simulation code as part of the design process. In the absence of fissile samples, operation of the gamma ray component of the system has been tested initially using calibrated Eu-152 and Ba-133 gamma ray sources and the use of archived spectra from fissile material measurements. Initial testing of the neutron system has employed a calibrated Cf-252 spontaneous fission neutron source. Neutron testing has been confirmed using further MCNP simulations. The results of test measurements with the system are presented and compared with the results of MCNP measurement simulations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Perform Design Support with MCNP for New Measurements

This report incorporates our work carried out during our 5-month internship at LANL under an internship agreement with EAMEA (École des Applications Militaires de l’Énergie Atomique). After outlining the context in which we worked, we present our work as aid to modeling and predicting the neutronic behavior of nuclear systems, with a view to carrying out criticality experiments qualifying the MCNP code as part of innovative projects. Fourth generation reactors will enable to tackle a lot of issues such as environmental crisis, affordable energy access, and nuclear waste management. They seem to be one of the keys for a sustainable future. Most of the projects that emerge nowadays include the use of HALEU (high assay low enriched uranium) or MOX recycled fuels. Our projects are part of this dynamic and addresses concrete scientific research issues in the nuclear field. Studies of HALEU package are essential to anticipate the need, therefore the Optimus L (OPTImal Modular Universal Shipping cask technology) designed by NAC (Nuclear Assurance Corporation) international but filled with 20 % enriched uranium dioxide (UO 2 ) will be studied to support safe transportation. However, it seems there is no benchmark with a high correlation with the combination of this fuel and this package to validate MCNP simulations. As such, the study will focus on the development of new criticality safety benchmarks for this case. On the other hand, there is a great need for critical benchmarks in the intermediate energy range with MOX fuel. An IER (Integral Experiment Request) has then be requested to answer it through a partnership between French institution IRSN and U.S. Department of Energy's Nuclear Criticality Safety Program (NCSP). This experience planned for 2025 requires to gather calculated data through a MCNP model to be realized safely. Finally, a presentation of our one-week experience at the DAF as part of our discovery of criticality experiments will be introduced in Appendix 1: Week at the DAF (Device Assembly Facility)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Experimental observation and integrated modelling of proton-beryllium fusion in He and D plasmas at JET

Validated integrated modelling of JET ITER-like wall experiments in which fusion performance is driven by reactions between fast ions and intrinsically present metal wall impurities is presented. A steady-state L-mode plasma with dominant proton-beryllium fusion and neutron yields of up to ≈ 6·10 13 s -1 is developed in He and D, via radiofrequency heating of a H minority. The fusion drive is unambiguously confirmed by the neutral particle analyser, fast ion loss detector, and γ-ray diagnostics. Experiments are analysed via an integrated modelling framework, developed to model the two-stage proton beryllium-fusion chain and produce high-fidelity fusion product source terms. The modelling chain comprises TRANSP and JETTO for plasma core modelling, LOCUST for full orbit product tracking and collisional slowing-down, DRESS to resolve two- and three-body fusion kinematics, and MCNP for neutron transport calculations. Modelling shows that the primary 9 Be(p,n) 9 B reaction is the dominant neutron emitter at naturally present concentrations of beryllium in these experiments. The yield contribution of secondary reactions between fusion products and beryllium, 9 Be(d,n) 10 B and 9 Be(α,n) 12 C, is found to be negligible. The proton-deuteron knock-on effect in D plasmas is modelled, which is calculated to contribute ≈ 25% to the total neutron yield. For both He and D discharges the total computed neutron rates match fission chamber (FC) measurements within the combined experimental and computational uncertainty, with an average discrepancy of ≈ ± 20%. Realistic proton-beryllium neutron sources are propagated through JET’s MCNP neutron transport model which shows that 235 U FCs’ response is sensitive to p–Be source changes, with up to ≈ 10% variation compared to a D–D neutron source. We show that the high-energy tail of the fast proton minority can be studied with multi-foil neutron activation. The framework is also applied to the study of interactions between fast protons and boron impurities, of relevance to ITER. We calculate that in JET conditions a significant alpha source with DT-like energies could be generated through 11 B(p, α)2α fusion, and detected via γ-emission in secondary interactions between fast alphas and boron. The work represents an important step towards validating predictive integrated modelling capabilities for non-standard fusion reactions.

JET↗

Whole-core comet solutions to a 3-dimensional PWR benchmark problem with gadolinium

A pressurized water reactor (PWR) benchmark problem with gadolinium was used to determine the accuracy and computational efficiency of the coarse mesh radiation transport method COMET. The benchmark problem contains 193 square fuel assemblies. The COMET solution (eigenvalue, assembly averaged and fuel pin averaged fission density distributions) was compared with those obtained from the corresponding Monte Carlo reference solution using the same 2-group material cross section library. The comparison showed that both the core eigenvalue and fission density distribution averaged over each assembly and fuel pin predicated by COMET agree very well with the corresponding MCNP reference solution if the incident flux response expansion used in COMET is truncated at 2nd order in the two spatial and the two angular variables. The benchmark calculations indicate that COMET has Monte Carlo accuracy. In, particular, the eigenvalue difference between the codes ranged from 17 pcm to 35 pcm, being within 2 standard deviations of the calculational uncertainty. The mean flux weighted relative differences in the assembly and fuel pin fission densities were 0.47% and 0.65%, respectively. It was also found that COMET's full (whole) core computational speed is 30,000 times faster than MCNP in which only 1/8 of the core is modeled. It is estimated that COMET would have been about over 6 orders of magnitude faster than MCNP if the full core were also modeled in MCNP. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High-purity germanium semiconductor modeling in the detector response function toolkit

In this study, we have extended the detector response function toolkit (DRiFT) to provide modeling capabilities of semiconductor sensors. DRiFT provides realistic nuclear instrumentation response by post-processing Monte-Carlo N-particle (MCNP®) radiation transport outputs. MCNP® is capable of modeling radiation transport in complex environments, but has limited detector physics and readout electronics modeling capabilities. Semiconductor detector response can be calculated with a high-fidelity for a flexible range of environments by utilizing MCNP® to simulate radiation interactions inside of detector volumes, and then using DRiFT to model charge transport and signal formation in the semiconductor, as well as the readout electronics. DRiFT models charge transport in the semiconductor, the preamplifier, shaping amplifier, pulse pile-up, and electronic noise to generate detector response. The semiconductor application in DRiFT can model a range of semiconductor materials, shapes, and sizes; and is demonstrated here for a large volume coaxial high-purity germanium (HPGe) detector. Here, we compare detector response functions of a coaxial HPGe detector with measurement of 60 Co, 133 Ba, and 137 Cs at varying count rates, and we conduct a parameter study to demonstrate the effect of changing parameters in the DRiFT simulation. The HPGe detector response function shows excellent agreement with measurements of difference sources with varying dead times and count rates.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

MCNPy

SAND2026-20425O MCNPy runs and analyzes simulations from MCNP, a software that models radiation transport of neutrons and gamma rays. MCNPy uses Python to start MCNP, retrieve event data files, and convert them into graph structures for detailed analysis. It offers visualization tools, including 2D views of particle histories, making complex simulation data easier to interpret for researchers and engineers. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

Nowack, Aaron [Sandia National Lab. (SNL-CA), Live↗

Review of the Oak Ridge National Laboratory (ORNL) Neutronic Calculations Regarding the Conversion of the High Flux Isotope Reactor (HFIR) to the Use of Low Enriched Uranium (LEU) Fuel

In awareness of the risk caused by the proliferation of nuclear materials, the international community initiated in 1978 the Reduced Enrichment for Research and Test Reactors (RERTR) Program. The goal of the RERTR program is to seek solutions to convert the reactors using High Enriched Uranium (HEU) fuel ( 235 U/U ≥ 20 wt. %) to the use of Low Enriched Uranium (LEU) Fuel ( 235 U/U < 20 wt. %). Among the 200 reactors worldwide currently in the scope of the program, about 70 have been converted or shut down prior to conversion. In the U.S., six high reactors remain to be converted. One of the U.S. reactors is the High Flux Isotope Reactor (HFIR) located at Oak Ridge National Laboratory (ORNL). This multipurpose reactor is mainly used for neutron scattering experiments and isotope production. It has two fuel elements made of involute-shaped fuel plates. The ORNL Research Reactor Division (RRD) has prepared during FY2011 two series of calculations, one neutronic and one thermal-hydraulic (TH), in support of the conversion activities for HFIR. The neutronic calculations cover mainly the evaluation of the matrix of performance (cycle length and flux), the distribution of power in the fuel elements, kinetics parameters and reactivity coefficients. As defined by the RRD procedures, the calculations have to be reviewed. Per internal decision, the RRD proposed an external reviewer, Argonne National Laboratory (ANL). The results of the Argonne neutronic review are presented in this report. The results of the TH review are presented in a separate report. The main conclusion of the ORNL reactor physics analyses is that the proposed reference HFIR LEU core (using the so-called monolithic UMo fuel) could maintain the current level of performance (obtained with a HEU core at 85MW) if it is operated at 100MW. However, to operate within the same safety limits as the current HEU core, the thermal hydraulic analyses show that the peaking of the power distribution occurring at the bottom of the fuel plates has to be reduced. The RRD has found a potential solution by reducing the fuel thickness at the bottom of the fuel plate. This would require addition of an axial grading to the already radially-graded fuel. Thus, the review has been focused on independently recalculating all of the ORNL calculations and making an as fair as possible comparison between the current 85MW HEU core and the reference 100MW LEU core design. ORNL and ANL have used the same computational tools to perform the work. The codes MCNP and ORIGEN have been used to perform the steady state and depletion calculations, respectively. The communication between the two codes has been carried out by the code VESTA. Even though the code MCNP is able to handle complex geometries, it cannot model the involute shape. Nonetheless, ANL has developed an innovative methodology to model the involute shape with MCNP with an approximation that can be made as accurate as desired. Thus, ANL has been able to perform the calculations using an explicit (or discrete) description of the fuel elements. The current modeling used by ORNL consists in homogenizing the different fuel element components (cladding, fuel, coolant…). Although this method is easier to implement, it may also affect the reactor physics which may have an impact on the calculated results. Even though some differences in the physics parameters (flux for instance) have been detected between the HEU explicit (ANL) and homogenized (ORNL) models, they are small enough to not significantly impact the performance parameters. However, this is not the case for the LEU core where the differences in the physics parameters lead to significant differences in the performance parameters, especially the cycle length (found to be considerably longer in the explicit model). Nonetheless, the others parameters calculated by ANL (power distribution, kinetics parameters…) are in good agreement with those calculated by ORNL. Some inconsistencies have been detected in the original ORNL inputs that have been corrected in the reference ones used by ANL for the review and, therefore, the cycle length difference between the LEU homogenized and explicit models has been amplified. Consequently, the ANL comparison of the HFIR HEU and LEU cores has given different results than the ORNL results. The ORNL matrix of performance for the 100MW LEU reference design shows that the reactor cycle length would be as long and the flux as high as the current 85MW HEU core. In contrast, the ANL matrix of performance shows that the 100MW LEU design would have a significantly longer cycle length (~33days) than the current 85MW HEU core (~25days). In addition, it has been shown that, if the reactor is operated for 33days, the neutron fluence would be high enough to increase the overall reactor performance by a significant factor. If the reactor is forced to shut down after 25days, calculations show that the neutron flux would be high enough to maintain the performance of the neutron scattering experiments but losses are predicted to occur in several irradiation facilities. Increasing the cycle length may not be possible for other reasons. For instance, a longer cycle length can increase the cumulative oxide layer thickness growth on the plate. The increase in the neutron fluence may affect the lifetime of the reactor vessel and, by reducing the time for maintenance, pose operational challenges. An assessment of these as well as other potential questions associated with a 33 day cycle length is needed. Since such an assessment has not been carried out, ANL cannot conclude that the ORNL reference LEU core should or could operate at 100MW. The conclusions of the ANL review are different than those of the ORNL investigation mainly because the explicit and homogenized methodologies used to model the core are different. In this report, an attempt is made to justify and explain as clearly as possible why the ANL methodology is thought to be better suited for HFIR analyses. It is acknowledged that because the ANL and ORNL LEU models have not been benchmarked against experiments, it is possible that doubts could be raised concerning the ability of the models to correctly capture all of the physics behavior of the reactor.

Bergeron, A.↗

Modeling and Analysis of the BUK/BES-5 Fast Reactor Using MCNP [Slides]

The Soviet Union launched a series of satellites designated RORSAT that were designed to essentially collect information on NATO and US Navy. 33 satellites were launched between 1967 and 1988, the majority (31) of which were powered by the BUK reactor design or minor variants of BUK, The motivation for this effort is to better understand this small space reactor design. In addition, there is much to be learned from the issues encountered while modeling this reactor concept and its agreement with the real reactor deployed in the space.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Optimization of a dry, mixed nuclear fuel storage array for nuclear criticality safety - 14339

A dry storage array of used nuclear fuel at the Idaho National Laboratory contains a mixture of more than twenty different types of research and test reactor fuels. Significant neutronic interaction occurs between fuel canisters within this large storage arrangement of 636 fuel storage canisters. Nuclear criticality safety is maintained primarily by moderation control and through control of neutron interaction and leakage provided by the credited geometry. Over the nearly forty years of operation the collected fuels are now arranged in a configuration that does not minimize the system neutron multiplication factor (k{sub eff}). This creates a concern since the entire remaining facility storage capacity cannot be utilized without exceeding the system minimum subcritical state (k{sub safe}) for ensuring nuclear criticality safety. This work determines a more optimal arrangement for the stored fuels with a goal to minimize the system k{sub eff}, but also with a minimum of potential fuel canister relocation movements. The ultimate goal is to establish a solution for this multiple-objective optimization problem by finding the Pareto Front for maximum k{sub eff} reduction versus the number of canister relocations. The process first required additional analysis of the array by Monte-Carlo neutron transport codes (specifically, MCNP) to gain a greater understanding of the interaction characteristics of the disparate fuel types. Some stored fuels are composed of primarily metallic and oxide components that provide very little inherent moderation, while other fuel types contain graphite or zirconium-hydride which can contribute a significant amount of moderation. These factors cause a complex relationship within a fuel arrangement with respect to thermal utilization and axial leakage probability. The new MCNP analysis was utilized to develop a simplified spreadsheet-based empirical model which simulates the array k{sub eff} result. This empirical model was benchmarked against additional MCNP calculations for various arrangements and then employed in the optimization analysis. A Tabu metaheuristic and stochastic approximation optimization was applied to the empirical model with the optimization function to minimize the array k{sub eff}. The optimization algorithm interrogated the empirical model to determine the optimal rearrangement for up to sixty canister movements. The Pareto Front was found, which shows that a limited number of canister relocation movements could allow the existing thirty-one empty array positions to be filled with common research reactor fuels while also yielding a reduction in the overall array k{sub eff} relative to the current arrangement. This will impart an improvement in the facility utilization while also offering an enhancement in the safety margin. The reduction in array reactivity has been confirmed by a final set of MCNP calculations with the proposed new arrangements. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Dynamic reactor modeling with applications to SPR and ZEDNA

A dynamic reactor model has been developed for pulse-type reactor applications. The model predicts reactor power, axial and radial fuel expansion, prompt and delayed neutron population, and prompt and delayed gamma population. All model predictions are made as a function of time. The model includes the reactivity effect of fuel expansion on a dynamic timescale as a feedback mechanism for reactor power. All inputs to the model are calculated from first principles, either directly by solving systems of equations, or indirectly from Monte Carlo N-Particle Transport Code (MCNP) derived results. The model does not include any empirical parameters that can be adjusted to match experimental data. Comparisons of model predictions to actual Sandia Pulse Reactor SPR-III pulses show very good agreement for a full range of pulse magnitudes. The model is also applied to Z-pinch externally driven neutron assembly (ZEDNA) type reactor designs to model both normal and off-normal ZEDNA operations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗