Scalar transport from a point source in flows over wavy walls

Simultaneous measurements of the velocity and concentration field in fully developed turbulent flows over a wavy wall are described. The concentration field originates from a low-momentum plume of a passive tracer. PLIF and digital particle image velocimetry are used to make spatially resolved measurements of the structure of the scalar distribution and the velocity. The measurements are performed at three different Reynolds numbers of Re _b = 5,600, Re _b = 11,200 and Re _b = 22,400, respectively, based on the bulk velocity u _b and the total channel height 2h. The velocity field and the scalar field are investigated in a water channel with an aspect ratio of 12:1, where the bottom wall of the test section consists of a train of sinusoidal waves. The wavy wall is characterized by the amplitude to wavelength ratio α = 0.05 and the ratio β between the wave amplitude and the half channel height where β = 0.1. The scalar is released from a point source at the wave crest. For the concentration measurements, Rhodamine B is used as tracer dye. At low to moderate Reynolds number, the flow field is characterized through a recirculation zone which develops after the wave crest. The recirculation zone induces high intensities of the fluctuations of the streamwise velocity and wall-normal velocity. Furthermore, large-scale structures are apparent in the flow field. In previous investigations it has been shown that these large-scale structures meander laterally in flows over wavy bottom walls. The investigations show a strong effect of the wavy bottom wall on the scalar mixing. In the vicinity of the source, the scalar is transported by packets of fluid with a high scalar concentration. As they move downstream, these packets disintegrate into filament-like structures which are subject to strong gradients between the filaments and the surrounding fluid. The lateral scale of the turbulent plume is smaller than the lateral scale of the large-scale structures in the flow field and the plume dispersion is dominated by the structures in the flow field. Due to the lateral meandering of the large-scale structures of the flow field, also the scalar plume meanders laterally. Compared to turbulent plumes in plane channel flows, the wavy bottom wall enhances the mixing effect of the turbulent flow and the spreading rate of the scalar plume is increased.     

Hydrodynamics study in a plane ultrafiltration module using an electrochemical method and particle image velocimetry visualization

 The wall shear stress is determined at the surface of a plane sheet of Plexiglas, taking the place of a membrane, using an electrochemical method. Several microelectrodes are mounted flush to this plane plate, and maps of shear stress are determined for two inlet and outlet configurations and three channel heights. The heterogeneity of the wall shear stress is observed for both configurations. Furthermore, the study of the turbulence features of the flow shows a decreasing fluctuating rate of velocity gradient when the channel height is decreased. The wall velocity gradients and turbulent intensity rates analysis are confirmed by a flow visualization using the particle image velocimetry method. Received: 25 September 2000 / Accepted: 23 April 2001     

Temporal Evolution of Vortical Structures in the Wall Region of Turbulent Channel Flow

Hairpin-like vortical structures that form in the wall region of turbulent channel flow are investigated. The analysis is performed by following a procedure in which the Navier-Stokes equations are first integrated by means of a computational code based on a mixed spectral-finite difference technique in the case of the flow in a plane channel. A DNS turbulent-flow database, representing the turbulent statistically steady state of the velocity field through 10 viscous time units, is computed and the vortex-detection method of the imaginary part of the complex eigenvalue pair of the velocity-gradient tensor is applied to the velocity field. As a result, hairpin-like vortical structures are educed. Flow visualizations are provided of the processes of evolution that characterize hairpin vortices in the wall region of turbulent channel flow. The relationship is investigated between vortex dynamics and 2nd- and 4th- quadrant events, showing that ejections and sweeps play a fundamental role in the way the morphological evolution of a hairpin vortex develops with time.     

Study of laminar–turbulent flow transition under pulsatile conditions in a constricted channel

