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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Stress Relaxation Cracking in 347H Stainless Steel Arc Welds: Susceptibility Evaluation of Heat-Affected Zone

Stress relaxation cracking (SRC) is considered one of the major failure mechanisms for 347H stainless steel welds at elevated service temperatures or during post weld heat treatment (PWHT), especially within the heat-affected zone (HAZ). This work focuses on the characterization of SRC susceptibility within 347H physically simulated arc welded HAZ at elevated temperatures. A four-step SRC thermomechanical test in combination with finite element modeling (FEM) of the welding and testing processes is developed to establish a susceptibility map for HAZ. The test first runs a thermal cycle with three different peak temperatures (1335, 1275, and 1150 °C) to duplicate representative HAZ subzone microstructures, followed by time-to-failure examination under a variety of pre-stress (260–600 MPa) and pre-strain conditions (0.03–0.19) as a function of reheat temperatures between 750 and 1050 °C. With the aid of FEM, SRC susceptibility maps are generated to identify the threshold stress, plastic strain, and creep strain as a function of test temperature. It was found out that HAZ subzone with a lower peak temperature (1150 °C) appears to be slightly less susceptible to SRC than the other two subzones that experienced higher peak temperatures. Generally, time-to-fracture reduces with increasing initially applied stress and strain for all test temperatures. The pre-stress thresholds decrease from about 500 to 330 MPa as the testing temperature increases from 800 to 1050 °C, while the corresponding initial plastic strain thresholds reduces from 0.15 to 0.06. The SRC susceptibility was also evaluated through the Larson–Miller Parameter (LMP) analysis as a function of plastic strain, initial stress and starting stress upon reaching the testing temperature, respectively. The 1050 °C test with a high pre-applied strain (0.1) exhibits an extremely short time to failure (t = 3 s) that lies outside the general trend in LMP analysis. Additionally, it was identified that a plastic strain above 0.07 is identified to significantly reduce the bulk creep strain tolerance to fracture and therefore increases SRC susceptibility. Hardness measurement and fractography analysis indicated that the strain aging of niobium carbonitrides and other potential phases in conjunction with intergranular precipitates contributes to an increase in microhardness and increased intergranular cracking susceptibility.

14 SOLAR ENERGY↗
DOE OSTI2022

Morphological-informed Thermal Property Prediction in the Engineering Domain

Using first-principle, atomic scale methods to predict thermal carrier (electron and phonon) behavior at nanoscale yields highly resolved thermal properties. Density functional theory simulations can model myriad effects on the transport of phonons and electrons (e.g., phonon and electron coupling, electron correlation, defect presence, carrier scattering), and yield thermal properties for small geometric domains such as thermal and electrical conductivity, or heat capacity. While these thermal properties are accurate in atomistic and nanoscopic systems, they must be scaled up to the microstructural domain to be a useful predictor for experimental basis, as the microscale is where myriad changes and physical phenomenon occur within a material (e.g., grain boundaries, precipitate aggregations, interfaces, defect clusters). My talk will discuss the necessity of thermal property predictions for nuclear fuels applications and more. Additionally I will discuss the necessity in developing multi-physics, multi-scale methods for future applications in predictions for not only nuclear fuels, but other materials applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phase-field modeling of radiation-induced segregation for multicomponent alloys

Structural alloys under irradiation are known to undergo radiation-induced solute redistribution (RIS) at grain boundaries, leading to detrimental effects such as intergranular corrosion and stress-assisted cracking. To better understand the phenomenon, improved models of RIS applicable to concentrated, multicomponent alloys, and mesoscale microstructures are needed. In this talk, we present a novel grand-potential-based phase-field model to account for the complete set of multicomponent kinetic and thermodynamic couplings between atoms and point defects in the Onsager transport equations. We demonstrate multiscale modeling capability by deriving the Onsager coefficient matrix from atomistic-based Kinetic Monte Carlo simulations. Model predictions and validations of RIS and the effect of defect production, grain boundary sink strength and density are demonstrated for a model FCC FeCrNi system. We also demonstrate the novel capability to describe RIS in the presence of equilibrium segregation described using a density-based CALPHAD thermodynamic framework. This modeling approach overcomes certain limitations in conventional RIS models and provides a step closer towards multiscale modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Co-Extrusion of Dissimilar Aluminum Alloys via Shear-Assisted Processing and Extrusion

