3.1.1. Classification of the flow fields

First, the effect of the free surface on the wake flow regime is investigated. Various types of wake interference for two tandem circular cylinders in steady flow have been identified by analysing the overshooting, reattachment and rolling-up behaviours of the shear layers separating from the upstream cylinder (Igarashi Reference Igarashi1981; Zdravkovich Reference Zdravkovich1987; Zhou & Yiu Reference Zhou and Yiu2006). However, with the presence of gravity and a deformable free surface, the wake dynamics can be more intricate (González-Gutierrez, Gimenez & Ferrer Reference González-Gutierrez, Gimenez and Ferrer2019), and is sensitive to both the geometric configuration and the Froude number. Furthermore, the present study is aimed to reveal the underlying physics responsible for the flow influenced by free-surface effects. Therefore, a supplementary classification approach, based on the behaviours of the shear layer originating from the free surface, is employed. The contours of normalized vorticity $\omega _{z}D/U_{\infty }$ in the mid-section of the cylinders are shown in figure 4 at $Fr=0.2$ and figure 5 at $Fr=0.8$, representing low and high Froude numbers, respectively. These cases are selected to illustrate the free-surface effects on the instantaneous flow fields especially at small gap ratios. Through careful observation and classification of the wake flow fields of all cases in the parameter space, eight distinct wake patterns have been identified, as sketched in figure 6, which can be categorized into two basic types: those characterized by the presence of the free-surface shear layer and those without it. The former type is denoted by a prefix ‘F’ to emphasize this distinction. The characteristics of those flow patterns are summarized as follows.

Figure 4.

Instantaneous normalized vorticity $(\omega _{z}D/U_{\infty })$ contours in the mid-section of the cylinders, illustrating the different flow regimes observed in the flow around two circular cylinders in tandem near a free surface at $Fr=0.2$. The green lines indicate the gas–fluid interfaces.

Figure 5.

Instantaneous normalized vorticity $(\omega _{z}D/U_{\infty })$ contours in the mid-section of the cylinders, illustrating the different flow regimes observed in the flow around two circular cylinders in tandem near a free surface at $Fr=0.8$.

Figure 6.

Flow patterns identified in this study for two tandem circular cylinders near a free surface: (a) TW-O; (b) SW-SR; (c) SW-AR; (d) CW-S; (e) F-TW-O; ( f) F-SW-SR; (g) F-SW-AR; (h) F-CW-S. The blue arrows indicate the jet-like flow.

Pattern TW-O: trailing wake with the shear layer overshooting the downstream cylinder, as sketched in figure 6(a). At extremely small gap ratios ($G=0.2$) and low Froude numbers ($Fr \leqslant 0.4$), the lower shear layer separating from the upstream cylinder overshoots the downstream cylinder, while the upper shear layer reattaches onto the downstream cylinder due to the confinement imposed by the presence of the free surface, as shown in figure 4(a,d,g). The stretched lower shear layer resembles a long tail trailing behind the upstream cylinder, hence the name ‘trailing wake’. The primary instabilities of the two shear layers separating from the upstream cylinder are inhibited, possibly due to their significant distance in the cross-flow direction. The presence of the downstream cylinder interrupts the communication between the upper and lower shear layers, effectively stabilizing the wake of the upstream cylinder. Zhao et al. (Reference Zhao, Wang, Zhu, Cao, Bao, Zhou and Han2022) made similar observations in the flow around a single circular cylinder near a free surface at small gap ratios.

Pattern SW-SR: single vortex shedding wake behind the downstream cylinder with steady reattachment of shear layers in the gap region, as sketched in figure 6(b). As the submergence depth increases, the free-surface confinement effect on the upper shear layers weakens. This allows for the vortex shedding to gradually resume, and a single Kármán vortex street can be observed in the wake of the downstream cylinder. The pattern SW-SR actually corresponds to the ‘extended-body’ regime in the classification by Zdravkovich (Reference Zdravkovich1987) and Sumner (Reference Sumner2010). The shear layers separating from the upstream cylinder may continuously reattach onto the downstream cylinder and quasi-stationary vortices are formed in the gap between the cylinders, as shown in figure 4(c).

Pattern SW-AR: single vortex shedding wake behind the downstream cylinder with alternating reattachment of shear layers in the gap region, as sketched in figure 6(c). As the separation distance between the cylinders increases, the shear layers separating from the upstream cylinder may flap up and down, alternately reattaching onto the downstream cylinder, corresponding to the Zdravkovich (Reference Zdravkovich1987) ‘alternate reattachment’ regime. This flow pattern principally involves the reattachment of the shear layers from the upstream cylinder, along with the formation and shedding of eddies in the gap region between the two cylinders, as exemplified in figure 4(b,e,f).

Pattern CW-S: co-shedding wake with vortex shedding in the gap region, as sketched in figure 6(d). At large gap ratios and spacing ratios, a distinct vortex street forms behind each of the two cylinders in this ‘co-shedding’ flow regime (Zhou & Yiu Reference Zhou and Yiu2006). The downstream cylinder is now located outside the vortex formation region of the upstream cylinder and is subjected to the periodic impingement of shed vortices from the upstream cylinder, as illustrated in figure 4(h,i).

Pattern F-TW-O: similar to pattern TW-O, but with a jet-like flow forming in the gap between the tandem cylinders, as sketched in figure 6(e). In this flow pattern, a distinct shear layer separating from the free surface and two shear layers separating from the upstream cylinder can be identified, as shown in figure 5(a,f,k,l). At small gap ratios, the presence of the free-surface shear layer is attributed to the formation of a jet-like flow penetrating through the gap between the cylinders. The occurrence of the jet-like flow has also been observed in the flow around a single circular cylinder near a free surface both in experimental and numerical studies (Sheridan et al. Reference Sheridan, Lin and Rockwell1997; Reichl et al. Reference Reichl, Hourigan and Thompson2005; Zhao et al. Reference Zhao, Wang, Zhu, Ping, Bao, Zhou, Cao and Cui2021).

Pattern F-SW-SR: similar to pattern SW-SR, but with the addition of a free-surface shear layer, as sketched in figure 6( f). In this flow pattern, the reattached shear layer in the gap deflects downward slightly due to the push of the accelerating flow between the free surface and the upstream cylinder, as shown in figure 5(b–d,g–i). The reattachment of the upper shear layer onto the front surface of the downstream cylinder has also been observed by Sumner et al. (Reference Sumner, Price and Paidoussis2000) in the flow around two staggered circular cylinders at very small pitch ratios and small angles of incidence (the shear layer reattachment flow pattern denoted therein as ‘SLR’). The deviation angles of the stagnation point of the upstream cylinder of pattern F-SW-SR exactly fall within the incidence angle range for the flow pattern ‘SLR’ as it turns out in the pressure measurement (see figure 17 in § 3.3.1). The lower shear layer is significantly stretched, deflecting downward and overshooting the downstream cylinder without reattaching to its surface. Large-scale quasi-stationary vortices (Igarashi Reference Igarashi1981), enclosed by the reattached shear layer on the upper side and the free shear layer on the lower side, occupy the gap between the cylinders, effectively blocking the flow through the gap. Under the right conditions, the oblique jet-like flow impinges onto downstream cylinder, as depicted in figure 5(g). After the impingement, the jet-like flow splits into two curved wall jets along the convex surface of the downstream cylinder in both clockwise and anticlockwise directions. The split wall jets are attached to the cylinder wall due to the Coand$\check{\rm a}$ effect (Wille & Fernholz Reference Wille and Fernholz1965). The shear layer separating from the upper surface of the upstream cylinder is compressed by the jet-like flow and captured in the gap between the cylinders. On the other hand, the shear layer separating from the lower side of the upstream cylinder is pushed away without any reattachment to the downstream cylinder. The conclusion regarding the impingement of the jet-like flow is further supported by the visualization of the mean streamlines in figure 13 (see § 3.2.1).

