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 397 records · Page 22

Using Simulated Micrometeoroid Impacts to Understand the Progressive Space Weathering of the Surface of Mercury

The surfaces of airless bodies such as Mercury are continually modified by space weathering, which is driven by micrometeoroid impacts and solar wind irradiation. Space weathering alters the chemical composition, microstructure, and spectral properties of surface regolith. In lunar and ordinarychondritic style space weathering, these processes affect the reflectance properties by darkening (lowering of reflectance), reddening (increasing reflectance with increasing wavelength), and attenuation of characteristic absorption features. These optical changes are driven by the production of nanophase Febearing particles (npFe). While our understanding of these alteration processes has largely been based on data from the Moon and near-Earth S-type asteroids, the space weathering environment at Mercury is much more extreme. The surface of Mercury experiences a more intense solar wind flux and higher velocity micrometeoroid impacts than its planetary counterparts at 1 AU. Additionally, the composition of Mercury’s surface varies significantly from that of the Moon. Most notably, a very low albedo unit has been identified on Mercury’s surface, known as the low reflectance material (LRM). This unit is enriched with up to 4 wt.% carbon, likely in the form of graphite, over the local mean. In addition, the surface concentration of Fe across Mercury’s surface is low (<2 wt.%) compared to the Moon. Our understanding of how these low-Fe and carbon phases are altered as a result of space weathering processes is limited. Since Fe plays a critical role in the development of space weathering features on other airless surfaces (e.g., npFe), its limited availability on Mercury may strongly affect the space weathering features in surface materials. In order to understand how space weathering affects the chemical, microstructural, and optical properties of the surface of Mercury, we can simulate these processes in the laboratory [7]. Here we used pulsed laser irradiation to simulate the short duration, high temperature events associated with micrometeoroid impacts. We used forsteritic olivine, likely present on the Mercurian surface, with varying FeO contents, each mixed with graphite, in our experiments. We then performed reflectance spectroscopy and electron microscopy to investigate the spectral, chemical, and microstructural changes in these samples.

Thompson, Michelle S.↗

Temperature-dependent mechanical properties and crystal plasticity parameters for additively manufactured Haynes-214 alloy: Experiments and numerical modeling

Our experimental mechanical testing data demonstrated that the additively manufactured (AM) laser powder bed fusion (L-PBF) Haynes-214 alloy exhibits non-linear mechanical properties as the temperature rises from ambient to 870 °C. Crystal plasticity (CP) simulations provide an effective approach to gaining deeper insights into microstructure-property linkages under thermomechanical loading. This method can reduce the need for costly high-temperature mechanical testing while accounting for the effects of crystallographic texture and grain morphology on the mechanical behavior of AM materials. However, calibrating a CP model is time-consuming because individual simulations are computationally expensive and hundreds (or more) of iterations over parameter sets may be required. To address this issue, we have designed a machine learning-differential evolution (ML-DE) CP framework that can accurately interpolate the tensile properties of AM L-PBF Haynes-214 alloy across a wide temperature range from ambient to 870 °C, with minimal reliance on experimental data. The framework uses electron backscatter diffraction (EBSD) measurements to generate statistically equivalent microstructural volume elements to serve as inputs to the CP modeling framework. Stress–strain curves were generated from 1000 CP simulations, which serve as the training data set for the three ML regression algorithms explored: linear, extra-trees, and multi-layer perceptron. These three regression models were independently evaluated to compare their efficiency and identify the most suitable algorithm for the given problem. Results revealed that the extra-trees ML regressor outperforms the other models in both qualitative and quantitative aspects with an R 2 of 0.98. Subsequently, the differential evolution optimization approach is employed to calibrate the ML-based CP material parameters with experimental results obtained at various temperatures. Finally, temperature-dependent CP material parameters are formulated. The effectiveness and efficiency of the designed framework are validated through comparison with experimental results, demonstrating a high degree of agreement. These calibrated parametric constitutive equations enable further use of the CP model to study the deformation behavior of this alloy under a wide range of thermo-mechanical loading conditions.

36 MATERIALS SCIENCE↗

Discrete-Element and Material-Point Method (DEM and MPM) Based Solvers for Sustainable Technologies

