5.1. Overview and terminology

Depicted in figure 1, the FTLE field extracted at a spanwise location of $z/\delta _{99,0}=12.57$ is compared with the corresponding magnitude of vorticity $\omega$ in a $xy$ plane through the flow field. We have chosen a grey-scale map to be consistent with earlier representations in literature. Thus, the brightest areas correspond to minima and the darkest areas to maxima in the depicted data (cf. grey-scale bar). To demonstrate that sufficient individual structures have been collected to ensure statistical convergence, the whole subdomain used in the quantitative analyses later on is depicted. Essentially, this figure illustrates the well-known picture: towards the wall, the turbulent structures become finer and finer, and laminar inclusions from the undisturbed external flow are present in the outer, intermittent region of the boundary layer. The two (arbitrarily) chosen close-ups provide a more detailed comparison between the levels of detail provided by the FTLE-based and vorticity-based visualisations.

Figure 1.

Comparison of the FTLE field (a) and vorticity magnitude (b) in a $xy$ plane.

Qualitatively, both panels (a) and (b) share great similarities with representations available in the literature in which FTLE maxima or ridge lines were associated with vortices in the flow field; see, e.g. Green et al. (Reference Green, Rowley and Haller2007) and Pan et al. (Reference Pan, Wang and Zhang2009), to mention only two. A comparison between panels (a) and (b) illustrates the already reported phenomenological similarities between the FTLE and vorticity structures, although most structures appear less sharply defined for the vorticity in panel (b) compared with the FTLE in panel (a). However, it is not readily evident from the literature how far-reaching these benefits could be. Note that the vorticity in panel (b) has the advantage of being much faster to compute and being Galilean invariant as well, as long as the reference frame is not rotating. However, shear layers and vortices cannot be distinguished, which led (among others) to the development of other methods based on the Eulerian view ($Q$, $\lambda _2$, swirling strength, etc.; see Chong, Perry & Cantwell Reference Chong, Perry and Cantwell1990; Jeong & Hussain Reference Jeong and Hussain1995; Chakraborty, Balachandar & Adrian Reference Chakraborty, Balachandar and Adrian2005).

To complement the insight given in figure 1 and allow for an assessment of the data's three dimensionality, figure 2 shows three wall-parallel slices at different wall distances $y=\mathrm {const}$. As the wall distance expressed in inner units $y^+$ changes over the streamwise direction, it is given as an averaged value $\overline {y^+}$ over the entire $x$ range for each slice. Three $\overline {y^+}$ values are chosen, one close to the wall at $\overline {y^+}=21$, one at about $50\,\%$ of the averaged boundary-layer thickness at $\overline {y^+}=315$ and one at about $112\,\%$ of the averaged boundary-layer thickness at $\overline {y^+}=706$. Again, close-ups are given for a more detailed comparison: panel (d) focuses on an ‘eruption’ event near the boundary-layer edge, panels (e) and (f) focus on regions containing local FTLE minima.

Figure 2.

The FTLE fields in different wall-parallel planes: (a) $\overline {y^+} \approx 706$, (b) $\overline {y^+} \approx 315$, (c) $\overline {y^+} \approx 21$, (d), (e), (f) zoom in corresponding to (a), (b), (c), respectively.

In contrast to other studies in the literature, it should be noted that figures 1 and 2 show a much larger data range, especially in the streamwise direction. Moreover, due to the strongly three-dimensional nature of the flow field considered (which contains mass transfer perpendicular to the cutting planes), the present visualisations substantially differ from studies of two-dimensional isotropic turbulence, e.g. Farazmand & Haller (Reference Farazmand and Haller2012). It is further mentioned that – owing to the fact that three eigenvalues occur in the three-dimensional case compared with only two in the two-dimensional case – methods and findings developed for two-dimensional flows usually cannot be simply adopted and applied.

As motivated in § 2, it is a primary goal of this paper to investigate the usually less discussed bright regions of the FTLE fields in figures 1 and 2 (in contrast to the dark regions there), to line out their physical relevance and to associate them with known processes within the turbulent boundary layer. According to the definition of the FTLE, these regions correspond to regions of very little stretching/compression. It is recalled that all three eigenvalues $\lambda _{min}$, $\lambda _{int}$ and $\lambda _{max}$ tend towards one when the maximum eigenvalue tends towards one, implying that the corresponding FTLE values computed for $\lambda _{min}$, $\lambda _{int}$ and $\lambda _{max}$ tend towards zero.

