TRPIV investigation of space-time correlation in turbulent flows over flat and wavy walls

The present experimental work is devoted to in- vestigate a new space-time correlation model for the turbulent boundary layer over a flat and a wavy walls. A turbulent boundary layer flow at Reo = 2460 is measured by tomographic time-resolved particle image velocimetry （Tomo-TRPIV）. The space-time correlations of instantaneous streamwise fluctuation velocity are calculated at 3 different wall-normal locations in logarithmic layer. It is found that the scales of coherent structure increase with moving far away from the wall. The growth of scales is a manifestation of the growth of prevalent coherent structures in the turbulent boundary layer like hairpin vortex or hairpin packets when they lift up. The resulting contours of the space-time correlation exhibit elliptic-like shapes rather than straight lines. It is suggested that, instead of Taylor hypothesis, the elliptic model of the space-time correlation is valid for the wallbounded turbulent flow over either a flat wall or a wavy wall. The elliptic iso-correlation curves have a uniform preferred orientation whose slope is determined by the convection velocity. The convection velocity derived from the space-time correlation represents the velocity at which the large-scale eddies carry small-scale eddies. The sweep velocity rep- resents the distortions of the small-scale eddies and is intimately associated with the fluctuation velocity in the logarithmic layer of turbulent boundary layers. The nondimensionalized correlation curves confirm that the elliptic model is more proper for approximating the space-time correlation than Taylor hypothesis, because the latter can not embody the small-scale motions which have non-negligible distortions. A second flow over a wavy wall is also recorded using TRPIV. Due to the combined effect of shear layers and the adverse pressure gradient, the space-time correlation does not show an elliptic-like shape at some specific heights over the wavy wall, but in the outer region of the wavy wallbounded flow, the elliptic model remai     

Numerical study of turbulent flows over wavy boundaries

The fully developed turbulent flows over wavy boundaries are investigated by means of thek-ε model. Predicted flow characteristics over rigid wavy walls are in good agreement with the vailable experimental data. Moreover drag reduction has been found in a 2-dimensional channel with periodical wavy walls. The energy input from turbulent wind to regular waves is also studied in the paper by the same turbulence model with carefully posed boundary conditions at wind-wave interface. Better agreement has been obtained in the predication of the growth rates of wind waves as compared with the previous theoretical and numerical results. The project supported by the National Natural Science Foundation of China.     

Numerical simulation of fully-developed compressible flows over wavy surfaces

Rough surfaces are common on high-speed vehicles, for example on heat shields, but compressibility is not usually taken into account in the flow modelling other than through the mean density. In the present study, supersonic fully-developed turbulent rough wall channel flows are simulated using direct numerical simulation to investigate whether strong compressibility effects significantly alter the mean flow and turbulence properties across the channel. The simulations were run for three different Mach numbers M = 0.3, 1.5 and 3.0 over a range of wall amplitude-to-wavelength ratios from 0.01 to 0.08, corresponding to transitionally and fully rough cases respectively. The velocity deficit values are found to decrease with increasing Mach number. It is also found that at Mach 3.0 significant differences occur in the mean flow and turbulence statistics throughout the channel and not just in a roughness sublayer. These differences are found to be due to the presence of strong shock waves created by the peaks of the roughness elements.     

Adverse pressure gradient turbulent flows over rough walls

An experimental study of a fully developed turbulent channel flow and an adverse pressure gradient (APG) turbulent channel flow over smooth and rough walls has been performed using a particle image velocimetry (PIV) technique. The rough walls comprised two-dimensional square ribs of nominal height, k = 3 mm and pitch, p = 2k, 4k and 8k. It was observed that rib roughness enhanced the drag characteristics, and the degree of enhancement increased with increasing pitch. Similarly, rib roughness significantly increased the level of turbulence production, Reynolds stresses and wall-normal transport of turbulence kinetic energy and Reynolds shear stress well beyond the roughness sublayer. On the contrary, the distributions of the eddy viscosity, mixing length and streamwise transport of turbulence kinetic energy and Reynolds shear stress were reduced by wall roughness, especially in the outer layer. Adverse pressure gradient produced a further reduction in the mean velocity (in comparison to the results obtained in the parallel section) but increased the wall-normal extent across which the mean flow above the ribs is spatially inhomogeneous in the streamwise direction. APG also reinforced wall roughness in augmenting the equivalent sand grain roughness height. The combination of wall roughness and APG significantly increased turbulence production and Reynolds stresses except in the immediate vicinity of the rough walls. The transport velocities of the turbulence kinetic energy and Reynolds shear stress were also augmented by APG across most part of the rough-wall boundary layer. Further, APG enhanced the distributions of the eddy viscosity across most of the boundary layer but reduced the mixing length outside the roughness sublayer.     

Scattering of acoustic waves from a point source over an impedance wedge

