Open access peer-reviewed chapter - ONLINE FIRST

Leading Edge Receptivity at Subsonic and Moderately Supersonic Mach Numbers

By Marvin E. Goldstein and Pierre Ricco

Submitted: October 11th 2018Reviewed: December 19th 2018Published: January 30th 2019

DOI: 10.5772/intechopen.83672

Downloaded: 138

Abstract

This chapter is a review of the receptivity and resulting global instability of boundary layers due to free-stream vortical and acoustic disturbances at subsonic and moderately supersonic Mach numbers. The vortical disturbances produce an unsteady boundary layer flow that develops into oblique instability waves with a viscous triple-deck structure in the downstream region. The acoustic disturbances (which have phase speeds that are small compared to the free stream velocity) produce boundary layer fluctuations that evolve into oblique normal modes downstream of the viscous triple-deck region. Asymptotic methods are used to show that both the vortically and acoustically-generated disturbances ultimately develop into modified Rayleigh modes that can exhibit spatial growth or decay depending on the nature of the receptivity process.

Keywords

  • boundary layer
  • boundary layer receptivity
  • compressible boundary layers
  • global instability

1. Introduction

This chapter is concerned with the effect of unsteady free-stream disturbances on laminar to turbulent transition in boundary layer flows. The exact mechanism depends on the nature and intensity of the disturbances. Transition at high disturbance levels (say >1%) usually begins with the excitation of low frequency streaks in the boundary layer flow that eventually break down into turbulent spots. This phenomena was initially studied by Dryden [1] and much later for compressible flows by Marensi et al. [2]. But the focus of this chapter is on low free steam disturbances levels (say less than 1%) where the transition usually results from a series of events beginning with the generation of spatially growing instability waves by acoustic and/or vortical disturbances in the free-stream. This so-called receptivity phenomenon results in a boundary value problem and therefore differs from classical instability theory which results in an eigenvalue problem for the Rayleigh or Orr-Sommerfeld equations that only apply when the mean flow can be treated as being nearly parallel (see, for example, Reshotko, [3]). The relevant boundary conditions cannot be imposed on the Orr-Sommerfeld or Rayleigh equations in the infinite Reynolds number limit being considered here but the free-stream disturbances can produce unsteady boundary layer perturbations in regions of rapidly changing mean flow that eventually produce unstable Rayleigh or Orr-Sommerfeld equation eigensolutions further downstream. These regions of nonparallel flow can result from surface roughness elements [4, 5], blowing or suction effects [6] or from the nonparallel mean flow that occurs near the boundary layer leading edge [7, 8].

The mechanism is similar in all cases but the simplest and arguably the most fundamental of these is the one resulting from the nonparallel leading edge flow and the focus here is, therefore, on that case. The initial studies were carried out for two dimensional incompressible flows. Ref. [7] used a low frequency parameter matched asymptotic expansion to show that there is an overlap domain where appropriate asymptotic solutions to the forced boundary layer equations (which apply near the edge) match onto the so-called Tollmien-Schlichting waves that satisfy the Orr-Sommerfeld equation in a region that lies somewhat further downstream. The coupling to the free-stream disturbances turns out to be fairly weak for the two dimensional incompressible flow considered in [7] due to the relatively large decay of boundary layer disturbances upstream of the Tollmien-Schlichting wave region where the Orr-Sommerfeld equation applies.

But there can be a much stronger coupling in supersonic flows which can support a number of different instabilities [9]. The coupling mechanism can be either viscous or inviscid and the instability can either be of the viscous Tollmien-Schlichting type or can be purely inviscid when the mean boundary layer flow has a generalized inflection point. The inviscid coupling, which was first analyzed in [10], tends to be dominant when the obliqueness angle θof the disturbance differs from the critical angle, θccos11/M, where the Mis the free-stream Mach number, by an O1amount. Figure 1 shows that the theoretical results of Ref. [10] are in good agreement with experimental data when Δθθcθ=O1but the agreement breaks down when θθc[12] and a new rescaled analysis was carried out in Ref. [11] to deal with this case.

Figure 1.

Comparison of the Fedorov/Khokhlov solution with experiment [12].

