Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “microstructure simulation”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Ferroelectric phase-transition frustration near a tricritical composition point

Phase transition describes a mutational behavior of matter states at a critical transition temperature or external field. Despite the phase-transition orders are well sorted by classic thermodynamic theory, ambiguous situations interposed between the first- and second-order transitions were exposed one after another. Here, we report discovery of phase-transition frustration near a tricritical composition point in ferroelectric Pb(Zr 1-x Ti x )O 3 . Our multi-scale transmission electron microscopy characterization reveals a number of geometrically frustrated microstructure features such as self-assembled hierarchical domain structure, degeneracy of mesoscale domain tetragonality and decoupled polarization-strain relationship. Associated with deviation from the classic mean-field theory, dielectric critical exponent anomalies and temperature dependent birefringence data unveil that the frustrated transition order stems from intricate competition of short-range polar orders and their decoupling to long-range lattice deformation. With supports from effective Hamiltonian Monte Carlo simulations, our findings point out a potentially universal mechanism to comprehend the abnormal critical phenomena occurring in phase-transition materials.

42 ENGINEERING↗
DOE OSTI2023

Parallel simulation via SPPARKS of on-lattice kinetic and Metropolis Monte Carlo models for materials processing

Abstract SPPARKS is an open-source parallel simulation code for developing and running various kinds of on-lattice Monte Carlo models at the atomic or meso scales. It can be used to study the properties of solid-state materials as well as model their dynamic evolution during processing. The modular nature of the code allows new models and diagnostic computations to be added without modification to its core functionality, including its parallel algorithms. A variety of models for microstructural evolution (grain growth), solid-state diffusion, thin film deposition, and additive manufacturing (AM) processes are included in the code. SPPARKS can also be used to implement grid-based algorithms such as phase field or cellular automata models, to run either in tandem with a Monte Carlo method or independently. For very large systems such as AM applications, the Stitch I/O library is included, which enables only a small portion of a huge system to be resident in memory. In this paper we describe SPPARKS and its parallel algorithms and performance, explain how new Monte Carlo models can be added, and highlight a variety of applications which have been developed within the code.

36 MATERIALS SCIENCE↗
NASA NTRS1985

A computer simulation of thin film nucleation and growth: The Volmer-Weber case

The computer simulation of thin film nucleation and growth, which was previously performed for the case of single monolayer, was modified to include multilayer growth via the Volmer-Weber mechanism. The simulation results show that: (1) the kinetics of multilayer film growth is nearly identical to that of monolayer growth; (2) when no reevaporation takes place, the cluster density resulting from multilayer growth is higher at high coverage than that resulting from monolayer growth; (3) when reevaporation does take place, the cluster density resulting from multilayer growth is nearly identical to that resulting from monolayer growth. This is not due, however, to similarity in microstructure.

Salik, J.↗
NASA NTRS1989

Mechanical and SEM analysis of artificial comet nucleus samples

Since 1987 experiments dealing with comet nucleus phenomena have been carried out in the DFVLR space simulation chambers. The main objective of these experiments is a better understanding of thermal behavior, surface phenomena and especially the gas dust interaction. As a function of different sample compositions and exposure to solar irradiation (xenon-bulbs) crusts of different hardness and thickness were measured. The measuring device consists of a motor driven pressure foot (5 mm diameter), which is pressed into the sample. The applied compressive force is electronically monitored. The microstructure of the crust and dust residuals is investigated by scanning electron microscopy (SEM) techniques. Stress-depth profiles of an unirradiated and an irradiated model comet are given.

Thiel, K.↗
DOE OSTI2021

A crystal plasticity model with an atomistically informed description of grain boundary sliding for improved predictions of deformation fields

