3. Preparatory computational tests

3.1. Mesh and domain

In our study, we investigated bumps with heights as low as 5 % of the local displacement thickness. For the TS wave of $F = 90 \times 10^{-6}$, a 5 % high bump with length $L=18.2$ centred at ${\textit {Re}} _{\delta ^*} = 950$ ($x_R=320.8$) can increase the TS wave amplitude by 2.5 % relative to that for the smooth plate. To investigate this effect, a compromise between accuracy and computational cost had to be reached. We considered the results independent of mesh if the maximum variation of amplitude were approximately 0.5 %. For this study, the TS wave was excited by wall suction and blowing. The disturbance was centred at ${\textit {Re}} _{\delta ^*} = 700$ ($x_D=174.0$) and had a slot length of $L_D = 18.0$. The meshes ‘0’ to ‘2’ used in the tests are discussed with the help of figure 3. The properties of all meshes are summarised in table 3. For a smooth plate, the grid requirements are not so severe as shown in our previous study (Himeno et al. Reference Himeno, Carvalho and Medeiros2023). The minimum mesh spacing of $\Delta y_{min} = 0.05$ along the $y$-direction is sufficient. The typical mesh spacing in the $x$-direction ($\Delta x_{typ}$) is better defined in terms of the number of nodes per TS wavelength. A value of $\Delta x_{typ}$ corresponding to approximately 13 grid points per TS wavelength is necessary to provide the accuracy required, which led to a value of $1.5$. The ‘mesh 0’ in figure 3(a,b) denotes the mesh that produced independent solutions for the smooth case. For a plate containing a bump, the most critical mesh parameter is the refinement at the bump corners. Figure 3(c) presents the solutions for a monochromatic TS wave interacting with a 5 % high bump conducted with different levels of refinement at the bump region, indicated in figure 3(a,b) as ‘mesh 1’ and ‘mesh 2.’ As discussed in § 2.2, $\Delta x_{min}$ was made equal to $\Delta y_{min}$ also to diminish errors related to the pressure boundary condition at the corners. As shown by the relative amplitude error in figure 3(d), the mesh with $\Delta x_{min} = \Delta y_{min} = 0.01$ satisfies the criterion of a relative error around 0.5 %. Despite being a very small grid element, the highest refinement level region is very localised and the mesh is rapidly stretched along both directions up to the refinement necessary to simulate the smooth plate TS wave evolution.

Figure 3.

Mesh independence analysis for a TS wave interacting with a 5 %$\times \delta ^*$ high bump. The bump is located at $x_R = 320.82$ ($Re_{\delta ^*} = 950$) and stretches along a length of $L = 18.18$, indicated by the shaded area. The figure shows the mesh properties along (a) the x-direction and (b) the y-direction. The plot (c) shows the TS amplitude evolution along the plate for each mesh and (d) the respective absolute change with respect to the finest grid (mesh 2). Zoomed views are also provided on areas indicated by the magenta boxes. Here, $F = 90 \times 10^{-6}$.

Table 3.

Properties of the meshes used in the present simulations. The last column (observations) provides some details of the simulations conducted with each mesh; ‘S’: smooth plate, ‘R’: rough plate, ‘M’: monochromatic disturbance, ‘W’: wave packet disturbance. The last column also includes the sections the meshes are used in the paper.

In the validation test section (§ 3.3), we compare our results with those from Wörner et al. (Reference Wörner, Rist and Wagner2003), who studied the evolution of a TS wave of $F = 49.34 \times 10^{-6}$. For that, meshes ‘3’ and ‘4’ were used. The typical $x$ spacing of $\Delta x _{typ} = 2.8$ is consistent with the criterion of 13 grid points per TS wavelength for that parameter space. The minimum grid sizes in the $x$- and $y$-directions are of the order of 0.01. The value of $\Delta x_{min} = \Delta y_{min}=0.0117$ fits the grid to the authors’ roughness height.

In § 4, the simulations were conducted with meshes ‘5’ and ‘6.’ This section employed wave packet simulations. The typical $x-$ grid used was $\Delta x _{typ} = 1.25$ to satisfy the criterion of 13 grid points per TS wavelength for the lowest frequency of interest in the packet.

In §§ 5 and 6, besides meshes ‘5’ and ‘6’, meshes ‘8’ and ‘9’ were used to account for the higher bump-location Reynolds numbers. A corresponding smooth plate mesh was also used. Figure 4 shows every 4 nodes along both the $x$- and $y$-directions. This plot illustrates how the mesh nodes are distributed along the domain relative to the dimensions of the bump of $h=0.40$ at ${\textit {Re}} _{\delta ^*} = 950$.

Figure 4.

‘Mesh 6’ grid showing one every 4 nodes along both directions. A detail is provided for the region (I) containing the bump of $h=0.40$ at ${\textit {Re}} _{\delta ^*} = 950$ (blue). A detailed view is given of region (II) at the bump leading edge. Buffer zones are indicated by the magenta areas.

Besides the grid refinement, it is also necessary to verify that the results are independent of the computational domain. As will be seen later, this becomes more critical for larger bumps as the Mach number approaches 1. Figure 5 compares the results of mesh ‘6’ and mesh ‘7’ which is twice as long in the vertical direction (see table 3). The figure indicates that mesh ‘6’ is not domain independent as $M_\infty$ and $h$ increase. In the compressibility effect studies (§ 8), mesh ‘7’ was used for simulations with $M_\infty \geq 0.7$ and $h \geq 0.30$. Other than that, mesh ‘6’ was used. Results suggest that simulations with mesh ‘7’ are sufficiently accurate to support the conclusions.

Figure 5.

Domain independence analysis. (a) Amplitude evolution relative to the smooth plate for the case with a bump of $h=0.40$ at ${\textit {Re}} _{\delta ^*} = 950$, $M_\infty = 0.90$ and selected frequencies, considering meshes ‘6’ and ‘7’ with different domains. (b) The amplitude of the dominant TS wave at $x=500.0$ relative to the smooth case for bumps with $h=0.30$ and $h=0.40$, and $M_\infty \geq 0.7$, considering different domains.

3.2. Linearity of TS wave and compressibility issues