Fedorov and Khokhlov [10] analyzed the generation of inviscid instabilities in a supersonic flat plate boundary layer by fast and slow acoustic disturbances in the free stream. They showed that the slow acoustic mode propagates downstream/upstream when the obliqueness angle θof the acoustic disturbances is smaller/larger than the critical angle θcand that downstream propagating slow acoustic modes with Δθ>0generate unsteady boundary layer disturbances that match onto the inviscid 1st Mack mode instability without undergoing any significant decay. The focus of that reference was on hypersonic flows while the interest here is in the moderately supersonic regime (Mach number less than 4), where the so called 1st Mack mode is the dominant instability, but (as shown in Section 6) emerges much too far downstream to be of practical interest when generated by the inviscid mechanism analyzed in [7]. The instability produced by the small Δθanalysis of Ref. [11] can, however, occur much further upstream when Δθis sufficiently small. But there is a smallest value of Δθfor which the instability wave coupling can occur.

Smith [13] showed that viscous instabilities, which exhibit the same triple-deck structure as the subsonic Tollmien-Schlichting waves, can also occur at supersonic speeds when the obliqueness angles θis greater than the critical angle θc. Their phase speeds are very small and they must therefore be produced by a viscous wall layer mechanism similar to the one identified in [7].

The analysis of Ref. [7] was extended to compressible subsonic and supersonic flat plate boundary layer flows by Ricco and Wu [14] who showed that highly oblique vortical disturbances can generate a limiting form of the Smith instability [13]. They found that the instability wave lower branch lies further upstream at supersonic speeds than the subsonic lower branch and much further upstream than the incompressible lower branch considered in [7], which means that the instability wave/free-stream disturbance coupling is much greater at supersonic speeds than it is in the incompressible flow considered in [7]. Goldstein and Ricco [11] show that the instability does not possess an upper branch in this case and matches onto a low frequency (short streamwise wavenumber) Rayleigh instability (that can be identified with the 1st Mack mode) when the downstream distance is slightly smaller than the downstream distance where acoustically generated instability corresponding to the smallest possible Δθemerges. It therefore makes sense to consider both of these receptivity mechanisms simultaneously.

As noted above, the present chapter is concerned with the unsteady flow in a flat plate boundary layer generated by mildly oblique vortical disturbance and small Δθacoustic disturbances in a moderately supersonic Mach number free stream. The results are expected to be relevant to transition in the straight wing boundary layers on supersonic aircraft such as the low-sweep Aerion AS2 Bizjet, shown in Figure 2.

Figure 2.

Low-sweep Aerion AS2 supersonic Bizjet. M∞≤1.5. Posted by Tim Brown on the Manufacturer Newsletter.

2. Imposed free-stream disturbances

Since the boundary layer is believed to be convectively unstable, the receptivity phenomena are best illustrated by considering a small amplitude harmonic distortion with angular frequency ωsuperimposed on a subsonic or moderately low Mach number supersonic flow of an ideal gas past an infinitely thin flat plate with uniform free-stream velocity U, temperature T, dynamic viscosity μand density ρ. The velocities, pressure fluctuations, temperature and dynamic viscosity are normalized by U, ρU2, Tand μ, respectively. The time tis normalized by ωand the Cartesian coordinates, say xyz, are normalized by LU/ωwith the coordinate ybeing normal to the plate.

As noted above the phenomenon is analyzed by requiring the Reynolds number Re=ρUL/μto be large, or equivalently requiring the frequency parameter F1/Reto be small, and using asymptotic theory to explain how the imposed harmonic distortion generates oblique instabilities at large downstream distances in the viscous boundary layer that forms on the surface of the plate. The natural expansion parameter turns out to be

εF1/6.E1

The free-steam disturbances will be inviscid at the lowest order of approximation and, as is well known [15], can be decomposed into an acoustic component that carries no vorticity, and vortical and entropic components that produce no pressure fluctuations. But only the first two will be considered here.

The vortical disturbance uvis given

uv=uvvvwv=δ̂uvwexpixt+γy+βz,E2

where δ̂1is a common scale factor and u,v,wsatisfy the continuity condition

u+γv+βw=0E3

but are otherwise arbitrary constants while the acoustic component is governed by the linear wave equation which has a fundamental plane wave solution

uapa=uavawapa=δ̂1ααγβ1αeiαx+γy+βzt,E4

for the velocity and pressure perturbation uapawhere

γ=M21αα1αα2,α1,2=M2±M2+β2M21M21E5

and, as noted in Section 1, Mdenotes the free-stream Mach number.

