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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effect of void positioning on the detonation sensitivity of a heterogeneous energetic material

We show although it is well-established that voids profoundly influence the initiation and reaction behaviors of heterogeneous energetic materials such as polymer-bonded explosives (PBX) and propellants, there has been little study of how void location in different constituents in the microstructures of such materials affect the macroscale behavior. Here, we use three-dimensional (3D) mesoscale simulations to study how void placement within the reactive grains versus the polymer binder influences the shock-to-detonation transition (SDT) in a polymer-bonded explosive. The material studied here has a microstructure comprised of 75% PETN (pentaerythritol tetranitrate) grains and 25% HTPB (hydroxyl-terminated polybutadiene) polymer binder by volume. Porosities up to 10% in the form of spherical voids distributed in both the grains and polymer are considered. An Arrhenius reactive burn relation is used to model the chemical kinetics of the PETN grains under shock loading, thereby resolving the heterogeneous detonation behavior of the PBX. The influence of void location on the shock initiation sensitivity of the material is quantitatively ranked by comparing the predicted run distance to detonation (RDD) for each sample. The analysis includes inherent quantification of uncertainties arising from the stochastic variations in the microstructure morphologies and void distributions by using statistically equivalent microstructure sample sets (SEMSS), leading to probabilistic formulations for the RDD as a function of shock pressure. The calculations reveal that the location of voids in the composite microstructure significantly affects the RDD. Specifically, voids exclusively within the grains cause the PBX to be more sensitive (having shorter RDD) than voids in the polymer binder. Unique probabilistic relationships are derived to map the probability of observing RDD for each void location material case, allowing for prediction of initiation behavior anywhere in the shock pressure – RDD space. These findings agree with trends reported in the literature.

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

Microstructural dependence of defect formation in iron-oxide thin films

In this report passivating iron-oxide films are grown atop iron films simulating the corrosion process in a nuclear reactor environment. Two oxide films grown via physical vapor deposition at 600 °C and room temperature exhibited dense-epitactic and columnar-polycrystalline, microstructures respectively. A third oxide film grown in open air at 600 °C exhibited an eqiuaxed, porous morphology. Cubic maghemite and magnetite phases in each oxide film were identified via grazing incidence X-ray diffraction. Positron annihilation spectroscopy was used to characterize point defects and measure their depth and size distributions in each oxide layer and showed a range of average positron lifetimes from 0.23 ns in the high temperature, vapor deposited film, 0.35 ns in the room temperature-grown film, and 0.31 ns in the thermally grown oxide. These data indicate that the film morphology, which varies greatly in these films, leads to very different defect content. Finally, four-dimensional scanning transmission electron microscopy was used to measure the internal stress of each film and was correlated to the strain state presented in the X-ray diffraction spectra. The defect formation in each film is reasoned through using a thin film growth model.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phase-field simulations of fission gas bubble growth and interconnection in U-(Pu)-Zr nuclear fuel

Abstract The growth and interconnection of fission gas bubbles in the hotter central regions of U-(Pu)-Zr nuclear fuel has been simulated with a phase-field model. The Cahn-Hilliard equation was used to represent the two-phase microstructure, with a single defect species. The volume fraction of the bubble phase and surface area of the bubble-matrix interface were determined during growth and interconnection. Surface area increased rapidly during the initial stages of growth, then slowed and finally decreased as bubble interconnection began and coarsening acted to reduce surface area. The fraction of the bubbles vented to a simulation domain boundary, f V , was quantified as a measure of the microstructure’s interconnectivity and plotted as a function of porosity p . The defect species diffusivity was varied; although changes in diffusivity significantly affected the microstructure, the plots of f V vs. p did not change significantly. The percolation threshold p c was calculated to be approximately 0.26, depending on the assumed diffusivity and using an initial bubble number density based on experimental observations. This is slightly smaller than the percolation threshold for continuum percolation of overlapping 3D spheres. The simulation results were used to parameterize two different engineering-scale swelling models for U-(Pu)-Zr in the nuclear fuel performance code BISON.

Aagesen, Larry K. (ORCID:000000034936676X)↗
DOE OSTI2024

