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 325 records · Page 18

Stress Distribution During Deformation of Polycrystalline Aluminum by Molecular-Dynamics and Finite-Element Modeling

In this paper, a multiscale modelling strategy is used to study the effect of grain-boundary sliding on stress localization in a polycrystalline microstructure with an uneven distribution of grain size. The development of the molecular dynamics (MD) analysis used to interrogate idealized grain microstructures with various types of grain boundaries and the multiscale modelling strategies for modelling large systems of grains is discussed. Both molecular-dynamics and finite-element (FE) simulations for idealized polycrystalline models of identical geometry are presented with the purpose of demonstrating the effectiveness of the adapted finite-element method using cohesive zone models to reproduce grain-boundary sliding and its effect on the stress distribution in a polycrystalline metal. The yield properties of the grain-boundary interface, used in the FE simulations, are extracted from a MD simulation on a bicrystal. The models allow for the study of the load transfer between adjacent grains of very different size through grain-boundary sliding during deformation. A large-scale FE simulation of 100 grains of a typical microstructure is then presented to reveal that the stress distribution due to grain-boundary sliding during uniform tensile strain can lead to stress localization of two to three times the background stress, thus suggesting a significant effect on the failure properties of the metal.

Yamakov, V.↗

Investigating Sulfur-Rich Mercury Analogs Exposed to Simulated Micrometeoroid Bombardment in the Laboratory

Space weathering (SW) continually alters the spectral, microstructural, and chemical characteristics of the surface of airless bodies across the solar system. The effects of SW vary depending on the heliocentric distance and the initial composition of the target surface. While SW on the Moon and S-type asteroids is well documented, our understanding of how this process affects Mercury is at an early stage. Mercury’s interplanetary environment is harsh, with the surface of the planet experiencing an intense solar wind flux as well as a higher flux and velocity of micrometeoroid impactors compared to the Moon and S-type asteroids. In addition, Mercury is also a geochemical endmember, with a surface composition low in Fe (<2 wt.%) and enriched in volatile components, such as sulfur (up to 4 wt.% in the low reflectance material (LRM)). These volatile components are thought to play a major role in the formation of hollows via their sublimation. Sulfur has been hypothesized to occur at the surface of Mercury as sulfide minerals (MgS, CaS) based on its correlation with Mg and Ca in remote sensing data. However, recent observations of chaotic terrains in the north polar area of Mercury and of glacier-like features at lower latitudes indicated that octasulfur (S 8 , elemental sulfur) is another likely constituent of a volatile-rich layer in Mercury’s crust. Its behavior on Mercury may result in a complex cycle of enrichment and depletion. While S is often depleted on small body surfaces, Mercury’s gravity could result in ejected S subsequently returning to the surface and coating regolith grains. Further, the reaction of reduced S-rich gas with glasses of a Mercury-like composition also produced S-rich coatings. However, the precise behavior and evolution of S-rich species exposed to the harsh SW on Mercury remains poorly understood and needs to be further investigated in the laboratory. Here, we present the results of our analyses of the spectral, microstructural, and chemical characteristics of S-rich Mercury analogs irradiated by pulsed laser to simulate the short duration, high temperature events associated with micrometeoroid impacts.

N Bott↗

An approach to combine neutron and ion irradiation data to accelerate material qualification for nuclear reactors

Next-generation nuclear power plants are generally characterized by higher operating temperatures, increased neutron fluences and energies, and distinct corrosive coolant environments versus the existing light water reactor fleet. Whether using existing materials in new environments, newly developed materials tailored for these environments, or new manufacturing methods, the traditional decades-long approach for materials qualification does not facilitate rapid deployment. Ion irradiation has demonstrated success in reproducing material microstructure and select property evolution resulting from neutron irradiation with three to four orders of magnitude reduction in time and cost, making it an ideal candidate for accelerated irradiation testing. Because microstructure has a large impact on bulk material properties, limited neutron irradiation data at lower damage levels can in principle be combined with accelerated ion testing results and modeling and simulation to form an accurate prediction of microstructure evolution and select properties under different neutron irradiation conditions and at higher damage levels. The objective of this work is to present a conceptual framework of specific steps to fulfill several technical challenges associated with qualifying materials for performance in radiation environments on an accelerated time frame. A brief review of the regulatory landscape for materials in nuclear environments is presented, followed by additional overviews to understand the current state of the art for correlation of materials properties across radiation environments using experimental and computational methodologies. Finally, the roles of academia, national laboratories, and industry in the advancement of this accelerated materials qualification framework are discussed as a path forward, with possible case studies presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Multi-scale characterization and simulation of impact welding between immiscible Mg/steel alloys

