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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Measurement of Creep Properties of Ultra-High-Temperature Materials by a Novel Non-Contact Technique

A non-contact technique for measuring the creep properties of materials has been developed and validated as part of a collaboration among the University of Massachusetts, NASA Marshall Space Flight Center Electrostatic Levitation Facility (ESL), and the University of Tennessee. This novel method has several advantages over conventional creep testing. The sample is deformed by the centripetal acceleration from the rapid rotation, and the deformed shapes are analyzed to determine the strain. Since there is no contact with grips, there is no theoretical maximum temperature and no concern about chemical compatibility. Materials may be tested at the service temperature even for extreme environments such as rocket nozzles, or above the service temperature for accelerated testing of materials for applications such as jet engines or turbopumps for liquid-fueled engines. The creep measurements have been demonstrated to 2400 C with niobium, while the test facility, the NASA MSFC ESL, has processed materials up to 3400 C. Furthermore, the ESL creep method employs a distribution of stress to determine the stress exponent from a single test, versus the many tests required by conventional methods. Determination of the stress exponent from the ESL creep tests requires very precise measurement of the surface shape of the deformed sample for comparison to deformations predicted by finite element models for different stress exponents. An error analysis shows that the stress exponent can be determined to about 1% accuracy with the current methods and apparatus. The creep properties of single-crystal niobium at 1985 C showed excellent agreement with conventional tests performed according to ASTM Standard E-139. Tests on other metals, ceramics, and composites relevant to rocket propulsion and turbine engines are underway.

Hyers, Robert W.↗
DOE OSTI2025

Improved pressure-gradient sensor for the prediction of separation onset in RANS models

Here, we improve upon two key aspects of the Menter shear stress transport (SST) turbulence model: (1) We propose a more robust adverse pressure gradient sensor based on the strength of the pressure gradient in the direction of the local mean flow; (2) We propose two alternative eddy viscosity models to be used in the adverse pressure gradient regions identified by our sensor. Direct numerical simulations of the Boeing Gaussian bump are used to identify the terms in the baseline SST model that need correction, and a posteriori Reynolds-averaged Navier-Stokes calculations are used to calibrate coefficient values, leading to a model that is both physics driven and data informed. The two sensor-equipped models are applied to two thick airfoils representative of modern wind turbine applications, the FFA-W3-301 and the DU00-W-212. The proposed models improve the prediction of stall (onset of separation) with respect to the prediction of the baseline SST model.

17 WIND ENERGY↗
DOE OSTI2022

Microstructure-Sensitive Uncertainty Quantification for Crystal Plasticity Finite Element Constitutive Models Using Stochastic Collocation Methods

Uncertainty quantification (UQ) plays a major role in verification and validation for computational engineering models and simulations, and establishes trust in the predictive capability of computational models. In the materials science and engineering context, where the process-structure-property-performance linkage is well known to be the only road mapping from manufacturing to engineering performance, numerous integrated computational materials engineering (ICME) models have been developed across a wide spectrum of length-scales and time-scales to relieve the burden of resource-intensive experiments. Within the structure-property linkage, crystal plasticity finite element method (CPFEM) models have been widely used since they are one of a few ICME toolboxes that allows numerical predictions, providing the bridge from microstructure to materials properties and performances. Several constitutive models have been proposed in the last few decades to capture the mechanics and plasticity behavior of materials. While some UQ studies have been performed, the robustness and uncertainty of these constitutive models have not been rigorously established. In this work, we apply a stochastic collocation (SC) method, which is mathematically rigorous and has been widely used in the field of UQ, to quantify the uncertainty of three most commonly used constitutive models in CPFEM, namely phenomenological models (with and without twinning), and dislocation-density-based constitutive models, for three different types of crystal structures, namely face-centered cubic (fcc) copper (Cu), body-centered cubic (bcc) tungsten (W), and hexagonal close packing (hcp) magnesium (Mg). Our numerical results not only quantify the uncertainty of these constitutive models in stress-strain curve, but also analyze the global sensitivity of the underlying constitutive parameters with respect to the initial yield behavior, which may be helpful for robust constitutive model calibration works in the future.

36 MATERIALS SCIENCE↗
NASA NTRS1995