Crystal plasticity (CP) is a powerful meso-scale technique for deformation modeling in polycrystalline materials. Crystal plasticity models typically do not include an explicit description of grain boundary (GB) sliding, which could lead to inaccurate predictions of strain distributions, especially in the vicinity of GBs. In the present study, a CP model is developed that includes GB sliding as an additional deformation mechanism and models the interaction between slip and GB sliding. Atomistic simulations are used to formulate the constitutive model for GB sliding and its interaction with incident slip. The deformation fields and structure of the GBs are obtained directly from experimental characterization and are faithfully reproduced as bicrystal systems for molecular dynamics simulations. The GBs are modeled as random-type (non-coincidence, asymmetrical) boundaries as observed from experimental data. The underlying atomistic-scale mechanism for pure sliding of random GBs was found to be analogous to fluid-flow. The interaction between incident slip and a sliding GB caused local increases in stress concentration, which further led to a momentary increase in the local sliding rate. The displacement profiles at the sliding GBs computed from the CP model with sliding accommodation are 38% more accurate than the baseline CP model as quantified by the mean squared error between the simulations and experiment. Here these results help improve the sophistication and accuracy of deformation modeling by including physics-based descriptions for GB sliding and its interaction with slip, eventually leading to more reliable predictions of micromechanical quantities.

36 MATERIALS SCIENCE↗
NASA NTRS2026

Mechanical Properties of Carbon Fiber Reinforced Composites Exposed to Cryogenic Conditions and Space Radiation via Simulation and Testing

As NASA missions extend beyond low Earth orbit, increasing reliance is placed on carbon fiber reinforced polymer (CFRP) composites for spacecraft structures where mass efficiency, durability, and long-term reliability are critical. In service, these materials are subjected to a combination of ultraviolet radiation, vacuum, ionizing radiation, atomic oxygen, and extreme thermal excursions under sustained mechanical loading. Flight systems such as the Boeing Starliner and SpaceX Dragon employ external composite structures that will experience these environments for extended durations. Although prior spaceflight and ground studies have reported limited changes in bulk mechanical properties, the synergistic effects of these environments on composite microstructure, particularly at the fiber matrix interphase, remain insufficiently characterized and represent a potential qualification and reliability risk. This study investigates the effects of short-term cryogenic exposure on a radiation shielding carbon epoxy composite, SC2020, as a ground-based analog for space relevant thermal extremes. The SC2020 material system has previously flown on the International Space Station under the Materials International Space Station Experiment (MISSE) program. Composite specimens were exposed to liquid nitrogen for 6 and 24 hours and evaluated using a multiscale characterization framework that combined ASTM D3039 tensile testing, Atomic Force Microscopy (AFM) based interphase analysis, and helium gas permeability measurements. Tensile testing showed no statistically significant or permanent degradation in global strength or modulus following cryogenic exposure. In contrast, AFM measurements revealed reductions in interphase modulus, weakened adhesion, and increased nanoscale heterogeneity, indicating localized degradation mechanisms not captured by conventional bulk testing. Gas permeability measurements showed a progressive increase in helium diffusion with exposure duration, consistent with micro-void formation or partial interfacial debonding. The results indicate that cryogenic exposure initiates degradation at the fiber matrix interphase while leaving global mechanical properties largely unchanged over short durations. These findings underscore the importance of multiscale diagnostics for identifying early-stage damage mechanisms that may influence long term performance and qualification margins for spaceflight composite structures. The data presented establish a cryogenic baseline for comparison with forthcoming MISSE flight exposure results and support ongoing NASA Established Program to Stimulate Competitive Research (EPSCoR) efforts aimed at improving composite qualification methodologies, risk assessment, and reliability prediction for space environments.

composite reliability↗
DOE OSTI2020

Deformation behavior and phase transformation of nanotwinned Al/Ti multilayers

We report nanotwinned Al/Ti multilayers have exhibited size-dependent microstructure evolution and high strength. However, their deformation mechanisms are less well understood. In this work, we investigated the deformation mechanisms of nanotwinned Al/Ti multilayers with FCC/HCP layer interfaces by using in situ micropillar compression tests. Nanotwinned Al/Ti multilayers exhibit compressive strength up to 2.4 GPa and good work hardening capability. Post-compression TEM analyses reveal high-density stacking faults and the HCP-to-FCC phase transformations in Ti. Molecular dynamics simulations elucidate the mechanisms of deformation induced phase transformation in Ti and the influence of collective movement of partial dislocations on the deformability of Al/Ti multilayers.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Physics-informed machine learning exploration of Na storage mechanisms in disordered carbon

