Search

The paper reports on an alternative approach to modelling the turbulent scalar fluxes that arise from time averaging the transport equation for a scalar. In this approach, a functional relationship between these fluxes and various tensor quantities is constructed with guidance from the exact equations governing the transport of fluxes. Results from tensor representation theory are then used to obtain an explicit relationship between the fluxes and the terms in the assumed functional relationship. Where turbulence length– and time–scales are implied, these are determined from two scalar quantities: the turbulence kinetic energy and its rate of dissipation by viscous action. The general representation is then reduced by certain justifiable assumptions to yield a practical model for the turbulent scalar fluxes that is explicit and algebraic in these quantities and one that correctly reflects their dependence on the gradients of mean velocity and on the details of the turbulence. Examination of alternative algebraic models shows most to be subsets of the present proposal. The new model is calibrated using results from large–eddy simulations (LESs) of homogeneous turbulence with passive scalars and then assessed by reference to benchmark data from heated turbulent shear flows. The results obtained show the model to correctly predict the anisotropy of the turbulent diffusivity tensor. The asymmetric nature of this tensor is also recovered, but only qualitatively, there being significant quantitative differences between the model predictions and the LES results. Finally, comparisons with data from benchmark two–dimensional free shear flows show the new model to yield distinct improvements over other algebraic scalar–flux closures.

This paper presents a rational approach to modelling the triple velocity correlations that appear in the transport equations for the Reynolds stresses. All existing models of these correlations have largely been formulated on phenomenological grounds and are defective in one important aspect: They all neglect to allow for the dependence of these correlations on the local gradients of mean velocity. The mathematical necessity for this dependence will be domonstrated in the paper. The present contribution lies in the novel use of Group Representation Theory to determine the most general tensorial form of these correlations in terms of all the second- and third-order tensor quantities that appear in the exact equations that govern their evolution. The requisite representation did not exist in the literature and therefore had to be developed specifically for this purpose by Professor G. F. Smith. The outcome of this work is a mathematical framework for the construction of algebraic, explicit, and rational models for the triple velocity correlations that are theoretically consistent and include all the correct dependencies. Previous models are reviewed, and all are shown to be an incomplete subset of this new representation, even to lowest order.

An important consistency condition for non-linear algebraic Reynolds stress models involving their dependence on rotational strains is discussed from a basic theoretical standpoint. A variety of recently proposed models are demonstrated to violate this condition which leads to physically inconsistent results and unrealizable solutions in rotating flows — an undesirable state of affairs that was avoided in the majority of earlier models. It is shown that when non-linear algebraic models are systematically derived from the Reynolds stress transport equation, such inconsistencies are automatically avoided. The implications of these results for turbulence modeling are demonstrated quantitatively by the illustrative example of isotropic turbulence in a rotating frame.

A critical comparison of explicit vs traditional algebraic stress models of turbulence is made in an effort to clear up the confusion that appears to have been generated by the recently published literature on the subject, in which disparate approaches are adopted. The only way that general second-order closures can formally lead to fully explicit algebraic stress models, in a global sense, is in the limit of equilibrium homogeneous turbulence. When these fully explicit models are then applled to turbulent flows that are far from equilibrium, a singularity can arise, which can be removed by a systematic regularization. When solved explicitly either analytically or numerically, the traditional, implicit algebraic stress models have either multiple solutions or singularities, which tends to explain why they have had problems in applications to complex flows. Thus, it is argued that traditional algebraic stress models are intrinsically ill-behaved and should be abandoned in future applications in favor of regularized, explicit algebraic stress models. It is furthermore argued that these should be based on the homogeneous equilibrium hypothesis, which allows for more general second-order closures to be used to obtain single-valued models.

We present an explicit algebraic stress model (EASM) based on the more current pressure-strain model of Speziale et al. (the SSG model). The SSG model represents a recent upgrade of the Launder et al. model that does not require wall reflection terms. The ability of the proposed model to predict nonequilibrium flows over an airfoil and a wing will be established. The ISAAC Navier-Stokes code is used to simulate the flows considered.

During the past few years, explicit algebraic stress models have been developped that are formally consistent with full second-order closures in the limit of homogeneous turbulence in equilibrium. These models allow for the solution of complex turbulent flows with a substantially reduced level of computation compared with full second-order closures, since they constitute two-equation models.The purpose of the present note is to show how the effects of anisotropic dissipation can be systematically incorporated into these explicit algebraic stress models by a simple readjustmentof the coefficients.

