Search NASA⌕ Search

SEARCH · Search NASA

Results for “TURBULENT”

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 289 records · Page 16

Investigations of turbulent scalar fields using probability density function approach

Scalar fields undergoing random advection have attracted much attention from researchers in both the theoretical and practical sectors. Research interest spans from the study of the small scale structures of turbulent scalar fields to the modeling and simulations of turbulent reacting flows. The probability density function (PDF) method is an effective tool in the study of turbulent scalar fields, especially for those which involve chemical reactions. It has been argued that a one-point, joint PDF approach is the one to choose from among many simulation and closure methods for turbulent combustion and chemically reacting flows based on its practical feasibility in the foreseeable future for multiple reactants. Instead of the multi-point PDF, the joint PDF of a scalar and its gradient which represents the roles of both scalar and scalar diffusion is introduced. A proper closure model for the molecular diffusion term in the PDF equation is investigated. Another direction in this research is to study the mapping closure method that has been recently proposed to deal with the PDF's in turbulent fields. This method seems to have captured the physics correctly when applied to diffusion problems. However, if the turbulent stretching is included, the amplitude mapping has to be supplemented by either adjusting the parameters representing turbulent stretching at each time step or by introducing the coordinate mapping. This technique is still under development and seems to be quite promising. The final objective of this project is to understand some fundamental properties of the turbulent scalar fields and to develop practical numerical schemes that are capable of handling turbulent reacting flows.

Gao, Feng↗

Research activities at the Center for Modeling of Turbulence and Transition

The main research activities at the Center for Modeling of Turbulence and Transition (CMOTT) are described. The research objective of CMOTT is to improve and/or develop turbulence and transition models for propulsion systems. The flows of interest in propulsion systems can be both compressible and incompressible, three dimensional, bounded by complex wall geometries, chemically reacting, and involve 'bypass' transition. The most relevant turbulence and transition models for the above flows are one- and two-equation eddy viscosity models, Reynolds stress algebraic- and transport-equation models, pdf models, and multiple-scale models. All these models are classified as one-point closure schemes since only one-point (in time and space) turbulent correlations, such as second moments (Reynolds stresses and turbulent heat fluxes) and third moments, are involved. In computational fluid dynamics, all turbulent quantities are one-point correlations. Therefore, the study of one-point turbulent closure schemes is the focus of our turbulence research. However, other research, such as the renormalization group theory, the direct interaction approximation method, and numerical simulations are also pursued to support the development of turbulence modeling.

Shih, Tsan-Hsing↗

A new approach to the study of the turbulence layer

Turbulent thermal convection is of considerable importance in fluid dynamical transport phenomena occurring, for example, in the planetary boundary layer of the Earth, the interiors of stars, and accretion disks. In particular, during a significant portion of the evolutionary phase of many stars having convectively unstable cores or outer envelopes, a substantial fraction of energy is transported from the central layers to the outer layers by thermal convection. Moreover, as much of the interior of a star is in highly turbulent motion, a complete theory of stellar structure and evolution requires the explicit consideration of turbulence in order to have expressions for the turbulent quantities arising in the stellar structure equations, and particularly, the turbulent fluxes that appear in the total flux conservation equation, such as the convective flux, kinetic energy flux, etc. A reliable quantification of these fluxes continues to present a challenge in astrophysical fluid dynamics, primarily because astrophysical turbulence is almost always fully-developed and nearly inviscid, and therefore governed by strong nonlinear interactions that distribute the energy among a very wide spectrum of eddies with scales ranging from the characteristic dimension of the flow to those sufficiently small to be affected by viscous dissipation. Furthermore, astrophysical flows are invariably compressible, anisotropic, and inhomogeneous, which requires the consideration of the dynamics of longitudinal modes and their interaction with the transverse modes, as well as complicated boundary conditions. In order to reach a compromise between analytical and numerical tractability and the basic physics of turbulent convection, we have constructed a model of stationary turbulent convection that yields various turbulence statistics, including the convective flux, that are required in stellar evolution models.

Levenson, Ruth A.↗

An application of a two-equation model of turbulence to three-dimensional chemically reacting flows

