Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRESS TENSOR”

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 127 records · Page 7

Turbulence Modeling for Shock Wave/Turbulent Boundary Layer Interactions

Accurate aerodynamic computational predictions are essential for the safety of space vehicles, but these computations are of limited accuracy when large pressure gradients are present in the flow. The goal of the current project is to improve the state of compressible turbulence modeling for high speed flows with shock wave / turbulent boundary layer interactions (SWTBLI). Emphasis will be placed on models that can accurately predict the separated region caused by the SWTBLI. These flows are classified as nonequilibrium boundary layers because of the very large and variable adverse pressure gradients caused by the shock waves. The lag model was designed to model these nonequilibrium flows by incorporating history effects. Standard one- and two-equation models (Spalart Allmaras and SST) and the lag model will be run and compared to a new lag model. This new model, the Reynolds stress tensor lag model (lagRST), will be assessed against multiple wind tunnel tests and correlations. The basis of the lag and lagRST models are to preserve the accuracy of the standard turbulence models in equilibrium turbulence, when the Reynolds stresses are linearly related to the mean strain rates, but create a lag between mean strain rate effects and turbulence when nonequilibrium effects become important, such as in large pressure gradients. The affect this lag has on the results for SWBLI and massively separated flows will be determined. These computations will be done with a modified version of the OVERFLOW code. This code solves the RANS equations on overset grids. It was used for this study for its ability to input very complex geometries into the flow solver, such as the Space Shuttle in the full stack configuration. The model was successfully implemented within two versions of the OVERFLOW code. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWBLI assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.↗

A Non-Invasive Approach for Elucidating the Spatial Distribution of In Situ Stress in Deep Subsurface Geologic Formations Considered for CO 2 Storage (Task 5 Report Field Scale Stress Modeling)

This report describes research accomplishments achieved with funding provided through Department of Energy (DOE) Contract DE-FE0031686, for the project “A Non-Invasive Approach for Elucidating the Spatial Distribution of In Situ Stress in Deep Subsurface Geologic Formations Considered for CO 2 Storage.” The purpose of this DOE funding is to develop non-invasive methods to obtain important information about subsurface stresses. The overall research goal of the project was to develop and improve methods for determining the spatial distribution of stresses, including magnitude and orientation of the three principal stress components in the subsurface, based on new methods that extract stress information from seismic data combined with well measurements, extended with well tests and logs, and unified in a numerical model that permits computation of the full stress tensor throughout the domains under investigation.

58 GEOSCIENCES↗

Behavior of a turbulent boundary layer subjected to sudden transverse strain

Data from two experiments on the development of the components of the Reynolds stress tensor after a sudden application of transverse strain are compared. Computations were based on four different turbulence models: a first-order mixing length model, a second-order two-equation eddy viscosity model, and two second-order Reynolds stress models. The second-order models do not produce dramatic improvements over the simple mixing length model. The Reynolds stress models still need development to represent the physics of shear-strained turbulence well. Finally, it is demonstrated that the assumption of a scalar eddy viscosity, often used in engineering calculations of three-dimensional boundary layers, is quite reasonable.

Higuchi, H.↗

Second moment closure analysis of the backstep flow database

A Second Moment Closure computation (SMC) is compared in detail with the Direct Numerical Simulation (DNS) data of Le and Moin for the backstep flow at Re = 5,000 in an attempt to understand why the intensity of the backflow and, consequently, the friction coefficient in the recirculation bubble are severely underestimated. The data show that this recirculation bubble is far from being laminar except in the very near wall layer. A novel 'differential a priori' procedure was used, in which the full transport equation for one isolated component of the Reynolds stress tensor was solved using DNS data as input. Conclusions are then different from what would have been deduced by comparing a full simulation to a DNS. One cause of discrepancy was traced back to insufficient transfer of energy to the normal stress by pressure strain, but was not cured. A significant finding, confirmed by the DNS data in the core region of a channel flow, is that the coefficient that controls destruction of dissipation, C epsilon(sub 2), should be decreased by a factor of 2 when production is vanishing. This is also the case in the recirculation bubble, and a new formulation has cured 25% of the backflow discrepancy.

Parneix, S.↗

Organized kinetic energy backscatter in the hurricane boundary layer from radar measurements

