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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Turbulence Measurements in the Near Field of a Wingtip Vortex

The roll-up of a wingtip vortex, at Reynolds number based on chord of 4.6 million was studied with an emphasis on suction side and near wake measurements. The research was conducted in a 32 in. x 48 in. low-speed wind tunnel. The half-wing model had a semi-span of 36 in. a chord of 48 in. and a rounded tip. Seven-hole pressure probe measurements of the velocity field surrounding the wingtip showed that a large axial velocity of up to 1.77 U(sub infinity) developed in the vortex core. This level of axial velocity has not been previously measured. Triple-wire probes have been used to measure all components of the Reynolds stress tensor. It was determined from correlation measurements that meandering of the vortex was small and did not appreciably contribute to the turbulence measurements. The flow was found to be turbulent in the near-field (as high as 24 percent RMS w - velocity on the edge of the core) and the turbulence decayed quickly with streamwise distance because of the nearly solid body rotation of the vortex core mean flow. A streamwise variation of the location of peak levels of turbulence, relative to the core centerline, was also found. Close to the trailing edge of the wing, the peak shear stress levels were found at the edge of the vortex core, whereas in the most downstream wake planes they occurred at a radius roughly equal to one-third of the vortex core radius. The Reynolds shear stresses were not aligned with the mean strain rate, indicating that an isotropic-eddy-viscosity based prediction method cannot accurately model the turbulence in the cortex. In cylindrical coordinates, with the origin at the vortex centerline, the radial normal stress was found to be larger than the circumferential.

Chow, Jim↗
DOE OSTI2022

Microstructural changes of proton irradiated Hastelloy-N and in situ micropillar compression testing of one single grain at different local damage levels

Here, in situ micropillar compression was used to study the deformation of proton-irradiated Hastelloy-N at different damage levels. Multiple pillars were prepared from a single domain along the cross-section of 2.5 MeV proton-irradiated Hastelloy-N. Depending on the location of micropillars, the critical resolved shear stress was obtained as a function of local damage levels. Such an approach eliminates the variation of yield stress due to the difference in the Schmid factor. Microstructural characterization showed that irradiation hardening is contributed by complicated defect structures, including (a) dislocation loops with many in corduroy-like alignments, (2) dislocations pile up, and (3) Cottrell atmosphere. Silicon atoms are found to segregate at dislocation lines and loops and form complicated Cottrell cloud patterns. The critical resolved shear stress, $t_{crss}$, and hardening exponents were both extracted as a function of displacements per atom values up to 2.3. There was a 60% increase in $t_{crss}$ at the highest damage level.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1976

A normal shock-wave turbulent boundary-layer interaction at transonic speeds

Experimental results, including surveys of the mean and fluctuating flow, and measurements of surface pressure, skin friction, and separation length, are compared with solutions to the Navier-Stokes equations utilizing various algebraic eddy viscosity models to describe the Reynolds shear stresses. The experimental data, obtained at a free-stream Mach number of 1.5 and Reynolds numbers between 10 million and 80 million, show that a separated zone forms near the foot of the shock and that its length is proportional to the initial boundary-layer thickness; that a supersonic region forms downstream of the shock; and that the shear stress increases significantly through the interaction and subsequently decays downstream. The computations adequately represent the qualitative features of the flow field throughout the interaction but quantitatively underpredict the extent of separation and the downstream level of skin friction.

Mateer, G. G.↗
NASA NTRS2023

Exploring the Langtry-Menter Transition Model for High Speed Applications Using FUN3D