A numerical study of three dimensional chemically reacting and non-reacting flowfields is conducted using a two-equation model of turbulence. A generalized flow solver using an implicit Lower-Upper (LU) diagonal decomposition numerical technique and finite-rate chemistry has been coupled with a low-Reynolds number two-equation model of turbulence. This flow solver is then used to study chemically reacting turbulent supersonic flows inside combustors with synergetic fuel injectors. The reacting and non-reacting turbulent combustor solutions obtained are compared with zero-equation turbulence model solutions and with available experimental data. The hydrogen-air chemistry is modeled using a nine-species/eighteen reaction model. A low-Reynolds number k-epsilon model was used to model the effect of turbulence because, in general, the low-Reynolds number k-epsilon models are easier to implement numerically and are far more general than algebraic models. However, low-Reynolds number k-epsilon models require a much finer near-wall grid resolution than high-Reynolds number models to resolve accurately the near-wall physics. This is especially true in complex flowfields, where the stiff nature of the near-wall turbulence must be resolved. Therefore, the limitations imposed by the near-wall characteristics and compressible model corrections need to be evaluated further. The gradient-diffusion hypothesis is used to model the effects of turbulence on the mass diffusion process. The influence of this low-Reynolds number turbulence model on the reacting flowfield predictions was studied parametrically.

Lee, J.↗

A new path for the electron bulk energization in solar flares: Fermi acceleration by magnetohydrodynamic turbulence in reconnection outflows

We recently proposed that a magnetohydrodynamic (MHD) turbulent cascade produces the bulk energization of electrons to approximately 25 keV in the impulsive phase of solar flares (LaRosa & Moore 1993). In that scenario, (1) the cascading MHD turbulence is fed by shear-unstable Alfvenic outflows from sites of strongly driven reconnection in the low corona, and (2) the electrons are energized by absorbing the energy that flows down through the cascade. We did not specify the physical mechanism by which the cascading energy is ultimately transferred to the electrons. Here we propose that Fermi acceleration is this mechanism, the process by which the electrons are energized and by which the cascading MHD turbulence is dissipated. We point out that in the expected cascade MHD fluctuations of scale 1 km can Fermi-accelerate electrons from 0.1 keV to approximately 25 keV on the subsecond timescales observed in impulsive flares, provided there is sufficient trapping and scattering of electrons in the MHD turbulence. We show that these same fluctuations provide the required trapping; they confine the electrons within the turbulent region until the turbulence eis dissipated. This results in the energization of all of the lectrons in each large-scale (5 x 10(exp 7)cm) turbulent eddy to 25 keV. The Fermi process also requires efficient scattering so that the pitch-angle distribution of the accelerating electrons remains isotropic. We propose that the electrons undergo resonant scattering by high-frequency plasma R-waves that, as suggested by others (Hamilton & Petrosian 1992), are generated by the reconnection. Ions are not scattered by R-waves. Provided that there is negligible generation of ion-scattering plasma turbulence (e.g., L-waves) by the reconnection or the MHD turbulence, the ions will not Fermi-accelerate and the cascading energy is transferred only to the electrons. We conclude that, given this situation, electron Fermi acceleration can plausibly account for the electron bulk energization in impulsive solar flares.

Larosa, T. N.↗

Effects of pressure gradients on turbulent premixed flames

The influence of a constant acceleration on a turbulent premixed flame is studied by direct numerical simulation. This acceleration induces a mean pressure gradient across the flame brush, leading to a modification of the turbulent flame structure due to differential buoyancy mechanisms between heavy cold fresh and light hot burnt gases. Such a pressure gradient may be encountered in practical applications in ducted flames. A favorable pressure gradient, i.e. the pressure decreases from unburnt to burnt gases, is found to decrease the flame wrinkling, the flame brush thickness, and the turbulent flame speed. A favorable pressure gradient also promotes counter-gradient turbulent transport. On the other hand, adverse pressure gradients tend to increase the flame brush thickness and turbulent flame speed, and promote classical gradient turbulent transport. The balance equation for the turbulent flux of the Favre averaged progress variable is also analyzed. The first results show that the fluctuating pressure term, cannot be neglected as generally assumed in models. Simple models assuming that a high mean pressure gradient may only be balanced by the cross-dissipation term seem too approximate. This analysis has to be continued to compare simulation data and closure schemes proposed for the transport equation. The analysis developed by Veynante et al.(1995) has been extended to imposed acceleration and mean pressure gradients. A simple model for the turbulent flux is proposed and validated from simulation data. Then, a modified criterion is derived to delineate between counter-gradient and gradient turbulent diffusion. In fact, counter-gradient diffusion may occur in most practical applications, especially for ducted flames.

Veynante, D.↗

Turbulence Modeling and Computation of Turbine Aerodynamics and Heat Transfer

