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,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Program Model Checking as a New Trend

This paper introduces a special section of STTT (International Journal on Software Tools for Technology Transfer) containing a selection of papers that were presented at the 7th International SPIN workshop, Stanford, August 30 - September 1, 2000. The workshop was named SPIN Model Checking and Software Verification, with an emphasis on model checking of programs. The paper outlines the motivation for stressing software verification, rather than only design and model verification, by presenting the work done in the Automated Software Engineering group at NASA Ames Research Center within the last 5 years. This includes work in software model checking, testing like technologies and static analysis.

Havelund, Klaus↗
NASA NTRS1998

Thermal Effects Modeling Developed for Smart Structures

Applying smart materials in aeropropulsion systems may improve the performance of aircraft engines through a variety of vibration, noise, and shape-control applications. To facilitate the experimental characterization of these smart structures, researchers have been focusing on developing analytical models to account for the coupled mechanical, electrical, and thermal response of these materials. One focus of current research efforts has been directed toward incorporating a comprehensive thermal analysis modeling capability. Typically, temperature affects the behavior of smart materials by three distinct mechanisms: Induction of thermal strains because of coefficient of thermal expansion mismatch 1. Pyroelectric effects on the piezoelectric elements; 2. Temperature-dependent changes in material properties; and 3. Previous analytical models only investigated the first two thermal effects mechanisms. However, since the material properties of piezoelectric materials generally vary greatly with temperature (see the graph), incorporating temperature-dependent material properties will significantly affect the structural deflections, sensory voltages, and stresses. Thus, the current analytical model captures thermal effects arising from all three mechanisms through thermopiezoelectric constitutive equations. These constitutive equations were incorporated into a layerwise laminate theory with the inherent capability to model both the active and sensory response of smart structures in thermal environments. Corresponding finite element equations were formulated and implemented for both the beam and plate elements to provide a comprehensive thermal effects modeling capability.

Lee, Ho-Jun↗
DOE OSTI2022

Film Cooling Performance Prediction for Air and Supercritical CO 2

The current study explores the possibility of cooling the vanes and blades of a direct-fired sCO 2 turbine using film cooling. The operating conditions of a direct-fired sCO 2 cycle and thermophysical properties of the fluid at those conditions can alter the flow field characteristics of the coolant jet and its mixing with the mainstream. Very little information is present in the literature regarding the performance of film cooling geometries employing supercritical CO 2 . The objective of this study is to estimate the resulting film cooling effectiveness while also capturing the effects of the crossflow-to-mainstream velocity ratio on the coolant jet. Herein, a computational fluid dynamic model is used to study the coolant jet exiting a cylindrical hole located on a flat plate, with the coolant fed by an internal channel. Steady-state Reynolds-averaged Navier–Stokes equations were solved along with the (shear-stress transport) SST k–ω model to provide the turbulence closure. The operating conditions for the direct-fired sCO 2 turbine are obtained using an in-house Cooled Turbine Model. Numerical predictions revealed that the crossflow effects and jet lift-off were more pronounced in the case of sCO 2 when compared to air. Spatial distribution of flow field and cooling effectiveness are presented at different operating conditions.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Burst Pressure Solutions of Thin and Thick-Walled Cylindrical Vessels

Pressure vessels (PVs) are widely used in the energy industry. Accurate burst pressure is critical to structural design and safe operation for both thin and thick-walled PVs. The traditional strength theories utilized a single-parameter material property, such as the yield stress or the ultimate tensile stress (UTS) to develop failure models for determining the yield or ultimate pressure carrying capacity in the PV design. The UTS-based Barlow formula is a typical burst pressure model developed from the Tresca strength theory that provides the basis for developing regulation rules and failure models for different industry design codes, such as ASME BPVC, ASME B31.3, and ASME B31G, among others. In order to reduce the conservatism of the Tresca strength model, ASME BPVC recently adapted failure models developed from the von Mises strength theory for the PV design and analysis. It has been commonly accepted that the burst pressure of pipelines depends on the UTS and strain hardening exponent, n, of the pipeline steel. An average shear stress yield theory was thus developed, and the Zhu-Leis solution of burst pressure was obtained as a function of UTS and n for thin-walled line pipes. Experiments showed that the Zhu-Leis solution provides an accurate, reliable prediction of burst pressure for defect-free thin-walled pipes. In order to extend the Zhu–Leis solution to thick-walled cylindrical PVs, this paper defined three new flow stresses, modified the traditional strength theories, and obtained three new burst pressure solutions that are valid for both thin and thick-walled cylindrical vessels. The proposed flow stresses are able to describe the tensile strength and the plastic flow response of PVs for a strain hardening steel. The associated strength theories were then developed in terms of the Tresca, von Mises, and Zhu-Leis yield criteria. From these new strength theories, three burst pressure solutions were obtained for thick-walled cylinders, where the von Mises solution is an upper bound prediction, the Tresca solution is a lower bound prediction, and the Zhu-Leis solution is an intermediate prediction of burst pressure for thick-walled cylinders. Lastly, the proposed burst pressure solutions were evaluated and validated by two large datasets of full-scale burst tests for thick-walled tubes and for thin-walled pipes.

