Search NASA⌕ Search

SEARCH · Search NASA

Results for “COUPLED STRESS”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A parallel discrete dislocation dynamics/kinetic Monte Carlo method to study non-conservative plastic processes

Non-conservative processes play a fundamental role in plasticity and are behind important macroscopic phenomena such as creep, dynamic strain aging, loop raft formation, etc. In the most general case, vacancy-induced dislocation climb is the operating unit mechanism. While dislocation/vacancy interactions have been modeled in the literature using a variety of methods, the approaches developed rely on continuum descriptions of both the vacancy population and its fluxes. However, there are numerous situations in physics where point defect populations display heterogeneous concentrations and/or non-smooth kinetics. Here, a kinetic Monte Carlo (kMC) approach for modeling vacancy transport in response to arbitrary stress fields is used. Vacancies are treated as point particles and are coupled to the dislocation substructure representing a deformed material via an advection term defined by the local stress gradients. The stress fields and the dislocation substructure are evolved using a discrete dislocation dynamics (DDD) module. To extend the coupled model to the treatment of large systems, we have implemented it in the massively-parallel DDD code ParaDiS. To avoid numerical incompatibilities associated with merging deterministic (DDD) and stochastic (kMC) integration algorithms, we cast the entire elasto-plastic-diffusive problem within a single stochastic framework, taking advantage of a parallel kMC algorithm to evolve the system as a single event-driven process. The large-scale implementation enables the study of the evolution of a variety of dislocation-defect scenarios governed by non-conservative transport kinetics. After carrying out an exhaustive numerical and computational analysis of our parallel algorithm, we show results that emphasize situations where inhomogeneous vacancy dynamics are of relevance, and compare discrete kinetics to continuum solutions for several cases.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Acoustoelastic characterization of plates using zero group velocity Lamb modes

Acoustoelasticity, a characteristic of material anharmonicity, gives rise to a link between wave propagation velocity and the stress state in materials. Ultrasonic techniques to monitor this coupling, particularly with high sensitivity and in a noncontact manner, can have widespread application both in the quantification of applied and residual stress and in the characterization of nonlinear material behavior through measurement of higher order elastic constants. Here, we use a laser ultrasonic technique to excite and detect zero group velocity (ZGV) Lamb wave resonances in aluminum plates under uniaxial loading. A laser line source is used to excite these resonances at different orientations with respect to the applied load, and the signals are detected using an interferometer. The effects of stress and source orientation on ZGV resonance frequencies are validated using the theory of acoustoelastic Lamb wave propagation. In addition, a model-based inversion technique is used to extract Murnaghan's third-order elastic constants from measurements of the stress dependence of the first two ZGV modes generated parallel and perpendicular to the applied load. Laser generation and detection of ZGV resonances is shown to be an effective and powerful approach for the noncontact and nondestructive acoustoelastic characterization of elastic waveguides.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A micromorphic filter for determining stress and deformation measures from direct numerical simulations of lower length scale behavior

A micromorphic filter is presented for the extraction of stress and deformation measures from underlying direct numerical simulation (DNS) of microstructural mechanical response. The filter is consistent with the micromorphic continuum theory of Eringen and Suhubi (1964), and allows the interrogation of a DNS to generate higher order constitutive models. The filter adopts aspects of the overlap coupling techniques for atomistic-continuum scale-bridging in order to update the nodal degrees of freedom of the finite-element-based filter. Here, two methods for the extraction of quantities of interest from the underlying DNS are presented. The first is a least squares fitting technique which is appropriate for finite-element-style macro-scale domains, and the second is a variationally based approach which generalizes the method to any macro-scale domain which has a variational form. No restrictions on the form of the micro-scale simulation are inherent to the approach, provided that the required quantities (stress, position, etc.) can be provided. The filter is applied to a simple homogeneous medium undergoing homogeneous deformation to demonstrate the recovery of a classical continuum response, what we call the null case. We then apply the framework to the same homogeneous medium with a central and edge penny-shaped crack to show micromorphic effects. We demonstrate that varying the size of the filter causes the micromorphic behaviors to arise and then dissipate as the filter becomes smaller than the crack's length scale. Additionally we show the effect of filter size on a chiral element network along with a plate perforated with evenly spaced and perturbed holes through the thickness. The micromorphic effects are shown to be present for the chiral elements under any loading but become most pronounced when a single element is located within each averaging domain. Micromorphic effects become pronounced for the plate upon the onset of localization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Final Report for Closing the Loop Between In Situ Stress Complexity and EGS Fracture Complexity