Pattern F-SW-AR: similar to pattern SW-AR, but with the addition of a shear layer separating from the free surface, as sketched in figure 6(g). Examples of flow visualizations for this flow pattern are shown in figure 5(e,j).

Pattern F-CW-S: similar to pattern CW-R, but with the addition of free-surface shear layers, as sketched in figure 6(h). At relatively large gap ratios, the free-surface shear layer may not significantly alter the wake, as shown in figure 5(o). However, at relatively small gap ratios, the near wake of the upstream cylinder deflects downward due to the push of a jet-like flow, allowing the spatial development of shear layers to be no longer constrained by the downstream cylinder. The resulted co-shedding from the two tandem cylinders is similar to that observed in the flow around two staggered cylinders. Examples of this flow pattern are shown in figure 5(m,n). This unique flow pattern will be discussed in detail in the next section. As the gap ratio further increases, the wake transitions from staggered to parallel co-shedding, as shown in figure 5(o).

It should be noted that the absence of the free-surface shear layer in a flow pattern does not imply that the flow is unaffected by the presence of the free surface. Even in the absence of the free-surface shear layer, the flow patterns can still exhibit differences compared with those in a single-phase flow unless the cylinders are sufficiently submerged, as shown in figure 4. In the flow past an isolated cylinder near a free surface, more types of flow patterns have been identified at higher Froude numbers by Bouscasse et al. (Reference Bouscasse, Colagrossi, Marrone and Souto-Iglesias2017) due to violent free-surface distortions.

3.1.2. Staggered flow pattern in two tandem cylinders

With the presence of the free surface, the flow passage above the cylinders is restricted and the oncoming flow is prone to slightly deflect to the bottom side of the upstream cylinder where there is no restriction. Due to the deflection of the oncoming flow, the two tandem cylinders are equivalently exposed to the free stream at small angles of incidence. Therefore, the flow around two tandem circular cylinders of equal diameter may exhibit flow patterns similar to the flow around two staggered circular cylinders due to the effects of the free surface. At small gap ratios, as shown in figure 5(b), the flow fields may show some similarity to the ‘SLR’ flow pattern at very small pitch ratios (Sumner et al. Reference Sumner, Price and Paidoussis2000). At intermediate gap ratios, as shown in figure 5(m,n), the flow field resembles the vortex pairing and enveloping (‘VPE’) flow pattern (Sumner et al. Reference Sumner, Price and Paidoussis2000) and ‘mode S-II’ (Hu & Zhou Reference Hu and Zhou2008; Tong et al. Reference Tong, Cheng and Zhao2015) for two staggered circular cylinders at intermediate incidence angle. To highlight this finding, we present instantaneous flow fields of the case at $(Fr, G, S)=(0.6, 0.6, 3.0)$ at six time instants during a vortex shedding period in figure 7. The specific time instants are labelled with letters (a–f) and marked with red points on the lift force coefficient curve of the upstream cylinder in the bottom inset of figure 7. A movie showing the temporal development of the wake is available as supplementary movie 1 at https://doi.org/10.1017/jfm.2024.1066. The fluid dynamics and key physics of the flow around the two cylinders can be summarized as follows.

Figure 7.

Vorticity contours with streamlines in the mid-section of the cylinders at $(Fr, G, S)=(0.6, 0.6, 3.0)$. The time instants (a–f) are indicated by the red points in the inset of the time history of the lift force coefficient ($C_L$) of the upstream cylinder at the bottom of the figure. The internal captions identify the major vorticity behaviours discussed in the text. Contours vary from blue to red in the range $-3 \leqslant \omega _{z}D/U_{\infty } \leqslant 3$. The dashed lines represent negative vorticity. The green lines indicate the gas–fluid interfaces.

In (a) a hydraulic jump can be observed in the overtopping region of the upstream cylinder. The flow in this region is accelerated as the flow passage narrows due to the presence of the free surface. The accelerated flow forms a jet-like flow, penetrating into the gap region between the cylinders. As will be demonstrated subsequently, the formation of the jet-like flow triggers a sequence of intricate vortex interactions that dominate the wake evolution. The jet-like flow obliquely dives into the gap region along the curved free surface. The deflection of the jet-like flow leaves a slow recirculating region between the jet and the free surface, which leads to the onset of a shear layer of positive vorticity from the free surface. The shear layers separating from the upstream cylinder are pushed away from the free surface and deflect downward by an angle of more than $20$ degrees. The deflected shear layers, previously constrained by the blockage of the downstream cylinder, now have sufficient space to roll up and shed vortices. Meanwhile, an incipient negative vortex begins to develop in the shear layer separating from the upper side of the downstream cylinder.

In (b), with the rolling up process of the shear layer separating from the top side of the upstream cylinder, the jet-like flow approaches closer to the downstream cylinder. As a result, the slow recirculating zone between the jet and the free surface becomes larger, facilitating the continuous development of the positive shear layer from the free surface. The shed vortex with positive vorticity from the bottom side of the upstream cylinder now impinges onto and sweeps over the downstream cylinder. Obviously, the upstream cylinder undergoes an inclined vortex shedding process in the gap between the cylinders. Meanwhile, an incipient negative vortex has formed in the shear layer on the top side of the downstream cylinder. On the other hand, the bottom side of the downstream cylinder lies in the path of the vortex street generated by the upstream cylinder and the shear layer separating from that side is weak.

In (c), as the vortex with negative vorticity is shed from the upper side of the upstream cylinder, the jet-like flow progressively impinges onto the downstream cylinder. The vorticity fed to the positive shear layer separating from the free surface is cut off and that provokes the shedding of a vortex in front of the downstream cylinder.

In (d), the vortex with negative vorticity in the gap between the cylinders catches up the positive vortex below the downstream cylinder after impingement. A relatively small vortex shed from the lower side of the downstream cylinder is trapped between the shed vortices from the upstream cylinder, progressively merging with the large-scale positive vortex. The vortex shedding from the free-surface shear layer is pushed by the jet-like flow and is entrained, passing through the gap between the downstream cylinder and the free surface.

In (e), below the downstream cylinder, the merged positive vortex now pairs with the negative one, translating further downstream along with the free stream. The free-surface vortex decays rapidly when interacting with the shear layer from the upper side of the downstream cylinder. The vorticity of the shear layer separating from the upper side of the downstream cylinder grows stronger as the oncoming flow passage opens again.

In ( f), the image at this instant is essentially similar to figure 7(a), corresponding to the start of the next cycle of wake evolution depicted in figure 7.

It is clear from the cycle of the wake evolution in figure 7 that three vortex shedding processes are present in this flow pattern. The upstream cylinder generates an angled vortex sheet, sweeping over the lower side of the downstream cylinder. This vortex shedding process is highly regular and easily distinguishable in the wake, with alternative vortex shedding from the top and bottom sides of the upstream cylinder. In term of the downstream cylinder, the shear layer separating from the top side experiences a relatively regular roll-up and shedding procedure, while the vortex shedding from the bottom side is rather weak and highly dependent on the local flow, as it lies in the path of the vortex street generated by the upstream cylinder. The presence of the upstream vortex street can disrupt the shedding process and result in a less organized pattern of vortex shedding from the bottom side of the downstream cylinder. Furthermore, the third vortex shedding process is associated with the shear layer separating from the free surface. The free-surface distortion and wave breaking can introduce a substantial quantity of positive vorticity to the vicinity of the free surface. The formation of the jet-like flow above the upstream cylinder facilitates the entry of this surface vorticity into the bulk flow. The dynamic motion of the jet-like flow allows the surface vorticity to be entrained and transported into the gap between the downstream cylinder and the free surface, leading to the shedding process of the shear layer separating from the free surface. In addition, it is the formation of the jet-like flow that alters the incidence angle of the oncoming flow, leading to staggered flow patterns for two tandem cylinders.