The fluid mechanics of hurricanes strongly depends on boundary layer energetics due to the warm-core nature of the system with peak velocities located at lower levels. One barrier that has inhibited a more complete characterization of energy transfer in the boundary layer is a lack of observations that resolve large, turbulent eddies. In particular, the occurrence and structure of upscale energy transfer (backscatter) in the hurricane boundary layer as well as the effects of backscatter on the vortex intensity are unknown. The analysis presented here of very high-resolution, three-dimensional wind observations from Hurricane Rita (2005) at peak intensity reveals large regions of organized backscatter in the boundary layer associated with coherent, turbulent eddies. Strong forwardscatter is also found next to the backscatter regions due to the interaction between adjacent eddies. Two components of the stress tensor are primarily responsible for this alternating scatter structure, as shown by large correlation coefficients between the fields: the radial–vertical component (τ13) and the azimuthal–vertical component (τ23) with average correlations of 79% and 49 %, respectively. The Leonard, Reynolds and cross-term stress components are also provided. The impact of the sub-filter-scale energy transfer is estimated by computing the kinetic energy budget for the resolved-scale and eddy-scale motions. The results show that the sub-filter-scale energy transfer term is of the same order as the other terms in the eddy-scale budgets, contributing between 16% and 40% to the local time tendency with an average contribution of approximately 30 %. These results indicate that the coherent turbulent eddies can affect the vortex dynamics through wave–wave nonlinear interactions, which can subsequently influence the wave–mean flow interactions. This is the first study to examine the full sub-filter-scale

Sydney Sroka↗

The Turbulent Stress Spectrum in the Inertial and Subinertial Ranges

For velocity and magnetic fields, the turbulent pressure and, more generally, the squared fields such as the components of the turbulent stress tensor, play important roles in astrophysics. For both one and three dimensions, we derive the equations relating the energy spectra of the fields to the spectra of their squares. We solve the resulting integrals numerically and show that for turbulent energy spectra of Kolmogorov type, the spectral slope of the stress spectrum is also of Kolmogorov type. For shallower turbulence spectra, the slope of the stress spectrum quickly approaches that of white noise, regardless of how blue the spectrum of the field is. For fully helical fields, the stress spectrum is elevated by about a factor of two in the subinertial range, while that in the inertial range remains unchanged. We discuss possible implications for understanding the spectrum of primordial gravitational waves from causally generated magnetic fields during cosmological phase transitions in the early universe. We also discuss potential diagnostic applications to the interstellar medium, where polarization and scintillation measurements characterize the square of the magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗

3-D laser anemometer measurements in a labyrinth seal

The flow field inside a seven cavity labyrinth seal with a 0.00127 m clearance was measured using a 3-D laser Doppler anemometer system. Through the use of this system, the mean velocity vector and the entire Reynolds stress tensor distributions were measured for the first, third, fifth, and seventh cavities of the seal. There was one large recirculation region present in the cavity for the flow condition tested, Re = 28,000 and Ta = 7,000. The axial and radial mean velocities as well as all of the Reynolds stress term became cavity independent by the third cavity. The azimuthal mean velocity varied from cavity to cavity with its magnitude increasing as the flow progressed downstream.

Morrison, G. L.↗

Enhanced Microseismicity During Production Pumping Cessation at the San Emidio Geothermal Field (Nevada, USA) in December 2016

Abstract Tectonic activity, geothermal fluids, and microseismic events (MSEs) tend to occur in similar locations as a result of spatiotemporal changes in the subsurface stress state. To quantify this association, we analyze data from a dense seismic array deployed at the San Emidio geothermal field, Nevada for 1 week in December 2016 to coincide with a 19.45‐hr shutdown of all injection and production pumping operations. 123 MSEs were detected, of which 101 occurred during the shutdown. The spatial association of the MSEs with the production wells suggests a causal relationship between the production cessation and the MSEs. Here we performed a detailed analysis to investigate reservoir material properties, distribution of seismically activated faults, and local stress state. We determined the hypocenters, magnitudes, and focal mechanisms for the MSEs, P‐wave tomographic velocity model, and local stress tensor. The results show that most MSEs occurred near the production wells. Magnitudes fall between −2.2 and 0.0 with larger events located closer to the production wells. Most MSEs occurred within a westward‐dipping normal fault zone in the reservoir associated with anomalously low P‐wave velocity values. The focal mechanism and stress inversion results show predominantly normal faulting with the maximum horizontal stress oriented north‐south. We suggest that the MSEs during shutdown were triggered on pre‐existing, small‐scale, critically stressed fault patches in the reservoir as the pore pressure increased around the production wells when the production pumping ceased. We interpret the larger MSE magnitudes closer to the production wells as a result of higher pore pressure increase.

15 GEOTHERMAL ENERGY↗

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.↗

Polarized-Cathodoluminescence Study Of Stress In GaAs/Si

