Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two-equation turbulence models”

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 37 records · Page 2

Numerical techniques for the solution of the compressible Navier-Stokes equations and implementation of turbulence models

The time-splitting explicit numerical method of MacCormack is applied to separated turbulent boundary layer flow problems. Modifications of this basic method are developed to counter difficulties associated with complicated geometry and severe numerical resolution requirements of turbulence model equations. The accuracy of solutions is investigated by comparison with exact solutions for several simple cases. Procedures are developed for modifying the basic method to improve the accuracy. Numerical solutions of high-Reynolds-number separated flows over an airfoil and shock-separated flows over a flat plate are obtained. A simple mixing length model of turbulence is used for the transonic flow past an airfoil. A nonorthogonal mesh of arbitrary configuration facilitates the description of the flow field. For the simpler geometry associated with the flat plate, a rectangular mesh is used, and solutions are obtained based on a two-equation differential model of turbulence.

Baldwin, B. S.↗

Non-equilibrium radiation from viscous chemically reacting two-phase exhaust plumes

A knowledge of the structure of the rocket exhaust plumes is necessary to solve problems involving plume signatures, base heating, plume/surface interactions, etc. An algorithm is presented which treats the viscous flow of multiphase chemically reacting fluids in a two-dimensional or axisymmetric supersonic flow field. The gas-particle flow solution is fully coupled with the chemical kinetics calculated using an implicit scheme to calculate chemical production rates. Viscous effects include chemical species diffusion with the viscosity coefficient calculated using a two-equation turbulent kinetic energy model.

Penny, M. M.↗

A reattaching free shear layer in compressible turbulent flow - A comparison of numerical and experimental results

An investigation of a two-dimensional, free turbulent shear layer reattaching on an inclined surface at Mach 2.92 and at a high Reynolds number is described. The test geometry is specifically designed to isolate the reattachment process of a high-speed separated flow. A numerical solution of the time-dependent, Reynolds-averaged, Navier-Stokes equations for the entire flow field, employing a two-equation eddy viscosity turbulence model, is presented. Detailed comparisons of prediction and experiment are made in the free shear layer, at reattachment, and in the developing boundary layer downstream. These comparisons include mean surface quantities as well as mean and fluctuating flowfield quantities. Although the overall features of this complex flow field are predicted, there are several deficiencies in the numerical solution, particularly in the region downstream of reattachment. Modifications of the turbulence model to correct these deficiencies are discussed.

Horstman, C. C.↗

Atmospheric flow over two-dimensional bluff surface obstructions

The phenomenon of atmospheric flow over a two-dimensional surface obstruction, such as a building (modeled as a rectangular block, a fence or a forward-facing step), is analyzed by three methods: (1) an inviscid free streamline approach, (2) a turbulent boundary layer approach using an eddy viscosity turbulence model and a horizontal pressure gradient determined by the inviscid model, and (3) an approach using the full Navier-Stokes equations with three turbulence models; i.e., an eddy viscosity model, a turbulence kinetic-energy model and a two-equation model with an additional transport equation for the turbulence length scale. A comparison of the performance of the different turbulence models is given, indicating that only the two-equation model adequately accounts for the convective character of turbulence. Turbulence flow property predictions obtained from the turbulence kinetic-energy model with prescribed length scale are only insignificantly better than those obtained from the eddy viscosity model. A parametric study includes the effects of the variation of the characteristics parameters of the assumed logarithmic approach velocity profile. For the case of the forward-facing step, it is shown that in the downstream flow region an increase of the surface roughness gives rise to higher turbulence levels in the shear layer originating from the step corner.

Bitte, J.↗

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

Effect of Mach number and Reynolds number on a normal shock-wave/turbulent boundary-layer interaction

Experimental values of wall pressure and skin friction have been obtained for Mach numbers from 1.32 to 1.48 and for Reynolds numbers from 8.5 x 10 to the 6th to 225 x 10 to the 6th. Increasing the Mach number reduces the wall shear and promotes incipient separation. Reynolds number variations have little effect except at the lowest value where the flow abruptly approaches separation. Comparisons are made with solutions to the time-dependent, Reynolds-averaged, Navier-Stokes equations incorporating a two-equation, Wilcox-Rubesin turbulence model. The computations are in agreement with the experimental results. Additional numerical results indicated that the wind-tunnel walls constrained the flow and suppressed the formation of a separation bubble at the shock wave.