3.1.3. Three-dimensional wake

The 3-D wake structures of the cylinders under the free-surface effects are further investigated by showing the isosurfaces of the $Q$-criterion, which represents the second invariant of the velocity gradient tensor. Figure 8 shows the isosurfaces of $QD^2/U_{\infty }^2=0.1$ under the influence of spacing distance between the two tandem cylinders at $(Fr, G)=(0.8, 0.6)$. The blue transparent surface in the figure represents the instantaneous free surface. The 3-D wake pattern at $S=1.5$ is characterized by a pure mode B structure with $6$ streamwise vortex pairs (figure 8a), resulting in a spanwise wavelength of approximately $1.3D$, which is slightly larger than the typical mode B wavelength (around $1D$) for a single circular cylinder (Barkley & Henderson Reference Barkley and Henderson1996; Williamson Reference Williamson1996). According to Williamson (Reference Williamson1996) and Barkley & Henderson (Reference Barkley and Henderson1996), the Reynolds-number range for the transition to mode B instability in a single circular cylinder is $230$ to $260$, which is significantly higher than $Re=180$ in the present study. The early onset of mode B instability may be attributed to the fact that the flow pattern falls within the ‘extended-body’ regime (F-SW-SR pattern in the present study), where the extended characteristic scale of the structure effectively raises the Reynolds number beyond $180$. An arrangement of vortex structures with the same periodicity is also observed in the free-surface shear layer. As seen in figure 8(b–e), when $S>1.5$, the spanwise vorticity in the wake is dominated by the primary vortices with spanwise-periodic deformation at various wavelengths. These 3-D flow structures strongly resemble mode A for a single circular cylinder, although the spanwise wavelengths can be different from the typical value of $4D$ obtained from linear stability analysis by Barkley & Henderson (Reference Barkley and Henderson1996). For example, the spanwise wavelength is approximately $4D$ at $S=2.0$ and $2.5$, $1.6D$ at $S=3.0$ and $2.7D$ at $S=3.5$. In addition, the spanwise periodicity is less regular, which may be attributed to the influence of the free surface.

Figure 8.

Instantaneous flow structures of the tandem cylinders at different spacing ratios with $(Fr, G)=(0.8, 0.6)$, represented by isosurfaces of $QD^2/U_{\infty }^2=0.1$ coloured by normalized streamwise vorticity $\omega _{x}D/U_{\infty }$. Panels show (a) $S=1.5$; (b) $S=2.0$; (c) $S=2.5$; (d) $S=3.0$; (e) $S=3.5$.

Figure 9 shows the isosurfaces of $QD^2/U_{\infty }^2=0.1$ under the influence of submergence depths at $(Fr, L)=(0.8, 4.0)$. At shallow submergence depths, a jet-like flow forms in the gap region between the two cylinders, with the shear layer separating from the free surface dominating the wake, as shown in figure 9(a,b). Hairpin vortex structures are formed in the free-surface shear layer, characterized by an array of large-scale streamwise vortices. Similar hairpin vortex structures, although on a smaller scale, are also identified at much higher Reynolds numbers (see Appendix A). As shown in figure 9(c,d), at intermediate submergence depths, large-scale spanwise vortex structures shedding from the free-surface shear layer can be observed. These free-surface vortices subsequently break down into streamwise structures in the overtopping region of the downstream cylinder. The primary vortex cores exhibit wavy deformation with spanwise wavelengths similar to those of the free-surface vortex structures. When deeply submerged, the wavy deformation of the primary vortex cores disappears, and the flow becomes two-dimensional, as shown in figure 9(e,f).

Figure 9.

Instantaneous flow structures of the tandem cylinders at different gap ratios with $(Fr, S)=(0.8, 4.0)$, represented by isosurfaces of $QD^2/U_{\infty }^2=0.1$ coloured by $\omega _{x}D/U_{\infty }$. Panels show (a) $G=0.2$; (b) $G=0.4$; (c) $G=0.6$; (d) $G=0.8$; (e) $G=1.2$; ( f) $G=2.0$.

Figure 10 shows the isosurfaces of $QD^2/U_{\infty }^2=0.1$ under the influence of different Froude numbers at $(G, L)=(0.6, 4.0)$. The 3-D wake patterns at $Fr=0.2$ and $0.4$ are mode A with spanwise wavelengths of approximately $4D$, as shown in figure 10(a,b). However, as the Froude number increases to $0.6$ and $0.8$, a large quantity of irregular vortex structures appear near the free surface, leading to more pronounced spanwise inhomogeneities in the wake of the cylinders. This indicates that the three-dimensionality of the flow around two tandem cylinders near a free surface becomes more pronounced with increasing Froude number.

Figure 10.

Instantaneous flow structures of the tandem cylinders at different Froude numbers with $(G, S)=(0.6, 4.0)$, represented by isosurfaces of $QD^2/U_{\infty }^2=0.1$ coloured by $\omega _{x}D/U_{\infty }$. Panels show (a) $Fr=0.2$; (b) $Fr=0.4$; (c) $Fr=0.6$; (d) $Fr=0.8$.

3.1.4. Summarized dependence of flow patterns on $G$ and $S$

The dependence of flow patterns on spacing ratio and gap ratio is summarized in figure 11. This map may contribute to the prediction of the flow structure in the wake of two tandem circular cylinders near a free surface. At different Froude numbers, the distribution of wake patterns varies on the $G$–$S$ plane, as shown by the lines in figure 12.

Figure 11.

Dependence of flow patterns on $G$ and $S$ in the wake of two tandem circular cylinders. To differentiate the flow patterns, the following symbols have been used: $\star$, TW-O; $\diamond$, SW-SR; $\circ$, SW-AR; $\triangle$, CW-S; $\triangledown$, F-TW-O; $\triangleleft$, F-SW-SR; $\square$, F-SW-AR; $\triangleright$, F-CW-S. Four flow patterns at the four Froude numbers are plotted at each coordinate grid point, with the light green symbol representing $Fr=0.2$, the light red symbol representing $Fr=0.4$, the black symbol representing $Fr=0.6$ and the light blue symbol representing $Fr=0.8$. Boundaries between various flow behaviours at $Fr=0.6$ are superimposed for reference. The internal captions of different colours correspond to the boundaries of the same colour.

Figure 12.

Variation of the flow pattern distribution as a function of the Froude number on the $G$–$S$ plane. (a) The influence of the free-surface shear layer, the dashed lines indicate the upper limit where the free-surface shear layer is present in the wake, the dash-dotted lines of the same colour denote the region where a jet-like flow forms. (b) Boundaries of the reattachment regime and the co-shedding regime.

