SEARCH · Search NASA

Results for “strain 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Evolution of the shock front and turbulence structures in the shock/turbulence interaction

The interaction of a weak shock front with isotropic turbulence has been investigated using Direct Numerical Simulation (DNS). Two problems were considered: the ability of the field equation (the equation for a propagating surface) to model the shock; and a quantitative study of the evolution of turbulence structure using the database generated by Lee et al. Field equation model predictions for front shape have been compared with DNS results; good agreement is found for shock wave interaction with 2D turbulence and for a single steady vorticity wave. In the interaction of 3D isotropic turbulence with a normal shock, strong alignment of vorticity with the intermediate eigenvector of the rate of strain tensor (S(sup *)(sub ij) = S(sub ij) - (1/3)(delta(sub ij))(S(sub kk))) is seen to develop upstream of the shock and to be further amplified on passage through the shock. Vorticity tends to align at 90 deg to the largest eigenvector, but there is no preferred alignment with the smallest eigenvector. Upstream of the shock, the alignments continue to develop even after the velocity derivative skewness saturates. There is a significant tendency, which increases with time throughout the computational domain, for velocity to align with vorticity. The alignment between velocity and vorticity is strongest in eddy regions and weakest in convergence regions.

Kevlahan, N.↗

NASA NTRS1993

Wrinkle ridges, reverse faulting, and the depth penetration of lithospheric stress in lunae planum, Mars

Tectonic features on a planetary surface are commonly used as constraints on models to determine the state of stress at the time the features formed. Quantitative global stress models applied to understand the formation of the Tharsis province on Mars constrained by observed tectonics have calculated stresses at the surface of a thin elastic shell and have neglected the role of vertical structure in influencing the predicted pattern of surface deformation. Wrinkle ridges in the Lunae Planum region of Mars form a conentric pattern of regularly spaced features in the eastern and southeastern part of Tharsis; they are formed due to compressional stresses related to the response of the Martian lithosphere to the Tharsis bulge. As observed in the exposures of valley walls in areas such as the Kasei Valles, the surface plains unit is underlain by an unconsolidated impact-generated megaregolith that grades with depth into structurally competent lithospheric basement. The ridges have alternatively been hypothesized to reflect deformation restricted to the surface plains unit ('thin skinned deformation') and deformation that includes the surface unit, megaregolith and basement lithosphere ('thick skinned deformation'). We have adopted a finite element approach to quantify the nature of deformation associated with the development of wrinkle ridges in a vertically stratified elastic lithosphere. We used the program TECTON, which contains a slippery node capability that allowed us to explicitly take into account the presence of reverse faults believed to be associated with the ridges. In this study we focused on the strain field in the vicinity of a single ridge when slip occurs along the fault. We considered two initial model geometries. In the first, the reverse fault was assumed to be in the surface plains unit, and in the second the initial fault was located in lithospheric basement, immediately beneath the weak megaregolith. We are interested in the conditions underwhich strain in the surface layer and basement either penetrates or fails to penetrate through the megaregolith. We thus address the conditions required for an initial basement fault to propagate through the megaregolith to the surface, as well as the effect of the megareolith on the strain tensor in the vicinity of a fault that nucleates in the surface plains unit.

Zuber, M. T.↗

NASA NTRS1993

The terrestrial gravitational wave environment from known sources

The objective of this project was to produce a gravitational wave spectral line list of all known binary stars producing expected strain amplitudes at Earth in excess of h = 10 (exp -21), or gravitational wave fluxes in excess of F = 10 (exp -12) erg cm(exp -2) s(exp -1). These strain and flux limits lie above the anticipated detection thresholds for space-borne laser interferometers capable of detecting gravitational radiation in the 10 micron Hz to 1 Hz frequency range. The source list was intended to provide frequency (including each harmonic), amplitude and phase (for each polarization and harmonic), and celestial coordinates for each system, lacking only the orientation of the principal polarization axis with respect to the pole of the coordinate system, and the sign of the source phase and frequency (or, equivalently, of the sense of rotation of the strain tensor with time) from providing a complete source description. Such a spectral line list would lay essential groundwork for high-sensitivity, low-frequency searches for gravitational radiation.

Webbink, Ronald F.↗

