Search NASA⌕ Search

SEARCH · Search NASA

Results for “local stress difference”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Controlling Deformation in Al/Ti: How Interface Roughness and Orientation Drive Bimetal Mechanics

The microstructural characteristics and morphology of interfaces in metals can be crucial in governing the initiation of plasticity and early deformation mechanisms under extreme conditions. During high-strain-rate deformation, these interfaces significantly affect dislocation nucleation, twinning, and other mechanisms that directly impact material strength and failure. Despite their importance, a substantial knowledge gap remains between the observed macroscopic material behavior and the underlying role of bimetal interfaces in plasticity initiation. Here, to address this gap, large-scale molecular dynamics simulations are performed on Al/Ti bimetal structures to examine the effect of interface characteristics on the onset of plasticity under uniaxial compression. Specifically, this study investigates how interface roughness (flat vs waveform interfaces) modifies the initiation of plastic events in the microstructure. Atomistic simulations indicate that interface roughness (a microscopic behavior) reduces the stress required for dislocation nucleation, thereby reducing the peak stress relative to a flat interface. For the square interface, the step height strongly influences plasticity initiation by setting the separation of locally flat regions. Varying the interface rotation relative to the loading direction (macroscopic behavior) reveals that the peak stress for both flat and waveform interfaces initially decreases and then increases with rotation, accompanied by a shift in dominant mechanisms—from Al twinning (flat) to interface sliding and Ti-dominated phase transformation and twinning at lower angles. The change in the stress-strain slope during initial compression reflects how rotation alters the resolved shear stress and activates different slip systems. Overall, the study provides valuable insights into the role of bimetal interfaces in controlling plasticity initiation and early deformation pathways for Al/Ti under extreme loading conditions.

36 MATERIALS SCIENCE↗
NASA NTRS2006

User-Defined Material Model for Progressive Failure Analysis

An overview of different types of composite material system architectures and a brief review of progressive failure material modeling methods used for structural analysis including failure initiation and material degradation are presented. Different failure initiation criteria and material degradation models are described that define progressive failure formulations. These progressive failure formulations are implemented in a user-defined material model (or UMAT) for use with the ABAQUS/Standard1 nonlinear finite element analysis tool. The failure initiation criteria include the maximum stress criteria, maximum strain criteria, the Tsai-Wu failure polynomial, and the Hashin criteria. The material degradation model is based on the ply-discounting approach where the local material constitutive coefficients are degraded. Applications and extensions of the progressive failure analysis material model address two-dimensional plate and shell finite elements and three-dimensional solid finite elements. Implementation details and use of the UMAT subroutine are described in the present paper. Parametric studies for composite structures are discussed to illustrate the features of the progressive failure modeling methods that have been implemented.

Knight, Norman F. Jr.↗
NASA NTRS2026

LAI Assimilation Schedules to Constrain Uncertain Cultivars and Soils in CERES-Maize

In anticipation of large-domain crop model applications where precise local configuration and calibration is not possible, we describe benefits and potential drawbacks of employing a crop pest module to achieve leaf area index (LAI) assimilation into a high performing CERES-Maize crop model configuration at a field experiment site in Perry, Iowa. Simulation experiments explore the use of MODIS satellite-derived LAI to constrain and adjust LAI to counter imprecise cultivar and soil configurations often occurring in the absence of high-quality local information. Simulations using single-day, window, and continuous LAI replacement across 5 cultivars and 2 soil calibration approaches for nine corn rotation years from 2004 to 2020 led to different yield outcomes and reverberations throughout the field environment. Evaluating variance and mean bias, results indicate minimal interventions in early vegetative and grain-filling stages were more beneficial than use of continuous LAI adjustments, as they minimized disruptions to the internal resource balances governing plant stresses and grain production. LAI adjustment was particularly helpful in constraining growth related to uncertain thermal unit requirements and leaf tip appearance rates (P1 and PHINT cultivar parameters, respectively). Findings underscore the need to assimilate additional state variables to ensure internal biophysical coherence. This approach shows promise for applications spanning wider domains with prediction time pressure where detailed configuration, more complex assimilation methods, or recalibration of crop model parameters may not be practical.

Phenology↗
NASA NTRS1990

Development of impact design methods for ceramic gas turbine components

Impact damage prediction methods are being developed to aid in the design of ceramic gas turbine engine components with improved impact resistance. Two impact damage modes were characterized: local, near the impact site, and structural, usually fast fracture away from the impact site. Local damage to Si3N4 impacted by Si3N4 spherical projectiles consists of ring and/or radial cracks around the impact point. In a mechanistic model being developed, impact damage is characterized as microcrack nucleation and propagation. The extent of damage is measured as volume fraction of microcracks. Model capability is demonstrated by simulating late impact tests. Structural failure is caused by tensile stress during impact exceeding material strength. The EPIC3 code was successfully used to predict blade structural failures in different size particle impacts on radial and axial blades.

