Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 199 records · Page 11

Turbulent solutions of equations of fluid motion

Some turbulent solutions of the unaveraged Navier-Stokes equations (equations of fluid motion) are reviewed. Those equations are solved numerically in order to study the nonlinear physics of incompressible turbulent flow. The three components of the mean-square velocity fluctuations are initially equal for the conditions chosen. The resulting solutions show characteristics of turbulence, such as the linear and nonlinear excitation of small-scale fluctuations. For the stronger fluctuations the initially nonrandom flow develops into an apparently random turbulence. The cases considered include turbulence that is statistically homogeneous or inhomogeneous and isotropic or anisotropic. A statistically steady-state turbulence is obtained by using a spatially periodic body force. Various turbulence processes, including the transfer of energy between eddy sizes and between directional components and the production, dissipation, and spatial diffusion of turbulence, are considered. It is concluded that the physical processes occurring in turbulence can be profitably studied numerically.

Deissler, R. G.↗

Investigation of the validity of Reynolds averaged turbulence models at the frequencies that occur in turbomachinery

Turbulent flows subjected to various kinds of unsteady disturbances were simulated using a large-eddy-simulation computer code for flow in a channel. The disturbances were: a normal velocity expressed as a traveling wave on one wall of the channel; staggered blowing and suction distributions on the opposite walls of the channel; and oscillations of the mean flow through the channel. The wall boundary conditions were designed to simulate the effects of wakes of a stator stage passing through a rotor channel in a turbine. The oscillating flow simulated the effects of a pressure pulse moving over the rotor blade boundary layer. The objective of the simulations was to provide better understanding of the effects of time-dependent disturbances on the turbulence of a boundary layer and of the underlying physical phenomena regarding the basic interaction between the turbulence and external disturbances of the type found in turbomachinery. Results showed that turbulence is sensitive to certain ranges of frequencies of disturbances. However, no direct connection was found between the frequency of imposed disturbances and characteristic burst frequency of turbulence. New insight into the nature of turbulence at high frequencies was found. The viscous phenomena near solid walls was found to be the dominant influence for high frequency perturbations. At high frequencies, the turbulence was found to be undisturbed, remaining the same as for the steady mean flow. A transition range exists between the high frequency range and the low, or quasi-steady, range in which the turbulence is not predictable by either quasi-steady models or the steady flow model. The limiting lowest frequency for use of the steady flow turbulence model is that for which the viscous Stokes layer based on the blade passing frequency is thicker than the laminar sublayer.

Kuhn, Gary D.↗

An investigation of turbulence models

The accuracy to which a turbulent boundary layer or wake can be predicted numerically depends on the validity of the turbulence closure model used. The modeling of turbulence physics is one of the most difficult problems in computational fluid dynamics (CFD). In fact, it is one of the pacing factors in the development of CFD. In general, there are three main approaches to the description of trubulence physics. First is turbulence modeling in which the Reynolds averaged Navier-Stokes equations are used and some closure approximation is made for the Reynolds stresses. A second approach to turbulence is large eddy simulation (LES) in which the computational mesh is taken to be fine enough that the large scale structure of the turbulence can be calculated directly. An empirical assumption must be made for the small scale sub-grid turbulence. The third approach is direct simulation. In this technique the Navier-Strokes equations are solved directly on a mesh which if fine enough to resolve the smallest length scale of the turbulence. The Reynolds averaged equations are not used and no closure assumption is required. These last two approaches require extensive computer resources and as such are not engineering tools. The purpose of the work was to investigate the various engineering turbulence models for accuracy and ease of programming. This involved comparison of the models with each other and with experimental data.

Doria, Michael L.↗

Calculation of reattaching shear layers in divergent channel with a multiple-time-scale turbulence model

Numerical calculations of turbulent reattaching shear layers in a divergent channel are presented. The turbulence is described by a multiple-time-scale turbulence model. The turbulent flow equations are solved by a control-volume based finite difference method. The computational results are compared with those obtained using k-epsilon turbulence models and algebraic Reynolds stress turbulence models. It is shown that the multiple-time-scale turbulence model yields significantly improved computational results than the other turbulence models in the region where the turbulence is in a strongly inequilibrium state.

Kim, S.-W.↗

Calculation of reattaching shear layers in divergent channel with a multiple-time-scale turbulence model

Numerical calculations of turbulent reattaching shear layers in a divergent channel are presented. The turbulence is described by a multiple-time-scale turbulence model. The turbulent flow equations are solved by a control-volume based finite difference method. The computational results are compared with those obtained using k-epsilon turbulence models and algebraic Reynolds stress turbulence models. It is shown that the multiple-time-scale turbulence model yields significantly improved computational results than the other turbulence models in the region where the turbulence is in a strongly inequilibrium state.

Kim, S.-W.↗

Investigation of turbulent transport in an axisymmetric sudden expansion

Simultaneous two-component laser velocimeter measurements were made in the incompressible turbulent flowfield following an axisymmetric sudden expansion. Mean velocities, Reynolds stresses, and triple products were measured and are presented at axial positions ranging from x/H = 0.2-14. A balance of the turbulent kinetic energy in the flow was performed. The production, convection, and diffusion of turbulent kinetic energy were computed directly from the experimental data using central differencing. A specially designed correction lens was employed to correct for optical aberrations introduced by the circular tube. This lens system allowed the accurate simultaneous measurement of axial and radial velocities in the test section. The experimental measurements were compared to predictions generated by a code that employed the k-epsilon turbulence model. Agreement was good for mean axial velocities, turbulent kinetic energy, and turbulent shear stresses. However, the modeled turbulent normal stresses where in poor agreement with the measured values. The modeled diffusion of turbulent kinetic energy was underpredicted in the region between the shear layer and the centerline of the flow giving lower values of turbulent kinetic energy downstream of the potential core than measured.

