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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

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↗
NASA NTRS1994

Coronal heating by the resonant absorption of Alfven waves: The effect of viscous stress tensor

The time-dependent linearized magnetohydrodynamics (MHD) equations for a fully compressible, low-beta, viscoresistive plasma are solved numerically using an implicit integration scheme. The full viscosity stress tensor (Braginskii 1965) is included with the five parameters eta(sub i) i = 0 to 4. In agreement with previous studies, the numerical simulations demonstrate that the dissipation on inhomogeneities in the background Alfven speed occurs in a narrow resonant layer. For an active region in the solar corona the values of eta(sub i) are eta(sub o) = 0.65 g/cm/s, eta(sub 1) = 3.7 x 10(exp -12) g/cm/s, eta(sub 2) = 4 eta(sub 1), eta(sub 3) = 1.4 x 10(exp -6) g/cm/s, eta(sub 4) = 2 eta(sub 3), with n = 10(exp 10)/cu cm, T = 2 x 10(exp 6) K, and B = 100 G. When the Lundquist number S = 10(exp 4) and R(sub 1) much greater than S (where R(sub 1) is the dimensionless shear viscous number) the width of the resistive dissipation layer d(sub r) is 0.22a (where a is the density gradient length scale) and d(sub r) approximately S(exp -1/3). When S much greater than R(sub 1) the shear viscous dissipation layer width d(sub r) scales as R(sub 1)(exp -1/3). The shear viscous and the resistive dissipation occurs in an overlapping narrow region, and the total heating rate is independent of the value of the dissipation parameters in agreement with previous studies. Consequently, the maximum values of the perpendicular velocity and perpendicular magnetic field scale as R(sub 1)(exp -1/3). It is evident from the simulations that for solar parameters the heating due to the compressive viscosity (R(sub 0) = 560) is negligible compared to the resistive and the shear viscous (R(sub 1)) dissipation and it occurs in a broad layer of order a in width. In the solar corona with S approximately equals 10(exp 4) and R(sub 1) approximately equals 10(exp 14) (as calculated from the Braginskii expressions), the shear viscous resonant heating is of comparable magnitude to the resistive resonant heating.

Ofman, L.↗
NASA NTRS1985

Altimetry data and the elastic stress tensor of subduction zones

The stress field in the lithosphere caused by the distribution of density anomalies associated to the geoidal undulations observed by the GEOS-3 and SEASAT Earth satellites in the Tonga region was studied. Different models of the lithosphere were generated with different assumptions on the density distribution and geometry, all generating a geoid profile almost identical to the observed one. The first model is the Airy isostatic hypothesis which consists of a crust of density 2.85 laying on a lithosphere of density 3.35. The models obtained with different compensation depths give residual shortwavelength anomalies of the order of several tens of mgal and several tens of meters geoidal undulations. It indicates that there is no isostasy of the Airy type in the Tonga region because the observed geoid has very smooth undulation of about 25 m over a distance of 2000 km. The Pratt isostatic hypothesis is used in a model consisting of a crust of variable density laying on a lithosphere of higher density. This model gives smaller residual anomalies but still shows that there is no isostasy of the Pratt type in the Tonga region because the observed geoidal undulation are much smaller and smoother than the residual undulations associated to the Pratt model of isostasy.

Caputo, M.↗
NASA NTRS1987

Altimetry data and the elastic stress tensor of subduction zones

The maximum shear stress (mss) field due to mass anomalies is estimated in the Apennines, the Kermadec-Tonga Trench, and the Rio Grande Rift areas and the results for each area are compared to observed seismicity. A maximum mss of 420 bar was calculated in the Kermadec-Tonga Trench region at a depth of 28 km. Two additional zones with more than 300 bar mss were also observed in the Kermadec-Tonga Trench study. Comparison of the calculated mss field with the observed seismicity in the Kermadec-Tonga showed two zones of well correlated activity. The Rio Grande Rift results showed a maximum mss of 700 bar occurring east of the rift and at a depth of 6 km. Recorded seismicity in the region was primarily constrained to a depth of approximately 5 km, correlating well to the results of the stress calculations. Two areas of high mss are found in the Apennine region: 120 bar at a depth of 55 km, and 149 bar at the surface. Seismic events observed in the Apennine area compare favorably with the mss field calculated, exhibiting two zones of activity. The case of loading by seamounts and icecaps are also simulated. Results for this study show that the mss reaches a maximum of about 1/3 that of the applied surface stress for both cases, and is located at a depth related to the diameter of the surface mass anomaly.