A new approach for the development of second-order closure models suitable for nonequilibrium turbulent shear flows is presented. The central idea is based on the implementation of a relaxation time approximation to the nonequilibrium extension of an explicit algebraic stress model. Here, the algebraic stress model is systematically derived from the full Speziale et al. 1991 (SSG) second-order closure in the equilibrium limit of homogeneous turbulence. It is then extended to nonequilibrium turbulent flows by means of a Pade approximation, whereby approximate consistency with the rapid distortion theory (RDT) solution for homogeneous shear flow is established. The resulting model is tested in homogeneous shear flow turbulence under a wide range of shear rates. Substantially improved results are obtained in comparison to conventional second-order closure models that are based on the direct modeling of terms in the Reynolds stress transport equation.

The consistency of second‐order closure models with results from hydrodynamic stability theory is analyzed for the simplified case of homogeneous turbulence. In a recent study, Speziale, Gatski, and Mac Giolla Mhuiris [Phys. Fluids A 2, 1678 (1990)] showed that second‐order closures are capable of yielding results that are consistent with linear stability theory for the case of homogeneous shear flow in a rotating frame. It is demonstrated in this paper that this success is due to the fact that the stability boundaries for rotating homogeneous shear flow are not dependent on the details of the spatial structure of the disturbances. For those instances where they are—such as in the case of elliptical flows where the instability mechanism is more subtle—the results are not so favorable. The origins and extent of this modeling problem are examined in detail along with a possible resolution based on Rapid Distortion Theory (RDT) and its implications for turbulence modeling.

The explicit algebraic stress model of Gatski and Speziale is considered. This constitutes a two-equation model with an anisotropic eddy viscosity that is systematically derived from the SSG second-order model via the algebraic stress approximation for equilibrium turbulent flows.

The restricted Euler equation captures many important features of the behavior of the velocity gradient tensor observed in direct numerical simulations (DNS) of isotropic turbulence. However, in slightly more complex flows the agreement is not good, especially in regions of low dissipation. In this paper, it is demonstrated that the Reynolds-averaged restricted Euler equation violates the balance of mean momentum for virtually all homogeneous turbulent flows with only two major exceptions: isotropic and homogeneously-sheared turbulence. A new model equation which overcomes this shortcoming and is more widely applicable is suggested. This model is derived from the Navier-Stokes equation with a restricted Euler type approximation made on the fluctuating velocity gradient field. Analytical solutions of the proposed modified restricted Euler equation appear to be difficult to obtain. Hence, a strategy for numerically calculating the velocity gradient tensor is developed. Preliminary calculations tend to indicate that the modified restricted Euler equation captures many important aspects of the behavior of the fluctuating velocity gradients in anisotropic homogeneous turbulence.

Dimensional similarity arguments proposed by Millikan are used with the Morkovin hypothesis to deduce logarithmic laws for compressible turbulent boundary layers as an alternative to the traditional van Driest analysis. It is shown that an overlap exists between the wall layer and the defect layer, and this leads to logarithmic behavior in the overlap region. The von Karman constant is found to depend parametrically on the Mach number based on the friction velocity, the dimensionless total heat flux, and the specific heat ratio. Even though it remains constant at approximately 0.41 for a freestream Mach number range of 0 to 4.544 with adiabatic wall boundary conditions, it rises sharply as the Mach number increases significantly beyond 4.544. The intercept of the logarithmic law of the wall is found to depend on the Mach number based on the friction velocity, the dimensionless total heat flux, the Prandtl number evaluated at the wall, and the specific heat ratio. On the other hand, the intercept of the logarithmic defect law is parametric in the pressure gradient parameter and all of the aforementioned dimensionless variables except the Prandtl number. A skin friction law is also deduced for compressible boundary layers. The skin friction coefficientmore » is shown to depend on the momentum thickness Reynolds number, the wall temperature ratio, and all of the other parameters already mentioned. 26 refs.« less

Explicit algebraic stress models that are valid for three-dimensional turbulent flows in non-inertial frames are systematically derived from a hierarchy of second-order closure models. This represents a generalization of the model derived by Pope (1975) who based his analysis on the Launder, Reece & Rodi model restricted to two-dimensional turbulent flows in an inertial frame. The relationship between the new models and traditional algebraic stress models – as well as anisotropic eddy viscosity models – is theoretically established. A need for regularization is demonstrated in an effort to explain why traditional algebraic stress models have failed in complex flows. It is also shown that these explicit algebraic stress models can shed new light on what second-order closure models predict for the equilibrium states of homogeneous turbulent flows and can serve as a useful alternative in practical computations.