Location-Specific Microstructure Characterization Within AM Bench 2022 Nickel Alloy 718 3D Builds

Abstract The Additive Manufacturing Benchmark Test Series (AM Bench) is a broad effort to produce rigorous measurement datasets for validating AM computer simulations across the range of processing, structure, and properties, for many additive manufacturing (AM) build methods and material classes. Here, the microstructures of nickel alloy 718 AM Bench 2022 test artifacts produced using laser-based powder bed fusion (PBF-LB), in both as-built and fully heat-treated conditions, are examined. Cross sections are primarily characterized using large area scanning electron microscopy (SEM) electron backscatter diffraction (EBSD) and example analyses of the crystallographic textures are described. These data are part of a large set of in situ and ex situ measurements from both three-dimensional builds and laser tracks on bare plates. All the measurement data are available online with download links at www.nist.gov/ambench .

Levine, L. E. (ORCID:0000000334484229)↗
DOE OSTI2025

Optimizing Bendability and Hardness of Age-Hardenable Aluminum Sheets through Local Thermo-Mechanical Processing

This paper introduces a novel thermo-mechanical process to modify the local mechanical properties of 6xxx aluminum (Al-Mg-Si-Cu) alloy sheets. In this process, two pairs of rollers travel along the length of a sheet while locally bending and unbending it so that the final shape and thickness of the sheet remains unchanged. Room temperature bending/unbending (B/U) produces a deformation gradient through thickness and local hardening in both T4 and T6 sheets (with 41% and 18% greater Vickers hardness, respectively). High temperature bending/unbending performed at ~ 500°C via induction heating, produces significant improvements in local formability as measured by bend testing. Formability levels equivalent to the as-received T4 temper are achieved within the processed zones of T6 sheets with bend angles of ~ 150°, without disturbing the T6 temper in the remainder of the sheet. Finite element simulations and characterization explained the mechanical property improvements and the heterogenous microstructure with weakened Cube texture. Due to its similarity with the roller hemming process, and its compact dimensions, the apparatus developed for B/U has the potential for seamless integration with roller hemming robots to selectively process high strength aluminum sheet materials within mass-production settings.

Efe, Mert [BATTELLE (PACIFIC NW LAB)]↗
DOE OSTI2020

The passivity of low-temperature carburized austenitic stainless steel AISI-316L in a simulated boiling-water-reactor environment

The passive film on low-temperature-carburized AISI-316L austenitic stainless steel was studied in an oxygenated, simulated boiling-water-reactor environment with normal water chemistry. The microstructure was studied using scanning electron microscopy, grazing-incidence X-ray diffractometry, Auger electron spectroscopy, X-ray energy-dispersive spectroscopy, and transmission electron microscopy. For non-surface engineered AISI-316L, three distinct layers were observed after exposure: an outer layer of larger nickel-oxide enriched particles with the spinel structure, an intermediate layer of small Cr-oxide-enriched hematite particles, and a compact inner layer of Cr-enriched oxide also having the spinel structure. Exposure of low-temperature-carburized AISI-316L, in contrast, resulted in only two distinct layers. These layers did not include a compact Cr-rich inner layer. Instead, the innermost layer consisted of an oxide with the spinel structure that was a mixture of Fe, Cr and Ni oxide. The outer layer on the low-temperature-carburized specimen consisted of large loosely stacked particles with the spinel structure of mixed composition. In addition to passive-film characterization, the corrosion rate was studied by monitoring mass loss. Finally, after 500 h of immersion, the mass changes measured were relatively low for all specimens indicating good adhesion and low oxidation rates.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Molecular-dynamics study of diffusional creep in uranium mononitride

Uranium mononitride (UN) is a promising advanced nuclear fuel due to its high thermal conductivity and high fissile density. Yet, many aspects of its mechanical behavior and microstructural features are currently unknown. In this paper, molecular dynamics (MD) simulations are used to study UN's diffusional creep. Nanometer-sized polycrystals are used to simulate diffusional creep and to calculate an effective GB width. It is found that Nabarro-Herring creep is not dominant in the temperature range of 1700–2000 K and that the dominant diffusional creep mechanism is Coble creep with an activation energy of 2.28 ± 0.09 eV. A method is proposed to calculate the diffusional GB width and its temperature dependence in polycrystals. The effective GB width of UN is calculated as 2.69 ± 0.08 nm. This value fits very well with the prefactor of the phenomenological Coble creep formula. It is demonstrated that the most comprehensive thermal creep model for UN can be represented as the combination of our Coble creep model and the dislocation creep model proposed by Hayes et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Mechanisms of Asymmetric Membrane Formation in Nonsolvent-Induced Phase Separation