Report describes experiments in which linearly polarized cathodoluminesce scanning electron microscopy, in both imaging and spectroscopic forms, used to analyze spatial variations in stress tensors in specimens of GaAs/Si. Specimens fabricated by atmospheric-pressure metallo-organic chemical vapor deposition of GaAs to thickness of 4 micrometers on Si substrates. Combination of mismatch between Si and GaAs crystal lattices and differential thermal contraction of materials from deposition temperature (700 degrees C) to room temperature gives rise to stresses in GaAs layer. Experiments conducted to contribute to understanding of effects of these stresses and ramifications thereof for fabrication of electronic devices in GaAs/Si system.

Rich, Daniel H.↗

Evaluating the Taylor hardening model in polycrystalline Ti using high energy X-ray diffraction microscopy

High energy X-ray diffraction microscopy (HEDM) was employed to index nearly 1000 grains in a polycrystalline Ti specimen and characterize their deformation during an incremental tensile test. For each grain, the positions of its associated diffraction peaks were used for analyzing its evolving crystal orientation and stress tensor. The azimuthal breadth of each peak at different load steps was measured, allowing measurement of initial yielding and providing an estimate of the dislocation density evolution in each grain. Furthermore, the Taylor hardening model was evaluated grain by grain using the above data, indicating that it can represent hardening in most grains. The yield stress and the strain hardening coefficient for all grains were statistically analyzed with respect to grain orientation, grain size, slip systems, and the surrounding neighborhood to examine how the yield stress and hardening are correlated.

36 MATERIALS SCIENCE↗

Crystallographic Texture, Structure, and Stress Transmission in Nugget Sandstone Examined With X‐Ray Tomography and Diffraction Microscopy

Subsurface processes in sandstones are controlled by porosity, permeability, and deformation mechanisms, all of which are controlled by a complex interplay of crystallographic rock texture, structure, and micromechanics. Texture, structure, and micromechanics have historically been studied using optical and electron microscopy of thin-sections. Here, we employed a new combination of in situ X-ray tomography and ray diffraction microscopy to study crystallographic texture, structure, and grain stresses in 3D. We examined these features in a sample of Nugget sandstone, a sandstone constituting hydrocarbon reservoirs across the American West. Our aims are threefold. First, we demonstrate the utility of X-ray diffraction microscopy probes for revealing texture, structure, and stress transmission in 3D. Second, we apply these techniques to Nugget sandstone and discuss findings in the context of prior work. Third, we study grain stress tensor evolution during mechanical compression to examine whether their heterogeneity and orientation evolution reflect that of inter-particle forces in granular materials. Our results show: (a) larger grains featured higher intra-granular misorientations, possibly from an increased prevalence of cements; (b) pores closed parallel to the loading direction and opened normal to loading; (c) grain stresses featured heterogeneity and orientations similar to inter-particle forces in non-cohesive granular materials; (d) grains featured compressive stresses in the loading direction and tensile stresses orthogonal to the loading direction, the latter resisting sample dilation and grain separation. Our work demonstrates the first known application of multi-modal X-ray tomography and diffraction microscopy to sandstone, providing new 3D insight into the nature of quartz cement and stress evolution.

high energy diffraction microscopy↗

A critique of some recent second-order turbulence closure models for compressible boundary layers

Computations based on two recently developed second-order turbulence closure models are compared with a series of boundary-layer experiments and with predictions of these experiments using an algebraic mixing length model. One of the models employs an eddy viscosity, whereas the other evaluates components of the Reynolds stress tensor. For flat plates, the computations are compared with the van Driest skin-friction transformation to assess the handling of compressibility. For boundary layers in pressure gradients, four experiments at Mach 4 and one at Mach 6.7 are used as the bases for comparison. In general, both models represent mean velocities and skin friction reasonably well, but represent the turbulence shear stress less accurately.

Rubesin, M. W.↗

The Acoustic Analogy: A Powerful Tool in Aeroacoustics with Emphasis on Jet Noise Prediction

The acoustic analogy introduced by Lighthill to study jet noise is now over 50 years old. In the present paper, Lighthill s Acoustic Analogy is revisited together with a brief evaluation of the state-of-the-art of the subject and an exploration of the possibility of further improvements in jet noise prediction from analytical methods, computational fluid dynamics (CFD) predictions, and measurement techniques. Experimental Particle Image Velocimetry (PIV) data is used both to evaluate turbulent statistics from Reynolds-averaged Navier-Stokes (RANS) CFD and to propose correlation models for the Lighthill stress tensor. The NASA Langley Jet3D code is used to study the effect of these models on jet noise prediction. From the analytical investigation, a retarded time correction is shown that improves, by approximately 8 dB, the over-prediction of aft-arc jet noise by Jet3D. In experimental investigation, the PIV data agree well with the CFD mean flow predictions, with room for improvement in Reynolds stress predictions. Initial modifications, suggested by the PIV data, to the form of the Jet3D correlation model showed no noticeable improvements in jet noise prediction.

