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 19 records

A low-Reynolds-number two-equation turbulence model for predicting heat transfer on turbine blades

A modified form of the Lam-Bremhorst low-Reynolds number kappa-epsilon turbulence model was developed for predicting transitional boundary layer flows under conditions characteristic of gas turbine blades. The application of the model to flows with pressure gradients is described. Tests against a number of turbine blade cascade data sets are included. Some additional refinements of the model that were made in recent months are explained.

Patankar, Suhas V.↗

Critical examination of two-equation turbulence closure models

Comparison of the Jones-Launder, Ng-Spalding, Saffman-Wilcox, and Wilcox-Traci two-equation turbulence models has been conducted. It was shown that the Saffman-Wilcox and Wilcox-Traci dissipation-rate formulations admit straightforward integration through the viscous sublayer, whereas integration through the viscous sublayer is a more difficult issue with the Jones-Launder dissipation-function and the Ng-Spalding length-scale formulations. Numerical computations were conducted in which the models were applied to four equilibrium boundary layer flows including adverse, zero, and favorable pressure gradients. Computations of zero pressure gradient flow over a convex wall composed the final part of the comparison.

Chambers, T. L.↗

Modeling of three-dimensional mixing and reacting ducted flows

A computer code, based upon a finite element solution algorithm, was developed to solve the governing equations for three-dimensional, reacting boundary region, and constant area ducted flow fields. Effective diffusion coefficients are employed to allow analyses of turbulent, transitional or laminar flows. The code was used to investigate mixing and reacting hydrogen jets injected from multiple orifices, transverse and parallel to a supersonic air stream. Computational results provide a three-dimensional description of velocity, temperature, and species-concentration fields downstream of injection. Experimental data for eight cases covering different injection conditions and geometries were modeled using mixing length theory (MLT). These results were used as a baseline for examining the relative merits of other mixing models. Calculations were made using a two-equation turbulence model (k+d) and comparisons were made between experiment and mixing length theory predictions. The k+d model shows only a slight improvement in predictive capability over MLT. Results of an examination of the effect of tensorial transport coefficients on mass and momentum field distribution are also presented. Solutions demonstrating the ability of the code to model ducted flows and parallel strut injection are presented and discussed.

Zelazny, S. W.↗

Analysis of neutrally stable atmospheric flow over a two-dimensional forward facing step

A mathematical model describing ground-wind induced flow fields around surface obstructions such as buildings, bridges or other man-made structures is presented. The Navier-Stokes equations with a two-equation turbulence model are used to analyze the flow over a two-dimensional forward facing step. A study of the assumed logarithmic approach velocity profile shows that an increase of surface roughness produces a rise to higher turbulence levels in the shear layer originating from the step corner. This rise leads to higher Reynolds stress and to faster reattachment of the separated flow.

Bitte, J.↗

Computation of turbulent boundary layers on curved surfaces, 1 June 1975 - 31 January 1976

An accurate method was developed for predicting effects of streamline curvature and coordinate system rotation on turbulent boundary layers. A new two-equation model of turbulence was developed which serves as the basis of the study. In developing the new model, physical reasoning is combined with singular perturbation methods to develop a rational, physically-based set of equations which are, on the one hand, as accurate as mixing-length theory for equilibrium boundary layers and, on the other hand, suitable for computing effects of curvature and rotation. The equations are solved numerically for several boundary layer flows over plane and curved surfaces. For incompressible boundary layers, results of the computations are generally within 10% of corresponding experimental data. Somewhat larger discrepancies are noted for compressible applications.

Wilcox, D. C.↗

Application of a general boundary layer analysis to turbulent boundary layers subjected to strong favorable pressure gradients

Theoretical predictions of turbulent boundary layer development under the influence of strong favorable pressure gradients made using a finite-difference calculation procedure are compared to experimental data. Comparisons are presented for low speed flows with and without wall heat transfer as well as for supersonic flows with adiabatic walls. The turbulence model used is governed by an integral form of the turbulence kinetic energy equation and the results are compared with predictions made using a conventional equilibrium turbulence model based upon Prandtl's mixing length, a Clauser-type eddy viscosity model used by Cebecci and Mosinskis, and a two-equation turbulence energy model of Launder and Jones.

