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 163 records · Page 9

Traction free finite elements with the assumed stress hybrid model

An effective approach in the finite element analysis of the stress field at the traction free boundary of a solid continuum was studied. Conventional displacement and assumed stress finite elements were used in the determination of stress concentrations around circular and elliptical holes. Specialized hybrid elements were then developed to improve the satisfaction of prescribed traction boundary conditions. Results of the stress analysis indicated that finite elements which exactly satisfy the free stress boundary conditions are the most accurate and efficient in such problems. A general approach for hybrid finite elements which incorporate traction free boundaries of arbitrary geometry was formulated.

Kafie, Kurosh↗

Experimental study of a three-dimensional, shear-driven, turbulent boundary layer using a three-dimensional laser Doppler velocimeter

Experimental data have been obtained in a three-dimensional, shear-driven, turbulent boundary layer subjected to sudden transverse strain. Measurements made with a newly developed, three-dimensional laser Doppler velocimeter include all components of mean flow, turbulent Reynolds stresses, and triple-product correlations. Eddy viscosities, production, convection, turbulent diffusion, and dissipation (balance of kinetic energy equation) are extracted from the data. It was shown experimentally that eddy viscosity is strongly direction-dependent (not a scalar), indicating a need for full Reynolds stress modeling. Turbulence models for eddy viscosity, stress diffusion, pressure rate-of-strain, and anisotropy of Reynolds stresses are tested using the data. Calculations using a Reynolds stress equation model are compared with the data.

Driver, D. M.↗

Calculations of Supersonic and Hypersonic Flows using Compressible Wall Functions

The present paper presents a numerical procedure to calculate supersonic and hypersonic flows using the compressible law of the wall. The turbulence models under consideration include the Launder-Reece-Rodi-Gibson Reynolds-stress model and the k-epsilon model. The models coupled with the proposed wall function technique have been tested in both separated and unseparated flows. The flows include (1) an insulated flat plate flow over a range of Mach numbers, (2) a Mach 5 flat plate flow with cold wall conditions, (3) a two dimensional supersonic compression corner flow, (4) a hypersonic flow over an axisymmetric flare, and (5) a hypersonic flow over a 2-D compression corner. Results indicate that the wall function technique gives improved predictions of skin friction and heat transfer in separated flows compared with models using wall dampers. Predictions of the extent of separation are not improved over the wall damper models except with the Reynolds-stress model for the supersonic compression corner flow case.

Huang, P. G.↗

Utah FORGE 2-2439: A Multi-Component Approach to Characterizing In-Situ Stress: Laboratory, Modeling and Field Measurement - Workshop Presentation

This is a presentation on A Multi-Component Approach to Characterizing In-Situ Stress at the U.S DOE FORGE EGS Site: Laboratory, Modeling and Field Measurement project by Battelle [Columbus, OH], presented by Mark Kelley. The project's objective was to characterize stress in the Utah FORGE EGS reservoir using three methods: a laboratory rock-core stress estimation combined with a Machine Learning approach for estimation of in-situ stress from field sonic-log data, a field based in-situ measurement (min-frac) approach, and a modeling approach. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 7, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗

Turbulence Modeling Workshop

Advances in turbulence modeling are needed in order to calculate high Reynolds number flows near the onset of separation and beyond. To this end, the participants in this workshop made the following recommendations. (1) A national/international database and standards for turbulence modeling assessment should be established. Existing experimental data sets should be reviewed and categorized. Advantage should be taken of other efforts already under-way, such as that of the European Research Community on Flow, Turbulence, and Combustion (ERCOFTAC) consortium. Carefully selected "unit" experiments will be needed, as well as advances in instrumentation, to fill the gaps in existing datasets. A high priority should be given to document existing turbulence model capabilities in a standard form, including numerical implementation issues such as grid quality and resolution. (2) NASA should support long-term research on Algebraic Stress Models and Reynolds Stress Models. The emphasis should be placed on improving the length-scale equation, since it is the least understood and is a key component of two-equation and higher models. Second priority should be given to the development of improved near-wall models. Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) would provide valuable guidance in developing and validating new Reynolds-averaged Navier-Stokes (RANS) models. Although not the focus of this workshop, DNS, LES, and hybrid methods currently represent viable approaches for analysis on a limited basis. Therefore, although computer limitations require the use of RANS methods for realistic configurations at high Reynolds number in the foreseeable future, a balanced effort in turbulence modeling development, validation, and implementation should include these approaches as well.

Rubinstein, R.↗

