Search NASA⌕ Search

SEARCH · Search NASA

Results for “stress equilibrium”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Imposing equilibrium on experimental 3-D stress fields using Hodge decomposition and FFT-based optimization

Here, we present a methodology to impose micromechanical constraints, i.e. stress equilibrium at grain and sub-grain scale, to an arbitrary (non-equilibrated) voxelized stress field obtained, for example, by means of synchrotron X-ray diffraction techniques. The method consists in finding the equilibrated stress field closest (in L 2 -norm sense) to the measured non-equilibrated stress field, via the solution of an optimization problem. The extraction of the divergence-free (equilibrated) part of a general (non-equilibrated) field is performed using the Hodge decomposition of a symmetric matrix field, which is the generalization of the Helmholtz decomposition of a vector field into the sum of an irrotational field and a solenoidal field. The combination of: a) the Euler–Lagrange equations that solve the optimization problem, and b) the Hodge decomposition, gives a differential expression that contains the bi-harmonic operator and two times the curl operator acting on the experimental stress field. These high-order derivatives can be efficiently performed in Fourier space. The method is applied to filter the non-equilibrated parts of a synthetic piecewise constant stress fields with a known ground truth, and stress fields in Gum Metal, a beta-Ti-based alloy measured in-situ using Diffraction Contrast Tomography (DCT). In both cases, the largest corrections were obtained near grain boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Self-equilibrated backstresses induce compensation between hardening and softening: Micromechanical and microstructural features

For crystal plasticity formulation neglecting the strain gradients at the dislocation scale, backstresses are necessary to satisfy the mesoscopic stress equilibrium and strain compatibility requirements. Backstresses in a material point and those between neighboring points are treated separately. Physically backstresses can dynamically induce either hardening or softening effects to adjust the dislocation hardening of slip systems, whereas the sum of the stresses resulting in backstresses in the representative volume element (RVE) is zero. Backstress hardening suppresses certain slip activities, while backstress softening triggers excess dislocations. A two-phase brass is studied. Backstresses make the hard phase deform smaller and the soft phase deform more than those without employing backstresses in a model. In the former, more homogeneous lattice strain and higher yield stress are obtained. Backstress has limited effects on the evolution of the critical resolved shear stress and crystallographic texture. Furthermore, it is a major cause of the banded structures featuring sharp crystallographic contrasts which are not parallel to the slip planes. Tests indicate that backstress induces a serrated behavior at the early stages of deformation and induces the geometrically necessary bands (GNBs), which are parallel to the macroscopic axes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electrochemical-Mechanical Coupling Strongly Affects the Performance of Nanopore, Thin-Film, and Solid-State Batteries

Here, the present study focuses on the electrochemical-mechanical (ECM) coupling effects of a thin-film, solid-state battery with only stiff, ceramic materials, in contrast to prior investigations that focus on individual active material particles or aggregations of particles. We model the impacts of ECM couplings including stress-transport and stress-equilibrium potential on the full-cell performance and potential mechanical failure modes of a thin-film battery conformally deposited in a nanopore scaffold, which is an experimentally achievable device. Model results indicate electrode volume changes due to lithium insertion or removal, along with mechanical boundary conditions, result in stress gradients that alter the lithium-ion flux, reduce lithium concentration gradients, and improve cell rate capability. However, the high stress levels in the cell can also lead to mechanics-related failure such as the separation of cell layers. For the parameter set in this work, stress-transport coupling has a much greater influence on rate capability than stress-potential coupling. Optimization of thin-film batteries to harness the benefits of ECM coupling effects requires leveraging geometric design and material selection. The current work underscores the need for further theoretical and experimental investigation into ECM coupling effects in thin-film batteries to enhance their understanding and design optimization.

25 ENERGY STORAGE↗
DOE OSTI2023

A study of explosive-induced fracture in polymethyl methacrylate (PMMA)

