Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “microstructure simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Embedding Sensors in 3D Printed Metal Structures

The Transformational Challenge Reactor (TCR) program is leveraging recent advances in modeling and simulation, materials, and additive manufacturing (AM) technologies to design a modern nuclear reactor. Some of the main TCR technologies include in situ monitoring and the integration of sensors during the manufacturing of quality-significant nuclear reactor components. This report describes the general procedure and process optimization for embedding sensors within generic stainless steel 316 (SS316) components using laser powder bed fusion (LPBF). A more detailed, quality-significant test plan and supporting procedures are available upon request (ORNL/TM-2021/2127). LPBF involves the use of a scanning laser to selectively melt regions of a powder bed, additively building a part layer by layer. This report describes the LPBF processing technique and discusses the effects of LPBF processing parameters on the success of the sensor embedding process. Experiments used machined cavities in the form of channels in an SS316 base for the sensors to lay in while material is additively built over the top, thereby embedding them in an SS316 matrix. A preliminary investigation involved using empty SS316 sheaths as surrogates to explore the effects of various LPBF processing parameters and the dimensional requirements of the machined channels. Microstructural investigations showed that a smaller channel width/depth combination closer to the sensor’s diameter was best for the embedding process. After the desired parameters were selected, Type-K thermocouples were embedded and evaluated post-embedding using nondestructive thermal testing, as well as destructive sectioning and microscopy. Post-build characterization showed that the thermocouples were well-bonded to the SS316 matrix and were fully functional after embedding. During thermal testing to temperatures up to 500 °C, the embedded thermocouples read consistently with one another and deviated only slightly from the readings of a nonembedded thermocouple located within the furnace. This slight discrepancy was most likely due to differences in the thermal time constants for a nonembedded thermocouple vs. a thermocouple embedded in a solid SS316 block. The results presented in this report will serve as the foundation for future work that will focus on embedding sensors in relevant TCR reactor components and eventually testing those components under neutron irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Development of Novel Ferritic-Martensitic Steels with Superior Creep Properties for Power Plant Applications

A creep resistant martensitic steel, CPJ-7, was developed with an operating temperature approaching 650°C. Subsequently, another creep resistant martensitic steel, JMP, was designed with the potential to operate at, or slightly above, 650°C. This report describes the development of these alloys from early iterations of CPJ steels to the final CPJ-7 formulation as well as the more heat resistant JMP steel formulation. The design originated from computational modeling for phase stability and precipitate strengthening using fifteen constituent elements. Approximately forty heats of CPJ and ten of JMP, each weighing ~7 kg, were vacuum induction melted. A computationally optimized heat treatment schedule was developed to homogenize the ingots prior to hot forging and rolling prior to final normalization and tempering. Overall, wrought and cast versions of CPJ-7 present superior creep properties when compared to wrought and cast versions of COST alloys for steam turbine and wrought and cast versions of P91/92 for boiler applications. For instance, the Larson Miller Parameter curve for CPJ-7 at 650°C almost coincides with that of COST E at 620°C. The prolonged creep life was attributed to slowing down the process of the destabilization of the MX and M 23 C 6 precipitates at 650°C. On the other hand, the cast version of CPJ-7 also revealed superior mechanical performance, well above commercially available cast 9% Cr martensitic steel or derivatives, especially for creep life at 650°C. The casting process employed slow cooling to simulate the conditions of a thick wall full-size steam turbine casing but utilized a separate homogenization step prior to final normalization and tempering. To advance the development of CPJ-7 for commercial applications, a process was used to scale up the production of the alloy using vacuum induction melting (VIM) and electroslag remelting (ESR), which underlined the importance of melt processing control of intentionally designed minor and trace elements in these advanced alloys. Following the work on CPJ-7, the JMP steels were designed with higher Co for increased solid solution strengthening, Si for oxidation resistance and increased W (with low Mo content) for matrix strength and stability as well as solid solution strengthening. The JMP steels showed increases in creep life compared to CPJ-7 between 118 to 150% at 650°C for testing at various stresses between 138 MPa and 207 MPa. On a Larson-Miller plot, the performance of the JMP steels surpasses that of state-of-the-art MARBN and other MARBN-type steels. The influence of various elements within the composition of the alloys on the microstructure and mechanical properties are discussed. This report presents approximately 420,000 h of in-house creep testing, the equivalent of almost 50 years of cumulative creep tests.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Extending the Takagi-Taupin equations for x-ray nanobeam Bragg coherent diffraction