We present the use of discrete element method (DEM) and material point method (MPM) in three relevant green technology applications that include biomass feedstock handling, lithium-ion battery manufacturing, and high-pressure reverse osmosis. Our open-source DEM and MPM solvers are developed using performance portable grid and particle management library, AMReX, thus enabling superior performance on NVIDIA and AMD GPUs with > 100 million particles. Our DEM solver resolves the motion of individual particles in a granular system and includes a bonded sphere method for modeling non-spherical particles along with Hertzian and liquid bridge-based contact models. We simulate highly variable biomass feedstock flows in large-scale hoppers for biofuel production and electrode calendering in battery manufacturing using DEM. Our simulations predict flow blockage in large scale biomass hoppers and electrode microstructure variations, thus providing valuable information for biofuel and battery manufacturers, respectively. The second half of the talk will be on MPM and its application towards pore resolved simulations of reverse osmosis membranes under compressive loads. We present a validation study of our MPM simulations with membrane microscopy imaging thus providing useful insights on membrane stability under high pressure conditions. We also present a spectral stability analysis of using linear hat, quadratic and cubic spline basis in MPM indicating regions of numerical stability.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗

Multiscale Modeling of Structurally-Graded Materials Using Discrete Dislocation Plasticity Models and Continuum Crystal Plasticity Models

A multiscale modeling methodology that combines the predictive capability of discrete dislocation plasticity and the computational efficiency of continuum crystal plasticity is developed. Single crystal configurations of different grain sizes modeled with periodic boundary conditions are analyzed using discrete dislocation plasticity (DD) to obtain grain size-dependent stress-strain predictions. These relationships are mapped into crystal plasticity parameters to develop a multiscale DD/CP model for continuum level simulations. A polycrystal model of a structurally-graded microstructure is developed, analyzed and used as a benchmark for comparison between the multiscale DD/CP model and the DD predictions. The multiscale DD/CP model follows the DD predictions closely up to an initial peak stress and then follows a strain hardening path that is parallel but somewhat offset from the DD predictions. The difference is believed to be from a combination of the strain rate in the DD simulation and the inability of the DD/CP model to represent non-monotonic material response.

Saether, Erik↗

Designed Iron Single Atom Catalysts for Highly Efficient Oxygen Reduction Reaction in Alkaline and Acid Media

Abstract Single atom catalysts (SACs) have attracted much attentions due to their advantages of high catalysis efficiency and excellent selectivity. However, for industrial applications, synthesis of SACs in large and practical quantities is very important. The challenge is to develop synthesis methods with controllability and scalability. Herein, a well‐characterized and scalable method is demonstrated to synthesize atomically dispersed iron atoms coordinated with nitrogen on graphene, SAFe @ NG, with high atomic loading (≈4.6 wt%) through a one‐pot pyrolysis process. The method is scalable for the fabrication of Fe SACs with high quantities. The Fe–N–G catalyst exhibits high intrinsic oxygen reduction reaction (ORR) performance, reaching half potential of 0.876 and 0.702 V in alkaline and acidic solutions, respectively, with excellent microstructure stability. Furthermore, the density functional theory (DFT) simulation confirms that the Fe atoms in coordination with four nitrogen atoms, FeN4, in graphene is the active center for the 4‐electron ORR process. This work demonstrates an efficient design pathway for single atom catalysts as highly active and stable electrocatalysts for high‐performance ORR applications.

Zhao, Shiyong↗

An OpenMP GPU-offload implementation of a non-equilibrium solidification cellular automata model for additive manufacturing