A comparison between DNS and linear incompressible parabolised stability equations (PSE) results for the TS wave amplitude evolution along the smooth plate is presented in figure 6 for some selected frequencies. The DNS simulations used a wave packet disturbance. The PSE code was obtained from Juniper, Hanifi & Theofilis (Reference Juniper, Hanifi and Theofilis2014). The amplitude of each mode obtained by the PSE code was matched to the DNS at $x=x_R=320.82$. Two measures of amplitude were used, namely, energy (based on $E = \int _0 ^{y_f} \tfrac {1}{2}( |\widehat {u^\prime }|^2 + |\widehat {v^\prime }|^2)\,\mathrm {d}y$) and the maximum amplitude of the fluctuating streamwise velocity along the wall-normal direction ($|\widehat {u^\prime }|_{max}$). The good agreement confirms that the packet evolution is linear and accurate and that our $M_\infty = 0.1$ simulations are effectively incompressible.

Figure 6.

Amplitude evolution for TS waves of selected frequencies along the smooth plate obtained by the present DNS code at $M_\infty = 0.1$ and by the incompressible linear PSE code. (a) Amplitude is based on the integrated fluctuation energy and (b) based on the maximum streamwise-velocity fluctuation. Here, ${\textit {Re}} _{\delta ^*} = 950$ at the bump location.

3.3. Comparison with literature

Results are compared with those of Wörner et al. (Reference Wörner, Rist and Wagner2003) that investigated the effect of rectangular two-dimensional bumps on the evolution of the TS wave of $F=49.34 \times 10^{-6}$. The reference presented incompressible results, hence we used $M_\infty =0.1$. The Reynolds number based on the conditions of the smooth plate at the roughness centre is ${\textit {Re}} _{\delta ^*} = 1088$. The bumps had a length of $L=20$ and the heights were $h=0.235$ and $h=0.470$. The TS wave of $F=49.34 \times 10^{-6}$ ($\omega = 5.37 \times 10^{-2}$) was included by using the monochromatic excitation disturbance (2.8a,b) with the amplitude of $A_D = 1 \times 10^{-5}$. The excitation slot was centred at $x_D = 254.8$ (${\textit {Re}} _{\delta ^*} = 906$) and extended along the length $L_D = 33.6$.

Figure 7(a) shows the evolution of the TS wave amplitude we obtained for these bumps. Wörner et al. (Reference Wörner, Rist and Wagner2003) report that the bumps of $h=0.235$ and $h=0.470$ were responsible for an amplitude increase of 1.4 and 3.2 compared with the smooth plate. They report that their simulations overestimated the amplitude obtained in their corresponding experiments by 7 %. Figure 7(b) shows the results we found for the far downstream amplification factor relative to smooth and the results from Wörner et al. (Reference Wörner, Rist and Wagner2003). Perfect agreement with their experimental results is obtained.

Figure 7.

(a) Current DNS results of the TS wave ($F=49.34 \times 10^{-6}$) amplitude evolution for selected cases investigated by Wörner et al. (Reference Wörner, Rist and Wagner2003). (b) Far downstream (${\textit {Re}} _{\delta ^*} = 1400$) amplification factor relative to the smooth case compared with that from Wörner et al. (Reference Wörner, Rist and Wagner2003). Bumps are centred at ${\textit {Re}} _{\delta ^*} = 1088$ and have a length of $L = 20$. $M_{\infty } = 0.1$.

4. Base flow and TS wave evolution

4.1. Base flow

In this section, we investigate how the incompressible boundary layer is deformed by the inclusion of a bump. The bumps are located at ${\textit {Re}}_{\delta ^*} = 950$ ($x_R=320.8$) for a physical domain extending up to ${\textit {Re}} _{\delta ^*} = 1210$ ($x=520.0$). This value of ${\textit {Re}}_{\delta ^*}$ is the same as that of de Paula, Würz & Medeiros (Reference de Paula, Würz and Medeiros2008); de Paula et al. (Reference de Paula, Würz, Mendonça and Medeiros2017). The base flow shows important features that can affect the evolution of TS waves. Figure 8 presents the streamlines in the bump region. A region of recirculating flow (red streamlines) is formed immediately behind the bump. For each case, the reattachment point (${\partial {u}}/{\partial {y}}| _{wall} = 0$) is indicated by a red marker. Even for $h=0.05$ a very small recirculating flow region is formed, however, as it reattaches on the bump lateral edge, for this case we considered $x_{reattach} = x_{R_{TE}}= x_R + L/2$. Figure 9 shows how the recirculation length, $L_{recirc}$ scales with $h$. For heights above $h=0.20$ a parabola fits the data very well.

Figure 8.

Streamlines for cases with bumps of different heights. The red streamlines indicate the recirculating flow region that reattaches on the plate at $x_{reattach}$ marked with red symbols. In all cases, the bumps are located at ${\textit {Re}} _{\delta ^*} = 950$ and have $L=18.2$. Here, $M_{\infty } = 0.1$.

Figure 9.

Length of the recirculating flow region vs bump height squared. The dashed line denotes a quadratic function fitting the length given by $h=0.30$ and $h=0.40$. Bumps centred at ${\textit {Re}} _{\delta ^*} = 950$ with length $L=18.2$. Here, $M_{\infty } = 0.1$.

Figure 10(a,b) shows the deviation caused by the bumps on both the streamwise and wall-normal velocity profiles, at selected $x$-stations downstream of the bump. Hereinafter, we refer to this as profile deviation. For each station, $y$ is normalised by the local Blasius displacement thickness, $\delta ^* _B$. The maximum deviations for $h=0.40$ are indicated for reference. For both velocity components, the deviation essentially decreases in magnitude with the maximum located progressively farther from the wall. At each $x$-station, the profile deviations are very similar for all bumps, except for the recirculating flow region. This exception is noted particularly for the $v$-profile deviation at the two initial stations for $h=0.40$.

Figure 10.

Velocity deviation profiles caused by bumps with different heights at selected stations in the bump wake. (a) Streamwise component ($u-u_{B}$) and (b) wall-normal component ($v-v_B$). Bumps are located at ${\textit {Re}} _{\delta ^*} = 950$ and have $L=18.2$. Here, $M_{\infty } = 0.1$.

The maximum deviation of $(u-u_B)$ is fairly large for the bumps. As a reference, in order to cause a deviation of $-0.175$ ($h=0.40$) the shape factor ($m$) of the Falkner–Skan profile family would have to be $m=-0.06$. Moreover, the maximum deviation of the Falkner–Skan profiles occurs at approximately $y/{\delta ^* _B}=1.2$, while for the bump wake profiles in figure 10(a), they occur, in general, much closer to the wall. This is important because the closer to the wall the deviation is, the stronger its impact on the stability of the boundary layer.

