Three-dimensional receptivity of boundary layers

The paper presents a review of results of some recent (mainly experimental) studies devoted to a quantitative investigation of the problem of receptivity of the 2D and 3D boundary layers with respect to various 3D (in general) external perturbations. The paper concentrates on the mechanisms of excitation and development of stationary and travelling instability modes in a 3D boundary layer on a swept wing, as well as in 2D boundary layers including the Blasius flow and a self-similar boundary layer with an adverse pressure gradient. In particular, the following problems of the boundary-layer receptivity are discussed: (i) receptivity to localized 3D surface roughness, (ii) receptivity to localized 3D surface vibrations, (iii) acoustic receptivity in presence of 3D surface roughness, and (iv) acoustic receptivity in the presence of 3D surface vibrations. All experiments described in the paper were conducted using controlled disturbance conditions with the help of simulation of the stationary and non-stationary perturbations by means of several disturbance generators. This approach gives us the possibility to obtain quantitative results which are independent of any uncontrolled background perturbations of the flow and the experimental model. In contrast to the data obtained at “natural” environmental conditions these results can be directly compared with calculations without any significant assumptions about the physical nature of the disturbances under investigation. The complex (amplitude and phase) coefficients of the boundary-layer receptivity to external perturbations, obtained as functions of the disturbance frequency and the spanwise wavenumber (or the wave propagation angle), represent the main results of the experiments described. These results can be used for the evaluation of the initial amplitudes and phases of the instability modes generated by various external perturbations, as well as for quantitative verification of linear receptivity theories. Several examples of the comparison of experimental results with calculations are also presented in this paper. A brief analysis of the state-of-art in the field is performed and some general properties of different receptivity mechanisms are discussed.     

Non-Parallel Acoustic Receptivity of a Blasius Boundary Layer Using an Adjoint Approach

Localized and non-localized acoustic receptivity for a Blasius boundary layer is investigated using the adjoint Parabolized Stability Equations. The scattering of an acoustic wave onto a hump, a rectangular roughness or a wall steady blowing and suction is described. Comparisons with local approaches, triple deck theory, direct numerical simulations and experiments are successfully shown. Non-parallel effects are discussed. For comparable parameters, the non-localized receptivity problem produces amplitudes one order of magnitude larger than for the case of localized receptivity. This revised version was published online in July 2006 with corrections to the Cover Date.     

Boundary-layer receptivity due to distributed surface imperfections of a deterministic or random nature

Acoustic receptivity of a Blasius boundary layer in the presence of distributed, two-dimensional surface irregularities is investigated analytically. It is shown that, out of the entire spatial spectrum of the surface irregularities, only a narrow band of Fourier components can lead to an efficient conversion of the acoustic input at any given frequency to an unstable eigenmode of the boundary-layer flow. The location and the width of this most receptive band of wave numbers is fixed by the requirement of a relative detuning of O(R _inf1.b. ^sup–3/8 ) or less with respect to the instability wave number at the lower-branch station for the frequency under consideration. Surface imperfections in the form of discrete-mode waviness in this range of wave numbers then lead to initial instability amplitudes which are larger by a factor of O(R _inf1.b. ^sup3/8 ) than the amplitudes resulting from a single, isolated roughness element of streamwise extent comparable with the instability wavelength at the lower-branch location. In contrast, random irregularities which are spatially homogeneous in nature, and also possess a continuous spectrum in the streamwise direction, lead to instability amplitudes that are intermediate to those caused by the periodic and isolated irregularities, respectively, being, in fact, of the same order as the geometrical mean of the amplitudes in the latter two cases. A physical explanation for these asymptotic scalings is given, in addition to providing an analytical expression for the expected value of the instability amplitude for an ensemble of statistically irregular surfaces with random phase distributions. The duality between the localized and distributed receptivity analyses is also discussed.Financial support for this work was provided by the Theoretical Flow Physics Branch, Fluid Mechanics Division, NASA Langley Research Center, Hampton, VA, under contract NAS1-19299.     

