Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomistic 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.

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

At least 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A critical review of irradiation-induced changes in reactor pressure vessel steels

A large body of work has been conducted to investigate the embrittlement and degradation of reactor pressure vessel (RPV) steels. This includes experiments on alloys with different compositions, performed in research and test reactors and ion accelerators that span various temperatures, fluxes, and fluences. In this paper, we perform a critical review of the published experimental data and compile experimentally reported values for dislocation loop size/density, precipitate size/density and yield stress into an easily downloadable format that can be used by both experimentalists and modelers. This thorough experimental review is complemented by a brief review of simulation efforts at atomistic and mesoscopic length scales. Finally, this paper highlights key aspects of the behavior of RPV steels under irradiation, identifies gaps or discrepancies in current understanding, and identifies priority future research directions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Identifying chemically similar multiphase nanoprecipitates in compositionally complex non-equilibrium oxides via machine learning

Abstract Characterizing oxide nuclear fuels is difficult due to complex fission products, which result from time-evolving system chemistry and extreme operating environments. Here, we report a machine learning-enhanced approach that accelerates the characterization of spent nuclear fuels and improves the accuracy of identifying nanophase fission products and bubbles. We apply this approach to commercial, high-burnup, irradiated light-water reactor fuels, demonstrating relationships between fission product precipitates and gases. We also gain understanding of the fission versus decay pathways of precipitates across the radius of a fuel pellet. An algorithm is provided for quantifying the chemical segregation of the fission products with respect to the high-burnup structure, which enhances our ability to process large amounts of microscopy data, including approaching the atomistic-scale. This may provide a faster route for achieving physics-based fuel performance modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Theory of Defect-Induced Crystal Field Perturbations in Rare-Earth Magnets

We present a theory describing the single-ion anisotropy of rare-earth (RE) magnets in the presence of point defects. Taking the RE-lean 1 ∶ 12 magnet class as a prototype, we use first-principles calculations to show how the introduction of Ti substitutions into SmFe 12 perturbs the crystal field, generating new coefficients due to the lower symmetry of the RE environment. We then demonstrate that these perturbations can be described extremely efficiently using a screened point charge model. We provide analytical expressions for the anisotropy energy that can be straightforwardly implemented in atomistic spin dynamics simulations, meaning that such simulations can be carried out for an arbitrary arrangement of point defects. The significant crystal field perturbations calculated here demonstrate that a sample that is single phase from a structural point of view can nonetheless have a dramatically varying anisotropy profile at the atomistic level if there is compositional disorder, which may influence localized magnetic objects like domain walls or skyrmions. Published by the American Physical Society 2024

25 ENERGY STORAGE↗
DOE OSTI2021

Advancements in modeling fuel pulverization and cladding behavior during a LOCA

During a loss-of-coolant accident (LOCA), it is possible for nuclear fuel rods to undergo a three-step process known as fuel fragmentation, relocation, and dispersal (FFRD). The chance of FFRD occurring increases as the fuel burnup increases. To support the nuclear industry's desire to increase the discharge burnup of nuclear fuels in light-water reactors (LWRs), it is imperative to understand the mechanisms driving the evolution of FFRD. In this work, a multiscale modeling approach is used to garner insight into underlying mechanisms leading to the fine fragmentation (also known as pulverization) of nuclear fuel during a LOCA. This report includes a summary of the atomistic and phase-field studies to develop a new pulverization criterion for use in the engineering-scale Bison fuel performance code. Details are also provided on cladding model improvements related to hydrogen/hydride embrittlement and damage and anisotropic thermal creep. The new models are used on the existing integral and separate effects LOCA validation cases available in Bison.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Atomistic Simulations for Thermophysical Properties of Uranium-Containing Halide Molten Salts