The distributions of the displacement thickness ($\delta ^*$), momentum thickness ($\theta$) and shape factor ($H$) are plotted in figure 11(a,b). The deviation from the Blasius profile increases with $h$. The value of $\delta ^*$ presents a maximum deviation just upstream of the bump trailing edge. The $\theta$ deviation oscillates in the bump region. There is a global maximum just upstream of the bump leading edge, a global minimum right at the bump leading edge and a local maximum downstream of the bump. In spite of that, these plots indicate that the smooth case Blasius profile is asymptotically recovered downstream. Figure 12(a,b) shows the evolution of $H-H_B$ normalised by $h$ as a function of the distance from the leading edge and from the reattachment point for each bump, respectively. Along the bump and the recovery region, the shape factor deviation, $H-H_B$, scales linearly only up to $h=0.10$ as shown by the collapse of the curves in red and blue. Above it, the shape factor deviation increases nonlinearly with the bump height. A higher bump naturally produces a higher deviation. However, higher bumps also are exposed to a higher flow speed, which should further distort the profiles. We believe this is a possible explanation for the nonlinear scaling observed. Along the recovery region, the deviation evolves similarly in particular if each evolution is shifted by the respective reattachment point. nevertheless, the recovery region extension is not substantially affected by the bump height. For all heights $(H-H_B)/h < 0.01$ occurs nearly at the same location, $x \approx 420.0$ ($\approx 90.0$ relative to the bump).

Figure 11.

Bump effect on the boundary layer (a) displacement and momentum thicknesses, $\delta ^*$ and $\theta$, respectively, and (b) shape factor ($H$). Bumps at ${\textit {Re}} _{\delta ^*} = 950$ with $L=18.2$. Here, $M_{\infty } = 0.1$.

Figure 12.

Evolution of the shape factor deviation normalised by $h$ as a function of distance from (a) the plate leading edge and (b) the reattachment location ($x-x_{reattach}$). Bumps at ${\textit {Re}} _{\delta ^*} = 950$ with $L=18.2$. Here, $M_{\infty } = 0.1$.

4.2. The TS wave

Next, we evaluate the impact of the bump on the evolution of a TS wave. The simulations used wave packet disturbances, but in this section the analysis was performed for a single frequency. This is adequate since our analysis is restricted to the linear regime, as shown in § 3.2. The following sections analyse the full packet. The disturbance was excited along a slot centred at ${\textit {Re}} _{\delta ^*} = 700$ ($x_D=173.7$). In the present section, the TS wave had the frequency $F = 90 \times 10^{-6}$, which roughly represents the dominant TS wave at the end of the physical domain considered. The bump length is $L = 18.2$. Figure 13 shows the evolution of the TS wave maximum amplitude. As the TS wave travels, it destabilises along the region upstream of the bump, then it stabilises along the bump element and, finally, destabilises in the wake. These three stages of the TS wave evolution under the influence of the bump were also observed by Nayfeh et al. (Reference Nayfeh, Ragab and Al-Maaitah1988) and Wörner et al. (Reference Wörner, Rist and Wagner2003). The fact that the TS wave reduces in amplitude along the bump should not be taken as surprising. The TS instability depends on a particular arrangement of velocity profiles close to the wall. This balance is affected by the profile distortions caused by the bump. It would be rather surprising if, instead, a large increase in the TS wave was seen on the bump. Figure 13 shows that, despite the relatively large fluctuations along the bump element, up to the bump trailing edge, the net impact is negligible. This is shown by the fact that, at the trailing edge of the bump, the curves for all the bump heights converged to the smooth plate curve. The dominant and enduring effect occurs in the wake region where the TS wave is massively amplified as the bump height increases. Far downstream of the bump, the TS wave amplitude evolves almost parallel to the smooth case, which indicates the growth rate associated with the Blasius profile is recovered. The TS wave recovery region (as coined by Klebanoff & Tidstrom Reference Klebanoff and Tidstrom1972) ends around $x \approx 420.0$ (figure 12), at the same location where the distortion caused by the bump on the Blasius profile becomes negligible.

Figure 13.

Amplitude evolution of the TS wave of $F = 90 \times 10^{-6}$ for different bump heights (normalisation considers the TS amplitude at $x=200.0$). The bumps are located at ${\textit {Re}} _{\delta ^*} = 950$ and have $L=18.2$. Here, $M_{\infty } = 0.1$.

The evolution of the growth rate deviation from Blasius ($\sigma - \sigma _B$, with $\sigma = (\textrm {d}A_{TS}/{\textrm {d}\kern 0.06em x})/A_{TS}$), normalised by the bump height is shown in figure 14. This plot shows that $\sigma - \sigma _B$ scales linearly with the bump height up to $h=0.10$. Since the deviation is normalised by $h$ the linear scaling up to $h=0.1$ is shown by the collapse of red and blue curves. Above $h=0.10$, the growth rate deviation varies nonlinearly with $h$, most of the nonlinear activity taking place very close to the bump. Consistent with the analysis above (figures 10–12), in this $h-$ normalised scale, the location where the smooth plate behaviour is recovered is around $x \approx 420.0$, independent of bump height.

Figure 14.

Evolution of the growth rate deviation normalised by $h$ for the TS wave of $F = 90 \times 10^{-6}$. The bumps are located at ${\textit {Re}} _{\delta ^*} = 950$ and have $L=18.2$. Here, $M_{\infty } = 0.1$.

The fact that the growth rate deviation scales linearly with bump height is noteworthy. This deviation is associated with the profile deviation, which up to $h=0.10$, scales linearly with $h$. So the results in figure 14 are consistent. In turn, where the growth rate deviation scales linearly with $h$, the TS wave amplification factor associated with the bump should scale exponentially. Figure 15 presents the amplification factor vs $h$ at $x=490.0$ (indicated by the red markers in figure 13) but any $x$-station downstream of $x=420.0$ (end of the recovery region) provides the same results. Consistent with the arguments above, up to $h=0.10$ the amplification factor associated with the bump should scale as $\exp {(c_1 \times h)}$. This expression was fit to the two smaller bump results, including the smooth case, and provided $c_1=0.5$. According to figure 14, above $h=0.10$ the amplification factor should scale super-exponentially. As shown in figure 15, a good fit was obtained for $\exp {(6.0 \times h^2)}$.