Song, J.↗
NASA NTRS1974

Calculation of turbulent shear stress in supersonic boundary layer flows

An analysis of turbulent boundary layer flow characteristics and the computational procedure used are discussed. The integrated mass and momentum flux profiles and differentials of the integral quantities are used in the computations so that local evaluation of the streamwise velocity gradient is not necessary. The computed results are compared with measured shear stress data obtained by using hot wire anemometer and laser velocimeter techniques. The flow measurements were made upstream and downstream of an adiabatic unseparated interaction of an oblique shock wave with the turbulent boundary layer on the flat wall of a two dimensional wind tunnel. A comparison of the numerical analysis and actual measurements is made and the effects of small differences in mean flow profiles on the computed shear stress distributions are discussed.

Sun, C. C.↗
NASA NTRS1990

Iosipescu shear characterization of polymeric and metal matrix composites

The paper outlines the results of an experimental/analytical investigation of deformation and stress fields in the test section of graphite/epoxy and boron/aluminum Iosipescu specimens using high-resolution moire interferometry and finite-element analysis. Very good correlation between experiment and analysis is demonstrated for the graphite/epoxy specimens. This, in turn, justifies the use of correction factors employed in the previous investigations to eliminate differences in the initial shear moduli obtained from 0-deg and 90-deg Iosipescu specimens. The significant discrepancies observed in the case of boron/aluminum specimens are attributed to localized yielding at the notch tips caused by the clamping action, the associated in-plane bending effects and changing boundary conditions due to slippage in the grips.

Pindera, M.-J.↗
DOE OSTI2024

Historical Subsurface Cooling in the Tropical Pacific and Its Dynamics

Abstract Understanding how the tropical Pacific responds to rising greenhouse gases in recent decades is of paramount importance given its central role in global climate systems. Extensive research has explored the long-term trends of tropical Pacific sea surface temperatures (SSTs) and the overlying atmosphere, yet the historical change in the upper ocean has received far less attention. Here, we present compelling evidence of a prominent subsurface cooling pattern along the thermocline in the central-to-eastern tropical Pacific since 1958. This subsurface cooling has been argued to be contributing to the observed cooling or lack of warming of the equatorial cold tongue SST. We further demonstrate that different mechanisms are responsible for different parts of the subsurface cooling. In the central-to-eastern equatorial Pacific and the southeastern off-equatorial Pacific, where zonal wind stress strengthens, a pronounced subsurface cooling trend emerges just above the thermocline that is closely tied to increased Ekman pumping. In the eastern equatorial Pacific where zonal wind stress weakens, the westward surface current and eastward Equatorial Undercurrent weaken as well, resulting in reduced vertical current shear and increased ocean stability, which suppresses vertical mixing and leads to local cooling. We conclude that the historical subsurface cooling is primarily linked to dynamical adjustments of ocean currents to tropical surface wind stress changes.

Jiang, Feng↗
NASA NTRS2024

Hawai’i Climate: Utilizing NASA Earth Observations to Assess Thermal Stress Impacts on Coastal Hawaiian Fishponds

Hawaiian fishponds (loko i’a), ancient aquaculture systems created over 1500 years ago, historically served to harvest prized native fish species. Colonization-induced changes in land management have reduced active loko i’a from over 500 to ~20, destabilizing nearshore ecosystems and diminishing native Hawaiian cultural traditions. In recent decades, efforts to restore loko i’a have intensified as part of a broader emphasis on preserving Hawaiian cultural identity. Anthropogenic-induced climate change threatens these systems, particularly through the increase in frequency and magnitude of marine heatwaves. Though largely unexplored, these events are effectively studied utilizing remote sensing via Earth observations. Leveraging this technique enables autonomous extraction of water quality parameters to inform fishpond management. This study examined seven loko i’a across the main Hawaiian Islands (Hawai’i Island, Moloka’i, Maui, O’ahu, Kaua’i). Sea surface temperature and normalized difference chlorophyll index measurements were derived in fishponds utilizing Landsat 8 Thermal Infrared Sensor Thermal Infrared Sensor and Sentinel-2 MultiSpectral Instrument between 2013 to 2023, while regional patterns in temperature and chlorophyll-a were observed via Aqua’s Moderate Resolution Imaging Spectroradiometer (2002-2021). Employing these sensors, we focused analysis around the 2015 and 2019 marine heatwaves in the Pacific. We observed a heat stress gradient across the archipelago, creating thermally sheltered and unsheltered fishponds, yet identified no notable localized relationships between temperature and chlorophyll-a. Additionally, we applied pre-determined thermal maxima for culturally significant fish species (Striped Mullet: Mugil cephalus; Milkfish: Chanos Chanos) to assess historical habitat suitability within fishponds. While the direct impacts of marine heatwaves on loko i’a ecology are unknown, temperature patterns observed in this study are cause for concern.

