The generation of Tollmien–Schlichting waves by free-stream turbulence in transonic flow
Part of a collection:
Practical Asymptotics VII
A. I. Ruban ORCID: orcid.org/0000-0003-1028-88561 &
H. Broadley2
Journal of Engineering Mathematics volume 129, Article number: 1 (2021) Cite this article
This paper studies the generation of Tollmien–Schlichting waves by free-stream turbulence in transonic flow over a half-infinite flat plate with a roughness element using an asymptotic approach. It is assumed that the Reynolds number (denoted Re) is large, and that the free-stream turbulence is uniform so it can be modelled as vorticity waves. Close to the plate, a Blasius boundary layer forms at a thickness of \(O(\mathrm{{Re}}^{-{1}/{2}})\), and a vorticity deformation layer is also present with thickness \(O(\mathrm{{Re}}^{-{1}/{4}})\). The report shows that there is no mechanism by which the vorticity waves can penetrate from the vorticity deformation layer into the classical boundary layer; therefore, a transitional layer is introduced between them in order to prevent a discontinuity in vorticity. The flow in the interaction region in the vicinity of the roughness element is then analysed using the triple-deck model for transonic flow. A novel asymptotic expansion is used to analyse the upper deck, which enables a viscous–inviscid interaction problem to be derived. In order to make analytical progress, the height of the roughness element is assumed to be small, and from this, we find an explicit formula for the receptivity coefficient of the Tollmien–Schlichting wave far downstream of the roughness.
Laminar–turbulent transition is a phenomenon which has been studied extensively by researchers for over a century, beginning with the pioneering experiments of Reynolds in the late nineteenth century. It is an extremely complicated process, which begins when small perturbations from the free-stream flow penetrate into the boundary layer and interact with roughness on a body surface, producing Tollmien–Schlichting waves. Tollmien–Schlichting waves are instability modes that travel downstream of a roughness element; however, they must often travel far downstream of this roughness before they cause transition to turbulent flow. In quiet flows (such as those seen in aerodynamic applications), the amplitude of these waves is at first too small to cause transition from laminar to turbulent flow. However, when the Tollmien–Schlichting waves grow in amplitude in the boundary layer, they eventually produce nonlinear effects which are quickly followed by a transition to turbulent flow. The study of laminar–turbulent transition met a major obstacle when it was found that the transition would not necessarily take place at the same point in two experiments which used the same principal similarity parameters. For example, a model of an aircraft wing tested in different wind tunnels, with the same Mach number and Reynolds number of flow, would have different points at which the flow changed from laminar to turbulent.
Receptivity theory was introduced in order to account for this discrepancy, which is mostly due to the difference in perturbations present in the oncoming flow. Its aims are to describe the generation and propagation of Tollmien–Schlichting waves in the boundary layer, as well as to predict the point at which the amplitude becomes large enough to cause laminar–turbulent transition. In order to do this, it is necessary to identify the external perturbations which penetrate into the boundary layer and the external perturbations that do not. The objective of experimental works such as those performed by [1] was to study different types of perturbations and find which respective category each fell into. Examples of external perturbations that do penetrate into the boundary layer include acoustic waves; free-stream turbulence; and entropy waves. In this paper, Tollmien–Schlichting waves generated by a combination of free-stream turbulence and wall roughness are studied.
The effect of external acoustic noise on Tollmien–Schlichting wave generation in the boundary layer has been extensively studied for the past 40 years, beginning with [2, 3]. However, neither of these papers predicted that Tollmien–Schlichting waves could be effectively produced by acoustic waves. Their work was extended independently by [4, 5], who demonstrated that the strict resonance conditions, which have to be satisfied in order for perturbations to penetrate into the boundary layer and, therefore, generate Tollmien–Schlichting waves, of [6] can be applied to receptivity theory. In simple mechanical systems, resonance conditions are often easy to satisfy; for example, for a simple spring system, the only resonance condition is that the frequency of external forcing is the same as the natural frequency. For Tollmien–Schlichting waves, however, we require both the frequency and the wave number of external perturbations to match those of the natural oscillations in the boundary layer. References [4, 5] considered incompressible flow, and found that in order for both resonance conditions to be satisfied in the presence of acoustic waves, the frequency of the waves had to be of O(\(\mathrm{{Re}}^{-{1}/{4}}\)) and that wall roughness had to be introduced; otherwise, the wavelength of the perturbation was too large. Under these conditions, the amplitude of Tollmien–Schlichting waves was calculated. More recently, [7] studied the generation of Tollmien–Schlichting waves by acoustic noise in compressible transonic flows, using the triple-deck model of [8, 9].
The effect of free-stream turbulence on instabilities arising in the boundary layer has also been the subject of investigation. Reference [10] considered how uniform free-stream turbulence, taking the form of vorticity waves, could penetrate into the boundary layer in incompressible flow. The vorticity waves that form the uniform turbulence can also be thought of as convective gusts. These do not automatically satisfy the impermeability condition, which led Hunt and Graham to introduce a new 'vorticity deformation layer' of thickness O(\(\mathrm{{Re}}^{-{1}/{4}}\)). Further work by [11] analysed the effect of free-stream turbulence on the boundary layer, and concluded that vorticity waves do not penetrate far enough into the boundary layer to generate Tollmien–Schlichting waves. However, [12] demonstrated that this is not the case, despite the fact that pressure perturbations do not penetrate into the boundary layer. The authors introduced vorticity waves to incompressible flow over a half-infinite flat plate with roughness and solved the resulting viscous–inviscid interaction problem to find the amplitude of Tollmien–Schlichting waves produced.
The scope of the present study is to extend the analysis of [12] to transonic flow. Free-stream turbulence will be assumed to be uniform and, therefore, modelled by sinusoidal vorticity waves. Solutions to the Navier–Stokes equations will first be sought upstream of the roughness element where the main regions to be studied are the unperturbed flow region, the vorticity deformation region, and the classical boundary layer. Due to the lack of pressure perturbations outside of the boundary layer, it is necessary to introduce a transition layer in which vorticity waves decay in order to prevent a discontinuity in vorticity. Following this, the flow in the interaction region forming in the vicinity of the roughness will be studied. To accurately model the flow in the neighbourhood of the roughness, the triple-deck model with scalings given by [8, 9] will be used. The main difficulty is in the upper tier of the triple-deck structure, where the leading order perturbations equations appear to be degenerate, and in order to deduce the viscous-interaction problem for pressure higher-order terms have to be considered. To make analytic progress, we then assume that the roughness height is small, which is justified physically when considering the vertical extent of roughness elements found on aircraft wings. Solving the resulting linearised equations will enable an expression to be found for the amplitude of the Tollmien–Schlichting waves generated far downstream of the roughness element, at which point the objective of the current paper will have been met.
Problem formulation
Consider an ideal gas flow over a flat plate aligned with the free stream velocity. We assume that in the flow, there exist small amplitude vorticity waves (uniform free-stream turbulence). In addition, we shall have a roughness element on the flat plate a distant L from the leading edge of said plate. Placing the origin of our coordinate axes at the leading edge of the plate, it is natural to use (\(\hat{x}\), \(\hat{y}\)) as our coordinates. A schematic of the flow regime can be found in Fig. 1. We shall assume that the flow is two dimensional and that the body force is negligible. We define \((\hat{u}, \hat{v})\) as the velocity vector, \(\hat{p}\) the pressure, \(\hat{\rho }\) the density, \(\hat{h}\) the enthalpy, and \(\hat{\mu }\) the viscosity coefficient. We denote the free-stream velocity as \(U_\infty \), and the free-stream pressure, density, and viscosity coefficient as \(p_\infty \), \(\rho _\infty \), and \(\mu _\infty \), respectively. In order to non-dimensionalise the compressible Navier–Stokes equations, we introduce the following variables:
Schematic of flow regime
Substituting these new variables back into our dimensional Navier–Stokes equations gives us the following non-dimensional set of equations:
$$\begin{aligned}&\rho \left( \frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y}\right) =-\frac{\partial p}{\partial x}+\frac{1}{\mathrm{{Re}}}\left( \frac{\partial }{\partial x}\left[ \mu \left( \frac{4}{3} \frac{\partial u}{\partial x}-\frac{2}{3}\frac{\partial v}{\partial y}\right) \right] +\frac{\partial }{\partial y} \left[ \mu \left( \frac{\partial u}{\partial y}+\frac{\partial v}{\partial x}\right) \right] \right) , \end{aligned}$$
(2.2a)
$$\begin{aligned}&\rho \left( \frac{\partial v}{\partial t} + u \frac{\partial v}{\partial x} + v \frac{\partial v}{\partial y}\right) =-\frac{\partial p}{\partial y}+\frac{1}{\mathrm{{Re}}}\left( \frac{\partial }{\partial y}\left[ \mu \left( \frac{4}{3} \frac{\partial v}{\partial y}-\frac{2}{3}\frac{\partial u}{\partial x}\right) \right] +\frac{\partial }{\partial x}\left[ \mu \left( \frac{\partial u}{\partial y}+\frac{\partial v}{\partial x}\right) \right] \right) , \end{aligned}$$
(2.2b)
$$\begin{aligned}&\begin{aligned} \rho \left( \frac{\partial h}{\partial t} + u \frac{\partial h}{\partial x} + v \frac{\partial h}{\partial y}\right)&= \frac{\partial p}{\partial t} + u \frac{\partial p}{\partial x} + v \frac{\partial p}{\partial y}+\frac{1}{\mathrm{{Re}}} \Bigg [\frac{1}{\mathrm{{Pr}}}\left( \frac{\partial }{\partial x} \left( \mu \frac{\partial h}{\partial x}\right) +\frac{\partial }{\partial y}\left( \mu \frac{\partial h}{\partial y}\right) \right) \\&\quad + \mu \left( \frac{4}{3} \frac{\partial u}{\partial x}-\frac{2}{3}\frac{\partial v}{\partial y}\right) \frac{\partial u}{\partial x}+\mu \left( \frac{4}{3} \frac{\partial v}{\partial y}-\frac{2}{3}\frac{\partial u}{\partial x}\right) \frac{\partial v}{\partial y}+ \mu \left( \frac{ \partial u}{ \partial y}+\frac{ \partial v}{\partial x}\right) ^2 \Bigg ], \end{aligned} \end{aligned}$$
(2.2c)
$$\begin{aligned}&\frac{\partial \rho }{\partial t}+\frac{\partial (\rho u)}{\partial x}+\frac{\partial (\rho v)}{\partial y}=0, \end{aligned}$$
(2.2d)
$$\begin{aligned}&h=\frac{1}{(\gamma -1)M_{\infty }^2}\frac{1}{\rho }+\frac{\gamma }{\gamma -1}\frac{p}{\rho }. \end{aligned}$$
(2.2e)
Here, we have defined the Reynolds number (Re) and the free-stream Mach number (\(M_{\infty }\)) as follows:
$$\begin{aligned} \mathrm{{Re}}=\frac{\rho _{\infty } U_{\infty } L}{\mu _{\infty } }, \quad M_{\infty }=\frac{U_{\infty }}{a_{\infty }}, \end{aligned}$$
respectively, where \(a_{\infty }=\sqrt{(\gamma p_\infty )/\rho _\infty }\) is the speed of sound in the undisturbed flow upstream of the plate. In the analysis that follows, we will take the limit \(\mathrm{{Re}} \rightarrow \infty \) and keep the Mach number as a constant which is close to 1.
Flow upstream of the roughness element
Free-stream flow (region 1)
First, we shall consider the region far away from the plate (and upstream of the roughness). We look for solutions to the Navier–Stokes equations of the form
$$\begin{aligned} \begin{aligned} u&=1+\varepsilon u_1(x,y,t), \quad v=\varepsilon v_1(x,y,t), \quad \rho =1+\varepsilon \rho _1(x,y,t), \\ p&=\varepsilon p_1(x,y,t),\quad h=\frac{1}{(\gamma -1)M_{\infty }^2}+\varepsilon h_1(x,y,t), \end{aligned} \end{aligned}$$
where \(\varepsilon \) is a small parameter. Substituting these asymptotic expansions into the Navier–Stokes equations, we obtain
$$\begin{aligned}&\frac{\partial u_1}{\partial t} + \frac{\partial u_1}{\partial x} =-\frac{\partial p_1}{\partial x}, \end{aligned}$$
$$\begin{aligned}&\frac{\partial v_1}{\partial t} + \frac{\partial v_1}{\partial x}=-\frac{\partial p_1}{\partial y}, \end{aligned}$$
$$\begin{aligned}&\frac{\partial h_1}{\partial t} + \frac{\partial h_1}{\partial x} =\frac{\partial p_1}{\partial t} + \frac{\partial p_1}{\partial x}, \end{aligned}$$
$$\begin{aligned}&\frac{\partial \rho _1}{\partial t}+\frac{\partial \rho _1}{\partial x}+\frac{\partial u_1}{\partial x} + \frac{\partial v_1}{\partial y}=0, \end{aligned}$$
$$\begin{aligned}&h_1=-\frac{1}{(\gamma -1)M_{\infty }^2}\rho _1+\frac{\gamma }{\gamma -1}p_1. \end{aligned}$$
It is clear that we can immediately eliminate \(h_1\) from (3.2c) using (3.2e) giving
$$\begin{aligned} \frac{\partial \rho _1}{\partial t}+\frac{\partial \rho _1}{\partial x}=M_\infty ^2 \left( \frac{\partial p_1}{\partial t}+\frac{\partial p_1}{\partial x} \right) . \end{aligned}$$
We can look for wavelike solutions to this equation and (3.2a, 3.2b, 3.2d) of the form \(u_1(\xi ,y)\), \(v_1(\xi ,y)\), \(\rho _1(\xi ,y)\), \(p_1(\xi ,y)\) where \(\xi =x-t\). Consequently, we write
$$\begin{aligned} u_1=U\mathrm{{e}}^{\mathrm{{i}}(\alpha \xi + \beta y)},\quad v_1= V\mathrm{{e}}^{\mathrm{{i}}(\alpha \xi + \beta y)},\quad \rho _1=R\mathrm{{e}}^{\mathrm{{i}}(\alpha \xi +\beta y)},\quad p_1=P\mathrm{{e}}^{\mathrm{{i}}(\alpha \xi + \beta y)}, \end{aligned}$$
where \(\alpha \) and \(\beta \) are the wave numbers in the x and y directions, respectively, and U, V, P, and R are constants. Substituting these representations into our equations results in a eigenvalue problem for c. In order for us to obtain non-trivial solutions, we require either (i) \(U=V=P=0\), R arbitrary or (ii) \(P=R=0\). Vorticity waves correspond to the case where \(R=P=0\), so this is the case we shall be concerned with. We will consider the case of one harmonic, so we define vorticity as follows:
$$\begin{aligned} \omega =\frac{\partial v_1}{\partial x}-\frac{\partial u_1}{\partial y}=\omega _0\mathrm{{e}}^{\mathrm{{i}}(\alpha \xi +\beta y)} \end{aligned}$$
where \(\omega _0\) the amplitude of the wave. Substituting \(R=P=0\) into our governing equations, we see that the only equation not automatically satisfied is
$$\begin{aligned} \frac{\partial u_1}{\partial x} + \frac{\partial v_1}{\partial y}=0. \end{aligned}$$
Combining this with our definition of vorticity, and without loss of generality setting \(\omega _0=1\), we obtain the following solutions for the velocity components:
$$\begin{aligned} u_1&=\frac{\mathrm{{i}} \beta }{\alpha ^2 +\beta ^2}\mathrm{{e}}^{\mathrm{{i}}(\alpha \xi +\beta y)}, \end{aligned}$$
$$\begin{aligned} v_1&=-\frac{\mathrm{{i}} \alpha }{\alpha ^2 +\beta ^2}\mathrm{{e}}^{\mathrm{{i}}(\alpha \xi +\beta y)}. \end{aligned}$$
Of course due to Eq. (3.2e), we have that \(h_1\equiv 0\).
Vorticity deformation layer (region 2)
The work of [10] showed that introducing free-stream turbulence to flow over a flat plate causes a new layer of flow to be formed at \(y \sim O(\mathrm{{Re}}^{-{1}/{4}})\). With that in mind, we begin conducting analysis of this region by introducing the scaled variables:
$$\begin{aligned} \overline{x}=\mathrm{{Re}}^{-{1}/{4}}x, \quad \overline{y}=\mathrm{{Re}}^{-{1}/{4}}y, \quad \overline{t}=\mathrm{{Re}}^{-{1}/{4}}t. \end{aligned}$$
We look for perturbed solutions to the Euler equations of the form:
$$\begin{aligned} \begin{aligned}&u=1+\varepsilon \overline{u}_2(\overline{x},\overline{y},\overline{t}), \quad v= \varepsilon \overline{v}_2(\overline{x},\overline{y},\overline{t}), \quad \rho =1+\varepsilon \overline{\rho }_2(\overline{x},\overline{y},\overline{t}),\\&p=\varepsilon \overline{p}_2 (\overline{x},\overline{y},\overline{t}), \quad h=\frac{1}{(\gamma -1)M_{\infty }^2}+\varepsilon \overline{h}_2(\overline{x},\overline{y},\overline{t}). \end{aligned} \end{aligned}$$
Eliminating pressure we find that
$$\begin{aligned} \frac{\partial }{\partial \overline{t}} \left( \frac{\partial \overline{v}_2}{\partial \overline{x}}- \frac{\partial \overline{u}_2}{\partial \overline{y}} \right) +\frac{\partial }{\partial \overline{x}} \left( \frac{\partial \overline{v}_2}{\partial \overline{x}}- \frac{\partial \overline{u}_2}{\partial \overline{y}} \right) =0. \end{aligned}$$
Using this in conjunction with the definition of vorticity, we can see that vorticity does not change but is merely transported downstream with the mean free-stream velocity. As the vorticity upstream of the plate is given by (3.5), we can write
$$\begin{aligned} \frac{\partial \overline{v}_2}{\partial \overline{x}}-\frac{\partial \overline{u}_2}{\partial \overline{y}}=\mathrm{{e}}^{\mathrm{{i}}(\alpha \overline{\xi }+\beta \overline{y})}, \end{aligned}$$
where \(\overline{\xi }=\overline{x}-\overline{t}\). This equation suggests that \(\overline{u}_2\) and \(\overline{v}_2\) are functions of \(\overline{\xi }\) and \(\overline{y}\) only. Looking for solutions to the Euler equations of this form, we find
$$\begin{aligned} \frac{\partial \overline{p}_2}{\partial \overline{y}}=0. \end{aligned}$$
Taking into account that there are no pressure perturbations in the free stream, we see that \(\overline{p}_2=0\). Consequently, the energy equation becomes
$$\begin{aligned} \frac{\partial \overline{h}_2}{\partial \overline{t}} + \frac{\partial \overline{h}_2}{\partial \overline{x}} =0. \end{aligned}$$
In the previous section, we found that upstream of the leading edge of the plate \(h_1=0\). By the above equation, we must, therefore, have \(\overline{h}_2 \equiv 0\) in the entire vorticity deformation layer, from which we can deduce that \(\overline{\rho }_2\equiv 0\). The continuity equation is then reduced to
$$\begin{aligned} \frac{\partial \overline{u}_2}{\partial \overline{x}} + \frac{\partial \overline{v}_2}{\partial \overline{y}}=0. \end{aligned}$$
Using this, we can introduce a stream function \(\overline{\psi }_2\) which is found to be
$$\begin{aligned} \overline{\psi }_2=\mathrm{{e}}^{\mathrm{{i}}\alpha \overline{\xi }}(C\mathrm{{e}}^{\alpha \overline{y}} +D\mathrm{{e}}^{-\alpha \overline{y}})+\frac{1}{\alpha ^2+\beta ^2}\mathrm{{e}}^{\mathrm{{i}}(\alpha \overline{\xi }+ \beta \overline{y})}, \end{aligned}$$
where C and D are constants of integration. Therefore, we have
$$\begin{aligned} \overline{u}_2&=\alpha \mathrm{{e}}^{\mathrm{{i}}\alpha \overline{\xi }}\left( C\mathrm{{e}}^{\alpha \overline{y}} -D\mathrm{{e}}^{-\alpha \overline{y}}\right) +\frac{\mathrm{{i}} \beta }{\alpha ^2+\beta ^2}\mathrm{{e}}^{\mathrm{{i}}(\alpha \overline{\xi }+ \beta \overline{y})}, \end{aligned}$$
(3.16a)
$$\begin{aligned} \overline{v}_2&=- \mathrm{{i}}\alpha \mathrm{{e}}^{\mathrm{{i}}\alpha \overline{\xi }}\left( C\mathrm{{e}}^{\alpha \overline{y}} +D\mathrm{{e}}^{-\alpha \overline{y}}\right) -\frac{\mathrm{{i}} \alpha }{\alpha ^2+\beta ^2}\mathrm{{e}}^{\mathrm{{i}}(\alpha \overline{\xi }+ \beta \overline{y})}. \end{aligned}$$
(3.16b)
Clearly, we require \(C=0\) otherwise \(\overline{u}_2\) will grow exponentially with \(\overline{y}\). We can calculate D using the impermeability condition at the plate from which we can deduce that
$$\begin{aligned} D=-\frac{1}{\alpha ^2 + \beta ^2}. \end{aligned}$$
Steady boundary layer flow (region 3)
When the y coordinate becomes of order \((\mathrm{{Re}}^{-{1}/{2}})\), a viscous boundary layer forms close to the plate. As there are no leading order pressure perturbations in the vorticity deformation layer, free-stream turbulence does not penetrate into the boundary layer. Therefore, the flow in the boundary layer is steady. Assuming that that x remains an O(1) quantity, and that \(Y=\mathrm{{Re}}^{{1}/{2}}y\) is an O(1) quantity as \(\mathrm{{Re}} \rightarrow \infty \), we look for solutions to the compressible Navier–Stokes equations of the form:
$$\begin{aligned} \begin{aligned} u=U_0(x,Y)+\cdots , \quad v=\mathrm{{Re}}^{-\frac{1}{2}} V_0(x,Y)+\cdots , \quad \rho =\rho _0(x,Y)+\cdots ,&\quad \\p=\mathrm{{Re}}^{-\frac{1}{2}}P_1(x,Y)+\cdots , \quad h=h_0(x,Y)+\cdots , \quad \mu =\mu _0(x,Y)+\cdots . \end{aligned} \end{aligned}$$
Substituting these representations into the governing equations and letting \(\mathrm{{Re}} \rightarrow \infty \), we obtain (to leading order)
$$\begin{aligned}&\rho _0 U_0 \frac{\partial U_0}{\partial x}+\rho _0 V_0 \frac{\partial U_0}{\partial y}=\frac{\partial }{\partial Y}\left( \mu _0 \frac{\partial U_0}{\partial Y}\right) , \end{aligned}$$
$$\begin{aligned}&\rho _0 U_0 \frac{\partial h_0}{\partial x}+\rho _0 V_0 \frac{\partial h_0}{\partial y}=\frac{1}{\mathrm{{Pr}}}\frac{\partial }{\partial Y}\left( \mu _0 \frac{\partial h_0}{\partial Y}\right) +\mu _0\left( \frac{\partial U_0}{\partial Y}\right) ^2, \end{aligned}$$
$$\begin{aligned}&\frac{\partial }{\partial x}(\rho _0 U_0)+\frac{\partial }{\partial Y}(\rho _0 V_0)=0, \end{aligned}$$
(3.19c)
$$\begin{aligned}&h_0=\frac{1}{(\gamma -1)M_{\infty }^2}\frac{1}{\rho _0}. \end{aligned}$$
(3.19d)
The boundary conditions of the system are as follows. At the leading edge of the plate, we have
$$\begin{aligned} U_0=1, \quad h_0=\frac{1}{(\gamma -1)M_{\infty }^2} \quad \mathrm {at} \quad x=0, \ Y \in [0, \infty ). \end{aligned}$$
On the surface of the plate, we have the no-slip and impermeability conditions, which are expressed as follows:
$$\begin{aligned} U_0=V_0=0 \quad \mathrm {at} \quad Y=0, \ x \in [0,1]. \end{aligned}$$
The matching condition at the outer edge of the boundary is
$$\begin{aligned} U_0=1, \quad h_0=\frac{1}{(\gamma -1)M_{\infty }^2} \quad \mathrm {as} \quad Y \rightarrow \infty , \ x \in [0,1]. \end{aligned}$$
We will also choose a thermally insulated wall, which has a constant temperature. We write this condition as follows:
$$\begin{aligned} h_0=h_w \quad \mathrm {at} \quad Y=0, \ x \in [0,1]. \end{aligned}$$
Self similar solutions to this system of the variable \(\eta ={Y}/{\sqrt{x}}\) exist. We seek solutions of the form
$$\begin{aligned} \begin{aligned}&U_0(x,Y)=U(\eta ), \quad V_0(x.Y)=\frac{1}{\sqrt{x}}V(\eta ), \quad \rho _0(x,Y)=\rho (\eta ),\\&h_0(x,Y)=h(\eta ), \quad \mu _0(x,Y)= \mu (\eta ). \end{aligned} \end{aligned}$$
In the limit \(\eta \rightarrow \infty \), we have the solutions:
$$\begin{aligned} U=1-\frac{B}{2(\eta +A)}\mathrm{{e}}^{-\frac{1}{4}(\eta +A)^2} ,\quad \rho =1, \quad V=-\frac{1}{2}A,\quad h=\frac{1}{(\gamma -1)M_\infty ^2}+\frac{C}{2(\eta +A)}\mathrm{{e}}^{-\frac{\mathrm{{Pr}}}{4}(\eta +A)^2} . \end{aligned}$$
where A, B, and C are constants.
In the limit \(\eta \rightarrow 0\), we consider the flow in the neighbourhood of the roughness element and Taylor expands our solutions in this limit around \(x=1\). We find to leading order
$$\begin{aligned} u=\lambda Y, \quad V=0, \quad \rho = \rho _w, \quad h=h_w,\quad \mu =\mu _w, \end{aligned}$$
where \(\lambda \) is the usual Blasius constant.
Transition layer (region 4)
Outside of the boundary layer, we can see that there are no pressure perturbations to leading order. As a consequence, there is no mechanism by which our vorticity waves can penetrate into the boundary layer. However, we cannot have a discontinuity in vorticity at the outer edge of the boundary layer. This leads us to introduce a transition region between the classical boundary layer and the vorticity deformation layer, where the vorticity waves decay. We shall seek the solution to the Navier–Stokes equations in this region in the form:
$$\begin{aligned} u&=U_4(x,Y)+\varepsilon u_4(x, \overline{\xi }, \overline{Y})+ \cdots , \end{aligned}$$
$$\begin{aligned} v&=\mathrm{{Re}}^{-{1}/{2}}V_4(x, Y)+ \varepsilon v_4(x, \overline{\xi }, \overline{Y})+\cdots , \end{aligned}$$
$$\begin{aligned} h&=H_4(x,Y)+\sigma h_4(x, \overline{\xi }, \overline{Y})+\cdots , \end{aligned}$$
$$\begin{aligned} \rho&=R_4(x, Y)+\sigma \rho _4(x, \overline{\xi }, \overline{Y})+\cdots , \end{aligned}$$
$$\begin{aligned} \mu&=\mu (x,Y)+\sigma \mu _4(x, \overline{\xi }, \overline{Y})+\cdots , \end{aligned}$$
(3.27e)
where we assume that the perturbations of pressure are small enough to be neglected. The perturbations of u and v are assumed to be of the same order as those in the vorticity deformation layer, and the order of \(\sigma \) is to be found. The coordinates (x, Y) are the same as those in the boundary layer, and \(\overline{\xi }\) is as in the vorticity deformation layer. We define \(\overline{Y}\) through
$$\begin{aligned} y=\mathrm{{Re}}^{-{1}/{2}}\left[ \Delta (\mathrm{{Re}})\sqrt{x}-A\sqrt{x}+\delta (\mathrm{{Re}})\overline{Y}\right] , \end{aligned}$$
where we have that as \(\mathrm{{Re}} \rightarrow \infty \), \(\Delta (\mathrm{{Re}}) \rightarrow \infty , \ \delta (\mathrm{{Re}}) \rightarrow 0\). Rearranging for \(\overline{Y}\), we have
$$\begin{aligned} \overline{Y}=\frac{y\,\mathrm{{Re}}^{{1}/{2}}+A\sqrt{x}- \Delta (\mathrm{{Re}}) \sqrt{x}}{\delta (\mathrm{{Re}})}. \end{aligned}$$
Substituting the expansions into the Navier–Stokes equations, we find that
$$\begin{aligned} U_4-1=-\frac{B}{2\Delta }\mathrm{{e}}^{-\frac{\Delta ^2}{4}}\mathrm{{e}}^{-\frac{\Delta \delta \overline{Y}}{2 \sqrt{x}}}. \end{aligned}$$
If this is constant, then the periodic solution in \(\overline{\xi }\) would be a superposition of exponential functions to the power of \(\overline{Y}\). This would make it impossible to satisfy both boundary conditions for vorticity at each edge of the transition layer, so in order to prevent this, we set
$$\begin{aligned} \Delta \delta =2. \end{aligned}$$
In what follows we will take our leading order expansions \(U_4\), \(V_4\), \(\rho _4\), and \(\mu _4\) from the outer edge of the Blasius boundary layer studied in the previous section. To avoid degeneration of the first term in the energy equation, we will set
$$\begin{aligned} \mathrm{{Re}}^{\frac{1}{4}}\frac{B}{2 \Delta } \mathrm{{e}}^{-\frac{\Delta ^2}{4}}=\frac{1}{\delta ^2}. \end{aligned}$$
This reduces our equations to
$$\begin{aligned}&-\frac{\sigma }{\delta ^2}\mathrm{{e}}^{-\frac{\overline{Y}}{\sqrt{x}}} \frac{\partial \rho _4}{\partial \overline{\xi }}+ \frac{\varepsilon }{\delta } \mathrm{{Re}}^{\frac{1}{2}} \frac{\partial v_4}{\partial \overline{Y}}- \frac{\varepsilon }{\delta ^2 \sqrt{x}} \frac{\partial u_4}{\partial \overline{Y}}- \frac{\sigma }{\delta ^2 \sqrt{x}} \frac{\partial \rho _4}{\partial \overline{Y}}=0, \end{aligned}$$
$$\begin{aligned}&-\frac{\sigma }{\delta ^2}\mathrm{{e}}^{-\frac{\overline{Y}}{\sqrt{x}}} \frac{\partial h_4}{\partial \overline{\xi }}+ \varepsilon v_4\, \mathrm{{Re}}^{\frac{1}{2}} \frac{\partial h_4}{\partial Y}- \frac{\sigma }{\delta ^2 \sqrt{x}} \frac{\partial h_4}{\partial \overline{Y}}= \frac{1}{\mathrm{{Pr}}} \frac{\sigma }{\delta ^2} \frac{\partial ^2 h_4}{\partial \overline{Y}^2}+ \frac{\varepsilon ^2}{\delta ^2} \left[ \left( \frac{\partial u_4}{\partial \overline{Y}}\right) ^2 +\frac{4}{3} \left( \frac{\partial v_4}{\partial \overline{Y}}\right) ^2 \right] , \end{aligned}$$
$$\begin{aligned}&-\frac{\varepsilon }{\delta ^2} \mathrm{{e}}^{-\frac{\overline{Y}}{\sqrt{x}}} \frac{\partial u_4}{\partial \overline{\xi }} -\frac{\varepsilon }{\delta ^2 \sqrt{x}} \frac{\partial u_4}{\partial \overline{Y}}+ \varepsilon v_4\, \mathrm{{Re}}^{\frac{1}{2}} \frac{\partial U_4}{\partial \overline{Y}}= \frac{\varepsilon }{\delta ^2} \frac{\partial ^2 u_4}{\partial \overline{Y}^2}, \end{aligned}$$
$$\begin{aligned}&h_4=-\frac{\rho _4}{(\gamma -1)M_{\infty }^2}. \end{aligned}$$
By the principle of least degeneration, we want to keep as many of the terms in these equations as possible. To this end, we set
$$\begin{aligned} \delta =\mathrm{{Re}}^{-\frac{1}{2}}, \end{aligned}$$
which is consistent with \(\delta \rightarrow 0\) as \(\mathrm{{Re}} \rightarrow \infty \), as well setting \(\sigma =\varepsilon \). Assuming that the solutions are periodic in \(\overline{\xi }\) and matching with the perturbations in both the Blasius boundary layer and vorticity deformation layer, we find
$$\begin{aligned} u_4(\overline{\xi }, \overline{Y}; x)=\mathrm{{e}}^{\mathrm{{i}} \alpha \overline{\xi }} \frac{1}{\alpha - \mathrm{{i}} \beta } J_0 \left( 2 \sqrt{\alpha x \mathrm{{i}}}\ \mathrm{{e}}^{-\frac{\overline{Y}}{2 \sqrt{x}}}\right) ,\quad h_4 \equiv 0,\quad \rho _4 \equiv 0, \end{aligned}$$
where \(J_0\) is the zeroth-order Bessel function of the first kind.
Flow in the interaction region
In order to analyse the flow around the roughness element, we shall use the same triple-deck structure as [8, 9]. We require a time scale of \(t=O(\mathrm{{Re}}^{-{2}/{9}})\) in order for the resonance conditions with Tollmien–Schlichting waves to be satisfied, and we shall further assume that the Mach number is of the form:
$$\begin{aligned} M_\infty =1+\mathrm{{Re}}^{-{1}/{9}}Q_\infty \end{aligned}$$
where \(Q_\infty \) is the Kármán–Guderley parameter.
Lower deck (region 5)
The horizontal extent of the roughness element is assumed to be of the order \(x=O(\mathrm{{Re}}^{-{1}/{3}})\), and the vertical extent of the lower deck is \(y=O(\mathrm{{Re}}^{-{11}/{18}})\). We shall, therefore, assume that the shape of the roughness is of the form:
$$\begin{aligned} y=\mathrm{{Re}}^{-{11}/{18}}F\left( \frac{x-1}{\mathrm{{Re}}^{-{1}/{3}}}\right) . \end{aligned}$$
It is natural to introduce re-scaled coordinates (\(x_{*}\),\(y_5\)) and time \(t_*\) as follows:
$$\begin{aligned} x=1+\mathrm{{Re}}^{-{1}/{3}}x_*, \quad y=\mathrm{{Re}}^{-{11}/{18}}y_5, \quad t=\mathrm{{Re}}^{-{2}/{9}}t_*. \end{aligned}$$
We seek solutions to the Navier–Stokes equations in this region in the form:
$$\begin{aligned} u&=\mathrm{{Re}}^{-\frac{1}{9}}u_5(x_*,y_5)+\cdots + \varepsilon \tilde{u}_5(x_*,y_5, t_*)+\cdots , \end{aligned}$$
$$\begin{aligned} v&=\mathrm{{Re}}^{-\frac{7}{18}}v_5( x_*,y_5)+\cdots + \varepsilon \, \mathrm{{Re}}^{-\frac{5}{18}}\tilde{v}_5(x_*,y_5, t_*)+\cdots , \end{aligned}$$
$$\begin{aligned} p&=\mathrm{{Re}}^{-\frac{2}{9}}p_5(x_*,y_5)+ \cdots + \varepsilon \, \mathrm{{Re}}^{-\frac{1}{9}}\tilde{p}_5(x_*,y_5, t_*)+\cdots , \end{aligned}$$
$$\begin{aligned} h&=h_w+\cdots , \end{aligned}$$
$$\begin{aligned} \rho&=\rho _w+ \cdots , \end{aligned}$$
$$\begin{aligned} \mu&=\mu _w+ \cdots , \end{aligned}$$
(4.4f)
where we assume that close to the roughness element enthalpy, viscosity, and density are constants, and \(\varepsilon<< \mathrm{{Re}}^{-{1}/{9}}\) to prevent the expansion becoming disordered. The first terms in these expansions correspond to the steady flow solution, whereas those that are O(\(\varepsilon \)) are due to time-dependent perturbations. Substituting these expansions into the non-dimensional Navier–Stokes equations yields to O(1):
$$\begin{aligned}&\rho _w\left( u_5 \frac{\partial u_5}{\partial x_*}+ v_5\frac{\partial u_5}{\partial y_5}\right) = -\frac{\partial p_5}{\partial x_*}+\mu _w \frac{\partial ^2 u_5}{\partial y_5^2}, \end{aligned}$$
$$\begin{aligned}&\frac{\partial p_5}{\partial y_5}=0, \end{aligned}$$
$$\begin{aligned}&\frac{\partial u_5}{\partial x_*}+ \frac{\partial v_5}{\partial y_5}=0. \end{aligned}$$
Our boundary conditions for these equations come from our no-slip conditions on the roughness element, and matching u with flow upstream of the roughness as \(x_* \rightarrow - \infty \). These are
$$\begin{aligned}&u_5=v_5=0 \quad \mathrm {at} \quad y_5=F(x_*), \end{aligned}$$
$$\begin{aligned}&u_5=\lambda y_5 \quad \mathrm {as} \quad x_* \rightarrow -\infty . \end{aligned}$$
Clearly, we can see that the pressure \(p_5\) does not change across the lower deck with respect to the y-direction. Using the continuity equation we find that in the limit \(y_5\rightarrow \infty \)
$$\begin{aligned} u_5=\lambda y_5+A_*,\quad v_5=-\frac{\mathrm{{d}}A_*}{\mathrm{{d}}x_*}y_5, \end{aligned}$$
where the function \(A_*\) is still to be found. We can perform a similar analysis on the O(\(\varepsilon \)) terms in our expansions. The boundary conditions for this new problem are
$$\begin{aligned}&\tilde{u}_5=\tilde{v}_5=0 \quad \mathrm {at} \quad y_5=F(x_*), \end{aligned}$$
$$\begin{aligned}&\tilde{u}_5=0 \quad \mathrm {as} \quad x_* \rightarrow - \infty . \end{aligned}$$
Substituting into the Navier–Stokes equations, we obtain the O(\(\varepsilon \)) equations:
$$\begin{aligned}&\rho _w \left( \frac{\partial \tilde{u}_5}{\partial t_*}+ u_5\frac{\partial \tilde{u}_5}{\partial x_*}+ \tilde{u}_5 \frac{\partial u_5}{\partial x_*}+ v_5 \frac{\partial \tilde{u}_5}{\partial y_5}+ \tilde{v}_5 \frac{\partial u_5}{\partial y_5}\right) =- \frac{\partial \tilde{p}_5}{\partial x_*}+\mu _w \frac{\partial ^2 \tilde{u}_5}{\partial y_5^2}, \end{aligned}$$
$$\begin{aligned}&\frac{\partial \tilde{p}_5}{\partial y_5}=0, \end{aligned}$$
$$\begin{aligned}&\frac{\partial \tilde{u}_5}{\partial x_*}+\frac{\partial \tilde{v}_5}{\partial y_5}=0. \end{aligned}$$
In this case, we find that as \(y_5 \rightarrow \infty \)
$$\begin{aligned} \tilde{u}_5=\tilde{A}_*(x_*, t_*), \quad \tilde{v}_5=-\frac{\partial \tilde{A}_*}{\partial x_*}y_5. \end{aligned}$$
In order to understand, the physical meanings of \(A_*\) and \(\tilde{A}_*\) consider the streamline slope angle \(\theta \) at the outer edge of region 5. Taking the limit as \(y_5 \rightarrow \infty \), we have
$$\begin{aligned} \theta =\mathrm {tan}^{-1}\frac{v}{u}=\mathrm {tan}^{-1}\left( -\frac{\mathrm{{Re}}^{-\frac{7}{18}}\frac{\mathrm{{d}} A_*}{\mathrm{{d}}x_*}y_5+ \varepsilon \,\mathrm{{Re}}^{-\frac{5}{18}}\frac{\partial \tilde{A}_*}{\partial x_*}y_5}{\mathrm{{Re}}^{-\frac{1}{9}}(\lambda y_5+A_*)+\varepsilon \tilde{A}_*} \right) \approx -\frac{1}{\lambda }\left( \mathrm{{Re}}^{-\frac{5}{18}} \frac{\mathrm{{d}}A_*}{\mathrm{{d}}x_*}+ \varepsilon \,\mathrm{{Re}}^{-\frac{1}{6}} \frac{\partial \tilde{A}_*}{\partial x_*} \right) . \end{aligned}$$
Due to this relationship between \(\theta \), \(A_*\), and \(\tilde{A}_*\) we call \(A_*\) the O(1) displacement function and \(\tilde{A}_*\) the O(\(\varepsilon \)) displacement function.
Middle deck (region 6)
The middle deck region is a continuation of the boundary layer upstream of the roughness. As a consequence, it has vertical extent \(y=O(\mathrm{{Re}}^{-{1}/{2}})\) which leads us to re-scale our y-coordinate as follows:
$$\begin{aligned} y=\mathrm{{Re}}^{-\frac{1}{2}}y_6. \end{aligned}$$
Relating this to the solution for the lower deck by considering the matching conditions for u and v as \(y_6 \rightarrow 0\), we can see that in this limit the solution in the middle deck is
$$\begin{aligned} u&= \lambda y_6 + \mathrm{{Re}}^{-\frac{1}{9}}A_*(x_*)+\cdots + \varepsilon \tilde{A}_*(x_*, t_*)+\cdots , \end{aligned}$$
$$\begin{aligned} v&=\mathrm{{Re}}^{-\frac{5}{18}} y_6 \frac{\mathrm{{d}}A_*}{\mathrm{{d}}x_*}+\cdots + \varepsilon \, \mathrm{{Re}}^{-\frac{1}{6}} y_6 \frac{\partial \tilde{A}_*}{\partial x_*}+\cdots \ . \end{aligned}$$
Consequently, we shall seek solutions in this region of the form:
$$\begin{aligned} u&=u_{6}(y_6)+\mathrm{{Re}}^{-\frac{1}{9}}u_{6}^*(x_*, y_6)+\cdots + \varepsilon \tilde{u}_6(x_*, y_6, t_*)+\cdots , \end{aligned}$$
$$\begin{aligned} v&=\mathrm{{Re}}^{-\frac{5}{18}}v_6(x_*, y_6)+\cdots \varepsilon \, \mathrm{{Re}}^{-\frac{1}{6}} \tilde{v}_6(x_*, y_6, t_*)+\cdots , \end{aligned}$$
$$\begin{aligned} p&=\mathrm{{Re}}^{-\frac{2}{9}}p_6(x_*, y_6)+\cdots \varepsilon \, \mathrm{{Re}}^{-\frac{1}{9}} \tilde{p}_6(x_*, y_6, t_*)+\cdots , \end{aligned}$$
$$\begin{aligned} \rho&= \rho _{60}(y_6)+\mathrm{{Re}}^{-\frac{1}{9}} \rho _{61}(x_*, y_6)+\cdots + \varepsilon \tilde{\rho }_6(x_*, y_6, t_*)+\cdots , \end{aligned}$$
$$\begin{aligned} h&=h_{60}(y_6)+\mathrm{{Re}}^{-\frac{1}{9}} h_{61}(x_*, y_6)+\cdots + \varepsilon \tilde{h}_6(x_*, y_6, t_*)+\cdots , \end{aligned}$$
$$\begin{aligned} \mu&=\mu _{60}(y_6)+\mathrm{{Re}}^{-\frac{1}{9}}\mu _{61}(x_*, y_6)+\cdots + \varepsilon \tilde{\mu }_6(x_*, y_6, t_*)+\cdots \ . \end{aligned}$$
(4.14f)
Substituting these expansions into the Navier–Stokes equations and using elimination to find u and v, we find the streamline slope angle to be
$$\begin{aligned} \theta =\mathrm{{Re}}^{-\frac{5}{18}}K(x_*)+ \varepsilon \, \mathrm{{Re}}^{-\frac{1}{6}} \tilde{K}(x_*, t_*), \end{aligned}$$
where \(K(x_*)\), \(\tilde{K}(x_*, t_*)\) are unknown functions. As \(\theta \) has no dependence on \(y_6\) we can match this with streamline slope angle from the lower deck yielding
$$\begin{aligned} \theta =-\frac{1}{\lambda }\left( \mathrm{{Re}}^{-\frac{5}{18}} \frac{\mathrm{{d}}A_*}{\mathrm{{d}} x_*} +\varepsilon \, \mathrm{{Re}}^{-\frac{1}{6}} \frac{\partial \tilde{A}_*}{\partial x_*}\right) . \end{aligned}$$
Upper deck (region 7)
In order to analyse the upper deck, we are required to re-scale y again. As we have that the vertical extent of the upper deck is \(y=O(\mathrm{{Re}}^{-{5}/{18}})\), we shall introduce a new variable \(y_7\) defined as follows:
$$\begin{aligned} y=\mathrm{{Re}}^{-\frac{5}{18}}y_7. \end{aligned}$$
We then look for solutions to the Navier–Stokes equations in this region of the form:
$$\begin{aligned} u&=1+\mathrm{{Re}}^{-{2}/{9}}u_{st}(x_*, y_7)\cdots +\varepsilon u_\alpha (t_*)+\cdots +\varepsilon \mathrm{{Re}}^{-{1}/{9}}u_*(t_*, x_*, y_7)+ \varepsilon \,\mathrm{{Re}}^{-{2}/{9}} \tilde{u_*}(t_*, x_*, y_7)+\cdots , \end{aligned}$$
$$\begin{aligned} v&=\mathrm{{Re}}^{-{5}/{18}}v_{st}(x_*, y_7)+\cdots +\varepsilon \mathrm{{Re}} ^{-{1}/{6}}v_*(t_*, x_*, y_7)+ \varepsilon \, \mathrm{{Re}}^{-{5}/{18}} \tilde{v_*}(t_*, x_*, y_7)+\cdots , \end{aligned}$$
$$\begin{aligned} p&=\mathrm{{Re}}^{-{2}/{9}}p_{st}(x_*, y_7)+\cdots +\varepsilon Re^{-{1}/{9}}p_* (t_*, x_*, y_7)+ \varepsilon \, \mathrm{{Re}}^{-{2}/{9}} \tilde{p_*}(t_*, x_*, y_7)+\cdots , \end{aligned}$$
$$\begin{aligned} \rho&= 1+\mathrm{{Re}}^{-{2}/{9}} \rho _{st}(x_*, y_7)+\cdots + \varepsilon \mathrm{{Re}}^{-{1}/{9}} \rho _*(t_*, x_*, y_7)+ \varepsilon \mathrm{{Re}}^{-{2}/{9}} \tilde{\rho _*}(t_*, x_*, y_7)+\cdots , \end{aligned}$$
$$\begin{aligned} h&=\frac{1}{(\gamma -1)M_\infty ^2}+\mathrm{{Re}}^{-{2}/{9}}h_{st}(x_*, y_7)+ \cdots + \varepsilon \mathrm{{Re}}^{-{1}/{9}} h_*(t_*, x_*, y_7) + \varepsilon \mathrm{{Re}}^{-{2}/{9}} \tilde{h_*}(t_*, x_*, y_*)+\cdots , \end{aligned}$$
where \(u_\alpha \) is a time-dependent perturbation due to vorticity which is known to be of the form \(e^{i \alpha t_*}\), and \(u_{st}\), \(v_{st}\), \(p_{st}\), \(\rho _{st}\), and \(h_{st}\) are the solutions in the steady flow. Derivatives of \(u_\alpha \) will be cancelled out by next order expansions of the perturbation terms, which have been omitted in the asymptotic expansion of u. Substituting these expansions into the Navier–Stokes equations, we obtain to O(1):
$$\begin{aligned}&\frac{\partial u_{st}}{\partial x_*}=-\frac{\partial p_{st}}{\partial x_*}, \end{aligned}$$
$$\begin{aligned}&\frac{\partial v_{st}}{\partial x_*}=-\frac{\partial p_{st}}{\partial y_7}, \end{aligned}$$
$$\begin{aligned}&\frac{\partial \rho _{st}}{\partial x_*}=\frac{\partial p_{st}}{\partial x_*}, \end{aligned}$$
$$\begin{aligned}&\frac{\partial u_{st}}{\partial x_*}+\frac{\partial \rho _{st}}{\partial x_*}=0. \end{aligned}$$
Using (4.19c) in (4.19d) and comparing this to (4.19a), we can see that this system is degenerate. It is, therefore, necessary to retain the next order terms, from which we obtain the equation
$$\begin{aligned} 2 Q_\infty \frac{\partial p_{st}}{\partial x_*} +\frac{\partial v_{st}}{\partial y_7}=0. \end{aligned}$$
We, therefore, have that the governing equation for \(p_{st}\) is given by
$$\begin{aligned} 2Q_\infty \frac{\partial ^2 p_{st}}{\partial x_*^2}-\frac{\partial ^2 p_{st}}{\partial y_7^2}=0. \end{aligned}$$
This is accompanied by the boundary condition that the perturbations are expected to decay far away from the roughness
$$\begin{aligned} u_{st}, v_{st}, p_{st}, \rho _{st} \rightarrow 0 \quad \mathrm {as} \quad x_*^2+y_7^2 \rightarrow \infty . \end{aligned}$$
We now turn to the unsteady terms in our expansions (4.18). Our O(\(\varepsilon \)) equations are
$$\begin{aligned}&\frac{\partial u_*}{\partial x_*} =-\frac{\partial p_*}{\partial x_*}, \end{aligned}$$
$$\begin{aligned}&\frac{\partial v_*}{\partial x_*}=- \frac{\partial p_*}{\partial y_7}, \end{aligned}$$
$$\begin{aligned}&\frac{\partial p_*}{\partial x_*} - \frac{\partial \rho _*}{\partial x_*} = 0, \end{aligned}$$
$$\begin{aligned}&\frac{\partial u_*}{\partial x_*}+ \frac{\partial \rho _*}{\partial x_*} =0. \end{aligned}$$
This system is again degenerate, so it is necessary to consider the next order equations which are
$$\begin{aligned}&u_\alpha \frac{\partial u_{st}}{\partial x_*}+\frac{\partial u_*}{\partial t_*}+ \frac{\partial \tilde{u}_*}{\partial x_*} =- \frac{\partial \tilde{p}_*}{\partial x_*}, \end{aligned}$$
$$\begin{aligned}&u_\alpha \frac{\partial v_{st}}{\partial x_*}+ \frac{\partial v_*}{\partial t_*}+\frac{\partial \tilde{v_*}}{\partial x_*} = -\frac{\partial \tilde{p_*}}{\partial y_7}, \end{aligned}$$
$$\begin{aligned}&u_\alpha \frac{\partial p_{st}}{\partial x_*}- u_\alpha \frac{\partial \rho _{st}}{\partial x_*}+ \frac{\partial p_*}{\partial t_*} - \frac{\partial \rho _*}{\partial t_*} +\frac{\partial \tilde{p}_*}{\partial x_*} - \frac{\partial \tilde{\rho }_*}{\partial x_*}+2Q_\infty \frac{\partial \rho _*}{\partial x_*}+2Q_\infty u_\alpha \frac{\partial \rho _{st}}{\partial x_*} = 0, \end{aligned}$$
$$\begin{aligned}&u_\alpha \frac{\partial \rho _{st}}{\partial x_*}+\frac{\partial \rho _*}{\partial t_*}+ \frac{\partial \tilde{\rho _*}}{\partial x_*}+ \frac{\partial \tilde{u_*}}{\partial x_*}+ \frac{\partial v_*}{\partial y_7} =0. \end{aligned}$$
Making use of Eqs. (4.19) and (4.23), this problem can be reduced to the following equation for the unsteady pressure perturbation in the upper tier
$$\begin{aligned} 2Q_\infty \frac{\partial ^2 p_*}{\partial x_*^2} + 2\frac{\partial ^2 p_*}{\partial t_* \partial x_*}-\frac{\partial ^2 p_*}{\partial y_7^2} =-2 u_\alpha \frac{\partial ^2 p_{st}}{\partial x_*^2} \end{aligned}$$
and of course, we expect a similar decay of these perturbations as in the O(1) case.
In order to close our viscous-inviscid interaction problem for \(p_*\), we require an additional boundary condition at \(y_7=0\) for both steady and unsteady flow. To do this, we consider our streamline slope angle \(\theta \) for region 6. So far, we have not considered the continuation of our transition layer into the triple-deck model. However, using a similar argument to [12], we can see that this is not necessary. The displacement thickness of the boundary layer has only been affected by the viscous lower deck, where the velocity was of small magnitude in comparison with other regions. We expect the velocity to be continuous across the transition layer of the triple deck, so it would be natural for it to be of the same order as the velocity in regions 6 and 7. Therefore, we can reasonably assume that the streamline slope angle will not be affected by this layer, as it is unaffected by region 6 which has the same magnitude of velocity and is much thicker. Consequently, the streamline slope angle at the bottom of region 7 is the same as it is in region 6 so we have
$$\begin{aligned} \theta =-\mathrm{{Re}}^{-\frac{5}{18}}\frac{1}{\lambda }\left( \mathrm{{Re}}^{-\frac{5}{18}} \frac{\mathrm{{d}}A_*}{\mathrm{{d}}x_*} +\varepsilon \, \mathrm{{Re}}^{-\frac{1}{6}} \frac{\partial \tilde{A}_*}{\partial x_*}\right) \quad \mathrm {at} \quad y_7=0. \end{aligned}$$
From this, we can derive the boundary condition on \(p_{st}\) and \(p_*\) at \(y_7=0\) which are
$$\begin{aligned}&\frac{\partial p_{st}}{\partial y_7}=\frac{1}{\lambda }\frac{\mathrm{{d}}^2 A_* }{\partial x_*^2} , \end{aligned}$$
$$\begin{aligned}&\mathrm {and} \quad \frac{\partial p_*}{\partial y_7}= \frac{1}{\lambda }\frac{\partial ^2 \tilde{A}_*}{\partial x_*^2}\quad \mathrm {at} \quad y_7=0, \end{aligned}$$
respectively. For steady flow, we, therefore, have
$$\begin{aligned}&2Q_\infty \frac{\partial ^2 p_{st}}{\partial x_*^2}-\frac{\partial ^2 p_{st}}{\partial y_7^2}=0, \end{aligned}$$
$$\begin{aligned}&p_{st}\rightarrow 0 \quad \mathrm {as} \quad x_*^2+y_7^2 \rightarrow \infty , \end{aligned}$$
$$\begin{aligned}&\frac{\partial p_{st}}{\partial y_7}=\frac{1}{\lambda }\frac{\mathrm{{d}}^2 A_* }{\mathrm{{d}}x_*^2} \quad \mathrm {at} \quad y_7=0. \end{aligned}$$
For unsteady perturbations, we have
$$\begin{aligned}&2Q_\infty \frac{\partial ^2 p_*}{\partial x_*^2} + 2\frac{\partial ^2 p_*}{\partial t_* \partial x_*}-\frac{\partial ^2 p_*}{\partial y_7^2} =-2 u_\alpha \frac{\partial ^2 p_{st}}{\partial x_*^2}., \end{aligned}$$
$$\begin{aligned}&p_* \rightarrow 0 \quad \mathrm {as} \quad x_*^2+y_7^2 \rightarrow \infty , \end{aligned}$$
$$\begin{aligned}&\frac{\partial p_*}{\partial y_7}= \frac{1}{\lambda } \frac{\partial ^2 \tilde{A}_*}{\partial x_*^2}\quad \mathrm {at} \quad y_7=0. \end{aligned}$$
Linear receptivity
In order to make analytic progress with our systems of equations (both steady and non-steady), we assume that the height of the roughness is small in comparison to the lower deck. We, therefore, write the equation of the roughness to be
$$\begin{aligned} y=\delta f(x) \quad \mathrm {where} \quad \delta \ll 1. \end{aligned}$$
Using this assumption, we can linearise our equations for the lower deck. In this region, we, therefore, write the expansions for our velocities terms and pressure as follows:
$$\begin{aligned} u&=\lambda y_5+\delta u_5'(x_*, y_5)+ \cdots + \delta \sigma \tilde{u}_5'(x_*, y_5, t_*)+\cdots , \end{aligned}$$
$$\begin{aligned} v&= \delta v_5'(x_*, y_5)+ \cdots + \delta \sigma \tilde{v}_5'(x_*, y_5, t_*)+ \cdots , \end{aligned}$$
$$\begin{aligned} p&= h P_5'(x_*, y_5)+ \cdots +\delta \sigma \tilde{P}_5'(x_*, y_5, t_*)+ \cdots , \end{aligned}$$
where \(\sigma \) is a small parameter. We shall also perform similar expansions for the displacement function A, and upper deck pressure p of the form:
$$\begin{aligned} A&=\delta A_*'(x_*)+\cdots +\delta \sigma \tilde{A}_*'(x_*, t_*)+ \cdots , \end{aligned}$$
$$\begin{aligned} p&=\delta p_{st}'(x_*, y_7)+\cdots +\delta \sigma p_*'(x_*, t_*)+ \cdots . \end{aligned}$$
We now substitute these representations into our equations for the lower deck. First, we shall consider the problem for steady flow.
Steady problem
Our linearised equations for steady flow to \(O(\delta )\) read
$$\begin{aligned}&\rho _w \left( \lambda y_5 \frac{\partial u_5'}{\partial x_*} + \lambda v_5' \right) =-\frac{\partial P_5'}{\partial x_*}+ \mu _w \frac{\partial ^2 u_5'}{\partial y_5^2}, \end{aligned}$$
$$\begin{aligned}&\frac{\partial P_5'}{\partial y_5}=0, \end{aligned}$$
$$\begin{aligned}&\frac{\partial u_5'}{\partial x_*}+\frac{\partial v_5'}{\partial y_5}=0. \end{aligned}$$
We also have the boundary conditions:
$$\begin{aligned} u_5'&=-\lambda f(x_*) \quad \mathrm {at} \quad y_5=0, \end{aligned}$$
$$\begin{aligned} v_5'&=0 \quad \mathrm {at} \quad y_5=0, \end{aligned}$$
$$\begin{aligned} u_5'&=A_*' \quad \mathrm {as} \quad y_5\rightarrow \infty , \end{aligned}$$
$$\begin{aligned} u_5'&=0 \quad \mathrm {as} \quad x_* \rightarrow - \infty . \end{aligned}$$
We also need to relate \(P_5'\) to \(p_{st}\). As the pressure distribution does not change across the middle deck, we can relate \(P_5'\) back to the pressure in the lower deck by the equation:
$$\begin{aligned} p_{st}|_{y_7=0}=P_5'. \end{aligned}$$
This system of equations will be solved in the usual way of transforming our variables into Fourier space. We define the Fourier transforms as follows:
$$\begin{aligned} (\overline{P}_5, \overline{p}_{st})=\int _{-\infty }^\infty \mathrm{{e}}^{-\mathrm{{i}}kx_*}(P_5', p_{st} )\mathrm{{d}}x_*. \end{aligned}$$
Defining a new variable \(z=(\mathrm{{i}} k)^{\frac{1}{3}} y_5\), we find the solution for subsonic flow (\(Q_\infty <0\)) to be
$$\begin{aligned} \overline{P}_5&=\frac{ |k|}{ \sqrt{2 |Q_\infty |}}\frac{ \overline{f}(k)q^{\frac{2}{3}} Ai'(0)}{q^{\frac{2}{3}} Ai'(0)- \frac{(i k)^{\frac{1}{3}} |k|}{\lambda \mu _w \sqrt{2|Q_\infty |}} \int _0^\infty Ai(q^{\frac{2}{3}} z') \mathrm{{d}}z' }, \end{aligned}$$
$$\begin{aligned} \overline{p}_{st}&=\frac{ |k|}{ \sqrt{2 |Q_\infty |}}\frac{ \overline{f}(k)q^{\frac{2}{3}} Ai'(0)}{q^{\frac{2}{3}} Ai'(0)- \frac{(i k)^{\frac{1}{3}} |k|}{\lambda \mu _w \sqrt{2|Q_\infty |}} \int _0^\infty Ai(q^{\frac{2}{3}} z') \mathrm{{d}}z' } \mathrm{{e}}^{-|k| \sqrt{2|Q_\infty |} y_7}, \end{aligned}$$
where \(\overline{p}_{st}\) is the Fourier transform of \(p_{st}\), and \(\overline{f}\) is the Fourier transform of the roughness shape. To find the solution for supersonic flow (\(Q_\infty >0\)), we have to take into account the causality condition, which yields the solutions
$$\begin{aligned} \overline{P}_5&=\frac{\mathrm{{i}} k}{ \sqrt{2 Q_\infty }}\frac{ \overline{f}(k)q^{\frac{2}{3}} Ai'(0)}{q^{\frac{2}{3}} Ai'(0)+\mathrm{{i}} \frac{(\mathrm{{i}} k)^{\frac{1}{3}} k}{\lambda \mu _w \sqrt{2Q_\infty }} \int _0^\infty Ai(q^{\frac{2}{3}} z') \mathrm{{d}}z' } , \end{aligned}$$
$$\begin{aligned} \overline{p}_{st}&=\frac{i k}{ \sqrt{2 Q_\infty }}\frac{ \overline{f}(k)q^{\frac{2}{3}} Ai'(0)}{q^{\frac{2}{3}} Ai'(0)+\mathrm{{i}} \frac{(\mathrm{{i}} k)^{\frac{1}{3}} k}{\lambda \mu _w \sqrt{2Q_\infty }} \int _0^\infty Ai(q^{\frac{2}{3}} z') \mathrm{{d}}z' } \mathrm{{e}}^{-\mathrm{{ik}} \sqrt{2Q_\infty } y_7}. \end{aligned}$$
At this point, we note that the branch cut of \((\mathrm{{i}} k)^{\frac{1}{3}}\) in the complex k-plane is chosen in the usual way to be the positive imaginary semi-axis.
Unsteady problem
In order to solve our unsteady problem, we consider the terms of \(O(\sigma \delta )\) in our linearised unsteady lower deck equations which give
$$\begin{aligned}&\rho _w \left( \frac{\partial \tilde{u}_5'}{\partial t_*}+\lambda y_5 \frac{\partial \tilde{u}_5'}{\partial x_*}+\lambda \tilde{v}_5'\right) =-\frac{\partial \tilde{P}_5'}{\partial x_*}+\mu _w \frac{\partial ^2 \tilde{u}_5'}{\partial y_7^2}, \end{aligned}$$
$$\begin{aligned}&\frac{\partial \tilde{P}_5'}{\partial y_5}=0, \end{aligned}$$
$$\begin{aligned}&\frac{\partial \tilde{u}_5'}{\partial x_*}+\frac{\partial \tilde{v}_5'}{\partial y_5}=0. \end{aligned}$$
The boundary conditions for this problem are
$$\begin{aligned} \tilde{u}_5'&=\tilde{v}_5'=0 \quad \mathrm {at} \quad y_5=0, \end{aligned}$$
$$\begin{aligned} \tilde{u}_5'&=\tilde{A}_*' \quad \mathrm {as} \quad y_5\rightarrow \infty , \end{aligned}$$
$$\begin{aligned} \tilde{u}_5'&=0 \quad \mathrm {as} \quad x_* \rightarrow - \infty . \end{aligned}$$
In order to make analytical progress in finding solutions to these equations, we shall assume that our solutions are periodic in time, with wavenumber \(\alpha \) corresponding to the wavenumber of \(u_\alpha (t)\). This means we can write them in the form:
$$\begin{aligned} \left( \tilde{u}_5', \tilde{v}_5', \tilde{P}_5', p_*' , \tilde{A}_*' \right) = \mathrm{{e}}^{\mathrm{{i}} \alpha t} \left( u_u, v_u, P_u, p_u , A_u \right) +\mathrm{{c}} \cdot \mathrm{{c}} \end{aligned}$$
where the functions on the right-hand side have no time dependence, and \(\mathrm{{c}}\cdot \mathrm{{c}}\) denotes the complex conjugate. Substituting these expressions into (5.10) and (5.11) yields
$$\begin{aligned}&\rho _w \left( \mathrm{{i}} \alpha u_u + \lambda y_5 \frac{\partial u_u}{\partial x_*}+\lambda v_u \right) = - \frac{\partial P_u}{\partial x_*}+ \mu _w \frac{\partial ^2 u_u}{\partial y_7^2}, \end{aligned}$$
$$\begin{aligned}&\frac{\partial P_u}{\partial y_5}=0, \end{aligned}$$
$$\begin{aligned}&\frac{\partial u_u}{\partial x_*}+ \frac{\partial v_u}{\partial y_5}=0, \end{aligned}$$
$$\begin{aligned}&u_u=v_u=0 \quad \mathrm {at} \quad y_5=0, \end{aligned}$$
$$\begin{aligned}&u_u= A_u \quad \mathrm {as} \quad y_5\rightarrow \infty , \end{aligned}$$
$$\begin{aligned}&u_u=0 \quad \mathrm {as} \quad x_* \rightarrow - \infty . \end{aligned}$$
Performing the same substitution in the boundary value problem (4.29) derived for the pressure in the upper deck \(p_*\) gives
$$\begin{aligned}&2\mathrm{{i}} \alpha \frac{\partial p_u}{\partial x_*}+ 2 Q_\infty \frac{\partial ^2 p_u}{\partial x_*^2}- \frac{\partial ^2 p_u}{\partial y_7^2}=-2 \frac{\partial ^2 p_{st}}{\partial x_*^2}, \end{aligned}$$
$$\begin{aligned}&p_u \rightarrow 0 \quad \mathrm {as} \quad x_*^2+y_7^2 \rightarrow \infty , \end{aligned}$$
$$\begin{aligned}&\frac{\partial p_u}{\partial y_7}=\frac{1}{\lambda }\frac{\mathrm{{d}}^2 A_u}{\mathrm{{d}} x_*^2} \quad \mathrm {at} \quad y_7=0. \end{aligned}$$
In similar fashion to how we approached the steady problem, we begin by Fourier transforming our variables in \(x_*\), so that the above problem for \(p_u\) becomes
$$\begin{aligned}&\frac{\mathrm{{d}}^2 \overline{p}_u}{\mathrm{{d}}y_7^2}+ \overline{p}_u\left( 2k\alpha +2Q_\infty k^2 \right) =-2k^2 \overline{p}_{st}, \end{aligned}$$
$$\begin{aligned}&\overline{p}_u \rightarrow 0 \quad \mathrm {as} \quad y_7 \rightarrow \infty , \end{aligned}$$
$$\begin{aligned}&\frac{\mathrm{{d}}\overline{p}_u}{\mathrm{{d}}y_7}=-\frac{k^2}{\lambda }\overline{A}_u \quad \mathrm {at} \quad y_7=0, \end{aligned}$$
where \(\overline{p}_u\) and \(\overline{A}_u\) are the Fourier transforms of \(p_u\) and \(A_u\), respectively, and \(\overline{p}_{st}\) was found in our analysis of the steady problem. The solution to this problem of course depends on the nature of flow (i.e. whether it is subsonic or supersonic) due to the fact that the forcing on the right-hand side of (5.15a) depends on the sign of \(Q_\infty \). In addition, as we need to use an inverse Fourier transform to return to physical space, we require a solution for all k. Therefore, for each sign of \(Q_\infty \), we have to split our solution into 3 different regions. For \(Q_\infty <0\), we have
\(k<0\), \(\quad \quad \quad \quad \quad 2k \alpha +2Q_\infty k^2<0\);
\(k \in \left( 0, -\frac{\alpha }{Q_\infty } \right) \), \(\quad 2k\alpha +2Q_\infty k^2>0\);
\(k>-\frac{\alpha }{Q_\infty }\), \(\quad \quad \quad 2k\alpha +2Q_\infty k^2<0\).
For \(Q_\infty >0\), we have
\(k<-\frac{\alpha }{Q_\infty }\), \(\quad \quad \quad 2k \alpha +2Q_\infty k^2>0\);
\(k \in \left( -\frac{\alpha }{Q_\infty },0 \right) \), \(\quad 2k\alpha +2Q_\infty k^2<0\);
\(k>0\), \(\quad \quad \quad \quad \quad 2k\alpha +2Q_\infty k^2>0\).
We begin by considering the case of subsonic flow.
Subsonic flow
In the case of subsonic flow (\(Q_\infty <0\)), the solution to the system (5.15) is
$$\begin{aligned} \overline{p}_u=\frac{\frac{k^2}{\lambda } \overline{A}_u+ \frac{k}{\alpha }K_1|k| \sqrt{2 |Q_\infty |}}{\Omega } \mathrm{{e}}^{-\Omega y_7}-\frac{k K_1}{\alpha } \mathrm{{e}}^{-|k| \sqrt{2 |Q_\infty |}y_7}, \end{aligned}$$
where \(\Omega \) is
$$\begin{aligned} \Omega ={\left\{ \begin{array}{ll} \mathrm{{i}} \sqrt{2 k^2 |Q_\infty |+2 |k| \alpha } \ \quad k<0, \\ \sqrt{2 k \alpha -2k^2 |Q_\infty |} \quad \ \ k \in \left( 0, -\frac{\alpha }{Q_\infty } \right) ,\\ \mathrm{{i}} \sqrt{2 k^2 |Q_\infty |-2k \alpha } \quad \quad k> -\frac{\alpha }{Q_\infty }, \end{array}\right. } \end{aligned}$$
and \(K_1\) is defined through
$$\begin{aligned} K_1=\frac{ |k|}{ \sqrt{2 |Q_\infty |}}\frac{ \overline{f}(k)q^{\frac{2}{3}} Ai'(0)}{q^{\frac{2}{3}} Ai'(0)- \frac{(\mathrm{{i}} k)^{\frac{1}{3}} |k|}{\lambda \mu _w \sqrt{2|Q_\infty |}} \int _0^\infty Ai(q^{\frac{2}{3}} z') \mathrm{{d}}z' }. \end{aligned}$$
Relating \(\overline{P}_u\) to \(\overline{p}_u\) using the condition (5.6), we have that
$$\begin{aligned} \overline{P}_u=\frac{\frac{k^2}{\lambda } \overline{A}_u+ \frac{k}{\alpha }K_1|k| \sqrt{2 |Q_\infty |}}{\Omega }- \frac{k}{\alpha }K_1. \end{aligned}$$
We now return to our linearised lower deck equations. After performing Fourier transforms and eliminating \(\overline{v}_u\) (the Fourier transform of \(v_u\)) in a similar fashion to the analysis performed for the steady problem, we obtain the boundary value problem for the linearised lower deck variables in Fourier space
$$\begin{aligned}&\mathrm{{i}}\rho _w \frac{\mathrm{{d}}\overline{u}_u}{\mathrm{{d}}y_5}(\alpha + \lambda y_5 k)= \mu _w \frac{\mathrm{{d}}^3 \overline{u}_u}{\mathrm{{d}}y_5^3}, \end{aligned}$$
$$\begin{aligned}&\overline{u}_u=0 \quad \mathrm {at} \quad y_5=0, \end{aligned}$$
$$\begin{aligned}&\overline{u}_u= \overline{A}_u \quad \mathrm {as} \quad y_5 \rightarrow \infty , \end{aligned}$$
$$\begin{aligned}&\frac{\mathrm{{d}}^2 \overline{u}_u}{\mathrm{{d}}y_5^2}= \frac{\mathrm{{i}} k}{\mu _w} \left( \frac{\frac{k^2}{\lambda } \overline{A}_u+ \frac{k}{\alpha }K_1|k| \sqrt{2 |Q_\infty |}}{\Omega }- \frac{k}{\alpha }K_1 \right) \quad \mathrm {at} \quad y_5=0. \end{aligned}$$
We now make the transformations
$$\begin{aligned} \zeta =\lambda (\mathrm{{i}} k)^{\frac{1}{3}} y_5+\zeta _0, \quad \zeta _0= \frac{\mathrm{{i}} \alpha }{(\mathrm{{i}}k)^{\frac{2}{3}}}. \end{aligned}$$
At this point, we note that we are required to introduce a branch cut of \((\mathrm{{i}}k)^{\frac{1}{3}}\), which we will again place along the positive imaginary semi-axis in the k-plane. Our boundary problem for \(\overline{u}_u\) is, therefore,
$$\begin{aligned}&\frac{\mathrm{{d}}^3 \overline{u}_u}{\mathrm{{d}} \zeta ^3}-\frac{\rho _w}{\lambda \mu _w}\zeta \frac{\mathrm{{d}}\overline{u}_u}{\mathrm{{d}}\zeta }=0, \end{aligned}$$
$$\begin{aligned}&\frac{\mathrm{{d}}^2 \overline{u}_u}{\mathrm{{d}}\zeta ^2}=\frac{(\mathrm{{i}} k)^{\frac{1}{3}}}{\mu _w \lambda }\left( \frac{\frac{k^2}{\lambda } \overline{A}_u+ \frac{k}{\alpha }K_1|k| \sqrt{2 |Q_\infty |}}{\Omega }- \frac{k}{\alpha }K_1 \right) \quad \mathrm {at} \quad \zeta =\zeta _0, \end{aligned}$$
$$\begin{aligned}&\overline{u}_u=\overline{A}_u \quad \mathrm {as} \quad \zeta \rightarrow \infty , \end{aligned}$$
$$\begin{aligned}&\overline{u}_u=0 \quad \mathrm {at} \quad \zeta = \zeta _0. \end{aligned}$$
The general form of the solution to the ODE above is
$$\begin{aligned} \frac{\mathrm{{d}}\overline{u}_u}{\mathrm{{d}}\zeta }=C Ai\left( q_1^{\frac{2}{3}} \zeta \right) + D Bi\left( q_1^{\frac{2}{3}} \zeta \right) \end{aligned}$$
where \(q_1=\sqrt{\frac{\rho _w}{\lambda \mu _w}}\), and Ai and Bi are the Airy functions. In order to satisfy (5.22c), we require \(D=0\); otherwise, the derivative of \(\overline{u}_u\) will be infinite while \(\overline{u}_u\) is tending to a finite value which is a contradiction. Differentiating this solution once and matching it with (5.22b), we have
$$\begin{aligned} q_1^{\frac{2}{3}}C Ai'\left( q_1^{\frac{2}{3}} \zeta _0\right) =\frac{(\mathrm{{i}} k)^{\frac{1}{3}}}{\mu _w \lambda }\left( \frac{\frac{k^2}{\lambda } \overline{A}_u+ \frac{k}{\alpha }K_1|k| \sqrt{2 |Q_\infty |}}{\Omega }- \frac{k}{\alpha }K_1 \right) . \end{aligned}$$
In addition, we can integrate (5.23) and considering (5.22d), we obtain
$$\begin{aligned} \overline{u}_u=C \int _{\zeta _0}^{\zeta } Ai\left( q_1^{\frac{2}{3}}\zeta '\right) \mathrm{{d}} \zeta '. \end{aligned}$$
In the limit \(\zeta \rightarrow \infty \), we can deduce from (5.22c) that
$$\begin{aligned} \overline{A}_u=C \int _{\zeta _0}^\infty Ai\left( q_1^{\frac{2}{3}}\zeta '\right) \mathrm{{d}} \zeta '. \end{aligned}$$
Eliminating C using (5.24), we have that \(\overline{A}_u\) can be expressed in the following form:
$$\begin{aligned} \overline{A}_u=\frac{(\mathrm{{i}} k)^{\frac{1}{3}}}{\alpha } \frac{\lambda k K_1 R_1 I_1}{\lambda ^2 \mu _w q_1^{\frac{2}{3}}Ai'\left( q_1^{\frac{2}{3}} \zeta _0\right) \Omega -(\mathrm{{i}} k)^{\frac{1}{3}} k^2I_1 }, \end{aligned}$$
where we have defined the parameters:
$$\begin{aligned} R_1=|k| \sqrt{2 |Q_\infty |}-\Omega , \quad I_1= \int _{\zeta _0}^\infty Ai\left( q_1^{\frac{2}{3}}\zeta '\right) \mathrm{{d}}\zeta ' . \end{aligned}$$
Substituting this into (5.19), we find
It is now necessary to return to physical space. Performing an inverse Fourier transform for \(\overline{P}_u\) yields
$$\begin{aligned} \tilde{P}_5'(x_*, t)= \frac{\mathrm{{e}}^{\mathrm{{i}} \alpha t}}{2 \pi }\int _{-\infty }^\infty \frac{\lambda ^2 \mu _w k K_1 R_1 q_1^{\frac{2}{3}}Ai'\left( q_1^{\frac{2}{3}} \zeta _0\right) }{\alpha \left( \lambda ^2 \mu _w q_1^{\frac{2}{3}}Ai'(q_1^{\frac{2}{3}} \zeta _0)\Omega -(\mathrm{{i}}k)^{\frac{1}{3}}k^2 I_1 \right) } \mathrm{{e}}^{\mathrm{{i}} k x_*} \mathrm{{d}}k. \end{aligned}$$
In order to evaluate this integral, we are going to deform the integration contour into the complex plane. In preparation for this, we shall evaluate the poles of the integrand. These occur when we have
$$\begin{aligned} \mu _w \lambda ^2 q_1^{\frac{2}{3}}Ai'\left( q_1^{\frac{2}{3}} \zeta _0\right) \Omega -(\mathrm{{i}}k)^{\frac{1}{3}} k^2 \int _{\zeta _0}^\infty Ai\left( q_1^{\frac{2}{3}}\zeta '\right) \mathrm{{d}}\zeta '=0. \end{aligned}$$
This is a dispersion relation for k. In the limit of \(k\rightarrow 0\), if we assume that \(\zeta _0\) remains finite (i.e. \(\alpha \sim k^{{2}/{3}}\)), then this equation becomes
$$\begin{aligned} Ai'\left( q_1^{\frac{2}{3}} \zeta _0\right) =0. \end{aligned}$$
It is known that this equation has an infinite number of roots, all lying on the negative real semi-axis of the \(\zeta _0\)-plane. This suggests that the dispersion relation has an infinite number of roots as well. In order to study the behaviour of the first few roots of this equation, it is only necessary to consider the behaviour of roots in the limit \(k\rightarrow \infty \). Assuming that \(\zeta _0\) again remains finite, this is equivalent to solving
$$\begin{aligned} \int _{\zeta _0}^\infty Ai \left( q_1^{\frac{2}{3}}\zeta '\right) \mathrm{{d}}\zeta '=0. \end{aligned}$$
Similar to Eq. (5.31), this equation has infinite roots. It is known from properties of the Airy function that these come in complex conjugate pairs and all lie in the negative real half of the \(\zeta _0\)-plane. It can be seen that the first root of (5.31) behaves differently in the \(\zeta _0\)-plane. As \(k \rightarrow \infty \), instead of tending to a finite point which is a root of (5.32), it instead tends to infinity. This root also behaves differently in the k-plane. As \(\alpha \) goes from zero to infinity, all other roots remain in the second quadrant of the k-plane for all \(\alpha \). The first root, however, starts in the second quadrant and crosses the negative real semi-axis at a critical point \(k_*<0\) which has a corresponding critical frequency \(\alpha _*\). For all \(\alpha >\alpha _*\), this root remains in the third quadrant. This is the root which corresponds to the wavenumber of the Tollmien–Schlichting wave. Trajectories of the first five roots in both the k and \(\zeta _0\) planes can be found in Fig. 2, where the root corresponding to the Tollmien–Schlichting wave is in bold.
Behaviour of the first five roots of (5.30) with \(Q=-1\) in both the k-plane and the \(\zeta _0\)-plane
The process used to calculate the behaviour of the roots as \(\alpha \) increased was as follows. First, Eq. (5.30) is written in terms of \(\zeta _0\) only, using (5.21). Then, Eq. (5.31) is used to find the position of each root at \(\alpha =0\). For each subsequent value of \(\alpha \), the previous value of \(\zeta _0\) was used as an initial guess, and the dispersion relation was solved using Newton Iterations.
Returning to (5.29), we deform our contour of integration in the k-plane as shown in Fig. 3. We have had to split our contours into the regions of \(k<0\) and \(k>0\) due to the branch cut we introduced along the positive imaginary semi-axis. We, therefore, split our integral up in the following way:
$$\begin{aligned} \tilde{P}_5'(x_*, t)= \frac{\mathrm{{e}}^{\mathrm{{i}} \alpha t}}{2 \pi } \left( \int _{-\infty }^0+\int _0^{-\frac{\alpha }{2Q_\infty }}+\int _{-\frac{\alpha }{2 Q_\infty }}^\infty \right) \frac{\lambda ^2 \mu _w k K_1 R_1 q_1^{\frac{2}{3}}Ai'\left( q_1^{\frac{2}{3}} \zeta _0\right) }{\alpha \left( \lambda ^2 \mu _w q_1^{\frac{2}{3}}Ai'\left( q_1^{\frac{2}{3}} \zeta _0\right) \Omega -(\mathrm{{i}}k)^{\frac{1}{3}}k^2 I_1 \right) } \mathrm{{e}}^{\mathrm{{i}} k x_*} \mathrm{{d}}k. \end{aligned}$$
Deformed integration contour for subsonic flow
If we let the radii of the segments shown in Fig. 3 go to infinity then by Jordan's lemma, we can disregard the integrals along the arcs. Using Watson's lemma, we find that
$$\begin{aligned}&\int _{C_-}\frac{\lambda ^2 \mu _w k K_1 R_1 q_1^{\frac{2}{3}}Ai'\left( q_1^{\frac{2}{3}} \zeta _0\right) }{\alpha \left( \lambda ^2 \mu _w q_1^{\frac{2}{3}}Ai'\left( q_1^{\frac{2}{3}} \zeta _0\right) \Omega -(\mathrm{{i}}k)^{\frac{1}{3}}k^2 I_1 \right) } \mathrm{{e}}^{\mathrm{{i}} k x_*} \mathrm{{d}}k=-\frac{2\mathrm{{i}}\overline{f}(0)}{\alpha \sqrt{2|Q_\infty |}}\frac{1}{x_*^3}, \end{aligned}$$
$$\begin{aligned}&\int _{C_+}\frac{\lambda ^2 \mu _w k K_1 R_1 q_1^{\frac{2}{3}}Ai'\left( q_1^{\frac{2}{3}} \zeta _0\right) \mathrm{{e}}^{\mathrm{{i}} k x_*}}{\alpha \left( \lambda ^2 \mu _w q_1^{\frac{2}{3}}Ai'\left( q_1^{\frac{2}{3}} \zeta _0\right) \Omega -(\mathrm{{i}}k)^{\frac{1}{3}}k^2 I_1 \right) } \mathrm{{d}}k=\frac{2\mathrm{{i}} \overline{f}(0)}{\alpha \sqrt{2|Q_\infty |}}\frac{1}{x_*^3}. \end{aligned}$$
As the integrals along \(C_a\) and \(C_b\) cancel out, by the residue theorem, we find that at large values x, the Tollmien-Schlichting wave emerges from the perturbation field with
$$\begin{aligned} P_{TS}(t, x_*)=\Lambda (\alpha , Q_\infty )\overline{f}(k_1) \mathrm{{e}}^{\mathrm{{i}}(\alpha t +k_1 x)}+\cdots \end{aligned}$$
where \(\Lambda \) is the receptivity coefficient found to be
$$\begin{aligned} \Lambda (\alpha , Q_\infty )= \frac{-\frac{3}{\alpha }\mathrm{{i}}\lambda ^2 \mu _w k_1 K_{\Lambda } R_{\Lambda } q_1^{\frac{2}{3}}Ai'\left( q_1^{\frac{2}{3}} \zeta _{\Lambda }\right) }{ 2 Ai\left( q_1^{\frac{2}{3}} \zeta _{\Lambda }\right) \frac{\zeta _{\Lambda }}{k_1} C_\Lambda +7(\mathrm{{i}}k_1)^{\frac{1}{3}} k_1 I_{\Lambda }-\frac{3(2|Q_\infty |k_1+\alpha )}{\Omega _{\Lambda }} \lambda ^2 \mu _w q_1^{\frac{2}{3}} Ai'\left( q_1^{\frac{2}{3}} \zeta _{\Lambda }\right) } \end{aligned}$$
with the parameters with subscript \(\Lambda \) defined as follows:
$$\begin{aligned} \Omega _\Lambda&= \sqrt{-2k_1^2 |Q_\infty |-2|k_1|\alpha },\quad \quad R_{\Lambda }=|k_1|\sqrt{2|Q_\infty |}-\Omega _{\Lambda }, \end{aligned}$$
$$\begin{aligned} K_{\Lambda }&=\frac{ |k_1|}{ \sqrt{2 |Q_\infty |}}\frac{ q^{\frac{2}{3}} Ai'(0)}{q^{\frac{2}{3}} Ai'(0)- \frac{(\mathrm{{i}} k_1)^{\frac{1}{3}} |k_1|}{\lambda \mu _w \sqrt{2|Q_\infty |}} \int _0^\infty Ai(q^{\frac{2}{3}} z') \mathrm{{d}}z' }, \end{aligned}$$
$$\begin{aligned} \zeta _{\Lambda }&=\frac{(\mathrm{{i}}\alpha )}{(\mathrm{{i}}k_1)^{\frac{2}{3}}}, \quad \quad I_{\Lambda }=\int _{\zeta _{\Lambda }}^\infty Ai(q_1^{\frac{2}{3}}\zeta ')\mathrm{{d}}\zeta ',\quad \quad C_\Lambda \Big ( 2(\mathrm{{i}}k_1)^{\frac{1}{3}} k_1^2-\lambda ^2 \mu _w \Omega _{\Lambda } q^{\frac{7}{3}}\zeta _{\Lambda } \Big ). \end{aligned}$$
Here, it is important to note that the receptivity coefficient is a function of the Kármán–Guderley parameter \(Q_\infty \) and frequency \(\alpha \) only, and does not depend on the roughness shape. Numerical calculations of the receptivity coefficient can be found in Fig. 4.
Subsonic receptivity coefficient modulus and argument
Supersonic flow
For supersonic flow, the analysis proceeds in much the same way. The solution to the problem (5.15) in this case is
$$\begin{aligned} \overline{p}_u=\frac{\frac{k^2}{\lambda } \overline{A}_u+ \mathrm{{i}} \frac{k^2}{\alpha }K_2 \sqrt{2 |Q_\infty |}}{\Omega }\mathrm{{e}}^{-\Omega y_7}-\frac{k}{\alpha }K_2 \mathrm{{e}}^{-\mathrm{{i}}k\sqrt{2Q_\infty }y_7}, \end{aligned}$$
where \(\Omega \) is now defined to be
$$\begin{aligned} \Omega = {\left\{ \begin{array}{ll} \sqrt{2k^2 Q_\infty -2|k|\alpha } \qquad \ k<-\frac{\alpha }{Q_\infty }, \\ \mathrm{{i}}\sqrt{2|k|\alpha -2 k^2Q_\infty } \qquad k \in \left( -\frac{\alpha }{Q_\infty }, 0 \right) ,\\ \sqrt{2k\alpha + 2k^2Q_\infty } \qquad k>0, \end{array}\right. } \end{aligned}$$
and \(K_2\) is
$$\begin{aligned} K_2=\frac{\mathrm{{i}} k}{ \sqrt{2 Q_\infty }}\frac{ \overline{f}(k)q^{\frac{2}{3}} Ai'(0)}{q^{\frac{2}{3}} Ai'(0)+i \frac{(\mathrm{{i}} k)^{\frac{1}{3}} k}{\lambda \mu _w \sqrt{2Q_\infty }} \int _0^\infty Ai(q^{\frac{2}{3}} z') \mathrm{{d}}z' }. \end{aligned}$$
Performing an identical analysis to the subsonic case, we find that for \(Q_\infty >0\)
$$\begin{aligned} \overline{A}_u&=\frac{\lambda (\mathrm{{i}}k)^{\frac{1}{3}}kK_2R_2I_1}{\alpha \left( \lambda ^2 \mu _w q_1^{\frac{2}{3}}Ai'\left( q_1^{\frac{2}{3}}\zeta _0\right) \Omega -(\mathrm{{i}}k)^{\frac{1}{3}}k^2I_1\right) }, \end{aligned}$$
$$\begin{aligned} \overline{P}_u&=\frac{q_1^{\frac{2}{3}} \lambda ^2 \mu _w Ai'(q_1^{\frac{2}{3}}\zeta _0)k K_2 R_2}{\alpha \left( \lambda ^2 \mu _w q_1^{\frac{2}{3}}Ai'\left( q_1^{\frac{2}{3}}\zeta _0\right) \Omega -(\mathrm{{i}}k)^{\frac{1}{3}}k^2I_1\right) }, \end{aligned}$$
where \(q_1\), \(\zeta _0\), and \(I_1\) are the same as in the subsonic case and \(R_2\) is defined through
$$\begin{aligned} R_2=\mathrm{{i}}k\sqrt{2 Q_\infty }-\Omega . \end{aligned}$$
Performing an inverse Fourier transform in order to return back to physical space, we have
$$\begin{aligned} \tilde{P}_5'(t, x_*)=\frac{\mathrm{{e}}^{\mathrm{{i}}\alpha t}}{2\pi } \int _{-\infty }^\infty \frac{q_1^{\frac{2}{3}} \lambda ^2 \mu _w Ai'(q_1^{\frac{2}{3}}\zeta _0)k K_2 R_2}{\alpha \left( \lambda ^2 \mu _w q_1^{\frac{2}{3}}Ai'\left( q_1^{\frac{2}{3}}\zeta _0\right) \Omega -(\mathrm{{i}}k)^{\frac{1}{3}}k^2I_1\right) }\mathrm{{e}}^{\mathrm{{i}}kx_*}\mathrm{{d}}k. \end{aligned}$$
This yields the dispersion relation for supersonic flow
$$\begin{aligned} \lambda ^2 \mu _w q_1^{\frac{2}{3}}Ai'\left( q_1^{\frac{2}{3}}\zeta _0\right) \Omega -(ik)^{\frac{1}{3}}k^2 \int _{\zeta _0}^\infty Ai\left( q_1^{\frac{2}{3}}\zeta '\right) \mathrm{{d}}\zeta '=0. \end{aligned}$$
This is nearly identical to the dispersion relation for subsonic flow (see 5.30), except for change in definition of \(\Omega \). The roots of this equation behave in a similar manner in both the k and \(\zeta _0\) planes as in the subsonic case. In order to analyse the integral (5.46), we again split the integral into three regions and deform each line into a contour in the complex k-plane. The deformed integration contour can be seen in Fig. 5. The behaviour of poles of the integrand is the same as in the subsonic case. We note that it does not matter which contour contains the root \(k_1\). Again, we can disregard the integral along both arcs in view of Jordan's Lemma, and making use of Watson's Lemma, we find the integrals along \(C_1\), \(C_2\), and \(C_3\) sum to zero at large values of \(x_*\). Therefore, using residue calculus, we have that for supersonic flow, the Tollmien–Schlichting wave emerges from the perturbation field such that
Deformed integration contour in the k-plane for supersonic flow
$$\begin{aligned} P_{TS}=\Lambda (\alpha , Q_\infty )\overline{f}(k_1) \mathrm{{e}}^{\mathrm{{i}}(\alpha t +k_1 x_*)}+\cdots , \end{aligned}$$
where in this case, the receptivity coefficient \(\Lambda \) is
$$\begin{aligned} \Lambda (\alpha , Q_\infty )=\frac{\frac{3}{\alpha }\mathrm{{i}}\lambda ^2 \mu _w k_1 K_{\Lambda } R_{\Lambda } q_1^{\frac{2}{3}}Ai'(q_1^{\frac{2}{3}} \zeta _{\Lambda })}{ 2 Ai(q_1^{\frac{2}{3}} \zeta _{\Lambda })\frac{\zeta _{\Lambda }}{k_1}C_{\Lambda } +7(\mathrm{{i}}k_1)^{\frac{1}{3}} k_1 I_{\Lambda }-\frac{3(2 Q_\infty k_1+\alpha )}{\Omega _{\Lambda }} \lambda ^2 \mu _w q_1^{\frac{2}{3}} Ai'(q_1^{\frac{2}{3}} \zeta _{\Lambda }) } \end{aligned}$$
with the coefficients having the subscript \(\Lambda \) defined through
$$\begin{aligned} \Omega _{\Lambda }&= \mathrm{{i}}\sqrt{2|k_1|\alpha -2k_1^2 Q_\infty },\quad \quad R_{\Lambda }=\mathrm{{i}} k_1 \sqrt{2Q_\infty }-\Omega _{\Lambda }, \end{aligned}$$
$$\begin{aligned} K_{\Lambda }&=\frac{\mathrm{{i}} k_1}{ \sqrt{2 Q_\infty }}\frac{ q^{\frac{2}{3}} Ai'(0)}{q^{\frac{2}{3}} Ai'(0)- i\frac{(\mathrm{{i}} k_1)^{\frac{1}{3}} k_1}{\lambda \mu _w \sqrt{2Q_\infty }} \int _0^\infty Ai\left( q^{\frac{2}{3}} z'\right) \mathrm{{d}}z' }, \end{aligned}$$
$$\begin{aligned} \zeta _{\Lambda }&=\frac{(\mathrm{{i}}\alpha )}{(ik_1)^{\frac{2}{3}}}, \quad \quad I_{\Lambda }=\int _{\zeta _{\Lambda }}^\infty Ai\left( q_1^{\frac{2}{3}}\zeta '\right) \mathrm{{d}}\zeta ', \quad \quad C_\Lambda =\left( 2(\mathrm{{i}}k_1)^{\frac{1}{3}} k_1^2-\lambda ^2 \mu _w \Omega _{\Lambda } q^{\frac{7}{3}}\zeta _{\Lambda } \right) . \end{aligned}$$
While similar to the subsonic case, we note that the receptivity coefficient for supersonic flow is not identical to (5.36). The results of the numerical calculations for receptivity coefficient in supersonic flow can be found in Fig. 6.
Supersonic receptivity coefficient modulus and argument
This paper studied the receptivity of the Blasius boundary layer to free-stream turbulence in transonic flow. As the receptivity analysis is local, it is, therefore, applicable to an arbitrary boundary layer, for example, a wing surface. We have seen that, unlike acoustic waves, vorticity waves do not produce pressure perturbations to leading order in the free stream or the vorticity deformation layer. Therefore, the vorticity waves cannot penetrate into the boundary layer, so the process by which Tollmien–Schlichting waves are generated differs greatly from the acoustic case. In the present case, Tollmien–Schlichting waves arise due to the interaction between the free-stream turbulence and a roughness element on the plate surface. In order for this interaction to take place, we required that the strict double resonance condition associated with Tollmien–Schlichting waves be followed, i.e. the frequency of the vorticity wave is close to the frequency of the Tollmien–Schlichting wave, and in the Fourier spectrum of the perturbations produced by the roughness element is a wave number matching the wave number of the Tollmien–Schlichting wave. Unlike in the case of the acoustic wave, the interaction of the vorticity wave with steady perturbations (produced by the roughness element) takes place in the upper tier of the triple-deck structure.
After first non-dimensionalising the governing equations, solutions in the flow upstream of the roughness element were derived. Due to the lack of leading order pressure perturbations, it was necessary to introduce a transition layer at the outer edge of the boundary layer in which the velocity perturbations caused by the free-stream turbulence decayed. The analysis of the boundary layer was, therefore, simple, as there was no need to introduce a Stokes layer close to the plate as is necessary in the case of acoustic waves.
To analyse the flow around the roughness, we used the scalings corresponding to transonic viscous–invisicid interaction theory, which enabled equations governing the velocity and pressure distributions in the three regions of the triple-deck structure to be derived. We note at this point that while the analysis performed on the lower and middle decks is similar to previous research, the inclusion of free-stream turbulence meant that the asymptotic expansions of the upper deck were novel. The external perturbations present in this deck interacted with the perturbations caused by the roughness element to produce the Tollmien–Schlichting wave.
When finding the amplitude of the Tollmien–Schlichting wave, it was necessary to assume that the height of the roughness element was small compared to \(L \mathrm{{Re}}^{-\frac{11}{18}}\) in order to make analytic progress. This assumption allowed the equations for the lower deck of the triple-deck model to be linearised, so the resulting viscous–inviscid interaction problems could be solved analytically. The behaviour of Tollmien–Schlichting waves far downstream of the roughness element is of primary interest when modelling laminar–turbulent transition as noted in the Introduction. By first solving the steady problem, and then assuming the unsteady solutions were periodic in time before making use of Fourier transforms, elementary results from complex analysis, and Watson's Lemma, a formula for the amplitude of Tollmien–Schlichting waves in this limit was found, depending only on the frequency of the vorticity wave and the Kármán–Guderley parameter. It is interesting that unlike the acoustic case studied by [7], the receptivity coefficient did depend on the sign of the Kármán–Guderley parameter.
To conclude, in this paper, we have derived an explicit equation for the amplitude of Tollmien–Schlichting waves generated by free-stream turbulence in transonic flow far downstream of roughness. This first required solving the Navier–Stokes equations in the various regions upstream of the roughness, before introducing the triple-deck model of [8, 9], the use of which allowed the equations governing the flow around the roughness element to be derived. The receptivity coefficient for both subsonic and supersonic flow were found and calculated numerically, and depended only on the frequency of the vorticity wave and the Kármán–Guderley parameter.
Of course, the generation of the Tollmien–Schlichting waves can be also studied by means of the direct numerical simulations of the Navier–Stokes equations. However, a clear advantage of the asymptotic theory presented here is that it provides an explicit dependence of the amplitude of the Tollmien–Schlichting wave on the shape of the wall roughness. We see that in order to delay the laminar–turbulent transition, one needs to minimise the Fourier transform of the wall roughness function at a particular value of the wave number that corresponds to the Tollmien–Schlichting wave. This result may be directly used to improve the wing design. For example, of particular interest is a small gap between the nasal retractable section of the wing with the main body of the wing. While it cannot be avoided without increasing the landing speed of an aircraft, the effect of this gap can be minimised by simply changing the gap width.
Schubauer GB, Skramstad HK (1948) Laminar boundary oscillations and transition on a flat plate. NACA Rep 38:251–292
Fedorov AV (1982) Generation of instability waves in a boundary-layer flow of compressbile gas exposed to an acoustic field. Num Meth Continuum Mech (Novosibirsk) 13:106
Goldstein ME (1983) The evolution of Tollmien–Schlichting waves near a leading edge. J Fluid Mech 127:59–81
Article MathSciNet Google Scholar
Goldstein ME (1985) Scattering of acoustic waves into Tollmien–Schlichting waves by small streamwise variations in surface geometry. J Fluid Mech 154:509–529
Ruban AI (1984) On the generation of Tollmien–Schlichting waves by sound. Fluid Dyn 19:709–717
Kachanov YS, Kozlov VV, Levchenko VY (1982) Beginning of turbulence in the boundary layer. Nauka, Novosibirsk
Ruban AI, Bernots T, Kravtsova MA (2016) Linear and nonlinear receptivity of the boundary layer in transonic flows. J Fluid Mech 786:154–189
Bowles RI, Smith FT (1993) On boundary-layer transition in transonic flow. J Eng Math 27:309–342
Timoshin S (1990) Asymptotic form of the lower branch of the neutral curve in transonic flow. Uch Zap TsAGI 21(6):50–57
Hunt JCR, Graham JMR (1978) Free-stream vorticity near plate boundaries. J Fluid Mech 84:209–235
Guliaev AN, Kozlov VE, Kuznetzov VR, Mineev BI, Secundov AN (1989) Interaction of a laminar boundary layer with external turbulence. Fluid Dyn 24:700–710
Duck PW, Ruban AI, Zhikharev CN (1996) The generation of Tollmien–Schlichting waves by free-stream turbulence. J Fluid Mech 312:341–371
Department of Mathematics, Imperial College London, 180 Queen's Gate, London, SW7 2BZ, UK
A. I. Ruban
School of Mathematics, University of Manchester, Manchester, M13 9PL, UK
H. Broadley
Correspondence to A. I. Ruban.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Ruban, A.I., Broadley, H. The generation of Tollmien–Schlichting waves by free-stream turbulence in transonic flow. J Eng Math 129, 1 (2021). https://doi.org/10.1007/s10665-021-10138-y
Boundary layer receptivity
Free-stream turbulence
Transonic flow