A series of Reynolds-averaged Navier-Stokes (RANS) simulations were performed using the FUN3D flow solver to explore the capabilities of the Langtry-Menter Shear-Stress Transport (LM-SST) transition model for predicting transition for aircraft inlet applications. Two geometries were simulated: a zeropressure-gradient flat plate and an axisymmetric cone exposed to hypersonic flow. In addition to the transition-sensitized LM-SST model investigations, simulations were run with the one-equation SpalartAllmaras (SA) and the two-equation Menter Shear-Stress Transport (SST-V) RANS models in fully turbulent mode to identify the natural RANS model transition behavior as a function of Mach number when executed in fully turbulent mode. The flat plate simulations showed that (1) the transition model was able to predict rapid transition at a freestream Mach number of 0.2, which is expected but (2) the predicted transition location moved downstream as the freestream Mach number was increased for the simulations that used the SST-V turbulence model. The latter is significant as it is usually assumed that one- and two-equation turbulence models will produce fully turbulent flow very near the boundary layer origin. The flat plate simulation freestream Mach number trend was confirmed with simulations using the Wind-US code, which also saw a similar trend when employing the SA turbulence model. For the axisymmetric cone simulations, the transition location was highly sensitive to the inflow turbulence levels. This is significant as the prediction of the transition location is crucial when trying to predict inlet performance, especially for hypersonic vehicle applications. It was also noted that the predicted transition location for the cone when using the SST-V turbulence model agreed well with the predicted transition location from the equivalent zero-pressure-gradient cold wall flat plate case.

Transition Model↗
NASA NTRS2012

Cubic Zig-Zag Enrichment of the Classical Kirchhoff Kinematics for Laminated and Sandwich Plates

A detailed anaylsis and examples are presented that show how to enrich the kinematics of classical Kirchhoff plate theory by appending them with a set of continuous piecewise-cubic functions. This analysis is used to obtain functions that contain the effects of laminate heterogeneity and asymmetry on the variations of the inplane displacements and transverse shearing stresses, for use with a {3, 0} plate theory in which these distributions are specified apriori. The functions used for the enrichment are based on the improved zig-zag plate theory presented recently by Tessler, Di Scuva, and Gherlone. With the approach presented herein, the inplane displacements are represented by a set of continuous piecewise-cubic functions, and the transverse shearing stresses and strains are represented by a set of piecewise-quadratic functions that are discontinuous at the ply interfaces.

Nemeth, Michael P.↗
NASA NTRS1976

Oceanic lithosphere and asthenosphere - Thermal and mechanical structure

A coupled thermomechanical subsolidus model of the oceanic lithosphere and asthenosphere is developed which includes vertical heat conduction, a temperature-dependent thermal conductivity, heat advection by a horizontal and vertical mass flow that depends on depth and age, contributions of viscous dissipation or shear heating, a linear or nonlinear deformation law relating shear stress and strain rate, as well as a temperature- and pressure-dependent viscosity. The model requires a constant horizontal velocity and temperature at the surface, but zero horizontal velocity and constant temperature at great depths. The depth- and age-dependent temperature, horizontal and vertical velocities, and viscosity structure of the lithosphere and asthenosphere are determined along with the age-dependent shear stress in those two zones. The ocean-floor topography, oceanic heat flow, and lithosphere thickness are deduced as functions of ocean-floor age; seismic velocity profiles which exhibit a marked low-velocity zone are constructed from the age-dependent geotherms and assumed values of the elastic parameters. It is found that simple boundary-layer cooling determines the thermal structure at young ages, while effects of viscous dissipation become more important at older ages.

Schubert, G.↗
NASA NTRSDate not supplied

Exploring the Langtry-Menter Transition Model for High Speed Applications Using FUN3D

