Search NASA⌕ Search

SEARCH · Search NASA

Results for “strain 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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Meso-scale modeling and damage analysis of carbon/epoxy woven fabric composite under in-plane tension and compression loadings

The mechanical properties and damage behaviors of carbon/epoxy woven fabric composite under in-plane tension and compression are studied at the meso-scale level through experiment and simulation. An efficient representative volume element (RVE) modeling method with consistent mesh, high yarn volume fraction and realistic geometry is proposed. The material constitutive laws with plasticity, tension-compression asymmetry and damage evolution are established for the three components - yarn, matrix and interface, respectively. Significantly different mechanical properties and damage evolutions are observed depending on loading conditions and initial geometry characteristics. It shows a non-linear stress-strain curve with clear transition region and intensive damage in tension, while a quasi-linear behavior up to facture is observed in compression with little damage prior to final fracture. Furthermore, compared to the constant Poisson's ratio with straining in compression, a dramatic increase in Poisson's ratio appears in tension. Simulation shows damage mechanisms including transverse damage, matrix damage and delamination, which all play critical roles in the property evolution. In particular, the rapid damage accumulation after elastic deformation destroys the strong bonds and causes the easy deformation of transverse yarns which results in the transition region and large Poisson's ratio in tension. All the mechanical behaviors and damage evolutions are well captured and explained with the current RVE model.

42 ENGINEERING↗
DOE OSTI2020

Deciphering atomistic mechanisms of the gas-solid interfacial reaction during alloy oxidation

Gas-solid interfacial reaction is critical to many technological applications from heterogeneous catalysis to stress corrosion cracking. A prominent question that remains unclear is how gas and solid interact beyond chemisorption to form a stable interphase for bridging subsequent gas-solid reactions. Here, we report real-time atomic-scale observations of Ni-Al alloy oxidation reaction from initial surface adsorption to interfacial reaction into the bulk. We found distinct atomistic mechanisms for oxide growth in O 2 and H 2 O vapor, featuring a “step-edge” mechanism with severe interfacial strain in O 2 , and a “subsurface” one in H 2 O. Ab initio density functional theory simulations rationalize the H 2 O dissociation to favor the formation of a disordered oxide, which promotes ion diffusion to the oxide-metal interface and leads to an eased interfacial strain, therefore enhancing inward oxidation. Our findings depict a complete pathway for the Ni-Al surface oxidation reaction and delineate the delicate coupling of chemomechanical effect on gas-solid interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Experimental Investigation of Elastodynamic Nonlinear Response of Dry Intact, Fractured and Saturated Rock

Nonlinear elastodynamic response of fractured rocks carries crucial information on fracture features that can be exploited to forecast flow properties, friction constitutive behavior and poromechanical response. Well-controlled laboratory experiments are designed to measure the nonlinear elastodynamic response of Westerly granite in three states: dry intact, dry fractured and saturated fractured. We study the effect of fracturing and saturation in modifying the elastodynamic response of the rock. Each sample is tested at a normal stress level of 15 MPa. We measure the elastodynamic response of an intact L-shaped sample of Westerly Granite subjected to normal stress oscillations of prescribed amplitudes (0.2–1 MPa) and frequencies (0.1, 1, 10 Hz). Ultrasonic waves transmitted across the sample are used to monitor the evolution of wave velocity before, during and after dynamic stressing. Furthermore, the nonlinearity of the elastodynamic response is measured in terms of: (1) the offset in normalized wave velocity; (2) the amplitude of wave velocity fluctuation during the oscillations; and (3) recovery rate of the wave velocity post-oscillation. We observe that the three nonlinearity parameters show a similar trend. Irrespective of the parameter, the nonlinearity measures higher for sample in dry-intact condition than that for dry-fractured and the saturated-fractured sample exhibits smaller nonlinearity than the dry-fractured sample. As expected, the saturated sample exhibits less nonlinearity than the dry intact and fractured samples due to the presence of interstitial fluid and the resulting increased interface stiffness. Conversely, the dry intact rock shows a higher nonlinearity than the dry fractured. We use numerical simulations to show that the presence of fracture significantly alters the strain distribution across the bulk of the sample and only the contacting asperities are highly strained, thus resulting in a decrease in the measured elastodynamic nonlinearity.