Figure 15.

Amplification factor vs (a) $h$ and (b) $h^2$ at the $x$-station $x=490.0$ indicated in figure 13. The lines are best-fit curves. The bumps are located at ${\textit {Re}} _{\delta ^*} = 950$ and have $L=18.2$. Here, $M_{\infty } = 0.1$.

5. Effect of Reynolds number on the base flow

This section investigates the effect of the Reynolds number on the base flow deviation caused by a bump. In addition to the ${\textit {Re}} _{\delta ^*} = 950$ case, we considered bumps placed at ${\textit {Re}} _{\delta ^*} = 1500$ and ${\textit {Re}} _{\delta ^*} = 2000$. The bump heights were fixed to $h=0.40$ relative to the local displacement thickness of each Reynolds number. In all cases, the non-dimensional bump length was fixed to $L=18.2$, which is the same as for ${\textit {Re}} _{\delta ^*} = 950$. The Mach number was set to $M_\infty = 0.1$ (incompressible).

Figure 16 compares the streamlines for the bump placed at different Reynolds numbers. Consistent with the observations of Klebanoff & Tidstrom (Reference Klebanoff and Tidstrom1972), the length of the recirculating flow region increases with ${\textit {Re}} _{\delta ^*}$. For a fixed bump height, the recirculation region length roughly scales linearly with Reynolds number, figure 17.

Figure 16.

Streamlines for cases with bump $h=0.40$ placed at different ${\textit {Re}} _{\delta ^*}$. The red streamlines indicate the recirculating flow region that reattaches to the plate at $x_{reattach}$ marked with red symbols. Here, $M_{\infty } = 0.1$ and $L=18.2$.

Figure 17.

Length of the recirculating flow region as a function of ${\textit {Re}} _{\delta ^*}$, for the bump with $h=0.40$ located at ${\textit {Re}} _{\delta ^*} = 950$. Here, $M_{\infty } = 0.1$ and $L=18.2$.

Figure 18 shows the effect of the Reynolds number on the shape factor ($H$) evolution. The curves indicate that the Reynolds number influences mainly the bump wake. It does not only increase the extent of the wake but also the magnitude of the shape factor deviation. The deviation velocity profiles are shown in figure 19. The Reynolds number affects the deviation magnitude, but not its distribution. It is known that the Reynolds number does not affect the shape of the laminar wake (Schlichting & Gersten Reference Schlichting and Gersten2000). Our observations of the bump wake are consistent with these results. As the bumps have the same height relative to the local Blasius profile $\delta ^* _B$, it is expected that the Reynolds number would not affect the shape of the laminar wake. The process by which the deviation profile vanishes is governed by diffusion. As the Reynolds number increases, the diffusion reduces relative to the convection in the boundary layer. Therefore, the wake requires a longer extension to return to the Blasius profile as the Reynolds number increases. This explains part of the destabilising influence of the Reynolds number on the effect of a bump on a TS wave evolution.

Figure 18.

Effect of the bump ${\textit {Re}} _{\delta ^*}$ on the boundary layer shape factor ($H$), for $h=0.40$. Here, $M_{\infty } = 0.1$ and $L=18.2$.

Figure 19.

Streamwise-velocity deviation profiles, $u-u_{B}$, caused by a bump of $h=0.40$ at different ${\textit {Re}} _{\delta ^*}$. Here, $M_{\infty } = 0.1$ and $L=18.2$.

6. Effect of frequency and Reynolds number on the TS wave–bump interaction

Here, we extend the analysis to account for the effect of Reynolds number on the TS wave–bump interaction. As will be seen, it is intertwined with an effect of TS wave frequency. Figure 20 shows contours of growth rate as a function TS wave frequency along the streamwise direction, for bumps located at ${\textit {Re}} _{\delta ^*}=950$ with heights $h=0.10$, $0.20$ and $0.40$. These results were obtained with simulations of the evolution of wave packets. As observed in figure 20, the frequency range of the packets ensured that all relevant modes were excited. These plots extend in frequency the analysis from § 4.2. We observe two distinct regions characterised by an increase in the growth rate that was not clearly apparent when focusing solely on the frequency of $F=90\times 10^{-6}$ (see figure 14). The first region is that corresponding to the recirculation flow (see figure 16) where the growth rate is notably high, albeit extending only for a short distance from the bump. The second region is in the far wake where the increase in the growth rate is comparatively lower, but persists over a much longer distance. As the latter exerts a more sustained influence, the effect along the far wake tends to be dominant. We are not confident on the stabilising effect observed downstream of $x=420.0$ around the lower branch, more clearly noticed for $h=0.40$. The evolution of TS waves with such low frequencies becomes progressively more compromised because they may not have been fully established before reaching the bump. In any case, these frequencies are so low that do not contribute to the $N$-factor curves discussed below and hence do not affect the conclusions of the paper. The contour plots in figure 14 are consistent with those presented by Xu et al. (Reference Xu, Mughal, Gowree, Atkin and Sherwin2017), Thomas, Mughal & Ashworth (Reference Thomas, Mughal and Ashworth2017) and Thomas et al. (Reference Thomas, Mughal, Roland, Ashworth and Martinez-Cava2018). For instance, their results show an increase in the growth rate right at an indentation while ours show a reduction in the growth rate right at the bump.

Figure 20.

Contours of growth rate for bumps with $h=$ (a) $0.10$, (b) $0.20$ and (c) $0.40$. Black-dashed lines denote iso-levels of $\sigma = [0.001;0.004;0.006]$ for a smooth plate case, the dash-dotted lines are the neutral curves. Here, ${\textit {Re}} _{\delta ^*} = 950$, $M_{\infty } = 0.1$ and $L=18.2$.