A series of Reynolds-averaged Navier-Stokes (RANS) simulations were performed using the FUN3D flow solver to explore the capabilities of the Langtry-Menter Shear-Stress Transport (LM-SST) transition model for predicting transition for aircraft inlet applications. Two geometries were simulated: a zero-pressure-gradient flat plate and an axisymmetric cone exposed to hypersonic flow. In addition to the transition-sensitized LM-SST model investigations, simulations were run with the one-equation Spalart-Allmaras (SA) and the two-equation Menter Shear-Stress Transport (SST-V) RANS models in fully turbulent mode to identify the natural RANS model transition behavior as a function of Mach number when executed in fully turbulent mode. The flat plate simulations showed that (1) the transition model was able to predict rapid transition at a freestream Mach number of 0.2, which is expected but (2) the predicted transition location moved downstream as the freestream Mach number was increased for the simulations that used the SST-V turbulence model. The latter is significant as it is usually assumed that one- and two-equation turbulence models will produce fully turbulent flow very near the boundary layer origin. The flat plate simulation freestream Mach number trend was confirmed with simulations using the Wind-US code, which also saw a similar trend when employing the SA turbulence model. For the axisymmetric cone simulations, the transition location was highly sensitive to the inflow turbulence levels. This is significant as the prediction of the transition location is crucial when trying to predict inlet performance, especially for hypersonic vehicle applications. It was also noted that the predicted transition location for the cone when using the SST-V turbulence model agreed well with the predicted transition location from the equivalent zero-pressure-gradient flat plate case.

Transition Model↗
NASA NTRS1974

Three-dimensional turbulent boundary layer in a rotating helical channel

Analytical and experimental investigation of the characteristics of a three-dimensional turbulent boundary layer in a rotating helical channel. Expressions are developed for the velocity profiles in the inner layer, where viscous effects dominate, and the outer layer, where viscous effects are small. The velocity profiles, wall shear stress, and limiting streamline angles are measured inside the passages of a flat-plate inducer at various radial and chordwise locations using rotating probes. Flow near the blade tip is found to be highly complex, due to interaction of blade boundary layers and annulus wall, resulting in appreciable radial inward flow, as well as a defect in mainstream velocity near the mid-passage. A wall shear stress correlation, which includes the effect of both Reynolds number and rotation parameter, is derived from the measured data.

Anand, A. K.↗
NASA NTRS2015

Evaluation of Industry Standard Turbulence Models on an Axisymmetric Supersonic Compression Corner

Reynolds-averaged Navier-Stokes computations of a shock-wave/boundary-layer interaction (SWBLI) created by a Mach 2.85 flow over an axisymmetric 30-degree compression corner were carried out. The objectives were to evaluate four turbulence models commonly used in industry, for SWBLIs, and to evaluate the suitability of this test case for use in further turbulence model benchmarking. The Spalart-Allmaras model, Menter's Baseline and Shear Stress Transport models, and a low-Reynolds number k- model were evaluated. Results indicate that the models do not accurately predict the separation location; with the SST model predicting the separation onset too early and the other models predicting the onset too late. Overall the Spalart-Allmaras model did the best job in matching the experimental data. However there is significant room for improvement, most notably in the prediction of the turbulent shear stress. Density data showed that the simulations did not accurately predict the thermal boundary layer upstream of the SWBLI. The effect of turbulent Prandtl number and wall temperature were studied in an attempt to improve this prediction and understand their effects on the interaction. The data showed that both parameters can significantly affect the separation size and location, but did not improve the agreement with the experiment. This case proved challenging to compute and should provide a good test for future turbulence modeling work.

Turbulent Boundary Layer↗
NASA NTRS2003

Analysis of Plasma-Sprayed Thermal Barrier Coatings With Homogeneous and Heterogeneous Bond Coats Under Spatially Uniform Cyclic Thermal Loading

This report summarizes the results of a numerical investigation into the spallation mechanism in plasma-sprayed thermal barrier coatings observed under spatially-uniform cyclic thermal loading. The analysis focuses on the evolution of local stress and inelastic strain fields in the vicinity of the rough top/bond coat interface during thermal cycling, and how these fields are influenced by the presence of an oxide film and spatially uniform and graded distributions of alumina particles in the metallic bond coat aimed at reducing the top/bond coat thermal expansion mismatch. The impact of these factors on the potential growth of a local horizontal delamination at the rough interface's crest is included. The analysis is conducted using the Higher-Order Theory for Functionally Graded Materials with creep/relaxation constituent modeling capabilities. For two-phase bond coat microstructures, both the actual and homogenized properties are employed in the analysis. The results reveal the important contributions of both the normal and shear stress components to the delamination growth potential in the presence of an oxide film, and suggest mixed-mode crack propagation. The use of bond coats with uniform or graded microstructures is shown to increase the potential for delamination growth by increasing the magnitude of the crack-tip shear stress component.