Receptivity mechanisms for Görtler vortex modes

The receptivity problem for Görtler vortices induced by wall roughness or freestream disturbances is reviewed. To date, receptivity studies for this problem have been exclusively linear in character and here we show how the roughness and freestream disturbance mechanisms can each play dominant or inconsequential roles as possible routes to transition. The importance of each process is dependent on the exact situation at hand. For example, distributed wall-roughness elements tend to be more important in the generation of O(1) wave-number vortices than are isolated roughness patches whilst variations in the freestream velocity can easily provoke high wave-number disturbances on which roughness distributions typically have negligible effect.It has only been in recent times that the influence a spanwise component of the underlying basic boundary-layer flow may have on the Görtler mechanism has come to be appreciated. In some new computations we show that the imposition of such a spanwise component can lead to an increase in the coupling coefficient (a measure of the efficiency of a generating process) for modes provoked by wall roughness. However, such crossflow tends to reduce the global amplification rate of the most unstable mode so has the overall effect of restricting vortex growth downstream of any roughness element. We also demonstrate how the nonparallelism of Görtler vortices implies that conclusions concerning vortex receptivity properties can only be drawn after taking full account of upstream conditions and the precise form of the generating mechanism but it appears that for a large class of flows distributed wall forcing is more important in the provocation of modes than are either isolated roughness or freestream disturbances.This work was completed whilst APB was on study leave at the School of Mathematics, University of New South Wales, Sydney. He is indebted to the Royal Society of London and the Australian Research Council without whose grants (the latter to Dr. Peter Blennerhassett) his visit would not have been possible. In addition he is grateful to the staff and students of New College, UNSW for their provision of a Visiting Fellowship and to the School of Mathematics for their hospitality.     

On the Secondary Instability of Three-Dimensional Boundary Layers

One of the possible transition scenarios in three-dimensional boundary layers, the saturation of stationary crossflow vortices and their secondary instability to high-frequency disturbances, is studied using the Parabolized Stability Equations (PSE) and Floquet theory. Starting from nonlinear PSE solutions, we investigate the region where a purely stationary crossflow disturbance saturates for its secondary instability characteristics utilizing global and local eigenvalue solvers that are based on the Implicitly Restarted Arnoldi Method and a Newton–Raphson technique, respectively. Results are presented for swept Hiemenz flow and the DLR swept flat plate experiment. The main focuses of this study are on the existence of multiple roots in the eigenvalue spectrum that could explain experimental observations of time-dependent occurrences of an explosive growth of traveling disturbances, on the origin of high-frequency disturbances, as well as on gaining more information about threshold amplitudes of primary disturbances necessary for the growth of secondary disturbances. Received 13 July 1998 and accepted 7 July 2000     

波涡相互作用研究的某些进展(Ⅱ)

<正> 5 波涡共振 从第3节的感受性问题再前进一步,自然要问在什么条件下入射波激发起涡中之波的最大响应.这就导致了入射波与层状或轴状涡中受激波之间共振的概念,简称波涡共振.一般说来,在流体内部若有两个或多个波相会,它们将互相穿透而沿原来的方向离去.但若它们的波矢量和频率满足一定的关系(参见Craik 1985),就会在相会点产生新的波.2阶扰动的振幅可达到1阶扰动振幅的量级,而且流场中会出现一些重要的独特性质.这就是流体内部波共振,波涡共振是其一类情形.      

Linear stability of three-dimensional boundary layers