The leading edge interaction will produce large scattered fields for O1values of the incidence angles tan1va/ua=tan1γ/αand tan1vv/uvof the acoustic and vortical disturbances, respectively. And, in order to focus on the fundamental mechanisms, we assume that the incidence angles of the vortical disturbances are small and that the incidence angles of the acoustic disturbances are zero, which requires that

v/u1E6

for the former disturbances and that

α=α=Mcosθ/Mcosθ1,θtan1β/α,E7

for the latter, where the subscripts −/+ refer to the slow/fast acoustic modes. Eq. (7) shows that the slow mode wavenumber becomes infinite when the obliqueness angle is equal to the critical angle referred to in the introduction.

3. Boundary layer disturbances

As indicated above our interest here is in explaining how the incident harmonic distortions generate oblique instabilities at large downstream distances in the viscous boundary layer that forms on the surface of the plate. We begin by considering the fluctuations imposed on this flow by the free-stream vortical disturbance (2).

3.1 Boundary layer disturbances generated by the free-stream vorticity

As noted in the introduction, these disturbances will generate oblique Tollmien-Schlichting instability waves which are known to exhibit a triple-deck structure in the vicinity of their lower branch which lies at an Oε2distance downstream [13] of the leading edge in the high Reynolds number flow being considered here. The Tollmien-Schlichting waves will have Oε1spanwise wavenumbers and we therefore require that

β¯εβ=O1E8

since the spanwise wavenumber must remain constant as the disturbances propagate downstream.

The continuity condition (3) and the obliqueness restriction (6) will be satisfied if we put

w¯w/ε=O1,v¯v/ε=O1,γ¯εγ=O1.E9

The vortical velocity (2) will then interact with the plate to produce an inviscid velocity field [12] that generates a slip velocity at the surface of the plate which must be brought to zero in a thin viscous boundary layer whose temperature, density and streamwise velocity, say Tη,ρη,Uη, respectively, are assumed to be functions of the Dorodnitsyn-Howarth variable

η1ε32x0yρxy˜dy˜E10

and are determined from the similarity equations given in Stewartson [16] and Ref. [14].

We begin by considering the flow in the vicinity of the leading edge where the streamwise length scale is x=O1. Since the inviscid velocity field can only depend on the streamwise coordinate through this relatively long streamwise length scale the solution for the velocity and temperature perturbation uuvwϑin this region is given by [14], [17]

u=δ̂uu¯v¯0ϑ¯+β¯w¯+iv¯u¯0v¯0w¯0ϑ¯0eiβ¯z/εt,E11

where u¯0xηv¯0xηw¯0xηϑ¯0xηsatisfies the three dimensional compressible linearized boundary layer equations (with unit spanwise wavenumber) subject to the boundary conditions [14]

u¯0,ϑ¯00,w¯0eix,asη,E12

while u¯xηv¯xη0ϑ¯xηexpiβ¯z/εtis a quasi-two dimensional solution that satisfies the two dimensional linearized boundary layer equations subject to the boundary conditions

u¯eix,w¯,ϑ¯0asη.E13

The lowest order triple-deck solution will match onto the quasi-two dimensional solution u¯v¯0ϑ¯expiβ¯z/εtof the two dimensional boundary layer equations, where the spanwise dependence only enters parametrically through the exponential factor in(11) .

Prandtl [18], Glauert [19] and Lam and Rott [20] showed that

u¯xη=BxUηT2x,ϑ¯xy=BxTηTη2x,E14
v¯xη=iBx+dBdxUηBxUηηc2x,E15

where

ηc1Tη0ηTη˜dη˜E16

is an exact eigensolution of the two-dimensional linearized unsteady boundary layer equations that satisfies the homogeneous boundary conditions u¯xη,w¯xη,ϑ¯xη0asηfor all Bx, but does not necessarily satisfy the no-slip condition at the wall.

Lam and Rott [20], [21] analyzed the two dimensional flat plate boundary layer and showed that the linearized equations possess asymptotic eigensolutions that satisfy a no-slip condition at the wall when xbecomes large. These solutions exhibit a two-layer structure consisting of an outer region that encompasses the main part of the boundary layer and a thin viscous region near the wall. The outer solution is given by (14) and (15) with the arbitrary function Bxdetermined by matching with the viscous wall layer flow.

Ref. [14] showed that the Lam and Rott [20, 21] analysis also applies to compressible flows when the full compressible solution (14) and (15) is used in the outer region and the viscous wall layer solution is slightly modified to account for the temperature and viscosity variations. The function Bxis then given by