Arnold, Steven M.↗
NASA NTRS1997

Microgravity Foam Structure and Rheology

To exploit rheological and multiple-light scattering techniques, and ultimately microgravity conditions, in order to quantify and elucidate the unusual elastic character of foams in terms of their underlying microscopic structure and dynamics. Special interest is in determining how this elastic character vanishes, i.e. how the foam melts into a simple viscous liquid, as a function of both increasing liquid content and shear strain rate. The unusual elastic character of foams will be quantified macroscopically by measurement of the shear stress as a function of static shear strain, shear strain rate, and time following a step strain; such data will be analyzed in terms of a yield stress, a static shear modulus, and dynamical time scales. Microscopic information about bubble packing and rearrangement dynamics, from which these macroscopic non-Newtonian properties presumably arise, will be obtained non-invasively by novel multiple-light scattering diagnostics such as Diffusing-Wave Spectroscopy (DWS). Quantitative trends with materials parameters, such as average bubble size, and liquid content, will be sought in order to elucidate the fundamental connection between the microscopic structure and dynamics and the macroscopic rheology.

Durian, Douglas J.↗
NASA NTRS2003

Numerical Simulation of High Drag Reduction in a Turbulent Channel Flow with Polymer Additives

The addition of small amounts of long chain polymer molecules to wall-bounded flows can lead to dramatic drag reduction. Although this phenomenon has been known for about fifty years, the action of the polymers and its effect on turbulent structures are still unclear. Detailed experiments have characterized two distinct regimes (Warholic et al. 1999), which are referred to as low drag reduction (LDR) and high drag reduction (HDR). The first regime exhibits similar statistical trends as Newtonian flow: the log-law region of the mean velocity profile remains parallel to that of the Newtonian ow but its lower bound moves away from the wall and the upward shift of the log-region is a function of drag reduction, DR. Although streamwise fluctuations are increased and transverse ones are reduced, the shape of the rms velocity profiles is not qualitatively modified. At higher drag reductions, of the order of 40-50%, the ow enters the HDR regime for which the slope of the log-law is dramatically augmented and the Reynolds shear stress is small (Warholic et al. 1999; Ptasinski et al. 2001). The drag reduction is eventually bounded by a maximum drag reduction (MDR) (Virk & Mickley 1970) which is a function of the Reynolds number. While several experiments report mean velocity profiles very close to the empirical profile of Virk & Mickley (1970) for MDR conditions, the observations regarding the structure of turbulence can differ significantly. For instance, Warholic et al. (1999) measured a near-zero Reynolds shear stress, whereas a recent experiment (Ptasinski et al. 2001) shows evidence of non-negligible Reynolds stress in their MDR flow. To the knowledge of the authors, only the LDR regime has been documented in numerical simulations (Sureshkumar et al. 1997; Dimitropoulos et al. 1998; Min et al. 2001; Dubief & Lele 2001; Sibilla & Baron 2002). This paper discusses the simulation of polymer drag reduced channel ow at HDR using the FENE-P (Finite Elastic non-linear extensibility-Peterlin) model which was used for the first LDR simulation by Sureshkumar et al. (1997). Flow and polymer parameters are close to realistic polymer drag reducing conditions. High drag reductions are achieved by using finite differences and a robust time stepping technique. A minimal channel flow is also used as a numerical experiment to investigate the effect of the outer region turbulent structures on the overall drag at HDR. The drag reducing action of the model is finally studied through the structure of energy transfers from the polymers to the velocity components. This investigation sheds some light on the details of polymer drag reduction.

Dubief, Yves↗
NASA NTRS1984