Kreskovsky, J. P.↗

Theoretical study of reactive and nonreactive turbulent coaxial jets

The hydrodynamic properties and the reaction kinetics of axisymmetric coaxial turbulent jets having steady mean quantities are investigated. From the analysis, limited to free turbulent boundary layer mixing of such jets, it is found that the two-equation model of turbulence is adequate for most nonreactive flows. For the reactive flows, where an allowance must be made for second order correlations of concentration fluctuations in the finite rate chemistry for initially inhomogeneous mixture, an equation similar to the concentration fluctuation equation of a related model is suggested. For diffusion limited reactions, the eddy breakup model based on concentration fluctuations is found satisfactory and simple to use. The theoretical results obtained from these various models are compared with some of the available experimental data.

Gupta, R. N.↗

Turbulent atmospheric flow over a backward-facing step

The phenomenon of atmospheric shear layer separation over a man-made structure such as a building (modeled as a backward-facing step) has been analyzed theoretically by (1) solving the two-dimensional equations of motion in the two variables, stream function and vorticity, and by (2) employing an approximate integral technique. Boundary conditions for the undisturbed flow are that of the turbulent atmospheric shear flow over a rough terrain. In the first approach a two-equation model of turbulence was used. In the second approach an approximate technique was utilized in an attempt to describe the details of the flow in the recirculation zone behind the step. The results predict velocity profiles in sufficient detail that the presence of the corner eddy in the region of negative surface pressure gradient is evident. The magnitude of the reversed flow velocity in the recirculation eddy has been found to agree with that found from experiments. Also, a surface eddy viscosity distribution has been an outgrowth of the method which realistically follows the magnitude of the surface pressure gradient distribution as found experimentally.

Kaul, U. K.↗

A complete model of turbulence

A set of constitutive equations suitable for a priori computation of turbulent shear flows has been developed. Since no properties of a given turbulent flow need be known in advance in order to obtain a solution, the equations comprise a complete model of turbulence. Perturbation analysis shows that the model predicts a composite five-layer structure for an incompressible turbulent boundary layer, viz, a defect layer, a law-of-the-wall layer, a viscous sublayer, a near-surface roughness layer, and a viscous superlayer at the boundary-layer edge. Analysis of the defect layer demonstrates the key improvement of the model over its predecessor, the Saffman-Wilcox two-equation model of turbulence. Examination of model-predicted sublayer structure yields model-parameter boundary conditions appropriate for surfaces with roughness and mass injection. Results of numerical computations of compressible and incompressible equilibrium boundary layers show that, for such flows, the model is as accurate as mixing-length theory. Applications to transitional boundary layers and to nonequilibrium relaxation of a boundary layer passing from a rough to a smooth surface indicate that the model's applicability extends far beyond that of mixing-length theory's.

Wilcox, D. C.↗

Recent Improvements in Pronghorn for Advanced Reactor Modeling

Pronghorn is a thermal-hydraulics computational tool developed using the Idaho National Laboratory's Multiphysics Object-Oriented Simulation Environment (MOOSE). It is designed to support Computational Fluid Dynamics (CFD) modeling, ranging from subchannel and porous media analysis to Reynolds Averaged Navier-Stokes (RANS) turbulence modeling. As an integral part of the MOOSE-based suite of tools, Pronghorn seamlessly couples with other MOOSE-based applications to simulate a variety of physical phenomena. This article highlights recent significant enhancements to Pronghorn's CFD modeling capabilities and demonstrates their application to advanced nuclear reactor designs. The recent improvements in Pronghorn primarily focus on modifications to its turbulence modeling capabilities, near-wall corrections and numerical schemes. In terms of turbulence modeling, the two-equation $k-\epsilon$ and $k-\omega$ SST models have been implemented and validated with both equilibrium and non-equilibrium wall treatments. Additionally, corrections for wall roughness, and curvature, and wall-channeling in pebble beds have been introduced in the near-wall modeling. These developments enable more accurate simulations of advanced nuclear reactors. Two case studies are presented in this work: a pool-type Molten Chloride Reactor and a salt-cooled Pebble-Bed High Temperature Reactor. In both cases, the previous models in Pronghorn are compared with the new implementations, demonstrating the improved accuracy achieved with the updated models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

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

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

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

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

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