Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational materials science”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Impact of the plastic deformation microstructure in metals on the kinetics of recrystallization: A phase-field study

The sensitivity of recrystallization kinetics in metals to the heterogeneity of microstructure and deformation history is a widely accepted experimental fact. However, most of the available recrystallization models employ either a mean field approach or use grain-averaged parameters, and thus neglecting the mesoscopic heterogeneity induced by prior deformation. Here in the present study, we investigate the impact of deformation-induced dislocation (subgrain) structure on the kinetics of recrystallization in metals using the phase-field approach. The primary focus here is upon the role of dislocation cell boundaries. The free energy formulation of the phase-field model accounts for the heterogeneity of the microstructure by assigning localized energy to the resulting dislocation microstructure realizations generated from experimental data. These microstructure realizations are created using the universal scaling laws for the spacing and the misorientation angles of both the geometrically necessary and incidental dislocation boundaries. The resulting free energy is used into an Allen-Cahn based model of recrystallization kinetics, which are solved using the finite element method. The solutions thus obtained shed light on the critical role of the spatial heterogeneity of deformation in the non-smooth growth of recrystallization nuclei and on the final grain structure. The results showed that, in agreement with experiment, the morphology of recrystallization front exhibits protrusions and retrusions. By resolving the subgrain structure, the presented algorithm paves the way for developing predictive kinetic models that fully account for the deformed state of recrystallizing metals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Scaling microstructural processes in the sintering of ionic ceramics

A multi-scale framework, combining a multiphase field formulation and large deformation mechanics, was developed as a stepping stone to perform the data analytics of the microstructural level kinetics of a sintering solid. Relevant microstructural information from this framework, such as grain, stress, and porosity statistics, was scaled up to describe the macroscopic level sintering kinetics. Here, the developed formulation was applied to describe the electric field assisted sintering of Y 2 O 3 . Microstructural inhomogeneities in a multi-granular solid result in the formation of a field of compressive stress networks, which interleave with low compression and weakly tensile regions, defining a scaffolding for sintering concentration regions to develop. A Poisson effect-induced lateral stress network is also naturally self-induced as a result of the mechanical constraints imposed by the sintering apparatus. For long sintering times, localized shear stresses enhancing mass flow along grain boundaries and internal surfaces develop. Three-sided pores are removed by either vacancy transport to the surrounding pores, or move towards the external surfaces through grain boundary diffusion. Four- and higher order-sided pores stabilize because an equal amount of vacancies are gained and lost through the connecting grain boundaries. Grain dewetting contributes to pore coalescence, suggesting that pore kinetics and grain growth are coupled and should be analyzed in concert. The combined sintering and grain growth kinetics define six regimes of sintering behavior: (1) T, the transient regime; (2) E$_Υ$, the surface energy dominated, early sintering regime, where the grain growth exponent, p = 1, and the stress concentration factor, $f$ ~ $1/\hat{ρ}^{4.6}$; (3) E S , the stress dominated, early sintering regime, where p = 1 and $f$ ~ $1/\hat{ρ}^{4}$; (4) I$_Υ$, the surface energy dominated, intermediate sintering regime, where p = 2 and $f$ ~ $1/\hat{ρ}^{4.6}$; (5) I S , the stress dominated, intermediate sintering regime, where p = 2 and $f$ ~ $1/\hat{ρ}^{4}$; and (6) L, the late sintering regime, where p = 3 and $f$ ~ 1. At the macroscopic level, the rapid densification and suppression of grain growth observed in the electric field assisted sintering process is a consequence of the compounding effects of the underlying stress-, transport-, and interfacial-energy-induced energy minimization kinetics, as predicted by the multi-scale framework.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Modeling brittle fracture due to anisotropic thermal expansion in polycrystalline materials

Here, this work investigated brittle fracture of polycrystalline materials due to thermal stresses arising from anisotropic thermal expansion. We used phase-field fracture simulations with the properties of alpha-uranium (α-U) and assumed a linear elastic mechanical response. Three-dimensional simulations were used to predict fracture for various conditions and crystallographic textures. We found that fracture was more pronounced during cooling than during heating because the anisotropy increased with temperature. We also found that the total crack surface area increased with increasing average misorientation, while the net shape change of the material decreased with increasing misorientation. Two-dimensional simulations in which one crystallographic coefficient of thermal expansion (CTE) was set to zero indicated that the expansion behavior in the crystallographic direction with the smallest CTE was the primary cause of fracture.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