The objective of the present research is to develop improved turbulence models for the computation of complex flows through turbomachinery passages, including the effects of streamline curvature, heat transfer and secondary flows. Advanced turbulence models are crucial for accurate prediction of rocket engine flows, due to existance of very large extra strain rates, such as strong streamline curvature. Numerical simulation of the turbulent flows in strongly curved ducts, including two 180-deg ducts, one 90-deg duct and a strongly concave curved turbulent boundary layer have been carried out with Reynolds stress models (RSM) and algebraic Reynolds stress models (ARSM). An improved near-wall pressure-strain correlation has been developed for capturing the anisotropy of turbulence in the concave region. A comparative study of two modes of transition in gas turbine, the by-pass transition and the separation-induced transition, has been carried out with several representative low-Reynolds number (LRN) k-epsilon models. Effects of blade surface pressure gradient, freestream turbulence and Reynolds number on the blade boundary layer development, and particularly the inception of transition are examined in detail. The present study indicates that the turbine blade transition, in the presence of high freestream turbulence, is predicted well with LRN k-epsilon models employed. The three-dimensional Navier-Stokes procedure developed by the present authors has been used to compute the three-dimensional viscous flow through the turbine nozzle passage of a single stage turbine. A low Reynolds number k-epsilon model and a zonal k-epsilon/ARSM (algebraic Reynolds stress model) are utilized for turbulence closure. An assessment of the performance of the turbulence models has been carried out. The two models are found to provide similar predictions for the mean flow parameters, although slight improvement in the prediction of some secondary flow quantities has been obtained by the ARSM model. It's found that the wake profiles inside the endwall boundary layers are predicted better than those near the mid-span.

Lakshminarayana, B.↗

An intercomparison of plasma turbulence at three comets: Grigg-Skjellerup, Giacobini-Zinner, and Halley

We examine and intercompare the LF plasma wave turbulence at three comets: Grigg-Skjellerup (GS), Giacobini-Zinner (GZ), and Halley (H). All three have power spectral peaks at the local ion cyclotron frequency (the pump wave) at approx. 10(exp -2) Hz, and a power-law fall-off at higher frequencies that suggest the development of turbulent cascades. The power laws for the three comets are approximately f(exp -1.9), f(exp -1.9) and f(exp -2.1), respectively. However, other than the similarities in the power spectra, we find the magnetic field turbulence is considerably different at the three comets. Phase steepening is demonstrated to occur at the trailing edges of the GS waves. This is probably due to nonlinear steepening plus dispersion of the left-hand mode components. A coherency analysis of GZ turbulence indicates that it is primarily composed of righthanded mode components, i.e., the turbulence is 'whistlermode.' This too can be explained by nonlinear steepening plus dispersion of the magnetosonic waves. At the level of GS and GZ turbulence development when the spacecraft measurements were made, classical three-wave processes, such as the decay or modulation instabilities do not appear to play important roles. It is most likely that the nonlinear steepening and dispersive time scales are more rapid than three-wave processes, and the latter had not had time to develop for the relatively 'new' turbulence. The wave turbulence at Halley is linearly polarized. The exact nature of this turbulence is still not well understood at this time. Several possibilities are suggested, based on our preliminary analyses.

Tsurutani, Bruce T.↗

Toward a Turbulence Constitutive Relation for Rotating Flows

In rapidly rotating turbulent flows the largest scales of the motion are in approximate geostrophic balance. Single-point turbulence closures, in general, cannot attain a geostrophic balance. This article addresses and resolves the possibility of constitutive relation procedures for single-point second order closures for a specific class of rotating or stratified flows. Physical situations in which the geostrophic balance is attained are described. Closely related issues of frame-indifference, horizontal nondivergence, Taylor-Proudman theorem and two-dimensionality are, in the context of both the instantaneous and averaged equations, discussed. It is shown, in the absence of vortex stretching along the axis of rotation, that turbulence is frame-indifferent. A derivation and discussion of a geostrophic constraint which the prognostic equations for second-order statistics must satisfy for turbulence approaching a frame-indifferent limit is given. These flow situations, which include rotating and nonrotating stratified flows, are slowly evolving flows in which the constitutive relation procedures are useful. A nonlinear non-constant coefficient representation for the rapid-pressure strain covariance appearing in the Reynolds stress and heat flux equations consistent with the geostrophic balance is described. The rapid-pressure strain model coefficients are not constants determined by numerical optimization but are functions of the state of the turbulence as parameterized by the Reynolds stresses and the turbulent heat fluxes. The functions are valid for all states of the turbulence attaining their limiting values only when a limit state is achieved. These issues are relevant to strongly vortical flows as well as flows such as the planetary boundary layers, in which there is a transition from a three-dimensional shear driven turbulence to a geostrophic or horizontal turbulence.