Reynolds stress calculations of homogeneous turbulent shear flow are conducted with a second-order closure model modified to account for non-equilibrium vortex stretching in the dissipation rate transport equation, as recently proposed by Bernard and Speziale. As with the earlier reported k-epsilon model calculations incorporating this vortex stretching effect, a production-equals-dissipation equilibrium is obtained with bounded turbulent kinetic energy and dissipation. However, this equilibrium is not achieved until the dimensionless time greater than 60, an elapsed time that is at least twice as large as any of those considered in previous numerical and physical experiments on homogeneous shear flow. Direct quantitative comparisons between the model predictions and the results of experiments are quite favorable. In particular, it is shown that the inclusion of this non-equilibrium vortex stretching effect has the capability of explaining the significant range of production to dissipation ratios observed in experiments.

Turbulent channel flow and homogeneous shear flow have served as basic building block flows for the testing and calibration of Reynolds stress models. A direct theoretical connection is made between homogeneous shear flow in equilibrium and the log-layer of fully-developed turbulent channel flow. It is shown that if a second-order closure model is calibrated to yield good equilibrium values for homogeneous shear flow it will also yield good results for the log-layer of channel flow provided that the Rotta coefficient is not too far removed from one. Most of the commonly used second-order closure models introduce an ad hoc wall reflection term in order to mask deficient predictions for the log-layer of channel flow that arise either from an inaccurate calibration of homogeneous shear flow or from the use of a Rotta coefficient that is too large. Illustrative model calculations are presented to demonstrate this point which has important implications for turbulence modeling.

Equations governing the motion of a specific class of singularities of the Euler equation in the extended complex spatial domain are derived. Under some assumptions, it is shown how this motion is dictated by the smooth part of the complex velocity at a singular point in the unphysical domain. These results are used to relate the motion of complex singularities to the stability of steady solutions of the Euler equation. A sufficient condition for instability is conjectured. Several examples are presented to demonstrate the efficacy of this sufficient condition which include the class of elliptical flows and the Kelvin–Stuart cat’s eye.

A near-wall two-equation turbulence model of the K - epsilon type is developed for the description of high-speed compressible flows. The Favre-averaged equations of motion are solved in conjunction with modeled transport equations for the turbulent kinetic energy and solenoidal dissipation wherein a variable density extension of the asymptotically consistent near-wall model of So and co-workers is supplemented with new dilatational models. The resulting compressible two-equation model is tested in the supersonic flat plate boundary layer - with an adiabatic wall and with wall cooling - for Mach numbers as large as 10. Direct comparisons of the predictions of the new model with raw experimental data and with results from the K - omega model indicate that it performs well for a wide range of Mach numbers. The surprising finding is that the Morkovin hypothesis, where turbulent dilatational terms are neglected, works well at high Mach numbers, provided that the near wall model is asymptotically consistent. Instances where the model predictions deviate from the experiments appear to be attributable to the assumption of constant turbulent Prandtl number - a deficiency that will be addressed in a future paper.

The assumption of self-preservation permits an analytical determination of the energy decay in isotropic turbulence. Batchelor (1948), who was the first to carry out a detailed study of this problem, based his analysis on the assumption that the Loitsianskii integral is a dynamic invariant – a widely accepted hypothesis that was later discovered to be invalid. Nonetheless, it appears that the self-preserving isotropic decay problem has never been reinvestigated in depth subsequent to this earlier work. In the present paper such an analysis is carried out, yielding a much more complete picture of self-preserving isotropic turbulence. It is proven rigorously that complete self-preserving isotropic turbulence admits two general types of asymptotic solutions: one where the turbulent kinetic energy K ∼ t−1 and one where K ∼ t−α with an exponent α > 1 that is determined explicitly by the initial conditions. By a fixed-point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K ∼ t−1 power law decay is the asymptotically consistent high-Reynolds-number solution; the K ∼ t−α decay law is only achieved in the limit as t → ∞ and the turbulence Reynolds number Rt vanishes. Arguments are provided which indicate that a t−1 power law decay is the asymptotic state toward which a complete self-preserving isotropic turbulence is driven at high Reynolds numbers in order to resolve an O(R1½) imbalance between vortex stretching and viscous diffusion. Unlike in previous studies, the asymptotic approach to a complete self-preserving state is investigated which uncovers some surprising results.

