Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CRYSTAL MODELS”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Crystal Plasticity Model of Reactor Pressure Vessel Embrittlement in GRIZZLY

The integrity of reactor pressure vessels (RPVs) is of utmost importance to ensure safe operation of nuclear reactors under extended lifetime. Microstructure-scale models at various length and time scales, coupled concurrently or through homogenization methods, can play a crucial role in understanding and quantifying irradiation-induced defect production, growth and their influence on mechanical behavior of RPV steels. A multi-scale approach, involving atomistic, meso- and engineering-scale models, is currently being pursued within the GRIZZLY project to understand and quantify irradiation-induced embrittlement of RPV steels. Within this framework, a dislocation-density based crystal plasticity model has been developed in GRIZZLY that captures the effect of irradiation-induced defects on the flow stress behavior and is presented in this report. The present formulation accounts for the interaction between self-interstitial loops and matrix dislocations. The model predictions have been validated with experiments and dislocation dynamics simulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Crystal Plasticity Model Validation Using Combined High-Energy Diffraction Microscopy Data for a Ti-7Al Specimen

High-Energy Diffraction Microscopy (HEDM) is a 3-d x-ray characterization method that is uniquely suited to measuring the evolving micromechanical state and microstructure of polycrystalline materials during in situ processing. The near-field and far-field configurations provide complementary information; orientation maps computed from the near-field measurements provide grain morphologies, while the high angular resolution of the far-field measurements provide intergranular strain tensors. The ability to measure these data during deformation in situ makes HEDM an ideal tool for validating micro-mechanical deformation models that make their predictions at the scale of individual grains. Crystal Plasticity Finite Element Models (CPFEM) are one such class of micro-mechanical models. While there have been extensive studies validating homogenized CPFEM response at a macroscopic level, a lack of detailed data measured at the level of the microstructure has hindered more stringent model validation efforts. Here, we utilize an HEDM dataset from an alphatitanium alloy (Ti-7Al), collected at the Advanced Photon Source, Argonne National Laboratory, under in situ tensile deformation. The initial microstructure of the central slab of the gage section, measured via near-field HEDM, is used to inform a CPFEM model. The predicted intergranular stresses for 39 internal grains are then directly compared to data from 4 far-field measurements taken between ~4% and ~80% of the macroscopic yield strength. In conclusion, the intergranular stresses from the CPFEM model and far-field HEDM measurements up to incipient yield are shown to be in good agreement, and implications for application of such an integrated computational/experimental approach to phenomena such as fatigue and crack propagation is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A crystal plasticity model incorporating the effects of precipitates in superalloys: Application to tensile, compressive, and cyclic deformation of Inconel 718

An elasto-plastic polycrystal plasticity model is developed and applied to an Inconel 718 (IN718) superalloy that was produced by additive manufacturing (AM). The model takes into account the contributions of solid solution, precipitates shearing, and grain size and shape effects into the initial slip resistance. Non-Schmid effects and backstress are also included in the crystal plasticity model for activating slip. The hardening law for the critical resolved shear stress is based on the evolution of dislocation density. In using the same set of material and physical parameters, the model is compared against a suite of compression, tension, and large-strain cyclic mechanical test data applied in different AM build directions. We demonstrate that the model is capable of predicting the particularities of both monotonic and cyclic deformation to large strains of the alloy, including decreasing hardening rate during monotonic loading, the non-linear unloading upon the load reversal, the Bauschinger effect, the hardening rate change during loading in the reverse direction as well as plastic anisotropy and the concomitant microstructure evolution. It is anticipated that the general model developed here can be applied to other multiphase alloys containing precipitates.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development of a MOOSE-based crystal plasticity model with irradiation defect evolution for irradiation creep in 316

Irradiation creep and irradiation swelling are two of the lifetime limiting factors for structural materials in nuclear reactors. These mechanical effects are driven by irradiation defect evolution and the interaction of those defects with dislocations in the microstructure. We present here a coupled cluster dynamics – crystal plasticity approach to model irradiation swelling and creep behavior in 316 SS. The time-dependent evolution of irradiation defects is calculated with a cluster dynamics approach and passed to the crystal plasticity model to compute the dislocation evolution. We show the impact of the irradiation defect evolution on the stress state in the material, which drives inelastic deformation through dislocation motion. The inelastic deformation in the 316 SS is dependent on the dose rate, where the inelastic deformation driven by the early-stage irradiation defect evolution determines the mechanical behavior of the 316 SS.