Raz Wachtel↗
DOE OSTI2025

10-th order of accuracy for numerical solution of 3-D elasticity equations for heterogeneous materials on unfitted Cartesian meshes

We have developed the Optimal Local Truncation Error Method (OLTEM) with 10-th order of accuracy on unfitted Cartesian meshes for a system of 3-D elasticity equations with smooth irregular interfaces. 5 x 5 x 5 = 125-point stencils (similar to those for quadratic finite elements) for elastic heterogeneous materials are used for OLTEM. There are no unknowns at the interface points between different materials; the structure of the global discrete equations is the same for homogeneous and heterogeneous materials. The calculation of unknown stencil coefficients is based on the minimization of the local truncation error of the stencil equations and yields the optimal 10-th order of accuracy for OLTEM on unfitted Cartesian meshes, i.e., the increase by 7 orders in accuracy compared to quadratic finite elements on conformal meshes. A new post-processing procedure provides the 9-th order of accuracy for stresses in the 3-D case. Similar to basic computations it uses OLTEM with the 125-point stencils, the interface conditions and the elasticity equations. It was shown that the use of the elasticity equations for post-processing improves the accuracy of 0.1% stresses by 6 orders compared to post-processing without the use of PDEs. At an accuracy of for stresses, OLTEM with the new post-processing procedure reduces the number of degrees of freedom by 360 - 8000 times compared to quadratic finite elements with similar stencils. OLTEM with the 125-point stencils yields even more accurate results than high-order finite elements with much wider stencils. OLTEM provides accurate numerical results for compressible and nearly incompressible materials.

elasticity equations↗
NASA NTRS2012

An Empirical Temperature Variance Source Model in Heated Jets

An acoustic analogy approach is implemented that models the sources of jet noise in heated jets. The equivalent sources of turbulent mixing noise are recognized as the differences between the fluctuating and Favre-averaged Reynolds stresses and enthalpy fluxes. While in a conventional acoustic analogy only Reynolds stress components are scrutinized for their noise generation properties, it is now accepted that a comprehensive source model should include the additional entropy source term. Following Goldstein s generalized acoustic analogy, the set of Euler equations are divided into two sets of equations that govern a non-radiating base flow plus its residual components. When the base flow is considered as a locally parallel mean flow, the residual equations may be rearranged to form an inhomogeneous third-order wave equation. A general solution is written subsequently using a Green s function method while all non-linear terms are treated as the equivalent sources of aerodynamic sound and are modeled accordingly. In a previous study, a specialized Reynolds-averaged Navier-Stokes (RANS) solver was implemented to compute the variance of thermal fluctuations that determine the enthalpy flux source strength. The main objective here is to present an empirical model capable of providing a reasonable estimate of the stagnation temperature variance in a jet. Such a model is parameterized as a function of the mean stagnation temperature gradient in the jet, and is evaluated using commonly available RANS solvers. The ensuing thermal source distribution is compared with measurements as well as computational result from a dedicated RANS solver that employs an enthalpy variance and dissipation rate model. Turbulent mixing noise predictions are presented for a wide range of jet temperature ratios from 1.0 to 3.20.

Khavaran, Abbas↗
DOE OSTI2026

A Data-Driven Method for Synthetic Extreme Weather Generation and Solar Impact Assessment: Preprint

High-resolution, high-fidelity weather datasets are essential for testing and evaluating the resilience of power systems, particularly under extreme weather conditions. However, existing extreme weather datasets are typically derived from historical events that are localized and may lack the spatial and temporal resolution or scenario diversity needed to test largescale power systems. In this work, we propose a synthetic extreme weather simulation approach capable of generating targeted extreme events, such as hurricanes, using publicly available data sources. Preliminary results demonstrate the impact of a simulated Category 1 hurricane on renewable generation and critical infrastructure in California. The work aims to provide a flexible approach for creating multiple types of extreme weather scenarios across different regions, enabling comprehensive system stress testing, training, and resilience assessment.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

In situ quantification of fracture slip induced by hydraulic injections in a deep borehole: A comparison of two different borehole techniques