Characterizing in situ stress is essential in Enhanced Geothermal System (EGS) for both risk assessment and operation design. Standard stress measurement techniques such as diagnostic fracture injection tests (DFITs) require creating hydraulic fractures and apply analytical methods such as the G -function analysis to estimate the stress. These classical approaches rely on several assumptions, including but not limited to: (i) the induced fracture is planar and persistent, and (ii) the fracture plane is normal to the minimum horizontal stress. However, these assumptions may not hold in EGS environments like Utah FORGE, where thermo-hydro-mechanical (THM) coupling, wellbore deviation, and rock heterogeneity can alter the local stress condition and consequently lead to fracture complexities in both near-well and far-field regions.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Estimating perturbed stress from 3-D borehole displacements induced by fluid injection in fractured or faulted shales

SUMMARY Hydrofracturing stress measurements in fractured and anisotropic shales are notoriously difficult, because opening of existing geological features tends to prevent the creation of a pure hydraulic fracture perpendicular to the least compressive principal stress. Here we show how adding 3-D borehole-displacement measurements while conducting the hydraulic injection test helps to better constrain the principal stress orientations and magnitudes. We developed a 3-D fully coupled hydromechanical numerical model to analyse the displacement, fluid pressure and injection flow-rate data measured during an injection pressure-step-rate test conducted to activate a faulted borehole interval in the Mont Terri Opalinus Clay (Switzerland). We find that injected fluids can only penetrate the fault when it is at or above the Coulomb failure pressure. Borehole displacement orientations are sensitive to a ∼15° variation in the stress–tensor orientation and a 1 MPa stress magnitude variation. Although some dispersion occurs while rupture is propagating along the fault plane ∼4 m away from the borehole, the maximum density of displacement orientations consistently informs about the stress orientation. Thus, an extended injection step-rate approach coupled with an accurate in situ measurement of the borehole wall displacements can be used to better constrain the local stress field perturbations in fractured shales and in heterogeneous rock in general.

Guglielmi, Yves↗
DOE OSTI2021

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗
DOE OSTI2020

Residual Stress Modeling with Phase Transformation for Wire Arc Additive Manufacturing of B91 Steel

Wire arc additive manufacturing (WAAM) is an energy-efficient manufacturing technique used for near-net-shape production of functional industrial components. However, heat accumulation during deposition and the associated mechanical and metallurgical changes result in complex residual stress profiles across the cross section of the fabricated components. These residual stresses are detrimental to the service life of the components. In this study, sequentially coupled thermomechanical analysis of WAAM B91 steel is conducted to quantify the residual stress variation across the component. The thermomechanical analysis includes a transient heat transfer model and a static stress model that incorporates the transformation-induced plasticity due to martensitic phase transformation. Here, the experimentally calibrated heat transfer model mirrors the temperature variation of the system during the deposition. The results from the stress model are validated via x-ray diffraction measurements, and the numerical results are in good agreement with the experimental data.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2026

Evaluating ecosystem water use efficiency under drought stress: a case study of the Helan Mountain region, northwest China