Mateer, G. G.↗

An evaluation of several compressible turbulent boundary-layer models Effects of pressure gradient and Reynolds number

Computations, employing several turbulence models, are compared with a series of attached supersonic turbulent boundary-layer experiments over an extensive range of Reynolds numbers (11.7 x 10 to the 6th to 314 x 10 to the 6th). These experiments included measurements of surface pressure and skin friction for adverse pressure gradients ranging in magnitude from those of previous investigations to an order of magnitude greater. The turbulence models evaluated include algebraic and two-equation eddy-viscosity models and two full Reynolds stress models. In general all the models tested performed well independently of the magnitude of the pressure gradient or Reynolds number and could predict the measured skin friction for most cases with sufficient accuracy for engineering purposes.

Horstman, C. C.↗

Prediction of swirling reacting flow in ramjet combustors

Numerical computations have been undertaken for a basic two-dimensional axisymmetric flowfield which is similar to that found in conventional gas turbine and ramjet combustors. A swirling flow enters a larger chamber via a sudden or gradual expansion. The calculation method involves a staggered grid system for axial and radial velocities, a line relaxation procedure for efficient solution of the equations, a two-equation turbulence energy-turbulence dissipation rate turbulence model, a stairstep boundary representation of the expansion flow, and realistic accommodation of swirl effects. The results include recirculation zone characterization and predicted mean streamline patterns. Predictions with and without chemical reaction are obtained. An associated isothermal experimental flow study is providing a useful data base. Successful outcomes of the work can be incorporated into the more combustion- and hardware-oriented activities of industrial concerns.

Lilley, D. G.↗

Computation of turbulent flows

The essential ingredients and the effectiveness of several levels of turbulent-flow partial differential equation models are considered. Zero-equation models use only the partial differential equation (pde) for the mean velocity field and do not employ turbulence pde's. One-equation models make use of an additional pde relating to the turbulence velocity scale. Attention is also given to two-equation models, stress-equation models, and large eddy simulations. Large-eddy simulations are concerned with a three-dimensional time-dependent numerical computation of the large-scale turbulence.

Reynolds, W. C.↗

A critical evaluation of invariant second-order closure models for subsonic boundary layers

Computations based on several second-order turbulence models, including full Reynolds stress and two-equation models, are compared with a number of boundary-layer experiments. In general, the models represent the data reasonably well, with skin friction tending to be somewhat overpredicted in the far downstream region of the adverse pressure gradient experiments. A discussion of the behavior of the ARAP full Reynolds stress model in predicting the components of the Reynolds stress tensor is given. It is concluded that compatibility at the wall may necessitate the use of more than one length scale.

Rubesin, M. W.↗

Progress in turbulence modeling for complex flow fields including effects of compressibility

Two second-order-closure turbulence models were devised that are suitable for predicting properties of complex turbulent flow fields in both incompressible and compressible fluids. One model is of the "two-equation" variety in which closure is accomplished by introducing an eddy viscosity which depends on both a turbulent mixing energy and a dissipation rate per unit energy, that is, a specific dissipation rate. The other model is a "Reynolds stress equation" (RSE) formulation in which all components of the Reynolds stress tensor and turbulent heat-flux vector are computed directly and are scaled by the specific dissipation rate. Computations based on these models are compared with measurements for the following flow fields: (a) low speed, high Reynolds number channel flows with plane strain or uniform shear; (b) equilibrium turbulent boundary layers with and without pressure gradients or effects of compressibility; and (c) flow over a convex surface with and without a pressure gradient.

Wilcox, D. C.↗

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

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

Effects of turbulence model selection on the prediction of complex aerodynamic flows

