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 811 records · Page 45

Interaction between a shock wave and a turbulent boundary layer in transonic flow

Interaction between a shock wave and an unseparated turbulent boundary layer is considered. The method of matched asymptotic expansions is used, with solutions valid in the double limit as Reynolds number tends to infinity and Mach number tends to unity. The shock is weak enough that interaction effects can be considered as perturbations to the undisturbed flow; the case considered is that where the sonic line is near the outer edge of the boundary layer. It is shown that, with order estimates for Reynolds stress perturbations, the induced wall pressure distribution can be calculated using only the two outer interaction regions, independent of a specific closure condition and that this solution is in fact a turbulent free interaction solution. A detailed analysis of the inner regions, for which an eddy viscosity model for the Reynolds shear stress is used, provides a description of the variations in velocity, temperature and density near and at the wall.

Adamson, T. C., Jr.↗

A streamwise vortex embedded in a plane mixing layer

Experimental results are presented on the interaction of a single streamwise vortex with a plane turbulent mixing layer. In the present setup, the vortex is generated by a half-delta wing mounted in the settling chamber of a blower-driven wind tunnel. Initially, the vortex is shown to ride just below the mixing layer, locally distorting the mean velocity and turbulence Reynolds stress distributions in the lower part of the mixing layer. Once the vortex becomes embedded within the mixing layer, the distortions are found to spread through the whole width of the layer, and the induced extra strain rates in the mixing layer result in the generation of additional, relatively large, Reynolds normal and shear stresses, mainly in the regions above the vortex.

Bell, James H.↗

The influence of preexisting host rock damage on fault network localization

The transition from stable to unstable fracture propagation occurs when fractures begin to interact and link. Thus, fracture network coalescence controls how rocks and engineered structures fail. Here, to constrain the factors that influence localization in shear zones under brittle conditions, we build discrete element method models with a rough fault embedded in a shear zone. We add varying numbers of diffuse, randomly-placed weaknesses to examine the influence of diffuse damage on fracture network localization. The number of weaknesses controls the localization behavior of the fault network and the final fault geometry. We quantify localization using the Gini coefficient of the fracture volume, which measures the nonuniformity in a population. Each model generally increases in localization toward failure. However, models with more diffuse damage experience delocalization phases that are superimposed on the overall trend of increasing localization. The observed link between delocalization and host rock damage may help explain the varying localization of low magnitude seismicity in southern California. Models with more diffuse damage produce more complex fault geometries comprised of several parallel strands of wing cracks. The propagation of these wing cracks reduces the shear stress acting on the model boundaries, indicating that this fracture development increases the mechanical efficiency of the system.

58 GEOSCIENCES↗

Two-equation turbulence modeling of oscillatory boundary layers

Three time-dependent turbulent boundary layer problems are treated using a k-epsilon turbulence model. The flows are characterized by large pressure gradients and either reversing or rotating wall shear stress vectors. Two existing and one new near-wall formulations are tested. The computed results are compared with data from Direct Numerical Simulations of the same flows. In one case previously unpublished DNS-data are given. It is found that the k-epsilon model can provide good predictions for the wall shear as well as profiles of velocities, turbulent kinetic energy, and dissipation. The crucial part is the choice of near-wall terms.

Justesen, Peter↗

Cellular Responses to Mechanical Stress Selected Contribution: A Three-Dimensional Model for Assessment of in Vitro Toxicity in Balaena Mysticetus Renal Tissue

