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 is a screening spreadsheet and is used as a tool to calculate the source term for annual reporting, a new proposed project (APAD) or to calculate a periodic confirmatory measurement (PCM). NASA is doing a study of possible locations for a nuclear thermal propulsion ground test site and would like to use the spreadsheet to refine their study and provide guidance for future planning.
It has been established that incremental improvements in beyond design basis accident performance can be achieved through accident tolerant fuel (ATF). However, they have the potential to recover margin with respect to conventional fuel and therefore enhance plant economics through uprate or cycle length increase. To realize this potential, it is necessary to quantify the reduction in source term due to use of ATF, and correspondingly how this is affected by increasing the cycle length and/or burnup. This requires development of a risk informed analysis methodology for ATF under high burnup conditions, which is being developed within LWRS. To this end, a MELCOR simulation of FeCrAl ATF in a recovered Large Break LOCA (LBLOCA) scenario has been developed, using a model based on the Zion Pressurized Water Reactor (PWR). The new user defined material capability, along with the inclusion of detailed neutronics- and depletion-derived parameters such as core power distribution, decay heat behavior, and fission product inventories, allows a more detailed simulation of FeCrAl clad material properties and behavior than has previously been possible using MELCOR. These detailed FeCrAl results were compared with a zircaloy clad model to investigate the differences between the two clad types and quantify the benefits of FeCrAl ATF with an 18-month cycle. Next, the fuel cycle with ATF was extended to 24 months, to determine whether any additional safety margin provided by FeCrAl ATF could be leveraged to implement high burnup FeCrAl-clad fuel while retaining the same operating and safety limits as current zircaloy-clad fuel. For the particular scenario analyzed, the delay to fuel failure from using ATF was of the same order as the LOCA recovery time, and hence significant in reducing fission product release, with the 24-month FeCrAl core performing better than the 18-month zircaloy core. It is noted that for other transients, the reduction in release due to using FeCrAl may be less significant. Furthermore, the material model developed here can be used in such further studies in support of determining the overall source term.
Transpiration cooling has renewed interest of study as a renewable system of thermal protection for atmospheric entry systems. The presence of coolant within the void space of a porous material changes its effective porosity and theorized to lend the external heatshield as a potential means of porous surface control. Effective porosity is explored in two forms: saturation of void space where exiting coolant creates ‘blowing’ over the surface, and partial evacuation of the void space where withdrawal of the coolant creates ‘suction’ on the surface. Thus this study evaluates using transpiration cooling as a as means of active aerodynamic control. The feasibility of the system is evaluated numerically using Reynolds-averaged Navier Stokes based computational fluid dynamics for the aerodynamic assessments by coding source terms of transpiration cooling products at a heat shield surface boundary at the wind side shoulder, creating an asymmetric scheme. Study results show induced moments about the pitch axis incurred following the source terms of a cooling fluid at a near-shoulder location on the heat shield. Lift and drag saw increasing modulation with increase in mass flux rate. Pitch moment was altered a total of 388 N-m between lowest and highest ‘suction’ rates and changed 194 N-m maximum between ‘blowing’ rates.
A modification is developed for application to the source terms used in combustion modeling. The modification accounts for the error of the finite difference scheme in regions where chain-branching chemical reactions produce exponential growth of species densities. The modification is first applied to a one-dimensional scalar model problem. It is then generalized to multiple chemical species, and used in quasi-one-dimensional computations of shock-induced combustion in a channel. Grid refinement studies demonstrate the improved accuracy of the method using this modification. The algorithm is applied in two spatial dimensions and used in simulations of steady and unsteady shock-induced combustion. Comparisons with ballistic range experiments give confidence in the numerical technique and the 9-species hydrogen-air chemistry model.
A comprehensive review of the current state of knowledge surrounding molten salt aerosol behavior and bubble transport in molten salt systems relevant to mechanistic source term (MST) analyses of molten salt reactors (MSR) and fluoride high temperature reactors (FHR) was provided. Although relatively little experimental work has been carried out on investigating these phenomena in molten salts compared to that for other nuclear fluids of interest (e.g., aqueous, sodium), the amount of work that has been reported still provides a baseline understanding to inform and design future experimental campaigns on these topics. The bulk of this data is derived from the Molten Salt Reactor Program (MSRP) which was carried out at Oak Ridge National Laboratory (ORNL) nearly 50 years ago, leading to the successful operation of the Molten Salt Reactor Experiment (MSRE). Many other related research efforts are also highlighted in this review, such as recent efforts under the DOE:NE Advanced Reactor Technology (ART) MSR Campaign. The goal of this report is to provide a foundation for modelers to base assumptions on, as well as to make recommendations to experimentalists attempting to close information gaps. As such, Section 4 contains recommendations for future experimental work that may be needed to close current information gaps related to MST modeling of MSRs. Additionally, given the flowing fuel nature of MSRs, the relevance of aerosol and bubble transport behavior in molten salts goes beyond just MST modeling, evidenced by experiences with the pump bowl and off-gas system of the MSRE.
This report summarizes FY23 activities to improve mechanistic source term modeling for MSR concepts. Relevant MELCOR capability enhancements made during FY23 are summarized including development of a flexible python-based EOS generator (MELEOS), porous domain modeling capabilities for validation applications, and development of a MELCOR model for the LSTL facility in anticipation of upcoming molten salt experiments.
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 include non-LWR reactor concepts. Current MELCOR development efforts include providing the NRC with the analytical capabilities to support regulatory readiness for licensing non-LWR technologies 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 enhancements 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 fluoride-salt-cooled high-temperature reactor (FHR) safety evaluations. The generic plant model is based on publicly-available FHR design information. For plant aspects (e.g., reactor building leak rate and details of the cover-gas system) that are not described in the FHR references, the analysts made assumptions needed to construct a MELCOR full-plant model. The FHR model uses a TRi-structural ISOtropic (TRISO)-particle fuel pebble-bed reactor with a primary system rejecting heat to two coiled tube air heat ex changers. Three passive direct reactor auxiliary cooling systems provide heat removal to supplement or replace the emergency secondary system heat removal during accident conditions. Surrounding the reactor vessel is a low volume reactor cavity that insulates the reactor with fire bricks and thick concrete walls. A refractory reactor liner system provides water cooling to reduce the concrete wall temperatures. 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 reactivity control, the primary system decay heat removal and also a piping leak of the line to the coolant drain tank.
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 over the past decades have expanded its application scope to include non-LWR reactor concepts. Current MELCOR development efforts include providing the NRC with the analytical capabilities to support regulatory readiness for licensing non-LWR technologies under Strategy 2 of the NRC's near-term Implementation Action Plans. Beginning with the Next Generation Nuclear Project (NGNP), MELCOR ha s undergone a range of enhancements 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 high-temperature gas reactor (HTGR) safety evaluations. The generic plant model is based on publicly available PMBR-400 design information. For plant aspects (e.g., reactor building size and leak rate) that are not described in the PBMR-400 references, the analysts made assumptions needed to construct a MELCOR full-plant model. The HTGR model uses a TRi-structural ISOtropic (TRISO)-particle fuel pebble-bed reactor with a primary system rejecting heat to a recuperative heat exchange r. Surrounding the reactor vessel is a reactor cavity contained within a confinement room cooled by the Reactor Cavity Cooling System (RCCS). 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 reactivity control, primary system heat removal and reactor vessel decay heat removal, and confinement cooling.
To assist both the advanced reactor industry and U.S. Nuclear Regulatory Commission (NRC) in the pursuit of reactor design and licensing, the U.S. Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has established a mechanistic source term (MST) research project under the Multiphysics Applications technical area. As part of this effort, an MST mod/sim development pathway was developed in FY21 under the NEAMS program, which outlines the high-level objectives and near-term tasks necessary to achieve the project objectives. One such recommended task was the assessment of a previously developed thermodynamic database to model radionuclide vaporization from liquid sodium pools (specifically in the context of MST analyses), which is the focus of this topic report. As such, a thermodynamic database was created in FactSage for liquid sodium pools containing oxygen and relevant elements which represent radionuclides pertinent to MST analyses. The database was demonstrated in calculating vapor fraction curves for each element in the system at concentrations relevant to SFR MST.
The purpose of this document is to provide source term estimates and the associated technical basis for Strategic Capabilities Office (SCO) microreactor designs that use TRISO fuel. This estimate can be used for relevant environmental and safety analyses that will be done as part of the project for all of the selected suppliers for a safety design strategy. This is a combined quantitative and qualitative analysis.
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).
Two approaches are used to extend the essentially non-oscillatory (ENO) schemes to treat conservation laws with stiff source terms. One approach is the application of the Strang time-splitting method. Here the basic ENO scheme and the Harten modification using subcell resolution (SR), ENO/SR scheme, are extended this way. The other approach is a direct method and a modification of the ENO/SR. Here the technique of ENO reconstruction with subcell resolution is used to locate the discontinuity within a cell and the time evolution is then accomplished by solving the differential equation along characteristics locally and advancing in the characteristic direction. This scheme is denoted ENO/SRCD (subcell resolution - characteristic direction). All the schemes are tested on the equation of LeVeque and Yee (NASA-TM-100075, 1988) modeling reacting flow problems. Numerical results show that these schemes handle this intriguing model problem very well, especially with ENO/SRCD which produces perfect resolution at the discontinuity.
Safety analyses for a nuclear power plants need to consider postulated accidents that results in a risk of accidental release of radiation. The regulations require plant specific safety analysis reports to include an evaluation of the requirements of 10CFR50.67. Such a safety analysis report mandates limits such that calculated radiological consequences at certain locations do not exceed established limits in 10CFR50.67 following a postulated release of radioactivity. Analyses should demonstrate, with reasonable assurance, that these prescribed limits are complied with. Conventional methodologies utilize conservative approaches to address lack of uncertainty quantification in the utilized approaches, methods, and/or inputs. These built-in excess conservatisms often result in with compounding effects and hence overly conservative results in the estimated radiological consequences, which leads to inaccurate margin evaluation for operation and accident mitigation. Therefore, evaluating accurate dose consequences is needed for both operational and safety reasons. This paper illustrates the use of the best estimate source term (BEAST) methodology to the case of the loss of coolant accident (LOCA) dose analysis.
The objectives of this paper are to investigate how various numerical treatments of the nonlinear source term in a model reaction-convection equation can affect the stability of steady-state numerical solutions and to show under what conditions the conventional linearized analysis breaks down. The underlying goal is to provide part of the basic building blocks toward the ultimate goal of constructing suitable numerical schemes for hypersonic reacting flows, combustions and certain turbulence models in compressible Navier-Stokes computations. It can be shown that nonlinear analysis uncovers much of the nonlinear phenomena which linearized analysis is not capable of predicting in a model reaction-convection equation.
Safety analyses for a nuclear power plants need to consider postulated accidents that results in at risk of accidental release of radiation. The regulations require plant specific safety analysis reports to include an evaluation of the requirements of 10CFR50.67. Such a safety analysis report mandates limits such that calculated radiological consequences relative to certain dose locations do not exceed a total effective dose equivalent (TEDE) limits following a postulated release of radioactivity. Calculations are performed to estimate the radiological consequences, in terms of dose, to people and equipment to ensure the estimated doses are within the prescribed limits. Analysis should demonstrate, with reasonable assurance, that these prescribed limits are complied with. Conventional methodologies utilize conservative approaches to address lack of uncertainty quantification in the utilized approaches, methods, and/or inputs. These built-in excess conservatisms often result in compounding effects and hence overly conservative results in the estimated radiological consequences, which leads to inaccurate margin evaluation for operation and accident mitigation. Therefore, evaluating accurate dose consequences is needed for both operational and safety reasons. The research documented in this paper provides an overview of Constellation's Best Estimate Alternate Source Term (BEAST) Methodology. BEAST methodology relies upon the use of realistic yet bounding input distributions for key analysis parameters, which replaces use of conservative deterministic singular inputs that bound overall analysis domain. This approach enables evaluating a more accurate accident analysis response, while enabling a bounding licensing basis envelope via more accurate quantification of uncertainty in the application. Therefore, built-in margin for a given scenario is more accurately evaluated.
There have been several attempts to introduce approximations into the exact form of Lilley's equation in order to express the source term as the sum of a quadrupole whose strength is quadratic in the fluctuating velocities and a dipole whose strength is proportional to the temperature fluctuations. The purpose of this note is to show that it is possible to choose the dependent (i.e., the pressure) variable so that this type of result can be derived directly from the Euler equations without introducing any additional approximations.
The industry-standard MGB approach to predicting the noise generated by a given aerodynamic flow field requires that the turbulence velocity correlation be specified so that the source terms in the Lighthill acoustic analogy may be computed. The velocity correlation traditionally used in MGB Computations is inconsistent with a number of basic qualitative properties of turbulent flows. In the present investigation the effect on noise prediction of using two alternative velocity correlations is examined.
An overview of the hardware and software developed for the Source Term Analysis of Xenon (STAX) project is presented which includes the data collection from two stack monitoring systems installed at medical isotope production facilities, infrastructure to transfer data to a central repository, and methods for sharing data from the repository with users. STAX is an experiment to collect radioxenon emission data from industrial nuclear facilities with the goal of developing a better understanding of the global radioxenon background and the effect industrial radioxenon releases have on nuclear explosion monitoring. The final goal of this work is to utilize collected data along with atmospheric transport modeling to calculate the contribution of a peak or set of peaks detected by the International Monitoring System (IMS) to provide desired discriminating information to the International Data Centre (IDC) and National Data Centers (NDCs). Types of data received from the STAX equipment are shown and collected data was used for a case study to predict radioxenon concentrations at two IMS stations closest to the Institute for RadioElements (IRE) in Belgium. The initial evaluation of results indicate that the data is very valuable to the nuclear explosion monitoring community.