Sodium-ion batteries are a cost-effective, sustainable alternative to lithium-ion systems for large-scale energy storage. However, optimizing sodium storage in carbon-based anodes with microstructural complexity and atomic disorder remains a major challenge. The intrinsic inhomogeneity of these materials produces diverse local environments, making it difficult for conventional methods to predict and control ion dynamics. Hard carbon (HC) anodes, composed of ranges of ordered-to-disordered graphitic and amorphous nanodomains, offer tunable ion storage and rate capacity, yet rationale design remains a challenge due to poorly understood correlation between local atomic feature and ion transport mechanism. Here, to address this challenge, we introduce a data-driven framework that integrates validated machine-learned interatomic potentials, large-scale molecular dynamics simulations, and machine learning to elucidate sodium transport mechanisms as a function of carbon and sodium loading densities. By computing per-ion structural descriptors and applying unsupervised learning, we identify distinct diffusion modes governed by microscopic features. Supervised analysis and correlation mapping then establish quantitative links between these transport regimes and processing variables such as bulk carbon density and sodium content. This physics-informed approach establishes quantitative structure–transport relationships and offers actionable design principles for engineering high-performance HC anodes.

Data-driven framework↗
DOE OSTI2022

Corrosion behavior of a series of combinatorial physical vapor deposition coatings on SiC in a simulated boiling water reactor environment

Here, the simulated boiling water condition corrosion properties of multilayer and combinatorial physical vapor deposition metal coatings (Zr, Cr, and/or Ti) deposited on a SiC f -SiC composite were investigated. Various compositions within the ternary system were corroded to the most suitable composition preventing underlying SiC f -SiC from further attack. Mass measurements and visual observations were conducted to identify which coatings formed protective oxides and which coatings spalled. Transmission electron microscopy was conducted to observe the microstructure of the physical vapor deposition metal layer before and after corrosion as well as identify the oxide which were forming.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Understanding the fundamental driver of semiconductor radiation tolerance with experiment and theory

Space is the operating environment of a multitude of systems that our society is heavily reliant on today, however, maintaining operability there necessitates special consideration of the electronic systems' tolerance of space radiation. Electronic systems are critically dependent on the electronic properties of their semiconductor components, which are modified by space radiation with an adverse impact on the space system performance. What innate property allows some semiconductors to sustain little damage while others accumulate defects rapidly with dose is poorly understood, which limits the extent to which radiation tolerance can be implemented as a design criterion. Here, to gain insight into what properties are drivers of semiconductor radiation tolerance, the first step is to generate a dataset of the relative radiation tolerance of a broad sampling of semiconductors. To accomplish this, Rutherford backscatter channeling experiments are used to compare the displaced lattice atom buildup in InAs, InP, GaP, GaN, ZnO, MgO, and Si as a function of stepwise alpha particle dose. With this experimental information on radiation-induced incorporation of interstitial defects in hand, hybrid density functional theory electron densities (and their derived quantities) are calculated and their gradient and Laplacian are evaluated to obtain key fundamental information about the interactions in each material. It is shown that simple, undifferentiated values (which are typically used to describe bond strength) are insufficient to predict radiation tolerance. Instead, the curvature of the electron density at bond critical points provides a measure of radiation tolerance consistent with the experimental results obtained. This curvature and associated forces surrounding bond critical points have the potential to disfavor the localization of displaced lattice atoms at these points, favoring their diffusion toward perfect lattice positions. With this criterion to predict radiation tolerance, simple density functional theory simulations can be conducted on potential new materials to gain insight into how they may operate in demanding high radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hardness Prediction by Incorporating Heat Transfer and Molten Pool Fluid Flow in a Multi-pass, Multi-layer Weld for Onsite Repair of Grade 91 steel