Gould, Richard D.↗

RNG in turbulence and modeling of bypass transition

Two projects are considered: the Renormalization Group (RNG) analysis of turbulence modeling, and the calculation of bypass transition through turbulence modeling. RNG is a process which eliminates small scales on the uneliminated large scales as the change in the transport properties. It is because of this property of RNG that it was previously suggested that RNG could be used as a model builder in turbulence modeling. The possibility is studied of constructing RNG based turbulence models, and to try to proceed to do the modeling through RNG in parallel with the classical approach. The numerical predictions made by RNG models and by classical models is compared against data from Direct Numerical Simulation. While in an environment with freestream turbulence, the transition initiated by the instability of the laminar boundary layer to Tollmien-Schlichting waves is found to be a bypass one in which turbulent spots are formed without T-S wave amplification. The formation is a random process, and flow within a turbulent spot is almost fully turbulent. This suggests the possibility of using turbulence modeling to describe and predict the bypass transition.

Yang, Zhigang↗

Some issues on modeling atmospheric turbulence experienced by helicopter rotor blades

The atmospheric turbulence velocities seen by nonrotating aircraft components and rotating blades can be substantially different. The differences are due to the spatial motion of the rotor blades, which move fore and aft through the gust waves. Body-fixed atmospheric turbulence refers to the actual atmospheric turbulence experienced by a point fixed on a nonrotating aircraft component such as the aircraft's center of gravity or the rotor hub, while blade-fixed atmospheric turbulence refers to the atmospheric turbulence experienced by an element of the rotating rotor blade. An example is presented, which, though overly simplified, shows important differences between blade- and body-fixed rotorcraft atmospheric turbulence models. All of the information necessary to develop the dynamic equations describing the atmospheric turbulence velocity field experienced by an aircraft is contained in the atmospheric turbulence velocity correlation matrix. It is for this reason that a generalized formulation of the correlation matrix describing atmospheric turbulence that a rotating blade encounters is developed. From this correlation matrix, earlier treated cases restricted to a rotor flying straight and level directly into the mean wind can be recovered as special cases.

Costello, Mark↗

Second order closure modeling of turbulent buoyant wall plumes

Non-intrusive measurements of scalar and momentum transport in turbulent wall plumes, using a combined technique of laser Doppler anemometry and laser-induced fluorescence, has shown some interesting features not present in the free jet or plumes. First, buoyancy-generation of turbulence is shown to be important throughout the flow field. Combined with low-Reynolds-number turbulence and near-wall effect, this may raise the anisotropic turbulence structure beyond the prediction of eddy-viscosity models. Second, the transverse scalar fluxes do not correspond only to the mean scalar gradients, as would be expected from gradient-diffusion modeling. Third, higher-order velocity-scalar correlations which describe turbulent transport phenomena could not be predicted using simple turbulence models. A second-order closure simulation of turbulent adiabatic wall plumes, taking into account the recent progress in scalar transport, near-wall effect and buoyancy, is reported in the current study to compare with the non-intrusive measurements. In spite of the small velocity scale of the wall plumes, the results showed that low-Reynolds-number correction is not critically important to predict the adiabatic cases tested and cannot be applied beyond the maximum velocity location. The mean and turbulent velocity profiles are very closely predicted by the second-order closure models. but the scalar field is less satisfactory, with the scalar fluctuation level underpredicted. Strong intermittency of the low-Reynolds-number flow field is suspected of these discrepancies. The trends in second- and third-order velocity-scalar correlations, which describe turbulent transport phenomena, are also predicted in general, with the cross-streamwise correlations better than the streamwise one. Buoyancy terms modeling the pressure-correlation are shown to improve the prediction slightly. The effects of equilibrium time-scale ratio and boundary condition are also discussed.

Zhu, Gang↗

Large-eddy simulation of compressible turbulence

The increase in the range of length scales with increasing Reynolds number limits the direct simulation of turbulent flows to relatively simple geometries and low Reynolds numbers. However, since most flows of engineering interest occur at much higher Reynolds number than is currently within the capabilities of full simulation, prediction of these flow fields can only be obtained by solving some suitably-averaged set of governing equations. In the traditional Reynolds-averaged approach, the Navier-Stokes equations are averaged over time. This in turn yields correlations between various turbulence fluctuations. It is these terms, e.g. the Reynolds stresses, for which a turbulence model must be derived. Turbulence modeling of incompressible flows has received a great amount of attention in the literature. An area of research that has received comparatively less attention is the modeling of compressible turbulent flows. An approach to simulating compressible turbulence at high Reynolds numbers is through the use of Large-Eddy Simulation (LES). In LES the dependent variables are decomposed into a large-scale (resolved) component and a sub-grid scale component. It is the small-scale components of the velocity field which are presumably more homogeneous than the large scales and, therefore, more easily modeled. Thus, it seems plausible that simpler models, which should be more universal in character than those employed in second-order closure schemes, may be developed for LES of compressible turbulence. The objective of the present research, therefore, is to explore models for the Large-Eddy Simulation of compressible turbulent flows. Given the recent successes of Zeman in second order closure modeling of compressible turbulence, model development was guided by principals employed in second-order closures.

Squires, Kyle D.↗

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