Three-dimensional CTOA and constraint effects during stable tearing in a thin-sheet material

A small strain theory, three-dimensional elastic-plastic finite element analysis was used to simulate fracture in thin sheet 2024-T3 aluminum alloy in the T-L orientation. Both straight and tunneled cracks were modeled. The tunneled crack front shapes as a function of applied stress were obtained from the fracture surface of tested specimens. The stable crack growth behavior was measured at the specimen surface as a function of applied stress. The fracture simulation modeled the crack tunneling and extension as a function of applied stress. The results indicated that the global constraint factor, alpha(sub g), initially dropped during stable crack growth. After peak applied stress was achieved, alpha(sub g) began to increase slightly. The effect of crack front shape on alpha(sub g) was small, but the crack front shape did greatly influence the local constraint and through-thickness crack-tip opening angle (CTOA) behavior. The surface values of CTOA for the tunneled crack front model agreed well with experimental measurements, showing the same initial decrease from high values during the initial 3mm of crack growth at the specimen's surface. At the same time, the interior CTOA values increased from low angles. After the initial stable tearing region, the CTOA was constant through the thickness. The three-dimensional analysis appears to confirm the potential of CTOA as a two-dimensional fracture criterion.

Dawicke, D. S.↗
NASA NTRS1988

A closure model of diffusion transport of the Reynolds-stress equations and its application to a turbulent step flow

A study is made for the development and computations of the separating and reattaching shear flows. The highlight of the study is an attempt to predict the third-moments of turbulent velocity which is responsible for the diffusion transport of the Reynolds stresses. The present computations show that the third-moments obtained by employing the low-Reynolds number model of transport equations improve the prediction of the third-moments.

Amano, R. S.↗
NASA NTRS2020

Soil Evaporation Stress Determines Soil Moisture‐Evapotranspiration Coupling Strength in Land Surface Modeling

Model‐based estimates of soil moisture (SM)‐evapotranspiration (ET) coupling strength(ρ) vary widely and are prone to bias. Here we apply numerical modeling and remote sensing to identify the process‐level source of modeledρbias with the goal of improving the fidelity of current Earth system models. Results illustrate that modeledρis most strongly determined by soil evaporation (E) stress, and (generally positive)ρmodeling bias is attributable to the oversimplification of soil texture impacts on E stress. Based on new remotely sensed estimates ofρ, we demonstrate that removingρbias via a single optimized E stress parameter leads to improved ET accuracy and resolves a well‐known modeling bias in the partitioning of ET into E and T. As such, we highlight the importance of the stress function relating E and SM and its central role in regulating land‐atmosphere coupling processes impacting local climate.

Jianzhi Dong↗
DOE OSTI2025

Computational materials reliability assessment of hydrogen fueled gas turbine power generation engines

The use of blended fuel sources in land based gas turbine engines drives variations in the resulting operational profile (temperatures and pressures) which can impact engine reliability. Furthermore, variability in the manufacture of components affects the resulting microstructure which directly impacts material performance and reliability. Currently, data-driven models are typically used for maintaining and inspecting fleets of engines. Without explicitly capturing material and operational sources of variability conservatism must be used in developing component-level reliability models. Therefore, there exists an opportunity to use information from materials-scale physics models to better inform reliability modeling and reduce conservatism; the impact is more cost-efficient operation and maintenance of current and future fleets. Specifically, this work establishes a computational framework for evaluating the probabilistic high temperature creep performance of hot-section Ni-based superalloys where uncertainty comes from both microstructural and operational variability. A novel high-fidelity physics model which phenomenologically captures grain-boundary sensitive phenomena has been established. A probabilistic calibration procedure was used to calibrate the model and capture uncertainty in the parameterized model coefficients. A design of experiments methodology was established for identifying informative microstructural digital representations for suitable for forward model evaluation. Results show that training a machine-learning surrogate using this design criteria outperforms random selection of microstructural representations. Finally, two surrogate models were developed: (1) a deterministic surrogate model which predicts the local field response given microstructure, constitutive model parameters, and operating conditions (stress, temperature) and (2) a probabilistic model, where uncertainty comes from constitutive law uncertainty, built using denoising diffusion probabilistic models which samples responses given (1) microstructure and (2) operating conditions. These surrogate models enable partner Siemens Energy to rapidly perform UQ analysis specific to creep deformation across a range of microstructures and operating conditions. The impact is that these ML and physics codes can be used to establish more advanced reliability models for the inspection, servicing, and maintenance of land based gas turbine engines.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Estimating the daily course of Konza Prairie latent heat fluxes using a modified Tergra model