We report the first simulations of nonsolvent-induced phase separation (NIPS) that predict membrane microstructures with graded asymmetric pore size distribution. In NIPS, a polymer solution film is immersed in a nonsolvent bath, enriching the film in nonsolvent, and leading to phase separation that forms a solid polymer-rich membrane matrix and polymer-poor membrane pores. We demonstrate how mass-transfer-induced spinodal decomposition, thermal fluctuations, and glass-transition dynamics implemented with mobility contrast between the polymer-rich and polymer-poor phases are essential to the formation of asymmetric membrane microstructures. Specifically, we show that the competition between the propagation of the phase-separation and glass-transition fronts determines the degree of pore-size asymmetry. In conclusion, we also explore the sensitivity of these microstructures to the initial film composition, and compare their formation in 2D and 3D.

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

Reactive boride infusion stabilizes Ni-rich cathodes for lithium-ion batteries

Engineered polycrystalline electrodes are critical to the cycling stability and safety of lithium-ion batteries, yet it is challenging to construct high-quality coatings at both the primary- and secondary-particle levels. Here, we present a room-temperature synthesis route to achieve full surface coverage of secondary particles and facile infusion into grain boundaries, thus offering a complete “coating-plus-infusion” strategy. Cobalt boride metallic glass is successfully applied to Ni-rich layered cathode LiNi 0.8 Co 0.1 Mn 0.1 O 2 . Here, it dramatically improves the rate capability and cycling stability, including under high-discharge-rate and elevated-temperature conditions and in pouch full cells. The superior performance originates from simultaneous suppression of microstructural degradation of intergranular cracking and side reactions with electrolyte. Atomistic simulations identified the critical role of strong selective interfacial bonding, which offers not only a large chemical driving force to ensure uniform reactive wetting and facile infusion but also lowered the surface/interface oxygen activity, contributing to the exceptional mechanical and electrochemical stabilities of the infused electrode.

25 ENERGY STORAGE↗
DOE OSTI2024

Role of dislocations on martensitic transformation temperatures and microstructure: A molecular dynamics study

Microstructure and defects strongly affect martensitic transformations in metallic alloys. Significant progress has been made in understanding the atomic-level processes that control the role of grain boundaries and precipitates in these solid-to-solid phase transformations. Yet, the role of dislocations and their structures on martensitic transformation temperature and the resulting microstructure remains unclear. Therefore, we used large-scale molecular dynamics simulations to study the forward and reverse transformation of a martensitic material modeled after Ni63Al37 under cyclic thermal loading. The simulations reveal that dislocations in the austenite phase act as one-dimensional seeds for the martensite phase, which is present at temperatures significantly above the martensite start value. We find a reduction in the dislocation density during cyclic thermal loading, which results in the increase in martensite and austenite transition temperatures, in agreement with experiments. Importantly, we extracted a critical martensitic nuclei size for developing stable domains and found that relatively low dislocation densities are needed to grow independent martensitic variants resulting in a multi-domain structure.

Physics↗
DOE OSTI2023

In-Situ Transmission Electron Microscopy Study of the Evolution of Extended Defects in Oxide Nuclear Fuels

