Three-dimensionality of grooved channel flows at intermediate Reynolds numbers

 This paper describes the three-dimensional flow structure in grooved channels with different cavity lengths at intermediate Reynolds numbers. For steady flow, the three-dimensional effects are dominant near the side walls of the channel. However, after the onset of self-sustained oscillatory flow due to Tollmien–Schlichting waves as the primary instability, a secondary instability produces a three-dimensional flow with Taylor–Geortler-like vortical structure, at the bottom of the groove. This trend becomes more significant as the cavity length increases. Furthermore, the reason for three-dimensional flow is discussed using additional numerical analysis, and it is confirmed that the source of three-dimensional instability is the groove vortices due to the presence of side walls, rather than the channel traveling wave. Received: 7 September 1999/Accepted: 11 November 2000     

Experimental visualization of temperature fields and study of heat transfer enhancement in oscillatory flow in a grooved channel

An experimental study was conducted of incompressible, moderate Reynolds number flow of air over heated rectangular blocks in a two-dimensional, horizontal channel. Holographic interferometry combined with high-speed cinematography was used to visualize the unsteady temperature fields in self- sustained oscillatory flow. Experiments were conducted in the laminar, transitional and turbulent flow regimes for Reynolds numbers in the range from Re = 520 to Re = 6600. Interferometric measurements were obtained in the thermally and fluiddynamically periodically fully developed flow region on the ninth heated block. Flow oscillations were first observed between Re = 1054 and Re = 1318. The period of oscillations, wavelength and propagation speed of the Tollmien–Schlichting waves in the main channel were measured at two characteristic flow velocities, Re = 1580 and Re = 2370. For these Reynolds numbers it was observed that two to three waves span one geometric periodicity length. At Re = 1580 the dominant oscillation frequency was found to be around 26 Hz and at Re = 2370 the frequency distribution formed a band around 125 Hz. Results regarding heat transfer and pressure drop are presented as a function of the Reynolds number, in terms of the block-average Nusselt number and the local Nusselt number as well as the friction factor. Measurements of the local Nusselt number together with visual observations indicate that the lateral mixing caused by flow instabilities is most pronounced along the upstream vertical wall of the heated block in the groove region, and it is accompanied by high heat transfer coefficients. At Reynolds numbers beyond the onset of oscillations the heat transfer in the grooved channel exceeds the performance of the reference geometry, the asymmetrically heated parallel plate channel. Received on 26 April 2000     

Linear stability analysis in compressible, flat-plate boundary-layers

The stability problem of two-dimensional compressible flat-plate boundary layers is handled using the linear stability theory. The stability equations obtained from three-dimensional compressible Navier–Stokes equations are solved simultaneously with two-dimensional mean flow equations, using an efficient shoot-search technique for adiabatic wall condition. In the analysis, a wide range of Mach numbers extending well into the hypersonic range are considered for the mean flow, whereas both two- and three-dimensional disturbances are taken into account for the perturbation flow. All fluid properties, including the Prandtl number, are taken as temperature-dependent. The results of the analysis ascertain the presence of the second mode of instability (Mack mode), in addition to the first mode related to the Tollmien–Schlichting mode present in incompressible flows. The effect of reference temperature on stability characteristics is also studied. The results of the analysis reveal that the stability characteristics remain almost unchanged for the most unstable wave direction for Mach numbers above 4.0. The obtained results are compared with existing numerical and experimental data in the literature, yielding encouraging agreement both qualitatively and quantitatively.      

Linear stability of two-dimensional steady flow in wavy-walled channels

Linear stability of fully developed two-dimensional periodic steady flows in sinusoidal wavy-walled channels is investigated numerically. Two types of channels are considered: the geometry of wavy walls is identical and the location of the crest of the lower and upper walls coincides (symmetric channel) or the crest of the lower wall corresponds to the furrow of the upper wall (sinuous channel). It is found that the critical Reynolds number is substantially lower than that for plane channel flow and that when the non-dimensionalized wall variation amplitude is smaller than a critical value (about 0.26 for symmetric channel, 0.28 for sinuous channel), critical modes are three-dimensional stationary and for larger , two-dimensional oscillatory instabilities set in. Critical Reynolds numbers of sinuous channel flows are smaller for three-dimensional disturbances and larger for two-dimensional disturbances than those of symmetric channel flows. The disturbance velocity distribution obtained by the linear stability analysis suggests that the three-dimensional stationary instability is mainly caused by local concavity of basic flows near the reattachment point, while the critical two-dimensional mode resembles closely the Tollmien–Schlichting wave for plane Poiseuille flow.     