Experimental tests of the Tergra-2 model are based on data collected under moderately wet to wet and very dry soil moisture conditions. Further testing of the model under intermediate soil moisture conditions is required and additional testing under very dry conditions may lead to modifications that make the model more suitable to water-stressed conditions. Combining the Tergra model with a soil evaporation routine should enhance the accuracy of the model and allow it to be employed in situations where vapor fluxes are not almost solely attributable to transpiration.

Hope, Allen S.↗
NASA NTRS2003

Damage Accumulation in SiC/SiC Composites with 3D Architectures

The formation and propagation of multiple matrix cracks in relatively dense ceramic matrix composites when subjected to increasing tensile stress is necessary for high strength and tough composites. However, the occurrence of matrix cracks at low stresses may limit the usefulness of some non-oxide composite systems when subjected to oxidizing environments for long times at stresses sufficient to cause matrix cracking. For SiC fiber-reinforced composites with two-dimensional woven architectures and chemically vapor infiltrated (CVI) SiC matrix and melt-infiltrated (MI) Si/SiC matrix composites, the matrix cracking behavior has been fairly well characterized for different fiber-types and woven architectures. It was found that the occurrence, degree, and growth of matrix cracks depends on the material properties of the composite constituents as well as other physical properties of the composite or architecture, e.g., matrix porosity and size of the fiber bundle. In this study, matrix cracking in SiC fiber reinforced, melt-infiltrated SiC composites with a 3D orthogonal architecture was determined for specimens tested in tension at room temperature. Acoustic emission (AE) was used to monitor the matrix cracking activity, which was later confirmed by microscopic examination of specimens that had failed. The determination of the exact location of AE demonstrated that initial cracking occurred in the matrix rich regions when a large z-direction fiber bundle was used. For specimens with large z-direction fiber tows, the earliest matrix cracking could occur at half the stress for standard 2D woven composites with similar constituents. Damage accumulation in 3D architecture composites will be compared to damage accumulation in 2D architecture composites and discussed with respect to modeling composite stress-strain behavior and use of these composites at elevated temperatures.

Morscher, Gregory N.↗
NASA NTRS2008

Fatigue Behavior of a Third Generation PM Disk Superalloy

The fatigue behavior of a 3rd generation PM disk alloy, LSHR, was studied at 1300 F. Tensile, creep, and fatigue tests were run on smooth and notched (Kt = 2) bars under a variety of conditions. Analysis of smooth bar fatigue data, run under strain and load control with R ratios of 0 and -1, showed that a stress based Smith-Watson-Topper approach could collapse the data set. While the tensile and creep data showed substantial notch strengthening at 1300 F, the fatigue data showed a life deficit for the notch specimens. A viscoplastic finite element model, which accounted for stress relaxation at the notch tip, provided the best correlation between the notched and smooth bar behavior, although the fatigue data was not fully rationalized based on this simplified viscoplastic model of the stresses at the notch tip.Inclusion of a 90 sec dwell at peak load was found to dramatically decrease notch fatigue life. This result was shown to be consistent with a simple linear creep-fatigue damage rule, where creep damage dominated at low stresses and fatigue damage was more prevalent at higher stresses.

Gayda, John↗
NASA NTRS1991

North American plate dynamics

Deformation within the North American plate in response to various tectonic processes is modeled using an elastic finite element analysis. The tectonic processes considered in the modeling include ridge forces associated with the normal thermal evolution of oceanic lithosphere, shear and normal stresses transmitted across transforms, normal stresses transmitted across convergent boundaries, stresses due to horizontal density contrasts within the continent, and shear tractions applied along the base of the plate. Model stresses are calculated with respect to a lithostatic reference stress state. Shear stresses transmitted across transform boundaries along the San Andreas and Caribbean are small, of the order of 5-10 MPa. Also, compressive stresses of the order of 5-10 MPa transmitted across the major transforms improve the fit to the data. Compressive stresses across convergent margins along the Aleutians and the Middle America trench are important.

