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

Defect annihilation in heavy ion irradiated polycrystalline gold

In this study, we explore the interaction of electron wind force (EWF) with defects originating from ion irradiation in-situ inside a transmission electron microscope. Nanocrystalline gold specimens were self-ion irradiated to a dose of 5 × 10 15 ions/cm 2 (45 displacement per atom) to generate a high density of displacement damage. We also developed a molecular dynamics simulation model to understand the associated atomic scale mechanisms. Both experiments and simulations show that the EWF can impart significant defect mobility even at low temperatures, resulting in the migration and elimination of defects in a few minutes. Overall, we propose that the EWF interacts with defects to create highly glissile Shockley partial dislocations, which makes the fast and low temperature defect annihilation possible.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Predicting structural material degradation in advanced nuclear reactors with ion irradiation

Swelling associated with the formation and growth of cavities is among the most damaging of radiation-induced degradation modes for structural materials in advanced nuclear reactor concepts. Ion irradiation has emerged as the only practical option to rapidly assess swelling in candidate materials. For decades, researchers have tried to simulate the harsh environment in a nuclear reactor in the laboratory at an accelerated rate. Here we present the first case in which swelling in a candidate alloy irradiated ~ 2 years in a nuclear reactor was replicated using dual ion irradiation in ~ 1 day with precise control over damage rate, helium injection rate, and temperature and utilize physical models to predict the effects of radiation in reactors. The capability to predict and replicate the complex processes surrounding cavity nucleation and growth across many decades of radiation dose rate highlights the potential of accelerated radiation damage experiments. More importantly, it demonstrates the capability to predict the swelling evolution and the possibility to predict other features of the irradiated microstructure evolution that control material property degradation required to accelerate the development of new, radiation-tolerant materials.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Structure and dynamics of $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquids close to Earth's liquid core conditions

Using an artificial neural-network machine learning interatomic potential, we have performed molecular dynamics simulations to study the structure and dynamics of $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquid close to the Earth's liquid core conditions. Here, the simulation results reveal that the short-range structural order (SRO) in the $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquid is very strong. About 80% of the atoms are arranged in crystallinelike SRO motifs. In particular, ~ $70$% of Fe-centered clusters can be classified as either hexagonal-close-packed-like or icosahedrallike SRO motifs. The SRO clusters centered on Fe, Si, or O atoms are strongly intermixed and homogenously distributed throughout the liquid. The atomic structure of the liquid and the fractions of dominant SRO clusters are not sensitive to pressure/temperature used in the simulations except that the SRO of the O-centered clusters is enhanced close to inner core pressures. The O-diffusion coefficient is about two to three times larger than the Fe and Si ions and increases more rapidly in the deeper core regions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Toward Addressing the Challenge to Predict the Heat Capacities of RDX and HMX Energetic Materials

Availability of heat capacity as function of pressure and temperature is an essential prerequisite for development of a computational multiscale strategy capable to address the evolution of microstructure and energy release in advanced high energy density materials. In the case of 1,3,5-trinitro-1,3,5-triazinane (RDX) and 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) systems as two of the most studied energetic materials, there are substantial gaps in experimental data, with available heat capacities values distributed only in a region close to standard ambient conditions. In this study we demonstrate how these major experimental limitations can be addressed in the case of the RDX and HMX systems based on the combined use of classical and quantum mechanical calculations. We show that by considering ideal gas properties evaluated using quantum mechanical methods, and residual properties obtained from molecular simulations using fully flexible atomistic force field models, excellent agreement can be obtained for the predicted heat capacities to the most recent experimental values. An important advantage of the current computational methodology is that it allows evaluation of both constant-volume and constant-pressure heat capacities for a broad interval of temperatures and pressures, which encompasses solid and liquid phases conditions. In the case of the solid α and γ phases of RDX and the β phase of HMX, the predicted results follow closely both the available experimental data at standard ambient conditions and the results obtained using density functional theory calculations at high pressures, a regime where experimental data are not available. A perspective to expand the current methodology is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Irradiation-induced gas production in REBCO-based magnet materials used for future compact fusion reactors