Characterizing the thermophysical properties in both fuel and coolant salts are critical in modeling, developing, process optimizing and utilizing molten salt reactors (MSRs), as these properties directly relate to operation metrics and can inform on the selection of candidate salts. The demand for consistent, accurate and publicly available thermophysical property data has become more apparent in recent years as interests have increased from molten salt reactor developers. There are a number of challenges in experimentally measuring properties such as thermal conductivity, viscosity, density and heat capacity , which have led to sparse and often times conflicting data points or molten salts in general. Additionally, there are a number of hazards to consider when synthesizing, storing, using, treating and disposing of molten salts. With the advances in computational capabilities over the last 10 years, the use of atomistic simulations can be implemented to support these efforts. The primary objective of this work is characterize the thermophysical transport properties in a number of molten chloride salts, and in particular NaCl-UCl 3 using ab-initio molecular dynamic (AIMD) simulations. In this binary salt the UCl 3 acts as the primary fissile material and NaCl acts as a carrier salt due with its’ high solubility for actinides A number of studies on the thermophysical properties of NaCl-UCl 3 have been published but there is not a vast amount of viscosity data for this system. In 1975, Desyatnik, et al published a study reporting dynamic viscosities that were calculated from kinematic viscosity measurements, and using the coefficients provided the viscosity in a 70:30 NaCl:UCl 3 mixture is 2.29 cP and 2.88 for a 60:40 mixture. Termini et al. recently reported viscosities in the range of 2.75 – 3 cP for the 63:37 NaCl-UCl 3 mixture in the same temperature range using rolling ball viscosity measurements. Computational viscosity of a similar mixture (64:36) can be obtained from the work Andersson et al. using the reported diffusion coefficients, and the hydrodynamic radius from the pair-radial distribution functions (RDFs). Using Eq (1) (vida infra), the viscosity would be 2.50 cP at 1100K. This is not to say that these values are incorrect due to the varying reported values, but aims to highlight the necessity of this work. The data reported in this ongoing work are computations on a 64:36 mixture of NaCl-UCl 3 at 987K. This work is likely to be expanded into varying concentrations of this mixture along with the inclusion of other salt candidate mixtures.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Multi-Scale Modeling and Prototype Development for Electrochemical CO2 Reduction (CRADA Final Report)

In this CRADA project, Lawrence Livermore National Laboratory, Stanford University, SLAC National Laboratory, and TotalEnergies collaboratively executed a multidisciplinary investigation of electrochemical reduction of CO2 to produce sustainable fuels and chemicals. Overall, the project led to an increased understanding of the fundamental processes involved in CO2 electrolysis, from the atomistic scale to the full electrolyzer device scale, ultimately leading to design guidelines for CO2 electrolyzers that will help in their future commercialization. As the model systems, Ag- and Cu-based catalysts were investigated in various forms depending on the electrochemical platform that was utilized to study the activity, selectivity, and durability towards electrochemical CO2 reduction. By employing experimental, theoretical, and computational techniques, the project team experimentally validated multi-physics models, evaluated the experimental levers that lead to increased electrolyzer reaction selectivity and energy efficiency, and used computational optimization to design higher performance electrodes. The learnings of this project were extensively documented in publicly available peer-reviewed journal publications and conference presentations, which serve as a foundation for further work to build from.

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

Predicting and Controlling Corrosion (Abbreviated Final Report)

Corrosion accumulates significant replacement costs in the transportation, utility, manufacturing, and infrastructure sectors. It also underpins several of LLNL’s core missions in stockpile stewardship, defense, and energy security. Whereas most corrosion and aging models are empirically parameterized to describe well-defined conditions late in the reaction, the factors that determine the early stages of corrosion—during which mitigation could be most impactful— are poorly understood. This project addressed the critical need for new approaches to predict the kinetics of corrosion initiation based on firm physical and chemical understanding. The activities encompassed degradation of relevant metals in both hydrogen-rich and environmental corrosion scenarios. The team integrated state-of-the-art multiscale simulation, in situ characterization, and data science within three technical thrusts: hydriding of Ti alloys; aqueous corrosion of Al and Ni-Cr alloys; and degradation of additively manufactured 316L stainless steel. In each case, novel capabilities were developed to identify and track the impacts of key atomistic, compositional, and microstructural features on the metal systems. For hydriding, protocols were developed to tightly integrate multiscale models, advanced multimodal characterization, and machine learning to determine how hydrogen interacts with native passivating surface oxides and nucleates new undesired phases, shedding new light on the critical role of grain boundaries, interfaces, and atomically disordered regions. For aqueous corrosion, the project demonstrated methods to predict dissolution rates of metal surfaces in corrosive solutions, to measure and understand microstructural and grain orientation effects on corrosion susceptibility, and to investigate competing growth and dissolution kinetics of surface oxides. For additively manufactured metals, analysis using state-of-the-art microscopy techniques revealed the role of specific heterogeneities invoked during laser processing, including cellular structure, dislocations, and precipitates, on corrosion susceptibility. In addition to new capabilities and understanding, the project provided an avenue for workforce development, as well as key partnerships with stakeholders in corrosion science.

