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,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Factors Influencing Dwell Fatigue Cracking in Notches of Powder Metallurgy Superalloys

The influences of heat treatment and cyclic dwells on the notch fatigue resistance of powder metallurgy disk superalloys were investigated for low solvus high refractory (LSHR) and ME3 disk alloys. Disks were processed to produce material conditions with varied microstructures and associated mechanical properties. Notched specimens were first subjected to baseline dwell fatigue cycles having a dwell at maximum load, as well as tensile, stress relaxation, creep rupture, and dwell fatigue crack growth tests at 704 C. Several material heat treatments displayed a bimodal distribution of fatigue life with the lives varying by two orders-of-magnitude, while others had more consistent fatigue lives. This response was compared to other mechanical properties, in search of correlations. The wide scatter in baseline dwell fatigue life was observed only for material conditions resistant to stress relaxation. For selected materials and conditions, additional tests were then performed with the dwells shifted in part or in total to minimum tensile load. The tests performed with dwells at minimum load exhibited lower fatigue lives than max dwell tests, and also exhibited early crack initiation and a substantial increase in the number of initiation sites. These results could be explained in part by modeling evolution of peak stresses in the notch with continued dwell fatigue cycling. Fatigue-environment interactions were determined to limit life for the fatigue cycles with dwells.

Gabb, T. P.↗
NASA NTRS1992

Scanning electron acoustic microscopy of residual stresses in ceramics - Theory and experiment

The paper presents a three-dimensional mathematical model of signal generation and contrast in brittle materials and uses the model to simulate the effect of residual stress fields on the scanning electron acoustic microscopy (SEAM)-generated electron acoustic signal. According to the model, a positive (tensile) strain produces an increase in the output signal, whereas a negative (compressive) strain produces a decrease in the ouput signal. Dark field contrast conditions occur at a chopping frequency at which V2 - V1 is greater than 0 (where V2 = V is the SEAM output in a region of residual stresses, and V1 is the output in a stress-free region of the sample). Under ideal conditions (maximum contrast) V1 approaches zero. It was found that tensile strains of the order 0.2-0.3 percent, possible in brittle materials, would produce a variation of the acoustic output signal of the order 10 nV (about 1 percent), well within the image contrast and signal processing capability of the SEAM electronics.

Cantrell, John H.↗
NASA NTRS2001

Computational Analyses of Propulsion Aeroacoustics for Mixed Flow Nozzle Pylon Installation at Takeoff

A CFD analyses is presented for a set of baseline and noise suppression mixed flow nozzles with and without a pylon installation. The five model configurations are as follows; a baselinecore/fan dual-stream nozzle with an external plug, a chevron mixer nozzle with a peak on the symmetry plane with external plug, both of the above nozzles with an installed bifurcatingpylon and lastly a clocked chevron mixer nozzle such that a trough is aligned with the center of the pylon. The fluid flow is simulated by solving the asymptotically steady, compressible, Reynolds-averaged Navier-Stokes equations using an implicit, up-wind, flux-difference splitting finite volume scheme and standard two equation k-epsilon turbulence model with a linear stress representation. All computations are performed using the multiblock, parallel, structuredcode PAB3D. Results indicate that the clocked chevron with pylon case achieves the most optimal levels of average and peak turbulence kinetic energy and vorticity and therefore is expected to be the quietest of the five configurations tested. Further study is required to refine expressions which are indicative of noise and mate these with rigorous noise prediction models.

Massey, Steven J.↗
DOE OSTI2021

Embedded training of neural-network subgrid-scale turbulence models

We report that the weights of a deep neural-network model are optimized in conjunction with the governing flow equations to provide a model for subgrid-scale stresses in a temporally developing plane turbulent jet at Reynolds number Re 0 = 6000 . The objective function for training is first based on the instantaneous filtered velocity fields from a corresponding direct numerical simulation, and the training is by a stochastic gradient descent method, which uses the adjoint Navier-Stokes equations to provide the end-to-end sensitivities of the model weights to the velocity fields. In-sample and out-of-sample testing on multiple dual-jet configurations show that its required mesh density in each coordinate direction for prediction of mean flow, Reynolds stresses, and spectra is half that needed by the dynamic Smagorinsky model for comparable accuracy. The same neural-network model trained directly to match filtered subgrid-scale stresses, without the constraint of being embedded within the flow equations during the training, fails to provide a qualitatively correct prediction. The coupled formulation is generalized to train based only on mean-flow and Reynolds stresses, which are more readily available in experiments. The mean-flow training provides a robust model, which is important, though a somewhat less accurate prediction for the same coarse meshes, as might be anticipated due to the reduced information available for training in this case. The anticipated advantage of the formulation is that the inclusion of resolved physics in the training increases its capacity to extrapolate. This is assessed for the case of passive scalar transport, for which it outperforms established models due to improved mixing predictions.