NASA NTRSDate not supplied

The Influence of Adaptive Mesh Refinement on the Prediction of Vortex Interactions about a Generic Missile Airframe

The complex interaction of forebody and wing vortices significantly impacts missile aerodynamics. The formation of these vortices involves smooth regions of the geometry or geometric discontinuities like leading edges, trailing edges, tips, and corners. Regions of supersonic flow and complex shock topologies interact with boundary layers and vortices. Smooth-body separation and 3D viscous effects strain current Reynolds-averaged Navier-Stokes (RANS) techniques. The quantification and control of discretization error is critical to obtaining reliable simulation results and often turbulence model assessments are made in the presence of unquantified (and potentially large) discretization errors. Two mesh adaptation schemes are applied to steady RANS simulations. Multiscale unstructured mesh adaptation is applied to control interpolation error estimates of the Mach field, which resolves boundary layers, vortices, and shocks. A dual-mesh approach with overset communication is applied between an expert-crafted near-body unstructured mesh and an adaptive off-body Cartesian mesh refined with Q-criterion scaled by the strain tensor magnitude. A generic missile configuration is examined in a supersonic flow field to show the interaction of mesh adaptation and turbulence model. Turbulence model modifications for rotational correction and a quadratic constitutive relationship show a strong influence on adaptive mesh refinement and predicted rolling moment.

slender body aerodynamics↗

NASA NTRS1997

Time Evolution of Modeled Reynolds Stresses in Planar Homogeneous Flows

The analytic expression of the time evolution of the Reynolds stress anisotropy tensor in all planar homogeneous flows is obtained by exact integration of the modeled differential Reynolds stress equations. The procedure is based on results of tensor representation theory, is applicable for general pressure-strain correlation tensors, and can account for any additional turbulence anisotropy effects included in the closure. An explicit solution of the resulting system of scalar ordinary differential equations is obtained for the case of a linear pressure-strain correlation tensor. The properties of this solution are discussed, and the dynamic behavior of the Reynolds stresses is studied, including limit cycles and sensitivity to initial anisotropies.

Jongen, T.↗

NASA NTRS1990

A turbulent two-phase flow model for nebula flows

A new and very efficient turbulent two-phase flow numericaly model is described to analyze the environment of a protoplanetary nebula at a stage prior to the formation of planets. Focus is on settling processes of dust particles in flattened gaseous nebulae. The model employs a perturbation technique to improve the accuracy of the numerical simulations of such flows where small variations of physical quantities occur over large distance ranges. The particles are allowed to be diffused by gas turbulence in addition to settling under gravity. Their diffusion coefficients is related to the gas turbulent viscosity by the non-dimensional Schmidt number. The gas turbulent viscosity is determined by the means of the eddy viscosity hypothesis that assumes the Reynolds stress tensor proportional to the mean strain rate tensor. Zero- and two-equation turbulence models are employed. Modeling assumptions are detailed and discussed. The numerical model is shown to reproduce an existing analytical solution for the settling process of particles in an inviscid nebula. Results of nebula flows are presented taking into account turbulence effects of nebula flows. Diffusion processes are found to control the settling of particles.

Champney, Joelle M.↗

NASA NTRS1973

A finite element tensor approach to plate buckling and postbuckling.

A practical finite element method is presented for geometrically nonlinear (von Karman) plate problems. Symmetric strain energy tensors provide an efficient formulation and solution, and allow the effects of initial imperfections to be considered. An unbalanced force iteration is employed, with the advantage that the time consuming Gauss decomposition needs to be performed only once. The convergence is controlled by an automatically computed relaxation factor, and this makes the method applicable in advanced stages of non-linearity. Several different triangular elements are used in explicitly deriving the required tensors.

Vos, R. G.↗

NASA NTRS1981

Numerical simulation of a compressible homogeneous, turbulent shear flow

