SEARCH · Engineering Papers
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.
Presentation slides on Atomistically informed 3D phase field dislocation dynamics (PFDD) model for nanomaterials; Nanocrystalline Metals; Atomistically Informed PFDD; Phase Field Dislocations; DFT-Informed PFDD; Grain Size Effects; and Interface-driven Plasticity.
This study focused on the analysis of geometrically necessary dislocation (GND) densities for five selected fine-grained magnesium samples. Among the samples, three were tested under different fatigue-loading conditions at 0 °C, one experienced quasi-static tensile loading at 0 °C, and one represented the as-rolled state. The fatigue-tested samples were chosen according to the relationship between the maximum loading stress of a test and the material’s yield strength. This study provides new insights on the deformation mechanism of fine-grained magnesium at 0 °C. It is observed that the average GND densities were increased by 95~111% for the tested samples when compared with the as-rolled sample. It is especially interesting that there is a significant increase in the average GND density for the sample that experienced the fatigue loading with a low-maximum applied stress, and the maximum applied stress was lower than the material’s yield strength. This observation implies that the grain boundary mediated the dislocation-emission mechanism.
Efforts are under way to develop a theoretical, multi-scale model for the prediction of fracture toughness of ferritic steels in the ductile-to-brittle transition temperature (DBTT) region that accounts for temperature, irradiation, strain rate, and material condition (chemistry and heat treatment) effects. This new model is intended to address difficulties associated with existing empirically-derived models of the DBTT region that cannot be extrapolated to conditions for which data are unavailable. Dislocation distribution equations, derived from the theories of Yokobori et al., are incorporated to account for the local stress state prior to and following initiation of a microcrack from a second-phase particle. The new model is the basis for the DISlocation-based FRACture (DISFRAC) computer code being developed at the Oak Ridge National Laboratory (ORNL). The purpose of this code is to permit fracture safety assessments of ferritic structures with only tensile properties required as input. The primary motivation for the code is to assist in the prediction of radiation effects on nuclear reactor pressure vessels, in parallel with the EURATOM PERFORM 60 project.
A crystalline material structure with reduced dislocation density and method of producing same is provided. The crystalline material structure is annealed at temperatures above the brittle-to-ductile transition temperature of the crystalline material structure. One or more stress elements are formed on the crystalline material structure so as to annihilate dislocations or to move them into less harmful locations.
Material defects play a large role in material response under shock loading, yet the process of how these defects initiate, propagate, and annihilate is not well understood. Using a newly developed diagnostic, dark field X-ray microscopy (DFXM), we can now visualize the behavior of dislocation defects in materials at the mesoscale under varying conditions. Using data from the European Synchrotron Radiation Facility, we apply a variety of image processing techniques to capture relevant features to locate and characterize size and orientation of dislocation defects in bulk aluminum DFXM images. Beyond simple visualization, this analysis drives statistical characterization of the defects and their dynamic response to further improve the relevant physics models.
Using a partitioned-energy thermodynamic framework which assigns energy to that of atomic configurational stored energy of cold work and kinetic-vibrational, in this study we derive an important constraint on the Taylor-Quinney coefficient, which quantifies the fraction of plastic work that is converted into heat during plastic deformation. Associated with the two energy contributions are two separate temperatures – the ordinary temperature for the thermal energy and the effective temperature for the configurational energy. We show that the Taylor-Quinney coefficient is a function of the thermodynamically defined effective temperature that measures the atomic configurational disorder in the material. Finite-element analysis of recently published experiments on the aluminum alloy 6016-T4 [1], using the thermodynamic dislocation theory (TDT), shows good agreement between theory and experiment for both stress-strain behavior and temporal evolution of the temperature. The simulations include both conductive and convective thermal energy loss during the experiments, and significant thermal gradients exist within the simulation results. Computed values of the differential Taylor-Quinney coefficient are also presented and suggest a value which differs between materials and increases with increasing strain.
Continuum-scale material deformation models, such as crystal plasticity (CP), can significantly enhance their predictive accuracy by incorporating input from lower-scale (i.e. mesoscale) models. The procedure to generate and extract the relevant information is however typically complex and ad hoc, involving decision and intervention by domain experts, leading to long development times. In this study, we develop a principled approach for calibration of continuum-scale models using lower scale information by representing a CP flow rule as a Gaussian process model. This representation allows for efficient parameter space exploration, guided by the uncertainty embedded in the model through a process known as Bayesian optimization (BO). We demonstrate a semi-autonomous BO loop which instantiates discrete dislocation dynamics simulations whose initial conditions are automatically chosen to optimize the uncertainty of a model CP flow rule. Our self-guided computational pipeline efficiently generated a dataset and corresponding model whose error, uncertainty, and physical feature sensitivities were validated with comparison to an independent dataset four times larger, demonstrating a valuable and efficient active learning implementation readily transferable to similar material systems.
Abstract Single crystals fabricated in glass by localized heating can develop uniquely deformed lattices stabilized by the surrounding amorphous medium. The development of lattice curvature appears to be intrinsic to the crystal growth process in some systems, while the result of the locally changing crystallography in others. In this work, a model laser‐fabricated rotating lattice Sb 2 S 3 crystal grown in stoichiometric glass is used to demonstrate fabrication of novel macroperiodic metastructures that utilize intrinsic lattice curvature superimposed with subtle crystallographic influences. The limited availability of slip systems drives the lattice curvature magnitude to vary with crystal growth direction, maximizing for lattices aligned with the predominant Burgers vector along with corresponding increases in dislocation density. Misaligned lattice orientations form smaller secondary lattice curvatures arising from misaligned Burgers vectors with further elastic contributions. Over extended crystal growth, these secondary components align the lattice to rotate about either the <001> or <010> crystal axes forming repeating metastructures of lattice orientation with periodicity 20–160 microns in length. The mechanistic approach used in this work may be expanded to other systems with known slip systems to better understand and design macroperiodic metastructures.
In this study, the fatigue-induced microstructure produced in a nickel-based polycrystalline superalloy that was subjected to cyclic loading was characterized by polychromatic x-ray microdiffraction (PXM) together with in-situ neutron diffraction and transmission-electron microscopy (TEM). In-situ neutron-diffraction measurements reveal two distinct stages of the fatigue damage: cyclic hardening followed by cyclic softening. Three-dimensional spatially resolved PXM micro-Laue measurements find an increase in the density of geometrically necessary dislocations near the grain boundaries, which is accompanied by lattice rotations and grain subdivisions. The PXM results are in agreement with the in-situ neutron-diffraction and TEM results.
Misfit dislocation arrays (MDAs) at semi-coherent heterogeneous interface have been demonstrated effectiveness to trap and then outgas helium (He) for alleviating void swelling in metals, but limited by atomic-level thin film deposition that is difficult to scale-up for industrial production. In this work, we proposed an innovative strategy that could generate semi-coherent homophase interfaces in a same grain by varying surface energy density assisted by high-temperature graphene (Gr) deposition. This strategy transforms strain energy and surface energy, two originally distinctive processes into successive processes, and demonstrates experimentally and theoretically by introducing {111} copper (Cu) in a {100} Cu grain. Subsequent He implantation and electrical resistivity measurement results suggest that high-density MDAs originated at this {111}/{100} interface could form stable He nanochannels near Cu surface with Gr assistance, manifesting by ultra-low increased ratio of the electrical resistivity (i.e. ~3800% of Cu vs. ~90% of Gr/Cu) after He implantation. Finally, the present findings may provide a new strategy for efficiently managing He and achieving surface protection of advanced metals, merely by depositing a super-rigid material on metal surface to greatly reduce surface energy density.
Here, in this work, the intragranular distributions of the densities of geometrically necessary dislocations (GNDs) generated in a polycrystalline tantalum (Ta) sample loaded in shock compression to spall are investigated. The sample studied is a highly resolved three-dimensional (3D) map of the orientations, grain boundary morphology, and voids across a 0.276mm 3 (with sample dimensions of 688.5µm x 669µm x 600µm) region of the microstructure about the spall plane obtained by TriBeam tomography. We present a formulation for calculating the crystallographic screw and edge character GND densities at a material point within body centered cubic (BCC) crystals when provided any number of potential slip systems or modes. Combining the 3D characterization, GND formulation, and a sample containing ≈ 6000 grains, we sought statistically relevant correlations between GND density per grain with grain size, grain orientation, void formation, and proximity from the spall plane. We show that GND density per grain increases with decreasing distance from the spall plane. Most grains had heterogeneously distributed GND densities but were free of any detectable voids, at least within the ≈ 1.5µm resolution of the characterization technique. It is found that the GND density for most grains is insensitive to grain size and grain orientation. Only the fraction of grains with the highest GND densities tended to have a [001] parallel to shock-direction orientation. Intragranular 3D maps of GNDs reveals that grains that contain voids have high GND density concentrations in the intragranular region immediately surrounding the void. The local intragranular GND mapping also showed voids that are connected to grain boundaries preferred to grow into the grain with the larger size, regardless of orientation.
Here we introduce Interfacial Line Defect Analysis (ILDA), a method for identifying and extracting interfacial dislocations and disconnections with little-to-no input from the user on the nature of the interface. By simply providing the orientations and coherency strains for the crystals in a coherent reference state of the interface, ILDA provides exact Burgers vectors. Alternatively, these orientations and strains can be estimated using local atomic deformation gradients, making ILDA fully automated and providing estimated Burgers vectors. ILDA also determines the step height associated with each defect line segment, in case the associated defect is a disconnection. The heart of the method is the identification of atoms residing at co-incidence sites between the two crystals and the construction of a surface mesh connecting these sites that is used to compose Burgers circuits and insert defect line segments. We demonstrate the performance of ILDA in two test cases: a twist grain boundary and a phase boundary.
Models for grain size dependence of flow stress typically ignore contributions from geometrically necessary dislocations (GNDs) in the grain boundary regions. Ashby had proposed the following equation for flow stress (τ) as a function of grain size (d) and effective plastic strain ($\bar{ε}$): $τ = τ_{0} + C'G\sqrt{^{b\bar{ε}}_{d}}$ using an unknown empirical coefficient, C’, In this study, the coefficient, C’, and equation for the grain size dependence of flow stress were investigated in Mg-4Al using nanoindentation. C’ ranges from 1.4 to 7.11 that is related to the hardness difference (∆H0) between the adjacent grains as $C' + 0.62\sqrt{ΔH_{0}}$. Ashby equation was used to predict uniaxial stress-strain curves in grain sizes of 55 µm, 187 µm, and 333 µm. The results showed lower flow stress and lower hardening rate compared to experimentally measured tensile stress-strain response. The discrepancy seems to be a result of hardening effect of SSDs. Here, we propose a modified Ashby equation that takes into account the hardening effect of SSDs and GNDs separately, resulting in a better agreement between calculated and experimentally measured values. This analysis demonstrates estimation of uniaxial stress-strain response using nanoindentation measurements that captures effects of grain size, and strain hardening from SSDs and GNDs. The combination of the revised equation and nanoindentation method will enable rapid exploration of stress-strain response of a diverse range of single-phase polycrystalline alloys.
The influence of microstructural anisotropy on shear response of high-purity titanium was studied using the compact forced-simple-shear specimen (CFSS) loaded under quasi-static loading conditions. Post-mortem characterization reveals significant difference in shear response of different directions in the same material due to material crystallographic texture anisotropy. Shear bands are narrower in specimens in which the shear zone is aligned along the direction with a strong {0001} basal texture. Twinning was identified as an active mechanism to accommodate strains in the shear region in both orientations. This paper confirms the applicability of the CFSS design for the investigation of differences in the shear response of materials as a function of process-induced crystallographic texture. A detailed, systematic approach to quantifying shear band evolution by evaluating geometrically necessary dislocations (GND) associated with crystallographic anisotropy is presented. Finally, the results show that: i) line average GND density profiles, for Ti samples that possess a uniform equiaxed-grain structure, but with strong crystallographic anisotropy, exhibit significant differences in GND density close to the shear band center; ii) GND profiles decrease steadily away from the shear band as the plastic strain diminishes, in agreement with Ashby's theory of work hardening, where the higher GND density in the through-thickness (TT) orientation is a result of restricted < a > type slip in the shear band compared with in-plane (IP) samples; iii) the anisotropy in deformation response is derived from initial crystallographic texture of the materials, where GND density of < a > GNDs are higher adjacent to the shear band in the through-thickness sample oriented away from easy slip, but the density of < c+a > type GNDs are very similar in these two samples; and iv) the increase in grain average GND density was determined to have strong correlation to an increase in the Euler Φ angle of the grain average orientation, indicating an increased misorientation angle evolution.
Using the recently developed Bond Order Potential for Cr, we investigate the mobility of isolated 1/2[111] screw dislocations under stress at 0 K. These studies provide the dependence of the critical resolved shear stress (i.e., the Peierls stress) on the orientation of the maximum resolved shear stress plane and on the shear stress perpendicular to the slip direction. Similar studies are made for loading in tension and compression to show that the onset of yielding in Cr single crystals can be described using a combination of the shear stresses parallel and perpendicular to the slip direction. Based on these results, we develop a single crystal yield criterion for Cr that includes all non-Schmid contributions to the Peierls stress. This criterion is used to predict the shape of the yield surface and compare it with the regular Tresca's hexagon arising from the Schmid law. For all directions of uniaxial load, the criterion is used to quantify the character and magnitude of the tension-compression asymmetry.
Here, in this paper, we investigate plasticity in irradiated FeCrAl nanopillars using discrete dislocation dynamics simulations (DDD), with comparisons to transmission electron microscopic (TEM) in situ tensile tests of ion- and neutron-irradiated commercial C35M FeCrAl alloy. The effects of irradiation-induced defects, such as a/2<111> and a<100> type loops and composition fluctuations representative of phase separation in irradiated FeCrAl alloys, are investigated separately as well as superposed together in simulations. We explore the effects of defects on the stress-strain behavior, specifically yield strength and hardening response, of FeCrAl nanopillars. Our simulations confirm the widely accepted fact that irradiated alloys exhibit a stress-strain response with higher yield strength as compared to unirradiated alloys. However, our DDD calculations reveal an atypical superposition of the hardening contributions due to composition inhomogeneity and irradiation loops wherein composition inhomogeneity annihilates the hardening due to irradiation loops at small scales. As a result, we observe that the yield strength in irradiated alloys, after taking into consideration the effects of both composition inhomogeneity and irradiation loops, is smaller than the yield strength of the alloys with only irradiation loops and is approximately same for the alloy with composition inhomogeneity alone. This is referred to as “destructive interference” between the hardening contributions due to composition fluctuations and irradiation loops in the paper. We also identify this destructive interference in the superposition in our parallel TEM in situ tensile tests on unirradiated, ion-irradiated, and neutron-irradiated C35M FeCrAl alloy. This destructive interference in the hardening contributions contrasts with the dispersed barrier hardening (DBH) models widely utilized by the experimental community to model the hardening contributions due to different irradiation induced defects. The effects of the loading orientations on the yield strength and hardening are investigated and the mechanisms for the hardening in irradiated FeCrAl alloys are also reported.
In situ tensile testing in a scanning electron microscope (SEM) in conjunction with high-resolution electron backscatter diffraction (HR-EBSD) under load was used to characterize the evolution of geometrically necessary dislocation (GND) densities at individual grain boundaries as a function of applied strain in a polycrystalline Mg—4Al alloy. The increase in GND density was investigated at plastic strains of 0 %, 0.6 %, 2.2 %, 3.3 % from the area including 76 grains and correlated with (i) geometric compatibility between slip systems across grain boundaries, and (ii) plastic incompatibility. We develop expressions for the grain boundary GND density evolution as a function of plastic strain and plastic incompatibility, from which uniaxial tensile stress-strain response of polycrystalline Mg—4Al are computed and compared with experimental measurement. The findings in this study contribute to understanding the mechanisms governing the strain hardening response of single-phase polycrystalline alloys and more reliable prediction of mechanical behaviors in diverse microstructures.
We introduced nonlocal elastic constitutive laws for crystals containing defects such as dislocations and disclinations. Additionally, the pointwise elastic moduli tensors adequately reflect the elastic response of defect-free regions by relating stresses to strains and couple-stresses to curvatures, elastic cross-moduli tensors relating strains to couple-stresses and curvatures to stresses within convolution integrals are derived from a nonlocal analysis of strains and curvatures in the defects cores. Sufficient conditions are derived for positive-definiteness of the resulting free energy, and stability of elastic solutions is ensured. The elastic stress/couple stress fields associated with prescribed dislocation/disclination density distributions and solving the momentum and moment of momentum balance equations in periodic media are determined by using a Fast Fourier Transform spectral method. Here, the convoluted cross-moduli bring the following results: (i) Nonlocal stresses and couple stresses oppose their local counterparts in the defects core regions, playing the role of restoring forces and possibly ensuring spatio-temporal stability of the simulated defects, (ii) The couple stress fields are strongly affected by nonlocality. Such effects favor the stability of the simulated grain boundaries and allow investigating their elastic interactions with extrinsic defects, (iii) Driving forces inducing grain growth or refinement derive from the self-stress and couple stress fields of grain boundaries in nanocrystalline configurations.