Vaporizing foil actuator spot welding method is used in this paper to join magnesium alloy AZ31 and uncoated high-strength steel DP590, which are typically considered as un-weldable due to their high physical property disparities, low mutual solubility, and the lack of any intermetallic phases. Characterization results from scanning electron microscopy (SEM) and high-resolution transmission electron microscopy (HRTEM) of the weld interface indicate that the impact creates an Mg nanocrystalline interlayer with abundant Fe particles. The interlayer exhibits intact bonding with both DP590 and AZ31 substrates. To investigate the fundamental bond formation mechanisms at the interface, a finite element (FE)-based process simulation is first performed to calculate the local temperature and deformation at the interface under the given macroscopic experimental condition. Finally, taking the FE results at the interface as inputs, molecular dynamics (MD) simulations are conducted to study the interlayer formation at the Mg/Fe interface during the impact and cooling. The results found a high velocity shearing-induced mechanical mixing mechanism that mixes Mg/Fe atoms at the interface and creates the interlayer, leading to the metallurgical bond between Mg/steel alloys.

36 MATERIALS SCIENCE↗

Enhancing fatigue life of aluminum alloy castings through cavitation water jet peening: Experiments and simulations

This study presents an investigation into the enhancement of the fatigue life of aluminum castings through the application of cavitation water-jet peening (CWJP). CWJP harnesses the impacts of water cavitation to induce surface compressive residual stress within metallic materials. In this work, CWJP was applied to a high pressure die-cast (HPDC) Al–Si alloy A380 with three different water-jet traverse velocities. The fatigue-life improvement, evaluated in a 4-point bending configuration (stress ratio R = 0.1), was found to vary with applied stress level and ranges from 1.6 to 10 times that of the parent alloy. The data also shows that decreasing the traverse velocity results in greater compressive residual stresses within the surface layer and a concurrent increase in surface roughness. This residual stress layer extends to a depth of 400 μm below the surface, as confirmed by through-thickness residual stress and microhardness measurements. CWJP treatment effectively slows down fatigue crack propagation, as evidenced by microstructural observations of narrower striation spacing. Simulations reveal that compressive residual stresses, in addition to surface hardening during CWJP, are key to improving fatigue life. A 20% increase in surface hardness and 150 MPa compressive residual stress imposed by CWJP process provides an average 5-fold enhancement of fatigue life across different stress levels. This study demonstrates the potential of CWJP as an effective surface treatment to enhance the fatigue life of aluminum castings, such as HPDC components for automotive applications.

Al casting↗

A Computational Study to Investigate the Effect of Defect Geometries on the Fatigue Crack Driving Forces in Powder-Bed AM Materials

Powder-bed additive manufacturing (AM) processes are associated with the formation of multiple types of process-specific pores, including but not limited to lack-of-fusion (LoF) and keyhole pores. The performance of an AM component is dependent on the type of pores, their density and their proximity to the free surface, and other heterogeneities in the microstructure. In order to characterize the influence of porosity on the mechanical behavior of AM materials, it is imperative to quantitatively analyze the heterogeneous strain accumulation in the vicinity of porosity. Process-specific microstructure models are generated using SPPARKS, an open-source process simulation code. Spherical keyhole or irregular LoF pores are embedded into the microstructure models, which are meshed and input into a finite element code, ScIFEN, to solve for the heterogeneous strain localization in the vicinity of the pores. Given the non-smooth geometries of LoF pores, they readily promote strain accumulation in their vicinity thereby increasing the propensity of initiating fatigue cracks.