We present an approach for simulating x-ray nanobeam Bragg coherent diffraction patterns based on the Takagi-Taupin equations. Compared to conventional methods, the current approach can be universally applied to any weakly strained system including semi-infinite crystals that diffract dynamically. It addresses issues such as the curved wave front and redivergence of the focused incident beam. In conclusion, we show excellent agreement with experimental data for a strained La 0.7 ⁢Sr 0.3⁢ MnO 3 thin film on a SrTiO 3 substrate and a path to extracting physical information using automatic differentiation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Oil-Pressure Based Apparatus for In-Situ High-Energy Synchrotron X-Ray Diffraction Studies During Biaxial Deformation

Background: Understanding biaxial loading response at the microstructural level is crucial in helping better design sheet manufacturing processes and calibrate/validate material deformation models. Objective: The objective of this work was to develop a low-cost testing apparatus to probe, with sufficient spatial resolution, the micro-mechanical response of a sheet material in-situ under biaxial loading conditions. Methods: The testing apparatus fabricated as a part of this study operates in a similar fashion to a standard bulge test and uses oil pressure to generate biaxial loading conditions. This biaxial testing apparatus was operated within a synchrotron beamline to characterize the mechanical response of a flash-processed steel sheet using in-situ high-energy X-ray diffraction (XRD) measurements. Further, the GSAS-II package was utilized to develop a workflow for the analysis of the large volume of diffraction data acquired. The workflow was then used to extract the peak position, width, and integrated intensity of the XRD peaks corresponding to the major body-centered cubic phase. Results: The equi-biaxial nature of the loading in the measured area was independently corroborated using experimental (XRD) and simulation (finite element analysis) methods. Furthermore, we discuss the evolution of elastic strain in the major body-centered cubic phase as a function of applied oil pressure and location on the steel sheet. Conclusions: A key advantage of the biaxial apparatus fabricated in this synchrotron study is demonstrated using the results obtained for the flash-processed steel sheet – i.e., mapping the lattice plane-dependent response to biaxial loading for a relatively large sample area in a spatially resolved manner.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A First Principles Approach to Spectral Phonon Transport in Heterostructures

Understanding thermal transport across interfaces which give rise to a thermal resistance (also known as Kapitza resistance) is a critical issue affecting the development of nanotechnologies. Much modern and emergent nanotechnology consist of adjacent materials, and phonon mediated heat transfer governs thermal behavior across internal interfaces in these devices. The physics of thermal transport in solids are governed both by phenomena occurring at the atomic scale and interactions with the material's microstructure. The forecasting of fundamental quantities such as temperature, heat flux and thermal conductivity typically employs the semi-classical Boltzmann transport equation to predict the macroscopic behavior of materials in terms of the microscopic dynamics of its heat carriers. Kapitza resistance was first discovered in liquid helium experiments and has led to a fundamental research thrust in micro and nano-scale heat transport, the behavior of thermal carriers across internal interfaces. Thermal interfacial resistance (TIR) is a widely studied phenomenon, first engaged by Swartz and Pohl through their development of the acoustic and diffuse mismatch methods, then continued through myriad efforts with varying methods and approaches in an attempt to resolve carrier behavior at thermal interfaces. Many of the fundamental approaches to TIR have been at the nanoscale, and research is conducted with molecular dynamics (MD) and density functional theory (DFT) methods. The limitations of these methods is system size; atomistic methods tend to be limited to system sizes of 100,000 atoms or less. Larger length-scale methods have also been pursued, based on the principles of acoustic or diffuse mismatch, but not all include simulation of TIR using a full phonon band spectrum, or temperature dependent methods. Our approach to enabling phonon transport in layered materials draws upon our previous work of demonstrating spectrally coupled phonon transport in homogeneous and heterogeneous materials. We use a semi-analytical approach in which the Bose-Einstein (B-E) statistics set the strength of the phonon radiance in a frequency group, but the B-E statistics are informed with information from the transport system. The B-E statistics in a single frequency group feels the influence of all the groups through the spatial temperature. We also include a new field term which is an indicator of the amount of non-equilibrium behavior of the phonon spectrum---this is added to the phonon source term in all groups to ensure closure and conservation of energy, as the phonon groups in the transport system and the analytical systems are coupled. This work builds upon our previous approach by adding a phonon coupling term at an internal interface, using the principles of the DMM through transmission and reflection coefficients. In this work, the coefficients are determined through computing a common temperature at the interface, influenced by the phonon band structure of both materials, in effect, providing mixing between the two material systems and using the common temperature to set the strength of the phonon radiance at the boundaries on either side of the interface. Our approach uses material properties computed along various crystallographic orientations, and while some isotropy is built into the interface condition, the material properties weight the phonon distributions in the proper crystalline direction. Greater resolution of phonon behavior in proximity to an interface, and more accurate predictions of TIR are obtained. While it is true the assumption of diffuse mismatch can yield inconsistent results compared to experiment especially at low temperatures, this work focuses on room temperature and beyond effects, for future applications in nuclear fuel, or thermoelectric devices; a modified mismatch approach may be feasible if applied properly. Additionally, our methods focus on bridging mesoscale to engineering scale

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of Electrode/Electrolyte Coupling on Birnessite (δ-MnO 2 ) Mechanical Response and Degradation