This study established two- and three-dimensional renal proximal tubular cell cultures of the endangered species bowhead whale (Balaena mysticetus), developed SV40-transfected cultures, and cloned the 61-amino acid open reading frame for the metallothionein protein, the primary binding site for heavy metal contamination in mammals. Microgravity research, modulations in mechanical culture conditions (modeled microgravity), and shear stress have spawned innovative approaches to understanding the dynamics of cellular interactions, gene expression, and differentiation in several cellular systems. These investigations have led to the creation of ex vivo tissue models capable of serving as physiological research analogs for three-dimensional cellular interactions. These models are enabling studies in immune function, tissue modeling for basic research, and neoplasia. Three-dimensional cellular models emulate aspects of in vivo cellular architecture and physiology and may facilitate environmental toxicological studies aimed at elucidating biological functions and responses at the cellular level. Marine mammals occupy a significant ecological niche (72% of the Earth's surface is water) in terms of the potential for information on bioaccumulation and transport of terrestrial and marine environmental toxins in high-order vertebrates. Few ex vivo models of marine mammal physiology exist in vitro to accomplish the aforementioned studies. Techniques developed in this investigation, based on previous tissue modeling successes, may serve to facilitate similar research in other marine mammals.

Goodwin, T. J.↗

A scheme for computing surface fluxes from mean flow observations

A computational method for estimating turbulent surface stress, sensible heat flux, and humidity flux from mean velocity, temperature, and humidity at a single height in the atmospheric surface layer is presented. The method is based on coupling a Monin-Obukhov similarity profile to a 'force-restore' formulation for the evolution of surface soil temperature to yield the local values of shear stress, heat flux, and surface temperature. In addition to reference-level mean flow properties, the parameters required to implement this method are the thermal heat capacity of the soil, the surface aerodynamic roughness, latitude, solar declination, surface albedo, surface emissivity, and atmospheric transmissivity. Sample calculations are presented for (1) constant atmospheric forcing at the reference level, and (2) variable atmospheric forcing corresponding to Kahle's (1977) measurements of windspeed, air temperature, and radiometer soil surface temperature under dry vegetatively sparse conditions in the Mohave Desert, Calif.

Hoffert, M. I.↗

The strength of Miranda's lithosphere

In attempting to understand the endogenic processes which have shaped the surface of an icy satellite, it is desirable to quantify the failure strength of the satellite's lithosphere. In a crust that is fractured on a large scale, frictional sliding along pre-existing fractures occurs in response to lower differential stresses than required to initiate fracture of pristine rock, thus governing failure of a brittle lithosphere. Failure is predicted along favorably oriented fracture planes; if fractures of all orientations are assumed to be present in the crust (as is expected of a heavily cratered lithosphere), frictional failure relations are directly applicable. The Coulomb criterion predicts that the shear stress (sigma sub t) and normal stress (sigma sub n) components on a fracture plane at failure are related as sigma sub t = mu-sigma sub n + S sub o, where S sub o is the cohesion and mu is the coefficient of friction. At moderate to high pressures, the frictional sliding strength of most materials is found to be sigma sub t = 0.85 sigma sub n.

Pappalardo, Robert↗

Dislocation Patterning in Deforming Crystals: Theory, Computational Predictions and Validation (Final Technical Report)

This project was awarded for an initial period of three years, followed by a no-cost extension for one year and a funded one-year renewal. During the first four years, the project focused on investigating the role of dislocation reactions and dislocation correlation in dislocation patterning in FCC metals. During the fifth year, the scope of research was expanded to investigate the effects of composition inhomogeneity on the mesoscale plastic response of Body-Centered Cubic (BCC) alloys. In addition to its impact on metal hardening during deformation, dislocation patterning provides the microstructure information required to understand phenomena like fracture and recrystallization in metals. The Continuum Dislocation Dynamics (CDD) framework was the methodology used during the first four years, followed by the use of Discrete Dislocation Dynamics (DDD) during the fifth year. CDD is a density-based formalism of dislocation dynamics in which the plastic deformation of crystals is predicted together with the mesoscale dislocation patterns by tracking the space-time evolution of dislocations as driven by the applied stress, short-range and long-range interactions of dislocations, and cross slip. CDD is a crystal mechanics approach in which the plastic constitutive part is replaced with the equations of dislocation dynamics, which is driven by the internal stress via a mobility laws, while the evolution of the density gives the eigenstrain required to update the internal stress itself. The governing equations are thus those of crystal mechanics cast as an eigenstrain problem and those of dislocation transport and reactions. Our investigations during the first four years were driven by the hypothesis that that dislocation patterning is triggered by spatiotemporal dislocation density and internal stress fluctuations, and that such fluctuations influence the collective dislocation dynamic through their effects on the short-range reaction rates and the collective dislocation mobility. This hypothesis was tested by modeling the influence of the dislocation reactions and dislocation correlations on collective dislocation dynamics. The main research components during the first four years were: 1) Reformulation of the CDD framework to integrate the dislocation correlations and dislocation reactions. 2) Simulation of the dislocation correlations and quantifying their contribution to the long-range stress of complex dislocation systems. 3) Computational implementation of the updated CDD model for FCC single crystals and development of an efficient CDD code. 4) Investigation of dislocation patterning in FCC crystals based on the updated CDD model. Some key findings from these investigations were: 1) The patterning of dislocations is initiated by cross slip and internal stress fluctuations, and the refinement of the pattern results from junction formation. 2) Patterning is more prominent under high stress. 3) The correlation stress was found to be a significant part of the mean field stress in continuum representations of dislocation dynamics. During the last year of the project, we have established a method for performing dislocation dynamics in inhomogeneous alloys and demonstrated the impact of composition inhomogeneity on the characteristics of initial yielding and dislocation character in BCC alloys. We also tested this method using an irradiated ferritic alloy in which the composition and dislocation loops effects are both present. In doing so, we have considered two aspects of the role of inhomogeneity, the creation of internal coherency stress due to the dependence of the lattice parameter on the local composition, and the dependence of the dislocation mobility on the local composition. These two aspects resulted in an unexpected behavior of the dislocation system, and, in turn, the yielding behavior of BCC alloys. Key findings include: 1) The composition inhomogeneity in BCC alloys alters the well-known role of screw dislocations by forcing a level of waviness on the average dislocation character, thus destroying the screw character dominance in the yielding observed in pure or homogeneous BCC metals. 2) In the case of irradiated alloys, while composition inhomogeneity by itself and dislocation loops by themselves result in some hardening, the superposition of the two mechanisms does not result in the sum of the two contributions via any known rule. The latter result contradicts the classical works related hardening via multiple mechanisms, which state that various hardening mechanisms can be added linearly or using some Pythagorean additive rule. We were able to rationalize this unexpected result by a closer loop at the origin of the hardening mechanisms themselves, which both involved long-range elastic interactions with dislocations. We reached the conclusion that the hardening resulting from the total stress associated with these two mechanisms (current work) differs from the sum of the effects of the individual stress fields (classical literature). We have also found a major flaw in the classical theory of spinodal hardening. The latter theory never considered the impact of composition undulation on the solute hardening. We have built a new general theory accounting for the effect of solute friction together with the composition undulations on the dislocation configuration at Critical Resolved Shear Stress (CRSS), which is yielding a new definition of the spinodal hardening and new different results. This work is the first to recognize the lack of the role of solute in spinodal hardening theory of alloys.

