Evidence that superstructures comprise self-similar coherent motions in high Reynolds number boundary layers

We present experimental evidence that the superstructures in turbulent boundary layers comprise smaller, geometrically self-similar coherent motions. The evidence comes from identifying and analysing instantaneous superstructures from large-scale particle image velocimetry datasets acquired at high Reynolds numbers, capable of capturing streamwise elongated motions extending up to 12 times the boundary layer thickness. Given the challenge in identifying the constituent motions of the superstructures based on streamwise velocity signatures, a new approach is adopted that analyses the wall-normal velocity ﬂuctuations within these very long motions, which reveals the constituent motions unambiguously. The conditional streamwise energy spectra of the Reynolds shear stress and the wall-normal ﬂuctuations, corresponding exclusively to the superstructure region, are found to exhibit the well-known distance-from-the-wall scaling in the intermediate-scale range. It suggests that geometrically self-similar motions are the constituent motions of these very-large-scale structures. Investigation of the spatial organization of the wall-normal momentum-carrying eddies, within the superstructures, also lends empirical support to the concatenation hypothesis for the formation of these structures. The association between the superstructures and self-similar motions is reafﬁrmed on comparing the vertical coherence of the Reynolds-shear-stress-carrying motions, by computing conditionally averaged two-point correlations, which are found to match with the mean correlations. The mean vertical coherence of these motions, investigated for the log region across three decades of Reynolds numbers, exhibits a unique distance-from-the-wall scaling invariant with Reynolds number. The ﬁndings support modelling of these dynamically signiﬁcant motions via data-driven coherent structure-based models.


