A systematic approach to identify initiating events and its relationship to Probabilistic Risk Assessment: Demonstrated on the Molten Salt Reactor Experiment
Not Available
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.
Not Available
Abstract Characterizing oxide nuclear fuels is difficult due to complex fission products, which result from time-evolving system chemistry and extreme operating environments. Here, we report a machine learning-enhanced approach that accelerates the characterization of spent nuclear fuels and improves the accuracy of identifying nanophase fission products and bubbles. We apply this approach to commercial, high-burnup, irradiated light-water reactor fuels, demonstrating relationships between fission product precipitates and gases. We also gain understanding of the fission versus decay pathways of precipitates across the radius of a fuel pellet. An algorithm is provided for quantifying the chemical segregation of the fission products with respect to the high-burnup structure, which enhances our ability to process large amounts of microscopy data, including approaching the atomistic-scale. This may provide a faster route for achieving physics-based fuel performance modeling.
Abstract Internal kink (IK) instability is investigated for European demonstration fusion reactor (EU DEMO) plasmas in both negative triangularity (NT) and positive triangularity (PT) configurations. For NT plasmas, the IK becomes more unstable as an ideal conformal wall moves away from the plasma boundary, with the mode growth rate saturating at the wall radial location of about b / a = 1.5 , where a is the plasma minor radius and b the wall radial location. The plasma resistivity destabilizes the IK mode. The effect of sub-sonic toroidal plasma flow is sufficiently weak and can thus be ignored for these EU DEMO equilibria. These results are consistent with those for PT plasmas, albeit with larger mode growth rate in the NT configuration. Both perturbative and self-consistent magneto-hydrodynamic (MHD)-kinetic hybrid calculations predict (partial) stabilization of the IK modes in both NT and PT configurations, with inclusion of various kinetic contributions. Precessional drift motion of trapped fusion-born alphas in EU DEMO produces weak stabilization to the IK mode. Stronger stabilization occurs with the toroidal precession of trapped thermal particles (ions and electrons) and the bounce-transit motion of thermal ions. The stabilization is similar between the NT and PT configurations, due to the similarity of the mode eigenfunction (occupying a nearly circular region in the plasma core) despite the sign difference in the triangularity. The non-perturbative MHD-kinetic hybrid model predicts much less stabilization of the mode than the perturbative model, primarily due to the self-consistent determination of the mode eigenvalue in the former. Generally, no significant difference in the IK mode stability is found between the NT and PT plasmas in EU DEMO.
The understanding of the physics underlying the L-H transition has strong implications for ITER experimental reactor and demonstration power plant (DEMO). In many tokamaks, including JET, it has been observed that, at a particular plasma density, n e,min , the power necessary to access H-mode P L-H is minimum. In the present work, L-H transitions of JET deuterium plasmas heated by neutral beam injection (NBI) are studied for the first time by means of a power balance analysis to characterize the main contributions in the transition, through integrated transport modelling. In the pulses analysed, we do observe a minimum of the L-H power threshold in density, indicating the presence of density branches and of n e,min . Electron and ion heat fluxes at the transition are estimated separately. The electron/ion equipartition power results in favour of the ions, as shown by QuaLiKiz quasilinear gyrokinetic simulations, which predict a larger ion transport that causes T e > T i . The resulting edge ion heat flux also shows a clear change of slope below n e,min , similarly to ASDEX-Upgrade (AUG) NBI pulses (Ryter et al 2014 Nucl. Fusion 54 083003). JET NBI data are compared to radio-frequency heated AUG and Alcator C-mod pulses (Schmidtmayr et al 2018 Nucl. Fusion 58 056003), showing a different trend of the power, coupled to ions at the L-H transition with respect to the linearity observed in the radio-frequency heated plasmas. Finally, the presence of n e,min and the role of the ion heat flux is discussed in the paper, although it seems it is not possible to explain the presence of a P L-H minimum in density by a critical ion heat flux and by the equipartition power for the JET NBI-heated plasmas analysed.
An optimized plasma current ramp-down strategy is critical for safe and fast termination of plasma discharges in a tokamak demonstration fusion reactor (DEMO), both in planned and emergency scenarios, avoiding plasma disruptions and excessive heat loads to the first wall. Plasma stability limits and machine-specific technical requirements constrain the stable envelope through which the plasma must be navigated. Large amounts of auxiliary heating are required throughout the ramp-down phase, to avoid a radiative collapse in the presence of intrinsic tungsten and seeded xenon impurities, as quantitatively estimated in this work. As the plasma current is reduced, the current density becomes increasingly peaked, reflected by a growing value of the internal inductance $\ell$ $i3$ , resulting in reduced controllability of the vertical position of the plasma. The feasibility of different plasma current ramp-down rates is tested by applying an automated optimization framework embedding the RAPTOR core transport solver. Optimal time traces for plasma current $I$ $p$ $(t)$ and plasma elongation $κ(t)$ are proposed, to satisfy an $I$ $p$ -dependent upper limit on the plasma internal inductance, as obtained from vertical stability studies using the CREATE-NL code, as well as a constraint on the time evolution of $q$ 95 , to avoid an ideal MHD mode. A negative current density near the plasma edge is observed in our simulations, even for the most conservative $I$ $p$ ramp-down rate, indicating significant transient dynamics due to a large resistive time.
The Korean Superconducting Tokamak Advanced Research has been focused on exploring the key physics and engineering issues for future fusion reactors by demonstrating the long pulse operation of high beta steady-state discharge. Advanced scenarios are being developed with the goal for steady-state operation, and significant progress has been made in high ℓ i , hybrid and high beta scenarios with β N of 3. In the new operation scenario called fast ion regulated enhanced (FIRE), fast ions play an essential role in confinement enhancement. GK simulations show a significant reduction of the thermal energy flux when the thermal ion fraction decreases and the main ion density gradient is reversed by the fast ions in FIRE mode. Optimization of 3D magnetic field techniques, including adaptive control and real-time machine learning control algorithm, enabled long-pulse operation and high-performance ELM-suppressed discharge. Symmetric multiple shattered pellet injections (SPIs) and real-time disruption event characterization and forecasting are being performed to mitigate and avoid the disruptions associated with high-performance, long-pulse ITER-like scenarios. Finally, the near-term research plan will be addressed with the actively cooled tungsten divertor, a major upgrade of the NBI and helicon current drive heating, and transition to a full metallic wall.
As several advanced reactor concepts are maturing, test beds are needed to accelerate the demonstration and deployment of these advanced nuclear technologies. The National Reactor Innovation Center (NRIC) is building new or enhancing existing US Department of Energy infrastructure to support testing of components and systems. Demonstration of Microreactor Experiments (DOME) will utilize the Experimental Breeder Reactor-II (EBR-II) dome containment structure to host reactor demonstrations. A reactor supplemental shielding is needed so that DOME dose requirements are met. To accelerate the confirmatory analysis required for the reactor demonstration, the NRIC Virtual Test Bed (VTB) is developing a virtual model of the DOME shield that will be made available on the VTB public repository. This will allow developers to leverage advanced modeling and simulation tools to ensure their reactor demonstration concept will meet dose requirement and the limit concrete temperature in the shield during steady state and transient operation conditions. This paper presents the model developed for the DOME shield using open-source tools: MOOSE heat transfer module, Monte Carlo code OpenMC, and Cardinal to calculate the DOME shield temperature distribution during steady state
The Department of Energy (DOE) Office of Nuclear Energy National Reactor Innovation Center accelerates the deployment of novel reactor concepts by establishing both physical and virtual spaces for building and testing various components, systems, and complete pilot plants. The Virtual Test Bed represents the virtual arm of the National Reactor Innovation Center and is a joint effort with the DOE Nuclear Energy Advanced Modeling and Simulation Program. The Virtual Test Bed mission is to accelerate the deployment of advanced reactors by facilitating the adoption of cutting-edge DOE advanced modeling and simulation tools to design, evaluate, and license reactors. This is primarily achieved by storing example challenge problems in an externally available repository and by developing models to fill the M&S gaps needed for potential demonstrators. Activities conducted this fiscal year focused on developing of a Demonstration of Microreactor Experiments shield model to help accelerate the confirmatory analysis required for the reactor demonstration. This model and workflow will allow developers to leverage advanced modeling and simulation tools to ensure their reactor demonstration concept will meet dose requirements and that the surrounding shield will stay within concrete temperature limits during steady-state and transient operation conditions. An initial model has been developed to evaluate the temperature distribution in the concrete shield during steady-state operation, including neutron and gamma heating effects. Various modeling strategies have been examined to understand their applicability and limitations with different reactor designs to make the workflow as reactor-agnostic as possible and computationally effective to maximize its usability. In addition to describing the Demonstration of Microreactor Experiments shield model and associated results, this report summarizes other accomplishments regarding repository maintenance and improvement and new external models hosted on the repository.
In support of the US Nuclear Regulatory Commission non-light-water reactor fuel cycle demonstration project, the capabilities of SCALE 6.3.1 for radionuclide characterization, criticality, and shielding were demonstrated through several scenarios in a molten salt reactor (MSR) nuclear fuel cycle. Three scenarios were selected for this study, with the molten salt breeder reactor (MSBR) serving as the reference design. The original thorium-based fuel salt used in MSBR was replaced with a 235 U-enriched fuel salt to better reflect anticipated future MSR concepts. The first scenario focused on the fuel salt preparation stage. Criticality analyses, using SCALE/CSAS-Shift and SCALE/SAMPLER, examined potential critical configurations in the fresh fuel salt container based on variables such as 235 U enrichment, UF 4 molar fraction, temperature, and container geometry. The second scenario explored the release of fission products during reactor operation in which the tritium buildup in the primary fuel salt was quantified through depletion calculations using SCALE/TRITON with continuous fuel makeup and fractional fission product removal. This scenario also assessed radioactivity levels in the off-gas system (OGS), quantifying radiotoxic nuclides that could be released into the environment. SCALE/SAMPLER was employed to perform uncertainty quantification to identify key operating parameters influencing the buildup of radiotoxic nuclides in the OGS. Additionally, radiation dose rates were calculated using SCALE/MAVRIC to estimate the maximum exposure levels within the reactor cell complex during operation. The final scenario evaluated the decay heat and radioactivity of the primary fuel salt collected in the drain tank, as well as the radiation dose rate outside the building from the presence of the drain tank. Overall, the application of SCALE’s various capabilities for radionuclide inventory generation, criticality analysis, and shielding in MSR fuel cycles was successfully demonstrated across these scenarios.
MELCOR is an integrated thermal hydraulics, accident progression, and source term code for reactor safety analysis that has been developed at Sandia National Laboratories for the United States Nuclear Regulatory Commission (NRC) since the early 1980s. Though MELCOR originated as a light water reactor (LWR) code, development and modernization efforts have expanded its application scope to includ e non-LWR reactor concepts. Current MELCOR development efforts include providing the NRC with the analytical capabilities to support regulatory readiness for licensing non-LWR techno logies under Strategy 2 of the NRC?s near- term Implementation Action Plans. Beginning with the Next Generation Nuclear Project (NGNP), MELCOR has undergone a range of enha ncements to provide analytical capabilities for modeling the spectrum of advanced non-LWR concepts. This report describes the generic plant model developed to demonstrate MELCOR capabilities to perform heat pipe reactor (HPR) safety evaluations. The generic plant mode l is based on a publicly-available Los Alamos National Laboratory (LANL) Megapower design as modified in the Idaho National Laboratory (INL) Design A description. For plant aspects (e.g., reactor building size and leak rate) that are not described in the LANL and INL references , the analysts made assumptions needed to construct a MELCOR full-plant model. The HP R uses high assay, low-enrichment uranium (HALEU) fuel with steel cladding that uses heat pipes to transfer heat to a secondary Brayton air cycle. The core region is surrounded by a stainless-steel shroud, alumina reflector, core barrel and boron carbide neutron shield. The reactor is secured inside a below-grade cavity, with the operating floor located above the cavity. Example calculations are performed to show the plant response and MELCOR capabilities to characterize a range of accident conditions. The accidents selected for evaluation consider a range of degraded and failed modes of operation for key safety functions providing re activity control, the primary and secondary system heat removal, and the effectiveness of th e confinement natural circulation flow into the reactor cavity (i.e., a flow blockage).
In support of the US Nuclear Regulatory Commission non–light-water reactor fuel cycle demonstration project, SCALE 6.3.1 capabilities for radiation source term and shielding calculations are demonstrated for scenarios in the sodium-cooled fast reactor (SFR) fuel cycle. A postulated accident scenario, which consists of a seismic event causing the refueling machine to fall and release a spent fuel assembly inside the containment building (CB), is analyzed in this paper. Radiation source terms were generated for a U/TRU-10Zr metal fuel assembly with a 16.5% initial transuranic waste content and a discharge burnup of approximately 95 GWd/tHM; source terms were also generated for a high-assay low-enriched uranium metal fuel assembly (U-10Zr) with a 16.5% initial enrichment and a discharge burnup of 149.74 GWd/tHM. These radiation source terms were then used to determine the dose rate inside the CB and near the outer surface of the CB for a range of spent fuel assembly cooling times. The dose rate produced by the analyzed SFR assemblies is similar to that produced by a typical pressurized water reactor assembly with a discharge burnup of 50 GWd/MTU. Ultimately, the validation of the source terms predicted for SFRs with SCALE will need to be demonstrated via the use of assay measurements.
Solar and other renewable energy driven gas-solid thermochemical energy storage (TCES) technology is a promising solution for the next generation energy storage systems due to its high operating temperature, efficient energy conversion, ultra-long storage duration, and potential high energy density. Experimental and theoretical studies suggest that the respective gravimetric and volumetric TCES energy storage densities vary from 200 to 3000 kJ kg –1 and 1–3 GJ m –3 . Solar radiation or heat generated from electric furnaces powered by renewable electricity can be stored in the form of chemical energy through endothermic reactions, while the stored chemical energy can be converted to thermal energy via an exothermic reaction when needed. The design of highly effective reactors requires a deep understanding of materials, thermodynamics, chemical kinetics, and transport phenomena. At time of writing, TCES reactors are yet to be deployed at commercially relevant scales, leaving a substantial gap between development efforts and commercial feasibility. Therefore, this review aims to examine the state-of-the-art design and performance of particle-based TCES reactors with different reactive materials. Fundamentals related to TCES reactive materials, reaction conditions, thermodynamics and kinetics, and transport phenomena are reviewed in detail to provide a comprehensive understanding of the reactor design and operation. Five major types of TCES reactors have been comprehensively reviewed and compared, including fixed, moving, rotary, fluidized, and entrained bed reactors. Most reported prototype reactors in the literature operate at lab scale with thermal inputs below 40 kW, and scaled TCES reactors (e.g., at megawatt level) are yet to be demonstrated. The nominal reactor operating temperatures range from 300 to 1500 °C, depending on the selected chemistry, reactive material, and heat sources. To evaluate their designs, the reactors are assessed in aspects of performance, cost, and durability. Discrepancies in performance indicators of energy storage density, extent of reaction, and various energy efficiencies are highlighted. The scale-up of reactors and power block integration, which hold the key to the successful commercialization of TCES systems, are critically analyzed. Furthermore, advanced materials (both reactive materials and ceramic reactor housing materials), effective particle flow control, advanced modeling tools, and novel system design may bring significant improvement to the energy efficiency, storage density and cost competitiveness of particle-based TCES reactors.
We report the precision reactor oscillation and spectrum experiment, PROSPECT, has made world-leading measurements of reactor antineutrinos at short baselines. In its first phase, conducted at the high flux isotope reactor (HFIR) at Oak Ridge National Laboratory, PROSPECT produced some of the strongest limits on eV-scale sterile neutrinos, made a precision measurement of the reactor antineutrino spectrum from 235 U, and demonstrated the observation of reactor antineutrinos in an aboveground detector with good energy resolution and well-controlled backgrounds. The PROSPECT collaboration is now preparing an upgraded detector, PROSPECT-II, to probe yet unexplored parameter space for sterile neutrinos and contribute to a full resolution of the reactor antineutrino anomaly, a longstanding puzzle in neutrino physics. By pressing forward on the world’s most precise measurement of the 235 U antineutrino spectrum and measuring the absolute flux of antineutrinos from 235 U, PROSPECT-II will sharpen a tool with potential value for basic neutrino science, nuclear data validation, and nuclear security applications. Following a two-year deployment at HFIR, an additional PROSPECT-II deployment at a low enriched uranium reactor could make complementary measurements of the neutrino yield from other fission isotopes. PROSPECT-II provides a unique opportunity to continue the study of reactor antineutrinos at short baselines, taking advantage of demonstrated elements of the original PROSPECT design and close access to a highly enriched uranium reactor core.
The MOOSE-based reactor physics code Griffin was assessed on a heterogeneous pin-resolved model of a prototype lead-cooled fast reactor assembly. This model was developed in preparation for future use in MOOSE-based multiphysics calculations for computing hot channel factors. Heterogeneous multigroup cross sections were prepared using the fast reactor multi-group cross section processing code MC{sup 2}-3 using a two-step method. Griffin simulations were performed using the DFEM-SN solver on 576 cores on Argonne's LCRC cluster. Diffusion-based acceleration methods were applied (NDA and CMFD). Reference solutions were generated with continuous energy MCNP and the hybrid MOC/finite element solver PROTEUS-MOC for code-to-code comparison. Space-angle convergence studies were conducted to observe convergence in k-eigenvalue and axial pin power distributions. The fully resolved Griffin calculation was within 68 pcm of the MCNP eigenvalue and exhibited max 1.2% relative error in the axial pin power distribution. Griffin produced nearly identical results to PROTEUS-MOC when using the same 9-group multigroup cross-section set. Griffin demonstrated favorable scaling in wall-clock time and memory usage when using diffusion-based acceleration methods. Griffin is capable of simulating the pin-resolved heterogeneous LFR assembly with good accuracy and performance, and is suitable for future use in coupled high-fidelity hot channel factor simulations. (authors)
The SHARP toolkit is a high-fidelity reactor simulation tool developed under the U.S. Department of Energy, Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Campaign. SHARP toolkit is comprised of the neutronics module PROTEUS thermal hydraulics module Nek5000, and structural mechanics module Diablo. During FY17 and FY18, the PROTEUS and Nek5000 components of SHARP were applied to solve challenging sodium-cooled fast reactor (SFR) problems. In particular, selected hot channel factors (HCF) for a prototype metal-fueled SFR design (the AFR-100) were analyzed in high fidelity, and the “SHARP zooming capability” for SFRs was developed and demonstrated to reduce computational expense for full core problems in cases where detailed data is needed in selected fuel assemblies. After the previous success applying SHARP to challenging SFR problems, the focus in FY19 and FY20 expanded to additional fast reactor applications including lead cooled fast reactors (LFR) and sodium cooled fast reactors (SFR). The specific technical tasks were (1) assessment of hot channel factors for LFR, for which no data currently exists, and (2) demonstration of zooming capability in assemblies of the Versatile Test Reactor (VTR). First-of-a-kind hot channel factor (HCF) estimation for LFR with high fidelity codes (PROTEUS/Nek5000) was successfully demonstrated in this study which began in FY19 and continued in FY20. Selected HCF were computed and compared with SFR data (AFR-100, EBR-II). The findings confirm that different reactor types, design parameters and uncertainties lead to different HCFs. Careful estimation of HCF for a specific design is necessary to obtain appropriate HCFs. In addition to improvement in HCF accuracy, high fidelity tools generate data to help the designer better understand the mechanism of the impact from these uncertainties. For example, the impact of cladding thickness manufacturing tolerance resulted in non-intuitive effects in the corner pins of the LFR assembly. This procedure of computing HCF using high fidelity models shows promise and flexibility for being repeated for any arbitrary reactor of choice. Along with the application on SFR and LFR, the capability of the tools has also been matured to deal with different reactor types and designs. Progress was made towards extending the previously demonstrated SHARP zooming capability to non-fueled SFR assemblies. In particular, in FY19 a gamma transport capability was implemented in both high fidelity PROTEUS solvers in order to accurately account for heat deposition caused by gamma particles, which accounts for ~10% of total core power. Neutronics verification cases were carried out for a candidate Versatile Test Reactor (VTR) design using the new gamma transport capability in PROTEUS. Comparisons were made with continuous energy MCNP calculations and shown to agree well. The models for the full core design with heterogeneous control and fuel assemblies is in progress for PROTEUS-SN and completed with MCNP. The MCNP power distributions were transferred to Nek5000 to perform thermal hydraulic calculations of the control and fuel assembly.
Advanced reactor concepts currently under development in the United States [e.g., the High-Temperature Gas Reactor (HTGR) of X-Energy and the Molten Chloride Fast Reactor (MCFR) of TerraPower] are striving to improve the economics of nuclear power primarily by using passively and ‘walk-away’ safe technologies in the reactor. The recently published report, ‘Future of Nuclear Energy in a Carbon-Constrained World’ by the Massachusetts Institute of Technology Energy Institute (MIT EI) has shown that a significant contributor to the capital cost of NPPs is the ‘civil works’ i.e., the construction of the balance of plant (BoP: all SSCs except those involved in power generation, e.g., reactor vessel) such as the buildings, containment dome, foundation, etc., which amount to almost half of the overnight capital cost. The report suggests that using advanced construction technologies such as seismic isolation and high-performance concrete to significantly reduce the capital costs of advanced reactors. This report demonstrates the application of two cost-cutting solutions on advanced reactor designs that are currently under development. These solutions include, (1) seismic base isolation, and (2) risk- and cost-based seismic design optimization including seismic isolation of individual components.
The U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program develops an integrated suite of advanced reactor physics tools built upon the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Each code generally requires an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The meshing process is often burdensome for the complex geometries present in reactors due to lack of easy-to-use, open-source meshing tools. To address the bottleneck associated with meshing complex geometries found in nuclear reactors, NEAMS initiated the development of the MOOSE Reactor Module starting in FY21. The Reactor Module builds off the existing MOOSE Mesh System to include targeted meshing capabilities such as the ability to generate hexagonal pin cells, assemblies with ducts, rotating control drums, cores, peripheral zones around a core, as well as the automatic labeling (“reporting IDs”) of pin, assembly, and planar zones to simplify post-processing of results. As a Physics Module in MOOSE, the Reactor Module is open-source, available with any MOOSE installation, directly compatible with MOOSE-based tools, and can be invoked from MOOSEbased applications to generate meshes. Functionality from the Reactor Module has been applied to several advanced reactor concepts to demonstrate user workflow improvements and accuracy. The primary objective of the Reactor Module is to improve useability of MOOSEbased tools by streamlining mesh generation and output inspection processes. During FY22, the functionality of the Reactor Module (and accompanying Mesh System) has been expanded based on user needs. First, the Reactor Geometry Mesh Builder capability developed primarily in FY21 has been refactored and merged to the public MOOSE repository. This capability wraps underlying Reactor Module mesh generators into a “Pin – Assembly – Core” workflow appropriate for conventional Cartesian and hexagonal geometries, and notably assigns material IDs during mesh generation stage and generates only the minimal number of blocks needed in order to reduce computational burden. Biasing and boundary layer options have been added to the base mesh generators as required by thermal hydraulics solvers. The reporting ID functionality has been expanded to differentiate ring-wise and azimuthal sectors within a pin for use with depletion algorithms, and VectorPostProcessor and Reporter objects are now available to integrate solution variables across zones based on ID combinations. Functionality to trim hexagonal meshes along the center or periphery has been developed so users may leverage symmetry and reflective boundary conditions to reduce the mesh size. A flexible and powerful tool to fill the space between two sidesets has been introduced to the framework and can be used for transition layers such as stitching two assemblies together with different numbers of pins, or for complex geometries which do not follow conventional Cartesian/hexagonal patterns. Finally, additional verification problems were performed with NEAMS physics tools in complement with existing NEAMS work.