The Essential Role of Westerly Wind Bursts in ENSO Dynamics and Extreme Events Quantified in Model “Wind Stress Shaving” Experiments

Abstract Westerly wind bursts (WWBs)—brief but strong westerly wind anomalies in the equatorial Pacific—are believed to play an important role in El Niño–Southern Oscillation (ENSO) dynamics, but quantifying their effects is challenging. Here, we investigate the cumulative effects of WWBs on ENSO characteristics, including the occurrence of extreme El Niño events, via modified coupled model experiments within Community Earth System Model (CESM1) in which we progressively reduce the impacts of wind stress anomalies associated with model-generated WWBs. In these “wind stress shaving” experiments we limit momentum transfer from the atmosphere to the ocean above a preset threshold, thus “shaving off” wind bursts. To reduce the tropical Pacific mean state drift, both westerly and easterly wind bursts are removed, although the changes are dominated by WWB reduction. As we impose progressively stronger thresholds, both ENSO amplitude and the frequency of extreme El Niño decrease, and ENSO becomes less asymmetric. The warming center of El Niño shifts westward, indicating less frequent and weaker eastern Pacific (EP) El Niño events. Removing most wind burst–related wind stress anomalies reduces ENSO mean amplitude by 22%. The essential role of WWBs in the development of extreme El Niño events is highlighted by the suppressed eastward migration of the western Pacific warm pool and hence a weaker Bjerknes feedback under wind shaving. Overall, our results reaffirm the importance of WWBs in shaping the characteristics of ENSO and its extreme events and imply that WWB changes with global warming could influence future ENSO.

Meteorology & Atmospheric Sciences↗

Finite element modeling of stress in the Nazca plate - Driving forces and plate boundary earthquakes

The state of stress within the Nazca plate due to plate driving forces and large plate boundary earthquakes has been analyzed by applying a finite element method using the wave front solution technique to models of the intraplate stress field in a single plate using a refined grid. Although only static elastic models have been explicitly calculated, certain limiting cases of an elastic plate over a viscous asthenosphere were also treated. A state of nearly east-west compression inferred from the source mechanism of thrust earthquakes in the interior of the plate requires ridge pushing forces. The net pulling force on the oceanic plate by the subducted slab has a maximum value comparable to pushing forces. The estimated horizontal deviatoric stress in intraplate regions, based on potential forces associated with the ridge, is on the order of a few hundred bars. The intraplate stress field in the region of the 1960 earthquake may change by a few tens of bars at most once the asthenosphere has relaxed, with changes on the order of one bar occurring at greater distances into the plate. The changes in the intraplate stress field are probably not noticeable unless the lithosphere is near failure.

Richardson, R. M.↗

Numerical Solutions for the CAWAPI Configuration on Structured Grids at NASA LaRC, United States

In this chapter numerical simulations of the flow around F-16XL are performed as a contribution to the Cranked Arrow Wing Aerodynamic Project International (CAWAPI) using the PAB3D CFD code. Two turbulence models are used in the calculations: a standard k-epsilon model, and the Shih-Zhu-Lumley (SZL) algebraic stress model. Seven flight conditions are simulated for the flow around the F-16XL where the free stream Mach number varies from 0.242 to 0.97. The range of angles of attack varies from 0 deg to 20 deg. Computational results, surface static pressure, boundary layer velocity profiles, and skin friction are presented and compared with flight data. Numerical results are generally in good agreement with flight data, considering that only one grid resolution is utilized for the different flight conditions simulated in this study. The Algebraic Stress Model (ASM) results are closer to the flight data than the k-epsilon model results. The ASM predicted a stronger primary vortex, however, the origin of the vortex and footprint is approximately the same as in the k-epsilon predictions.

Elmiligui, Alaa A.↗

Reynolds stress closure modeling in wall-bounded flows

This report describes two projects. Firstly, a Reynolds stress closure for near-wall turbulence is described. It was motivated by the simpler k-epsilon-(v-bar(exp 2)) model described in last year's annual research brief. Direct Numerical Simulation of three-dimensional channel flow shows a curious decrease of the turbulent kinetic energy. The second topic of this report is a model which reproduces this effect. That model is described and used to discuss the relevance of the three dimensional channel flow simulation to swept wing boundary layers.

Durbin, Paul A.↗

Elliptic Relaxation of a Tensor Representation for the Redistribution Terms in a Reynolds Stress Turbulence Model

