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 577 records · Page 32

Steady state stresses in ribbon parachute canopies

An experimental study of the steady state stresses in model ribbon parachute canopies is presented. The distribution of circumferential stress was measured in the horizontal ribbons of two parachutes using Omega sensors. Canopy pressure distributions and overall drag were also measured. Testing was conducted in the University of Minnesota Low-Speed Wind Tunnel at dynamic pressures ranging from 1.0 to 1.5 inches of water. The stresses in the parachute canopies were calculated using the parachute structural analysis code, CANO. It was found that the general shape of the measured and calculated stress distributions was fairly similar; however, the measured stresses were somewhat less than the calculated stresses.

Garrard, W. L.↗

Overview and highlights of Early Warning and Crop Condition Assessment project

Work of the Early Warning and Crop Condition Assessment (EW/CCA) project, one of eight projects in the Agriculture and Resources Inventory Surveys Through Aerospace Remote Sensing (AgRISTARS), is reviewed. Its mission, to develop and test remote sensing techniques that enhance operational methodologies for crop condition assessment, was in response to initiatives issued by the Secretary of Agriculture. Meteorologically driven crop stress indicator models have been developed or modified for wheat, maize, grain sorghum, and soybeans. These models provide early warning alerts of potential or actual crop stresses due to water deficits, adverse temperatures, and water excess that could delay planting or harvesting operations. Recommendations are given for future research involving vegetative index numbers and the NOAA and Landsat satellites.

Boatwright, G. O.↗

Vertical movement in mare basins - Relation to mare emplacement, basin tectonics, and lunar thermal history

A hypothesis is developed that the spatial and temporal distributions of linear rilles and mare ridges in the mare regions of the moon are the products of two superposed stress systems: a local stress due to lithospheric loading by basalt fill in the mare basin, and a global thermal stress due to the thermal history of the lunar interior. Quantitative models for stress in the Serenitatis basin area, including the global thermal stress associated with lunar thermal history, are presented and used to account for the spatial distribution, the orientations, and the formation times of rilles and ridges in that region. The models provide a simple explanation for the localization of the most recent eruption sites of mare basalt magmas to mare basin edges.

Solomon, S. C.↗

Current capabilities for simulating the extreme distortion of thin structures subjected to severe impacts

The explicit transient dynamics technology in use today for simulating the impact and subsequent transient dynamic response of a structure has its origins in the 'hydrocodes' dating back to the late 1940's. The growth in capability in explicit transient dynamics technology parallels the growth in speed and size of digital computers. Computer software for simulating the explicit transient dynamic response of a structure is characterized by algorithms that use a large number of small steps. In explicit transient dynamics software there is a significant emphasis on speed and simplicity. The finite element technology used to generate the spatial discretization of a structure is based on a compromise between completeness of the representation for the physical processes modelled and speed in execution. That is, since it is expected in every calculation that the deformation will be finite and the material will be strained beyond the elastic range, the geometry and the associated gradient operators must be reconstructed, as well as complex stress-strain models evaluated at every time step. As a result, finite elements derived for explicit transient dynamics software use the simplest and barest constructions possible for computational efficiency while retaining an essential representation of the physical behavior. The best example of this technology is the four-node bending quadrilateral derived by Belytschko, Lin and Tsay. Today, the speed, memory capacity and availability of computer hardware allows a number of the previously used algorithms to be 'improved.' That is, it is possible with today's computing hardware to modify many of the standard algorithms to improve their representation of the physical process at the expense of added complexity and computational effort. The purpose is to review a number of these algorithms and identify the improvements possible. In many instances, both the older, faster version of the algorithm and the improved and somewhat slower version of the algorithm are found implemented together in software. Specifically, the following seven algorithmic items are examined: the invariant time derivatives of stress used in material models expressed in rate form; incremental objectivity and strain used in the numerical integration of the material models; the use of one-point element integration versus mean quadrature; shell elements used to represent the behavior of thin structural components; beam elements based on stress-resultant plasticity versus cross-section integration; the fidelity of elastic-plastic material models in their representation of ductile metals; and the use of Courant subcycling to reduce computational effort.

Key, Samuel W.↗

Modeling of Turbulent Swirling Flows

