Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stress Modeling”

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 1,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Micromechanical predictions of crack propagation and fracture energy in a single fiber boron/aluminum model composite

The axisymmetric finite element model and associated computer program developed for the analysis of crack propagation in a composite consisting of a single broken fiber in an annular sheath of matrix material was extended to include a constant displacement boundary condition during an increment of crack propagation. The constant displacement condition permits the growth of a stable crack, as opposed to the catastropic failure in an earlier version. The finite element model was refined to respond more accurately to the high stresses and steep stress gradients near the broken fiber end. The accuracy and effectiveness of the conventional constant strain axisymmetric element for crack problems was established by solving the classical problem of a penny-shaped crack in a thick cylindrical rod under axial tension. The stress intensity factors predicted by the present finite element model are compared with existing continuum results.

Adams, D. F.↗
NASA NTRS1984

On the Influence of Turbulent Motions on Non-radial Oscillations

The effect of turbulent motions on oscillations is studied, considering only the coupling between turbulent and oscillatory velocities. In this case, the turbulence affects the oscillations through the Reynolds stresses in the momentum equation for the pulsations. A simple model of turbulence is adopted to evaluate these Reynolds stresses and the perturbed eigenfrequencies are expressed as a function of certain averages of the turbulent velocities.

Durney, B. R.↗
NASA NTRS1995

System Identification of Damped Truss-Like Space Structures

A spacecraft payload flown on a launch vehicle experiences dynamic loads. The dynamic loads are caused by various phenomena ranging from the start-up of the launch vehicle engine to wind gusts. A spacecraft payload should be designed to meet launch vehicle dynamic loads. One of the major steps taken towards determining the dynamic loads is to correlate the finite element model of the spacecraft with the test results of a modal survey test. A test-verified finite element model of the spacecraft should possess the same spatial properties (stiffness, mass, and damping) and modal properties (frequencies and mode shapes) as the test hardware representing the spacecraft. The test-verified and correlated finite element model of the spacecraft is then coupled with the finite element model of the launch vehicle for analysis of loads and stress. Modal survey testing, verification of a finite element model, and modification of the finite element model to match the modal survey test results can easily be accomplished if the spacecraft structure is simple. However, this is rarely the case. A simple structure here is defined as a structure where the influence of nonlinearity between force and displacement (uncertainty in a test, for example, with errors in input and output), and the influence of damping (structural, coulomb, and viscous) are not pronounced. The objective of this study is to develop system identification and correlation methods with the focus on the structural systems that possess nonproportional damping. Two approaches to correct the nonproportional damping matrix of a truss structure were studied, and have been implemented on truss-like structures such as the National Aeronautics and Space Administration's space station truss. The results of this study showed nearly 100 percent improvement of the correlated eigensystem over the analytical eigensystem. The first method showed excellent results with up to three modes used in the system identification process. The second method could handle more modes, but required more computer usage time, and the results were less accurate than those of the first method.

Armand, Sasan↗
NASA NTRS2000

Modeling of Horns for Sonic/Ultrasonic Applications

JPL has a requirement for telerobotic tools for planetary sample acquisition, which require low power and have the ability to work in harsh environments. We are currently investigating the possibility of using ultrasonic horns to develop a family of ultrasonic tools for these environments. In an effort to determine control parameters a one-dimensional Mason's model for a stepped ultrasonic horn assembly was developed which includes the effects of mechanical and electrical losses in the piezoelectric material and acoustic elements. The model is separated into three regions; the piezoelectric stack including stress bolt the backing layer and the horn. The model is found to predict the impedance data of the horn assembly very accurately up to the first coupled (radial) resonance. The model also allows for the calculation of the velocity and force and power delivered to each acoustic element. FEM modeling and accelerometer data from the horn tip were used to corroborate the model. The difficulties associated with modeling the load impedance of various devices will be discussed and current directions noted.

Sherrit, Stewart↗
DOE OSTI2026

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗
DOE OSTI2020

