Numerical simulation of wall-bounded turbulent shear flows

Developments in three dimensional, time dependent numerical simulation of turbulent flows bounded by a wall are reviewed. Both direct and large eddy simulation techniques are considered within the same computational framework. The computational spatial grid requirements as dictated by the known structure of turbulent boundary layers are presented. The numerical methods currently in use are reviewed and some of the features of these algorithms, including spatial differencing and accuracy, time advancement, and data management are discussed. A selection of the results of the recent calculations of turbulent channel flow, including the effects of system rotation and transpiration on the flow are included

Large eddy simulation of turbulent channel flow: ILLIAC 4 calculation

The three-dimensional time dependent equations of motion were numerically integrated for fully-developed turbulent channel flow. A large scale flow field was obtained directly from the solution of these equations, and small scale field motions were simulated through an eddy viscosity model. The calculations were carried out on the ILLIAC 4 computer. The computed flow patterns show that the wall layer consists of coherent structures of low speed and high speed streaks alternating in the spanwise direction. These structures were absent in the regions away from the wall. Hot spots, small localized regions of very large turbulent shear stress, were frequently observed. The profiles of the pressure velocity-gradient correlations show a significant transfer of energy from the normal to the spanwise component of turbulent kinetic energy in the immediate neighborhood of the wall ('the splatting effect')

Numerical investigation of turbulent channel flow

Fully developed turbulent channel flow was simulated numerically at Reynolds number 13800, based on centerline velocity and channel halt width. The large-scale flow field was obtained by directly integrating the filtered, three dimensional, time dependent, Navier-Stokes equations. The small-scale field motions were simulated through an eddy viscosity model. The calculations were carried out on the ILLIAC IV computer with up to 516,096 grid points. The computed flow field was used to study the statistical properties of the flow as well as its time dependent features. The agreement of the computed mean velocity profile, turbulence statistics, and detailed flow structures with experimental data is good. The resolvable portion of the statistical correlations appearing in the Reynolds stress equations are calculated. Particular attention is given to the examination of the flow structure in the vicinity of the wall

The structure of the vorticity field in turbulent channel flow. Part 1: Analysis of instantaneous fields and statistical correlations

An investigation into the existence of hairpin vortices in turbulent channel flow is conducted using a database generated by the large eddy simulation technique. It is shown that away from the wall the distribution of the inclination angle of vorticity vector attains its maximum at about 45 deg to the wall. Two point correlations of velocity and vorticity fluctuations strongly support a flow model consisting of vortical structures inclined at 45 deg to the wall. The instantaneous vorticity vectors plotted in planes inclined at 45 deg show that the flow contains an appreciable number of hairpins. Vortex lines are used to display the three dimensional structure of hairpins, which are shown to be generated from deformation of transverse vortex filaments

The structure of the vorticity field in turbulent channel flow. Part 2: Study of ensemble-averaged fields

Several conditional sampling techniques are applied to a data base generated by large-eddy simulation of turbulent channel flow. It is shown that the bursting process is associated with well-organized horseshoe vortices inclined at about 45 deg. to the wall. These vortical structures are identified by examining the vortex lines of three-dimensional, ensemble averaged vorticity fields. Two distinct horseshoe-shaped vortices corresponding to the sweep and ejection events are detected. These vortices are associated with high Reynolds shear stress and hence make a significant contribution to turbulent energy production. The dependency of the ensemble averaged vortical structures on the detection criteria, and the question of whether this ensemble-averaged structure is an artifact of the ensemble averaging process are examined. The ensemble-averaged pattern of these vortical structures that emerge from the analysis could provide the basis for a hypothetical model of the organized structures of wall-bounded shear flows

Direct numerical simulation of curved turbulent channel flow

Low Reynolds number, mildly curved, turbulent channel flow has been simulated numerically without subgrid scale models. A new spectral numerical method developed for this problem was used, and the computations were performed with 2 million degrees of freedom. A variety of statistical and structural information has been extracted from the computed flow fields. These include mean velocity, turbulence stresses, velocity skewness, and flatness factors, space time correlations and spectra, all the terms in the Reynolds stress balance equations, and contour and vector plots of instantaneous velocity fields. The effects of curvature on this flow were determined by comparing the concave and convex sides of the channel. The observed effects are consistent with experimental observations for mild curvature. The most significant difference in the turbulence statistics between the concave and convex sides was in the Reynolds shear stress. This was accompanied by significant differences in the terms of the Reynolds shear stress balance equations. In addition, it was found that stationary Taylor-Gortler vortices were present and that they had a significant effect on the flow by contributing to the mean Reynolds shear stress, and by affecting the underlying turbulence

Evolution of a curved vortex filament into a vortex ring

The deformation of a hairpin-shaped vortex filament under self-induction and in the presence of shear is studied numerically using the Biot-Savart law. It is shown that the tip region of an elongated hairpin vortex evolves into a vortex ring and that the presence of mean shear impedes the process. Evolution of a finite-thickness vortex sheet under self-induction is also investigated using the Navier-Stokes equations. The layer evolves into a hairpin vortex which in turn produces a vortex ring of high Reynolds stress content. These results indicate a mechanism for the generation of ring vortices in turbulent shear flows, and a link between the experimental and numerical observation of hairpin vortices and the observation of ring vortices in the outer regions of turbulent boundary layers

Reynolds-stress and dissipation rate budgets in a turbulent channel flow

The budgets for the Reynolds stresses and for the dissipation rate of the turbulence kinetic energy are computed using direct simulation data of a turbulent channel flow. The budget data reveal that all the terms in the budget become important close to the wall. For inhomogeneous pressure boundary conditions, the pressure-strain term is split into a return term, a rapid term, and a Stokes term. The Stokes term is important close to the wall. The rapid and return terms play different roles depending on the component of the term. A split of the velocity pressure-gradient term into a redistributive term and a diffusion term is proposed, which should be simpler to model. The budget data is used to test existing closure models for the pressure-strain term, the dissipation rate, and the transport rate. In general, further work is needed to improve the models

Computation of turbulent flows over backward-facing step

A numerical method for computing incompressible turbulent flows is presented. The method is tested by calculating laminar recirculating flows and is applied in conjunction with a modified Kappa-epsilon model to compute the flow over a backward-facing step. In the laminar regime, the computational results are in good agreement with the experimental data. The turbulent flow study shows that the reattachment length is underpredicted by the standard Kappa-epsilon model. The addition of a term to the standard model that accounts for the effects of rotation on turbulent flow improves the results in the recirculation region and increases the computed reattachment length

Sampling inhomogeneous turbulent fields

The reconstruction of an inhomogeneous random process from a finite number of discrete samples can be performed in terms of the Karhunen-Loeve (KL) expansion for that process. The n(th) eigenfunction has n - 1 zero crossings which are the sampling points for the inhomogeneous process. The rapid variation of the KL eigenfunctions makes it unnecessary to have a high density of sampling (or grid points) near the wall. However, this result should not be construed to indicate that with spectral simulations significantly fewer grid points are required with the KL expansion as compared to other orthogonal expansions. Moin and Moser (1989) have shown that the advantage of the KL expansion over Chebychev expansion rapidly diminishes when high percentage (say 90 percent) energy recovery is demanded