Stability of compressible three-dimensional boundary layers on a swept wing model is studied within the framework of the linear theory. The analysis based on the approximation of local self-similarity of the mean flow was performed within the Falkner-Skan-Cooke solution extended to compressible flows. The calculated characteristics of stability for a subsonic boundary layer are found to agree well with the measured results. In the case of a supersonic boundary layer, the results calculated for a Mach number M = 2 are also in good agreement with the measured spanwise scales of nonstationary vortices of the secondary flow. The calculated growth rates of disturbances, however, are substantially different from the measured values. This difference can be attributed to a high initial amplitude of disturbances generated in the experiment, which does not allow the linear stability theory to be applied. The evolution of natural disturbances with moderate amplitudes is fairly well predicted by the theory. The effect of compressibility on crossflow instability modes is demonstrated to be insignificant. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 2, pp. 3–14, March–April, 2008.     

Vortex-speed selection within the boundary-layer flow over a rotating sphere placed in an enforced axial flow

This paper furthers existing work into the instability mechanisms within the boundary-layer flow over a rotating sphere through the study of amplification rates within the convectively-unstable region. The onset of convective instability is associated with the experimentally observed onset of spiral vortices reported in the literature. Axial flow is found to stabilize the boundary layer by both delaying the onset of convective instability at all latitudes and also by significantly reducing the spatial amplification rates. We find that the type II (streamline curvature) mode becomes increasingly amplified with respect to the type I (crossflow) mode and is therefore likely to be selected in practice for sufficiently high axial flow rates. Furthermore, in experiments where special care is taken to remove all surface roughness, we predict that vortices will rotate at around 75% of the local surface speed. This is consistent with the experimental observations of Kobayashi & Arai who note a speed of around 76% under particular experimental conditions. These predictions are entirely consistent with related work on the rotating-disk and cone boundary layers.     

Stabilization of cylinder wakes in shallow water flows by means of roughness elements: an experimental study

Shallow wakes that occur in a wide range of natural flows tend to generate instabilities that develop into large, 2D coherent structures (2DCS). We present the results of an experimental study to stabilize shallow wakes by local, enhanced bottom roughness. Two successful stabilization strategies are compared to a base case of an unsteady bubble wake. First, localized bed roughness is placed in the lateral shear layers near the shoulders of the cylinder. Second, a local roughness element is placed at the end of the recirculation bubble, in the downstream region where large-scale vortices would normally shed. Dye visualization is used to assess the qualitative behavior of the wake, and two-component laser Doppler velocimetry (LDV) measurements are made to measure the Reynolds stress distributions and time-averaged velocity profiles. In both stabilization cases, a minimum patch size of the enhanced roughness elements is required for stabilization, which depends on the momentum thickness of the shear layers and the locations of enhanced Reynolds shear stresses. The main effect of the wake stabilization is a reduction in momentum exchange with the ambient flow due to damping of the large 2DCS. This reduction in eddy diffusivity results in a narrower wake and a slower decay of the centerline velocity deficit with downstream distance compared to the base case of an unsteady bubble wake.     

Effect of roughness on the velocity profile of a laminar boundary layer

Flow past a rough wall is examined. Calculations are made to find the roughness-induced mean velocity which is expressed in an integral form in terms of the spectral density of the roughness and an influence function. Values of the influence function are tabulated using the known values of the modified Hankel functions of order 1/3 and their integrals. The first order change in lower critical Reynolds number due to the roughness-induced change in profile is calculated; the stability of the profile is increased due to the presence of roughness. This work was supported by U.S. Naval Ordnance Systems Command under Contract NOw 65-0123-d through the Garfield Thomas Water Tunnel, Ordnance Research Laboratory.     

Disturbance evolution in the leading-edge region of a blunt flat plate in supersonic flow

The spatio-temporal dynamics of small disturbances in viscous supersonic flow over a blunt flat plate at freestream Mach number M_∞=2.5 is numerically simulated using a spectral approximation to the Navier–Stokes equations. The unsteady solutions are computed by imposing weak acoustic waves onto the steady base flow. In addition, the unsteady response of the flow to velocity perturbations introduced by local suction and blowing through a slot in the body surface is investigated. The results indicate distinct disturbance/shock-wave interactions in the subsonic region around the leading edge for both types of forcing. While the disturbance amplitudes on the wall retain a constant level for the acoustic perturbation, those generated by local suction and blowing experience a strong decay downstream of the slot. Furthermore, the results prove the importance of the shock in the distribution of perturbations, which have their origin in the leading-edge region. These disturbance waves may enter the boundary layer further downstream to excite instability modes.     