Aircraft engine combustors generally involve turbulent swirling flows in order to enhance fuel-air mixing and flame stabilization. It has long been recognized that eddy viscosity turbulence models are unable to appropriately model swirling flows. Therefore, it has been suggested that, for the modeling of these flows, a second order closure scheme should be considered because of its ability in the modeling of rotational and curvature effects. However, this scheme will require solution of many complicated second moment transport equations (six Reynolds stresses plus other scalar fluxes and variances), which is a difficult task for any CFD implementations. Also, this scheme will require a large amount of computer resources for a general combustor swirling flow. This report is devoted to the development of a cubic Reynolds stress-strain model for turbulent swirling flows, and was inspired by the work of Launder's group at UMIST. Using this type of model, one only needs to solve two turbulence equations, one for the turbulent kinetic energy k and the other for the dissipation rate epsilon. The cubic model developed in this report is based on a general Reynolds stress-strain relationship. Two flows have been chosen for model evaluation. One is a fully developed rotating pipe flow, and the other is a more complex flow with swirl and recirculation.

Shih, Tsan-Hsing↗

The present state and future direction of second order closure models for compressible flows

The topics are presented in viewgraph form and include: (1) Reynolds stress closure models; (2) Favre averages and governing equations; (3) the model for the deviatoric part of the pressure-strain rate correlation; (4) the SSG pressure-strain correlation model; (5) a compressible turbulent dissipation rate model; (6) variable viscosity effects; (7) near-wall stiffness problems; (8) models of the Reynolds mass and heat flux; and (9) a numerical solution of the compressible turbulent transport equation.

Gatski, Thomas B.↗

Crack-surface displacements for two symmetric cracks emanating from a circular hole in an infinite plate

An analytical model is presented for the surface displacement associated with two symmetric cracks radiating from a circular hole in a plate undergoing remote uniform stress. The model equations account for the stress intensity factor (Kic), the boundary conditions enveloping the hole, and the effects of symmetric forces concentrated around the hole. The displacement(s) of the surface are calculated at a given distance from the hole along the crack(s).

Mall, S.↗

Improved two-equation k-omega turbulence models for aerodynamic flows

Two new versions of the k-omega two-equation turbulence model will be presented. The new Baseline (BSL) model is designed to give results similar to those of the original k-omega model of Wilcox, but without its strong dependency on arbitrary freestream values. The BSL model is identical to the Wilcox model in the inner 50 percent of the boundary-layer but changes gradually to the high Reynolds number Jones-Launder k-epsilon model (in a k-omega formulation) towards the boundary-layer edge. The new model is also virtually identical to the Jones-Lauder model for free shear layers. The second version of the model is called Shear-Stress Transport (SST) model. It is based on the BSL model, but has the additional ability to account for the transport of the principal shear stress in adverse pressure gradient boundary-layers. The model is based on Bradshaw's assumption that the principal shear stress is proportional to the turbulent kinetic energy, which is introduced into the definition of the eddy-viscosity. Both models are tested for a large number of different flowfields. The results of the BSL model are similar to those of the original k-omega model, but without the undesirable freestream dependency. The predictions of the SST model are also independent of the freestream values and show excellent agreement with experimental data for adverse pressure gradient boundary-layer flows.

Menter, Florian R.↗

The effect of entrance flow development on vortex formation and wall shear stress in a curved artery model

In this work, we numerically investigate the effect of entrance condition on the spatial and temporal evolution of multiple three-dimensional vortex pairs and the wall shear stress distribution in a curved artery model. We perform this study using a Newtonian blood-analog fluid subjected to a pulsatile flow with two inflow conditions. The first flow condition is fully developed while the second condition is undeveloped (i.e., uniform). We discuss the connection along the axial direction between regions of organized vorticity observed at various cross sections of the model and compare results between the different entrance conditions. We model a human artery with a simple, rigid 180° curved pipe with a circular cross section and constant curvature, neglecting the effects of taper, torsion, and elasticity. Numerical results are computed from a discontinuous high-order spectral element flow solver. The flow rate used in this study is physiological. We observe differences in secondary flow patterns, especially during the deceleration phase of the physiological waveform where multiple vortical structures of both Dean-type and Lyne-type coexist. The results indicate that decreased axial velocities under an undeveloped condition produce smaller secondary flows that ultimately inhibit growth of any interior flow vortices. We highlight the effect of the entrance condition on the formation of these structures and subsequent appearance of abnormal inner wall shear stresses, which suggest there may be a lower prevalence of cardiovascular disease in curved arteries where the flow is rather undeveloped—a potentially physiologically significant result to help understand the influence of blood flow development on disease.

60 APPLIED LIFE SCIENCES↗

Zonal Two Equation Kappa-Omega Turbulence Models for Aerodynamic Flows