Understanding the deformation of energy storage electrodes at a local scale and its correlation to electrochemical performance is crucial for designing effective electrode architectures. In this work, the effect of electrolyte cation and electrode morphology on birnessite (δ-MnO 2 ) deformation during charge storage in aqueous electrolytes was investigated using a mechanical cyclic voltammetry approach via operando atomic force microscopy (AFM) and molecular dynamics (MD) simulation. In both K 2 SO 4 and Li 2 SO 4 electrolytes, the δ-MnO 2 host electrode underwent expansion during cation intercalation, but with different potential dependencies. When intercalating Li + , the δ-MnO 2 electrode presents a nonlinear correlation between electrode deformation and electrode height, which is morphologically dependent. Further, these results suggest that the stronger cation–birnessite interaction is the reason for higher local stress heterogeneity when cycling in Li 2 SO 4 electrolyte, which might be the origin of the pronounced electrode degradation in this electrolyte.

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

Invariant discovery of features across multiple length scales: Applications in microscopy and autonomous materials characterization

Physical imaging is a foundational characterization method in areas from condensed matter physics and chemistry to astronomy and spans length scales from atomic to universe. Images encapsulate crucial data regarding atomic bonding, materials microstructures, and dynamic phenomena such as microstructural evolution and turbulence, among other phenomena. The challenge lies in effectively extracting and interpreting this information. Variational Autoencoders (VAEs) have emerged as powerful tools for identifying the underlying factors of variation in image data, providing a systematic approach to distilling meaningful patterns from complex data sets. However, a significant hurdle in their application is the definition and selection of appropriate descriptors reflecting local structures. Here, we introduce the scale-invariant VAE approach (SI-VAE) based on the progressive training of the VAE with the descriptors sampled at different length scales. The SI-VAE allows the discovery of the length scale-dependent factors of variation in the system. Here, we illustrate this approach using the ferroelectric domain images and generalize it to the movies of the electron-beam induced phenomena in graphene and topography evolution across combinatorial libraries. This approach can further be used to initialize the decision making in automated experiments including structure–property discovery and can be applied across a broad range of imaging methods. This approach is universal and can be applied to any spatially resolved data including both experimental imaging studies and simulations, and can be particularly useful for exploration of phenomena such as turbulence and scale-invariant transformation fronts.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mechanically Robust Bismuth-Embedded Carbon Microspheres for Ultrafast Charging and Ultrastable Sodium-Ion Batteries