The fracture response of geologic materials is of interest for applications, including geothermal energy harnessing and containment of underground explosions. To better understand the explosively induced fracture response of geomaterials, polymethyl methacrylate (PMMA) was used as a transparent rock surrogate to allow imaging of internal shock propagation and fracture growth processes. Experiments were conducted using high-speed shadowgraphy and photon Doppler velocimetry (PDV), which were compared to numerical simulations. Experiments measured fractures produced in 304.8mm × 304.8mm × 304.8mm PMMA cubes with two simultaneously initiated detonators. The cubes were subjected to varying amounts and directions of externally applied uniaxial stresses, including no stress, 2 MPa stress, and 20 MPa stress. The fracture radius as a function of time was extracted from the high-speed videos. Post-test images of the PMMA cubes aided in the determination of three-dimensional effects not directly imaged by the cameras. The surface velocity history and the shock response captured in PDV and the high-speed videos were compared to the simulated explosive-induced shock response. The simulation results indicate that the shock drives the fracture for the first 20 μs corresponding to a fracture radius of approximately 15 mm in the experiments. The gas-driven fracture extent was estimated analytically using an equilibrium stress distribution calculated after the shock wave propagation through the sample. Reduction in the gas pressure due to the leakage of the explosive products through the crack as a function of time was accounted for. In conclusion, the estimated fracture lengths were in agreement with the experimentally observed fracture lengths.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

A FFT-based mesoscale continuum dislocation mechanics with defect energy: Applications to composites and polycrystals

A crystal plasticity elastoviscoplastic FFT (fast Fourier transform) formulation with a mesoscale continuum field dislocation mechanics model is presented, which incorporates a defect energy density that depends on GND densities and an associated material length scale. This allows to thermodynamically derive internal length scale dependent intra-crystalline backstress and Peach–Koehler force acting on GND densities. The model considers GND density evolution through a filtered numerical spectral approach, which is coupled with stress equilibrium through the elastoviscoplastic FFT algorithm. The discrete Fourier transform (DFT) method together with finite difference (FD) schemes is applied to solve both the backstress tensor and the Fourier–Green operator. Numerical results are first reported for two-phase laminate composites with plastic single crystal channels and elastic precipitates for shear loadings. Channel size effects are simulated and analyzed on the overall and local hardening behaviors during monotonous loadings. In addition, the evolutions of GND densities and the role of their associated backstress on size effects are examined during reversible shear loading. In a second part, the role of the defect energy internal length scale on polycrystal’s hardening during tension–compression is discussed. The results are compared to those obtained using FFT-based continuum field dislocation mechanics without defect energy.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Evaluating Novel Polymer-Composite Resins for the Ability to Derive Elastic Properties for Explicit Formulations, or to Develop Nonlinear Models Thereof

This project evaluates the ability to derive elastic material properties for novel additively manufactured carbon-polymer matrix materials. Initially stress-strain equilibrium equations were thought to be sufficient to characterize the elastic region of material performance. However, initial test data shows the cured material behaves nonlinearly in the elastic region. To properly determine the elastic performance, additional strain data must be collected in future testing. This paper discusses the initial test data, methods and models used to describe nonlinear elasticity, and the determination of mean elastic material properties, such as Lame’s constants, using the stress-strain equilibrium equations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phase field modeling of dislocations and obstacles in InSb

We present a phase-field dislocation dynamics (PFDD) model informed by first-principle calculations to elucidate the competitive dislocation nucleation and propagation between the glide and shuffle sets in InSb diamond cubic crystal. The calculations are directly informed with generalized stacking fault energy curves on the (111) slip plane for both the “glide set,” with the smaller interplanar spacing, and the “shuffle set,” with the larger interplanar spacing. The formulation also includes elastic anisotropy and the gradient term associated with the dislocation core. The PFDD calculations show that under no stress the equilibrium structure of screw glide set dislocations dissociates into Shockley partials, while those of the shuffle set dislocations do not dissociate, remaining compact. The calculated dislocation core widths of these InSb dislocations agree well with the measured values for other semiconductor materials, such as Si and GaN. We find that a shuffle set dislocation emits from a dislocation source at an applied stress about three times smaller than that needed to emit leading and trailing partials successively on the glide set plane. Once the partial dislocations in the glide set are emitted, they propagate faster than the shuffle set perfect dislocation at the same stress level.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Tuning nonequilibrium phase transitions with inertia