On the receptivity and nonparallel stability of traveling disturbances in rotating-disk flow

The generation and evolution of small amplitude long wavelength traveling disturbances in rotating-disk flow is the subject of this paper. The steady rotational speed of the disk is perturbed so as to introduce high-frequency oscillations in the flow field. Secondly, we introduce surface imperfections on the disk such as roughness elements. The interaction of these two disturbances will generate the instability waves whose evolution is governed by parabolic partial differential equations which are solved numerically. It is found that, for the class of disturbances considered here (wavelength on the order of Reynolds number), eigensolutions exist which decay or grow algebraically in the radial direction. However, these solutions grow only for frequencies larger than 4.58 times the steady rotational speed of the disk. The computed receptivity coefficient shows that there is an optimum size of roughness for which these modes are preferentially excited. The width of these roughness elements in the radial direction is about 0.1r ₀ ^* where r ₀ ^* is the radial location of the roughness. It is also found that the receptivity coefficient is larger for a negative spanwise wave number than for a positive one. The cumulative wave pattern produced from the roughness site shows that the typical wave angles for these disturbances are about –26° with about seven waves around the circumference. This is in contrast with the wave angles of 10°–14° found for the 30 or so inviscid cross-flow vortices.This work was sponsored by NASA Langley Research Center under Contracts NAS1-18240 (P.B. and M.R.M.) and NAS1-18605 (P.H.).     

Boundary-layer receptivity to external disturbances using multiple scales

Non-homogeneous multiple scales are introduced to solve the resonant problem of non-parallel boundary-layer receptivity originating from the quadratic mixing of environmental disturbances. The resulting algorithm is computationally inexpensive and can be efficiently included in industrial codes for transition prediction. The mutual interactions between acoustic wave, vorticity wave, wall vibration and wall roughness are discussed in detail and the receptivity coefficient, which relates the amplitude of the excited wave to the amplitude of the exciting sources, is computed. The largest effect is found for the interaction between acoustic waves and wall roughness perturbations. Other coupling mechanisms are less effective. By comparing parallel and non-parallel results, it is found that flow non-parallelism can play a non-negligible role even in Blasius’ boundary layer, although the largest effects are evident for the three-dimensional boundary layer over an infinite swept wing. For the particular case of wall roughness—wall vibration mixing, the velocity disturbance is shown to be exactly equal to the velocity perturbation induced by wall roughness alone on a wall vibrating in the normal direction.     

Traveling disturbance appearing in boundary layer transition in a Yawed cylinder

Boundary layers that develop over a body in fluid flow are in most cases three-dimensional owing to the spin, yaw, or surface curvature of the body. Therefore, the study of three-dimensional (3D) boundary-layer transition is essential to work in practical aerodynamics. The present investigation is concerned with the problem of 3D boundary layers over a yawed body. A yawed cylinder model that represents the leading edge portion of a swept wing and the mechanism of crossflow instability are investigated in detail using hot-wire velocimetry and a flow visualization technique. As a result, traveling disturbances having frequencies f₁ and f₂, which differ by about one order of magnitude, are detected in the transition region. The phase velocities and directions of travel of those disturbances are measured. Results for the low-frequency disturbance f₁ show qualitative coincidence with results numerically predicted for a crossflow unsteady disturbance. Nameley, F₁ travels nearly spanwise to the yawed cylinder and very close to the cylinder wall. The results for the high-frequency disturbance f₂ good agreement with the existing experimental results. The ₂ disturbance is found to be the high-frequency inflectional secondary instability that appears in 3D boundary layer transition in general. A two-stage transition process, where stationary crossflow vortices appear as the primary instability and a traveling inflectional disturbance is generated as a secondary instability, was observed. Secondary instability seems to play a major role in turbulent transition.     

