No abstract

A spatially developing supersonic adiabatic flat plate boundary layer flow (at M-infinity=2.25 and Re(theta)approximate to4000) is analyzed by means of direct numerical simulation. The numerical algorithm is based on a mixed weighted essentially nonoscillatory compact-difference method for the three-dimensional Navier-Stokes equations. The main objectives are to assess the validity of Morkovin's hypothesis and Reynolds analogies, and to analyze the controlling mechanisms for turbulence production, dissipation, and transport. The results show that the essential dynamics of the investigated turbulent supersonic boundary layer flow closely resembles the incompressible pattern. The Van Driest transformed mean velocity obeys the incompressible law-of-the-wall, and the mean static temperature field exhibits a quadratic dependency upon the mean velocity, as predicted by the Crocco-Busemann relation. The total temperature has been found not to be precisely uniform, and total temperature fluctuations are found to be non-negligible. Consistently, the turbulent Prandtl number is not unity, and it varies between 0.7 and 0.8 in the outer part of the boundary layer. Nonetheless, a modified strong Reynolds analogy is still verified. In agreement with the low Mach number results, the streamwise velocity component and the temperature are only weakly anti-correlated. The turbulent kinetic energy budget also shows similarities with the incompressible case provided all terms of the equation are properly scaled; indeed, the leading compressibility contributions are negligible throughout the boundary layer. (C) 2004 American Institute of Physics

We study the organization of turbulence in supersonic boundary layers through large-scale direct numerical simulations (DNS) at M-infinity = 2, and momentum-thickness Reynolds number up to Re-delta 2 approximate to 3900 (corresponding to Re-tau approximate to 1120) which significantly extend the current envelope of DNS in the supersonic regime. The numerical strategy relies on high-order, non-dissipative discretization of the convective terms in the Navier Stokes equations, and it implements a recycling/rescaling strategy to stimulate the inflow turbulence. Comparison of the velocity statistics up to fourth order shows nearly exact agreement with reference incompressible data, provided the momentum-thickness Reynolds number is matched, and provided the mean velocity and the velocity fluctuations are scaled to incorporate the effects of mean density variation, as postulated by Morkovin's hypothesis. As also found in the incompressible regime, we observe quite a different behaviour of the second-order flow statistics at sufficiently large Reynolds number, most of which show the onset of a range with logarithmic variation, typical of 'attached' variables, whereas the wall-normal velocity exhibits a plateau away from the wall, which is typical of 'detached' variables. The modifications of the structure of the flow field that underlie this change of behaviour are highlighted through visualizations of the velocity and temperature fields, which substantiate the formation of large jet-like and wake-like motions in the outer part of the boundary layer. It is found that the typical size of the attached eddies roughly scales with the local mean velocity gradient, rather than being proportional to the wall distance, as happens for the wall-detached variables. The interactions of the large eddies in the outer layer with the near-wall region are quantified through a two-point amplitude modulation covariance, which characterizes the modulating action of energetic outer-layer eddies

We study passive scalars in turbulent plane channels at computationally high Reynolds number, thus allowing us to observe previously unnoticed effects. The mean scalar profiles are found to obey a generalized logarithmic law which includes a linear correction term in the whole lower half-channel, and they follow a universal parabolic defect profile in the core region. This is consistent with recent findings regarding the mean velocity profiles in channel flow. The scalar variances also exhibit a near universal parabolic distribution in the core flow and hints of a sizeable log layer, unlike the velocity variances. The energy spectra highlight the formation of large scalar-bearing eddies with size proportional to the channel height which are caused by a local production excess over dissipation, and which are clearly visible in the flow visualizations. Close correspondence of the momentum and scalar eddies is observed, with the main difference being that the latter tend to form sharper gradients, which translates into higher scalar dissipation. Another notable Reynolds number effect is the decreased correlation of the passive scalar field with the vertical velocity field, which is traced to the reduced effectiveness of ejection events.

We review numerical methods for direct numerical simulation (DNS) and large-eddy simulation (LES) of turbulent compressible flow in the presence of shock waves. Ideal numerical methods should be accurate and free from numerical dissipation in smooth parts of the flow, and at the same time they must robustly capture shock waves without significant Gibbs ringing, which may lead to nonlinear instability. Adapting to these conflicting goals leads to the design of strongly nonlinear numerical schemes that depend on the geometrical properties of the solution. For low-dissipation methods for smooth flows, numerical stability can be based on physical conservation principles for kinetic energy and/or entropy. Shock-capturing requires the addition of artificial dissipation, in more or less explicit form, as a surrogate for physical viscosity, to obtain nonoscillatory transitions. Methods suitable for both smooth and shocked flows are discussed, and the potential for hybridization is highlighted. Examples of the application of advanced algorithms to DNS/LES of turbulent, compressible flows are presented