SEARCH · Engineering Papers

Results for “conditionally specified models”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Numerical Study on the Effect of Methane Doping in Hydrogen-Air Rotating Detonation Engines for Various Temperatures and Pressures

Rotating detonation engines (RDEs) have gained attention as a promising technology for future aviation engines. However, the numerical studies of these systems pose severe challenges due to the broad range of spatial and temporal scales. In this study, we use an adaptive mesh refinement based compressible, reactive solver PeleC to resolve the broad range of scales and accurately capture shock and detonation waves using high-resolution numerical schemes. Multi- species transport along with compressible Navier-Stokes equations are solved in the model along with a finite-rate based chemistry model. Embedded boundary method is used to model the complex geometry consisting of discrete fuel nozzles and the combustion chamber. The fuel consisting predominantly of hydrogen is doped with varying levels of methane while air is used as the oxidizer. For a specified total pressure and temperature, the number of stable detonation waves is found to decrease with increasing methane concentration in the fuel mixture. Additionally, no stable detonation solutions are observed for methane composition higher than 20% by volume for the range of operating conditions studied (300-900K, 10 Atm). The increased presence of high temperature zones is also indicative of higher thermal NOx emissions at low methane concentrations. The effect of fuel-air mixture composition and temperature on the detonability, detonation wave structure, mode transitions and their stability are analyzed in this study in addition to their implication on NOx emission.

automatic mesh refinement↗

DOE OSTI2021

Development of Mission Profiles for Humidity Models in the Reliability Testing of PV Inverters: Preprint

To understand the impacts of humidity on PV inverters, mission profiles were developed to accurately describe the different processes and rates based on the environmental factors of temperature, relative humidity, and irradiance. The operating environment of the device introduces stress factors dependent on component temperature, relative humidity, and component voltage. The development of profiles started from looking at historical climate conditions at Cocoa, Florida, a hot and humid location. The profile developed from field location was compared with existing humidity test profiles in industrial and military standards to determine acceleration factors. Prior to the evaluation of the mission profiles, critical components in the PV inverter were identified based on vendor discussions along with failure modes of interest. The profiles considered in this study was demonstrated using a 4-kVA string inverter. The measurements from the devices were substituted in device manufacturer specified equation to determine the severity of each test. The severity of each test was compared against the severity of the 24-hour result to determine the acceleration factor of each test. The acceleration factors were determined to be over a range of 1.2-6 for the existing humidity tests from standards, when compared against the mission profile developed from field data. This optimal humidity/thermal profile will be passed to the IEC standards development for inverter qualification testing.

humidity models↗

DOE OSTI2020

GADRAS Batch Inject Tool User Guide

Gamma Detector Response and Analysis Software (GADRAS) is used by the radiation detection and emergency response community to perform modeling and spectral analysis for gamma detector systems. Built into GADRAS is the ability to define a detector, geometry, background characteristics and source composition to generate synthetic spectra for drills and exercises (injects). Consequence Management is currently in development of a sample result data simulator tool in which a deposition model is probed for source conditions at moments in time and locations in space. These values are used to generate realistic sample results for use in drills and exercises. In addition to sample results, there is a need to simulate the actual spectra that would be observed in the field by downlooking HPGe instruments given a deposition activity. This way, the FRMAC Gamma Spectroscopist can practice their process of generating quantified results from spectra on realistic data as well. Recognizing the decades of work done in GADRAS to accurately generate synthetic spectra, this team decided to build a link between the new simulator and GADRAS to generate these spectra quickly and easily. The simulator tool will generate a file that specifies the name of the spectra, its location, date/time of measurement, duration of measurement, height off the ground, and the deposition activity and age for every radionuclide in the simulation. Then, a new tool within the Inject Tab of GADRAS was developed to read in this file given a detector selection and generate In-Situ spectra for each row in the file in any file format the user chooses. This way, simulation cell staff can take these files and then upload them to the appropriate data system (RAMS or RadResponder) for use during drills and exercises. An advanced feature of this tool allows for generating any spectra given an appropriate model and mapping of source to model layer in the batch inject tool. This way, spectra from field sample counts, mobile laboratories, or even fixed laboratories can be generated in bulk given an estimate of the radioactivity concentration or total radioactivity in an import file. This expands the capabilities of this tool a great deal and will make it a more useful tool for CM and others to help estimate detector response for nearly any situation. This user guide will explain the steps needed to perform a batch inject file generation.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2020