42 ENGINEERING↗
NASA NTRS2000

Physical Activity and Bone Density in Women

A mathematical model of bone density regulation as a function of the daily tissue "effective" stress has been derived. Using the model, the influence of daily activity in the form of a daily loading history has been related to bone density of the calcaneus. The theory incorporates a stress exponent m to account for differences in the importance of magnitude and number of load cycles experienced during daily activity. We have derived a parameter from the model, the "Bone Density Index" (BDI). We have developed a method of collecting daily habitual loading histories using an insole force sensor interfaced to a portable digital data logger carried in a fanny pack. Our goal for this study was to determine a stress exponent, m, relating GRFz history to Calcaneal Bone Mineral Density (CBMD).

Bowley, Susan M.↗
NASA NTRS2004

Modeling the Nonlinear, Strain Rate Dependent Deformation of Shuttle Leading Edge Materials with Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites, such as the Reinforced Carbon Carbon (RCC) material used on the leading edges of the Space Shuttle. In the developed model, the differences in the tension and compression deformation behaviors have also been accounted for. State variable viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear response is independent of the stiffness in the normal directions. The developed equations have been implemented into LS-DYNA through the use of user defined subroutines (UMATs). Several sample qualitative calculations have been conducted, which demonstrate the ability of the model to qualitatively capture the features of the deformation response present in woven ceramic matrix composites.

Goldberg, Robert K.↗
DOE OSTI2024

Shape Matters: Understanding the Effect of Electrode Geometry on Cell Resistance and Chemo-Mechanical Stress

Rechargeable batteries that incorporate shaped three-dimensional electrodes have been shown to have increased power and energy densities when compared to a conventional geometry, i.e. a planar cathode and anode that sandwich an electrolyte. Electrodes can be shaped to enable a higher active material loading, while keeping ion transport distances small. However, the relationship between electrical and mechanical performance of shaped electrodes remains poorly understood. Many electrode designs have been explored, where the electrodes are individually shaped or intertwined, and advances in manufacturing and shape/topology optimization have made such designs a reality. Here, we explore sinusoidal half cells and interdigitated full cells. First, we use a simple electrostatics model to understand the cell resistance as a function of shape. Here we focus on low-temperature conditions, where the electrolyte conductivity decreases relative to that of the electrode; here, LiPF 6 EC:DMC electrolyte and MnO 2 electrode are considered. Next, we use a chemo-mechanics model to examine the stress that arises due to intercalation-driven volume expansion. We show that shaped electrodes provide a significant reduction in resistance in low-temperature conditions, however, they exhibit unfavorable stress concentrations. Overall, we find that the fully interdigitated electrodes may provide the best balance with respect to this resistance-stress trade-off.

25 ENERGY STORAGE↗
NASA NTRS2005

Navier-Stokes Computations With One-Equation Turbulence Model for Flows Along Concave Wall Surfaces

This report presents the use of a time-marching three-dimensional compressible Navier-Stokes equation numerical solver with a one-equation turbulence model to simulate the flow fields developed along concave wall surfaces without and with a downstream extension flat wall surface. The 3-D Navier- Stokes numerical solver came from the NASA Glenn-HT code. The one-equation turbulence model was derived from the Spalart and Allmaras model. The computational approach was first calibrated with the computations of the velocity and Reynolds shear stress profiles of a steady flat plate boundary layer flow. The computational approach was then used to simulate developing boundary layer flows along concave wall surfaces without and with a downstream extension wall. The author investigated the computational results of surface friction factors, near surface velocity components, near wall temperatures, and a turbulent shear stress component in terms of turbulence modeling, computational mesh configurations, inlet turbulence level, and time iteration step. The computational results were compared with existing measurements of skin friction factors, velocity components, and shear stresses of the developing boundary layer flows. With a fine computational mesh and a one-equation model, the computational approach could predict accurately the skin friction factors, near surface velocity and temperature, and shear stress within the flows. The computed velocity components and shear stresses also showed the vortices effect on the velocity variations over a concave wall. The computed eddy viscosities at the near wall locations were also compared with the results from a two equation turbulence modeling technique. The inlet turbulence length scale was found to have little effect on the eddy viscosities at locations near the concave wall surface. The eddy viscosities, from the one-equation and two-equation modeling, were comparable at most stream-wise stations. The present one-equation turbulence model is an effective approach for turbulence modeling in the near solid wall surface region of flow over a concave wall.