58 GEOSCIENCES↗
DOE OSTI2021

Tool-workpiece stick-slip conditions and their effects on torque and heat generation rate in the friction stir welding

Friction stir welding (FSW) has found increased applications in automotive and aerospace industries due to its advantages of solid-state bonding, no fusion and melting, and versatility in various working conditions and material combinations. However, the relationship among processing parameters, material properties, and bonding extent and fidelity remains largely empirical, primarily because of the lack of the mechanistic understanding of the tool-workpiece frictional behavior that affects our subsequent understanding of microstructural evolution and interface bonding formation. While the tool-workpiece stick-slip condition is believed to dictate the resulting torque and heat generation rate during the welding process, it remains rare and elusive to conduct a quantitative experimental measurement of such interfacial field. On the other hand, numerical simulations based on Computational Fluid Dynamics (CFD) rely on ad hoc assumptions of interfacial pressure and shear-stress conditions, but predictions can only be validated via the medium- and far-range temperature field which is known to be insensitive to the interfacial frictional behavior. This work first presents a comparison among two CFD-based simulation methodologies and the Coupled Eulerian Lagrangian (CEL) model in finite element method, the last of which uses the Coulomb friction so that the stick-slip is naturally developed. Based on the Hill-Bower similarity relationship in the contact analysis, an analytical model is developed here to prove why a constant stick-slip fraction will be developed in the steady state, to correlate the stick-slip fraction to processing parameters such as the tool spin rate, and further to derive dimensionless functions for torque and heat-generation-rate predictions. Pros and cons of various numerical approaches in predicting stick-slip are discussed, and our analytical model has been found to agree well with our numerical simulation and literature experimental results. These analyses provide the critical strain-rate and temperature fields that are needed for the bonding analysis in our future work.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Comparative analysis of plasticity-based GND density estimation methods in crystal plasticity finite element models

In crystal plasticity finite element (CPFE) simulations, accurately quantifying geometrically necessary dislocations (GNDs) is critical for capturing strain gradients in polycrystals. We compare different methods for quantifying GNDs, all of which originate from the Nye tensor, which is computed as the curl of the plastic deformation gradient. The projection technique directly decomposes the Nye tensor onto individual screw and edge dislocation components to compute GNDs. This approach requires converting a nine-component Nye tensor into densities for a larger number of dislocation systems, a fundamentally underdetermined (non-unique) process, which is resolved using L2 minimization. In contrast, when employing CPFE analysis, one could directly compute dislocation densities on each slip system using shear gradients. Projection and slip gradient methods are compared with respect to their prediction of GNDs with changing grain size, strain, and grain neighborhoods, including multigrain junctions. Although these techniques match analytical GND densities for single slip, single crystal deformation, and are consistent with anticipated overall GND trends, we find that the GND densities from projection techniques are significantly lower than those predicted from CPFE-based slip gradients in polycrystals. A suggested improvement of only using the active dislocation systems in the projection technique almost entirely resolved this mismatch.

Crystal plasticity↗
DOE OSTI2021

Simulations of crack extensions in arc-shaped tension specimens of uncharged and hydrogen-charged 21-6-9 austenitic stainless steels using cohesive zone modeling with varying cohesive parameters

Crack extensions in uncharged and hydrogen-charged side-grooved A(T) specimens of conventionally forged 21-6-9 austenitic stainless steels are simulated using the cohesive zone modeling approach. Two-dimensional plane strain finite element analyses with fixed cohesive parameters are first conducted to fit the experimental load-displacement curves. Similar analyses using varying cohesive parameters as functions of the crack extension are then conducted to fit the experimental load-crack extension and crack extension-displacement curves. The computational results with varying cohesive parameters can fit very well the experimental data. The computational results also indicate that the average cohesive energy for the hydrogen-charged A(T) specimen is lower than that for the uncharged A(T) specimen.

42 ENGINEERING↗
DOE OSTI2020

Shear Softening in a Metallic Glass: First-Principles Local-Stress Analysis