Numerical Simulation of Riblet Controlled Spatial Transition in a Zero-Pressure-Gradient Boundary Layer

To analyze the fundamental physical mechanism which determines the damping effect of a riblet surface on three-dimensional transition several numerical simulations of spatial transition in a flat plate zero-pressure-gradient boundary layer above a riblet wall are performed in this study. Two types of forced transition scenarios are investigated. The first type of transition is defined by K-type transition induced by a dominant two-dimensional Tollmien–Schlichting (TS) wave and a weak spanwise disturbance. The second type of transition is purely excited by two oblique waves. By a qualitative analysis of the occurring maximum wall-normal and spanwise velocity components and the Fourier modes of the disturbances the two-dimensional TS waves are found to be amplified by riblets, whereas three-dimensional structures, i.e., Λ-, hairpin, and streamwisely aligned vortices, are damped. At oblique transition the breakdown to turbulence is delayed by the riblets compared to transition on a clean surface. The investigation of the near wall flow structure reveals secondary flows induced by the riblets and reduced wall normal ejections as well as a reduced downwash.     

Transition on a slender cone in hypervelocity flow

 An experimental exploratory study of laminar, transitional and turbulent boundary layers on a slender, sharp cone at zero incidence in hypervelocity flow was conducted in the shock tunnel T5, whose high operating pressure permits achieving the necessary high Reynolds numbers at high total enthalpy. Surface heat transfer rate and resonantly enhanced flow visualization were used to determine transition location and flow structure. Results suggest that the transition mechanism at the conditions tested is via the Tollmien–Schlichting instability, but since no detailed information about the noise spectrum in the facility is available, a definite conclusion is not possible. A strong, gas-dependent increase of the transition Reynolds number evaluated at the reference temperature on specific total enthalpy was observed. This effect increases monotonically with decreasing the dissociation energy of the gas. Received: 1 February 1996/Accepted: 28 June 1996     

Active cancellation of artificially introduced Tollmien–Schlichting waves using plasma actuators

In the present work artificially excited Tollmien–Schlichting (TS) waves were cancelled using plasma actuators operated in pulsed mode. In order to achieve this a vibrating surface driven by an electromagnetic turbulator was flush mounted in a flat plate to excite the TS waves. These were amplified by an adverse pressure gradient induced by an insert on the upper wall of the test section. A control plasma actuator positioned downstream of the excitation actuator attenuates the waves by imparting an unsteady force into the boundary layer to counteract the oscillation. As a result the amplitude of the velocity fluctuations at the excitation frequency is reduced significantly depending on the distance from the wall. A parameter study was performed to identify the influence of several operation parameters of the control actuator.     

Application of the Discrete Fourier Transform to Study the Dynamics of Wave Packets

A method for solving equations that describe the dynamics of wave packets of the Tollmien–Schlichting waves in the boundary layer is proposed. The method of splitting the initial problem into the linear and nonlinear parts at each time step is used. The linear part is resolved by using an equation for spectral components of the wave packet with a subsequent Fourier transform from the space of wavenumbers to the physical space. A system of ordinary differential equations is solved in the physical space. The Fourier transform is performed by means of the library procedure of the fast Fourier transform. As examples, the problems solved were the linear dynamics of the wave packet concentrated in the vicinity of the instability region (i.e., a set of wave vectors in the space of wavenumbers for which the imaginary part of the eigenfrequency of the Tollmien–Schlichting waves is positive) and the nonlinear dynamics of the wave packet overlapping the instability region.     

On Tollmien–Schlichting-like waves in streaky boundary layers

