Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shear Layers”

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 1,549 records · Page 86

Comparison of AGE and Spectral Methods for the Simulation of Far-Wakes

Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This report demonstrates that a good compromise is possible with the Advected Grid Explicit (AGE) method. AGE has proven to be both efficient and accurate for simulating turbulent free-shear flows, including planar mixing layers and planar jets. Its efficiency results from its localized fully explicit finite difference formulation (Bisset 1998a,b) that is very straightforward to compute, outweighing the need for a fairly small timestep. Also, most of the successful simulations were slightly under-resolved, and therefore they were, in effect, large-eddy simulations (LES) without a sub-grid-scale (SGS) model, rather than direct numerical simulations (DNS). The principle is that the role of the smallest scales of turbulent motion (when the Reynolds number is not too low) is to dissipate turbulent energy, and therefore they do not have to be simulated when the numerical method is inherently dissipative at its resolution limits. Such simulations are termed 'auto-LES' (LES with automatic SGS modeling) in this report.

Bisset, D. K.↗

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

More accurate predictions with transonic Navier-Stokes methods through improved turbulence modeling

Significant improvements in predictive accuracies for off-design conditions are achievable through better turbulence modeling; and, without necessarily adding any significant complication to the numerics. One well established fact about turbulence is it is slow to respond to changes in the mean strain field. With the 'equilibrium' algebraic turbulence models no attempt is made to model this characteristic and as a consequence these turbulence models exaggerate the turbulent boundary layer's ability to produce turbulent Reynolds shear stresses in regions of adverse pressure gradient. As a consequence, too little momentum loss within the boundary layer is predicted in the region of the shock wave and along the aft part of the airfoil where the surface pressure undergoes further increases. Recently, a 'nonequilibrium' algebraic turbulence model was formulated which attempts to capture this important characteristic of turbulence. This 'nonequilibrium' algebraic model employs an ordinary differential equation to model the slow response of the turbulence to changes in local flow conditions. In its original form, there was some question as to whether this 'nonequilibrium' model performed as well as the 'equilibrium' models for weak interaction cases. However, this turbulence model has since been further improved wherein it now appears that this turbulence model performs at least as well as the 'equilibrium' models for weak interaction cases and for strong interaction cases represents a very significant improvement. The performance of this turbulence model relative to popular 'equilibrium' models is illustrated for three airfoil test cases of the 1987 AIAA Viscous Transonic Airfoil Workshop, Reno, Nevada. A form of this 'nonequilibrium' turbulence model is currently being applied to wing flows for which similar improvements in predictive accuracy are being realized.

Johnson, Dennis A.↗

Near-wall model for compressible turbulent boundary layers based on an inverse velocity transformation

In this work, a near-wall model, which couples the inverse of a recently developed compressible velocity transformation (Griffin et al., Proc. Natl Acad. Sci., vol. 118, 2021, p. 34) and an algebraic temperature–velocity relation, is developed for high-speed turbulent boundary layers. As input, the model requires the mean flow state at one wall-normal height in the inner layer of the boundary layer and at the boundary-layer edge. As output, the model can predict mean temperature and velocity profiles across the entire inner layer, as well as the wall shear stress and heat flux. The model is tested in an a priori sense using a wide database of direct numerical simulation high-Mach-number turbulent channel flows, pipe flows and boundary layers (48 cases, with edge Mach numbers in the range 0.77–11, and semi-local friction Reynolds numbers in the range 170–5700). The present model is significantly more accurate than the classical ordinary differential equation (ODE) model for all cases tested. The model is deployed as a wall model for large-eddy simulations in channel flows with bulk Mach numbers in the range 0.7–4 and friction Reynolds numbers in the range 320–1800. When compared to the classical framework, in the a posteriori sense, the present method greatly improves the predicted heat flux, wall stress, and temperature and velocity profiles, especially in cases with strong heat transfer. In addition, the present model solves one ODE instead of two, and has a computational cost and implementation complexity similar to that of the commonly used ODE model.

