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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Strain-Induced Growth of Twisted Bilayers during the Coalescence of Monolayer MoS 2 Crystals

Tailoring the grain boundaries (GBs) and twist angles between two-dimensional (2D) crystals are two crucial synthetic challenges to deterministically enable envisioned applications such as moiré excitons, emerging magnetism, or single-photon emission. Here, we reveal how twisted 2D bilayers can be synthesized from the collision and coalescence of two growing monolayer MoS 2 crystals during chemical vapor deposition. The twisted bilayer (TB) moiré angles are found to preserve the misorientation angle (θ) of the colliding crystals. The shapes of the TB regions are rationalized by a kink propagation model that predicts the GB formed by the coalescing crystals. Optical spectroscopy measurements reveal a θ-dependent long-range strain in crystals with stitched grain boundaries and a sharp (θ > 20°) threshold for the appearance of TBs, which relieves this strain. Reactive molecular dynamics simulations explain this strain from the continued growth of the crystals during coalescence due to the insertion of atoms at unsaturated defects along the GB, a process that self-terminates when the defects become saturated. The simulations also reproduce atomic-resolution electron microscopy observations of faceting along the GB, which is shown to arise from the growth-induced long-range strain. These facets align with the axes of the colliding crystals to provide favorable nucleation sites for second-layer growth of a TB with twist angles that preserve the misorientation angle θ. This interplay between strain generation and aligned nucleation provides a synthetic pathway for the growth of TBs with deterministic angles.

2D materials↗
NASA NTRS1996

The Effect of Strain Rate on the Evolution of Plane Wakes Subjected to Irrotational Strains

Direct numerical simulations of time-evolving turbulent plane wakes developing in the presence of irrotational plane strain applied at three different strain rates have been generated. The strain geometry is such that the flow is compressed in the streamwise direction and expanded in the cross-stream direction with the spanwise direction being unstrained. This geometry is the temporally evolving analogue of a spatially evolving wake in an adverse pressure gradient. A pseudospectral numerical method with up to 16 million modes is used to solve the equations in a reference frame moving with the irrotational strain. The initial condition for each simulation is taken from a previous turbulent self-similar plane wake direct numerical simulation at a velocity deficit Reynolds number, Re, of about 2,000. Although the evolutions of many statistics are nearly collapsed when plotted against total strain, there are some differences owing to the different strain rate histories. The impact of strain-rate on the wake spreading rate, the peak velocity deficit, the Reynolds stress profiles, and the flow structure is examined.

Rogers, Michael M.↗
NASA NTRS1997

Reynolds Stress Balance in Plane Wakes Subjected to Irrotational Strains

Direct numerical simulations of time-evolving turbulent plane wakes developing in the presence of various irrotational plane strains have been generated. A pseudospectral numerical method with up to 25 million modes is used to solve the equations in a reference frame moving with the irrotational strain. The initial condition for each simulation is taken from a previous turbulent self-similar plane wake direct numerical simulation at a velocity deficit Reynolds number, R(sub e), of about 2,000. All the terms in the equations governing the evolution of the Reynolds stresses have been calculated. The relative importance of the various terms is examined for the different strain geometries and the behavior of the individual terms is used to better assess whether the strained wakes are evolving self-similarly.

Rogers, Miichael M.↗
DOE OSTI2022

X-ray free-electron laser based dark-field X-ray microscopy: a simulation-based study

Dark-field X-ray microscopy (DFXM) is a nondestructive full-field imaging technique providing three-dimensional mapping of microstructure and local strain fields in deeply embedded crystalline elements. This is achieved by placing an objective lens in the diffracted beam, giving a magnified projection image. So far, the method has been applied with a time resolution of milliseconds to hours. In this work, the feasibility of DFXM at the picosecond time scale using an X-ray free-electron laser source and a pump–probe scheme is considered. Thermomechanical strain-wave simulations are combined with geometrical optics and wavefront propagation optics to simulate DFXM images of phonon dynamics in a diamond single crystal. Using the specifications of the XCS instrument at the Linac Coherent Light Source as an example results in simulated DFXM images clearly showing the propagation of a strain wave.

47 OTHER INSTRUMENTATION↗
NASA NTRS2016

Applications of Displacement Transfer Functions to Deformed Shape Predictions of the G-III Swept-Wing Structure