In-situ transmission electron microscopy (TEM) with simulated extreme environments is an effective tool for understanding and tracking microstructural changes down to the atomic scale. The objective of this research is to study the effect of temperature on the evolution of extended defects in ThO2 and UO2. Here, we present in-situ TEM isochronal thermal annealing experiments using a micro-electro-mechanical-system (MEMS)-based heating holder. ThO2 and UO2 single crystal specimens were grown inside an inert silver ampoule using hydrothermal synthesis. Both samples were irradiated using 2 MeV protons at 600oC up to 0.1 dpa at Texas A&M University’s Accelerator Laboratory. Fig. 1a shows a weak-beam dark field TEM image of ThO2 after irradiation, indicating the presence of faulted 1/3 <111> type dislocation loops. This presentation will discuss the effect of annealing temperature (600oC, 800oC, 1000oC, 1100oC, etc.) on dislocation loop density, size and distribution, loop nature (interstitial/vacancies) and Burgers vector, as well as the formation of voids in both ThO2 and UO2. This work will also discuss the interaction between defects during annealing. The in-situ TEM annealing cycle is shown in Fig 1b. This research significantly improves the understanding of defect behavior in oxide nuclear fuels with temperature and will aid computational modeling efforts. This work was supported as part of the Center for Thermal Energy Transport under Irradiation (TETI) Energy Frontier Research Center, funded by the U.S. Department of Energy Office of Science.

36 MATERIALS SCIENCE↗
NASA NTRS1981

Further Study of near Solidus Intergranular Cracking in Inconel 718

A series of tests, performed to determine the strain necessary to initiate intergranular cracking in Inconel 718 as a function of temperature, contained enough scatter near the melting temperature that questions remained as to the best curve of curves to fit to the data. Fracture surface analysis showed that the scatter was due to incipient melting in the grain boundary region. The melting contributed to low fracture strain but had only a small on the incipient cracking strain. Gleeble tests, which could be interrupted by water quenching, were used to study the incipient intergranular melting of Inconel 718. This modified weld simulation test provided a sufficiently rapid quench to preserve the intergranular microstructure created during incipient melting. This structure was studied both microscopically and with energy dispensive X-ray analysis. The implications of incipient melting and low-strain incipient cracking on the development of microfissuring envelopes are discussed.

Thompson, R. G.↗
NASA NTRS2016

Multiscale Modeling of Carbon/Phenolic Composite Thermal Protection Materials: Atomistic to Effective Properties

Next generation ablative thermal protection systems are expected to consist of 3D woven composite architectures. It is well known that composites can be tailored to achieve desired mechanical and thermal properties in various directions and thus can be made fit-for-purpose if the proper combination of constituent materials and microstructures can be realized. In the present work, the first, multiscale, atomistically-informed, computational analysis of mechanical and thermal properties of a present day - Carbon/Phenolic composite Thermal Protection System (TPS) material is conducted. Model results are compared to measured in-plane and out-of-plane mechanical and thermal properties to validate the computational approach. Results indicate that given sufficient microstructural fidelity, along with lowerscale, constituent properties derived from molecular dynamics simulations, accurate composite level (effective) thermo-elastic properties can be obtained. This suggests that next generation TPS properties can be accurately estimated via atomistically informed multiscale analysis.

Materials Science↗
DOE OSTI2024

Modeling Ni Coarsening under Humid Atmosphere in the Electrode of Solid Oxide Cells

Ni coarsening is an important degradation mechanism in solid oxide cells. It is reported that Ni coarsening is faster under a humid atmosphere experimentally, but the underline mechanism is not well understood. In this work, Ni coarsening through Ni(OH)x surface diffusion was investigated by a combination of density-functional theory (DFT) and phase-field modeling (PFM). DFT is used to evaluate the surface coverage and diffusivity of Ni(OH)x on Ni (111) surface and the results are used as an input to the PFM to simulate Ni coarsening under humid atmospheres on both reconstructed and synthetic microstructures. It is found that steam partial pressure and microstructure effect alone cannot explain the faster Ni coarsening under a humid atmosphere observed experimentally. Large local overpotential in fuel cell mode is needed to explain the fast Ni coarsening.

Lei, Yinkai↗
DOE OSTI2021

Coupling of Spark Plasma Sintering with Advanced Modeling to Enable Process Scale-Up: LDRD Final Project Poster

Electric-field assisted sintering (EFAS) manufacturing techniques reduce energy usage requirements, fabricate near-net-shape parts, and produce bulk nanostructured materials through rapid heating and cooling rates achievable through Joule heating. Process parameters of pulsed DC voltage and constant pressure are applied through stainless steel rams and graphite tooling to the ceramic yttrium oxide powder to be sintered. Our high-fidelity modeling code application advances manufacturing for extreme environments by enabling process-informed design and improving microstructure consistency through multiphysics multiscale process-structure-property-performance correlation simulations.