42 ENGINEERING↗
NASA NTRSDate not supplied

Effect of Boundary Conditions on Process-Induced Stresses in a Plain Weave Unit Cell

Woven polymer matrix composites (PMCs) are leveraged in aerospace applications for their desirable specific properties, yet they are vulnerable to high residual stresses during manufacturing and their complex geometry makes experimental results difficult to observe. Process modeling is needed to characterize the effects of the curing and predict end stress states. Finite element software can be used to model woven architectures, however accurate representation of processing conditions remains a challenge when it comes to selecting boundary conditions. The effect of BCs on process-induced stress within woven PMCs is studied. The commercial Finite Element Analysis (FEA) software Abaqus is coupled with user-written subroutines in a process modeling framework. A two-dimensionally (2D) woven PMC repeating unit cell (RUC) is modeled with TexGen and Abaqus. Virtual curing is imposed on the bulk matrix. The BC study is conducted with Free, Periodic, Flat, and Flat-Free configurations. Results show that the end stress state is sensitive to the boundary condition assumptions. Flat BC results show great agreement with Periodic BCs. Residual stress results from process modeling are then compared with a linear-elastic thermal cooldown analysis in Abaqus. Cooldown results indicate an overestimation in matrix stresses compared with process modeling.

micromechanics↗
NASA NTRS1992

Towards a structure-based turbulence model

The author is developing a new type of turbulence model in which a new one-point quantity, the eddy structure tensor, carries information about the two-point structure of the turbulence. The model was motivated by the observation that conventional one-point turbulence models based only on the turbulent stresses do not predict the rapid changes in state that are found when anisotropic homogeneous turbulence is subjected to mean rotation, and hence are fundamentally incorrect for rotation. The model appears to give topologically correct predictions for the changes in stress state and structure state under all types of rapid distortions of homogeneous turbulence.

Reynolds, W. C.↗
DOE OSTI2022

Model calibration of the liquid mercury spallation target using evolutionary neural networks and sparse polynomial expansions

The mercury constitutive model predicting the strain and stress in the target vessel plays a central role in improving the lifetime prediction and future target designs of the mercury targets at the Spallation Neutron Source. We leverage the experiment strain data collected over multiple years to improve the mercury constitutive model through a combination of large scale simulations of the target behavior and the use of machine learning tools for parameter estimation. We present two interdisciplinary approaches for surrogate-based model calibration of expensive simulations using evolutionary neural networks and sparse polynomial expansions. The newly calibrated simulations achieve 7% average improvement on the prediction accuracy and 8% reduction in mean absolute error compared to previously reported reference parameters, with some individual sensors experiencing up to 30% improvement. The calibrated simulations can aid in fatigue analysis to estimate the mercury target lifetime, which reduces abrupt failure and saves tremendous amount of costs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2004

Thermostructural Analysis of Carbon Cloth Phenolics "Ply Lifting" and Correlation to LHMEL Test Results

This paper provides a discussion of the history of Carbon Cloth Phenolic (CCP) ply lifting in the Redesigned Solid Rocket Motor (RSRM) Program, a brief presentation of theoretical methods used for analytical evaluation, and results of parametric analyses of CCP material subject to test conditions of the Laser Hardened Material Evaluation Laboratory. CCP ply lift can occur in regions of the RSRM nozzle where ply angle to flame surface is generally less than about 20 degrees. There is a heat rate dependence on likelihood and severity of the condition with the higher heating rates generally producing more ply lift. The event occurs in-depth, near the heated surface, where the load necessary to mechanically separate the CCP plies is produced by the initial stages of pyrolysis gas generation due to the thermal decomposition of the phenolic resin matrix. Due to the shallow lay-up angle of the composite, normal components of the indepth mechanical load, due to "pore pressure", are imparted primarily as a cross-ply tensile force on the interlaminar ply boundaries. Tensile capability in the cross-ply (out of plane) direction is solely determined by the matrix material capability. The elevated temperature matrix material capabilities are overcome by pressure induced mechanical normal stress and ply-lift occurs. A theoretical model used for CCP in-depth temperature, pressure, and normal stress prediction, based on first principles, is briefly discussed followed by a parametric evaluation of response variables subject to boundary conditions typical of on-going test programs at the LHMEL facility. Model response demonstrates general trends observed in test and provides insight into the interactivity of material properties and constitutive relationships.