Bx=x3/2Bnexp23/2e/43λςn3/2Twμw1/2x3/2+.E17

where TwT0, μwμT0, λU0and ζnis a root of

Aiςn=0,forn=0,1,2,3.E18

The only difference from the Lam-Rott result is the Tw/μw1/2factor in the exponent. The asymptotic solution to the full inhomogeneous boundary value problem can now be expressed as the sum of a Stokes layer solution plus a number of these asymptotic eigensolutions. The first few Bnwere determined from numerical solutions to the boundary layer problem in Ref. [8]. But we are primarily concerned with the lowest order n=0mode because that is the only one that matches onto a spatially growing oblique Tollmien-Schlicting wave further downstream [11]. The receptivity problem can then be solved by combining the numerical computations with appropriate matched asymptotic expansions to relate the instability wave amplitude to that of the free-stream disturbance. But we will analyze the boundary layer disturbances generated by the free-stream acoustic disturbances before considering these expansions.

3.2 Boundary layer disturbances generated by the Fedorov/Khokhlov mechanism for obliqueness angles close to critical angle

Fedorov and Khokhlov [10] used matched asymptotic expansions to analyze the generation of Mack mode instabilities by oblique acoustic waves of the form (4) where the wavenumbers αand βsatisfy the dispersion relation (7) when the incidence angle γis equal to zero, which, as noted, above is the case being considered here. Their focus was on hypersonic flows where the most rapidly growing disturbances are usually two dimensional 2nd Mack modes, while, as noted in the introduction, the focus of the present chapter is on the relatively low supersonic Mach number regime (say, less than about 4) where the most rapidly growing instability waves are highly oblique 1st Mack modes. Numerical results [9] show that the obliqueness angle of the most rapidly growing 1st mode lies between 50 and 70 degrees at Mach numbers between 2 and 6.

Ref. [10] shows that the boundary layer disturbance produced by diffraction of the slow acoustic wave by the nonparallel mean flow in the region where x=oε3can be matched onto a 1st Mack mode instability in the downstream region where x=Oε6when the deviation

ΔθθcθE19

of the obliqueness angle θfrom the critical angle

cosθc1/ME20

takes on O1positive values. The diffraction region has a double layer structure which consists of a region that fills the mean boundary layer and an outer diffraction region of thickness O1/ε3/2. (The purely passive Stokes layer near the wall does not play a role in the diffraction process and can be ignored).

The instability emerges from the downstream limit of the solution in this region. But as noted in the introduction this occurs too far downstream to be of practical interest when scaled up to actual flight conditions if Δθ=O1[14] at the moderately supersonic Mach numbers being considered here. It will however emerge much further upstream when θis close to the critical angle θc, i.e., when Δθ1. But the solution in Ref. [10] does not apply when Δθ1and a new analysis was developed in Ref. [11] to extend their result into the small -Δθregime.

It follows from (7) that

α=α˜/Δθ+α˜1+,β=β1=β˜/ΔθE21

where

α˜1/tanθc,β˜1,α˜11/sin2θcE22

when Δθ1since tanθcΔθ=tanθcΔθ/cos2θc+OΔθ2in that case.

This shows that αalso becomes large when Δθ1and that αwill expand in powers of Δθas indicated in (21) if βis fixed at the indicated value to all orders in Δθ(which we now assume to be the case).

The spanwise wavenumber will equal the vortical spanwise wavenumber (8) when Δθ=Oεand as in that case the diffraction wave solution will eventually develop a triple-deck structure but the resulting solution will (as shown in [11]) not decay at large wall normal distances and is therefore invalid. This means that the diffraction region solution cannot be continued downstream for Δθ=Oε.

Figure 3.

Structure of diffraction region for Δθ=Oε2/3.

Ref. [11] shows that the smallest value of Δθis Δθ=Oε2/3and the diffraction region will then occur at an Oε4/3distance downstream. The relevant solution will have the triple-deck structure shown in Figure 3: a main boundary layer region that fills the mean boundary layer (region 1), an outer diffraction region of thickness Oε1/3(region 2) and an Oε3thick viscous wall layer in which the unsteady, convective and viscous terms all balance.

The pressure in region 2 is of the form

p=1+δ̂p2x2y2eiα˜/Δθ+α˜1x+β˜z/Δθt,E23

where

x2xε4/3=O1,y2yε1/3=O1E24

