1. Introduction
The wake phenomenon is widespread in the nature, and its prevalent examples are environmental flows and flows occurring in many engineering applications. The wake is the highly turbulent region behind the bodies laid on the flat surfaces formed due to flow separation. The turbulent properties of the cylinder wake are a well-researched topic in fluid dynamics. Many investigations are available in the past literature [
1,
2,
3,
4,
5,
6,
7] to study the cylinder wake. Djeridi et al. [
1] conducted experimental research to explore the turbulent characteristics in the flow field of a circular cylinder in the near-bed upstream and near-wake region for a high value of Reynolds number equal to 14,000. They also analyzed the mean and turbulent flow field for a moderate blockage and an aspect ratio. Konstantinidis et al. [
2] experimentally determined the mean and fluctuating velocity fields in the near wake of a circular cylinder under periodic velocity fluctuations incident on the mean flow. They described the wake in terms of the recirculation region, vortex formation region, the maximum intensity of velocity fluctuations, and vortex street’s wavelength. Braza et al. [
3] studied the flow field of a circular cylinder with a high blockage and a low aspect ratio. The flow field was explored in the starting point of the critical flow regime for a Reynolds number of 14,000. Akoz [
4] experimentally estimated the characteristics of flow structures behind the cylinder, instantaneous and time-averaged velocity field, vorticity contours, streamline topology, and Reynolds stress concentrations for the cylinder Reynolds number,
ReD between 1000 and 7000. They dictated separation point from the cylinder surface and the variation of length of the primary and secondary downstream separation regions with Reynolds numbers. Akoz and Kirkgoz [
5] investigated the turbulent flow field of a bed-mounted horizontal cylinder for Reynolds numbers between 1000 and 7000. Devi and Hanmaiahgari [
6] experimentally analyzed a wall-mounted cylinder’s mean and turbulent characteristics at a downstream location. Devi et al. [
7] conducted a comparative study for the flow field of a wall-mounted cylinder to determine the effect of gravel-bed roughness compared to sand bed roughness. Despite numerous research in this field, many issues related to the cylinder wake are still unresolved, especially in shallow flows.
In the shallow flows, the horizontal length scale is generally larger than the wall-normal length scale [
8]. Destabilizing shallow flows are caused by an abrupt change in topology, such as a bluff body over a flat surface, resulted in shallow wakes. Many researchers [
9,
10,
11] have studied the nature of shallow wakes. Akilli and Rockwell [
11] used flow visualization techniques in addition to particle image velocimetry (PIV) to study the shallow flows in the near wake region of a circular cylinder. They discussed the variation of mean velocity, Reynolds stress, vorticity, and streamline topology at three different wall-normal distances. The channel bed can significantly affect the flow characteristics of the shallow wake, and bed effects are attenuated while moving away from it towards the free surface. In the near-bed region of the flow, the bed restricts the development of the wake, diminishes the interaction between shear layers, and stabilizes the flow [
12].
Self-preservation of turbulent wakes downstream of a cylinder at high Reynolds numbers was initially reported by Townsend [
13,
14] in his elaborate investigations. Out of these wakes, the one which frequently occurs in the environment is the shallow wake. The concept of self-preservation of turbulent flows is also known as local similarity or self-similarity. If some or all statistical characteristics of a turbulent flow are a function of only scales, the flow is called self-preserving in its behaviour. In other words, a turbulent flow in its self-preserving condition is the only function of its length and velocity scales. The far-field of a turbulent free shear layer such as a wake and a jet are classic examples of self-preserving flows. Maji et al. [
15] reviewed the occurrence of von Kármán vortex streets as a function of solid volume fraction in the wake region of a group of wall-mounted vertical cylinders. Sadeque et al. [
16] analyzed the flow patterns of the near-wakes of the wall-mounted vertical cylinder in the turbulent shallow flow. They found that flow fields in the region away from the wall are similar. They also perceived that the turbulent kinetic energy and the primary Reynolds stress distributions displayed the self-preserving behaviour in moderate to deeply submerged cylinders. Balachandar et al. [
8] and Tachie and Balachandar [
17] deduced that even in a near-wake region of the flow, the self-preserving characteristics can be established in terms of the time-averaged velocity distributions by using appropriate velocity and length scales. Maji et al. [
18] analyzed the self-preserving behaviour in the interior and the wake region of the vegetation patch regarding the mean flow velocities in the streamwise and wall-normal directions and turbulence intensities, RSS, and turbulent kinetic energy (TKE). They observed that the time-averaged streamwise velocity distributions depict self-preserving characteristics for both emergent and sparse vegetation patches in their outer, interior, and wake regions when scaled by their depth-averaged value. The self-preserving behaviour of the time-averaged streamwise velocity distributions was independent of the aspect ratios and the Reynolds numbers. Similar observations were also made for the RSS profiles. In the past literature [
8,
16,
17], different velocity and length scales were used for velocity, RSS, and turbulence intensities distributions. Balachandar et al. [
8] and Tachie and Balachandar [
17] used
as a velocity scale and
as a length scale for the velocity distributions to establish the self-preserving characteristics of the wake region of a sharp-edged plate. Here,
and
are the maximum and the minimum velocities, respectively, in the wake region at a particular streamwise distance,
, and
is the wall-normal distance at which
is equal to one half of the
. Sadeque et al. [
16] used
as the velocity scale for velocity profiles for the bed-mounted cylinder. In their study, the
is equivalent to
of Balachandar et al. [
8] and Tachie and Balachandar [
17]. They also observed an excellent collapse of the velocity profiles in the near wake region of the cylinder. The velocity deficit profiles collapse on a narrow band for this scaling, representing that the wall-wake region of the pipe is self-preserving in nature. Despite that much research in this area, there are very few studies that explore the self-preserving characteristics of the wall-wake region of a bed-mounted horizontal circular pipe in the shallow flow with proper scaling for the velocity, RSS, and turbulence intensities.
Therefore, the main objective of this study is to assess the self-preserving characteristics of the wall-wake flow of a bed-mounted circular pipe in terms of streamwise velocity, RSS, and turbulence intensities for different inflow conditions. Further, the variation of the wall-normal distributions of third-order correlations along the streamwise direction in the wake region was also evaluated. For this purpose, an ADV probe was used to measure instantaneous three-dimensional (3D) velocity data. After post-processing the velocity data, it was used to determine turbulence statistics of the flow and their deficits and maximum values. Henceforth, how the wall-normal profiles of the time-averaged streamwise velocity , RSS , streamwise turbulence intensity , wall-normal turbulence intensity in non-dimensional forms and third-order correlations change in the wake region of the pipe with streamwise distance, were also documented. Additionally, the wall-normal distributions of non-dimensional velocity deficit , non-dimensional RSS deficit , non-dimensional streamwise turbulence intensity deficit and non-dimensional wall-normal turbulence intensity deficit were investigated to determine the self-preserving behaviour of velocity, RSS, and turbulence intensities, respectively. Besides this, the decay of maximum velocity deficit , maximum RSS deficit , maximum streamwise turbulence intensities deficit and maximum wall-normal turbulence intensities deficit with streamwise distance, were also evaluated. This was followed by the analysis of the streamwise growth of half velocity profile width, and half RSS profile width, .
3. Results and Discussion
The approach flow shear velocity,
, was used as the velocity scale to non-dimensionalize the turbulent characteristics and the flow depth,
, was used as the wall-normal length scale for the wall-normal distance,
, as
. The wall-normal profiles of the non-dimensional time-averaged streamwise velocity,
, in the wake region of the pipe at a different streamwise distance,
, for Run 3 and Run 4 are plotted in
Figure 3. It is clear from
Figure 3 that the streamwise velocity,
, has a negative sign in the near-bed region because of flow separation while it becomes positive in the region away from the bed. The velocity in the free-surface region is also affected by the presence of the pipe, and its magnitude is different from the corresponding approach flow velocity in the free surface region. The velocity in the free-surface region of the recirculation region is greater than the undisturbed upstream velocity at the same elevation because of lower near-bed velocities in the recirculation region which shows mass conservation. In addition, this result reveals that the flow was diverted upward and accelerated in the near-surface region due to the separation of the flow from the pipe. The maximum magnitude of the streamwise velocity is occurring near the free surface for all streamwise locations,
. The maximum magnitude of the streamwise velocity is increasing along with streamwise distance up to
and then decreasing with a further increase in
. It is noteworthy that the negative velocity is fully recovered at a streamwise distance of
8, in both the runs (Run 3 and Run 4). The velocity,
-, distributions in the wake region display an upward concavity in their profiles, and this asymmetric velocity distribution occurs due to the wall-wake effects.
The points of inflection () in the individual profiles are observed near the top level of the pipe, . In the free-surface region, the velocity - distributions are attaining almost the same velocity gradient, . The extrapolation of the streamwise velocity profiles up to the channel bed verified the preservation of the no-slip condition.
Figure 4 depicts the schematic diagram of a typical wall-normal distribution of streamwise velocity deficit,
, in the wall-wake region of the pipe along with the definition of the half velocity profile width,
. The velocity deficit,
, is the difference between the upstream undisturbed time-averaged streamwise velocity at particular elevation
and the downstream wall-wake time-averaged streamwise velocity at the same elevation. The velocity deficit,
(
z), at a specific wall-normal distance,
, can be written as
. Here,
is the upstream undisturbed time-averaged streamwise velocity at an elevation,
, while
is the downstream wall-wake time-averaged streamwise velocity at the same elevation. The half velocity profile width,
, is the wall-normal elevation from the bed at which the velocity deficit is equal to one-half of the maximum velocity deficit,
.
Figure 5a displays the variation of
against
for Run 1, Run 2, Run 3, and Run 4 in the wall-wake region of the pipe, depicting the self-preserving characteristics of the velocity deficit,
, distributions. The maximum velocity deficit,
, is used as the velocity scale, and the half velocity profile width,
, is taken as the wall-normal length scale. The velocity deficit profile is exhibiting a single peak in its distribution at
. The self-preserving characteristics of the velocity deficit profiles are found to be independent of approach flow conditions and diameter of the cylinder,
. All of the streamwise velocity profiles are showing a good collapse over a narrow band for the wall-normal domain of the study,
2, and for the complete streamwise domain of the research. Above the wall-normal domain,
2 velocity profiles are not self-preserving in nature; therefore, they are excluded from the analysis. It can also be observed that the value of non-dimensional velocity deficit is equal to unity at
, which simply implied that
is the wall-normal location of the maximum velocity deficit,
.
Figure 5b portrays the streamwise growth of non-dimensional half velocity profile width,
for Run 1, Run 2, Run 3, and Run 4 in the wall-wake region of the pipe. The diameter of pipe
is used as the wall-normal length scale. It can be observed from
Figure 5b that the half velocity profile width,
, grows consistently and gradually except for starting two points. This variation in the half velocity profile width is resulting from an increase followed by a decrease in the thickness of the wake region downstream of the cylinder. A similar observation was also made by Tachie and Balachandar [
17], with a slightly different value of streamwise distance as
, the half velocity profile width grows constantly and gradually with
. Contrary to this, Sadeque et al. [
16] found a constant decline in the half velocity profile width with streamwise distance for a bed-mounted vertical cylinder.
The streamwise decay of non-dimensional maximum velocity deficit,
, for Run 1, Run 2, Run 3, and Run 4 are demonstrated in
Figure 5c. The approach flow shear velocity,
, is taken as the velocity scale to normalize the maximum velocity deficit,
. It is perceived from
Figure 5c that the maximum velocity deficit,
, decays in the streamwise direction with
. Its decay indicates the recovery of the wall-wake flow velocity similar to its undisturbed upstream values. The flow is uniformly and quickly recovering to its undisturbed upstream values. Schmeeckle and Nelson [
16] also made similar observations.
The wall-normal profiles of the non-dimensional RSS,
in the wake region of the pipe at a different streamwise distance,
, for Run 3 and Run 4 are shown in
Figure 6. Here, Reynolds shear stress (RSS),
is expressed as
. It is visible from
Figure 6 that RSS starts with a small positive value at the bed and increases sharply with an increase in the wall-normal distance,
. It attained its peak at the top level of the pipe (
z/D ≈ 1) and decreased quickly with a further rise in
. Henceforth, the value of RSS becomes very small and almost invariant of the wall-normal distance,
in the free-surface region. The location of the peak of the RSS corresponds to the point of inflection (
) in the velocity profiles. It is understood that the peaks in RSS at the top level of the pipe are caused by the enhanced mixing of turbulence at that level.
Figure 7 illustrates the diagram of the typical wall-normal distribution of RSS deficit,
, in the wall-wake region of the pipe along with the definition of the half RSS profile width,
. The RSS deficit,
, is the difference between the downstream wall-wake RSS at an elevation,
, and upstream undisturbed RSS at the same elevation. The RSS deficit,
, at a particular elevation,
, can be written as
. Here,
is the upstream undisturbed RSS at an elevation,
, and
is the downstream wall-wake RSS at the same elevation. The half RSS profile width,
, was estimated similarly to the half velocity profile width,
. The half RSS profile width
is the wall-normal height from the bed at which the RSS deficit is equal to
.
Figure 8a expresses the variation of
against
for Run 1, Run 2, Run 3, and Run 4 in the wall-wake region of the pipe, by depicting the self-preserving characteristics in the RSS deficit (
) profiles. The self-preservation of the wall-wake is evaluated based upon the collapse of the RSS deficit profiles on a narrow band. The maximum RSS deficit,
is used as velocity scale and the half RSS profile width,
is taken as the wall-normal length scale. The RSS deficit profiles are collapsing on a narrow band for this scaling, which demonstrates that the wall-wake region of the pipe is self-preserving in nature. The self-preserving characteristics of the RSS deficit profiles are found to be independent of the approach flow conditions similar to the velocity deficit profiles. Similar to velocity profiles, all RSS profiles show good collapse with the streamwise distance,
, for the wall-normal domain of study,
2. It is noteworthy that the RSS deficit profiles display a single peak, and it occurs at
. It can also be observed from
Figure 8a that the value of non-dimensional RSS deficit
was equal to unity for
. It simply reveals the peak value of the RSS deficit occurs at
i.e.,
.
Figure 8b displays the streamwise growth of non-dimensional half RSS profile width,
, for Run 1, Run 2, Run 3, and Run 4 in the wall-wake region of the pipe. The diameter of the pipe,
D, is used as the wall-normal length scale. It can be observed from
Figure 8b that half RSS profile width,
, grows in the streamwise direction with an increase in the streamwise distance,
.
The streamwise decay of non-dimensional maximum RSS deficit,
for Run 1, Run 2, Run 3, and Run 4 are depicted in
Figure 8c. The negative sign is used before the maximum RSS deficit
, since
is greater than
. The approach flow shear velocity,
, is taken as the velocity scale to normalized the maximum RSS deficit,
.
Figure 8c demonstrates that the maximum RSS deficit,
, decays in the streamwise direction with
, which indicates the recovery of the flow with an increase in
. Similar results were also revealed by Sadeque et al. [
16], although their case is of a bed-mounted vertical cylinder. These patterns in the variation of
in the streamwise direction are happening due to the dampening of the turbulence mixing with an increase in
and recovering the flow similar to its undisturbed upstream form.
The turbulence intensity downstream of the pipe exists due to the fluctuations in the instantaneous flow velocity influenced by the pipe and bed roughness. Non-dimensional turbulence intensities in the streamwise,
and wall-normal directions,
are given as
and
, respectively. The streamwise,
and wall-normal,
(
) turbulence intensities are the root mean square (RMS) of the fluctuating velocity components in respective directions, representing the strength of turbulence. The wall-normal profiles of the non-dimensional streamwise,
and wall-normal turbulence intensities,
in the wake region of the pipe at a different streamwise distance,
for Run 3 and Run 4 are depicted in
Figure 9 and
Figure 10. It is observed from
Figure 9 and
Figure 10 that the wall-normal distributions of turbulence intensities are akin to RSS distribution with smaller magnitudes. Similar to RSS distributions, at the bed,
and
for both runs, are starting with small positive values, and increasing sharply with an increase in
until they attain their peaks at the top-level of the pipe (
), which indicates the corresponding high-velocity fluctuations due to enhanced turbulence mixing at that level. Turbulence intensities are declining with a further rise in
for
. The peaks of turbulence intensities show increasing trends in the near-wake region and decreasing patterns in the far-wake region along the streamwise direction,
. These patterns are attributed to the enhanced turbulence in the former region and the attenuated turbulence in the latter region. Below the top level of the pipe,
, the turbulence intensities dampened with a fall in
due to a reduction in the turbulence mixing with a decline in
. In the free-surface region,
, turbulence intensities attained very small magnitudes and became approximately constant against
. This result indicates that the free-surface region is far away from the influence of the bed and the cylinder.
Figure 11a,b display the variation of
and
, respectively, against
for Run 1, Run 2, Run 3, and Run 4 in the wall-wake region of the circular pipe, depicting the self-preserving characteristics of the turbulence intensities distributions. The self-preservation of the wall-wake is also evaluated based upon the collapse of the turbulence intensities deficit profiles. The maximum turbulence intensities deficits,
and
, are used as velocity scales, respectively, for
and
, and the half RSS profile width,
, is taken as the wall-normal length scale. The turbulence intensities deficits,
and
, are calculated as the difference between the downstream wall-wake turbulence intensities at a particular elevation,
, and upstream undisturbed turbulence intensities at the same elevation. The turbulence intensities deficits,
and
, at a particular elevation,
, can be written as
and
, respectively. Here,
and
are the upstream undisturbed turbulence intensities in streamwise and wall-normal directions, respectively, at an elevation,
, while
and
are turbulence intensities in streamwise and wall-normal directions, respectively, in the wall-wake region of the pipe at particular
at the same elevation. The turbulence intensities deficit profiles collapse on a narrow band when scaled by these characteristic scales of velocity and length, representing that the wall-wake region of the pipe is self-preserving in nature. The self-preserving behaviour of the turbulence intensities distributions is found to be independent of approach flow conditions. All of the turbulence intensities profiles show good collapse with the streamwise distance,
, for the wall-normal domain of the study,
2.
The streamwise decay of non-dimensional maximum turbulence intensities deficits,
and
for Run 1, Run 2, Run 3, and Run 4 are illustrated in
Figure 11c,d. The approach flow shear velocity,
, is taken as the velocity scale to normalize the maximum turbulence intensities deficits,
and
. It is observed from
Figure 11c that the maximum streamwise turbulence intensity deficit decreases slowly with
. However, the maximum wall-normal turbulence intensity deficit increases in the beginning and attaining maximum value at
and after that, it falls with a further rise in the streamwise distance, as is visible from
Figure 11d.
The third-order correlations give valuable information about the fluxes and diffusion of the Reynolds normal stresses (RNSs) [
25]. To be precise, the effective response of the bursting events can also be determined by the third-order correlations [
24]. The third-order correlations (
) are expressed as
where
= 3,
, and
. The third-order correlations can be categorized into four categories, namely: flux of streamwise RNS in the flow direction,
; flux of the wall-normal RNS in the wall-normal direction,
; diffusion of streamwise RNS in wall-normal direction,
; and diffusion of wall-normal RNS in the streamwise direction,
. The wall-normal profiles of the third-order correlations,
, in the wake region of the pipe at different
for Run 3 and Run 4 are illustrated in
Figure 12.
A closer look at
Figure 12 demonstrates that at the upstream undisturbed location (
),
and
begin with an insignificantly small negative value in the near-bed region and then consistently increased (negative magnitude) with a rise in
for both runs which are typical of turbulent shear flows over rigid surfaces [
22]. Although, in the wake region,
and
start with small positive values,
and
decrease with a rise in the wall-normal distance from the bed and change their signs near the top level of the cylinder for both the runs. With a further rise in
, their negative magnitude increase sharply and attain their negative peaks at
0.4. After that, negative magnitudes of
and
decrease gradually with
and become almost invariant of
in the free surface zone. It is noteworthy that the variation of
and
occur very similarly to
and
, however with the opposite signs. Unlike
and
distributions,
and
begin at bed with small negative values in the wake region. The negative magnitudes of
and
decrease with a further rise in
, and change their signs at
. After that, the value of
and
increase gradually with a rise in the vertical distance from the bed and attain peaks and decrease with a further rise in
, and become almost constant and invariant of
in the near free surface region. The positive and negative values of
and
, respectively, reveal that the downstream flux of streamwise RNS was associated with a downward flux of wall-normal RNS in the near-bed region, while the positive and negative values of
and
indicate the downstream advection of streamwise RNS and downward advection of wall-normal RNS in the near-bed region. Therefore, it can be understood that the dominance of sweep events in the near-bed region makes
and
to be positive and
and
to be negative. The negative values of
and
and positive values of
and
away from the boundary indicates the influence of the ejection events in the disturbed outer layer region. The sweep event represents a high-speed fluid parcel rushed towards the bed and is mainly responsible for bed mobility. From overall observations, it is clear that the enhanced contribution of
(positive) and
(negative) in the downstream region of the pipe reveal an intense momentum exchange in the near-bed region.