In its physical association, minimal divergence or compression over a period of time means that fluid regions involved remain together over the period of time under consideration, implying this motion to be ‘coherent’. To avoid confusion with the variety of structures already associated with the term of a ‘coherent structure’ (e.g. vortices, shear layers, LCS, etc.), it seems reasonable to distinguish these structures from others by calling them ‘LAMS’ in the following. ‘Lagrangian’ as they are obtained from a Lagrangian point of view, and ‘areas of minimal stretching’ due to the mathematical background of their detection method. Furthermore, for a clear distinction from the term of a ‘coherent structure’, the related structures will be termed to be ‘cohesive’. By definition, the LAMS lie between the ridges of the FTLE field and, thus, between neighbouring attracting and repelling LCS; whether each FTLE ridge represents an LCS on closer analysis or not; compare § 1 and Farazmand & Haller (Reference Farazmand and Haller2012). A single LAMS thus could be described as a cohesive, quiet island in the turbulent flow, which is embedded in the turbulent flow similar to a colour blob transported with the flow without being rapidly diffused.

Note that the LAMS term proposed is a kind of superordinate term with respect to the classification of LCS made by Haller (Reference Haller2015) and others, as LAMS can contain both elliptic and parabolic LCS. In comparison with the mathematical criteria that must be fulfilled to unambiguously identify elliptic or parabolic LCS, LAMS are less precisely defined, but much easier to identify visually.

5.2. Detection and interpretation of LAMS in the flow field

In the following, the general behaviour of LAMS is characterised. As already seen in figure 2, these findings depend on the wall-normal position investigated due to the three dimensionality of the flow problem considered. In the plane near the boundary-layer edge at $\overline {y^+}=706$ ($y/\bar {\delta }_{99}\approx 1.12$, $y/\delta _{99,0}\approx 8.68$) in figure 2(a), the FTLE field contains large regions of laminar external flow with interspersed ‘islands’ of increased turbulent activity rising locally upward into the laminar outer flow; see also figure 1. In the plane at $\overline {y^+} \approx 315$ ($y/\bar {\delta }_{99}\approx 0.50$, $y/\delta _{99,0}\approx 3.88$) in figure 2(b), the proportion of bright regions (LAMS) in the overall image is much smaller. Nevertheless, connected areas with small FTLE values (LAMS) still clearly stand out as isolated bright regions. The minimal FTLE values are close to zero and located in shallow ‘valleys’ surrounded by steeply rising FTLE ‘flanks’, compare the detail in figure 2(e). A comparison with the longitudinal section in figure 1(a) indicates that most of the affected regions have clear connections to the laminar external flow. Therefore, the majority of bright regions (LAMS) are expected to be simply associated with laminar regions of the outer flow sinking downward in compensation for the turbulent zones rising from below at this wall-normal position. In the lowest wall-parallel plane at $\overline {y^+} \approx 21$, the arrangement of LAMS becomes rather streaky in the flow direction, matching the structures known from visualisations of the velocity or vorticity field. Due to the complexity of the three-dimensional FTLE field, it is by far not obvious anymore whether these regions still continue to have clear connections to the laminar external flow or not. However, the underlying flow map identifies LAMS that descend a significant distance toward the wall over the time period considered. This indication is further supported by the fact that the prevailing wall-normal velocity component is negative for a large number of LAMS until near the wall; a more detailed discussion will be provided in § 5.4.

To complete the overview over the obtained three-dimensional FTLE field, the wall-normal variation of the FTLE values is presented in the form of histograms and wall-normal distributions of characteristic data values in figure 3. For this purpose, histograms of the FTLE distributions are determined and plotted for every available wall-parallel slice as probability density distributions $P$. With $400 \times 256$ grid points per slice and time step, each curve in figure 3(a) contains a total of $102\ 400 \times 100$ data points. The associated histograms are predominantly in the form of skewed Gaussian distributions evolving toward a Gaussian normal distribution with decreasing wall distance. Near the boundary-layer edge at $\overline {y^+} \approx 706$, very low FTLE values dominate the distribution due to the laminar free-stream flow that exhibits minimal stretching. In the next four curves towards the wall, the probability density function (PDF) maximum successively shifts to larger FTLE values, the distribution flattens out and less data points with small FTLE values near zero occur. No second maximum occurs with small FTLEs, which would be a welcome criterion for the detection of LAMS, because a local minimum between low and high FTLE values could be used as a threshold. In the direct vicinity of the wall at $\overline {y^+}=2$, the maximum shifts to about $0.25$ and virtually no data points with FTLE values smaller than $0.15$ are detected anymore.

Figure 3.