Bimetallic tubes are used when a component requires more than one performance requirement, for instance, strength/creep resistance and oxidation/corrosion resistance. Shear-assisted processing and extrusion (ShAPE) can be used to fabricate extrudates that are comparable or superior in performance relative to conventionally manufactured extrudates. For the first time, ShAPE has been successfully employed in producing bimetallic Al tubing consisting of 6061 and 7075 Al alloys. Both light and electron microscopy techniques were used to investigate the integrity of the tubes, especially the interface between the core and cladding of the bimetallic tubes. Void-free bimetallic tubes were produced using ShAPE. Quantification of tube integrity was carried out with tensile testing under as-extruded and post-aged conditions. All the bimetallic tubes in the as-extruded samples exhibited uniform elongation above 5% with good tensile strength. Key insights such as material flow during bimetallic tube extrusion were obtained from the characterization of remnant billet and simulation results from a smoothed particle hydrodynamics model.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Supersonic hot jet ablative testing and analysis of boron nitride nanotube hybrid composites

Boron nitride nanotubes (BNNTs) are high-strength, high-modulus nanotubes with high thermal and oxidative stabilities. Two hybrid composites were prepared with satin weave carbon fiber (CF) and resole-type phenolic resin: one with surface layers of BNNTs and one with alternating interlayers of BNNTs. The samples were subjected to hot jet tests that simulate realistic high-pressure-temperature conditions to understand the behavior of BNNTs under high-pressure erosion. Adding BNNTs to CF/phenolic laminates enhanced the ablation resistance by reinforcing the char material and mitigated localized thermal damage. Hybrid laminates exhibited up to 14% lower weight loss, 55% increase in flexural modulus, higher thermal diffusivity, and improved char yield and microstructure compared to CF/phenolic samples. The surface layer hybrid had many surviving nanotubes reinforcing the char and crystalline oxide structures that could mitigate further oxygen diffusion. Further, various characterization methods were used to deduce possible mechanisms and their products, indicating that BNNTs could serve as growth templates for direct crystalline boron oxide formation. Overall, hybrid BNNT/CF/phenolic laminates displayed better ablation resistance and favorable microstructure evolution under high-pressure conditions.

42 ENGINEERING↗
DOE OSTI2020

Beyond Idealized Models of Nanoscale Metal Hydrides for Hydrogen Storage

Metal hydrides are attractive for compact, low-pressure hydrogen storage, yet a foundational understanding of factors governing their thermodynamics and kinetics is still lacking. Predictive modeling from the atomic to the microstructural scale plays a critical role in addressing these gaps, particularly for nanoscale materials, which promise improved performance but are difficult to probe. In this paper, we summarize strategies being developed within the Hydrogen Materials—Advanced Research Consortium (HyMARC) for going beyond conventional models to incorporate more complex physics, more realistic structures, and better approximation of operation conditions in simulations of nanoscale metal hydrides. We highlight four beyond-ideal factors that influence predicted performance: (1) surface anharmonic dynamics, (2) interface and surface energy penalties, (3) mechanical stress under confinement, and (4) the presence of native surface oxide. Approaches for addressing these factors are demonstrated on model materials representative of high-capacity hydrogen storage systems, and implications for understanding performance under operating conditions are discussed.

08 HYDROGEN↗
DOE OSTI2020

The Sealing Mechanisms of a Fracture in Opalinus Clay as Revealed by in situ Synchrotron X-Ray Micro-Tomography