36 MATERIALS SCIENCE↗
NASA NTRS2022

Simulation and Scaling of the Turbulent Vertical Heat Transport and Deep-Cycle Turbulence Across the Equatorial Pacific Cold Tongue

Microstructure observations in the Pacific cold tongue reveal that strong turbulence often penetrates into the thermocline producing hundreds of W/m2 of downward heat transport during nighttime and early morning. However, virtually all observations of this deep-cycle turbulence (DCT) are from 0◦N,140◦W. Here, a hierarchy of ocean process simulations including submesoscale-permitting regional models and turbulence-permitting large eddy simulations (LES) embedded in a regional model provide insight into DCT at and beyond 0◦N,140◦W. A regional hindcast quantifies the spatio-temporal variability of subsurface turbulent heat fluxes throughout the cold tongue from 1999-2016. Climatologically, the mean heat flux is strongest (∼100W/m2) within2◦ of the equator. The heat flux is slightly (∼10 W/m2) stronger in the northern than southern hemisphere throughout the cold tongue. The seasonal cycle of the heat flux, which is not everywhere aligned with the surface flux, ranges from 150 W/m2 near the equator to 30 W/m2 and 10 W/m2 at 4◦N and S respectively. Aseasonal variability of turbulent heat fluxes is logarithmically distributed and temporally uncorrelated with surface heat fluxes. The aseasonal variability of turbulence is highlighted in 34-day LES of Boreal autumn at 0◦N and 3◦N,140◦W. Intense DCT occurs frequently above the undercurrent at 0◦N and intermittently at 3◦N. Daily-mean heat fluxes scale with the bulk vertical shear and the wind stress, which together explain ∼90% of the daily variance across both LES. Observational validation of the scaling at 0◦N,140◦W is encouraging, but observations beyond 0◦N,140◦W can facilitate needed refinement of mixing parameterization in ocean models.

Simulation↗
NASA NTRS1994

Direct numerical simulation of sheared turbulent flow

The summer assignment to study sheared turbulent flow was divided into three phases which were: (1) literature survey, (2) computational familiarization, and (3) pilot computational studies. The governing equations of fluid dynamics or Navier-Stokes equations describe the velocity, pressure, and density as functions of position and time. In principle, when combined with conservation equations for mass, energy, and thermodynamic state of the fluid a determinate system could be obtained. In practice the Navier-Stokes equations have not been solved due to the nonlinear nature and complexity of these equations. Consequently, the importance of experiments in gaining insight for understanding the physics of the problem has been an ongoing process. Reasonable computer simulations of the problem have occured as the computational speed and storage of computers has evolved. The importance of the microstructure of the turbulence dictates the need for high resolution grids in extracting solutions which contain the physical mechanisms which are essential to a successful simulation. The recognized breakthrough occurred as a result of the pioneering work of Orzag and Patterson in which the Navier-Stokes equations were solved numerically utilizing a time saving toggling technique between physical and wave space, known as a spectral method. An equally analytically unsolvable problem, containing the same quasi-chaotic nature as turbulence, is known as the three body problem which was studied computationally as a first step this summer. This study was followed by computations of a two dimensional (2D) free shear layer.

Harris, Vascar G.↗
DOE OSTI2023

Modeling and Simulation of Advanced Manufacturing Techniques using MOOSE and MALAMUTE

Advanced manufacturing techniques offer increased geometry complexity, energy and material usage efficiency improvements, and an expanded palette of materials as compared to conventional manufacturing approaches. Advanced-manufacturing-produced parts can experience wide variations in the final microstructure, and these microstructure variations significantly impact the parts’ performance. In this chapter, we present recent code developments within Multiphysics Object-Oriented Simulation Environment (MOOSE) and in the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE). Here we demonstrate applying these modeling and simulation codes to two advanced manufacturing process types: advanced sintering techniques and laser-based additive manufacturing techniques. The multiphysics and multiscale capabilities of these codes enable prediction of the microstructure evolution resulting from variations in the Advanced manufacturing process parameters.

36 MATERIALS SCIENCE↗