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,387 records · Page 77

Results of the third Marine Ice Sheet Model Intercomparison Project (MISMIP+)

We present the result of the third Marine Ice Sheet Model Intercomparison Project, MISMIP+. MISMIP+ is intended to be a benchmark for ice-flow models which include fast sliding marine ice streams and floating ice shelves and in particular a treatment of viscous stress that is sufficient to model buttressing, where upstream ice flow is restrained by a downstream ice shelf. A set of idealized experiments first tests that models are able to maintain a steady state with the grounding line located on a retrograde slope due to buttressing and then explore scenarios where a reduction in that buttressing causes ice stream acceleration, thinning, and grounding line retreat. The majority of participating models passed the first test and then produced similar responses to the loss of buttressing. We find that the most important distinction between models in this particular type of simulation is in the treatment of sliding at the bed, with other distinctions – notably the difference between the simpler and more complete treatments of englacial stress but also the differences between numerical methods – taking a secondary role.

58 GEOSCIENCES↗

The prediction of two- and three-dimensional asymmetrical turbulent wakes - A comparision of the performance of three turbulence models for the effects of streamline curvature and rotation

Two- and three-dimensional turbulent wakes which develop under the influence of streamline curvature and rotation were calculated with three turbulence closure models. The first model comprised of transport equations for the turbulent kinetic energy and the rate of energy dissipation. The second model comprised of equations for the rate of turbulent kinetic energy dissipation and Reynolds stresses, but the effects of the convection and diffusion in the Reynolds stress transport equation were handled collectively. The third model utilizes equations of turbulent kinetic energy dissipation and Reynolds stresses in nearly exact form. All of the three models were modified for the effects of streamline curvature. The second and third models include the effects of rotation through the rotation-originated redistribution term in the transport equation of Reynolds stresses. The numerical results for the wakes of airfoils, cascades and turbomachinery rotor blades demonstrate that the second and third models provide accurate predictions, but computer time and storage can be considerably saved with the second model.

Hah, C.↗

Exploring the Langtry-Menter Transition Model for High Speed Applications Using FUN3D

A series of Reynolds-averaged Navier-Stokes (RANS) simulations were performed using the FUN3D flow solver to explore the capabilities of the Langtry-Menter Shear-Stress Transport (LM-SST) transition model for predicting transition for aircraft inlet applications. Two geometries were simulated: a zeropressure-gradient flat plate and an axisymmetric cone exposed to hypersonic flow. In addition to the transition-sensitized LM-SST model investigations, simulations were run with the one-equation SpalartAllmaras (SA) and the two-equation Menter Shear-Stress Transport (SST-V) RANS models in fully turbulent mode to identify the natural RANS model transition behavior as a function of Mach number when executed in fully turbulent mode. The flat plate simulations showed that (1) the transition model was able to predict rapid transition at a freestream Mach number of 0.2, which is expected but (2) the predicted transition location moved downstream as the freestream Mach number was increased for the simulations that used the SST-V turbulence model. The latter is significant as it is usually assumed that one- and two-equation turbulence models will produce fully turbulent flow very near the boundary layer origin. The flat plate simulation freestream Mach number trend was confirmed with simulations using the Wind-US code, which also saw a similar trend when employing the SA turbulence model. For the axisymmetric cone simulations, the transition location was highly sensitive to the inflow turbulence levels. This is significant as the prediction of the transition location is crucial when trying to predict inlet performance, especially for hypersonic vehicle applications. It was also noted that the predicted transition location for the cone when using the SST-V turbulence model agreed well with the predicted transition location from the equivalent zero-pressure-gradient cold wall flat plate case.

Transition Model↗

A Non-Invasive Approach for Elucidating the Spatial Distribution of in-situ Stress in Deep Subsurface Geologic Formations Considered for CO 2 Storage (Task 4 (Field Testing) Report)