Context Water use efficiency (WUE) is a fundamental ecological indicator links carbon assimilation and water loss in terrestrial ecosystems. Understanding its responses to drought stress is essential for adaptive ecosystem management, particularly in climate-sensitive mountain landscapes. Objectives This study aimed to investigate drought-driven variations in WUE across major vegetation types in the Helan Mountain region of Northwest China. Specifically, we sought to identify dominant ecological drivers of WUE variability and to disentangle their relative importance and causal pathways. Methods We quantified WUE using the Moderate Resolution Imaging Spectroradiometer (MODIS) products and the Drought Severity Index (DSI) data from 2001 to 2020. To examine WUE – drought relationships across contrasting vegetation types, we employed a spatially explicit analytical framework integrating Random Forest (RF) modeling, partial correlation analysis, and structural equation modeling (SEM). Results Regional WUE exhibited relatively stable interannual dynamics, yet pronounced spatial heterogeneity that was strongly modulated by drought conditions. Vegetation properties, particularly Leaf Area Index (LAI) and Normalized Difference Vegetation Index (NDVI), emerged as the dominant determinants of WUE, with NDVI alone explaining over 20% of its spatial variance in forest and grassland during non-drought periods. SEM analyses revealed that climate forcing influenced WUE mainly through indirect pathways mediated by soil moisture availability and vegetation structural dynamics, rather than through direct climatic controls. Among all regulating factors, LAI acted as the central control node governing ecosystem carbon–water coupling. In contrast, short-term climatic stress, especially atmospheric demand and drought duration, exerted weak or negative direct effects on WUE. Ecosystem-specific responses were observed, with croplands mainly regulated by soil water availability, whereas forests and grasslands showed more sensitive to atmospheric drought stress. Together, these results reveal a hierarchical control framework where soil–vegetation interactions mediate climate impacts on WUE, driving strong spatial heterogeneity in drought responses across mountain landscapes. Conclusions Our findings highlight the pivotal role of indirect drought effects mediated by vegetation and soil processes in shaping ecosystem WUE. The identified soil–vegetation–climate regulatory hierarchy provides mechanistic insight into landscape–scale drought sensitivity and supports integrated modeling approaches for evaluating ecosystem resilience and sustainable management in arid mountain regions.

China↗
DOE OSTI2021

Predicting Cement-Sheath Integrity with Consideration of Initial State of Stress and Thermoporoelastic Effects

In conventional wellbore-integrity analysis, the cement sheath’s initial state of stress and transient thermoporoelastic effects are often neglected. However, the initial state of stress is prerequisite information for accurately predicting the safe operating conditions that prevent a cemented well from being damaged. In addition, transient thermoporoelastic effects can have a profound effect on when damage will occur. Here in this paper, we propose a model that includes these effects to predict the safe operating pressures and temperatures that will prevent cement-sheath failure. For the initial state of stress, we proposed an empirical model using measurements. Subsequent stress changes are evaluated by a fully coupled transient thermoporoelastic model to analyze the mechanical behavior of the cement sheath. We predict the safe operating envelope (SOE) for shear, tensile, and debonding cement-sheath failures caused by pressure and temperature perturbations after the cement sets. Our model predicts that pore pressure is a key factor for cement failure, especially for rapid temperature changes. If the formation is low permeability, the transient pore pressures are amplified, causing the risk of damage to increase. Compared with conventional thermoelastic models, the thermoporoelastic model predicts a smaller SOE when heating the internal casing fluid and a larger envelope when cooling the internal casing fluid. Finally, the heating rate was considered with respect to field applications. The heating rate was also considered, and slower heating/cooling rates can prevent damage to the cement sheath. Finally, the thermoporoelastic model was applied to explain several laboratory and field experiments and achieved good matches.

02 PETROLEUM↗
DOE OSTI2024

Sequential Dual Alignments Introduce Synergistic Effect on Hexagonal Boron Nitride Platelets for Superior Thermal Performance

Planarly aligning 2D platelets is challenging due to their additional orientational freedom compared to 1D materials. Here, this study reports a sequential dual-alignment approach, employing an extrusion-printing-induced shear force and rotating-magnetic-field-induced force couple for platelet planarly alignment in a yield-stress support bath. It is hypothesized that the partial alignment induced by a directional shear force facilitates subsequent axial rotation of the platelets for planar alignment under an external force couple, resulting in a synergistic alignment effect. This sequential dual-alignment approach achieves better planar alignment of 2D modified hexagonal boron nitride (mhBN). Specifically, the thermal conductivity of the 40 wt% mhBN/epoxy composite is significantly higher (692%) than that of unaligned composites, surpassing the cumulative effect of individual methods (only 133%) with a 5 times more synergistic effect. For 30, 40, and 50 wt% mhBN composites, the thermal conductivity values (5.9, 9.5, and 13.8 W m –1 K –1 ) show considerable improvement compared to the previously reported highest values (5.3, 6.6, and 8.6 W m –1 K –1 ). Additionally, a 3D mhBN/epoxy heat sink is printed and evaluated to demonstrate the feasibility of device fabrication. The approach enables the planar alignment of electrically or thermally conducting 2D fillers during 3D fabrication.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Data‐Driven Safety Risk Prediction of Lithium‐Ion Battery