and the surface pressure p2x20is related to the up-wash velocity v1x2limηv1x2ηat the outer edge of the boundary layer by

p2x20=p1x2=1x22πiα˜M2101σ1σiα˜v1x2σx2σ,E25

where p1x2denotes the pressure in the boundary layer region 1 (which is independent of the wall normal direction) and the wall normal velocity v1x2is given in terms of

ξ2i1/32x2/α˜λ2/3Tw/μw1/3E26

and the integral and the derivative of the Airy function Aiξby

v1x22x2=ip1x2α˜2+β˜2Tw2ξ2λAiξ2ξ2Aiξ,E27

which behaves like

v1x2/2x2ip1x2α˜2+β˜2Tw2/λE28

as x2since ([22], pp. 446–447)

Aiξ/ξAiqdqξasξ.E29

Inserting (28) and (27) into (25) shows that

p1x2=1γ0x201σ1σp1σx2,asx2E30

where

γ0α˜2+β˜2α˜1/2Tw2λ2πiM21,E31

which is formally the same as the equation considered in [10] who showed that the solution behaves like

p1x2expγ02πx22asx2.E32

The acoustically and vortically generated boundary layer disturbances considered in this section will eventually evolve into propagating eigensolutions in regions that lie further downstream. The resulting flow will have a triple-deck structure of the type considered in [13], [23] and [14] in the former (i.e., vortically generated) case. But the acoustically generated disturbance will only develop an eigensolution structure much further downstream. The minimum distance occurs when Δθ=Oε2/3. We begin by considering the triple-deck region.

4. The viscous triple-deck region

Refs. [13, 14, 23] show that the linearized Navier-Stokes equations possess an eigensolution of the form

uvwp=δ̂Πyεei1ε30x1κx1εdx1+β¯z¯tE33

in the triple-deck region where δ̂1is the common scale factor introduced at the beginning of Section 2,

Πyε=A¯x1UηTηiκ0A¯x1Uη2x1,ε2β¯P¯Tηκ0Uηε2P¯E34

in the main boundary layer where η=O1,

x1ε2x=O1E35

and

z¯z/ε=zω/εUE36

is a scaled transverse coordinate. The complex wavenumber κhas the expansion [11].

κx1ε=κ0x1+εκ1x1+ε2κ2x1+.,E37

where the lowest order term in this expansion satisfies the following dispersion relation ([13, 14, 23])

κ02+β¯2=1iκ01/3λ2x15/3μwTw71/3β¯2M21κ021/2Aiξ0ξ0AiqdqE38

and

ξ0=i1/32x1κ0λ2/3Tw/μw1/3E39

whose solution must satisfy the inequality

Reβ¯2M21κ021/20E40

in order to insure that the eigensolution does not exhibit unphysical wall normal growth.

This requirement will be satisfied for all M<1but will only be satisfied at supersonic Mach numbers when the obliqueness angle θis greater than the critical angle θccos11/M[11, 13]. The dispersion relation (38) and (39) reduces to the dispersion relation given by Eqs. (4.52), (5.2) and (5.3) of [7] when β¯and Mare set to zero.

4.1 Matching with the Lam-Rott solution

The dispersion relation (38) and (39) will be satisfied at small values of x1if κ0x1and ξ0ζn,forn=0,1,2 asx10,where ςnis the nth root of the Lam-Rott dispersion relation (18). Inserting this into (38) shows that

κ01λςn3/22Twx1iμw1/2asx10.E41

The cross flow velocity wdrops out of (33) as x10and the flow in the main deck is therefore compatible with the quasi-two dimensional Lam-Rott solution (14)(17).

4.2 Numerical results

The dispersion relation (38) is expected to have at least one root corresponding to each of the infinitely many roots of (18). But only the lowest order n=0root can produce the spatially growing modes of (38). The wall temperature Twand viscosity μwcan be scaled out of this equation by introducing the rescaled variables.

κ0=κ0Tw1/2μw1/6,x1=x1Tw2/μw2/3,β¯=β¯Tw1/2μw1/6.E42

The real and negative imaginary parts of κ0calculated from (38) together with the n=0Lam-Rott initial condition (41) are plotted as a function of the scaled streamwise coordinate x¯1in Figures 4 and 5 for three values of the frequency scaled transverse wavenumber β¯2. The insets are included to more clearly show the changes at small x¯1. The dashed curves in the insets denote the real and imaginary parts of the small-x¯1asymptotic formula(41) .