316 Stainless Steels↗
DOE OSTI2021

Crystal plasticity model of residual stress in additive manufacturing using the element elimination and reactivation method

Selective laser melting is receiving increasing interest as an additive manufacturing technique. Residual stresses induced by the large temperature gradients and inhomogeneous cooling process can favour the generation of cracks. In this work, a crystal plasticity finite element model is developed to simulate the formation of residual stresses and to understand the correlation between plastic deformation, grain orientation and residual stresses in the additive manufacturing process. The temperature profile and grain structure from thermal-fluid flow and grain growth simulations are implemented into the crystal plasticity model. An element elimination and reactivation method is proposed to model the melting and solidification and to reinitialize state variables, such as the plastic deformation, in the reactivated elements. The accuracy of this method is judged against previous method based on the stiffness degradation of liquid regions by comparing the plastic deformation as a function of time induced by thermal stresses. The method is used to investigate residual stresses parallel and perpendicular to the laser scan direction, and the correlation with the maximum Schmid factor of the grains along those directions. The magnitude of the residual stress can be predicted as a function of the depth, grain orientation and position with respect to the molten pool. The simulation results are directly comparable to X-ray diffraction experiments and stress–strain curves.

316 stainless steel↗
DOE OSTI2022

Critical Shape for the Growth of Grain Boundary Twin Embryos in Mg and Mg Alloys: Crystal Plasticity Modeling

Application of polycrystalline hexagonal close packed (HCP) metals in engineering designs has been constrained by their anisotropic responses due to twinning and limited plasticity. In deformation, twins most often initiate at grain boundaries (GBs), and thicken and propagate across the grain. In this work, the GB twin embryos in Mg and Mg alloys, and the conditions that influence their propagation are investigated. Using a micromechanical crystal plasticity model, the role of embryo shape on the driving forces prevailing at the embryo boundaries that could support its expansion is studied. The modeled embryos are either planar, extending more in the shear direction than normal to the twin plane, or equiaxed. Results show that the thinner the embryo, the greater the driving forces for both thickening and forward propagation. Alloys with low prismatic-to-basal critical resolved shear stress (CRSS) ratios promote embryo thickening and large CRSS values for the slip mode that primarily accommodates the twin shear encourage propagation. The neighboring grains with orientations that enable local accommodation of the embryo twin shear by pyramidal slip promote forward propagation but have little effect on thickening. When two like embryos lie along the same GB, their paired interaction promotes forward propagation but hinders thickening.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Investigating the Difference in Mechanical Stability of Retained Austenite in Bainitic and Martensitic High-Carbon Bearing Steels using in situ Neutron Diffraction and Crystal Plasticity Modeling

In situ neutron diffraction of the uniaxial tension test was used to study the effect of the surrounding matrix microstructure on the mechanical stability of retained austenite in high-carbon bearing steels. Comparing the samples with bainitic microstructures to those with martensitic ones, it was found that the retained austenite in a bainitic matrix starts transforming into martensite at a lower strain compared to that within a martensitic matrix. On the other hand, the rate of transformation of the austenite was found to be higher within a martensitic microstructure. Crystal plasticity modeling was used to analyze the transformation phenomenon in these two microstructures and determine the effect of the surrounding microstructure on elastic, plastic, and transformation components of the strain. The results showed that the predominant difference in the deformation accumulated was from the transformation strain and the critical transformation driving force within the two microstructures. The retained austenite was more stable for identical loading conditions in case of martensitic matrix compared to the bainitic one. It was also observed that the initial volume fraction of retained austenite within the bainitic matrix would alter the onset of transformation to martensite, but not the rate of transformation.

bainite↗
DOE OSTI2025

Understanding the limitations and potential of micro tensile testing of tungsten and needs for crystal plasticity modelling