Technical Performance of Refractory Liners for Molten Chloride Salt Thermal Energy Storage Systems

A chloride-based molten-salt system that uses a ternary blend of MgCl2/KCl/NaCl is investigated to provide higher temperature thermal energy storage capability. Despite higher thermal stability, molten chlorides present several unique challenges, including the design of internal refractory-ceramic liners to prevent the corrosion and thermal stress of alloy tank shells. This work discusses issues and potential solutions related to containment of molten chloride salt, specifically the optimization of the refractory material at the molten salt interface (hot face). The down-selected hot face candidate limits permeation of salt through the material and forms a highly stable secondary surface phase in equilibrium with the molten salt. A mortar is created using the corrosion resistant hot face brick. Brick and mortar composite are subjected to mechanical stress/strain analysis, in order to calculate composite material properties and better inform thermomechanical models. The U.S. Department of Energy Generation 3 (DOE Gen3) program seeks to develop higher efficiency CSP plants that can provide cost-competitive, flexible power in the U.S. electric grid. The proposed Gen3 Liquid Pathway CSP plant closely resembles the configuration of current nitrate salt power towers with two-tank storage (Gen2). The differences between Gen2 and Gen3 include the types of compatible materials used in salt storage tank construction. Stainless steel loses strength at Gen3 temperatures, and although nickel superalloys would be capable of withstanding sustained high temperatures, these materials are prohibitively expensive at scale. Uninsulated tank shells also pose a significant risk as common steels are highly susceptible to chemical attack from molten chloride salt. To address these concerns, refractory-ceramic based containment materials are proposed to line the inside of the hot and cold storage tanks. In doing so, stainless or carbon steel shells may be used in construction depending on the level of insulation provided. The composition of the internal liner requires careful consideration to maximize the efficacy of multiple parameters including corrosion resistance, strength at operating temperature, durability, and cost. This is particularly true for the material at the interface with the salt, known as the "hot face", which is responsible for protecting the insulating layers between the tank shell and the hot face brick layer. The molten salt in this system is superheated over 300 °C above its freeze temperature. Therefore, unlike other industrial processes which use refractory-lined vessels, it is not expected that a freeze plane will develop in the hot face. Therefore, the hot face must be designed to withstand chemical corrosion and inhibit permeation of molten salt into the insulating layers. A down selection was performed to identify a hot face candidate best equipped to maintain thermal, mechanical, and chemical integrity when exposed to molten salt over extended periods of time. Long-duration chemical capability experiments were conducted with the down selected hot face refractory fully immersed in molten chloride salt for up to 3000 hours. The average salt penetration does not exceed 100 microns. When extrapolated to 20 and 30 years of continuous exposure, the expected salt penetration depth is approximately 2.0 mm and 2.9 mm, respectively. A magnesium-rich secondary phases develops at the salt/refractory interface. X-ray diffraction identifies the material as forsterite (Mg2SiO4), which is reported to form synthetically in molten chloride salt solutions. These results suggest the selected hot face will adequately inhibit salt permeation. While there is optimism that the hot face brick will inhibit salt permeation, mortar joints are typically the weakest point of a refractory brick liner. From a thermochemical perspective, differing thermal expansion coefficients may result in the mortar and brick to grow independent of each other, creating gaps through which molten salt can penetrate. To address this issue, NREL has developed an in-house mortar composed of the down selected hot face brick that has been shown to be compatible with the salt. Compressive stress/strain analyses have been performed on the brick/ mortar composites to generate stress/strain curves. Modulus of elasticity and Poisson's ratio of the composite may be calculated from the stress/strain curves, in order to provide more representative data to finite element mechanical models for accurate approximation of stress on the tank shell and the amount of thermal expansion expected within the tank liner.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2023

A constitutive model for sheared dense fiber suspensions