The impact of anisotropic thermal expansion on the isothermal annealing of polycrystalline α -uranium

Although grain growth impacts microstructural evolution in a wide variety of materials systems, the effect of anisotropic thermal expansion on grain boundary mobility and texture evolution has not been widely studied. Anisotropic thermal expansion occurs in multiple non-cubic metals, and the thermomechanical processing behavior of these materials can be better understood with further study into the impact of thermal expansion on grain boundary mobility and texture evolution. In this work, we develop a mesoscale phase field model of grain growth that includes the effect of anisotropic thermal expansion, which is applied to study polycrystalline α-uranium, a highly anisotropic metal. Three-dimensional simulations on polycrystalline a-uranium with and without thermal expansion eigenstrains are performed to study the grain boundary mobility and texture evolution as a function of temperature. A strain-free temperature of 933 K is selected, and the system is studied at within the range of 873 – 933 K at intervals of ten degrees, resulting in increasing thermal eigenstrain with decreasing temperature. We also estimate a grain boundary mobility prefactor and activation energy based on existing experimental data of isothermal annealing of α-uranium. The grain boundary mobility is found to display significant deviation from Arrhenius behavior with the inclusion of thermal expansion eigenstrain as the amount of thermal eigenstrain (and thus elastic strain energy within the system) increases. This result explains an experimentally observed grain boundary mobility deviation from Arrhenius behavior. Furthermore, the texture evolution is affected, such that the grain orientations become less random with increasing thermal eigenstrain, which could explain experimentally observed texture behavior. These results indicate that the effect of thermal expansion should be considered when predicting the thermomechanical processing behavior of α-uranium and other materials with anisotropic thermal expansion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Impact of amorphous pockets on displacement damage evolution in silicon

Silicon has long been known to exhibit amorphization in response to heavy particle bombardment. For doses below the total amorphization threshold, partial amorphization is observed in the form of scattered amorphous pockets. While extensive research has gone into modeling the formation and evolution of amorphous pockets in response to irradiation, no studies yet investigate their impact on the evolution of other damage such as interstitial supersaturation and clustering. In this study, we survey the impact of amorphous pockets on defect evolution in silicon when treated as static sinks. MD is first used to show that amorphous pockets provide energetically favorable sites for point defects relative to the crystalline bulk, supporting the hypothesis that they act as sinks. A 0-D cluster dynamics model is then constructed, taking an interstitial clustering model from the literature and including amorphous pockets as a sink species. We conduct our survey for temperatures between 30 and 400 °C and sink strengths between 1 to 6 x 10 10 cm −2 . Both implantation- and radiation-induced damage states are investigated using interstitial and vacancy concentrations as initial condition variables. We find that, due to the differing migration rates of the interstitial and the vacancy, amorphous pockets have a non-monotonic impact on the final damage state depending on the effective sink strength of the amorphous pockets, resulting in increased damage formation in regimes of intermediate amorphization. In conclusion, this result emphasizes the important role of amorphous pockets in governing the evolution of damage in partially amorphized crystalline materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Bayesian discovery of optimal reduced order models from mechanistic and experimental data: A case study of Pd penetration in TRISO fuels using BISON

TRistructural ISOtropic (TRISO) particles rely on a silicon carbide (SiC) layer as the primary structural material and barrier to metallic fission products (FPs) release. Accurate prediction of palladium (Pd) transport and penetration is therefore critical for qualifying TRISO fuels for advanced reactors. The empirical correlation for Pd penetration in BISON is derived from historical particle-fuel data, but cannot explain the large scatter in the experimental data that arises from varying experimental conditions. To aid fuel qualification, we previously developed a mechanistic reduced order model (ROM) using BISON that resolves these dependencies. Here, in this work we build on that mechanistic ROM and perform validation and quantify its uncertainty using Bayesian uncertainty quantification (UQ). calibration against a suite of in-pile and out-of-pile experiments spanning particle compositions, geometries, and operating conditions, and we benchmark it against the empirical correlation. Bayesian UQ identifies influential parameters, calibrates them to data, and yields predictive intervals. Results show that while the empirical correlation can be tuned to fit a single experiment type, it transfers poorly; the mechanistic ROM sustains accuracy with credible uncertainty across disparate conditions. This demonstrates a practical path—via Bayesian UQ applied to mechanistic ROMs—to leverage single-effect experiments for inferring in-reactor behavior and supporting TRISO fuel qualification.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Incorporating corrosion design constraints in desalination process optimization: A case study in mechanical vapor compression