Caputo, Michele↗
NASA NTRS1986

Viscosity and the Chew-Goldberger-Low equations in the solar corona

A general discussion of the dominant terms in the stress tensor in a magnetized plasma such as the solar corona is presented. The importance of dissipative terms such as electrical resistivity, heat conduction, and interspecies collisions is assessed. For average coronal conditions, the proton stress tensor is found to reduce to the dominant terms in the classical expression for the viscous stress. The classical expression can fail in the transition region, however. In the diffusion region of reconnection, classical viscosity will be appropriate if the resistivity is very large, so that the diffusion region is broad, but in that case the viscous heating is small compared to the resistive heating. On the other hand, the more general expression for the stress tensor is required if the diffusion region is thin; the stress tensor will be important in this case. The electron stress tensor is also considered, and it is shown how the classical expression for electron viscosity can fail in the transition region and lower corona.

Hollweg, J. V.↗
NASA NTRS1991

Turbulence modeling

The performance of existing two-equation eddy viscosity models was examined. An effort was made to develop better models for near-wall turbulence using direct numerical simulations of plane channel and boundary layer flows. The asymptotic near-wall behavior of turbulence was used to examine the problems of current second order closure models and develop new models with the correct near-wall behavior. Rapid Distortion Theory was used to analytically study the effects of mean deformation on turbulence, obtain analytical solutions for the spectrum tensor, Reynolds stress tensor, anisotropy tensor and its invariants, which can be used in the turbulence model development. The potential of the renormalization group theory in turbulence modeling was studied, as well as compressible turbulent flows, and modeling of bypass transition.

Shih, Tsan-Hsing↗
NASA NTRS1984

A convective model for turbulent mixing in rotating convection zones

The effects of rotation are included in an analytical model for the convective motions in a plane-parallel layer of an ideal fluid. The turbulent stress tensor, formed by taking products and averages of the various velocity components, is calculated for an arbitrary eddy size and shape. Heuristic formulae presented for determining the size and shape of the dominant eddy then give a fully specified stress tensor. Applications for this stress tensor in problems of stellar internal dynamics, heat flow, scalar diffusion, and dynamo theory are suggested. The resultant stresses tend to produce differential rotation profiles with rapidly rotating equators and interiors. The dynamo activity associated with these convective motions tends to occur near the lower boundary of the convection zone.

Hathaway, D. H.↗
NASA NTRS1996

Element Library for Three-Dimensional Stress Analysis by the Integrated Force Method

The Integrated Force Method, a recently developed method for analyzing structures, is extended in this paper to three-dimensional structural analysis. First, a general formulation is developed to generate the stress interpolation matrix in terms of complete polynomials of the required order. The formulation is based on definitions of the stress tensor components in term of stress functions. The stress functions are written as complete polynomials and substituted into expressions for stress components. Then elimination of the dependent coefficients leaves the stress components expressed as complete polynomials whose coefficients are defined as generalized independent forces. Such derived components of the stress tensor identically satisfy homogenous Navier equations of equilibrium. The resulting element matrices are invariant with respect to coordinate transformation and are free of spurious zero-energy modes. The formulation provides a rational way to calculate the exact number of independent forces necessary to arrive at an approximation of the required order for complete polynomials. The influence of reducing the number of independent forces on the accuracy of the response is also analyzed. The stress fields derived are used to develop a comprehensive finite element library for three-dimensional structural analysis by the Integrated Force Method. Both tetrahedral- and hexahedral-shaped elements capable of modeling arbitrary geometric configurations are developed. A number of examples with known analytical solutions are solved by using the developments presented herein. The results are in good agreement with the analytical solutions. The responses obtained with the Integrated Force Method are also compared with those generated by the standard displacement method. In most cases, the performance of the Integrated Force Method is better overall.

Kaljevic, Igor↗
NASA NTRS1977

Aerodynamic sound in a relaxing medium

A theory of aerodynamic sound propagation, when inhomogeneities characterized by a relaxation process are present in both the source and propagation region, is formulated. The details of the relaxation process need not be specified at the outset, although the relaxation process is characterized by a relaxation time and by an equilibrium and a frozen sound speed in a propagation region which is otherwise in equilibrium. Propagation is described in terms of a D'Alembertian characterized by the frozen sound speed relaxing toward one characterized by the equilibrium sound speed, while the source is interpreted in terms of a frozen Lighthill stress tensor relaxing toward the equilibrium stress tensor. An appropriate Green's function for the three-dimensional relaxing wave propagation operator is used to construct an exact integral for the aerodynamic sound. The sound generated far from the source is then estimated in terms of the aerodynamic sound source.