In this paper, direct numerical simulation is performed to investigate a pulsatile flow in a constricted channel to gain physical insights into laminar–turbulent–laminar flow transitions. An in-house computer code is used to conduct numerical simulations based on available high-performance shared memory parallel computing facilities. The Womersley number tested is fixed to 10.5 and the Reynolds number varies from 500 to 2000. The influences of the degree of stenosis and pulsatile conditions on flow transitions and structures are investigated. In the region upstream of the stenosis, the flow pattern is primarily laminar. Immediately after the stenosis, the flow recirculates under an adverse streamwise pressure gradient, and the flow pattern transitions from laminar to turbulent. In the region far downstream of the stenosis, the flow becomes re-laminarised. The physical characteristics of the flow field have been thoroughly analysed in terms of the mean streamwise velocity, turbulence kinetic energy, viscous wall shear stresses, wall pressure and turbulence kinetic energy spectra.     

Reconstruction of a quasi-instantaneous image of coherent structures from hotwire signals obtained by a multi-point simultaneous measurement system

The turbulent velocity components (u, v) at 11 points in a reciprocating oscillatory turbulent flow have been measured simultaneously by a set of eleven X-type hotwire probes located in a plane perpendicular to the mean flow. Using a conditional sampling technique and a new method of data analysis for the inverse estimation of flow fields called the “virtual plate/load and MASCON model”, a quasi-instantaneous three-dimensional image of coherent structures of turbulence was first reconstructed directly from the experimental velocity data. The quasi-instantaneous image was expressed in terms of the velocity components u, v, w and the vorticity components ω_x ω_y, ω_z and we found that the large-scale coherent structure was composed of a pair of counter-rotating fluid motions with asymmetry which was quite different from that of the ensemble-averaged one. Flow patterns induced by the large-scale structure have been clarified by perspective representations visualized by computer simulations that produce timelines and streaklines of fluid particle traces. Results showed that the new experimental method was applicable for investigating the three-dimensional feature of coherent structures including asymmetry.     

Scanning liquid-crystal thermometry and stereo velocimetry for simultaneous three-dimensional measurement of temperature and velocity field in a turbulent Rayleigh-Bérnard convection

In this paper, we report on an experimental technique for the simultaneous measurement of temperature and three components of velocity in a three-dimensional thermal flow using scanning liquid-crystal thermometry and stereo velocimetry. The temperature is measured by the color image analysis of the liquid-crystal particles suspended in a fluid, while the three velocity components are measured by stereo particle image velocimetry (stereo PIV) with the aid of tracer particles. The measurement is carried out by scanning the light-sheet plane while capturing the sequential color images of the liquid crystals and tracer particles. This measurement allows the reconstruction of the three-dimensional distribution of temperature and full velocity field simultaneously. The present experimental technique is applied to the horizontal fluid layer of a turbulent Rayleigh-Bérnard convection and the three-dimensional structures of thermal plumes are evaluated. The experimental results indicate that the structures of plumes are often correlated with the vertical velocity of the fluid, but they behave randomly in space, influenced by the large-scale turbulence evident in the middle of the fluid layer.     

Quantitative visualization of the near-wall structures in a turbulent pipe flow by image correlation velocimetry

Laser-induced fluorescent dye visualization and image correlation velocimetry were employed to delineate near-wall turbulent structures in a pipe flow. The sweeping and ejection events near the wall and the downstream evolution of a large-scale eddy structure rotating in a counter-clockwise direction were clearly reflected in the instantaneous fluctuating velocity fields. This eddy structure was found to form mostly in the logarithmic region and to dominate the flow structures there, while the ejection and sweeping events in the log layer were greatly influenced by the existence of the large-scale eddy structure. Received: 29 January 2001 / Accepted: 22 October 2001     

Measurement of velocity and velocity derivatives based on pattern tracking in 3D LIF images

Pattern tracking in consecutive 3D LIF images based on least squares matching (LSM) of grey levels has been developed recently for velocity and velocity gradient measurements. The shortcomings of this method are clearly shown. The present article presents an improvement on this method by introducing a local multi-patch (LMP) technique through the LSM approach. The method is validated using the flow field of a turbulent channel flow obtained by direct numerical simulation (DNS) and a synthetic image with grey-level patterns. The results show that LMP matching allows the determination of the velocity and the velocity gradient fields with high accuracy including the second derivatives. Measurements of a round non-buoyant jet are presented which demonstrate the good performance of the method when applied under laboratory conditions. This method can also be applied on two-dimensional images provided that the flow is strictly two-dimensional. Received: 7 July 1999/Accepted: 13 November 1999     

