Abstract

Numerical simulations are carried out to study the flow around two tandem square cylinders (SC) under the effect of spacing ratio(g/D) and splitter plate length (l/D) for a fixed Reynolds number (Re) = 100. The g/D is varied from 0 to 10 and l/D is varied from 0.5 to 10. The splitter plate length is found to have strong effect on vortex shedding and fluid forces. The maximum reduction in mean drag coefficient is observed at l/D = 8, that is 15% and 78% for upstream and downstream cylinders, respectively. The maximum reduction in root-mean-square value of lift coefficient is found at l/D = 10, that is 99%. The flow pattern at both of these points is steady flow. There is 100% vortex shedding suppression for l/D > 5. The observed flow patterns for flow past tandem cylinders without splitter plate are; single bluff body (SBB), steady flow (SF), quasi-steady flow (QSF), fully developed flow (FDF) and fully developed two-row vortex street flow (FDTRVS) regimes. SBB, QSF and SF regimes were observed in presence of splitter plate.

1. Introduction

The study of fluid forces and vorticity dynamics through the bluff body has been greatly improved in the last few decades. Most of these studies have been taken for single cylinder. However, much less attention has been paid to multiple cylinders flow compared to single cylinder. The addition of another bluff body can change the different aspects of flow such as fluid forces and transition threshold. To study the flow interference different models such as side-by-side, tandem and staggered models have been used. Among these models, tandem is the simplest one. A number of experimental and numerical studies have been conducted for flow past tandem cylinders. Kondo and Matsukuma [1] numerically analyzed the salient flow features through two tandem circular cylinders at Re = 1000. They found that the forces exerted on the last cylinder are substantially higher than those exerted on the first cylinder. While in case of tandem square cylinders, the upstream cylinder forces are larger than those of the downstream cylinder [2]. Azuma et al. [3] experimentally measured forces acting on square structures both side-by-side and in tandem. They observed that the drag forces of two tandem cylinders changed significantly when the separation ratio was greater than g/D = 3 (where is the spacing between cylinders and D is cylinder diameter) as compared to side-by-side arrangement. Abbasi et al. [4] numerically analyzed the flow regime transitions over the flow around three tandem cylinders for the g/D between 0.5 and 16 at Re = 200 and found five different types of flow regimes. In addition, they noted that the g/D is the important analysis parameter. Abbasi et al. [5] analyzed the flow regime transitions around two-, three- and four-tandem cylinders for Re ranging from 1 to 130 at g/D = 2 and 5. They observed two separate trends in terms of flow transitions for g/D = 2 and 5. The more reduction in terms of drag is observed at g/D = 2. Mittal et al. [6] numerically analyzed the flow characteristics through two cylinders for Re = 100 and 1000 and found that the shedding starts at Re = 100. Abbasi et al. [7] numerically found three various types of flow features for two tandem cylinders. Vikram et al. [8] also analyzed the salient flow features past tandem cylinders. They observed that compared to downstream cylinder the upstream cylinder experience higher lift force. Sharman et al. [9] numerically analyzed the flow through two inline cylinders under the influence of g/D varying from 2 to 10 for Re = 100. They identified g/D = 3.75–4 as critical spacing, where the flow characteristics changes abruptly. Liu and Jerry [10] experimentally analyzed flow on two cylinders in tandem with Re = 2000–16000 and g/D = 1.5–9. They reported hysteresis with two discontinuous jumps associated with different flow regimes. Sohankar and Etminan [11] analyzed flow and heat transfer across cylinders arranged inline for Re = 1–200. They found that vortex shedding formation begins in the range 35 ≤ Re ≤ 40.

The addition of an eddy promoter improved the heat transfer of the downstream cylinder and controls the CD. Bearman [12, 13] used the splitter plates to increase the formation length and reduce the drag. One of the most important studies was carried out by Zdravkovich [14], which divided the vortex suppression into three different branches, namely: (i) surface protrusions, (ii) shrouds and (iii) near-wake stabilizers. In literature both active and passive methods exist. Active control procedures generally require continuous external energy to reduce the drag forces and improve the lift [15, 16]. While for the passive methods there is no need internal energy, these include tripping rods [17], control cylinder [1820], splitter plate [21], T-shaped plate [22] etc.

Anderson and Szewczyk [23] experimentally analyzed the impacts of low length splitter plates for the flow through circular cylinder and noticed that less than D/8 splitter length appreciably affect the wake zone for Re = 2700 to 46,000. You et al. [24] observed that significant noise reduction can take place for flow around circular cylinder once the splitter plate length is identical to the cylinder diameter. Uffinger et al. [25] analyzed the effects of different shapes on flow control behind a cylinder and found that the wedge shape behind the square body significantly control the flow. Ali et al. [26] computationally analyzed the influence of splitter plate length on a cylinder at Re = 150.