Large Eddy Simulation of boundary layer transition over an isolated ramp-type micro roughness element

Boundary layer transition over an isolated surface roughness element is investigated by means of numerical simulation. Large Eddy Simulation (LES) flow-modeling approach is employed to study flow characteristics and transition phenomenon past a roughness element immersed within an incoming developing boundary layer, at a height-based Reynolds number of 1170. LES numerical results are compared to experimental data from literature showing the time-averaged velocity distribution, the velocity fluctuation statistics and the instantaneous flow topology.Despite slight difference in the intensity of streamwise velocity fluctuations, the present LES results and experimental data show very good agreement. The mean flow visualization shows streamwise counter-rotating vortices pairs formation downstream of the obstacle. The primary pair induces an upwash motion and a momentum deficit that creates a Kelvin-Helmholtz type flow instability. The instantaneous flow topology reveals the formation of coherent K-H vortices downstream that produce turbulent fluctuations in the wake of the roughness element. These vortices are streched and lifted up when moving downstream. The velocity fluctuations results show that the onset of the turbulence is dominated by the energy transfer of large-scale vortices.     

Stability of the laminar boundary layer flow encountering a row of roughness elements: Biglobal stability approach and DNS

The stability of the laminar boundary layer developing on a flat plate in the presence of a periodic row of roughness elements is investigated. A Direct Numerical Simulation is performed to compute the steady flow downstream of the roughness elements, which contains a pair of two counter-rotating streamwise vortices per element, which can be considered as a “pre-streaky” structure. The linear stability of this base flow is analyzed by means of the so-called “biglobal” stability approach. Three-dimensional eigenmodes are found, which are shown to be the continuation of the Tollmien–Schlichting waves present in the case of an unperturbed boundary layer. Moreover, a stabilizing effect due to the roughness-induced vortices is found. A Direct Numerical Simulation of the interaction between a two-dimensional Tollmien–Schlichting wave and the roughness array is also performed. The computed perturbation traveling downstream of the roughness elements is shown to be a linear combination of the biglobal eigenmodes.     

On the etiology of shear layer vortices

In an effort to isolate the mechanism by which streamwise structures form in turbulent wall layers, evolution equations were derived for the streamwise velocity and vorticity perturbations about a mean turbulent fully developed channel flow. The stability of these equations, which take their most concise form when derived from the Generalized Lagrangian mean equations of Andrews and McIntyre, are studied assuming normal modes and infinitesimal disturbances. The resulting stability diagram yields, inter alia, the spanwise periodicity of the resulting structures, which we term shear layer vortices. If streaks are thought of as the footprints of these vortices, we then have a formal way of determining the spacing of streaks. The first three modes of instability are determined; at the first not just two vortices form per period, but four. It is also evident that an intense local shear layer forms about the plane in which the convection velocity equals the mean Eulerian velocity.Dedicated to Professor J.L. Lumley on the occasion of his 60th birthday.This work was supported in part by the U.S. Office of Naval Research under SRO IV Grant No. N00014-85-K-0172 and in part by the National Science Foundation Grant CTS-9008477.     

Görtler Vortices with System Rotation

