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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,765 records · Page 98

All sources NASA NTRS DOE OSTI

DOE OSTI2024

A macro-micro approach for identifying crystal plasticity parameters for necking and failure in nickel-based alloy haynes 282

Here, this work develops a two-scales macro-micro approach to address the challenge in calibrating crystal plasticity microstructural models when samples undergo necking prior to fracture. The crystal plasticity models are crucial for predicting the materials’ plastic deformation and failure at the microstructure level, identifying the materials’ intrinsic properties as well as investigating the microstructure-properties relationships. However, after necking occurs, the experimentally measured stress-strain curves fail to reflect the materials ‘true’ stress-strain behavior and cannot be directly fitted into crystal plasticity models. The proposed macro-micro approach employs a top-down strategy to address this challenge, which has been studied with experimental tests on precipitation-strengthened Ni-based superalloy Haynes® 282®. In this approach, a macro rate-dependent anisotropic plasticity model with Voce-type hardening and Rice-Tracey damage law is first utilized to model the deformation and failure of the tensile bar, and calibrated by matching the stress-strain curves, necking strain, and reduction of area. Especially, to match the testing results under different applied strain rates, the rate-sensitivity parameter m and saturation stress in the elasticity model are modified to incorporate dependence on the local strain rate. Then, the ‘true’ stress-strain behaviors are extracted from the necking zone of the macro-model, which are used to calibrate a micro-model with explicit microstructures and governed by an extended crystal plasticity law. The consistency between the micro-model and macro-model are enforced during calibration. The calibration outcomes from the crystal plasticity model elucidate the materials intrinsic properties for slip, hardening, and failure, which is vital for further investigations into the microstructure-properties relationship and for accurate prediction of the material behavior under various test and service conditions.

36 MATERIALS SCIENCE↗

NASA NTRS1987

Stress history of the Tharsis Region, Mars

The Tharsis topographic rise of Mars is roughly 5000 km wide and 10 km high and is believed to have originated more than 3.5 BY ago. Within its boundaries lie the four largest volcanoes on the planet. It is also the locus of a series of fracture traces which extend over approximately a hemisphere. The events leading to the formation of the Tharsis region continue to generate debate. Three geophysical models of the formation of Tharsis are now in general contention and each of these models has been used to predict a characteristic stress-field. These models are: the volcanic construct model, the isostatic compensation model, and the lithospheric uplift model. Each has been used by its proponents to predict some of the features observed in the Tharsis region but none accurately accounts for all of the fracture features observed. This is due, in part, to the use of fractures too young to be directly related to the origin of Tharsis. To constrain the origin of Tharsis, as opposed to its later history, one should look for the oldest fractures related to Tharsis and compare these to the predictions made by the models. Mapping of old terrains in and around the Tharsis rise has revealed 175 hitherto unknown old fracture features.

Francis, Robert A.↗

NASA NTRS2008

Ares I-X Upper Stage Simulator Structural Analyses Supporting the NESC Critical Initial Flaw Size Assessment

The structural analyses described in the present report were performed in support of the NASA Engineering and Safety Center (NESC) Critical Initial Flaw Size (CIFS) assessment for the ARES I-X Upper Stage Simulator (USS) common shell segment. The structural analysis effort for the NESC assessment had three thrusts: shell buckling analyses, detailed stress analyses of the single-bolt joint test; and stress analyses of two-segment 10 degree-wedge models for the peak axial tensile running load. Elasto-plastic, large-deformation simulations were performed. Stress analysis results indicated that the stress levels were well below the material yield stress for the bounding axial tensile design load. This report also summarizes the analyses and results from parametric studies on modeling the shell-to-gusset weld, flange-surface mismatch, bolt preload, and washer-bearing-surface modeling. These analyses models were used to generate the stress levels specified for the fatigue crack growth assessment using the design load with a factor of safety.

Knight, Norman F., Jr.↗

DOE OSTI2021

Dislocation drag and its influence on elastic precursor decay