42 ENGINEERING↗

Lubricating and wear mechanisms for a hemisphere sliding on polyimide-bonded graphite fluoride film

Friction, wear life, rider wear, and film wear for a 440 C high-temperature-stainless-steel, hemispherically tipped rider sliding against polyimide-bonded graphite fluoride films were evaluated in a moist-air atmosphere at 25 C. Optical microscope and surface profilometry observations were made at various sliding intervals to determine how film thickness affected the lubricating and failure mechanisms of the films. Two lubrication regimes operated for the same load. In the first, the film supported the load and the lubricating mechanism consisted of the shear (plastic flow) of a thin layer of the lubricant between the metallic rider and the film surface. In the second, the film did not support the load (it was worn away) and the lubricating mechanism consisted of the shear of very thin lubricant films between flat areas generated on the rider and on sandblasted metallic asperities in the film wear track. Lubricant was supplied from the valleys between the asperities or from the sides of the wear track. With thicker films, wear life increased since a greater lubricant supply was available from the sides of the wear track.

Fusaro, R. L.↗

Nonlinear adhesive behavior effects in a cracked orthotropic sheet stiffened by a semi-infinite orthotropic sheet

The stress-intensity factors are determined for a cracked orthotropic sheet adhesively bonded to an orthotropic stringer where the adhesive layer is modeled with a nonlinear stress-strain curve. By the use of Green's functions and the complex variable theory of orthotropic elasticity, a set of integral equations is obtained. The integral equations are replaced by an equivalent set of algebraic equations, which are solved to obtain the shear stress distribution in the adhesive layer, with which the crack-tip stress-intensity factors are found. When the adhesive was modeled with a nonlinear stress-strain curve, the peak shear stresses in the adhesive were considerably reduced in comparison to the solution for the linear elastic adhesive. This resulted in increases in the stress-intensity factors for the nonlinear adhesive solution compared to the linear adhesive solution. The nonlinear adhesive has no significant effect on the stress-intensity factor unless the near crack tip is beneath the stringer. It is assumed that the adhesive bond remains intact and it is predicted that onset of adhesive failure occurs at decreasing levels of applied stress as the crack propagates beneath the stringer.

Bigelow, C. A.↗

Nonlinear adhesive behavior effects in a cracked orthotropic sheet stiffened by a semi-infinite orthotropic sheet

The stress-intensity factors are determined for a cracked orthotropic sheet adhesively bonded to an orthotropic stringer where the adhesive layer is modeled with a nonlinear stress-strain curve. By the use of Green's functions and the complex variable theory of orthotropic elasticity, a set of integral equations is obtained. The integral equations are replaced by an equivalent set of algebraic equations, which are solved to obtain the shear stress distribution in the adhesive layer, with which the crack-tip stress-intensity factors are found. When the adhesive was modeled with a nonlinear stress-strain curve, the peak shear stresses in the adhesive were considerably reduced in comparison to the solution for the linear elastic adhesive. This resulted in increases in the stress-intensity factors for the nonlinear adhesive solution compared to the linear adhesive solution. The nonlinear adhesive has no significant effect on the stress-intensity factor unless the near crack tip is beneath the stringer. It is assumed that the adhesive bond remains intact and it is predicted that onset of adhesive failure occurs at decreasing levels of applied stress as the crack propagates beneath the stringer.

Bigelow, C. A.↗

Wind shear modeling for aircraft hazard definition

Mathematical models of wind profiles were developed for use in fast time and manned flight simulation studies aimed at defining and eliminating these wind shear hazards. A set of wind profiles and associated wind shear characteristics for stable and neutral boundary layers, thunderstorms, and frontal winds potentially encounterable by aircraft in the terminal area are given. Engineering models of wind shear for direct hazard analysis are presented in mathematical formulae, graphs, tables, and computer lookup routines. The wind profile data utilized to establish the models are described as to location, how obtained, time of observation and number of data points up to 500 m. Recommendations, engineering interpretations and guidelines for use of the data are given and the range of applicability of the wind shear models is described.