A direct, low Reynolds number, numerical simulation was performed on a homogeneous turbulent shear flow. The full compressible Navier-Stokes equations were used in a simulation on the ILLIAC IV computer with a 64,000 mesh. The flow fields generated by the code are used as an experimental data base, to examine the behavior of the Reynols stresses in this simple, compressible flow. The variation of the structure of the stresses and their dynamic equations as the character of the flow changed is emphasized. The structure of the tress tensor is more heavily dependent on the shear number and less on the fluctuating Mach number. The pressure-strain correlation tensor in the dynamic uations is directly calculated in this simulation. These correlations are decomposed into several parts, as contrasted with the traditional incompressible decomposition into two parts. The performance of existing models for the conventional terms is examined, and a model is proposed for the 'mean fluctuating' part.

Feiereisen, W. J.↗

NASA NTRS2017

Effects of Adiabatic Heating on the High Strain Rate Deformation of Polymer Matrix Composites

Polymer matrix composites (PMCs) are increasingly being used in aerospace structures that are expected to experience complex dynamic loading conditions throughout their lifetime. As such, a detailed understanding of the high strain rate behavior of the constituents, particularly the strain rate, temperature, and pressure dependent polymer matrix, is paramount. In this paper, preliminary efforts in modeling experimentally observed temperature rises due to plastic deformation in PMCs subjected to dynamic loading are presented. To this end, an existing isothermal viscoplastic polymer constitutive formulation is extended to model adiabatic conditions by incorporating temperature dependent elastic properties and modifying the components of the inelastic strain rate tensor to explicitly depend on temperature. It is demonstrated that the modified polymer constitutive model is capable of capturing strain rate and temperature dependent yield as well as thermal softening associated with the conversion of plastic work to heat at high rates of strain. The modified constitutive model is then embedded within a strength of materials based micromechanics framework to investigate the manifestation of matrix thermal softening, due to the conversion of plastic work to heat, on the high strain rate response of a T700Epon 862 (T700E862) unidirectional composite. Adiabatic model predictions for high strain rate composite longitudinal tensile, transverse tensile, and in-plane shear loading are presented. Results show a substantial deviation from isothermal conditions; significant thermal softening is observed for matrix dominated deformation modes (transverse tension and in-plane shear), highlighting the importance of accounting for the conversion of plastic work to heat in the polymer matrix in the high strain rate analysis of PMC structures.

polymer matrix composites↗

NASA NTRS2003

Implementation of an Associative Flow Rule Including Hydrostatic Stress Effects Into the High Strain Rate Deformation Analysis of Polymer Matrix Composites

A previously developed analytical formulation has been modified in order to more accurately account for the effects of hydrostatic stresses on the nonlinear, strain rate dependent deformation of polymer matrix composites. State variable constitutive equations originally developed for metals have been modified in order to model the nonlinear, strain rate dependent deformation of polymeric materials. To account for the effects of hydrostatic stresses, which are significant in polymers, the classical J2 plasticity theory definitions of effective stress and effective inelastic strain, along with the equations used to compute the components of the inelastic strain rate tensor, are appropriately modified. To verify the revised formulation, the shear and tensile deformation of two representative polymers are computed across a wide range of strain rates. Results computed using the developed constitutive equations correlate well with experimental data. The polymer constitutive equations are implemented within a strength of materials based micromechanics method to predict the nonlinear, strain rate dependent deformation of polymer matrix composites. The composite mechanics are verified by analyzing the deformation of a representative polymer matrix composite for several fiber orientation angles across a variety of strain rates. The computed values compare well to experimentally obtained results.

Goldberg, Robert K.↗

NASA NTRS1986

On fracture of fibrous composites

The problem of an existing through-the-thickness center notch in a unidirectional fibrous composite is investigated using a homogeneous anisotropic elasticity analysis in conjunction with various crack growth criteria. Tensor polynomial, strain energy density, and normal stress ratio criteria are applied as crack growth direction criteria and compared to experimental observations. For the cases and criteria studied, only the normal stress ratio criterion exhibits an ability to correctly predict crack growth direction. The normal stress ratio criterion is subsequently applied to the same cases as a failure criterion. Comparison with experiment indicates that these predictions exhibit reasonable correlation with observed critical stresses. Finally, the macroscopic approach to the problem used in this study is compared to a microscopic shear lag analysis of experimental results in the literature. The comparison indicates that both approaches yield essentially identical predictions. The comparisons made with other criteria and approaches indicate that application of the normal stress ratio within a macroscopic level analysis has great potential as an approach to modeling composite crack growth.