36 MATERIALS SCIENCE↗

Release adiabat measurements on minerals: The effect of viscosity

The current inversion of pressure-particle velocity data for release from a high pressure shock state to a pressure-density path is analyzed. It is assumed that the release process is isentropic. It was shown that for geological materials below stresses of 150 GPa, the effective viscosity must be 1000 kg/m/s in order that the viscous (irreversible) work carried out on the material in the shock state remains small compared to the mechanical work recovered upon adiabatic rarefaction. The available data pertaining to the offset of the Rayleigh line from the Hugoniot for minerals, the magnitude of the shear stress in the high pressure shock state for minerals, and the direct measurements of the viscosities of several engineering materials shocked to pressures below 150 GPa yield effective viscosities of 1000 kg/m/s or less. An inferance that this indicates that the conditions for isentropic release of minerals from shock states are achieved, and a conclusion that the application of the Riemann integral to obtain pressure-density states along the release adiabats of minerals in shock experiments is valid are made.

Jeanloz, R.↗

Slip Analysis in a Ni-base Superalloy

A Ni-base superalloy single crystal with Gamma/Gamma' structure was tested at room temperature along the , <110> and <111> directions. Consistent with previously reported investigations, this alloy did not obey the Schmid law and the CRSS (critical resolved shear stress) was noticeably lower for the <111>-oriented samples. Furthermore, the strain hardening rate decreased and the degree of deformation localization increased in the order of <111>, and <110> orientations. The appearance and orientation of deformation traces were found to depend on the loading orientation as well as the amount of strain. In general, when Gamma'-particles were sheared, the traces followed the expected octahedral shear planes. It is demonstrated that the wavy deformation traces that do not follow the {111} planes are associated with changes in the gamma-channels width and the falling off of the gamma-particles. In this paper the evolution of deformation bands are discussed in terms of plastic localization at microscopic, mesoscopic and macroscopic levels.