To the best of ours’ knowledge, for the proposed problem there is no numerical and experimental study reported in the literature. Hopefully, this work will help engineer’s working on fluid forces reduction. The main motivation is to examine the splitter plate effect on flow regimes. The other main objective is to stabilize the critical gap spacing. Furthermore, we want to systematically examine the change over behavior of the mechanism of wake structure behind the downstream cylinder and amplitude variation of fluid forces.

This paper is ordered as follows: The proposed problem, numerical method, initial and boundary conditions are discussed in Section 2. The effect of computational domain and code validation is done in section 3. The computed results are presented in section 4. Finally, the conclusion is presented in section 5.

2. Problem Description and Numerical Details

2.1. Problem Description

A schematic representation of the flow through two fixed tandem SCs in presence of attached splitter plate is shown in Figure 1. D, l and are the width of the cylinder, length of the splitter plate and spacing between the cylinders, respectively. C1 and C2 are the first and second cylinders, respectively. Lu, Ld, Lx, and Ly are the upstream location, downstream location, length, and height of the domain, respectively. Parabolic velocity profile is adopted at the entrance position of the domain, while the convective boundary condition is adopted at the exit of the domain [27]. No-slip boundary condition is applied to the surface of the cylinders and splitter plate [28]. Free-slip boundary condition is adopted at the top and bottom lateral boundaries of the channel [29]. The method of momentum exchange [30] is used for the calculation of forces.

2.2. Lattice Boltzmann Method

We have established the 2-D code for the flow through two tandem SCs with and without attached splitter at  = (0–10, 0.5–10) for Re = 100. The lattice Boltzmann method (LBM) is applied for the numerical results obtained in this paper. As compared to well known numerical techniques, the LBM is relatively new technique [28]. It was developed by Frisch et al. [31].

Here we will present a short overview of this method. In this study we applied Q2N9 (where Q is the space dimension and N is the number of particles) lattice model. In this model, each computational node consists of a rest particle (i = 0) along with eight moving particles (i = 1–8) (see Figure 2).

The density evolution equation is given by;where, and are the particle distribution function, the corresponding equilibrium distribution function, dimensionless time, position of particles, velocity directions, and single relaxation time, respectively.

The corresponding equilibrium distribution function is:

Here, , and are the fluid density, velocity and corresponding weighting functions ( for i = 0,  = 1/9 for i = 1–4, and  = 1/36 for i = 5–8), respectively. Equation (1) can be solved by two steps: collision which use a Bhatnagar-Gross-Krook (BGK) operator [32, 33] and propagation. These steps can be defined as:

Collision:

Streaming:

Equations (5) and (6) are used to calculate the density and velocity at each computational node

The boundary conditions are applied after the streaming step (equation (4)), and the entire process is solved iteratively. The pressure can be calculated from the equation of state: ; where cs is the speed of sound and its value is .

3. Computational Domain Study and Code Validation

In this section we will present the grid independence and code validation.

3.1. Grid Independence Study

Table 1 presents the grid independence study. We have taken the error with respect to D = 30. From the results it is observed that as the size of D is increased from 20 to 30, the maximum variation in mean drag coefficient (CDmean) and Strouhal number (, where is the vortex shedding frequency) are 0.336% and 0.9%, respectively. Therefore, for the rest of simulation we used D = 20.

3.2. Computational Domain

The computational domain size strongly affects characteristics of flow around bluff bodies. The small change in computational domain changes the vortex shedding procedure, which affects the physical parameters. We have calculated the CDmean and St for different values of Lu, Ld and Ly, so that we can choose a suitable domain. The data in Table 2 shows that there is a slight difference in between the resulted values. So we can consider any case given in Table 2. We will simulate the given problem by using Lu = 8D; Ld = 25D and Ly = 10D. The numbers 1 and 2 in subscripts represents the physical parameters of upstream and downstream cylinders.

3.3. Code Validation

The validation of code is done for the flow past single cylinder for Re = 100. For the validation of the code we have used only two parameters that is CDmean and St. For comparison the numerical results of Dutta et al. [34] as well as experimental data of Norberg et al. [35] and Okajima [36] are given in Table 3. The results shows that the present calculations of CDmean and Stare in good agreement with the experimental data of Norberg et al. [35] and Okajima [36]and numerical data of Dutta et al. [34].

4. Results and Discussion

A two-dimensional (2D) numerical study is carried out to study the effect of /D and l/D at Re = 100. The splitter plate is attached in the middle of the two cylinders. In this study we will discuss the vorticity contours visualization, time analysis of CD and CL and power spectra of CL. During this study the /D will be varied from 0 to 10 and l/D will be varied from 0.5 to 10.