Nuclear fusion is an enticing alternative to current sources of energy, with multilayered Rare-Earth Barium Copper Oxide (REBCO) coated conductors deemed pivotal in the race toward fully realized, commercially viable, and magnetic confinement fusion reactors. In this study, we simulated the ion spectrum expected to evolve from REBCO's nickel-based Hastelloy C-276 substrate and copper stabilizer in an affordable robust compact-like reactor. We then emulated this gas production through helium implantation to investigate changes in materials and superconducting properties. Our results revealed that the substrate and stabilizer are capable of producing protons energetic enough to recoil throughout the tape thickness in appreciable doses, and alphas energetic enough to deposit 7.54 × 1014 ions/cm2 or 50.1 helium appm in the superconducting layer over a 30-year reactor lifetime. The superconducting layer of SuperPower® tapes exhibited at least double the swelling rate of the other major layers, and both SuperPower and Fujikura Ltd. tapes displayed microstructural changes in the REBCO layer not observed in isotropic metals. For the estimated lifetime fluence, the Fujikura tapes showed a ∼1 K reduction in critical temperature and a 32% degradation in critical current for compact reactor-relevant conditions (16 T, 20 K). Nuclear transmutation, low-temperature solder implantations, gas-ion evolution, the influence of gas production on vortex dynamics, and other related considerations are also discussed.

Reis, Chris↗
DOE OSTI2020

Data-driven learning of nonlocal models: from high-fidelity simulations to constitutive laws

We show that machine learning can improve the accuracy of simulations of stress waves in one-dimensional composite materials. We propose a data-driven technique to learn nonlocal constitutive laws for stress wave propagation models. The method is an optimization-based technique in which the nonlocal kernel function is approximated via Bernstein polynomials. The kernel, including both its functional form and parameters, is derived so that when used in a nonlocal solver, it generates solutions that closely match high-fidelity data. The optimal kernel therefore acts as a homogenized nonlocal continuum model that accurately reproduces wave motion in a smaller-scale, more detailed model that can include multiple materials. We apply this technique to wave propagation within a heterogeneous bar with a periodic microstructure. Several one-dimensional numerical tests illustrate the accuracy of our algorithm. The optimal kernel is demonstrated to reproduce high-fidelity data for a composite material in applications that are substantially different from the problems used as training data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Microstructurally informed synchrotron x-ray analysis revealing helium defect transitions in ultrafine grained tungsten

The formation of insoluble gaseous defects in materials due to nuclear transmutation or ion implantation involves the diffusion of impurity atoms to form atomic defect clusters that coalesce into bubbles or cavities and ultimately degrade the material properties. Transmission electron microscopy (TEM) is limited in its ability to resolve sub-nanometer gas clusters whereas X-ray diffraction (XRD) provides information pertaining to local atomic changes. Here, in this study, helium (He) implanted ultrafine grained tungsten is explored through a multimodal defect characterization campaign combining TEM-informed Small Angle X-ray Scattering (SAXS) analysis, XRD lattice parameter measurements, and nanoscale He cluster quantification from a region of reciprocal space accessible via Wide Angle X-ray Scattering (WAXS). Moderate elevated temperature implantations are shown to produce high concentrations of sub-nanoscale He clusters and small, homogeneously distributed cavities, which collectively are linked to lattice expansion and further substantiated through complementary atomistic simulations. Increased implantation temperatures encourage the diffusion of these defects to the grain boundaries (GBs), leading to lattice relaxation and the growth of large GB cavities manifesting as bimodal size distributions in the SAXS analysis. Overall, our results demonstrate the utility of multimodal synchrotron X-ray analysis in bridging the gap between microscale He cavity quantification and atomic-scale defect analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Documentation, assessment, and improvements of fission gas modeling capabilities in BISON

