Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shear stress”

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 307 records · Page 17

The elliptical crack subjected to nonuniform shear loading.

The equations of elasticity are solved for the problem of a flat elliptical crack which has nonuniform shear stresses applied to its surfaces. The shear stresses are prescribed independently in two directions on the crack surface and are expressed in the form of a third-degree polynomial. Mode-two and mode-three stress-intensity factors are presented in analytical and graphical form as functions of position along the crack border.

Smith, F. W.↗

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

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

Menter, Florian R.↗

Shear-Sensitive Liquid Crystal Coating Method: Surface-Inclination Effects on Shear Vector Measurements

The shear-sensitive liquid crystal coating (SSLCC) method is an image-based technique for both visualizing dynamic surface-flow phenomena, such as transition and separation, and for measuring the continuous shear-stress vector distribution acting on an aerodynamic surface. Under proper lighting and viewing conditions (discussed below), the coating changes color in response to an applied aerodynamic shear. This color-change response is continuous and reversible, with a response time of milliseconds, and is a function of both the shear magnitude and the shear vector orientation relative to the observer. The liquid crystal phase of matter is a weakly-ordered, viscous, non-Newtonian fluid state that exists between the nonuniform liquid phase and the ordered solid phase of certain organic compounds. Cholesteric liquid crystal compounds possess a helical molecular arrangement that selectively scatters white light, incident along the helical axis, as a three-dimensional spectrum. This property is linked to the helical pitch length, which is within the range of wavelengths in the visible spectrum. The pitch length, and hence the wavelength of the scattered light, is influenced by shear stress normal to the helical axis. This unique optical property produces a measurable color change in response to an applied shearing force. The full-surface shear stress vector measurement method, developed at NASA-Ames, is schematically illustrated. As with the visualization method, the coated test surface is illuminated from the normal direction with white light and the camera is positioned at an above-plane view angle of approximately 30 deg. Experiments have been initiated at NASA Ames to begin the process of quantifying surface-inclination (surface-curvature) effects on shear vector measurement accuracy. In preliminary experiments, surface-inclination angles theta(sub x), theta(sub y) of 0, +/-5, +/-10, and +/-15 deg were employed. In this arrangement, white-light illumination was positioned normal to the untilted test surface, and the camera above-plane view angle was set at 30 deg relative to the untilted test surface. As can be seen, vector-aligned lambda(sub d) values showed no dependence on theta(sub x) or theta(sub y) for absolute values of these tilt angles is less than or equal to 15 deg. Acquisition and analyses of full-surface color images are presently underway to definitively document the insensitivity limits of the shear vector measurement methodology to surface-slope variations.

Reda, Daniel C.↗

Zonal Two Equation Kappa-Omega Turbulence Models for Aerodynamic Flows

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

Menter, Florian R.↗

Surface temperatures and glassy state investigations in tribology

Measurements were made of the limiting shear stress for two naphthenic oils of differing molecular weight and three blends of the lower molecular weight oil and polymers of differing molecular weight. All reached the same limiting shear stress for the same temperature and pressure; although the polymer solutions reduced the limiting shear stress by about fifteen percent. A falling body viscometer was constructed to operate to 230 C and to 0.6 GPa and another was constructed to extend the pressure range to 1.1 GPa. A concentrated contact simulator was developed which allows recording of the traction force while the slide-roll ratio is continuously varied and the rolling speed is maintained essentially constant. Measurement of lubricant minimum film thickness of elliptical EHD contacts of various aspect ratios were made by optical interferometry. The experimental data were thirty percent greater than that predicted by the Hamrock and Dowson model. Preliminary development of the application of a scanning infrared radiation system to a tribological system was completed.

Bair, S.↗

Application of fiber bridging models to fatigue crack growth in unidirectional titanium matrix composites

Several fiber bridging models were reviewed and applied to study the matrix fatigue crack growth behavior in center notched (0)(sub 8) SCS-6/Ti-15-3 and (0)(sub 4) SCS-6/Ti-6Al-4V laminates. Observations revealed that fatigue damage consisted primarily of matrix cracks and fiber matrix interfacial failure in the (0)(sub 8) SCS-6/Ti-15-3 laminates. Fiber-matrix interface failure included fracture of the brittle reaction zone and cracking between the two carbon rich fiber coatings. Intact fibers in the wake of the matrix cracks reduce the stress intensity factor range. Thus, an applied stress intensity factor range is inappropriate to characterize matrix crack growth behavior. Fiber bridging models were used to determine the matrix stress intensity factor range in titanium metal matrix composites. In these models, the fibers in the wake of the crack are idealized as a closure pressure. An unknown constant frictional shear stress is assumed to act along the debond or slip length of the bridging fibers. The frictional shear stress was used as a curve fitting parameter to available data (crack growth data, crack opening displacement data, and debond length data). Large variations in the frictional shear stress required to fit the experimental data indicate that the fiber bridging models in their present form lack predictive capabilities. However, these models provide an efficient and relatively simple engineering method for conducting parametric studies of the matrix growth behavior based on constituent properties.