The detailed mechanisms of the sealing of a single fracture, from hydration to almost complete closure by increase of confining pressure, as monitored from in situ synchrotron X-ray microtomography during the flow of carbonated water, is here shown for the first time. Different mechanisms play the key role at different stages in the evolution of the fracture. Hydration mechanically weakens the surfaces of the fracture and induces a first closure due to microcracking at the asperity contacts, increasing their size and creating choke points. Increase in confining stress promptly hydraulically seals the fracture by closing the main choke point, with a relative small deformation of the sample. Finally, the more pervasive mechanical deformation observed at higher stresses almost completely seals the whole fracture. The evolution of the sample has been described and quantified using 4D image processing, focusing on the evolution of aperture and digital volume correlation. Hydraulic properties of the sample at different stages have been modeled via Stokes flow simulation, and the results compared to the morphometric analysis, finding positive correlations with the average fracture aperture variation along the flowpath in function of confining pressure. Opalinus Clay is found to be a rock markedly prone to sealing in case of flow with carbonated water; this behavior is the result of the large fraction of clays and of its microstructure, lacking both cementing phases and large stiff particles. CO2 in this sample has no evident role in the evolution of the fracture; chemically induced weathering on the surface has not been detected, in contrast with the behavior observed in samples with carbonates as cementing phase.

58 GEOSCIENCES↗
DOE OSTI2025

Quantifying Impacts of Biomass Pelletization on Fast Pyrolysis Using a Single-Particle Reactor, X-ray Computed Tomography, and Computational Modeling

The pore structure and density of lignocellulosic feedstocks dictate intraparticle transport phenomena and thereby play an important role in thermochemical conversion processes such as fast pyrolysis for biofuel and biochemical production. Variations in microstructure are inherent from different biomass species and can be introduced by preprocessing techniques such as cutting and pelletization. Morphological changes also occur during conversion and lead to vastly different pore structures and behavior during pyrolysis, which impact required conversion times and product distributions. The current work presents a comprehensive comparison of fast pyrolysis of neat and pelletized pine feedstocks, which includes single-particle experiments, modeling, and 3D imaging by X-ray computed tomography (XCT). The particle-scale model included anisotropic heat and mass transport in a shrinking particle with pyrolysis reactions based on the CRECK mechanism with boundary conditions informed by reactor-scale simulations of the single-particle reactor. The models were validated by measurements of the temperature and mass loss from single-particle pyrolysis experiments of neat and pelletized pine. Quantitative analysis of XCT geometries revealed that pyrolytic conversion yielded chars with increased porosity and permeability compared to the unpyrolyzed materials, along with decreased tortuosity and anisotropy. Pelletization of the pine feedstock resulted in a much denser, less permeable material, which converted slower and produced more residual char after pyrolysis compared to neat pine. The results from particle modeling revealed that accounting for the dynamic and anisotropic heat and mass transport caused by differences in pore structure is critical to achieving agreement with experimental results. Overall, this study highlights the dramatic differences in conversion behavior imparted by pelletization and the importance of capturing microstructural attributes in computational models to guide the design and optimization of pyrolysis processes for specific biomass feedstocks.

09 BIOMASS FUELS↗
DOE OSTI2021

Interfacial characteristics of P91 steel - Inconel 740H bimetallic structure fabricated using wire-arc additive manufacturing

Although bimetallic structures can offer unique solutions to engineering problems with varying functionality, joining two dissimilar alloys using additive manufacturing involves challenges such as cracking and segregation. In this work, bimetallic structures of Inconel 740H superalloy and P91 steel are processed using wire-arc additive manufacturing (WAAM) for the first time. P91 steel is directly deposited over 740H with a cooling time of 20 minutes in between the two depositions. During microstructure characterization, a large gradient zone with coarse grains spanning for ~2200 μm is observed. Moreover, the microhardness of the gradient zone is the least in comparison with the P91 and 740H deposits. No deleterious phase formation is found; however, intergranular cracks are observed along the gradient zone. In order to deduce the reason for cracking, non-equilibrium solidification simulations are performed for the mixed compositions of P91 and 740H in the gradient zone. Here, the reasons for crack formation in the gradient zone are the sudden change in the volumetric coefficient of thermal expansion due to the formation of MC carbides as well as the development of local strains from the thermal residual stresses. Therefore, altering the deposition sequence and introducing gradient layers with mixed composition between the P91 and 740H deposits will be beneficial in achieving defect-free builds. Through the integration of experiments and simulation, this work gains valuable insights into the processability of dissimilar alloy manufacturing using the WAAM technique.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

