Numerical research for the Burgers' model turbulence based on the characteristic functional formalism
Investigation of one dimensional Burgers model of turbulence
Hosokawa, I.↗
SEARCH · Search NASA
Results for “turbulence model”
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.
Computational simulation of corotational vortex merger using 0, 1, and 2 equation turbulence models
A computational study has been conducted to gain a better understanding of the role of turbulence in the merger of corotating vortices and in the decay of a single vortex. The flow field is simulated by using an inviscid flow model and zero-, one-, and two-equation turbulent-flow models. An implicit finite-difference procedure is employed to integrate the unsteady, two-dimensional, incompressible flow equations in a cross-plane. In the zero- and one-equation models, a spatially varying mixing-length model that simulates the streamline curvature effect is used to define an eddy viscosity. The computational merging distance criterion is in good agreement with the experimental one. In the case of the two-equation model, the turbulent kinetic energy and a modified rate of dissipation equation incorporating a curvature correction term are solved. A comparative study of the computational results corresponding to the various flow models and the available experimental data is presented.
Raj, P.↗
Three velocity component, nonhomogeneous atmospheric boundary layer turbulence modeling
The vertical nonhomogeneous character of turbulence in the atmospheric boundary layer results in a non-stationary turbulence process relative to an aircraft during takeoff and landing despite the fact that the turbulence statistics can be horizontally homogeneous. The simulation of the three components of the turbulent winds which include the nonstationary aspect of atmospheric turbulence is the subject of this paper. A procedure is developed and demonstrated to generate the three components of a turbulence ramdom process field, u sub i(x,z) where x and z denote horizontal and vertical coordinates and u sub i, i = 1,2,3 are the three orthogonal components of the turbulent random field. This field satisfies any desired one point auto spectra as well as two point statistics (interlevel correlations). By use of Taylors frozen eddy hypothesis we can transform the turbulent random field into the time domain and obtain the random turbulence along an aircraft trajectory.
Perlmutter, M.↗
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.↗
Kolmogorov Behavior of Near-Wall Turbulence and Its Application in Turbulence Modeling
The near-wall behavior of turbulence is re-examined in a way different from that proposed by Hanjalic and Launder and followers. It is shown that at a certain distance from the wall, all energetic large eddies will reduce to Kolmogorov eddies (the smallest eddies in turbulence). All the important wall parameters, such as friction velocity, viscous length scale, and mean strain rate at the wall, are characterized by Kolmogorov microscales. According to this Kolmogorov behavior of near-wall turbulence, the turbulence quantities, such as turbulent kinetic energy, dissipation rate, etc. at the location where the large eddies become Kolmogorov eddies, can be estimated by using both direct numerical simulation (DNS) data and asymptotic analysis of near-wall turbulence. This information will provide useful boundary conditions for the turbulent transport equations. As an example, the concept is incorporated in the standard k-epsilon model which is then applied to channel and boundary flows. Using appropriate boundary conditions (based on Kolmogorov behavior of near-wall turbulence), there is no need for any wall-modification to the k-epsilon equations (including model constants). Results compare very well with the DNS and experimental data.
Shih, Tsan-Hsing↗
Separated transonic airfoil flow calculations with a nonequilibrium turbulence model
Navier-Stokes transonic airfoil calculations based on a recently developed nonequilibrium, turbulence closure model are presented for a supercritical airfoil section at transonic cruise conditions and for a conventional airfoil section at shock-induced stall conditions. Comparisons with experimental data are presented which show that this nonequilibrium closure model performs significantly better than the popular Baldwin-Lomax and Cebeci-Smith equilibrium algebraic models when there is boundary-layer separation that results from the inviscid-viscous interactions.
King, L. S.↗
Modeling turbulent boundary layers in adverse pressure gradients
Many of the turbulent layers encountered in practical flows develop in adverse pressure gradients; hence, the dynamics of the thickening and possible separation of the boundary layer has important implications for design practices. What are the key physical processes that govern how a turbulent boundary layer responds to an adverse pressure gradient, and how should these processes be modeled? Despite the ubiquity of such flows in engineering and nature, these equations remain largely unanswered. The turbulence closure models presently used to describe these flows commonly use 'wall functions' that have ad hoc corrections for the effects of pressure gradients. There is, therefore, a practical and theoretical need to examine the effects of adverse pressure gradients on wall bounded turbulent flows in order to develop models based on sound physical principle. The evolution of a turbulent boundary layer on a flat wall with an externally imposed pressure gradient is studied.
Belcher, Stephen E.↗
Measurement of terms and parameters in turbulent models
Experimental measurements of the mean and turbulent velocity field in a water flow, turn-around-duct is documented. The small radius of curvature duct experiments were made over a range of Reynolds numbers (based on a duct height of 10 cm) from 70,000 to 500,000. For this particular channel, the flow is dominated by the inertia forces. Use of the local bulk velocity to non-dimensionalize the local velocity was found to limit Reynolds number effects to the regions very close to the wall. Only secondary effects on the flow field were observed when the inlet or exit boundary conditions were altered. The flow over the central two-thirds of the channel was two-dimensional. Mean tangetial and radial velocities, streamlines, pressure distributions, surface shear stress; tangential, radial and lateral turbulent velocities and the Reynolds turbulent shear values are tabulated in other reports. It is evident from the experimental study that a complex numerical modeling technique must be developed to predict the flow in the turn-around-duct. The model must be able to predict relaminarization along the inner-convex-wall. It must also allow for the major increase in turbulence produced by the outer-concave-wall.
Sandborn, Virgil A.↗
A critical comparison of several low Reynolds number k-epsilon turbulence models for flow over a backward facing step
Turbulent backward-facing step flow was examined using four low turbulent Reynolds number k-epsilon models and one standard high Reynolds number technique. A tunnel configuration of 1:9 (step height: exit tunnel height) was used. The models tested include: the original Jones and Launder; Chien; Launder and Sharma; and the recent Shih and Lumley formulation. The experimental reference of Driver and Seegmiller was used to make detailed comparisons between reattachment length, velocity, pressure, turbulent kinetic energy, Reynolds shear stress, and skin friction predictions. The results indicated that the use of a wall function for the standard k-epsilon technique did not reduce the calculation accuracy for this separated flow when compared to the low turbulent Reynolds number techniques.
Steffen, C. J., Jr.↗
Computation of unsteady turbulent boundary layers with flow reversal and evaluation of two separate turbulence models
A procedure which solves the governing boundary layer equations within Keller's box method was developed for calculating unsteady laminar flows with flow reversal. This method is extended to turbulent boundary layers with flow reversal. Test cases are used to investigate the proposition that unsteady turbulent boundary layers also remain free of singularities. Turbulent flow calculations are performed. The governing equations for both models are solved. As in laminar flows, the unsteady turbulent boundary layers are free from singularities, but there is a clear indication of rapid thickening of the boundary layer with increasing flow reversal. Predictions of both turbulence models are the same for all practical purposes.
Cebeci, T.↗
A critical comparison of several low Reynolds number k-epsilon turbulence models for flow over a backward-facing step
Turbulent backward-facing step flow was examined using four low turbulent Reynolds number k-epsilon models and one standard high Reynolds number technique. A tunnel configuration of 1:9 (step height: exit tunnel height) was used. The models tested include: the original Jones and Launder; Chien; Launder and Sharma; and the recent Shih and Lumley formulation. The experimental reference of Driver and Seegmiller was used to make detailed comparisons between reattachment length, velocity, pressure, turbulent kinetic energy, Reynolds shear stress, and skin friction predictions. The results indicated that the use of a wall function for the standard k-epsilon technique did not reduce the calculation accuracy for this separated flow when compared to the low turbulent Reynolds number techniques.
Steffen, Christopher J., Jr.↗
On local approximations of the pressure-strain term in turbulence models
The results of numerical simulations of turbulent channel flows were used to examine the validity of the local approximation of the pressure-strain term in the Reynolds stress transport equation. Outside of the viscous sublayer the local approximation compares very well with the exact pressure strain. This agreement is due, at least in part, to the high correlation between the rapid pressure and its Laplacian, which suggests that only the near parts of the flow contribute to the rapid pressure at a point. In the viscous sublayer the distance over which the mean shear can be considered constant is comparable to the length scale in the normal direction of the correlations of velocity gradients, leading to failure of the local approximation.
Bradshaw, P.↗
Advancements in engineering turbulence modeling
Some new developments in two-equation models and second order closure models are presented. Two-equation models (k-epsilon models) have been widely used in computational fluid dynamics (CFD) for engineering problems. Most of low-Reynolds number two-equation models contain some wall-distance damping functions to account for the effect of wall on turbulence. However, this often causes the confusion and difficulties in computing flows with complex geometry and also needs an ad hoc treatment near the separation and reattachment points. A set of modified two-equation models is proposed to remove the aforementioned shortcomings. The calculations using various two-equation models are compared with direct numerical simulations of channel flow and flat boundary layers. Development of a second order closure model is also discussed with emphasis on the modeling of pressure related correlation terms and dissipation rates in the second moment equations. All the existing models poorly predict the normal stresses near the wall and fail to predict the 3-D effect of mean flow on the turbulence (e.g. decrease in the shear stress caused by the cross flow in the boundary layer). The newly developed second order near-wall turbulence model is described and is capable of capturing the near-wall behavior of turbulence as well as the effect of 3-D mean flow on the turbulence.
Shih, T.-H.↗
Modeling turbulent compressible flows - The mass fluctuating velocity and squared density
This paper deals with single-point closure theory for compressible turbulent flow, including the effects of compressibility on the turbulence. In particular, the combination of the pressure dilatation and the dilatation dissipation, terms which appear on the turbulent kinetic energy equation, are modeled. Model parameters in these transport equations are determined by comparing predictions with boundary layer measurements. Finally, predictions with a k-epsilon model, including the new formulations, are presented for the compressible shear layer.
Taulbee, D.↗
Influence of freestream values on k-omega turbulence model predictions
The Wilcox (1988, 1991) k-omega model for eddy-viscosity turbulence does not require damping functions in the viscous sublayer, and its equations are less stiff near the wall. It has been designed to predict the requisite wake length in equilibrium, adverse pressure-gradient boundary-layer flows. When applied to free shear layers, however, a strong dependency of its results on the freestream value of omega has been noted. This feature is presently investigated via the self-similar equations for incompressible equilibrium boundary layers.
Menter, F. R.↗
Turbulence modeling for non-equilibrium flow
The work performed during this year has involved further assessment and extension of the k-epsilon-v(exp 2) model, and initiation of work on scalar transport. The latter is introduced by the contribution of Y. Shabany to this volume. Flexible, computationally tractable models are needed for engineering CFD. As computational technology has progressed, the ability and need to use elaborate turbulence closure models has increased. The objective of our work is to explore and develop new analytical frameworks that might extend the applicability of the modeling techniques. In past years the development of a method for near-wall modeling was described. The method has been implemented into a CFD code and its viability has been demonstrated by various test cases. Further tests are reported herein. Non-equilibrium near-wall models are needed for some heat transfer applications. Scalar transport seems generally to be more sensitive to non-equilibrium effects than is momentum transport. For some applications turbulence anisotropy plays a role and an estimate of the full Reynolds stress tensor is needed. We have begun work on scalar transport per se, but in this brief I will only report on an extension of the k-epsilon-v(exp 2) model to predict the Reynolds stress tensor.
Durbin, P. A.↗
Modelling turbulent flame ignition and blowout
A statistical mixing model incorporating an overall rate equation to describe the fuel oxidation process was developed for studies of ignition and blowout in a combustor primary zone. This zone is treated as a partially stirred reactor whose composition is described by a statistical ensemble of equal mass fluid elements. This ensemble experiences mixing interactions, which represent the turbulent mixing process, at time intervals governed by an empirically determined mixing frequency. Each mixing interaction is computed by ramdomly selecting two different elements which are then allowed to mix completely so that they reach a mean composition depending on their thermodynamic states prior to mixing. The two elements then separate, and the chemical kinetics proceed depending on their new composition and temperature.
Radhakrishnan, K.↗
A Lagrangian dynamic subgrid-scale model turbulence
A new formulation of the dynamic subgrid-scale model is tested in which the error associated with the Germano identity is minimized over flow pathlines rather than over directions of statistical homogeneity. This procedure allows the application of the dynamic model with averaging to flows in complex geometries that do not possess homogeneous directions. The characteristic Lagrangian time scale over which the averaging is performed is chosen such that the model is purely dissipative, guaranteeing numerical stability when coupled with the Smagorinsky model. The formulation is tested successfully in forced and decaying isotropic turbulence and in fully developed and transitional channel flow. In homogeneous flows, the results are similar to those of the volume-averaged dynamic model, while in channel flow, the predictions are superior to those of the plane-averaged dynamic model. The relationship between the averaged terms in the model and vortical structures (worms) that appear in the LES is investigated. Computational overhead is kept small (about 10 percent above the CPU requirements of the volume or plane-averaged dynamic model) by using an approximate scheme to advance the Lagrangian tracking through first-order Euler time integration and linear interpolation in space.
Meneveau, C.↗