Turbulence measurements in a complex flowfield using a crossed hot-wire

An X-wire probe was used to measure the time-mean and fluctuating velocities and shear stress in nonswirling nonreacting confined jet flows. Data were taken from an axisymmetric confined jet with an expansion ratio of 2 and an expansion angle of 90 deg, and from the same segment with a contraction nozzle. Velocity profiles developed faster in the confined jet than in the free jet, with the former experiencing higher turbulence levels and larger time-mean velocities. The X-wire is concluded to furnish more accurate results for the turbulent shear stress than a multioriented single-wire technique.

Lilley, D. G.↗
NASA NTRS1985

Possible role of crustal flexure in the initial detachment of extensional allochthons

The existence of low-angle normal faults indicates that the ratio of shear stress (TAU) to normal stress (sigma sub N) needed to cause slip on faults is substantially less than would be predicted based on experimental data. Because the tensional strength of rock at a large scale is exceedingly low, the upper plate of a low-angle normal fault cannot be pulled down the fault ramp, but must be driven down it by its own weight. The active or recently active Sevier Desert detachment fault in western Utah dips regionally at 12 deg. The radio of shear stress to normal stress due to the weight of the upper plate on a 12 deg dipping fault surface is 0.2. In contrast, laboratory experiments indicate that slip on fracture surfaces occurs with almost all rock types when (tau/sigma) reaches values of 0.6 to 0.85, corresponding to normal-fault dips of 30 deg to 40 deg. Seismological data indicate that low deviatoric stresses are associated with movement on faults of other geometrics and are not unique to low-angle normal faults. It thus appears that approximately planar fault zones with surface areas of hundreds to thousands of square kilometers have different mechanical properties than would be predicted based on laboratory studies of fractured rock. Modeling of stresses indicates that deviatoric stresses greater than 5kb exist and are sustainable in oceanic lithosphere, and that failure occurs when tau/sigma sub N approaches 0.6 to 0.85.

Spencer, J. E.↗
NASA NTRS2019

The Intermediate Wake of a Thin Flat Plate with a Circular Trailing Edge

The intermediate wakes of thin flat plates with circular trailing edges (TEs) are investigated here with direct numerical simulations (DNSs). The separating boundary layers are turbulent in all cases. The near wake in two thin-plate cases (IN & NS), with a focus on the vortex shedding process, was explored in a recent article. Intermittent shedding was observed in Case IN. Case NS, with half the TE diameter of Case IN, was an essentially non-shedding case. A third case (ST) with a sharp trailing edge was also investigated and found to exhibit an intermittent wake instability. The objectives of the present study are twofold. The first is to determine if the wake instability found in Case ST exists in Cases IN and NS as well. The second is to provide the distributions of the turbulent normal intensities and shear stress in the wake and to understand these distributions via the budget terms in the corresponding transport equations. The results show that both Cases IN & NS exhibit a wake instability in the intermediate wake region, that is similar to that found earlier in Case ST. We note that in Case IN, the presence of an intermediate-wake instability results in the co-existence of two different types of instability within a single wake. The distributions of the turbulent normal intensities and shear stress, and the budget terms for the streamwise intensity are included and discussed here. All the budget terms contribute appreciably to the overall budget in the transport equation for streamwise normal intensity.

Rai, Man Mohan↗
DOE OSTI2020

Evaluating Implementations of the Immersed Boundary Method in the Weather Research and Forecasting Model