Introduction and motivation
Over the past two decades, the study of high Reynolds number (Re τ O(10 4 )) wallbounded flows has become synonymous with very-large-scale motions (VLSMs), also known as 'superstructures', which play a predominant role in the dynamics and spatial organization of wall turbulence.Here, Re τ = δU τ /ν, where δ is the boundary layer thickness, ν is the kinematic viscosity and U τ is the skin-friction velocity, with the latter two used to normalize the statistics in viscous units (indicated by superscript '+').The superstructures can extend beyond 20δ in the streamwise direction (Kim & Adrian 1999;Hutchins & Marusic 2007) and also exhibit 'meandering' when viewed on a wall-parallel plane (de Silva et al. 2015), particularly in the logarithmic region of the flow.Such a large spatial footprint permits these motions to carry significant proportions of the total turbulent kinetic energy and the Reynolds shear stresses of the flow (Liu et al. 2001;Guala et al. 2006;Balakumar & Adrian 2007).Given that the shear stress is responsible for the wall-normal momentum transfer, this suggests that the VLSMs/superstructures also contribute significantly to the high Re τ turbulent skin-friction drag (Deck et al. 2014).Hence, an improved understanding of the origin of these VLSMs/superstructures, towards which this study is directed, stands to advance our knowledge in both a fundamental and an applied perspective.Hutchins & Marusic (2007) used the terminology 'superstructures' when referring to the spectrogram of the streamwise velocity fluctuations (u) from a high Re τ boundary layer, as shown in figure 1.The spectrogram presents the premultiplied u-energy spectra as a function of the viscous-scaled streamwise wavelengths (λ + x = λ x U τ /ν) and wall-normal distance (z + = zU τ /ν), with λ x = 2π/k x , where k x is the streamwise wavenumber.The high Re τ u-spectrogram is seen to have two prominent peaks.One is located in the inner-region synonymous with the well-documented near-wall cycle (Kline et al. 1967), consisting of high and low-speed viscous-scaled streaks (λ + x ≈ 1000), which are responsible for intense local production of turbulent kinetic energy.The second peak is in the outer region of the flow (typically in the logarithmic/inertial region), and corresponds to the superstructures, which have a spectral signature at very long wavelengths (λ x ∼ 6δ) and also extend down to the wall (Hutchins & Marusic 2007).It is worth noting here that this second peak is only visible for Re τ 2000, owing to the insufficient separation of scales and weaker energy of the superstructures at lower Re τ (Hutchins & Marusic 2007).Between the innerand outer-peaks, a nominal plateau is seen in the spectrogram which corresponds to the distance-from-the-wall (z)-scaled eddies coexisting in the log-region; these eddies make up the increased range of scales with increasing Re τ .In the literature, these intermediate scaled eddies have been described by various structures or motions, including the large-scale motions (LSMs; Kim & Adrian 1999, Adrian et al. 2000), uniform momentum zones (UMZs; Meinhart & Adrian 1995, de Silva et al. 2016), attached eddies (Baars et al. 2017;Marusic & Monty 2019;Hu et al. 2020;Deshpande et al. 2021a) and so forth.In the remainder of this section, for simplicity, we will refer to these motions as LSMs.It should also be noted that the terminology 'VLSMs' and 'superstructures' have been conventionally associated with the very-large-scale motions in internal (Kim & Adrian 1999) and external wall-bounded flows (Hutchins & Marusic 2007), respectively.Considering this study focuses solely on zero-pressure gradient turbulent boundary layers, we henceforth refer to either of these structures simply as superstructures.
To date, several studies have investigated the probable mechanisms responsible for the formation of superstructures, with two theories hypothesized most often: (i) the formation of superstructures via concatenation of the LSMs (Kim & Adrian 1999;Adrian et al. 2000;Lee & Sung 2011;Dennis & Nickels 2011), or (ii) the emergence of superstructures due to a linear instability mechanism (Del Alamo & Jimenez 2006;McKeon & Sharma 2010;Hwang & Cossu 2010).The present study does not focus on comparing and contrasting the likelihood of one mechanism over the other.Rather, it builds upon recent compelling evidence in support of the concatenation mechanism (Wu et al. 2012;Baltzer et al. 2013;Lee et al. 2014Lee et al. , 2019)), to investigate the characteristics of the constituent motions forming the superstructures.The formation of superstructures via streamwise concatenation of the relatively smaller motions has been confirmed by several studies conducted across all canonical wall-bounded flows (turbulent boundary layers, channels, pipes), through: (i) investigation of the time evolution of instantaneous flow fields (Lee & Sung 2011;Dennis & Nickels 2011;Wu et al. 2012;Lee et al. 2019), (ii) statistical analysis of the superstructure formation frequency/population density (Lee et al. 2014) and (iii) spatial correlations of the low-pass filtered velocity fields (Baltzer et al. 2013;Lee et al. 2019).In comparison, few studies have presented similar statistical arguments in favour of the linear instability mechanism.For instance, Bailey et al. (2008) supported the linear instability argument by noting different spanwise widths of the superstructures and LSMs in the inertial region of a turbulent pipe flow.Their estimates, however, were limited to two-point velocity correlations reconstructed in a particular wallparallel plane, which cannot be uniquely associated with the LSMs responsible for the superstructure formation (Deshpande et al. 2020).Considering that superstructures extend down from the log-region to the wall, Deshpande et al. (2021b) reconstructed two-point velocity correlations across two wall-parallel planes located in the near-wall and the log-region.These statistics, which are purely representative of the large 'wallcoherent' motions, revealed similar spanwise extents of the coexisting superstructures and LSMs for all canonical wall flows, thereby favouring the concatenation argument.
Despite substantial support for the concatenation argument, several unanswered questions are still associated with this mechanism.For instance, there is no universal agreement on what facilitates the streamwise concatenation of LSMs to form superstructures.While few studies have associated this with the spanwise alternate positioning of low and high momentum LSMs (Lee et al. 2014), others have conjectured the role played by secondary roll cells (Baltzer et al. 2013;Lee et al. 2019) in favourably organizing the relatively smaller motions.Progress in this regard has been hindered by the lack of understanding of the constituent motions forming the superstructures; for instance, are superstructures purely composed of the inertial δ-scaled motions corresponding to the extreme right end of region II in figure 1? Or do they also comprise of the geometrically self-similar, i.e. z-scaled hierarchy of eddies encompassing the entirety of region II?The present study aims to answer these questions by analyzing the characteristics of the constituent motions.
In the past, clarifying such information on the constituent motions has not been possible due to the low to moderate Re τ ( 2000) of the experiments/simulations analyzing the concatenation argument, which severely constricts the extent of region (II) in figure 1.This prevents an unambiguous delineation between the δ-scaled and z-scaled inertial motions coexisting in region II.Further, the statistical signature of the superstructures is also very weak at these Re τ (Hutchins & Marusic 2007), making it challenging to identify and isolate them from the other motions in the flow.However, increased access to high Re τ data over the past decade has substantially increased our knowledge of these inertial eddies coexisting in the log and outer regions (Marusic et al. 2015;Baidya et al. 2017;Deshpande et al. 2021a).This has also led to growing acceptance of the existence of the geometrically selfsimilar attached eddy hierarchy in the inertial region (de Silva et al. 2016;Baars et al. 2017;Hwang & Sung 2018;Hu et al. 2020;Deshpande et al. 2020Deshpande et al. , 2021a)), which can be modelled conceptually (Marusic & Monty 2019).These advancements make it compelling to investigate whether these self-similar inertial motions are associated with the formation of superstructures, a conjecture that has previously shown promising results when tested for low Re τ channel flows (Lozano-Durán et al. 2012), and when implemented in coherent structure-based models (Deshpande et al. 2021b).If this conjecture is proven true, then the preferred streamwise alignment of this energy-containing hierarchy of motions (to form superstructures) would have implications on Townsend's attached eddy hypothesis, which otherwise assumes a random distribution of attached eddies in the flow field (Townsend 1976;Marusic & Monty 2019).The investigation can also help answer the long-standing contradiction (Guala et al. 2006;Balakumar & Adrian 2007;Wu et al. 2012) between: (i) the attached eddy hypothesis, which classifies turbulent superstructures to be 'inactive' (Deshpande et al. 2021a), and (ii) instantaneous flow field observations, per which Focus on Fluids articles must not exceed this page length these streamwise elongated motions carry significant Reynolds shear stresses (and hence behave as 'active' motions).
To this end, the present study investigates the geometric scalings exhibited by the constituent motions of the superstructures.Experimental data is employed from a moderate to high Re τ turbulent boundary layer (2500 Re τ 7500), which is an order of magnitude higher than the simulation studies reported previously, to ensure coexistence of a broad range of inertial scales (region II).The dataset comprises of sufficiently resolved large-scale velocity fluctuations acquired in a physically thick boundary layer via unique, large field-of-view (LFOV) particle image velocimetry (PIV), capturing instantaneous flow fields with an extent of 12δ in the streamwise direction (x).In contrast to most studies to date, which have investigated the superstructures by analyzing the large-scale u-fluctuations, here we adopt a unique strategy to investigate the wall-normal (w) velocity fluctuations within the superstructure region.This is because deciphering smaller constituent u-motions from within a larger u-motion can be inconclusive, as can be noted from a sample DNS flow field shown in figures 2(a,b).On the other hand, the wfluctuations can bring out the individual constituent motions more distinctly, which is evident from figure 2(c) and will be analyzed here by computing conditional statistics.It can be noted from figures 2(a-c) that the individual w-eddies within the region associated with a long u-motion are much smaller in streamwise extent (than u), and exhibit sort of a clustered/packed organization plausibly leading to the appearance of a u-superstructure.This scenario is recreated in figure 2(d,e), using an idealized distribution of prograde vortices, which suggests the possibility of strong u-as well as w-correlations extending across large streamwise separations.Such a flow organization, which adds further credibility to the streamwise concatenation hypothesis, will be investigated here via conditional statistics from high Re τ data.3a).Terminology has been defined in §2.∆x + and ∆z + indicate viscous-scaled spatial resolution along x and z directions, respectively.
It is important to note that in the present study, any reference to concatenation henceforth refers to the spatial organization of constituent motions over extended streamwise distances, such as in figures 2(d,e).Given the experimental limitations, the study cannot directly comment on the dynamics/mechanism behind how this spatial organization comes into existence.Also, the terminology 'attached eddies' is used here to refer to any eddies/motions scaling with their distance from the wall, and hence is not limited to the eddies physically extending to the wall.