(a) Probability density functions (PDFs) of FTLE at different wall distances and (b) mean value, median of the FTLE and $\arg \max$ of the PDF of the FTLE versus wall distance.

Figure 3(b) quantifies the dependence of FTLE maximum, median and mean values on wall-normal position. Note that all these values differ only marginally. As expected, the mean value is slightly larger than the median. Qualitatively, this increase of FTLE towards the wall is in line with the visualisations in figure 2, where bright FTLE regions (LAMS) become smaller and rarer as one moves closer to the wall. Another plausible cause is the strongly increasing streamwise shear in the mean flow that leads to increasingly stronger stretching of the fluid near the wall.

5.3. Association with the flowmap

Although figures 1 and 2 allow for a first qualitative assessment of the FTLE fields discussed, their interpretability is somewhat limited due to the scalar nature of the FTLE values depicted as well as their two-dimensional representation. Thus, to more clearly illustrate the property of Lagrangian regions of minimal stretching, a snapshot view on particles that have been used to compute the FTLE field is depicted in figure 4. Therein, particle positions are depicted after fourteen (backward) integration steps, which initially have been uniformly distributed on the numerical grid points at $t_0$ and $\overline {y^+}\approx 315$ with constant distances ${\rm \Delta} x$ and ${\rm \Delta} z$ (${\rm \Delta} x > {\rm \Delta} z$). The top view of the integration domain is depicted in panel (d), three representative close-ups are depicted in panels (a–c), and the respective side views illustrating the wall-normal position of the particles are depicted in panels (e–f). As can be seen, the amount of integration steps is high enough to allow for a significant redistribution of the particles initially set, qualitatively revealing regions of high particle density in the flow map where particles are accumulated. Note, however, that not all of these regions are real attractors, as also regions where particles are displaced in the vertical direction such that they ‘pile up’ above each other appear as regions of high particle density in the top view. To allow for a connection with the previous illustrations of the FTLE field, particles with low FTLE and low stretching (${\rm FTLE}< 0.053$, corresponding to $20\,\%$ of the data points) are coloured red, while particles with high FTLE and high stretching (${\rm FTLE}> 0.104$, corresponding to $30\,\%$ of the data points) are coloured blue. Conceptually, at least three different types of structures should be observable in figure 4, all satisfying the criteria of less stretched regions: (I) regions of laminar external flow sinking to the wall, (II) flow regions that have erected over the integration time containing particles remaining at an almost constant distance from each other, (III) the cores of vortices with solid-body-like rotation.

Figure 4.

Particle positions in a wall-parallel plane at $\overline {y^+} \approx 315$ (d). (a–c) The $x-z$ cutaways corresponding to (d). (e–g) The corresponding $x-y$ side views of the cutaways in (d).

Visually, the flow field depicted in figure 4 contains large connected regions of red particles (see, e.g. panel a), that can easily be attributed to the first mentioned structures (I). These regions clearly exhibit a time-line structure in the top view, implying that their relative distance to each other has only slightly been changed over the considered time interval. In consequence, these regions have to be characterised by an almost homogeneous velocity, which causes these areas to move in their entirety as a ‘cohesive’ patch of fluid; see the red plateau in, e.g. panel (e). As will be elaborated in more detail below, most of these regions move upward when tracked in backward time (see, e.g. panel e), implying that they are characterised by a negative wall-normal velocity component when tracked in forward time; for the blue coloured particles, the opposite holds. Thus, for the wall-normal distance of $\overline {y^+}\approx 315$ selected, it is more likely that the large red coloured cohesive regions of minimal stretching can be associated with regions connected to the laminar free-stream flow.

However, a large number of red particles are also located in regions of high particle concentration and in close proximity to blue particles with high stretching when viewed in top view. These particles can be associated with the second type of structures (II) mentioned and are located in flow regions that have erected nearly perpendicular to the wall over the integration time while maintaining a constant distance from each other. Thus, these regions form parts of flanks in the flow field only appearing as pseudo attractors in the present backward integration. An example are the red particles in the middle left region of panel (b), which are almost horizontally stacked in the side view in panel (f). Note that the third type of structure (III) mentioned above will be discussed later in connection with a spanwise cut.

To graphically underline these statements, figure 5 presents a three-dimensional detailed view on the LAMS already selected in figure 4(b) along with its surrounding instantaneous vortices; the latter have been identified by an isosurface $\lambda _2= {\rm const.}$, see Jeong & Hussain (Reference Jeong and Hussain1995), and coloured by wall distance. In figure 5 the red time-line blob ‘floats’ in an area with little vortical activity. On the contrary, looking at areas of maximal stretching, blue particles appear much closer to the vortices.

