Search NASA⌕ Search

SEARCH · Search NASA

Results for “shear layer turbulence”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

An alternative assessment of second-order closure models in turbulent shear flows

The performance of three recently proposed second-order closure models is tested in benchmark turbulent shear flows. Both homogeneous shear flow and the log-layer of an equilibrium turbulent boundary layer are considered for this purpose. An objective analysis of the results leads to an assessment of these models that stands in contrast to that recently published by other authors. A variety of pitfalls in the formulation and testing of second-order closure models are uncovered by this analysis.

Speziale, Charles G.↗

Reynolds shear stress development in pressure-driven three-dimensional turbulent boundary-layers

The development of the Reynolds stresses has been examined experimentally in an initially two-dimensional boundary layer which is driven to three dimensionality by a spanwise pressure gradient. The pressure field was imposed by an upstream-facing wedge. Two different wedge angles were used in order to vary the level of boundary layer skewing. Bradshaw's Al parameter was found to decrease with the rate of decrease being dependent on the level of skewing between the freestream and the wall flow. It was also concluded that the ratio of the cross-stream to streamwise shear stress components was governed by the rate of freestream turning.

Anderson, S. D.↗

Visualization of Capsule Reentry Vehicle Heat Shield Ablation Using Naphthalene PLIF

The Orion Multi-Purpose Crew Vehicle (MPCV) will use an ablative heat shield and improved understanding of the ablation process would be beneficial for design purposes. Given that ablation is a multi-physics process involving heat and mass transfer, codes aiming to predict heat shield ablation are in need of experimental data pertaining to the turbulent transport of ablation products for validation. At The University of Texas at Austin, a technique is being developed that uses planar laser-induced fluorescence (PLIF) of a low-temperature sublimating ablator (naphthalene) to visualize the transport of ablation products in a supersonic flow. Since ablation at reentry temperatures can be difficult to recreate in a laboratory setting it is desirable to create a limited physics problem and simulate the ablation process at relatively low temperature conditions using naphthalene. A scaled Orion MPCV model with a solid naphthalene heat shield has been tested in a Mach 5 wind tunnel at various angles of attack in the current work. PLIF imaging reveals the distribution of the ablation products as they are transported into the heat-shield boundary layer and over the capsule shoulders into the separated shear layer and backshell recirculation region. Visualizations of the capsule shear layer using both naphthalene PLIF and Schlieren imaging compared favorably. High concentrations of naphthalene in the capsule separated flow region, intermittent turbulent structures on the heat shield surface, and interesting details of the capsule shear layer structure were observed using the naphthalene PLIF technique. The capsule shear layer was also shown to generally appear to be more turbulent at lower angles of attack. Furthermore, the PLIF signal increased steadily over the course of a run indicating that during a wind tunnel run the model heated up and the rate of naphthalene ablation increased. The shear layer showed increasing signs of turbulence over the course of a wind tunnel run as well, likely because of the combination of increased surface roughness and surface blowing rate. PLIF imaging also detected regions with a relatively low concentration of naphthalene in the capsule backshell recirculation region that are most likely the result of cross-flow-induced vortices on the capsule afterbody.

Combs, Christopher S.↗

The diurnal cycle of entrainment and detrainment in LES of the Southern Ocean driven by observed surface fluxes and waves

Empirical rules for both entrainment and detrainment are developed from LES of the Southern Ocean boundary layer when the turbulence, stratification, and shear cannot be assumed to be in equilibrium with diurnal variability in surface flux and wave (Stokes drift) forcing. A major consequence is the failure of downgradient eddy viscosity, which becomes more serious with Stokes drift and is overcome by relating the angle between the stress and shear vectors to the orientations of Lagrangian shear to the surface and of local Eulerian shear over 5 m. Thus, the momentum flux can be parameterized as a stress magnitude and this empirical direction. In addition, the response of a deep boundary layer to sufficiently strong diurnal heating includes boundary layer collapse and the subsequent growth of a morning boundary layer, whose depth is empirically related to the time history of the forcing, as are both morning detrainment and afternoon entrainment into weak diurnal stratification. Below the boundary layer, detrainment rules give the maximum buoyancy flux and its depth, as well a specific stress direction. Another rule relates both afternoon and nighttime entrainment depth and buoyancy flux to surface layer turbulent kinetic energy production integrals. These empirical relationships are combined with rules for boundary layer transport to formulate two parameterizations; one based on eddy diffusivity and viscosity profiles and another on flux profiles of buoyancy and of stress magnitude. Evaluations against LES fluxes show the flux profiles to be more representative of the diurnal cycle, especially with Stokes drift.

