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 English/Portuguese brochure describes the Industrial Technologies Program Save Energy Now model and provides information on tools and resources to help manufacturing facilities reduce industrial energy intensity.
The Multiphysics Applications technical area of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is tasked with assessing, demonstrating, and applying NEAMS tools in solving multiphysics problems of nuclear reactors, such as liquid-fueled molten salt reactors (MSRs), which are the focus of this report. MSRs are actively being pursued as a potential candidate for power and heat generation by the nuclear industry. However, the inherent multiphysics nature of liquid-fueled MSRs, stemming from the strong interrelationship of neutronics, thermal fluids, and chemistry phenomena, provides unique challenges in modeling and simulation (mod/sim). Therefore, it will be important to develop mod/sim tools with varying types of multiphysics coupling that depend on the problem. The current work is focused on the initiation and development of MSR chemistry modeling benchmarks useful for validating current and potential future NEAMS tools. Development of such benchmarks include the following actions: 1) Summarize available chemistry data from the Oak Ridge National Laboratory (ORNL) MSR Program (MSRP) including operation of the Molten Salt Reactor Experiment (MSRE) and design studies for the Molten Salt Breeder Reactor (MSBR) concept; 2) Recommend simulation problems in MSR chemistry mod/sim based on (1); 3) Assess the current state of NEAMS tools that may support (2); 4) Demonstrate and validate the capabilities of the NEAMS tools in (3) while providing iterative feedback on future code development activities. The objective of this report is to initiate this effort by completing actions (1) and (2), which are discussed in Section 2. An overview of the relevant available MSRE experimental data is provided with examples of how this data may be useful in chemistry mod/sim problems, with the caveat that most of this data is over 50 years old, therefore some problem details as well as uncertainty estimates are often not provided. In Section 3, the current state of NEAMS tools is assessed for potential use in these mod/sim problems, with considerations for future code development activities, supporting action (3). Future work supporting this project under NEAMS may include a deeper dive into the specific phenomena discussed here including tasks such as the compilation of additional available data, updates in code development activities, and ultimately the demonstration and validation of these tools, supporting action (4).
We present a computational framework for the shape optimization of a Horizontal-Axis Wind Turbine (HAWT) rotor blade under uncertainty. Our framework integrates aerodynamic simulations based on the blade element method which utilizes reduced order models of the blade structure and wind load with design sensitivity analysis and nonlinear programming. The wind velocity is modeled as a stochastic process to account for variations in time and space. An additional stochastic process accounts for uncertainties in the structural material properties. The uncertainty propagation is based on a non-intrusive polynomial chaos expansion that allows accurate estimation of stochastic performance metrics such as generated power and structural compliance. Sensitivities of cost and constraint functions with respect to shape parameters namely twist angles are computed with an efficient scheme to enable gradient based optimization. To demonstrate the effect of uncertainty, designs obtained from optimization under uncertainty are compared to those obtained from deterministic optimization.
No abstract prepared.
Since the events at the Fukushima-Daiichi nuclear power plant in March 2011 significant research has unfolded at national laboratories, universities and other institutions into alternative materials that have potential enhanced ac- cident tolerance when compared to traditional UO2 fuel zircaloy clad fuel rods. One of the potential replacement claddings are iron-chromium-alunimum (FeCrAl) alloys due to their increased oxidation resistance [1–4] and higher strength [1, 2]. While the oxidation characteristics of FeCrAl are a benefit for accident tolerance, the thermal neu- tron absorption cross section of FeCrAl is about ten times that of Zircaloy. This neutronic penalty necessitates thinner cladding. This allows for slightly larger pellets to give the same cold gap width in the rod. However, the slight increase in pellet diameter is not sufficient to compensate for the neutronic penalty and enriching the fuel beyond the current 5% limit appears to be necessary [5]. Current estimates indicate that this neutronic penalty will impose an increase in fuel cost of 15-35% [1, 2]. In addition to the neutronic disadvantage, it is anticipated that tritium release to the coolant will be larger because the permeability of hydrogen in FeCrAl is about 100 times higher than in Zircaloy [6]. Also, radiation-induced hardening and embrittlement of FeCrAl need to be fully characterized experimentally [7]. Due to the aggressive development schedule for inserting some of the potential materials into lead test assemblies or rods by 2022 [8] multiscale multiphysics modeling approaches have been used to provide insight into these the use of FeCrAl as a cladding material. The purpose of this letter report is to highlight the multiscale modeling effort for iron-chromium-alunimum (FeCrAl) cladding alloys as part of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program through its Accident Tolerant Fuel (ATF) High Impact Problem (HIP). The approach taken throughout the HIP is to utilize lower length scale approaches (e.g., density functional theory, cluster dynamics, rate theory, phase field, and Visco-Plastic- Self-Consistent (VPSC)) to develop more physically informed models at the engineering scale for use in the BISON [9] fuel performance code.
The Single Primary Heat Extraction and Removal Emulator (SPHERE) and Microreactor Applications Research Validation and Evaluation (MARVEL) were utilized to demonstrate thermal hydraulics codes of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. A recent long-duration test of a sodium heat pipe conducted at the SPHERE facility was utilized for the validation of the heat pipe code Sockeye. The first 10 hours of operation featured a startup of the heat pipe from room temperature; this portion was modeled with Sockeye and compared to the experimental data. Measurements included several thermocouple measurements across the length of the heat pipe surface. Two of Sockeye’s heat pipe models were applied to experiment, as well as several models of the cooling at the condenser end of the heat pipe. Sockeye results gave a decent match to the experimental data, with notable discrepancy in the condenser section, where experimental data indicated a significant inactive length of the heat pipe. The MARVEL microreactor was modeled through a coupled simulation of System Analysis Module (SAM) and SubChannel Module (SCM). This multiscale coupled simulation was performed as an exercise to demonstrate the compatibility and applicability of these codes for thermal-hydraulic analysis of microreactors. The coupling between the SCM and SAM codes for multiscale modeling of MARVEL, was achieved through a domain overlapping approach. In this example a transient SAM simulation provided boundary conditions to SCM, and SCM provided pressure drop to SAM.
This English/Russian brochure describes the Industrial Technologies Program Save Energy Now model and provides information on tools and resources to help manufacturing facilities reduce industrial energy intensity.
The U.S. Department of Energy (DOE)’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering-scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. Fuel behavior in a nuclear reactor is governed by a complex network of mechanisms interacting with various other physics aspects in the reactor system. Any model developed to represent the fuel behavior will likely be idealized resulting in uncertainties in their predictions compared to the observed data. As such, this report was motivated by the need to identify the sources of uncertainties and quantify and propagate them through the fuel model outputs. Such quantification of uncertainties will establish a level of model trustworthiness, identify approaches to improve the model trustworthiness, and even guide optimal experiment design for maximal information gain. To accomplish the uncertainty quantification for computational models, this report has relied on the Bayesian framework which provides probabilistic treatment of models their inputs and outputs. The current state-of-the-art on performing Bayesian Uncertainty Quantification (UQ) for nuclear engineering models using High Performance Computing (HPC) resources have been reviewed. Implementation of capabilities for massively parallel Bayesian UQ in Multiphysics Object-Oriented Simulation Environment (MOOSE) is discussed. Several verification cases are discussed to verify the accuracy of the quantified uncertainties using the developed computational capabilities in MOOSE. Then, the problem of quantifying the uncertainties in TRI-Structural isOtropic (TRISO) fuel silver release is addressed. For the first time, the uncertainties arising from the TRISO Fission Gas Release (FGR) model due to model inadequacy and experimental noise are quantified. Also, the Bayesian capabilities are applied to the calibration of the MATPRO creep model, a widely used model in several fuel assessment cases. The impact of the prediction uncertainties in the MATPRO model on the fuel cladding behavior as part of the TRIBULATION assessment case (which is an integral effects case) is investigated. This report concludes with a discussion on the future work for the UQ for computational models.
The U.S. Department of Energy (DOE)’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel-cycle systems. This program has been providing engineering-scale support for the continued development of BISON, a high-fidelity, high-resolution fuel performance tool. Fuel behavior in nuclear reactors is governed by a complex network of mechanisms that interact with various other physics aspects in the reactor system. Any model developed to represent fuel behavior will likely be idealized, resulting in uncertainties when comparing their predictions against the observed data. In Fiscal Year (FY)-23, we initiated the Uncertainty Quantification (UQ) work by using Bayesian methods to establish a level of model trustworthiness and further improve it, with a particular emphasis on TRI-Structural isOtropic (TRISO) nuclear fuel. This year, we further expanded on that UQ work by investigating an approach to quantifying model inadequacy and accounting for lower-length scale (LLS) effects in TRISO silver (Ag) release modeling. Furthermore, we are implementing parallel active learning capabilities to reduce the computational cost (i.e., required computational resources and elapsed time) of performing UQ. Specifically, we utilized The Kennedy O’Hagan framework for Bayesian uncertainty quantification (KOH) to account for model inadequacy in TRISO Ag release predictions made by BISON. The KOH framework represents an improvement over the standard Bayesian framework used in FY-23. Explicitly accounting for model inadequacy in the Bayesian framework helps establish the level of experimental noise uncertainty in the Advanced Gas Reactor (AGR) data. We compared the inverse UQ results obtained from both the standard Bayesian and KOH frameworks in light of the AGR-2/3/4 data, and also compared the predictive UQ results obtained from these two frameworks in light of the AGR-1 data. Next, we investigated the impact of considering LLS effects in the Ag release simulations. We developed an expanded database of LLS simulated effective diffusivities for Ag, covering a wide range of microstructures and temperatures. Using this database, we developed a framework for incorporating LLS effects into the engineering-scale Ag release UQ. We developed both parametric and non-parametric approaches for bridging the length scales. We then investigated the inverse UQ results in light of the AGR-2/3/4 data and the predictive UQ results in light of the AGR-1 data, and compared the LLS-informed approach and the Arrhenius equation, which does not include microstructure information. Finally, we discussed implementing parallel active learning capabilities in the Multiphysics Object Oriented Simulation Environment (MOOSE)/BISON to reduce the computational cost (i.e., computational resources and elapsed time) of Bayesian UQ. For verification purposes, we first tested these new capabil ities on a species interaction problem. We then demonstrated them on the TRISO Ag release application, showing that parallel active learning capabilities can enhance the accuracy of UQ while also substantially reducing the computational cost in comparison to the reference methods developed in FY-23.
NASA’s supply of radioisotopes for Radioisotope Heat Units (RHU) and Radioisotope Thermoelectric Generator (RTG) power sources is facing a crisis due to shortages of Pu-238 for future missions. Am-241 is a possible replacement for Pu-238 since its stockpile from the nuclear weapons program has remained relatively intact. The purpose of this project was to assess the safety of Am-241 and the risks associated with it as compared with PU-238. The project has resulted in a journal article, F. De-La-Torre-Aguilar, N. White, M. Prelas, R. V. Tompson, and S. K. Loyalka "Space nuclear power system accidents: Doses from Pu-238 and Am- 241 inhalation," Progress in Nuclear Energy vol.100C (2017), pp.171-182, which showed that the radiation dosage resulting from Am-241 inhalation are about 1/2.5 of the dosage resulting from Pu-238 inhalation per KW of Pu-238 power. Considerable progress was also made on the Direct Simulation Monte Carlo Method (DSMC) for simulation of aerosol evolution during accidents resulting in three publications: 1.I. Saldivar, F. De-La-Torre-Aguilar, M. Boraas, and S. K. Loyalka “Benchmark problems in aerosol evolution: Comparison of some exact and DSMC results” in Annals of Nuclear Energy vol.117 (2018), pp.213-222. 2. I. Saldivar and S. K. Loyalka, “Evolution of Aerosols Coupled to Environments: Verification of Direct Simulation Monte Carlo (DSMC)”, Nuclear Technology 204 (2018), pp.172-183. 3. M. Boraas and S. K. Loyalka, “Mesh Free Simulation of Spatially Inhomogeneous aerosols in Arbitrary Geometries” Nuclear Science and Engineering (2018). Further, we have computations on Am-241/Pu-238 vapor and aerosol interactions, and found that the vapor condensation rate is relatively insensitive to the mass effects, confirming some previous results. We have also completed computations of charge effects on evolution of different species aerosols and a manuscript is in preparation. Additionally, we have compared atmospheric dispersion and deposition of Pu-238 and Am-241 using the HYSPLIT code. We have also used a real lung model and the CFD program FLUENT to compute PuO2 and AmO2 deposition to explore improvements upon the ICPR-66 model. We have found that the deposition efficiencies of PuO2 and AmO2 are similar. These results (four manuscripts) will be submitted for publication in the post project period. We have made progress using framework of Mars 2020 mission nuclear risk assessment model to integrate all results, but our progress has been limited as the improved consequence modeling itself has been quite involved and challenging. We will integrate and complete the more extensive risk modeling computations and publish the results in the post-project period. In addition to the above, the project has resulted in research training of four PhD students
Tsunami is a series of wave trains which travels with high speed on the sea surface. This traveling wave is caused by the displacement of a large volume of water after the occurrence of an underwater earthquake or volcano eruptions. The speed of tsunami decreases when it reaches the sea shore along with the increase of its amplitudes. Two large tsunamis had occurred in the last decades in Indonesia with huge casualties and large damages. Indonesian Tsunami Early Warning System has been installed along the west coast of Sumatra. This early warning system will give about 10-15 minutes to evacuate people from high risk regions to the safe areas. Here in this paper, a mathematical model for Tsunami evacuation is presented with the city of Padang as a study case. In the model, the safe areas are chosen from the existing and selected high rise buildings, low risk region with relatively high altitude and (proposed to be built) a flyover ring road. Each gathering points are located in the radius of approximately 1 km from the ring road. The model is formulated as an optimization problem with the total normalized evacuation time as the objective function. The constraints consist of maximum allowable evacuation time in each route, maximum capacity of each safe area, and the number of people to be evacuated. The optimization problem is solved numerically using linear programming method with Matlab. Numerical results are shown for various evacuation scenarios for the city of Padang.
Many evaluation, assessment, and modeling tasks at nuclear power plants require spatial in-formation; this often requires physical visits to locations within the facility because the 2D or 3D schematics and current models do not contain enough detail or do not capture as-built and real-world conditions. These visits require extensive manual labor for not only the requesting party, but also support groups, such as security. Light Detection and Ranging (LiDAR) mapping is trying to solve that problem by providing very detailed 3D models for low costs. However, the use of these models can be very limited because either component reference information is missing and too costly to add or there is no way to extract specific spatial data needed for other tools. This report outlines two main efforts. First, to reduce the effort of “Tagging” data in large 3D models, a general Application Programming Interface (API) was developed to import a variety of existing plant database information into a 3D-visualization engine. Filters allow the user to have only zone-specific items listed; then, they can simply click and assign the information to a specific spot or component in the 3D model. The second part of the work is the development of an interface for importing pieces from the3D-LiDAR model into other systems needed for modeling and simulation, outlined around fire modeling. This interface allows for the retrieval of item location and boundaries, enabling the auto generation of models for varying tools.
This report summarizes the mesoscale modeling work performed in fiscal year 2024 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to capture the microstructural evolution and restructuring observed in the dark regions of high burnup UO 2 nuclear fuel. This is the first attempt to realistically simulate the restructuring behavior observed in different region of a high burnup fuel. We employ a grand-potential based phase-field model to concurrently evaluate the formation of subgrains and growth of fission bubbles within the fuel. A energy-based subgrain formation criteria is introduced to simulate the restructuring process. Effect of different initial conditions and different modeling parameters are studies systematically to capture how each of these parameters influence the characteristics of the restructured fuel. It is observed that the subgrain formation begins around existing fission gas bubbles and then proceeds towards triple junctions, grain boundaries and grain interiors. It is demonstrated that restructuring is influenced by a combination of initial dislocation densities, subgrain formation rate, and temperature. Rate of restructuring increases with increase in fuel temperature. A restructuring bias is observed within the microstructure due to variation in defect accumulation among different grains. Furthermore, bubble sizes and distribution does not have a significant effect on rate of restructuring. The predicted microstructures resembles the characteristics of the restructured regions as observed in experiments. Finally, a correlation is presented that demonstrates the evolution of the restructuring volume fraction as a function of local effective burnup. This work provides a first of its kind restructuring model for darkzone that can be used by BISON for performance prediction of high burnup UO 2 fuel.
Pool-type Molten Salt Reactors (MSRs) are a promising advanced nuclear reactor concept relying on passive physical behavior to provide increased safety characteristics. The high heat capacity, high boiling point, unpressurized liquid salt contains the nuclear fuel and as it passes through the vessel achieves criticality and produces nuclear power. The salt then carries away the heat produced along with the delayed neutron precursors towards heat exchangers and primary pumps before entering the core. This forms a highly coupled multiphysics problem between neutronics and fluid flow, which makes modeling these reactors challenging. The safety/licensing case heavily relies on modeling and simulation to prove their safety, and significant effort has been expanded by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to develop the appropriate simulation tools. Their safety/licensing will require a thorough analysis of their behavior during earthquake-based transients. Earthquakes are a transient condition commonly experienced in many regions, and are safely survived by dozens of reactors every year. Nonetheless, significant conservatisms were introduced with limited predictive modeling and simulation available. The first-of-a-kind capability developed by this research will allow for engineers to analyze ranges of design-basis and beyond-design basis accident scenarios to assess the integrity of the reactor vessel and surrounding system, possibly allowing for limiting the costly conservatisms in the design. This projects uses MASTODON, Griffin, and Pronghorn in the NEAMS ecosystem to model the coupled multiphysics problem posed by an earthquake in a molten salt reactor. It marches through a number of number of coupling schemes, from uncoupled simulations to tightly two-way coupled simulations. The appropriate level of coupling for these simulations will then be determined based on the accuracy in the quantities of interest (QoI) and computational effort involved in each scheme. The QoIs chosen for this paper include the mass flow rate across the heat exchanger (for Pronghorn simulations), power output of the reactor (for Griffin simulations), and the maximum Von Mises stress (for MASTODON simulations). This summary paper reports on the current progress of this project, with 2D tightly coupled multiphysics results.
The authors briefly summarize the development of, and test experience with, an initial data interface between the structural mechanics code Diablo and the SHARP reactor simulation system data hub MOAB. That interface has been exercised both to write MOAB databases from Diablo, and then also to use such a database to read in part of a simulation definition for a subsequent Diablo execution. All enhancements are integrated into the central Diablo source repository. The SHARP software system for advanced simulation of nuclear reactors and power plant systems is sponsored by DOE's Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. SHARP has been architected as a federation of single-physics simulation tools to permit flexibility in programming languages and leveraging of past and on-going investments. Solution of multi-physics problems will be coordinated by, and data passed through, a central 'hub'. SHARP's hub implementation is utilizing MOAB: a Mesh-Oriented datABase. This same data hub approach is also intended to enable multi-resolution simulations, e.g, lower-dimension plant-scale simulations can be informed by high-fidelity 3D models of particular critical components.
This report summarizes the lower length scale modeling work performed in the fiscal year 2023 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to capture the microstrutural evolution and associated pulverization criteria in the high burnup regions. This year, the resolution mechanisms within the cluster dynamics code has been updated and it influence on the bubble growth in HBS regions at the mesoscale is studied. To improve the accuracy of phase-field models of fission gas bubble growth, a new Helmholtz free energy for high-density Xe gas is derived from a virial equation of state. Two previously used strategies for representing net vacancy production in phase-field models of fission gas bubble growth are compared with each other and with analytical models of bubble growth. The vacancy source only model is found to be more convenient to parameterize realistically compared with the vacancy source+sink model. Furthermore, the phase-field-fracture simulation have been performed using MD- informed failure stress values and realistic HBS structure obtained from the phase-field simulations. We also present the uncertainty bands on prediction of the critical stress to account for the effect of the lower length scale variabilities on the failure criteria at the mesoscale. It is observed that pulverization may occur in partially restructured regions with restructuring fraction as low as 17%. In light of this, an update to the BISON’s pulverization criteria is recommended.
An effective model for a safety assessment of a conceptual trench repository system, in which low-and intermediate level radioactive waste (LILW) is disposed of has been developed. The repository system is assumed to be located in the surface area of a mountain hill near the seashore. An easy-to-use computer program based on this model, developed as a GoldSim template, has been extended based on the previous one in consideration of the rather possible complex pathways in and around the trench for its performance assessment, which is able to evaluate a nuclide release from the repository and farther transport into the geosphere and biosphere under various normal, disruptive natural, and manmade events, and scenarios that can occur after a failure of waste drums with associated uncertainty. To demonstrate the nuclide release and transport behavior through these various possible pathways in the near- and far-fields of the repository system under certain scenarios, some illustrative evaluations and a simple sensitivity investigation with C-14 for the concrete barriers of the trench were made. (authors)
This report summarizes work performed during Fiscal Year (FY) 2021 at Idaho National Laboratory (INL) for the U.S. Department of Energy’s Nuclear Engineering Advanced Modeling and Simulation (NEAMS) program for the Structural Materials and Chemistry Technical Area in the work package entitled “MS-21IN050106 - Structural Materials - INL.” This effort mainly focused on development and application of capabilities for engineering-scale analysis of nuclear reactor structural components in the Grizzly and BlackBear codes. These efforts include performance improvements and application of models for high-temperature component response, cluster dynamics modeling of precipitation in light water reactor pressure vessel steel, engineering-scale capability development in the areas of concrete simulation and reactor pressure vessel analysis, and preparation for and conducting an independent assessment of the adherence of these codes to software quality standards.