Figure 5.

Comparison of particle positions from enlargement (a) in figure 4 with instantaneous vortices ($\lambda _2={\rm const.}$).

5.4. Association with the Eulerian representation of the flow field

Previous investigations, like Green et al. (Reference Green, Rowley and Haller2007), Pan et al. (Reference Pan, Wang and Zhang2009) or Wilson et al. (Reference Wilson, Tutkun and Cal2013) for instance, have already shown that detailed comparisons of FTLE with Eulerian representations of instantaneous flow fields have several beneficial effects: corroboration and consolidation of observations and interpretations, as well as additional insight. Figure 6 serves this purpose in addition to providing a cross-stream cut through the data to complement the previous illustrations. A superposition of FTLE (using a red pseudo-colour scale, see legend), the in-plane velocity directions in grey (using the line integral convolution (LIC) algorithm, see Cabral & Leedom Reference Cabral and Leedom1993), and green semi-transparent areas with negative $\lambda _2$, see Jeong & Hussain (Reference Jeong and Hussain1995), are shown. To separate high from low FTLE areas (LAMS), the median isoline of the FTLE field is included as a grey contour line. The dashed horizontal line is at $y^+ = 20$ for reference. Exemplary hand-drawn vectors illustrate the flow direction belonging to the LIC field, such that upward and downward motions can be identified. The latter is especially useful to confirm the presence of vortices in case that a vortex has been cut normal to its axis. Two such examples appear at $y/\delta _{99,0}\approx 7$, $z/\delta _{99,0}\approx 10$ and $y/\delta _{99,0}\approx 5$, $z/\delta _{99,0}\approx 14$. Due to the complicated three-dimensional nature of the data (cf. figure 5) vortices are only occasionally cut normal to their axes. However, in the lowest level around $y^+ \approx 20$ many streamwise oriented small vortices are intersected. These vortices belong to the well-known ‘legs’ of the lowest layer of hairpins in the log layer. A majority of these contain small FTLE values, i.e. LAMS. With respect to the expectation (III) of observing LAMS in the cores of vortices formulated above, this figure proves that the cores of several vortices indeed contain LAMS. Thus, as expected, LAMS can be viewed as a superordinate term for parabolic and elliptic LCS.

Figure 6.

Superposition of different fields in the $yz$ plane at $x/\delta _{99,0}\approx 210$: (a) FTLE and LIC, (b) FTLE and $\lambda _2$. Dashed horizontal line is at $y^+ \approx 20$; arrows indicate flow direction.

In addition, figure 6 yields insight into the dynamics of the turbulent boundary layer. Mushroom-like eruptions near the boundary-layer edge contain areas of high FTLE and vortices with local upward and downward motions. Interspersed in the spanwise direction are quieter areas with LAMS character sinking towards the wall and discernible down to $y/\delta _{99,0} \approx 5$, especially at $z/\delta _{99,0} \approx 5 \to 9$ and $z/\delta _{99,0} \approx 13$. The whole ‘middle layer’ of the boundary layer contains a mixture of high and low FTLE regions that appear entangled by the many vortices. It is clear that the view of LAMS as mere laminar inclusions coming from the free stream is not viable in such an environment. Nevertheless, it seems evident that the vortices act like a three-dimensional ‘gear’ that transports and mixes the flow. In this process new areas of high and low FTLE can be created due to shear layers in between neighbouring vortices. Thus, the effects observed in figure 6 can be understood as a mix of two effects: (i) engulfment and transport of laminar patches (LAMS) from the free stream towards the wall, and (ii) generation of new LAMS by vortex interactions. Newly generated LAMS will have a finite minimal strength when they are formed in an area of high three-dimensional stretching. Therefore, the interpretation of LAMS as quasi-laminar patches coming from the free stream must be successively extended to an interpretation also considering newly produced LAMS the smaller the distance to the wall becomes. Especially, the latter effect should be dominant also in the outer layer in a channel flow or a pipe flow where no laminar free stream exists. Despite this interpretation LAMS still catch the eye very close to the wall until $y^+ \approx 20$. Only below this value, the FTLE minima increase to discernibly higher (redder) values; see figure 3. On closer inspection these are entrained by the streamwise vortices next to the wall.