54 ENVIRONMENTAL SCIENCES↗

Flow Phenomena in the Very Near Wake of a Flat Plate with a Circular Trailing Edge

The very near wake of a flat plate with a circular trailing edge, exhibiting pronounced shedding of wake vortices, is investigated with data from a direct numerical simulation. The separating boundary layers are turbulent and statistically identical thus resulting in a wake that is symmetric in the mean. The focus here is on the instability of the detached shear layers, the evolution of rib-vortex induced localized regions of reverse flow that detach from the main body of reverse flow in the trailing edge region and convect downstream, and phaseaveraged velocity statistics in the very near wake. The detached shear layers are found to exhibit unstable behavior intermittently, including the development of shear layer vortices as in earlier cylinder flow investigations with laminar separating boundary layers. Only a small fraction of the separated turbulent boundary layers undergo this instability, and form the initial shed vortices. Pressure spectra within the shear layers show a broadband peak at a multiple of shedding frequency. Phase-averaged intensity and shear stress distributions of the randomly fluctuating component of velocity are compared with those obtained in the near wake. The distributions of the production terms in the transport equations for the turbulent stresses are also provided.

Turbulence↗

Turbulence structures associated with the bursting event

Turbulence structures in a wall bounded shear layer during the bursting event detected by a conditional sampling technique are investigated using data obtained from large eddy simulation of turbulent channel flow. Streamlines are constructed from the ensemble-averaged velocity field to illustrate the flow patterns associated with the bursting event. The spatting motions during the sweep event and the existence of a pair of counterrotating streamwise vortices during the ejection process are exhibited.

Kim, J.↗

Turbulence structures associated with the bursting event

Turbulence structures in a wall bounded shear layer during the bursting event detected by a conditional sampling technique are investigated using data obtained from large eddy simulation of turbulent channel flow. Streamlines are constructed from the ensemble-averaged velocity field to illustrate the flow patterns associated with the bursting event. The splatting motions during the sweep event and the existence of a pair of counterrotating streamwise vortices during the ejection process are exhibited.

Kim, J.↗

An integral turbulent kinetic energy analysis of free shear flows

Mixing of coaxial streams is analyzed by application of integral techniques. An integrated turbulent kinetic energy (TKE) equation is solved simultaneously with the integral equations for the mean flow. Normalized TKE profile shapes are obtained from incompressible jet and shear layer experiments and are assumed to be applicable to all free turbulent flows. The shear stress at the midpoint of the mixing zone is assumed to be directly proportional to the local TKE, and dissipation is treated with a generalization of the model developed for isotropic turbulence. Although the analysis was developed for ducted flows, constant-pressure flows were approximated with the duct much larger than the jet. The axisymmetric flows under consideration were predicted with reasonable accuracy. Fairly good results were also obtained for the fully developed two-dimensional shear layers, which were computed as thin layers at the boundary of a large circular jet.

Peters, C. E.↗

Time Resolved Particle Image Velocimetry in Nozzle Flow

Temporally resolved particle image velocimetry (TR-PIV) is the newest and most exciting tool recently developed to support our continuing efforts to characterize and improve our understanding of the decay of turbulence in jet flows—a critical element for understanding the acoustic properties of the flow. A new TR-PIV system has been developed at the NASA Glenn Research Center which is capable of acquiring planar PIV image frame pairs at up to 25 kHz. The data reported here were collected at Mach numbers of 0.5 and 0.9 and at temperature ratios of 0.89 and 1.76. The field of view of the TR-PIV system covered six nozzle diameters along the lip line of the 50.8 mm diameter jet. The cold flow data at Mach 0.5 were compared with hotwire anemometry measurements in order to validate the new TR-PIV technique. The axial turbulence profiles measured across the shear layer using TR-PIV were thinner than those measured using hotwire anemometry and remained centred along the nozzle lip line. The collected TR-PIV data illustrate the differences in the single point statistical flow properties of cold and hot jet flows. The planar, time-resolved velocity records were then used to compute two-point space–time correlations of the flow at the Mach 0.9 flow condition. The TR-PIV results show that there are differences in the convective velocity and growth rate of the turbulent structures between cold and hot flows at the same Mach number.