Wang, Chi R.↗
DOE OSTI2021

The Impala's horn applied to posterior samples of Ti-6Al-4V strength model parameters

We have generated strength model parameters for the Preston-Tonks-Wallace (PTW) model of plastic deformation based on a collection of Ti-6Al-4V alloy data and a Bayesian framework. This data collected from the literature includes different chemistries and processes. The experiments include quasistatic compression, Split Hopkinson Pressure Bar (SHPB) compression, and Taylor cylinders. Here we use posterior parameter distributions applied to a hypothetical impact referred to as the Impala's horn. The Impala's horn is like a Taylor cylinder, but modified to have a conical pro le that is truncated before coming to a point. The smaller diameter at the impact end results in higher calculated strain rates and strains than the same velocity impact with the traditional Taylor profile. The modified geometry exercises the behavior of the strength model outside the range of conditions found in the calibration data. Here we sort posterior strength models based on their flow stress at characteristic values for this hypothetical deformation. Then we compare the performance of these posterior models based on two simple observable final state properties of the Impala's horn.

36 MATERIALS SCIENCE↗
NASA NTRS2004

Modeling the Nonlinear, Strain Rate Dependent Deformation of Woven Ceramic Matrix Composites With Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites with a plain weave fiber architecture. In the developed model, the differences in the tension and compression response have also been considered. State variable based viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear stiffness is independent of the stiffness in the normal directions. The developed equations have been implemented into a commercially available transient dynamic finite element code, LS-DYNA, through the use of user defined subroutines (UMATs). The tensile, compressive, and shear deformation of a representative plain weave woven ceramic matrix composite are computed and compared to experimental results. The computed values correlate well to the experimental data, demonstrating the ability of the model to accurately compute the deformation response of woven ceramic matrix composites.

Goldberg, Robert K.↗
DOE OSTI2025

Linking Pressure to Electrochemical Evolution in Solid-State Conversion Cathode Composites

Conversion-type cathodes, such as sulfur, FeS 2 , and FeF 3 , offer high theoretical capacities in solid-state lithium batteries but are hindered by substantial volume changes during cycling, leading to interfacial contact loss, crack formation, and microstructural degradation. Here, we investigate the relationships between electrochemical, mechanical, and structural evolution in solid-state electrode composites with these three active materials. Using real-time stack-pressure monitoring, synchrotron X-ray absorption spectroscopy, and electrokinetic modeling, we elucidate how stress evolution is linked to reversible and irreversible redox reactions. Nonlinear stack pressure evolution in cells with sulfur, FeS 2 , and FeF 3 electrode composites is found to arise from material-specific volume changes, the balance of volume change between the working and counter electrode, and the formation of distinct reaction intermediates. The three materials exhibit distinct stack pressure evolution, which is closely related to the different reaction processes in the materials, as demonstrated with X-ray absorption spectroscopy measurements. Through mesoscale modeling, we relate the experimental measurements to species evolution at the particle scale and track the dynamic coexistence of intermediate phases. Our findings highlight the importance of designing for volume changes of a given active material in solid-state battery systems.

batteries↗
DOE OSTI2020

Tensile testing and fracture mechanism analysis of polyvinyl alcohol nanofibrous webs

Abstract A tensile properties testing study was conducted to understand the influence of thickness, cross‐head speed (speed of testing), gauge length (GL; specimen test length), and sample shape on important tensile properties of polyvinyl alcohol (PVA) nanofiber webs. The effects of each testing parameter on load at break, extension at break, Young's modulus, and tensile stress–strain curve of PVA nanofiber webs are analyzed. The Welch two sample t ‐tests show the significant difference among tested data. Using interaction plots, two‐way analysis of variance, and margin mean plots, the interaction effects among testing parameters have been analyzed. Of all the factors, cross‐head speed, the interaction among GL, and sample thickness (GL: Thickness) and the interaction among GL, testing speed and sample thickness (GL: Speed: Thickness) have significant influence on the tensile properties of PVA nanofiber webs. Moreover, the hypothesized model of mechanism of tensile strain–stress curve of PVA nanofiber webs has been proposed. Based on the model, the tensile strain–stress curve can be split into three stages: linear elastic, partial break up, and complete breakage. This study will provide a better understanding of tensile testing parameters' effects and their interaction effects on the tensile properties of nanowebs.

Lou, Lihua↗
NASA NTRS1981

Large area sheet task: Advanced dendritic web growth development