"This document presents a summary of TASK 4 (Field Testing) of the project entitled “A Non-Invasive Approach for Elucidating the Spatial Distribution of in-situ Stress in Deep Subsurface Geologic Formations Considered for CO 2 Storage (FE0031686)”. This three-year project is part of the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL) Carbon Storage program to research and address gaps that affect the economics of commercial CCS projects. One of these key gaps is the lack of certainty in predicting the geomechanical impacts of pressure migration due to CO 2 injection into a storage complex. This report documents field testing that was conducted under TASK 4 to obtain input data needed to develop the site-specific geomechanical model in TASK 5 and to verify the model results (i.e. calculated stresses) for this site. The field work entailed collecting geophysical logs and core samples and conducting geomechanical stress tests in the Core Energy LLC State Otsego Lake (SOL 8-15A) well, located in Otsego County, Michigan."

58 GEOSCIENCES↗

Stress-constrained topology optimization of lattice-like structures using component-wise reduced order models

We report lattice-like structures can provide a combination of high stiffness with light weight that is useful in many applications, but a resolved finite element mesh of such structures results in a computationally expensive discretization. This computational expense may be particularly burdensome in many-query applications, such as optimization. We develop a stress-constrained topology optimization method for lattice-like structures that uses component-wise reduced order models as a cheap surrogate, providing accurate computation of stress fields while greatly reducing run time relative to a full order model. We demonstrate the ability of our method to produce large reductions in mass while respecting a constraint on the maximum stress in a pair of test problems. The ROM methodology provides a speedup of about 150x in forward solves compared to full order static condensation and provides a relative error of less than 5% in the relaxed stress.

97 MATHEMATICS AND COMPUTING↗

Thermal Stresses in Silicon-Web Growth

Mathematical model represents thermomechanical behavior of solid web of rectangular cross section as pulled from melt of same material. Intended for application to growth of silicon webs or ribbons, model general enough to be applicable to other materials as well. Shows suitable temperature distribution results in stress-free ribbon. Conceptual ribbon of rectangular cross section pulled from melt. Simplified configuration used to calculate thermal conditions for stress-free ribbon.

Wade, Ben K.↗

A Genomic Catalog of Stress Response Genes in Anaerobic Fungi for Applications in Bioproduction

Anaerobic fungi are a potential biotechnology platform to produce biomass-degrading enzymes. Unlike model fungi such as yeasts, stress responses that are relevant during bioprocessing have not yet been established for anaerobic fungi. In this work, we characterize both the heat shock and unfolded protein responses of four strains of anaerobic fungi ( Anaeromyces robustus, Caecomyces churrovis, Neocallimastix californiae , and Piromyces finnis ). The inositol-requiring 1 (Ire1) stress sensor, which typically initiates the fungal UPR, was conserved in all four genomes. However, these genomes also encode putative transmembrane kinases with catalytic domains that are similar to the metazoan stress-sensing enzyme PKR-like endoplasmic reticulum kinase (PERK), although whether they function in the UPR of anaerobic fungi remains unclear. Furthermore, we characterized the global transcriptional responses of Anaeromyces robustus and Neocallimastix californiae to a transient heat shock. Both fungi exhibited the hallmarks of ER stress, including upregulation of genes with functions in protein folding, ER-associated degradation, and intracellular protein trafficking. Relative to other fungi, the genomes of Neocallimastigomycetes contained the greatest gene percentage of HSP20 and HSP70 chaperones, which may serve to stabilize their asparagine-rich genomes. Taken together, these results delineate the unique stress response of anaerobic fungi, which is an important step toward their development as a biotechnology platform to produce enzymes and valuable biomolecules.

Swift, Candice L.↗

Second-order closure models for supersonic turbulent flows

Recent work on the development of a second-order closure model for high-speed compressible flows is reviewed. This turbulent closure is based on the solution of modeled transport equations for the Favre-averaged Reynolds stress tensor and the solenoidal part of the turbulent dissipation rate. A new model for the compressible dissipation is used along with traditional gradient transport models for the Reynolds heat flux and mass flux terms. Consistent with simple asymptotic analyses, the deviatoric part of the remaining higher-order correlations in the Reynolds stress transport equations are modeled by a variable density extension of the newest incompressible models. The resulting second-order closure model is tested in a variety of compressible turbulent flows which include the decay of isotropic turbulence, homogeneous shear flow, the supersonic mixing layer, and the supersonic flat-plate turbulent boundary layer. Comparisons between the model predictions and the results of physical and numerical experiments are quite encouraging.