In striking contrast to equilibrium systems, inertia can profoundly alter the structure of active systems. Here, we demonstrate that driven systems can exhibit effective equilibrium-like states with increasing particle inertia, despite rigorously violating the fluctuation–dissipation theorem. Increasing inertia progressively eliminates motility-induced phase separation and restores equilibrium crystallization for active Brownian spheres. This effect appears to be general for a wide class of active systems, including those driven by deterministic time-dependent external fields, whose nonequilibrium patterns ultimately disappear with increasing inertia. The path to this effective equilibrium limit can be complex, with finite inertia sometimes acting to accentuate nonequilibrium transitions. The restoration of near equilibrium statistics can be understood through the conversion of active momentum sources to passive-like stresses. Unlike truly equilibrium systems, the effective temperature is now density dependent, the only remnant of the nonequilibrium dynamics. Further, this density-dependent temperature can in principle introduce departures from equilibrium expectations, particularly in response to strong gradients. Our results provide additional insight into the effective temperature ansatz while revealing a mechanism to tune nonequilibrium phase transitions.

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

Bottlebrush Polymers at Liquid Interfaces: Assembly Dynamics, Mechanical Properties, and All-Liquid Printed Constructs

Bottlebrush polymer surfactants (BPSs), formed by the interfacial interactions between bottlebrush polymers (BPs) with poly(acrylic acid) side chains dissolved in an aqueous phase and amine-functionalized ligands dissolved in the oil phase, assemble and bind strongly to the fluid-fluid interface. The ratio between N BB (backbone degree of polymerization) and N SC (side chain degree of polymerization) defines the initial assembly kinetics, interface packing efficiency, and stress relaxation. The equilibrium interfacial tension (γ) increases when N BB < N SC , but decreases when N BB $\gg$ N SC , correlating to a pronounced change in the effective shape of the BPs from being spherical to worm-like structures. The apparent surface coverage (ASC), i.e., the interfacial packing efficiency, decreases as N BB increases. The dripping-to-jetting transition of an injected polymer solution, as well as fluorescence recovery after photobleaching experiments, revealed faster initial assembly kinetics for BPs with higher N BB . Furhter, Euler buckling of BPS assemblies with different N BB values was used to characterize the stress relaxation behavior and bending modulus. The stress relaxation behavior was directly related to the ASC, reflecting the strong influence of macromolecular shape on packing efficiency. The bending modulus of BPSs decreases for N BB < N SC , but increased when N BB $\gg$ N SC , showing the effect of molecular architecture and multisite anchoring. All-liquid printed constructs with lower N BB BPs yielded more stable structured liquids, underscoring the importance of macromolecular packing efficiency at fluid interfaces. Overall, this work elucidates fundamental relationships between nanoscopic structures and macroscopic properties associated with various bottlebrush polymer architectures, which translate to the stabilization of all-fluidic printed constructs.

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

Experimental Investigation of Low-Frequency Distributed Acoustic Sensor Responses to Two Parallel Propagating Fractures

Low-frequency distributed acoustic sensing (LF-DAS) is a diagnostic tool for hydraulic fracture propagation with far-field monitoring using fiber optic sensors. LF-DAS senses strain rate variation caused by stress field change due to fracture propagation. Fiber optic sensors are installed in the monitoring wells in the vicinity of a fractured well. From the strain responses, fracture propagation can be evaluated. To understand subsurface conditions with multiple propagating fractures, a laboratory-scale hydraulic fracture experiment was performed simulating the LF-DAS response to fracture propagation with embedded distributed optical fiber strain sensors under these conditions. The experiment was performed using a transparent cube of epoxy with two parallel radial initial flaws centered in the cube. Fluid was injected into the sample to generate fractures along the initial flaws. The experiment used distributed high-definition fiber optic strain sensors with tight spatial resolutions. The sensors were embedded at two different locations on opposite sides of the initial flaws, serving as observation/monitoring locations. We also employed finite element modeling to numerically solve the linear elastic equations of equilibrium continuity and stress–strain relationships. The measured strains from the experiment were compared to simulation results from the finite element model. The experimentally derived strain and strain-rate waterfall plots from this study show the responses to both fractures propagating, while the fracture at the lower position took most of the fluid during the experiment. Interestingly, a fracture first began propagating from the upper flaw of the two flaws, but once the lower fracture was initiated, it grew much faster than the upper fracture. Both fibers were intercepted by the lower fracture, further verifying the strain signature as a fracture is approaching and intersecting an offset fiber.

Chemistry↗
DOE OSTI2024

Stress-hybrid virtual element method on six-noded triangular meshes for compressible and nearly-incompressible linear elasticity