4.1. The Flow Over Two Tandem Square Cylinders without Splitter Plate

Figures 3(a)–3(f) presents the instantaneous flow visualization of flow field through two SCs in tandem arrangement without splitter plate. The different /D values in Figures 3(a)–3(f) represented the flow filed characteristics is SBB (/D = 0.5 and 1), SF regime (/D = 2 and 3), FDF regime (/D = 5) and FDTRVS regime (/D = 9), respectively. The pattern of vortex was a single row in SBB flow regime (Figure 3(a)). The phenomenon was analogous to the flow past a single SC. When the value of /D is increased up to 1, the shear layers detached from the C1 reattach on the upper surface of the C2 (Figure 3(b)). When the /D is increased to 2 and 3, no vortex shedding occurs behind the second cylinder (Figures 3(c), and 3(d)). When /D = 5, the fully developed vortex shedding can be clearly seen behind the upstream square cylinder. Under these conditions, the flow behind the downstream SC was almost matching to that around single cylinder, and there was a little affect by the upstream SC (Figure 3(e)). On the other hand, the downstream SC was still affected by the upstream SC, appearing to be a two-rows vortex street at /D = 9 (Figure 3(f)).

In the FDF regime, the two inline cylinders separately shed vortices at the same frequency. Igarashi [37] also reported that the two inline circular cylinders shed vortices at the same frequency at reasonably large spacing ratio. It is observed that the FDF or co shedding flow regime occurs at g/D = 5 where the vortices shed separately from the two cylinders but with same frequency. The negative vortex shedding from C1 hits the front face of the C2 and at the same time the negative and positive shed vortices develop from the cylinder C2. These different shed vortices interact with each other and results in a FDTRVS regime with multi-frequency variation in the spectra of the lift coefficient of C2 (ref. spectrum in Figure 4(j)).

Instantaneous streamline visualization shows the alternate generation of negative and positive shed vortices presented in Figures 5(a)–5(f). It is found that the wake width and vortex formation length are considerably dependent on spacing ratio. In SBB regime, the shear layers separated from the front edges of C1 quickly reattached on the lateral surface of C2, forming two small recirculation regions within the cylinders (Figure 5(b)). Due to strong suction a negative drag value is observed for C2. At g/D = 5, both the cylinders (C1 and C2) shed vortices and the flow characteristics is called a FDF regime.

Time variation CD and CL of both cylinders are shown in Figures 6(a)–6(l) for different spacing ratios. As spacing ratio increases, the amplitude of drag coefficient also increases for FDF regime and FDTRVS regime (Figures 6(i), and 6(k)). It is observed that the CD of C2 is modulated. At g/D = 0.5 and 1, the CL become sinusoidal, and the first cylinder and second cylinder shed vortices result in an inphase mode (Figures 6(b), and 6(d)). At /D = 2, the CL of both cylinders becomes constant (Figure 6(f)).

The spectra graph shows a single dominant peak corresponding to the primary frequency; no secondary frequency was observed in the spectra in case of single bluff body flow regime (Figures 4(a)–4(d)). The small peak in the power spectra (Figures 4(h) and 4(j)) suggests that the second cylinder flow is still affected by the first cylinder.

4.2. The Flow Around Two Tandem Square Cylinders with Attached Splitter Plate

Figures 7(a)–7(e) presents the vorticity contours visualization of flow field around through inline SCs. The different splitter plate length within the spacing between the two cylinders represented the flow filed characteristics in reattachment flow regime (/D = l/D = 0.5 and 1) and steady flow regime (/D = l/D = 2.5, 5 and 10). In shear layers reattachment flow (SLR) regime, the shear layers separated from the C1 reattach on the upper surface of the C2 (Figures 7(a) and 7(b)). When /D = l/D ≥ 2.5, the steady flow regime can be clearly seen behind both cylinders (Figures 7(c)–7(e)). Under these circumstances, the flow behind the upstream and downstream SCs was almost similar, and it was hardly affected by the downstream SC. Inspection of Figures 8(a)–8(e) further confirms that the flow through the downstream SC becomes steady only in existence of the splitter plate of l/D ≥ 2.5.

The higher upward mean lift force is observed for the downstream SC when the two cylinders are closely spaced as the viscous effect is more towards the lower face of the downstream square cylinder. From Figures 9(a), 9(c), 9(e) and 9(g) it can be seen that at l/D = 0.5, 1, 2.5 and 10 the CD1 and CD2 are both constants. Furthermore, the CD2 having negative drag coefficient at l/D = 0.5, 1 and 2.5. On the other hand, the lift force has periodic behaviour with decreasing amplitude as the value of l/D increased. Furthermore, at l/D = 10, the lift force is a straight line due to steady flow behaviour behind the cylinders.