The growth of silicon dendritic web for photovoltaic applications was investigated. The application of a thermal model for calculating buckling stresses as a function of temperature profile in the web is discussed. Lid and shield concepts were evaluated to provide the data base for enhancing growth velocity. An experimental web growth machine which embodies in one unit the mechanical and electronic features developed in previous work was developed. In addition, evaluation of a melt level control system was begun, along with preliminary tests of an elongated crucible design. The economic analysis was also updated to incorporate some minor cost changes. The initial applications of the thermal model to a specific configuration gave results consistent with experimental observation in terms of the initiation of buckling vs. width for a given crystal thickness.

Duncan, C. S.↗
NASA NTRS2011

Investigating Some Technical Issues on Cohesive Zone Modeling of Fracture

This study investigates some technical issues related to the use of cohesive zone models (CZMs) in modeling fracture processes. These issues include: why cohesive laws of different shapes can produce similar fracture predictions; under what conditions CZM predictions have a high degree of agreement with linear elastic fracture mechanics (LEFM) analysis results; when the shape of cohesive laws becomes important in the fracture predictions; and why the opening profile along the cohesive zone length needs to be accurately predicted. Two cohesive models were used in this study to address these technical issues. They are the linear softening cohesive model and the Dugdale perfectly plastic cohesive model. Each cohesive model constitutes five cohesive laws of different maximum tractions. All cohesive laws have the same cohesive work rate (CWR) which is defined by the area under the traction-separation curve. The effects of the maximum traction on the cohesive zone length and the critical remote applied stress are investigated for both models. For a CZM to predict a fracture load similar to that obtained by an LEFM analysis, the cohesive zone length needs to be much smaller than the crack length, which reflects the small scale yielding condition requirement for LEFM analysis to be valid. For large-scale cohesive zone cases, the predicted critical remote applied stresses depend on the shape of cohesive models used and can significantly deviate from LEFM results. Furthermore, this study also reveals the importance of accurately predicting the cohesive zone profile in determining the critical remote applied load.

Wang, John T.↗
DOE OSTI2026

Low-stress creep behavior of multiple salt caverns under cyclic operations

Here, the use of multiple salt caverns for large-scale underground hydrogen storage offers a strategic solution to the growing demand for efficient and flexible energy storage, due to the inherent low permeability, high mechanical strength, and self-healing properties of the salt formation. Despite these advantages, the long-term geomechanical response of salt caverns under cyclic injection and withdrawal remains a significant challenge for ensuring operational stability. This mechanistic study highlights the critical role of low-stress creep, also known as pressure solution creep, in governing deformation of salt caverns during cyclic operations. Unlike conventional creep mechanisms that dominate under high-stress conditions, low-stress creep becomes particularly relevant under the moderate stress changes induced by repeated injection and withdrawal, which provides a previously underexplored control on cavern closure. Our three-dimensional computational model integrates the low-stress creep mechanism into the Munson–Dawson creep equation across a spectrum of operational scenarios including variations in cavern pressure cycles, injection-withdrawal frequencies, and caprock mechanical property. The low-stress creep components accelerate the cavern creep closure, which is amplified with shorter cycle frequency and larger pressure difference. Sensitivity analyses capture the complex interactions between multiple caverns and surround formations at the field scale that allows for a quantitative assessment of deformation patterns and closure rates over time and thereby informs the design and operation of sustainable underground hydrogen storage system.

Cyclic operation↗
NASA NTRSDate not supplied

Improvements in Simulating a Mach 0.80 Transonic Truss-Braced Wing Configuration using the Spalart-Allmaras and k-ω SST Turbulence

Teams from the NASA Ames Research Centers (ARC) and Langley Research Center (LaRC) have been working on validating their computational fluid dynamics (CFD) results for the Boeing Mach 0.80 Transonic Truss-Braced Wing (TTBW) configuration. Experimental data used for the validation were gathered from a test conducted of a 4.5% scale Mach 0.80 TTBW model in the NASA Ames Research Center 11-by 11-Foot Transonic Wind Tunnel. The CFD simulations were initially run with both LAVA and USM3D Mixed Element solvers utilizing the Spalart-Allmaras (SA) turbulence model. A discrepancy was observed between CFD and experimental loads and moments ranging on average from 0.047 to 0.063 for lift coefficient, 16.5 to 27 in drag counts, and -0.02 to -0.032 for pitching moment, varying with the solver used. With introducing the refactored version of LAVA these values dropped to 0.034 to 0.047 for lift coefficient, 7.74 to 16.5 in drag counts, and lastly -0.012 to -0.02 for pitching moment. Based on findings from Boeing, who observed an improved comparison to experimental data when using the k-ω Shear Stress Transport (SST) turbulence model, the NASA teams conducted simulations with SST to investigate turbulence modeling sensitivities. OVERFLOW and USM3D V6 solvers were used for the comparison of SST and SA simulations. CFD results using the SST turbulence model demonstrate an improvement in matching with experimental CL values, reducing the discrepancies seen by 0.021-0.055 (62-93% reduction in discrepancy respective to the OVERFLOW and USM3D V6). The SST model has varying effects on CD, based on the solver / grid paradigm; for USM3D this prediction is only improved at higher angles of attack (above the cruise design point) resulting in a reduced discrepancy of 8.8–25 drag counts However, at the mid-to-low angles of attack, SST increases the CD discrepancy by 4-55 counts, increasing at lower angles of attack. For OVERFLOW this discrepancy only exists at the low angles of attack and maxes out at 10 counts of drag. Above α = 1.5° OVERFLOW SST reduces discrepancy by 9.3 drag counts on average.