Introduction: In the current fleet of fossil-fired power plants, creep strength enhanced ferritic steels (CSEF) are used to sustain the harsh service conditions. Enhanced properties of Grade 91 steel result from tempered martensite with a fine distribution of MX and M23C6 carbides. Grade 91 steel is subjected to onsite welding repair to remedy their degradation due to extreme service condition. Knowledge of weld repairability of these steels, such as as-welded hardness distribution, is essential to establishing sound repair procedures. Experimental trial and error tests can consume a lot of time as many welding variables need to be studied. For numerical modelling, most of the multi-pass multi-layer models are based on finite element method, which are limited to solve the heat conduction equation and ignore convective heat transfer due to melt flow. Moreover, the mesh has to be pre-built based on a known or assumed weld cross-section geometry. These finite element based models thus have limited predictive capability as defects are not considered and nugget size are pre-assumed. This research aims at developing a thermal and microstructure evolution model incorporating molten pool dynamics in a multi-pass multi-layer material deposition to predict the as-welded hardness distribution. Technical Approach: All the thermal, physical, and metallurgical properties of Grade 91 as a function of temperature are collected from the literature and inputted into the thermo-fluid model based on Flow-3D, a computational fluid dynamics software. A multi-pass, multi-layer material deposition is simulated where the melting of filler wire into the molten pool is directly considered based on the volume of fluid (VOF) method. The flow behaviour of the molten pool is used to understand the formation of deposition geometry and defects. The temperature profiles during the multi-pass, multi-layer welding are calculated. The results computed using the new model are compared against the experimental data of fusion zone geometry and thermal cycles. Hardness prediction in the heat-affected zone (HAZ) are made using Johnson-Mehl-Avrami (JMA) equation for solid-state phase transformation kinetics. The JMA parameters are extracted from the experimental data available in the literature. For comparison, a standard finite element heat conduction model is also developed to predict the thermal cycles and hardness distribution in the multi-pass, multi-layer weld. Expected Result: Results obtained using the molten pool dynamic simulation versus the finite element heat conduction model are compared. Specifically, the effects of convective heat transfer on the accuracy of the calculated thermal history, bead shape and size, and HAZ hardness distribution are examined.

Modeling, Grade 91 steel, Multi-pass multi-layer d↗
DOE OSTI2025

Computational Modeling of Molten Salt Infiltration and Oxidation in Nuclear Graphite

Graphite is utilized as a moderator and reflector in advanced nuclear reactor designs due to its high thermal conductivity, neutron moderation properties, and resistance to radiation damage. However, its longterm performance and reliability are challenged by degradation mechanisms such as molten salt infiltration in molten salt reactors (MSRs) and oxidation in gas-cooled reactors (GCRs). These mechanisms can compromise the structural integrity and operational lifetime of graphite components, necessitating a more detailed assessment of their physical behavior. This report focuses on the development of computational models for molten salt infiltration and oxidation of graphite to aid the design and performance analysis of graphite components. For molten salt infiltration, a computational framework is developed that couples incompressible Navier-Stokes and phase-field model to simulate the penetration of molten salt into graphite?s interconnected pore structure. Initial model verification is performed using two-phase flows in two dimensions, demonstrating the models ability to capture fundamental physical behavior and agree with analytical solution. This framework is then applied to a realistic IG110 nuclear graphite , where a computed tomography extracted pore geometry is used to analyse the infiltration behavior of FLiNaK molten salt. This model provides insights into how the microstructure and other relevant parameters influence the transport pathways of molten salt into graphite, potentially offering a means to rapidly evaluate a graphite grade?s resistance to infiltration. For oxidation, the report details pore-scale mass and heat transport models, describing the diffusion of gases, reaction kinetics, and thermal effects. Additionally, this report highlights inconsistencies in the existing volume-averaged macroscopic model, particularly in upscaling of reaction kinetics and flux terms, and surface to volume transformations. These inconsistencies suggest that current formulations may not accurately capture the experimentally observed graphite oxidation process, highlighting the need for improved model development. This work advances the development of physics-based computational models for graphite degradation, contributing to improved predictive models for next-generation nuclear reactor designs. Future efforts will focus on refining the infiltration model to address non-physical behaviors and enhance its robustness. Additionally, for oxidation, further studies will employ the principles of volume averaging to rigorously derive the upscaled equations, potentially in collaboration with subject matter experts.