Abstract Inevitable safety issues have pushed battery engineers to become more conservative in battery system design; however, battery‐involved accidents still frequently are reported in headlines. Identifying, understanding, and predicting safety risks have become priorities to further accelerate technology and industry development. However, diverse loading scenarios, significantly varied stress‐induced short circuit mechanisms, and highly coupled mechanical–electrochemical safety behaviors have remained grand challenges. Herein, the safety risk is termed as the probability of the mechanical triggering of an internal short circuit, to reflect the safety related behaviors of lithium‐ion batteries. Based on a mechanical model and experimental results, a sufficient dataset is generated consisting of strain states and their corresponding safety risks, covering both cylindrical and pouch cells, various states of charges, and loading conditions. Machine‐learning tools combined with the established finite element mechanical model are applied to predict the safety risks of the cells. The results achieve a high level of accuracy on the test data (the relative error of the average short circuit prediction deviation is less than 6.2%.). This work underpins the safety risk concept and highlights the promise of physics combined with data‐driven modeling methodology to predict the safety behaviors of energy storage systems.

Jia, Yikai↗
DOE OSTI2024

A mechanistic interpretation of Nelson curves for PVP failures under high temperature hydrogen attack

As an empirically established design criterion, Nelson curves that relate the service temperature and the allowable hydrogen partial pressure have been developed and utilized for more than sixty years in pressure vessels and piping (PVP) safety design. Despite a relatively clear thermodynamic understanding of the high-temperature-hydrogen-attack (HTHA) problem, the detailed fracture process on the microstructural length scales, however, remains elusive, and a quantitative assessment of the PVP lifetime under HTHA from the available creep fracture dataset is still not possible. This work develops a microstructure-informed and micromechanics-based model by incorporating a synergy between hydrogen transport and intergranular-cavity-based fracture process. Based on the available creep lifetime data of C-0.5Mo steels, we are able to calibrate material constitutive parameters, and then conduct nonlinear finite element simulations that reveal a real-time stress-induced hydrogen diffusional transport along grain boundaries, coupled with a microstructure-explicit failure process, from which Nelson curves can be computed. Such failure analyses allow us to delineate two distinct regimes on the Nelson curves, i.e., dislocation-creep-controlled or grain boundary diffusion-assisted cavity growth. More importantly, we found that a small change of the pipe thickness and applied stresses can significantly shift these lifetime curves. However, these two parameters are usually not provided in Nelson curves, thus limiting their usage in material selection and safety design. In conclusion, this discrepancy can clearly be mitigated by extensive parametric studies from our micromechanical modeling/simulation framework.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Corrosive Influence of Carbon Dioxide on Crack Initiation in Quartz: Comparison With Liquid Water and Vacuum Environments

The stimulation of crack growth in quartz and siliceous materials by injecting carbon dioxide (CO 2 ) represents a key technology in long-term carbon storage and in the development of natural gas wells. While this technology is widely used, the molecular impact of CO 2 interactions on the solid matrix is only incompletely understood. In this work, we employ reactive molecular dynamics simulations to study how the CO 2 fluid environment affects the mechanical properties of pre-cracked single-crystal quartz. The thermodynamic conditions of interest are those relevant to subsurface reservoirs. We report how structural properties of quartz—bond length distribution and crack tip shape—evolve upon introduction of a fluid. These properties are directly related to macroscopic quantities of the global stress-strain curves, thus reaffirming the inherent coupling across multiple scales for fluid-solid interactions in the subsurface. Here, we find that CO 2 reduces the fracture toughness of quartz by 12.1% compared to that of quartz in vacuum, thereby promoting crack growth and enhancing fluid transport in the subsurface.

58 GEOSCIENCES↗
DOE OSTI2024

Revisiting shear stress tensor evolution: Nonresistive magnetohydrodynamics with momentum-dependent relaxation time

This study aims to develop second-order relativistic viscous magnetohydrodynamics (MHD) derived from kinetic theory within an extended relaxation time approximation (momentum/energy dependent) for the collision kernel. The investigation involves a detailed examination of shear stress tensor evolution equations and associated transport coefficients. The Boltzmann equation is solved using a Chapman-Enskog-like gradient expansion for a charge-conserved conformal system, incorporating a momentum-dependent relaxation time. The derived relativistic nonresistive, viscous second-order MHD equations for the shear stress tensor reveal significant modifications in the coupling with dissipative charge current and magnetic field due to the momentum dependence of the relaxation time. By utilizing a power law parametrization to quantify the momentum dependence of the relaxation time, the anisotropic magnetic field-dependent shear coefficients in the Navier-Stokes limit have been investigated. The resulting viscous coefficients are seen to be sensitive to the momentum dependence of the relaxation time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