Frost, W.↗

Electric Field Observations of Plasma Convection, Shear, Alfven Waves, and other Phenomena Observed on Sounding Rockets in the Cusp and Boundary Layer

On December 14,2002, a NASA Black Brant X sounding rocket was launched equatorward from Ny Alesund, Spitzbergen (79 N) into the dayside cusp and subsequently cut across the open/closed field line boundary, reaching an apogee of771 km. The launch occurred during Bz negative conditions with strong By negative that was changing during the flight. SuperDarn (CUTLASS) radar and subsequent model patterns reveal a strong westward/poleward convection, indicating that the rocket traversed a rotational reversal in the afternoon merging cell. The payload returned DC electric and magnetic fields, plasma waves, energetic particle, suprathermal electron and ion, and thermal plasma data. We provide an overview of the main observations and focus on the DC electric field results, comparing the measured E x B plasma drifts in detail with the CUTLASS radar observations of plasma drifts gathered simultaneously in the same volume. The in situ DC electric fields reveal steady poleward flows within the cusp with strong shears at the interface of the closed/open field lines and within the boundary layer. We use the observations to discuss ionospheric signatures of the open/closed character of the cusp/low latitude boundary layer as a function of the IMF. The electric field and plasma density data also reveal the presence of very strong plasma irregularities with a large range of scales (10 m to 10 km) that exist within the open field line cusp region yet disappear when the payload was equatorward of the cusp on closed field lines. These intense low frequency wave observations are consistent with strong scintillations observed on the ground at Ny Alesund during the flight. We present detailed wave characteristics and discuss them in terms of Alfven waves and static irregularities that pervade the cusp region at all altitudes.

Pfaff, R. F.↗

Kink band and shear band localization in anisotropic perfectly plastic solids

Shear-driven strain localization has been observed in a wide variety of materials and may take the form of shear bands or kink bands. Based on observations of kink bands in plastically anisotropic metallic nanolaminates and single crystal metals, here we posit that, for the specific case of isochoric deformation, the kinematics of kink band formation are indistinguishable from those of plane strain shear band formation. The only distinction between shear bands and kink bands in these systems would then be that kink bands ‘lock up’ at a particular value of material rotation while shear bands may progress to arbitrarily high strains. In order to investigate whether strong material anisotropy is sufficient to arrest shear localization at a geometry that matches the classic kink band geometry, we model the development of a band of simple shear within an anisotropic perfectly plastic material. The resulting analytical model provides the stress state needed to maintain the kinematics of simple shear as a function of material anisotropy, deformation band orientation, and shear strain (or equivalently, material rotation). It is found that plastic anisotropy can promote either kink band or shear band formation depending on the loading orientation. When the deviatoric stress is positive parallel to the plane of anisotropy, shear localization may progress without bound and a shear band is produced. When the deviatoric stress is negative parallel to the plane of anisotropy, shear localization is arrested after a certain material rotation, resulting in a kink band. Examination of the requisite applied stress state during kink band formation provides an explanation for the experimentally-observed ‘lock up’ geometry. Solutions for the band boundary inclination angle are obtained and used to provide bounds on permissible band angles for both shear bands and kink bands.

36 MATERIALS SCIENCE↗

Midlatitude sporadic-E layers

Rocket borne probes and incoherent scatter radar were demonstrated to be effective methods of studying the structure of midlatitude sporadic E layers. Layers are formed when metal ions are converged vertically in a wind shear to produce a local enhancement of electron density. Rocket and radar observations show that the layers may occassionally have complex structure produced by an unstable wind shear. The partial transparency to radio waves of sporadic E layers is shown to be due to localized regions of high electron density.

Smith, L. G.↗

