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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Ion-beam radiation-induced Eshelby transformations: The mean and variance in hydrostatic and shear residual stresses

Ion beam plays a pivotal role in ion implantations and the fabrication of nanostructures. However, there lacks a quantitative model to describe the residual stresses associated with the ion-beam radiation. Radiation-induced residual stress/transformation strain have been mostly recognized in the hydrostatic sub strain space. Here, in this study, we use molecular dynamics (MD) simulations to show that the response of a material to irradiation is generally anisotropic that depends on the ion-beam direction, and should be described using tensorial quantities. We demonstrate that accelerator-based ion beam irradiation, combined with the intrinsic lattice anisotropy and externally induced anisotropy (such as anisotropic mechanical loadings), causes radiation-actuated shear transformation strains in addition to hydrostatic expansion. We map out these complex correlations for several materials. Radiation-induced defects are shown to be responsible for residual shear stresses in the manner of Eshelby inclusion transformation. We propose such tensorial response model should be considered for accurate nanoscale fabrication using ion-beam irradiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Nonlinear convergence in contact mechanics: Immersed boundary finite volume

In this report we present an immersed boundary finite volume (IBM) method for simulating quasistatic contact mechanics of linearly elastic domains at small strains. In IBM, all external boundaries and internal contacts of an object are represented by embedded surfaces inside a Cartesian mesh, which need not conform to the grid lines. The contact constraints consist of the non-penetrability condition and Coulomb’s friction law, which are discretized using special interpolation stencils and enforced via penalty parameters. The resulting nonlinear system depends on displacement unknowns only. To solve it, we use the Newton method but find that it diverges frequently. To understand the divergence pattern, we analyze a simplified 2-cell problem and show that the global convergence of Newton cannot be ensured for any choice of penalty parameters. We thus propose a modified Newton solver, which guarantees convergence for the 2-cell problem and is numerically verified to converge for all the challenging simulations considered herein. While both 1 st - and 2 nd -order variants of IBM, in displacement unknowns, are proposed, the modified Newton solver applies only to the 1 st -order variant.

42 ENGINEERING↗
DOE OSTI2022

Development of integrated screening, cultivar optimization, and verification research (DISCOVR): A coordinated research-driven approach to improve microalgal productivity, composition, and culture stability for commercially viable biofuels production

To address major knowledge gaps and barriers to the commercial development of algal biomass for biofuels and co-products, a collaborative consortium, Development of Integrated Screening, Cultivar Optimization, and Verification Research (DISCOVR), was established in 2016. Funded by the U.S. Department of Energy (DOE) Bioenergy Technologies Office (BETO), this consortium constitutes a partnership between four DOE national laboratories - Pacific Northwest National Laboratory (PNNL), Los Alamos National Laboratory (LANL), the National Renewable Energy Laboratory (NREL), and Sandia National Laboratories (SNL) - and the Arizona Center for Algae Technology and Innovation (AzCATI) at Arizona State University. To address the barriers of strain selection for achieving high seasonal productivities with a suitable composition and culture resilience, a tiered strain down-selection pipeline is implemented. At Tier I, the temperature and salinity tolerance of strains is determined in flask cultures; at Tier II, the areal biomass productivity and composition is determined in climate-simulation photobioreactors; at Tier III, the productivity and culture stability are determined in outdoor raceways. The top performing strains move forward to long-term testing at the algae testbed site at AzCATI to generate annual biomass productivity data. Concurrent to the strain down-selection in the DISCOVR pipeline, hypotheses for increasing biomass productivity, shifting biomass composition to enhance intrinsic value, and improving culture stability and resistance to pests are also tested. Techno-economic analyses are carried out to determine whether promising findings from laboratory studies or proposed modifications in outdoor pond cultivation conditions translate into reductions in the minimum biomass selling price (MBSP). Notably, in the three years following the launch of DISCOVR, annual biomass productivity has increased from 11.7 to 17.6 g m -2 day -1 , resulting in an MBSP decrease from 824 to 611 $ ton -1 .

09 BIOMASS FUELS↗
DOE OSTI2023

Microstructure characterization and elastic-plastic self-consistent simulation studies of anisotropic deformation of β -tin