The contours in figure 20 also show that along the bump wake the impact on the growth rate is much stronger for the TS waves around the upper branch. This is what leads to the upward shift in the frequency of the dominant TS wave, observed in early studies (Klebanoff & Tidstrom Reference Klebanoff and Tidstrom1972; Nayfeh et al. Reference Nayfeh, Ragab and Al-Maaitah1988; Masad & Iyer Reference Masad and Iyer1994). Recall that the base flow deviation along the wake profiles (see figure 10a) also resembles that caused by the presence of an adverse pressure gradient. The adverse pressure gradient extends the unstable region of the instability diagram (Schlichting & Gersten Reference Schlichting and Gersten2000) mostly by moving the upper branch upwards. We refer to this phenomenon as the upper branch effect. Clearly, besides the effect of frequency, this also incorporates an effect of the Reynolds number at the bump location.

To better substantiate the effect of the Reynolds number we run simulations with a bump of $h=0.40$ placed on two other locations (${\textit {Re}}_{\delta ^*}=1500$ and $2000$). The results of the $n$- and the $N$-factor evolution are shown in figure 21(a). Here, the $n$-factor corresponds to a TS wave amplification for a single frequency, $n(x)=\ln (A_{TS}(x)/A_{TS,0})$, where $A_{TS,0}$ corresponds to the TS wave amplitude at the lower branch. The $N$-factor corresponds to the envelope of the $n$-factor evolution curves for all frequencies at all $x$-stations. Figure 21(b) shows the change on the $n$- and the $N$-factors for the bump cases relative to that from the smooth one (black lines in figure 21a). Both pictures display clearly that the higher the Reynolds number the larger the impact of a given bump. This result incorporates both the effect of the Reynolds number on the base flow and the upper branch effect.

Figure 21.

(a) The $n$-factor (thin lines) and $N$-factor (thick lines) evolution for a bump of $h=0.40$ placed at ${\textit {Re}}_{\delta ^*}=950$, ${\textit {Re}}_{\delta ^*}=1500$ and ${\textit {Re}}_{\delta ^*}=2000$; (b) $n$- and $N$-factor change relative to those from the smooth plate. For selected TS frequencies the $n$-factor curves are highlighted and indicated by different line styles. For these modes, results from PSE calculations (black lines) are also shown for reference. Here, $M_{\infty } = 0.1$ and $L = 18.2$ in all cases.

Figure 21 also illustrates another interesting aspect associated with the upper branch effect that affects the $\Delta N$ curves. The frequency of the TS wave at the transition location reduces as the transition moves downstream independently of the presence of the bump. However, the lower frequencies are less affected by the bump, as also seen by Klebanoff & Tidstrom (Reference Klebanoff and Tidstrom1972) and Wang & Gaster (Reference Wang and Gaster2005). Consequently, while the $\Delta n$ curves tend asymptotically to a constant value, the $\Delta N$ curves decrease even downstream of the recovery region. Because of that, the effective $\Delta N$ factor depends also on the magnitude of the TS wave content in the boundary layer. For example, from figure 21, if the environment is such that $N_{cr}=5$, for a bump at ${\textit {Re}}_{\delta ^*}=1500$ (red curve), the transition would occur around ${\textit {Re}}_{\delta ^*} \approx 1600$, with a $\Delta N = 1.8$. For the same bump in a reduced noise level environment such that $N_{cr}=7.5$, the transition would occur at ${\textit {Re}}_{\delta ^*}=2400$ with $\Delta N \approx 1.1$. An almost 40 % reduction in the $\Delta N$.

As an example, the overall effect of the Reynolds number can be quantified as follows. Suppose that transition occurs at ${\textit {Re}}_{\delta ^*}=2150$ for the bump at ${\textit {Re}}_{\delta ^*}=2000$ (yellow curve in figure 21). In this case, $\Delta N=2.3$ and $N_{cr}=7.6$. On the other hand, for the bump at ${\textit {Re}}_{\delta ^*}=1500$, the transition would occur at ${\textit {Re}}_{\delta ^*}=2400$ with $\Delta N = 1.1$. For the bump located at ${\textit {Re}}_{\delta ^*}=950$, $\Delta N$ would be approximately 0.5 or less, with an almost negligible impact on transition. A twofold variation in Reynolds number could lead to more than fourfold variation in $\Delta N$. Recall that these variations in $\Delta N$ incorporate both the effect of the Reynolds number on the base flow and the upper branch effect. We quantified these effects as ratios between two values of $\Delta N$ because these ratios should remain similar for bumps of other heights at these locations. These ratios are very significant and neglecting these effects should lead to inaccuracy in transition prediction.

Figure 21 also shows that the bump causes a frequency shift. It does not represent the difference in the TS frequency at the transition location for a bump case relative that of the smooth one. The frequency shift is on the dominant TS wave at each $x$-station (Himeno et al. Reference Himeno, Carvalho and Medeiros2023). It is associated with the difference between the dotted blue line, the dashed red line and the dashed-dotted yellow line and their respective envelopes. This shift depends on the bump height and Reynolds number, but in the current parameter space, it represents an error in $\Delta N$ below 0.1, which is not relevant for predictions.

There are some studies in the literature that present results for different bump-location Reynolds numbers. However, only two studies (Klebanoff & Tidstrom Reference Klebanoff and Tidstrom1972; Masad & Iyer Reference Masad and Iyer1994) were found that gave enough detail for the analysis below. Klebanoff & Tidstrom (Reference Klebanoff and Tidstrom1972) is a wind tunnel experiment that measured the impact of circular-section bumps (wires) on the transition location under natural conditions. Masad & Iyer (Reference Masad and Iyer1994) give a theoretical study of the impact of smooth-shaped bumps (humps) using the IBL method and linear stability analysis to obtain the amplification of the TS waves. In Masad & Iyer (Reference Masad and Iyer1994), the transition location is based on the $e^N$ method with $N_{cr} = 9$. The use of $N_{cr} = 9$ could be questioned, but the conclusions of the analysis do not depend on the value of $N_{cr}$. These data are interesting also because they extend the parameter space towards bumps substantially larger than ours. Klebanoff & Tidstrom (Reference Klebanoff and Tidstrom1972) used bumps with $h=0.72$ and $h=0.77$ and for the selected cases from Masad & Iyer (Reference Masad and Iyer1994), the heights are between $h=0.58$ and $h=1.17$. These studies also extend the analysis to include other methods, namely, experiment and theory, while ours is numerical.