An investigation of thermally stratified turbulent channel flow with temperature oscillation on the bottom wall by large eddy simulation

Thermally stratified shear turbulent channel flow with temperature oscillation on the bottom wall of the channel is calculated to investigate the behavior of turbulent flow and heat transfer by use of large eddy simulation (LES) approach coupled with dynamic subgrid-scale (SGS) models. The objective of this study is to deal with the effect of the temperature oscillation on turbulent behavior of thermally stratified turbulent channel flow and to examine the effectiveness of the LES technique for predicting statistically unsteady turbulent flow driven by time-varying buoyancy force. To validate the present calculation, thermally stratified shear turbulent channel flow is computed and compared with available data obtained by direct numerical simulation (DNS), which confirm that the present approach can be used to predict thermally stratified turbulent channel flow satisfactorily. Further, to illustrate the effect of the temperature oscillation with different Richardson numbers and periods of the oscillation on turbulence characteristics, the phase-averaged mean value and fluctuation of the resolved velocities and temperature, and instantaneous velocity fluctuation structures are analyzed.     

Passive Scalar in a Turbulent Channel Flow with Wall Velocity Disturbances

Direct Numerical Simulations (DNS) of a passive scalar in a turbulent channel flow with a normal velocity disturbance on the lower wall are presented for high and low Reynolds numbers. The aim is to reproduce the complex physics of turbulent rough flows without dealing with the geometric complexity. In addition, isothermal walls that cannot be easily assigned in an experiment, are considered. The paper explains the increase of heat transfer through the changes of the velocity and thermal structures. As in real rough flows, the transpiration produces an isotropization of the turbulence near the wall.     

Tracking of coherent thermal structures on a heated wall.

The temporal evolution of a thermal pattern observed on a heated wall by infrared camera is correlated with the propagation velocity of the thermal perturbations calculated by DNS. In the experiment the propagation velocity was measured by using PIV-based analysis of infrared images of the thermal pattern on the wall. To verify the experimental technique of image analysis, a sequence of synthetic images, simulating thermal patterns on the wall, was generated from the DNS solution, and the convective velocity was evaluated. It was found that the convective velocity of thermal structures obtained by PIV-based analysis of the experimental and synthetic images was in relatively good agreement with that calculated from the DNS solution. The present study confirmed that for a high Prandtl number fluid (water) the propagation velocity of the thermal perturbations is only about half of the convective velocity of the velocity perturbations. It was also found that the convection velocity observed for hot spots is distinctly lower than that for the cold spots.     

Transient forced convection heat transfer in a channel

A numerical analysis is made of incompressible transient turbulent flow heat transfer between two parallel plates when there is a step jump in space along the channel in wall heat flux or wall temperature. The variation of the fluid velocity and effective diffusivity over the channel cross section are accounted for. The fluid is assumed to have a fully-developed turbulent velocity profile throughout the length of the channel. The thermal responses of the system are obtained by solving energy equation for air by a digital computer. The results are presented in graphical forms. The stability of the finite difference solution is studied and condition for the stability of the difference solution is derived. A method is given to obtain velocity distributions from the distribution of turbulent eddy diffusivity of momentum. Variations of Nusselt numbers are obtained as a function of time and space. Steady-state values are also given and compared with the published results.     

The nature of near-wall convection velocity in turbulent channel flow

A novel notion of turbulent structure the local cascade structure-is introduced to study the convection phenomenon in a turbulent channel flow. A space-time cross-correlation method is used to calculate the convection velocity. It is found that there are two characteristic convection speeds near the wall, one associated with small-scale streaks of a lower speed and another with streamwise vortices and hairpin vortices of a higher speed. The new concept of turbulent structure is powerful to illustrate the dominant role of coherent structures in the near-wall convection, and to reveal also the nature of the convection-the propagation of patterns of velocity fluctuations-which is scale-dependent.     

Large-eddy simulation of turbulent channel flow at transcritical states