Observations of Laminar Surface Shear Stress Patterns on the Sidewall of a Rectangular, Mach 3 Quiet Tunnel Using Liquid Crystal Coatings

A shear-stress--sensitive liquid crystal coating (LCC) was used to visualize the surface shear stress distribution on the flat sidewall of the MSU quiet supersonic tunnel as a function of flow stagnation pressure. Under conditions of quiet operation, the LCC color-change response indicated the existence of a nonuniform surface shear stress distribution. This shear pattern was characterized by an elongated, down st ream-point ing triangular region of relatively low shear, with its apex on and its axis coincident with, the sidewall centerline. This low-shear zone was bounded symmetrically by two regions of relatively higher shear; these high-shear zones originated within the concave-curvature portion of the nozzle, in the corners between the flat sidewall and the contoured upper and lower nozzle surfaces. A 3-D Navier-Stokes code was used to compute the pressure and surface shear distributions on the sidewall. Flow-expansion-induced transverse pressure gradients on the nozzle sidewall generated symmetric inflows from the corners towards the sidewall centerline; these inflows caused a thickening of the sidewall boundary layer along the centerline, resulting in lower shear stresses consistent with the liquid crystal results. Peripherally nonuniform laminar boundary layer development, and the associated stability of such complex 3-D flows, must be considered in quiet-tunnel applications using rectangular nozzles. A color video will be shown.

Reda, D.↗

Navier-Stokes Computations With One-Equation Turbulence Model for Flows Along Concave Wall Surfaces

This report presents the use of a time-marching three-dimensional compressible Navier-Stokes equation numerical solver with a one-equation turbulence model to simulate the flow fields developed along concave wall surfaces without and with a downstream extension flat wall surface. The 3-D Navier- Stokes numerical solver came from the NASA Glenn-HT code. The one-equation turbulence model was derived from the Spalart and Allmaras model. The computational approach was first calibrated with the computations of the velocity and Reynolds shear stress profiles of a steady flat plate boundary layer flow. The computational approach was then used to simulate developing boundary layer flows along concave wall surfaces without and with a downstream extension wall. The author investigated the computational results of surface friction factors, near surface velocity components, near wall temperatures, and a turbulent shear stress component in terms of turbulence modeling, computational mesh configurations, inlet turbulence level, and time iteration step. The computational results were compared with existing measurements of skin friction factors, velocity components, and shear stresses of the developing boundary layer flows. With a fine computational mesh and a one-equation model, the computational approach could predict accurately the skin friction factors, near surface velocity and temperature, and shear stress within the flows. The computed velocity components and shear stresses also showed the vortices effect on the velocity variations over a concave wall. The computed eddy viscosities at the near wall locations were also compared with the results from a two equation turbulence modeling technique. The inlet turbulence length scale was found to have little effect on the eddy viscosities at locations near the concave wall surface. The eddy viscosities, from the one-equation and two-equation modeling, were comparable at most stream-wise stations. The present one-equation turbulence model is an effective approach for turbulence modeling in the near solid wall surface region of flow over a concave wall.

Wang, Chi R.↗

Microstructure and mechanical properties of ultrasonic spot welding of AA7075-T6 and A380 casting aluminum alloy

Ultrasonic spot welding (USW) was employed to join AA7075-T6 and casting Al alloy A380. Under a constant welding power, the influence of the welding time duration (from 1 s to 3 s) on peak temperature, microstructure, hardness, and mechanical strength of the joints were studied. It was found that the peak temperature increased with welding time. A sound solid-state metallurgical bonding was obtained at the joint interface. The cracks occurred at the edge of the joints were eventually eliminated as the welding time increased. A layer of refined grains formed in A380 adjacent to the bonding line. Further, hardness reduction occurred in AA7075 at the joint interface and the degree of hardness reduction and the volume increased with the peak temperature. By increasing welding time, lap shear tensile tests show that the peak load increased from 5.3 kN to 6.55 kN. The corresponding interfacial shear strength ranged from 96.29 MPa to 99.72 MPa, lower than the shear strength of both base materials (331 MPa for AA7075 and 185 MPa for A380). All joint samples exhibited interfacial fracture after lap shear tensile testing, with a thin layer of A380 remnant attached to AA7075 at the center of the fractured surface.