Saikumar R Yeratapally↗

Reversible Binding Interfaces Made of Microstructured Polymer Brushes

The polymer brush architecture of the end-tethered polymer molecules is one of the most widely used efficient methods to regulate interfacial interactions in colloidal systems found in live matter and manufactured materials. Emerging applications of polymer brush structures require solutions to new tasks in the control of interfacial interactions. The rapid development of live cell manufacturing relies on scalable and efficient cell harvesting methods. Stimuli-responsive surfaces made of surface-grafted poly(N-isopropylacrylamide) (PNIPAM) can bind and detach the adherent cell upon changes in temperature and have been used for cell growth and harvesting. The applications are limited by the requirement to satisfy a range of PNIPAM coating characteristics that depend on the dimensions of the integrin complex in the cell membrane and the basal surface. The analysis of the microstructured surfaces, when adhesive and disjoining functions of the microdomains are decoupled, shows that many limitations of PNIPAM one-component coatings can be avoided by using a much broader range of structural characteristics of the microstructured interfaces composed of alternating disjoining PNIPAM domains and adhesive polymeric domains with cell-affinity functional groups. Temperature-controlled reversible adhesion to such microstructured interfaces is studied here experimentally with model systems of solid spherical particles and by employing simulations for solid and soft membranes interacting with the microstructured surfaces to mimic interactions with soft and solid disk-like particles.

36 MATERIALS SCIENCE↗

Development of a numerical model for simulating stress corrosion cracking in spent nuclear fuel canisters

Abstract Prediction and detection of the chloride-induced stress corrosion cracking (CISCC) in Type 304 stainless steel spent nuclear fuel canisters are vital for the lifetime extension of dry storage canisters. This paper conducts a critical review that focuses on the numerical modeling and simulation on the research progress of the CISCC. The numerical models emphasizing the residual stress, susceptible microstructure, and corrosive environment are summarized individually. Meanwhile, the simulation studies on the role of hydrogen-assisted cracking are reviewed. Finally, a multi-physical numerical model, which combines the different fields is proposed based on our recent investigation.

Materials Science↗

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗

Modeling the 4D discharge of lithium-ion batteries with a multiscale time-dependent deep learning framework

The lithium-ion battery (LIB) field is moving towards the direction of investigating spatially resolved physical phenomena in the 3D porous microstructure of electrodes. These pore-scale simulations give new insights into the local dynamics of lithiation/de-lithiation and charge transport, Nevertheless, the computational time of these simulations limits the integration of these models in optimization workflows of cycling conditions or electrode manufacturing processes. Machine learning models present a way of assessing in real-time the performance of materials. While several successful techniques for replicating simulations with machine learning have been proposed, this case study presents a more demanding problem, due to the necessity of understanding the behavior of heterogeneous 3D local data, as it evolves in time: this poses both a scientific and a technical challenge. To this end, we propose an autoregressive multiscale convolutional neural network model to predict relevant quantities at the pore-scale in the solid phase: the lithium concentration (in the active material) and potential (in the active material and carbon binder). Here, these are ultimately used to reconstruct the battery discharge curve. 3D images of the electrode microstructures are the input to the network, trained with a dataset of finite element method simulations to predict the discharge behavior of the cathode side in lithium ion batteries. We propose this machine learning model as a proof-of-concept of the applicability of multiscale networks for time-dependent physics problems. The trained model exhibits very high accuracy (with errors lower than 2 %) in forecasting the discharge behavior of new unseen cathodes.

25 ENERGY STORAGE↗

Using simulation to study solid oxide cell degradation

Slides of the presentation given by Harry Abernathy from the National Energy Technology Laboratory at the 2023 International Conference and Expo on Advanced Ceramics and Composites. The presentation was an overview of NETL's solid oxide cell performance degradation modeling, including a demonstration of the impact of microstructural analysis methods on the resulting degradation simulations.

Abernathy, Harry↗

Integrated Process-Structure-Property Simulations for Additive Manufacturing Using the Open-Source Materialite Package