08 HYDROGEN↗
DOE OSTI2022

Maximum strength and dislocation patterning in multi–principal element alloys

Multi–principal element alloys (MPEAs) containing three or more components in high concentrations render a tunable chemical short-range order (SRO). Leveraging large-scale atomistic simulations, we probe the limit of Hall-Petch strengthening and deformation mechanisms in a model CrCoNi alloy and unravel chemical ordering effects. The presence of SRO appreciably increases the maximum strength and lowers the propensity for faulting and structure transformation, accompanied by intensification of planar slip and strain localization. Deformation grains exhibit notably different microstructures and dislocation patterns that prominently depend on their crystallographic orientation and the number of active slip planes. Grain of single-planar slip attains the highest volume fraction of deformation-induced structure transformation, and grain with double-slip planes develops the densest dislocation network. These results advancing the fundamental understanding of deformation mechanisms and dislocation patterning in MPEAs suggest a mechanistic strategy for tuning mechanical behavior through simultaneously tailoring grain texture and local chemical order.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multiscale thermal properties prediction in the Multiphysics Object Oriented Simulation Environment (MOOSE) via a general Boltzmann solver [Poster]

The project objectives were: (1) Engineering-scale fuel performance modeling relies on accurate thermal properties; (2) Thermal properties are inherently multi-scale, arising from atomistic processes and interactions with a material’s microstructure; (3) Heat transport in solids via conduction occurs through transport and scattering of electrons and phonons; (4) Use the Boltzmann transport equation (BTE) to predict the macroscopic behavior of a materials system in terms of the microscopic dynamics of its heat carriers; and (5) This project establishes a new MOOSE (Multiphysics Object Oriented Simulation Environment) module, Boltzmann, dedicated to phonon and thermal electron transport.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Active learning for robust, high-complexity reactive atomistic simulations

Machine learned reactive force fields based on polynomial expansions have been shown to be highly effective for describing simulations involving reactive materials. Nevertheless, the highly flexible nature of these models can give rise to a large number of candidate parameters for complicated systems. In these cases, reliable parameterization requires a well-formed training set, which can be difficult to achieve through standard iterative fitting methods. In this paper, we present an active learning approach based on cluster analysis and inspired by Shannon information theory to enable semi-automated generation of informative training sets and robust machine learned force fields. The use of this tool is demonstrated for development of a model based on linear combinations of Chebyshev polynomials explicitly describing up to four-body interactions, for a chemically and structurally diverse system of C/O under extreme conditions. We show that this flexible training database management approach enables development of models exhibiting excellent agreement with Kohn–Sham density functional theory in terms of structure, dynamics, and speciation.

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

Decoding defect statistics from diffractograms via machine learning

Abstract Diffraction techniques can powerfully and nondestructively probe materials while maintaining high resolution in both space and time. Unfortunately, these characterizations have been limited and sometimes even erroneous due to the difficulty of decoding the desired material information from features of the diffractograms. Currently, these features are identified non-comprehensively via human intuition, so the resulting models can only predict a subset of the available structural information. In the present work we show (i) how to compute machine-identified features that fully summarize a diffractogram and (ii) how to employ machine learning to reliably connect these features to an expanded set of structural statistics. To exemplify this framework, we assessed virtual electron diffractograms generated from atomistic simulations of irradiated copper. When based on machine-identified features rather than human-identified features, our machine-learning model not only predicted one-point statistics (i.e. density) but also a two-point statistic (i.e. spatial distribution) of the defect population. Hence, this work demonstrates that machine-learning models that input machine-identified features significantly advance the state of the art for accurately and robustly decoding diffractograms.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Diffusion, atomic transport, and ordering in Al-Zr alloys: FCC and liquid phases

Additive manufacturing of materials with controlled microstructure demands knowledge of atomic scale properties near the solid-liquid transition state. Many of these properties are not affordable by experimental techniques and computer modeling is the possible solution to the problem. In this paper, we present the results of an extended atomistic study of intrinsic atomic transport due to vacancy diffusion in FCC and L12 solid phases and diffusion in the liquid phase of Al-Zr alloys. A deceleration of the overall self-diffusion was observed when Zr was added to Al. The effect was stronger in the solid and weaker in the liquid. Additionally, the effect was strongly temperature dependent in the solid phases, but not in the liquid. Atomic transport was chemically biased: transport of Zr atoms was significantly slower than that of Al atoms, and this bias effect was stronger in the solid phases. The overall diffusion and chemical ordering processes in the liquid state were five to six orders in magnitude faster than in the solid. Chemical short-range order parameters in the liquid saturated at values close to those in the ordered L12 structure of Al3Zr. Chemical and structural ordering in the solid phases was negligible over the modeled microsecond time scale. Here, the results are discussed in view of optimizing additive manufacturing parameters for the controlled formation of metastable L1 2 precipitates.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Multiscale analysis of large twist ferroelectricity and swirling dislocations in bilayer hexagonal boron nitride