We present well-resolved large-eddy simulations (LES) of a channel flow solving the fully compressible Navier–Stokes equations in conservative form. An adaptive look-up table method is used for thermodynamic and transport properties. A physically consistent subgrid-scale turbulence model is incorporated, that is based on the Adaptive Local Deconvolution Method (ALDM) for implicit LES. The wall temperatures are set to enclose the pseudo-boiling temperature at a supercritical pressure, leading to strong property variations within the channel geometry. The hot wall at the top and the cold wall at the bottom produce asymmetric mean velocity and temperature profiles which result in different momentum and thermal boundary layer thicknesses. Different turbulent Prandtl number formulations and their components are discussed in context of strong property variations.     

Turbulent Channel Flow with an Artificial Two-Dimensional Wall Layer

Turbulent flow of an incompressible fluid in a plane channel with parallel walls is considered. The three-dimensional time-dependent Navier-Stokes equations are solved numerically using the spectral finite-difference method. An artificial force which completely suppresses lateral oscillations of the velocity is introduced in the near-wall zone (10 % of the channel half-width in the neighborhood of each wall). Thus, the three-dimensional flow zone, in which turbulent oscillations can develop, is separated from the wall by a fluid layer. It is found that the elimination of three-dimensionality in the neighborhood of the walls leads to a significant reduction in the drag. However, complete laminarization does not occur. The flow in the stream core remains turbulent and can be interpreted as a turbulent flow in a channel with walls located on the boundary of the two-dimensional layer and traveling at the local mean-flow velocity. The oscillations developing inside the two-dimensional layer, which have significant amplitude, distort the flow only in the adjacent zone. Beyond this zone the distributions of the mean characteristics and the structure of instantaneous fields completely correspond to ordinary turbulent flow in a channel with rigid walls. The results obtained confirm the hypothesis of the unimportance of the no-slip boundary conditions for the fluctuating velocity component in the mechanism of onset and self-maintenance of turbulence in wall flows.     

PIV study on a shock-induced separation in a transonic flow

A transonic interaction between a steady shock wave and a turbulent boundary layer in a Mach 1.4 channel flow is experimentally investigated by means of particle image velocimetry (PIV). In the test section, the lower wall is equipped with a contour profile shaped as a bump allowing flow separation. The transonic interaction, characterized by the existence in the outer flow of a lambda shock pattern, causes the separation of the boundary layer, and a low-speed recirculating bubble is observed downstream of the shock foot. Two-component PIV velocity measurements have been performed using an iterative gradient-based cross-correlation algorithm, providing high-speed and flexible calculations, instead of the classic multi-pass processing with FFT-based cross-correlation. The experiments are performed discussing all the hypotheses linked to the experimental set-up and the technique of investigation such as the two-dimensionality assumption of the flow, the particle response assessment, the seeding system, and the PIV correlation uncertainty. Mean velocity fields are presented for the whole interaction with particular attention for the recirculating bubble downstream of the detachment, especially in the mixing layer zone where the effects of the shear stress are most relevant. Turbulence is discussed in details, the results are compared to previous study, and new results are given for the turbulent production term and the return to isotropy mechanism. Finally, using different camera lens, a zoom in the vicinity of the wall presents mean and turbulent velocity fields for the incoming boundary layer.     

An experimental investigation of the flow structure over a corrugated waveform in a transpired air collector

An experimental investigation of the flow dynamics in a channel with a corrugated surface is presented. Particle image velocimetry was used to obtain two-dimensional velocity fields at three different locations along the channel length, over a range of Reynolds numbers. The results show a significant impact of the corrugation waveform on the mean and turbulent flow structure inside the channel. Strong bursting flow originating from the trough, sweeping flow from the bulk region and the vortex shedding off the crest were observed. Their interactions created a complex three-dimensional flow structure extended over almost the entire channel. The mean velocity profiles indicate a strong diffusion of shear. The profiles of various turbulent properties show the enhancement of turbulence in the vicinity of the waveform. It was found that the turbulence in the channel was almost entirely produced in this region above the corrugation trough. Significant momentum transfer from the corrugation wall by the turbulent velocity field was also observed. The mean and turbulent flow behaviour was found to be periodic with respect to the waveform over most of the channel length. The results show the presence of strong turbulence even at the Reynolds number that falls within the conventional laminar range.     