Metallic glasses deform elastically under stress. However, the atomic-level origin of elastic properties of metallic glasses remain unclear. In this Letter using ab initio molecular dynamics simulations of the Cu 5 0Zr 50 metallic glass under shear strain, we show that the heterogeneous stress relaxation results in the increased charge transfer from Zr to Cu atoms, enhancing the softening of the shear modulus. Changes in compositional short-range order and atomic position shifts due to the nonaffine deformation are discussed. It is shown that the Zr subsystem exhibits a stiff behavior, whereas the displacements of Cu atoms from their initial positions, induced by the strain, provide the stress drop and softening.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Physical Regime Sensitivity

Abstract This work presents a novel sensitivity approach that quantifies sensitivity to regimes of a model’s state variables rather than constitutive model parameters. This Physical Regime Sensitivity (PRS) determines which regimes of a model’s independent variables have the biggest influence on an experiment or application. PRS analysis is demonstrated on a strength model used in the simulation of a copper Taylor cylinder. In a series of simulations, the strength model was perturbed sequentially in local regimes of plastic strain, plastic strain rate, temperature and pressure, and then the prediction of cylinder shape was compared to unperturbed calculations. Results show, for example, that the deformed length of the cylinder was most sensitive to strength at a strain rate of 1.0 × 10 4 /sec., but the deformed footprint radius was most sensitive to strength at a strain rate of about 4.0 × 10 4 /sec. Compared to current histogram approaches, PRS can be used to design or interpret integrated experiments by identifying not just which regimes are accessed somewhere in the experiment but the causality question of which regimes actually affect the measured data. PRS should allow one to focus experimental and modeling efforts where they are most needed and to better interpret experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Neural networks for large eddy simulations of wall-bounded turbulence: numerical experiments and challenges

Here, we examine the application of neural network-based methods to improve the accuracy of large eddy simulations of incompressible turbulent flows. The networks are trained to learn a mapping between flow features and the subgrid scales, and applied locally and instantaneously—in the same way as traditional physics-based subgrid closures. Models that use only the local resolved strain rate are poorly correlated with the actual subgrid forces obtained from filtering direct numerical simulation data. We see that highly accurate models in a priori testing are inaccurate in forward calculations, owing to the preponderance of numerical errors in implicitly filtered large eddy simulations. A network that accounts for the discretization errors is trained and found to be unstable in a posteriori testing. We identify a number of challenges that the approach faces, including a distribution shift that affects networks that fail to account for numerical errors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Quantifying the Agyrotropy of Proton and Electron Heating in Turbulent Plasmas

An important aspect of energy dissipation in weakly collisional plasmas is that of energy partitioning between different species (e.g., protons and electrons) and between different energy channels. Here we analyse pressure–strain interaction to quantify the fractions of isotropic compressive, gyrotropic, and nongyrotropic heating for each species. An analysis of kinetic turbulence simulations is compared and contrasted with corresponding observational results from Magnetospheric Multiscale Mission data in the magnetosheath. In assessing how protons and electrons respond to different ingredients of the pressure–strain interaction, we find that compressive heating is stronger than incompressive heating in the magnetosheath for both electrons and protons, while incompressive heating is stronger in kinetic plasma turbulence simulations. Concerning incompressive heating, the gyrotropic contribution for electrons is dominant over the nongyrotropic contribution, while for protons nongyrotropic heating is enhanced in both simulations and observations. Variations with plasma β are also discussed, and protons tend to gain more heating with increasing β.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Off-Stoichiometric Restructuring and Sliding Dynamics of Hexagonal Boron Nitride Edges in Conditions of Oxidative Dehydrogenation of Propane