Bakuckas, J. G., Jr.↗

Turbulence stress measurements in a nonadiabatic hypersonic boundary layer

Turbulent shear stress and direct turbulent total heat-flux measurements have been made across a nonadiabatic, zero pressure gradient, hypersonic boundary layer by using specially designed hot-wire probes free of strain-gauging and wire oscillation. Heat-flux measurements were in reasonably good agreement with values obtained by integrating the energy equation using measured profiles of velocity and temperature. The shear-stress values deduced from the measurements, by assuming zero correlation of velocity and pressure fluctuations, were lower than the values obtained by integrating the momentum equation. Statistical properties of the cross-correlations are similar to corresponding incompressible measurements at approximately the same momentum-thickness Reynolds number.

Mikulla, V.↗

Tensile and compressive plastic deformation behavior of medium-entropy Cr-Co-Ni single crystals from cryogenic to elevated temperatures

The equiatomic Cr-Co-Ni medium-entropy alloy has the face-centered cubic (FCC) structure. Bulk single crystals of this alloy were grown and tested in tension and compression between 14 K and 1373 K with the loading axis parallel to [23]. At room temperature, the critical resolved shear stress (CRSS) for {111}<10> slip is 65 ± 5 MPa and does not exhibit a tension-compression asymmetry. It does, however, increase significantly as the test temperature decreases. A dulling of this temperature dependence occurs below 50 K, which may be due to the inertia effect. When the measured values above 50 K are extrapolated to lower temperatures, a value of 225 MPa is estimated for the CRSS at 0 K. This is larger than that (168 MPa) previously determined for the equiatomic Cr-Mn-Fe-Co-Ni high-entropy alloy using a similar procedure. The stacking fault energy of the present Cr-Co-Ni is estimated to be about 14 mJm-2, which is sufficiently low to account for deformation twinning both at 77 K and room temperature. Twinning at 77 K occurs on conjugate (1) planes at an onset shear stress of 482 MPa after primary slip and propagates in the form of Lüders deformation. At room temperature, twinning occurs uniformly throughout the gauge section on primary (111) planes at an onset shear stress of 381 MPa after primary and subsequent conjugate slip. Thin layers with the hexagonal close-packed stacking are observed in association with twinning both at 77 K and room temperature.

36 MATERIALS SCIENCE↗

Modeling the role of local crystallographic correlations in microstructures of Ti-6Al-4V using a correlated structure visco-plastic self-consistent polycrystal plasticity formulation

This paper presents a multi-level crystal plasticity-based simulation framework for modeling mechanical response and microstructure evolution of Ti-6Al-4V with α-lath/lamellar microstructures. Here, the model is a correlated structure visco-plastic self-consistent (CS-VPSC) formulation linking three scales: a single crystals micro-scale, a lath/lamellar colony meso-scale, and a lath/lamellar aggregate macro-scale. A selected hardening law for the evolution of critical resolved shear stress per slip system used in CS-VPSC is phenomenological. However, it adjusts the resistances of basal and prismatic slip systems based on the geometry of slip transfer between adjacent lamellae. Consistent with experimental evidences, the resolved shear stress on the pyramidal slip planes is dependent not only on the stress in the direction of slip but also on the two orthogonal shear stress components and the three normal stress components (non-Schmid effects). Electron backscatter diffraction (EBSD) data in conjunction with a procedure relying on α→β phase transformation is used to construct paired variants of α-laths/lamellae satisfying their local crystallographic correlations. The procedure fits volume fractions of individual laths/lamellae with the experimental EBSD data and selects a distribution of habit planes between adjacent variants with respect to the loading direction. The simulation framework is applied to interpret the deformation behavior in tension and compression along two sample directions of Ti-6Al-4V fabricated via laser powder bed fusion. Moreover, the model is used to simulate texture evolution during rolling of the material to large strains. It is demonstrated that the model is capable of predicting plastic anisotropy/asymmetry and the concomitant texture evolution. While the model reveals a significant effect of habit plane inclination with respect to the loading direction on yield stress, the comparison of the data and model predictions shows that a random distribution of habit planes fits the flow response. It is further inferred that the tension-compression asymmetry arises from the non-Schmid effects.

36 MATERIALS SCIENCE↗

Computer Modeling of Thermal Convection in Melts to Explain Glass Formation in Low Gravity and on Earth