Uniaxial compression tests of 99.9% pure polycrystalline $β$-tin (Sn) were conducted at various strain rates (10 -3 /s, 10 -1 /s and 1/s) and temperatures (294 K and 193 K) to understand the effects of these variables on the stress–strain response, microstructure evolution, and stress relaxation behavior. Multiple Sn specimens were subjected to complex compressive loading/unloading/reloading paths at different strain rates. Specimens initially compressed at higher strain rates showed strain rate-dependent texture evolution and more pronounced relaxation upon reloading, as compared to those pre-strained at lower rates at room temperature (RT). Compression tests conducted at low temperature (193 K) revealed increased strength, similar to the enhanced strength observed when the strain rate was increased at RT. Neutron diffraction was employed to characterize the initial and final bulk textures of the RT specimens. Electron backscatter diffraction was utilized to examine the crystallographic grain orientation and morphology, thereby identifying the signatures of dislocation-mediated deformation, grain refinement, recrystallization, and twinning behaviors. Elastic-plastic self-consistent simulations were performed to investigate the deformation modes responsible for the strain rate-dependent macroscopic stress–strain response and texture evolution. Effects of crystallographic orientation on stress relaxation behavior was also examined. In conclusion, the model predictions are in reasonable agreement with experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A finite viscoelastic constitutive model for low to high strain rate response of elastomers with application of strain rate-induced glass transition

Amorphous elastomers exhibit significant rate-stiffening and unique viscous flow characteristics across a wide range of strain rates, often undergoing glass transition above a strain rate threshold. We have developed a thermodynamically-consistent and micromechanically-inspired constitutive model for soft elastomeric materials to capture the rate-dependent stress-strain behavior and hysteresis when subjected to low to high strain rates. Here, our proposed constitutive model encapsulates the viscous flow of materials through molecular motion at low strain rates and local rearrangement and alignment of the molecules trying to overcome the intermolecular resistance at high strain rates, essentially covering the glass transition. We applied our constitutive model to uniaxial compression experiments performed at low and high strain rates for polyborosiloxane (PBS) to identify the material parameters, and subsequently, performed numerical simulations of single and multi-cycle compression, stress relaxation, and small amplitude oscillatory tension-compression. Our analyses indicate that the model predicts higher total energy dissipation with increasing strain rate; however, dissipation associated with molecular relaxation decreases (forming a cusp) because, beyond a crossover strain rate, intermolecular rearrangement and alignment become dominant, which is consistent with the onset of the glass transition. For cyclic loading-unloading, we observed that dissipation over a cycle remains constant at low strain rates but decreases non-monotonically at high strain rates before becoming constant, with the peak stress over the cycle becoming higher, which can be interpreted as more loading being carried elastically by the polymer network as the intermolecular rearrangement process occurs. Additionally, our model was able to predict the qualitative nature of the storage modulus and loss modulus in the limit of small strain over a wide range of frequency sweeps.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Anomalous temperature dependence of elastic limit in metallic glasses

Abstract Understanding the atomistic mechanisms of inelastic deformation in metallic glasses (MGs) remains challenging due to their amorphous structure, where local carriers of plasticity cannot be easily defined. Using molecular dynamics (MD) simulations, we analyzed the onset of inelastic deformation in CuZr MGs, specifically the temperature dependence of the elastic limit, in terms of localized shear transformation (ST) events. We find that although the ST events initiate at lower strain with increasing temperature, the elastic limit increases with temperature in certain temperature ranges. We explain this anomalous behavior through the framework of an energy-strain landscape (ESL) constructed from high-throughput strain-dependent energy barrier calculations for the ST events identified in the MD simulations. The ESL reveals that the anomalous behavior is caused by the transition of ST events from irreversible to reversible with increasing temperature. An analytical formulation is developed to predict this transition and the temperature dependence of the elastic limit.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Multiscale analysis of large twist ferroelectricity and swirling dislocations in bilayer hexagonal boron nitride