Beuth, J. L., Jr.↗

NASA NTRS1987

Direct numerical simulation of buoyantly driven turbulence

Numerical simulations of homogeneous turbulence subject to buoyant forcing were performed. The presence of a mean temperature gradient combined with a gravitational field results in a forcing term in the momentum equations. The development of the turbulence was studied and compared to the decay of similar fields in the absence of gravity. In the buoyantly driven field, the vorticity is preferentially aligned with the intermediate eigenvector of the strain-rate tensor and the local temperature gradient is more likely to be aligned with the most compressive eigenvector. These relationships are qualitatively similar to those observed in previous shear flow results studied by Ashurst (1987). A tensor diffusivity model for passive scalar transport developed from shear flow results in Rogers, Moin, and Reynolds (1986) also predicts this buoyant scalar transport, indicating that the relationship between the scalar flux and the Reynolds stress is similar in both flows.

Ashurst, William T.↗

NASA NTRS1990

A contribution toward rational modeling of the pressure-strain-rate correlation

A novel method of obtaining an analytical expression of the 'linear part' of the pressure-strain-rate tensor in terms of the anisotropy tensor of the Reynolds stresses has been developed, where the coefficients of the seven independent tensor terms are functions of the invariants of the Reynolds-stress anisotropy. The coefficients are evaluated up to fourth order in the anisotropy of the Reynolds stresses to provide guidance for development of a turbulence model.

Lee, Moon Joo↗

NASA NTRS1992

Pancakelike domes on Venus

The shape of seven large domes on the plains of Venus, with volumes between 100 and 1000 cu km, is compared with that of an axisymmetric gravity current spreading over a rigid horizontal surface. Both the altimetric profiles and the horizontal projection of the line of intersection of domes on the SAR images agree well with the theoretical similarity solution for a newtonian fluid, but not with the shape calculated for a rigid-plastic rheology, nor with that for a static model with a strong skin. As a viscous current spreads, it generates an isotropic strain rate tensor whose magnitude is independent of radius. Such a flow can account for the randomly oriented cracks that are uniformly distributed on the surface of the domes. The stress induced by the flow in the plains material below is obtained, and is probably large enough to produce the short radial cracks in the surface of the plains beyond the domes. The viscosity of the domes can be estimated from their thermal time constants if spreading is possible only when the fluid is hot, and lies between 10(exp 14) and 10(exp 17) Pa s. Laboratory experiments show that such viscosities correspond to temperatures of 610 - 690 C in dry rhyolitic magmas. These temperatures agree with laboratory measurements of the solidus temperature of wet rhyolite. These results show that the development of the domes can be understood using simple fluid dynamical ideas, and that the magmas involved can be produced by wet melting at depths below 10 km, followed by eruption and degassing.

Mckenzie, Dan↗

NASA NTRS1995

Global Dynamic Modeling of Space-Geodetic Data

The proposal had outlined a year for program conversion, a year for testing and debugging, and two years for numerical experiments. We kept to that schedule. In first (partial) year, author designed a finite element for isostatic thin-shell deformation on a sphere, derived all of its algebraic and stiffness properties, and embedded it in a new finite element code which derives its basic solution strategy (and some critical subroutines) from earlier flat-Earth codes. Also designed and programmed a new fault element to represent faults along plate boundaries. Wrote a preliminary version of a spherical graphics program for the display of output. Tested this new code for accuracy on individual model plates. Made estimates of the computer-time/cost efficiency of the code for whole-earth grids, which were reasonable. Finally, converted an interactive graphical grid-designer program from Cartesian to spherical geometry to permit the beginning of serious modeling. For reasons of cost efficiency, models are isostatic, and do not consider the local effects of unsupported loads or bending stresses. The requirements are: (1) ability to represent rigid rotation on a sphere; (2) ability to represent a spatially uniform strain-rate tensor in the limit of small elements; and (3) continuity of velocity across all element boundaries. Author designed a 3-node triangle shell element which has two different sets of basis functions to represent (vector) velocity and all other (scalar) variables. Such elements can be shown to converge to the formulas for plane triangles in the limit of small size, but can also applied to cover any area smaller than a hemisphere. The difficult volume integrals involved in computing the stiffness of such elements are performed numerically using 7 Gauss integration points on the surface of the sphere, beneath each of which a vertical integral is performed using about 100 points.