Tensile deformation behavior of twist grain boundaries in CoCrFeMnNi high entropy alloy bicrystals

High entropy alloys (HEA) are a class of materials that consist of multiple elemental species in similar concentrations. The use of elements in far from dilute concentrations introduces a multi-dimensional composition design space by which the properties of metallic systems can be tailored. While the mechanical behavior of HEAs has been the subject of active research recently, the role of grain boundaries (GBs) in their deformation behavior remains poorly understood. Motivated by recent experiments on HEAs demonstrating that GBs act as nucleation sites for deformation twins, herein, we leverage atomistic simulations to construct a series of equiatomic CoCrFeMnNi HEA bicrystals with $\langle 110 \rangle$ and $\langle 111 \rangle$ symmetric twist GBs and examine their tensile behavior and underlying deformation mechanisms at 77 K. Simulation results reveal that plastic deformation proceeds by the nucleation of partial dislocations from GBs, which then grow with further loading by bowing into the bulk crystals leaving behind stacking faults. Variations in the nucleation stress exist as function of GB character, defined in this work by the twist angle. Our results provide future avenues to explore GBs as a microstructure design tool to develop HEAs with tailored properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Transforming an Ionic Conductor into an Electronic Conductor via Crystallization: In Situ Evolution of Transference Numbers and Structure in (La,Sr)(Ga,Fe)O 3–x Perovskite Thin Films

Mixed-conducting perovskites are workhorse electrochemically active materials, but typical high-temperature processing compromises their catalytic activity and chemo-mechanical integrity. Low-temperature pulsed laser deposition of amorphous films plus mild thermal annealing is an emerging route to form homogeneous mixed conductors with exceptional catalytic activity, but little is known about the evolution of the oxide-ion transport and transference numbers during crystallization. Here the coupled evolution of ionic and electronic transport behavior and structure in room-temperature-grown amorphous (La,Sr)(Ga,Fe)O 3-x films as they crystallize is explored. In situ ac-impedance spectroscopy with and without blocking electrodes, simultaneous capturing-synchrotron-grazing-incidence X-ray diffraction, dc polarization, transmission electron microscopy, and molecular dynamics simulations to evaluate isothermal and non-isothermal crystallization effects and the role of grain boundaries on transference numbers is combined. Ionic conductivity increases by ≈2 orders of magnitude during crystallization, with even larger increases in electronic conductivity. Consequently, as crystallinity increases, LSGF transitions from a predominantly ionic conductor to a predominantly electronic conductor. The roles of evolving lattice structural order, microstructure, and defect chemistry are examined. Grain boundaries appear relatively nonblocking electronically but significantly blocking ionically. The results demonstrate that ionic transference numbers can be tailored over a wide range by tuning crystallinity and microstructure without having to change the cation composition.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Intermetallics Based on Sodium Chalcogenides Promote Stable Electrodeposition–Electrodissolution of Sodium Metal Anodes

Sodiophilic micro-composite films of sodium-chalcogenide intermetallics (Na 2 Te and Na 2 S) and Cu particles are fabricated onto commercial copper foam current collectors (Na 2 Te@CF and Na 2 S@CF). For the first time a controllable capacity thermal infusion process is demonstrated. Enhanced wetting by the metal electrodeposition leads to state-of-the-art electrochemical performance. For example, Na 2 Te@CF-based half-cells demonstrate stable cycling at 6 mA cm -2 and 6 mAh cm -2 , corresponding to 54 µm of Na electrodeposited/electrodissolved by geometric area. Sodium metal batteries with Na 3 V 2 (PO 4 ) 3 cathodes are stable at 30C (7 mA cm -2 ) and for 10 000 cycles at 5C and 10C. Cross-sectional cryogenic focused ion beam (cryo-FIB) microscopy details deposited and remnant dissolved microstructures. Sodium metal electrodeposition onto Na 2 Te@CF is dense, smooth, and free of dendrites or pores. On unmodified copper foam, sodium grows in a filament-like manner, not requiring cycling to achieve this geometry. Substrate–metal interaction critically affects the metal–electrolyte interface, namely the thickness and morphology of the solid electrolyte interphase. Density functional theory and mesoscale simulations provide insight into support-adatom energetics, nucleation response, and early-stage morphological evolution. On Na 2 Te sodium atomic dispersion is thermodynamically more stable than isolated clusters, leading to conformal adatom coverage of the surface.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A microstructural signature of the coesite-quartz transformation: New insights from high-pressure experiments and EBSD