Plastic deformation is mediated by the creation and movement of dislocations, and at high stress the latter is dominated by dislocation drag from phonon wind. Here, by simulating a 1-D shock impact problem we analyze the importance of accurately modeling dislocation drag and dislocation density evolution in the high stress regime. Dislocation drag is modeled according to a first-principles derivation as a function of stress and dislocation character, and its temperature and density dependence are approximated to the extent currently known. Much less is known about dislocation density evolution, leading to far greater uncertainty in these model parameters. In studying anisotropic fcc metals with character dependent dislocations, the present work generalizes similar earlier studies by other authors.

36 MATERIALS SCIENCE↗

NASA NTRS1993

Bounds on internal state variables in viscoplasticity

A typical viscoplastic model will introduce up to three types of internal state variables in order to properly describe transient material behavior; they are as follows: the back stress, the yield stress, and the drag strength. Different models employ different combinations of these internal variables--their selection and description of evolution being largely dependent on application and material selection. Under steady-state conditions, the internal variables cease to evolve and therefore become related to the external variables (stress and temperature) through simple functional relationships. A physically motivated hypothesis is presented that links the kinetic equation of viscoplasticity with that of creep under steady-state conditions. From this hypothesis one determines how the internal variables relate to one another at steady state, but most importantly, one obtains bounds on the magnitudes of stress and back stress, and on the yield stress and drag strength.

Freed, Alan D.↗

DOE OSTI2020

Lifetime estimation of a BWR core shroud in terms of IGSCC

The continued operation of aging Boiling Water Reactors (BWRs) worldwide requires gradually increasing cost of inspection, maintenance, and repair. Intergranular Stress Corrosion Cracking (IGSCC) in sensitized austenitic stainless steel piping first became a major issue for BWRs in the 1980s, resulting in recognition of the susceptibility of reactor internals to IGSCC. Shroud cracking identified in 1993–1994 confirmed that IGSCC of internals is a significant issue for BWRs. IGSCC is a time-dependent, material degradation process, which is caused and accelerated by the presence of residual stresses, material sensitization, irradiation, cold work, elevated temperature, and corrosive environments. This paper emphasizes the importance of accounting for corrosive environments, or more exactly, electrochemical phenomena in modelling IGSCC and predicting the service life of BWR in-vessel components, and stresses the necessity of performing such modeling not just for a single state point under full power conditions, but for the whole operating history of the reactor, including startups and shutdowns. Overall, this paper demonstrates that ignoring electrochemical considerations may result in underestimating component lifetime and lead to unnecessary expenses for inspection and repair.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Raman spectroscopic characterization of SiO 2 phase transformation and Si substrate stress relevant to EBC performance

To accurately model the long-term durability of environmental barrier coatings (EBCs), a more complete understanding of the phase composition and transformations of the thermally grown oxide SiO 2 (TGO) is desired. For the TGO formed during thermal cycling in steam, cristobalite formation and the subsequent β- to α-cristobalite transformation has been identified as a potentially life-limiting mechanism. In this study, Raman micro-spectroscopy was used to quantify the cristobalite transformation on a polycrystalline Si coupon that was exposed to steam at 1350°C for 100 h. The phase transformation was mapped at 200–260°C on the TGO surface at different ramp rates using a heating stage and a micro-positioning stage. The stress in the Si substrate was also determined using Raman spectroscopy by measuring the stress induced peak shift. The α→β phase transformation produced a 300–500 MPa tensile stress in the Si substrate, which compared well to the stress predicted from the volumetric expansion of the cristobalite. In conclusion, quantifying the phase transformation and residual stress are critical tools in developing the next generation of high performance EBCs.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Stress Intensity Thresholds for Development of Reliable Brittle Materials