Speziale, Charles G.↗

Second-order closure models for supersonic turbulent flows

Recent work by the authors on the development of a second-order closure model for high-speed compressible flows is reviewed. This turbulence closure is based on the solution of modeled transport equations for the Favre-averaged Reynolds stress tensor and the solenoidal part of the turbulent dissipation rate. A new model for the compressible dissipation is used along with traditional gradient transport models for the Reynolds heat flux and mass flux terms. Consistent with simple asymptotic analyses, the deviatoric part of the remaining higher-order correlations in the Reynolds stress transport equation are modeled by a variable density extension of the newest incompressible models. The resulting second-order closure model is tested in a variety of compressible turbulent flows which include the decay of isotropic turbulence, homogeneous shear flow, the supersonic mixing layer, and the supersonic flat-plate turbulent boundary layer. Comparisons between the model predictions and the results of physical and numerical experiments are quite encouraging.

Speziale, Charles G.↗

Exploring the Langtry-Menter Transition Model for High Speed Applications Using FUN3D

A series of Reynolds-averaged Navier-Stokes (RANS) simulations were performed using the FUN3D flow solver to explore the capabilities of the Langtry-Menter Shear-Stress Transport (LM-SST) transition model for predicting transition for aircraft inlet applications. Two geometries were simulated: a zero-pressure-gradient flat plate and an axisymmetric cone exposed to hypersonic flow. In addition to the transition-sensitized LM-SST model investigations, simulations were run with the one-equation Spalart-Allmaras (SA) and the two-equation Menter Shear-Stress Transport (SST-V) RANS models in fully turbulent mode to identify the natural RANS model transition behavior as a function of Mach number when executed in fully turbulent mode. The flat plate simulations showed that (1) the transition model was able to predict rapid transition at a freestream Mach number of 0.2, which is expected but (2) the predicted transition location moved downstream as the freestream Mach number was increased for the simulations that used the SST-V turbulence model. The latter is significant as it is usually assumed that one- and two-equation turbulence models will produce fully turbulent flow very near the boundary layer origin. The flat plate simulation freestream Mach number trend was confirmed with simulations using the Wind-US code, which also saw a similar trend when employing the SA turbulence model. For the axisymmetric cone simulations, the transition location was highly sensitive to the inflow turbulence levels. This is significant as the prediction of the transition location is crucial when trying to predict inlet performance, especially for hypersonic vehicle applications. It was also noted that the predicted transition location for the cone when using the SST-V turbulence model agreed well with the predicted transition location from the equivalent zero-pressure-gradient flat plate case.

Transition Model↗

Investigation of turbulent transport in an axisymmetric sudden expansion

Simultaneous two-component laser velocimeter measurements were made in the incompressible turbulent flowfield following an axisymmetric sudden expansion. Mean velocities, Reynolds stresses, and triple products were measured and are presented at axial positions ranging from x/H = 0.2-14. A balance of the turbulent kinetic energy in the flow was performed. The production, convection, and diffusion of turbulent kinetic energy were computed directly from the experimental data using central differencing. A specially designed correction lens was employed to correct for optical aberrations introduced by the circular tube. This lens system allowed the accurate simultaneous measurement of axial and radial velocities in the test section. The experimental measurements were compared to predictions generated by a code that employed the k-epsilon turbulence model. Agreement was good for mean axial velocities, turbulent kinetic energy, and turbulent shear stresses. However, the modeled turbulent normal stresses where in poor agreement with the measured values. The modeled diffusion of turbulent kinetic energy was underpredicted in the region between the shear layer and the centerline of the flow giving lower values of turbulent kinetic energy downstream of the potential core than measured.

Gould, Richard D.↗

Toughening of Boron Carbide Composites by Hierarchical Microstructuring