On the one hand, as shown in figure 12(a), the horizontal borders indicate the varying degrees of influence of the free surface on the wake flow patterns, with a predominant dependence on the $y$-direction variable $G$. As discussed previously, a shear layer separating from the free surface is introduced into the wake due to the distortion of the free surface and the blocking effect of the cylinders at small and intermediate gap ratios. It is clear that higher Froude numbers can extend the effect of the free surface to deeper submergence depths, as indicated by the dashed lines in figure 12(a). At the minimum Froude number ($Fr=0.2$) considered in this study, no free-surface shear layer forms. However, the free-surface shear layer occurs at $G \leqslant 0.8$ at $Fr=0.4$, extending significantly to $G \leqslant 1.6$ at $Fr=0.8$. It can be seen that it is easier for the free-surface shear layer to form large spacing ratios. In particular, the free-surface shear layer occurs only when $S \geqslant 3.00$ at $Fr=0.4$, as depicted by the region enclosed by the red dashed lines figure 12(a). The strength and the dynamic behaviours of the free-surface shear layer play a decisive role in the flow state in the gap between the two tandem cylinders. As can be observed in figure 5, the formation of the jet-like flow in the gap between the cylinder is consistently accompanied by a strong free-surface shear layer. The jet-like flow typically forms at high Froude numbers ($Fr \geqslant 0.6$), with higher $Fr$ corresponding to a higher upper limit of $G$, as indicated by the dash-dotted lines denoting the ‘jet-like flow’ border in figure 12(a). The jet-like flow exhibits two forms: one penetrates the gap region between the two cylinders, and the other impinges on the downstream cylinder, as indicated by the blue arrows in figures 6(e) and 6( f), respectively. Regarding the penetrating jet-like flow, the roll-up of the shear layers separating from the upstream cylinder is often interrupted, resulting in the cessation of the regular vortex shedding process, as shown in figure 5( f,k). On the other hand, the impinging jet-like flow pushes downward the shear layer separating from the upper side of the upstream cylinder, causing it to reattach on the front surface of the downstream cylinder, as shown in figure 5(b,g). The angles of the impingement of the jet-like flow varies with $Fr$, which will be discussed later in § 3.3.1. At large gap ratios, the formation of the shear layer separating from the free surface is no longer observed, which means that the influence of the free surface on the flow pattern is significantly diminished.

On the other hand, the vertical borders illustrate different behaviours of shear layers separating from the upstream cylinder, with a higher dependence on the variable $S$ in the $x$ direction, as shown in figure 12(b). It can be observed that the critical spacing ratios at which the reattachment regime transitions to the co-shedding regime are different at varying submergence depths and Froude numbers. At a low Froude number of $Fr=0.2$, the critical gap ratio remains at $3.25$ when the cylinders are located close to the free surface ($G \leqslant 0.6$), but decreases to $3.00$ when $G \geqslant 0.8$. A similar trend is observed for $Fr=0.4$, except that the critical spacing ratio increases to $3.25$ when $G \geqslant 1.6$. As the Froude number increases to $0.6$, the critical spacing ratio reaches a significant low value of $2.75$ at $G=0.6$, and then remains at $3.25$ over a large range of gap ratios ($0.8 \leqslant G \leqslant 1.8$). The advance of the critical spacing ratio is attributed to the formation of an inclined jet-like flow in the gap between the cylinders, which suppresses vortex shedding from the upstream cylinder. However, this advance does not occur for $Fr=0.8$, since the greater deflection angle of the jet-like flow at higher Froude numbers allows the formation of staggered vortex shedding, as discussed previously in § 3.1.2. At $Fr=0.8$, a slightly larger critical spacing ratio of $3.5$ is observed at intermediate submergence depths ($1.0 \leqslant G \leqslant 1.2$). The critical spacing ratio recovers to $3.5$ as the cylinders are deeply submerged ($G \geqslant 2.0$). In general, the critical spacing ratio tends to be larger at higher Froude numbers, although the jet-like flow with small deflection angle formed at shallow submergence depths may cause the critical value to advance.

The different distribution of the flow patterns at varying $Fr$ is mainly due to the different amplitudes of free-surface distortions and the local flow state in the vicinity of the free surface (Reichl et al. Reference Reichl, Hourigan and Thompson2005; Zhao et al. Reference Zhao, Wang, Zhu, Cao, Bao, Zhou and Han2022). At a low Froude number ($Fr=0.2$), only a negligible amount of vorticity resides at the free surface due to the small deformation of the free surface. As a result, the flow patterns without the free-surface shear layer are primarily observed at low Froude numbers, as indicated with light blue symbols in figure 11. However, at higher Froude numbers ($Fr \geqslant 0.4$), in addition to the substantial quantity of vorticity produced due to free-surface sharpening and wave breaking, the curved profile of the free surface effectively directs this free-surface vorticity into the bulk flow, significantly interfering with the wake evolution.

3.2.1. Mean flow topology

It was established in the preceding section that the instantaneous flow features are significantly influenced by the presence of the free surface. We further present the free-surface effects on the mean flow topology in this section. Figure 13 shows several representative mean flow fields, revealing a significant asymmetry in the mean wake of the two circular cylinders due to the presence of the free surface. The mean flow topology can be correlated with the transient flow fields discussed in the previous section, offering a deeper insight into the underlying physical mechanisms involved.

Figure 13.

Time-averaged streamlines with mean streamwise velocity ($u/U_\infty$) in the wake of two tandem circular cylinders at $Fr=0.6$. The green lines indicate the time-averaged gas–fluid interfaces.

For the averaged wake of the upstream cylinder, the recirculation bubble in the near wake is narrowed and deflects toward the bottom side. In cases with small spacing ratios (figure 13a–c), the wake of the upstream cylinder is tightly confined within the narrow gap between the two tandem cylinders and pushed downward by the accelerated flow above the cylinders. There is a wall jet attached onto the front surface of the downstream cylinder, penetrating through the narrow gap between the cylinders. This is in contrast to the stagnant or quasi-steady fluid observed for two tandem circular cylinders without the free surface at small spacing ratios (Igarashi Reference Igarashi1981). The difference is due to the loss of symmetry in the system configuration with the presence of the free surface and the jet impingement, as discussed previously. This penetrating gap flow was also observed by Sumner et al. (Reference Sumner, Price and Paidoussis2000) for two staggered circular cylinders at small angles of incidence ($10^\circ$ to $30^\circ$). At large spacing ratios, the wake of the upstream cylinder is driven downward and biased away from the centreline by the inclined jet-like flow formed above the upstream cylinder, as observed in figure 13(e,f). A comparison of figure 13(e,f) reveals that the jet-like flow has a smaller inclination angle at larger gap ratio, resulting in less confinement to the near wake of the upstream cylinder. The reduced confinement and the angled wake enable the vortex shedding process to occur for the upstream cylinder, leading to the pattern F-CW-S discussed previously. As the gap ratio and spacing ratio increase, the wake of the upstream cylinder progressively restores its symmetry, as exemplified in figure 13(h).

In terms of the averaged wake of the downstream cylinder, at small spacing ratios, the flow above the cylinders is effectively blocked and the entrainment requirement of the wake is not sufficiently met. Consequently, a subsequent low-velocity recirculation zone forms behind the cylinders, as shown in figure 13(a). In particular, when the separation spacing between the two cylinders is large enough for the jet-like flow to penetrate through the gap between the cylinders without impinging onto the downstream cylinder, the downstream cylinder will be completely immersed within the large-scale recirculation region formed behind the jet-like flow near the free surface, as illustrated in figure 13(e). It is evident that the length of the recirculation region decreases with deeper submergence depths. The extremely short formation region behind the downstream cylinder is attributed to the impingement of the upstream wake, as suggested by Carmo et al. (Reference Carmo, Meneghini and Sherwin2010b).

3.2.2. Mean mass flux ratio