Tungsten is the leading plasma facing material candidate due to its exceptional properties. Understanding the response to neutron irradiation is crucial for the lifetime evaluation of tungsten. The limited space in nuclear reactors and the high levels of radioactivity of the specimens after irradiation are significant barriers to accurately measuring the mechanical properties after neutron irradiation. Testing of micro tensile specimens is one approach to reduce the total amount of irradiated material needed for a set of mechanical testing experiments. However, micro samples are not necessarily measuring the bulk material properties, as size effects produce higher measured mechanical properties than is observed in engineering size counterparts. Determining the minimal sample size for reliable bulk property measurement, combined with modeling efforts, is essential. We conducted room temperature tensile tests on tungsten micro tensile specimens fabricated with focused ion beam, plasma focused ion beam and femto-second laser ablation system to dimension of 2×2×7?µm³, 7×7×18?µm³, and 80×100×233?µm³ (width×thickness×gauge-length), respectively. Elevated tensile testing was performed up to 475°C on 5×5×18?µm³ specimens fabricated with focused ion beam. The smallest specimens exhibited a high degree of ductility and strength, whereas the largest specimens demonstrated behavior akin to bulk tungsten. To investigate the effects introduced by the micro specimen fabrication processes and to obtain the necessary bulk dislocation density for the crystal plasticity model, we employed X-ray diffraction depth profiling. This measurement was performed using different X-ray sources with varying penetration depths. Finally, the potential and limitations of micro mechanical tests for tungsten will be discussed.

36 - MATERIALS SCIENCE↗
DOE OSTI2021

A stochastic solver based on the residence time algorithm for crystal plasticity models

Abstract The deformation of crystalline materials by dislocation motion takes place in discrete amounts determined by the Burgers vector. Dislocations may move individually or in bundles, potentially giving rise to intermittent slip. This confers plastic deformation with a certain degree of variability that can be interpreted as being caused by stochastic fluctuations in dislocation behavior. However, crystal plasticity (CP) models are almost always formulated in a continuum sense, assuming that fluctuations average out over large material volumes and/or cancel out due to multi-slip contributions. Nevertheless, plastic fluctuations are known to be important in confined volumes at or below the micron scale, at high temperatures, and under low strain rate/stress deformation conditions. Here, we develop a stochastic solver for CP models based on the residence-time algorithm that naturally captures plastic fluctuations by sampling among the set of active slip systems in the crystal. The method solves the evolution equations of explicit CP formulations, which are recast as stochastic ordinary differential equations and integrated discretely in time. The stochastic CP model is numerically stable by design and naturally breaks the symmetry of plastic slip by sampling among the active plastic shear rates with the correct probability. This can lead to phenomena such as intermittent slip or plastic localization without adding external symmetry-breaking operations to the model. The method is applied to body-centered cubic tungsten single crystals under a variety of temperatures, loading orientations, and imposed strain rates.

42 ENGINEERING↗
NASA NTRS2019

Analysis of Experimental Ice Accretion Data and Assessment of a Thermodynamic Model During Ice Crystal Icing

This paper evaluates a thermodynamic ice crystal icing model that has been previously presented to describe the possible mechanisms of icing within the core of a turbofan jet engine. The model functions between two distinct ice accretions based on a surface energy balance: freeze-dominated icing and melt-dominated icing. Freeze-dominated icing occurs when liquid water (from melted ice crystals) freezes and accretes on a surface along with the existing ice of the impinging water and ice mass. This freeze-dominated icing is characterized as having strong adhesion to the surface. The amount of ice accretion is partially dictated by a freeze fraction, which is the fraction of impinging liquid water that freezes. Melt-dominated icing occurs as unmelted ice on a surface accumulates. This melt-dominated icing is characterized by weakly bonded surface adhesion. The amount of ice accumulation is partially dictated by a melt fraction, which is the fraction of impinging ice crystals that melts. Experimentally observed ice growth rates suggest that only a small fraction of the impinging ice remains on the surface, implying a mass loss mechanism such as splash, runback, bounce, or erosion. The fraction of mass loss must be determined in conjunction with the fraction of freezing liquid water or fraction of melting ice on an icing surface for a given ice growth rate. This mass loss parameter, however, along with the freeze fraction and melt fraction, are the only experimental parameters that are currently not measured directly. Using icing growth rates from ice crystal icing experiments, a methodology that has been previously proposed is used to determine these unknown parameters. This work takes ice accretion data from tests conducted by the National Aeronautics and Space Administration (NASA) at the Glenn Research Center in 2018 that examined the fundamental physics of ice crystal icing. This paper continues evaluation of the thermodynamic model from a previous effort, with additions to the model that account for sub-freezing temperatures that have been observed at the leading edge of the airfoil during icing. The predicted temperatures were generally in good agreement with measured temperatures. Other key findings include the total wet-bulb temperature being a good first order indicator of whether icing is freeze-dominated (sub-freezing values) or melt-dominated (above freezing). Maximum sticking efficiency values, the fraction of impinging mass that adheres to a surface, was calculated to be about 0.2, and retained this maximum value for a range of melt ratios (0.3 to 0.65 and possibly higher), which is defined as the ratio of liquid water content to total water content. Higher air velocities reduced the maximum sticking efficiency and shifted the icing regime to higher melt ratio values. Finally, the leading edge ice accretion angle was found to be related to ice growth (lower growth rates for smaller angles) and melt ratio (smaller melt ratios resulted in smaller angles, likely due to erosion effects).