Advancements in the development of fast-charging and long-lasting microstructured alloying anodes with high volumetric capacities are essential for enhancing the operational efficiency of sodium-ion batteries (SIBs). These anodes, however, face challenges such as declined cyclability and rate capability, primarily due to mechanical degradation reduced by significant volumetric changes (over 252%) and slow kinetics of sodium-ion storage. Herein, we introduce a novel anode design featuring densely packed bismuth (Bi) embedded within highly conductive carbon microspheres to overcome the aforementioned challenges. Remarkably, the high loading Bi anode within carbon microspheres with a high tap density of 2.59 g cm -3 possesses significant mechanical strength exceeding 590 MPa and limits volume swelling of only 10.9% post-sodiation. This anode demonstrates a high volumetric capacity (908.3 mAh cm -3 ), ultrafast chargeability (200 A g -1 , full charge/discharge in just 5.5 s), and outstanding cyclability over 12,000 cycles and maintains exceptional cycling stability even at -30 °C. The full cell paired with a Na 3 V 2 (PO 4 ) 3 cathode retains over 80% capacity after 600 cycles at 36 C, demonstrating a remarkable rate capability of 126 C (full charge/discharge in 28.6 s). Our comprehensive experimental evaluations and chemo-mechanical simulations shed light on the mechanisms underpinning the anode's superior performance. In conclusion, this development marks a significant advancement in the design of durable and fast-charging anodes for high-performance SIBs.

25 ENERGY STORAGE↗
DOE OSTI2020

Synthesis of metallic glass nanoparticles by inert gas condensation

Inert Gas Condensation (IGC) is a versatile method for the synthesis of ultrafine nanoparticles. While it has been mostly used in the synthesis of nanocrystalline particles it has also been employed in the synthesis of metallic glass nanoparticles, which can be further used in the generation of nanostructured metallic glasses, known as nanoglasses. In this work, we use molecular dynamics simulations to investigate the IGC processes leading to the synthesis of Cu 64 Zr 36 amorphous nanoparticles. We investigate the growth of nanoparticles up to 7.3 nm diameter using long simulations of up to 130 ns, under relatively low pressures in the range of 4 to 10 bars. Results indicate a significant effect of pressure on the final structure of the amorphous nanoparticles, yet no clear influence on the segregation of Cu atoms to its surface. Cu segregation occurring at all pressures results in a thin Cu rich layer, about 2.5 Å thick, that coats all nanoparticles. Here, statistics of atomic Voronoi polyhedra shows that the fraction of full icosahedra in the nanoparticles produced at 4 bars is equivalent to that of bulk metallic glasses produced by melt and quench at 10 9 K/s. Higher pressures result in a lower fraction of full icosahedra.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tailoring composition and deformation modes at the microstructural level for next generation low-cost high-strength austenitic stainless steels

The objective of this project is to enable deliberate development of cost-effective, hydrogen resistant alloys by establishing detailed relationships specific to the effects of alloy composition, short-range order (SRO), and microsegregation in the presence of hydrogen on the transition between homogeneous deformation and localized plasticity in shear bands. In collaboration with the International Institute for Carbon-Neutral Energy Research, I2CNER, at Kyushu University in Japan, we conceptualized, designed, and manufactured four austenitic alloys that maintain corrosion resistance and ensure lower cost relative to baseline commercial alloys. The mechanical properties and deformation modes of the novel alloys (KU alloys) were assessed in the presence of hydrogen (H). Correlations between composition and performance revealed that two of the KU alloys are suitable replacements for 316 steel, while another is a viable replacement for 304 steel at room temperature. We found that, in the presence of other austenite stabilizing elements namely Mn and N, replacing Ni with Cu does not lead to martensite formation as has been previously reported.1–3 Furthermore, we found that the addition of Cu leads to an earlier onset of multiple slip resulting in an relative earlier onset of a higher work hardening rate (WHR). Greater understanding of the relationships between alloy composition and SRO required the development of a novel advanced electron diffraction methodology to characterize SRO in complex FCC alloys. This innovative approach, which combines fluctuation and correlation analyses of diffuse-scattering signals, successfully differentiated between SRO and long-range ordering (LRO). Further investigations into annealed austenitic stainless steels could provide insights into manipulating SRO and its effects on material properties. Atomistic simulations provided understanding of SRO behavior that was difficult to capture experimentally. This project created the first spin cluster expansion model that is able to capture and describe SRO effects in Fe-Ni-Cr FCC alloys, accounting for the non-negligible effects of magnetism. An automated computational workflow was established to provide reliable predictions of SRO in Fe-Ni-Cr austenitic alloys, both with and without the presence of H atoms. Analysis of the propensity for SRO in Fe-Ni-Cr alloys revealed that H tends to cluster with specific, well-defined SRO domains. The computational framework is general purpose and can be extended to realistic stainless steels across diverse composition ranges. With confidence that SRO is possible in austenitic stainless steels, we developed a discrete dislocation finite element code to understand the interaction of dislocations with SRO in the presence of H. By incorporating H effects on the dislocation emission and SRO stress field we show that the critical stress for the dislocation pileup to breakthrough the SRO domain decreases in the presence of H, which directly contributes localized deformation at the macroscale. Through the simulation of a uniaxial tension test, we demonstrated that H-induced weakening of SRO stress field and H-enhanced dislocation emission can lead to the onset of shear localization at lower macroscopic strains. As a whole, this project identified three novel alloys that show improvements in performance and cost efficiency for H-facing applications by studying correlations between alloy chemistry and deformation behavior. We also made significant advancements to experimental and computational methodologies necessary to study the chemistry and distribution of SRO across a range of alloys, which in turn allowed us to demonstrate how deformation mechanisms change due to the contributions of SRO in austenitic alloys in the presence of H. The combined advancements in fundamental understanding with novel alloy development in this project has increased the viability of next generation H-technologies for the broader public through accessible low-cost alloys and accelerated development towards future H-infrastructure.