The power spectra of CL of the two cylinders with different splitter plate length are shown in Figures 10(a)–10(f). At l/D = 0.5, the lift force oscillation of the upstream cylinder was stable, which ensures that the interference is weak between the cylinders. When the splitter plate length is varied from 0.5 to 1, the power spectra not changed. With the further increase in splitter plate length, the interference among the two tandem cylinders slowly disappeared, and the oscillation became stable again. For the downstream cylinder, its lift force oscillation had almost the same behaviour as the upstream square cylinder. The difference was observed when the splitter plate length was varying from 5 to 10. Some minor peaks in the spectra of the upstream square cylinder were observed. This is due to the vortex generated by the C1 directly acts on the C2.

The graphical representation of the different observed flow regimes is given in Figure 11. It can be observed that at small value of g/D and splitter plate length the single bluff body flow regime occurs. At intermediate values of /D and l/D the shear layers reattachment flow regime occurs. The FDF regime and FDTRVS regime can be seen only for the two tandem square cylinders without splitter plate. It can also be observed from the flow diagram that the splitter plate considerably suppressed the vortex shedding behind the cylinders and completely convert it to steady flow regime.

4.3. Analysis of Force Statistics

In Figure 12(a), compared with the CDmean of a single SC, the CDmean of C1 was slightly reduced up to /D = 4.5 without splitter plate. As the g/D increased, the CDmean of the C1 remained unchanged and was almost equal to an isolated cylinder. The CDmean of the C2 was considerably smaller than that of single cylinder due to the influence by the wake of the upstream SC. The value of CDmean of the C2 was negative when g/D is varied from 0.5 to 4. When the value of g/D increased to 4, the CDmean of C2 changed from negative to positive. As the g/D continued to increase (g/D > 4), the CDmean of the downstream cylinder first increased and then slowly starts to decrease. The transition from shear layers reattachment flow regime to two-row vortex street flow regime was actually a positive-negative conversion of the CD2. It is found that in presence with the splitter plate, the value of CDmean of the C2 is considerably smaller than in the cases without splitter plate. The splitter plate length has marginal role on drag reduction for /D > 5.

CDrms1 and CDrms2 are the root-mean-square drag values of the first cylinder and second cylinder, respectively (Figure 12(b)). The discontinuity is observed in CDrms2 without splitter plate in Figure 12(b) for the case of /D = 4.5 due to transition of flow regime from shear layers reattachment flow regime (/D = 4) to fully developed flow regime (/D = 4.5). For a splitter plate length of l/D 1 to 10, the CDrms of both cylinders is nearly equal to that of the isolated cylinder. However, it is found without splitter plate that the CDrms2 is irregular.

From Figure 13(a), it is seen that the CL of the cylinder C2 is more affected as compared to the cylinder C1 due to the change in /D. The dependency of CLrms2 is more pronounced at higher spacing ratios. On the increase of length of the splitter plate, the CLrms2 apparently lower than the value of an isolated SC. The shedding frequency is one of the fluidic parameters which can affect the C2and alter wake dynamics. Similar to the numerical results of Sharman et al. [30], the St values of the C1 and C2 are the same with all g/D values considered in this study.

Figure 14 presents the percentage reduction of CDmean of flow past two tandem cylinders. Regarding the CDmean2, the control effectiveness increases with increasing the splitter plate length: i.e., 49, 53.4, 57.2, 56.1, 62, 63.1, 67.1, 60, 53.1, 59.4, 61, 65.2, 79.1, 76.1, and 74.2% for l/D = 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, 6, 7, 8, 9, and 10, respectively. The larger splitter plate length will lead to a better control of fluid forces.

5. Conclusions

In this article, the numerical simulation for flow past two tandem cylinders in presence of attached splitter plate was carried out at Re = 100. This study investigates the drag reduction and the interference among the cylinders at various splitter plate length of l/D = 0.5 to 10. The observed flow regimes were divided into five different regimes: single-bluff body, shear layers reattachment, steady, fully developed two-row vortex shedding and fully developed flow regimes. The CDmean2 was negative when g/D is varied from 0.5 to 4. We not found fully developed shed vortices between the cylinders in case of single-bluff body and shear layers reattachment flow regimes. The discontinuity is observed in CDrms2 without splitter plate for the case of /D = 4.5 due to transition of flow regime from shear layers reattachment flow regime (/D = 4) to fully developed flow regime (/D = 4.5). The dependency of root-mean-square value of lift coefficient of the downstream cylinder (CLrms2) is more pronounced at higher spacing ratios. Regarding the mean drag coefficient of the downstream square cylinder, the maximum reduction about 79.1% is observed for l/D = 8.

Data Availability

Data will be provided on demand.

Conflicts of Interest

The authors declare that they have no conflicts of interest.