Fission gas behavior critically affects nuclear fuel performance, influencing thermal conductivity, microstructure, swelling, and fuel-cladding interactions. The Simple Integrated Fission Gas Release and Swelling (Sifgrs) model, implemented within the BISON fuel performance code, provides advanced capabilities for modeling these phenomena. This manuscript reviews and assesses Sifgrs’ modular structure and submodels describing intragranular, intergranular, and release mechanisms. The modular implementation enables rapid improvements as our understanding advances and new fuel forms emerge. After an overview of capabilities, we present a comprehensive validation suite, comparing predictions to experimental data from tests under normal and transient conditions, including reproduction of the Vitanza curve. Results demonstrate Sifgrs’ accuracy, robustness, and applicability. Current limitations and ongoing efforts are discussed, including mechanistic, multiscale modeling, uncertainty quantification, and validation within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, supporting the safe and efficient use of nuclear fuels.

BISON↗
DOE OSTI2023

Combining synchrotron X-ray diffraction, mechanistic modeling and machine learning for in situ subsurface temperature quantification during laser melting

Laser melting, such as that encountered during additive manufacturing, produces extreme gradients of temperature in both space and time, which in turn influence microstructural development in the material. Qualification and model validation of the process itself and the resulting material necessitate the ability to characterize these temperature fields. However, well established means to directly probe the material temperature below the surface of an alloy while it is being processed are limited. To address this gap in characterization capabilities, a novel means is presented to extract subsurface temperature-distribution metrics, with uncertainty, from in situ synchrotron X-ray diffraction measurements to provide quantitative temperature evolution data during laser melting. Temperature-distribution metrics are determined using Gaussian process regression supervised machine-learning surrogate models trained with a combination of mechanistic modeling (heat transfer and fluid flow) and X-ray diffraction simulation. The trained surrogate model uncertainties are found to range from 5 to 15% depending on the metric and current temperature. The surrogate models are then applied to experimental data to extract temperature metrics from an Inconel 625 nickel superalloy wall specimen during laser melting. The maximum temperatures of the solid phase in the diffraction volume through melting and cooling are found to reach the solidus temperature as expected, with the mean and minimum temperatures found to be several hundred degrees less. The extracted temperature metrics near melting are determined to be more accurate because of the lower relative levels of mechanical elastic strains. However, uncertainties for temperature metrics during cooling are increased due to the effects of thermomechanical stress.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Ultrasonic Nondestructive Characterization of Adhesive Bonds

Adhesives and adhesive joints are widely used in various industrial applications to reduce weight and costs, and to increase reliability. For example, advances in aerospace technology have been made possible, in part, through the use of lightweight materials and weight-saving structural designs. Joints, in particular, have been and continue to be areas in which weight can be trimmed from an airframe through the use of novel attachment techniques. In order to save weight over traditional riveted designs, to avoid the introduction of stress concentrations associated with rivet holes, and to take full advantage of advanced composite materials, engineers and designers have been specifying an ever-increasing number of adhesively bonded joints for use on airframes. Nondestructive characterization for quality control and remaining life prediction has been a key enabling technology for the effective use of adhesive joints. Conventional linear ultrasonic techniques generally can only detect flaws (delamination, cracks, voids, etc) in the joint assembly. However, more important to structural reliability is the bond strength. Although strength, in principle, cannot be measured nondestructively, a slight change in material nonlinearity may indicate the onset of failure. Furthermore, microstructural variations due to aging or under-curing may also cause changes in the third order elastic constants, which are related to the ultrasonic nonlinear parameter of the polymer adhesive. It is therefore reasonable to anticipate a correlation between changes in the ultrasonic nonlinear acoustic parameter and the remaining bond strength. It has been observed that higher harmonics of the fundamental frequency are generated when an ultrasonic wave passes through a nonlinear material. It seems that such nonlinearity can be effectively used to characterize bond strength. Several theories have been developed to model this nonlinear effect. Based on a microscopic description of the nonlinear interface binding force, a quantitative method was presented. Recently, a comparison between the experimental and simulated results based on a similar theoretical model was presented. A through-transmission setup for water immersion mode-converted shear waves was used to analyze the ultrasonic nonlinear parameter of an adhesive bond. In addition, ultrasonic guided waves have been used to analyze adhesive or diffusion bonded joints. In this paper, the ultrasonic nonlinear parameter is used to characterize the curing state of a polymer/aluminum adhesive joint. Ultrasonic through-transmission tests were conducted on samples cured under various conditions. The magnitude of the second order harmonic was measured and the corresponding ultrasonic nonlinear parameter was evaluated. A fairly good correlation between the curing condition and the nonlinear parameter is observed. The results show that the nonlinear parameter might be used as a good indicator of the cure state for adhesive joints.