Here, in this paper, performance strategies on GPU-based HPC platforms of a cellular automata (CA) simulation code for non-equilibrium solidification, including nucleation, grain growth, solute partitioning and transport for the metal additive manufacturing (AM) process are investigated using OpenMP 4.5. To accurately report the speed-up for multicore CPUs and GPUs, a rigorous performance analysis employed optimizations appropriate for both CPU-only code (baseline) and GPU offload codes for an isothermal test problem. The performance results on Summit at the Oak Ridge Leadership Computing Facility indicate that using a precomputed list of interface cells significantly decreased the wall-clock time on GPUs. The speedup due to GPU acceleration was evaluated for a full Summit node and measured to be 1.8X when comparing a 6 MPI tasks run with 6 GPUs versus 36 MPI tasks on the CPU only. That speed-up was found to be 7.9X when comparing 6 MPI tasks with 6 GPUs versus the 6 MPI tasks running on the CPU only. Performance measurements showed that system total time is almost constant for runs with more than 96 MPI tasks (or GPUs), indicating that the GPU-accelerated code showed an excellent weak scaling performance. Finally, a rapid directional solidification problem was considered to demonstrate the CA code capability on Summit. It was found that a mesh size of at least 0.05 μm is recommended for the AM-like simulations in order to obtain accurate elongated grain microstructure and elongated subgrain features, which are in qualitative good agreement with experimental data. The results presented in this study indicate that the performance strategies on GPU-based HPC platforms for the CA code are appropriate for novel HPC exascale platforms.

36 MATERIALS SCIENCE↗

A parallel and performance portable implementation of a full-field crystal plasticity model

We have developed a parallel implementation of an Elasto-Viscoplastic Fast Fourier Transform-based (EVPFFT) micromechanical solver to enable computationally efficient crystal plasticity modeling for polycrystalline materials. Our primary focus lies in achieving performance portability, allowing a single EVPFFT implementation to run optimally on various homogeneous architectures, including multi-core Central Processing Units (CPUs), as well as on heterogeneous computer architectures comprising multi-core CPUs and Graphics Processing Units (GPUs) from different vendors. To accomplish this goal, we have leveraged MATAR, a C++ software library that simplifies the creation and utilization of multidimensional dense or sparse matrix and array data structures. These data structures are designed to be portable across diverse architectures through the use of Kokkos, a performance-portable library. Additionally, we have employed the Message Passing Interface (MPI) to efficiently distribute the computational workload among processors. The heFFTe (Highly Efficient FFT for Exascale) library is used to facilitate the performance portability of the fast Fourier transforms (FFTs) computation. The computational performance of EVPFFT is evaluated and presented in terms of parallel scalability and simulation runtime on different high-performance computing (HPC) architectures. As a result, the utility of the developed framework to efficiently simulate the micro-mechanical fields in polycrystalline microstructures in engineering applications is discussed.

36 MATERIALS SCIENCE↗

Microstructure and mechanical properties of high Mn-containing ferritic-martensitic alloys exposed to cyclical thermal treatment

Substantial residual tensile stress tends to accumulate in currently available high-Cr ferritic martensitic steels that are subjected to cyclical heat treatment, which leads to premature brittle fracture. By tailoring the alloy composition, this thermal cycling can be exploited to induce a high number density of nanoprecipitates and phase transformations countering residual tensile stresses. In this work, three new variants of ferritic-martensitic steels have been designed with computational thermodynamics to meet the goals of mitigating residual tensile stresses by lowering martensite start temperatures and of enhancing mechanical strength and irradiation sink strength by increasing the number density of nanoprecipitates. Furthermore, ctenast materials were subjected to cyclical heat treatment. The thermally cycled samples were evaluated with mechanical testing and microstructural analysis to identify the optimal composition in which figures of merit include low residual stress and a high density of nanoscale MX (M = metal, X = C/N) precipitates, leading to high yield strength with reasonable ductility. The noticeably higher density of nanoprecipitates in the optimal alloy favor its higher yield strength, which is supported by the microstructure-derived yield strength calculation and precipitation kinetics simulation.

36 MATERIALS SCIENCE↗

Considering interplay between multiple physical phenomena to elucidate single crystal-like texture, phase transformations, and mechanical behavior of directed energy deposited SS316L