Numerical simulations of viscous transonic flow over a circular-arc airfoil and in a diffuser are described. The simulations are made with a new computer program designed to serve as a tool in the development of improved turbulence models for complex flows. The program incorporates zero-, one-, and two-equation eddy viscosity models and includes a variety of subsonic and supersonic boundary conditions. The airfoil flow contains a shock-separated boundary-layer interaction that has resisted previous attempts at simulation. The diffuser flow also contains a shock-boundary-layer interaction, which has not been simulated previously. Calculations using standard turbulence models, developed originally for incompressible unseparated flows, are described. Results indicate that although there are interesting differences in predictions between the various models, none of them predict the flows accurately. Suggestions for improved turbulence models are discussed.

Coakley, T. J.↗

Behavior of a turbulent boundary layer subjected to sudden transverse strain

Data from two experiments on the development of the components of the Reynolds stress tensor after a sudden application of transverse strain are compared. Computations were based on four different turbulence models: a first-order mixing length model, a second-order two-equation eddy viscosity model, and two second-order Reynolds stress models. The second-order models do not produce dramatic improvements over the simple mixing length model. The Reynolds stress models still need development to represent the physics of shear-strained turbulence well. Finally, it is demonstrated that the assumption of a scalar eddy viscosity, often used in engineering calculations of three-dimensional boundary layers, is quite reasonable.

Higuchi, H.↗

Influence of external disturbances and compressibility on free turbulent mixing

It is shown that disturbances in external flow can significantly affect, by as much as an order of magnitude, the turbulent mixing rate in free shear layers and that the length scale of the external flow disturbances is as important as the amplitude. The difference between the effect of wide-band and narrow-band disturbances is stressed. The model for pressure fluctuation term in the kinetic energy equation is included in a two-equation model. The reduced spreading rate in high Mach number, high Reynolds number, adiabatic, free turbulent shear layers is predicted.

Oh, Y. H.↗

Time-History Statistics of Soot Formation in A Model Gas Turbine Combustor

Soot formation is a complex dynamic and intermittent process determined by properties of the fuel, combustor design, and combustor operation. Although the major steps in soot formation (i.e., formation of precursors, inception, growth and evolution) are similar for a variety of carbonaceous fuels, applications, and operating conditions, it remains unclear when the temporal transition between these steps occurs. An engineering prediction tool coupled with computational fluid physics (CFD), therefore needs to accurately model all these complex steps. To develop such a model, we propose the time-history concept for understanding the time dependency of soot formation as a function of local properties (i.e., temperature, velocity, local fuel air ratio, etc.). We continue our previous work with modeling the DLR aero-combustor [1] with our updated in-house CFD code, Open National Combustion Code (OpenNCC), that now includes a Multiple Time-Scale Flamelet Progress Variable approach and a the semi-empirical two-equation soot model. We injected massless tracer particles upstream of the injector region of the combustor to collect time-history statistics of the solution variables. The correlations between the collected statistics with respect to the experimental soot volume fraction data showed that time-history effect of certain flow variables, including turbulent kinetic energy (TKE), and multiple species is indeed important for soot formation. We then conducted a time-history based correlation analysis to determine the key species and the concentration ranges critical for soot formation (C6H5-based nucleation, acetylene-based surface growth, and oxidation with OH and O2). Based on the time-history correlation coefficient (THCC) analysis, we propose possible modifications to improve the current two-equation model.

LES↗

Effect of turbulence models on criticality conditions in swirling flows

The critical state of vortex cores downstream of vortex breakdown has been studied. Base vortical flows were computed using the Reynolds-averaged, axisymmetric Navier-Stokes equations. Standard K - epsilon, RNG and second-order Reynolds stress models were employed. Results indicate that the return to supercriticality is highly dependent on the turbulence model. The K - epsilon model predicted a rapid return of the vortex to supercritical conditions, the location of which showed little sensitivity to changes in the swirl ratio. The Reynolds stress model predicted that the vortex remains subcritical to the end of the domain for each of the swirl ratios employed, and provided results in qualitative agreement with experimental work. The RNG model produced intermediate results, with a downstream movement in the critical location with increasing swirl. Calculations for which area reductions were introduced at the exit in a subcritical flow were also performed using the Reynolds stress model. The structure of the resulting recirculation zone was altered significantly. However, when area reductions were employed within supercritical flows as predicted using the two-equation models, no significant influence on the recirculation zone was noted.

Gatski, Thomas B.↗