In this work we study diffraction of a spherical acoustic wave due to a point source, by an impedance wedge In the exterior of the wedge the acoustic pressure satisfies the stationary wave (Helmholtz) equation and classical impedance boundary conditions on two faces of the wedge, as well as Meixner’s condition at the edge and the radiation conditions at infinity. Solution of the boundary value problem is represented by a Weyl type integral and its asymptotic behavior is discussed. On this way, we derive various components in the far field interpreting them accordingly and discussing their physical meaning.     

Relationship between wall shear stresses and streamwise vortices in turbulent flows over wavy boundaries

The relationship between wall shear stresses and near-wall streamwise vortices is investigated via a direct numerical simulation(DNS) of turbulent flows over a wavy boundary with traveling-wave motion. The results indicate that the wall shear stresses are still closely related to the near-wall streamwise vortices in the presence of a wave. The wave age and wave phase significantly affect the distribution of a two-point correlation coefficient between the wall shear stresses and streamwise vorticity. For the slow wave case of c/Um = 0.14, the correlation is attenuated above the leeward side while the distribution of correlation function is more elongated and also exhibits a larger vertical extent above the crest. With respect to the fast wave case of c/U_m=1.4, the distribution of the correlation function is recovered in a manner similar to that in the flat-wall case. In this case, the maximum correlation coefficient exhibits only slight differences at different wave phases while the vertical distribution of the correlation function depends on the wave phase.     

Turbulence production in flow over a wavy wall

Measurements of the spatial and time variation of two components of the velocity have been made over a sinusoidal solid wavy boundary with a height to length ratio of 2a/λ = 0.10 and with a dimensionless wave number of α⁺ = (2π/λ)(v/u ^?) = 0.02. For these conditions, both intermittent and time-mean flow reversals are observed near the troughs of the waves. Statistical quantities that are determined are the mean streamwise and normal velocities, the root-meansquare of the fluctuations of the streamwise and normal velocities, and the Reynolds shear stresses. Turbulence production is calculated from these measurements. The flow is characterized by an outer flow and by an inner flow extending to a distance of about α^?1 from the mean level of the surface. Turbulence production in the inner region is fundamentally different from flow over a flat surface in that it is mainly associated with a shear layer that separates from the back of the wave. Flow close to the surface is best described by an interaction between the shear layer and the wall, which produces a retarded zone and a boundary-layer with large wall shear stresses. Measurements of the outer flow compare favorably with measurements over a flat wall if velocities are made dimensionless by a friction velocity defined with a shear stress obtained by extrapolating measurements of the Reynolds stress to the mean levels of the surface (rather than from the drag on the wall).     

Supersonic gas flow over a wavy surface

Exact numerical solutions of the problem are obtained for a compressible viscous heat-conducting gas flowing over a slightly wavy surface. The effect of the free-stream parameters on the stresses, the temperature fields and the heat fluxes at the wavy interface is studied in detail.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 15–22, November–December, 1988.     

Method of calculating surface disturbances over a point source and a dipole

Formulas for calculating the deviations of the free surface of a heavy fluid in steady flow past a point source and a point dipole are derived. Examples of numerical calculations made in accordance with these formulas are presented.     

Unsteady natural convection in an enclosure with vertical wavy walls

In this paper, unsteady heat transfer and fluid flow characteristics in an enclosure are investigated. The enclosure consists of two vertical wavy and two horizontal straight walls. The top and the bottom walls are considered adiabatic. Two wavy walls are kept isothermal and their boundaries are approximated by a cosine function. Governing equations including continuity, momentum and energy were discretized using the finite-volume method and solved by SIMPLE method in curvilinear coordinate. Simulation was carried out for a range of Grashof number Gr = 10³–10⁶, Prandtl number Pr = 0.5–4.0, wave ratio A (defined by amplitude/wavelength) 0.0–0.35 and aspect ratio W (defined by average width/wavelength) 0.5–1.0. Streamlines and isothermal lines are presented to corresponding flow and thermal fields. Local and average Nusselt number distributions are presented. The obtained results are in good agreement with available numerical and experimental data.     

Scalar interfaces in digital images of turbulent flows

A scalar interface is defined as the surface separating the scalar-marked regions of a turbulent flow from the rest. The problem of determining the two-dimensional intersections of scalar interfaces is examined, taking as a specific example digital images of an axisymmetric jet visualized by laser-induced fluorescence. The usefulness of gradient and Laplacian techniques for this purpose is assessed, and it is shown that setting a proper threshold on the pixel intensity works well if the signal/noise ratio is high. Two methods of determining the proper threshold are presented, and the results are discussed. As one application of the technique, the fractal dimension of the scalar interface is calculated.     

Scalar transport in plane mixing layers