As discussed above, the flow state within the gap significantly influences the wake evolution. To avoid ambiguity, the term ‘gap’ in this subsection specifically refers to the gap between the upper edge of the cylinders and the free surface. Here, two mean mass flux ratios ($R_1$ and $R_2$) are defined to quantify the gap flow state, as illustrated in figure 14. Jiang et al. (Reference Jiang, Cheng, Draper and An2017) defined a similar quantity to characterize the flow within the gap between a circular cylinder and a moving wall. The mean mass flux ratio $R_1$ is defined as the ratio of the time-averaged mass flux through the gap at the upstream cylinder to the time-averaged mass flux of the free-stream flow, which is measured between the free surface and the cylinder centreline ($\kern0.7pt y=0$) at $10D$ in front of the upstream cylinder. The mean mass flux ratio $R_2$ is defined as the ratio of the time-averaged mass flux through the gap at the downstream cylinder to the time-averaged mass flux through the gap at the upstream cylinder. The definitions are given as

(3.1a,b)$$\begin{gather} R_1=\frac{F_1}{F_0},\quad R_2=\frac{F_2}{F_1}, \end{gather}$$

(3.2)$$\begin{gather}F=\frac{1}{T} \int_{t}^{t+T} \iint_A \rho \boldsymbol{U} \boldsymbol{\cdot} \mathrm{d} \boldsymbol{A} \,\mathrm{d} t, \end{gather}$$

where $F$ represents the time-averaged mass flux through the planar surface $A$ over the time duration $T$, while $F_0$, $F_1$ and $F_2$ denote the mean mass flux at the respective locations as marked in figure 14.

Figure 14.

Sketch of the definition of the mean mass flux ratios.

Figure 15 shows the variation of the mean mass flux ratios as a function of spacing ratio and gap ratio at different Froude numbers. It can be observed that the variation of $R_1$ with $G$ exhibits a consistent trend across different spacing ratios and Froude numbers, which suggests that the gap flow at the upstream cylinder is minimally influenced by the downstream cylinder. As the submergence depth increases, $R_1$ gradually increases from a value less than $0.4$ to nearly $1.0$. With the presence of the free surface, the mean mass flux ratios at the upstream cylinder stay below $1.0$ due to the fact that the flow passage above the cylinder is restricted by the free surface and it is easier for the flow to deflect to the bottom side of the cylinders where there is less restriction (Zhao et al. Reference Zhao, Wang, Zhu, Cao, Bao, Zhou and Han2022). Moreover, as the gap ratio increases, the restriction imposed by the free surface becomes progressively weaker. When the cylinders are placed extremely far below the free surface, it is expected that the mean mass flux ratio $R_1$ would approach a value of $1.0$ due to the restoration of symmetry in the system.

Figure 15.

Variation of the mean mass flux ratio with $G$ and $S$ at (a) $Fr=0.2$, (b) $Fr=0.4$, (c) $Fr=0.6$ and (d) $Fr=0.8$. The solid and dashed lines denote $R_1$ and $R_2$, respectively.

In contrast, the gap flow at the downstream cylinder is significantly affected by the upstream cylinder and the distortion of the free surface. As shown in figure 15(a), at a low Froude number ($Fr=0.2$), the mean mass flux ratio $R_2$ keeps above approximately $0.5$, indicating that a substantial portion of the fluid entering the gap at upstream cylinder is transferred to the gap at the downstream cylinder. As the Froude number increases, the formation of the jet-like flow at small gap ratios redistributes the flow from the gap at the upstream cylinder. At small spacing ratios, the jet-like flow splits into two parts along the surface of the downstream cylinder in both clockwise and anticlockwise directions due to the impingement onto the downstream cylinder. The anticlockwise one diverts part of the flow, resulting in a decrease in the mean mass flux ratio $R_2$, which can be inferred from the flow visualizations in figures 5(a) and 13(a). When the jet-like flow penetrates through the region between the cylinders without impinging onto the downstream cylinder, a large-scale recirculation region forms behind the jet-like flow near the free surface, as visualized in figures 5( f) and 13(e). The downstream cylinder is completely immersed in this low-velocity recirculation region. Consequently, the mean mass flux ratio $R_2$ has an extremely small value, approaching zero, as shown in figure 15(d). As the gap ratio increases, the jet-like flow disappears and the mean mass flux ratio $R_2$ gradually approaches $1.0$. The representative topologies of the mean flow field corresponding to the variations of the mean mass flux ratios are shown in the insets in figure 15(d).

3.3.1. Pressure distribution

The free-surface effects on the flow fields discussed so far are, of course, also reflected in the pressure coefficient, $2\bar {p}_{rgh}/\rho U_{\infty }^{2}$. The variations of the mean pressure coefficients around the two cylinders with the circumferential angle $\theta$ at $Fr=0.6$ are shown in Cartesian and polar coordinates in figure 16, where the time-averaged pressure at the middle cross-section ($z=4D$) is illustrated. Note that the contribution of hydrostatic pressure to the pressure coefficient has been excluded. Here, $\theta$ is the angle on the circumference taken from the rear base point on each cylinder, as illustrated in figure 16(b).

Figure 16.

Mean pressure coefficient distributions around two tandem cylinders at various submergence depths at $Fr=0.6$ for $S=1.50$ (a,b), $S=2.50$ (c,d) and $S=3.50$ (e,f). In Cartesian plots, the dashed and solid lines represent the pressure distribution around the upstream and downstream cylinders, respectively. The grey shaded region indicates the part of the cylinder surface facing the free surface. In polar plots, the direction of the arrows indicates the sign of the pressure, with positive values pointing inward towards the cylinder and negative values pointing outward away from the cylinder.

For the upstream cylinder, it is obvious that the mean pressure distribution around the cylinder surface is asymmetrical with respect to the cylinder centreline. The stagnation point deviates from the leading edge ($\theta =180^\circ$) towards the upper side facing the free surface. The angle of deviation increases rapidly as the cylinder approaches the free surface, as depicted by the dashed arrow in figure 16(a). The maximum of deviation angle is approximately $22.9^\circ$ at $G=0.20)$ for all the gap ratios, and the minimum of deviation angle is approximately $3.5^\circ$ at $G=2.00)$, as shown in figure 16(a). A comparison between figures 16(a), 16(c) and 16(e) reveals that the separation distance between the two tandem cylinders hardly changes the deviation angle, which indicates that the downstream cylinder has a limited influence on the flow in front of the upstream cylinder. The pressure on the upper windward side is lower at small gap ratios than at large gap ratios, as shown in figure 16(b,d,f). This is due to the narrowing of the flow passage above the cylinders when they are placed close to the free surface. As a result, the flow is accelerated in that region, leading to a subsequent reduction in pressure. The decrease in pressure on the base side with increasing gap ratio indicates the enhancement of the vortex shedding from the upstream cylinder (Tong et al. Reference Tong, Cheng and Zhao2015).

In contrast, the mean pressure distribution around the downstream cylinder is significantly influenced by the presence of the upstream cylinder and the free surface. At small gap ratios, the inclined jet-like flow would impinge onto the downstream cylinder when it is placed close to the upstream cylinder, resulting in a region with positive pressure, as shown in figure 16(b,d). The location of the maximum pressure coefficient, corresponding to the jet impingement point, shifts upward as the gap ratio increases, as depicted by the solid arrow in figure 16(a). Additionally, the magnitude of the maximum pressure attenuates as the gap ratio increases due to the weakening of the jet-like flow. In the absence of the jet-like flow or the jet impingement, the pressure coefficient around the downstream cylinder is consistently negative, as indicated by the outward arrows in figure 16(b–f). Furthermore, figure 16(d,f) shows that, for the downstream cylinder, the red lines at small gap ratios are nearly enclosed within the blue lines at large gap ratios, indicating that the pressure magnitude decreases as the cylinder approaches the free surface. The decline of pressure magnitude implies a more energetic flow around the downstream cylinder due to the agitation of the distorted free surface. It can be observed that the mean pressure distribution around the downstream cylinder is characterized by one, two and one evident local maximum on the windward side at small, intermediate and large gap ratios, as shown in figure 16(a,c,e). This corresponds to the transition in the behaviour of the shear layers in the gap between the two tandem cylinders from impingement to reattachment to shedding.