In support of the Adaptive Compliant Trailing Edge [ACTE] project at the NASA Armstrong Flight Research Center, displacement transfer functions were applied to the swept wing of a Gulfstream G-III airplane (Gulfstream Aerospace Corporation, Savannah, Georgia) to obtain deformed shape predictions. Four strainsensing lines (two on the lower surface, two on the upper surface) were used to calculate the deformed shape of the G III wing under bending and torsion. There being an insufficient number of surface strain sensors, the existing G III wing box finite element model was used to generate simulated surface strains for input to the displacement transfer functions. The resulting predicted deflections have good correlation with the finite-element generated deflections as well as the measured deflections from the ground load calibration test. The convergence study showed that the displacement prediction error at the G III wing tip can be reduced by increasing the number of strain stations (for each strain-sensing line) down to a minimum error of l.6 percent at 17 strain stations; using more than 17 strain stations yielded no benefit because the error slightly increased to 1.9% when 32 strain stations were used.

displacement theory↗
DOE OSTI2020

Effects of particle inclusions on cracking in ultrathin barrier films

Here, the effects of particle inclusions with different mechanical properties on the development of cracks in ultra-thin barrier films have been studied. Two types of particles, alumina and poly(styrene-b-methyl methacrylate) with large differences in their elastic modulus were used, representing both hard and soft inclusions. Particles were seeded onto polyethylene naphthalate substrates and subsequently coated with atomic layer deposited Al 2 O 3 as a barrier film. Ca corrosion tests used to determine the performance of the barrier films revealed that soft polymer particle inclusions severely affect the development of cracking in barrier films as compared to stiff ceramic particles. Tensile tests also confirmed that the barrier films with polymer particle inclusions have the lowest value of critical onset strain. Parametric simulations, conducted by varying the elastic modulus of the particle inclusions over a wide range, resulted in a larger strain energy release rate near particle inclusions which promotes crack formation and negatively impacts the performance of barrier films with softer particle inclusions. Thus, the combination of high residual stresses in the barrier films with the presence of soft particle inclusions is found the most detrimental for the development of high-performance ultrathin barrier films.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Deformation, dislocation evolution and the non-Schmid effect in body-centered-cubic single- and polycrystal tantalum

A physically-informed continuum crystal plasticity model is presented to elucidate deformation mechanisms, dislocation evolution and the non-Schmid effect in body-centered-cubic (bcc) tantalum widely used as a key structural material for mechanical and thermal extremes. We show the unified structural modeling framework informed by mesoscopic dislocation dynamics simulations is capable of capturing salient features of the large inelastic behavior of tantalum at quasi-static (10 -3 s -1 ) to extreme strain rates (5000 s -1 ) and at low (77 K) to high temperatures (873 K) at both single- and polycrystal levels. Here we also present predictive capabilities of the model for microstructural evolution in the material. To this end, we investigate the effects of dislocation interactions on slip activities, instability and the non-Schmid behavior at the single crystal level. Furthermore, ex situ measurements on crystallographic texture evolution and dislocation density growth are carried out for polycrystal tantalum specimens at increasing strains. Numerical simulation results also support that the modeling framework is capable of capturing the main features of the polycrystal behavior over a wide range of strains, strain rates and temperatures. The theoretical, experimental and numerical results at both single- and polycrystal levels provide critical insight into the underlying physical pictures for micro- and macroscopic responses and their relations in this important class of refractory bcc materials undergoing large inelastic deformations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simulation driven adaptive sampling for neutron-diffraction based strain mapping of additively manufactured parts

Neutron diffraction based strain mapping is a useful technique for measuring residual strains in additively manufactured (AM) metal parts. The measurement is traditionally done by scanning the sample in a point-wise raster pattern to extract the strain at each position. Since the overall scan can span several hours, adaptive sampling approaches using Bayesian optimization based on Gaussian process (BO-GP) regression have been introduced—demonstrating that even with a fraction of the typically made measurements the dominant strain patterns in the sample can be reconstructed. However, the parameters of the BO-GP algorithm have to be carefully chosen for best performance, and the movement time between arbitrary points can offset the time savings from a reduced number of measurement locations. In this paper, we propose algorithms to refine the BO-GP based methods by using simulations of strain patterns in AM parts based on the materials and the process used to print them. We demonstrate that the simulated strain patterns can be used to help choose better parameters for the BO-GP based framework—leading to low reconstruction error for the final strain pattern. Furthermore, we show that the strain mapping experiment can be initialized with a sampling pattern learnt from the simulation data and ordered to reduce movement time, dramatically enabling reduction in the overall time required to run the baseline BO-GP method.