With its atomically thin structure and intrinsic ferroelectric properties, heterodeformed bilayer hexagonal boron nitride (hBN) has gained prominence in next-generation non-volatile memory applications. However, studies to date have focused almost exclusively on small-twist bilayer hBN, leaving the question of whether ferroelectricity can persist under small heterostrain and large heterodeformation entirely unexplored. In this work, we establish the crystallographic origin of ferroelectricity in bilayer hBN configurations heterodeformed relative to high-symmetry configurations such as AA-stacking and 21.786789° twisted configurations (Σ7), using Smith normal form bicrystallography. We then demonstrate out-of-plane ferroelectricity in bilayer hBN across configurations vicinal to both the AA and Σ7 stackings. Atomistic simulations reveal that AA-vicinal systems support ferroelectricity under both small twist and small strain, with polarization switching in the latter governed by the deformation of swirling dislocations rather than the straight interface dislocations seen in the former. For Σ7-vicinal systems, where existing interatomic potentials underperform particularly under extreme out-of-plane compression, we develop a density-functional-theory-informed continuum framework—the bicrystallography-informed frame-invariant multiscale (BFIM) model, which captures out-of-plane ferroelectricity in heterodeformed configurations vicinal to Σ7 stacking. Interface dislocations in these large heterodeformed bilayer configurations exhibit markedly smaller Burgers vectors compared to interface dislocations in small-twist and small-strain bilayer hBN. The BFIM model reproduces experimental results and provides a powerful, computationally efficient framework for predicting ferroelectricity in large-unit-cell heterostructures where atomistic simulations are prohibitively expensive.

