Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulent length scale”

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 55 records · Page 3

Temperature and velocity measurements in premixed turbulent flames

Turbulent flame speed data for premixed flames of methane-air, propane-air and ethylene-air mixtures stabilized in grid turbulence are reported and discussed. It is shown that turbulence effects on flame speed cannot be fully correlated by the turbulence length scale and r.m.s. velocity in the cold flow. Rather there appear to be significant flame-flow-turbulence interactions affecting both turbulence level in the reaction zone and measured flame speeds. Results of detailed velocity measurements, including autocorrelations, by laser velocimetry are used to elucidate the nature of these interactions. It is concluded that flame speed experiments must be designed and conducted to provide sufficient information (e.g., boundary conditions) to allow for reconstruction of the flow field and these interactions by modelers if the data are to be of value in turbulent combustion model development and evaluation.

Dandekar, K. V.↗

A dynamic subgrid-scale model for LES of the G-equation

Turbulent combustion is a difficult subject as it must deal with all of the issues found in both turbulence and combustion. (We consider only premixed flames in this paper, but some of the ideas can be applied to the non-premixed case.) As in many other fields, there are two limiting cases that are easier to deal with than the general case. These are the situations in which the chemical time scale is either much shorter or much longer than the time scale associated with the turbulence. We deal with the former case. In this limit, the flame is thin compared to the turbulence length scales and can be idealized as an infinitely thin sheet. This is commonly called the flamelet regime; it has been the subject of many papers and the basis for many models (see, e.g., Linan & Williams 1993). In the flamelet model, the local flame structure is assumed to be identical to the laminar flame structure; thus the flame propagates normal to itself at the laminar flame speed, S(sub L). This allows the use of simple approximations. For example, one expects the rate of consumption of fuel to be proportional to the area of the flame surface. This idea allowed Damkohler (1940) to propose that the wrinkled flame could be replaced by a smooth one which travels at the turbulent flame speed, S(sub T), defined by S(sub T)/S(sub L) = A(sub L) /A(sub P) where A(sub L) is the total flame surface area and AP is the area projected onto the mean direction of propagation. This relation can be expected to be valid when the flame structure is modified only slightly by the turbulence. More recent approaches have attempted to relate the turbulent flame speed to turbulence intensity, u(sub '), which presumably, characterizes the wrinkling of the flame.

Bourlioux, A.↗

Example on how to model and simulate turbulence for flight simulators

Analytical developments relative to gust response are discussed. Turbulence length scale, spectral functions, zero crossing values, gust loads analysis, power spectral techniques for analyzing the response of aircraft in turbulence, the spectrum of the rolling moment coefficient, and the spectrum correction factor are among the issues considered.

Houbolt, John C.↗

Some practical turbulence modeling options for Reynolds-averaged full Navier-Stokes calculations of three-dimensional flows

New turbulence modeling options recently implemented for the 3D version of Proteus, a Reynolds-averaged compressible Navier-Stokes code, are described. The implemented turbulence models include: the Baldwin-Lomax algebraic model, the Baldwin-Barth one-equation model, the Chien k-epsilon model, and the Launder-Sharma k-epsilon model. Features of this turbulence modeling package include: well documented and easy to use turbulence modeling options, uniform integration of turbulence models from different classes, automatic initialization of turbulence variables for calculations using one- or two-equation turbulence models, multiple solid boundaries treatment, and fully vectorized L-U solver for one- and two-equation models. Good agreements are obtained between the computational results and experimental data. Sensitivity of the compressible turbulent solutions with the method of y(+) computation, the turbulent length scale correction, and some compressibility corrections are examined in detail. Test cases show that the highly optimized one- and two-equation turbulence models can be used in routine 3D Navier-Stokes computations with no significant increase in CPU time as compared with the Baldwin-Lomax algebraic model.

Bui, Trong T.↗

Evaluation of turbulence models for three primary types of shock-separated boundary layers

Zero-equation (algebraic), one-equation (kinetic energy), and two-equation (kinetic energy plus length scale) turbulence eddy viscosity models were used in computing three basic types of shock-separated boundary-layer flows. The three basic types of shock boundary-layer interaction discussed are: (1) a normal shock wave at transonic speeds, (2) a compression corner shock at supersonic speeds, and (3) an incident oblique shock at hypersonic speeds. The models tested are simple, unmodified models used extensively for incompressible, unseparated flows. A comparison of computed and measured results for the compressible, separated flows described herein indicates that model performance is dependent on flow configuration with no distinct superiority of one model over the other for all three flow configurations.

Coakley, T. J.↗

Evaluation of Full Reynolds Stress Turbulence Models in FUN3D for Axisymmetric Jet Flow

Two full seven-equation turbulence models available in the FUN3D code are evaluated for their ability to improve the computation of challenging mixing flows encountered in aerospace propulsion. These models are the SSG/LRR and Wilcox full second-moment Reynolds stress models. They solve equations for the six components of the Reynolds stress and a seventh equation for the turbulent length scale. Two standard eddy viscosity models are also evaluated for comparison, the Spalart-Allmaras (SA) one-equation model and the Menter Shear Stress Transport (SST-V) two-equation turbulence model. Flow through an axisymmetric reference nozzle is examined at three flow conditions: subsonic unheated, subsonic heated, and near sonic unheated. Centerline profiles of velocity and turbulent kinetic energy and radial profiles of velocity, turbulent kinetic energy and turbulent stresses are examined. Results showed that the SA model did well at predicting the jet potential core length, but over-mixed the downstream flow, whereas the SST-V model over-predicted the potential core length. The Wilcox-model significantly over-predicted the potential core length and under-predicted the mixing and was not well-suited for the jet flows evaluated, however the SSG/LRR Reynolds stress model did well at predicting the mixing rate and mean velocity for all cases examined.

mixing layers↗

Some practical turbulence modeling options for Reynolds-averaged full Navier-Stokes calculations of three-dimensional flows

New turbulence modeling options recently implemented for the 3-D version of Proteus, a Reynolds-averaged compressible Navier-Stokes code, are described. The implemented turbulence models include: the Baldwin-Lomax algebraic model, the Baldwin-Barth one-equation model, the Chien k-epsilon model, and the Launder-Sharma k-epsilon model. Features of this turbulence modeling package include: well documented and easy to use turbulence modeling options, uniform integration of turbulence models from different classes, automatic initialization of turbulence variables for calculations using one- or two-equation turbulence models, multiple solid boundaries treatment, and fully vectorized L-U solver for one- and two-equation models. Validation test cases include the incompressible and compressible flat plate turbulent boundary layers, turbulent developing S-duct flow, and glancing shock wave/turbulent boundary layer interaction. Good agreement is obtained between the computational results and experimental data. Sensitivity of the compressible turbulent solutions with the method of y(sup +) computation, the turbulent length scale correction, and some compressibility corrections are examined in detail. The test cases show that the highly optimized one-and two-equation turbulence models can be used in routine 3-D Navier-Stokes computations with no significant increase in CPU time as compared with the Baldwin-Lomax algebraic model.

Bui, Trong T.↗

In Situ Wind and Turbulence Measurements in a Field of Full-Size Parabolic Trough Solar Collectors

Concentrated Solar Power (CSP) is a promising method for using Solar power for electricity generation with thermal energy storage. One of the primary drivers of structural design costs of CSP collector structures is wind loading. To date, the design of these structures has relied on data from wind tunnels that do not adequately capture the dynamic effects observed at scale. NREL initiated a field measurement campaign at the operational Nevada Solar One (NSO) powerplant that uses parabolic troughs as solar collectors. The aim of the project is a detailed characterization of prevailing wind and turbulence conditions and resulting operational loads on parabolic troughs, providing insights on structural dynamic response, and generating a first-of-a-kind, comprehensive, high-resolution wind-loading dataset available for validating simulations of wind loading on collector structures. The measurements at NSO consist of Sonic anemometers on masts at different heights to characterize the incoming flow and conditions within and above the trough field, in combination with a Doppler Lidar scanning the horizontal plane above the troughs. The wind measurements at NSO have been continuously operating since October 2021 and provide a year-long dataset characterizing wind and turbulence conditions. The structural load measurements start in November 2022 and will complement the wind measurements. In this poster, we present first results of the wind measurement campaign by highlighting days with different atmospheric flow conditions. We identified three main factors altering the flow over the parabolic troughs: Wind speed, wind direction, and the angle of the parabolic troughs. The highest loads are expected when the wind blows perpendicular to the trough rows. In this case, the first rows experience the highest loads and block the subsequent rows, creating conditions with decreased wind speed and enhanced turbulence within the trough field. Also, turbulent length scales change. These affected conditions produce unique load cases on the structures that will be captured by the load measurements. If the wind blows along the trough rows, wind and turbulence conditions are less influenced by the troughs.

atmospheric turbulence↗

In Situ Wind and Turbulence Measurements in a Field of Full-Size Parabolic Trough Collectors

Concentrated Solar Power (CSP) is a promising method for using Solar power for electricity generation with thermal energy storage. One of the primary drivers of structural design costs of CSP collector structures is wind loading. To date, the design of these structures has relied on data from wind tunnels that do not adequately capture the dynamic effects observed at scale. NREL initiated a field measurement campaign at the operational Nevada Solar One (NSO) powerplant that uses parabolic troughs as solar collectors. The aim of the project is a detailed characterization of prevailing wind and turbulence conditions and resulting operational loads on parabolic troughs, providing insights on structural dynamic response, and generating a first-of-a-kind, comprehensive, high-resolution wind-loading dataset available for validating simulations of wind loading on collector structures. The measurements at NSO consist of Sonic anemometers on masts at different heights to characterize the incoming flow and conditions within and above the trough field, in combination with a Doppler Lidar scanning the horizontal plane above the troughs. The wind measurements at NSO have been continuously operating since October 2021 and provide a year-long dataset characterizing wind and turbulence conditions. The structural load measurements start in November 2022 and will complement the wind measurements. In this poster, we present first results of the wind measurement campaign by highlighting days with different atmospheric flow conditions. We identified three main factors altering the flow over the parabolic troughs: Wind speed, wind direction, and the angle of the parabolic troughs. The highest loads are expected when the wind blows perpendicular to the trough rows. In this case, the first rows experience the highest loads and block the subsequent rows, creating conditions with decreased wind speed and enhanced turbulence within the trough field. Also, turbulent length scales change. These affected conditions produce unique load cases on the structures that will be captured by the load measurements. If the wind blows along the trough rows, wind and turbulence conditions are less influenced by the troughs.

concentrated solar power↗

A study of the effects of Rotating Frame Turbulence (RFT) on helicopter flight mechanics

The turbulence actually experienced by a helicopter blade-element significantly differs from the space-fixed free atmospheric turbulence. The turbulence in the rotor disk requires a rotationally sampled description in a rotating frame of reference. It is referred to as the rotating frame turbulence or RFT which exhibits a striking phenomenon. The RFT spectral density versus frequency shows high peak values at 1P,2P, or 3P, frequencies. The energy increase at these peaks is balanced by an energy decrease primarily at the lower-than-1P frequency range. Particularly for low altitude flight regimes of pure helicopters, such as the nap-of-the-earth maneuvers, the conventional space-fixed description of turbulence is not a good approximation, since the turbulence scale length can have values comparable to the rotor radius. Accordingly the flight mechanics characteristics with RFT description are compared with those based on the conventional space-fixed turbulence description. The results demonstrate that the RFT qualitatively and quantitatively affects the prediction of helicopter flight mechanics characteristics in turbulence. Such comparisons should play an important role in the new development of handling qualities specifications for helicopters.

Schrage, D. P.↗

Turbulence modeling of shock separated boundary-layer flows

Computations of transonic and hypersonic shock-separated boundary-layer flows using zero-equation (algebraic), one-equation (kinetic energy), and two-equation (kinetic energy plus length scale) turbulence eddy viscosity models are described and compared with measurements. The computations make use of a new Navier-Stokes computer algorithm that has reduced computing times by one to two orders of magnitude. The algorithm, and how the turbulence models are incorporated into it, are described. Results for the transonic flow show that the unmodified one-equation model is superior to the zero-equation model in skin-friction predictions. For the hypersonic flow, a highly modified one-equation model that accurately predicts surface pressure and heat transfer is described. Preliminary two-equation model results are also presented.

Coakley, T. J.↗

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

Numerical solutions of atmospheric flow over semielliptical simulated hills

Atmospheric motion over obstacles on plane surfaces to compute simulated wind fields over terrain features was studied. Semielliptical, two dimensional geometry and numerical simulation of flow over rectangular geometries is also discussed. The partial differential equations for the vorticity, stream function, turbulence kinetic energy, and turbulence length scale were solved by a finite difference technique. The mechanism of flow separation induced by a semiellipse is the same as flow over a gradually sloping surface for which the flow separation is caused by the interaction between the viscous force, the pressure force, and the turbulence level. For flow over bluff bodies, a downstream recirculation bubble is created which increases the aspect ratio and/or the turbulence level results in flow reattachment close behind the obstacle.

Shieh, C. F.↗

Effects of Freestream Turbulence on Turbine Blade Heat Transfer

Experiments have shown that moderate turbulence levels can nearly double turbine blade stagnation region heat transfer. Data have also shown that heat transfer is strongly affected by the scale of turbulence as well as its level. In addition to the stagnation region, turbulence is often seen to increase pressure surface heat transfer. This is especially evident at low to moderate Reynolds numbers. Vane and rotor stagnation region, and vane pressure surface heat transfer augmentation is often seen in a pre-transition environment. Accurate predictions of transition and relaminarization are critical to accurately predicting blade surface heat transfer. An approach is described which incorporates the effects of both turbulence level and scale into a CFD analysis. The model is derived from experimental data for cylindrical and elliptical leadng edges. Results using this model are compared to experimental data for both vane and rotor geometries. The comparisons are made to illustrate that using a model which includes the effects of turbulence length scale improves agreement with data, and to illustrate where improvements in the modeling are needed.

Boyle, Robert J.↗

Cloud Simulations in Response to Turbulence Parameterizations in the GISS Model E GCM

The response of cloud simulations to turbulence parameterizations is studied systematically using the GISS general circulation model (GCM) E2 employed in the Intergovernmental Panel on Climate Change's (IPCC) Fifth Assessment Report (AR5).Without the turbulence parameterization, the relative humidity (RH) and the low cloud cover peak unrealistically close to the surface; with the dry convection or with only the local turbulence parameterization, these two quantities improve their vertical structures, but the vertical transport of water vapor is still weak in the planetary boundary layers (PBLs); with both local and nonlocal turbulence parameterizations, the RH and low cloud cover have better vertical structures in all latitudes due to more significant vertical transport of water vapor in the PBL. The study also compares the cloud and radiation climatologies obtained from an experiment using a newer version of turbulence parameterization being developed at GISS with those obtained from the AR5 version. This newer scheme differs from the AR5 version in computing nonlocal transports, turbulent length scale, and PBL height and shows significant improvements in cloud and radiation simulations, especially over the subtropical eastern oceans and the southern oceans. The diagnosed PBL heights appear to correlate well with the low cloud distribution over oceans. This suggests that a cloud-producing scheme needs to be constructed in a framework that also takes the turbulence into consideration.

Yao, Mao-Sung↗

“Gray Zone” Simulations Using a Three-Dimensional Planetary Boundary Layer Parameterization in the Weather Research and Forecasting Model

Abstract Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δ x ) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δ x ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation. Significance Statement The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.