Computational Modeling of Molten Salt Infiltration↗
DOE OSTI2021

A Continuum Dislocation Dynamics Crystal Plasticity Approach to Irradiated Body-Centered Cubic α-Iron

Radiation-induced embrittlement of reactor pressure vessel (RPV) steels can potentially limit the operating life of nuclear power plants. Over extended exposure to radiation doses, these body-centered cubic (BCC) irons demonstrate irradiation damage. Here, we present a continuum dislocation density (CDD) crystal plasticity model to capture the interaction among dislocations and self-interstitial atom (SIA) loops in α-iron. We demonstrate the importance of modeling cross slip using a combined stochastic Monte Carlo approach and the role of slip system strength anisotropy in capturing stochastic cross slip interactions. Through these captured interactions, the CDD crystal plasticity model can capture both the stress response and the physical evolution of dislocations on different slip system planes. Single-crystal verification experiments are used to calibrate the CDD crystal plasticity model, and a set of simplified polycrystalline simulations demonstrates the model’s ability to capture the stress response from tensile experiments on α-iron.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Controlling solute channel formation using magnetic fields

Solute channel formation introduces compositional and microstructural variations in a range of processes, from metallic alloy solidification, to salt fingers in ocean and water reservoir flows. Applying an external magnetic field interacts with thermoelectric currents at solid/liquid interfaces generating additional flow fields. This thermoelectric (TE) magnetohydrodynamic (TEMHD) effect can impact on solute channel formation, via a mechanism recently drawing increasing attention. To investigate this phenomenon, we combined in situ synchrotron X-ray imaging and Parallel-Cellular-Automata-Lattice-Boltzmann based numerical simulations to study the characteristics of flow and solute transport under TEMHD. Observations suggest the macroscopic TEMHD flow appearing ahead of the solidification front, coupled with the microscopic TEMHD flow arising within the mushy zone are the primary mechanisms controlling plume migration and channel bias. Two TE regimes were revealed, each with distinctive mechanisms that dominate the flow. Further, we show that grain orientation modifies solute flow through anisotropic permeability. These insights led to a proposed strategy for producing solute channel-free solidification using a time-modulated magnetic field.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modes of strain accommodation in Cu-Nb multilayered thin film on indentation and cyclic shear

Two-phase layered thin films with a high density of semi-coherent interfaces exhibit excellent mechanical properties and thermal stability. Here, in this study, a magnetron-sputtered Cu-Nb dual-layered thin film (~500 nm for Cu and ~150 nm for Nb) having an amorphous interface between Cu and Nb with a high density of aligned growth twins in Cu is subjected to severe surface deformation. The material is loaded using indentation and cyclic shear under tribological testing. The strain accommodation in the subsurface microstructure after deformation varies based on the local structure and deformation mode. Grain refinement and crack formations in the stressed region of the Nb layer and localized crystallization of the amorphous interface are observed after indentation and scratch testing. Pronounced detwinning of growth twins in the Cu layer under the cyclic shear strain leaves large dislocations sites and loops which are observed both by high-resolution transmission electron microscopy and experiment-guided molecular dynamic (MD) simulations. Our simulations provided insights into understanding the pathway for the detwinning process under cyclic shear loading.

36 MATERIALS SCIENCE↗
DOE OSTI2021

All-Atom Molecular Dynamics Simulations of the Temperature Response of Densely Grafted Polyelectrolyte Brushes