Corrosion is an expensive and complex challenge for desalination, yet current design approaches do not explicitly account for corrosion mechanisms in process modeling and technoeconomic analysis. Here, to address this gap, we present a workflow for incorporating corrosion design constraints directly into desalination process optimization models. We develop surrogate models for general and localized corrosion metrics as functions of temperature, pH, salinity, dissolved oxygen, and material using data from OLI Systems’ Corrosion Analyzer. We then integrate these surrogates as corrosion design constraints in a cost-optimization MVC model that minimizes the levelized cost of water (LCOW). For a case study of mechanical vapor compression (MVC) treating seawater across a range of recoveries, we find dissolved oxygen (DO) is the dominant driver of localized corrosion, and thus of cost-optimal material choice and operating conditions. Reducing the DO from 8 mg/L to 0.5 mg/L reduces the LCOW by 15-35%, informing the breakeven costs for implementing DO removal or selecting highly corrosion-resistant alloys. This framework is broadly applicable across corrosion types, materials, and components and enables desalination process design that minimizes capital costs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Automatic information extraction from neutron radiography imaging to estimate axial fuel expansion in $\mathrm{EBR-II}$

Approximately 130,000 metal fuel pins were irradiated in the Experimental Breeder Reactor II (EBR-II) during its 30 years of operation to develop and characterize existing and prospective fuels. For many of the metal fuel irradiation experiments, neutron radiography imaging was performed to characterize fuel behavior, like fuel swelling. However, due to the lack of technology or resources, many of the images have not been processed or were processed manually through visual examination. This paper represents first-attempt to develop an image processing algorithm capable of automatically extracting information regarding the degree of fuel swelling from neutron radiography imaging. The algorithm was applied to 120 images of three different metallic fuel pin compositions—U-10Zr, U-8Pu-10Zr, and U-19Pu-10Zr. The algorithm performs operations of image intensity adjustment, image binarization, region finding, and labeling to extract information about fuel swelling. The average growth for U-10Zr was found to be 8.49% with 95% Confidence Interval (CI) [8.33 –8.66%], for U-8Pu-10Zr — 7.50% with 95% CI [7.23 – 7.78%], and for U-19Pu-10Zr — 3.15%, with 95% CI [2.40 – 3.91%]. The results obtained by applying this automatic image processing algorithm are consistent with previously reported studies of the same types of fuels. The automatic image processing algorithm will be expanded to include thousands of available neutron radiography images and different types of fuels to investigate empirical dependencies of fuel swelling, which can be subsequently applied to advanced fuel modeling. Results from this study can later be compared to BISON simulations to further benchmark modeling efforts and develop assessment cases.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

A mechanistic, multiscale model for predicting Pd penetration in TRISO fuels using BISON

TRistructural ISOtropic (TRISO) particles use silicon carbide (SiC) as the primary structural member and barrier against metallic fission product (FP) release. palladiums (PDs), produced by fission in the fuel kernel, can diffuse to and chemically interact with the SiC layer, degrading its structural integrity and ability to contain radioactive FPs. Existing temperature-dependent correlations for predicting Pd penetration rely on experimental data with significant scatter due to varying conditions, potentially complicating ongoing fuel qualification and licensing efforts for advanced reactors that would subject TRISO fuels to operating conditions outside of those examined in the experiments. A mechanistic model of Pd production, transport, and reaction is developed in this work to better understand and predict PD attack of SiC in TRISO particles. molecular dynamicss (MDs) simulations are utilized to calculate the Pd diffusivity in SiC grain bulk and grain boundaries. A mesoscale phase-field diffusion model, informed by the MD diffusivities, is used to develop a reduced order model (ROM) for the effect of SiC microstructure on the Pd penetration rate. The engineering scale BISON model calculates Pd production and transport, utilizing the ROM to predict penetration rates consistent with experimental data. This novel mechanistic ROM captures the effect of temperature, microstructure, and irradiation history on Pd penetration. In conclusion, these new capabilities are expected to support ongoing qualification and licensing efforts associated with near-term TRISO-fueled reactor applications.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Multiphysics simulation of TRISO fuel compacts and the effects of homogenization on silver release predictions