Experimental datasets and methodology
2.1.Description of the experimental datasets Five multipoint datasets are used from previously published high Re τ experiments (table 1).Four of these are acquired via two-dimensional (2-D) two-component PIV in the Melbourne wind tunnel (HRNBLWT; Marusic et al. 2015) and span the Re τ range ∼ 2500-14500.The test section of this wind tunnel has a cross-section of 0.92 m × 1.89 m, and has a large streamwise development length of ∼27 m, with maximum possible free-stream speeds (U ∞ ) of up to 45 ms −1 .Such a large-scale facility permits the generation of a sufficiently high Re τ canonical boundary layer flow facilitated by substantial increment in its boundary layer thickness, along its long streamwise fetch.This capability is leveraged in the four PIV datasets employed in the present study, which will be described next.
Three of the PIV datasets comprise snapshots of very large streamwise wall-normal flow fields of a turbulent boundary layer (x × z ∼ 12δ × 1.2δ), and are thus henceforth referred to as the large field-of-view (LFOV) PIV datasets (de Silva et al. 2015(de Silva et al. , 2020)).To the best of the authors' knowledge, this is the only published lab-based dataset giving access to sufficient LFOV instantaneous flow fields at Re τ 5000 (to achieve statistical convergence), thereby making the analysis presented in this paper unique as well as ideally-suited for investigating turbulent superstructures.The LFOV is made possible by stitching the imaged flow fields from eight high-resolution 14 bit PCO 4000 PIV cameras, each with a sensor resolution of 4008 × 2672 pixels.Figure 3(a) shows a schematic of the experimental setup for the LFOV PIV, where the region shaded in orange indicates the individual FOVs combined from the eight cameras.These measurements were conducted at the upstream end of the test section, with the LFOV starting at x ≈ 4.5 m from the start of the test section.The experiments were conducted at three free-stream speeds (U ∞ ≈ 10, 20 and 30 ms −1 ), which led to a corresponding variation in Re τ of 2500, 5000 and 7500, respectively.Here, U τ and δ used to estimate the flow Re τ , were computed at the middle of the LFOV, using the method outlined in Chauhan et al. (2009).The boundary layer thickness is nominally δ ≈ 0.11 m for all three Re τ cases.
Considering the focus of the experiment was on a LFOV, a homogeneous seeding density was ensured across the entire test section of the tunnel for these measurements, and the particles were illuminated by a Big Sky Nd-YAG double pulse laser (∼1 ṁm thickness), delivering 120 mJ/pulse.The last optical mirror to direct this laser sheet was tactically placed within the test section (figure 3a), for ensuring adequate laser illumination levels across the LFOV.This optic arrangement, however, was sufficiently downstream of the PIV flow field and introduced no adverse effects (such as blockage, etc) on the measurement (de Silva et al. 2015).Figures 3(b,c) gives an example of the viscous-scaled u-and w-fluctuations estimated from the LFOV PIV experiment at Re τ ≈ 2500, which successfully captures a turbulent superstructure (of length L x ), as highlighted by a dashed green box in the u-field.Analysis on such a dataset not only avoids uncertainties due to Taylor's hypothesis approximation (Dennis & Nickels 2008;del Álamo & Jiménez 2009;Wu et al. 2012), but also permits identification of these superstructures directly from an instantaneous flow field of a high Re τ boundary layer (where superstructures are statistically significant).The latter represents another unique feature of the present study, and overcomes the limitations experienced by past experimental studies (Liu et al. 2001;Guala et al. 2006;Balakumar & Adrian 2007), which were restricted to isolating superstructure characteristics based on Fourier-filtering, or Proper orthogonal decomposition (POD)-based decomposition of ensemble/time-averaged statistics.The accuracy of these LFOV PIV datasets have been firmly established in appendix 1 ( §5), which compares the premultiplied 1-D spectra obtained from the present data, with those acquired via multiwire anemometry published previously (Morrill-Winter et al. 2015;Baidya et al. 2017).Readers can also refer to the same appendix section for details associated with the computation of the velocity spectra from PIV flow fields, which is relevant to the analysis presented ahead in the paper.
The fourth and final PIV dataset comprises of relatively smaller flow fields in the x-z plane (in terms of δ-scaling), and is hence referred to as simply the PIV dataset.This was acquired at U ∞ ≈ 20 ms −1 , close to the downstream end of the test section (x ≈ 21 m from the trip), where δ ≈ 0.3 m, yielding a high Re τ ≈ 14500.The full velocity field captured in this experiment was also made possible by using the same eight PCO 4000 cameras, arranged in two vertical rows of four cameras each, to capture the significantly thicker boundary layer (refer to figures 1-2 of de Silva et al. 2014).This limits the streamwise extent of the flow field to x ∼ 2δ in this case, and is hence not used for identifying the turbulent superstructures in instantaneous fields, but rather used to compute the two-point correlations of u-and w-fluctuations along the z-direction (limited to the inner-region).It is owing to this reason that only a part of the full flow field (x × z ∼ 2δ × 0.4δ), from this dataset, has been considered in the present study.The image pairs from all four PIV datasets were processed via an in-house PIV package developed by the Melbourne group (de Silva et al. 2014), with the final window sizes (∆x + ,∆z + ) used for processing given in table 1. Interested readers may refer to the cited references for further details about the experimental setup and methodology adopted for acquiring these datasets.
The fifth dataset, which is at the highest Re τ ∼ O(10 6 ), was acquired at the Surface Layer Turbulence and Environmental Science Test (SLTEST) facility in the salt flats of western Utah.The data is acquired from a spanwise and wall-normal array of 18 sonic anemometers (Campbell Scientific CSAT3) arranged in an 'L'shaped configuration (refer to figure 1 of Hutchins et al. 2012).While the full dataset comprises of continuous measurements of all three velocity components as well as the temperature at the SLTEST site over a duration of nine days, here we limit our attention solely to one hour of data associated with near-neutral (i.e.near canonical) atmospheric boundary layer conditions (Hutchins et al. 2012).These conditions were confirmed based on estimation of the Monin-Obukhov similarity parameter, determined on averaging across the 10 sonic anemometers placed along the spanwise array, at a fixed distance from the wall (z ≈ 2.14 m).For the present analysis, we are solely interested in the u-and w-fluctuations measured synchronously by the nine sonic anemometers on the wall-normal array, which were placed between 1.42 m ≤ z ≤ 25.69 m with logarithmic spacing.Mean streamwise velocity measurements reported by Hutchins et al. (2012) confirm that all these z-locations fall within the log-region of the atmospheric boundary layer.This data is also used here to compute the two-point correlations of u-and w-fluctuations along the z-direction, for comparison with those obtained from the PIV datasets acquired in the laboratory.