In this paper, we present a first-order Stress-Hybrid Virtual Element Method (SH-VEM) on six-noded triangular meshes for linear plane elasticity. Here, we adopt the Hellinger–Reissner variational principle to construct a weak equilibrium condition and a stress based projection operator. In each element, the stress projection operator is expressed in terms of the nodal displacements, which leads to a displacement based formulation. This stress-hybrid approach assumes a globally continuous displacement field while the stress field is discontinuous across each element. The stress field is initially represented by divergence-free tensor polynomials based on Airy stress functions, but we also present a formulation that uses a penalty term to enforce the element equilibrium conditions, referred to as the Penalty Stress-Hybrid Virtual Element Method (PSH-VEM). Numerical results are presented for PSH-VEM and SH-VEM, and we compare their convergence to the composite triangle FEM and B-bar VEM on benchmark problems in linear elasticity. The SH-VEM converges optimally in the L 2 norm of the displacement, energy seminorm, and the L 2 norm of hydrostatic stress. Furthermore, the results reveal that PSH-VEM converges in most cases at a faster rate than the expected optimal rate, but it requires the selection of a suitably chosen penalty parameter.

42 ENGINEERING↗
DOE OSTI2022

Disclination-dislocation based model for grain boundary stress field evolution due to slip transmission history and influence on subsequent dislocation transmission

This work demonstrates how the structure of a grain boundary (GB) and its evolution due to slip transmission history influences subsequent dislocation transmission. First, a model for the evolution of stress fields within grain boundaries that accounts for the effects of coherent dislocation transmission is introduced. Starting with a disclination-based construct of GBs at minimum energy (equilibrium), the model describes the evolution of the GB stress field to a state characteristic of excess energy (non-equilibrium) due to the incorporation of residual Burgers vector content following sequential slip transmission events. Several essential features of this model are verified via molecular dynamics simulations of lattice dislocation absorption. Second, this model is implemented into a discrete dislocation dynamics (DDD) code and simulations are performed to understand the influence of Non-equilibrium GB stress fields, conditioned by the slip transmission history, on subsequent dislocation transmission. DDD simulations reveal that the critical resolved stress necessary for slip propagation can be reduced by with continued absorption of residual dislocation content. Moreover, DDD simulations prove that a comprehensive consideration of both the binding and driving stresses, and the evolution of the transmission configuration is necessary to quantify the influence of the mechanical state of the GB on slip propagation. In general, this work provides important insights into the role of GB structure evolution, conditioned by prior deformation history, on intergranular plasticity.

42 ENGINEERING↗
DOE OSTI2024

Iterative Stress Reconstruction Algorithm to Estimate Three-Dimensional Residual Stress Fields in Manufactured Components

Residual stress (RS) significantly impacts the mechanical performance of components. Measurement of RS often provides incomplete data in terms of components of stress and spatial density. Employing such fields in finite element simulations results in significant modification of the field to achieve equilibrium and compatibility among strains. To overcome this, an iterative stress reconstruction algorithm (ISRA) is developed to estimate 3D RS fields that satisfy equilibrium, are stress component-wise complete, and represent the characterized data sampled. An Al 7075-T651 plate and an additively manufactured (AM) A36 steel wall are considered for RS reconstruction using measurement data from the literature. A maximum variation of ~2.5 MPa in the Al plate, and ~10 MPa in the steel wall are observed between the reconstructed and measured stresses. Furthermore, unknown stress components emerge and reach significant magnitudes (upto ~2.3 MPa in the Al plate and ~45 MPa in the AM wall) during ISRA. Indeed, it is found that minor errors in measurement or data processing are eliminated through the physical requirements during ISRA. Employing a reconstructed RS field is hence not just more accurate given its compatibility, but it additionally corrects for minor errors in measurement. Furthermore, it is found that spatially dense measurement data result in convergence with fewer iterations. Finally, although ISRA yields a nonunique solution dependent on boundary conditions, measurement errors, fitting errors, and mesh density, it accommodates for uncertainties and inaccuracies in measurement, as opposed to failing to reach a physically realistic converged solution.

42 ENGINEERING↗
DOE OSTI2026

Measurements of the inert Hugoniot and observation of mechanical ignition in Ni(V) + Al reactive multilayers via x-ray diffraction