Farassat, F.↗

Wave-driven torques to drive current and rotation

In the classic Landau damping initial value problem, where a planar electrostatic wave transfers energy and momentum to resonant electrons, a recoil reaction occurs in the nonresonant particles to ensure momentum conservation. To explain how net current can be driven in spite of this conservation, the literature often appeals to mechanisms that transfer this nonresonant recoil momentum to ions, which carry negligible current. However, this explanation does not allow the transport of net charge across magnetic field lines, precluding E X B rotation drive. In this study, we show that in steady state, this picture of current drive is incomplete. Using a simple Fresnel model of the plasma, we show that for lower hybrid waves, the electromagnetic energy flux (Poynting vector) and momentum flux (Maxwell stress tensor) associated with the evanescent vacuum wave become the Minkowski energy flux and momentum flux in the plasma and are ultimately transferred to resonant particles. Thus, the torque delivered to the resonant particles is ultimately supplied by the electromagnetic torque from the antenna, allowing the nonresonant recoil response to vanish and rotation to be driven. We present a warm fluid model that explains how this momentum conservation works out locally, via a Reynolds stress that does not appear in the one-dimensional initial value problem. This model is the simplest that can capture both the nonresonant recoil reaction in the initial-value problem, and the absence of a nonresonant recoil in the steady-state boundary value problem, thus forbidding rotation drive in the former while allowing it in the latter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Experimental Investigation of a Wing-Fuselage Junction Model in the NASA Langley 14- by 22-Foot Subsonic Tunnel

Current turbulence models, such as those employed in Reynolds-averaged Navier-Stokes CFD, are unable to reliably predict the onset and extent of the three-dimensional separated flow that typically occurs in wing-fuselage junctions. To critically assess, as well as to improve upon, existing turbulence models, experimental validation-quality flow-field data in the junction region is needed. In this report, we present an overview of experimental measurements on a wing-fuselage junction model that addresses this need. The experimental measurements were performed in the NASA Langley 14- by 22-Foot Subsonic Tunnel. The model was a full-span wing-fuselage body that was configured with truncated DLR-F6 wings, both with and without leading-edge extensions at the wing root. The model was tested at a fixed chord Reynolds number of 2.4 million, and angles-of-attack ranging from -10 degrees to +10 degrees were considered. Flow-field measurements were performed with a pair of miniature laser Doppler velocimetry (LDV) probes that were housed inside the model and attached to three-axis traverse systems. One LDV probe was used to measure the separated flow field in the trailing-edge junction region. The other LDV probe was alternately used to measure the flow field in the leading-edge region of the wing and to measure the incoming fuselage boundary layer well upstream of the leading edge. Both LDV probes provided measurements from which all three mean velocity components, all six independent components of the Reynolds-stress tensor, and all ten independent components of the velocity triple products were calculated. In addition to the flow-field measurements, static and dynamic pressures were measured at selected locations on the wings and fuselage of the model, infrared imaging was used to characterize boundary-layer transition, oil-flow visualization was used to visualize the separated flow in the leading- and trailing-edge regions of the wing, and unsteady shear stress was measured at limited locations using capacitive shear-stress sensors. Sample results from the measurement techniques employed during the test are presented and discussed.

Kegerise, Michael A.↗

A finite-strain rate- and pressure-dependent constitutive framework for analyzing shock compression behavior of cemented tungsten carbides to 100 GPa

In the present study a thermodynamically-consistent finite-strain rate-and-pressure-dependent constitutive framework is implemented to analyze the shock-compression behavior of cemented tungsten carbides to 100 GPa. Central to this framework is the use of logarithmic strain with a set of invariant basis that allow the Cauchy stress tensor to be expressed as a sum of three response terms that are mutually orthogonal, thus permitting a complete separation of the deviatoric and volumetric (pressure) response. An overstress viscoplasticity model that includes strain and strain rate hardening along with thermal softening is used to represent the deviatoric response, while a complete Mie-Grüneisen equation of state (EoS) is used to obtain the pressure response. Using this formulation, the shock-induced compression behavior of cemented tungsten carbide - obtained from planar plate impact experiments using a 30 mm powder gun to peak stresses of up to ~100 GPa - is analyzed to better understand the structure of the measured shock wave profiles and the associated in-material shock quantities. Of particular interest is the evolution of material inelasticity and strength, and temperature in the tungsten carbide samples during the shock compression process.

Cemented tungsten carbide↗