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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Strength and deformation properties of basaltic lava flows on planetary surfaces

Basaltic rocks are thought to constitute a volumetrically significant rock type on the Moon, Mercury, Mars, and Venus, in addition to the Earth. Spacecraft images of surfaces with known or suspected basaltic composition on these bodies, particularly on Venus, indicate that these rocks have been deformed in the brittle regime to form faults and perhaps dilatant cracks, in addition to folding and more distributed types of deformation. Predictions of brittle fracture or other types of deformation are made by comparing calculated stresses from a tectonic model to some criterion for rock strength. Common strength criteria used in the planetary science literature for near-surface deformation include a Griffith tensile-strength criterion for intact rock, a Mohr envelope for intact basalt, and a brittle strength envelope based on Byerlee's law of rock frictional resistance. However, planetary terrains of basaltic composition consist of much more than just intact basaltic rock. The aggregate basaltic material, termed the 'rock mass,' consists of both the intact rock and the associated fracture, faults, lithologic contacts, and other discontinuous surfaces. A basaltic rock mass is the relevant material for which strength properties must be defined and calculated model stresses must be compared to in order to more accurately predict brittle deformation. For example, the various strengths of a rock mass are less than that of intact material of the same composition. This means that tectonic models which compare stresses to intact failure strengths overestimate the stresses required for fracture and so underestimate the extent and magnitude of brittle deformation predicted in these models. On the other hand, rock mass shear strength can be greater than that predicted from Byerlee's law. The concept of rock mass strength is central to many engineering design studies in which calculated stresses are used to predict brittle fracture, and this experience indicates that brittle strength envelopes which assume properties for intact rock (Griffith parabolas) or sliding along a single, continuous surface (Byerlee's law) inadequately characterize the tensile, compressive, and shear strengths of rock masses. The criterion adopted here to relate stresses to rock mass fracture is based on a Griffith-type curve for tensile normal stress and a concave downward curve for compressive normal stress. It is the only available criterion that explicitly considers the weakening effects of discontinuities within the rock mass on the stress state required for fracture.

Schultz, Richard A.↗
NASA NTRS2005

Representation of Vegetation and Other Nonerodible Elements in Aeolian Shear Stress Partitioning Models for Predicting Transport Threshold

The presence of nonerodible elements is well understood to be a reducing factor for soil erosion by wind, but the limits of its protection of the surface and erosion threshold prediction are complicated by the varying geometry, spatial organization, and density of the elements. The predictive capabilities of the most recent models for estimating wind driven particle fluxes are reduced because of the poor representation of the effectiveness of vegetation to reduce wind erosion. Two approaches have been taken to account for roughness effects on sediment transport thresholds. Marticorena and Bergametti (1995) in their dust emission model parameterize the effect of roughness on threshold with the assumption that there is a relationship between roughness density and the aerodynamic roughness length of a surface. Raupach et al. (1993) offer a different approach based on physical modeling of wake development behind individual roughness elements and the partition of the surface stress and the total stress over a roughened surface. A comparison between the models shows the partitioning approach to be a good framework to explain the effect of roughness on entrainment of sediment by wind. Both models provided very good agreement for wind tunnel experiments using solid objects on a nonerodible surface. However, the Marticorena and Bergametti (1995) approach displays a scaling dependency when the difference between the roughness length of the surface and the overall roughness length is too great, while the Raupach et al. (1993) model's predictions perform better owing to the incorporation of the roughness geometry and the alterations to the flow they can cause.

King, James↗
NASA NTRS1993

Critical assessment of Reynolds stress turbulence models using homogeneous flows

In modeling the rapid part of the pressure correlation term in the Reynolds stress transport equations, extensive use has been made of its exact properties. These, for example, have been employed in obtaining the widely used LRR model. Some recent proposals have dropped one of these properties to obtain new models. We demonstrate, by computing some simple homogeneous flows, that doing so does not lead to any significant improvements over the LRR model and it is not the right direction in improving the performance of existing models.

Shabbir, Aamir↗
NASA NTRS2001

Isolating Curvature Effects in Computing Wall-Bounded Turbulent Flows

The flow over the zero-pressure-gradient So-Mellor convex curved wall is simulated using the Navier-Stokes equations. An inviscid effective outer wall shape, undocumented in the experiment, is obtained by using an adjoint optimization method with the desired pressure distribution on the inner wall as the cost function. Using this wall shape with a Navier-Stokes method, the abilities of various turbulence models to simulate the effects of curvature without the complicating factor of streamwise pressure gradient can be evaluated. The one-equation Spalart-Allmaras turbulence model overpredicts eddy viscosity, and its boundary layer profiles are too full. A curvature-corrected version of this model improves results, which are sensitive to the choice of a particular constant. An explicit algebraic stress model does a reasonable job predicting this flow field. However, results can be slightly improved by modifying the assumption on anisotropy equilibrium in the model's derivation. The resulting curvature-corrected explicit algebraic stress model possesses no heuristic functions or additional constants. It lowers slightly the computed skin friction coefficient and the turbulent stress levels for this case (in better agreement with experiment), but the effect on computed velocity profiles is very small.

Rumsey, Christopher L.↗
NASA NTRS1989

A numerical and experimental study of coaxial jets

An algebraic stress model and the standard k-epsilon model is applied to predict the mean and turbulence quantities for axisymmetric, nonswirling coaxial jets without confinement. To investigate the effects of numerical (false) diffusion on the predicted results, three different discretization schemes, namely, hybrid, power-law, and the flux-spline, are employed. In addition, an experimental study is conducted to provide data of good quality, especially near the inlet, for model assessment. The results show that the use of the algebraic stress model leads to better agreement between the numerical results and experimental data.

Nikjooy, M.↗
NASA NTRS1992

Critical assessment of Reynolds stress turbulence models using homogeneous flows

In modeling the rapid part of the pressure correlation term in the Reynolds stress transport equations, extensive use has been made of its exact properties which were first suggested by Rotta. These, for example, have been employed in obtaining the widely used Launder, Reece and Rodi (LRR) model. Some recent proposals have dropped one of these properties to obtain new models. We demonstrate, by computing some simple homogeneous flows, that doing so does not lead to any significant improvements over the LRR model and it is not the right direction in improving the performance of existing models. The reason for this, in our opinion, is that violation of one of the exact properties can not bring in any new physics into the model. We compute thirteen homogeneous flows using LRR (with a recalibrated rapid term constant), IP and SSG models. The flows computed include the flow through axisymmetric contraction; axisymmetric expansion; distortion by plane strain; and homogeneous shear flows with and without rotation. Results show that for most general representation for a model linear in the anisotropic tensor, performs either better or as good as the other two models of the same level.

Shabbir, Aamir↗
NASA NTRS1991

Creep and stress relaxation modeling of polycrystalline ceramic fibers

A variety of high performance polycrystalline ceramic fibers are currently being considered as reinforcement for high temperature ceramic matrix composites. However, under mechanical loading above 800 C, these fibers display creep-related instabilities which can result in detrimental changes in composite dimensions, strength, and internal stress distributions. As a first step toward understanding these effects, this study examines the validity of mechanistic-based empirical model which describes primary stage tensile creep and stress relaxation of polycrystalline ceramic fibers as independent functions of time, temperature, and applied stress or strain. To verify these functional dependencies, a simple bend test is used to measure stress relaxation for four types of commercial ceramic fibers for which direct tensile creep data are available. These fibers include both nonoxide (SCS-6, Nicalon) and oxide (PRD-166, FP) compositions. The results of the bend stress relaxation (BSR) test not only confirm the stress, time, and temperature dependencies predicted by the model but also allow measurement of model empirical parameters for the four fiber types. In addition, comparison of model predictions and BSR test results with the literature tensile creep data show good agreement, supporting both the predictive capability of the model and the use of the BSR test as a simple method for parameter determination for other fibers.

Dicarlo, James A.↗
NASA NTRS1994

Creep and stress relaxation modeling of polycrystalline ceramic fibers

A variety of high performance polycrystalline ceramic fibers are currently being considered as reinforcement for high temperature ceramic matrix composites. However, under mechanical loading about 800 C, these fibers display creep related instabilities which can result in detrimental changes in composite dimensions, strength, and internal stress distributions. As a first step toward understanding these effects, this study examines the validity of a mechanism-based empirical model which describes primary stage tensile creep and stress relaxation of polycrystalline ceramic fibers as independent functions of time, temperature, and applied stress or strain. To verify these functional dependencies, a simple bend test is used to measure stress relaxation for four types of commercial ceramic fibers for which direct tensile creep data are available. These fibers include both nonoxide (SCS-6, Nicalon) and oxide (PRD-166, FP) compositions. The results of the Bend Stress Relaxation (BSR) test not only confirm the stress, time, and temperature dependencies predicted by the model, but also allow measurement of model empirical parameters for the four fiber types. In addition, comparison of model tensile creep predictions based on the BSR test results with the literature data show good agreement, supporting both the predictive capability of the model and the use of the BSR text as a simple method for parameter determination for other fibers.

Dicarlo, James A.↗
NASA NTRS1993

A comprehensive comparison of turbulence models in the far wake

In the present study, the far wake was examined numerically using an implicit, upwind, finite-volume, compressible Navier-Stokes code. The numerical grid started at 500 equivalent circular cylinder diameters in the wave, and extended to 4000 equivalent diameters. By concentrating only on the far wake, the numerical difficulties and fine mesh requirements near the wake-generating body were eliminated. At the time of this writing, results for the K-epsilon and K-omega turbulence models at low Mach number have been completed and show excellent agreement with previous incompressible results and far-wake similarity solutions. The code is presently being used to compare the performance of various other turbulence models, including Reynolds stress models and the new anisotropic two-equation turbulence models being developed at NASA Langley. By increasing our physical understanding of the deficiencies and limits of these models, it is hoped that improvements to the universality of the models can be made. Future plans include examination of two-dimensional momentumless wakes as well.

Cimbala, John M.↗
NASA NTRS1995

Turbulence program for propulsion systems

Program goals at the Center for Modeling of Turbulence and Transition (CMOTT), NASA Lewis Research Center, are (1) to develop reliable turbulence (including bypass transition) and combustion models for complex flows in propulsion systems and (2) to integrate developed models into deliverable CFD tools for propulsion systems in collaboration with industry. This viewgraph presentation covers the following topics: development of turbulence and combustion models; collaboration with industry and technology transfer; isotropic eddy viscosity models; algebraic Reynolds stress models; scalar turbulence models; second order closure models; multiple scale k-epsilon models; and PDF modeling of turbulent reacting flows.

Shih, Tsan-Hsing↗
NASA NTRS1991

Software Models Impact Stresses

Generalized Impact Stress Software designed to assist engineers in predicting stresses caused by variety of impacts. Program straightforward, simple to implement on personal computers, "user friendly", and handles variety of boundary conditions applied to struck body being analyzed. Applications include mathematical modeling of motions and transient stresses of spacecraft, analysis of slamming of piston, of fast valve shutoffs, and play of rotating bearing assembly. Provides fast and inexpensive analytical tool for analysis of stresses and reduces dependency on expensive impact tests. Written in FORTRAN 77. Requires use of commercial software package PLOT88.

Hanshaw, Timothy C.↗
NASA NTRS2006

Modeling the Stress Strain Behavior of Woven Ceramic Matrix Composites

Woven SiC fiber reinforced SiC matrix composites represent one of the most mature composite systems to date. Future components fabricated out of these woven ceramic matrix composites are expected to vary in shape, curvature, architecture, and thickness. The design of future components using woven ceramic matrix composites necessitates a modeling approach that can account for these variations which are physically controlled by local constituent contents and architecture. Research over the years supported primarily by NASA Glenn Research Center has led to the development of simple mechanistic-based models that can describe the entire stress-strain curve for composite systems fabricated with chemical vapor infiltrated matrices and melt-infiltrated matrices for a wide range of constituent content and architecture. Several examples will be presented that demonstrate the approach to modeling which incorporates a thorough understanding of the stress-dependent matrix cracking properties of the composite system.

Morscher, Gregory N.↗
NASA NTRS2018

Plasticity and Damage Modeling of Stress Asymmetry and Dynamic Behavior of AFS Additive Manufactured Aluminum Alloy 2219

The Solid State Additive Manufacturing (AM) process referred as MELD that fabricated the samples in this study, provides a new path for repairing, coating, joining and additive manufacturing metals and metal matrix composites. This research will be the first application of a physics-based microstructure dependent internal state variable (ISV) plasticity and damage material model to capture the mechanical response of an AM Aluminum Alloy (AA) 2219 via the MELD process. In this research, a microstructure-based internal state variable (ISV) plasticity-damage model was used to capture the mechanical behavior of AFS 2219 aluminum alloy. Aeroprobe Corporation, creator and patent holder for the MELD process, fabricated the material by pushing a solid filler rod of AA2219-T861 material through a hollow rotating tool onto an AA2219 T851 plate substrate. As feedstock, solid or powder precursor metals are pushed through a nonconsumable rotating cylindrical tool. Herein, added layers are deposited and metallurgically bonded to substrate material or previously deposited layers by the heat generated from the rotating tool through plastic deformation of the filler material. Once a layer has been added, the tool height increases, and starts the deposition of the next layer. This process results in beneficial properties such as grain refinement, homogenization and reduced porosity (fully dense). This process will experience temperatures similar to those in the weld nugget zone (WNZ) in friction stir welding (FSW), ranging from 0.6-0.9 Tm, with Tm being the melting point of the material. MELD is highly scalable with AA deposition rates reaching over 1000 cm3/hr, which allows for MELD being used for repairs, coatings, and building components. A motivating factor driving the research for physics-based history dependent material modeling of MELD components is the ability to accurately capture the stress-state and strain rate dependence in the material caused by variations in material microstructure from the MELD processing of new or repaired components. The ISV model incorporates microstructural content and is consistent with continuum level kinematics, kinetics, and thermodynamics. These features allow the ISV model to capture large deformations at the structural scale using the kinematic and isotropic hardening, while microscale damage is obtained from the microstructural features. The benefits of the ISV model arise from the inclusion of structure-property relationships identified from microstructural characterization and experimentation. The Bauschinger effect (BE) is an important concept, vital in the accurate prediction of cyclic stress-strain response of ductile materials such as metals. The ISV model has been successfully used to capture the behavior and damage, and the BE of different aluminum alloys and steels. The ISV model uses kinematic and isotropic hardening to help capture deformations of the material at the macro scale. To understand this hardening relationship, calculating the kinematic and isotropic hardening relationship in the material is warranted for a high-fidelity model. Electron Backscattered Diffraction (EBSD) was used to characterize the as-fabricated microstructure, where a fully-dense equiaxed grain morphology with average grain size of 2.5 μm was observed. Microhardness mapping of the as-built structures, monotonic tension and compression experiments at both quasi-static (0.001/s) strain rates, tension-followed-by-compression and compression-followed-by-tension experiments were performed to obtain the set of plasticity and damage constants necessary to capture strain rate and stress state behavior of this additive material. To calibrate the plasticity-damage model, a single set of constants were determined to capture the different stress states the MELD AA2219. One set of the constants was determined from experimental true stress-strain curves for the tension and compression data. Additionally, microstructural information and data from the open literature were used as the other model constants. This research is a first of its kind for AFS AA2219, includes correlating the ISV model to the monotonic experimental results that capture the isotropic and kinematic plasticity mechanical response.

Rivera, O. G.↗
NASA NTRS1994

Computation and modeling of aero-thermal fields in turbine cascades and strongly curved ducts

Advanced turbulence models are crucial for accurate prediction of rocket engine flows, due to existence of very large extra strain rates, such as strong streamline curvature. Numerical simulation of the turbulent flow in a strongly curved turn-around duct (TAD) has been carried out with a Reynolds stress model (RSM), an algebraic Reynolds stress model (ARSM) and a kappa-epsilon model. the RSM model and the ARSM model are found to capture the turbulence damping due to the convex curvature, but underpredict the turbulence enhancement caused by the concave curvature. To capture the concave curvature effects, it is necessary to modify the epsilon-equation. The modification of episilon-equation suggested by Launder, et.al, provides the correct trend, but over-corrects the curvature effects. A comparative study of two modes of transition in gas turbine, the by-pass transition and the separation-induced transition, has been carried out with several low-Reynolds-number (LRN) kappa-epsilon models. Effects of blade surface pressure gradient, freestream turbulence and Reynolds number on the blade boundary layer development, and particularly the inception of transition are examined in detail. The present study indicates that the turbine blade transition, in the presence of high freestream turbulence, is predicted well with LRN kappa-epsilon models employed.

Luo, J.↗
NASA NTRS1990

Large-eddy simulation of transitional channel flow

A large-eddy simulation (LES) of transition in plane channel flow was carried out. The LES results were compared with those of a fine direct numerical simulation (DNS), and with those of a coarse DNS that uses the same mesh as the LES, but does not use a residual stress model. While at the early stages of transition, LES and coarse DNS give the same results: the presence of the residual stress model was found to be necessary to predict accurately mean velocity and Reynolds stress profiles during the late stages of transition (after the second spike stage). The evolution of single Fourier modes is also predicted more accurately by the LES than by the DNS. As small scales are generated, the dissipative character of the residual stress starts to reproduce correctly the energy cascade. As transition progresses, the flow approaches its fully developed turbulent state, the subgrid scales tend towards equilibrium, and the model becomes more accurate.

Piomelli, Ugo↗
NASA NTRS2017

Shallow Lunar Seismic Activity and the Current Stress State of the Moon

A vast, global network of more than 3200 lobate thrust fault scarps has been revealed in high resolution Lunar Reconnaissance Orbiter Camera (LROC) images. The fault scarps are very young, less than 50 Ma, based on their small scale and crisp appearance, crosscutting relations with small-diameter impact craters, and rates of infilling of associated small, shallow graben and may be actively forming today. The population of young thrust fault scarps provides a window into the recent stress state of the Moon and offers insight into the origin of global lunar stresses. The distribution of orientations of the fault scarps is non-random, inconsistent with isotropic stresses from late-stage global contraction as the sole source of stress. Modeling shows that tidal stresses contribute significantly to the current stress state of the lunar crust. Tidal stresses (orbital recession and diurnal tides) superimposed on stresses from global contraction result in non-isotropic compressional stress and may produce thrust faults consistent with lobate scarp orientations. At any particular point on the lunar surface, peak compressive stress will be reached at a certain time in the diurnal cycle. Coseismic slip events on currently active thrust faults are expected to be triggered when peak stresses are reached. Analysis of the timing of the 28 the shallow moonquakes recorded by the Apollo seismic network shows that 19 indeed occur when the Moon is closer to apogee, while only 9 shallow events occur when the Moon is closer to perigee. Here we report efforts to refine the model for the current stress state of the Moon by investigating the contribution of polar wander. Progress on relocating the epicentral locations of the shallow moonquakes using an algorithm designed for sparse networks is also reported.

structure↗
NASA NTRS1990

Near-wall modelling of compressible turbulent flows

Work was carried out to formulate near-wall models for the equations governing the transport of the temperature-variance and its dissipation rate. With these equations properly modeled, a foundation is laid for their extension together with the heat-flux equations to compressible flows. This extension is carried out in a manner similar to that used to extend the incompressible near-wall Reynolds-stress models to compressible flows. The methodology used to accomplish the extension of the near-wall Reynolds-stress models is examined and the actual extension of the models for the Reynolds-stress equations and the near-wall dissipation-rate equation to compressible flows is given. Then the formulation of the near-wall models for the equations governing the transport of the temperature variance and its dissipation rate is discussed. Finally, a sample calculation of a flat plate compressible turbulent boundary-layer flow with adiabatic wall boundary condition and a free-stream Mach number of 2.5 using a two-equation near-wall closure is presented. The results show that the near-wall two-equation closure formulated for compressible flows is quite valid and the calculated properties are in good agreement with measurements. Furthermore, the near-wall behavior of the turbulence statistics and structure parameters is consistent with that found in incompressible flows.

So, Ronald M. C.↗
NASA NTRS1985

Burner liner thermal/structural load modeling: TRANCITS program user's manual

Transfer Analysis Code to Interface Thermal/Structural Problems (TRANCITS) is discussed. The TRANCITS code satisfies all the objectives for transferring thermal data between heat transfer and structural models of combustor liners and it can be used as a generic thermal translator between heat transfer and stress models of any component, regardless of the geometry. The TRANCITS can accurately and efficiently convert the temperature distributions predicted by the heat transfer programs to those required by the stress codes. It can be used for both linear and nonlinear structural codes and can produce nodal temperatures, elemental centroid temperatures, or elemental Gauss point temperatures. The thermal output of both the MARC and SINDA heat transfer codes can be interfaced directly with TRANCITS, and it will automatically produce stress model codes formatted for NASTRAN and MARC. Any thermal program and structural program can be interfaced by using the neutral input and output forms supported by TRANCITS.

Maffeo, R.↗