08 HYDROGEN↗
DOE OSTI2025

Design and Development of Novel Equiatomic Refractory Multi Principal elemental Alloys Based on MoNbTi System for Use in Irradiation Environments

Multi Principal Elemental Alloys (MPEA) have emerged as promising materials for next-generation nuclear reactors due to their exceptional irradiation resistance. Eight equiatomic MPEA based on the MoNbTi system, comprising of elements with low thermal neutron absorption cross-sections were explored using a combined approach employing empirical parameter estimations, and CALPHAD simulations by which the phases and elemental segregation observed in all the alloys in the as-cast state were predicted. Solution heat treatment at 1500°C transformed five alloys into single-phase matrix materials, enhancing homogeneity and reducing hardness. The densities of the alloys ranged between 6.47 to 7.68 g/cm3, hardness between 472 and 656 VHN, Young’s modulus between 142 GPa to 169 GPa, shear modulus between 54 GPa and 62 GPa, bulk modulus between 117 GPa to 194 GPa and Poisson’s ratio between 0.3 to 0.35. The in-situ high temperature Xray diffraction results, differential scanning calorimetry and dilatometry results up to 1000°C suggested the high temperature phase stability of the MPEA. Subsequent ageing heat treatment at 800 and 1000 oC for 96 hours revealed significant secondary phase precipitation in MoNbTiZr, MoNbTiZrV, and MoNbTiCrA. Oxidation studies at 800 oC for 24 hours in air revealed superior oxidation resistance and cubic rate law dependence in Cr containing MPEA especially in MoNbTiCrAl, while severe mass gain resulting in total disintegration and exfoliation in MoNbTiZr and MoNbTiZrV. This comprehensive study underscores the potential of novel MPEA as promising materials for advanced nuclear reactor applications, shedding light on their microstructural control, mechanical properties, thermal stability, and oxidation resistance.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Towards understanding particle rigid-body motion during solid-state sintering

A quantitative understanding of particle rigid body (RB) motion that inherently accompanies grain boundary (GB) diffusion is highly desirable to understand and control the dynamic interplay between coarsening and densification during solid state sintering. In this work, by computer simulation using a multi-phase-field approach, we analyze systematically the roles played by each of these processes at different stages of the shrinkage of the internal pore in a three-particle green body as a function of particle size as well as thermodynamic and kinetic factors of interfaces. We demonstrate that particle RB translation promotes both neck growth, and pore rounding and shrinkage. Moreover, the forces acting at GBs and pulling neighboring particles towards one another dynamically evolve as particles fuse. In contrast, particle RB rotation has no contribution to pore shrinkage. The translational force acting on an individual particle varies with not only its size, but also the number and sizes of its neighboring particles.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of sintering temperature on adhesion of spray-on piezoelectric transducers