The microstructure and properties of additively manufactured (AM) metals are strongly dependent on process conditions. Therefore, process-structure-property (PSP) simulations are a useful tool for exploring process parameter space, studying process variations, and quantifying uncertainty in material properties. However, integrating process-structure and structure-property simulations often involves connecting multiple software packages. Each package may use unique data structures and require substantial domain knowledge. This presentation demonstrates PSP simulation capabilities of Materialite, an open-source package developed at NASA Langley Research Center. Materialite simplifies model linkages by using a common data structure and model interface, enabling straightforward simulation across a PSP model chain. Physics-based models, including kinetic Monte Carlo and crystal plasticity, are implemented within the package. The model interface is also intended to simplify implementation of new models and enable integration with external simulation tools. Example use cases include uncertainty quantification with PSP models and GPU-accelerated powder bed fusion AM process models.

additive manufacturing↗

Fluid Dynamics Effects on Microstructure Prediction in Single-Laser Tracks for Additive Manufacturing of IN625

Single-track laser fusion were simulated using a heat-transfer-solidification-only (HTS) model and its extension with fluid dynamics (HTS_FD) model using a parallel open-source code, which included laminar fluid dynamics, flat-free surface of the molten alloy, heat transfer, phase-change, evaporation, and surface tension phenomena. The results illustrate that the fluid dynamics affects the solidification and ensuing microstructure. For the HTS_FD simulations, thermal gradient, G was found to exhibit a maximum at the extremity of the solidified pool ( i.e. , at the free surface), while for HTS simulations, G exhibited a maximum around the entire edge of the solidified pool. HTS_FD simulations predicted a wider range of cooling rates than the HTS simulations, exhibited an increased spread in the solidification speed, V variation within the melt-pool with respect to the HTS model results. Primary dendrite arm spacing (PDAS) were evaluated based on power law correlations and marginal stability theory models using the ( G , V ) from HTS and HTS_FD simulations to quantify the effect of the fluid dynamics on the microstructure. At low-laser powers and low-scan speeds, the PDAS obtained with the fluid dynamics model (HTS_FD) was larger by more than 30 pct with respect to the PDAS calculated with the simple HTS model. A new PDAS correlation, i.e. , \( \lambda_{1} \left[ {\mu {\text{m}}} \right] = 832\;G\left[ {\text{K/m}} \right]^{ - 0.5} V\left[ {\text{m/s}} \right]^{ - 0.25} \) , which uses the ( G , V ) results from the HTS_FD model was developed and validated against experimental results.

36 MATERIALS SCIENCE↗

An integrated modeling framework with open architecture for phase field simulation of multi-component alloys

An integrated modeling framework (PanPhaseField) has been developed, which enables a direct and fast coupling between CALPHAD calculations and large-scale phase field simulations for multi-component alloys. Further, it adopts an open architecture allowing for integration of user-defined phase field models in a plug-and-play manner by taking full advantage of the user-friendly graphical interface of Pandat software. The developed modeling platform becomes an enabling tool that can be used to simulate the evolution of spatially varying microstructures of industrial complex alloys for various engineering applications.

36 MATERIALS SCIENCE↗

Phase evolution and mechanical performance of an ettringite-based binder during hydrothermal aging

Little is known about the performance of ettringite-based binders in hydrothermal conditions. This investigation aims to gain insights into the phase evolution and corresponding mechanical performance of an ettringite-based binder considering crystallization pressure caused by late-reaction products. Additionally, the role of fiber reinforcement on the strength retention of the binder was investigated. When aged at an elevated temperature under water-saturated conditions, hard-burned MgO hydrated to form brucite. The precipitation and growth of the brucite crystals led to a crystallization pressure of approximately 200 MPa calculated using thermodynamic modelling. Damage was observed after 4 months of aging with cracks in the microstructure and eventually a failure at the macro scale. Ettringite remained stable at 60 °C due to the water-saturated conditions. Polypropylene fiber delayed crack propagation and thus reduced the damage caused by crystallization pressure. The fiber improved the flexural performance of composite attaining deflection-hardening behavior regardless of aging conditions.

36 MATERIALS SCIENCE↗