The terrain-following coordinate system used by many atmospheric models can cause numerical instabilities due to discretization errors as resolved terrain slopes increase and the grid becomes highly skewed. The immersed boundary (IB) method, which does not require the grid to conform to the terrain, has been shown to alleviate these errors, and has been used successfully for high-resolution atmospheric simulations over steep terrain, including vertical building surfaces. Since many previous applications of IB methods to atmospheric models have used very fine grid resolution (5 m or less), the present study seeks to evaluate IB method performance over a range of grid resolutions and aspect ratios. Two classes of IB algorithms, velocity reconstruction and shear stress reconstruction, are tested within the common framework of the Weather Research and Forecasting (WRF) Model. Performance is evaluated in two test cases, one with flat terrain and the other with the topography of Askervein Hill, both under neutrally stratified conditions. WRF-IB results are compared to similarity theory, observations, and native WRF results. Despite sensitivity to the location at which the IB intersects the model grid, the velocity reconstruction IB method shows consistent performance when used with a hybrid RANS/LES surface scheme. The shear stress reconstruction IB method is not sensitive to the grid intersection, but is less consistent and near-surface velocity errors can occur at coarse resolutions. This study represents an initial investigation of IB method variability across grid resolutions in WRF. Future work will focus on improving IB method performance at intermediate to coarse resolutions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Scaling patch analysis of turbulent planar plume

Proper scaling in turbulent planar plumes is investigated in this study using a scaling patch approach. Based on the scaled boundary conditions, a proper velocity scale for the mean axial flow is the plume centerline velocity U ref = U ctr , and a proper temperature scale for the temperature excess is Θ ref = T ctr – T ∞ , where T ctr is the plume centerline temperature and T ∞ is the ambient fluid temperature. By seeking an admissible scaling, a key concept in the scaling patch approach, for the mean continuity, mean momentum, and mean energy equations, respectively, the following is found: (1) a proper scale for the mean transverse flow is V ref = ( d δ / d x ) U ctr , where d δ / d x is the growth rate of the plume width. (2) A proper scale for the Reynolds shear stress is R vu , ref = U ctr V ref = ( d δ / d x ) U ctr 2 , a mix of the scales for the mean axial and transverse flows. (3) A proper scale for the turbulent heat flux is R v θ , ref = V ref Θ ctr , a mix of the scales for the mean transverse flow and mean temperature excess. The mean transverse flow thus plays a critical role in the scaling of turbulent planar plumes. Approximate functions are developed for the scaled mean transverse flow, Reynolds shear stress, and turbulent temperature flux, and are found to agree favorably with experimental and numerical simulation data. The integral analysis of the mean momentum equation yields a Richardson number Ri, which remains invariant in the axial direction. The Richardson number is defined as Ri = def g β Θ ctr δ t / ( U ctr V ref ) ≈ 1 / 2 , where g is the gravitational acceleration, β is the thermal expansion coefficient, and δ t is the plume half-width based on the mean temperature profile. This Richardson number arises directly from the scaling patch analysis of the mean momentum equation, including both the streamwise and transverse velocity scales.

42 ENGINEERING↗
NASA NTRS1995

A preliminary attempt to use neural networks for turbulent eddy classification

This note describes an attempt to use standard neural network tools to fashion a means of detecting eddy patterns in the wall region of a turbulent flow. The research was motivated by the desire to formulate a means to use only flow parameters that can be sensed on the wall to describe the passing eddy structure. If a simple formulation can be obtained, it could conceivably be utilized to control actuators embedded in the wall. Such actuators have been developed by Jacobson and Reynolds (1993a), Blackwelder and Liu (1994), Tung et al. (1995), and others. These actuators have the common characteristics that they are small and are typically flush with the wall when not deployed. When they are activated, it is assumed that they will be able to interact constructively with the turbulent eddies near their location to either decrease the wall shear stress, enhance or reduce the mixing, etc. At present, there is only a nascent understanding of the interaction dynamics between the actuators and the eddies in the flow. Nevertheless, for such interaction to succeed, methods to couple the actuators to the oncoming flow must be obtained. General methods must be found that will detect the space and temporal location of the desired structure. In particular, it will be necessary to know when the eddies will arrive at the location of the actuator. This research attempted to use the shear stress measurements on the wall in the vicinity of an actuator location to predict when a particular eddy pattern would arrive and/or occur at the designated location. In this work the eddy pattern to be detected was identified by its velocity signature only.

Blackwelder, Ron F.↗