To connect the LAMS observed in the FTLE field with the unsteady velocity field of the streamwise velocity component, the $u^\prime$ field of the wall-parallel section located at $\overline {y^+}=315$ is depicted in figure 7(a) first. Regions with positive $u'$ are coloured red, with negative $u'$ blue and isolines with $u^\prime =0$ are depicted as black solid lines. To allow for a simpler assessment of the correlation between the LAMS and $u^\prime$, the isolines from figure 7(a) are overlayed in panel (b) with the FTLE field repeated; figure 7(c) depicts the $u'$ field from panel (a) filtered by ${\rm FTLE}<0.06$. By comparison between panels (a) and (b), a strikingly high correlation between areas of minimal stretching (${\rm FTLE} \rightarrow 0$) and the red coloured high-speed streaks can be observed. This impression is confirmed in panel (c), where at least the largest (or more contiguous) LAMS are characterised by positive $u^\prime$ values and, thus, travel downstream faster than the mean flow. This is consistent with the findings from figure 4(e–g), where a majority of LAMS have been identified to sink down towards the wall and, thus, carry fluid being faster than the mean flow.

Figure 7.

Snapshots in a wall-parallel plane at $\overline {y^+} \approx 315$ of (a) streamwise velocity fluctuations $u^\prime$, (b) the FTLE field corresponding to (a), (c) streamwise velocity fluctuations $u^\prime$ filtered by ${\rm FTLE}<0.06$ corresponding to (a) and (b).

When extending this study to other wall-parallel planes, the exact choice of the threshold value $\epsilon _{{\rm FTLE}}$ for extracting areas belonging to LAMS becomes crucial. While this could be chosen rather arbitrarily so far, the interpretation of trends in the wall-normal direction is only reasonable if it takes into account the dependence of the FTLE value on the wall distance already noted in figure 3, where all FTLE values, including the minimum, increase towards the wall. While the median (50 % percentiles) depicted in figure 3(b) could be a meaningful choice to separate the data points into two halves – those with small FTLE values containing LAMS and the rest – smaller percentiles of the data have been shown to be more advantageous as they more clearly extract the ‘cores’ of the LAMS in the unsteady field. Further note that the usage of percentiles makes the explicit definition of an FTLE threshold obsolete. Thus, only data points which belong to a specific percentile of all available ones (per $y^+$ location) are evaluated in the following.

Depicted in figure 8 and in direct extension to figure 7, panel (a) gives the ratio of data points of a respective percentile with positive $u^\prime$ components to those with negative components, $R_{(u'>0)/(u'<0)}$. To allow for an assessment of the robustness and, thus, the reliability of the trends observed, the percentiles were successively reduced from $50\,\%$ to $5\,\%$ of all the total data points, resulting in $y$-independent ratios between filtered and unfiltered data points; see figure 8(c).

Figure 8.

(a) Ratio of positive data points to negative data points of the fluctuations of the streamwise velocity component $R_{(u^\prime >0)/(u^\prime <0)}$ for different FTLE thresholds. (b) Ratio of positive data points to negative data points of the fluctuations of the wall-normal velocity component $R_{(v^\prime >0)/(v^\prime <0)}$ for different FTLE thresholds. (c) Ratio of filtered data points to all data points of the fluctuations of the streamwise velocity component $R_{(u^\prime _f)/(u^\prime _{all})}$ for different FTLE thresholds.

From panel (a), it is clearly confirmed that the data points with positive $u^\prime$ dominate over those with negative $u^\prime$ for all percentiles between $y^+\approx 10$ and $600$, which verifies the red shift already observed in figure 7(c). In general, the ratio $R_{(u'>0)/(u'<0)}$ shows two peaks, one at about $y^+ \approx 10$ to $20$ and one near the boundary-layer edge ($\overline {y^+}\approx Re_\tau$). For wall distances greater than about $y^+ \approx 5$ to $10$, depending on the percentile considered, the ratio $R_{(u'>0)/(u'<0)}$ is always greater than one, below, it becomes smaller than one. Although the dominance of faster regions increases sharply as the threshold value is reduced (especially for the two maxima), the qualitative results stay the same for all percentiles. As smaller threshold values emphasise more the centres of the LAMS, this implies a stronger dominance of positive $u$ fluctuations in the centres of the LAMS compared with their edges, further implying that the regions identified as LAMS are predominantly faster than the mean flow and, therefore, belong to high-speed streaks. The peak at the boundary-layer edge stems from the fact that the large sinking regions of largely undisturbed far-field flow are initially very uniform without being penetrated by the turbulence surrounding them. Below this peak, a process of ‘consuming’ the LAMS then begins, whereby they become increasingly smaller and slower. In this process, the ratio of fast to slow subregions decreases, but always remains greater than one, however. Interestingly, the ratio rises again to strikingly high values in the range $y^+ \approx 10\unicode{x2013}20$, before falling to values less than one in the direct vicinity of the wall. As will be discussed in more detail further down, the latter means that the wall-closest LAMS are dominated by low-speed streaks; compare also the blueish figure 10(c) depicted further down. Note that the minimum smaller than one in the uppermost data level (i.e. outside the boundary layer) should be taken with caution, as regions of minimal stretching are very large and extend over more than the specified percentile of the data points. If the FTLE values generally or predominantly approach zero, then restricting them to a specific percentile will lead to results that are no longer representative at the top data level.

