SEARCH · Engineering Papers
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.
This software allows for a lidar-generated point cloud to be converted to a water-tight geometry, and finally an MCNP radiation transport input deck automatically. This includes the automatic registration of detectors and placement of pre-generated radiation transport models of these detectors in the scene.
The purpose of the COMMAND code is to provide a flexible, scalable tool for use in developing, integrating, and testing the technologies necessary for achieving autonomous operations of advanced nuclear reactors. The code enables users to efficiently implement custom simulations and experiments by combining key methods from different software modules. These modules are focused on: modeling and simulation tools, such as nuclear simulation tools used for high-fidelity modeling (e.g., Reactor Excursion and Leak Analysis Program [RELAP5-3D] and Monte Carlo N-Particle [MCNP]); machine learning and optimization tools (e.g., anomaly detection and data-driven modeling techniques); advanced control in its digital, high-performance, and supervisory control forms (e.g., proportional integral derivative (PID) control and model predictive control (MPC); and integration with hardware through industrial communication protocols. To ensure flexibility and scalability, COMMAND was designed to be both modular—the software “pieces” all inherit from generic building blocks and can be combined and connected to create complicated simulations—and high performing—designed for parallel processing, enabling simulations and experiments to take advantage of multi-core computers, servers, and nodes. The code is written in the Python programming language due to the language's popularity, active community, and open-source and cross-platform nature. Maintaining consistency with other simulation tools used within the nuclear energy community, users implement simulations and experiments through text input files, which define components, parameters, connections, etc., through lines of text. Given that COMMAND is written in Python, these input files are native Python scripts, and so use the standard Python structure and formatting. This also enables users to take advantage of Python's extensive package library to develop custom capabilities for their specific use cases.
It is a geometry optimizer, written for MCNP geometries but also applicable to other simulation programs. It uses NSGA-II optimization algorithm to change the specified input parameters to iteratively improve the requested tally values over generations.
A tool for converting SCALE input to an equivalent MCNP input.
The Nuclear Material and Composition (NuMaCo) server is a web-based platform for storing, managing, and serving material composition data. It includes the SCALE Standard Composition Library and the PNNL Compendium, along with their associated references and uncertainty data. NuMaCo also provides tools to generate, compare, and validate material input definitions for both SCALE and MCNP.
Modeling and simulation can be used to estimate the detection of different types and strengths of radioactive sources. Simulations of shielded gamma-ray sources in complex environments, such as searches in urban areas using vehicle-mounted detectors, can demonstrate what source strengths and shielding combinations should be able to be detected. For faster calculations, gamma-emitting sources that could have different shielding layers can be simulated ahead of time, resulting in equivalent point source emission distributions that can be used in larger simulation models. Larger-scale models will then run faster using equivalent point sources because the shielding layers around the sources have already been accounted for. Accompanying this paper are 260 equivalent point source emission files for 72 sources with different shielding layers, formatted for three different Monte Carlo codes. Users wanting to request versions of the SCALE and/or MCNP software can find details here: https://rsicc.ornl.gov/
Plasma Material Interaction (PMI) is a major concern in fusion reactor design and analysis. The Material-Plasma Exposure eXperiment (MPEX) will explore PMI under fusion reactor plasma conditions. Samples with accumulated displacements per atom (DPA) damage produced by fast neutron irradiations in the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (ORNL) will be studied in the MPEX facility. This paper presents assessments of the calculated induced radioactivity and resulting radiation dose rates of a variety of potential fusion reactor plasma-facing materials (such as tungsten). The scientific code packages MCNP and SCALE were used to simulate irradiation of the samples in HFIR including the generation and depletion of nuclides in the material and the subsequent composition, activity levels, gamma radiation fields, and resultant dose rates as a function of cooling time. A challenge of the MPEX project is to minimize the radioactive inventory in the preparation of the samples and the sample dose rates for inclusion in the MPEX facility.
During the ignition experimental campaign, the National Ignition Facility (NIF) is expected to perform shots with varying fusion yield (up to 20 MJ or 7.1 x 10 18 neutrons per shot) and a maximum annual yield of 1200 MJ. A detailed MCNP model of the Target Bay (TB) and the two switchyards (SY) has been developed to estimate the post-shot radiation environment inside the facility. During D-T shots, a pulse of 14.1 MeV neutrons streaming outside the Target Chamber (TC) will activate the air present inside the TB and the argon gas inside the laser tubes. Smaller levels of activity are also generated in the SY air and in the argon portion of the SY laser beam path. The activated TB air will be mixed with fresh air from the Operations Support Building (OSB) before release through the stack. Flow of activated air from the Target Bay is controlled by the heating, ventilating, and air conditioning (HVAC) system. 16 N (T 1/2 = 7.13 s) dominates the radiation levels during the first minute following the shot. It is expected that 16 N will decay away during the confinement time before releasing the TB air through the stack. The other major contributors are 13 N (T 1/2 = 9.97 min) and 41 Ar (T 1/2 = 1.83 h). In general a low dose rate of < 1 μSv/h is expected near the stack during the first few hours following a 20 MJ shot. Here, the amount of activated Target Bay air released through the stack is very small and does not pose significant hazard to personnel or the environment. In the mean time, due to a very small leakage rate out of the laser tubes, the activated argon gas decays within the tubes and any resulting release to the environment is insignificant.
In this work, we present a method to map the spectral and spatial distributions of radioactive sources using a limited number of detectors. Locating and identifying radioactive materials is important for border monitoring, in accounting for special nuclear material in processing facilities, and in cleanup operations following a radioactive material spill. Most methods to analyze these types of problems make restrictive assumptions about the distribution of the source. In contrast, the source mapping method presented here allows an arbitrary three-dimensional distribution in space and a gamma peak distribution in energy. To apply the method, the problem is cast as an inverse problem where the system’s geometry and material composition are known and fixed, while the radiation source distribution is sought. A probabilistic Bayesian approach is used to solve the resulting inverse problem since the system of equations is ill-posed. The posterior is maximized with a Newton optimization method. The probabilistic approach also provides estimates of the confidence in the final source map prediction. A set of adjoint, discrete ordinates flux solutions, obtained in this work by the Denovo code, is required to efficiently compute detector responses from a candidate source distribution. These adjoint fluxes form the linear mapping from the state space to the response space. The test of the method’s success is simultaneously locating a set of 137 Cs and 60 Co gamma sources in a room. This test problem is solved using experimental measurements that we collected for this purpose. Because of the weak sources available for use in the experiment, some of the expected photopeaks were not distinguishable from the Compton continuum. However, by supplanting 14 flawed measurements (out of a total of 69) with synthetic responses computed by MCNP, the proof-of-principle source mapping was successful. The locations of the sources were predicted within 25 cm for two of the sources and 90 cm for the third, in a room with an ~4-x 4-m floor plan. Finally, the predicted source intensities were within a factor of ten of their true value.
An integrated nondestructive assay (NDA) system combining active (neutron generator) and passive neutron detection and passive gamma (PG) detection is being analyzed in order to estimate the amount of plutonium, verify initial enrichment, burnup, and cooling time, and detect partial defects in a spent fuel assembly (SFA). Active signals are measured using the differential die-away (DDA), delayed neutron (DN), and delayed gamma (DG) techniques. Passive signals are measured using total neutron (TN) counts and both gross and spectral resolved gamma counts. To quantify how a system of several NDA techniques is expected to perform, all of the relevant NDA techniques listed above were simulated as a function of various reactor conditions such as initial enrichment, burnup, cooling time, assembly shuffling pattern, reactor operating conditions (including temperature, pressure, and the presence of burnable poisons) by simulating the NDA response for five sets of light water reactor assemblies. In this study, we compare the performance of several exploratory model-fitting options (including neural networks, adaptive regression with splines, iterative bias reduction smoothing, projection pursuit regression, and regression with quadratic terms and interaction terms) to relate data simulated with measurement and model error effects from various subsets of the NDA techniques to the total Pu mass. Isotope masses for SFAs and expected detector responses (DRs) for several NDA techniques are simulated using MCNP, and the DRs become inputs to the fitting process. Such responses include eight signals from DDA, one from DN, one from TN, and up to seven from PG; the DG signal will be examined separately. Results are summarized using the root-mean-squared estimation error for plutonium mass in held-out subsets of the data for a range of model and measurement error variances. Different simulation assumptions lead to different spent fuel libraries relating DRs to Pu mass. Some results for training with one library and testing with another library are also given.
In this paper, a correction-based resonance self-shielding method is developed that allows annular subdivision of the fuel rod. The method performs the conventional iteration of the embedded self-shielding method (ESSM) without subdivision of the fuel to capture the interpin shielding effect. The resultant self-shielded cross sections are modified by correction factors incorporating the intrapin effects of radial variation of the shielded cross section, radial temperature distribution, and resonance interference. A quasi–one-dimensional slowing-down equation is developed to calculate such correction factors. The method is implemented in the DeCART code and compared with the conventional ESSM and subgroup method with benchmark MCNP results. The new method yields substantially improved results for both spatially dependent reaction rates and eigenvalues for typical pressurized water reactor pin cell cases with uniform and nonuniform fuel temperature profiles. Finally, the new method is also proved effective in treating assembly heterogeneity and complex material composition such as mixed oxide fuel, where resonance interference is much more intense.
The particle transport code MCNP has been used to produce a background radiation data file on a worldwide grid that can easily be sampled as a source in the code. Location-dependent cosmic showers were modeled by Monte Carlo methods to produce the resulting neutron and photon background flux at 2054 locations around Earth. An improved galactic-cosmic-ray feature was used to model the source term as well as data from multiple sources to model the transport environment through atmosphere, soil, and seawater. A new elevation scaling feature was also added to the code to increase the accuracy of the cosmic neutron background for user locations with off-grid elevations. Furthermore, benchmarking has shown the neutron integral flux values to be within experimental error.
This work is devoted to the analysis of the neutronic characterization of the ENEA-Casaccia center RSV TAPIRO fast neutrons source reactor performed during the 80's in the frame of a joint ENEA-SCK-CEN Mol collaboration. The Belgian scientist Albert Fabry, who was available to investigate a fast research reactor to perform a full characterization of the system neutron flux, led the TAPIRO experimental campaign. In particular, the paper is focused on the comparison between some 'recommended' Fabry experimental reaction rates traverses and nowadays MCNP Monte Carlo and ERANOS deterministic calculation results. It appears that the quality of the results depends on the penetration depth into the copper reflector. The influence of copper nuclear data are more evident in the ERANOS results. Massimo Salvatores has provided a great contribution to the interpretation of the TAPIRO neutronics: this humble work is dedicated to him.
Ensuring system safety in the design, licensing, and operation phases is a priority in the nuclear industry. Performing extensive, full-scale reactor safety experiments is often prohibitive due to the large associated costs. Computational simulations offer an alternative safety analysis method, typically with significant cost reductions. Recent high-level developments in technology and increased availability of computational resources have allowed the development of high-fidelity, high-resolution multi-physics coupled codes. Such developments may be used to generate reference models for deterministic core calculations. Under this framework, the high-fidelity continuous energy Monte Carlo-based neutron transport code MCNP6.2 was coupled externally with the state-of-the-art thermal-hydraulics subchannel code, CTF4.0 for VVER (Water-Water Energetic Reactor) applications. A VVER-1000 fuel assembly model was used to demonstrate the capability of the coupled code. The converged coupled solution is compared to initial results, consisting of the first MCNP evaluation and first CTF evaluation after initialization. The VVER-1000 type assembly results are compared to other evaluations of the same assembly model. The findings indicate good agreement with expectations and reference cases, where available. The initial results of the coupled MCNP6.2/CTF4.0 calculations at the assembly level for steady-state calculations are presented in this study, which may support future work toward high-fidelity coupled full-core modeling capabilities. The multi-physics model may be further improved, enhanced, and expanded for both cycle depletion and transient applications. Such a multi-physics system will also be applicable to the VVER-1200 and other triangular lattice designs and support their deployment and operation safely and economically. (authors)
The importance of high-quality thermal neutron scattering kernels coupled with a lack of accurate sub-thermal total cross section data prompted researchers at Rensselaer Polytechnic Institute (RPI) to design a cold moderation system to enhance sub-thermal neutron flux (defined as flux below 5 meV), for nuclear data measurements at the RPI Gaerttner LINAC facility. Using a series of neutronic and heat transfer simulations obtained with MCNP and COMSOL, an optimized final design was produced. The final design exceeded current target sub-thermal neutron producing capabilities by a factor of 8 in the sub-thermal region and is easily configurable with existing neutron producing targets.
A parametric study was simulated in MCNP to investigate the impact of spontaneous fission neutron energy on neutron detector response. A total of 880 configurations were simulated by varying 3 He pressure and polyethylene reflector thickness for 8 different isotopes. Results showed that deviations between the simulated isotopes become more prevalent at higher energies with 238 U showing the greatest deviation. It was also observed that configurations comparable to most detection setups, 10 atm and 6 - 8 cm polyethylene, were closer to the upper bound of deviations from a 252 Cf. 240 Pu, a common spontaneous fission source of interest, ranged from 2% - 3% within this space. This suggests that further investigation is required to determine exact correction factors when calibrating neutron detectors with a 252 Cf. Although demonstrated here for a 4π-detection system, future plans include simulating a more detailed model of the MC-15 detection setup to determine exact correction factors for future experiments.
IRSN (France), LLNL (USA) and ORNL (USA) began a long-term collaboration effort in 2015 to update the nuclear criticality Slide Rule for the emergency response to a nuclear criticality accident. The Slide Rule permits the estimation of neutron and gamma dose rates and integrated doses based upon estimated fission yields, as a function of distance from the fission source, and time after criticality accidents for different critical systems. This paper presents results from the fourth phase of the current update of the Slide Rule project, in which delayed fission-product gamma (DFG) dose rates of unreflected plutonium critical systems were compared by several modern 3D radiation transport codes (MCNP, COG, SCALE), using updated flux-to-dose conversion factors. Dose rates are calculated for fissile material at five moderation ratios (H/Pu), at 1 m above the ground as a function of distance (between 30 cm and 1,200 m) from the external surface of the source to the center of the detector, and for periods between 1 s and 1,000 min after the critical instantaneous event. Further efforts have been devoted to the determination of the delayed gamma source, by comparing the time-dependent energy spectra obtained from several methods. Overall, DFG dose rates calculated by each participant led to consistent results. Extra effort is under way to identify the cause of the remaining differences, by comparing precisely the gamma source, and particularly nuclides inventories.
MuSim is a new user-friendly program designed to interface to many different particle simulation codes, regardless of their data formats or geometry descriptions. It presents the user with a compelling graphical user interface that includes a flexible 3-D view of the simulated world plus powerful editing and drag-and-drop capabilities. All aspects of the design can be parametrized so that parameter scans and optimizations are easy. It is simple to create plots and display events in the 3-D viewer (with a slider to vary the transparency of solids), allowing for an effortless comparison of different simulation codes. Simulation codes: G4beamline, MAD-X, and MCNP; more coming. Many accelerator design tools and beam optics codes were written long ago, with primitive user interfaces by today's standards. MuSim is specifically designed to make it easy to interface to such codes, providing a common user experience for all, and permitting the construction and exploration of models with very little overhead. For today's technology-driven students, graphical interfaces meet their expectations far better than text-based tools, and education in accelerator physics is one of our primary goals.