Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model initialization”

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.

At least 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Thermo-Mechanical Phase-Field Modeling of Fracture in High-Burnup UO2 Fuels Under Transient Conditions

This study presents a novel multiphysics phase-field fracture model to analyze high-burnup uranium dioxide (UO2) fuel behavior under transient reactor conditions. Fracture is treated as a stochastic phase transition, which inherently accounts for the random microstructural effects that lead to variations in the value of fracture strength. Moreover, the model takes into consideration the effects of temperature and burnup on thermal conductivity. Therefore, the model is able to predict crack initiation, propagation, and complex morphologies in response to thermal gradients and stress distributions. Several simulations were conducted to investigate the effects of operational and transient conditions on fracture behavior and the resulting cracking patterns. High-burnup fuels exhibit reduced thermal conductivity, elevating temperature gradients and resulting in extensive radial and circumferential cracks. Transient heating rates and temperatures significantly affect fracture patterns, with higher heating rates generating steeper gradients and more irregular crack trajectories. This approach provides critical insights into fuel integrity during accident scenarios and supports the safety evaluation of extended burnup limits.

Chemistry↗
DOE OSTI2024

Prediction of Dielectric Constant in Series of Polymers by Quantitative Structure-Property Relationship (QSPR)

This work is devoted to the investigation of dielectric permittivity which is influenced by electronic, ionic, and dipolar polarization mechanisms, contributing to the material’s capacity to store electrical energy. In this study, an extended dataset of 86 polymers was analyzed, and two quantitative structure–property relationship (QSPR) models were developed to predict dielectric permittivity. From an initial set of 1273 descriptors, the most relevant ones were selected using a genetic algorithm, and machine learning models were built using the Gradient Boosting Regressor (GBR). In contrast to Multiple Linear Regression (MLR)- and Partial Least Squares (PLS)-based models, the gradient boosting models excel in handling nonlinear relationships and multicollinearity, iteratively optimizing decision trees to improve accuracy without overfitting. The developed GBR models showed high R2 coefficients of 0.938 and 0.822, for the training and test sets, respectively. An Accumulated Local Effect (ALE) technique was applied to assess the relationship between the selected descriptors—eight for the GB_A model and six for the GB_B model, and their impact on target property. ALE analysis revealed that descriptors such as TDB09m had a strong positive effect on permittivity, while MLOGP2 showed a negative effect. These results highlight the effectiveness of the GBR approach in predicting the dielectric properties of polymers, offering improved accuracy and interpretability.

Ascencio-Medina, Estefania↗
DOE OSTI2024

"Source Term Modeling for Advanced Gas Micro-Reactors"

Maintaining the safety of the public, environment, and operating personnel is the most important factor in designing, operating, maintaining, and decommissioning nuclear reactors. In recent years, there has been a growing interest in the development of micro-reactors employing TRi-structural ISOtropic (TRISO)-coated particle fuel. In gas reactors, TRISO fuel plays an important role in the safety case for high temperature reactors because of the fission product retention properties of the fuel. This ability enables the use of a functional containment strategy for the reactor where multiple barriers are used to prevent fission product release to the environment. Part of the safety analysis of these advanced reactors is the assessment of radionuclide releases under normal and accident conditions through the multiple credited safety barriers. Using conservative assumptions, a mechanistic analysis can be performed to quantify these releases that combines the probabilistic assessment of failure with analytic solutions to radionuclide transport equations. Source term modeling for TRISO fuel has been performed for previous reactor designs; however, these models are outdated, in many cases proprietary, and need updates to be applied to the current state of TRISO fuel technology and alternative gas reactor core configurations [1]. Currently, the only publicly available source term assessment for gas reactors is an expert-based Monte Carlo simulation based on the effectiveness of the fuel kernel, coating layers, and graphite block in a modular high temperature gas reactor [2]. Thus, there is a need to develop a simple, versatile, and mechanistic model of fission product release and transport in gas reactor cores that could be applied to a variety of reactors through user inputs and reactor-specific radionuclide inventories. The release is calculated by the diffusion of the key safety important fission products through the kernel, silicon carbide (SiC), graphite for both intact and defective TRISO particles based on fuel and graphite temperatures in the reactor under normal operation. These releases from the fuel enter the coolant where they can plate-out on cooler surfaces. A clean-up model is included for designs with a coolant purification system to remove fission gases. This initial distribution of fission products in the reactor serves as an initial condition for potential releases under postulated accident conditions. The model then can calculate the fission product release for any transient temperature profile and fission product releases can then be used to assess radiological dose to the workers and the public using conventional dose tools. Data on the diffusion of fission products is based on historic German TRISO experiments and the more current Department of Energy (DOE) Advanced Gas Reactor (AGR) TRISO fuel development program. The model is coded in python with inputs and outputs in excel spreadsheets, as well as python plotting utilities to aid in the interpretation of the results. References: [1] INL, NGNP Mechanistic Source Term White Paper, INL-10-17997, July 2010. [2] David A. Petti, Richard R. Hobbins, Peter Lowry, Hans Gougar, “Representative Source Terms and The Influence of Reactor Attributes on Functional Containment in Modular High Temperature Gas-cooled Reactors,” Nuclear Technology, Vol. 184, p. 181-197, Nov. 2013.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2019