Temporally Resolved Particle Image Velocimetry↗

Time-Resolved PIV for Space-Time Correlations in Hot Jets

Temporally Resolved Particle Image Velocimetry (TR-PIV) is being used to characterize the decay of turbulence in jet flows a critical element for understanding the acoustic properties of the flow. A TR-PIV system, developed in-house at the NASA Glenn Research Center, is capable of acquiring planar PIV image frame pairs at up to 10 kHz. The data reported here were collected at Mach numbers of 0.5 and 0.9 and at temperature ratios of 0.89 and 1.76. The field of view of the TR-PIV system covered 6 nozzle diameters along the lip line of the 50.8 mm diameter jet. The cold flow data at Mach 0.5 were compared with hotwire anemometry measurements in order to validate the new TR-PIV technique. The axial turbulence profiles measured across the shear layer using TR-PIV were thinner than those measured using hotwire anemometry and remained centered along the nozzle lip line. The collected TR-PIV data illustrate the differences in the single point statistical flow properties of cold and hot jet flows. The planar, time-resolved velocity records were then used to compute two-point space-time correlations of the flow at the Mach 0.9 flow condition. The TR-PIV results show that there are differences in the convective velocity and growth rate of the turbulent structures between cold and hot flows at the same Mach number.

Wernet, Mark P.↗

Time Resolved PIV for Space-Time Correlations in Hot Jets

Temporally Resolved Particle Image Velocimetry (TR-PIV) is the newest and most exciting tool recently developed to support our continuing efforts to characterize and improve our understanding of the decay of turbulence in jet flows -- a critical element for understanding the acoustic properties of the flow. A new TR-PIV system has been developed at the NASA Glenn Research Center which is capable of acquiring planar PIV image frame pairs at up to 25 kHz. The data reported here were collected at Mach numbers of 0.5 and 0.9 and at temperature ratios of 0.89 and 1.76. The field of view of the TR-PIV system covered 6 nozzle diameters along the lip line of the 50.8 mm diameter jet. The cold flow data at Mach 0.5 were compared with hotwire anemometry measurements in order to validate the new TR-PIV technique. The axial turbulence profiles measured across the shear layer using TR-PIV were thinner than those measured using hotwire anemometry and remained centered along the nozzle lip line. The collected TR-PIV data illustrate the differences in the single point statistical flow properties of cold and hot jet flows. The planar, time-resolved velocity records were then used to compute two-point space-time correlations of the flow at the Mach 0.9 flow condition. The TR-PIV results show that there are differences in the convective velocity and growth rate of the turbulent structures between cold and hot flows at the same Mach number

Wernet, Mark P.↗

Toward large eddy simulation of turbulent flow over an airfoil

The flow field over an airfoil contains several distinct flow characteristics, e.g. laminar, transitional, turbulent boundary layer flow, flow separation, unstable free shear layers, and a wake. This diversity of flow regimes taxes the presently available Reynolds averaged turbulence models. Such models are generally tuned to predict a particular flow regime, and adjustments are necessary for the prediction of a different flow regime. Similar difficulties are likely to emerge when the large eddy simulation technique is applied with the widely used Smagorinsky model. This model has not been successful in correctly representing different turbulent flow fields with a single universal constant and has an incorrect near-wall behavior. Germano et al. (1991) and Ghosal, Lund & Moin have developed a new subgrid-scale model, the dynamic model, which is very promising in alleviating many of the persistent inadequacies of the Smagorinsky model: the model coefficient is computed dynamically as the calculation progresses rather than input a priori. The model has been remarkably successful in prediction of several turbulent and transitional flows. We plan to simulate turbulent flow over a '2D' airfoil using the large eddy simulation technique. Our primary objective is to assess the performance of the newly developed dynamic subgrid-scale model for computation of complex flows about aircraft components and to compare the results with those obtained using the Reynolds average approach and experiments. The present computation represents the first application of large eddy simulation to a flow of aeronautical interest and a key demonstration of the capabilities of the large eddy simulation technique.

Choi, Haecheon↗

Turbulence characteristics inside a turbulent spot in plane Poiseuille flow