On the Solution of ℓ 0 -Constrained Sparse Inverse Covariance Estimation Problems

The sparse inverse covariance matrix is used to model conditional dependencies between variables in a graphical model to fit a multivariate Gaussian distribution. Estimating the matrix from data are well known to be computationally expensive for large-scale problems. Sparsity is employed to handle noise in the data and to promote interpretability of a learning model. Although the use of a convex ℓ 1 regularizer to encourage sparsity is common practice, the combinatorial ℓ 0 penalty often has more favorable statistical properties. In this paper, we directly constrain sparsity by specifying a maximally allowable number of nonzeros, in other words, by imposing an ℓ 0 constraint. Here, we introduce an efficient approximate Newton algorithm using warm starts for solving the nonconvex ℓ 0 -constrained inverse covariance learning problem. Numerical experiments on standard data sets show that the performance of the proposed algorithm is competitive with state-of-the-art methods.

$\ell_0$-Constrained↗

DOE OSTI2024

Sensitivity of thermodynamic profiles retrieved from ground-based microwave and infrared observations to additional input data from active remote sensing instruments and numerical weather prediction models

Accurate and continuous estimates of the thermodynamic structure of the lower atmosphere are highly beneficial to meteorological process understanding and its applications, such as weather forecasting. In this study, the Tropospheric Remotely Observed Profiling via Optimal Estimation (TROPoe) physical retrieval is used to retrieve temperature and humidity profiles from various combinations of input data collected by passive and active remote sensing instruments, in situ surface platforms, and numerical weather prediction models. Among the employed instruments are microwave radiometers (MWRs), infrared spectrometers (IRSs), radio acoustic sounding systems (RASSs), ceilometers, and surface sensors. TROPoe uses brightness temperatures and/or radiances from MWRs and IRSs, as well as other observational inputs (virtual temperature from the RASS, cloud-base height from the ceilometer, pressure, temperature, and humidity from the surface sensors) in a physical iterative retrieval approach. This starts from a climatologically reasonable profile of temperature and water vapor, with the radiative transfer model iteratively adjusting the assumed temperature and humidity profiles until the derived brightness temperatures and radiances match those observed by the MWR and/or IRS instruments within a specified uncertainty, as well as within the uncertainties of the other observations, if used as input. In this study, due to the uniqueness of the dataset that includes all the above-mentioned sensors, TROPoe is tested with different observational input combinations, some of which also include information higher than 4 km above ground level (a.g.l.) from the operational Rapid Refresh numerical weather prediction model. These temperature and humidity retrievals are assessed against independent collocated radiosonde profiles under non-cloudy conditions to assess the sensitivity of the TROPoe retrievals to different input combinations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Cost Impacts of Risk-Based Methods for Defining AoR and PISC Duration of a CO 2 Storage Project Using NRAP Tools and FE/NETL CO 2 Saline Storage Cost Model