A formulation to include the effects of wall proximity in a second-moment closure model that utilizes a tensor representation for the redistribution terms in the Reynolds stress equations is presented. The wall-proximity effects are modeled through an elliptic relaxation process of the tensor expansion coefficients that properly accounts for both correlation length and time scales as the wall is approached. Direct numerical simulation data and Reynolds stress solutions using a full differential approach are compared for the case of fully developed channel flow.

Carlson, J. R.↗

Turbulence modelling of flow fields in thrust chambers

Following the consensus of a workshop in Turbulence Modelling for Liquid Rocket Thrust Chambers, the current effort was undertaken to study the effects of second-order closure on the predictions of thermochemical flow fields. To reduce the instability and computational intensity of the full second-order Reynolds Stress Model, an Algebraic Stress Model (ASM) coupled with a two-layer near wall treatment was developed. Various test problems, including the compressible boundary layer with adiabatic and cooled walls, recirculating flows, swirling flows, and the entire SSME nozzle flow were studied to assess the performance of the current model. Detailed calculations for the SSME exit wall flow around the nozzle manifold were executed. As to the overall flow predictions, the ASM removes another assumption for appropriate comparison with experimental data to account for the non-isotropic turbulence effects.

Chen, C. P.↗

Development of an algebraic stress/two-layer model for calculating thrust chamber flow fields

Following the consensus of a workshop in Turbulence Modeling for Liquid Rocket Thrust Chambers, the current effort was undertaken to study the effects of second-order closure on the predictions of thermochemical flow fields. To reduce the instability and computational intensity of the full second-order Reynolds Stress Model, an Algebraic Stress Model (ASM) coupled with a two-layer near wall treatment was developed. Various test problems, including the compressible boundary layer with adiabatic and cooled walls, recirculating flows, swirling flows and the entire SSME nozzle flow were studied to assess the performance of the current model. Detailed calculations for the SSME exit wall flow around the nozzle manifold were executed. As to the overall flow predictions, the ASM removes another assumption for appropriate comparison with experimental data, to account for the non-isotropic turbulence effects.

Chen, C. P.↗

Examination of various turbulence models for application in liquid rocket thrust chambers

There is a large variety of turbulence models available. These models include direct numerical simulation, large eddy simulation, Reynolds stress/flux model, zero equation model, one equation model, two equation k-epsilon model, multiple-scale model, etc. Each turbulence model contains different physical assumptions and requirements. The natures of turbulence are randomness, irregularity, diffusivity and dissipation. The capabilities of the turbulence models, including physical strength, weakness, limitations, as well as numerical and computational considerations, are reviewed. Recommendations are made for the potential application of a turbulence model in thrust chamber and performance prediction programs. The full Reynolds stress model is recommended. In a workshop, specifically called for the assessment of turbulence models for applications in liquid rocket thrust chambers, most of the experts present were also in favor of the recommendation of the Reynolds stress model.

Hung, R. J.↗

Research on complex turbulent flows at the UW (overall effort)

Development of a pressure-strain model, an algebraic stress model, and wall functions appropriate for flows with spanwise variations in the local wall shear stress are accomplished. Furthermore, a hot-wire measurement technique was also developed for determining the local mean velocity and Reynolds stresses in a complex flow. Experiments were performed on supersonic and subsonic turbulent flow in a square duct, flow about a strut-endwall, flow within a transition duct, and on co-flowing annular jets with swirl. All results are presented in a viewgraph format.

Gessner, F. B.↗

Reduced Stress Tensor and Dissipation and the Transport of Lamb Vector

We develop a methodology to ensure that the stress tensor, regardless of its number of independent components, can be reduced to an exactly equivalent one which has the same number of independent components as the surface force. It is applicable to the momentum balance if the shear viscosity is constant. A direct application of this method to the energy balance also leads to a reduction of the dissipation rate of kinetic energy. Following this procedure, significant saving in analysis and computation may be achieved. For turbulent flows, this strategy immediately implies that a given Reynolds stress model can always be replaced by a reduced one before putting it into computation. Furthermore, we show how the modeling of Reynolds stress tensor can be reduced to that of the mean turbulent Lamb vector alone, which is much simpler. As a first step of this alternative modeling development, we derive the governing equations for the Lamb vector and its square. These equations form a basis of new second-order closure schemes and, we believe, should be favorably compared to that of traditional Reynolds stress transport equation.

Wu, Jie-Zhi↗

Strength and deformation properties of basaltic lava flows on planetary surfaces

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

Schultz, Richard A.↗

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

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

King, James↗