Brittle material failure in high consequence systems can appear random and unpredictable at subcritical stresses. Gaps in our understanding of how structural flaws and environmental factors (humidity, temperature) impact fracture propagation need to be addressed to circumvent this issue. A combined experimental and computational approach composed of molecular dynamics (MD) simulations, numerical modeling, and atomic force microscopy (AFM) has been undertaken to identify mechanisms of slow crack growth in silicate glasses. AFM characterization of crack growth as slow as 10 -13 m/s was observed, with some stepwise crack growth. MD simulations have identified the critical role of inelastic relaxation in crack propagation, including evolution of the structure during relaxation. A numerical model for the existence of a stress intensity threshold, a stress intensity below which a fracture will not propagate, was developed. This transferrable model for predicting slow crack growth is being incorporated into mission-based programs.

36 MATERIALS SCIENCE↗

NASA NTRS1989

SSME LOX post flow analysis/fluid structure interaction. Volume 1: Flow analysis

The realization of measures to improve the performance of the Space Shuttle is, to a large extent, dependent on an improved understanding of the fluid flow phenomena occurring in the main engine. The overall arrangement of the primary components of the Space Shuttle Main Engine (SSME) are presented. The impingement of the hot gases from the transfer ducts onto the LOX posts causes them to be subjected to severe thermal and gas dynamic loads, which in the past have resulted in the occasional breakage of some elements of the outer row of posts during test firings of the engine, particularly at higher power levels. Large velocities in the gap between the LOX posts may also be a contributing factor in causing dynamic stability problems. The deforming structural response of the posts to the pressure loading may likely affect the gas flowfield by producing a moving flowfield boundary, thereby creating a dynamically coupled unsteady fluid-structure system. The objective was to investigate the three-dimensional, turbulent flow around a simplified SSME LOX post array using an existing Reynolds averaged Navier-Stokes flow solver and a suitable turbulence model to parameterize the turbulent shear stresses. Numerical computations were performed to analyze the effect on the flowfield of varying the spacing between the LOX posts, which were modeled as rigid, three-dimensional circular cylinders. The methodology used in the computations is described. Results are presented.

Burke, Roger W.↗

NASA NTRS1990

Modeling of near-wall turbulence

An improved k-epsilon model and a second order closure model is presented for low Reynolds number turbulence near a wall. For the k-epsilon model, a modified form of the eddy viscosity having correct asymptotic near wall behavior is suggested, and a model for the pressure diffusion term in the turbulent kinetic energy equation is proposed. For the second order closure model, the existing models are modified for the Reynolds stress equations to have proper near wall behavior. A dissipation rate equation for the turbulent kinetic energy is also reformulated. The proposed models satisfy realizability and will not produce unphysical behavior. Fully developed channel flows are used for model testing. The calculations are compared with direct numerical simulations. It is shown that the present models, both the k-epsilon model and the second order closure model, perform well in predicting the behavior of the near wall turbulence. Significant improvements over previous models are obtained.

Shih, T. H.↗

NASA NTRS1991

The dynamical evolution of the protosolar nebula

Evolutionary models for protostellar nebulae are calculated under the hypothesis that the only source for the turbulent viscosity is thermal convection. The viscous stress is approximated by an 'alpha' model, and the constant alpha is calculated in terms of the properties of turbulent thermal convection. A relatively sensitive dependence of the Rosseland mean opacity on temperature is needed for the vertical temperature gradient of the nebula to become convectively unstable. However, this requires that the vertical optical depth in the nebula is relatively sensitive to the total surface density, and as the disk is depleted of matter by protostellar accretion, the Rosseland mean optical depth drops and the disk must become convectively unstable. This limits the amount of mass that the nebula can accretionally process before the convective turbulence ceases. The resulting disk evolutionary properties are calculated and comparisons with the solar system are made.

Ruden, Steven P.↗

DOE OSTI2021

Finite Element Analysis (FEA) for Water-Foam Fracturing of Granite Rock