Conventionally sol-gel spray-on transducers require a high-temperature (> 700 ◦C) sintering process; however, this process can affect the microstructure of the substrate material. For mechanical elbows and valves utilized for fluid transport in the energy sector, the components are designed to have a specific microstructure, and deviations from these specifications can create weak points in the system. For this reason it is important to investigate how the temperature of the deposition process affects the substrate. This paper investigates the effect of high-temperature and low-temperature (< 150 ◦C) processing conditions on the surface composition of the substrate. Furthermore, the resultant transducers from high- and low-temperature fabrication processes are compared to determine if a low-temperature processing method is feasible. For these studies a sol-gel spray-on process is employed to deposit piezoelectric ceramics onto a stainless-steel 316L substrate. Energy-dispersive X-ray spectroscopy is utilized to determine the composition of the substrate surface before and after transducer deposition. Results indicate that the high-temperature processing conditions may alter the surface composition of the metal due to a diffusion of the metal into the ceramic, which results in a metal surface that is bonded to the ceramic. Furthermore, it is shown that low-temperature processing of spray-on transducers is a viable method for transducer fabrication where the resultant transducers meet the industry minimum requirement of 30 dB signalto-noise ratio. In parallel simulation calculations, finite-element method (FEM) studies were performed to model the adhesive strength of the low-temperature processed transducer to the substrate surface. Comparisons between the simulations and experiments suggest that the bond strength is much greater than the commercial gel bonds and closer to hardened epoxy glue bonds. These results indicate that spray-on transducers fabricated under lowtemperature processing conditions are a viable solution for leave-in-place monitoring of structures.

M. Sinding, Kyle↗
DOE OSTI2020

All-atom molecular dynamics simulations of weak polyionic brushes: influence of charge density on the properties of polyelectrolyte chains, brush-supported counterions, and water molecules

All atom molecular dynamics (MD) simulations of planar Na+-counterion-neutralized polyacrylic acid (PAA) brushes are performed for varying degrees of ionization (and thereby varying charge density) and varying grafting density. Variation in the PE charge density (or degree of ionization) and grafting density leads to massive changes of the properties of the PE molecules (quantified by the changes in the height and the mobility of the PE brushes) as well as the local arrangement and distribution of the brush-supported counterions and water molecules within the brushes. The effect on the counterions is manifested by the corresponding variation of the counterion mobility, counterion concentration, extent of counterion binding to the charged site of the PE brushes, water-in-salt-like structure formation, and counterion-water–oxygen radial distribution function within the PE brushes. On the other hand, the effect on water molecules is manifested by the corresponding variation of water–oxygen–water–oxygen RDF, local water density, water–water and water–PE functional group hydrogen bond networks, static dielectric constant of water molecules, orientational tetrahedral order parameter, and water mobility. Enforcing such varying degree of ionization of weak polyelectrolytes is possible by changing the pH of the surrounding medium. Thus, our results provide insights into the changes in microstructure (at the atomistic level) of weak polyionic brushes at varying pH. Finally, we anticipate that this knowledge will prove to be vital for the efficient design of several nano-scale systems employing PE brushes such as nanomechanical gates, current rectifiers, etc.

42 ENGINEERING↗
DOE OSTI2025

Multi-objective surrogate-assisted calibration of CPFEM models using macroscopic response and in situ EBSD measurements of grain reorientation trajectories