Turbulence models are developed by supplementing the renormalization group (RNG) approach of Yakhot and Orszag [J. Sci. Comput. 1, 3 (1986)] with scale expansions for the Reynolds stress and production of dissipation terms. The additional expansion parameter (η≡SK/■) is the ratio of the turbulent to mean strain time scale. While low‐order expansions appear to provide an adequate description for the Reynolds stress, no finite truncation of the expansion for the production of dissipation term in powers of η suffices−terms of all orders must be retained. Based on these ideas, a new two‐equation model and Reynolds stress transport model are developed for turbulent shear flows. The models are tested for homogeneous shear flow and flow over a backward facing step. Comparisons between the model predictions and experimental data are excellent.

The ability of two-equation models to accurately predict separated flows is analyzed from a combined theoretical and computational standpoint. Turbulent flow past a backward facing step is chosen as a test case in an effort to resolve the variety of conflicting results that were published during the past decade concerning the performance of two-equation models. It is found that the errors in the reported predictions of the k-epsilon model have two major origins: (1) numerical problems arising from inadequate resolution, and (2) inaccurate predictions for normal Reynolds stress differences arising from the use of an isotropic eddy viscosity. Inadequacies in near wall modelling play a substantially smaller role. Detailed calculations are presented which strongly indicate the standard k-epsilon model - when modified with an independently calibrated anisotropic eddy viscosity - can yield surprisingly good predictions for the backstep problem.

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a963grid using Fourier collocation methods. Extensive comparisons between the direct and modelled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence – conducted on a coarse323grid – is shown to yield results that are in excellent agreement with the fine-grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.

The equilibrium structure of homogeneous turbulent shear flow is investigated from a theoretical standpoint. Existing turbulence models, in apparent agreement with physical and numerical experiments, predict an unbounded exponential time growth of the turbulent kinetic energy and dissipation rate; only the anisotropy tensor and turbulent time scale reach a structural equilibrium. It is shown that if a residual vortex stretching term is maintained in the dissipation rate transport equation, then there can exist equilibrium solutions, with bounded energy states, where the turbulence production is balanced by its dissipation. Illustrative calculations are presented for a k–ε model modified to account for net vortex stretching. The calculations indicate an initial exponential time growth of the turbulent kinetic energy and dissipation rate for elapsed times that are as large as those considered in any of the previously conducted physical or numerical experiments on homogeneous shear flow. However, vortex stretching eventually takes over and forces a production-equals-dissipation equilibrium with bounded energy states. The plausibility of this result is further supported by independent calculations of isotropic turbulence which show that when this vortex stretching effect is accounted for, a much more complete physical description of isotropic decay is obtained. It is thus argued that the inclusion of vortex stretching as an identifiable process may have greater significance in turbulence modeling than has previously been thought and that the generally accepted structural equilibrium for homogeneous shear flow, with unbounded energy growth, could be in need of re-examination.

A variety of two-equation turbulence models,including several versions of the K-epsilon model as well as the K-omega model, are analyzed critically for near wall turbulent flows from a theoretical and computational standpoint. It is shown that the K-epsilon model has two major problems associated with it: the lack of natural boundary conditions for the dissipation rate and the appearance of higher-order correlations in the balance of terms for the dissipation rate at the wall. In so far as the former problem is concerned, either physically inconsistent boundary conditions have been used or the boundary conditions for the dissipation rate have been tied to higher-order derivatives of the turbulent kinetic energy which leads to numerical stiffness. The K-omega model can alleviate these problems since the asymptotic behavior of omega is known in more detail and since its near wall balance involves only exact viscous terms. However, the modeled form of the omega equation that is used in the literature is incomplete-an exact viscous term is missing which causes the model to behave in an asymptotically inconsistent manner. By including this viscous term and by introducing new wall damping functions with improved asymptotic behavior, a new K-tau model (where tau is identical with 1/omega is turbulent time scale) is developed. It is demonstrated that this new model is computationally robust and yields improved predictions for turbulent boundary layers.