Richardson, Randall M.↗
DOE OSTI2026

On the roles of welding residual stresses in determination of fracture toughness in austenitic stainless steel SUS 304 pipeline girth welds

Welding residual stresses especially the high tensile stresses are proved to have negative impacts on the fatigue and fracture behaviors of welded structures. In this study, a virtual fabrication of test specimens from welding process to specimen preparation was carried out by numerical simulation. An austenitic stainless steel multi-pass pipe welding was simulated by transient thermal–mechanical finite element analysis, the residual stresses were then mapped into the test specimen to evaluate fracture toughness. The findings in this study confirmed that, residual stress can be high in a sub-sized compact tensile specimen, which may accelerate or hinder the crack propagation during actual fatigue and fracture tests as reported in recent years. The influence of the cutting location and orientation of the specimen on fracture performance was investigated systematically to provide a fundamental understanding of welding residual stress and necessary insights into the specimen preparation procedure. Considering the limitation of measuring techniques and the complexity of the stress distribution, the developed numerical model can be a very useful tool to elucidate the stress evolution and quantify the effect of remaining welding stress on fracture toughness.

Fracture behavior↗
DOE OSTI2023

3-D Simulations of earthquakes rupture jumps: 1. Homogeneous pre-stress conditions

SUMMARY Observational and modelling studies indicate that earthquake ruptures can jump between fault sections as large as ∼3 and ∼5 km for compressional and extensional offsets, respectively. Here, we compare characteristics of the rupture jump process on parallel but offset fault sections from traditional 3-D dynamic rupture simulations governed by slip weakening friction using the finite element code, FaultMod, to those from quasi-dynamic simulations governed by rate- and state-dependent friction (rate-state friction) using the code RSQSim. These simulations use spatially uniform initial stresses. For a variety of measures the rupture renucleation position on the offset fault, the rate-state friction and slip weakening friction models produce very similar results. The principal difference is the additional occurrence of delayed rupture jumps that arise from the time- and stress-dependent nucleation that is characteristic of rate-state friction. For immediate rupture jumps, models with slip weakening friction span greater offsets than those with rate-state friction. However the jump distances are nearly identical when delayed rupture jumps are included in the comparisons. We propose that delayed rupture jumps are the likely mechanism for adjacent large-earthquake pairs and clusters. Based on the similarity of renucleation positions with both dynamic and quasi-dynamic models, we conclude that the renucleation positions for rupture initiation on the receiver fault (separated by less than ∼3 km from the source fault) are primarily controlled by static stress changes induced by slip on the initiating fault. However, in light of the slightly greater maximum jump distances (>3 km) seen with the dynamic slip weakening friction model, dynamic stress changes from seismic waves play an increasingly important role as offset distances increase.

, RSQSim↗
NASA NTRS1995

Dynamics of the Antarctic Circumpolar Current. Evidence for Topographic Effects from Altimeter Data and Numerical Model Output

Geosat altimeter data and numerical model output are used to examine the circulation and dynamics of the Antarctic Circumpolar Current (ACC). The mean sea surface height across the ACC has been reconstructed from height variability measured by the altimeter, without assuming prior knowledge of the geoid. The results indicate locations for the Subantarctic and Polar Fronts which are consistent with in situ observations and indicate that the fronts are substantially steered by bathymetry. Detailed examination of spatial and temporal variability indicates a spatial decorrelation scale of 85 km and a temporal e-folding scale of 34 days. Empirical Orthogonal Function analysis suggests that the scales of motion are relatively short, occuring on 1000 km length-scales rather than basin or global scales. The momentum balance of the ACC has been investigated using output from the high resolution primitive equation model in combination with altimeter data. In the Semtner-Chervin quarter-degree general circulation model topographic form stress is the dominant process balancing the surface wind forcing. In stream coordinates, the dominant effect transporting momentum across the ACC is bibarmonic friction. Potential vorticity is considered on Montgomery streamlines in the model output and along surface streamlines in model and altimeter data. (AN)

Gille, Sarah T.↗
NASA NTRS1991

Design Strategy for a Formally Verified Reliable Computing Platform