The results are plotted in figure 22, on top of the smooth plate instability diagram. The squares represent the Klebanoff & Tidstrom (Reference Klebanoff and Tidstrom1972) results and the circles those of Masad & Iyer (Reference Masad and Iyer1994). The vertical green line indicates the transition Reynolds number for the smooth case from Masad & Iyer (Reference Masad and Iyer1994). The red markers indicate the Reynolds number for the transition location ($N=9$) with the bump. The blue markers indicate the Reynolds number at the bump location. The distance between the red marker and the green line gives the movement of the transition location caused by the bump. In the plot, for all markers the frequency was the dominant one at the transition location for the case with the bump. The bump shifts the dominant frequency upwards. The larger the transition point movement, the higher the frequency shift. There is a bump, located shortly above ${\textit {Re}} _{\delta ^*}=2000$, that causes the largest transition location movement and higher-frequency shift. For bumps upstream of this, the impact of the bump is given by the distance between the green line and the red markers, the displacement of the transition point. For bumps downstream of this location, the transition point movement is necessarily smaller than the maximum, because the bump will not move the transition to a point upstream of the bump. The transition point is then so close to the bump that $N_{cr}=9$ is reached upstream of the point of its maximum $\Delta N$. In other words, the bump does contribute to transition with its full potential. For bumps at these positions, the magnitude of the bump effect is not measured by the distance between the smooth plate transition point (green line) and the bump case transition point (red markers). It is given by the distance between the bump position (blue marker) and the transition point for the bump case (red markers). Downstream of the critical position, the higher the bump impact the closer to the bump the transition occurs. There is clearly a bump position that moves the transition point the most, which is that for which the critical $N$ is reached right where $\Delta N$ is maximum. For reference, we called it the critical position. This position would depend on the background noise and probably on the bump shape.

Figure 22.

Data extracted from Klebanoff & Tidstrom (Reference Klebanoff and Tidstrom1972) (squares) and Masad & Iyer (Reference Masad and Iyer1994) (circles) plotted over the instability diagram. The Reynolds numbers at the bump are indicated by light- ($h<1.0$) and dark-blue ($h \geq 1.0$) symbols. The Reynolds number at the transition locations are shown by the red symbols. For each case, the frequency corresponds to that which is critical at the transition location.

In figure 22, for each of the two references (Klebanoff & Tidstrom Reference Klebanoff and Tidstrom1972; Masad & Iyer Reference Masad and Iyer1994), the size of the markers (square or circle) represents the relative magnitude of the bump height. The data were not meant to investigate exclusively the effect of either bump height or Reynolds number and for every bump-location Reynolds number, a different height relative to the local displacement thickness was used. Yet, useful information on the different effects can be obtained. By the way in which these tests were performed, in general, the larger bumps were placed at a higher Reynolds number, two effects that tend to produce a larger impact on transition, either by causing a larger transition movement (for bump upstream of the critical position) or by moving transition closer to the bump (for bumps downstream of the critical position). However, there are a few roughness pairs (indicated by the numbers 1, 2 and 3) for which a smaller bump was tested downstream of a larger one and had a higher impact. The bump impact is assessed by either comparing the length of the transition movement or by comparing the transition point proximity with the bump. There is also a case (indicated by 4) in which a smaller bump had almost the same impact as a larger one. In general the bump height has a higher impact, but for these 4 cases, the effect of Reynolds number is manifested. This phenomenon is apparent downstream of the maximum growth rates of the TS wave on a smooth plate, consistent with the arguments raised above about the effect of proximity to the upper branch. In fact, in the figure, it is noticeable that, under the influence of the bump, the transition occurs shortly upstream of the upper branch. This is also possibly an underlying manifestation of the upper branch effect.

In summary, the upper branch effect dictates that the bump impact increases the closer the bump is to transition. Under these circumstances, the TS wave leading to transition interacts with the bump closer to its maximum $\Delta N$ location. On the other hand, if the bump is located far upstream of transition, the TS wave that leads to transition interacts with the bump far from the point where it crosses the upper branch and the bump has a smaller impact on the TS wave. It is important to add that the results analysed in this section include both the effect of the Reynolds number on the base flow as well as the Reynolds number effect associated with the proximity to the upper branch. Both of them dictate that the bump effect is larger as the bump moves downstream, closer to the transition location. The upper branch effect also dictates that for a fixed bump location, if transition moves upstream (caused for example by an increase in free stream turbulence) the bump effect increases.

7. The combined effects of the bump height and local Reynolds number

The impact of a bump on the transition location has been evaluated in terms of several dimensionless parameters (Tani et al. Reference Tani, Hama and Mituisi1940; Fage & Preston Reference Fage and Preston1941; Dryden Reference Dryden1953; Carmichael Reference Carmichael1959). One parameter widely used is a Reynolds number defined as ${\textit {Re}}_{hh} = u_h ^{\star } h^{\star } / \nu ^{\star }$, where $u_h ^\star$ denotes the streamwise velocity evaluated at the bump height for the smooth surface boundary layer. Considering the smooth case is a Blasius flow, for a bump with non-dimensional height $h = h^\star / {\delta ^*}^\star$ placed at ${\textit {Re}}_{\delta ^*}$, we can write ${\textit {Re}}_{hh} = 0.571~ {\textit {Re}} _{\delta ^*} ~ h^2$ (Dryden Reference Dryden1953; Morkovin Reference Morkovin1990). As the ${\textit {Re}}_{hh}$ parameter encompasses both the $h$ and ${\textit {Re}}_{\delta ^*}$ effects, this parameter was not convenient to discuss the individual contributions of the bump height and of the local Reynolds number. However, the combined effects analysed in terms of the ${\textit {Re}}_{hh}$ give further valuable insight into understanding of the phenomenon.

First, recall some experimental results from Tani (Reference Tani1961) that investigated the impact of circular-section bumps (wires) on the transition location of a flat plate boundary layer. For each bump with fixed dimensional height ($h^\star$) and location ($x_R ^\star$), Tani (Reference Tani1961) measured the transition location ($x_{tr}^\star$) as the free-stream velocity ($U_\infty ^\star$) increased. As expected, for a fixed dimensional bump case, $x_{tr}^\star$ moves towards $x_{R}^\star$ as $U_\infty ^\star$ increases. Tani (Reference Tani1961) also shows that, for each run, there exists a $U_\infty ^\star$ such that the Reynolds number at the transition location, ${\textit {Re}} _{x_{tr}} = U_\infty ^\star x_{tr}^\star / \nu ^\star$, reaches a minimum value, $({\textit {Re}} _{x_{tr}})_{min}$. For this particular condition, the transition does not occur at the roughness, but at some distance downstream. Tani (Reference Tani1961) also shows that the $({\textit {Re}} _{x_{tr}})_{min}$ condition was not related to a constant ${\textit {Re}}_{hh}$, as previously stated by Dryden (Reference Dryden1953), but varies along a range of ${\textit {Re}}_{hh}$.