Ultra-high pressure (UHP) metamorphism is difficult to identify in continental crust as few petrological barometers are suitable for dominantly felsic lithologies. In such cases, burial to extreme depths is commonly identified through the preservation of coesite, a high-pressure polymorph of SiO 2 that typically forms at depths exceeding ∼ 100 km (i.e., > 2 GPa pressure). Unfortunately, coesite readily transforms to quartz upon exhumation, meaning that UHP terranes may often be overlooked. While some studies have suggested that quartz may inherit an orientation signature indicative of former coesite, both the specific nature of this signature and the conditions favouring its development remain uncertain. Here, to address this problem, we combine electron backscatter diffraction analysis of natural and experimental samples to explore microstructural evolution across the coesite-quartz phase transformation. We demonstrate that neighbouring domains of quartz commonly feature an 84 ± 4° rotation of [c] axes around the pole of a common {m} plane. This orientation relationship is a product of epitaxy, whereby the {$11\bar{2}2$} Japan twin plane in quartz nucleates on the (010) plane in coesite. In supercell simulations, the nucleation of Japan twins can be explained by the energetically favourable alignment of quartz tetrahedra on parental coesite tetrahedra. Through experiments, we demonstrate that this signature emerges over a broad range of conditions, regardless of the availability of nucleation sites (e.g., grain boundaries) or the density of crystal lattice defects (e.g., dislocations). Overall, our work provides a quantitative and unambiguous tool for identifying UHP terranes from quartz in isolation.

Coesite↗
DOE OSTI2023

Incorporating dynamic recrystallization into a crystal plasticity model for high-temperature deformation of Ti-6Al-4V

During hot deformation of (α + β) titanium alloys, the simultaneous action of strain and temperature in the (α + β) regime facilitates dynamic recovery, dynamic recrystallization (DRX), and phase transformations via non-equilibrium paths. DRX is manifested in the form of fine recrystallized α or β grains. Here, in the present study, a two-phase crystal plasticity finite element framework (CP-DRX) has been developed which incorporates DRX kinetics into the crystal plasticity (CP) model to predict the flow characteristics of Ti-6Al-4V alloys during thermo-mechanical processing. The CP slip system parameters, as well as elastic properties from both α and β phases of Ti, are calibrated for different strain rate conditions. An EBSD-informed two-phase microstructure representation has been utilized in the CP-DRX framework to explore the time-dependent evolution of DRX microstructure and crystal orientation for different strain rate conditions. The proposed CP-DRX can capture the evolution of crystal orientation and plastic flow stress–strain response of polycrystalline Ti-6Al-4V during the deformation process. Furthermore, the proposed model is able to capture the softening behavior, observed in average stress–strain response from experiments performed using a Gleeble thermomechanical simulator and predict the recrystallization texture.

(α+ β)-Titanium (Ti) Alloys↗
DOE OSTI2024

Benchmarking machine learning strategies for phase-field problems

Abstract We present a comprehensive benchmarking framework for evaluating machine-learning approaches applied to phase-field problems. This framework focuses on four key analysis areas crucial for assessing the performance of such approaches in a systematic and structured way. Firstly, interpolation tasks are examined to identify trends in prediction accuracy and accumulation of error over simulation time. Secondly, extrapolation tasks are also evaluated according to the same metrics. Thirdly, the relationship between model performance and data requirements is investigated to understand the impact on predictions and robustness of these approaches. Finally, systematic errors are analyzed to identify specific events or inadvertent rare events triggering high errors. Quantitative metrics evaluating the local and global description of the microstructure evolution, along with other scalar metrics representative of phase-field problems, are used across these four analysis areas. This benchmarking framework provides a path to evaluate the effectiveness and limitations of machine-learning strategies applied to phase-field problems, ultimately facilitating their practical application.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Characterization and Modeling of Deformation Induced Damage in Titanium Alloys

