Search NASA⌕ Search

Engineering topics

Wilcox, D. C.

Publications and source records attributed to Wilcox, D. C..

Transition region modeling for compressible flow

The ongoing effort at NASA Langley Research Center aimed at developing transitional zone models for Reynolds-averaged Navier-Stokes calculations is described. Results from three models that we have examined - two zero-equation models and a two-equation model of the k-omega type - are presented for high-speed flows. The use of analytical and computational tools to provide guidance and to provide the transitional flow-field data bases necessary to develop such models is discussed. Some results from validation of these tools are given; qualitative comparisons between modeled and theoretical/computational representations of the transitional zone are provided.

Chang, C.-L.↗

Algorithm for rapid integration of turbulence model equations on parabolic regions

While developing a three-dimensional boundary layer program using a standard parabolic matching scheme, the author has found computing time with the Wilcox-Rubesin (1979) two-equation turbulence model to be very lengthy. The long computing time occurs because converged solutions are possible only when very small streamwise steps are taken. The proposed remedy reduces computing time by increasing the maximum permissible step size.

Wilcox, D. C.↗

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

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

Wilcox, D. C.↗

Transitional Boundary-Layer Solutions Using a Mixing-Length and a Two-Equation Turbulence Model

Boundary-layer solutions were obtained using the conventional two-layer mixing-length turbulence model and the Wilcox-Traci two-equation model of turbulence. Both flatplate and blunt-body geometries were considered. The most significant result of the study is development of approximations for the two-equation model which permit streamwise stepsize comparable to that used in mixing-length computations. Additionally, a set of model-equation boundary conditions derived which apply equally well to both flat-plate and blunt-body geometries. Solutions obtained with the two-equations turbulence model are compared with experimental data and/or corresponding solutions obtained using the mixing-length model. Agreement is satisfactory for flat-plate boundary layers but not for blunt body boundary layers.

Anderson, E. C.↗

Vorticity interaction effects on blunt bodies

Numerical solutions of the viscous shock layer equations governing laminar and turbulent flows of a perfect gas and radiating and nonradiating mixtures of perfect gases in chemical equilibrium are presented for hypersonic flow over spherically blunted cones and hyperboloids. Turbulent properties are described in terms of the classical mixing length. Results are compared with boundary layer and inviscid flowfield solutions; agreement with inviscid flowfield data is satisfactory. Agreement with boundary layer solutions is good except in regions of strong vorticity interaction; in these flow regions, the viscous shock layer solutions appear to be more satisfactory than the boundary layer solutions. Boundary conditions suitable for hypersonic viscous shock layers are devised for an advanced turbulence theory.

Anderson, E. C.↗

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

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

Streamline curvature effects on turbulent boundary layers

A theoretical tool has been developed for predicting, in a nonempirical manner, effects of streamline curvature and coordinate-system rotation on turbulent boundary layers. The second-order closure scheme developed by Wilcox and Traci has been generalized for curved streamline flow and for flow in a rotating coordinate system. A physically based straightforward argument shows that curvature/rotation primarily affects the turbulent mixing energy; the argument yields suitable curvature/rotation terms which are added to the mixing-energy equation. Singular-perturbation solutions valid in the wall layer of a curved-wall boundary layer and a fully developed rotating channel flow demonstrate that, with the curvature/rotation terms, the model predicts the curved-wall and the rotating coordinate system laws of the wall. Results of numerical computations of curved-wall boundary layers and of rotating channel flow show that curvature/rotation effects can be computed accurately with second-order closure.

Wilcox, D. C.↗

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