A self-sustaining process model of inertial layer dynamics in high Reynolds number turbulent wall flows
Abstract
Field observations and laboratory experiments suggest that at high Reynolds numbers Re the outer region of turbulent boundary layers self-organizes into quasi-uniform momentum zones (UMZs) separated by internal shear layers termed ‘vortical fissures’ (VFs). Motivated by this emergent structure, a conceptual model is proposed with dynamical components that collectively have the potential to generate a self-sustaining interaction between a single VF and adjacent UMZs. A large-Re asymptotic analysis of the governing incompressible Navier–Stokes equation is performed to derive reduced equation sets for the streamwise-averaged and streamwise-fluctuating flow within the VF and UMZs. The simplified equations reveal the dominant physics within—and isolate possible coupling mechanisms among—these different regions of the flow.
This article is part of the themed issue ‘Toward the development of high-fidelity models of wall turbulence at large Reynolds number’.
1. Introduction
Turbulent wall flows, including the canonical flat-plate turbulent boundary layer (BL), exhibit quasi-coherent flow structures on a wide range of spatio-temporal scales. In the near-wall region, i.e. the viscous sublayer and buffer layer, the interaction among space-filling streamwise vortices and streaks has been shown to give rise to a regeneration cycle capable of locally sustaining the turbulence on scales of approximately 100 viscous units ν/u_{τ} in extent, where ν is the kinematic viscosity and u_{τ} is the wall friction velocity [1]. There is a large body of literature documenting the occurrence and properties of these near-wall coherent structures in laboratory experiments [2] and direct numerical simulations (DNS) [3] and establishing quantitative theories for their self-sustenance [4,5,6,7,8,9,10]. At large values of the friction Reynolds number Re_{τ}≡u_{τ}h/ν, where h is an appropriate outer length scale such as the BL thickness, quasi-coherent flow structures also arise in the outer region of the BL (i.e. outside the buffer layer), where the long-time mean dynamics are inertially dominated. These structures, too, have been observed and studied in laboratory experiments and DNS. Moreover, Hwang et al. [11,12] recently have employed over-damped large-eddy simulations strategically designed to isolate large-scale structures away from the wall for Re_{τ}≈100–1000. To date, however, quantitative theories—systematically developed from the governing incompressible Navier–Stokes (NS) equation—that describe the essential nonlinear dynamics by which self-organized flow features in the inertial layer are maintained have not been developed to the extent that is true for their near-wall counterparts. One exception is the large Reynolds number asymptotic theory developed by Deguchi & Hall [13], who demonstrate that coherent vortical structures (nearly) in the free stream of laminar BL flows that approach their free-stream speeds exponentially with distance from the wall can be self-sustained independently of near-wall excitation. (In fact, these authors show that the free-stream structures can drive strong streaky flows adjacent to the wall.) Nevertheless, their theory does not specifically account for flow features observed within the inertial layer, the subject of interest here.
Such quantitative theories are desirable for at least two reasons. First, they provide a framework for understanding the mechanisms controlling the formation and evolution of flow structures within the inertial layer, which ultimately may facilitate high-fidelity numerical simulations of wall-bounded flows at extreme values of Re_{τ} (>10^{4}) by guiding the development of customized, e.g. multi-scale, numerical algorithms. Second, elucidation of the intrinsic nonlinear dynamics and structure of the inertial layer is a necessary prerequisite for increased understanding of ‘inner–outer’ interactions, understanding that may be leveraged for the design of improved flow control strategies. As an initial step towards the development of a first-principles quantitative theory, we propose in this investigation a simplified model for a self-sustaining process (SSP) that may support quasi-coherent structures away from the wall in turbulent shear flows at extreme values of the Reynolds number. The model, which provides a conceptual framework for interpreting the emergence of characteristic inertial layer structures (described below), is made quantitative in the physically relevant limit using multi-scale asymptotic analysis.
Various quasi-coherent flow structures have been observed in the outer region of turbulent wall flows. These structures have length scales much larger than 100 ν/u_{τ} penetrate into or exist wholly within the inertial region, and include so-called large-scale and very-large-scale motions (LSMs and VLSMs, respectively) and superstructures [14], which exhibit spatial coherence over increasing streamwise length scales. In particular, inertial-layer superstructures, emergent flow features carrying positive and negative streamwise velocity fluctuations, have been shown to extend for roughly 5–15 h in the streamwise direction, to carry a significant fraction of the total turbulent kinetic energy of the flow, and to impress a modulational signature in the near-wall region [15]. Despite their likely significant impact on turbulent transport and inner–outer interactions, the origin of these various large-scale structures is not wholly clear. For example, it has been argued that LSMs and VLSMs arise from the spontaneous organization of attached and/or detached hairpin vortices and hairpin vortex packets [14]. Here, we investigate a complementary possibility, namely that certain large-scale quasi-coherent flow structures in the outer region of turbulent wall flows arise directly from a self-sustaining, multiple space and time-scale process, in loose analogy with the way in which near-wall streaks and vortices have been argued to arise from a single-scale (uniformly space-filling) SSP [1]. It should be emphasized that a similar thesis has been advanced recently by Hwang and co-workers [11,12,16], but here we focus specifically on developing a semi-analytical, mechanistic description that has the potential to explain the formation and maintenance of uniform momentum zones (UMZs) and interlaced ‘vortical fissures’ (VFs), arguably the primal coherent structures in the outer part of turbulent wall flows at extreme Reynolds number.
As first documented by Meinhart & Adrian [17], the wall-normal structure of the instantaneous streamwise velocity in a turbulent BL exhibits a staircase-like variation, with relatively uniform regions (UMZs) segregated by VFs across which there are discernible jumps in the streamwise velocity; that is, the VFs essentially are internal shear layers. A spate of more recent investigations [18,19] has confirmed and further quantified the basic conception of UMZs advanced in the pioneering study by Meinhart & Adrian [17]; see figure 1a,b, which is adapted from de Silva et al. [19]. These latter authors also show that the number of UMZs (and VFs) increases logarithmically with increasing Reynolds number. Klewicki [20,21] has proposed a kinematic description of the spanwise vorticity associated with this staircase UMZ/VF structure that is consistent both with the mean momentum equation and with the intriguing notion that the turbulent BL comprises logarithmically many viscous internal layers containing most of the vorticity and, hence, dissipation. In a related study, Cuevas et al. [22] demonstrate that the logarithmic profile and other statistical features associated with the long-time mean streamwise flow can be recovered by ensemble averaging staircase-like streamwise velocity profiles with semi-empirically determined distributions for the spacing between and magnitude of the velocity jumps (or ‘steps’).
Motivated by observations of the preponderance of and likely fundamental dynamical role played by UMZ/VF structures in the inertial region of turbulent wall flows, we propose, in §2, a conceptual model of a single (geometrically) idealized VF and adjacent UMZs. A large Reynolds number asymptotic analysis is presented in §3 to render this model quantitative. The analysis exploits the emergence of distinct dominant balances of forces in different subdomains of the flow and yields simplified (‘reduced’) equations that govern the flow in each region, which then must be coupled via suitable matching conditions. This rather technical section is followed in §4 by a summary of the resulting reduced model and a description of a suitable iterative solution algorithm; here, the global structure of the nonlinear couplings characterizing the underlying SSP is laid bare. (Some readers may prefer initially to skim §3 and subsequently to return to this section after reading §4.) We conclude in §5 with a brief discussion of the mathematical and physical attributes of the present SSP model, including its connection to observed BL features.
2. Conceptual model
Perhaps the most fundamental question to be addressed is why, in the first instance, should regions of quasi-uniform momentum appear in a sheared flow? To this end, we note that UMZs are observed not only in flat-plate turbulent BLs (for which, in the absence of turbulent motions, the flow indeed would be quite uniform except within an exquisitely thin laminar BL adjacent to the wall) but also in plane Poiseuille and other canonical wall flows whose laminar states are sheared on O(h) wall-normal length scales. The innovative large Reynolds number analysis by Deguchi & Hall [13] discussed in §1 exploits a presumed pre-existing shear-layer or laminar-BL structure but does not yield insights into physical processes that would be capable of inducing sharply varying velocity profiles were those profiles not supported by external agencies. The hypothesis we propose in this investigation is based on our interpretation that large-scale structures in the outer region consist of streamwise (x-directed) vortices and streaks, much like the streamwise vortices and streaks that occur in the near-wall region, with the streaks being associated with spanwise anomalies in streamwise flow speed between laterally adjacent UMZs and the streamwise vortices accounting for a comparably weak rotary flow in wall-normal/spanwise (y/z) or ‘perpendicular’ (⊥) planes within the UMZs. However, we postulate a number of crucial differences between near-wall and inertial-layer streaks and vortices, as discussed next.
First, even at asymptotically large Reynolds number, the streaks associated with lower- and also upper-branch exact coherent states (ECS) arising in the SSP theory of Waleffe and co-workers [5,23] and in the large Reynolds number limit of this theory independently developed by Hall and co-workers [6,7] (who employ the terminology ‘vortex–wave interaction’ or VWI theory; see Chini [24]), the x-averaged streaky streamwise flow varies smoothly throughout the entire spatial domain. The absence of anything akin to a UMZ is consistent with the effective Reynolds number for this x-mean flow being O(1) [7,8]. Although it has long been understood that comparably weak roll motions in the presence of O(1) dimensionless shear can induce O(1) streaks through the lift-up effect associated with roll advection (indeed, this is a cornerstone of VWI theory), a primary facet of the model we propose is that the streamwise vortices—while still weaker than the streaks—must be sufficiently strong that the effective Reynolds number for the x-averaged flow is large.
Consistent with this requirement, we demonstrate in §3a that the large-Re_{τ} limiting form of the x-averaged streamwise momentum equation reduces to a two-dimensional advection–diffusion equation (i.e. in the ⊥-plane) in the weak diffusion limit for the averaged streamwise flow, with a prescribed advecting velocity field generated by the streamwise roll motions, i.e. the streaky flow essentially behaves as a passive scalar field at large effective ‘Peclet’ number. Crucially, numerous studies of two-dimensional scalar advection by a counter-rotating cellular velocity field in the weak diffusion (large Peclet number) limit confirm that the scalar field is strongly homogenized [25,26,27], thereby providing a possible mechanism for the observed quasi-uniformity of the streamwise momentum within UMZs. Empirical support for this proposition is provided by the DNS results of Papavassiliou & Hanratty [28], which confirm that low-momentum regions of large-scale (inertial layer) structures in turbulent plane Couette flow are separated by the vortex cores of nearly inviscid streamwise roll modes, implying a large effective Peclet number for the advected streamwise flow.
Second, the kinematic configuration of vortices and streaks in the near-wall region differs from that which we propose in our simplified model of inertial-layer dynamics. Motivated by the conceptual picture of UMZs and VFs illustrated in the schematic (adapted from Priyadarshana et al. [29]) shown in figure 2a, we postulate that away from the wall a spanwise array of streamwise-aligned, counter-rotating roll vortices stacked in the wall-normal direction and exposed to a background (e.g. laminar) wall-normal shear flow will—at large effective Reynolds number—naturally induce UMZs separated by VFs (figure 2b). If realized, this configuration would thereby provide a mechanism for the formation and sustenance of the internal shear layers exhibited by the streamwise velocity field. In the absence of the straining provided by the vortices, the shear layers (and, hence, VFs) ultimately would disintegrate owing to instability and viscous diffusion. To accord with observations, the stacked vortices must maintain the VFs through an advective/diffusive balance in which the jump in the streamwise-averaged streamwise velocity is O(u_{τ}) (figure 1). Owing to this kinematic configuration and the dynamics by which it is sustained, the SSP underlying our simplified model of inertial layer motions inherently involves multiple spatial scales, also in contrast with the near-wall SSP.
Finally, the feedback mechanism by which the streamwise rolls are sustained in the near-wall SSP differs from that conceived here for their inertial-layer counterparts. For the SSP operative in the near-wall region, the rolls are maintained by nonlinear interactions among x-varying Rayleigh instabilities arising from spanwise inflections of the streaky streamwise flow [5,7]. In contrast, in the model contemplated here, the rolls within the UMZs are triggered by instabilities within the VFs, i.e. outside the UMZs, because the fissures are the sites of strong wall-normal inflections and thus are prone to Rayleigh instabilities with disturbances predominantly varying in x/y planes. Although linear stability considerations indicate that the fastest-growing disturbances have streamwise wavelengths comparable to the (dimensional) VF thickness Δ_{f}, these modes cannot directly couple with motions in the UMZs because their wall-normal extent similarly scales with Δ_{f}, implying that such disturbances are confined to the fissure. In contrast, disturbances with streamwise wavelengths long compared with the fissure thickness induce pressure fluctuations that penetrate well outside the VF. Moreover, the near-marginal instabilities most relevant to an SSP are, in fact, long-wavelength modes. As shown in §3, the x-varying streamwise velocity component associated with these long-wavelength near-marginal modes is amplified within the VF, enabling the fluctuation dynamics to nonlinearly couple with the streamwise-averaged flow inside the fissure.
3. Large Reynolds number analysis
In this section, we perform a large Reynolds number asymptotic analysis of the NS equation to derive simplified equations governing the flow configuration shown in figure 2b. In the idealization considered here, the domain is imagined to have been excised from a large-Re_{τ} flow sufficiently far from the wall that the x-averaged streamwise velocity component varies by only O(u_{τ}) across an outer wall-normal distance l_{y}<h comparable to the wall-normal spacing between adjacent VFs, except near the fissures, where the x-averaged streamwise velocity also varies by O(u_{τ}) but across a distance comparable to the fissure width Δ_{f}.
Scaling all velocity components by u_{τ} and lengths by l_{y}, the dimensionless NS equation can be expressed as
We choose a local Cartesian coordinate system moving with the presumed planar VF (a geometrical configuration perhaps more likely to be relevant to VFs in fully developed internal flows, such as plane Poiseuille flow; see Eisma et al. [18]), which is taken to be coincident with the plane y=0. The fissure, or internal shear layer, has O(Δ) dimensionless thickness, where Δ(Re)→0 as . It transpires that in addition to the internal shear layer an even thinner ‘inner’ region, a critical layer (CL), emerges with thickness O(δ), where δ=o(Δ) as (figure 3). We separately analyse the flow within the UMZs above (y>0, denoted ‘+’) and below (y<0, denoted ‘−’) the fissure, within the internal shear layer and within the CL in the limit of large Reynolds number.
(a) Uniform momentum zones
First, we analyse the flow within UMZ^{+} and UMZ^{−}. As for near-wall streaks and vortices, we anticipate that at large Re a comparably weak vortex flow of size , where as , exposed to an O(1) dimensionless wall-normal shear can induce an O(1) streaky flow: a streamwise-averaged streamwise flow that varies in the spanwise (z) direction. Thus, we posit the following expansions for the velocity and pressure fields:
field | UMZ [y=O(1)] | VF [y=O(Δ)=O(Re^{−1/5})] | CL [y=O(δ)=O(Re^{−2/5})] |
---|---|---|---|
1 | 1 | Re^{−1/5} | |
Re^{−3/5} | Re^{−4/5} | Re^{−1} | |
Re^{−3/5} | Re^{−3/5} | Re^{−3/5} | |
Re^{−6/5} | Re^{−6/5} | Re^{−6/5} | |
u′ | Re^{−7/10} | Re^{−1/2} | Re^{−1/2} |
v′ | Re^{−7/10} | Re^{−7/10} | Re^{−9/10} |
w′ | Re^{−7/10} | Re^{−7/10} | Re^{−1/2} |
p′ | Re^{−7/10} | Re^{−7/10} | Re^{−7/10} |
Substituting these expansions into the NS equation (3.1) and the incompressibility condition and performing an x/t average yields the leading-order equations for the streamwise-averaged flow within the UMZs:
Subtracting equations (3.3) from the instantaneous NS and continuity equations yields the leading-order fluctuating momentum and continuity equations within the UMZs:
(b) Internal shear layer (VF)
Given smoothly varying roll motions within the UMZs, equation (3.1) indicates that a simple advective/wall-normal-diffusive balance is attained when wall-normal gradients in the x-averaged streamwise velocity occur over a wall-normal length scale . To analyse the flow within this internal shear layer (or VF), we introduce a rescaled wall-normal coordinate , so that , and posit the following asymptotic expansions for the various fields:
Substituting expansions (3.7) into the x-component of the NS equation (3.1), collecting terms at leading order and averaging yields
The leading-order fluctuation equations within the internal shear layer are
(c) Critical layer
Presuming as (as shown to be self-consistent below), then , where δ=o(Δ) is the thickness of the CL. Matching with the mean fields in the internal shear layer requires and , while incompressibility then yields . Similarly, given the hypothesized behaviour of as , inspection of equations (3.11) shows that within the CL the fluctuation fields must have the following magnitudes: u′=O(a′/Δ), v′=O(δa′/Δ) and w′=O(Δa′/δ), while p′=O(a′). Because w′ is amplified within the CL, whereas u′ is not, three-dimensional incompressibility can be restored within this layer by requiring δ=Δ^{2}.
The thicknesses of both the critical and internal shear layers can now be determined by considering the dominant balance of terms in the fluctuating spanwise momentum equation. To heal the singularity in as , wall-normal diffusion must be comparable to x-advection within the CL, yielding the relation Δ=δ^{3} Re. Substituting δ=Δ^{2} then gives the thickness of the internal shear layer Δ=Re^{−1/5} and, thence, the thickness of the CL δ=Re^{−2/5}. Recalling that the mean advective/diffusive balance in the internal shear layer requires , we deduce that the roll amplitude . Finally, the amplification of w′ within the CL renders the RSD component sufficiently large to drive an x-averaged flow in the spanwise direction provided that the size a′ of the fluctuations is tuned appropriately. The key observation is that because the CL is internal to the flow (i.e. away from any walls) matching of the x-mean x-vorticity between the VF and the CL must be enforced. Consequently, within the CL, is a factor of δ/Δ smaller than the expected estimate of . Using this refined scaling, it is then straightforward to show that the RSD component balances wall-normal diffusion of mean spanwise momentum within the CL if a′=Re^{−7/10}. Thus, , confirming that the RSD driving the rolls is concentrated within the CL. Table 1 summarizes the implied scalings of the various fields in each of the three subdomains.
Denoting the inner fields with capital letters, e.g. u=U(x,Y,z,τ,T;Re), where Y ≡Re^{2/5}y and the intermediate time scale τ≡Re^{−1/5}t, the leading-order mean equations within the CL reduce to
The leading-order fluctuation equations within the CL simplify to
Restricting attention to equilibrium ECS, we Fourier transform in x the steady version of the third of equations (3.14) to obtain an inhomogeneous Airy-like equation for the Fourier-transformed spanwise fluctuation velocity component ,
4. Summary of reduced model
The large-Re analysis described in §3 suggests the possibility of an inertial-layer SSP in which rolls within the UMZs are driven by the RSD acting in a CL of thickness δ=Re^{−2/5} that itself is embedded within an internal shear layer (i.e. a VF) of thickness Δ=Re^{−1/5}. The O(Re^{−3/5}) counter-rotating rolls, stacked in the wall-normal direction, are sufficiently strong to differentially homogenize the O(1) ambient shear and maintain the VF. Once the mean streamwise momentum is homogenized within the UMZs, the x-varying fluctuations may be expected to be largely irrotational (as the generalized Rayleigh equation, equation (3.6), reduces to Laplace’s equation for spatially uniform ). It should also be noted that, in accord with the Prandtl–Batchelor theorem for steady two-dimensional cellular flows at large (effective) Reynolds number [32], the x-mean x-directed vorticity also may be expected to be largely uniform within the UMZs. The key elements of the proposed SSP are summarized in figure 6.
To demonstrate the closure of this feedback loop, we next outline an algorithm to compute the global flow structure. We begin by discussing the determination of the (x-Fourier-transformed) fluctuation pressure within UMZ^{+} (say) via the solution of the generalized Rayleigh equation, equation (3.6), subject to appropriate boundary conditions. We seek solutions for the fluctuating pressure in the UMZ that are 2π/β-periodic in the z direction and that match with the fluctuating pressure in the internal shear layer. The matching condition on the pressure fluctuation can be enforced indirectly by matching the O(Re^{−7/10}) fluctuating wall-normal velocity components. Evaluating the y-component of the fluctuation momentum equations in UMZ^{+} as y→0^{+}, we obtain (in steady state)
From boundary condition (4.1), the magnitude of the Fourier-transformed pressure fluctuation is proportional to some scalar (i.e. z-independent) measure of the amplitude function A(z) arising in equations (3.11), which we take to be
As noted previously, the mean streamwise flow will be strongly homogenized (see the discussion in §2) except within asymptotically thin, spanwise-localized jet-like regions of positive and negative anomolies in the mean streamwise speed (and, of course, within the fissures). Thus, we simplify the following discussion of the global solution algorithm by setting (in UMZ^{+}), i.e. a constant known from symmetry considerations—a more refined analysis would properly treat the variability of within the UMZs. For brevity of notation, we omit the subscripts from the field variables in the following algorithm.
For a given spanwise domain length 2π/β and scalar fluctuation amplitude :
(i) Generate an initial iterate for λ(z), the streak-induced shear within the CL.
(ii) Obtain an integral representation of the solution to Laplace’s equation
to which the generalized Rayleigh equation (equation (3.6)) reduces for uniform , for within the UMZ^{+} domain y∈(0^{+},1^{−}) and z∈[0,2π/β], subject to the boundary conditions (see the discussion surrounding equation (4.1))Denote this solution, which is a functional of A(z) and depends on α, by .(iii) Solve the ordinary differential eigenvalue problem for A(z) and α obtained by setting
and by imposing the normalization constraint given by equation (4.2), where is determined via solution of equation (4.3) and is a functional of A(z) through , subject to a 2π/β-periodicity requirement in z. ( can be related to A(z) via further manipulation and partial solution of equations (3.14).)(iv) Time-advance to steady state the mean equations for the ⊥-flow within the VF. These equations (equation (3.9), noting that the RSD term vanishes; incompressibility; and the -independence of the mean pressure) can be cast in stream function/vorticity form,
where and , and must be solved subject to the symmetry and matching conditionsrespectively. Note that the solution for in UMZ^{+} is readily obtained up to a multiplicative constant via the solution of a Poisson equation for the UMZ stream function ψ(y,z), , where the right-hand side is constant owing to the Prandtl–Batchelor constraint, subject to ψ→0 around the periphery of each cell. We leave for future investigation the determination of this constant value of the core vorticity but anticipate that the imposition of global constraints, including a global energy balance and/or constancy of the wall-normal flux of mean streamwise momentum, may be useful [27].(v) Time-advance to steady state the mean equation for the streaky flow within the VF (equation (3.8)),
subject to the symmetry and matching conditions at and as .(vi) Evaluate the streak shear within the CL, , and return to step (ii) until convergence is achieved to within a desired tolerance.
5. Conclusion
Motivated by the occurrence of zones of quasi-uniform streamwise momentum (UMZs) separated by thin layers of concentrated spanwise vorticity (VFs) in the outer region of turbulent wall flows at extreme values of Re_{τ}, we have performed a large Reynolds number asymptotic analysis of the incompressible NS equation to investigate whether an SSP can directly support such flow structures. This possibility should be contrasted with the notion that large-scale, outer-region structures result from the spontaneous concatenation of smaller-scale flow features. In this work, we have taken an important step towards confirming the former scenario by deriving from the instantaneous NS equation a self-consistent reduced model of a single planar VF coupled to adjacent UMZs above and below the fissure. The analysis suggests a three-region asymptotic structure in the wall-normal direction, with: outer UMZs in which the streamwise-averaged (‘mean’) and streamwise-varying (‘fluctuation’) dynamics are largely inviscid; an internal shear layer across which the mean streamwise flow exhibits a sharp, step-like variation; and an even thinner CL (within the shear layer) in which the spanwise fluctuating velocity component is amplified. The wall-parallel gradient of the streamwise self-correlation of this velocity component drives a mean spanwise flow within the CL; this spanwise flow is viscously transmitted to the internal shear layer and, by continuity (and further viscous action), ultimately transformed into roll motions outside the shear layer. The induced rolls are weak compared with the background ambient shear (i.e. in dimensional terms, roll velocities are small compared with the friction velocity u_{τ}) but, unlike prior large Reynolds number asymptotic studies of ECS and the associated SSPs, the rolls are sufficiently strong to differentially homogenize the resulting streaky streamwise flow. The straining provided by the counter-rotating rolls, which are stacked in the wall-normal direction, reinforces the VF and prevents its diffusive disintegration. Of course, to verify the proposed SSP scenario, solutions of the reduced system of equations must be numerically constructed following the algorithm outlined in §4 (and, ideally, compared with full NS solutions). We hope to report on these computations in a subsequent publication.
Detailed comparisons of the proposed theory with measurements of turbulent flows from laboratory experiments are not yet feasible. Nevertheless, we can compare the predicted scaling of the VF thickness with that extracted from laboratory and field data. In particular, Klewicki [20] estimates that the dimensional fissure thickness Δ_{f} normalized by the outer length scale (labelled h in this study) varies in direct proportion to using data that span more than two decades in Re_{τ}. By contrast, our theory indicates that the corresponding dimensionless thickness of the internal shear layer Δ=O(Re^{−1/5}), where Re=(l_{y}/h) Re_{τ}; i.e. the theory evidently predicts a thicker fissure than has been observed. One potential explanation for this discrepancy is the necessity of arbitrarily choosing thresholds for the quantification of Δ_{f} from raw empirical data. It is also plausible that distinct large-Re reduced models, based on different scalings, may yield a fissure thickness that accords with the measured Δ_{f}/h. Even in that circumstance, we conjecture that elements of the global flow structure predicted by the present analysis still may be realized (e.g. the homogenizing action of the rolls). It is also worth noting that the present theory does, in fact, predict an inner region of thickness O(Re^{−2/5}), but the streamwise-averaged streamwise velocity is not sheared over this length scale. Finally, the conclusion reached above regarding the apparent discrepancy in the predicted and measured fissure thicknesses tacitly presumes the constancy of the length-scale ratio l_{y}/h when, in fact, this ratio appears to be a weak function of Re_{τ}.
More qualitatively, we comment that the UMZ/VF flow structure analysed here may constitute a ‘straight’ section of a longer superstructure that slowly meanders in the streamwise direction, and that the coupling mechanisms elucidated herein have the potential both to complement those described by the modelling framework of Sharma et al. [35] and to connect to the features observed in the numerical and physical experiments of Cossu & Hwang [36] and Baars et al. [37], respectively. In contrast to the recent work by Hwang et al. [11,12,16] discussed in §1, this study captures only a single interacting UMZ/VF ‘unit’. Nevertheless, based on the self-similar hierarchical rescaling of VWI states performed by Blackburn et al. [38] and of ‘resolvent modes’ by Moarref et al. [39], we speculate that the flow structures investigated here may constitute one member of a family of similar motions on a hierarchy—the varying spacing between VFs with distance from the wall will naturally generate a hierarchy of scales.
We conclude by commenting again on the relationship between the more well-studied near-wall and the newly proposed outer-region SSPs and associated ECS. In the near-wall region, the streamwise-averaged flow structures, i.e. rolls and streaks, are space filling and single scale even as . By contrast, the proposed SSP supports a multi-scale flow structure away from the wall in that limit. Continuing with this dichotomy, it is helpful to recall that Hall & Sherwin [7] describe the near-wall SSP they study as arising from a novel variant of ‘steady streaming’ within a CL, i.e. the single-wavenumber fluctuation/fluctuation nonlinearity within the CL drives a mean spanwise component of the roll flow. Because the effective Reynolds number felt by the mean flow is O(1), the streaming could be described more precisely as being of the Rayleigh–Nyborg–Westervelt variety [40,41,42]. The proposed SSP also may be viewed as arising from a novel form of steady streaming—but in the limit in which the effective Reynolds number is large. This so-called ‘Stuart’ streaming [43] necessarily involves the emergence of nested boundary or, here, internal layers. In summary, the self-sustaining interaction investigated here clearly has important dynamical connections both with the VWI exhibited by specific ECS in Couette flow and analysed by Hall & Sherwin [7] and with a vortex/Tollmien–Schlichting wave interaction arising in plane Poiseuille flow and analysed by Dempsey et al. [34]. Indeed, the present self-sustaining process theory blends in a novel fashion elements of each of these analyses, while also including an irreducible and essentially inviscid ‘core’ dynamics that accords with the behaviour observed in the outer region of turbulent wall flows at extreme values of Re_{τ}.
Authors' contributions
G.P.C. conceived of and with B.M. carried out the asymptotic analysis and computations and drafted the manuscript. J.K. conceived of the study and with C.M.W. participated in the model development. All authors read and approved the manuscript.
Competing interests
The authors declare that they have no competing interests.
Funding
This material is based upon work supported by the National Science Foundation under grant no. CBET-1545564. J.K. also acknowledges the support of the Australian Research Council under grant no. DP150102593.
Acknowledgements
G.P.C. acknowledges useful discussions with W. Young and N. Balmforth at an early stage of this work.
References
- 1
Hamilton JM, Kim J, Waleffe F . 1995 Regeneration mechanisms of near-wall turbulence structures. J. Fluid Mech. 287, 317–348. (doi:10.1017/S0022112095000978) Crossref, Web of Science, Google Scholar - 2
Kline SJ, Reynolds WC, Schraub FA, Rundstadler PW . 1967 The structure of turbulent boundary layers. J. Fluid Mech. 30, 741–773. (doi:10.1017/S0022112067001740) Crossref, Web of Science, Google Scholar - 3
Jiménez J, Moin P . 1991 The minimal flow unit in near-wall turbulence. J. Fluid Mech. 225, 213–240. (doi:10.1017/S0022112091002033) Crossref, Web of Science, Google Scholar - 4
Nagata M . 1990 Three-dimensional finite-amplitude solutions in plane Couette flow: bifurcation from infinity. J. Fluid Mech. 217, 519–527. (doi:10.1017/S0022112090000829) Crossref, Web of Science, Google Scholar - 5
Waleffe F . 1997 On a self-sustaining process in shear flows. Phys. Fluids 9, 883–900. (doi:10.1063/1.869185) Crossref, Web of Science, Google Scholar - 6
Hall P, Smith F . 1991 On strongly nonlinear vortex/wave interaction in boundary layer transition. J. Fluid Mech. 227, 641–666. (doi:10.1017/S0022112091000289) Crossref, Web of Science, Google Scholar - 7
Hall P, Sherwin S . 2010 Streamwise vortices in shear flows: harbingers of transition and the skeleton of coherent structures. J. Fluid Mech. 661, 178–205. (doi:10.1017/S0022112010002892) Crossref, Web of Science, Google Scholar - 8
Beaume C, Chini GP, Julien K, Knobloch E . 2015 Reduced description of exact coherent states in parallel shear flows. Phys. Rev. E 91, 043010. (doi:10.1103/PhysRevE.91.043010) Crossref, Web of Science, Google Scholar - 9
Schoppa W, Hussain F . 2002 Coherent structure generation in near-wall turbulence. J. Fluid Mech. 453, 57–108. (doi:10.1017/S002211200100667X) Crossref, Web of Science, Google Scholar - 10
Farrell BF, Ioannou PJ . 2012 Dynamics of streamwise rolls and streaks in wall-bounded shear flow. J. Fluid Mech. 708, 149–196. (doi:10.1017/jfm.2012.300) Crossref, Web of Science, Google Scholar - 11
Hwang Y, Bengana Y . 2016 Self-sustaining process of minimal attached eddies in turbulent channel flow. J. Fluid Mech. 795, 708–738. (doi:10.1017/jfm.2016.226) Crossref, Web of Science, Google Scholar - 12
Hwang Y, Willis AP, Cossu C . 2016 Invariant solutions of minimal large-scale structures in turbulent channel flow for Re_{τ} up to 1000. J. Fluid Mech. 802R1, 1–13. (doi:10.1017/jfm.2016.470) Crossref, Web of Science, Google Scholar - 13
Deguchi K, Hall P . 2014 Free-stream coherent structures in parallel boundary-layer flows. J. Fluid Mech. 752, 602–625. (doi:10.1017/jfm.2014.282) Crossref, Web of Science, Google Scholar - 14
Kim KC, Adrian RJ . 1999 Very large-scale motion in the outer layer. Phys. Fluids 11, 417–422. (doi:10.1063/1.869889) Crossref, Web of Science, Google Scholar - 15
Hutchins N, Marusic I . 2007 Evidence of very long meandering features in the logarithmic region of turbulent boundary layers. J. Fluid Mech. 579, 1–28. (doi:10.1017/S0022112006003946) Crossref, Web of Science, Google Scholar - 16
Hwang Y, Cossu C . 2010 Self-sustained process at large scales in turbulent channel flow. Phys. Rev. Lett. 105, 044505. (doi:10.1103/PhysRevLett.105.044505) Crossref, PubMed, Web of Science, Google Scholar - 17
Meinhart CD, Adrian RJ . 1995 On the existence of uniform momentum zones in a turbulent boundary layer. Phys. Fluids 7, 694. (doi:10.1063/1.868594) Crossref, Web of Science, Google Scholar - 18
Eisma J, Westerweel J, Ooms G, Elsinga G . 2015 Interfaces and internal layers in a turbulent boundary layer. Phys. Fluids 27, 055103. (doi:10.1063/1.4919909) Crossref, Web of Science, Google Scholar - 19
de Silva CM, Hutchins N, Marusic I . 2016 Uniform momentum zones in turbulent boundary layers. J. Fluid Mech. 786, 309–331. (doi:10.1017/jfm.2015.672) Crossref, Web of Science, Google Scholar - 20
Klewicki JC . 2013 A description of turbulent wall-flow vorticity consistent with mean dynamics. J. Fluid Mech. 737, 176–204. (doi:10.1017/jfm.2013.565) Crossref, Web of Science, Google Scholar - 21
Klewicki JC . 2013 On the singular nature of turbulent boundary layers. Proc. IUTAM 9, 69–78. (doi:10.1016/j.piutam.2013.09.007) Crossref, Google Scholar - 22
Cuevas J, Ebadi A, White C, Chini GP, Klewicki JC . 2016 A simple model of inertial layer dynamics in turbulent boundary layers. In Bull. 69th Annu. Meeting of the APS Division of Fluid Dynamics, Portland, OR, 20–22 November 2016, vol. 61, no. 20. College Park, MD: American Physical Society. Google Scholar - 23
Wang J, Gibson J, Waleffe F . 2007 Lower branch coherent states in shear flow: transition and control. Phys. Rev. Lett. 98, 204501. (doi:10.1103/PhysRevLett.98.204501) Crossref, PubMed, Web of Science, Google Scholar - 24
Chini GP . 2016 Exact coherent structures at extreme Reynolds number. J. Fluid Mech. 794, 1–4. (doi:10.1017/jfm.2016.154) Crossref, Web of Science, Google Scholar - 25
Childress S, Soward AM . 1989 Scalar transport and alpha-effect for a family of cat’s-eye flows. J. Fluid Mech. 205, 99–133. (doi:10.1017/S0022112089001965) Crossref, Web of Science, Google Scholar - 26
Rhines PB, Young WR . 1983 How rapidly is a passive scalar mixed within closed streamlines? J. Fluid Mech. 133, 133–145. (doi:10.1017/S0022112083001822) Crossref, Web of Science, Google Scholar - 27
Chini GP, Cox SM . 2009 Large Rayleigh number thermal convection: heat-flux predictions and strongly nonlinear solutions. Phys. Fluids 21, 083603. (doi:10.1063/1.3210777) Crossref, Web of Science, Google Scholar - 28
Papavassiliou DV, Hanratty TJ . 1997 Interpretation of large-scale structures observed in a turbulent plane Couette flow. Int. J. Heat Fluid Fl. 18, 55–69. (doi:10.1016/S0142-727X(96)00138-5) Crossref, Web of Science, Google Scholar - 29
Priyadarshana PJA, Klewicki JC, Treat S, Foss JF . 2007 Statistical structure of turbulent-boundary layer velocity–vorticity products at high and low Reynolds numbers. J. Fluid Mech. 570, 307–346. (doi:10.1017/S0022112006002771) Crossref, Web of Science, Google Scholar - 30
Hall P, Horseman NJ . 1991 The linear inviscid secondary instability of longitudinal vortex structures in boundary layers. J. Fluid Mech. 232, 357–375. (doi:10.1017/S0022112091003725) Crossref, Web of Science, Google Scholar - 31
Balmforth N, Del-Castillo-Negrete D, Young WR . 1997 Dynamics of vorticity defects in shear. J. Fluid Mech. 333, 197–230. (doi:10.1017/S0022112096004119) Crossref, Web of Science, Google Scholar - 32
Batchelor GK . 1956 On steady laminar flow with closed streamlines at large Reynolds number. J. Fluid Mech. 1, 177–190. (doi:10.1017/S0022112056000123) Crossref, Web of Science, Google Scholar - 33
Mantic-Lugo V, Arratia C, Gallaire F . 2014 Self-consistent mean flow description of the nonlinear saturation of the vortex shedding in the cylinder wake. Phys. Rev. Lett. 113, 084501. (doi:10.1103/PhysRevLett.113.084501) Crossref, PubMed, Web of Science, Google Scholar - 34
Dempsey LJ, Deguchi K, Hall P, Walton AG . 2016 Localized vortex/Tollmien–Schlichting wave interaction states in plane Poiseuille flow. J. Fluid Mech. 791, 97–121. (doi:10.1017/jfm.2016.50) Crossref, Web of Science, Google Scholar - 35
Sharma AS, Moarref R, McKeon BJ . 2017 Scaling and interaction of self-similar modes in models of high Reynolds number wall turbulence. Phil. Trans. R. Soc. A 375, 20160089. (doi:10.1098/rsta.2016.0089) Link, Web of Science, Google Scholar - 36
Cossu C, Hwang Y . 2017 Self-sustaining processes at all scales in wall-bounded turbulent shear flows. Phil. Trans. R. Soc. A 375, 20160088. (doi:10.1098/rsta.2016.0088) Link, Web of Science, Google Scholar - 37
Baars WJ, Hutchins N, Marusic I . 2017 Reynolds number trend of hierarchies and scale interactions in turbulent boundary layers. Phil. Trans. R. Soc. A 375, 20160077. (doi:10.1098/rsta.2016.0077) Link, Web of Science, Google Scholar - 38
Blackburn HM, Hall P, Sherwin SJ . 2013 Lower branch equilibria in Couette flow: the emergence of canonical states for arbitrary shear flows. J. Fluid Mech. 726R2, 1–12. (doi:10.1017/jfm.2013.254) Google Scholar - 39
Moarref R, Sharma AS, Tropp JA, McKeon BJ . 2013 Model-based scaling of the streamwise energy density in high-Reynolds-number turbulent channels. J. Fluid Mech. 734, 275–316. (doi:10.1017/jfm.2013.457) Crossref, Web of Science, Google Scholar - 40
Rayleigh L . 1883 On the circulation of air observed in Kundt’s tubes and some allied acoustical problems. Phil. Trans. R. Soc. 75, 1–21. (doi:10.1098/rspl.1883.0075) Google Scholar - 41
Nyborg WL . 1953 Acoustic streaming due to attenuated plane waves. J. Acoust. Soc. Am. 25, 68–75. (doi:10.1121/1.1907010) Crossref, Web of Science, Google Scholar - 42
Westervelt PJ . 1953 The theory of steady rotational flow generated by a sound field. J. Acoust. Soc. Am. 25, 60–67. (doi:10.1121/1.1907009) Crossref, Web of Science, Google Scholar - 43
Stuart JT . 1966 Double boundary layers in oscillatory viscous flow. J. Fluid Mech. 24, 673–687. (doi:10.1017/S0022112066000910) Crossref, Web of Science, Google Scholar