Gaussian process regression↗
NASA NTRS1995

The Evolution of Plane Wakes Subjected to Irrotational Strains

Three direct numerical simulations of time-evolving turbulent plane wakes developing in the presence of various irrotational plane strains have been generated. A pseudospectral numerical method with up to 26 million modes is used to solve the equations in a reference frame moving with the irrotational strain. The initial condition for each simulation is taken from a previous turbulent self-similar plane wake direct numerical simulation at a velocity deficit Reynolds number, Re, of about 2000. Three different plane strains (of the same magnitude) are imposed. In the first two simulations the strain is in a plane normal to the streamwise wake direction (the two cases having strain of opposite sign); in the third the wake is compressed in the streamwise direction and stretched in the inhomogeneous cross-stream direction. The two flows that are stretched in the cross-stream direction experience an exponential increase of Re; flow visualization indicates many small-scale vortices with little or no organized large-scale structure. In the flow that is compressed in the cross-stream direction Re decays exponentially and the layer appears to be relaminarizing. The evolution of several turbulence statistics in each of these flows is examined.

Rogers, M. R.↗
DOE OSTI2026

Strain age cracking in a simulated heat affected zone of Inconel 740H laser-powder bed fusion components

Inconel® 740H components produced via laser-powder-bed fusion (L-PBF) additive manufacturing are ideal for supercritical CO2 primary heat exchangers. Arc welding is often needed, and subsequent post-weld heat treatment aging at 790–840 °C is required to improve strength via gamma prime (γ’) precipitation; however, strain age cracking (SAC) can occur in the heat affected zone (HAZ) during this process. This study uses stress relaxation testing at 800 °C on simulated HAZ specimens from vertically and horizontally built L-PBF IN740H to assess SAC susceptibility across heating rates of 40–3480 °C/h and weld-induced strains of 3–10 %. Vertical builds require greater mechanical energy input and exhibit longer times to fracture than horizontal builds, and time to fracture occurs sooner at slower heating rates. Creep voids were observed in γ’-denuded regions along grain boundaries, and cracking propagated along migrated grain boundaries and at interfaces between elongated secondary phases (γ’ or carbides) and the γ-matrix.

14 SOLAR ENERGY↗
NASA NTRSDate not supplied

Enabling Nuclear and Solar Thermal Propulsion with the Computational Laboratory: Faster and Cheaper Materials Characterization

Nuclear and solar thermal propulsion offer enhanced efficiencies for in-space travel which may open up the outer Solar System to both crewed and automated craft. However, these advanced propulsion systems require cutting-edge materials which are often not easily tested in the laboratory. Ab initio thermodynamics and other computational methods offer a way to test materials at high pressures and temperatures more quickly, and with reduced cost. Here, several case studies are presented which highlight the utility of modern simulation techniques. First, an extensive thermodynamic study was performed at the ab initio level to elucidate the role of hot hydrogen flow on the erosion of four refractory carbides to examine their appropriateness for use as channel coatings. A detailed analysis of reaction products informs estimates of erosion which are validated in comparison with heritage NERVA data. Specific material suggestions are made for both nuclear and solar designs. Secondly, we examined the mechanical behavior and chemical stability of various components of an idealized nuclear cermet core design using techniques spanning from the atomic to microscale. Next, the mechanical and chemical behavior of tungsten (a potential matrix material) and the stability of uranium mononitride (a potential ceramic fuel element), are characterized in the presence of hydrogen and at elevated temperatures with ab initio thermodynamics. Finally, the mechanical response and failure of tungsten under tensile strain is simulated at the micron scale with a combined finite element and dislocation dynamics approach. The resultant stress-strain data agrees well with experiments and leads to a workflow which can incorporate ab initio thermodynamic data into micron-scale simulations and thus provide real-world relevant estimates of erosion and stress-strain behavior. The power of this framework and our plans to use it in the future will be discussed.

William C Tucker↗
NASA NTRS2011

Dynamic Shape Reconstruction of Three-Dimensional Frame Structures Using the Inverse Finite Element Method