Qu, Jianmin↗
DOE OSTI2026

Reduced‐Order Modeling of Energetic Materials Using Physics‐Aware Recurrent Convolutional Neural Networks in a Latent Space (LatentPARC)

Physics-aware deep learning (PADL) has gained popularity for use in spatiotemporal dynamics simulations, such as those in computational modeling of energetic materials (EM). We show that the challenge PADL methods face while learning complex field evolution problems can be simplified and accelerated by decoupling it into two tasks: learning complex geometric features in evolving fields and modeling dynamics over these features in a lower-dimensional feature space. We build upon our previous work on physics-aware recurrent convolutional neural networks (PARC). PARC embeds knowledge of underlying physics into its neural network architecture for more robust and accurate prediction of evolving physical fields. PARC was shown to effectively learn complex nonlinear features such as the formation of hotspots and coupled shock fronts in various initiation scenarios of EMs, as a function of microstructures, serving effectively as a microstructure-aware burn model. Here, we further accelerate PARC and reduce its computational cost by projecting the original dynamics onto a lower-dimensional invariant manifold, or “latent space.” The projected latent representation encodes the complex geometry of evolving fields (e.g., temperature and pressure) in a set of data-driven features. The reduced dimension of this latent space allows us to learn the dynamics during the initiation of EM with a lighter and more efficient model. We observe a significant decrease in training and inference time while maintaining results comparable to PARC at inference. This work takes steps towards enabling rapid prediction of EM thermomechanics at larger scales and characterization of EM structure–property–performance linkages at a full application scale.

Mathematics and Computing↗
DOE OSTI2021

Hydrogen permeability and stability of palladium films in a high pressure reactor

Palladium (Pd) thin films were fabricated by colloidal spray deposition (CSD) on substrates of porous yttrium stabilized zirconium oxide (YSZ). The films were sintered at 1140–1160 °C for 2–3 h in argon and had thickness of 5–11 μm. The hydrogen (H 2 ) permeation flux for the Pd film was measured with simulated feed gases (90% H 2 /He, syn-gas/steam, or syn-gas/steam with H 2 S) flowing on the film side and sweep gas (Ar) flowing on the substrate side. Here, the hydrogen flux with the Pd film was obtained by adjusting the flow rate and pressure of the gas. Hydrogen flux values of 121–173 cm 3 min -1 cm -2 were measured at 388–596 °C (feed: 90% H 2 /He at 200 mL min -1 ; sweep: Ar at 1400 mL min -1 ; and ΔP = ~170 psi), and the permeability was calculated and compared with literature values. Because Pd is known to be sensitive to some gases, the stability of Pd films was evaluated by examining the microstructures of samples after their exposure to gas mixtures containing CO, CO 2 , H 2 O, and H 2 S. From the results of stability tests in H 2 S-containing atmospheres, the Pd/Pd 4 S phase boundary was determined, and this measured boundary was compared to the Pd/Pd 4 S boundary calculated from the thermodynamic data in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Demonstration of neutron time-of-flight diffraction with an event-mode imaging detector

