Engineering PapersSearch

Engineering topics

Reynolds, W. C.

Publications and source records attributed to Reynolds, W. C..

At least 19 records

Nonlinear Stability and Structure of Compressible Reacting Mixing Layers

The parabolized stability equations (PSE) are used to investigate issues of nonlinear flow development and mixing in compressible reacting shear layers. Particular interest is placed on investigating the change in flow structure that occurs when compressibility and heat release are added to the flow. These conditions allow the 'outer' instability modes- one associated with each of the fast and slow streams-to dominate over the 'central', Kelvin-Helmholtz mode that unaccompanied in incompressible nonreacting mixing layers. Analysis of scalar probability density functions in flows with dominant outer modes demonstrates the ineffective, one-sided nature of mixing that accompany these flow structures. Colayer conditions, where two modes have equal growth rate and the mixing layer is formed by two sets of vortices, offer some opportunity for mixing enhancement. Their extent, however, is found to be limited in the mixing layer's parameter space. Extensive validation of the PSE technique also provides a unique perspective on central- mode vortex pairing, further supporting the view that pairing is primarily governed perspective sheds insight on how linear stability theory is able to provide such an accurate prediction of experimentally-observed, fully nonlinear flow phenomenon.

Day, M. J.

Structure of the Compressible Reacting Mixing Layer: A Linear Stability Analysis

Previous investigations have demonstrated that a mixing layer in compressible reacting conditions can develop two instability modes in addition to the more common central mode, which exists unaccompanied in incompressible non-reacting flows. These 'outer' modes were shown in the investigations of Planche and Reynolds to have a detrimental effect on mixing and combustion efficiency and therefore their presence has important consequences for applications of supersonic combustion. Results from a parametric study of the compressible reacting mixing layer's regime space using a linear stability technique will be discussed. The purpose is to develop a more complete understanding behind the effects of compressibility, heat release and the density, equivalence and velocity ratios on the amplification rate and structure of each instability mode.

Day, M. J.

A Particle Representation Model for the Deformation of Homogeneous Turbulence

In simple flows, where the mean deformation rates are mild and the turbulence has time to come to equilibrium with the mean flow, the Reynolds stresses are determined by the applied strain rate. Hence in these flows, it is often adequate to use an eddy-viscosity representation. The modern family of kappa-epsilon models has been very useful in predicting near equilibrium turbulent flows, where the rms deformation rate S is small compared to the reciprocal time scale of the turbulence (epsilon/kappa). In modern engineering applications, turbulence models are quite often required to predict flows with very rapid deformations (large S kappa/epsilon). In these flows, the structure takes some time to respond and eddy viscosity models are inadequate. The response of turbulence to rapid deformations is given by rapid distortion theory (RDT). Under RDT the nonlinear effects due to turbulence-turbulence interactions are neglected in the governing equations, but even when linearized in this fashion, the governing equations are unclosed at the one-point level due to the non-locality of the pressure fluctuations.

Kassinos, S. C.

An extended structure-based model based on a stochastic eddy-axis evolution equation

We have proposed and implemented an extension of the structure-based model for weak deformations. It was shown that the extended model will correctly reduce to the form of standard k-e models for the case of equilibrium under weak mean strain. The realizability of the extended model is guaranteed by the method of its construction. The predictions of the proposed model were very good for rotating homogeneous shear flows and for irrotational axisymmetric contraction, but were seriously deficient in the case of plane strain and axisymmetric expansion. We have concluded that the problem behind these difficulties lies in the algebraic constitutive equation relating the Reynolds stresses to the structure parameters rather than in the slow model developed here. In its present form, this equation assumes that under irrotational strain the principal axes of the Reynolds stresses remain locked onto those of the eddy-axis tensor. This is correct in the RDT limit, but inappropriate under weaker mean strains, when the non-linear eddy-eddy interactions tend to misalign the two sets of principal axes and create some non-zero theta and gamma.

Kassinos, S. C.

Modeling the two-point correlation of the vector stream function