Westbrooke, Eboni F.↗

An uncertainty visualization framework for large-scale cardiovascular flow simulations: A case study on aortic stenosis

We present a generalizable uncertainty quantification (UQ) and visualization framework for lattice Boltzmann method simulations of high Reynolds number vascular flows, demonstrated on a patient-specific stenosed aorta. The framework combines EasyVVUQ for parameter sampling with large-eddy simulation turbulence modeling in HemeLB, and executes ensembles on the Frontier exascale supercomputer. Spatially resolved metrics, including entropy and isosurface-crossing probability, are used to map uncertainty in pressure and wall shear stress fields directly onto vascular geometries. Two sources of model variability are examined: inlet peak velocity and the Smagorinsky constant. Inlet velocity variation produces high uncertainty downstream of the stenosis where turbulence develops, while upstream regions remain stable. Smagorinsky constant variation has little effect on the large-scale pressure field but increases WSS uncertainty in localized high-shear regions. In both cases, the stenotic throat manifests low entropy, indicative of robust identification of elevated WSS. By linking quantitative UQ measures to three-dimensional anatomy, the framework improves interpretability over conventional 1D UQ plots and supports clinically relevant decision-making, with broad applicability to vascular flow problems requiring both accuracy and spatial insight.

Hemodynamics↗

Effect of adhesive debond on stress-intensity factors in bonded composite panels

Stress-intensity factors are calculated for a cracked infinite sheet adhesively bonded to a stringer, and debonding of the adhesive layer is predicted. The stringer is modeled as a semiinfinite sheet. Adhesive nonlinearity is also included. Both the sheet and stringer are treated as homogeneous, orthotropic materials, a set of integral equations is formulated and solved to obtain the adhesive shear stresses and crack-tip stress-intensity factors. Adhesive debonding is predicted using a rupture criterion based on the combined adhesive stresses. A through-the-thickness crack is located in the infinite sheet perpendicular to the edge of the stringer. When the crack is not under the stringer, the debond extends along the edge of the stringer. When the crack tip is beneath the stringer, the debond extends to the crack tip, then along the edge of the stringer. Stress levels required for debond initiation decrease as the crack tip is moved beneath the stringer. With a nonlinear adhesive, the debond initiates at higher applied stress levels than in linear adhesive cases. Compared with the linear adhesive solution, modeling a nonlinear adhesive causes the stress-intensity factor to decrease when debonding is included.

Bigelow, C. A.↗

Particle Image Velocimetry Measurements of a Two/Three-dimensional Separating/Reattaching Boundary Layer Downstream of an Axisymmetric Backward-facing Step

Planar Particle Image Velocimetry measurements were obtained in the separating/reattaching flow region downstream of an axisymmetric backward-facing step. Data were acquired for a two-dimensional (2D) separating boundary layer at five different Reynolds numbers based on step height (Re(sub h)), spanning 5900-33000, and for a three-dimensional (3D) separating boundary layer at Re(sub h) = 5980 and 8081. Reynolds number effects were investigated in the 2D cases using mean-velocity field, streamwise and wall-normal turbulent velocity, and Reynolds stress statistics. Results show that both the reattachment length (x(sub r)) and the secondary separation point are Reynolds number dependent. The reattachment length increased with rising Re(sub h) while the secondary recirculation region decreased in size. These and other Re(sub h) effects were interpreted in terms of changes in the separating boundary layer thickness and wall-shear stress. On the other hand, in the 3D case, it was found that the imposed cross-flow component was relatively weak in comparison to the streamwise component. As a result, the primary influences of three dimensionality only affected the near-separation region rather than the entire separation bubble.

Hudy, Laura M.↗

Effects of edge-seal design on the mechanical and thermal performance of vacuum-insulated glazing