This study investigates the inert Hugoniot response, mechanical ignition, and reaction dynamics of Ni(V)+Al multilayers during longitudinal, laser-driven shock compression experiments. Ni(V)+Al multilayers, known for their self-propagating exothermic reactions, were subjected to longitudinal stresses exceeding 50 GPa using the laser shock facility within the Dynamic Compression Sector (DCS) at the Advanced Photon Source (APS). In situ x-ray diffraction (XRD) revealed that Ni(V) and Al were not in equilibrium during compression, with stress discrepancies attributed to twinning, grain structure effects, and/or dislocation density. However, the measured inert Hugoniot closely matched prior experimental and computational studies, confirming the utility of XRD for measuring the equation of state of thin, complex materials. Additionally, reaction was observed at significantly higher stresses than reported previously using laser-launched flyers. This discrepancy suggests a strong influence of externally imposed shear stress on reaction thresholds, which likely arose from deviations in flyer planarity during past experiments. Full reaction of the multilayer occurred within 40 ns after shock-wave passage, evidenced by complete melting of the constituents. Eulerian hydrocode simulations replicated experimental conditions, providing insights into equilibrium dynamics and experimental artifacts. The results highlight how even small shear forces facilitate ignition in Ni(V)+Al multilayers at lower stresses.

Crystallography↗
DOE OSTI2024

Initial Testing of an In Situ Load Retention Aging Vessel

A thermal aging vessel instrumented with load cells was fabricated. The primary function of the vessel is to continuously monitor the in situ load retention of up to three compressed polymer coupons undergoing thermally accelerated aging under nitrogen. A secondary function is to enable gas sampling of the vessel headspace during thermal aging. Heating of the vessel is achieved using a custom heater jacket. To improve upon our conventional aging study methods which require periodic interruption of aging to perform load testing in an Instron machine at room temperature, this technology aims to automate/facilitate data acquisition/analysis, improve data quality, and enable uninterrupted compression of the polymer which represents the service condition. As an example case to assess functionality of the in situ vessel, the load retention of a siloxane elastomer material additively manufactured by direct-ink-writing (DIW) was measured at three different isothermal aging temperatures for ~1 month. Initial compression of the coupons while near the aging temperature was achieved by temporarily opening the heated vessel to access the interior chamber and manually tightening four nuts to drive the heated compression plate down onto the heated coupons. Initial testing demonstrated achievement of the primary load retention monitoring function. Unfortunately, the vessel leaked which prevented gas sampling; an active purge was used to maintain a nitrogen atmosphere. Welded or otherwise sealed joints, which could be implemented in a future design, would likely eliminate leak paths. To apply time-temperature superposition (TTS), a technique used to provide long-term prediction of the load retention from short-term isothermal data, the load retention needed to be calculated relative to the load at an estimated “equilibrium” time, after most of the transient viscoelastic physical relaxation occurred. The peak load immediately after compression could not be used as the load retention basis for two reasons: (1) age-related changes must be isolated from non-age-related physical relaxation before applying TTS and (2) the manual mechanism used to compress the specimens at the aging temperature was neither smooth nor repeatable which affected the peak load value. To better understand the effect of the mode of initial compression on the measured load, and possibly better estimate “equilibrium” physical relaxation times, systematic stress relaxation experiments were performed using an Instron machine with a thermal chamber. At a given temperature, the DIW polymer was compressed to a fixed strain in either a stepped or continuous manner at two different rates, then held at that strain for 24 hrs. The results indicated that, at a given temperature, the different stress relaxation curves appeared to converge to the same curve at some “equilibrium” time when the non-age-related physical relaxation was mostly complete. Though this observation suggests that the discontinuous manual compression employed by the vessel is feasible, a compression mechanism that is rapid, smooth, and repeatable would enhance its use.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Transferable Water Potentials Using Equivariant Neural Networks

Machine learning interatomic potentials (MLIPs) have emerged as a technique that promises quantum theory accuracy for reduced cost. It has been proposed [J. Chem. Phys. 2023, 158, 084111] that MLIPs trained on solely liquid water data cannot accurately transfer to the vapor–liquid equilibrium while recovering the many-body decomposition (MBD) analysis of gas-phase water clusters. This suggests that MLIPs do not directly learn the physically correct interactions of water molecules, limiting transferability. In this work, we show that MLIPs using equivariant architecture and trained on 3200 liquid water structures reproduces liquid-phase water properties (e.g., density within 0.003 g/cm 3 between 230 and 365 K), vapor–liquid equilibrium properties up to 550 K, the MBD analysis of gas-phase water cluster up to six-body interactions, and the relative energy and the vibrational density of states of ice phases. We show that potentials developed using equivariant MLIPs allow transferability for arbitrary phases of water that remain stable in nanosecond long simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