Methodology employed to identify and extract turbulent superstructures
In the present study, we are interested in computing conditional statistics of the velocity fluctuations associated with the superstructures, identified from the individual flow fields in the LFOV PIV dataset.To identify these structures, we need to first define what we mean by a superstructure, for which we draw inspiration from past studies that have investigated these motions based on 3-D instantaneous flow fields (Hutchins & Marusic 2007;Dennis & Nickels 2011;Lee & Sung 2011).Those studies, as noted by Smits et al. (2011), refer to superstructures as "very long, meandering, features consisting of narrow regions of low-streamwise-momentum fluid flanked by regions of higher-momentum fluid", that "have also been observed in the logarithmic and wake regions of wall flows."Here, for the purpose of analyzing 2-D velocity fields, we define superstructures as very large-scale motions that persist spatially with coherent regions of streamwise velocity, and account for a significant fraction of the streamwise turbulent kinetic energy.Identifying these structures from the PIV field, hence, requires establishing logical thresholds to the geometric and kinematic properties of the fluctuating u-field (Hwang & Sung 2018;de Silva et al. 2020).For this, we consider previous findings and adopt the following thresholds: (Liu et al. 2001), where u 2 (z) is the root-mean-square of the u-fluctuations at z.
(iii) wall-normal extent should at least span across 2.6 √ Re τ z + 0.5Re τ (Guala et al. 2006;Hutchins & Marusic 2007;Balakumar & Adrian 2007;Deshpande et al. 2021b).In the process of identifying a superstructure, the threshold associated with the streamwise turbulent kinetic energy (i.e.(i)) is considered first before applying thresholds associated with the geometric extent ((ii) and (iii)).With regards to criteria (ii), we acknowledge that past studies investigating 3-D instantaneous flow fields (Hutchins & Marusic 2007;Lee & Sung 2011;Dennis & Nickels 2011) have found superstructures to be as long as 10-20δ.However, statistical analysis based on 1-D one-/two-point correlations (Guala et al. 2006;Hutchins & Marusic 2007;Balakumar & Adrian 2007;Deshpande et al. 2021b) suggests these structures have relatively modest lengths (on average), between 3-6δ.Considering that the present analysis is also limited to 2-D flow fields, we adapt the estimates from past statistical analyses and consider u-structures with streamwise extent, L x > 3δ as superstructures.Figure 3(d) gives an example of a -u superstructure identified and extracted by the algorithm (u| SS ), based on the aforementioned thresholds from the full flow field depicted in figure 3(b) (highlighted by the dashed green box).Streamwise extent/length of the identified structures (L x ) is judged based on the length of a rectangular bounding box (along x) that fully encompasses the identified structure.Our superstructure identification algorithm extracts the rectangular 2-D flow field within this box to conduct further conditional analysis associated with the superstructures.Although the choice of a rectangular box inevitably also brings in some part of the flow not associated with a superstructure, it only forms a minor part (∼20%) of the bounding box, suggesting conditional statistics can be predominantly associated with the superstructures.Interested readers are referred to the supplementary document provided along with this manuscript, which provides a step-by-step description of the superstructure identification and extraction procedure from a 2-D PIV flow field.
Besides identifying a superstructure, which is indicated by a dashed green box in figures 3(b,c), the algorithm also identifies a region of same length×height as the green box but not associated with a superstructure (u| noSS ).The u| noSS flow field region is allocated by the algorithm in the same wall-normal range as u| SS , but in a different streamwise location within the PIV image that does not satisfy criteria (i-ii) defined above, thereby ensuring it doesn't overlap with u| SS .This practice of extracting u| noSS , from the same PIV fields used to extract u| SS and of the same size as that of u| SS , is conducted across all three LFOV datasets to form a set of u| noSS and u| SS of equal ensembles.Conditional statistics are computed and compared from both u| SS and u| noSS , with the latter considered to confirm that the trends depicted by the former are not an artefact of aliasing or insufficient ensembling/noise.The superstructure extraction algorithm described above, identified superstructures of both +u and -u signatures, of varying lengths, from the three LFOV PIV datasets.A summary of their streamwise extents is presented in the form of a probability distribution function (pdf ) plot in figure 4. The plot is obtained by sorting the identified u-motions into bins of width 0.5δ (between 3.0:0.5:6.0), based on their respective lengths (L x ).The population associated with each bin is then normalized by the total number of −u and +u superstructures identified by the algorithm (for L x > 3δ), which is then plotted in the figure.It can be noted from the plots that the pdf s do not change significantly with Re τ for structures of lengths, L x < 5δ.It is only when L x is increased significantly (> 5δ) that notable differences appear for different Re τ .For example, no −u or +u-structures are identified in certain PIV datasets while in others, the probability is low.Further, the logarithmic scaling of the vertical axis of the plots reveals that the population density decreases near exponentially as the criteria (ii) to identify a superstructure (i.e.minimum length, L x ) is increased.The effect of increasing the minimum streamwise extent of a ustructure to qualify as a superstructure, on the conditionally averaged statistics, has been documented in figure 15 in Appendix 2 ( §5).Given that an increase in L x does not change the scaling behaviour, but significantly reduces the convergence of the x , computed from the LFOV PIV dataset at various Reτ .Dash-dotted golden and magenta lines represent the relationships λx ≈ 2z and λx ≈ 15z, respectively following Baidya et al. (2017).(c) Schematic of representative w and uw-carrying eddies centred at various distances from the wall (zr) in the log region, with light to dark shading used to suggest an increase in zr.Rww(z/zr) and Ruw(z/zr) respectively represent the vertical coherence of the w-and uw-carrying eddies centred at zr. conditioned statistics (due to fewer ensembles), reinforces the choice of L x 3δ in criteria (ii) discussed above.