A robust and efficient computational method for reconstructing the three-dimensional displacement field of truss, beam, and frame structures, using measured surface-strain data, is presented. Known as shape sensing , this inverse problem has important implications for real-time actuation and control of smart structures, and for monitoring of structural integrity. The present formulation, based on the inverse Finite Element Method (iFEM), uses a least-squares variational principle involving strain measures of Timoshenko theory for stretching, torsion, bending, and transverse shear. Two inverse-frame finite elements are derived using interdependent interpolations whose interior degrees-of-freedom are condensed out at the element level. In addition, relationships between the order of kinematic-element interpolations and the number of required strain gauges are established. As an example problem, a thin-walled, circular cross-section cantilevered beam subjected to harmonic excitations in the presence of structural damping is modeled using iFEM; where, to simulate strain-gauge values and to provide reference displacements, a high-fidelity MSC/NASTRAN shell finite element model is used. Examples of low and high-frequency dynamic motion are analyzed and the solution accuracy examined with respect to various levels of discretization and the number of strain gauges.

Gherlone, Marco↗
DOE OSTI2023

Molecular Dynamics Simulation of Hygroscopic Aging Effects in Epoxy Polymer

The automobile industry is incorporating more lightweight content in car designs to boost fuel-economy. New structural adhesives are needed to mitigate the corrosion and thermal expansion issues associated with joining dissimilar lightweight materials, but adhesive developers lack a fundamental understanding of the chemistry that occurs in the adhesive as the joint ages. In this study, we developed structural adhesive molecular models and applied classical molecular dynamics simulations and density functional theory calculations to gain molecular insights into the influence of water molecules on the properties of epoxy-based adhesives (DGEBA + Jeffamine (JD230)). The simulations were complemented by experimental synthesis and characterization. Our work underscores the impact of water molecules on the local structure of the epoxy network as well as resulting mechanical properties. Water molecules were mainly coordinated with hydroxyls, primary amines and secondary amines, but also weakly coordinated with ether linkages, which were found most probable to be labile. Simulated stress–strain data indicates that increasing the water content deteriorates the mechanical properties. The Young’s modulus decreased by ~ 30% when the water content increased to 3 wt%. We conclude, this integration of molecular-level chemical insights with mechanical property simulations of the hydrated epoxy system and experimental validation holds the promise to advance lightweight joint technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mapping the structural–mechanical landscape of amorphous carbon with ReaxFF molecular dynamics

We use ReaxFF molecular dynamics (MD) to investigate the relationship between structural and mechanical properties in bulk and nanostructured amorphous carbon (a-C). The liquid-quench MD method is used to generate isotropic bulk samples with mass densities ranging from 0.96 to 3.29 g/cm3. Structural analysis identifies two types of structures with distinct short- and medium-range order: lower-density sp2-dominated a-C, which is characterized by a bimodal ring-size distribution, and higher-density sp3-dominated tetrahedral amorphous carbon (ta-C), exhibiting a unimodal ring-size distribution. Stress–strain MD simulations and analysis reveal how an atomistic structure impacts elastic properties and post-yield atomic rearrangements. All stretched structures demonstrate elastic isotropy and plasticity driven by a ring-size expansion mechanism reflected in changes in ring statistics. The plastic region is substantially larger in ta-C than in a-C due to the post-yield shift from sp3 to sp2 C dominant bonding. In both a-C and ta-C, ultimate failure occurs when a reactive crack, traversed by long sp chains, forms and propagates predominantly perpendicular to the direction of the applied strain. Oxygen infiltration into the fractured region significantly reduces stress resistance, primarily through the early rupture of long sp chains. MD simulations and analysis are extended to a-C slabs, a-C nanotubes, and partially a-C nanotubes. The latter nanostructure highlights the differences between the elastically isotropic a-C walls, which develop circumferential cracking, and the crystalline walls, which tear along crystallographic directions. These results provide a strong foundation for further computational characterization of a-C materials.

Dernov, A. (ORCID:0009000004220973)↗
NASA NTRS2004

Effect of anharmonicity of interatomic potential on strain distribution in semiconductor nanostructures

Experiments and theory have shown that the energy spectrum of nanostructures is extremely sensitive to the built-in strain. Knowledge of the strain distribution is therefore Experiments and theory have shown that the energy spectrum of nanostructures is extremely sensitive to the built-in strain. Knowledge of the strain distribution is therefore of utmost importance for the design of optical devices with prescribed light emission spectrum.