In addition to the foam data that were obtained from literature and that were collected from the current study, simulation data was also generated from finite element analysis (FEA) conducted in this study using COMSOL Multiphysics software. The FEA models were built to simulate the experiments conducted at Oak Ridge National Laboratory (ORNL) on cement and granite samples. In these FEA models, temperature was kept at ambient while the pressure profile resembled the loading conditions during the ORNL experiments, where pressure was either monotonically increased or applied cyclically. The cement material was used as a model material and was used to study Von Mises stress and tensile stress distribution for different bore hole length geometry using a parametric sweep with water as fracturing fluid using solid-fluid interaction module. For the granite material, FEA models were developed for stress analysis of cylindrical samples with water or foam fluids. The solid mechanics module in COMSOL was implemented to solve for Von Mises stress and tensile stress. The fluid-structure interaction module was implemented to solve for water-foam interaction on granite cylinder with addition of fluid-loading on structure, i.e., large deformation in solid mechanics with no impact on fluid deformation. Foam was considered as a pseudo single-phase compressible fluid for which material properties were calculated from water and gas (nitrogen) phases. The density of foam is calculated as a function of the densities of water and nitrogen, while viscosity is a function of temperature. Four types of FEA analyses were modelled: 1. Monotonic injection with water 2. Monotonic injection with foam 3. Cyclic injection with water 4. Cyclic injection with foam All the COMSOL files are converted to a zip file which is save in .mph.

15 GEOTHERMAL ENERGY↗

NASA NTRS1983

Large-area sheet task: Advanced dendritic-web-growth development

Thermally generated stresses in the growing web crystal were reduced. These stresses, which if too high cause the ribbon to degenerate, were reduced by a factor of three, resulting in the demonstrated growth of high-quality web crystals to widths of 5.4 cm. This progress was brought about chiefly by the application of thermal models to the development of low-stress growth configurations. A new temperature model was developed which can analyze the thermal effects of much more complex lid and top shield configurations than was possible with the old lumped shield model. Growth experiments which supplied input data such as actual shield temperature and melt levels were used to verify the modeling results. Desirable modifications in the melt level-sensing circuitry were made in the new experimental web growth furnace, and this furnace has been used to carry out growth experiments under steady-state conditions. New growth configurations were tested in long growth runs at Westinghouse AESD which produced wider, lower stress and higher quality web crystals than designs previously used.

Duncan, C. S.↗

DOE OSTI2024

Assessment of Machine Learning Wall Modeling Approaches for Large Eddy Simulation of Gas Turbine Film Cooling Flows: An a Priori Study

Here, in this work, a priori analysis of machine learning (ML) strategies is carried out with the goal of data-driven wall modeling for large eddy simulation (LES) of gas turbine film cooling flows. High-fidelity flow datasets are extracted from wall-resolved LES (WRLES) of flow over a flat plate interacting with the coolant flow supplied by a single row of 7-7-7 shaped cooling holes inclined at 30 degrees with the flat plate at different blowing ratios (BR). The WRLES are performed using the high-order Nek5000 spectral element computational fluid dynamics (CFD) solver. Light gradient boosting machine (LightGBM) is employed as the ML algorithm for the data-driven wall model. Parametric tests are conducted to systematically assess the influence of a wide range of input flow features (velocity components, velocity gradients, pressure gradients, and fluid properties) on the accuracy of ML wall model with respect to prediction of wall shear stress. In addition, the use of spatial stencil and time delay is also explored within the ML wall modeling framework. It is shown that features associated with gradients of the streamwise and spanwise velocity components have a major impact on the prediction fidelity of wall model, while the effect of gradients of wall-normal velocity component is found to be negligible. Moreover, adding flow feature information from an x-y-z spatial stencil significantly improves the ML model accuracy and generalizability compared to just using local flow features from the matching location. Overall, highest prediction accuracy is achieved when both spatial stencil and time delay features are incorporated within the data-driven wall modeling paradigm.

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS1992

A creep cavity growth model for creep-fatigue life prediction of a unidirectional W/Cu composite