The condition which causes the lowest ${\textit {Re}} _{x_{tr}}$ corresponds to the maximum bump impact relative to the smooth surface. Then, in the context of the $e^N$-method, such a $({\textit {Re}} _{x_{tr}})_{min}$ would correspond to a maximum $\Delta N$ condition, $(\Delta N)_{max}$. As shown in figure 21(b), the $\Delta N$ curves reach a maximum value not far downstream of the bump, consistent with the Tani (Reference Tani1961) results.

Figure 23 shows $(\Delta N)_{max}$ vs ${\textit {Re}} _{hh}$ for the cases presented in §§ 4–6. The value of $(\Delta N)_{max}$ scales almost linearly with ${\textit {Re}}_{hh}$ as indicated by the dashed line. This linear relation may indicate why $Re_{hh}$ was taken as a good parameter to assess the effect of a bump on transition. The threshold value of ${\textit {Re}} _{hh} = 25$, often assumed as the upper limit for the effectively smooth surface (Morkovin Reference Morkovin1990), is indicated for reference. We note that there is not a ${\textit {Re}} _{hh}$ below which the $(\Delta N)_{max}$ is actually null. Nevertheless, for ${\textit {Re}} _{hh} < 25$, the plot indicates a $(\Delta N)_{max}<0.3$, which is likely to be within the uncertainty of measurements of the transition location in wind tunnels (Crouch & Kosorygin Reference Crouch and Kosorygin2020), hence there is consistency in the results.

Figure 23.

Value of $(\Delta N)_{max}$ vs the ${\textit {Re}}_{hh}$ for bumps of different $h$ and placed at different ${\textit {Re}}_{\delta ^*}$. The blue markers correspond to bumps placed ${\textit {Re}} _{\delta ^*} = 950$, the red marker at ${\textit {Re}} _{\delta ^*} =1500$ and the yellow marker at ${\textit {Re}} _{\delta ^*} =2000$. For all cases, $M_{\infty } = 0.1$ and $L = 18.2$.

Figure 24 shows the recirculation length vs the ${\textit {Re}}_{hh}$ for the same cases shown in figure 23. The value of $L_{recirc.}$ also scales linearly with ${\textit {Re}}_{hh}$, as indicated by the dashed line. In fact, since ${\textit {Re}} _{hh} = 0.571Re_{\delta ^*}h^2$, such agreement was expected from figures 9 and 17, where we showed that $L_{recirc.}$ varies almost quadratically with $h$ and linearly with ${\textit {Re}} _{\delta ^*}$. Moreover, the results in figures 23 and 24 also suggest that $(\Delta N)_{max}$ is associated with the recirculating bubble in the bump wake.

Figure 24.

Value of $L_{recirc.}$ vs ${\textit {Re}}_{hh}$ for bumps of different $h$ and placed at different ${\textit {Re}}_{\delta ^*}$. The blue markers correspond to bumps placed ${\textit {Re}} _{\delta ^*} = 950$ the red marker at ${\textit {Re}} _{\delta ^*} =1500$ and the yellow marker at ${\textit {Re}} _{\delta ^*} =2000$. For all cases, $M_{\infty } = 0.1$ and $L = 18.2$.

8. The destabilising effect close to sonic speeds

In this section, we investigate the effect of compressibility for Mach numbers from $0.1$ to $0.9$, for bumps placed at ${\textit {Re}} _{\delta ^*} = 950$. In all cases, the TS waves were excited by the wave packet described by (2.9) with the parameters $A_{0}$, $\omega _0$ and $n_k$ shown in table 2. For all cases discussed in this section, the simulation set-up was the same as that for the $M_\infty =0.1$ cases. Since the compressibility thickens the boundary layer, the effective non-dimensional bump heights ($h = h^\star / {\delta ^*}^\star$) reduced with $M_\infty$ compared with their nominal values ($h = 0.05$, $0.10$, $0.20$, $0.30$ and $0.40$) defined for $M_\infty = 0.1$. This reduction is generally small, but as it becomes significant for the higher heights and higher Mach numbers, the effective height will be considered.

The disturbance slot was also the same as that for $M_\infty =0.1$ cases. Hence, the effective ${\textit {Re}} _{\delta ^*}$ at the disturbance position was different from the nominal one. However, this variation does not affect the conclusions because only the packet evolution far from excitation is relevant for the current analysis. In the far field, owing to the large amplification rates of the boundary layer, small variations in the excitation position do not significantly affect the packet evolution in the linear regime.

Figure 25(a) shows, for selected nominal bump heights, the effect of $M_\infty$ on the amplitude of the dominant TS wave at $x = 500.0$. The corresponding frequencies are shown in figure 25(b). The results are qualitatively similar to those of Masad & Iyer (Reference Masad and Iyer1994) for both the amplitude and frequency for $M_\infty \leq 0.8$, the Mach number range they investigated. The results show a stabilising effect of $M_{\infty }$.

Figure 25.

Effect of Mach number on the amplitude and frequency of the dominant TS wave at $x=500.0$. Bumps at ${\textit {Re}} _{\delta ^*} = 950$ ($x_R = 320.8$) with length $L=18.2$.

As anticipated, we have to consider that the boundary layer thickness increases with $M_{\infty }$. The results presented in figure 25 are actually for lower values of bump height non-dimensionalised by the smooth plate boundary layer displacement thickness at the bump location, $\delta ^* _{h=0}(x_R)$. The nominal height results suggest an apparent further stabilising effect of the Mach number. Figure 26 shows the variation of $\Delta n$ with respect to the effective bump height, $h=h^\star / {{\delta ^* _{h=0}}^\star (x_R)}$. Note that $h$ reduces as $M_\infty$ increases. For each Mach number, $\Delta n$ was calculated for the TS wave that is dominant at $x=500.0$ for the smooth plate, i.e. for the frequencies shown by the black curve in figure 25(b). The $\Delta n$ curves in figure 26 essentially collapse up to $M_\infty =0.8$, but for $M_\infty =0.9$, the curve diverges substantially from the others. The value of $h$ reduces by approximately 5 %. As occurs in the incompressible regime (see figure 15), $\Delta n$ varies almost quadratically with $h$, as indicated by the best-fit parabolas shown in grey in figure 26. We then used the quadratic fitting curves to estimate the $\Delta n$ for the effective bump heights of $h=0.05$, $0.10$, $0.20$, $0.30$ and $0.40$ at every $M_\infty$. The resulting curves are shown in figure 27.