The purpose of this analysis is to assess the fiscal impact of using risk-based methods to define the area of review (AoR) and post-injection site care (PISC) duration for a geologic carbon dioxide (CO 2 ) storage project. The Office of Fossil Energy (FE)/National Energy Technology Laboratory (NETL) CO 2 Saline Storage Cost Model (“cost model”) was utilized to evaluate CO 2 storage economics in this study. The cost model is a widely used tool for evaluating CO 2 storage costs in geologic settings. It estimates costs for storage operations conducted under compliance to the United States Environmental Protection Agency (EPA) Underground Injection Control (UIC) Class VI regulations. This analysis evaluates and compares the costs of storage (on a first-year break-even dollar per tonne [2018$] basis) at the proposed FutureGen 2.0 site as a case study storage location by implementing different approaches to determine AoR and PISC given the prevailing geologic conditions and injection volume and duration considerations. The approaches include the use of 1) risk-based methods to define AoR and PISC duration, 2) an EPA-approved AoR and PISC documented in the FutureGen 2.0 UIC Class VI permit applications, and 3) cost model default settings that utilize uncertainty multipliers to estimate CO 2 plume and pressure front extent as part of AoR determination, as well as a 50-year PISC default. From all three sources, the AoR and PISC specified in the Class VI permit applications for FutureGen 2.0 were the largest and longest relative to the other two sources.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Bayesian Approach to Estimation of Water Table Elevations Using Historical Rasters as Prior Information 2019 - 20430

In cases of complex but only partially known geology and a lack of spatial control in observation well locations, water table elevation estimation is very challenging. In some cases, auxiliary information, such as observations of the movement of tracers, operation of injection and extraction wells, and calibration of groundwater models against the historical elevation data, can be combined with expert judgement to estimate flow directions in areas of sparse data and to aid in the production of more reliable contour maps (and associated rasters) than could be produced by relying on sparse well elevation data alone. Given a historical sequence of these raster maps, the question arises how to automate, to the extent possible, the process of producing new raster maps to reflect data from previous times, the current data and the operation of expert judgement. One solution is to adopt a Bayesian point of view and to regard the historical well elevation data, auxiliary information and historical raster maps as prior information. The well elevations for water table wells, as well as those for injection/extraction wells and the data associated with other relevant variables, can be viewed as predictors for the raster surface. From this prior information, we can, conditional on the values of these predictors for a new time period, compute an expected value map and a standard deviation map for the new raster. These then can be taken to specify a prior predictive distribution for the pixels in the new raster map. Then we condition the pixels, corresponding to water level observation wells within the raster, on the observed values in those wells (which in general will differ from the regression estimate) for the new time period. Given the smoothness of the water table surface, we then smooth the surface of deviations from the mean surface, based on the variograms of the historical rasters, and add this smoothed surface to the regression mean surface. The error structure of the produced raster map is defined by the regression error structure and the error due to smoothing based on the estimated variograms. This methodology has been developed and is being further refined for groundwater monitoring and remediation at LANL. It is a very flexible method that can also be applied with a variety of other predictors applied to model the water level wells in the area of interest over the historical record. The smoothness of the spatial process and its possible evolution over time can then be estimated from the residuals from this regression. This can be augmented by expert hydrogeological opinion based on site topography and hydrogeology. (authors)

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Kinetics of simultaneous hydrodesulfurization and hydrodenitrogenation reactions using CoMoP/Al 2 O 3 and NiMoP/Al 2 O 3

Hydrodesulfurization (HDS) is a key reaction to achieve diesel production at the specified low sulfur levels and is highly affected by a competing reaction involving nitrogen removal through hydrodenitrogenation (HDN). This work evaluated kinetic parameters of simultaneous reactions of HDS of dibenzothiophene (DBT) and HDN of quinoline (Q) using CoMoP/Al 2 O 3 and NiMoP/Al 2 O 3 catalysts under operational conditions that allow a wide range of reagent conversions. Estimated parameters were evaluated using rigorous statistical 2 analysis. Good fits for the evaluated experimental data were provided by both power-law and Langmuir-Hinshelwood models. Further, turnover frequency values highlight adsorption and competition effects between nitrogen-containing compounds and sulfur-containing compounds. NiMoP catalyst showed higher hydrogenating power than CoMoP, with larger absolute value of the estimated adsorption enthalpy (-120 kJ.mol -1 for NiMoP and -75 kJ.mol -1 for CoMoP), suggesting strong adsorption of nitrogen compounds. A catalyst with more hydrogenating power is also more capable of performing both HDN and HDS reactions simultaneously.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Improved Fuel Cycle Capability of Griffin for Fast Reactor Applications

