This document summarily provides brief descriptions of the MELCOR code enhancement made between code revision number 18019and 21440. Revision 18019 represents the previous official code release; therefore, the modeling features described within this document are provided to assist users that update to the newest official MELCOR code release, 21440. Along with the newly updated MELCOR Users’ Guide [2] and Reference Manual [3], users are aware and able to assess the new capabilities for their modeling and analysis applications.
This report represents an assessment of the gaps in Mechanistic Source Term (MST) and consequence assessment modeling for Molten Salt Reactors (MSRs). The current capabilities for MELCOR and the MELCOR Accident Code System (MACCS) are discussed, along with updates needed in order to address specific needs for MSRs. A test plan developed by Argonne National Laboratories is discussed as addressing some of these gaps, while some will require additional attention. Further recommendations are made on addressing these gaps. This report satisfies the DOE NE Milestone M2RD-21SN0601061 to leverage MELCOR and MACCS to identify parameters of importance for source term assessments for salt spill experiments.
This report aids in the development of models to perform characterization studies of aerosol dispersal and deposition within a spent fuel cask system. Due to the complex geometry in a spent-fuel canister, direct simulation of buoyancy-driven flow through the fuel assemblies to model aerosol deposition within the fuel canister is computationally expensive. Identification of an effective permeability as given in this work for a nuclear fuel assembly greatly simplifies the requirements for thermal hydraulic computations. The results of computations performed using OpenFOAM® to solve the Navier-Stokes Equations for laminar flow are used to determine an effective permeability by applying Darcy's Law. The computations are validated against an analytical solution for the special case of an infinite array of pins for which the numerical and analytical solutions have excellent agreement. The effective permeability of a 1717 PWR nuclear fuel assembly in a basket without spacer grids is numerically determined to be 1.85010 -6 m 2 for the range of fluid viscosities and pressure drops expected in a spent fuel storage canister. However, the flow is not uniform on the scale of multiple pins. Instead, significantly higher velocities are attained in the space between the assembly and the basket walls compared to the flow between the fuel pins within the assembly. Comparison with an analytical solution for fully developed flow through an infinite array of pins shows that the larger spacing near the basket walls results in about a 20% larger permeability compared to the analytical solution which does not include the enhanced flow in the space between the assembly and basket wall, or entrance and exit effects. A preliminary assessment of turbulence effects shows that with a k-epsilon model, significantly higher flow velocities are attained between the fuel pins within the assembly compared to the flow velocity in the space between the assembly and the basket walls. This is the opposite of what is determined for laminar flow.
This report documents the progress in preparing the ANSYS/Fluent® and MELCOR models to perform characterization studies of aerosols dispersal and deposition within a select spent fuel cask system. Steady state thermal-hydraulic cask response is modeled with both codes, at present, while the MELCOR source code is being modified to allow imposed thermal-hydraulic conditions with aerosol physics calculations. This will allow the MELCOR model to assume the thermal-hydraulic calculation from the ANSYS/Fluent®, while only computing the aerosol physics. Detailed results are presented on the thermal-hydraulic analysis of the MAGNASTOR® cask for the current ANSYS/Fluent® model, with convergent conditions observed over two fidelities. While the MELCOR computation computes the steady state conditions, they differ sufficiently from that computed with ANSYS/Fluent®. The MELCOR analysis include a set of UO2 sources to investigate system response for interim reporting. The airborne concentrations and evolving distributions are presented. Model development is anticipated to continue as additional components impacting flow conditions or aerosol deposition within the MAGNASTOR® cask and Westinghouse fuel assemblies are identified.
This report represents completion of milestone deliverable M2SF-21SN010309012 “Annual Status Update for OWL and Waste Form Characteristics” that provides an annual update on status of fiscal year (FY 2020) activities for the work package SF-20SN01030901 and is due on January 29, 2021. The Online Waste Library (OWL) has been designed to contain information regarding United States (U.S.) Department of Energy (DOE)-managed (as) high-level waste (DHLW), spent nuclear fuel (SNF), and other wastes that are likely candidates for deep geologic disposal, with links to the current supporting documents for the data (when possible; note that no classified or official-use-only (OUO) data are planned to be included in OWL). There may be up to several hundred different DOE-managed wastes that are likely to require deep geologic disposal. This draft report contains versions of the OWL model architecture for vessel information (Appendix A) and an excerpt from the OWL User’s Guide (Appendix B and SNL 2020), which are for the current OWL Version 2.0 on the Sandia External Collaboration Network (ECN).
One of the objectives of the United States (U.S.) Department of Energy's (DOE) Office of Nuclear Energy's Spent Fuel and Waste Science and Technology Campaign is to better understand the technical basis, risks, and uncertainty associated with the safe and secure disposition of spent nuclear fuel (SNF) and high-level radioactive waste. Commercial nuclear power generation in the U.S. has resulted in thousands of metric tons of SNF, the disposal of which is the responsibility of the DOE (Nuclear Waste Policy Act 1982). Any repository licensed to dispose the SNF must meet requirements regarding the longterm performance of that repository. For an evaluation of the long-term performance of the repository, one of the events that may need to be considered is the SNF achieving a critical configuration. Of particular interest is the potential behavior of SNF in dual-purpose canisters (DPCs), which are currently being used to store and transport SNF but were not designed for permanent geologic disposal. A two-phase study has been initiated to begin examining the potential consequences, with respect to longterm repository performance, of criticality events that might occur during the postclosure period in a hypothetical repository containing DPCs. Phase I, a scoping phase, consisted of developing an approach intended to be a starting point for the development of the modeling tools and techniques that may eventually be required either to exclude criticality from or to include criticality in a performance assessment (PA) as appropriate; Phase I is documented in Price et al. (2019). The Phase I approach guided the analyses and simulations done in Phase II to further the development of these modeling tools and techniques as well as the overall knowledge base. The purpose of this report is to document the results of the analyses conducted during Phase II. The remainder of Section 1 presents the background, objective, and scope of this report, as well as the relevant key assumptions used in the Phase II analyses and simulations. Subsequent sections discuss the analyses that were conducted (Section 2), the results of those analyses (Section 3), and the summary and conclusions (Section 4). This report fulfills the Spent Fuel and Waste Science and Technology Campaign deliverable M2SF-20SN010305061.
This report represents completion of milestone deliverable M2SF-19SNO10309013 "Online Waste Library (OWL) and Waste Forms Characteristics Annual Report" that reports annual status on fiscal year (FY) 2019 activities for the work package SF-19SN01030901 and is due on August 2, 2019. The online waste library (OWL) has been designed to contain information regarding United States (U.S.) Department of Energy (DOE)-managed (as) high-level waste (DHLW), spent nuclear fuel (SNF), and other wastes that are likely candidates for deep geologic disposal, with links to the current supporting documents for the data (when possible; note that no classified or official-use-only (OUO) data are planned to be included in OWL). There may be up to several hundred different DOE-managed wastes that are likely to require deep geologic disposal. This annual report on FY2019 activities includes evaluations of waste form characteristics and waste form performance models, updates to the OWL development, and descriptions of the management processes for the OWL. Updates to the OWL include an updated user's guide, additions to the OWL database content for wastes and waste forms, results of the beta testing and changes implemented from it. Also added are descriptions of the management/control processes for the OWL development, version control, and archiving. These processes have been implemented as part of the full production release of OWL (i.e., OWL Version 1.0), which has been developed on, and will be hosted and managed on, Sandia National Laboratories (SNL) systems. The version control/update processes will be implemented for updates to the OWL in the future. Additionally, another process covering methods for interfacing with the DOE SNF Database (DOE 2007) at Idaho National Laboratory on the numerous entries for DOE-managed SNF (DSNF) has been pushed forward by defining data exchanges and is planned to be implemented sometime in FY2020. The INL database is also sometimes referred to as the Spent Fuel Database or the SFDB, which is the acronym that will be used in this report. Once fully implemented, this integration effort will serve as a template for interfacing with additional databases throughout the DOE complex.
The U.S. Department of Energy is conducting research and development on generic concepts for disposal of spent nuclear fuel and high-level radioactive waste in multiple lithologies, including salt, crystalline (granite/metamorphic), and argillaceous (clay/shale) host rock. These investigations benefit greatly from international experience gained in disposal programs in many countries around the world. The focus of this study is the post-closure degradation and radionuclide-release rates for tristructural-isotropic (TRISO) coated particle spent fuels for various generic geologic repository environments.1,2,3 The TRISO particle coatings provide safety features during and after reactor operations, with the SiC layer representing the primary barrier. Three mechanisms that may lead to release of radionuclides from the TRISO particles are: (1) helium pressure buildup4 that may eventually rupture the SiC layer, (2) diffusive transport through the layers (solid-state diffusion in reactor, aqueous diffusion in porous media at repository conditions), and (3) corrosion5 of the layers in groundwater/brine. For TRISO particles in a graphite fuel element, the degradation in an oxidizing geologic repository was concluded to be directly dependent on the oxidative corrosion rate of the graphite matrix4, which was analyzed as much slower than SiC layer corrosion processes. However, accumulated physical damage to the graphite fuel element may decrease its post-closure barrier capability more rapidly. Our initial performance model focuses on the TRISO particles and includes SiC failure from pressure increase via alpha-decay helium, as exacerbated by SiC layer corrosion5. This corrosion mechanism is found to be much faster than solid-state diffusion at repository temperatures but includes no benefit of protection by the other outer layers, which may prolong lifetime. Our current model enhancements include constraining the material properties of the layers for porous media diffusion analyses. In addition to evaluating the SiC layer porosity structure, this work focuses on the pyrolytic carbon layers (inner/outer-IPyC/OPyC) layers, and the graphite compact, which are to be analyzed with the SiC layer in two modes: (a) intact SiC barrier until corrosion failure and (b) SiC with porous media transport. Our detailed performance analyses will consider these processes together with uncertainties in the properties of the layers to assess radionuclide release from TRISO particles and their graphite compacts.
TRISO nuclear fuel particles that are less than 1 mm in diameter are designed with multiple barrier layers to retain fission products both during reactor operations and for long-term geological disposal. The primary barrier is a 35 lam thick silicon carbide (SiC-a highly impermeable semi-metal) layer for which data are available on the diffusion of short-lived fission products at high temperatures (> 1000 degC). However, for a geological repository, this layer may contact brine and hence corrode even at ambient temperatures. As an initial approach to assess the effectiveness of the SiC barrier for geological repositories, ranges of fission product diffusivities and corrosion rates for SiC are modeled concurrently with the simultaneous effect of radioactive decay. Using measured corrosion rates of SiC, if the diffusivity is more than about 10 -20 m 2 /s, fission product releases may occur before the SiC barrier has corroded to the point of breach. For diffusivities less than about 10 -21 m 2 /s there may not be significant diffusional releases prior to SiC barrier removal/breach by corrosion. This work shows the importance of estimating diffusivities in SiC at geological repository temperatures, and highlights the relevance of evaluating the porosity/permeability evolution of the SiC layer in a geologic environment.
This report represents completion of milestone deliverable M2SF-18SNO10309013 "Inventory and Waste Characterization Status Report and OWL Update that reports on FY2018 activities for the work package (WP) SF-18SNO1030901. This report provides the detailed final information for completed FY2018 work activities for WP SF-18SN01030901, and a summary of priorities for FY2019. This status report on FY2018 activities includes evaluations of waste form characteristics and waste form performance models, updates to the OWL development, and descriptions of the two planned management processes for the OWL. Updates to the OWL include an updated user's guide, additions to the OWL database content for wastes and waste forms, results of the Beta testing and changes implemented from it. There are two processes being planned in FY2018, which will be implemented in FY2019. One process covers methods for interfacing with the DOE SNF DB (DOE 2007) at INL on the numerous entries for DOE managed SNF, and the other process covers the management of updates to, and version control/archiving of, the OWL database. In FY2018, we have pursued three studies to evaluate/redefine waste form characteristics and/or performance models. First characteristic isotopic ratios for various waste forms included in postclosure performance studies are being evaluated to delineate isotope ratio tags that quantitatively identify each particular waste form. This evaluation arose due to questions regarding the relative contributions of radionuclides from disparate waste forms in GDSA results, particularly, radionuclide contributions of DOE-managed SNF vs HLW glass. In our second study we are evaluating the bases of glass waste degradation rate models to the HIP calcine waste form. The HIP calcine may likely be a ceramic matrix material, with multiple ceramic phases with/without a glass phase. The ceramic phases are likely to have different degradation performance from the glass portion. The distribution of radionuclides among those various phases may also be a factor in the radionuclide release rates. Additionally, we have an ongoing investigation of the performance behavior of TRISO particle fuels and are developing a stochastic model for the degradation of those fuels that accounts for simultaneous corrosion of the silicon carbide (SiC) layer and radionuclide diffusion through it. The detailed model of the TRISO particles themselves, will be merged with models of the degradation behavior(s) of the graphite matrix (either prismatic compacts or spherical "pebbles") containing the particles and the hexagonal graphite elements holding the compacts.
Safety basis analysts throughout the U.S. Department of Energy (DOE) complex rely heavily on the information provided in the DOE Handbook, DOE - HDBK - 3010, Airborne Release Fractions/Rates and Respirable Fractions for Nonreactor Nuclear Facilities, to determine radionuclide source terms. In calculating source terms, analysts tend to use the DOE Handbook's bounding values on airborne release fractions (ARFs) and respirable fractions (RFs) for various categories of insults (representing potential accident release categories). This is typically due to both time constraints and the avoidance of regulatory critique. Unfortunately, these bounding ARFs/RFs represent extremely conservative values. Moreover, they were derived from very limited small-scale bench/laboratory experiments and/or from engineered judgment. Thus, the basis for the data may not be representative of the actual unique accident conditions and configurations being evaluated. The goal of this research is to develop a more accurate and defensible method to determine bounding values for the DOE Handbook using state-of-art multi-physics-based computer codes. This enables us to better understand the fundamental physics and phenomena associated with the types of accidents in the handbook. In this year, this research included improvements of the high-fidelity codes to model particle resuspension and multi-component evaporation for fire scenarios. We also began to model ceramic fragmentation experiments, and to reanalyze the liquid fire and powder release experiments that were done last year. The results show that the added physics better describes the fragmentation phenomena. Thus, this work provides a low-cost method to establish physics-justified safety bounds by taking into account specific geometries and conditions that may not have been previously measured and/or are too costly to perform.
In this work we have presented a particle resuspension model implemented in the SNL code SIERRA/Fuego, which can be used to model particle dispersal and resuspension from surfaces. The method demonstrated is applicable to a class of particles, but would require additional parametric fits or physics models for extension to other applications, such as wetted particles or walls. We have demonstrated the importance of turbulent variations in the wall shear stress when considering resuspension, and implemented both shear stress variation models and stochastic resuspension models (not shown in this work). These models can be used in simulations with of physically realistic scenarios to augment lab-scale DOE Handbook data for airborne release fractions and respirable fractions in order to provide confidences for safety analysts and facility designers to apply in their analyses at DOE sites. Future work on this topic will involve validation of the presented model against experimental data and extension of the empirical models to be applicable to different classes of particles and surfaces.