Experimental study of a three-dimensional large-scale structure in a reciprocating oscillatory flow

An experimental investigation has been conducted on the large-scale structure of the reciprocating oscillatory turbulent flow in a rectangular duct. Using hot-wire anemometers and conditional sampling technique, we found that the Reynolds stress generated explosively in the decelerating phase, a characteristic feature of the flow, is caused by the three-dimensional large-scale structure. These structures are revealed to be a pair of counter-rotating fluid motions which change their form gradually as the phase-averaged velocity distribution changes. Phase differences of the turbulence intensity distribution in an oscillatory flow are found to take place in accordance with the generation-growth-decay process of the large-scale structure.     

Influence of wavy surfaces on coherent structures in a turbulent flow

We describe how outer flow turbulence phenomena depend on the interaction with the wall. We investigate coherent structures in turbulent flows over different wavy surfaces and specify the influence of the different surface geometries on the coherent structures. The most important contribution to the turbulent momentum transport is attributed to these structures, therefore this flow configuration is of large engineering interest. In order to achieve a homogeneous and inhomogeneous reference flow situation two different types of surface geometries are considered: (1) three sinusoidal bottom wall profiles with different amplitude-to-wavelength ratios of α = 2a/Λ = 0.2 (Λ = 30 mm), α = 0.2 (Λ = 15 mm), and α = 0.1 (Λ = 30 mm); and (2) a profile consisting of two superimposed sinusoidal waves with α = 0.1 (Λ = 30 mm). Measurements are carried out in a wide water channel facility (aspect ratio 12:1). Digital particle image velocimetry (PIV) is performed to examine the spatial variation of the streamwise, spanwise and wall-normal velocity components in three measurement planes. Measurements are performed at a Reynolds number of 11,200, defined with the half channel height h and the bulk velocity U _B. We apply the method of snapshots and perform a proper orthogonal decomposition (POD) of the streamwise, spanwise, and wall-normal velocity components to extract the most dominant flow structures. The structure of the most dominant eigenmode is related to counter-rotating, streamwise-oriented vortices. A qualitative comparison of the eigenfunctions for different sinusoidal wall profiles shows similar structures and comparable characteristic spanwise scales Λ _z = 1.5 H in the spanwise direction for each mode. The scale is observed to be slightly smaller for α = 0.2 (Λ = 15 mm) and slightly larger for α = 0.2 (Λ = 30 mm). This scaling for the flow over the basic wave geometries indicates that the size of the largest structures is neither directly linked to the solid wave amplitude, nor to the wavelength. The characteristic spanwise scale of the dominant eigenmode for the developed flow over the surface consisting of two superimposed waves reduces to 0.85 H. However, a scale in the order of 1.3 H is identified for the second mode. The eigenvalue spectra for the superimposed waves is much broader, more modes contribute to the energy-containing range. The turbulent flow with increased complexity of the bottom surface is characterized by an increased number of dominant large-scale structures with different spanwise scales.     

Numerical prediction of turbulent boundary layer noise from a sharp-edged flat plate

An efficient hybrid uncorrelated wall plane waves–boundary element method (UWPW-BEM) technique is proposed to predict the flow-induced noise from a structure in low Mach number turbulent flow. Reynolds-averaged Navier-Stokes equations are used to estimate the turbulent boundary layer parameters such as convective velocity, boundary layer thickness, and wall shear stress over the surface of the structure. The spectrum of the wall pressure fluctuations is evaluated from the turbulent boundary layer parameters and by using semi-empirical models from literature. The wall pressure field underneath the turbulent boundary layer is synthesized by realizations of uncorrelated wall plane waves (UWPW). An acoustic BEM solver is then employed to compute the acoustic pressure scattered by the structure from the synthesized wall pressure field. Finally, the acoustic response of the structure in turbulent flow is obtained as an ensemble average of the acoustic pressures due to all realizations of uncorrelated plane waves. To demonstrate the hybrid UWPW-BEM approach, the self-noise generated by a flat plate in turbulent flow with Reynolds number based on chord Re_c = 4.9 × 10⁵ is predicted. The results are compared with those obtained from a large eddy simulation (LES)-BEM technique as well as with experimental data from literature.