A microstructural model was developed to predict creep-fatigue life in a (0)(sub 4), 9 volume percent tungsten fiber-reinforced copper matrix composite at the temperature of 833 K. The mechanism of failure of the composite is assumed to be governed by the growth of quasi-equilibrium cavities in the copper matrix of the composite, based on the microscopically observed failure mechanisms. The methodology uses a cavity growth model developed for prediction of creep fracture. Instantaneous values of strain rate and stress in the copper matrix during fatigue cycles were calculated and incorporated in the model to predict cyclic life. The stress in the copper matrix was determined by use of a simple two-bar model for the fiber and matrix during cyclic loading. The model successfully predicted the composite creep-fatigue life under tension-tension cyclic loading through the use of this instantaneous matrix stress level. Inclusion of additional mechanisms such as cavity nucleation, grain boundary sliding, and the effect of fibers on matrix-stress level would result in more generalized predictions of creep-fatigue life.

Kim, Young-Suk↗

DOE OSTI2020

Phase field modeling of discontinuous dynamic recrystallization in hot deformation of magnesium alloys

We have developed a quantitative phase field model of discontinuous dynamic recrystallization (PF-DDRX) for investigating grain structure evolution and its impact on mechanical response in magnesium alloy during hot deformation. A set of isothermal compression tests were conducted by Gleeble thermo-mechanical simulator and used to determine PF-DDRX model parameters such as critical stresses of recrystallization, activation energies of dislocation, grain boundaries and recrystallization grain nucleation. In the simulations, the experimentally measured average grain size and distribution were used to generate initial microstructures. The results show that the typical ‘necklace’ microstructures observed in DDRX forms at the earlier stage of recrystallization. The inhomogeneity of microstructures first increases rapidly and then decreases gradually and reaches a constant with the increase of dynamic recrystallization volume fraction. Both predicted grain structure and mechanical response are in good agreement with experimental results. This demonstrates that the model parameters determined by experimental data are reliable for the developed PF-DDRX model to predict the grain structure evolution and flow behavior of AZ80 magnesium alloy under hot deformation. This work presents a path to combine experiments and modeling to give quantitative prediction and minimize experimental cost.

magnesium alloy, Phase field method, dynamic recry↗

NASA NTRS2006

FORTRAN Versions of Reformulated HFGMC Codes

Several FORTRAN codes have been written to implement the reformulated version of the high-fidelity generalized method of cells (HFGMC). Various aspects of the HFGMC and its predecessors were described in several prior NASA Tech Briefs articles, the most recent being HFGMC Enhancement of MAC/GMC (LEW-17818-1), NASA Tech Briefs, Vol. 30, No. 3 (March 2006), page 34. The HFGMC is a mathematical model of micromechanics for simulating stress and strain responses of fiber/matrix and other composite materials. The HFGMC overcomes a major limitation of a prior version of the GMC by accounting for coupling of shear and normal stresses and thereby affords greater accuracy, albeit at a large computational cost. In the reformulation of the HFGMC, the issue of computational efficiency was addressed: as a result, codes that implement the reformulated HFGMC complete their calculations about 10 times as fast as do those that implement the HFGMC. The present FORTRAN implementations of the reformulated HFGMC were written to satisfy a need for compatibility with other FORTRAN programs used to analyze structures and composite materials. The FORTRAN implementations also afford capabilities, beyond those of the basic HFGMC, for modeling inelasticity, fiber/matrix debonding, and coupled thermal, mechanical, piezo, and electromagnetic effects.

Arnold, Steven M.↗

NASA NTRS1976

Mathematical modelling of undrained clay behavior

The proposed general analytical model describes the anisotropic, elastoplastic, path-dependent, stress-strain properties of inviscid saturated clays under undrained conditions. Model parameters are determined by using results from strain-controlled simple shear tests on a saturated clay. The model's accuracy is evaluated by applying it to predict the results of other tests on the same clay, including monotonic and cyclic loading. The model explains the very anisotropic shear strength behavior observed for weak marine clays.

Prevost, J. H.↗