Two new versions of the kappa-omega two-equation turbulence model will be presented. The new Baseline (BSL) model is designed to give results similar to those of the original kappa-omega model of Wilcox, but without its strong dependency on arbitrary freestream values. The BSL model is identical to the Wilcox model in the inner 50% of the boundary-layer but changes gradually to the standard kappa-epsilon model (in a kappa- omega formulation) towards the boundary-layer edge. The free shear layers. The second version of the model is called Shear-Stress Transport (SST) model. It is a variation of the BSL model with the additional ability to account for the transport of the principal turbulent shear stress in adverse pressure gradient boundary-layers. The model is based on Bradshaw's assumption that the principal shear-stress is proportional to the turbulent kinetic energy, which is introduced into the definition of the eddy-viscosity. Both models are tested for a large number of different flowfields. The results of the BSL model are similar to those of the original kappa-omega model, but without the undesirable freestream dependency. The predictions of the SST model are also independent of the freestream values but show better agreement with experimental data for adverse pressure gradient boundary-layer flows.

Menter, Florian R.↗

Predicting Failure Using Deep Learning SAND Report

Accurate prediction of ductile failure is critical to Sandia’s NW mission, but the models are computationally heavy. The costs of including high-fidelity physics and mechanics that are germane to the failure mechanisms are often too burdensome for analysts either because of the person-hours it requires to input them or because of the additional computational time, or both. In an effort to deliver analysts a tool for representing these phenomena with minimal impact to their existing workflow, our project sought to develop modern data-driven methods that would add microstructural information to business-as-usual calculations and expedite failure predictions. The goal is a tool that receives as input a structural model with stress and strain fields, as well as a machine-learned model, and output predictions of structural response in time, including failure. As such, our project spent substantial time performing high-fidelity, three-dimensional experiments to elucidate materials mechanisms of void nucleation and evolution. We developed crystal-plasticity finite-element models from the experimental observations to enrich the findings with fields not readily measured. We developed engineering length-scale simulations of replicated test specimens to understand how the engineering fields evolve in the presence of fine-scale defects. Finally, we developed deep learning convolutional neural networks, and graph-based neural networks to encode the findings of the experiments and simulations and make forward predictions in time for structural performance. This project demonstrated the power of data-driven methods for model development, which have the potential to vastly increase both the accuracy and speed of failure predictions. These benefits and the methods necessary to develop them are highlighted in this report. However, many challenges remain to implementing these in real applications, and these are discussed along with potential methods for overcoming them.

97 MATHEMATICS AND COMPUTING↗

SiC/SiC Life Modeling at Intermediate Temperatures under NASA UEET

Within the NASA UEET program, physics-based residual strength and rupture life models are being developed for a variety of time-temperature-stress-environment conditions for SiC/SiC components such as combustor liners and vanes. During service, a tensile stress condition that could be life-limiting typically occurs on the 'cold side' of these components due to thermal gradients and attachment schemes. An intermediate temperature (600 to 1000 C) stress-rupture model has been developed based on the underlying physical processes, i.e., oxygen penetration through matrix cracks, oxidation of the BN interphase, and subsequent fiber-fiber fusion within crack-bridging tows that eventually fail in time due to time-dependent strength degradation of the weakest fibers. Modifications to the model will be discussed to incorporate different environmental conditions as well as changes in oxidation kinetics due to improved interphase concepts.

Morscher, Gregory N.↗

Alkali-silica reaction expansion model for confined concrete with stress-dependency and casting direction anisotropy

Alkali-silica reaction (ASR) is a deleterious chemical reaction between alkali hydroxyl ions and types of silica found in some aggregates of concrete. Owners and regulators of nuclear power plants aim to ensure the safety of the concrete structures with optimal maintenance strategies. A new model was developed for predicting the expansion of concrete structures affected by alkali-silica reaction. The model includes a novel combination of existing models as an alkali-silica reaction advancement model, a casting direction anisotropic expansion model, a stress-dependent anisotropic expansion model, and a material property evolution model dependent on the degree of ASR expansion. The model parameters were calibrated based on existing literature data and data generated by previous efforts of this study. The calibrated model was then validated with the experiments carried out in previous efforts of this study. The model was shown to accurately predict the ASR-expansion of large-scale reinforced concrete specimens with confinement.

36 MATERIALS SCIENCE↗

A Unified Framework to Reconcile Different Approaches of Modeling Transpiration Response to Water Stress: Plant Hydraulics, Supply Demand Balance, and Empirical Soil Water Stress Function