Juliano, Timothy W.↗

Comparison of multiequation turbulence models for several shock-separated boundary-layer interaction flows

Several multiequation eddy viscosity models of turbulence are used with the Navier-Stokes equations to compute three classes of experimentally documented shock-separated turbulent boundary-layer flows. The types of flow studied are: (1) a normal shock at transonic speeds in both a circular duct and a two-dimensional channel; (2) an incident oblique shock at supersonic speeds on a flat surface; and (3) a two-dimensional compression corner at supersonic speeds. Established zero-equation (algebraic), one-equation (kinetic energy), and two-equation (kinetic energy plus length scale) turbulence models are each utilized to describe the Reynolds shear stress for the three classes of flows. These models are assessed by comparing the calculated values of skin friction, wall pressure distribution, velocity, Mach number, and turbulent kinetic energy profiles with experimental measurements. Of the models tested the two-equation model results gave the best overall agreement with the data.

Viegas, J. R.↗

Modeling near wall effects in second moment closures by elliptic relaxation

The elliptic relaxation model of Durbin (1993) for modeling near-wall turbulence using second moment closures (SMC) is compared to DNS data for a channel flow at Re(sub t) = 395. The agreement for second order statistics and even the terms in their balance equation is quite satisfactory, confirming that very little viscous effects (via Kolmogoroff scales) need to be added to the high Reynolds versions of SMC for near-wall-turbulence. The essential near-wall feature is thus the kinematic blocking effect that a solid wall exerts on the turbulence through the fluctuating pressure, which is best modeled by an elliptic operator. Above the transition layer, the effect of the original elliptic operator decays rapidly, and it is suggested that the log-layer is better reproduced by adding a non-homogeneous reduction of the return to isotropy, the gradient of the turbulent length scale being used as a measure of the inhomogeneity of the log-layer. The elliptic operator was quite easily applied to the non-linear Craft & Launder pressure-strain model yielding an improved distinction between the spanwise and wall normal stresses, although at higher Reynolds number (Re) and away from the wall, the streamwise component is severely underpredicted, as well as the transition in the mean velocity from the log to the wake profiles. In this area a significant change of behavior was observed in the DNS pressure-strain term, entirely ignored in the models.

Laurence, D.↗