With its atomically thin structure and intrinsic ferroelectric properties, heterodeformed bilayer hexagonal boron nitride (hBN) has gained prominence in next-generation non-volatile memory applications. However, studies to date have focused almost exclusively on small-twist bilayer hBN, leaving the question of whether ferroelectricity can persist under small heterostrain and large heterodeformation entirely unexplored. In this work, we establish the crystallographic origin of ferroelectricity in bilayer hBN configurations heterodeformed relative to high-symmetry configurations such as AA-stacking and 21.786789° twisted configurations (Σ7), using Smith normal form bicrystallography. We then demonstrate out-of-plane ferroelectricity in bilayer hBN across configurations vicinal to both the AA and Σ7 stackings. Atomistic simulations reveal that AA-vicinal systems support ferroelectricity under both small twist and small strain, with polarization switching in the latter governed by the deformation of swirling dislocations rather than the straight interface dislocations seen in the former. For Σ7-vicinal systems, where existing interatomic potentials underperform particularly under extreme out-of-plane compression, we develop a density-functional-theory-informed continuum framework—the bicrystallography-informed frame-invariant multiscale (BFIM) model, which captures out-of-plane ferroelectricity in heterodeformed configurations vicinal to Σ7 stacking. Interface dislocations in these large heterodeformed bilayer configurations exhibit markedly smaller Burgers vectors compared to interface dislocations in small-twist and small-strain bilayer hBN. The BFIM model reproduces experimental results and provides a powerful, computationally efficient framework for predicting ferroelectricity in large-unit-cell heterostructures where atomistic simulations are prohibitively expensive.