To complement the picture drawn from panel (a) – that LAMS above roughly $\overline {y^+}\approx 10$ predominantly belong to (kind of laminar) high-speed streaks sinking down towards the wall (negative $R_{(v'>0)/(v'<0)}$) and below roughly $\overline {y^+}\approx 10$ predominantly belong to low-speed streaks lifting up from the wall (positive $R_{(v'>0)/(v'<0)}$) – the ratio of data points with positive $v^\prime$ components to those with negative components, $R_{(v'>0)/(v'<0)}$, is depicted in figure 8(b). Essentially, the trends in panel (b) are reversed to that in panel (a), but show a qualitatively identical behaviour, which confirms the above-mentioned interpretation. Note that a linear scaling is more appropriate in panel (b) compared with panel (a).

To further characterise the combined motion of LAMS in the streamwise and wall-normal direction, it is common to define ‘ejection’ and ‘sweep’ events by simultaneously considering the perturbation velocities in the flow direction and normal to the wall; see Willmarth & Lu (Reference Willmarth and Lu1972). If $u^\prime >0$ and $v^\prime <0$, then the region under consideration descends with simultaneous overspeed, which is usually called a ‘sweep event’. Since this state lies in the fourth quadrant of the $u^\prime$–$v^\prime$ plane, occasionally also the term Q4 event is used for it. Accordingly, ‘ejections’ with $u^\prime <0$ and $v^\prime >0$ are located in the second quadrant (Q2). Both criteria were applied to the present data after they were previously filtered according to the percentiles shown in figure 8. Subsequently, the ratio of Q4 to Q2 events was formed and plotted against the wall distance in figure 9(a), which results in a distribution almost identical to that of $R_{(u'>0)/(u'<0)}$ plotted in figure 8(a). However, a quantitative comparison shows that the ratio of sweeps to ejections $R_{(sw)/(ej)}$ is slightly larger than the ratio of high-speed to low-speed streaks $R_{(u'>0)/(u'<0)}$.

Figure 9.

(a) Ratio of sweep event data points to ejection event data points $R_{(sw)/(ej)}$ for different FTLE thresholds. (b) Ratio of filtered data points to all data points of the ejection and sweep events $R_{(ejsw_f)/(u^\prime _{all})}$ for different FTLE thresholds. (c) Ratio of data points of the ejection and sweep events to data points of the fluctuations of the streamwise velocity component $R_{(ejsw_f)/(u^\prime _f)}$ for different FTLE thresholds. (d) Ratio of sweep event data points to ejection event data points $R_{(sw)/(ej)}$ for different FTLE thresholds for the ZPG case.

As a control, in figure 9(b), comparable to figure 8(c), the ratio of the remaining data points to the total data points per level is plotted for different percentiles. In direct comparison with figure 8(c), one can see a reduction in the remaining data points with slight maxima at $y^+ \approx 300$ and $y^+ \approx 10$. The difference is due to data points where Q1 and Q3 events occur.

For figure 9(c), the ratio of the remaining data points from figure 9(b) and figure 8(c) was computed. Here it is noticeable that consistently $60$ to $80\,\%$ of the data points of each percentile belong to ejections and sweeps. Only at the boundary-layer edge and below $y^+ \approx 10$ a dependence on the FTLE threshold can be seen, in the way that the values increase when the threshold is reduced. The increase is monotonic. Since one gets closer to the core of the LAMS with reduction of the FTLE value, this observation means that the dominance of sweeps is higher in the centres of the LAMS than at their edges. The physical consequence of these results is that LAMS are to a large extent associated with sweep events. Only immediately at the wall, in the viscous sublayer, the ratio reverses and ejections increasingly dominate over sweeps.