Figure 4.

Reκ0† as a function of x1† calculated from (38) together with the initial condition (41) for M∞=2,3,4 (double dot dashed, dot dashed, and solid lines, respectively) and three values of β¯†≥2. The dashed curve in the main graph is the rescaled large-x1† asymptote (49).

Figure 5.

−Imκ0† as a function of x1† calculated from (38) together with the initial condition (41) for M∞=2,3,4 (double dot dashed, dot dashed and solid lines, respectively) and three values of the frequency scaled transverse wavenumber.

The triple-deck eigensolution (33) (which contains the Lam-Rott solution as an upstream limit) can undergo a significant amount of damping before it turns into a spatially growing instability wave at the lower branch of the neutral stability curve.

The exponential damping in Eq. (33) is proportional to Im0xLBκx1dx=ε2Im0x1LBκx1dx1, where x1LBand xLBdenote the scaled and unscaled streamwise location of the lower branch of the neutral stability, which implies that the total damping is proportional to the area under the growth rate curve between zero and the lower branch in Figure 5. The inset shows that the length Δx1=0.01of this upstream region is very short and therefore that the total amount of damping is relatively small.

5. The inviscid triple-deck region

As noted above the acoustically driven solution will only match onto an eigensolution in the downstream region when OΔθε3/2. This region will lie downstream of the viscous triple-deck region considered above and will be closest to that region when OΔθ=ε3/2. It will have an inviscid triple- deck structure and the relevant dispersion relation can be obtained by putting ε/Δθ=Oε1/3in (21), inserting the rescaled variables

β¯¯=β¯/ε1/3,κ¯0=κ0/ε1/3,x̂1=x1ε4/3E43

into (38), using (29), and taking the limit as ε0with β¯¯, κ¯0and x̂1held fixed, to show that the rescaled wavenumber κ¯0satisfies the inviscid dispersion relation

κ¯02+β¯¯2=λβ¯¯2M21κ¯021/2κ¯02x̂1Tw2E44

when the square root β¯¯2M21κ¯021/2is required to remain finite as ε0.

5.1 Matching with the small ΔθFedorov/Khokhlov solution

It can then be shown by direct substitution that the solution κ¯0behaves like

κ¯0β¯¯M211/2β¯¯5α̂02x̂1asx̂10,E45

where α̂0M2Tw2/M217/4λ. The square root β¯¯2M21κ¯021/2still satisfies the inequality (40) when x̂10and (44) therefore remains valid in this limit.

The pressure component of the resulting solution will then match onto the downstream limit (32) and (30) of the acoustically generated diffraction region solution when β¯¯=Oε2/3/Δθand x2is given by (24) since it follows from (8),(35),(43) and (45) that

1/ε30x1κ0x1dx1=1/ε40x̂1κ¯0x̂1dx̂1α˜/Δθxεβ¯5α̂02x2/2=α˜/Δθxβ5α̂02ε3x2/2=α˜/Δθxα̂02x22/2.E46

5.2 Numerical results

Figure 6 is a plot of the scaled lowest order wavenumber κ¯0/β¯¯=κ0/β¯as a function of the scaled streamwise coordinate β¯¯Tw4x̂1/λ2=β¯Tw4x1/λ2for various values of the free-stream Mach number Mcalculated from the inviscid triple-deck dispersion relation(44) together with the asymptotic initial condition (45) which is shown by the dashed curves in the figure. The lowest order wave number κ¯0is purely real which means that exponential growth (if it occurs) can only occur at higher order. This suggests that the acoustically generated instabilities will be less significant than the vortically-generated instabilities which appear upstream.

Figure 6.

Scaled wavenumber κ¯0/β¯¯=κ0/β¯ as a function of β¯¯Tw4x̂1/λ2=β¯Tw4x1/λ2 for various values of M∞. The solid lines represent the numerical solution. Dashed lines are the asymptotic solution(45) .

6. The next stage of evolution

6.1 Downstream behavior of the triple-deck solution

Eqs. (29), (38) and (39) show that

β¯1κ01/3λ2x15/31Tw22x1κ0λ2/3=λκ0Tw22x1E47

when x1and, therefore, that

κ0λβ¯Tw22x1,E48

when κ0is allowed to approach zero as x1.