Griffin is a MOOSE based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, the fuel cycle capability has been significantly extended by improving the assembly shuffling option to allow flexible fuel reloading in the multi-cycle depletion calculation and incorporating decay between cycles. An equilibrium core calculation capability was also implemented to find an equilibrium core. Additionally, an enrichment search capability was added to determine the enrichment condition that allows a core to reach an equilibrium cycle with the end-of-cycle k-effective meeting a user-specified target value. The updated fuel management capability has been extensively tested using the three-dimensional ABTR problem with different batch schemes, exhibiting reasonable solutions in terms of manual shuffling, equilibrium cycle, and enrichment search calculations. The cross section generation workflow capability for fast reactors was further verified to produce microscopic cross sections as well as Griffin core inputs for an ring-heterogeneous configuration. Rigorous verification tests using the ABTR problem demonstrated that the RH core calculations, with cross sections and Griffin inputs generated from the cross section workflow of Griffin, produced accurate solution for fast reactor problems. In addition, an option to convert delay neutron parameter data generated from MC 2 -3 in the DLAYXS format into XML format was added to support transient calculations using MC 2 -3-generated data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Influence of initial conditions on data-driven model identification and information entropy for ideal mhd problems

Data-driven methods of model identification are able to discern governing dynamics of a system from data. Such methods are well suited to help us learn about systems with unpredictable evolution or systems with ambiguous governing dynamics given our current understanding. Many plasma problems of interest fall into these categories as there are a wide range of models that exist, however each model is only useful in a certain regime and often limited by computational complexity. To ensure data-driven methods align with theory, they must be consistent and predictable when acting on data whose governing dynamics are known. Weak Sparse Identification of Nonlinear Dynamics (WSINDy) is a recently developed data-driven method that has shown promise in learning governing dynamics from data with high noise levels [1]. This work examines how WSINDy acts on ideal MHD test problems as the initial conditions are varied and specifies limiting requirements for successful equation identification. Furthermore, it is hard to recover the governing dynamics from data that emphasize a single dominant behavior. In these low information cases, Shannon information entropy is able to pick up on the redundancies in the data that affect recoverability.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Development of A Multi-Spectral Pyrometry Sensor for High-Speed Transient Surface-Temperature Measurements in Combustion-Relevant Harsh Environments

Accurate and high-speed transient surface-temperature measurements of combustion devices including internal combustion (IC) engines, gas turbines, etc., provide validation targets and boundary conditions for computational fluid dynamics models, and are broadly relevant to technology advancements such as performance improvement and emissions reduction. Development and demonstration of a multi-infrared-channel pyrometry-based optical instrument for high-speed surface-temperature measurement is described. The measurement principle is based on multi-spectral radiation thermometry (MRT) and uses surface thermal radiation at four discrete spectral regions and a corresponding emissivity model to obtain surface temperature via non-linear least squares (NLLS) optimization. Rules of thumb for specifying the spectral regions and considerations to avoid interference with common combustion products are developed; the impact of these along with linear and non-linear MRT analysis are assessed as a function of temperature and signal-to-noise ratio. A multi-start method to determine the MRT-solution global optimum is described and demonstrated. The resulting multi-channel transient pyrometry instrument is described along with practical considerations including optical-alignment drift, matching intra-channel transient response, and solution-confidence indicators. The instrument demonstrated excellent >97% accuracy and >99% 2-sigma precision over the 400–800 °C range, with ~20 µs (50 kHz, equivalent to 0.2 cad at 2000 RPM IC-engine operation) transient response in the bench validation.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Large eddy simulation of lean blow-off in swirl-stabilized flame with the subgrid dissipation concept