The linear stability of the boundary layer developing on a flat plate in the presence of finite-amplitude, steady and spanwise periodic streamwise streaks is investigated. The streak amplitudes considered here are below the threshold for onset of the inviscid inflectional instability of sinuous perturbations. It is found that, as the amplitude of the streaks is increased, the most unstable viscous waves evolve from two-dimensional Tollmien–Schlichting waves into three-dimensional varicose fundamental modes which compare well with early experimental findings. The analysis of the growth rates of these modes confirms the stabilising effect of the streaks on the viscous instability and that this stabilising effect increases with the streak amplitude. Varicose subharmonic modes are also found to be unstable but they have growth rates which typically are an order of magnitude lower than those of fundamental modes. The perturbation kinetic energy production associated with the spanwise shear of the streaky flow is found to play an essential role in the observed stabilisation. The possible relevance of the streak stabilising role for applications in boundary layer transition delay is discussed.     

Three-Dimensional Instability of Planar Flows

We study the stability of two-dimensional solutions of the three-dimensional Navier–Stokes equations, in the limit of small viscosity. We are interested in steady flows with locally closed streamlines. We consider the so-called elliptic and centrifugal instabilities, which correspond to the continuous spectrum of the underlying linearized Euler operator. Through the justification of highly oscillating Wentzel–Kramers–Brillouin expansions, we prove the nonlinear instability of such flows. The main difficulty is the control of nonoscillating and nonlocal perturbations issued from quadratic interactions.     

Localized Coherent Structures in the Boundary Layer

A Blasius laminar boundary layer and a steady turbulent boundary layer on a flat plate in an incompressible fluid are considered. The spectral characteristics of the Tollmien—Schlichting (TS) and Squire waves are numerically determined in a wide range of Reynolds numbers. Based on the spectral characteristics, relations determining the three–wave resonance of TS waves are studied. It is shown that the three–wave resonance is responsible for the appearance of a continuous low–frequency spectrum in the laminar region of the boundary layer. The spectral characteristics allow one to obtain quantities that enter the equations of dynamics of localized perturbations. By analogy with the laminar boundary layer, the three–wave resonance of TS waves in a turbulent boundary layer is considered.     

On two-dimensional linear waves in Blasius boundary layer over viscoelastic layers

This work concerns the direct numerical simulation of small-amplitude two-dimensional ribbon-excited waves in Blasius boundary layer over viscoelastic compliant layers of finite length. A vorticity-streamfunction formulation is used, which assures divergence-free solutions for the evolving flow fields. Waves in the compliant panels are governed by the viscoelastic Navier's equations. The study shows that Tollmien–Schlichting (TS) waves and compliance-induced flow instability (CIFI) waves that are predicted by linear stability theory frequently coexist on viscoelastic layers of finite length. In general, the behaviour of the waves is consistent with the predictions of linear stability theory. The edges of the compliant panels, where abrupt changes in wall property occur, are an important source of waves when they are subjected to periodic excitation by the flow. The numerical results indicate that the non-parallel effect of boundary-layer growth is destabilizing on the TS instability. It is further demonstrated that viscoelastic layers with suitable properties are able to reduce the amplification of the TS waves, and that high levels of material damping are effective in controlling the propagating CIFI.     

Direct Spatial Resonance in the Compressible Boundary Layer on a Rotating-disk

This paper focuses on a resonance mechanism that can lead to significant disturbance amplification at conditions which are sub-critical to nonlinear instabilities. Particularly, direct spatial resonance instability is investigated, which is present in the basic three-dimensional viscous compressible boundary-layer flow due to a rotating-disk. Within this purpose, the linearized system of stability equations is treated numerically making use of a spectral Chebyshev collocation method. The analysis provides critical resonant Reynolds numbers above which growth occurs. Amplitudes of the response of the degeneracies decaying rapidly due to their high damping rates are shown to exist for small enough Reynolds numbers while the flow is still in the laminar state. If the flow is restricted to the incompressible case, the results of Turkyilmazoglu and Gajjar (in Sadhana Acad P Engs 25:601–617, 2000) are completely reproduced. The influences of compressibility are then explored by means of varying the Mach and Prandtl numbers in the cases of heating/cooling the wall as well as the isothermal wall. In general, compressibility effects are found strongly in favor of stabilizing as the Mach number increases, while a strong destabilization is observed by lowering the critical values of Reynolds numbers in the cases of wall heating and insulation. The modal interaction and coalescence of the eigenmodes calculated here create local algebraic growth by rapid development of relatively large amplitudes which might then provide the onset of nonlinear effects followed by transition.     