R-parity violating decays of Bino neutralino LSPs at the LHC

The R-parity violating decays of Bino neutralino LSPs are analyzed within the context of the B - L MSSM “heterotic standard model”. These LSPs correspond to statistically determined initial soft supersymmetry breaking parameters which, when evolved using the renormalization group equations, lead to an effective theory satisfying all phenomenological requirements; including the observed electroweak vector boson masses and the Higgs mass. The explicit RPV decay channels of these LSPs into standard model particles, the analytic and numerical decay rates and the associated branching ratios are presented. The analysis of these quantities breaks into two separate calculations; first, for Bino neutralino LSPs with mass larger than MW± and, second, when the Bino neutralino mass is smaller than the electroweak scale. The RPV decay processes in both of these regions is analyzed in detail. The decay lengths of these RPV interactions are discussed. It is shown that for heavy Bino neutralino LSPs the vast majority of these decays are “prompt”, although a small, but calculable, number correspond to “displaced” decays of various lengths. The situation is reversed for light Bino LSPs, only a small number of which can RPV decay promptly. The relation of these results to the neutrino hierarchy–either normal or inverted–is discussed in detail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows

Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Creation of a prototype biomimetic fish to better understand impact trauma caused by hydropower turbine blade strikes

Biomimetic model organisms could be useful surrogates for live animals in many applications if the models have sufficient biofidelity. One such application is for use in field and laboratory tests of fish mortality associated with passage through hydropower turbines. Laboratory trials suggest that blade strikes are especially injurious and often causes mortality when fish are struck by thinner blades moving at higher velocities. Dose-response relationships have been created from these data, but the exact relationship between fish mortality and the actual forces enacted on fish during simulated blade strike testing remains unknown. Here, we describe the methods used to create a prototype biomimetic model fish composed of ballistic gelatin and covered with a surrogate skin to better approximate the biomechanical properties of a fish body. Frozen fish were scanned with high-fidelity laser scanners, and a 3D-printed, reusable mold was created from which to cast our gelatin model. Computed tomography scan data, imaged directly or taken from online data repositories, were also successfully used to create CAD models for use in additive manufacturing of molds. One 3-axis accelerometer was embedded into the gelatin to compare accelerometer data to dose-response data from previous laboratory research on live fish. The resulting model ( i.e. , Gelfish) had a statistically indistinguishable tissue durometer to that of real fish tissue and preliminary blade strike impact testing suggested its overall flexibility was similar to that of live fish. Gelfish was designed with biofidelity as its guiding principle and our results suggest initial experimentation was successful. Future research will include replication of initial Gelfish test results, quantitative measurement of model flexibility relative to real fish, and inclusion of surrogate skeletal structures to enhance biofidelity. Use of more sophisticated sensors would also better quantify the physical forces of blade strike impact and help determine how said forces correlate with rates of mortality observed during tests on live fish.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Proton Spin from Small-x With Constraints from the Valence Quark Model

We apply the valence quark model to constrain the non-perturbative initial condition for the small-x helicity evolution. The remaining free parameters are constrained by performing a global analysis to the available polarized small-x deep inelastic scattering data. A good description of the world data is obtained with only 8 free parameters. The model parameters are tightly constrained by the data, allowing us to predict the proton polarized structure-function to be negative at small x. Furthermore, we obtain the small-x quark and gluon spins, depending on the applied running coupling prescription.