Mean statistics
Before investigating the conditionally averaged statistics associated with the superstructures, it is worth revisiting the scaling behaviour of the mean statistics, against which the former would be compared.Here, the mean statistics have been obtained by averaging across all 3000 flow fields, and considering the entire 12δ long flow fields in case of the LFOV PIV datasets.In the present study, since we are primarily interested in the w-velocity behaviour associated with superstructures, we investigate the mean spatial coherence of the w-carrying eddies in the log-region of a high Re τ boundary layer.We look at the spatial coherence in both the streamwise (figure 5) as well as wall-normal direction (figure 6), for both the w-fluctuations and the Reynolds shear stress (uw).Previous investigations on the vertical coherence have been rare compared to the streamwise coherence, particularly for the log-region of a high Re τ boundary layer, owing to the lack of large-scale PIV experiments of the kind utilized here.This makes the present investigation (figure 6) unique by itself.
Figures 5(a,b) depict the iso-contours of the premultiplied spectrogram of the w-velocity and the Reynolds shear stress respectively, computed from the three LFOV PIV datasets.These are plotted as a function of λ + x and z + .The iso-contours for the w-velocity spectrograms can be seen centred around the linear (z-)scaling indicated by λ x = 2z for all Re τ , which is consistent with previous observations in the literature (Baidya et al. 2017).Similarly, the iso-contours for the Reynolds shear stress spectrograms also follow a linear scaling (λ x = 15z) for all Re τ , again consistent with the literature (Baidya et al. 2017).This analysis not only validates the spectra estimated from the LFOV PIV, but also assists with the construction of a simplified 2-D conceptual picture of the w-and uw-carrying eddies in the log-region of a high Re τ boundary layer (figure 5c).Here, based on the z-scaling exhibited by the data, the lengths (λ x ) of the w-and uw-carrying eddies have been defined as 2z r and 15z r respectively, where z r represents the distance of the eddy centre from the wall.This scaling confirms the association of these w-and uw-carrying eddies with Townsend's attached eddy hierarchy, according to which attached eddies scale with z r (Townsend 1976;Baidya et al. 2017;Deshpande et al. 2021a).
While both these linear scalings, which represent the streamwise coherence of the w-and uw-carrying eddies, are well accepted in the literature, not much is known about the vertical/wall-normal coherence of the same eddying motions at high Re τ .Several previous studies (Comte-Bellot 1963;Tritton 1967;Sabot et al. 1973;Hunt et al. 1987;Liu et al. 2001;Sillero et al. 2014) have investigated their vertical coherence in low Re τ canonical wall flows via traditional two-point correlations, providing interesting insights on their scaling.Here, we are inspired by one such interesting result reported in the seminal work of Hunt et al. (1988), based on high Re τ unstably stratified atmospheric boundary layer data, who found the two-point correlation coefficients given by: to be a function of (z/z r ).Here, z r acts as the reference wall-normal location fixed in the log-region, thereby making R ww (or R uw ) representative of the vertical coherence of the eddy centred at z r .It should be noted here that these correlation functions are different from the conventionally used two-point correlations (which consider normalization by the root-mean-square of velocity at both z and z r ), and hence their values aren't restricted between -1 and 1. Equations (3.1), however, are ideally suited for the present study, which tests the self-similarity (i.e.z-scaling) of the vertical coherence of the momentum carrying eddies.We compute these correlations for the four high Re τ boundary layer datasets considered and plot them in figure 6, for various z r restricted to the log-region.It can be clearly observed that R uw curves for varying z r and Re τ collapse over one another (represented by a line in teal colour based on least-squares fit), suggesting Re τ -invariance via z-scaling of the vertical coherence of uw-carrying motions.On the other hand, the collapse in the R ww curves is not as good for the relatively low Re τ cases (< 7500), but certainly gets better for the very high Re τ atmospheric data (figure 6d).This case has a significantly thicker log-region than the boundary layers generated in the lab, suggesting the influence of the wall behind the relatively poor collapse of R ww at low Re τ .Accordingly, the z-scaling of the R ww curves has been represented by the golden lines in figures 6(a-d) (obtained by a least-squares fit), which are consistent with R ww curves in figure 6(d), as well as R ww estimated farthest from the wall (z + r ≈ 0.2Re τ ) in figures 6(a-c).The analytical expressions associated with these golden and teal lines are: (3.2) The fact that R ww and R uw are solely a function of z/z r represents geometric self-similarity in the vertical coherence of the w-and uw-carrying inertial eddies, reaffirming their association with Townsend's attached eddies.The analytical forms in (3.2) can thus be used in data-driven coherent structure-based models (Deshpande et al. 2021b) to simulate high Re τ boundary layers (such as atmospheric surface layers).It is worth noting that the collapse in the R uw and R ww curves, observed in figure 6, does not exist for w-and uw-carrying eddies centred far outside the log-region of the boundary layer (i.e.z r > 0.2δ; not shown here), which may be due to the growing influence of the turbulent/non-turbulent interface in the outer-region (de Silva et al. 2014).Investigations for z r below the log-region, however, were not possible owing to insufficient data points captured by the LFOV PIV.