Experiments conducted up to this time on glass forming melts in the low gravity environment of space show that glasses prepared in low-g are more chemically homogeneous and more resistant to crystallization than the comparable glasses prepared at 1-g on Earth. This result is somewhat surprising and opposite to the accepted concept on glass formation for a melt. A hypothesis based on "shear thinning" of a melt, a decrease in viscosity with increasing shear stress, is proposed as an explanation for the observed low-gravity results. This paper describes detailed simulation procedures to test the role of thermal convection in introducing shear stress in glass forming melts, using a lithium disilcate melt as a model. The simulation system in its idealized version consists of a cylinder that is heated at one end and cooled at the other with gravity acting in a transverse direction to the thermal gradient. The side wall of the cylinder is assumed to be insulating. The governing equations of motion and energy are solved using variable properties for viscosity (Arrehenius and non-Arrehenius behaviors) and density (constant and temperature dependent). Other parametric variables in the calculations include gravity level and gravity vector orientation. The shear stress in the system are then computed as a function of gravity from the calculated values of maximum melt velocity, and its effect on melt viscosity (shear thinning) is predicted. Also included and discussed are the modeling efforts related to other potential convective processes in glass forming melts and their possible effects on melt viscosity.

Ray, Chandra S.↗

In search of ice-stream sticky spots

The form drag of large bedrock bumps sticking into the base of an ice stream can produce effective 'sticky spots' supporting large basal shear stress. Bedrock regions surrounded by lubricating till at the same topographic level can cause sticky spots, but tend to collect lubricating water and thus are unlikely to support a shear stress of more than a few tenths of a bar unless they contain abundant large bumps. Raised regions on the ice-air surface also can cause moderate increases in the shear stress supported on the bed beneath. Surveys of large-scale bedrock roughness, strain grids across the margins of ice-surface highs, and possibly, water-pressure measurements in regions of thin or zero till would help identify and characterize sticky spots.

Alley, Richard B.↗

Models of recurrent strike-slip earthquake cycles and the state of crustal stress

Numerical models of the strike-slip earthquake cycle, assuming a viscoelastic asthenosphere coupling model, are examined. The time-dependent simulations incorporate a stress-driven fault, which leads to tectonic stress fields and earthquake recurrence histories that are mutually consistent. Single-fault simulations with constant far-field plate motion lead to a nearly periodic earthquake cycle and a distinctive spatial distribution of crustal shear stress. The predicted stress distribution includes a local minimum in stress at depths less than typical seismogenic depths. The width of this stress 'trough' depends on the magnitude of crustal stress relative to asthenospheric drag stresses. The models further predict a local near-fault stress maximum at greater depths, sustained by the cyclic transfer of strain from the elastic crust to the ductile asthenosphere. Models incorporating both low-stress and high-stress fault strength assumptions are examined, under Newtonian and non-Newtonian rheology assumptions. Model results suggest a preference for low-stress (a shear stress level of about 10 MPa) fault models, in agreement with previous estimates based on heat flow measurements and other stress indicators.

Lyzenga, Gregory A.↗

Surface temperatures and glassy state investigations in tribology, part 4

Measurements were made of the limiting shear stress for two naphthenic oils of differing molecular weight and three blends of the lower molecular weight oil and polyalkylmethacrylate polymers of differing molecular weight. The two base oils reached the same limiting shear stress for the same temperature and pressure. This was also true for all the polymer solutions although the polymer reduced the limiting shear stress by about 15 percent. It is shown that limiting stress is more a function of material type than viscosity or molecular weight. A new falling body viscometer was constructed to operate to 230 C and 0.6 GPa. Another viscometer was constructed to extend the pressure range to 1.1 GPa. A concentrated contact simulator was developed which allows recording of the traction force while the slide-roll ratio is continuously varied and the rolling speed is maintained essentially constant by a single drive motor. The configuration is that of a crowned roller against a disk. Measurement of lubricant minimum film thickness of elliptical EHD contacts of various aspect ratios were made by optical interferometry. The data collected were used to evaluate the Hamrock and Dowson minimum film thickness model over a range of contract ellipticity ratio where the major axis of the contact ellipse was aligned both parallel and perpendicular to the direction of motion. A statistical analysis of the measured film thickness data showed that on the average the experimental data were 30 percent greater than the film thickness predicted by the model. Preliminary development of the application of a scanning infrared radiation system to a tribo-system was completed.

Bair, S. S.↗

Direct numerical simulation of curved turbulent channel flow