Although vacuum-insulated glazing (VIG) has been proposed as a promising solution towards developing energy-efficient buildings, VIGs have not become popular in the market due to several technical challenges including the complexity of the fabrication process. In particular, the edge-seal is a key component that significantly affects the thermal insulation and mechanical performance, and the development of edge-seal with adequate thermal insulation, mechanical strength, and reasonable processing cost is essential to overcome such technical issues in VIG. For this purpose, effects of edge-seal design parameters on the VIG performance should be identified. In this research, we analyzed the edge-seal for thermal transport as well as structural stresses to study the effects, and then identified and evaluated the material mixes for the edge-seal requirements. Here, the finite element simulations showed the significance of VIG corner calculation on overall thermal transmittance and the importance of seal conductivity below 1 W/m.K. The experiments with the flexible seals with different ratios of fine glass powder demonstrated that the measured shear strength values for the seal with less than 30% glass powder were more than 10 times larger than the calculated shear stress values. Based on these simulation and experimental results, a flexible sealant was developed using a proprietary mix of ceramic materials that meets the requirements of the designed VIG edge-seal, including structural as well as thermal stress resistance and a low conductivity. Moreover, the sealant is self-curing under atmospheric conditions, and thus it does not require costly inline process of laser curing or oven baking.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

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

Effect of debond growth on stress-intensity factors in a cracked orthotropic sheet stiffened by a semi-infinite orthotropic sheet

Stress-intensity factors are determined for a cracked infinite sheet adhesively bonded to a stringer, and debonding of the adhesive layer is predicted. The stringer is modeled as a semi-infinite sheet. Adhesive nonlinearity is also included. Both the sheet and stringer are treated as homogeneous, orthotropic materials. A set of integral equations is formulated and solved to obtain the adhesive shear stresses and crack-tip stress-intensity factors. Adhesive debonding is predicted using a rupture criterion based on the combined adhesive stresses. When the crack is not under the stringer, the debond extends along the edge of the stringer. When the crack tip is beneath the stringer, the debond grows to the end of the crack, then along the edge of the stringer. Stress levels required for debond initiation decrease as the crack tip is moved beneath the stringer. With a nonlinear adhesive, the debond initiates at higher applied stress levels than in linear adhesive cases. Compared with the linear adhesive solution, modeling a nonlinear adhesive causes the stress-intensity factor to increase when the bond is assumed to remain intact but causes the stress-intensity factor to decrease when debonding is included.

Bigelow, C. A.↗

A micromechanics-based strength prediction methodology for notched metal matrix composites

An analytical micromechanics based strength prediction methodology was developed to predict failure of notched metal matrix composites. The stress-strain behavior and notched strength of two metal matrix composites, boron/aluminum (B/Al) and silicon-carbide/titanium (SCS-6/Ti-15-3), were predicted. The prediction methodology combines analytical techniques ranging from a three dimensional finite element analysis of a notched specimen to a micromechanical model of a single fiber. In the B/Al laminates, a fiber failure criteria based on the axial and shear stress in the fiber accurately predicted laminate failure for a variety of layups and notch-length to specimen-width ratios with both circular holes and sharp notches when matrix plasticity was included in the analysis. For the SCS-6/Ti-15-3 laminates, a fiber failure based on the axial stress in the fiber correlated well with experimental results for static and post fatigue residual strengths when fiber matrix debonding and matrix cracking were included in the analysis. The micromechanics based strength prediction methodology offers a direct approach to strength prediction by modeling behavior and damage on a constituent level, thus, explicitly including matrix nonlinearity, fiber matrix debonding, and matrix cracking.

Bigelow, C. A.↗

A micromechanics-based strength prediction methodology for notched metal-matrix composites

An analytical micromechanics-based strength prediction methodology was developed to predict failure of notched metal matrix composites. The stress-strain behavior and notched strength of two metal matrix composites, boron/aluminum (B/Al) and silicon-carbide/titanium (SCS-6/Ti-15-3), were predicted. The prediction methodology combines analytical techniques ranging from a three-dimensional finite element analysis of a notched specimen to a micromechanical model of a single fiber. In the B/Al laminates, a fiber failure criteria based on the axial and shear stress in the fiber accurately predicted laminate failure for a variety of layups and notch-length to specimen-width ratios with both circular holes and sharp notches when matrix plasticity was included in the analysis. For the SCS-6/Ti-15-3 laminates, a fiber failure based on the axial stress in the fiber correlated well with experimental results for static and postfatigue residual strengths when fiber matrix debonding and matrix cracking were included in the analysis. The micromechanics-based strength prediction methodology offers a direct approach to strength prediction by modeling behavior and damage on a constituent level, thus, explicitly including matrix nonlinearity, fiber matrix debonding, and matrix cracking.

Bigelow, C. A.↗