A new model for the two-point vector stream function correlation has been developed using tensor invariant arguments and evaluated by the comparison of model predictions with DNS data for incompressible homogeneous turbulent shear flow. This two-point vector stream function model correlation can then be used to calculate the two-point velocity correlation function and other quantities useful in turbulence modeling. The model assumes that the two-point vector stream function correlation can be written in terms of the separation vector and a new tensor function that depends only on the magnitude of the separation vector. The model has a single free model coefficient, which has been chosen by comparison with the DNS data. The relative error of the model predictions of the two-point vector stream function correlation is only a few percent for a broad range of the model coefficient. Predictions of the derivatives of this correlation, which are of interest in turbulence modeling, may not be this accurate.

Oberlack, M.

Compressibility Effects on the Passive Scalar Flux Within Homogeneous Turbulence

Compressibility effects on turbulent transport of a passive scalar are studied within homogeneous turbulence using a kinematic decomposition of the velocity field into solenoidal and dilatational parts. It is found that the dilatational velocity does not produce a passive scalar flux, and that all of the passive scalar flux is due to the solenoidal velocity.

Blaisdell, G. A.

Compressibility Effects on the Growth and Structure of Homogeneous Turbulent Shear Flow

Compressibility effects within decaying isotropic turbulence and homogeneous turbulent shear flow have been studied using direct numerical simulation. The objective of this work is to increase our understanding of compressible turbulence and to aid the development of turbulence models for compressible flows. The numerical simulations of compressible isotropic turbulence show that compressibility effects are highly dependent on the initial conditions. The shear flow simulations, on the other hand, show that measures of compressibility evolve to become independent of their initial values and are parameterized by the root mean square Mach number. The growth rate of the turbulence in compressible homogeneous shear flow is reduced compared to that in the incompressible case. The reduced growth rate is the result of an increase in the dissipation rate and energy transfer to internal energy by the pressure-dilatation correlation. Examination of the structure of compressible homogeneous shear flow reveals the presence of eddy shocklets, which are important for the increased dissipation rate of compressible turbulence.

Blaisdell, G. A.

Experiments on an unsteady, three-dimensional separation

Unsteady, three-dimensional flow separation occurs in a variety of technical situations including turbomachinery and low-speed aircraft. An experimental program at Stanford in unsteady, three-dimensional, pressure-driven laminar separation has investigated the structure and time-scaling of these flows; of particular interest is the development, washout, and control of flow separation. Results reveal that a two-dimensional, laminar boundary layer passes through several stages on its way to a quasi-steady three-dimensional separation. The quasi-steady state of the separation embodies a complex, unsteady, vortical structure.

Henk, R. W.

On the Yakhot-Orszag renormalization group method for deriving turbulence statistics and models

An independent, comprehensive, critical review of the 'renormalization group' (RNG) theory of turbulence developed by Yakhot and Orszag (1986) is provided. Their basic theory for the Navier-Stokes equations is confirmed, and approximations in the scale removal procedure are discussed. The YO derivations of the velocity-derivative skewness and the transport equation for the energy dissipation rate are examined. An algebraic error in the derivation of the skewness is corrected. The corrected RNG skewness value of -0.59 is in agreement with experiments at moderate Reynolds numbers. Several problems are identified in the derivation of the energy dissipation rate equations which suggest that the derivation should be reformulated.

Smith, L. M.

Towards a structure-based turbulence model

The author is developing a new type of turbulence model in which a new one-point quantity, the eddy structure tensor, carries information about the two-point structure of the turbulence. The model was motivated by the observation that conventional one-point turbulence models based only on the turbulent stresses do not predict the rapid changes in state that are found when anisotropic homogeneous turbulence is subjected to mean rotation, and hence are fundamentally incorrect for rotation. The model appears to give topologically correct predictions for the changes in stress state and structure state under all types of rapid distortions of homogeneous turbulence.

Reynolds, W. C.

The dissipation-range spectrum and the velocity-derivative skewness in turbulent flows

A dynamical equation for the energy dissipation rate is used together with a model spectrum to predict the velocity derivative skewness at a high Reynolds number. The objective is to determine the best choice of the exponent, m, so that the resulting model spectrum can be used to estimate statistical properties of the fine scales of turbulence. Skewness data suggest that m = 2 is the best model for laboratory-scale flows and m = 1 is a poor model for a broad range of Reynolds numbers.

Smith, L. M.

Three-dimensional simulations of large eddies in the compressible mixing layer