Low Reynolds number, mildly curved, turbulent channel flow has been simulated numerically without subgrid scale models. A new spectral numerical method developed for this problem was used, and the computations were performed with 2 million degrees of freedom. A variety of statistical and structural information has been extracted from the computed flow fields. These include mean velocity, turbulence stresses, velocity skewness, and flatness factors, space time correlations and spectra, all the terms in the Reynolds stress balance equations, and contour and vector plots of instantaneous velocity fields. The effects of curvature on this flow were determined by comparing the concave and convex sides of the channel. The observed effects are consistent with experimental observations for mild curvature. The most significant difference in the turbulence statistics between the concave and convex sides was in the Reynolds shear stress. This was accompanied by significant differences in the terms of the Reynolds shear stress balance equations. In addition, it was found that stationary Taylor-Gortler vortices were present and that they had a significant effect on the flow by contributing to the mean Reynolds shear stress, and by affecting the underlying turbulence.

Moser, R. D.↗

Analysis of two-equation turbulence models for recirculating flows

The two-equation kappa-epsilon model is used to analyze turbulent separated flow past a backward-facing step. It is shown that if the model constraints are modified to be consistent with the accepted energy decay rate for isotropic turbulence, the dominant features of the flow field, namely the size of the separation bubble and the streamwise component of the mean velocity, can be accurately predicted. In addition, except in the vicinity of the step, very good predictions for the turbulent shear stress, the wall pressure, and the wall shear stress are obtained. The model is also shown to provide good predictions for the turbulence intensity in the region downstream of the reattachment point. Estimated long time growth rates for the turbulent kinetic energy and dissipation rate of homogeneous shear flow are utilized to develop an optimal set of constants for the two equation kappa-epsilon model. The physical implications of the model performance are also discussed.

Thangam, S.↗

Postbuckling failure of composite plates with central holes

The postbuckling failure of square composite plates with central holes is analyzed numerically and experimentally. The particular plates studies have stacking sequences of: (+ and - 45/0/90)(sub 2S); (+ and - 45/0(sub 2))(sub 2S); (+ and - 45/0(sub 6))(sub S); and (+ and - 45)(sub 4S). A simple plate geometry, one with a hole diameter to plate width ratio of 0.3 is compared. Failure load, failure mode, and failure location are predicted numerically by using the finite element method. Predictions are compared with experimental results. In numerical failure analysis the interlaminar shear stresses, as well as the inplane stresses are taken into account. An issue addressed in this study is the possible mode shape change of the plate during loading. It is predicted that the first three laminates fail due to excessive stresses in the fiber direction, and more importantly, that the load level is independent of whether the laminate is deformed in a one-half or two-half wave configuration. It is predicted that the fourth laminate fails due to excessive inplane shear stress. Interlaminar shear failure is not predicted for any laminates. For the first two laminates the experimental observations correlated well with the predictions. Experimentally, the third laminate failed along the side support due to interlaminar shear strength S(sub 23). The fourth experimental laminate failed due to inplane shear in the location predicted, however material softening resulted in a different failure load from predictions.

Lee, H. H.↗

Calculation of free turbulent mixing by interaction approach.

The applicability of Bradshaw's interaction hypothesis to two-dimensional free shear flows was investigated. According to it, flows with velocity extrema may be considered to consist of several interacting layers. The hypothesis leads to a new expression for the shear stress which removes the usual restriction that shear stress vanishes at the velocity extremum. The approach is based on kinetic energy and the length scale equations. The compressible flow equations are simplified by restriction to low Mach numbers, and the range of their applicability is discussed. The empirical functions of the turbulence model are found here to be correlated with the spreading rate of the shear layer. The analysis demonstrates that the interaction hypothesis is a workable concept.

Morel, T.↗

Scaling patch analysis of planar turbulent wakes

In this work, a scaling patch approach is used to investigate the proper scales in planar turbulent wakes. A proper scale for the mean axial flow is the well-known maximum velocity deficit U ref = U ∞ –U ctr , where U ∞ is the free stream velocity and U ctr is the mean axial velocity at the wake centerline. From an admissible scaling of the mean continuity equation, a proper scale for the mean transverse flow is found as V ref = (dδ/dx)U ref , where dδ/dx is the growth rate of the wake width. From an admissible scaling of the mean momentum equation, a proper scale for the kinematic Reynolds shear stress is found as R uv,ref = U ∞ V ref , which is a mixed scale of the free stream velocity and the mean transverse flow scale. Expressions are derived for the scaled mean transverse velocity and Reynolds shear stress in the far field of planar turbulent wakes. Using a Gaussian function for the mean axial velocity deficit, approximate functions for the scaled mean transverse velocity and Reynolds shear stress are developed and found to agree well with experimental and simulation data. This work reveals that the mean transverse flow, despite its small magnitude, plays an important role in the scaling and understanding of the planar turbulent wake.

42 ENGINEERING↗