Figure 9(d) presents results for the zero-pressure gradient (ZPG) case from Wenzel et al. (Reference Wenzel, Gibis, Kloker and Rist2019) for comparison with the APG case considered throughout this paper. It can be directly compared with figure 9(a). In contrast to the APG case the ZPG analysis has been performed for one arbitrarily chosen data time step only instead of hundreds. As the number of time steps had already no influence on the statistics in the APG case, we can consider this safe. Figure 9(d) proofs that the streamwise pressure gradient has no significant influence on the present observations such that we can say that the latter are not specific for the chosen pressure gradient.

It is interesting to note that Houra et al. (Reference Houra, Tsuji and Nagano2000), using quadrant splitting and trajectory analyses, reveal that in their strongest APG flow, marked transfers of high momentum fluid toward the wall (Q4) occur frequently with a short duration time. On the other hand, the contributions from other coherent motions, especially the ejection motions (Q2), significantly decrease and their duration becomes longer. Compared with ZPG, they also find that fluid motions from the outer to inner layers predominate in the APG flow. Thus, it seems that we observe similar effects in our FTLE analysis but not yet a similar influence of the APG, because their observations are taken from a case with a five times stronger Clauser parameter.

Figure 10 serves as a graphical illustration of the observations previously described. For this purpose, instantaneous data of the $u^\prime$-velocity component are filtered with the $P_{12.5}$ percentile of the FTLE distribution for three wall-parallel planes, such that only $12.5\,\%$ of all data points with the lowest FTLE values are shown in each plane. Whereas the two wall-distant planes in panels (a) and (b) exhibit a clear red shift, the wall-closest plane in panel (c) is characterised by almost only blue structures, which visually underlines the statements made above. Furthermore, the LAMS structures become smaller and smaller from $y^+ = 315$ to $y^+ = 124$. At $y^+ = 2$, nearly all of the LAMS appear as elongated blue streaks in the direction of flow, which implies that the interpretation/association of the LAMS structures detected has to significantly change for different wall-normal positions.

Figure 10.

Wall-parallel snapshots of the streamwise velocity component filtered by FTLE (12.5 % percentile) at different wall distances: (a) $\overline {y^+} \approx 315$, (b) $\overline {y^+} \approx 21$, (c) $\overline {y^+} \approx 2$.

5.5. About the peak at $y^+ \approx 20$

As observed in figures 8(a) and 9(a), the LAMS in the unsteady turbulent flow field are dominated by high-speed streaks or sweeps, in particular, at $y^+ = 10-20$, suggesting a closer look at the possible causes for this behaviour in that region. To this end, the distribution of the $\lambda _2$ minimum (Jeong & Hussain Reference Jeong and Hussain1995) is plotted versus $y$ in figure 11, which is known to correlate with vortex strength (compare, e.g. Rist, Muller & Wagner Reference Rist, Muller and Wagner1998). The smaller the $\lambda _2$ value, the stronger the associated vortex. As depicted in figure 11, the $\lambda _2$ minimum and, thus, a maximum of vortex activity is located at $y^+ \approx 20$, which correlates well with the lower maxima seen in figures 8 and 9. For the upper maximum near the boundary-layer edge, conversely, the intensity of the vortex activity does not peak and continuously tends towards zero. However, it should be noted that vortices occur intermittently at the boundary-layer edge, so that the magnitudes of the associated averaged $\overline {\lambda _2}$ values are reduced in magnitude by containing also vortex-free regions.

Figure 11.

Mean value of $\lambda _2$.

Referring to the lower minimum, the vortices located in the buffer layer at about $y^+ \approx 20$ are created by the fact that the ‘legs’ of all vortices further above are stretched here in the direction of the flow, as illustrated in the work by, e.g. Robinson (Reference Robinson1991) or Adrian (Reference Adrian2007). In this process these vortices must undergo very strong stretching in the streamwise direction due to the mean velocity gradient normal to the wall, which causes the $x$ vorticity to be continuously formed here. This interpretation can be confirmed by calculating the vortex stretching in the present data, for which the dominant $x$ component is shown in figure 12, also peaking at $y^+ \approx 20$. Below $y^+ \approx 10$, the vortex stretching has a damping effect. Thus, here the processes reverse compared with the region farther from the wall, consistent with the observations made above when considering LAMS-based ejections and sweeps.

Figure 12.

Mean of vortex stretching $\overline {{\rm D} \omega _x /{\rm D} t}$ in the streamwise direction.

