Predicting shear transmission across grain boundaries with an iterative stress relief model
Not provided.
SEARCH · Engineering Papers
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.
Not provided.
Conditionally sampled, ensemble-averaged velocity measurements, made with a laser velocimeter, were taken in the flowfield over the rear half of an 18% thick circular arc airfoil at zero incidence tested at M = 0.76 and at a Reynolds number based on chord of 11 x 10(exp 6). Data for one cycle of periodic unsteady flow having a reduced frequency f of 0.49 are analyzed. A series of compression waves, which develop in the early stages of the cycle, strengthen and coalesce into a strong shock wave that moves toward the airfoil leading edge. A thick shear layer forms downstream of the shock wave. The kinetic energy and shear stresses increase dramatically, reach a maximum when dissipation and diffusion of the turbulence exceed production, and then decrease substantially. The response lime of the turbulence to the changes brought about by the shock-wave passage upstream depends on the shock-wave strength and position in the boundary layer. The cycle completes itself when the shock wave passes the midchord, weakens, and the shear layer collapses. Remarkably good comparisons are found with computations that employ the time-dependent Reynolds averaged form of the Navier-Stokes equations using an algebraic eddy viscosity model, developed for steady flows.
Conditionally sampled, ensemble-averaged velocity measurements, made with a laser velocimeter, were taken in the flowfield over the rear half of an 18% thick circular arc airfoil at zero incidence tested at M = 0.76 and of a Reynolds number based on chord of 11 x 10(exp 6). Data for one cycle of periodic unsteady flow having a reduced frequency bar-f of 0.49 are analyzed. A series of compression waves, which develop in the early stages of the cycle, strengthen and coalesce into a strong shock wave that moves toward the airfoil leading edge. A thick shear layer forms downstream of the shock wave. The kinetic energy and shear stresses increase dramatically, reach a maximum when dissipation and diffusion of the turbulence exceed production, and then decrease substantially. The response time of the turbulence to the changes brought about by the shock-wave passage upstream depends on the shock-wave strength and position in the boundary layer. The cycle completes itself when the shock wave passes the midchord, weakens, and the shear layer collapses. Remarkably good comparisons are found with computations that employ the time-dependent Reynolds averaged form of the Navier-Stokes equations using an algebraic eddy viscosity model, developed for steady flows.
Our team has developed a novel approach to measuring creep at extremely high temperatures using electrostatic levitation (ESL). This method has been demonstrated on niobium up to 2300 C, while ESL has melted tungsten (3400 C). This method has been extended to lower temperatures and higher stresses and applied to new materials, including a niobium-based superalloy, MASC. High-precision machined spheres of the sample are levitated in the NASA MSFC ESL, a national user facility and heated with a laser. The samples are rotated with an induction motor at up to 30,000 revolutions per second. The rapid rotation loads the sample through centripetal acceleration, producing a shear stress of about 60 MPa at the center, causing the sample to deform. The deformation of the sample is captured on high-speed video, which is analyzed by machine-vision software from the University of Massachusetts. The deformations are compared to finite element models to determine the constitutive constants in the creep relation. Furthermore, the non-contact method exploits stress gradients within the sample to determine the stress exponent in a single test.
The unstructured tetrahedral grid cell-centered finite volume flow solver USM3D has been recently extended to handle mixed element grids composed of hexahedral, prismatic, pyramidal, and tetrahedral cells. Presently, two turbulence models, namely, baseline Spalart-Allmaras (SA) and Menter Shear Stress Transport (SST), support mixed element grids. This paper provides an overview of the various numerical discretization options available in the newly enhanced USM3D. Using the SA model, the flow solver extensions are verified on three two-dimensional test cases available on the Turbulence Modeling Resource website at the NASA Langley Research Center. The test cases are zero pressure gradient flat plate, planar shear, and bump-inchannel. The effect of cell topologies on the flow solution is also investigated using the planar shear case. Finally, the assessment of various cell and face gradient options is performed on the zero pressure gradient flat plate case.
The laser interferometric skin friction meter was used to measure wall shear stress distributions in two interactions of fin-generated swept shock waves with turbulent boundary layers. The basic research configuration was an unswept sharp-leading-edge fin of variable angle mounted on a flatplate. The results indicate that such measurements are practical in high-speed interacting flows, and that a repeatability of + or - 6 percent or better is possible. Marked increases in wall shear were observed in both swept interactions tested.
Laser Doppler velocimetry (LDV) measurements and hydrogen bubble flow visualization techniques were used to examine the near-wall flow structure of 2D and 3D turbulent boundary layers (TBLs) over a range of low Reynolds numbers. The goals of this research were (1) an increased understanding of the flow physics in the near wall region of turbulent boundary layers,(2) to observe and quantify differences between 2D and 3D TBL flow structures, and (3) to document Reynolds number effects for 3D TBLs. The LDV data have provided results detailing the turbulence structure of the 2D and 3D TBLs. These results include mean Reynolds stress distributions, flow skewing results, and U and V spectra. Effects of Reynolds number for the 3D flow were also examined. Comparison to results with the same 3D flow geometry but at a significantly higher Reynolds number provided unique insight into the structure of 3D TBLs. While the 3D mean and fluctuating velocities were found to be highly dependent on Reynolds number, a previously defined shear stress parameter was discovered to be invariant with Reynolds number. The hydrogen bubble technique was used as a flow visualization tool to examine the near-wall flow structure of 2D and 3D TBLs. Both the quantitative and qualitative results displayed larger turbulent fluctuations with more highly concentrated vorticity regions for the 2D flow.
Turbulence models, with increasing complexity, up to triple product terms, are applied to the flow in a rotating pipe. The rotating pipe is a challenging case for turbulence models as it contains significant rotational and curvature effects. The flow field starts with the classic fully developed pipe flow, with a stationary pipe wall. This well defined condition is then subjected to a section of pipe with a rotating wall. The rotating wall introduces a second velocity scale, and creates Reynolds shear stresses in the radial-circumferential and circumferential-axial planes. Furthermore, the wall rotation introduces a flow stabilization, and actually reduces the turbulent kinetic energy as the flow moves along the rotating wall section. It is shown in the present work that the Reynolds stress models are capable of predicting significant reduction in the turbulent kinetic energy, but triple product improves the predictions of the centerline turbulent kinetic energy, which is governed by convection, dissipation and transport terms, as the production terms vanish on the pipe axis.
Turbulence models, with increasing complexity, up to triple product terms, are applied to the flow in a rotating pipe. The rotating pipe is a challenging case for turbulence models as it contains significant rotational and curvature effects. The flowfield starts with the classic fully developed pipe flow, with a stationary pipe wall. This well defined condition is then subjected to a section of pipe with a rotating wall. The rotating wall introduces a second velocity scale, and creates Reynolds-shear stresses in the radial-circumferential and circumferential-axial planes. Furthermore, the wall rotation introduces a ow stabilization, and actually reduces the turbulent kinetic energy as the flow moves along the rotating wall section. It is shown in the present work that the Reynolds stress models are capable of predicting significant reduction in the turbulent kinetic energy, but triple product improves the predictions of the centerline turbulent kinetic energy, which is governed by convection, dissipation and transport terms, as the production terms vanish on the pipe axis.
A fatigue failure criterion applicable to composite materials is developed and applied to predict the fatigue behavior of graphite/epoxy laminates with particular emphasis on the influence of temperature. Tensile stress-strain curves and tension-tension fatigue curves for various unidirectional, angle-ply and symmetrically balanced laminates were developed at test temperatures of 25 C, 74 C, and 114 C. For most laminates a reduction in both static strength and fatigue strength is observed with increasing temperature. This reduction appeared more severe in fatigue loading than in static tensile loading and most severe where the shear stress in the lamina is the dominate failure mode. Through an analytical formulation of shifting functions for the influences of temperature, all fatigue data are shown to be capable of being reduced to a single reference curve at some temperature. Examples are given which demonstrate the capability of the fatigue failure criterion to predict failure of complex symmetrically balanced laminates from relevant parameters obtained from the observed behavior of unidirectional and angle-ply laminates.
The eddy viscosity hypothesis is a popular method in wind turbine wake modeling for estimating turbulent Reynolds stresses. We document the downstream evolution of eddy viscosity in the wake of a wind turbine from experimental and large-eddy-simulation data. Wake eddy viscosity is isolated from its surroundings by subtracting the inflow profile, and the driving forces are identified in each wake region. Eddy viscosity varies in response to changes in turbine geometry and nacelle misalignment with larger turbines generating stronger velocity gradients and shear stresses. We propose a model for eddy viscosity based on a Rayleigh distribution. Model parameters are obtained from scaling the eddy viscosity hypothesis and demonstrate satisfactory agreement with the reference data. The model is implemented in the curled wake formulation in the FLOw Redirection and Induction in Steady State (FLORIS) framework and assessed through comparisons with the previous formulation. Our approach produced more accurate flow field estimates with lower total error for the majority of cases.
We analyze the cycle of great earthquakes along the San Andreas fault with a finite element numerical model of deformation in a crust with a nonlinear viscoelastic rheology. The viscous component of deformation has an effective viscosity that depends exponentially on the inverse absolute temperature and nonlinearity on the shear stress; the elastic deformation is linear. Crustal thickness and temperature are constrained by seismic and heat flow data for California. The models are for anti plane strain in a 25-km-thick crustal layer having a very long, vertical strike-slip fault; the crustal block extends 250 km to either side of the fault. During the earthquake cycle that lasts 160 years, a constant plate velocity v(sub p)/2 = 17.5 mm yr is applied to the base of the crust and to the vertical end of the crustal block 250 km away from the fault. The upper half of the fault is locked during the interseismic period, while its lower half slips at the constant plate velocity. The locked part of the fault is moved abruptly 2.8 m every 160 years to simulate great earthquakes. The results are sensitive to crustal rheology. Models with quartzite-like rheology display profound transient stages in the velocity, displacement, and stress fields. The predicted transient zone extends about 3-4 times the crustal thickness on each side of the fault, significantly wider than the zone of deformation in elastic models. Models with diabase-like rheology behave similarly to elastic models and exhibit no transient stages. The model predictions are compared with geodetic observations of fault-parallel velocities in northern and central California and local rates of shear strain along the San Andreas fault. The observations are best fit by models which are 10-100 times less viscous than a quartzite-like rheology. Since the lower crust in California is composed of intermediate to mafic rocks, the present result suggests that the in situ viscosity of the crustal rock is orders of magnitude less the rock viscosity determined in the laboratory.
High Cycle Fatigue (HCF) induced failures in aircraft gas-turbine and rocket engine turbopump blades is a pervasive problem. Single crystal turbine blades are being utilized in rocket engine turbopumps and jet engines throughout industry and NASA because of their superior creep, stress rupture, melt resistance and thermomechanical fatigue capabilities over polycrystalline alloys. Single-crystal materials have highly orthotropic properties making the position of the crystal lattice relative to the pan geometry a significant factor in the overall analysis. The failure modes of single crystal turbine blades is complicated to predict due to the material orthotropy and variations in crystal orientations. A fatigue failure criteria based on the maximum shear stress amplitude [delta t max] on the 30 slip systems, is presented for single crystal nickel superalloys (FCC crystal). This criteria reduces the scatter in uniaxial LCF test data, for four different specimen orientations, for PWA 1484 at 1200 F in air, quite well. A power law curve fit of the failure parameter, delta t max, vs. cycles to failure is presented.
A fatigue failure criterion applicable to composite materials is developed and applied to predict the fatigue behavior of graphite/epoxy laminates with particular emphasis on the influence of temperature. Tensile stress-strain curves and tension-tension fatigue curves for various unidirectional, angle-ply and symmetrically balanced laminates were developed at test temperatures of 25 C, 74 C and 114 C. In general for most laminates a reduction in both static strength and fatigue strength is observed with increasing temperature. This reduction appeared more severe in fatigue loading than in static tensile loading and most severe where the shear stress in the lamina is the dominant failure mode. Through an analytical formulation of shifting functions for the influences of temperature, all fatigue data are shown to be capable of being reduced to a single reference curve at some temperature. Additionally, examples are given which demonstrate the capability of the fatigue failure criterion to predict failure of complex symmetrically balanced laminates from relevant parameters obtained from the observed behavior of unidirectional and angle-ply laminates.
The existence of strike-slip faults was recently documented in two locations on Mars. Two clear examples are reviewed located southeast of Valles Marineris and preliminary evidence is presented for more widespread strike-slip deformation elsewhere in Coprates. The first two examples show that strike-slip faulting occurred in a broad zone east of the Coprates Rise spanning approximately 400 km east-west by perhaps 1000 km north-south. The last example suggests that the growth of major wrinkle ridges throughout Coprates may have been influenced by horizontally directed shear stresses and that more than one generation of ridges was produced. Thus, 'compressional' deformation of ridged plains south of Valles Marineris was spatially heterogeneous and a temporal change in stress may have been involved.
For an infinitely long plate strip of constant curvature there is set up the buckling condition for combined shear and axial stresses. Assumptions are made of the vanishing of the moments acting on the longitudinal edges and of the vanishing of the tangential displacements u and the transverse stresses sigma(sub y) (longitudinal stiffeners nonresistant to bending in the circumferential direction). The buckling condition is evaluated for the case of pure shear and an approximate formula is derived for the critical stresses for large curvatures. In connection with the discussion of the methods and results there is defined the range of validity and the buckling condition set up.
Hydrogen production techniques based on solar-water splitting have emerged as carbon-free energy systems. Many researchers have developed highly efficient thin-film photoelectrochemical (PEC) devices made of low-cost and earth-abundant materials. However, solar water splitting systems suffer from short lifetimes due to catalyst instability that is attributed to both chemical dissolution and mechanical stress produced by hydrogen bubbles. A recent study found that the nanoporous hydrogel could prevent the structural degradation of the PEC devices. In this study, we investigate the protection mechanism of the hydrogel-based overlayer by engineering its porous structure using the cryogelation technique. Tests for cryogel overlayers with varied pore structures, such as disconnected micropores, interconnected micropores, and surface macropores, reveal that the hydrogen gas trapped in the cryogel protector reduce shear stress at the catalyst surface by providing bubble nucleation sites. The cryogelated overlayer effectively preserves the uniformly distributed platinum catalyst particles on the device surface for over 200 h. Our finding can help establish semi-permanent photoelectrochemical devices to realize a carbon-free society.
Several conditional sampling techniques are applied to a data base generated by large-eddy simulation of turbulent channel flow. It is shown that the bursting process is associated with well-organized horseshoe vortices inclined at about 45 deg. to the wall. These vortical structures are identified by examining the vortex lines of three-dimensional, ensemble averaged vorticity fields. Two distinct horseshoe-shaped vortices corresponding to the sweep and ejection events are detected. These vortices are associated with high Reynolds shear stress and hence make a significant contribution to turbulent energy production. The dependency of the ensemble averaged vortical structures on the detection criteria, and the question of whether this ensemble-averaged structure is an artifact of the ensemble averaging process are examined. The ensemble-averaged pattern of these vortical structures that emerge from the analysis could provide the basis for a hypothetical model of the organized structures of wall-bounded shear flows.