Due to a unique combination of properties including high hardness, low density, chemical and thermal stability, semi-conductivity, and high neutron absorption, boron carbide (B4C) is a potential candidate for various applications involving extreme environment. However, B4C’s current application is limited because of its low fracture toughness. In this study, a hierarchical microstructure design with features including TiB2 grains and graphite platelets was used to toughen B4C by simultaneously utilizing multiple toughening mechanisms including crack deflection, bridging, and micro-crack toughening. Using field-assisted sintering technology (FAST), B4C composites with dense and hierarchical microstructure were fabricated. Previously, the fracture toughness of fabricated B4C composites was measured at micro-scale using micro-indentation to have up to 56% improvement. In this work, the B4C composites’ fracture toughness was characterized at macro-scale using four-point bending methods and compared with previous results obtained at micro-scale. Micromechanics modeling of fracture behaviors for B4C-TiB2 composites was also performed to evaluate the contributions from experimentally observed toughening mechanisms. From four-point bending tests, B4C composites reinforced with both TiB2 grains (~15 vol%) and graphite platelets (~8.7 vol%) exhibited the highest fracture toughness enhancement from 2.38 to 3.65 MPa∙m^(1/2). The measured values were lower than those obtained using micro-indentation but maintained the general trends. The discrepancy between the indentation and four-point bending test results originated from the complex deformation behaviors triggered by the high contact load during indentation tests. Through micromechanics modeling, introduced thermal residual stress due to thermal expansion mismatch between B4C and TiB2, and weak interphases at B4C-TiB2 boundaries were identified as the main causes for experimentally observed toughness enhancement. These results proved the effectiveness of hierarchical microstructure designs for B4C toughening and can provide reference for the future design of B4C composites with optimized microstructures for further fracture toughness enhancement.

composites↗

Weld stresses beyond elastic limit: Materials discontinuity

When welded structures depend on properties beyond the elastic limit to qualify their ultimate safety factor, and weld-parent materials abruptly change at the interface, then stress discontinuity is inevitable. The stress concentration is mildly sensitive to material relative strain hardening and acutely sensitive to applied stress fields. Peak stresses occur on the weld surface, at the interface, and dissipate within a 0.01-inch band. When the stress is intense, the weld will always fracture at the interface. The analysis incorporates a classical mechanics model to more sharply define stress spikes within the bandwidth, and suggests a relative material index and Poisson's ratio related to strain hardening. Implications are discussed which are applicable to industries of high performance structures.

Verderaime, V.↗

Prediction of Size Effects in Notched Laminates Using Continuum Damage Mechanics

This paper examines the use of a continuum damage model to predict strength and size effects in notched carbon-epoxy laminates. The effects of size and the development of a fracture process zone before final failure are identified in an experimental program. The continuum damage model is described and the resulting predictions of size effects are compared with alternative approaches: the point stress and the inherent flaw models, the Linear-Elastic Fracture Mechanics approach, and the strength of materials approach. The results indicate that the continuum damage model is the most accurate technique to predict size effects in composites. Furthermore, the continuum damage model does not require any calibration and it is applicable to general geometries and boundary conditions.

Camanho, D. P.↗

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

Constituent properties, along with volume fraction, have a first order effect on the microscale fields within a composite material and influence the macroscopic response. Therefore, there is a need to assess the significance of stochastic variation in the constituent properties of composites at the higher scales. The effect of variability in the parameters controlling the time-dependent behavior, in a unidirectional SCS-6 SiC fiber-reinforced RBSN matrix composite lamina, on the residual stresses induced during processing is investigated numerically. The generalized method of cells micromechanics theory is utilized to model the ceramic matrix composite lamina using a repeating unit cell. The primary creep phases of the constituents are approximated using a Norton-Bailey, steady state, power law creep model. The effect of residual stresses on the proportional limit stress and strain to failure of the composite is demonstrated. Monte Carlo simulations were conducted using a normal distribution for the power law parameters and the resulting residual stress distributions were predicted.

residual stress↗

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

