In this study, we experimentally investigate the high stain rate and spall behavior of Cantor high-entropy alloy (HEA), CoCrFeMnNi. First, the Hugoniot equations of state (EOS) for the samples are determined using laser-driven CoCrFeMnNi flyers launched into known Lithium Fluoride (LiF) windows. Photon Doppler Velocimetry (PDV) recordings of the velocity profiles find the EOS coefficients using an impedance mismatch technique. Following this set of measurements, laser-driven aluminum flyer plates are accelerated to velocities of 0.5–1.0 km/s using a high-energy pulse laser. Upon impact with CoCrFeMnNi samples, the shock response is found through PDV measurements of the free surface velocities. From this second set of measurements, the spall strength of the alloy is found for pressures up to 5 GPa and strain rates in excess of 106 s-1. Further analysis of the failure mechanisms behind the spallation is conducted using fractography revealing the occurrence of ductile fracture at voids presumed to be caused by chromium oxide deposits created during the manufacturing process.
Phase-field modeling is an effective but computationally expensive method for capturing the mesoscale morphological and microstructure evolution in materials. Hence, fast and generalizable surrogate models are needed to alleviate the cost of computationally taxing processes such as in optimization and design of materials. The intrinsic discontinuous nature of the physical phenomena incurred by the presence of sharp phase boundaries makes the training of the surrogate model cumbersome. We develop a framework that integrates a convolutional autoencoder architecture with a deep neural operator (DeepONet) to learn the dynamic evolution of a two-phase mixture and accelerate time-to-solution in predicting the microstructure evolution. We utilize the convolutional autoencoder to provide a compact representation of the microstructure data in a low-dimensional latent space. After DeepONet is trained in the latent space, it can be used to replace the high-fidelity phase-field numerical solver in interpolation tasks or to accelerate the numerical solver in extrapolation tasks.
Embrittling potency is a thermodynamic metric that assesses the influence of solute segregation to a grain boundary (GB) on intergranular fracture. Historically, authors of studies have reported embrittling potency as a single scalar value, assuming a single segregation site of importance at a GB and a particular cleavage plane. However, the topography of intergranular fracture surfaces is not generally known a priori. Accordingly, we, in this study, present a statistical ensemble approach to compute embrittling potency, where many free surface (FS) permutations are systematically considered to model fracture of a GB. The result is a statistical description of the thermodynamics of GB embrittlement. As a specific example, embrittling potency distributions are presented for Cr segregation to sites at two Ni $\langle 111 \rangle$ symmetric tilt GBs using atomistic simulations. We show that the average embrittling potency for a particular GB site, considering an ensemble of FS permutations, is not equal to the embrittling potency computed using the lowest energy pair of FSs. A mean GB embrittlement is proposed, considering both the likelihood of formation of a particular FS and the probability of solute occupancy at each GB site, to compare the relative embrittling behavior of two distinct GBs.
Elmslie, Timothy A.; Startt, Jacob K.; Soto-Medina, Sujeily S.; Fen, Keke F.; Zappala, Emma Z.; Frandsen, Benjamin A.; Meisel,
M.; Dingreville, Remi P.; Hamlin, james H.
Magnetic, specific heat, and structural properties of the equiatomic Cantor alloy system are reported for temperatures between 5 and 300 K, and up to fields of 70 kOe. Here magnetization measurements performed on as-cast, annealed, and cold-worked samples reveal a strong processing history dependence and that high-temperature annealing after cold working does not restore the alloy to a “pristine” state. Measurements on known precipitates show that the two transitions, detected at 43 and 85 K, are intrinsic to the Cantor alloy and not the result of an impurity phase. Experimental and ab initio density functional theory computational results suggest that these transitions are a weak ferrimagnetic transition and a spin-glass-like transition, respectively, and magnetic and specific heat measurements provide evidence of significant Stoner enhancement and electron-electron interactions within the material.
Radiation-induced segregation is a phenomenon commonly observed in many alloys which consists of the redistribution of elements (solute or interstitial impurities) under irradiation. The onset and development of radiation-induced segregation can only occur when a sufficient flux of defects is sustained and defect sinks are present. Irradiation dose, dose rate, and particle types all affect defect flux. In this work, we employ a phase-field model to examine the effects of dose, dose rate, and type of incident particles on radiation-induced segregation behavior in a model binary alloy. The phase-field model takes into account the formation and evolution of point defects as well as defect clusters, the diffusion and clustering of alloy species, the presence of additional extrinsic defect sinks in the form of dislocations, and two different methods of radiation-damage insertion, which are intended to simulate either light-ion/electron irradiation via Frenkel pairs or heavy-ion irradiation in the form of cascades. Our results show a dose-rate and particle-type dependence on the amount of solute segregation. We show that the material systems exposed to higher dose rates are less subjected to solute segregation at equivalent doses. We also show that such dose-rate-dependence behavior is due to a delay of the incubation dose at which radiation-induced segregation effectively starts. Particle type and the presence of dislocations can accentuate this behavior. Our model predictions correlate with many experimental observations made over the years on radiation-induced segregation providing credence to the simulation results. The methodology presented in this study allows for a first-order prediction of the dose rate at which proxy irradiation experiments could be performed to approximate radiation-induced segregation behaviors seen in targeted irradiation conditions.
Metals subjected to irradiation environments undergo microstructural evolution and concomitant degradation, yet the nanoscale mechanisms for such evolution remain elusive. Here, we combine in situ heavy ion irradiation, atomic resolution microscopy, and atomistic simulation to elucidate how radiation damage and interfacial defects interplay to control grain boundary (GB) motion. While classical notions of boundary evolution under irradiation rest on simple ideas of curvature-driven motion, the reality is far more complex. Focusing on an ion-irradiated Pt Σ3 GB, we show how this boundary evolves by the motion of 120° facet junctions separating nanoscale {112} facets. Our analysis considers the short- and mid-range ion interactions, which roughen the facets and induce local motion, and longer-range interactions associated with interfacial disconnections, which accommodate the intergranular misorientation. We suggest how climb of these disconnections could drive coordinated facet junction motion. These findings emphasize that both local and longer-range, collective interactions are important to understanding irradiation-induced interfacial evolution.
In this work, we present a combination of machine-learned models that predicts the surface elastic properties of general free surfaces in face-centered cubic (FCC) metals. These models are built by combining a semi-analytical method based on atomistic simulations to calculate surface properties with the artificial neural network (ANN) method or the boosted regression tree (BRT) method. The latter is also used to link bulk properties and surface orientation to surface properties. The surface elastic properties are represented by their invariants considering plane elasticity within a polar method. The resulting models are shown to accurately predict the surface elastic properties of seven pure FCC metals (Cu, Ni, Ag, Au, Al, Pd, Pt). The BRT model reveals the correlations between bulk and corresponding surface properties in terms of invariants, which can be used to guide the design of complex nano-sized particles, wires and films. Finally, by expressing the surface excess energy density as a function of surface elastic invariants, fast predictions of surface energy as a function of in-plane deformations can be made from these model constructs.
Digital twins are emerging as powerful tools for supporting innovation as well as optimizing the in-service performance of a broad range of complex physical machines, devices, and components. A digital twin is generally designed to provide accurate in-silico representation of the form (i.e., appearance) and the functional response of a specified (unique) physical twin. This paper offers a new perspective on how the emerging concept of digital twins could be applied to accelerate materials innovation efforts. Specifically, it is argued that the material itself can be considered as a highly complex multiscale physical system whose form (i.e., details of the material structure over a hierarchy of material length) and function (i.e., response to external stimuli typically characterized through suitably defined material properties) can be captured suitably in a digital twin. Accordingly, the digital twin can represent the evolution of structure, process, and performance of the material over time, with regard to both process history and in-service environment. This paper establishes the foundational concepts and frameworks needed to formulate and continuously update both the form and function of the digital twin of a selected material physical twin. The form of the proposed material digital twin can be captured effectively using the broadly applicable framework of n-point spatial correlations, while its function at the different length scales can be captured using homogenization and localization process-structure-property surrogate models calibrated to collections of available experimental and physics-based simulation data.
Garnet-type, solid electrolytes, such as Li7La3Zr2O12 (LLZO), are a promising alternative to liquid electrolytes for lithium-metal batteries. However, such solid-electrolyte materials frequently exhibit undesirable lithium (Li) metal plating and fracture along grain boundaries. In this study, we employ atomistic simulations to investigate the mechanisms and key fracture properties associated with intergranular fracture along one such boundary. Our results show that, in the case of a Σ5(310) grain boundary, this boundary exhibits brittle fracture behavior, i.e. the absence of dislocation activity ahead of the propagating crack tip, accompanied with a decrease in work of separation, peak stress, and maximum stress intensity factor as the temperature increases from 300 K to 1500 K. As the crack propagates, we predict two temperature-dependent Li clustering regimes. For temperatures at or below 900 K, Li tends to cluster in the bulk region away from the crack plane driven by a void-coalescence mechanism concomitant a simultaneous cubic-to-tetragonal phase transition. The tetragonalization of LLZO in this temperature regime acts as an emerging toughening mechanism. At higher temperatures, this phase transition mechanism is suppressed leading to a more uniform distribution of Li throughout the grain-boundary system and lower fracture properties as compared to lower temperatures.
Structural alloys may experience corrosion when exposed to molten chloride salts due to selective dissolution of active alloying elements. One way to prevent this is to make the molten salt reducing. For the KCl + MgCl2 eutectic salt mixture, pure Mg can be added to achieve this. However, Mg can form intermetallic compounds with nickel at high temperatures, which may cause alloy embrittlement. This study shows that an optimum level of excess Mg could be added to the molten salt which will prevent corrosion of alloys like 316 H, while not forming any detectable Ni-Mg intermetallic phases on Ni-rich alloy surfaces.
Echeverria, Marco J.; Galitskiy, Sergey; Mishra, Avanish; Dingreville, Remi P.; Dongare, Avinash M.
A hybrid atomic-scale and continuum-modeling framework is used to study the microstructural evolution during the laser-induced shock deformation and failure (spallation) of copper microstructures. A continuum two-temperature model (TTM) is used to account for the interaction of Cu atoms with a laser in molecular dynamics (MD) simulations. The MD-TTM simulations study the effect of laser-loading conditions (laser fluence) on the microstructure (defects) evolution during various stages of shock wave propagation, reflection, and interaction in single-crystal (sc) Cu systems. In addition, the role of the microstructure is investigated by comparing the defect evolution and spall response of sc-Cu and nanocrystalline Cu systems. The defect (stacking faults and twin faults) evolution behavior in the metal at various times is further characterized using virtual in situ selected area electron diffraction and x-ray diffraction during various stages of evolution of microstructure. The simulations elucidate the uncertain relation between spall strength and strain-rate and the much stronger relation between the spall strength and the temperatures generated due to laser shock loading for the small Cu sample dimensions considered here.
Metamaterials are artificial structures that can manipulate and control sound waves in ways not possible with conventional materials. While much effort has been undertaken to widen the bandgaps produced by these materials through design of heterogeneities within unit cells, comparatively little work has considered the effect of engineering heterogeneities at the structural scale by combining different types of unit cells. In this paper, we use the relaxed micromorphic model to study wave propagation in heterogeneous metastructures composed of different unit cells. We first establish the efficacy of the relaxed micromorphic model for capturing the salient characteristics of dispersive wave propagation through comparisons with direct numerical simulations for two classes of metamaterial unit cells: namely phononic crystals and locally resonant metamaterials. We then use this model to demonstrate how spatially arranging multiple unit cells into metastructures can lead to tailored and unique properties such as spatially-dependent broadband wave attenuation, rainbow trapping, and pulse shaping. In the case of the broadband wave attenuation application, we show that by building layered metastructures from different metamaterial unit cells, we can slow down or stop wave packets in an enlarged frequency range, while letting other frequencies through. In the case of the rainbow-trapping application, we show that spatial arrangements of different unit cells can be designed to progressively slow down and eventually stop waves with different frequencies at different spatial locations. Finally, in the case of the pulse-shaping application, our results show that heterogeneous metastructures can be designed to tailor the spatial profile of a propagating wave packet. Collectively, these results show the versatility of the relaxed micromorphic model for effectively and accurately simulating wave propagation in heterogeneous metastructures, and how this model can be used to design heterogeneous metastructures with tailored wave propagation functionalities.
Designing next generation thin films, tailor-made for specific applications, relies on the availability of robust processing-structure-property relationships. Traditional structure zone diagrams are limited to low-dimensional mappings, with machine-learning methods only recently attempting to relate multiple processing parameters to the final microstructure. Despite this progress, structure-processing relationships are unknown for processing conditions that vary during thin-film deposition, limiting the range of microstructures and properties achievable. In this project, we employed a phase-field computational model combined with a genetic algorithm (GA) to identify and design time-dependent processing protocols that achieve tailor-made microstructures. We simulate the physical vapor deposition of a binary-alloy thin film by employing a phase-field model, where deposition rates and diffusivities are controlled via the genetic algorithm. Our GA-guided protocols achieve targeted microstructures with lateral and vertical concentration modulations, as well as more complex, hierarchical microstructures previously not described in simple structure zone diagrams. Our algorithm provides insight to experimentalists looking for additional avenues to design novel thin-film microstructures.
With the rapid proliferation of additive manufacturing and 3D printing technologies, architected cellular solids including truss-like 3D lattice topologies offer the opportunity to program the effective material response through topological design at the mesoscale. The present report summarizes several of the key findings from a 3-year Laboratory Directed Research and Development Program. The program set out to explore novel lattice topologies that can be designed to control, redirect, or dissipate energy from one or multiple insult environments relevant to Sandia missions, including crush, shock/impact, vibration, thermal, etc. In the first 4 sections, we document four novel lattice topologies stemming from this study: coulombic lattices, multi-morphology lattices, interpenetrating lattices, and pore-modified gyroid cellular solids, each with unique properties that had not been achieved by existing cellular/lattice metamaterials. The fifth section explores how unintentional lattice imperfections stemming from the manufacturing process, primarily sur face roughness in the case of laser powder bed fusion, serve to cause stochastic response but that in some cases such as elastic response the stochastic behavior is homogenized through the adoption of lattices. In the sixth section we explore a novel neural network screening process that allows such stocastic variability to be predicted. In the last three sections, we explore considerations of computational design of lattices. Specifically, in section 7 using a novel generative optimization scheme to design novel pareto-optimal lattices for multi-objective environments. In section 8, we use computational design to optimize a metallic lattice structure to absorb impact energy for a 1000 ft/s impact. And in section 9, we develop a modified micromorphic continuum model to solve wave propagation problems in lattices efficiently.
This project focused on providing a fundamental physico-chemical understanding of the coupling mechanisms of corrosion- and radiation-induced degradation at material-salt interfaces in Ni-based alloys operating in emulated Molten Salt Reactor(MSR) environments through the use of a unique suite of aging experiments, in-situ nanoscale characterization experiments on these materials, and multi-physics computational models. The technical basis and capabilities described in this report bring us a step closer to accelerate the deployment of MSRs by closing knowledge gaps related to materials degradation in harsh environments.
Twin boundaries play an important role in the thermodynamics, stability, and mechanical properties of nanocrystalline metals. Understanding their structure and chemistry at the atomic scale is key to guide strategies for fabricating nanocrystalline materials with improved properties. In this work, we report an unusual segregation phenomenon at gold-doped platinum twin boundaries, which is arbitrated by the presence of disconnections, a type of interfacial line defect. By using atomistic simulations, we show that disconnections containing a stacking fault can induce an unexpected transition in the interfacial-segregation structure at the atomic scale, from a bilayer, alternating-segregation structure to a trilayer, segregation-only structure. This behavior is found for faulted disconnections of varying step heights and dislocation characters. Supported by a structural analysis and the classical Langmuir–McLean segregation model, we reveal that this phenomenon is driven by a structurally induced drop of the local pressure across the faulted disconnection accompanied by an increase in the segregation volume.
During the various stages of shock loading, many transient modes of deformation can activate and deactivate to affect the final state of a material. In order to fundamentally understand and optimize a shock response, researchers seek the ability to probe these modes in real-time and measure the microstructural evolutions with nanoscale resolution. Neither post-mortem analysis on recovered samples nor continuum-based methods during shock testing meet both requirements. High-speed diffraction offers a solution, but the interpretation of diffractograms suffers numerous debates and uncertainties. By atomistically simulating the shock, X-ray diffraction, and electron diffraction of three representative BCC and FCC metallic systems, we systematically isolated the characteristic fingerprints of salient deformation modes, such as dislocation slip (stacking faults), deformation twinning, and phase transformation as observed in experimental diffractograms. This study demonstrates how to use simulated diffractograms to connect the contributions from concurrent deformation modes to the evolutions of both 1D line profiles and 2D patterns for diffractograms from single crystals. Harnessing these fingerprints alongside information on local pressures and plasticity contributions facilitate the interpretation of shock experiments with cutting-edge resolution in both space and time.