This work studies the impact of explicit and homogenized modeling approaches on the multiphysics simulation of TRistructural ISOtropic (TRISO) fuel compacts in prismatic High Temperature Gas Reactors (HTGRs) and the silver release predictions. TRISO fuel compacts exhibit complex double heterogeneity that significantly affects heat conduction, neutron transport, and fission products release. In this work, we use Cardinal, a multiphysics tool based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, to couple neutron transport and heat conduction. OpenMC is used for neutron transport and the MOOSE heat transfer module is used for the heat conduction. Then, we use BISON for the silver release predictions. Two different modeling approaches—explicit modeling of individual TRISO particles and homogenized representation using effective thermal properties—are compared at high TRISO packing fractions (20% and 40%) across varying power densities. Results demonstrate that homogenization significantly under-predicts peak temperatures, for the case of high power/TRISO, there is a difference of 83.73 K in the maximum temperature. Additionally, homogenization underestimates the silver release predictions compared to explicit modeling. Specifically, the temperature and power density differences lead to significant differences in silver release predictions. This work demonstrates the importance of accurately modeling heterogeneity of the TRISO particles to reliably predict fission product release and assess reactor safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Promoting regulatory acceptance of combined ion and neutron irradiation testing of nuclear reactor materials: Modeling and software considerations

As the needs for the nuclear energy industry continue to evolve in the 21st century, timely adoption of new technological solutions acceptable to regulatory agencies is critical. Quantitative prediction of radiation damage in materials and its impact on mechanical properties is a key component of licensing and regulatory decisions regarding nuclear power plants. Accelerated testing methodologies such as combined ion and neutron irradiation data sets are crucial for the development and deployment of new materials and new manufacturing methods (e.g., additive manufacturing). However, regulatory acceptance of accelerated testing methodologies is necessary for their adoption. Further, the present work discusses the fundamental basis for comparing ion- and neutron-induced material microstructures, the theory behind interpreting radiation damage across length and time scales and radiation types, and the codes, standards, and quality assurance concerns surrounding different modeling methods and software. In particular, recommendations are given as to the path forward that will enable national laboratories, academia, and industry to develop the modeling and software basis for regulatory acceptance of the combined use of ion and neutron irradiation for material performance evaluation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Anisotropic grain boundary area and energy distributions in tungsten

Describing microstructure evolution in tungsten requires a quantitative description of the anisotropic grain boundary energy. Here we present a grain boundary energy function for tungsten that specifies the energy of an arbitrary boundary given its five macroscopic crystallographic parameters. A comparison of measured grain boundary areas and the grain boundary energies given by the function at the Σ11, Σ17b, and Σ33a misorientations, which are problematic to determine by measurement or atomistic calculations, reveals inverse correlations that are similar to what have been observed in other metals.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Economic NMPC for a Reversible Solid Oxide Cell

Reversible solid oxide fuel cells (rSOCs) offer the flexibility to operate in tandem with the electric grid by switching between fuel cell and electrolysis modes based on real-time electricity prices. However, their complex, tightly coupled dynamic behavior poses significant challenges in determining optimal operating strategies. In this work, we present an economic nonlinear model predictive control (E-NMPC) framework to optimize the operation of rSOCs. The proposed E-NMPC is applied to a detailed rSOC flowsheet model that includes a utility scale rSOC module as well as balance-of-plant equipment necessary for thermal management. Our results demonstrate that in fuel cell mode, the E-NMPC strategy reduces hydrogen consumption compared to conventional set-point tracking NMPC, while maintaining the same level of electricity output. Also, in electrolysis mode, the E-NMPC yields a marginal improvement in hydrogen production. In addition, we explore the integration of a battery with the rSOC system to enhance flexibility in meeting electricity production and consumption targets.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