The lean blow-off mechanism of the premixed swirl flame is numerically investigated by large eddy simulation (LES) with the subgrid dissipation concept (SDC) combustion model. Three simulated cases cover stable, near blow-off, and transient conditions. Compared with the experiment, the LES-SDC approach captures the flow and combustion features for stable and near blow-off conditions. More importantly, the predictions of the blowoff procedure and duration agree satisfactorily with the experiment, indicating that the LES-SDC approach is a promising tool for predicting strong, unsteady turbulent combustion processes. Further, the numerical results are used to investigate the blow-off mechanism. Two stages in the blow-off procedure are specified. The first is the necking and extinction of the downstream flame surface, and the second is the shrinking of the upstream flame surface. The blow-off mechanism is well explained by the theory of stretched flame extinction. At the end of the recirculation zone, the large negative radial velocity pushes the flame to the central line. The combustion process here can be abstracted as the stretched counter-flame of the reactant-reactant configuration. The excessive flame stretch dominates the flame extinction and triggers the blow-off event. The upstream flame resists the intense stretch with the help of hot product recirculation, and the combustion here can be idealized as the counter-flame of the reactant-product configuration. The alignment of the temperature gradient and flow velocity, together with the excessive stretch, clearly indicates the tendency of flame local extinction. A Damköhler number-based blow-off criterion is raised from the mean flow strain rate and laminar flame bulk extinction strain rate.

42 ENGINEERING↗

DOE OSTI2023

What is flat ΛCDM, and may we choose it?

The Universe is neither homogeneous nor isotropic, but it is close enough that we can reasonably approximate it as such on suitably large scales. The inflationary-Λ-Cold Dark Matter (ΛCDM) concordance cosmology builds on these assumptions to describe the origin and evolution of fluctuations. With standard assumptions about stress-energy sources, this system is specified by just seven phenomenological parameters, whose precise relations to underlying fundamental theories are complicated and may depend on details of those fields. Nevertheless, it is common practice to set the parameter that characterizes the spatial curvature, Ω K , exactly to zero. This parameter-fixed ΛCDM is awarded distinguished status as separate model, "flat ΛCDM." Ipso facto this places the onus on proponents of "curved ΛCDM" to present sufficient evidence that Ω K ≠ 0, and is needed as a parameter. While certain inflationary model Lagrangians, with certain values of their parameters, and certain initial conditions, will lead to a present-day universe well-described as containing zero curvature, this does not justify distinguishing that subset of Lagrangians, parameters and initial conditions into a separate model. Absent any theoretical arguments, we cannot use observations that suggest small Ω K to enforce Ω K = 0. Our track record in picking inflationary models and their parameters a priori makes such a choice dubious, and concerns about tensions in cosmological parameters and large-angle cosmic-microwave-background anomalies strengthens arguments against this choice. We argue that Ω K must not be set to zero, and that ΛCDM remains a phenomenological model with at least 7 parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

SiC cladding R&D execution plan

To advance the development of accident-tolerant fuel cladding for light-water reactors, the Advanced Fuels Campaign actively addresses the R&D gap between the specialized technological progress in SiC-based composite materials, led by industrial partners, and the fundamental scientific understanding of material behavior under reactor operating conditions. The laboratory activities are systematically organized based on the Development Strategy document, which outlines critical barriers to the use of SiC composites in nuclear energy applications for the coming decade and beyond. To complement the strategy, this execution plan specifies a detailed year-by-year work package and establishes key milestones for assessing progress in technological advancement. This plan has been formulated in collaboration with principal investigators at national laboratories, stakeholders, and industry partners with the goal of closing technology gaps and readying SiC-based cladding technologies for full-scale commercial testing. The research conducted at the national laboratories is categorized into three primary subjects: (1) addressing challenges associated with the performance of SiC cladding under normal operating environments, (2) identifying the failure limits of SiC cladding under accident conditions, and (3) modeling fuel performance for SiC cladding systems. Comprehensive laboratory activities integrate the development of advanced modeling tools for SiC-based cladding systems with experimental validation via integrated testing of multiphysics phenomena. This document lays out the laboratory plan to support the development of SiC-based technologies, which have the potential for significant operational and safety benefits for the U.S. nuclear industry. This plan aligns with industry developer needs, relevant Nuclear Energy University Programs, and Nuclear Energy Advanced Modeling and Simulation Programs to ensure SiC cladding technologies are matured as efficiently as possible and will be maintained and updated to accommodate the dynamically evolving progress of development, capability needs, and resource availability.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Structural Analysis of the US ITER ECH Transmission Line System