Here we propose a constitutive model to predict the viscosity of fiber suspensions, which undergoes shear thinning, at various volume fractions, aspect ratios, and shear stresses/rates. We calibrate the model using the data from direct numerical simulation and prove the accuracy by predicting experimental measurements from the literature. We use a friction coefficient decreasing with the normal load between the fibers to quantitatively reproduce the experimentally observed shear thinning in fiber suspensions. In this model, the effective normal contact force, which is directly proportional to the bulk shear stress, determines the effective friction coefficient. A rise in the shear stress reduces the effective friction coefficient in the suspension. As a result, the jamming volume fraction increases with the shear stress, resulting in a shear thinning in the suspension viscosity. Moreover, we extend the model to quantify the effects of fiber volume fraction and aspect ratio in the suspension. We calibrate this model using the data from numerical simulations for the rate-controlled shear flow. Once calibrated, we show that the model can be used to predict the relative viscosity for different volume fractions, shear stresses, and aspect ratios. The model predictions are in excellent agreement with the available experimental measurements from the literature. The findings of this study can potentially be used to tune the fiber size and volume fraction for designing the suspension rheology in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Optimization of overhang printing in ceramic binder jetting additive manufacturing

Cylindrical-shaped horizontally-hanging overhangs with different heights on a solid cubic body were printed via binder jetting additive manufacturing with sililcon carbide powder feedstock. Cracking on the overhang neck and sagging on the overhang body were observed both pre- and post-curing, especially on printed design with a smaller particle size of 10 μm. Here, it was found that the width of cracking and depth of sagging were dependent on the height relative to the solid body in the powder bed, and elevated position of overhang with large volume of underneath loose powder led to increased extend of cracking and sagging. Overhang diameter and weight were varied, which were found with no apparent effects on its deformation extend. To mitigate cracking and sagging, contactless support was designed and tested, which could significantly improve the overhang stability and avoid any overhang deformations. Support gap and support distance were studied for their effects on support capability. The printing processes with and without contactless support were simulated by a finite element modeling approach to understand the stress distribution and evolution of the overhang structures, which revealed stress concentrations on the upper neck point of the designed overhang.

Binder jetting↗
DOE OSTI2023

STRESS AND FATIGUE ANALYSIS OF AIR-TO-REFRIGERANT HEAT EXCHANGERS WITH NON-ROUND TUBE SHAPES

Air-to-refrigerant heat exchangers are fundamental components in Heating, Ventilation, Air Conditioning, and Refrigeration (HVAC&R) systems. Recent research has shown that heat exchangers utilizing small hydraulic diameter, non-round, shape-optimized tubes to be promising next design update. Such tubes can improve performance, reduce material and lower refrigerant charge. However, such tube designs may be more susceptible to fatigue issues. Therefore, stress and fatigue analyses of such heat exchangers are critical to ensure structural integrity, reliability, and manufacturability. In the present work, a mechanical model is developed to obtain the stress distributions in heat exchangers with round, ellipse, and airfoil-shaped tubes. The results indicate that under 3.45 MPa internal pressure, the maximum stresses don't surpass 190 MPa.

Zhang, Mingkan↗
NASA NTRS1987

Dislocation dynamics and the viscoplastic buckling of dendritic web type silicon ribbon

The effect of dendrites (reinforced edges) on the residual stresses, dislocation densities and buckling behavior during growth of web type silicon ribbon is studied. A viscoplastic material response function (Haasen-Sumino model) is used to calculate the stresses and the disloction density at each point in the silicon ribbon. In addition, the role of dendrites on the viscoplastic buckling behavior of the ribbon is investigated. The critical thicknesses, the corresponding deflection shapes and lateral deflection speeds are calculated. These results are then compared with similar data obtained for flat plates.

Tsai, C. T.↗
NASA NTRS1974

Design specification for the whole-body algorithm