Clayton, Louie↗
NASA NTRS2001

Intermediate Temperature Strength Degradation in SiC/SiC Composites

Woven silicon carbide fiber-reinforced, silicon carbide matrix composites are leading candidate materials for an advanced jet engine combustor liner application. Although the use temperature in the hot region for this application is expected to exceed 1200 C, a potential life-limiting concern for this composite system exists at intermediate temperatures (800 +/- 200 C), where significant time-dependent strength degradation has been observed under stress-rupture loading. A number of factors control the degree of stress-rupture strength degradation, the major factor being the nature of the interphase separating the fiber and the matrix. BN interphases are superior to carbon interphases due to the slower oxidation kinetics of BN. A model for the intermediate temperature stress-rupture of SiC/BN/SiC composites is presented based on the observed mechanistic process that leads to strength degradation for the simple case of through-thickness matrix cracks. The approach taken has much in common with that used by Curtin and coworkers, for two different composite systems. The predictions of the model are in good agreement with the rupture data for stress-rupture of both precracked and as-produced composites. Also, three approaches that dramatically improve the intermediate temperature stress-rupture properties are described: Si-doped BN, fiber spreading, and 'outside debonding'.

Morscher, Gregory N.↗
DOE OSTI2020

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗
DOE OSTI2025

Developing a robust strength model using physically-informed genetic programming

The strength of materials is influenced by a range of external conditions, such as temperature and deformation rate. Consequently, materials that demonstrate substantial variations in their mechanical behavior due to fluctuations in temperature and strain rate require complex strength models to accurately predict material performance in real-world applications. To predict such complex behavior, a robust and flexible strength model is necessary. In this work, we utilize genetic programming-based symbolic regression (GPSR) to develop data-driven strength models that accurately represent the measured stress–strain responses of tin across a wide range of strain, strain rate and temperature regimes. The GPSR models are constrained by physically-informed conditions, which leads to significant improvement in extrapolation. The best model is integrated into a multi-physics code to perform Taylor impact simulations, validating the model’s accuracy and robustness. In conclusion, the model predictions showed excellent agreement with experimental results, particularly when compared to predictions using traditional strength models.

Genetic programming↗
NASA NTRS1984

An elastic potential for the nonlinear response of unidirectional graphite composites

An elastic potential is proposed that is capable of modeling the reversible portion of the observed nonlinear response of unidirectional graphite fiber composites. The model includes both the stiffening stress-strain behavior as well as the softening Poisson's response for loading in the fiber direction. The model is compared with experimental results for Celion 6000/PMR-15 graphite-polyimide.

Pindera, M.-J.↗
NASA NTRS1989

Treatment of singularities in cracked bodies