The electron cyclotron heating (ECH) and current drive system is one of the main plasma heating systems for ITER. It uses high-power microwave beams with the power deposition location steerable across the plasma cross section. Microwave power is conveyed via transmission lines (TLs) that run from the gyrotrons in the radio frequency building through the assembly hall and tokamak building to the ECH launchers within the tokamak vacuum vessel. The ECH system includes a vast array of interconnected TL waveguides, in-line components, and support structures.Finite element (FE) modeling provides an essential means of simulating the system, applying loads and determining deflections, rotations, forces, moments, and stresses in order to evaluate various structural and microwave transmission performance metrics. A representative FE model of the overall ECH TL system is developed in ANSYS®. This top-level model defines the centerline of the waveguide system. Waveguide segments are represented by line elements (beams and pipes) with equivalent section properties, and support structures are represented by boundary conditions. A systematic approach is used to model each ECH component with lumped masses and structurally equivalent stiffness matrices or ANSYS superelements.The top-level TL FE model is used to evaluate the various loads (thermal, vacuum, seismic, etc.) and operating scenarios. The top-level model directly calculates stresses in the straight aluminum waveguide segments. The model provides the forces and moments acting on the in-line components for detailed submodel assessments. Displacement results from the top-level analysis feed into a separate microwave performance model to help determine operational efficiency. All TL performance and thermal-structural requirements are met, as specified by the applicable codes and standards, and successfully documented in numerous technical reports and demonstrated at the final design review.

Shanmugasundaram, Aravind↗

DOE OSTI2021

Constitutive Model Development for Aging Polymer Encapsulants (ASC P&EM FY2021 L2 Milestone 7836)

This SAND report fulfills the completion requirements for the ASC Physics and Engineering Modeling Level 2 Milestone 7836 during Fiscal Year 2021. The Sandia Simplified potential energy clock (SPEC) non-linear viscoelastic constitutive model was developed to predict a whole host of polymer glass physical behaviors in order to provide a tool to assess the effects of stress on these materials over their lifecycle. Polymer glasses are used extensively in applications such as electronics packaging, where encapsulants and adhesives can be critical to device performance. In this work, the focus is on assessing the performance of the model in predicting material evolution associated with long-term physical aging, an area that the model has not been fully vetted in. These predictions are key to utilizing models to help demonstrate electronics packaging component reliability over decades long service lives, a task that is very costly and time consuming to execute experimentally. The initiating hypothesis for the work was that a model calibration process can be defined that enables confidence in physical aging predictions under ND relevant environments and timescales without sacrificing other predictive capabilities. To test the hypothesis, an extensive suite of calibration and aging data was assembled from a combination of prior work and collaborating projects (Aging and Lifetimes as well as the DoD Joint Munitions Program) for two mission relevant epoxy encapsulants, 828DGEBA/DEA and 828DGEBA/T403. Multiple model calibration processes were developed and evaluated against the entire set of data for each material. A qualitative assessment of each calibration's ability to predict the wide range of aging responses was key to ranking the calibrations against each other. During this evaluation, predictions that were identified as non-physical, i.e., demonstrated something that was qualitatively different than known material behavior, were heavily weighted against the calibration performance. Thus, unphysical predictions for one aspect of aging response could generate a lower overall rating for a calibration process even if that process generated better quantitative predictions for another aspect of aging response. This insurance that all predictions are qualitatively correct is important to the overall aim of utilizing the model to predict residual stress evolution, which will depend on the interplay amongst the different material aging responses. The DSC-focused calibration procedure generated the best all-around aging predictions for both materials, demonstrating material models that can qualitatively predict the whole host of different physical aging responses that have been measured. This step forward in predictive capability comes from an unanticipated source, utilization of calorimetry measurements to specify model parameters. The DSC-focused calibration technique performed better than compression-focused techniques that more heavily weigh measurements more closely related to the structural responses to be predicted. Indeed, the DSC-focused calibration procedure was only possible due to recent incorporation of the enthalpy and heat capacity features into SPEC that was newly verified during this L2 milestone. Fundamentally similar aspects of the two material model calibrations as well as parametric studies to assess sensitives of the aging predictions are discussed within the report. A perspective on the next steps to the overall goal of residual stress evolution predictions under stockpile conditions closes the report.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Multiphysics Pebble-Bed Reactor Control Rod Withdrawal Study