Ristorcelli, J. R.↗

Large Eddy Simulation of Aircraft Wake Vortices: Atmospheric Turbulence Effects

Crow instability can develop in most atmospheric turbulence levels, however, the ring vortices may not form in extremely strong turbulence cases due to strong dissipation of the vortices. It appears that strong turbulence tends to accelerate the occurrences of Crow instability. The wavelength of the most unstable mode is estimated to be about 5b(sub 0), which is less than the theoretical value of 8.6b(sub 0) (Crow, 1970) and may be due to limited domain size and highly nonlinear turbulent flow characteristics. Three-dimensional turbulence can decay wake vortices more rapidly. Axial velocity may be developed by vertical distortion of a vortex pair due to Crow instability or large turbulent eddy motion. More experiments with various non-dimensional turbulence levels are necessary to get useful statistics of wake vortex behavior due to turbulence. Need to investigate larger turbulence length scale effects by enlarging domain size or using grid nesting.

Han, Jongil↗

Studying Turbulence Using Numerical Simulation Databases: Proceedings of the Summer Program - No. 7

The Seventh Summer Program of the Center for Turbulence Research took place in the four-week period, July 5 to July 31, 1998. This was the largest CTR Summer Program to date, involving thirty-six participants from the U. S. and nine other countries. Thirty-one Stanford and NASA-Ames staff members facilitated and contributed to most of the Summer projects. A new feature, and perhaps a preview of the future programs, was that many of the projects were executed on non-NASA computers. These included supercomputers located in Europe as well as those operated by the Departments of Defense and Energy in the United States. In addition, several simulation programs developed by the visiting participants at their home institutions were used. Another new feature was the prevalence of lap-top personal computers which were used by several participants to carry out some of the work that in the past were performed on desk-top workstations. We expect these trends to continue as computing power is enhanced and as more researchers (many of whom CTR alumni) use numerical simulations to study turbulent flows. CTR's main role continues to be in providing a forum for the study of turbulence for engineering analysis and in facilitating intellectual exchange among the leading researchers in the field. Once again the combustion group was the largest. Turbulent combustion has enjoyed remarkable progress in using simulations to address increasingly complex and practically more relevant questions. The combustion group's studies included such challenging topics as fuel evaporation, soot chemistry, and thermonuclear reactions. The latter study was one of three projects related to the Department of Energy's ASCI Program (www.llnl.gov/asci); the other two (rocket propulsion and fire safety) were carried out in the turbulence modeling group. The flow control and acoustics group demonstrated a successful application of the so-called evolution algorithms which actually led to a previously unknown forcing strategy for jets yielding increased spreading rate. A very efficient algorithm for flow in complex geometries with moving boundaries based on the immersed boundary forcing technique was tested with very encouraging results. Also a new strategy for the destruction of aircraft trailing vortices was introduced and tested. The Reynolds Averaged Modeling (RANS) group demonstrated that the elliptic relaxation concept for RANS calculations is also applicable to transonic flows with shocks; however, prediction of laminar/turbulent transition remains an important pacing item. A large fraction of the LES effort was devoted to the development and testing of a new algorithmic procedure (as opposed to phenomenological model) for subgrid scale modeling based on regularized de-filtering of the flow variables. This appears to be a very promising approach, and a significant effort is currently underway to assess its robustness in high Reynolds number flows and in conjunction with numerical methods for complex flows. As part of the Summer Program two review tutorials were given on Turbulent structures in hydrocarbon pool fires (Sheldon Tieszen), and Turbulent combustion modeling: from RANS to LES via DNS (Luc Vervisch); and two seminars entitled Assessment of turbulence models for engineering applications (Paul Durbin) and Subgrid-scale modeling for non-premixed, turbulent reacting flows (James Riley) were presented. A number of colleagues from universities, government agencies, and industry attended the final presentations of the participants on July 31 and participated in the discussions. There are twenty-six papers in this volume grouped in five areas. Each group is preceded with an overview by its coordinator.

Source record↗

Large Eddy Simulation of Homogeneous Rotating Turbulence