Conditionally averaged statistics associated with superstructures
With the scaling behaviour of the mean statistics established in §3, we progress next towards analyzing the conditionally averaged statistics (spectra and correlations) associated with superstructures.Figure 7 plots the conditionally averaged, premultiplied u-spectra computed from the extracted flow fields associated with superstructures (figure 3d), from the three LFOV PIV datasets.The spectra are plotted for z + ≈ 2.6 √ Re τ and 0.5Re τ , and estimated individually from the extracted flow fields associated with low-momentum (k x φ + uu | −u ss ; in blue) and high-momentum superstructures (k x φ + uu | +u ss ; in red).Also plotted is the conditionally averaged spectra considering both -u ss and +u ss (k x φ + uu | −u ss ,+u ss ; in green), which is compared against the mean u-spectra shown in figures 14(a-c).A noteworthy observation from the conditionally averaged spectra (k x φ + uu | −u ss ,+u ss ) is the enhanced largescale energy (λ + x 10 4 ) seen for all three Re τ cases.These enhanced energy levels are due to the significant streamwise turbulent kinetic energy associated with the superstructures, which is captured in the extracted flow fields and averaged across fewer ensembles, than those used for obtaining the mean spectra.To confirm that  (c,f) 7500.Dashed black lines correspond to the mean spectra obtained by ensembling across 3000 PIV images of the full flow field.While, the solid blue and red lines represent conditional spectra computed from the u-flow field extracted based on identification of a -uss and +uss, respectively.The spectra in green is computed by ensembling across both -uss and +uss.The mean spectra estimated from the full flow field (in black lines) is ensembled across all 3000 fields.While, the conditional spectra corresponds to extracted flow fields (∼ 300) of the same length×height associated (in green) and not associated (in brown) with the superstructures.these trends are not an artefact of aliasing or ensembling, figure 8 compares the conditionally averaged spectra associated with superstructures (green boxes in figure 3c) with that not associated with the superstructures (brown boxes in figure 3c).Given that both the conditional spectra are estimated from the same number of extracted flow fields, of the same length×height, the enhanced energy in the largest scales for k x φ + uu | −u ss ,+u ss (compared to k x φ + uu | noSS ) can be unambiguously associated with the turbulent superstructures.These trends give us confidence regarding the efficacy of the superstructure extraction algorithm.Also, they indicate that the scalings observed from the conditionally averaged u-, w-statistics can be associated with the constituent motions of superstructures.This is one of the advantages of analyzing very-large-scale motions based on extraction of instantaneous flow fields (present study), as compared with the much simpler approach of Fourier filtering (past studies).
Another interesting observation from the conditional spectra for low-and highmomentum motions, k x φ + uu | −u ss and k x φ + uu | +u ss , is their starkly different behaviour in the lower portion of the log-region (figures 7a-c) and outside of it (figures 7d-e).
While ensemble-averaged co-spectra (k x φ + uw ) at various z + within the log-region.Similar to that noted for the w-spectra, the z-scaling observed in the ensemble-averaged cospectra (λ x = 15z; figure 5b) is also noted for the k x φ + uw | −u ss ,+u ss , confirming our claim that the self-similar motions coexist within the superstructure region.This comparison between k x φ + uw | −u ss ,+u ss and k x φ + uw also showcases the significance of analyzing the very-large-scale motions by extracting instantaneous flow fields, than using pure Fourier filtering.While the latter is simpler to execute, it doesn't present the 'full physical picture' associated with the very-large-scale motions.It is only after extraction of the instantaneous flow fields at high Re τ that the present study can confirm the z-scaling characteristics associated with the constituent motions of the superstructures (k x φ + uw | −u ss ,+u ss ).In figure 10, again, high energy levels can be noted in k x φ + uw | −u ss ,+u ss at very large λ x , the magnitude of which is much greater than the energy levels for k x φ + uw | noSS and k x φ + uw at the same λ x .Further analysis is presented in §4.1 to reaffirm that these peaks do not represent very-large-scale w-motions existing in the physical flow field.
While the conditional 1-D spectra brings out the geometric characteristics of the constituent motions along the streamwise direction, the former can be understood for the wall-normal direction by computing the two-point correlations (R ww ; (3.1)) for the extracted flow fields.Figure 11 plots R ww | −u ss ,+u ss , i.e. the two-point correlations computed from the w-fluctuations associated with both -u ss and +u ss , for z r limited to the log-region.These are estimated for all three LFOV PIV datasets and compared with the least-squares fit (given by (3.2)) estimated from the mean statistics (plotted with a golden line).Consistent with our observations based on the mean statistics in figure 6, the collapse in the R ww | −u ss ,+u ss curves is not very good at low Re τ but improves significantly at Re τ ≈ 7500.Interestingly, however, R ww | −u ss ,+u ss curves estimated at all Re τ are close to the empirically obtained least-squares fit.Hence, investigation of the vertical coherence of the w-carrying eddies (associated with superstructures) also indicates that the geometrically self-similar eddies coexist within superstructure region, in line with interpretations based on figures 9 and 10.While the present study lacks the analysis to investigate the spanwise coherence of the constituent motions, consideration of the present findings in light of the recent knowledge on the log region & Sung 2018; Deshpande et al. 2020Deshpande et al. , 2021a,b) ,b) suggests that they likely exhibit self-similar characteristics along span as well.Notably, Deshpande et al. (2020Deshpande et al. ( , 2021b) ) found that the spanwise extent of the wall-coherent, intermediate-scaled motions (λ x 4δ) varies self-similarly with respect to their streamwise extent, which directly corresponds to the scale-range associated with the constituent motions of the superstructures.

Physical interpretations and discussions on the conditionally averaged statistics
Here, we discuss the physical interpretation of the conditionally-averaged spectra presented in figures 7-10, and how it advances our understanding of the constituent motions forming the turbulent superstructures.Given the geometry of individual w-eddies does not physically conform with the very-large-scale peaks noted in k x φ + ww | −u ss ,+u ss and k x φ + uw | −u ss ,+u ss (discussed previously based on figures 2c,e), these peaks are likely an artefact of the preservation of the covariance tensor, which is a property of the Fourier transform.However, the non-zero correlation between u and w-fluctuations, at large λ x , has often been misinterpreted to be representative of instantaneous w-features physically as long as the superstructures (as also highlighted by Lozano-Durán et al. 2012 andSillero et al. 2014), especially when one analyzes it from the perspective of the structure parameter (∼O(1) for large λ x ).Here, we prove from our analysis that this interpretation is incorrect.If one observes w| ss plotted in figure 3(e), which is conditioned with respect to a −u superstructure, it is clear there are no long and energetic w-features extending beyond 3δ.To the best of the authors' knowledge, energetic w-features of such long streamwise extents have never been noted in instantaneous flow fields, and their absence can also be confirmed from the negligible energy in the 1-D w-spectra plotted in figure 5(a), or in the literature (Baidya et al. 2017).Absence of very-long ( 3δ) w-features also means there are no very-large-scale Reynolds shear stress-carrying motions in the instantaneous flow (Lozano-Durán et al. 2012;Sillero et al. 2014).Such misinterpretations are the source for the long-standing contradictions between the attached eddy hypothesis and past studies (Guala et al. 2006;Balakumar & Adrian 2007;Wu et al. 2012) investigating the Reynolds shear stress co-spectra (refer §1), which we attempt to clarify here.
To reaffirm that the very-large-scale peaks in k x φ + ww | −u ss ,+u ss do not correspond with very long and energetic w-features in the instantaneous flow, figure 12 analyzes the streamwise extent of w-eddies (L w x ) in the extracted w| SS and w| noSS fields.For this analysis, the same algorithm is deployed to identify and characterize the weddies, as used to identify and extract superstructures in the u-field (refer §2.2 and the supplementary document).Figures 12(a,b) represent the same w| SS and w| noSS fields as in figures 3(e,g), but only consider motions with strong fluctuations (i.e.|w SS |, ).This threshold is based on Dennis & Nickels (2011) and assists with identification and extraction of individual, energetic w-eddies.Figures 12(c,d) present the probability distribution functions of the streamwise extents of the w-eddies identified within w| SS and w| noSS flow fields, extracted across all three PIV datasets.The pdf s confirm that the streamwise extent of w-eddies is limited to 3δ across both w| SS and w| noSS .This wavelength range closely corresponds with the geometrically self-similar hierarchy of eddies exhibiting distance-from-thewall scaling in figure 9, reaffirming the key finding of this study, based on direct analysis of the physical flow field.Although not shown here, a similar analysis on the Reynolds shear stress-carrying eddies also yields the same conclusion, reinforcing our earlier statements on the interpretation of k x φ + uw | −u ss ,+u ss .The analysis also confirms that energetic w-eddies do not physically extend along x, as long as the superstructures (> 3δ), meaning the only possible way of observing a physically long w-feature is when the individual w-eddies align along the x-direction.Indeed, the w| SS flow field indicates a much more closely-packed/clustered organization of the individual w-eddies, compared to w| noSS , in the x-direction (figure 12a,b).However, since the present analysis uses snapshot 2-D PIV data, this study cannot definitively comment on the dynamics associated behind the formation of superstructures.But, the conditional analysis presented in this section does lend empirical support in favour of the formation of superstructures, via streamwise concatenation of the intermediate-scaled eddies (Adrian et al. 2000).Interested readers are referred to the supplementary document, where we have utilized a simplified coherent structurebased model (i.e. the attached eddy model), to demonstrate a statistically plausible scenario of self-similar eddies aligning in x to 'form' a superstructure.
The present results are also consistent with the conclusions of Lozano-Durán et al. (2012), who found large-scale Reynolds shear stress-carrying structures to be essentially a concatenation of smaller uw-carrying eddies, having lengths ∼3 times their height.This also clarifies the contradiction in the literature on the 'active'/'inactive' status of the very-large-scale u-motions (i.e.superstructures).Given there are no very-large-scale w-(and consequently uw-) features in the instantaneous flow, the superstructures are indeed inactive as per the definition of Townsend (Deshpande et al. 2021a).Present evidence indicates that the superstructures comprise of several z-scaled w-carrying (i.e.active) motions, which explains the past empirical observations of superstructures carrying significant Reynolds shear stress.