A widespread implementation of large scale additive manufacturing (AM) processes, such as wire arc-directed energy deposition (WA-DED) AM can transform the current manufacturing supply chain networks. Naturally, such implementation requires control of the microstructural attributes, such as texture and phase evolution in the processed alloys. Currently, the texture evolution in fusion-based AM (F-BAM) processes is majorly rationalized by the phenomena occurring only during solidification. However, such rationalization is insufficient for understanding the evolution of primary and secondary crystallographic orientations, and consequently, fails to offer a comprehensive understanding and control of overall texture in F-BAM processed alloys. To this end, we report a single crystal (SX)-like texture in WA-DED processed SS316L for the first time. Furthermore, we assess the physical phenomena that may lead to such unique microstructural evolution during WA-DED AM. Subsequently, using microstructural characterization spanning the build height and thermomechanical simulations we investigate the effect of competitive growth and epitaxial growth occurring during solidification and thermally induced plastic deformation occurring post solidification on the overall texture of WA-DED processed SS316L. A spatial variation in solidification pathway is also established and correlated with variation in undercoolings across the build. Tensile tests revealed a strong orientation dependence of deformation mechanisms with over 110% elongation to failure of specimens deformed along <011>. Such anisotropy is rationalized using Schmid's analysis of dislocation slip and deformation twinning. Importantly, overall, the mechanisms outlined in this work will facilitate an enhanced understanding and subsequent control of texture evolution, solidification behavior and mechanical behavior of WA-DED processed steels.

36 MATERIALS SCIENCE↗

Dynamics of Radiation Damage Buildup in Ultrathin Hexagonal Boron Nitride Films under Ion Bombardment

Two-dimensional hexagonal boron nitride (hBN) is attractive for several emerging applications. Ion bombardment can be used to modify the hBN properties. However, the understanding of radiation damage buildup in hBN remains limited. Here, we investigate the effects of the dose rate and ion mass on radiation damage buildup by studying 40 nm-thick hBN films bombarded at room temperature with 500 keV 4 He, 15 N, 40 Ar, and 129 Xe ions and comparing with results for ion bombardment of polycrystalline hBN ceramics. Raman spectroscopy is used to quantify damage buildup, and transmission electron microscopy is used for microstructural analysis. Experiments are complemented by molecular dynamics simulations of the formation and evolution of point defects. Lighter ions are found to be more efficient at disordering hBN than heavier ions. This observation points to a critical role of intracascade defect processes. In contrast, a negligible dose rate effect observed suggests limited intercascade defect dynamic annealing processes for these irradiation conditions. These findings provide a fundamental basis for hBN defect engineering.

2D materials↗

An Eulerian crystal plasticity framework for modeling large anisotropic deformations in energetic materials under shocks

Here, this paper demonstrates a novel Eulerian computational framework for modeling anisotropic elastoplastic deformations of organic crystalline energetic materials (EM) under shocks. While Eulerian formulations are advantageous for handling large deformations, constitutive laws in such formulations have been limited to isotropic elastoplastic models, which may not fully capture the shock response of crystalline EM. The present Eulerian framework for high-strain rates, large deformation material dynamics of EM incorporates anisotropic isochoric elasticity via a hypo-elastic constitutive law and visco-plastic single-crystal models. The calculations are validated against atomistic calculations and experimental data and benchmarked against Lagrangian (finite element) crystal plasticity computations for shock-propagation in a monoclinic organic crystal, octahydro-1,3,5,7-tetranitro-1,3,5,7 tetrazocine (β-HMX). The Cauchy stress components and the resolved shear stresses calculated using the present Eulerian approach are shown to be in good agreement with the Lagrangian computations for different crystal orientations. The Eulerian framework is then used for computations of shock-induced inert void collapse in β-HMX to study the effects of crystal orientations on hotspot formation under different loading intensities. The computations show that the hotspot temperature distributions and the collapse profiles are sensitive to the crystal orientations at lower impact velocities (viz., 500 m/s); when the impact velocity is increased to 1000 m/s, the collapse is predominantly hydrodynamic and the role of anisotropy is modest. The present methodology will be useful to simulate energy localization in shocked porous energetic material microstructures and other situations where large deformations of single and polycrystals govern the thermomechanical response.

42 ENGINEERING↗

Niobium hydride formation in superconducting qubit thin films