From this discussion, it becomes obvious that the LAMS detected in the direct vicinity of the wall are substantially different to those observed in the outer layer. Whereas the LAMS in the outer region have been found to be associated with large regions of free-stream flow slowly sinking down towards the wall while increasingly becoming smaller and more ‘turbulent’, the near-wall LAMS have to be much more associated with vortex-related structures. This observation is already contained in the cut through the instantaneous flow field in figure 6, where a dashed horizontal line indicates the wall distance of $y^+\approx 20$. The underlying LIC (grey) and FTLE fields (red) show the mixing of areas of high and low FTLE values. From this visualisation, a majority of LAMS can be attributed to streamwise oriented vortices that are cut approximately perpendicular to their axes in the chosen representation. Areas of low and high stretching thus either get attracted towards the vortex centres or expelled from these. In between the vortices, alternating areas of upward and downward motions occur, which partly also contain LAMS. It is interesting to note that these turn into the horizontal direction at almost exactly $y^+= 20$ due to the proximity of the wall.

5.6. Corroboration of the results by other means

In summary, the LAMS detected so far have been attributed to originally laminar areas sinking down from the outer flow towards the wall in the outer region, and to vortex-bound structures in the direct vicinity of the wall. To allow for an assessment, to which degree the LAMS of the outer layer actually stay together on their way to the wall and how close they can actually reach to the wall, essentially two methods are conceivable. Firstly, fluid areas of predefined LAMS could be marked to be traced in their temporal evolution; secondly, by weakly heating the wall, a passive scalar field could be introduced; see Warhaft (Reference Warhaft2000). As the present compressible turbulent boundary layer is treated with an adiabatic, and thus, compared with the free-stream temperature, slightly hotter wall, the temperature field can serve as a direct indicator to trace the LAMS motions. To this end, the analysis from figures 8 and 9 is repeated in figure 13 for the temperature fluctuation $T^\prime$. For obvious reasons, a high qualitative agreement with the result for the $v$ component is observed, however, under the additional advantage that cold spots can be undoubtedly interpreted as structures sunk down from further outside.

Figure 13.

Ratio of positive data points to negative data points of the fluctuations of the temperature $R_{(T^\prime >0)/(T^\prime <0)}$ for different FTLE thresholds.

As can be seen from figure 14 and known from the literature, the temperature field shows individual cold regions in each snapshot (associated with LAMS) that indeed succeed in descending from the outer region down to $y^+ \approx 20$. Prominent examples can be found in panel (a) at the left edge, in the middle and at $x/\delta _{99,0} \approx 255$. To more easily allow for an assessment about the time evolution of these structures, two representative regions are marked in all time instances.

Figure 14.

Snapshots of the instantaneous temperature field at different times: (a) $t=26802$, (b) $t=26809$, (c) $t=26816$, (d) $t=26823$, (e) $t=26830$.

To further quantify the strength of the LAMS's wall-normal movement, the wall-normal velocity component $v=\bar {v}+v'$ for all data points remaining for successively smaller percentiles of the FTLE distribution were determined in each level $y={\rm const.}$ and plotted in figure 15; this method is analogous to the procedure above, e.g. in figure 8. Besides, the time-averaged mean flow $\bar {v}$ is shown for the case where no filtering is applied for comparison. Note that in figure 8(b) merely the ratio between data points with negative to positive $v^\prime$ components has been evaluated. In figure 15, outside the boundary layer, the LAMS move away from the wall at the mean velocity. Below, a rapid drop to a minimum follows as the local downward motion is stronger than the mean upward motion. Between $y^+ \approx 200$ and $100$, the downward motion is reduced before a renewed strengthening takes place in the $y^+ \approx 20$ region. As observed earlier, a reduction of the percentiles brings one closer to the centres of the involved LAMS and the associated mean downward velocity intensifies, meaning that the centres of the LAMS are sinking strongest. The upward motion below $y^+ \approx 3$ remains very small compared with the rest. When compared with $\bar {v}$, the downward motion in the individual LAMS is many times larger than the mean-flow value. To illustrate this behaviour, exemplary histograms at $y^+ \approx 315$ are shown in figure 16 as the probability density distribution $P$ for the different thresholds from figure 15. The colour coding of the individual cases is identical in both figures. Here, a very clear shift of the PDF maximum to negative $v$ is present with decreasing FTLE threshold. However, about one third of all values are still positive, pulling the mean towards zero. It is also interesting to observe that the strongest downward movement reaches close to $-0.1$, which is about ten times larger than the mean value in figure 15. This indicates that the downward motion in individual LAMS can be very violent, matching the observations in the temperature field in figure 14.

Figure 15.

Percentile-averaged wall-normal velocity component versus wall distance.

Figure 16.

Probability density functions of the wall-normal velocity component $v$ filtered for different FTLE percentiles in a wall-parallel plane at $\overline {y^+} \approx 315$.