The dashed curves in the main plot of Figure 4 represent the re-scaled large-x¯1asymptote (48). It confirms that the numerical results are well approximated by the (appropriately rescaled) large-x1asymptote (48).

As noted in [11], the solution to the reduced dispersion relation (44) satisfies the rescaled version

κ¯0λβ¯¯Tw22x̂1asx̂1E49

of (48), which can be considered to be a special case of this result if we put

β¯¯=β¯/εr,κ¯0=κ0/εr,x̂1=x1ε4rE50

and allow r to be zero or 1/3.

The expansion (37) then generalizes to [11]

κ¯x1ε=κ¯0x̂1+ε1rκ¯1x̂1+ε21rκ¯2x̂1+.,E51

where

κ¯,κ¯1,κ¯2κ/εr,κ1,κ2εrE52

and x̂1is defined in (43).

6.2 Derivation of the governing equations

Eq. (49) shows, among other things, that the lowest order wave number and streamwise growth rate approach zero but do not become negative as the disturbance propagates downstream. The boundary layer thickness which is Oε3xcontinues to increase and the triple-deck scaling breaks down at the streamwise location

x¯1=xε4+2r=O1,E53

where it becomes of the order of the spanwise length scale, which remains constant at Oε1r. This region is located well upstream of the region where the unsteady flow is governed by the full Rayleigh equation considered in [9].

Eqs. (37), (43), (51) and (52) show that the Tollmien-Schlichting wave becomes more oblique and

expi1ε30x1κx1εdx1+β¯z¯t=expi1ε31+r0x̂1κ¯0x̂1εdx̂1+1ε2+4r0x̂1κ¯1x̂1εdx̂1+1ε1+5r0x̂1κ¯2x̂1εdx̂1+Oε4r+εrβ¯¯z¯teiε22+r0x¯1α¯x¯1εdx¯1+β¯¯z¯¯tE54

as x̂1,where α¯x¯1is an O1function of x¯1(given by (53)) and

z¯¯εrz¯=z/ε1r,E55

which means that the solution should be proportional to expiε4+2r0x¯1α¯x¯1εdx¯1+β¯¯z¯¯t, where α¯x¯1is an O1function of x¯1that behaves like

α¯λβ¯¯Tw22x¯1+asx¯10E56

in this stage of evolution. The solution should remain inviscid in the main boundary layer and the viscous wall layer (i.e., a Stokes layer) is expected to be completely passive.

The scaled variable

y¯y/ε1rE57

will be O1in the main boundary layer since its thickness is now of the order of the spanwise length scale, Oε1r. It therefore follows from (53) and (57) that the transverse pressure gradients will come into play and the solution in this region should expand like

uvwp=U000+δ̂Ax¯1u¯y¯x¯1ε1rv¯y¯x¯1,ε1rw¯y¯x¯1ε21rp¯y¯x¯1expi1ε4+2r0x¯1α¯x¯1εdx¯1+β¯¯z¯¯tE58

where Ax¯1is a function of the slow variable x¯1. Substituting (58) into the linearized Navier-Stokes equations shows that the wall normal velocity perturbation v¯is determined by the incompressible reduced Rayleigh equation

Tddy¯1Tdv¯dy¯+Tα¯1α¯Uddy¯dU/dy¯Tβ¯¯2v¯=0E59

whose solution must satisfy the following boundary conditions

v¯eβ¯¯y¯asy¯,v¯=0aty¯=0.E60

Matching with the upstream solution (33) and (37) requires that α¯x¯1satisfy the matching condition (56) as x¯10.

Inserting (10) and (57) into (59), using (60) and assuming the ideal gas law ρT=1shows that

d1T2dv¯+α¯1α¯UUT2β¯¯2x¯12v¯=Oε21r,E61
v¯=0atη=0,E62

which means that

α¯=fβ̂,E63

where

β̂β¯¯2x¯1.E64

6.3 Matching with the triple-deck solution

Eq. (64) clearly approaches zero when x¯10,which means that α¯will be consistent with the matching condition (54) if we require that it behave like

α¯=λ/Tw2β̂+α1+α2β̂+asx¯10E65

where α1,α2are (in general complex) constants such that

α1=limx̂1κ¯1x̂1,α2=limx̂1κ¯2x̂1/β¯¯2x̂1.E66

Ref. [11] proved that (60)(64) possess an asymptotic solution of the form v¯=Uη+β̂v1+β̂2v2+..as β̂0when α¯satisfies (65) which implies that their solutions are able to match onto the lowest order triple-deck solution upstream and are consistent with the higher order solutions in this region.