Near-wall modeling of the dissipation rate equation is investigated and its asymptotic behavior is studied in detail using a k-epsilon model. It is found that all existing modeled dissipation rate equations predict an incorrect behavior for the dissipation rate near a wall. An improvement is proposed and the resulting near-wall dissipation rate distribution is found to be similar to that given by numerical simulation data. To further validate the improved k-epsilon model, it is used to calculate flat-plate turbulent boundary-layer flows at high- as well as low-turbulence Reynolds numbers, and the results are compared with measurements, numerical simulation data, and the calculations of three different two-equation models. These comparisons show that all the models tested give essentially the same flow properties away from the wall; significant differences only occur in a region very close to the wall. In this region, the calculations of the improved k-epsilon model are in better agreement with measurements and numerical simulation data. In particular, the modeled distribution of the dissipation rate is significantly improved and a maximum is predicted at the wall instead of away from the wall. Furthermore, the improved k-epsilon model is found to be the most asymptotically consistent among the four different two-equation models examined.

Analytical methods for the development of Reynolds stress models in turbulence are reviewed in detail. Zero, one and two equation models are discussed along with second-order closures. A strong case is made for the superior predictive capabilities of second-order closure models in comparison to the simpler models. The central points are illustrated by examples from both homogeneous and inhomogeneous turbulence. A discussion of the author's views concerning the progress made in Reynolds stress modeling is also provided along with a brief history of the subject.

The pressure driven, fully developed turbulent flow of an incompressible viscous fluid in curved ducts of square cross-section is studied numerically by making use of a finite volume method. A nonlinear K -1 model is used to represent the turbulence. The results for both straight and curved ducts are presented. For the case of fully developed turbulent flow in straight ducts, the secondary flow is characterized by an eight-vortex structure for which the computed flowfield is shown to be in good agreement with available experimental data. The introduction of moderate curvature is shown to cause a substantial increase in the strength of the secondary flow and to change the secondary flow pattern to either a double-vortex or a four-vortex configuration.

The structure of the subgrid scale fields in plane channel flow has been studied at various stages of the transition process to turbulence. The residual stress and subgrid scale dissipation calculated using velocity fields generated by direct numerical simulations of the Navier-Stokes equations are significantly different from their counterparts in turbulent flows. The subgrid scale dissipation changes sign over extended areas of the channel, indicating energy flow from the small scales to the large scales. This reversed energy cascade becomes less pronounced at the later stages of transition. Standard residual stress models of the Smagorinsky type are excessively dissipative. Rescaling the model constant improves the prediction of the total (integrated) subgrid scale dissipation, but not that of the local one. Despite the somewhat excessive dissipation of the rescaled Smagorinsky model, the results of a large eddy simulation of transition on a flat-plate boundary layer compare quite well with those of a direct simulation, and require only a small fraction of the computational effort. The inclusion of non-dissipative models, which could lead to further improvements, is proposed.

A quadratic nonlinear generalization of the linear Rotta model for the slow pressure‐strain correlation of turbulence is developed for high Reynolds number flows. The model is shown to satisfy realizability and to give rise to no stable nonzero equilibrium solutions for the anisotropy tensor in the case of vanishing mean velocity gradients. In order for any model to predict a return to isotropy for all relaxational flows, it is necessary to ensure that there is no nonzero stable fixed point that attracts realizable initial conditions. Both the phase space dynamics and the temporal behavior of the model are examined and compared against experimental data for the return to isotropy problem. It is demonstrated that the quadratic model successfully captures the experimental trends which clearly exhibit nonlinear behavior. Comparisons are also made with the predictions of the linear Rotta model, the quasilinear Lumley model, and the nonlinear model of Shih, Mansour, and Moin. The simple quadratic model proposed in this study does better than the Rotta model as anticipated, and also compares quite favorably with the other more complicated nonlinear models.

The applicability of linear and quadratic models for the fluctuating pressure-strain correlations is extended to flows with a system rotation by using the intrinsic spin tensor in place of the coordinate-dependent vorticity tensor. The extension is tested here for flow in a channel rotated about its spanwise axis but is equally applicable to all other models of rotation

Reynolds stress models and traditional large-eddy simulations are reexamined with a view toward developing a combined methodology for the computation of complex turbulent flows. More specifically, an entirely new approach to time-dependent Reynolds-averaged Navier-Stokes (RANS) computations and very large-eddy simulations (VLES) is presented in which subgrid scale models are proposed that allow a direct numerical simulation (DNS) to go continuously to a RANS computation in the coarse mesh/infinite Reynolds number limit. In between these two limits, we have a large eddy simulation (LES) or VLES, depending on the level of resolution. The Reynolds stress model that is ultimately recovered in the coarse mesh/infinite Reynolds number limit has built in nonequilibrium features that make it suitable for time-dependent RANS. The fundamental technical issues associated with this new approach, which has the capability of bridging the gap between DNS, LES and RANS, are discussed in detail. Illustrative calculations are presented along with a discussion of the future implications of these results for the simulation of the turbulent flows of technological importance.