Over the lifetime of their DOE support, and specifically during the past funding cycle, the PIs have primarily focused on ambient temperature studies of the mechanisms of heterogeneous deformation and shear transfer at GBs. The PIs have developed and applied a broad array of experimental approaches complemented by advanced computational methods and simulations to enable new understanding of GB shear transfer, and a new approach for modeling this behavior at the mechanistic level has been established. The PIs have examined how slip behavior differs on the surface and interior of polycrystals using a variety of approaches. Some of the associated data has been used to generate 3-D models of experimentally-characterized microstructures that both support interpretation of deformation mechanisms as well as reveal where modeling needs to be improved to more closely agree with what is experimentally observed. The PIs progress in constitutive modeling of crystal plasticity has introduced explicit GB resistance to shear that depends on the mutual geometry of active slip systems. This past work places the PIs in a good position to take on the challenges associated with mechanistically understanding the interplay between GBS and slip interactions.

42 ENGINEERING↗
DOE OSTI2020

The Role of Local Inhomogeneities on Dendrite Growth in LLZO-Based Solid Electrolytes

The majority of the ceramic solid electrolytes (LLZO, LATP) demonstrate polycrystalline grain/grain-boundary (G/GB) microstructure. Higher lithium (Li) concentration and lower mechanical stiffness result in current focusing at the GBs. Growth of Li dendrites through local inhomogeneities and subsequent short circuit of the cell is a major concern. Recent studies have revealed that bulk Li metal is a viscoplastic material that has low (~0.3 MPa) and high (~1.0 MPa) yield strength during deformation at smaller and larger rates of strain, respectively. It has been argued that during deposition at smaller current densities, due to its lower yield strength, Li metal should demonstrate plastic flow against stiff ceramic electrolytes, and Li dendrites will be prevented from penetrating through solid electrolytes. In this manuscript, a multiscale modeling framework has been developed for predicting properties of GBs and the bulk of ceramic electrolytes using atomistic calculations for input to mesoscale models. Using the parameters obtained from the atomistic simulations, the mesoscale model reveals that, given enough time, even at low charge rates, lithium dendrites can grow through the GBs of LLZO. The present multiscale model results also provide information regarding the dendrite growth velocity through LLZO.

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

Coupled cluster and dislocation dynamics modeling of microstructure evolution in irradiated materials

We develop here a coupled cluster and dislocation dynamics framework to study the microstructure evolution of irradiated materials. The framework not only accounts for the three dimensional diffusion of radiation-generated clusters, but also their interaction with dislocation networks and the resultant climb motion of discrete dislocations within finite crystals. The framework is solved with a superposition solution scheme, and is applied to investigate the evolution of the irradiation-induced dislocation loops in zirconium (Zr), considering the effects of various bias factors including the diffusion anisotropy difference (DAD) of interstitials and interstitial clusters, the dislocation bias of defects to discrete dislocation segments, and the production bias of defects from the radiation cascade. We find that the DAD is the most critical factor influencing the kinetics of the loop evolution in Zr, while the recombination/interaction of mobile defects can induce a strong spatial dependence of the loop evolution together with the DAD. Here, the method is also adopted to study the evolution of interstitial $\langle$a$\rangle$ and vacancy $\langle$c$\rangle$ dislocation loop ensembles consistent with the microstructure observed during irradiation-induced growth of Zr. Our findings not only reveal the spatial dependence of the size and ellipticity of the dislocation loops, but also suggest a limit on the anisotropy factor of interstitials to reproduce the co-growth of $\langle$a$\rangle$ and $\langle$c$\rangle$ loops in zirconium, in good agreement with experimental observations and other simulation results.

Bias factors↗