Constituent properties, along with volume fraction, have a first order effect on the microscale fields within a composite material and influence the macroscopic response. Therefore, there is a need to assess the significance of stochastic variation in the constituent properties of composites at the higher scales. The effect of variability in the parameters controlling the time-dependent behavior, in a unidirectional SCS-6 SiC fiber-reinforced RBSN matrix composite lamina, on the residual stresses induced during processing is investigated numerically. The generalized method of cells micromechanics theory is utilized to model the ceramic matrix composite lamina using a repeating unit cell. The primary creep phases of the constituents are approximated using a Norton-Bailey, steady state, power law creep model. The effect of residual stresses on the proportional limit stress and strain to failure of the composite is demonstrated. Monte Carlo simulations were conducted using a normal distribution for the power law parameters and the resulting residual stress distributions were predicted.

processing↗

Traveling Weather Disturbances in Mars Southern Extratropics: Sway of the Great Impact Basins

As on Earth, between late autumn and early spring on Mars middle and high latitudes within its atmosphere support strong mean thermal contrasts between the equator and poles (i.e. "baroclinicity"). Data collected during the Viking era and observations from both the Mars Global Surveyor (MGS) and Mars Reconnaissance Orbiter (MRO) indicate that this strong baroclinicity supports vigorous, large-scale eastward traveling weather systems (i.e. transient synoptic-period waves). Within a rapidly rotating, differentially heated, shallow atmosphere such as on Earth and Mars, such large-scale, extratropical weather disturbances are critical components of the global circulation. These wave-like disturbances act as agents in the transport of heat and momentum, and moreover generalized tracer quantities (e.g., atmospheric dust, water vapor and water-ice clouds) between low and high latitudes of the planet. The character of large-scale, traveling extratropical synoptic-period disturbances in Mars' southern hemisphere during late winter through early spring is investigated using a high-resolution Mars global climate model (Mars GCM). This global circulation model imposes interactively lifted (and radiatively active) dust based on a threshold value of the instantaneous surface stress. Compared to observations, the model exhibits a reasonable "dust cycle" (i.e. globally averaged, a more dusty atmosphere during southern spring and summer occurs). In contrast to their northern-hemisphere counterparts, southern synoptic-period weather disturbances and accompanying frontal waves have smaller meridional and zonal scales, and are far less intense synoptically. Influences of the zonally asymmetric (i.e. east-west varying) topography on southern large-scale weather disturbances are examined. Simulations that adapt Mars' full topography compared to simulations that utilize synthetic topographies emulating essential large-scale features of the southern middle latitudes indicate that Mars' transient barotropic/baroclinic eddies are significantly influenced by the great impact basins of this hemisphere (e.g., Argyre and Hellas). In addition, the occurrence of a southern storm zone in late winter and early spring is keyed particularly to the western hemisphere via orographic influences arising from the Tharsis highlands, and the Argyre and Hellas impact basins. Geographically localized transient-wave activity diagnostics are constructed that illuminate fundamental differences amongst such simulations and these are described.

Weather disturbances↗

STRANAL-PMC Version 2.0

Version 2.0 of the Strain Rate Dependent Analysis of Polymer Matrix Composites (STRANAL-PMC) software has been released. A prior version was reported in Analyzing Loads and Strains in Polymer- Matrix Composites (LEW-17227), NASA Tech Briefs, Vol. 26, No. 11 (November 2002), page 36. To recapitulate: Modified versions of constitutive equations of viscoplasticity of metals are used to represent deformation of a polymeric matrix. The equations are applied in a micromechanical approach, proceeding upward from slices of unit cells, through the ply level, to the laminate level. The constitutive equations are integrated in time by a Runge- Kutta technique. To predict the ultimate strength of each composite ply, failure criteria are implemented within the micromechanics. The inputs to STRANAL-PMC are the laminate geometry, properties of the fiber and matrix materials, and applied stress or strain versus time. The outputs are time-dependent stresses and strains at the slice, ply, and laminate levels. The improvements in version 2.0 include more rigorous representation of hydrostatic- stress effects in the matrix, refinement and extension of ply failure models, and capabilities to analyze transverse shear stresses. Version 2.0 can be implemented as a material-model code within transient dynamic finite-element codes.

Goldberg, Robert↗