Ice-Crystal Icing↗
DOE OSTI2016

Quantifying deformation processes near grain boundaries in α titanium using nanoindentation and crystal plasticity modeling

In this work, the influence of grain boundaries on plastic deformation was studied by carrying out nanoindentation near grain boundaries (GBs). Surface topographies of indentations near grain boundaries were characterized using atomic force microscopy (AFM) and compared to corresponding single crystal indent topographies collected from indentations in grain interiors. Comparison of the single crystal indents to indents adjacent to low-angle boundaries shows that these boundaries have limited effect on the size and shape of the indent topography. Higher angle boundaries result in a decrease in the pile-up topography observed in the receiving grain, and in some cases increases in the topographic height in the indented grain, indicating deformation transfer across these boundaries is more difficult. A crystal plasticity finite element (CPFE) model of the indentation geometry was built to simulate both the single crystal and the near grain boundary indentation (bi-crystal indentation) deformation process. The accuracy of the model is evaluated by comparing the point-wise volumetric differences between simulated and experimentally measured topographies. Good agreement, in both single and bi-crystal cases, suggests that the crystal plasticity kinematics plays a dominant role in single crystal indentation deformation, and is also essential to bi-crystal indentation. Despite the good agreement, some differences between experimental and simulated topographies were observed. These discrepancies have been rationalized in terms of reverse plasticity and the inability of the model to capture the full resistance of the boundary to slip. Lastly, this is discussed in terms of dislocation nucleation versus glide in the model and in the physics of the slip transfer process.

36 MATERIALS SCIENCE↗
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 OSTI2014

Measuring the x-ray resolving power of bent potassium acid phthalate diffraction crystals

This report presents the results from measuring the X-ray resolving power of a curved potassium acid phthalate (KAP(001)) spectrometer crystal using two independent methods. It is part of a continuing effort to measure the fundamental diffraction properties of bent crystals that are used to study various characteristics of high temperature plasmas. Bent crystals like KAP(001) do not usually have the same diffraction properties as corresponding flat crystals. Models that do exist to calculate the effect of bending the crystal on the diffraction properties have simplifying assumptions and their accuracy limits have not been adequately determined. The type of crystals that we measured is being used in a spectrometer on the Z machine at Sandia National Laboratories in Albuquerque, New Mexico. The first technique for measuring the crystal resolving power measures the X-ray spectral line width of the characteristic lines from several metal anodes. The second method uses a diode X-ray source and a double crystal diffractometer arrangement to measure the reflectivity curve of the KAP(001) crystal. The width of that curve is inversely proportional to the crystal resolving power. The measurement results are analyzed and discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1993

Orthopyroxenes as recorders of diogenite petrogenesis: Trace element systematics

Eucrite, howardite, and diogenite members of the achondrites are considered, by many, to be genetically related. Therefore, each provide a piece of the puzzle for reconstructing magmatic processes on the eucrite parent body (EPB). The relationship between eucrites and diogenites can be viewed within the context of two distinctly different models: (1) fractional crystallization; and (2) partial melting. In fractional crystallization models, eucrites and diogenites represent a complementary continuum of planetary fractional crystallization products in which the diogenites represent crystal accumulations during the crystallization of eucritic magmas at shallow to deep levels in the EPB. Alternatively, experimental studies may be interpreted as indicating eucrites represent peritectic partial melts of a primitive, chondritic EPB mantle. Within this type of model, the diogenites are also generally considered to be cumulates; however, their petrogenetic relationship to the eucrites is less clear. Sack et al. proposed that the olivine diogenites represent residua from the partial melting events that produced eucritic liquids. Initial trace element studies of orthopyroxene (OPX) are consistent with this model. However, this trace element modeling of the olivine diogenites is nonunique. As a further test of these models, we did the following three things: (1) analyzed OPX from cumulate diogenites to compare with the olibine diogenite data; (2) improved ion microprobe analytical techniques for the analysis of elements critical to our interpretations; and (3) selected more relevant Kd's for OPX-eucritic melt.