6.4 Numerical results

The Rayleigh eigenvalues α¯are determined by the boundary value problem (60), (61) and (62). We assume in the following that the Prandtl number is equal to unity and that the viscosity μTsatisfies the simple linear relation μT=Tη.

Parts (a) and (b) of Figure 7 are plots of the real and imaginary parts respectively of these eigenvalues as a function of β̂. They show that the numerical solution for α¯will be consistent with the matching conditions (65)and(66) if the higher order terms in the triple-deck expansion(51) satisfy Imlimx̂1κ¯1x̂1=0and limx̂1κ¯2x̂1/β¯¯2x̂1=±iC, where the values of C are given in the caption of Figure 7. They also show that α¯is initially real and eventually becomes complex. But these eigenvalues must occur in complex conjugate pairs since the coefficients in (61) are all real. The computations show that Imα¯eventually goes to zero at some finite value of β̂which is consistent with the fact that U/T2is equal to zero at some finite value of ηand Eq. (61) therefore has a generalized inflection point there.

Figure 7.

(a) Reα¯ and (b) Imα¯ vs. β̂ calculated from the modified Rayleigh solution. The red dashed curves are calculated from the asymptotic formula(56) . The red dashed lines in the inset are Imα¯=Cβ̂, where C=36 for M∞=2C=129.4 for M∞=3 and C=340.1 for M∞=4.

7. Conclusions

This chapter uses high Reynolds number asymptotics to study the nonlocal behavior of boundary layer instabilities generated by small amplitude free-stream disturbances at subsonic and moderate supersonic Mach numbers. The appropriate small expansion parameter turns out to be ε=F1/6,where Fdenotes the frequency parameter. The oblique 1st Mack mode instabilities generated by free-stream acoustic disturbances are compared with those generated by elongated vortical disturbances. The focus is on explaining the relevant physics and not on obtaining accurate numerical predictions.

The free-stream vortical disturbances generate unsteady flows in the leading edge region that produce short spanwise wavelength instabilities in a viscous triple-deck region which lies at an Oε2distance downstream from the leading edge. The mechanism was first considered for two dimensional incompressible flows in Ref. [7], but the instability onset occurs much further upstream in the supersonic case and is, therefore, much more likely to be important at the higher Mach numbers considered in this chapter. The lowest order triple-deck solution does not possess an upper branch and evolves into an inviscid 1st Mack mode instability with short spanwise wavelength at an Oε4distance downstream.

Fedorov and Khokhlov [10] used asymptotic methods to study the generation of inviscid instabilities in supersonic boundary layers by fast and slow acoustic disturbances in the free stream whose obliqueness angle θdeviated from its critical value by an O1amount and showed that slow acoustic disturbances generate unsteady boundary layer disturbances that produce O1spanwise wavelength inviscid 1st Mack mode instabilities a much larger Oε6distance downstream. But the calculations in Ref. [11] show that the physical streamwise distance x=U3/ω2νcorresponding to this scaled downstream location is at least equal to about 7 m for the typical supersonic flight conditions at M=3U=888m/sν=0.000264m2/send an altitude of 20 km with an upper bound of 100 kHz for the characteristic frequency. This means that this instability occurs too far downstream to be of any practical interest at the moderately low supersonic Mach numbers considered in this chapter.

But, the inviscid instability, which first appears at an Oε4+2/3distance downstream when Δθis reduced to Oε2/3can be significant when scaled to flight conditions. It is therefore appropriate to compare the vortically-generated instabilities with the instabilities generated by oblique acoustic disturbances with obliqueness angles in this range as done in this chapter.

Acknowledgments

This research was sponsored by NASA’s by the Transformative Aeronautics Concepts Program of the Aeronautics Research Mission Directorate under the Transformational Tools and Technologies (TTT) Project. PR was supported by the Air Force Office of Scientific Research award number AFOSR Grant FA9550-15-1-0248. We would also like to thank Dr. Meelan Choudhari for bringing the photograph in Figure 2 to our attention.

Download

chapter PDF

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Marvin E. Goldstein and Pierre Ricco (January 30th 2019). Leading Edge Receptivity at Subsonic and Moderately Supersonic Mach Numbers [Online First], IntechOpen, DOI: 10.5772/intechopen.83672. Available from:

chapter statistics

138total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us