Neutron diffraction beamlines have traditionally relied on deploying large detector arrays of 3 He tubes or neutron-sensitive scintillators coupled with photomultipliers to efficiently probe crystallographic and microstructure information of a given material. Given the large upfront cost of custom-made data acquisition systems and the recent scarcity of 3 He, new diffraction beamlines or upgrades to existing ones demand innovative approaches. This paper introduces a novel Timepix3-based event-mode imaging neutron diffraction detector system as well as first results of a silicon powder diffraction measurement made at the HIPPO neutron powder diffractometer at the Los Alamos Neutron Science Center. Notably, these initial measurements were conducted simultaneously with the 3 He array on HIPPO, enabling direct comparison. Data reduction for this type of data was implemented in the MAUD code, enabling Rietveld analysis. Results from the Timepix3-based setup and HIPPO were benchmarked against McStas simulations, showing good agreement for peak resolution. With further development, systems such as the one presented here may substantially reduce the cost of detector systems for new neutron instrumentation as well as for upgrades of existing beamlines.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1978

Shear rupture of a directionally solidified eutectic gamma/gamma prime - alpha (Mo) alloy

Directionally solidified Mo alloys are evaluated to determine the shear rupture strength and to possibly improve it by microstructural and heat treatment variations. Bars of the alloy containing nominally 5.7% Al and 33.5% Mo by weight with balance Ni were directionally solidified at rates between 10 and 100 mm per hour in furnaces with thermal gradients at the liquid-solid interface of 250 or 100 C per cm. A limited number of longitudinal shear rupture tests were conducted at 760 C and 207 MPa in the as - solidified and in several heat treated conditions. It is shown that shear rupture failures are partly transgranular and that resistance to failure is prompted by good fiber alignment and a matrix structure consisting mainly of gamma prime. Well aligned as - solidified specimens sustained the shear stress for an average of 81 hours. A simulated coating heat treatment appeared to increase the transformation of gamma to gamma prime and raised the average shear life of aligned specimens to 111 hours. However, heat treatments at 1245 C and especially at 1190 C appeared to be detrimental by causing partial solutioning of the gamma prime, and reducing lives to 47 and 10 hours, respectively.

Harf, F. H.↗
NASA NTRSDate not supplied

Multiscale Modeling of Woven Ablative Thermal Protection System Materials

The NASA Entry Systems Modeling project maintains a portfolio of computational model and tool development activities focused on reducing performance uncertainties in ablative Thermal Protection System (TPS) materials for NASA missions. The development activities span material scale and strive to allow microstructural characterization of material structure and properties, mesoscale analyses of damage, and macroscale evaluation of heatshield performance and recession in a given aerothermodynamic environment. This talk will detail the application of developed capabilities at all three scales to the woven TPS material that the Agency has selected as the heatshield for the Mars Sample Return Earth Entry System (MSR-EES) mission – 3D Mid-Density Carbon Phenolic (3MDCP). Each of the applications focuses on driving down uncertainties in material performance and thus risk for MSR-EES and other future missions that may leverage woven TPS. At the microscale, machine learning techniques are used to characterize images from destructive microscopy and inform structural variability. At the mesoscale, Lagrangian techniques are used to simulate ballistic impact and interpret damage modes noted in experiments. At the macroscale, coupled flow-material response techniques are validated by Arc Jet testing to enable heatshield design for missions with massive ablation.

Justin B Haskins↗
DOE OSTI2022

Improved Understanding of Zircaloy Hydrogen Pickup Mechanisms in BWRs. NSUF Report (Final Report)