36 MATERIALS SCIENCE↗

An approximate semi-analytical method for prediction of interlaminar shear stresses in an arbitrarily laminated thick plate

An approximate semianalytical method for determination of interlaminar shear stress distribution through the thickness of an arbitrarily laminated thick plate has been presented. The method is based on the assumptions of transverse inextensibility and layerwise constant shear angle theory (LCST) and utilizes an assumed quadratic displacement potential energy based finite element method (FEM). Centroid of the triangular surface has been proved from a rigorous mathematical point of view (Aubin-Nitsche theory), to be the point of exceptional accuracy for the interlaminar shear stresses. Numerical results indicate close agreement with the available three-dimensional elasticity theory solutions. A comparison between the present theory and that due to an assumed stress hybrid FEM suggest that the (normal) traction-free-edge condition is not satisfied in the latter approach. Furthermore, the present paper is the first to present the results for interlaminar shear stresses in a two-layer thick square plate of balanced unsymmetric angle-ply construction. A comparison with the recently proposed Equilibrium Method (EM) indicates the superiority of the present method, because the latter assures faster convergence as well as simultaneous vanishing of the transverse shear stresses on both of the exposed surfaces of the laminate. Superiority of the present method over the EM, in the case of a symmetric laminate, is limited to faster convergence alone. It has also been demonstrated that the combination of the present method and the reduced (quadratic order) numerical integration scheme yields convergence of the interlaminar shear stresses almost as rapidly as that of the nodal displacements, in the case of a thin plate.

Chaudhuri, Reaz A.↗

An investigation of turbulent scatter from the mesosphere as observed by coherent-scatter radar

Turbulent scatter from he mesosphere is observed using the Urbana coherent-scatter radar. The variation in signal-to-noise ratio as a function of time-of-day is examined. The origin of scattering regions is investigated by comparing the variations in scattered power and Doppler velocity. Nighttime echoes are shown for periods of enhanced electron concentration. The spectrum of the returned signal is studied with a resolution of ten seconds. Spectral information is used to increase altitude resolution and observe the motion of scatterers. The expected variation in signal-to-noise ratio with solar flux is observed. It is found that variations in the scattered power generally do not correspond to the gravity waves which are simultaneously observed. Turbulent layers are observed at altitudes with high shear in the horizontal velocity and at altitudes with low shear. The ten-second resolution is necessary to distinguish meteor echoes from echoes produced by the advection of a scattering layer through the radar beam.

Gibbs, K. P.↗

A Study of Intermittent Turbulence in Stable Arctic Boundary Layers

Abstract Polar boundary layers are difficult to model due to the existence of intermittent turbulence within stable layers. Here we present a case study evaluation of coherent structures in a stable boundary layer observed during a series of flights with an uncrewed aircraft system (DataHawk2) on 19 October 2016 at Oliktok Point, Alaska as part of the ERASMUS (Evaluation of Routine Atmospheric Sounding Measurements using Unmanned Systems) field campaign. During a sequence of five flights over a nine-hour period, 57 profiles of atmospheric properties (0–400 m a.g.l) were collected. Turbulence was identified using derived Richardson Number, temperature structure function parameter, and turbulence kinetic energy dissipation. Throughout all flights on this strongly stable day, intermittent turbulence was observed. These turbulent layers showed well-mixed potential temperature profiles embedded within otherwise stable potential temperature profiles; often resulting in a characteristic staircase pattern. Turbulent layers ranged from 1 to 30 m deep, with most individual layers being 1–2 m deep. Vertical propagation velocities of layers in the lower atmosphere were on the order of a few cm s −1 , typical of non-convective environments. In different regions of the profile, turbulence was driven by a different balance of buoyancy and shear forces, with turbulence in the near surface environment driven by strong shear forces overcoming strong resistance to buoyancy, while turbulence in elevated layers characterized by weaker shear forces overcoming weaker resistance to buoyancy. We discuss the potential of such datasets for improving subgrid parameterizations of small-scale turbulence embedded within stable boundary layers.