In situ measurements of fracture deformation during fluid injection are rare, yet essential for understanding the mechanical response of fractured rock. In this study, we evaluate the reliability of two methods by comparing their slip vector estimates: high-resolution borehole acoustic televiewer images captured before and after injection tests, and displacement data from a three-component borehole deformation probe recorded during the injections. Acoustic televiewer images capture only final in-plane displacement, whereas three-component borehole deformation measurements provide full 3D, transient fracture movement. Four injection tests in a fractured granitic rock mass along an inclined borehole at the Bedretto Underground Laboratory (Switzerland), beneath more than 1100 m of overburden, were analyzed. The two methods yielded consistent kinematics and comparable slip magnitudes, typically in the range of 0.2 – 0.6 mm. Angular differences between estimated slip directions ranged from 0.2° to 28°. For two test intervals with pre-existing fractures optimally oriented for slip, the measured slip vectors agreed with forward-modeled slip vectors based on the regional in situ stress field, supporting the Wallace–Bott hypothesis. Measurable shear displacement was identified in two test intervals containing fractures oriented nearly perpendicular to the least principal stress, including a newly formed hydraulic fracture. This demonstrates that even fractures unfavorably oriented for slip can undergo shear due to local stress perturbations around an inclined borehole.

Acoustic televiewer↗
NASA NTRS2001

Statistical Ensemble of Large Eddy Simulations

A statistical ensemble of large eddy simulations (LES) is run simultaneously for the same flow. The information provided by the different large scale velocity fields is used to propose an ensemble averaged version of the dynamic model. This produces local model parameters that only depend on the statistical properties of the flow. An important property of the ensemble averaged dynamic procedure is that it does not require any spatial averaging and can thus be used in fully inhomogeneous flows. Also, the ensemble of LES's provides statistics of the large scale velocity that can be used for building new models for the subgrid-scale stress tensor. The ensemble averaged dynamic procedure has been implemented with various models for three flows: decaying isotropic turbulence, forced isotropic turbulence, and the time developing plane wake. It is found that the results are almost independent of the number of LES's in the statistical ensemble provided that the ensemble contains at least 16 realizations.

Carati, Daniele↗
NASA NTRS1951

Theory of thin-walled rods

Starting with the Love equations for bending of extensible shells, "principal stress states" are sought for a thin-walled rod of arbitrary but open cross section. Principal stress states exclude those local states arising from end conditions which damp out with distance from the ends. It is found that for rods of intermediate length, long enough to avoid local bending at a support, and short enough that elementary torsion and bending are not the most significant stress states, four principal states exist. Three of these states are associated with the planar distribution of axial stress and are equivalent to the engineering theory of extension and bending of solid sections. The fourth state resembles that which has been called in the literature "bending stress due to torsional", except that cross sections are permitted to bend and the shear along the center line of the cross section is permitted to differ from zero.

LOADS AND STRESSES, STRUCTURAL - TENSION↗
NASA NTRS2013

DNS, Enstrophy Balance, and the Dissipation Equation in a Separated Turbulent Channel Flow

The turbulent flows through a plane channel and a channel with a constriction (2-D hill) are numerically simulated using DNS and RANS calculations. The Navier-Stokes equations in the DNS are solved using a higher order kinetic energy preserving central schemes and a fifth order accurate upwind biased WENO scheme for the space discretization. RANS calculations are performed using the NASA code CFL3D with the komega SST two-equation model and a full Reynolds stress model. Using DNS, the magnitudes of different terms that appear in the enstrophy equation are evaluated. The results show that the dissipation and the diffusion terms reach large values at the wall. All the vortex stretching terms have similar magnitudes within the buffer region. Beyond that the triple correlation among the vorticity and strain rate fluctuations becomes the important kinematic term in the enstrophy equation. This term is balanced by the viscous dissipation. In the separated flow, the triple correlation term and the viscous dissipation term peak locally and balance each other near the separated shear layer region. These findings concur with the analysis of Tennekes and Lumley, confirming that the energy transfer terms associated with the small-scale dissipation and the fluctuations of the vortex stretching essentially cancel each other, leaving an equation for the dissipation that is governed by the large-scale motion.

Balakumar, Ponnampalam↗
DOE OSTI2024

A Scale‐Dependent Analysis of the Barotropic Vorticity Budget in a Global Ocean Simulation