Adamiak, Daniel [Thomas Jefferson National Acceler↗
DOE OSTI2022

Initializing BSQ with Open-Source ICCING

While it is well known that there is a significant amount of conserved charges in the initial state of nuclear collisions, the production of these due to gluon splitting has yet to be thoroughly investigated. The ICCING (Initial Conserved Charges in Nuclear Geometry) algorithm reconstructs these quark distributions, providing conserved strange, baryon, and electric charges, by sampling a given model for the g → qq¯ splitting function over the initial energy density, which is valid at top collider energies, even when µB = 0. The ICCING algorithm includes fluctuations in the gluon longitudinal momenta, a structure that supports the implementation of dynamical processes, and the c++ version is now open-source. A full analysis of parameter choices on the model has been done to quantify the effect these have on the underlying physics. We find there is a sustained difference across the different charges that indicates sensitivity to hot spot geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Determining the extent of potential fugitive fluid migration from geologic carbon storage in hydrocarbon-bearing reservoirs: Insights from one-dimensional numerical modeling

Numerical modeling of Geologic Carbon Sequestration in permeable reservoirs initially containing hydrocarbons is conducted using the multi-phase, multi-component thermohydrologic simulator TOGA (TOUGH Oil, Gas, Aqueous; TOUGH stands for Transport Of Unsaturated Groundwater and Heat), to determine how phase and composition of the original fluids influence the extent of the zone where upward fugitive fluid migration could potentially occur, denoted R f . The area within R f comprises regions of substantially elevated pressure and free-phase CO 2 saturation, where a breach in reservoir sealing capacity would lead to upward fugitive fluid migration. The model examines the conditions within the storage reservoir that could lead to fugitive flow, but does not model the fugitive flow itself. A one-dimensional radial model of the storage reservoir is used, and three initial phase conditions are considered: single-phase aqueous, two-phase gas-aqueous, and three-phase oil-gas-aqueous. Components that may be present are H 2 O, CO 2 , CH 4 , C 4 H 10 , and C 10 H 22 . The most important factors controlling Rf are (1) the initial gas-phase saturation within the reservoir, and (2) the lateral extent of multi-phase initial conditions, particularly CO 2 . The composition of liquid and gas phases has a secondary effect. The impact of reservoir depth, thickness, injection rate, and hydrologic properties are also briefly examined, with thickness (or equivalently injection rate) having the biggest effect. These results can help to understand important trends in potential response of CO 2 -EOR fields being considered for dedicated CO 2 storage.

CO₂ plume migration↗
DOE OSTI2019

Simulating the Fracture of Notched Mortar Beams through Extended Finite-Element Method and Peridynamics

This paper simulates fracture in notched mortar beams under three-point bending using extended finite element method (XFEM) and peridynamics. A three-phase microstructure (i.e., cement paste, aggregates, and paste-aggregate interface) is used for constitutive modeling of the mortar to obtain the elastic properties for simulation. In the XFEM approach, the simulated homogenized elastic modulus is used along with the total fracture energy of the cement mortar in a damage model to predict the fracture response of the mortar including crack propagation and its fracture parameters (Mode I stress intensity factor, K IC and critical crack tip opening displacement, CTOD C ). The damage model incorporates a maximum principal stress-based damage initiation criteria and a traction-separation law for damage evolution. In the peridynamics approach, a bond-based model involving a prototype microelastic brittle (PMB) material model is used. The elastic properties and fracture energy release rates are used as inputs in the PMB model, along with the choice of peridynamic horizon size. Comparison with experimental fracture properties ( K IC , CTOD C ) as well as crack propagation paths from digital image correlation show that both the approaches yield satisfactory results, particularly for K IC and crack extension. Furthermore, both these methods can be adopted for fracture simulation of cement-based materials.

42 ENGINEERING↗
DOE OSTI2019

A Simple Hydride Model for Cerium Ejecta Particles

Cerium ejecta particles form after a shock wave impacts a cerium plate with a known roughness or prescribed surface perturbation. When these particles are ejected into hydrogen or deuterium gas the material reacts exothermically and raises the temperature of the particle relative to the gas. A simple model is developed, assuming that the particle remains intact and spherical, to capture the heat and mass transfer that occurs during this process. Model performance is evaluated by comparing with cerium ejecta experiments where the particles are ejected into deuterium gas at initial pressures of 4 and 8 atm. Overall, the model is able to capture the approximately 400 K increase in particle temperature above the surrounding gas temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A simple hydride model for cerium ejecta particles

Cerium ejecta particles form after a shock wave impacts a cerium plate with a known roughness or prescribed surface perturbation. When these particles are ejected into hydrogen or deuterium gas the material reacts exothermically and raises the temperature of the particle relative to the gas. A simple model is developed, assuming that the particle remains intact and spherical, to capture the heat and mass transfer that occurs during this process. Model performance is evaluated by comparing with cerium ejecta experiments where the particles are ejected into deuterium gas at initial pressures of 4 and 8 atm. Overall, the model is able to capture the approximately 400 K increase in particle temperature above the surrounding gas temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A simple hydride model for cerium ejecta particles

Cerium ejecta particles form after a shock wave impacts a cerium plate with a known roughness or prescribed surface perturbation. When these particles are ejected into hydrogen or deuterium gas the material reacts exothermically and raises the temperature of the particle relative to the gas. A simple model is developed, assuming that the particle remains intact and spherical, to capture the heat and mass transfer that occurs during this process. Model performance is evaluated by comparing with cerium ejecta experiments where the particles are ejected into deuterium gas at initial pressures of 4 and 8 atm. Overall, the model is able to capture the approximately 400 K increase in particle temperature above the surrounding gas temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

SAS4A/SASSYS-1 Validation with the FFTF LOFWOS Tests #10-12

The FFTF liquid sodium cooled test reactor was built to assist the development and testing of advanced fuels and material for fast reactors. The goals of the FFTF Passive Safety Testing program included confirming the safety margins of FFTF as a liquid metal reactor, providing data for computer code validation, and demonstrating the inherent and passive safety benefits of its specific design features. The program included thirteen unprotected loss of flow without scram tests. The results of modeling and simulation of LOFWOS Tests #10, #11, and #12 are presented in this report to support validation of Argonne’s SAS4A/SASSYS-1 fast reactor safety analysis code. These three tests did not use the primary loop pony motors and thus transitioned to natural circulation flow rates from varying initial power levels. This validation activity leverages modeling and simulation efforts that originated under an International Atomic Energy Agency Coordinated Research Project (CRP) for benchmark analysis of the FFTF LOFWOS Test #13. The SAS4A/SASSYS-1 model originally developed for the CRP has been extended to include LOFWOS Tests #10-12 and modernized to utilize new code features and modeling practices. Predicted results from SAS4A/SASSYS-1 were compared to measured test data and differences were investigated. Overall, it was concluded that, for LOFWOS Tests #10-12, there was good agreement between the flow, power, and temperature predictions and the measured data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
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 OSTI2018

Dual-Purpose Canister Filling Demonstration Project Progress Report, Rev. 1

This report discusses the initial progress made at the Oak Ridge National Laboratory to support direct disposal of dual-purpose canisters (DPCs)using filler materials to demonstrate that the probability of criticality in DPCs during disposal to be below the probability for inclusion in a repository performance assessment. In the initial phase of a multi-phase effort that will result in a full-scale demonstration, a computational fluid dynamics (CFD) model was developed to gauge the filling process and to uncover any unforeseen issues. The initial filling simulations of the lower region (mouse holes) of a prototypic DPC show successful removal of the inner space voids and smooth,even progression of the liquid level. In the initial phase, flow through a pipe that is similar to the drain pipe in a DPC will be investigated separately to gain valuable insight of flow regime inside a pipe. The initial experimental setups for validating the computational filling model have been designed,and the various assembly parts are being procured. The experience gained from the initial experiments will be applied to the next steps toward a full-scale demonstration and to the validation of multiphysics filling simulation models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2019

Dual-Purpose Canister Filling Demonstration Project Progress Report

This report discusses the initial progress made at the Oak Ridge National Laboratory to support direct disposal of dual-purpose canisters (DPCs) using filler materials to demonstrate that the probability of criticality in DPCs during disposal to be below the probability for inclusion in a repository performance assessment. In the initial phase of a multi-phase effort that will result in a full-scale demonstration, a computational fluid dynamics (CFD) model was developed to gauge the filling process and to uncover any unforeseen issues. The initial filling simulations of the lower region (mouse holes) of a prototypic DPC show successful removal of the inner space voids and smooth, even progression of the liquid level. In the initial phase, flow through a pipe that is similar to the drain pipe in a DPC will be investigated separately to gain valuable insight of flow regime inside a pipe. The initial experimental setups for validating the computational filling model have been designed, and the various assembly parts are being procured. The experience gained from the initial experiments will be applied to the next steps toward a full-scale demonstration and to the validation of multiphysics filling simulation models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Dual-Purpose Canister Filling Demonstration Project Progress Report

This report discusses the initial progress made at the Oak Ridge National Laboratory to support direct disposal of dual-purpose canisters (DPCs) using filler materials to demonstrate that the probability of criticality in DPCs during disposal to be below the probability for inclusion in a repository performance assessment. In the initial phase of a multi-phase effort that will result in a full-scale demonstration, a computational fluid dynamics (CFD) model was developed to gauge the filling process and to uncover any unforeseen issues. The initial filling simulations of the lower region (mouse holes) of a prototypic DPC show successful removal of the inner space voids and smooth, even progression of the liquid level. In the initial phase, flow through a pipe that is similar to the drain pipe in a DPC will be investigated separately to gain valuable insight of flow regime inside a pipe. The initial experimental setups for validating the computational filling model have been designed, and the various assembly parts are being procured. The experience gained from the initial experiments will be applied to the next steps toward a full-scale demonstration and to the validation of multiphysics filling simulation models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