ALBERT

Description: In the present work, the high Reynolds number flow past an inclined plate with a splitter plate placed in its wake is considered numerically. A numerical conformal mapping technique is employed to transform the two-plate system into the same number of cylinders: the flow field is assumed to be two-dimensional. The vortex shedding from the inclined plate is modelled using the discrete vortex method. It is shown that the splitter plate has a profound effect on the development of the flow over a range of values of a suitably defined offset parameter and for a range of positions of the leading edge of the splitter plate. The acoustic field is also calculated and the spectrum reflects the flow results.

Description: Spectral local isotropy tests are applied to direct numerical simulation data, mainly at the centerline of a fully developed turbulent channel flow. Despite the small Reynolds number of the simulation, the high-wavenumber behavior of velocity and vorticity spectra is consistent with local isotropy. This consistency is verified by the relationship between streamwise wavenumber spectra and spanwise wavenumber spectra. The high-wavenumber behavior of the pressure spectrum is also consistent with local isotropy and compares favorably with the calculation of Batchelor (1951), which assumes isotropy and joint normality of the velocity field at two points in space. The latter assumption is validated by the shape but not the magnitude of the quadruple correlation of the streamwise velocity fluctuation at small separations. There is only partial support for local spectral isotropy away from the centerline as the magnitude of the mean strain rate increases.

Description: A Navier-Stokes equations solver based on a pressure correction method with a pressure-staggered mesh and calculations of separated three-dimensional flows are presented. It is shown that the velocity pressure decoupling, which occurs when various pressure correction algorithms are used for pressure-staggered meshes, is caused by the ill-conditioned discrete pressure correction equation. The use of a partial differential equation for the incremental pressure eliminates the velocity pressure decoupling mechanism by itself and yields accurate numerical results. Example flows considered are a three-dimensional lid driven cavity flow and a laminar flow through a 90 degree bend square duct. For the lid driven cavity flow, the present numerical results compare more favorably with the measured data than those obtained using a formally third order accurate quadratic upwind interpolation scheme. For the curved duct flow, the present numerical method yields a grid independent solution with a very small number of grid points. The calculated velocity profiles are in good agreement with the measured data.

Description: A database obtained from direct numerical simulation of a turbulent channel flow is analyzed to extract the streamwise component of the propagation velocity V of velocity, vorticity, and pressure fluctuations from their space-time correlations. A surprising result is that V is approximately the same as the local mean velocity for most of the channel, except for the near-wall region. For y(+) less than 15, V is virtually constant, implying that perturbations of all flow variables propagate like waves near the wall. In this region V is 55 percent of the centerline velocity Uc for velocity and vorticity perturbations and 75 percent of U sub c for pressure perturbations. This is equal to U at y(+) = 15 for velocity and vorticity perturbations, and equal to U at y(+) = 20 for pressure perturbations, indicating that the dynamics of the nearwall turbulence is controlled by turbulence structures present near y(+) about 15-20. Scale dependence of V is also examined by analyzing the bandpass-filtered flow fields. This paper contains comprehensive documentation on the propagation velocities, which should prove useful in the evaluation of Taylor's hypothesis.

Description: Direct numerical simulation data for the lateral velocity derivative delta(u)/delta(y) at the centerline of a fully developed turbulent channel flow provide reasonable support for Wyngaard's analysis of the error involved in measuring this quantity using parallel hot wires. Numerical data in the wall region of the channel flow also provide a useful indication of how to select the separation between the wires. Justification for this choice is obtained by comparing several measured statistics of delta(u)/delta(y) with the corresponding numerical data.

Description: One of the key technical elements in NASA's high speed research program is reducing the noise level to meet the federal noise regulation. The dominant noise source is associated with the supersonic jet discharged from the engine exhaust system. Whereas the turbulence mixing is largely responsible for the generation of the jet noise, a broadband shock-associated noise is also generated when the nozzle operates at conditions other than its design. For both mixing and shock noise components, because the source of the noise is embedded in the jet plume, one can expect that jet noise can be predicted from the jet flowfield computation. Mani et al. developed a unified aerodynamic/acoustic prediction scheme by applying an extension of Reichardt's aerodynamic model to compute turbulent shear stresses which are utilized in estimating the strength of the noise source. Although this method produces a fast and practical estimate of the jet noise, a modification by Khavaran et al. has led to an improvement in aerodynamic solution. The most notable feature in this work is that Reichardt's model is replaced with the computational fluid dynamics (CFD) solution of Reynolds-averaged Navier-Stokes equations. The major advantage of this work is that the essential, noise-related flow quantities such as turbulence intensity and shock strength can be better predicted. The predictions were limited to a shock-free design condition and the effect of shock structure on the jet mixing noise was not addressed. The present work is aimed at investigating this issue. Under imperfectly expanded conditions the existence of the shock cell structure and its interaction with the convecting turbulence structure may not only generate a broadband shock-associated noise but also change the turbulence structure, and thus the strength of the mixing noise source. Failure in capturing shock structures properly could lead to incorrect aeroacoustic predictions.

Description: Low-Reynolds-number effects are observed in the inner region of a fully developed turbulent channel flow, using data obtained either from experiments or by direct numerical simulations. The Reynolds-number influence is observed on the turbulence intensities and to a lesser degree on the average production and dissipation of the turbulent energy. In the near-wall region, the data confirm Wei and Willmarth's (1989) conclusion that the Reynolds stresses do not scale on wall variables. One of the reasons proposed to account for this behavior, namely, the 'geometry' effect or direct interaction between inner regions on opposite walls, was investigated in some detail by introducing temperature at one of the walls, both in experiment and simulation. Although the extent of penetration of thermal excursions into the opposite side of the channel can be significant at low Reynolds numbers, the contribution these excursions make to the Reynolds shear stress and the spanwise vorticity in the opposite wall region is negligible. In the inner region, spectra and cospectra of the velocity fluctuations u and v change rapidly with the Reynolds number, the variations being mainly confined to low wavenumbers in the u spectrum.

Description: The fine-scale structure of turbulence in a fully developed turbulent duct flow is examined by considering the 3D velocity derivative field obtained from direct numerical simulations at two relatively small Reynolds numbers. The magnitudes of all mean-square derivatives (normalized by wall variables) increase with the Reynolds number, the increase being largest at the wall. These magnitudes are not consistent with the assumption of local isotropy except perhaps near the duct center-line. When the assumption of local isotropy is relaxed to one of local axisymmetry, or invariance with respect to rotation about a coordinate axis (here chosen in the streamwise direction), satisfactory agreement is indicated by the data outside the wall region. Support for axisymmetry is demonstrated by anisotropy invariant maps of the dissipation and vorticity tensors.

Description: Numerical calculation of a three dimensional turbulent flow of a jet in a crossflow using a multiple time scale turbulence model is presented. The turbulence in the forward region of the jet is in a stronger inequilibrium state than that in the wake region of the jet, while the turbulence level in the wake region is higher than that in the front region. The calculated flow and the concentration fields are in very good agreement with the measured data, and it indicated that the turbulent transport of mass, concentration, and momentum is strongly governed by the inequilibrium turbulence. The capability of the multiple time scale turbulence model to resolve the inequilibrium turbulence field is also discussed.