Figure 26.

The $\Delta n$-factor vs the effective bump height ($h^\star / {\delta ^* _{h=0}}^\star (x_R)$) for each simulated Mach number and plate temperature condition. The colours of symbols indicate the nominal bump height. For each $M_\infty$, the TS wave is that which is dominant at $x=500.0$ considering the smooth case. Bumps at ${\textit {Re}} _{\delta ^*} = 950$ ($x_R = 320.8$) with length $L=18.2$.

Figure 27.

The $\Delta n$-factor caused by bumps of varying height vs the Mach number. For each $M_\infty$, the TS wave is that which is dominant at $x=500.0$ considering the smooth case. Bumps at ${\textit {Re}} _{\delta ^*} = 950$ ($x_R = 320.8$) with length $L=18.2$.

Figure 27 presents results normalised by the smooth plate case. This removes the strong stabilising effect of $M_\infty$ on the TS wave smooth plate evolution noticed in figure 25 and focuses on the influence of compressibility on the effect of a bump on a TS wave evolution. Note that, in figure 27, the curves correspond to the effective bump height since we resorted to the curves in figure 26. This more precise and clearer view indicates that the stabilising effect of compressibility prevails only up to $M_\infty = 0.7$. Above it, compressibility becomes progressively more destabilising, such that the magnitude of the bump impact at $M_\infty = 0.9$ is higher than that at $M_\infty = 0.1$.

In order to investigate the origin of the destabilising effect of compressibility we first investigated its impact on the base flow. Figure 28 presents the evolution of $\delta ^*$, $\theta$ and $H$, for selected bump heights and Mach numbers. The curves are normalised by the effective $h$. The maximum deviation of $H$ occurs on top of the bump and is not substantially affected by Mach number. Upstream of the bump, the curves diverge, but as discussed previously, this region does not significantly affect the overall impact of the bump on the TS wave evolution. Downstream of the bump, for $M_\infty =0.8$ and especially for $M_\infty =0.9$, the shape factor deviation extends over a longer region. The displacement thickness deviation also stretches further downstream with $M_\infty$, while the momentum thickness deviation is less sensitive to $M_\infty$.

Figure 28.

Deviation caused by selected bumps relative to the smooth case in the distribution of the (a) displacement thickness, (b) momentum thickness and (c) shape factor. Results are normalised by the effective dimensionless bump height, $h_e = h / \delta ^* _{h=0}(x_R)$. Bumps at ${\textit {Re}} _{\delta ^*} = 950$ ($x_R = 320.8$) with length $L=18.2$.

The pressure distribution, $C_p$, is shown in figure 29(a). The values were taken at the iso-velocity curve of $u/U_{\infty } = 99\,\%$. The magnitude of the ${C_p}_{u_{99}}$ distribution increases with Mach number but, on the bump's top, the curves collapse fairly well if the Prandtl–Glauert similarity rule is applied, figure 29(b). Recall that the Prandtl–Glauert rules depend on the aspect ratio of the body ($h^\star / L_R^\star$) and not on the boundary layer thickness. Downstream of the bump, the collapse is not as good. Clearly, the higher pressure gradients observed at higher Mach numbers retard the reattachment point. This is consistent with and potentialised by the longer tail of the boundary layer displacement thickness deviation (figure 28a) at the higher $M_{\infty }$.

Figure 29.

(a) Distribution of pressure coefficient, $C_p$, along the iso-velocity curve where $u/U_{\infty }=99\,\%$ for selected bump heights and Mach numbers. (b) Collapse of $C_p$ curves by considering the Prandtl–Glauert similarity rule for compressible flows, $\sqrt {1-M_{\infty }}$. Bumps at ${\textit {Re}} _{\delta ^*} = 950$ ($x_R = 320.8$) with length $L=18.2$.

Figure 30(a) shows that the recirculation length ($L_{recirc.}$) correlates well with the value of ${C_p}_{u_{99}}$ at the bump trailing edge. Figure 30(b) shows that, even though the pressure distribution does not conform to the Prandtl–Glauert rule along the bump wake (figure 29b), $L_{recirc.}$ correlates very well with this similarity rule. The value of $L_{recirc.}$ is not significantly increased in the low-subsonic regime, but the figure suggests that it is likely to be massively larger in the transonic regime. Masad & Iyer (Reference Masad and Iyer1994) had already pointed out that, below $M_{\infty }=0.8$, the destabilising effect caused by the bump in the compressible scenario is connected with the length of the recirculation, but they indicated that this effect was smaller than the stabilising effect of $M_{\infty }$ in the TS waves. Our results confirm this destabilising effect observed by Masad & Iyer (Reference Masad and Iyer1994), but we found that it becomes dominant above $M_{\infty }=0.7$. Our numerical procedures could not cope with the sonic flow, but the results in figure 27 suggest that, towards $M_{\infty }=1.0$, there is a substantial increment in comparison with that at $M_{\infty }=0.7$.

Figure 30.

Length of the recirculating flow caused by the bump of $h = 0.40$ vs (a) ${C_p}_{u_{99}}$ at the bump trailing edge and (b) considering the Prandtl–Glauert similarity rule.

At transonic speeds, the flow in the boundary layer heats up considerably and the use of either isothermal or adiabatic walls needs to be considered (Mengaldo et al. Reference Mengaldo, Kravtsova, Ruban and Sherwin2015). For an adiabatic wall, the higher temperatures thickens the boundary layer and the issue of the effective non-dimensional bump height becomes more relevant. This is shown in figure 26, where the curve for the adiabatic wall at $M_{\infty }=0.9$ further diverges from those of lower $M_{\infty }$. The results for the adiabatic wall, corrected via figure 26, are also plotted in 27. The adiabatic wall enhances the effect of the bump on the TS wave and the destabilising effect at transonic speeds. This is particularly relevant because this condition is likely to better represent that of cruise flight.