nanostructures↗
DOE OSTI2022

FY22 Progress Report on BISON Metallic Fuel Model Development and V&V Using EBR-II Legacy Data

In this report, the activities and achievements made by Argonne National Laboratory for the Nuclear Energy Advanced Modeling and Simulation (NEAMS) BISON code metallic fuel validation and verification project in FY2022 are summarized. The FIPD-BISON integration powered metallic fuel low-burnup swelling framework has been enhanced to cover both radial and axial swelling analyses. Axial-dependent as-fabricated fuel radius profiles have been used to improve the accuracy in evaluating radial swelling strain. The simulations of the IFR experiment X447, which was focused on FCCI/CCCI induced cladding degradation and failure, have been converted into a BISON metallic fuel assessment case with detailed documentation, which is also the first FIPD-BISON integration powered assessment case in BISON. A FIPD-BISON integration repository has been established to be an optional submodule of BISON to support FIPD-BISON integration powered assessment cases. Preliminary trials to simulate the in-cell out-of-pile transient experiment have been made after implementing some essential models such as cladding plasticity models and liquid phase penetration models. The initial results are promising and help identify gaps that need to be filled in FY2023 and beyond, to eventually achieve OPTD-BISON integration powered assessment cases.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Designing Ti-6Al-4V microstructure for strain delocalization using neural networks

Abstract The deformation behavior of Ti-6Al-4V titanium alloy is significantly influenced by slip localized within crystallographic slip bands. Experimental observations reveal that intense slip bands in Ti-6Al-4V form at strains well below the macroscopic yield strain and may serially propagate across grain boundaries, resulting in long-range localization that percolates through the microstructure. These connected, localized slip bands serve as potential sites for crack initiation. Although slip localization in Ti-6Al-4V is known to be influenced by various factors, an investigation of optimal microstructures that limit localization remains lacking. In this work, we develop a novel strategy that integrates an explicit slip band crystal plasticity technique, graph networks, and neural network models to identify Ti-6Al-4V microstructures that reduce the propensity for strain localization. Simulations are conducted on a dataset of 3D polycrystals, each represented as a graph to account for grain neighborhood and connectivity. The results are then used to train neural network surrogate models that accurately predict localization-based properties of a polycrystal, given its microstructure. These properties include the ratio of slip accumulated in the band to that in the matrix, fraction of total applied strain accommodated by slip bands, and spatial connectivity of slip bands throughout the microstructure. The initial dataset is enriched by synthetic data generated by the surrogate models, and a grid search optimization is subsequently performed to find optimal microstructures. Describing a 3D polycrystal with only a few features and a combination of graph and neural network models offer robustness compared to the alternative approaches without compromising accuracy. We show that while each material property is optimized through a unique microstructure solution, elongated grain shape emerges as a recurring feature among all optimal microstructures. This finding suggests that designing microstructures with elongated grains could potentially mitigate strain localization without compromising strength.

Ahmadikia, Behnam↗
DOE OSTI2026

Assessment of the hydromechanical higher-order MPM for the simulation of geotechnical problems

The Material Point Method (MPM) has been increasingly used to simulate large strain deformations. Linear interpolation functions are commonly used to perform the spatial integration. It is well-known that the discontinuities in the interpolation function derivatives induce shock-like artifacts known as ‘cell-crossing’ error. These errors compound with volumetric locking errors when used with hydromechanical formulations for porous media, where different velocity fields are used for each phase. The capabilities of higher-order MPM frameworks have not been explored for real-scale geotechnical problems. As such, this paper aims to assess, validate, and further discuss a higher-order B-spline MPM (BS-MPM) framework. First, the BS-MPM framework is verified against the large-strain oedometer consolidation problem. Second, the framework is validated against a real-scale slope failure experiment triggered by pore water pressure recharge. Landslide features that are captured using the higher-order framework are specifically highlighted, and results (e.g., pore water pressure and deformation) are validated with field measurements. A generally convergent numerical solution is observed when using cubic interpolation functions. Third, a footing penetration problem is simulated using the multi-patch BS-MPM. Trends are examined with respect to penetration velocity and variation in hydraulic conductivity. The BS-MPM framework ultimately presents a stabilized numerical solution that captures plausible hydromechanical interaction trends important in geotechnical engineering applications.

36 MATERIALS SCIENCE↗