Study of turbulent flows in rotating reference frames has proven to be one of the more challenging areas of turbulence research. The large number of theoretical, experimental, and computational studies performed over the years have demonstrated that the effect of solid-body rotation on turbulent flows is subtle and remains exceedingly difficult to predict. Because of the complexities associated with non-homogeneous turbulence, it is worthwhile to examine the effect of steady system rotation on the evolution of an initially isotropic turbulent flow. The assumption of statistical homogeneity considerably simplifies analysis and computation; calculation of homogeneous turbulence is further motivated since it possesses the essential physics found in more complex rotating flows. The principal objectives of the present study have therefore been to increase our fundamental understanding of turbulent flows in rotating reference frames through an examination of the asymptotic state of homogeneous rotating turbulence; particularly as to the existence of an asymptotic state which is self similar. Knowledge of an asymptotic similarity state permits prediction of the ultimate statistical evolution of the flow without requiring detailed knowledge of the complex, and not well understood, non-linear transfer processes. Aside from examination of possible similarity states in rotating turbulence, of further interest in this study has been an examination of the degree to which solid-body rotation induces a two-dimensional state in an initially isotropic flow.

Squires, Kyle D.↗

Interaction of Particles and Turbulence in the Solar Nebula

The most widely accepted theories for the formation of the Solar system claim that small solid particles continue to settle into a thin layer at the midplane of the Solar nebula until it becomes gravitationally unstable and collapses directly into km-sized planetesimals. This scenario has been challenged on at least two grounds: (1) due to turbulence, the particles may not settle into a thin layer, and (2) a thin layer may not be unstable. The Solar nebula contains at least three sources of turbulence: radial shear, vertical shear, and thermal convection. The first of these is small and probably negligible, while the last is poorly understood. However, the second contribution is likely to be substantial. The particle-rich layer rotates at nearly the Keplerian speed, but the surrounding gaseous nebula rotates slower because it is partly supported by pressure. The resulting shear generates a turbulent boundary layer which stirs the particles away from the midplane, and forestalls gravitational instability. Our previous work used a 'zero-equation' (Prandtl) model to predict the intensity of shear-generated turbulence, and enabled us to demonstrate numerically that settling of particles to the midplane is self-limiting. However, we neglected the possibility that mass loading by particles might damp the turbulence. To explore this, we have developed a more sophisticated 'one-equation' model which incorporates local generation, transport, and dissipation of turbulence, as well as explicit damping of turbulence by particles. We also include a background level of global turbulence to represent other sources. Our results indicate that damping flattens the distribution of particles somewhat, but that background turbulence thickens the particle layer.

Dacles-Mariani, Jennifer S.↗

Large Eddy Simulation of Aircraft Wake Vortices in a Homogeneous Atmospheric Turbulence: Vortex Decay and Descent

The effects of ambient turbulence on decay and descent of aircraft wake vortices are studied using a validated, three-dimensional: large-eddy simulation model. Numerical simulations are performed in order to isolate the effect of ambient turbulence on the wake vortex decay rate within a neutrally-stratified atmosphere. Simulations are conducted for a range of turbulence intensities, by injecting wake vortex pairs into an approximately homogeneous and isotropic turbulence field. The decay rate of the vortex circulation increases clearly with increasing ambient turbulence level, which is consistent with field observations. Based on the results from the numerical simulations, simple decay models are proposed as functions of dimensionless ambient turbulence intensity (eta) and dimensionless time (T) for the circulation averaged over a range of radial distances. With good agreement with the numerical results, a Gaussian type of vortex decay model is proposed for weak turbulence: while an exponential type of Tortex decay model can be applied for strong turbulence. A relationship for the vortex descent based on above vortex decay model is also proposed. Although the proposed models are based on simulations assuming neutral stratification, the model predictions are compared to Lidar vortex measurements observed during stable, neutral, and unstable atmospheric conditions. In the neutral and unstable atmosphere, the model predictions appear to be in reasonable agreement with the observational data, while in the stably-stratified atmosphere, they largely underestimate the observed circulation decay with consistent overestimation of the observed vortex descent. The underestimation of vortex decay during stably-stratified conditions suggests that stratification has an important influence on vortex decay when ambient levels of turbulence are weak.

Han, Jongil↗

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

Progress in Development of an Airborne Turbulence Detection System

