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

Dynamical interpretation of conditional patterns

While great progress is being made in characterizing the 3-D structure of organized turbulent motions using conditional averaging analysis, there is a lack of theoretical guidance regarding the interpretation and utilization of such information. Questions concerning the significance of the structures, their contributions to various transport properties, and their dynamics cannot be answered without recourse to appropriate dynamical governing equations. One approach which addresses some of these questions uses the conditional fields as initial conditions and calculates their evolution from the Navier-Stokes equations, yielding valuable information about stability, growth, and longevity of the mean structure. To interpret statistical aspects of the structures, a different type of theory which deals with the structures in the context of their contributions to the statistics of the flow is needed. As a first step toward this end, an effort was made to integrate the structural information from the study of organized structures with a suitable statistical theory. This is done by stochastically estimating the two-point conditional averages that appear in the equation for the one-point probability density function, and relating the structures to the conditional stresses. Salient features of the estimates are identified, and the structure of the one-point estimates in channel flow is defined

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

Development of a long life thermal cell Final report

Development of long life thermal cells, environmental and electrochemical performance tests to evaluate their capability and reproductibilit

Ejection mechanisms in the sublayer of a turbulent channel

A possible model for the inception of vorticity ejections in the viscous sublayer of a turbulent rectangular channel is presented. It was shown that this part of the flow is dominated by protruding strong shear layers of z-vorticity, and it was proposed as a mechanism for their maintenance and reproduction which is essentially equivalent to that responsible for the instability of 2-D Tollmien-Schlichting waves. The efforts to isolate computationally a single structure for its study have failed up to now, since it appears that single structures decay in the absence of external forcing, but a convenient computation model was identified in the form of a long and narrow periodic computational box containing at each moment only a few structures. Further work in the identification of better reduced systems is in progress

Self similarity of two point correlations in wall bounded turbulent flows

The structure of turbulence at a height y from a wall is affected by the local mean shear at y, by the direct effect of the wall on the eddies, and by the action of other eddies close to or far from the wall. Some researchers believe that a single one of these mechanisms is dominant, while others believe that these effects have to be considered together. It is important to understand the relative importance of these effects in order to develop closure models, for example for the dissipation or for the Reynolds stress equation, and to understand the eddy structure of cross correlation functions and other measures. The specific objective was to examine the two point correlation, R sub vv, of the normal velocity component v near the wall in a turbulent channel flow and in a turbulent boundary layer. The preliminary results show that even in the inhomogeneous turbulent boundary layer, the two-point correlation function may have self similar forms. The results also show that the effects of shear and of blocking are equally important in the form of correlation functions for spacing normal to the wall. But for spanwise spacing, it was found that the eddy structure is quire different in these near flows. So any theory for turbulent structure must take both these effects into account

Sensitivity of mixing layers to three-dimensional forcing

It is well known that turbulent mixing layers are dominated by large scale, fairly coherent structures, and that these structures are related to the stability characteristics of the flow. These facts have led researchers to attempt controlling such flows by selectively forcing certain unstable modes, which can in addition have the effect of suppressing other modes. Much of the work on controlling the mixing layer has relied on forcing 2-D instabilities. The results of forcing 3-D instabilities are addressed. The objectives of the work are twofold: to understand how a mixing layer responds to 3-D perturbations, and to test the validity of an amplitude expansion in predicting the mixing layer development. The amplitude expansion could be very useful in understanding and predicting the 3-D response of the flow to a variety of initial conditions

Dynamics of coherent structures in a plane mixing layer

An incompressible, time developing 3-D mixing layer with idealized initial conditions was simulated numerically. Consistent with the suggestions from experimental measurements, the braid region between the dominant spanwise vortices or rolls develops longitudinal vortices or ribs, which are aligned upstream and downstream of a roll and produce spanwise distortion of the rolls. The process by which this distortion occurs is explained by studying a variety of quantities of dynamic importance (e.g., production of enstrophy, vortex stretching). Other quantities of interest (dissipation, helicity density) are also computed and discussed. The currently available simulation only allows the study of the early evolution (before pairing) of the mixing layer. New simulations in progress will relieve this restriction

Hyperbolic outer billiards : a first example

We present the first example of a hyperbolic outer billiard. More precisely we construct a one parameter family of examples which in some sense correspond to the Bunimovich billiards.Comment: 11 pages, 8 figures, to appear in Nonlinearit