Pairwise nonlinear interactions of tollmien-schlichting waves in flows of the boundary-layer type

For a flow of the boundary-layer type, the method of a phase plane is used to investigate the joint evolution of a pair of Tollmien–Schlichting waves as a function of the R numbers and the frequencies.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 33–37, March–April, 1977.     

Experiments on the linear instability of flow in a wavy channel

The effects of wall corrugation on the stability of wall-bounded shear flows have been examined experimentally in plane channel flows. One of the channel walls has been modified by introduction of the wavy wall model with the amplitude of 4% of the channel half height and the wave number of 1.02. The experiment is focused on the two-dimensional travelling wave instability and the results are compared with the theory [J.M. Floryan, Two-dimensional instability of flow in a rough channel, Phys. Fluids 17 (2005) 044101 (also: Rept. ESFD-1/2003, Dept. of Mechanical and Materials Engineering, The University of Western Ontario, London, Ontario, Canada, 2003)]. It is shown that the flow is destabilized by the wall corrugation at subcritical Reynolds numbers below 5772, as predicted by the theory. For the present corrugation geometry, the critical Reynolds number is decreased down to about 4000. The spatial growth rates, the disturbance wave numbers and the distribution of disturbance amplitude measured over such wavy wall also agree well with the theoretical results.     

Characterization of noise amplifiers with global singular modes: the case of the leading-edge flat-plate boundary layer

This article deals with the linear dynamics of a transitional boundary layer subject to two-dimensional Tollmien–Schlichting instabilities. We consider a flat plate including the leading edge, characterized by a Reynolds number based on the length of the plate equal to Re = 6 × 10⁵, inducing a displacement thickness-based Reynolds number of 1,332 at the end of the plate. The global linearized Navier–Stokes equations only display stable eigenvalues, and the associated eigen-vectors are known to poorly represent the dynamics of such open flows. Therefore, we resort to an input–output approach by considering the singular value decomposition of the global resolvent. We then obtain a series of singular values, an associated orthonormal basis representing the forcing (the so-called optimal forcings) as well as an orthonormal basis representing the response (the so-called optimal responses). The objective of this paper is to analyze these spatial structures and to closely relate their spatial downstream evolution to the Orr and Tollmien–Schlichting mechanisms. Analysis of the spatio-frequential diagrams shows that the optimal forcings and responses are respectively localized, for all frequencies, near the upstream neutral point (branch I) and the downstream neutral point (branch II). It is also shown that the spatial growth of the dominant optimal response favorably compares with the e ^N method in regions where the dominant optimal forcing is small. Moreover, thanks to an energetic input–output approach, it is shown that this spatial growth is solely due to intrinsic amplifying mechanisms related to the Orr and Tollmien–Schlichting mechanisms, while the spatial growth due to the externally supplied power by the dominant optimal forcing is negligible even in regions where the dominant optimal forcing is strong. The energetic input–output approach also yields a general method to assess the strength of the instability in amplifier flows. It is based on a ratio comparing two quantities of same physical dimension, the mean-fluctuating kinetic energy flux of the dominant optimal response across some boundary and the supplied mean external power by the dominant optimal forcing. For the present boundary-layer flow, we have computed this amplification parameter for each frequency and discussed the results with respect to the Orr and Tollmien–Schlichting mechanisms.     

Small Reynolds number instabilities in two-layer Couette flow