Aircraft encounters with turbulence are the leading cause of in-flight injuries (Tyrvanas 2003) and have occasionally resulted in passenger and crew fatalities. Most of these injuries are caused by sudden and unexpected encounters with severe turbulence in and around convective activity (Kaplan et al 2005). To alleviate this problem, the Turbulence Prediction and Warning Systems (TPAWS) element of NASA s Aviation Safety program has investigated technologies to detect and warn of hazardous in-flight turbulence. This effort has required the numerical modeling of atmospheric convection: 1) for characterizing convectively induced turbulence (CIT) environments, 2) for defining turbulence hazard metrics, and 3) as a means of providing realistic three-dimensional data sets that can be used to test and evaluate turbulence detection sensors. The data sets are being made available to industry and the FAA for certification of future airborne turbulence-detection systems (ATDS) with warning capability. Early in the TPAWS project, a radar-based ATDS was installed and flight tested on NASA s research aircraft, a B-757. This ATDS utilized new algorithms and hazard metrics that were developed for use with existing airborne predictive windshear radars, thus avoiding the installation of new hardware. This system was designed to detect and warn of hazardous CIT even in regions with weak radar reflectivity (i.e. 5-15 dBz). Results from an initial flight test of the ATDS were discussed in Hamilton and Proctor (2002a; 2002b). In companion papers (Proctor et al 2002a; 2002b), a numerical simulation of the most significant encounter from that flight test was presented. Since the presentation of these papers a second flight test has been conducted providing additional cases for examination. In this paper, we will present results from NASA s flight test and a numerical model simulation of a turbulence environment encountered on 30 April 2002. Progress leading towards FAA certification of industry built ATDS will also be discussed.

Hamilton, David W.↗

Turbulence Model Behavior in Low Reynolds Number Regions of Aerodynamic Flowfields

The behaviors of the widely-used Spalart-Allmaras (SA) and Menter shear-stress transport (SST) turbulence models at low Reynolds numbers and under conditions conducive to relaminarization are documented. The flows used in the investigation include 2-D zero pressure gradient flow over a flat plate from subsonic to hypersonic Mach numbers, 2-D airfoil flow from subsonic to supersonic Mach numbers, 2-D subsonic sink-flow, and 3-D subsonic flow over an infinite swept wing (particularly its leading-edge region). Both models exhibit a range over which they behave transitionally in the sense that the flow is neither laminar nor fully turbulent, but these behaviors are different: the SST model typically has a well-defined transition location, whereas the SA model does not. Both models are predisposed to delayed activation of turbulence with increasing freestream Mach number. Also, both models can be made to achieve earlier activation of turbulence by increasing their freestream levels, but too high a level can disturb the turbulent solution behavior. The technique of maintaining freestream levels of turbulence without decay in the SST model, introduced elsewhere, is shown here to be useful in reducing grid-dependence of the model's transitional behavior. Both models are demonstrated to be incapable of predicting relaminarization; eddy viscosities remain weakly turbulent in accelerating or laterally-strained boundary layers for which experiment and direct simulations indicate turbulence suppression. The main conclusion is that these models are intended for fully turbulent high Reynolds number computations, and using them for transitional (e.g., low Reynolds number) or relaminarizing flows is not appropriate.

Rumsey, Christopher L.↗

Turbulence Model Behavior in Low Reynolds Number Regions of Aerodynamic Flowfields

The behaviors of the widely-used Spalart-Allmaras (SA) and Menter shear-stress transport (SST) turbulence models at low Reynolds numbers and under conditions conducive to relaminarization are documented. The flows used in the investigation include 2-D zero pressure gradient flow over a flat plate from subsonic to hypersonic Mach numbers, 2-D airfoil flow from subsonic to supersonic Mach numbers, 2-D subsonic sink-flow, and 3-D subsonic flow over an infinite swept wing (particularly its leading-edge region). Both models exhibit a range over which they behave 'transitionally' in the sense that the flow is neither laminar nor fully turbulent, but these behaviors are different: the SST model typically has a well-defined transition location, whereas the SA model does not. Both models are predisposed to delayed activation of turbulence with increasing freestream Mach number. Also, both models can be made to achieve earlier activation of turbulence by increasing their freestream levels, but too high a level can disturb the turbulent solution behavior. The technique of maintaining freestream levels of turbulence without decay in the SST model, introduced elsewhere, is shown here to be useful in reducing grid-dependence of the model's transitional behavior. Both models are demonstrated to be incapable of predicting relaminarization; eddy viscosities remain weakly turbulent in accelerating or laterally-strained boundary layers for which experiment and direct simulations indicate turbulence suppression. The main conclusion is that these models are intended for fully turbulent high Reynolds number computations, and using them for transitional (e.g., low Reynolds number) or relaminarizing flows is not appropriate.

Rumsey, Christopher L.↗