Crystal plasticity finite element method (CPFEM) models are widely used to simulate the deformation behaviour of polycrystalline materials, but their calibration is often limited by their high computational cost and the non-convexity of the optimisation landscape. Here, this study develops a multi-objective surrogate-assisted calibration workflow that couples a multi-objective genetic algorithm (MOGA) with an adaptively trained deep neural network (DNN) surrogate model to efficiently identify CPFEM parameters from experimental data. The workflow is demonstrated on three crystal plasticity (CP) formulations of increasing complexity — Voce hardening (VH), two-coefficient latent hardening (LH2), and six-coefficient latent hardening (LH6) — using in situ electron backscatter diffraction (EBSD) measurements of Alloy 617 under uniaxial tensile loading. The CPFEM models are calibrated against the experimentally observed stress–strain response and reorientation trajectories of eight grains, then validated against eight additional trajectories and overall texture evolution. Across the CP formulations, the macroscopic response was reproduced reliably, while differences emerged in the robustness and accuracy of the grain-scale predictions. Including grain reorientation trajectories in the multi-objective calibration improved texture evolution predictions and filtered out physically inconsistent parameter sets that can arise from calibrating against only the stress–strain data. The workflow also demonstrates good transferability of calibrated parameters from a low- to a high-fidelity microstructural model. These results provide practical guidance for integrating in situ microstructural data into CPFEM through efficient, repeatable, and physically meaningful multi-objective calibration.

Crystal plasticity finite element method↗
DOE OSTI2024

Comprehensive new insights on the potential use of SiC as plasma-facing materials in future fusion reactors

Abstract The performance of silicon carbide as an alternative plasma facing material (PFM) was studied at various irradiation conditions relevant to ion energies and fluxes of a fusion reactor. This analysis involves detailed modeling of subsurface plasma/material interactions, sputtered particle transport above the surface and redeposition, and related changes in material composition and microstructure induced by steady-state and Edge Localized Mode ion fluxes. Transition of a crystalline SiC surface to semi-crystalline and amorphous phases was analyzed based on advanced modeling of DIII-D tokamak experiments where SiC was irradiated in single- and multiple- L-mode and H-mode discharges. This analysis shows that displacement damage, particle deposition/redeposition, and D accumulation on the SiC divertor surface can lead to significant microstructural changes that result in enhanced sputtering erosion in comparison with the original crystalline material. However, the resulting total net erosion rate for a full-coverage, advanced tokamak, SiC coated divertor may well be acceptably low. Moreover, the C sputtering yield from the evolved SiC surface can be seven times lower than from a pure graphite surface; this would imply significantly reduced tritium co-deposition rates in a D-T tokamak reactor, compared with a pure carbon surface. It was also determined that chemical sputtering of both C and Si should not result in any noticeable effect on the net erosion, for attached plasma regimes. Our results thus show encouraging results overall for use of SiC as a PFM in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Crystal plasticity finite element simulation of lattice rotation and x-ray diffraction during laser shock compression of tantalum

Here we present a crystal plasticity model tailored for high-pressure, high-strain-rate conditions that uses a multiscale treatment of dislocation-based slip kinetics. We use this model to analyze the pronounced plasticity-induced lattice rotations observed in shock-compressed polycrystalline tantalum via in situ x-ray diffraction. By making direct comparisons between experimentally measured and simulated texture evolution, we can explain how the details of the underlying slip kinetics control the degree of lattice rotation that ensues. Specifically, we show that only the highly nonlinear kinetics caused by dislocation nucleation can explain the magnitude of the rotation observed under shock compression. We demonstrate a good fit between our crystal plasticity model and x-ray diffraction data and exploit the data to quantify the dislocation nucleation rates that are otherwise poorly constrained by experiment in the dynamic compression regime.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Calcium silicate hydrate colloid at different humidities: Microstructure, deformation mechanism, and mechanical properties

Calcium silicate hydrate (C-S-H) is the main hydration product of Portland cement, while fundamental questions persist regarding its colloidal structure and water content. In this work, we developed a full atomistic C-S-H colloid model and investigated the physical and mechanical properties at different gel water concentrations. Water confined in nanometer gel pores degrades the mechanical properties of C-S-H because of its hydrolysis effect. A “ductile to brittle” transition was observed in the C-S-H colloid model with increasing gel water contents, evolving from transgranular failure to intergranular fracture. The gel water confined in nanometer-sized pores displays glassy characteristics, and it acts as a lubricant and facilitates structure distortions during the deformation. We elucidated the topology effect of the hydrogen bond network on the evolution of the viscoelasticity of C-S-H at different humidities. In conclusion, this work sheds light on the role of water in C-S-H colloids from nanostructure perspectives.

36 MATERIALS SCIENCE↗