Concluding remarks
The present study analyzes large-scale PIV datasets, acquired in moderate to high Re τ turbulent boundary layers (2500 Re τ 7500), to investigate the constituent motions of the turbulent superstructures.Considering that superstructures are statistically significant only at Re τ 2000 (Hutchins & Marusic 2007), the present datasets (providing sufficient scale separation) are ideally suited to identify superstructures and analyze their constituent motions.These unique datasets accurately capture the inertia-dominated instantaneous u-and w-fluctuations across a large streamwise wall-normal plane, extending up to 12δ in the x-direction.This facilitates a comprehensive investigation of the horizontal (via 1-D spectra) as well as vertical coherence (via two-point correlations) of the Reynolds shear stress-carrying eddies coexisting in the log-region, which are responsible for the momentum transfer in a high Re τ boundary layer (Baidya et al. 2017;Deshpande et al. 2021a).The statistics bring out the geometric self-similarity of these energetically significant eddies, which complements the well-established knowledge on the self-similarity exhibited by the wall-parallel velocity components in a canonical flow (Baars et al. 2017;Hwang & Sung 2018;Deshpande et al. 2020).We note that this motivates undertaking similar investigations of the momentum and heat flux in thermally stratified wall-bounded flows at high Re τ (for example atmospheric boundary layers), which can likely assist with coherent structure-based modelling of these practically relevant flows.The empirically derived scaling behaviour observed from these mean statistics (spectra and correlations) provide a benchmark for comparing and contrasting with the conditionally averaged statistics, associated with the turbulent superstructures.Such conditional statistics are made possible by the large-scale PIV flow fields, which permit identification of the superstructures directly from instantaneous flow fields.These statistics present a comprehensive picture of the superstructures, in comparison to the limited information available based on modal decompositions, used often in past studies (such as Fourier filtering, etc.).Considering the ambiguity involved while interpreting the smaller constituent motions from a u-flow field, the present study adopts the approach of investigating the w-fluctuations within the superstructure region, to understand its constituent motions.Notably, the conditional streamwise wand uw-spectra exhibit the classical z-scaling (λ x = 2z; λ x = 15z) in the intermediate scale range (Baidya et al. 2017), clearly suggesting that geometrically self-similar eddies co-exist within the superstructure region (represented schematically in figure 13).The same conclusion is demonstrated through the conditional two-point wcorrelations, along the vertical direction, which also exhibit self-similar scaling similar to that noted for the mean flow.Investigations of these kinds are only possible on analyzing instantaneous flow fields, highlighting the uniqueness of the present largescale high Re τ PIV dataset.
The argument regarding the self-similar motions, as the likely constituent motions of the turbulent superstructures, is reaffirmed by analyzing the geometry and population of individual w-eddies associated with these very-large-scaled structures.The maximum streamwise extent of the energetic w-eddies was found limited to 3δ within the superstructures, similar to that noted outside a superstructure, and conforming to the self-similar hierarchy of scales.The same analysis also revealed the spatial organization of these constituent w-eddies within the superstructures, which is consistent with the streamwise concatenation argument of forming superstructures.This also helps clarify longstanding contradiction in the literature on the active/inactive behaviour of the superstructures (Guala et al. 2006;Balakumar & Adrian 2007;Wu et al. 2012).Since there are no very-large-scale w-and uw-features in the instantaneous flow, the superstructures are indeed inactive per the definition of Townsend (Deshpande et al. 2021a).However, the study finds that superstructures comprise of several z-scaled Reynolds shear stress-carrying (i.e.active) motions (Lozano-Durán et al. 2012), which explains the past empirical observations of these very-large-scaled motions carrying significant Reynolds shear stress.
The present study concludes that superstructures are an assemblage of the attached eddy hierarchy in the streamwise wall-normal plane, hinting at a well-defined spatial organization of the attached eddies.This contradicts the original hypothesis of Townsend (1976), per which attached eddies are randomly distributed in the flow domain, suggesting the need to revisit the hypothesis (this has also been tested based on synthetic flow fields and presented in the supplementary document).The present empirical findings, specifically the Re τ -invariance of the vertical coherence of inertial eddies (R ww , R uw ), can also be used to further improve coherent structure-based models, such as the attached eddy model (Marusic & Monty 2019).This is possible through extending the data-driven approach proposed recently in Deshpande et al. (2021b), by defining the geometry of the representative eddies based on the leastsquares fits presented in (3.2).The present findings would also benefit the attached eddy model by acting as empirical evidence, to model superstructures as clusters of self-similar (attached) eddies, organized along the streamwise direction.

Figure 1 :
Figure 1: Premultiplied spectra of the streamwise velocity (kxφ + uu ) plotted against viscous-scaled wavelength (λ + x ) and distance from the wall (z + ) for a turbulent boundary layer at Reτ ≈ 7300 (Hutchins & Marusic 2007).× and marked in the plot correspond to the 'inner' and 'outer' peaks of the u-spectrogram noted previously in the literature.Regions (I), (II) and (III) are used to indicate spectral signatures of various coherent motions observed in the literature.Region (I) corresponds to the near-wall cycle captured via flow visualization by Prof. S. J. Kline (photo shared by Prof. D. Coles).Region (II) corresponds to the LSMs (conceptual sketch by Adrian et al. 2000), UMZs (particle image velocimetry (PIV) by de Silva et al. 2016) and attached eddies (attached eddy simulations by de Silva et al. 2016).Region (III) corresponds to the VLSMs/superstructures, visualized via time resolved PIV by Dennis & Nickels (2011).

