The boundary layer equations for plane, incompressible, and steady flow are 

$$\matrix{ {u{{\partial u} \over {\partial x}} + v{{\partial u} \over {\partial y}} = - {1 \over \varrho }{{\partial p} \over {\partial x}} + v{{{\partial ^2}u} \over {\partial {y^2}}},} \cr {0 = {{\partial p} \over {\partial y}},} \cr {{{\partial u} \over {\partial x}} + {{\partial v} \over {\partial y}} = 0.} \cr }$$

Boundary Layer Theory

An introduction to heat transfer


 The spectral ratio between horizontal and vertical components (H/V ratio) of microtremors measured at the ground surface has been used to estimate fundamental periods and amplification factors of a site, although this technique lacks theoretical background. The aim of this article is to formulate the H/V technique in terms of the characteristics of Rayleigh and Love waves, and to contribute to improve the technique. The improvement includes use of not only peaks but also troughs in the H/V ratio for reliable estimation of the period and use of a newly proposed smoothing function for better estimation of the amplification factor. The formulation leads to a simple formula for the amplification factor expressed with the H/V ratio. With microtremor data measured at 546 junior high schools in 23 wards of Tokyo, the improved technique is applied to mapping site periods and amplification factors in the area.

Ground Motion Characteristics Estimated from Spectral Ratio between Horizontal and Verticcl Components of Mietremors.

The Structure of Turbulent Shear Flow By Dr. A. A. Townsend. Pp. xii + 315. 8¾ in. × 5½ in. (Cambridge: At the University Press.) 40s.

/pdf/the-structure-of-turbulent-shear-flow-2yf3szhaqb.pdf

The Structure of Turbulent Shear Flow

A direct numerical simulation of incompressible channel flow at a friction Reynolds number ( ) of 5186 has been performed, and the flow exhibits a number of the characteristics of high-Reynolds-number wall-bounded turbulent flows. For example, a region where the mean velocity has a logarithmic variation is observed, with von Karman constant . There is also a logarithmic dependence of the variance of the spanwise velocity component, though not the streamwise component. A distinct separation of scales exists between the large outer-layer structures and small inner-layer structures. At intermediate distances from the wall, the one-dimensional spectrum of the streamwise velocity fluctuation in both the streamwise and spanwise directions exhibits dependence over a short range in wavenumber . Further, consistent with previous experimental observations, when these spectra are multiplied by (premultiplied spectra), they have a bimodal structure with local peaks located at wavenumbers on either side of the range.

Direct numerical simulation of turbulent channel flow up to

A spatially developing supersonic adiabatic flat plate boundary layer flow (at M∞=2.25 and Reθ≈4000) 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. Consiste...

Direct numerical simulation and analysis of a spatially evolving supersonic turbulent boundary layer at M=2.25

Conservative Hybrid Compact-WENO Schemes for Shock-Turbulence Interaction

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.

Numerical Methods for High-Speed Flows

The interaction of a spatially developing adiabatic boundary layer flow at M∞=2.25 and Reθ=3725 with an impinging oblique shock wave (β=33.2°) is analyzed by means of direct numerical simulation of the compressible Navier-Stokes equations. Under the selected flow conditions the incoming boundary layer undergoes mild separation due to the adverse pressure gradient. Coherent structures are shed near the average separation point and the flow field exhibits large-scale low-frequency unsteadiness. The formation of the mixing layer is primarily responsible for the amplification of turbulence, which relaxes to an equilibrium state past the interaction. Complete equilibrium is attained in the inner part of the boundary layer, while in the outer region the relaxation process is incomplete. Far from the interaction zone, turbulence exhibits a universal behavior and it shows similarities with the incompressible case. The interaction of the coherent structures with the incident shock produces acoustic waves that prop...

Direct numerical simulation of impinging shock wave/turbulent boundary layer interaction at M=2.25

The high-Reynolds-number behaviour of the canonical incompressible turbulent channel flow is investigated through large-scale direct numerical simulation (DNS). A Reynolds number is achieved ( , where is the channel half-height, and is the viscous length scale) at which theory predicts the onset of phenomena typical of the asymptotic Reynolds number regime, namely a sensible layer with logarithmic variation of the mean velocity profile, and Kolmogorov scaling of the velocity spectra. Although higher Reynolds numbers can be achieved in experiments, the main advantage of the present DNS study is access to the full three-dimensional flow field. Consistent with refined overlap arguments (Afzal & Yajnik, J. Fluid Mech. vol. 61, 1973, pp. 23–31; Jimenez & Moser, Phil. Trans. R. Soc. Lond. A, vol. 365, 2007, pp. 715–732), our results suggest that the mean velocity profile never achieves a truly logarithmic profile, and the logarithmic diagnostic function instead exhibits a linear variation in the outer layer whose slope decreases with the Reynolds number. The extrapolated value of the von Karman constant is . A near logarithmic layer is observed in the spanwise velocity variance, as predicted by Townsend’s attached eddy hypothesis, whereas the streamwise variance seems to exhibit a shoulder, perhaps being still affected by low-Reynolds-number effects. Comparison with previous DNS data at lower Reynolds number suggests enhancement of the imprinting effect of outer-layer eddies onto the near-wall region. This mechanisms is associated with excess turbulence kinetic energy production in the outer layer, and it reflects in flow visualizations and in the streamwise velocity spectra, which exhibit sharp peaks in the outer layer. Associated with the outer energy production site, we find evidence of a Kolmogorov-like inertial range, limited to the spanwise spectral density of , whereas power laws with different exponents are found for the other spectra. Finally, arguments are given to explain the ‘odd’ scaling of the streamwise velocity variances, based on the analysis of the kinetic energy production term.

Sergio Pirozzoli

Papers

Direct numerical simulation and analysis of a spatially evolving supersonic turbulent boundary layer at M=2.25

Conservative Hybrid Compact-WENO Schemes for Shock-Turbulence Interaction

Numerical Methods for High-Speed Flows

Direct numerical simulation of impinging shock wave/turbulent boundary layer interaction at M=2.25

Velocity statistics in turbulent channel flow up to