Boron-containing materials, such as hexagonal boron nitride (h-BN), recently shown to be active and selective catalysts for the oxidative dehydrogenation of propane (ODHP), have been shown to undergo significant surface oxyfunctionalization and restructuring. Although experimental ex situ studies have probed the change in chemical environment on the surface, the structural evolution of it under varying reaction conditions has not been established. Herein, we perform global optimization structure search with a grand canonical genetic algorithm to explore the chemical space of off-stoichiometric restructuring of the h-BN surface under ambient as well as ODHP-relevant conditions. A grand canonical ensemble representation of the surface is established, and the predicted 11B solid-state NMR spectra are consistent with previous experimental reports. In addition, we investigated the relative sliding of h-BN sheets and how it influences the surface chemistry with ab initio molecular dynamics simulations. Furthermore, the B–O linkages on the edges are found to be significantly strained during the sliding, causing the metastable sliding configurations to have higher reactivity toward the activation of propane and water.

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

Strain phase equilibria and phase‐field method of ferroelectric polydomain: A case study of monoclinic K x Na 1 − x NbO 3 thin films

Abstract Knowledge of the thermodynamic equilibria and domain structures of ferroelectrics is critical to establishing their structure–property relationships that underpin their applications from piezoelectric devices to nonlinear optics. Here, we establish the strain condition for strain phase separation and polydomain formation and analytically predict the corresponding domain volume fractions and wall orientations of, relatively low symmetry and theoretically more challenging, monoclinic ferroelectric thin films by integrating thermodynamics of ferroelectrics, strain phase equilibria theory, microelasticity, and phase‐field method. Using monoclinic K x Na 1 − x NbO 3 (0.5 < x < 1.0) thin films as a model system, we establish the polydomain strain–strain phase diagrams, from which we identify two types of monoclinic polydomain structures. The analytically predicted strain conditions of formation, domain volume fractions, and domain wall orientations for the two polydomain structures are consistent with phase‐field simulations and in good agreement with experimental results in the literature. The present study demonstrates a general, powerful analytical theoretical framework to predict the strain phase equilibria and domain wall orientations of polydomain structures applicable to both high‐ and low‐symmetry ferroelectrics and provide fundamental insights into the equilibrium domain structures of ferroelectric K x Na 1 − x NbO 3 thin films that are of technology relevance for lead‐free dielectric and piezoelectric applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Direct Imaging of Hydrogen‐Driven Dislocation and Strain Field Evolution in a Stainless Steel Grain

Hydrogen embrittlement (HE) poses a significant challenge to the durability of materials used in hydrogen production and utilization. Disentangling the competing nanoscale mechanisms driving HE often relies on simulations and electron-transparent sample techniques, limiting experimental insights into hydrogen-induced dislocation behavior in bulk materials. This study employs in situ Bragg coherent X-ray diffraction imaging to track three-dimensional (3D) dislocation and strain field evolution during hydrogen charging in a bulk grain of austenitic 316 stainless steel. Tracking a single dislocation reveals hydrogen-enhanced mobility and relaxation, consistent with dislocation dynamics simulations. Subsequent observations reveal dislocation unpinning and climb processes, likely driven by osmotic forces. Additionally, nanoscale strain analysis around the dislocation core directly measures hydrogen-induced elastic shielding. These findings experimentally validate theoretical predictions and offer mechanistic insights into hydrogen-driven dislocation behavior. The quantified nanoscale phenomena serve as critical inputs for multiscale modeling frameworks to predict bulk material responses and accelerate the development of HE-resistant alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Dislocation density-based plasticity model from massive discrete dislocation dynamics database

In this work, we present a dislocation density-based strain hardening model for single crystal copper through a systematic coarse-graining analysis of more than 200 discrete dislocation dynamics (DDD) simulations of plastic deformation under uniaxial tension. The proposed constitutive model has two components: a generalized Taylor relation connecting resolved shear stresses to dislocation densities on individual slip systems, and a generalized Kocks-Mecking model for dislocation multiplication. The DDD data strongly suggests a logarithmic dependence of flow stress on the plastic shear strain rate on each slip system, and, equivalently, an exponential dependence of the plastic shear strain rate on the resolved shear stress. Hence the proposed generalized Taylor relation subsumes the Orowan relation for plastic flow. The DDD data also calls for a correction to the Kocks-Mecking model of dislocation multiplication to account for the increase of dislocation density on slip systems with negligible plastic shear strain rate. This is accomplished by allowing the multiplication rate on each slip system to include contributions from the plastic strain rates of the two coplanar slip systems. The resulting constitutive model successfully captures the strain hardening rate dependence on the loading orientation as predicted by the DDD simulations, which is also consistent with existing experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Crystal plasticity model of BCC metals from large-scale MD simulations