Three-dimensional finite-element analyses of middle-crack tension (M-T) and bend specimens subjected to mode I loadings were performed to study the stress singularity along the crack front. The specimen was modeled using 20-node isoparametric elements. The displacements and stresses from the analysis were used to estimate the power of singularities using a log-log regression analysis along the crack front. The analyses showed that finite-sized cracked bodies have two singular stress fields of the form rho = C sub o (theta, z) r to the -1/2 power + D sub o (theta, phi) R to the lambda rho power. The first term is the cylindrical singularity with the power -1/2 and is dominant over the middle 96 pct (for Poisson's ratio = 0.3) of the crack front and becomes nearly zero at the free surface. The second singularity is a vertex singularity with the vertex point located at the intersection of the crack front and the free surface. The second term is dominant at the free surface and becomes nearly zero away from the the boundary layer. The thickness of the boundary layer depends on Poisson's ratio of the material and is independent of the specimen type. The thickness of the boundary layer varied from 0 pct to about 5 pct of the total specimen thickness as Poisson's ratio varied from 0.0 to 0.45. Because there are two singular stress fields near the free surface, the strain energy release rate (G) is an appropriate parameter to measure the severity of the crack.

Shivakumar, K. N.↗
NASA NTRS1990

Treatment of singularities in cracked bodies

Three-dimensional finite-element analyses of middle-crack tension (M-T) and bend specimens subjected to mode I loadings were performed to study the stress singularity along the crack front. The specimen was modeled using 20-node isoparametric elements. The displacements and stresses from the analysis were used to estimate the power of singularities using a log-log regression analysis along the crack front. The analyses showed that finite-sized cracked bodies have two singular stress fields of the form rho = C sub o (theta, z) r to the -1/2 power + D sub o (theta, phi) R to the lambda rho power. The first term is the cylindrical singularity with the power -1/2 and is dominant over the middle 96 pct (for Poisson's ratio = 0.3) of the crack front and becomes nearly zero at the free surface. The second singularity is a vertex singularity with the vertex point located at the intersection of the crack front and the free surface. The second term is dominant at the free surface and becomes nearly zero away from the boundary layer. The thickness of the boundary layer depends on Poisson's ratio of the material and is independent of the specimen type. The thickness of the boundary layer varied from 0 pct to about 5 pct of the total specimen thickness as Poisson's ratio varied from 0.0 to 0.45. Because there are two singular stress fields near the free surface, the strain energy release rate (G) is an appropriate parameter to measure the severity of the crack.

Shivakumar, K. N.↗
NASA NTRS1992

Free-edge stress analysis of glass-epoxy laminates with matrix cracks

The effect of matrix cracks on the composite delamination and interlaminar stresses is investigated in (+15/90n/-15)s glass-epoxy laminates (with values of n equal to 0, 1, 2, or 3) subjected to monotonically increasing tension loads. Three-dimensional (3D) and quasi-3D (Q3D) finite-element analyses are used to model the free-edge stress states in the laminates with and without a matrix crack, respectively. The Q3D results show that in-plane transverse tensile stresses exist in the +15 deg plies near the free edges of all of the laminates used and that only the interlaminar shear stress is high at the +15/theta interface. The results of 3D analysis indicate that large tensile interlaminar normal as well as shear stresses develop at the intersection of the matrix crack and the free edge. This suggests that the interlaminar normal stress plays a significant role in the failure of these laminates.

Fish, John C.↗
NASA NTRS1995

Calculations of Diffuser Flows with an Anisotropic K-Epsilon Model

A newly developed anisotropic K-epsilon model is applied to calculate three axisymmetric diffuser flows with or without separation. The new model uses a quadratic stress-strain relation and satisfies the realizability conditions, i.e., it ensures both the positivity of the turbulent normal stresses and the Schwarz' inequality between any fluctuating velocities. Calculations are carried out with a finite-element method. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to ensure numerical credibility of the solutions. The standard K-epsilon model is also used in order to highlight the performance of the new model. Comparison with the experimental data shows that the anisotropic K-epsilon model performs consistently better than does the standard K-epsilon model in all of the three test cases.

Zhu, J.↗
NASA NTRS2015

Numerical Calculations of 3-D High-Lift Flows and Comparison with Experiment

Solutions were obtained with the Navier-Stokes CFD code TLNS3D to predict the flow about the NASA Trapezoidal Wing, a high-lift wing composed of three elements: the main-wing element, a deployed leading-edge slat, and a deployed trailing-edge flap. Turbulence was modeled by the Spalart-Allmaras one-equation turbulence model. One case with massive separation was repeated using Menter's two-equation SST (Menter's Shear Stress Transport) k-omega turbulence model in an attempt to improve the agreement with experiment. The investigation was conducted at a free stream Mach number of 0.2, and at angles of attack ranging from 10.004 degrees to 34.858 degrees. The Reynolds number based on the mean aerodynamic chord of the wing was 4.3 x 10 (sup 6). Compared to experiment, the numerical procedure predicted the surface pressures very well at angles of attack in the linear range of the lift. However, computed maximum lift was 5% low. Drag was mainly under predicted. The procedure correctly predicted several well-known trends and features of high-lift flows, such as off-body separation. The two turbulence models yielded significantly different solutions for the repeated case.

Compton, William B, III↗
NASA NTRS2020

Tensile Creep Behavior of SiC(f)/SiC Ceramic Matrix Minicomposites

Single fiber-tow minicomposites represent the major load-bearing element of woven and laminate ceramic matrix composites (CMCs). To understand the effects of fiber type, fiber content, and matrix cracking on tensile creep in SiCf/SiC CMCs, single-tow SiCf/SiC minicomposites with different fiber types and contents were investigated. The minicomposites studied contained either Hi-Nicalon™ or Hi-Nicalon™ Type S SiC fibers with a boron nitride (BN) interphase and a chemical-vapor-infiltrated-silicon-carbide (CVI-SiC) matrix. Tensile creep was performed at 1200 °C in air. A bottom-up creep modeling approach was applied where creep parameters of the fibers and matrix were obtained separately at 1200 °C. Next, a theoretical model based on the rule of mixtures was derived to model the fiber and matrix creep-time-dependent stress redistribution. Fiber and matrix creep parameters, load transfer model results, and numerical modeling were used to construct a creep strain model to predict creep damage evolution of minicomposites with different fiber types and contents.

SiC/SiC Minicomposites Creep↗