This work studied the responses of both a generic gas- and a fluoride-cooled PBR concept---the gPBR-200 and gFHR, respectively---during reactivity insertion accidents. Both models rely on 2-D axisymmetric simulations to solve the neutron flux distribution, nuclide concentrations, and temperature across the core---in addition to numerous representative pebble and TRISO particle simulations for determining fuel and moderator temperatures. This not only allows for computing maximum temperatures in the core---thus enabling estimation of how near the fuel is to peak operational and safety limits---as prescribed by specified acceptable fuel design limits, which are determined in such a way that fuel is not damaged during operational or anticipated abnormal occurrences---but also predicting how much of the core exceeds a given temperature limit, as well as determining the local energy deposition rate. These models consider both control rod withdrawal and ejection events. The former introduces a great deal more reactivity, as all the control rods are withdrawn (as opposed to a single one in the latter case), though at a much slower pace. In addition, for the gPBR-200, two limiting cases were considered: one with the core starting under hot full-power conditions and one with it starting under cold zero-power conditions. While the amount of reactivity added in the latter case is much higher (due to the far lower temperatures and the lack of neutron poisons such as Xe-135), the margin to temperature limits is also much more significant. Overall, for the design considered, none of the accidents resulted in the maximum fuel temperature reaching values close to the TRISO limit. However, the methodology presented herein could be very relevant if some designs consider reduced margins (e.g., higher temperatures) to achieve enhanced economics. Further model improvement is needed to better capture control rod worth, both in terms of cusping effects (as the rods are slowly withdrawn) and differential worth, especially as the tips of the rods near the upper cavity.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

First-principles free energies by hybrid thermodynamic integration for phase equilibria and fission product solubility in molten salts

The overarching goal of this research project is quantitative microscopic modeling of phase diagrams and solubility of actinides and lanthanide-fission products in molten salts for molten salt reactors. Application of molten salts as coolants and liquid fuels for next-generation reactors requires quantitative understanding of phase diagrams, equations of state and transport properties of binary and ternary molten salt mixtures, as well as solubility of fuel, fission products and reactor-wall corrosion products, but these remain beyond the reach of current first-principles computational techniques. We will address this fundamental capability gap using a novel hybrid framework for thermodynamic integration to directly predict absolute free energies of solid and molten salts at any specified conditions of composition, temperature and pressure. This framework will combine a hierarchy of first-principles, molecular dynamics and continuum techniques to deliver first-principles-level accuracy at a fraction of the computational cost of conventional ab initio molecular dynamics approaches. In particular, this technique is expected to provide first-principles accuracy, while off-loading most of the computational work on much less expensive classical molecular dynamics and continuum techniques.

36 MATERIALS SCIENCE↗