The influence of the Froude number, spacing ratio and gap ratio on the mean stagnation pressure ($C_{ps}$) and shift angles ($\alpha$) of mean stagnation points is presented in figure 17. The shift angle of the mean stagnation point is defined as the difference between the leading edge of the cylinder ($\theta =180^\circ$) and the angle of stagnation point, as illustrated in figure 17(b). As mentioned earlier, under certain conditions, the mean pressure distribution on the surface of the downstream cylinder exhibits two local maxima. In this case, only the larger maximum is considered as the stagnation point on the downstream cylinder. In terms of the upstream cylinder, the maximum shift angle always occurs at the smallest gap ratio, with a value of around $20^\circ$. The shift angle then decreases towards zero as the gap ratio increases, due to the progressive restoration of symmetry in the wake of the upstream cylinder, as discussed previously in § 3.2.1. The Froude number has a weak effect on the shift of the stagnation points on the upstream cylinder, although the shift angle increases slightly as the Froude number increases. For the downstream cylinder, when the spacing ratio is below the critical value, the shift angle rises from approximately $10^\circ$ to over $55^\circ$ as the gap ratio increases from $0.2$ to $2.0$. Simultaneously, the stagnation pressure decreases from positive values to negative values, as shown in figure 17(a–d). The deviation in shift angles between low and high Froude numbers reaches approximately $40^\circ$ at $G=0.2$, while the difference weakens as the gap ratio increases. Additionally, the mean stagnation pressure is significantly higher at larger $Fr$. This notable difference is primarily attributed to the strength and deflection angle of the jet-like flow formed near the free surface. However, when the flow enters the ‘co-shedding’ regime, the shift angle remains relatively stable, fluctuating within a small range of $20^\circ$ to $35^\circ$, and the mean stagnation pressure stays below $0.4$, as shown in figure 17(e,f).

Figure 17.

Variations of mean stagnation pressure and shift angles of mean stagnation points as a function of gap ratios at (a) $S=1.5$, (b) $S=2.0$, (c) $S=2.5$, (d) $S=3.0$, (e) $S=3.5$ and ( f) $S=4.0$. The size of the markers represents the magnitude of the stagnation pressure, with solid markers indicating positive values and hollow markers indicating negative values.

3.3.2. Force coefficients

Figure 18 shows time histories of the lift and drag coefficients and their fast Fourier transform (FFT) spectra for the upstream and downstream cylinders at $(Fr, S)=(0.8, 3.5)$ for three submergence depths. At $G=0.2$, the flow pattern F-TW-O is observed, as shown in figure 5(k). In this flow pattern, a jet-like flow penetrates through the gap between the cylinders and vortex shedding is completely suppressed. Consequently, the r.m.s. values of hydrodynamic coefficients are very small, with $C_D^\prime =0.016$ and $C_L^\prime =0.024$ on the upstream cylinder, and $C_D^\prime =0.041$ and $C_L^\prime =0.050$ on the downstream cylinders, as shown in figure 18(a,d). Only a minor peak at half the vortex shedding frequency appears in the spectrum, as shown in figure 18(g). As the gap ratio increases to $0.6$, the flow pattern changes to F-CW-S, where staggered vortex shedding occurs, as illustrated in figure 5(m). Figure 18(b,e) shows a significant increase in the r.m.s. values of all forces, particularly for the downstream cylinder, where $C_D^\prime$ increases to $0.270$ and $C_L^\prime$ to $0.266$, while the mean drag coefficient rises from $-0.057$ to $0.596$. Obviously, this discontinuity in hydrodynamic forces corresponds to the change in the flow pattern from F-TW-O to F-CW-S. In the spectrum (figure 18a), alongside the prominent peak at the vortex shedding frequency, small peaks are also observed at half and twice this frequency due to the influence of the free-surface shear layer. The presence of the second subharmonic and harmonic of the vortex shedding frequency results in irregular time histories of the hydrodynamic forces, rather than a sinusoidal form, as shown in figure 18(b,e). As the gap ratio further increases to $1.0$, the wake transitions from staggered to parallel co-shedding, as illustrated in figure 5(o). The r.m.s. values of lift coefficients exhibit significant increase, as shown in figure 18(c), with $C_L^\prime$ rising to $0.420$ on the upstream cylinder and to $0.875$ on the downstream cylinder. The second subharmonic of the vortex shedding frequency is no longer present in the spectrum due to the weakening of the free-surface shear layer.

Figure 18.

Time histories of lift and drag coefficients and their FFT for the upstream and downstream cylinders at $(Fr, S)=(0.8, 3.5)$. Panels show (a,d,g) $G=0.2$, pattern F-TW-O; (b,e,h) $G=0.6$, pattern F-CW-S (staggered co-shedding); (c,f,i) $G=1.0$, pattern F-CW-S (parallel co-shedding).

Figure 19 summarizes the dependence of mean hydrodynamic forces on the spacing ratio and gap ratio for the two tandem cylinders near a free surface. Note that the lift force here excludes the contribution of hydrostatic force (buoyancy). The contours of mean force coefficients also reflect the dependence of the flow patterns on $G$ and $S$ by comparing with figure 11.

Figure 19.

Dependence on $G$ and $S$ of mean drag coefficients (a,c,e,g) and mean lift coefficients (b,d,f,h) on the two tandem cylinders at $Fr=0.2$ (a,b), $Fr=0.4$ (c,d), $Fr=0.6$ (e,f) and $Fr=0.8$ (g,h). The white lines indicate the critical spacing ratios. In each panel, the left contours represent values for the upstream cylinder and the right contours represent values for the downstream cylinder. Contours vary from blue to white to red, representing values ranging from negative to zero to positive. The solid and dashed lines represent positive and negative values, respectively.

In terms of the mean drag force, the ‘drag inversion’ (Carmo et al. Reference Carmo, Meneghini and Sherwin2010b) is well reproduced at any Froude number. As shown in figure 19(a,c,e,g), the mean drag coefficient on the upstream cylinder jumps from a relatively small value to a higher value approaching that of a single isolated cylinder with increasing spacing ratio. Similarly, the mean drag coefficient on the downstream cylinder exhibits a discontinuous jump from negative values (a thrust force Zdravkovich Reference Zdravkovich1977) to positive values as the spacing ratio increases. The discontinuity in the values of the hydrodynamic forces when the separation distance between the two tandem cylinders varies continuously has been reported extensively in the literature (e.g. Arie et al. Reference Arie, Kiya, Moriya and Mori1983; Ljungkrona et al. Reference Ljungkrona, Norberg and Sundén1991; Alam et al. Reference Alam, Moriya, Takai and Sakamoto2003).

Although the ‘drag inversion’ occurs for any submergence depth, the mechanism responsible for small gap ratios (roughly $G \leqslant 0.4$) is fundamentally different from that for larger gap ratios. The latter is attributed to the transition of the flow from the reattachment regime to the co-shedding regime, which is the result of the wake interference between the two cylinders. However, at small gap ratios, the reduction of the mean drag force on the upstream cylinder results from the incidence angle of the oncoming flow induced by the presence of the free surface. The angled oncoming flow in front of the upstream cylinder produces a downward inclined force on the upstream cylinder. Only the streamwise component of this force contributes to the drag force, resulting in a smaller drag compared with that in a uniform free stream, while the transverse component contributes the lift force on the upstream cylinder, leading to a negative value of the lift force, as shown in figure 19(b,d,f,h).