Plant responses to water stress is a major uncertainty to predicting terrestrial ecosystem sensitivity to drought. Different approaches have been developed to represent plant water stress. Empirical approaches (the empirical soil water stress (or Beta) function and the supply-demand balance scheme) have been widely used for many decades; more mechanistic based approaches, that is, plant hydraulic models (PHMs), were increasingly adopted in the past decade. However, the relationships between them—and their underlying connections to physical processes—are not sufficiently understood. This limited understanding hinders informed decisions on the necessary complexities needed for different applications, with empirical approaches being mechanistically insufficient, and PHMs often being too complex to constrain. Here we introduce a unified framework for modeling transpiration responses to water stress, within which we demonstrate that empirical approaches are special cases of the full PHM, when the plant hydraulic parameters satisfy certain conditions. We further evaluate their response differences and identify the associated physical processes. Finally, we propose a methodology for assessing the necessity of added complexities of the PHM under various climatic conditions and ecosystem types, with case studies in three typical ecosystems: a humid Midwestern cropland, a semi-arid evergreen needleleaf forest, and an arid grassland. Notably, Beta function overestimates transpiration when VPD is high due to its lack of constraints from hydraulic transport and is therefore insufficient in high VPD environments. With the unified framework, we envision researchers can better understand the mechanistic bases of and the relationships between different approaches and make more informed choices.

54 ENVIRONMENTAL SCIENCES↗

Postseismic viscoelastic deformation and stress. Part 2: Stress theory and computation; dependence of displacement, strain, and stress on fault parameters

A viscoelastic model for deformation and stress associated with earthquakes is reported. The model consists of a rectangular dislocation (strike slip fault) in a viscoelastic layer (lithosphere) lying over a viscoelastic half space (asthenosphere). The time dependent surface stresses are analyzed. The model predicts that near the fault a significant fraction of the stress that was reduced during the earthquake is recovered by viscoelastic softening of the lithosphere. By contrast, the strain shows very little change near the fault. The model also predicts that the stress changes associated with asthenospheric flow extend over a broader region than those associated with lithospheric relaxation even though the peak value is less. The dependence of the displacements, stresses on fault parameters studied. Peak values of strain and stress drop increase with increasing fault height and decrease with fault depth. Under many circumstances postseismic strains and stresses show an increase with decreasing depth to the lithosphere-asthenosphere boundary. Values of the strain and stress at distant points from the fault increase with fault area but are relatively insensitive to fault depth.

Cohen, S. C.↗

High Temperature High Vacuum Mechanical Property Assessment of Zirconium Nuclear Fuel Cladding

This report presents the mechanical characterization of a specific Zry-4 cladding batch serving as the foundation for a diverse range of fuel performance research at Oak Ridge National Laboratory (ORNL). This effort supports research needs for the U.S. Department of Energy (DOE), particularly regarding evaluating accident tolerant fuel (ATF) cladding coating concepts, expanding understanding of cladding response to loss-of-coolant accidents (LOCA) transients, refining post-critical heat flux (CHF) limits (t@T), and upcoming irradiation campaigns. The central objective was to define the baseline performance of the substrate Zircaloy-4 (Zry-4) material leveraged across ORNL Advanced Fuel Campaign (AFC) efforts through controlled high-temperature vacuum tensile testing. This work begins to address gaps in existing models where implementation based on nominal heat-treatment labels, such as stress relief annealed (SRA), often fail to capture the interplay of recovery, recrystallization, and grain growth. To quantify this, data was benchmarked against the Pacific Northwest National Laboratory (PNNL) stress strain model to determine where this material falls in comparison to assumed values for materials in the same heat treatment regime. Analysis of the tensile data revealed that this specific SRA batch exhibits a transitional microstructural state best described by an effective cold-work (CW) parameter of 0.09, diverging from the previous estimation of 0.5 for SRA materials. Additionally, comparative testing of Cr coated specimens demonstrated no distinct difference in axial strength relative to the bare substrate. This suggests that the strengthening benefits of Cr coatings observed in burst scenarios are driven by residual stress mechanisms acting solely in the hoop direction.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Global analysis of ocean surface wind and wind stress using a general circulation model and Seasat scatterometer winds

Instantaneous and 15-day time-averaged fields of surface wind, wind stress, curl of the wind stress, and wind divergence are presented. These fields are derived from the Goddard Laboratory for Atmospheres four-dimensional analysis/forecast cycle, for the period September 6-30, 1978, using conventional data, satellite temperature soundings, cloud-track winds, and subjectively dealiased Seasat scatterometer winds.

Kalnay, E.↗

New generalized rheological model for lubrication of a ball spinning in a nonconforming groove

The elastohydrodynamic theory for predicting the spinning friction of a ball in a nonconforming groove was modified to incorporate a rheological model. The rheological model is based on the exponential pressure viscosity relation for low shear stresses, but at high shear rates and pressures, the relation is altered to one in which the shear stress is porportional to the normal stress. The model was fitted to experimental spinning torques for four different lubricants: a synthetic paraffinic lubricant, di-2-ethylhexyl sebacate, a super-refined naphthenic mineral oil, and a polyphenyl ether (5P4E). Good agreement between the model and experiment was found.

Allen, C. W.↗