Ahmed, Md Tusher [Univ. of Illinois at Urbana-Cham↗
DOE OSTI2025

Artificial Intelligence and Multiscale Modeling for Sustainable Biopolymers and Bioinspired Materials

Abstract Biopolymers and bioinspired materials contribute to the construction of intricate hierarchical structures that exhibit advanced properties. The remarkable toughness and damage tolerance of such multilevel materials are conferred through the hierarchical assembly of their multiscale (i.e., atomistic to macroscale) components and architectures. Here, the functionality and mechanisms of biopolymers and bio‐inspired materials at multilength scales are explored and summarized, focusing on biopolymer nanofibril configurations, biocompatible synthetic biopolymers, and bio‐inspired composites. Their modeling methods with theoretical basis at multiple lengths and time scales are reviewed for biopolymer applications. Additionally, the exploration of artificial intelligence‐powered methodologies is emphasized to realize improvements in these biopolymers from functionality, biodegradability, and sustainability to their characterization, fabrication process, and superior designs. Ultimately, a promising future for these versatile materials in the manufacturing of advanced materials across wider applications and greater lifecycle impacts is foreseen.

Wang, Xing Quan [Department of Mechanical Engineer↗
DOE OSTI2023

2023 Roadmap on molecular modelling of electrochemical energy materials

New materials for electrochemical energy storage and conversion are the key to the electrification and sustainable development of our modern societies. Molecular modelling based on the principles of quantum mechanics and statistical mechanics as well as empowered by machine learning techniques can help us to understand, control and design electrochemical energy materials at atomistic precision. Therefore, this roadmap, which is a collection of authoritative opinions, serves as a gateway for both the experts and the beginners to have a quick overview of the current status and corresponding challenges in molecular modelling of electrochemical energy materials for batteries, supercapacitors, CO 2 reduction reaction, and fuel cell applications.

25 ENERGY STORAGE↗
DOE OSTI2021

Size-dependent radiation damage mechanisms in nanowires and nanoporous structures

Nanostructures with a high density of interfaces, such as in nanoporous materials and nanowires, resist radiation damage by promoting the annihilation and migration of defects. This study details the size effect and origins of the radiation damage mechanisms in nanowires and nanoporous structures in model face-centered (gold) and body-centered (niobium) cubic nanostructures using accelerated multi-cascade atomistic simulations and in-situ ion irradiation experiments. Our results reveal three different size-dependent mechanisms of damage accumulation in irradiated nanowires and nanoporous structures: sputtering for very small nanowires and ligaments, the formation and accumulation of point defects and dislocation loops in larger nanowires, and a face-centered-cubic to hexagonal-close-packed phase transformation for a narrow range of wire diameters in the case of gold nanowires. Smaller nanowires and ligaments have a net effect of lowering the radiation damage as compared to larger wires that can be traced back to the fact that smaller nanowires transition from a rapid accumulation of defects to a saturation and annihilation mechanism at a lower dose than larger nanowires. These irradiation damage mechanisms are accompanied with radiation-induced surface roughening resulting from defect-surface interactions. Comparisons between nanowires and nanoporous structures show that the various mechanisms seen in nanowires provide adequate bounds for the defect accumulation mechanisms in nanoporous structures with the difference attributed to the role of nodes connecting ligaments in nanoporous structures. Altogether, our results shed light on the compounded, size-dependent mechanisms leading to the radiation resistance of nanowires and nanoporous structures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Point defect evolution under irradiation: Finite size effects and spatio-temporal correlations

The evolution of point defect concentrations under irradiation is controlled by their diffusion properties, and by their formation and elimination mechanisms. The latter includes the mutual recombination of vacancies and interstitials, and the elimination of point defects at sinks. Two models are traditionally used to predict the defect concentration evolution, the standard rate theory (SRT) and the atomistic kinetic Monte Carlo (AKMC). We show in this work a large discrepancy in the defect concentrations between both methods when the average number of defects in the AKMC simulation box is close or lower than one. The reason is that AKMC naturally captures strong space and time correlations between vacancies and interstitials generated by the finite size of the periodic simulation box. These correlations strongly affect the recombination rate and the point defect concentrations. SRT fails to predict such correlations, and the corresponding solution deviates from the more accurate solution given by the AKMC simulation under similar conditions. These finite size correlation effects are strong when the elimination of point defects occurs by recombination only, but can still be significant in the presence of sinks. In order to account for the spacio-temporal correlation in a continuum framework, we introduce a Correlated Pair Theory (CPT). This theory fully takes into account the correlations between vacancy and interstitial pairs and predicts point defect concentrations in good agreement with AKMC simulations. Inversely, if the goal is for the AKMC to reproduce defect concentrations in bulk using small simulation boxes, the naturally occurring correlations need to be corrected. We show here that the CPT can be used to modify the elimination rates in the AKMC simulations, so as to yield point defect concentrations in agreement with SRT.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

PLAT domain protein 1 (PLAT1/PLAFP) binds to the Arabidopsis thaliana plasma membrane and inserts a lipid

Robust agricultural yields depend on the plant's ability to fix carbon amid variable environmental conditions. Over seasonal and diurnal cycles, the plant must constantly adjust its metabolism according to available resources or external stressors. The metabolic changes that a plant undergoes in response to stress are well understood, but the long-distance signaling mechanisms that facilitate communication throughout the plant are less studied. The phloem is considered the predominant conduit for the bidirectional transport of these signals in the form of metabolites, nucleic acids, proteins, and lipids. Lipid trafficking through the phloem in particular attracted our attention due to its reliance on soluble lipid-binding proteins (LBP) that generate and solubilize otherwise membrane-associated lipids. The Phloem Lipid-Associated Family Protein (PLAFP) from Arabidopsis thaliana is generated in response to abiotic stress as is its lipid-ligand phosphatidic acid (PA). PLAFP is proposed to transport PA through the phloem in response to drought stress. To understand the interactions between PLAFP and PA, nearly 100 independent systems comprised of the protein and one PA, or a plasma membrane containing varying amounts of PA, were simulated using atomistic classical molecular dynamics methods. In these simulations, PLAFP is found to bind to plant plasma membrane models independent of the PA concentration. When bound to the membrane, PLAFP adopts a binding pose where W41 and R82 penetrate the membrane surface and anchor PLAFP. This triggers a separation of the two loop regions containing W41 and R82. Subsequent simulations indicate that PA insert into the β -sandwich of PLAFP, driven by interactions with multiple amino acids besides the W41 and R82 identified during the insertion process. Finally, fine-tuning the protein-membrane and protein-PA interface by mutating a selection of these amino acids may facilitate engineering plant signaling processes by modulating the binding response.

59 BASIC BIOLOGICAL SCIENCES↗