This paper describes the application of the Eulerian, single-point, single-time joint-scalar probability density function (PDF) equation for predicting the scalar transport in mixing layer with a high-speed and a low-speed stream. A finite-volume procedure is applied to obtain the velocity field with the k-ε closure being used to describe turbulent transport. The scalar field is represented through the modelled evolution equation for the scalar PDF and is solved using a Monte Carlo simulation. The PDF equation employs gradient transport modelling to represent the turbulent diffusion, and the molecular mixing term is modelled by the LMSE closure. There is no source term for chemical reaction as only an inert mixing layer is considered here. The experimental shear layer data published by Batt is used to validate the computational results despite the fact that comparisons between experiments and computational results are difficult because of the high sensitivity of the shear layer to initial conditions and free stream turbulence phenomena. However, the bimodal shape of the RMS scalar fluctuation as was measured by Batt can be reproduced with this model, whereas standard gradient diffusion calculations do not predict the dip in this profile. In this work for the first time an explanation is given for this phenomenon and the importance of a micromixing model is stressed. Also it is shown that the prediction of the PDF shape by the LMSE model is very satisfactory. Received on 27 October 1998     

Stability and nonlinear wavy regimes in downward film flows on a corrugated surface

The linear and nonlinear stability of downward viscous film flows on a corrugated surface to freesurface perturbations is analyzed theoretically. The study is performed with the use of an integral approach in ranges of parameters where the calculated results and the corresponding solutions of Navier-Stokes equations (downward wavy flow on a smooth wall and waveless flow along a corrugated surface) are in good agreement. It is demonstrated that, for moderate Reynolds numbers, there is a range of corrugation parameters (amplitude and period) where all linear perturbations of the free surface decay. For high Reynolds numbers, the waveless downward flow is unstable. Various nonlinear wavy regimes induced by varying the corrugation amplitude are determined. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 110–120, January–February, 2007.     

Entropy generation and natural convection in square cavities with wavy walls

The aim of the present work is to study the entropy generation in the natural convection process in square cavities with hot wavy walls through numerical simulations for different undulations and Rayleigh numbers, while keeping the Prandtl number constant. The results show that the hot wall geometry affects notably the heat transfer rate in the cavity. It has been found in the present numerical study that the mean Nusselt number in the case of heat transfer in a cavity with wavy walls is lower, as compared to heat transfer in a cavity without undulations. Based on the obtained dimensionless velocity and temperature values, the distributions of the local entropy generation due to heat transfer and fluid friction, the local Bejan number, and the local entropy generation are determined and plotted for different undulations and Rayleigh numbers. The study is performed for Rayleigh numbers 10³ < Ra < 10⁵, irreversibility coefficients 10^?4 < φ < 10^?2, and Prandtl numbers Pr = 0.71. The total entropy generation is found to increase with increasing undulation number.     

Stability of flows in a peristaltic transport

The stability problem related to the basic flows induced by the peristaltic waves propagating along the deformable walls is investigated numerically. The neutral stability boundary is obtained by solving the relevant Orr–Sommerfeld equation via a verified preconditioned complex-matrix solver. The critical Reynolds number becomes 577.25 when the ratio of the wave speed to the maximum speed of the basic flow (c/u_max) becomes 10.     

Heavy particle dispersion from a point source in turbulent pipe flow

Dispersion of heavy particles from a point source in high-Reynolds pipe flow was studied using large-eddy simulation, LES. A stochastic Langevin type Lagrangian model developed by Berrouk et al. was used to account for heavy particle transport by the sub-grid scale motion. In both the LES and in an experiment by Arnason, the larger particles dispersed more than the small ones. The change in diffusivity with particle size is interpreted in terms of the effect of inertia and cross-trajectory effects and qualitatively compared with the analysis of heavy particle dispersion in isotropic turbulence by Wang and Stock. Particle inertia has a much larger influence on the dispersion than the crossing-trajectories effects.     

Scalar concentration measurements in liquid-phase flows with pulsed lasers

The suitability of pulsed lasers for laser-induced-fluorescence (LIF) measurements of the local concentration of scalars in liquid-phase flows is investigated. Experiments were performed to measure the fluorescence intensity of aqueous solutions of rhodamine-6G chloride excited by a Q-switched, frequency-doubled, Nd:YAG laser. The fluorescence intensity is found to be linear with dye concentration, but not with illumination power density. The fluorescence intensity saturates at laser power densities easily exceeded by short-pulse-duration lasers. A procedure for calibrating and normalizing the raw image data is discussed which relies only on weak absorption and the linearity of fluorescence intensity with concentration. This procedure enables quantitative concentration measurement with pulsed-laser-induced fluorescence.     

Structure of modulated wavy vortical flows in the circular Couette system

Results of experimental studies are presented on the evolution of the circular Couette flow with the outer cylinder at rest. In particular, several modes of modulated wavy vortical flows have been identified in a consistent way both from flow visualization and from signals of scattered laser beam. The latter technique has been extended to measure the spatial correlation between fluctuations as a function of the azimuthal separation of the points. It has been found that the two-point angular correlation reflects the pattern of modulation, and falls off as (R_d — R)^σ when the ‘disorder Reynolds number’ R_d is approached. The value of the exponent σ has been estimated to be about 1/3 in the experiment.