Shearer, C. K.↗
DOE OSTI2017

Crystal plasticity modeling of irradiation growth in Zircaloy-2

A reaction-diffusion based mean field rate theory model is implemented in the viscoplastic self-consistent (VPSC) crystal plasticity framework to simulate irradiation growth in hcp Zr and its alloys. A novel scheme is proposed to model the evolution (both number density and radius) of irradiation-induced dislocation loops that can be informed directly from experimental data of dislocation density evolution during irradiation. This framework is used to predict the irradiation growth behavior of cold-worked Zircaloy-2 and trends compared to available experimental data. The role of internal stresses in inducing irradiation creep is discussed. Effects of grain size, texture, and external stress on the coupled irradiation growth and creep behavior are also studied.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Utilizing integrated neutron diffraction and elastoplastic self-consistent crystal plasticity model to quantitatively assess the strengthening mechanism in Al–12.5Ce and Al–12.5Ce–0.4Mg alloys

An integrated in-situ neutron diffraction and elastic plastic self-consistent crystal plasticity (EPSC-CP) modeling scheme is performed on a binary Al–12Ce alloy and a ternary Al–12Ce–0.4Mg alloys. Using this scheme, the constitutive parameters, i.e. elastic constants and slip system parameters of individual phases can be calibrated which can be used in microstructure-based CP models to predict materials performance. From this study, it is shown that the elastic constants of Al 11 Ce 3 intermetallics calculated from density function theory calculation in the literature are rather accurate. When applied to the EPSC-CP model, the lattice strains of both the binary and ternary alloys are correctly predicted as compared with experiments, and large lattice strain differences between Al (100) plane and Al 11 Ce 3 (010) plane are demonstrated. The slip system parameters calibrated by the scheme shows that the addition of 0.4 wt% Mg in the alloy has little influence on the critical resolved shear stress of initial dislocation glide in the Al matrix which caused plastic yield in the material. This can be explained by the very dilute Mg solute content in the Al solid solution, causing large spacing of Al–Mg lattice misfit sites and little impact on resistance of dislocation glide at initial yield. The 0.4 wt% Mg addition, on the other hand, has a large influence on the hardening term in the slip system parameters, indicating those Al–Mg misfit sites do help dislocation accumulation during the deformation. The impact of dilute Mg addition on the Al slip system parameters is also reflected in the flow behavior of the ternary alloy: little impact on the yield stress, but a large impact on working hardening and tensile strength of the materials which is consistent with the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Exact equations of state for nematics

We propose a novel approach to the solution of nematic Liquid Crystal models based on the derivation of a system of nonlinear wave equations for order parameters such that the occurrence of uniaxial and biaxial phase transitions can be interpreted as the propagation of a two-dimensional shock wave in the space of thermodynamic parameters. We obtain the exact equations of state for an integrable model of biaxial nematic liquid crystals and show that the classical transition from isotropic to uniaxial phase in absence of external fields is the result of a van der Waals type phase transition, where the jump in the order parameters is a classical shock generated from a gradient catastrophe at a non-zero isotropic field. The study of the equations of state provides the first analytical description of the rich structure of nematics phase diagrams in presence of external fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

Modeling the Crystallization of Martian Igneous Compositions

The martian crust is dominated by basalts and ultramafic rocks [1]. Understanding the crystallization sequence of martian igneous rocks is critical for constraining the petrologic and geochemical evolution of the crust and mantle. Experimental studies have investigated the parental liquid compositions and crystallization behavior using martian meteorite compositions [2-7], as well as compositions measured by rovers on the martian surface [8,9]. While experiments deliver crucial results, they are often limited by time and resources. Therefore, modeling approaches are an ideal tool to complement experimental results. The application of the MELTS software package [10, 11], to calculate the crystallization of martian compositions, reproduces important compositional details of magma crystallization but significant uncertainties regarding the pressure of the olivine-pyroxene multiple saturation point have been reported [12], and in the temperatures of initial crystallization [12, 13]. Recently, a new melting model, MAGMARS, has been shown to reproduce martian primary melts more accurately [14]. Here, we constrain the applicability of Gibbs free energy minimization software Perple_X [15] for martian compositions by comparing our modeled results to experimentally derived crystallization paths for two martian compositions.

J Semprich↗