This paper presents a high-level design for a reliable computing platform for real-time control applications. The design tradeoffs and analyses related to the development of a formally verified reliable computing platform are discussed. The design strategy advocated in this paper requires the use of techniques that can be completely characterized mathematically as opposed to more powerful or more flexible algorithms whose performance properties can only be analyzed by simulation and testing. The need for accurate reliability models that can be related to the behavior models is also stressed. Tradeoffs between reliability and voting complexity are explored. In particular, the transient recovery properties of the system are found to be fundamental to both the reliability analysis as well as the "correctness" models.

Butler, Ricky W.↗
DOE OSTI2020

P-model v1.0: an optimality-based light use efficiency model for simulating ecosystem gross primary production

Terrestrial photosynthesis is the basis for vegetation growth and drives the land carbon cycle. Accurately simulating gross primary production (GPP, ecosystem-level apparent photosynthesis) is key for satellite monitoring and Earth system model predictions under climate change. While robust models exist for describing leaf-level photosynthesis, predictions diverge due to uncertain photosynthetic traits and parameters which vary on multiple spatial and temporal scales. Here, we describe and evaluate a GPP (photosynthesis per unit ground area) model, the P-model, that combines the Farquhar–von Caemmerer–Berry model for C 3 photosynthesis with an optimality principle for the carbon assimilation–transpiration trade-off, and predicts a multi-day average light use efficiency (LUE) for any climate and C 3 vegetation type. The model builds on the theory developed in Prentice et al. (2014) and Wang et al. (2017a) and is extended to include low temperature effects on the intrinsic quantum yield and an empirical soil moisture stress factor. The model is forced with site-level data of the fraction of absorbed photosynthetically active radiation (fAPAR) and meteorological data and is evaluated against GPP estimates from a globally distributed network of ecosystem flux measurements. Although the P-model requires relatively few inputs, the R 2 for predicted versus observed GPP based on the full model setup is 0.75 (8 d mean, 126 sites) – similar to comparable satellite-data-driven GPP models but without predefined vegetation-type-specific parameters. The R 2 is reduced to 0.70 when not accounting for the reduction in quantum yield at low temperatures and effects of low soil moisture on LUE. The R 2 for the P-model-predicted LUE is 0.32 (means by site) and 0.48 (means by vegetation type). Applying this model for global-scale simulations yields a total global GPP of 106–122 Pg C yr –1 (mean of 2001–2011), depending on the fAPAR forcing data. The P-model provides a simple but powerful method for predicting – rather than prescribing – light use efficiency and simulating terrestrial photosynthesis across a wide range of conditions. The model is available as an R package (rpmodel).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1985

Features of a reattaching turbulent shear layer in divergent channel flow

Experimental data have been obtained in an incompressible turbulent flow over a rearward-facing step in a diverging channel flow. Mean velocities, Reynolds stresses, and triple products that were measured by a laser Doppler velocimeter are presented for two cases of tunnel wall divergence. Eddy viscosities, production, convection, turbulent diffusion, and dissipation (balance of kinetic energy equation) terms are extracted from the data. These data are compared with various eddy-viscosity turbulence models. Numerical calculations incorporating the k-epsilon and algebraic-stress turbulence models are compared with the data. When determining quantities of engineering interest, the modified algebraic-stress model (ASM) is a significant improvement over the unmodified ASM and the unmodified k-epsilon model; however, like the others, it dramatically overpredicts the experimentally determined dissipation rate.

Driver, D. M.↗
NASA NTRS1974

Compressive strength of fiber reinforced composite materials

Results of an experimental and analytical investigation of the compressive strength of unidirectional boron-epoxy composite material are presented. Observation of fiber coordinates in a boron-epoxy composite indicates that the fibers contain initial curvature. Combined axial compression and torsion tests were conducted on boron-epoxy tubes, and it was shown that the shear modulus is a function of axial compressive stress. An analytical model which includes initial curvature in the fibers and permits an estimate of the effect of curvature on compressive strength is proposed. Two modes of failure which may result from the application of axial compressive stress are analyzed, delamination and shear instability. Based on tests and analysis, failure of boron-epoxy under axial compressive load is due to shear instability.

Davis, J. G., Jr.↗