Figure 2 :
Figure 2: Colour contours of the instantaneous (a,b) streamwise, u and (c) wall-normal velocity fluctuations, w in a boundary layer at Reτ ≈ 2000.This data has been extracted from a particular 3-D time block of the publicly available DNS dataset of Sillero et al. (2013).In (a), u is plotted on a wall-parallel plane at z ≈ 0.05δ, as well as on the cross-planes at x ≈ 2δ and 4δ.(b) and (c) respectively plot the u and w fluctuations in the streamwise wall-normal plane shaded in grey in (a).Dashed black line in (b,c) traces the top part of a long −u ramp type structure.(d,e) respectively plot an idealized distribution of u and w flow field induced by multiple prograde vortices (in green) positioned along the ramp (Adrian et al. 2000; de Silva et al. 2016).

Figure 3 :
Figure 3: Schematic of the experimental setup used to conduct LFOV PIV experiments in the streamwise wall-normal plane (x-z) in the HRNBLWT.Green shading indicates flow illuminated by the laser while the orange shading indicates the flow field cumulatively captured by the PIV cameras (shown in the background).Dash-dotted black line represents the streamwise evolution of the boundary layer thickness, with δ defined at the centre of the full flow field.(b,c) Instantaneous (b) u + and (c) w + -fluctuations from the LFOV PIV dataset at Reτ ≈ 2500.The dashed green box in (b,c) identifies a low-momentum turbulent superstructure (−uss) of length Lx based on the superstructure extraction algorithm described in §2.2.(d,e) show an expanded view of the u-and w-fluctuations within −uss, as identified in (b,c), respectively.Alternatively, the dashed brown box in (b,c) represents flow field of the same length×height as the dashed green box, but not associated with a turbulent superstructure (noSS).(f,g) show an expanded view of the u-and w-fluctuations within the noSS region identified in (b,c), respectively.

Figure 4 :
Figure 4: Probability density function (pdf ) of the lengths of the large and intense, (a) low and (b) high streamwise momentum motions detected by the superstructure extraction algorithm in PIV flow fields of various Reτ .Background shading indicates the bin sizes used to estimate the pdf, for which the total number of detected superstructures (i.e.addition of +uss and −uss) was used for normalization.Empty symbols indicate zero probability for the respective bin.

Figure 5 :
Figure 5: Iso-contours of the premultiplied streamwise 1-D (a) energy spectra of w-fluctuations and (b) co-spectra of the Reynolds shear stress plotted against z + and λ +x , computed from the LFOV PIV dataset at various Reτ .Dash-dotted golden and magenta lines represent the relationships λx ≈ 2z and λx ≈ 15z, respectively followingBaidya et al. (2017).(c) Schematic of representative w and uw-carrying eddies centred at various distances from the wall (zr) in the log region, with light to dark shading used to suggest an increase in zr.Rww(z/zr) and Ruw(z/zr) respectively represent the vertical coherence of the w-and uw-carrying eddies centred at zr.

Figure 6 :
Figure6: (a-d) Cross-correlation of w-fluctuations measured at z and zr, normalized by w 2 (zr) for various zr the log-region.(e-h) Cross-correlation between u(zr) and w(z), normalized by uw(zr) for the same zr as in (a-d). .(a,b,e,f) are estimated from the LFOV PIV datasets while (c,g) have been computed from the PIV case of the SLTEST dataset in zr listed in the legend corresponds to the 9 th , 8 th , 6 th and 4 th sonic positioned from the ground.Dashed green line corresponds to the linear relationship, z/zr while dash-dotted golden and teal lines correspond to the least-squares fit R a ww and R a uw defined in (3.2), respectively.

Figure 7 :
Figure 7: (a-f) Premultipled 1-D spectra of the u-fluctuations plotted versus λx/δ at (a-c) z + ≈ 2.6√ Reτ and (d-f) z + ≈ 0.5Reτ for LFOV PIV data at Reτ ≈ (a,d) 2500, (b,e) 5000 and (c,f) 7500.Dashed black lines correspond to the mean spectra obtained by ensembling across 3000 PIV images of the full flow field.While, the solid blue and red lines represent conditional spectra computed from the u-flow field extracted based on identification of a -uss and +uss, respectively.The spectra in green is computed by ensembling across both -uss and +uss.

Figure 8 :
Figure 8: Premultiplied 1-D spectra of the u-fluctuations plotted versus λx/δ at z + ≈ 2.6√ Reτ for Reτ ≈ (a) 2500, (b) 5000 and (c) 7500.The mean spectra estimated from the full flow field (in black lines) is ensembled across all 3000 fields.While, the conditional spectra corresponds to extracted flow fields (∼ 300) of the same length×height associated (in green) and not associated (in brown) with the superstructures.

Figure 11 :
Figure 11: (a-f) Conditionally averaged correlations between w-fluctuations at z and zr, normalized by w 2 (zr) for various zr.The correlations have been computed from the extracted w-flow fields associated with both -uss and +uss.Dashed black line corresponds to the linear relationship, z/zr while dashed dotted golden line corresponds to R a ww defined in (3.2).

Figure 12 :
Figure 12: (a,b) Examples of intense w-fluctuations (|wSS|, |wnoSS| > 1.3 w 2 (z)) present within flow fields associated (a) with superstructures (w|SS) and (b) not associated with superstructures (w|noSS).The w|SS and w|noSS flow fields used as examples in (a,b) essentially correspond to the extracted fields shown in figures 3(e,g).(c,d) Probability density function (pdf ) of the lengths (L w x ) of intense (c) w|SS and (d) w|noSS motions extracted from the corresponding flow fields at various Reτ .Background shading indicates the bin size used to estimate the pdf , for which the total number of detected w|SS and w|noSS were used for normalization.Empty symbols indicate zero probability for the respective bin.

Figure 13 :
Figure 13: Conceptual representation of the main conclusion of this study: z-scaled eddies are likely the constituent motions forming the turbulent superstructures.

Table 1 :
Table summarizing details of datasets comprising synchronized measurements of u-and w-fluctuations at various wall-normal locations.Reτ for the various PIV datasets is based on δ estimated at the centre of the flow field (figure