54 ENVIRONMENTAL SCIENCES↗

Experimental Studies of Low-Pressure Turbine Flows and Flow Control. Streamwise Pressure Profiles and Velocity Profiles

This report summarizes research performed in support of the NASA Glenn Research Center (GRC) Low-Pressure Turbine (LPT) Flow Physics Program. The work was performed experimentally at the U.S. Naval Academy faculties. The geometry corresponded to "Pak B" LPT airfoil. The test section simulated LPT flow in a passage. Three experimental studies were performed: (a) Boundary layer measurements for ten baseline cases under high and low freestream turbulence conditions at five Reynolds numbers of 25,000, 50,000, 100,000, 200,000, and 300,000, based on passage exit velocity and suction surface wetted length; (b) Passive flow control studies with three thicknesses of two-dimensional bars, and two heights of three-dimensional circular cylinders with different spanwise separations, at same flow conditions as the 10 baseline cases; (c) Active flow control with oscillating synthetic (zero net mass flow) vortex generator jets, for one case with low freestream turbulence and a low Reynolds number of 25,000. The Passive flow control was successful at controlling the separation problem at low Reynolds numbers, with varying degrees of success from case to case and varying levels of impact at higher Reynolds numbers. The active flow control successfully eliminated the large separation problem for the low Reynolds number case. Very detailed data was acquired using hot-wire anemometry, including single and two velocity components, integral boundary layer quantities, turbulence statistics and spectra, turbulent shear stresses and their spectra, and intermittency, documenting transition, separation and reattachment. Models were constructed to correlate the results. The report includes a summary of the work performed and reprints of the publications describing the various studies.This report summarizes research performed in support of the NASA Glenn Research Center (GRC) Low-Pressure Turbine (LPT) Flow Physics Program. The work was performed experimentally at the U.S. Naval Academy faculties. The geometry corresponded to "Pak B" LPT airfoil. The test section simulated LPT flow in a passage. Three experimental studies were performed: (a) Boundary layer measurements for ten baseline cases under high and low freestream turbulence conditions at five Reynolds numbers of 25,000, 50,000, 100,000, 200,000, and 300,000, based on passage exit velocity and suction surface wetted length; (b) Passive flow control studies with three thicknesses of two-dimensional bars, and two heights of three-dimensional circular cylinders with different spanwise separations, at same flow conditions as the 10 baseline cases; (c) Active flow control with oscillating synthetic (zero net mass flow) vortex generator jets, for one case with low freestream turbulence and a low Reynolds number of 25,000. The Passive flow control was successful at controlling the separation problem at low Reynolds numbers, with varying degrees of success from case to case and varying levels of impact at higher Reynolds numbers. The active flow control successfully eliminated the large separation problem for the low Reynolds number case. Very detailed data was acquired using hot-wire anemometry, including single and two velocity components, integral boundary layer quantities, turbulence statistics and spectra, turbulent shear stresses and their spectra, and intermittency, documenting transition, separation and reattachment. Models were constructed to correlate the results. The report includes a summary of the work performed and reprints of the publications describing the various studies. The folders in this supplement contain processed data in ASCII format. Streamwise pressure profiles and velocity profiles are included. The velocity profiles were acquired using single sensor and cross sensor hot-wire probes which were traversed from the wall to the freestream at various streamwise locations. In some of the flow control cases (3D Trips and Jets) profiles were acquired at multiple spanwise locations.

Volino, Ralph↗