Regarding the mean lift force, as shown in figure 19(b,d,f,h), it is not surprising to see that the mean lift coefficients are not equal to zero due to the asymmetry of the mean wake topology discussed in § 3.2.1. It can be seen that the majority of the lift coefficient values across the $G$–$S$ plane are negative. This can be attributed to the presence of the jet-like flow or the accelerated flow in the narrow gap between the cylinders and the free surface. The mean lift coefficients on the downstream cylinder recover to zero at large spacing ratios and gap ratios, but not for the upstream cylinder. This means influence of the free surface on the upstream cylinder is more severe than on the downstream cylinder.

Figure 20 is a compilation of the Strouhal numbers for the two tandem cylinders across the range of gap ratios and Froude numbers examined in this study. It is evident that the vortex shedding frequency is consistent for both the upstream and downstream cylinders, except in certain cases at small gap ratios. This frequency inconsistency is attributed to the no-shedding flow pattern, specifically patterns TW-O and F-TW-O, where the downstream cylinder is immersed in the low-velocity region formed behind the upstream cylinder, as shown in figures 4(a) and 5( f). Due to the low velocity of the oncoming flow, the Strouhal numbers are significantly higher for the downstream cylinder compared with the upstream cylinder. As shown in figure 20(a,b), the Strouhal numbers increase with deeper submergence depths for all Froude numbers at small spacing ratios ($S=1.5$, 2.0). At large spacing ratios ($S=3.5$, 4.0), the Strouhal numbers exhibit different trends depending on the Froude number: they rise from lower values (approximately $0.15$) with increasing $G$ at low Froude numbers ($Fr = 0.2$ and $0.4$), while they decrease from higher values (around $0.225$) with increasing $G$ at high Froude numbers ($Fr = 0.6$ and $0.8$), as shown in figure 20(e,f). Eventually, at large gap ratios, the Strouhal numbers converge to the same value ($\approx$0.175). The small Strouhal numbers observed at small gap ratios and low Froude numbers are attributed to the delayed roll-up process of the upper shear layer caused by the confinement of the free surface, as shown in figure 4(h). In contrast, the high Strouhal numbers at small gap ratios and high Froude numbers are due to the presence of the free surface that causes the oncoming flow to deviate from the centreline of the cylinder, thus reducing the effective oncoming flow velocity, as shown in figure 5(m).

Figure 20.

Variations of Strouhal numbers for the two tandem cylinders as a function of gap ratios at (a) $S=1.5$, (b) $S=2.0$, (c) $S=2.5$, (d) $S=3.0$, (e) $S=3.5$ and ( f) $S=4.0$.

3.4.1. Mean surface elevation

Figure 22 shows the profiles of the time-averaged elevations of the free surface along the streamwise direction for three Froude numbers addressed in the present study. Yu & Tryggvason (Reference Yu and Tryggvason1990) has suggested that the Froude number is the decisive parameter governing the free-surface distortion. Although not shown here, at the lowest Froude number ($Fr=0.2$), the free surface exhibits negligible distortions at any gap ratio, with the maximum mean elevation not exceeding $2\,\%$ of the cylinder diameter. At a similar Froude number ($Fr\approx 0.24$), Miyata, Shikazono & Kanai (Reference Miyata, Shikazono and Kanai1990) observed no discernible free-surface distortions in the flow around a circular cylinder steadily advancing beneath a free surface. This can be attributed to the low effective Froude number, where the small inertial force compared with gravity prevents significant distortions of the free surface (Sarpkaya Reference Sarpkaya1996). However, as the Froude number increases, the inertial force becomes increasingly dominant compared with gravity, resulting in noticeable deformations of the free surface due to the stirring and mixing effects caused by the significant inertial force, as shown in figure 22 (note the difference in the limits of the $y$-coordinate). The maximum mean elevation reaches approximately $0.05D$, $0.15D$ and $0.25D$ at $Fr=0.4$, $0.6$ and $0.8$, respectively. Noticeable surface distortions have also been observed in the flow around a single circular cylinder close to a free surface at large Froude numbers (e.g. Reichl et al. Reference Reichl, Hourigan and Thompson2005; Moreira & Peregrine Reference Moreira and Peregrine2010; Zhao et al. Reference Zhao, Wang, Zhu, Ping, Bao, Zhou, Cao and Cui2021).

Figure 22.

Variation of profile of mean free-surface elevation ($\bar {\eta }/D$) with $G$ and $S$ at $Fr=0.4$ (a–c), $Fr=0.6$ (d–f) and $Fr=0.8$ (g–i). The dashed lines represent vertical centrelines of the upstream and downstream cylinders.

As shown in figure 22, a pronounced hydraulic jump induced by the obstruction beneath the free surface occurs immediately above the upstream cylinder. This leads to a rather abrupt rise with a constant height in the free surface followed by a widely extended trough. In general, a sink with varying depth appears immediately behind each cylinder. In the region behind the two tandem cylinders, the mean free-surface elevation remains close to the still water level, while significant elevation fluctuations can be identified due to wave generation, as discussed in the next subsection.

The occurrence of the hydraulic jump is attributed to the abrupt change in flow regime from locally supercritical to locally subcritical, as suggested by Chanson (Reference Chanson2004). This phenomenon and its underlying mechanism are virtually the same as those observed in the flow around a single circular cylinder near a free surface, as investigated in previous studies (Zhao et al. Reference Zhao, Wang, Zhu, Ping, Bao, Zhou, Cao and Cui2021, Reference Zhao, Wang, Zhu, Cao, Bao, Zhou and Han2022). It can be observed from figure 22 that the variation of the height of the hydraulic jump almost remains consistent across varying spacing distance between the cylinders. Figure 23 shows variations of the hydraulic jump height ($H_{Jump}$) as a function of gap ratios for all Froude numbers at $S=3.5$. The height of the hydraulic jump shows a clear tendency to increase with higher Froude number and decrease with deeper submergence depth. The jump height reaches its maximum at the smallest gap ratio ($G=0.2$): $0.017D$ at $Fr=0.2$, $0.072D$ at $Fr=0.4$, $0.160D$ at $Fr=0.6$ and $0.274D$ at $Fr=0.8$. It then decreases to $0.1D$ at $Fr=0.8$ and falls below $0.05D$ at $Fr \leqslant 0.6$ at the largest gap ratio $G=2.0$ examined in the present study.

Figure 23.

Variations of the hydraulic jump height as a function of gap ratios at $S=3.5$.

3.4.2. Free-surface waves

The gravity waves induced by a submerged cylinder in a uniform flow near a free surface has attracted extensive attention. Assuming a steady irrotational potential flow with a linearized free-surface boundary condition and replacing the cylinder by an equivalent doublet at its centre, Lamb (Reference Lamb1932) observed a local distortion immediately above the obstacle followed by a train of stationary sinusoidal gravity waves downstream of the cylinder. Furthermore, Tuck (Reference Tuck1965) considered the nonlinear effects from the free-surface boundary and suggested that the nonlinear solution would involve highly non-sinusoidal or even breaking waves. Other authors who have more recently investigated free-surface gravity flow around a submerged cylinder include Dagan (Reference Dagan1971), Moreira & Peregrine (Reference Moreira and Peregrine2010) and Semenov & Wu (Reference Semenov and Wu2020). However, the present vortex system involves intricate arrays of vortices shed from both the cylinders and the free surface, and these vortices are not fixed in position but convect and diffuse within the viscous flow. Therefore, the distortions of the free surface are a consequence of the mutual induction and interaction of all vorticity concentrations (Cetiner & Rockwell Reference Cetiner and Rockwell2001).