Instabilities in two-layer Couette flow are investigated from a small Reynolds number expansion of the eigenvalue problem governing linear stability. The wave velocity and growth rate are given explicitly, and previous results for long waves and short waves are retrieved as special cases. In addition to the inertial instability due to viscous stratification, the flow may be subject to the Rayleigh–Taylor instability. As a result of the competition of these two instabilities, inertia may completely stabilise a gravity-unstable flow above some finite critical Froude number, or conversely, for a gravity-stable flow, inertia may give rise to finite wavenumber instability above some finite critical Weber number. General conditions for these phenomena are given, as well as exact or approximate values of the critical numbers. The validity domain of the many asymptotic expansions is then investigated from comparison with the numerical solution. It appears that the small-Re expansion gives good results beyond Re = 1, with an error less that 1%. For Reynolds numbers of a few hundred, we show from the eigenfunctions and the energy equation that the nature of the instability changes: instability still arises from the interfacial mode (there is no mode crossing), but this mode takes all the features of a shear mode. The other modes correspond to the stable eigenmodes of the single-layer Couette flow, which are recovered when one fluid is rigidified by increasing its viscosity or surface tension.     

Laminar–turbulent transition in oscillating boundary layer: experimental and numerical analysis using continuous wavelet transform

The laminar–turbulent transition of a forced oscillating boundary layer with a varying pressure gradient is experimentally and numerically investigated for two Strouhal numbers. Time-dependent characterization of the natural instability modes is carried out using continuous wavelet analysis of velocity signals. The periodic evolution of the total growth rates of the most unstable disturbances are measured and compared to the results of the linear stability theory. The “Tollmien–Schlichting” (TS) and “convective” transition modes are identified. It is shown that they correspond to the extrema of opposite signs of the skewness factor of unstable wavelet transform of the most unstable frequencies.     

Unsteady PTV velocity field past an airfoil solved with DNS: Part 1. Algorithm of hybrid simulation and hybrid velocity field at <Emphasis Type="Italic">Re</Emphasis> ≈ 10<Superscript>3</Superscript>

We develop a hybrid unsteady-flow simulation technique combining direct numerical simulation (DNS) and particle tracking velocimetry (PTV) and demonstrate its capabilities by investigating flows past an airfoil. We rectify instantaneous PTV velocity fields in a least-squares sense so that they satisfy the equation of continuity, and feed them to the DNS by equating the computational time step with the frame rate of the time-resolved PTV system. As a result, we can reconstruct unsteady velocity fields that satisfy the governing equations based on experimental data, with the resolution comparable to numerical simulation. In addition, unsteady pressure distribution can be solved simultaneously. In this study, particle velocities are acquired on a laser-light sheet in a water tunnel, and unsteady flow fields are reconstructed with the hybrid algorithm solving the incompressible Navier–Stokes equations in two dimensions. By performing the hybrid simulation, we investigate nominally two-dimensional flows past the NACA0012 airfoil at low Reynolds numbers. In part 1, we introduce the algorithm of the proposed technique and discuss the characteristics of hybrid velocity fields. In particular, we focus on a vortex shedding phenomenon under a deep stall condition (α = 15°) at Reynolds numbers of Re = 1000 and 1300, and compare the hybrid velocity fields with those computed with two-dimensional DNS. In part 2, the extension to higher Reynolds numbers is considered. The accuracy of the hybrid simulation is evaluated by comparing with independent experimental results at various angles of attack and Reynolds numbers up to Re = 10⁴. The capabilities of the hybrid simulation are also compared with two-dimensional unsteady Reynolds-Averaged Navier–Stokes (URANS) solutions in part 2. In the first part of these twin papers, we demonstrate that the hybrid velocity field approaches the PTV velocity field over time. We find that intensive alternate vortex shedding past the airfoil, which is predicted by the two-dimensional DNS, is substantially suppressed in the hybrid simulation and the resultant flow field is similar to the PTV velocity field, which is projection of the three-dimensional velocity field on the streamwise plane. We attempt to identify the motion that originates three-dimensional flow patterns by highlighting the deviation of the PTV velocity field from the two-dimensional governing equations at each snapshot. The results indicate that the intensive spots of the deviation appear in the regions in which three-dimensional instabilities are induced in the shear layer separated from the pressure side.