The steady primary instability of Görtler vortices developing along a curved Blasius boundary layer subject to spanwise system rotation is analysed through linear and nonlinear approaches, to clarify issues of vortex growth and wavelength selection, and to pave the way to further secondary instability studies.A linear marching stability analysis is carried out for a range of rotation numbers, to yield the (predictable) result that positive rotation, that is rotation in the sense of the basic flow, enhances the vortex development, while negative rotation dampens the vortices. Comparisons are also made with local, nonparallel linear stability results (Zebib and Bottaro, 1993) to demonstrate how the local theory overestimates vortex growth. The linear marching code is then used as a tool to predict wavelength selection of vortices, based on a criterion of maximum linear amplification.Nonlinear finite volume numerical simulations are performed for a series of spanwise wave numbers and rotation numbers. It is shown that energy growths of linear marching solutions coincide with those of nonlinear spatially developing flows up to fairly large disturbance amplitudes. The perturbation energy saturates at some downstream position at a level which seems to be independent of rotation, but that increases with the spanwise wavelength. Nonlinear simulations performed in a long (along the span) cross section, under conditions of random inflow disturbances, demonstrate that: (i) vortices are randomly spaced and at different stages of growth in each cross section; (ii) upright vortices are the exception in a universe of irregular structures; (iii) the average nonlinear wavelengths for different inlet random noises are close to those of maximum growth from the linear theory; (iv) perturbation energies decrease initially in a linear filtering phase (which does not depend on rotation, but is a function of the inlet noise distribution) until coherent patches of vorticity near the wall emerge and can be amplified by the instability mechanism; (v) the linear filter represents the receptivity of the flow: any random noise, no matter how strong, organizes itself linearly before subsequent growth can take place; (vi) the Görtler number, by itself, is not sufficient to define the state of development of a vortical flow, but should be coupled to a receptivity parameter; (vii) randomly excited Görtler vortices resemble and scale like coherent structures of turbulent boundary layers.A.Z. has been supported, during his stay at EPFL, by an ERCOFTAC Visitor Grant. A.B. acknowledges the Swiss National Fund, Grant No. 21-36035.92, for travel support associated with this research. This work was also supported by the Swedish Board of Technical Development (NUTEK), the Swedish Technical-Scientific Council (TFR), and an ERCOFTAC Visitor Grant, through which the stay of B.G.B.K. at the EPFL was made possible. Cray-2 computing time for this research was generously provided by the Service Informatique Centrale of EPFL.     

Mechanism of three-dimensional boundary-layer receptivity

Boundary-layer receptivity is always a hot issue in laminar-turbulent transition. Most actual laminar-turbulent transitions belong to three-dimensional flows. An infinite back-swept flat-plate boundary layer is a typical three-dimensional flow. Study of its receptivity is important both in theory and applications. In this paper, a freestream turbulence model is established. A modified fourth-order Runge-Kutta scheme is used for time marching, and compact finite difference schemes are used for space discretization. On these bases, whether unsteady cross-flow vortices can be excited in the three-dimensional boundary layer(the infinite back-swept flat-plate boundary layer) by free-stream turbulence is studied numerically. If so, effects of the level and the direction of free-stream turbulence on the three-dimensional boundary-layer receptivity are further studied. Differences of the three-dimensional boundary-layer receptivity are then discussed by considering the non-parallel effect, influence of the leading-edge stagnation point of the flat plate, and variation of the back-swept angle separately. Intensive studies on the three-dimensional boundary-layer receptivity will benefit the development of the hydrodynamic stability theory, and provide a theoretical basis for prediction and control of laminar-turbulent transition.     

Local receptivity in non-parallel boundary layer

The research on boundary-layer receptivity is the key issue for the laminarturbulent transition prediction in fluid mechanics. Many of the previous studies for local receptivity are on the basis of the parallel flow assumption which cannot accurately reflect the real physics. To overcome this disadvantage, local receptivity in the non-parallel boundary layer is studied in this paper by the direct numerical simulation (DNS). The difference between the non-parallel and parallel boundary layers on local receptivity is investigated. In addition, the effects of the disturbance frequency, the roughness location, and the multiple roughness elements on receptivity are also determined. Besides, the relations of receptivity with the amplitude of free-stream turbulence (FST), with the roughness height, and with the roughness length are ascertained as well. The Tollmien- Schlichting (T-S) wave packets are excited in the non-parallel boundary layer under the interaction of the FST and the localized wall roughness. A group of T-S waves are separated by the fast Fourier transform. The obtained results are in accordance with Dietz’s measurements, Wu’s theoretical calculations, and the linear stability theory (LST).