Turbulence characteristics inside a turbulent spot in plane Poiseuille flow are investigated by analyzing a database obtained from a direct simulation. The spot area is divided into two distinct regions - a turbulent area and a wave area. It is found that the flow structures inside the turbulent area have strong resemblance to those found in the fully-developed turbulent channel flow. A suitably defined mean and rms fluctuations as well as the internal shear-layer structures are found to be similar to the turbulent counterpart. In the wave area the inflexional mean spanwise profiles cause a rapid growth of oblique waves, which break down to turbulence. The rms fluctuations and Reynolds stress are found to be higher in that area, and the shear-layer structures are similar to those observed in the secondary instability of two-dimensional Tollmien-Schlichting waves.

Henningson, Dan S.↗

Evaluation of Baldwin-Barth turbulence model with thin-layer Navier-Stokes computation of an axisymmetric afterbody-exhaust jet flowfield

Axisymmetric steady flowfield over a steep circular arc afterbody with a choked exhaust jet in a subsonic freestream is numerically investigated using overflow code, developed by NASA Ames. The accuracy of Baldwin-Barth one-equation turbulence model for the computation of separated flow and plume shear layers is assessed by a comparison of predictions with available measurements from NASA Langley. The results show that the Baldwin-Barth turbulence model performs better than the Baldwin-Lomax algebraic model in predicting the surface pressure distribution in the separated flow region of the afterbody. In the jet mixing layer, the Baldwin-Barth model considerably unpredicts the growth of the shear layer, as does the Baldwin-Lomax model. It is shown that the one-equation turbulence model needs improvement for the prediction of separated flow due to smooth adverse pressure gradient, and mixing layer growth.

Kandula, M.↗

Coherent motion induced fluctuations in the primary transition region of a plane shear layer

The naturally occurring large scale motions in a single stream shear layer (that is initiated from a fully turbulent boundary layer) are made evident by the induced velocities in the entrainment region beyond the active shear layer. The distinctive attributes of these induced motions are particularly evident in the Michigan State Univerity Free Shear Flow Facility since the total test section length (3m) is nominally the same as the location of the first, fully formed, coherent motion, ca/x theta (0) = 400 (or 2.5 m). Hence, detailed studies of the induced motions can be executed. Individual coherent motions are identified by the induced velocity signatures and conditional-ensemble statistics are used to represent the irrotational field properties. Clusters of such motions exist; some of their properties are substantially different from the unconditionally averaged values.

Foss, J. F.↗

Predicting equilibrium states with Reynolds stress closures in channel flow and homogeneous shear flow

Turbulent channel flow and homogeneous shear flow have served as basic building block flows for the testing and calibration of Reynolds stress models. A direct theoretical connection is made between homogeneous shear flow in equilibrium and the log-layer of fully-developed turbulent channel flow. It is shown that if a second-order closure model is calibrated to yield good equilibrium values for homogeneous shear flow it will also yield good results for the log-layer of channel flow provided that the Rotta coefficient is not too far removed from one. Most of the commonly used second-order closure models introduce an ad hoc wall reflection term in order to mask deficient predictions for the log-layer of channel flow that arise either from an inaccurate calibration of homogeneous shear flow or from the use of a Rotta coefficient that is too large. Illustrative model calculations are presented to demonstrate this point which has important implications for turbulence modeling.

Abid, R.↗

Predicting equilibrium states with Reynolds stress closures in channel flow and homogeneous shear flow

Turbulent channel flow and homogeneous shear flow have served as basic building block flows for the testing and calibration of Reynolds stress models. A direct theoretical connection is made between homogeneous shear flow in equilibrium and the log-layer of fully-developed turbulent channel flow. It is shown that if a second-order closure model is calibrated to yield good equilibrium values for homogeneous shear flow it will also yield good results for the log-layer of channel flow provided that the Rotta coefficient is not too far removed from one. Most of the commonly used second-order closure models introduce an ad hoc wall reflection term in order to mask deficient predictions for the log-layer of channel flow that arise either from an inaccurate calibration of homogeneous shear flow or from the use of a Rotta coefficient that is too large. Illustrative model calculations are presented to demonstrate this point which has important implications for turbulence modeling.

Abid, R.↗

Numerical Study on the Aerodynamics of an Iced Airfoil With Scale-Resolving Simulations

Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.

ARMD↗