Using an instrumented autoclave operating at 288°C and 1100 psi, electrical impedance spectroscopy (EIS) measurements were taken from three neutron irradiated fuel channel samples (Zircaloy-2 and Zircaloy-4). The three high priority fuel channel samples were tested at 288°C using a simulated PWR water with the AC frequency range extending down through 10 -4 Hz, with corresponding DC measurements to validate the resistance measurements. EIS measurements were also taken at room temperature using a simulated PWR water from four neutron irradiated fuel channel samples (Zircaloy-2 and Zircaloy-4). and four water rod samples (Zircaloy-2). EIS measurements were also performed on unirradiated Zircaloy-2 tubes that were previously oxidized in autoclaves at various temperatures for different time durations to obtain oxide layers of different thicknesses. These experiments proved to be more challenging than anticipated because of the small size of the samples, which required development of fixtures to allow valid testing of the materials, some of which are radioactive, under pressure and temperature. Equivalent circuits were used to model the data. Some samples resulted in Nyquist plots where the arcs were too close to each other resulting in the inability to separate the different oxide morphological thicknesses. Microstructural characterization in the form of optical microscopy, scanning electron microscopy and (scanning) transmission electron microscopy was performed to document the microstructural features of the oxidized Zr-2 and Zr-4 samples (oxide and metal/oxide interface) and to interpret the EIS results. At the start of this work, it was assumed that despite the laminar cracks and porosity visible in both optical and SEM, the oxide scale is dense near the metal, effectively preventing the water from contacting the metal directly. While it is generally accepted that the large laminar cracks pores allow the water to percolate into the oxide film, it is unknown just how deep it can go, and how interconnected is the porosity, which will dictate how much the porosity is connected electrically to the water within the oxide scale. The EIS results demonstrated that there is in fact a distribution of oxide thicknesses over the various samples, and within each material condition, that are communicating to the liquid. This is derived from the fact that Nyquist plots have an arc radius whose center point is often below the Z’ axis, implying that instead of a single oxide thickness, these oxides have a distribution of thicknesses. This distribution of thicknesses produces arcs that are not clearly separable in most of the samples.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Beneficial effect of copper on pitting resistance of Ni-Cr-Fe alloys

This study examines the effect of copper alloying on pitting resistance in a model solid solution FCC Ni-13%Cr-10%Fe alloy through potentiodynamic and potentiostatic polarization in 0.1 M NaCl in conjunction with an analysis using first-principles competitive electro-chemisorption modeling. The pitting potential increased with increasing Cu content in the alloy. Furthermore, the extent of metastable pit growth was suppressed and the incubation time for metastable to stable pit transition increased with Cu content. The first-principles competitive adsorption calculations suggested that Cu alloying suppresses chloride ion adsorption on the alloy surface in a simulated pit environment, which inhibits active dissolution at the pit bottom, enhances proton adsorption, and thereby increases the local pH at the pit bottom. In conclusion, we propose that these two effects of Cu in solid solution combine to reduce pit stability and may act in addition to the enrichment of Cu on the corroding pit surface.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Topological Defect Mediated Helical Phase Reorientation by Uniaxial Stress

Strain engineering enables precise, energy-efficient control of nanoscale magnetism. However, unlike well-studied strain-dislocation interactions in mechanical deformation, the spatial evolution of strain-induced spin rearrangement remains poorly understood. Using Lorentz transmission electron microscopy, we manipulate and observe helical domain reorientation under quantitatively applied uniaxial tensile stress. Our findings reveal striking similarity to plastic deformation in metals, where the critical stress for propagation vector ( Q ) reorientation depends on its angle with the stress direction. Magnetic defects mediate reorientation via “break-and-reconnect” or “dislocation gliding–annihilation” processes. Simulations confirm that strain-induced anisotropic Dzyaloshinskii-Moriya interaction may play a key role. These insights advance strain-driven magnetism and offer a promising route for energy-efficient magnetic nanophase control in next-generation information technology. Published by the American Physical Society 2025

36 MATERIALS SCIENCE↗