Ahmed, Md Tusher [Univ. of Illinois at Urbana-Cham↗
DOE OSTI2023

Ab-initio simulations of coherent phonon-induced pumping of carriers in zirconium pentatelluride

Abstract Laser-driven coherent phonons can act as modulated strain fields and modify the adiabatic ground state topology of quantum materials. Here we use time-dependent first-principles and effective model calculations to simulate the effect of the coherent phonon induced by strong terahertz electric field on electronic carriers in the topological insulator ZrTe 5 . We show that a coherent A 1g Raman mode modulation can effectively pump carriers across the band gap, even though the phonon energy is about an order of magnitude smaller than the equilibrium band gap. We reveal the microscopic mechanism of this effect which occurs via Landau–Zener-Stückelberg tunneling of Bloch electrons in a narrow region in the Brillouin zone center where the transient energy gap closes when the system switches from strong to weak topological insulator. The quantum dynamics simulation results are in excellent agreement with recent pump-probe experiments in ZrTe 5 at low temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Large-scale atomistic studies of sliding friction in polycrystalline aluminum interfaces

In this work, we discuss large-scale non-equilibrium molecular dynamics (NEMD) simulations of ductile metal sliding comprising up to 1.8 × 10 9 atoms over time scales of 100 ns. The results of these simulations have identified a variety of physical mechanisms that are important in determining the steady-state frictional force for a wide range of velocities at compressed metal–metal interfaces. These include grain growth and refinement, the evolution of large plastic strains and strain rates, material mixing, and melting. These phenomena can be included in a strain, strain rate, and grain size model that gives good agreement with the NEMD simulations and can be applied to macroscopic continua.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Heterogeneous solute segregation suppresses strain localization in nanocrystalline Ag-Ni alloys

Solute segregation to individual grain boundaries is used by design to produce strong and stable nanocrystalline metallic alloys. Grain-boundary segregation, however, is known to cause adverse embrittlement effects from a strain-localization failure mechanism that imposes significant material limitations for structural applications. Here, using atomistic simulations, it is discovered that heterogeneous Ni segregation in nanocrystalline Ni-mixed Ag alloys dramatically shuts down localized shear bands during plastic deformation, while simultaneously increasing the tensile strength. Nanocrystalline Cu-mixed Ag metals are predicted to exhibit standard homogeneous Cu segregation and a tensile strength that saturates above a solute concentration of 8 at.% due to glass-like shear localization induced by grain boundaries. By contrast, it is found that heterogeneous Ni segregation in nanocrystalline Ag-Ni alloys forms solute-rich clusters along interfaces leading to strain delocalization at high strain and continuous strengthening at high solute concentrations up to 15 at.%. As a result, this study reveals the importance of heterogeneous versus homogeneous segregation behaviors on strain localization and points to a fundamentally new strategy to design failure-resistant nanostructured materials through grain boundary segregation engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Binding Interactions between Receptor-Binding Domain of Spike Protein and Human Angiotensin Converting Enzyme-2 in Omicron Variant

The emergence of new SARS-CoV-2 Omicron variant of concern (OV) has exacerbated the COVID-19 pandemic because of a large number of mutations in the spike protein, particularly in the receptor-binding domain (RBD), resulting in highly contagious and/or vaccine-resistant strains. Herein, we present a systematic analysis based on detailed molecular dynamics (MD) simulations in order to understand how the OV RBD mutations affect the ACE2 binding. We show that the OV RBD binds to ACE2 more efficiently and tightly predominantly because of strong electrostatic interactions, thereby promoting increased infectivity and transmissibility compared to other strains. Some of the OV RBD mutations are predicted to affect the antibody neutralization either through their role in the S-protein conformational changes, such as S371L, S373P, and S375F, or through changing its surface charge distribution, such as G339D, N440K, T478K, and E484A. Other mutations, such as K417N, G446S, and Y505H, decrease the ACE2 binding, whereas S447N, Q493R, G496S, Q498R, and N501Y tend to increase it.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Defect-Enabled Phase Programming of Transition Metal Dichalcogenide Monolayers

The ability to tune the local electronic transport properties of group VI transition metal dichalcogenide (TMD) monolayers by strain-induced structural phase transformations (“phase programming”) has stimulated much interest in the potential applications of such layers as ultrathin programmable and dynamically switchable nanoelectronics components. In this manuscript, we propose a new approach toward controlling TMD monolayer phases by employing macroscopic in-plane strains to amplify heterogeneous strains arising from tailored, spatially extended defects within the monolayer. Furthermore, the efficacy of our proposed approach is demonstrated via numerical simulations of emerging domains localized around arrays of holes, grain boundaries, and compositional heterointerfaces. Quantitative relations between the macroscopic strains required, spatial resolution of domain patterns, and defect configurations are developed. In particular, the introduction of arrays of holes is identified as the most feasible phase programming route.

2D materials↗
DOE OSTI2025

Mechanically and electrically tunable Rashba-Edelstein effect in ferroelectric semiconductors, CsGe⁢𝑋 3 (𝑋 = I, Br, Cl)

The ability of materials to convert charge current into spin current is fundamental to many spintronics applications. One means of realizing this conversion is via Rashba-Edelstein effect (REE). Using density functional theory simulations, we predict that REE can be induced in the recently discovered family of semiconducting ferroelectrics, CsGe⁢𝑋 3 (𝑋 = I, Br, Cl). The effect is quantified through Rashba-Edelstein coefficients, 𝜒 𝑥⁢𝑦 =−𝜒 𝑦⁢𝑥 , which are nonzero in valence and conduction bands. The largest values, obtained for CsGeI 3 , are 3.45 × 10 10 $\frac{ℏ}{Λ cm}$ and 0.97 × 10 10 $\frac{ℏ}{Λ cm}$ in the conduction and valence bands, respectively. The values are comparable to, and sometimes exceed, those of other inorganic materials, although the maximal values occur away from the band edges. The coefficients' sign couples to the direction of spontaneous polarization, offering opportunities for nonvolatile spin current manipulation via external electric field. Furthermore, these coefficients are highly tunable through strain engineering owing to strain-induced variations in spin textures and energy isosurfaces. Specifically, 𝜒 in the valence band of CsGeI 3 is enhanced to 3.61 × 10 10 $\frac{ℏ}{Λ cm}$ under 5% biaxial strain. Given this potential, we synthesized CsGeI 3 to validate our simulation structure and found excellent agreement between experiment and simulation, thereby allowing extrapolation of our findings to practically significant temperatures. Furthermore, our study identifies promising materials for nonvolatile, multifunctional spintronic applications.

Electronic structure↗
DOE OSTI2020

Modeling the role of local crystallographic correlations in microstructures of Ti-6Al-4V using a correlated structure visco-plastic self-consistent polycrystal plasticity formulation

This paper presents a multi-level crystal plasticity-based simulation framework for modeling mechanical response and microstructure evolution of Ti-6Al-4V with α-lath/lamellar microstructures. Here, the model is a correlated structure visco-plastic self-consistent (CS-VPSC) formulation linking three scales: a single crystals micro-scale, a lath/lamellar colony meso-scale, and a lath/lamellar aggregate macro-scale. A selected hardening law for the evolution of critical resolved shear stress per slip system used in CS-VPSC is phenomenological. However, it adjusts the resistances of basal and prismatic slip systems based on the geometry of slip transfer between adjacent lamellae. Consistent with experimental evidences, the resolved shear stress on the pyramidal slip planes is dependent not only on the stress in the direction of slip but also on the two orthogonal shear stress components and the three normal stress components (non-Schmid effects). Electron backscatter diffraction (EBSD) data in conjunction with a procedure relying on α→β phase transformation is used to construct paired variants of α-laths/lamellae satisfying their local crystallographic correlations. The procedure fits volume fractions of individual laths/lamellae with the experimental EBSD data and selects a distribution of habit planes between adjacent variants with respect to the loading direction. The simulation framework is applied to interpret the deformation behavior in tension and compression along two sample directions of Ti-6Al-4V fabricated via laser powder bed fusion. Moreover, the model is used to simulate texture evolution during rolling of the material to large strains. It is demonstrated that the model is capable of predicting plastic anisotropy/asymmetry and the concomitant texture evolution. While the model reveals a significant effect of habit plane inclination with respect to the loading direction on yield stress, the comparison of the data and model predictions shows that a random distribution of habit planes fits the flow response. It is further inferred that the tension-compression asymmetry arises from the non-Schmid effects.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Understanding the plasticity contributions during laser-shock loading and spall failure of Cu microstructures at the atomic scales

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 (nc) 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. Here, 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.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Atomistic Mechanisms of Stress-Dependent Molten Salt Corrosion in NiCr Alloys

Ni-based structural alloys in molten salt environments often experience simultaneous mechanical loading and corrosive attack, yet the mechanisms governing stress-corrosion interactions remain unclear. Prior studies largely emphasize tensile stress, while the role of compressive stress has received limited attention. Here, reactive molecular dynamics simulations are used to investigate the coupled effects of applied strain and corrosion in Ni 0.75 Cr 0.25 exposed to molten FLiNaK at 800 °C. A Σ5(210) grain boundary model is subjected to tensile (+4%) to compressive (−4%) uniaxial strains, and corrosion behavior is evaluated through fluorine adsorption, charge redistribution, and grain boundary evolution. Tensile strain accelerates intergranular corrosion susceptibility by reducing local atomic packing through elastic dilation and increasing excess free volume at the grain boundary, which enhances atomic mobility and salt infiltration. In contrast, compressive strain can suppress corrosion by promoting the formation of a ridge-like surface layer along the grain boundary, limiting salt access to the underlying alloy. These results provide atomistic insight into how stress states influence grain boundary corrosion in molten salts.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Unveiling nano-scale crystal deformation using coherent X-ray dynamical diffraction

Visualization of internal deformation fields in crystalline materials helps bridge the gap between theoretical models and practical applications. Applying Bragg coherent diffraction imaging under X-ray dynamical diffraction conditions provides a promising approach to the longstanding challenge of investigating the deformation fields in micron-sized crystals. Here, we present an automatic differentiation-based reconstruction method that integrates dynamical scattering theory to accurately reconstruct deformation fields in large crystals. Using this forward model, our simulated and experimental results demonstrate that three-dimensional local strain information inside a large crystal can be accurately reconstructed under coherent X-ray dynamical diffraction conditions with Bragg coherent X-ray diffraction imaging. These findings open an avenue for extending the investigation of local deformation fields to microscale crystals while maintaining nanoscale resolution, leveraging the enhanced coherence and brightness of advanced X-ray sources.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Frontal Maintenance in Submesoscale Flows

Classic deformation theory includes parameters—divergence, total strain, and vorticity—that are invariant to changes in the coordinate system. However, these parameters are sometimes ambiguous with respect to characterizing how fronts are formed and maintained because the presence of a front imposes a reference coordinate system. To help remedy this ambiguity, we propose a framework in frontal coordinates based on along- and cross-front velocity gradients to better characterize frontal maintenance, which can also be used to define divergence and normal strain in frontal coordinates. The framework with these four parameters (along-, cross-front velocity gradients, divergence, and normal strain in frontal coordinate) defines eight characteristic flow types at a front, providing a complete characterization of the flow that strengthens or weakens a front. This framework highlights the importance of the “strain efficiency” concept, which unambiguously defines the contribution of total strain to frontogenesis. Two examples, one based on a realistic simulation of submesoscales in the northern Gulf of Mexico and the other based on an idealized model with similar flow characteristics, are provided to demonstrate how this framework can be used to enhance our understanding of frontal dynamics in submesoscale flows.

frontogenesis↗