Accurate crystal plasticity models that faithfully capture the behavior of single crystals under a wide range of loading conditions, such as loading direction, strain rate, and temperature, are still lacking. Here we introduce a novel approach in which a crystal plasticity (CP) model is informed directly from and calibrated to large-scale quantum-accurate MD simulations in which single crystal BCC Ta serves as a testbed material. By analyzing our large set of MD simulations several key insights are obtained leading us to modify constitutive assumptions in order to address deficiencies of existing CP models. Importantly, we observe that the standard notion of fixed slip systems – pairs of slip directions and slip planes – is inadequate for describing high-rate plasticity in BCC tantalum at room temperature. Instead, pencil glide defined as dislocation motion in the maximum resolved shear stress planes (MRSSP) of each Burgers vector is fully consistent with our MD simulation data while providing significant simplifications of the constitutive relations. Herein, our resulting new CP model closely matches the behavior of single crystals observed in high-rate MD simulations while being fully consistent with lower rate experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mesh objective stochastic simulations of quasibrittle fracture

Continuum finite element (FE) modeling of damage and failure of quasibrittle structures suffers from the spurious mesh sensitivity due to strain localization. Here this issue has been addressed for deterministic analysis through the development of localization limiters. Here this study proposes a mechanism-based model to mitigate the mesh sensitivity in stochastic FE simulations of quasibrittle fracture. The interest is placed on the analysis of large-size structures, where the mesh size is conveniently chosen to be larger than the width of the fracture process zone as well as the correlation length of the random fields of constitutive properties. The present model is formulated within the framework of continuum damage mechanics. Two localization parameters are introduced to describe the evolution of the damage pattern of each finite element. These parameters are used to guide the energy regularization of the constitutive law, as well as to formulate the mesh-dependent probability distributions of constitutive properties. Depending on the prevailing damage pattern, different energy regularization schemes and mesh dependence of the probability distribution functions are used in the constitutive law. The model is applied to simulate the stochastic failure behavior of quasibrittle structures of different geometries featuring different failure processes including damage initiation, localization, and propagation. It is shown that using fixed probability distribution functions of constitutive properties could lead to strong mesh dependence of the prediction of the mean and variance of the structural load capacity. The probability distribution functions of constitutive properties must be linked to the damage pattern, which may evolve during the failure process. Such a mechanism-based modeling of the probability distributions of constitutive properties is essential for mitigating the spurious mesh sensitivity in stochastic FE analysis of quasibrittle fracture.

42 ENGINEERING↗
DOE OSTI2024

Robustness of magnetocrystalline anisotropy and coercivity in Fe–Co–B

As a potential candidate for critical chemical element free permanent magnet, (Fe 0.7 Co 0.3 ) 2 B, shows promising intrinsic magnetic properties, including large magnetization, high Curie temperature and moderate magnetocrystalline anisotropy (MCA). Still, the coercivity that has been achieved is relatively low. To understand and improve coercivity, here we study the robustness of MCA related to the fluctuation of chemical composition and the structural distortion and its effect on coercivity in (Fe 0.7 Co 0.3 ) 2 B from first principles DFT (density functional theory) calculation and micromagnetic simulation. A fluctuation of cobalt content of 10%, or 1% isotropic lattice strain, or 2% tetragonal distortion reduces MCA energy up to 30%. By considering randomly distributed local MCA energy reduction up to 30%, the calculated coercivity of anisotropic bulk magnets reduced from 3.6 kOe to 2.2 kOe. Rational microstructure design, such as grain size refinement and development of strong crystal texture of MCA easy-axis (i.e., degree of alignment for (Fe 0.7 Co 0.3 ) 2 B grains), can enhance coercivity and mitigate the effect of local MCA reduction. In addition to acceptable intrinsic magnetic properties, the promising hard magnetic phase needs good MCA robustness to facilitate the development of practical magnets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