An explicit algebraic stress equation, developed by Gatski and Speziale, is used in the framework of K-epsilon formulation to predict complex aerodynamic turbulent flows. The nonequilibrium effects are modeled through coefficients that depend nonlinearly on both rotational and irrotational strains. The proposed model was implemented in the ISAAC Navier-Stokes code. Comparisons with the experimental data are presented which clearly demonstrate that explicit algebraic stress models can predict the correct response to nonequilibrium flow.

A variety of modifications to the modeled dissipation rate transport equation that have been proposed during the past two decades to account for rotational strains are examined. The models are subjected to two crucial test cases: the decay of isotropic turbulence in a rotating frame and homogeneous shear flow in a rotating frame. It is demonstrated that these modifications do not yield substantially improved predictions for these two test cases and in many instances give rise to unphysical behavior. An alternative proposal, based on the use of the tensor dissipation rate, is made for the development of improved models.

Recent work on the development of a second-order closure model for high-speed compressible flows is reviewed. This turbulent closure is based on the solution of modeled transport equations for the Favre-averaged Reynolds stress tensor and the solenoidal part of the turbulent dissipation rate. A new model for the compressible dissipation is used along with traditional gradient transport models for the Reynolds heat flux and mass flux terms. Consistent with simple asymptotic analyses, the deviatoric part of the remaining higher-order correlations in the Reynolds stress transport equations are modeled by a variable density extension of the newest incompressible models. The resulting second-order closure model is tested in a variety of compressible turbulent flows which include the decay of isotropic turbulence, homogeneous shear flow, the supersonic mixing layer, and the supersonic flat-plate turbulent boundary layer. Comparisons between the model predictions and the results of physical and numerical experiments are quite encouraging.

A recently developed K-tau model for near wall turbulent flows is applied to a variety of test cases. The turbulent flows considered include the incompressible flat plate boundary layer with adverse pressure gradients, incompressible flow past a backward facing step, and the supersonic flat plate boundary layer at zero pressure gradient. Calculations are performed for this two-equation model using an anisotropic as well as isotropic eddy-viscosity. The model predictions are shown to compare quite favorably with experimental data.

A two-equation turbulence model of the K-epsilon type was recently derived by using Renormalization Group (RNG) methods. It was later reported that this RNG based model yields substantially better predictions than the standard K-epsilon model for turbulent flow over a backward facing step - a standard test case used to benchmark the performance of turbulence models in separated flows. The improvements obtained from the RNG K-epsilon model were attributed to the better treatment of near wall turbulence effects. In contrast to these earlier claims, it is shown in this paper that the original version of the RNG K-epsilon model substantially underpredicts the reattachment point in the backstep problem. This is a deficiency that is traced to the modeling of the production of dissipation term. However, with the most recent improvements in the RNG K-epsilon model, excellent results for the backstep problem are now obtained.

The physical properties of the commonly used second-order closure models are examined theoretically for rotating turbulent flows. Comparisons are made with results which are a rigorous consequence of the Navier—Stokes equations for the problem of fully-developed turbulent channel flow in a rapidly rotating framework. It is demonstrated that all existing second-order closures yield spurious physical results for this test problem of rotating channel flow. In fact, the results obtained are shown to be substantially more unphysical than those obtained from the simpler K − ϵ K - \epsilon and K − l K - l models. Modifications in the basic structure of these second-order closure models are proposed which can alleviate this problem.