TEC (ThermoElastoCreep) [SWR-25-167]

This repository provides a Python script to simulate the thermomechanical behavior of Concentrating Solar Power (CSP) receivers under periodic thermal loading. The simulation couples elastic and creep strains to predict stress relaxation, fatigue, and rupture life of the receiver structure. The work was done under the project Light Trapping Planar Cavity Receiver (LTPCR) funded by U.S. Department of Energy. All other relevant work and publications from LTPCR project can be found at https://github.com/NatLabRockies/LTPCR-Public

He, Xin [National Laboratory of the Rockies (NLR),↗
DOE OSTI2020

Multi‐scale simulation of rock compaction through breakage models with microstructure evolution

Abstract. Regional subsidence due to fluid depletion includes the interaction among multiple physical processes. Specifically, rock compaction is governed by coupled hydro-mechanical feedbacks involving fluid flow, effective stress change and pore collapse. Although poroelastic models are often used to explain the delay between depletion and subsidence, recent evidence indicates that inelastic effects could alter the rock microstructure, thus exacerbating coupling effects. Here, a constitutive law built within the framework of Breakage Mechanics is proposed to account for the inherent connection between rock microstructure, hydraulic conductivity, and pore compaction. Furthermore, it is embedded into a 1-D hydromechanical coupled finite element analysis (FEA) to explore the effects of micro-structure rearrangement on the development of reservoir compaction. Numerical examples with the proposed model are compared with simulations under constant hydraulic conductivity to illustrate the model capability to capture the non-linear processes of reservoir compaction induced by fluid depletion.

Buscarnera, Giuseppe↗
DOE OSTI2022

Modeling Isoprene Emission Response to Drought and Heatwaves Within MEGAN Using Evapotranspiration Data and by Coupling With the Community Land Model

We introduce two new drought stress algorithms designed to simulate isoprene emission with the Model of Emissions of Gases and Aerosols from Nature (MEGAN) model. The two approaches include the representation of the impact of drought on isoprene emission with a simple empirical approach for offline MEGAN applications and a more process-based approach for online MEGAN in Community Land Model (CLM) simulations. The two versions differ in their implementation of leaf-temperature impacts of mild drought. For the online version of MEGAN that is coupled to CLM, the impact of drought on leaf temperature is simulated directly and the calculated leaf temperature is considered for the estimation of isoprene emission. For the offline version, we apply an empirical algorithm derived from whole-canopy flux measurements for simulating the impact of drought ranging from mild to severe stage. In addition, the offline approach adopts the ratio ($f$ PET ) of actual evapotranspiration to potential evapotranspiration to quantify the severity of drought instead of using soil moisture. We applied the two algorithms in the CLM-CAM-chem (the Community Atmosphere Model with Chemistry) model to simulate the impact of drought on isoprene emission and found that drought can decrease isoprene emission globally by 11% in 2012. We further compared the formaldehyde (HCHO) vertical column density simulated by CAM-chem to satellite HCHO observations. We found that the proposed drought algorithm can improve the match with the HCHO observations during droughts, but the performance of the drought algorithm is limited by the capacity of the model to capture the severity of drought.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Calibration technique of wall shear stress sensors using oscillating plate

Disclosed are a system and methods for calibrating wall shear stress sensors. The system includes an oscillating plate coupled to an actuator and mounted on a rolling elements, and one or more sensors coupled to a height adjusting device. The system can further comprise a height control rod coupled to a height control base and a sensor holder configured to house the one or more sensors and supported on a connector, the connector configured to be rotatably disposed about the height control rod. The system can be calibrated by causing the actuator to oscillate the oscillating plate at a frequency, sensing, using the one or more sensors, shear stress at a wall, the shear stress at the wall being associated with a velocity field, and determining a theoretical wall shear stress based on fluid properties, the frequency, and the height of the one or more sensors above the oscillating plate.

Freidkes, Brett R.↗