The formation of nonsuperconducting hydrides in 160–170-nm-thick films of niobium is examined. We identify six elastically distinct orientation relationships ɛ−Nb 4 ⁢H 3 takes within the matrix, solid-solution 𝛼−NbH 𝑥 phase. We employ a phase field model to assess the impact of elastic energy induced by the strain of phase transformation on the morphology and transformation dynamics of ɛ precipitates within a thin film. We consider the dimensions of the thin film, crystallographic growth direction, and diffusion rates to predict the timescale of hydride evolution. Leveraging the finite element method, we predict two-dimensional and three-dimensional equilibrium shapes of ɛ-hydrides within a bulk sample and in a thin film that has a traction-free surface. Our results suggest that niobium hydrides migrate to the free surface of the film. Precipitates which reach the free surface coarsen, while precipitates within the film dissolve. Precipitates in both two dimensions and three dimensions experience a repulsive interaction force at the free surface, that is attractive in the bulk, shown in experiment and theory of previous studies.

36 - MATERIALS SCIENCE↗

Latent Pitfalls in Microstructure-Based Modeling for Thermally Aged 9Cr-1Mo-V Steel (Grade 91)

A case study was conducted on a mechanistic model development that predicted tensile strength deterioration with thermal aging of 9Cr-1Mo-V steel in supporting the 60-year design life expected for advanced nuclear reactors. For property prediction beyond practical testing times, mechanistic modeling is highly desired, as it taps into the physics of structure–property relationships and therefore can generate reliable results for extrapolation. Meanwhile, as mechanistic models are often complicated, reflecting the intricacy of microstructure and strengthening mechanisms, pitfalls that are difficult to detect often exist. Here, this paper discusses latent pitfalls that are common in mechanistic modeling or specific in this 9Cr-1Mo-V case development through using the American Society of Mechanical Engineers verification and validation in computational solid mechanics (ASME V&V 10) standard for evaluating credibility of modeling in materials engineering. Suggestions are also made for enhancing reliability of microstructure-based modeling.

36 MATERIALS SCIENCE↗

Evaluation of Ceramic Foams and Simulated Reentry Environments

NASA and Ames is conducting ongoing research in lower density reusable TPS material for future reentry vehicles. Ceramic foams from pre ceramic polymer routes have potential for use in both acreage TPS and as tile leading edges for reentry vehicles. One of the key factors to investigate, when developing new materials for reentry applications, is their oxidation behavior in the appropriate reentry environment whrch can be simulated using ground based arc jet (plasma jet) testing. Studies have shown that oxidation rates of materials will differ when exposed to either monatomatic or molecular oxygen and the amount of monoatomic oxygen depends on the conditions particular to a given test situation. Monoatomic oxygen is always present in reentry environment, therefore arc jet testing is required to provide the appropriate conditions (stagnation pressures, heat fluxes, enthalpies, heat loads and atmospheres) encountered during flight. This preliminary work looks at the response of ceramic foams (Si systems) exposed to simulated reentry environments and investigates the influence of microstructure and composition on the materials response.

Stackpoole, Mairead↗

Ceramic Foams for TPS Applications

Ceramic foams have potential in many areas of Thermal Protection Systems (TPS) including acreage and tile leading edges as well as being suitable as a repair approach for re-entry vehicles. NASA Ames is conducting ongoing research in developing lower-density foams from pre-ceramic polymer routes. One of the key factors to investigate, when developing new materials for re-entry applications, is their oxidation behavior in the appropriate re-entry environment which can be simulated using ground based arc jet (plasma jet) testing. Arc jet testing is required to provide the appropriate conditions (stagnation pressures, heat fluxes, enthalpies, heat loads and atmospheres) encountered during flight. This work looks at the response of ceramic foams (Si systems) exposed to simulated reentry environments and investigates the influence of microstructure and composition on the material? response. Other foam properties (mechanical and thermal) will also be presented.

Stockpoole, Mairead↗

NASA's Nuclear Thermal Propulsion Project