Abstract The climatological mean barotropic vorticity budget is analyzed to investigate the relative importance of surface wind stress, topography, planetary vorticity advection, and nonlinear advection in dynamical balances in a global ocean simulation. In addition to a pronounced regional variability in vorticity balances, the relative magnitudes of vorticity budget terms strongly depend on the length‐scale of interest. To carry out a length‐scale dependent vorticity analysis in different ocean basins, vorticity budget terms are spatially coarse‐grained. At length‐scales greater than 1,000 km, the dynamics closely follow the Topographic‐Sverdrup balance in which bottom pressure torque, surface wind stress curl and planetary vorticity advection terms are in balance. In contrast, when including all length‐scales resolved by the model, bottom pressure torque and nonlinear advection terms dominate the vorticity budget (Topographic‐Nonlinear balance), which suggests a prominent role of oceanic eddies, which are of km in size, and the associated bottom pressure anomalies in local vorticity balances at length‐scales smaller than 1,000 km. Overall, there is a transition from the Topographic‐Nonlinear regime at scales smaller than 1,000 km to the Topographic‐Sverdrup regime at length‐scales greater than 1,000 km. These dynamical balances hold across all ocean basins; however, interpretations of the dominant vorticity balances depend on the level of spatial filtering or the effective model resolution. On the other hand, the contribution of bottom and lateral friction terms in the barotropic vorticity budget remains small and is significant only near sea‐land boundaries, where bottom stress and horizontal viscous friction generally peak.

Khatri, Hemant↗
DOE OSTI2024

Direct numerical simulation of open-channel flow over a heterogeneous particle bed at low relative submergence

In this article, this study investigates turbulent open-channel flows over beds of irregularly arranged particles, using direct numerical simulations at a friction Reynolds number of Re T = 300⁠. Two distinct cases are examined: a polydisperse bed (P800) composed of multiple layers of randomly distributed spheres of varying sizes, and a monodisperse bed (M1015) formed by a random distribution of uniform sized spheres, with a bottommost single layer of varied-sized particles to introduce realistic randomness. Our investigation unveils a rich network of low- and high-speed streaks within the flow field, exhibiting distinctive behaviors in different bed configurations. The P800 case presents a poorly organized flow pattern induced by the varied particle sizes and arrangements, while the M1015 case shows a more regular flow pattern, marked by larger streaks. We also observe that total wall shear stress is substantially influenced by surface roughness-induced drag, extending beyond the effects documented in existing studies of open-channel flows. The present study reveals intricate secondary flow patterns over irregular particle beds. Large-scale circulations are discerned around particle crests in the P800 case and localized circulations with increased turbulence in the M1015 case. Furthermore, analysis of Reynolds stress tensor components indicates that roughness disrupts coherent turbulent eddies, consequently mitigating peak stress. We quantify correlations between drag force and local fluid velocity fluctuations. Notably, a larger deviation in drag is observed in the P800 case compared to M1015, accentuating the influence of particle size and distribution on fluid–particle interactions.

42 ENGINEERING↗
NASA NTRS1995

Microcracking in composite laminates under thermal and mechanical loading

Composites used in space structures are exposed to both extremes in temperature and applied mechanical loads. Cracks in the matrix form, changing the laminate thermoelastic properties. The goal of the present investigation is to develop a predictive methodology to quantify microcracking in general composite laminates under both thermal and mechanical loading. This objective is successfully met through a combination of analytical modeling and experimental investigation. In the analysis, the stress and displacement distributions in the vicinity of a crack are determined using a shear lag model. These are incorporated into an energy based cracking criterion to determine the favorability of crack formation. A progressive damage algorithm allows the inclusion of material softening effects and temperature-dependent material properties. The analysis is implemented by a computer code which gives predicted crack density and degraded laminate properties as functions of any thermomechanical load history. Extensive experimentation provides verification of the analysis. AS4/3501-6 graphite/epoxy laminates are manufactured with three different layups to investigate ply thickness and orientation effects. Thermal specimens are cooled to progressively lower temperatures down to -184 C. After conditioning the specimens to each temperature, cracks are counted on their edges using optical microscopy and in their interiors by sanding to incremental depths. Tensile coupons are loaded monotonically to progressively higher loads until failure. Cracks are counted on the coupon edges after each loading. A data fit to all available results provides input parameters for the analysis and shows them to be material properties, independent of geometry and loading. Correlation between experiment and analysis is generally very good under both thermal and mechanical loading, showing the methodology to be a powerful, unified tool. Delayed crack initiation observed in a few cases is attributed to a lack of preexisting flaws assumed by the analysis. Some interactions between adjacent ply groups are attributed to local stress concentrations. These two effects are not captured by the analysis due to its global nature. The analysis is conservative in these cases and agrees well with data after the observed onset of cracking.

Maddocks, Jason R.↗