Liu, J. T. C.↗
NASA NTRS1997

Mechanical signals in plant development: a new method for single cell studies

Cell division, which is critical to plant development and morphology, requires the orchestration of hundreds of intracellular processes. In the end, however, cells must make critical decisions, based on a discrete set of mechanical signals such as stress, strain, and shear, to divide in such a way that they will survive the mechanical loads generated by turgor pressure and cell enlargement within the growing tissues. Here we report on a method whereby tobacco protoplasts swirled into a 1.5% agarose entrapment medium will survive and divide. The application of a controlled mechanical load to agarose blocks containing protoplasts orients the primary division plane of the embedded cells. Photoelastic analysis of the agarose entrapment medium can identify the lines of principal stress within the agarose, confirming the hypothesis that cells divide either parallel or perpendicular to the principal stress tensors. The coincidence between the orientation of the new division wall and the orientation of the principal stress tensors suggests that the perception of mechanical stress is a characteristic of individual plant cells. The ability of a cell to determine a shear-free orientation for a new partition wall may be related to the applied load through the deformation of the matrix material. In an isotropic matrix a uniaxial load will produce a rotationally symmetric strain field, which will define a shear-free plane. Where high stress intensities combine with the loading geometry to produce multiaxial loads there will be no axis of rotational symmetry and hence no shear free plane. This suggests that two mechanisms may be orienting the division plane, one a mechanism that works in rotationally symmetrical fields, yielding divisions perpendicular to the compressive tensor, parallel to the long axis of the cell, and one in asymmetric fields, yielding divisions parallel to the short axis of the cell and the compressive tensor.

NASA Discipline Number 40-50↗
NASA NTRS1983

Subsurface stress evaluations under rolling/sliding contacts

A computer model has been developed for evaluating the subsurface stresses incurred within rolling/sliding (elastic) contacts. The model involves first defining the stress tensor at any point (x, y, or z) beneath the surface in terms of the surface stresses. The stress tensors are analyzed to determine the maximum shear stresses and stress reversals. As a result of computations with the model, several observations were made. For example, the maximum reversing shear stresses are on the plane of the orthogonal shear stress. Further, the magnitude of these stresses is not altered by friction. However, under very high friction (typical of dry contact) surface stresses can dominate over subsurface stresses.

Kannel, J. W.↗
NASA NTRS1974

Computation of turbulent flows

All of the work deals with a model that incorporates dynamical equations for the mean velocity field, for all components of the Reynolds stress tensor, and for the 'isotropic dissipation'. The trace of the Reynolds stress tensor provides a velocity scale for the turbulence. For completeness, a turbulence length or time scale is also needed; many workers impose the length scales, but in a generalized model it must be evolved by the turbulence itself. In order, the following special classes of flows are considered: (1) decay of isotropic turbulence, (2) return to isotropy in homogeneous turbulence without strain, (3) homogeneous turbulence with strain, and (4) inhomogeneous turbulence with strain and turbulent transport.

Reynolds, W. C.↗
NASA NTRS1981

Principal stress formulation of continua motion

Principal stress coordinates are formulated for a material without tension, with the equations of motion projected onto a selected coordinate system. The translations are achieved by forming the local basis vector in correspondence with the local stress tensor. The principal stress coordinates are described as orthogonal and curvilinear, and the metric is defined. The kinematics and dynamics of the principal stress coordinates are analyzed, the latter in orthogonal coordinates, and an analytic formulation of the dynamics of a no-tension material is provided. Finally, elastic granular material such as sand, soil, etc., which are more complicated versions of no-tension material, are explored.

Zak, M.↗
NASA NTRSDate not supplied

Implicit Large-Eddy Simulations of Compressible Mixing Layers

Implicit large-eddy simulations of the self-similar regions of two compressible mixing layers at high Reynolds number and convective Mach numbers of 0.381 and 0.690 were carried out. Experimental data was used as the input into a synthetic eddy method turbulent inflow boundary condition to initiate simulations. Computational grids with nearly isotropic cells were used and the grid counts ranged from 24 to 270 million points. Agreement with experimental velocity and Reynolds stress data was good given the simplifications of the simulation. The decrease in mixing layer growth rate with increasing compressibility compares well with the literature. Budgets for the Reynolds stress transport equation were extracted from the simulations. It was demonstrated that the grids used did not resolve the dissipation term at the experimental Reynolds number. However, low Reynolds number simulations indicate that the numerical dissipation replaces the physical dissipation and the dissipation term in the transport equation can be represented by the summation term that quantifies the imbalance in the budget. The budgets for the high Reynolds number mixing layers agree well with previous low Reynolds number studies; indicating that the primary mechanism for the reduced growth rate and increased anistropy of the Reynolds stress tensor due to compressibility is a reduction in the pressure-strain and production terms. A new scaling based on the magnitude of the Reynolds stress tensor is proposed. This scaling provides better relative comparisons between the data and shows that the decrease in pressure-strain in the shear component is the primary driver in reducing the transverse normal stress.