The necessary requirements and guidelines for the construction of a computer program of the whole-body algorithm are presented. The minimum subsystem models required to effectively simulate the total body response to stresses of interest are (1) cardiovascular (exercise/LBNP/tilt); (2) respiratory (Grodin's model); (3) thermoregulatory (Stolwijk's model); and (4) long-term circulatory fluid and electrolyte (Guyton's model). The whole-body algorithm must be capable of simulating response to stresses from CO2 inhalation, hypoxia, thermal environmental exercise (sitting and supine), LBNP, and tilt (changing body angles in gravity).

Fitzjerrell, D. G.↗
NASA NTRS1995

The ocean response to operational westerly wind bursts during the 1991 - 1992 El Nino

Numerical simulations of the remotely forced ocean response to westerly wind bursts prior to and during the 1991 - 1992 El Nino are examined; the models are forced by wind stress from the U.S. Navy's atmospheric global operational analysis/forecast system. The study focuses on (1) the relative response of the first and second internal modes to a single episode of westerly bursts; (2) the role of westerly bursts in producing the eastern Pacific sea level variations from October 1990 to February 1992; and (3) the relative importance of the remotely forced sea level responses generated by central and western Pacific wind anomalies. The simulation use the Naval Research Laboratory global multilayer formulation; the suite of experiments includes hydrodynamic simulations that use both one- and three- active-layer reduced gravity models as well as an experiment that also includes thermodynamic effects. The models are executed from January 1, 1990, to March 1, 1992, a period that includes 10 significant westerly wind bursts or burst clusters and the 1991 - 1992 El Nino event. The numerical experiments reveal an ability to accurately simulate the eastern Pacific sea level variations during this period. In response to a single burst, the three-layer hydrodynamic simulation reveals that the second internal mode Kelvin wave yields a sea level change at the eastern boundary that is approximately one third that of the first mode and a surface velocity signature that is equivalent to the first mode. During the onset of the El Nino event the inclusion of higher modes also produces a more realistic representation of the observed eastern boundary sea level signal. Furthermore, by comparing the model response to particular wind bursts with the observed sea level at Baltra, Galapagos, a value of 2.5 - 2.6 m/s is suggested as the most appropriate mean speed for the first internal mode Kelvin wave during the onset phase. A scenario for the onset of the 1991 - 1992 El Nino is given.

Kindle, John C.↗
DOE OSTI2025

Ponderosa pine hydraulic stress predicts more extreme wildfire behavior under future conditions in Bandelier National Monument, New Mexico

Background Live fuel moisture contributes to wildfire spread and reflects plant stress and physiological traits. The anticipated change in live fuel moisture under future conditions is likely non-linear, owing to physiological plant thresholds in water hydraulics. We constructed a mechanistic model of live fuel moisture’s response to water stress to understand the impact of future climate on live fuel moisture. We first gathered data on plant physiology and live fuel moisture for Pinus ponderosa at Bandelier National Monument, NM, USA, and modeled their relationship. We then parameterized a mechanistic plant hydrodynamics model (FATES-HYDRO) to simulate changes in plant stress and a statistical model to simulate the resulting impact on live fuel moisture. We then simulated FATES-HYDRO under future climate anomalies (SSP2-4.5 and SSP5-8.5: 2080–2100) to understand the change in plant stress and estimate its impact on live fuel moisture. Results We found that the number of days below crucial thresholds of live fuel moisture (100% and 79%) increased from contemporary levels (< 100%: 72 days, < 79%: 1.4 days) under SSP2-4.5 (< 100%: 185 days, < 79%: 10.2 day) and increased exponentially under SSP5-8.5 (< 100%: 215 days, < 79%: 65 days). We found that gross primary productivity decreased under both future climate scenarios (contemporary: 336 g C m −2 , SSP2-4.5: 203 g C m −2 , SSP5-8.5: 243 g C m −2 ); however, spring productivity increased under SSP5-8.5, potentially altering fuel loading. We additionally see a potentially lethal loss of conductivity in hydraulic P. ponderosa under SSP5-8.5. Conclusions The overall increase in plant water stress (as represented by loss of hydraulic conductivity and leaf water potential) and lower live fuel moisture appear to be driven by reduced precipitation during late summer monsoons typical of the region, extending the fire season later in the year. We predict increasing variability in the P. ponderosa wildfire regime under both climate projections driven by changing productivity, rising mortality, and an overall decrease in live fuel moisture.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Progress Report on Model Development for Aerosol Transport through Divergent Cracks Paths

This report summarizes the progress in developing a phenomenological model of aerosol transport, deposition, and plugging through microchannels. The purpose of this effort is to introduce to a user community—involving researchers, regulators, and industry—a generic, reliable numerical model for the prediction of aerosol transport while accounting for potential deposition and plugging of the leak paths to model spent nuclear fuel (SNF) release from postulated Stress Corrosion Cracks in canisters under storage or transportation. Current work focuses on expanding the model to predict aerosol (and gas) flow through complex microchannel (nozzle) geometries in addition to the rectangular and cylindrical geometries, as we approach more realistic Stress Corrosion Crack conditions. In this regard, a divergent nozzle geometry, used at Sandia National Laboratories (SNL) for testing aerosol release and retention (Durbin et al. 2021), was added to the model. The model was then validated with blowdown data for the particular microchannel, from experiments conducted at SNL. The report also presents preliminary non-benchmarked aerosol penetration fraction and mass flow rate for a monodisperse 10-micron (AED) particle concentration of 1.7e-08 kg/m 3 being released through the divergent microchannel geometry from the SNL’s aerosol experimental tank setup.Future work will involve validating the aerosol model and updating the model’s Graphical User Interface (GUI) to include the divergent nozzle geometry. This is expected to help stakeholders perform quick and easy first principles calculations without needing to understand the underlying MATLAB script.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Tracking Seasonal and Interannual Variability in Photosynthetic Downregulation in Response to Water Stress at a Temperate Deciduous Forest

Abstract The understanding and modeling of photosynthetic dynamics affected by climate variability can be highly uncertain. In this paper, we examined a well‐characterized eddy covariance site in a drought‐prone temperate deciduous broadleaf forest combining tower measurements and satellite observations. We find that an increase in spring temperature usually leads to enhanced spring gross primary production (GPP), but a GPP reduction in late growing season due to water limitation. We evaluated how well a coupled fluorescence‐photosynthesis model (SCOPE) and satellite data sets track the interannual and seasonal variations of tower GPP from 2007 to 2016. In SCOPE, a simple stress factor scaling of Vcmax as a linear function of observed predawn leaf water potential ( ψ pd ) shows a good agreement between modeled and measured interannual variations in both GPP and solar‐induced chlorophyll fluorescence (SIF) from the Global Ozone Monitoring Experiment‐2 (GOME‐2). The modeled and satellite‐observed changes in SIF yield are ~30% smaller than corresponding changes in light use efficiency (LUE) under severe stress, for which a common linear SIF to GPP scaling would underestimate the stress reduction in GPP. Overall, GOME‐2 SIF tracks interannual tower GPP variations better than satellite vegetations indices (VIs) representing canopy “greenness.” However, it is still challenging to attribute observed SIF variations unequivocally to greenness or physiological changes due to large GOME‐2 footprint. Higher‐resolution SIF data sets (e.g., TROPOMI) already show the potential to well capture the downregulation of late‐season GPP and could pave the way to better disentangle canopy structural and physiological changes in the future.

He, Liyin↗
NASA NTRS1987

Toward the large-eddy simulations of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model (in terms of Favre-filtered fields) is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The three dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96 to the third power grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. Future applications of these compressible subgrid-scale models to the large-eddy simulation of supersonic aerodynamic flows are discussed briefly.

Erlebacher, G.↗
NASA NTRS1990

Toward the large-eddy simulation of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96(exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.

Erlebacher, G.↗