Space fission power systems can provide a power rich environment anywhere in the solar system, independent of available sunlight. Space fission propulsion offers the potential for enabling rapid, affordable access to any point in the solar system. One type of space fission propulsion is Nuclear Thermal Propulsion (NTP). NTP systems operate by using a fission reactor to heat hydrogen to very high temperature (>2500 K) and expanding the hot hydrogen through a supersonic nozzle. First generation NTP systems are designed to have an Isp of approximately 900 s. The high Isp of NTP enables rapid crew transfer to destinations such as Mars, and can also help reduce mission cost, improve logistics (fewer launches), and provide other benefits. However, for NTP systems to be utilized they must be affordable and viable to develop. NASA's Advanced Exploration Systems (AES) NTP project is a technology development project that will help assess the affordability and viability of NTP. Early work has included fabrication of representative graphite composite fuel element segments, coating of representative graphite composite fuel element segments, fabrication of representative cermet fuel element segments, and testing of fuel element segments in the Compact Fuel Element Environmental Tester (CFEET). Near-term activities will include testing approximately 16" fuel element segments in the Nuclear Thermal Rocket Element Environmental Simulator (NTREES), and ongoing research into improving fuel microstructure and coatings. In addition to recapturing fuels technology, affordable development, qualification, and utilization strategies must be devised. Options such as using low-enriched uranium (LEU) instead of highly-enriched uranium (HEU) are being assessed, although that option requires development of a key technology before it can be applied to NTP in the thrust range of interest. Ground test facilities will be required, especially if NTP is to be used in conjunction with high value or crewed missions. There are potential options for either modifying existing facilities or constructing new ground test facilities. At least three potential options exist for reducing (or eliminating) the release of radioactivity into the environment during ground testing. These include fully containing the NTP exhaust during the ground test, scrubbing the exhaust, or utilizing an existing borehole at the Nevada National Security Site (NNSS) to filter the exhaust. Finally, the project is considering the potential for an early flight demonstration of an engine very similar to one that could be used to support human Mars or other ambitious missions. The flight demonstration could be an important step towards the eventual utilization of NTP.

Houts, Michael↗

High-Throughput Ion Irradiation and Microstructural Characterization of Multi-Principal-Element Alloys

Multi-principal-element alloys (MPEAs) have garnered interest in the nuclear community for their apparent resistance to microstructural damage under irradiation. While both experimentation and simulation have demonstrated reduced hardening and void swelling in MPEAs, the relationship between composition and irradiation response is still poorly understood since a rich experimental database for MPEAs does not yet exist. In this work, to accelerate the generation of experimental irradiation-response data for MPEAs, high-throughput synthesis, irradiation, and characterization techniques have been employed. Specifically, additive manufacturing has been used to produce arrays of various MPEAs (1-cm2 coupons), which have been subsequently irradiated using 4-MeV Ni2+ ions to a peak damage of 200 dpa at 500 °C. Hardening and void swelling have been assessed using nanoindentation and trenching via plasma focused ion beam (PFIB), respectively. Trends in the experimental data relating hardening and void swelling behavior to composition have been analyzed using different machine learning approaches.

36 MATERIALS SCIENCE↗

Helium causing disappearance of a/2<111> dislocation loops in binary Fe-Cr ferritic alloys

In this work, single and dual-beam self-ion irradiations were performed at 500°C on ultra-high purity Fe14%Cr alloy to ~0.33 displacements-per-atom (dpa), and 0 or 3030 atomic-parts-per-million (appm) helium/dpa, respectively. Using transmission electron microscopy, we reveal that helium can drastically modify the dislocation loop Burgers vector in Fe-Cr alloys. Helium co-implantation caused complete disappearance of a/2<111> type dislocation loops, and the microstructure consisted of only a<100> loops. Conversely, a/2<111> type loops were predominant without He co-implantation. The total loop density remained largely unaffected. The results strikingly contrast literature asserting that helium stabilizes a/2<111> type loops in bcc Fe alloys, based on low temperature irradiations. Collectively analyzing the results with literature suggest that the small positive interaction between helium and self-interstitial atoms (SIA) in Fe predicted by atomistic simulations maybe insufficient to holistically explain the dislocation loop microstructure development in presence of helium. Helium-SIA positive binding inadvertently implies elevated a/2<111> loop fraction and higher loop densities that the present results contradict. Helium induced high cavity density causing a preferential loss of highly glissile <111> clusters, leaving the matrix saturated with <100> type clusters is proposed as a potential mechanism. Further, the in-situ irradiations combined with Burgers vector analysis strengthened the evidence of Cr-induced dislocation loop mobility reduction that appears to stabilize the a/2<111> type loops and causes higher loop densities in Fe-Cr alloys.

36 MATERIALS SCIENCE↗