turbulence↗
NASA NTRS1993

On the large eddy simulation of turbulent flows in complex geometry

Application of the method of Large Eddy Simulation (LES) to a turbulent flow consists of three separate steps. First, a filtering operation is performed on the Navier-Stokes equations to remove the small spatial scales. The resulting equations that describe the space time evolution of the 'large eddies' contain the subgrid-scale (sgs) stress tensor that describes the effect of the unresolved small scales on the resolved scales. The second step is the replacement of the sgs stress tensor by some expression involving the large scales - this is the problem of 'subgrid-scale modeling'. The final step is the numerical simulation of the resulting 'closed' equations for the large scale fields on a grid small enough to resolve the smallest of the large eddies, but still much larger than the fine scale structures at the Kolmogorov length. In dividing a turbulent flow field into 'large' and 'small' eddies, one presumes that a cut-off length delta can be sensibly chosen such that all fluctuations on a scale larger than delta are 'large eddies' and the remainder constitute the 'small scale' fluctuations. Typically, delta would be a length scale characterizing the smallest structures of interest in the flow. In an inhomogeneous flow, the 'sensible choice' for delta may vary significantly over the flow domain. For example, in a wall bounded turbulent flow, most statistical averages of interest vary much more rapidly with position near the wall than far away from it. Further, there are dynamically important organized structures near the wall on a scale much smaller than the boundary layer thickness. Therefore, the minimum size of eddies that need to be resolved is smaller near the wall. In general, for the LES of inhomogeneous flows, the width of the filtering kernel delta must be considered to be a function of position. If a filtering operation with a nonuniform filter width is performed on the Navier-Stokes equations, one does not in general get the standard large eddy equations. The complication is caused by the fact that a filtering operation with a nonuniform filter width in general does not commute with the operation of differentiation. This is one of the issues that we have looked at in detail as it is basic to any attempt at applying LES to complex geometry flows. Our principal findings are summarized.

Ghosal, Sandip↗
NASA NTRS2024

Analysis of faults and pit chains in Noctis Labyrinthus: Implications for early extension and possible magmatic plumbing

Noctis Labyrinthus has been a region of disputed origin due to its complexity and poor understanding of how various processes and mechanisms may have combined to form it. The surface is an integrated record of intensive tectonic activity expressed by a multiple extended sets of dip-slip faults oriented in different directions, and thought to have acted on this region over its history. These faults are always coalescent to pits and pit chains displaying a complicated geological history in the region. To understand this geological history, we mapped the surface features in Noctis Labyrinthus using the High-Resolution Stereo Camera (HRSC) onboard Mars Express ND2 nadir channel basemaps, and we adapted the Digital Terrain Map (DTM) from the Mission Experiment Gridded Data Record (MEGDR) of Mars Orbiter Laser Altimeter (MOLA) onboard Mars Global Surveyor (MGS) for the topography. We have investigated the spatial distribution and trend of fault systems, the pit chains' morphology, and the correlation between these two types of features. Our results show three fault systems: i) NS and NNE-SSW, ii) EW and ENE-WSW, and iii) NNW-SSE and NW. The analysis of the faults trending, cross-cutting correlation and the superimposition led to identify multiple intersections between these faults that have been alongside with the reactivations of some inherited faults. We interpreted the first system of fault to be related to coeval lateral extension, generated by regional stress tensor, which is probably related to the slight bending of Valles Marineris within two phases of bidirectional extension. The second system of faults has been generated by the radial oblate stress tensor related to the formation of the small shield volcanoes in Syria Planum. However, the third system is likely related to the external driving process, probably in the Tharsis province. We classified pits in four evolutionary stages based on their morphometric attributes. We believe that the formation of the pit chains in Noctis Labyrinthus is related to a surface collapse after a pressure drop related to the magma chamber deflation associated with Syria Planum volcanic province. We propose a deformational model based on early extension and magmatic plumbing as driving processes for the formation of Noctis Labyrinthus.

Mayssa El Yazidi↗