Figure 24 presents the temporal evolution of the free surface at a spacing ratio of $S=4.00$ for the three Froude numbers examined in the present study. The selected gap ratios of $G=0.40$, 0.80, 1.60 represent shallow, intermediate and deep submergence, respectively. Gravity waves are discernible on the downstream side, although the regularity and continuity of the waves vary depending on the specific parameter combinations. Free-surface undulations can be observed between the two tandem cylinders, which results from the large drop of the hydraulic fall. This phenomenon is similar to the trapped waves induced by two successive submerged obstacles investigated by many other authors, such as Dias & Vanden-Broeck (Reference Dias and Vanden-Broeck2004), Page & Părău (Reference Page and Părău2014) and Shelton & Trinh (Reference Shelton and Trinh2023). At the free surface between the cylinders, the waves are not as well organized as the waves far downstream, and the wavelengths are shorter and more irregular. Note also that some ripples may be formed in front of the upstream cylinder at small Froude numbers ($Fr \leqslant 0.4,$), while a stationary hydraulic jump with a constant height is formed at a large Froude number ($Fr \geqslant 0.6$). The reason may be that the relatively small upstream Froude numbers might potentially intersect with the linear dispersion relation (Wang et al. Reference Wang, Chai, Părău, Page and Wang2022), leading to a train of waves upstream. It can be observed that the induced waves downstream of the cylinders have higher phase velocities and shorter wavelengths at higher Froude numbers, as indicated by the steeper slope of the wave propagation trajectories. A summarized profile of the free-surface distortion is illustrated in figure 25.

Figure 24.

Temporal evolution of the free surface at $S=4.00$ for $Fr=0.2$ (a–c), $Fr=0.4$ (d–f), $Fr=0.6$ (g–i) and $Fr=0.8$ (j–l). The green dashed lines indicate the centres of the upstream and downstream cylinders.

Figure 25.

Schematic profile of the free-surface distortion in the flow around two tandem circular cylinders near a free surface.

Figure 26 shows variations of wave height ($H_{Wave}$) as a function of gap ratios for all Froude numbers at $S=3.5$. Here, wave height is defined as the r.m.s. value of free-surface elevation multiplied by $2\sqrt {2}$. Wave heights are measured at two streamwise locations: $x/D=-10$ (upstream of the cylinders) and $x/D=10$ (downstream of the cylinders). As shown in figure 26(a), observable wave heights occur only at $Fr=0.4$, while the wave heights stay well below $0.03D$ for the other Froude numbers. This observation is consistent with figure 24, which indicates that wave generation upstream of the cylinders only occurs at $Fr=0.4$. A comparison between figures 26(a) and 26(b) reveals that the wave heights upstream at $Fr=0.4$ are consistently higher than the wave heights downstream at any submergence depth. As shown in figure 26(b), the wave height downstream of the cylinder increases and then decreases with increasing submergence depth. The height of wave is related to the strength of the vortex structures and the regularity of the vortex streets beneath the free surface. At shallow submergence depths, the vortex shedding from the cylinders can be disrupted or completely suppressed by the free-surface shear layer. As a result, waves generated behind the cylinders are less energetic, with low wave heights at small gap ratios. As the submergence depth increases, the interference effects of the free-surface shear layer are reduced, allowing for a more regular and continuous vortex street, so that the fluid can be stirred up violently, thus generating waves with large wave heights. As the gap ratio further increases, the vortex streets are too far from the free surface to effectively transfer energy to the free surface, leading to a decrease in wave height. As shown in figure 26(b), the wave height reaches its maximum value of $0.10D$ at $G=0.8$ for $Fr=0.4$, $0.15D$ at $G=1.0$ for $Fr=0.6$ and $0.22D$ at $G=1.2$ for $Fr=0.8$. It can also be observed that the wave height increases with higher Froude number, and the gap ratio at which the wave height reaches it maximum also increases with higher Froude number. This is consistent with the conclusion from other authors (Reichl et al. Reference Reichl, Hourigan and Thompson2005; Bouscasse et al. Reference Bouscasse, Colagrossi, Marrone and Souto-Iglesias2017; Zhao et al. Reference Zhao, Wang, Zhu, Cao, Bao, Zhou and Han2022) that a larger Froude number leads to more drastic free-surface deformation.

Figure 26.

Variations of wave height as a function of gap ratios for $S=3.5$ at (a) $x/D=-10$ and at (b) $x/D=10$.

Figure 27 shows variations of wave height as a function of spacing ratios for all Froude numbers at $G=0.8$. It can be observed that the wave heights stay below approximately $0.05D$ before the critical spacing ratio, while they increase abruptly at the critical spacing ratio. This is attributed to the wake transition from the ‘reattachment’ regime to the ‘co-shedding’ regime, where two vortex streets behind the two tandem cylinders in the ‘co-shedding’ regime can induce waves with large wave heights.

Figure 27.

Variations of wave height as a function of spacing ratios for $G=0.8$ at $x=10$.

3.4.3. Characteristics of the induced waves

Furthermore, figure 28 shows the dispersion relations of the waves on the downstream side at $S=4.00$, $G=1.60$ for three Froude numbers examined in the present study, obtained by performing the space–time fast Fourier transform of surface elevations. The abscissa represents the wavenumber ($k$) normalized by the diameter of the cylinder ($D$). The frequency is normalized in the form of Strouhal frequency. The linear dispersion relation of waves on a mean current in deep water is given by (Phillips Reference Phillips1966)

(3.3)\begin{equation} (\omega-kU)^2=gk, \end{equation}

where $g$ is the acceleration due to gravity, $k$ is the wavenumber, $\omega$ is the angular frequency and $U$ is the constant velocity of the mean current. On the downstream side of the cylinders, the flow velocity gradually recovers to the that of the free stream, i.e. $U=U_\infty$ in the present study. The dispersion relations of waves downstream of the cylinders show good agreement with the linear dispersion relation at harmonics of the vortex shedding frequency.

Figure 28.

The dispersion relations of the waves on the downstream side at $S=4.0$, $G=1.6$ for: (a) $Fr=0.2$, (b) $Fr=0.4$, (c) $Fr=0.6$ and (d) $Fr=0.8$. The white curves represent the linear dispersion relations of waves on a mean current in deep water.

It can be observed that the frequencies of the waves cluster around the vortex shedding frequency, as shown in figure 28. Figure 29 shows the time histories of the wave heights and their power spectral density (PSD) along the streamwise locations at $x/D=-10$, $-5$, 0, 5, 10 for $(Fr, G, S) = (0.4, 0.8, 1.6)$, corresponding to the case shown in figure 24( f). The peak frequencies of the waves at all streamwise locations are consistent and align with the vortex shedding frequency, which further confirms the correlation between the vortex shedding and wave generation on the downstream side.

Figure 29.

Wave height at different streamwise locations at $(Fr, G, S) = (0.4, 0.8, 1.6)$. (a) Time histories of the wave height at $x/D=-10$, $-5$, 0, 5, 10; (b) PSD of the wave heights.

It should be noted that well-organized waves on the downstream side are only allowed to occur when the vortex street is stable and the vorticity is well concentrated in the wake of the cylinders. In cases where vorticity is dispersed due to unstable vortex shedding, vortex impingement or interference from the free-surface vorticity, the generation of free-surface waves is inhibited, as shown in figure 24(j). Compared with the wave generation by a single fixed submerged point vortex (Shelton & Trinh Reference Shelton and Trinh2023) or a single fixed submerged doublet (Moreira & Peregrine Reference Moreira and Peregrine2010), the generation of waves by moving vortex arrays in the present system can be more complex and challenging.