A comparison of several commonly used turbulence models (including the K–ε model and three second-order closures) is made for the test problem of homogeneous turbulent shear flow in a rotating frame. The time evolution of the turbulent kinetic energy and dissipation rate is calculated for these models and comparisons are made with previously published experiments and numerical simulations. Particular emphasis is placed on examining the ability of each model to predict equilibrium states accurately for a range of the parameter Ω/S (the ratio of the rotation rate to the shear rate). It is found that none of the commonly used second-order closure models yield substantially improved predictions for the time evolution of the turbulent kinetic energy and dissipation rate over the somewhat defective results obtained from the simpler K–ε model for the unstable flow regime. There is also a problem with the equilibrium states predicted by the various models. For example, the K–ε model erroneously yields equilibrium states that are independent of Ω/S while the Launder, Reece & Rodi model and the Shih-Lumley model predict a flow relaminarization when Ω/S > 0.39 - a result that is contrary to numerical simulations and linear spectral analyses, which indicate flow instability for at least the range 0 [les ] Ω/S [les ] 0.5. The physical implications of the results obtained from the various turbulence models considered herein are discussed in detail along with proposals to remedy the deficiencies based on a dynamical systems approach.

The scaling properties of plane homogeneous turbulent shear flows in a rotating frame are examined mathematically by a direct analysis of the Navier-Stokes equations. It is proved that two such shear flows are dynamically similar if and only if their initial dimensionless energy spectrum E star (k star, 0), initial dimensionless shear rate SK sub 0/epsilon sub 0, initial Reynolds number K squared sub 0/nu epsilon sub 0, and the ration of the rotation rate to the shear rate omega/S are identical. Consequently, if universal equilibrium states exist, at high Reynolds numbers, they will only depend on the single parameter omega/S. The commonly assumed dependence of such equilibrium states on omega/S through the Richardson number Ri=-2(omega/S)(1-2 omega/S) is proven to be inconsistent with the full Navier-Stokes equations and to constitute no more than a weak approximation. To be more specific, Richardson number similarity is shown to only rigorously apply to certain low-order truncations of the Navier-Stokes equations (i.e., to certain second-order closure models) wherein closure is achieved at the second-moment level by assuming that the higher-order moments are a small perturbation of their isotropic states. The physical dependence of rotating turbulent shear flows on omega/S is discussed in detail along with the implications for turbulence modeling.

The role that helicity fluctuations play in relationship to small-scale intermittency and coherent structures in turbulent flows is examined from a theoretical standpoint. It has been argued recently that regions of large turbulent activity are associated with regions of small local helicity and, hence, that such a correlation exists. In this paper, it will be proven that fluctuations in the local helicity are not Galilean invariant and, consequently, it is possible to have completely different local helicity fluctuations for a given fluctuating velocity field and its associated turbulence statistics. It is thus highly doubtful that local helicity fluctuations can be of any fundamental value in correlating small-scale intermittency or coherent structures.

A surface functional approximation to the general theory of nonlocal Stokesian fluids is obtained. The relationship between this theory and the existing phenomenological theories of turbulence is demonstrated. The theory is then applied to the problem of turbulent channel flow. An exact solution is obtained for the mean velocity and turbulent stresses which is shown to he consistent with existing experimental data.

The consistency of the problem of Reynolds stress closure in turbulence with the fundamental principles of modern continuum mechanics is examined. It is shown that Reynolds stress closure is completely consistent with the principles of determinism and material frame-indifference of modern continuum mechanics. The consequences that this has on turbulence modeling are briefly discussed.

The mathematical properties of the pressure-velocity and vorticity-velocity formulations of the equations of viscous flow are compared. It is shown that a vorticity-velocity formulation exists which has the interesting property that non-inertial effects only enter the problem through the implementation of initial and boundary conditions. This valuable characteristic, along with other advantages of the vorticity-velocity approach, are discussed in detail.

Recent conjectures concerning the correlation between regions of high local helicity and low dissipation are examined from a rigorous theoretical standpoint based on the Navier-Stokes equations. It is proven that only the solenoidal part of the Lamb vector w × u (which is directly tied to the nonlocal convection and stretching of vortex lines) contributes to the energy cascade in turbulence. Consequently, it is shown that regions of low dissipation can be associated with either low or high helicity—a result which disproves earlier speculations concerning this direct connection between helicity and the energy cascade. Some brief examples are given along with a discussion of the consistency of these results with the most recent computations of helicity fluctuations in incompressible turbulent flows.

A subgrid‐scale model recently derived by Yoshizawa [Phys. Fluids 29, 2152 (1986)] for use in the large‐eddy simulation of compressible turbulent flows is examined from a fundamental theoretical and computational standpoint. It is demonstrated that this model, which is only applicable to compressible turbulent flows in the limit of small density fluctuations, correlates somewhat poorly with the results of direct numerical simulations of compressible isotropic turbulence at low Mach numbers. An alternative model, based on Favre‐filtered fields, is suggested which appears to reduce these limitations.