Water acts a “good” solvent for most polyelectrolyte (PE) systems due to its high dielectric strength. Theoretical and experimental studies have highlighted some remarkably unconventional properties of strongly water-swollen PE brushes such as a reduction in brush height with increase in temperature. However, an understanding of the atomistic-scale response of the grafted PE chains as well as the local arrangement and organization of the brush-supported counterions and water molecules to changes in temperature is missing. Here, we conduct an all-atom molecular dynamics (MD) study to probe the influence of temperature on the “microstructure” of densely grafted, fully ionized polyacrylic acid (PAA) brushes. The atomistic insights obtained from our study shed light on prior experimental observations and elucidate the key role played by changes in hydrogen bonding network in dictating the hydrophilicity of the brushes. In addition to the temperature-mediated reduction in height of the PE brush layer, we study the influence of temperature on several properties of the brush-supported water molecules such as the number distribution and kinetics of water-water and water-PE hydrogen bonds (HBs), mass density, orientational tetrahedral order parameter, and the activation energy associated with the self-diffusion of water. The effect of temperature on the properties of the brush-trapped counterions is also quantified via changes in their solvation structure, counterion-water radial distribution function (RDF), and the translational mobility. Our findings unravel several hitherto unknown phenomena such as the enhanced propensity for attaining PE brush-induced “water-in-salt”-like scenarios at elevated temperatures, and an increase in the activation energy for the self-diffusion of water with increase in the degree of PE brush-induced nanoconfinement (from ‘bulk water’ to ‘sparsely grafted PE brushes’ to ‘densely grafted PE brushes’). Further, we anticipate that this work will be helpful in improving our predictive capabilities of the temperature response of water-swollen PE brushes, which would be instrumental for better design of several PE brush systems used in a multitude of applications.

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

Multiscale formulation of frictional contact mechanics at the pore scale

Direct numerical simulation (DNS) yields the highest fidelity predictions of mechanical deformation at the pore scale, but is prohibitively expensive for analyzing large or many samples. Discrete element methods (DEM) are an efficient alternative, but are limited to granular media and incapable of estimating or controlling prediction errors. We present a pore-level multiscale method (PLMM) that approximates DNS efficiently and with controllable accuracy. We focus on the linear elastic response of a consolidated geologic porous medium with arbitrary microstructure, heterogeneous mineralogy, containing cracks or defects. PLMM decomposes the solid phase into non-overlapping subdomains, on which local basis functions are constructed. The bases are then coupled with a global interface problem that accounts for slip or stick contact conditions between the subdomains. PLMM produces an initial, but accurate, approximation to DNS that can be iteratively improved. It is amenable to parallelism and allows for different mesh, models, and physics in each subdomain. An algebraic interpretation of PLMM as a preconditioner is also presented to allow non-intrusive implementation into existing solvers. Lastly, this work extends previous developments of PLMM in fluid dynamics to solid mechanics and enables future extensions towards modeling coupled flow and mechanics problems.

42 ENGINEERING↗
DOE OSTI2026

A robust alloy design (RAD) strategy for next-generation (IV) nuclear fission reactors

Next-generation nuclear reactors demand structural materials capable of withstanding extreme conditions, including high temperatures, intense neutron flux, and corrosive environments. Multi-Principal Element Alloys (MPEAs) have emerged as promising candidates due to their exceptional radiation tolerance, thermal stability, and compositional flexibility. This study introduces a versatile and customizable Robust Alloy Design (RAD) strategy for systematically designing MPEAs for GEN-IV reactor fuel cladding. The RAD framework integrates nuclear-relevant selection criteria, empirical parameter assessments, and high-throughput CALPHAD simulations to efficiently narrow compositional space and identify stable alloys. A unified RAD score developed for the first time, combines key performance metrics, including fuel-clad chemical interaction (FCCI), neutron absorption cross-section (NAC), valence electron configuration (VEC), and melting point factor (MPF), into a flexible ranking system adaptable to reactor-specific priorities. Among 724 candidates, V555(5Al–5Cr–5Fe–85V) emerged as the top alloy, validated experimentally with a homogeneous single-phase BCC microstructure and superior mechanical properties (nano-indentation: 3.389 ± 0.258 GPa; Vickers hardness: 240 ± 6.7 HV), significantly outperforming Zircaloy-4 and V-4Cr-4Ti. Importantly, the RAD strategy is not limited to nuclear applications; its customizable weighting system enables scalability to other extreme environments. This adaptability positions RAD strategy as a versatile tool for advanced materials design across multiple industries.

Alloy design↗