Bird, Peter↗

NASA NTRS1999

Scale-Similar Models for Large-Eddy Simulations

Scale-similar models employ multiple filtering operations to identify the smallest resolved scales, which have been shown to be the most active in the interaction with the unresolved subgrid scales. They do not assume that the principal axes of the strain-rate tensor are aligned with those of the subgrid-scale stress (SGS) tensor, and allow the explicit calculation of the SGS energy. They can provide backscatter in a numerically stable and physically realistic manner, and predict SGS stresses in regions that are well correlated with the locations where large Reynolds stress occurs. In this paper, eddy viscosity and mixed models, which include an eddy-viscosity part as well as a scale-similar contribution, are applied to the simulation of two flows, a high Reynolds number plane channel flow, and a three-dimensional, nonequilibrium flow. The results show that simulations without models or with the Smagorinsky model are unable to predict nonequilibrium effects. Dynamic models provide an improvement of the results: the adjustment of the coefficient results in more accurate prediction of the perturbation from equilibrium. The Lagrangian-ensemble approach [Meneveau et al., J. Fluid Mech. 319, 353 (1996)] is found to be very beneficial. Models that included a scale-similar term and a dissipative one, as well as the Lagrangian ensemble averaging, gave results in the best agreement with the direct simulation and experimental data.

Sarghini, F.↗

NASA NTRS2000

Turbulence Model Predictions of Strongly Curved Flow in a U-Duct

The ability of three types of turbulence models to accurately predict the effects of curvature on the flow in a U-duct is studied. An explicit algebraic stress model performs slightly better than one- or two-equation linear eddy viscosity models, although it is necessary to fully account for the variation of the production-to-dissipation-rate ratio in the algebraic stress model formulation. In their original formulations, none of these turbulence models fully captures the suppressed turbulence near the convex wall, whereas a full Reynolds stress model does. Some of the underlying assumptions used in the development of algebraic stress models are investigated and compared with the computed flowfield from the full Reynolds stress model. Through this analysis, the assumption of Reynolds stress anisotropy equilibrium used in the algebraic stress model formulation is found to be incorrect in regions of strong curvature. By the accounting for the local variation of the principal axes of the strain rate tensor, the explicit algebraic stress model correctly predicts the suppressed turbulence in the outer part of the boundary layer near the convex wall.

Rumsey, Christopher L.↗

NASA NTRS1997

Damage Initiation and Ultimate Tensile Strength of Scaled [0 deg n/90 deg n/0 deg n]sub T Graphite-Epoxy Coupons

Previous research on scaling effects in composite materials has demonstrated that the stress levels at first ply failure and ultimate failure of composite laminates are dependent on the size of the laminate. In particular, the thickness dimension has been shown to be the most influential parameter in strength scaling of composite coupons loaded in tension. Geometrically and constitutively scaled laminates exhibit decreasing strength with increasing specimen size, and the magnitude of the strength-size effect is a function of both material properties and laminate stacking sequence. Some of the commonly used failure criteria for composite materials such as maximum stress, maximum strain, and tensor polynomial (e.g., Tsai-Wu) cannot account for the strength-size effect. In this paper, three concepts are developed and evaluated for incorporating size dependency into failure criteria for composite materials. An experimental program of limited scope was performed to determine the first ply failure stress in scaled cross-ply laminates loaded in tension. Test specimens were fabricated of AS-4/3502 graphite-epoxy composite material with laminate stacking sequences of [0 deg n/90 deg n/o deg n]subT where n=1-6. Two experimental techniques were used to determine first ply failure, defined as a transverse matrix crack in the 90 deg ply: (1) step loading with dye penetrant x-ray of the specimen at each load interval, and (2) acoustic emission. The best correlation between first ply failure analysis and experimental data was obtained using a modified Weibull approach which incorporated the residual thermal stress and the outer ply constraint, as well as the ply thickness effect. Finally, a second set of experiments was performed to determine the tensile response and ultimate failure of the scaled cross-ply laminates. The results of these experiments indicated no influence of specimen size on tensile response or ultimate strength.

Jackson, Karen E.↗