Consideration is given to the effect of compressibility on a plane mixing layer that is a prototype free shear layer, amenable to study by numerical simulation and experiment. The full time-dependent compressible Navier-Stokes equations are solved numerically for a temporally evolving mixing layer employing a mixed spectral and high-order finite difference method. Simulations with random initial conditions confirm the prediction of linear stability theory that at high Mach numbers oblique waves grow faster than two-dimensional waves. Simulations are presented of the nonlinear temporal evolution of the most rapidly amplified linear instability waves.

Sandham, N. D.

Progress in understanding the renormalization group skewness and kappa-epsilon models

The immediate goal is to understand and validate the Yakhot-Orszag model of the velocity-derivative skewness and model equation for the rate of energy dissipation epsilon. A summary of a more detailed manuscript in preparation is presented. The purpose is to clarify some limitations of the theory by careful examination of key assumptions and approximations, and thereby to encourage its improvement.

Smith, Leslie M.

Compressibility effects on the growth and structure of homogeneous turbulent shear flow

Direct numerical simulations of compressible homogeneous turbulent shear flow are used to provide insight into compressibility effects on turbulence. The simulations show a reduction in the growth rate of the turbulence compared to the incompressible case. Examination of the turbulent kinetic energy budget shows that the reduced growth rate is due to an increase in the dissipation rate due to the divergence of the velocity and to the pressure-dilatation correlation which acts to transfer energy between internal energy and kinetic energy. The structure of the turbulence is also examined. Visualizations of the flowfields reveal the presence of eddy shocklets. These shock structures are important contributors to the increased dissipation of compressible turbulence. A mechanism for the generation of the shocks is suggested.

Blaisdell, G. A.

Effects of rotation on homogeneous turbulence

Turbulence models are known to have difficulty in flows with strong rotation. This paper examines the effect of rotation, using rapid distortion theory for homogeneous turbulence as a guide. It is shown that rotation significantly modifies the turbulent stress anisotropy in a way not predicted by current turbulence models. The reasons for this failure of models is argued to be a lack of information about the turbulence structure. A new tensor quantity, the structure tensor, is defined and shown to be of critical importance in the rotation problem. The physical reasons for the effects of rotation are explained, and a simple model that does display the effects is proposed.

Reynolds, W. C.

Scalar entrainment in the mixing layer

New definitions of entrainment and mixing based on the passive scalar field in the plane mixing layer are proposed. The definitions distinguish clearly between three fluid states: (1) unmixed fluid, (2) fluid engulfed in the mixing layer, trapped between two scalar contours, and (3) mixed fluid. The difference betwen (2) and (3) is the amount of fluid which has been engulfed during the pairing process, but has not yet mixed. Trends are identified from direct numerical simulations and extensions to high Reynolds number mixing layers are made in terms of the Broadwell-Breidenthal mixing model. In the limit of high Peclet number (Pe = ReSc) it is speculated that engulfed fluid rises in steps associated with pairings, introducing unmixed fluid into the large scale structures, where it is eventually mixed at the Kolmogorov scale. From this viewpoint, pairing is a prerequisite for mixing in the turbulent plane mixing layer.

Sandham, N. D.

Use of passive scalar tagging for the study of coherent structures in the plane mixing layer

Data obtained from the numerical simulation of a 2-D mixing layer were used to study the feasibility of using the instantaneous concentration of a passive scalar for detecting the typical coherent structures in the flow. The study showed that this technique works quite satisfactorily and yields results similar to those that can be obtained by using the instantaneous vorticity for structure detection. Using the coherent events educed by the scalar conditioning technique, the contribution of the coherent events to the total turbulent momentum and scalar transport was estimated. It is found that the contribution from the typical coherent events is of the same order as that of the time-mean value. However, the individual contributions become very large during the pairing of these structures. The increase is particularly spectacular in the case of the Reynolds shear stress.

Ramaprian, B. R.

Advances in turbulent transport modeling based on direct simulations of turbulence

Recent work using numerical simulations to study the modeling of scalar transport is outlined. A model developed by reference to simulations for homogeneous shear flow is presented. This model predicts cross-gradient transports that can give rise to the illusion of counter-gradient transport. The model is found to work remarkably well when compared to simulations for channel flow, but fails to predict important features of the scalar fluxes in mixing layers. Reasons for this failure are discussed and suggestions are made for future work.

Reynolds, W. C.