AATT↗
NASA NTRS2011

Orion Handling Qualities During ISS Rendezvous and Docking

The Orion spacecraft was designed to rendezvous with multiple vehicles in low earth orbit (LEO) and beyond. To perform the required rendezvous and docking task, Orion must provide enough control authority to perform coarse translational maneuvers while maintaining precision to perform the delicate docking corrections. While Orion has autonomous docking capabilities, it is expected that final approach and docking operations with the International Space Station (ISS) will initially be performed in a manual mode. A series of evaluations was conducted by NASA and Lockheed Martin at the Johnson Space Center to determine the handling qualities (HQ) of the Orion spacecraft during different docking and rendezvous conditions using the Cooper-Harper scale. This paper will address the specifics of the handling qualities methodology, vehicle configuration, scenarios flown, data collection tools, and subject ratings and comments. The initial Orion HQ assessment examined Orion docking to the ISS. This scenario demonstrates the Translational Hand Controller (THC) handling qualities of Orion. During this initial assessment, two different scenarios were evaluated. The first was a nominal docking approach to a stable ISS, with Orion initializing with relative position dispersions and a closing rate of approximately 0.1 ft/sec. The second docking scenario was identical to the first, except the attitude motion of the ISS was modeled to simulate a stress case ( 1 degree deadband per axis and 0.01 deg/sec rate deadband per axis). For both scenarios, subjects started each run on final approach at a docking port-to-port range of 20 ft. Subjects used the THC in pulse mode with cues from the docking camera image, window views, and range and range rate data displayed on the Orion display units. As in the actual design, the attitude of the Orion vehicle was held by the automated flight control system at 0.5 degree deadband per axis. Several error sources were modeled including Reaction Control System (RCS) jet angular and position misalignment, RCS thrust magnitude uncertainty, RCS jet force direction uncertainty due to self plume impingement, and Orion center of mass uncertainty.

Hart, Jeremy J.↗
NASA NTRS2011

Failure Modes in Capacitors When Tested Under a Time-Varying Stress

Power-on failure has been the prevalent failure mechanism for solid tantalum capacitors in decoupling applications. A surge step stress test (SSST) has been previously applied to identify the critical stress level of a capacitor batch to give some predictability to the power-on failure mechanism [1]. But SSST can also be viewed as an electrically destructive test under a time-varying stress (voltage). It consists of rapidly charging the capacitor with incremental voltage increases, through a low resistance in series, until the capacitor under test is electrically shorted. When the reliability of capacitors is evaluated, a highly accelerated life test (HALT) is usually adopted since it is a time-efficient method of determining the failure mechanism; however, a destructive test under a time-varying stress such as SSST is even more time efficient. It usually takes days or weeks to complete a HALT test, but it only takes minutes for a time-varying stress test to produce failures. The advantage of incorporating a specific time-varying stress profile into a statistical model is significant in providing an alternative life test method for quickly revealing the failure mechanism in capacitors. In this paper, a time-varying stress that mimics a typical SSST has been incorporated into the Weibull model to characterize the failure mechanism in different types of capacitors. The SSST circuit and transient conditions for correctly surge testing capacitors are discussed. Finally, the SSST was applied for testing Ta capacitors, polymer aluminum capacitors (PA capacitors), and multi-layer ceramic (MLC) capacitors with both precious metal electrodes (PME) and base metal electrodes (BME). The test results are found to be directly associated with the dielectric layer breakdown in Ta and PA capacitors and are independent of the capacitor values, the way the capacitors were built, and the capacitors manufacturers. The test results also show that MLC capacitors exhibit surge breakdown voltages much higher than the rated voltage and that the breakdown field is inversely proportional to the dielectric layer thickness. The SSST data can also be used to comparatively evaluate the voltage robustness of capacitors for decoupling applications.

Liu, David (Donhang)↗