# Application of Potential Theory to Steady Flow Past Two Cylinders in Tandem Arrangement

Authors: Yangyang Gao; Danielle S. Tan; Zhiyong Hao; Xikun Wang; Soon Keat Tan
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

The wake flow patterns associated with flow past a cylinder and a cylinder-pair in tandem configuration are revisited, compared, and evaluated with respect to the streamline patterns generated based on potential flow theory and superposition of various potential flow elements. The wakes, which are vortex shedding in the lee of the cylinder(s), are reproduced by placing pairs of equal but opposite circulation elements in the potential flow field. The strength of the circulation elements determines the size of the vortices produced. The streamline patterns of flow past a pair of unequal cylinders in tandem configuration provide an indirect means to establish the threshold condition for the wake transition from that of a single bluff body to alternating reattachment behavior. This threshold condition is found to be a function of the diameter ratio,d
/
D (diameters d and D, d
≤
D ), spacing ratio, L
/
D (centre-to-centre distance, L, to cylinder diameter, D), and equivalent incident flow speed, U. A unique functional relationship f (L
/
D, d
/
D, U) of this threshold condition is established.


---

## Body

## 1. Introduction

Flow past a cylinder and two circular cylinders in tandem configuration are common occurrences in many engineering applications, such as mooring lines, riser, and bundled cylinders at offshore installations, bridge piers, and tube bundles in heat exchangers. The wake flow patterns in the lee of a single cylinder are well understood, while there is still much unknown in the case of a cylinder-pair in tandem configuration. A cylinder-pair may be of different diameters, with equal diameters being the special case. Many attempts have been made to categorize the flow behavior behind a pair of tandem cylinders with equal diameter (downstream of the first cylinder and in the lee of the cylinder-pair) (Igarashi [1, 2]; Ljungkrona and Sundén [3]; Sumner et al. [4]; Lin et al. [5]; Xu and Zhou [6]; Zhou and Yiu [7]; Carmo et al. [8, 9]; Sumner [10]). Zdravkovich [11] identified three basic types of interference based on the range of flow Reynolds number (Re): (i) single bluff body (composite) behavior at small centre-to-centre spacing (type 1); (ii) shear layer reattachment at intermediate spacing (type 2); and (iii) Karman vortex shedding from each cylinder independently at larger spacing (type 3). There is but a handful of studies on unequal cylinder-pairs (Igarashi [12]; Mahbub Alam and Zhou [13]; Zhao et al.  [14, 15]; Gao et al. [16, 17]). Under certain cross-flow speed (U or Re) and depending on the spacing ratio L
/
D and diameter ratio d
/
D, these three types of flow interference can also be observed, as illustrated in Figure 1. Here, for a fixed incident flow speed U
(
=
1
) and d
/
D
=
2
/
3, the reattachment and alternating vortex shedding patterns in the wake of the upstream cylinder begin to appear at L
/
D
=
2.4, with the cylinder-pair behaving like a single bluff body for L
/
D
<
2.4. This observation indicates that the value L
/
D
=
2.4 corresponds to a threshold condition that demarcates the end of single bluff body behavior and the beginning of alternating vortices in the region between the two cylinders. For L
/
D
=
2.4 and d
/
D
=
2
/
3, all three types of interference can be observed consecutively with increasing incident flow speed (Reynolds number), as shown in Figure 2, which indicates that the Reynolds number has significant effect on determination of the threshold condition.Instantaneous flow patterns around two tandem cylinders withd
/
D
=
2
/
3 at different L
/
D, for a fixed incident flow speed (Re). (a) Laboratory observation based on particle image velocimetry (PIV) and (b) schematic flow structure.
(a)
(b)Instantaneous flow patterns around two tandem cylinders atL
/
D
=
2.4 for a fixed d
/
D but different Re.
(a)
(b)
(c)Of the three types of interference, most engineering applications would prefer the bundled cylinders to behave as a single bluff body. Therefore, it is desirable to establish the threshold condition below which the single bluff body behavior is observed. The authors are not aware of any analytical expression for this threshold value ofL
/
D, or said expression, as a function of d
/
D and speed of the cross-flow U, or Re. Reported spacing ratios L
/
D at which type 1 and type 2 interference take place are listed in Table 1. It can be seen from this table that the upper limit of the spacing ratio for type 1 interference is sensitive to the Reynolds number and is consistent with Igarashi’s [1] suggestion.Table 1
Summary of spacing ratio ranges for single bluff body and reattachment flow patterns behind two tandem circular cylinders.
Author
Cylinder diameter:d, D
Re
(
L
/
D
)
I
(
L
/
D
)
II
Igarashi [1]
d
=
D
8700
–
5.2
×
10
4
1.03–1.91
2.09–3.53
Zdravkovich [11]
d
=
D
Not given
1–1.2~1.8
1.2~1.8–3.4~3.8
Ljungkrona and Sundén [3]
d
=
D
(
1.0
–
1.2
)
×
10
4
1.25
2.0, 4.0
Xu and Zhou [6]
d
=
D
800
–
4.2
×
10
4
1.0–2.0
2.0–3.0
Zhou and Yiu [7]
d
=
D
7000
1.3
2.5, 4.0
Gao et al. [16]
d
<
D
1200–4800
1.2
2.4
(
L
/
D
)
I: centre-to-centre spacing ratio range for single bluff body pattern.
(
L
/
D
)
II: centre-to-centre spacing ratio range for reattachment flow pattern.The authors noticed that there are uncanny similarities between the streamline patterns generated using potential theory for a pair of doublets placed at different centre-to-centre spacing ratios and the flow patterns shown in Figures1 and 2. It appears that one could draw an analogy between potential flow about a pair of doublet and flow past a cylinder-pair in tandem arrangement. Using potential flow theory as a tool to aid in the derivation of analytical/numerical solutions for high Reynolds number flows is not new. This technique has been applied in many studies on surface waves, the wake around a ship hull, and ship-to-ship interactions, where the fluid is assumed to be steady, incompressible, inviscid, and irrotational in the far field (Millward et al. [18]; Shahjada Tarafder and Suzuki [19]; Ahmed and Guedes Soares [20]). Potential flow around multiple cylinders has also been investigated by many researchers to derive analytical solutions for potential flow past multiple cylinders, as in Johnson and McDonald [21], Burton et al. [22], and Crowdy [23]. Considering the analogous patterns of potential flow streamlines and the flow patterns behind cylinder-pair in tandems, the obvious question arises as to whether one could deduce from these patterns a precise threshold condition analytically and account for the onset of type 2 and type 3 interferences.For this purpose, the authors drew on the potential flow theory approach to generate streamline patterns for the following cases:(1)
wake flow in the lee of a single cylinder;
(2)
flow patterns in the lee of a cylinder-pair, as a function of spacing ratio and diameter ratio, with the aim of determining the threshold spacing ratio for the onset of streamline patterns similar to that of type 2 and type 3 interference.


## 2. Derivation of Streamline Patterns Based on Potential Flow Theory

### 2.1. Flow Past a Single Cylinder

When an ideal fluid flows past a single circular cylinder, the streamlines and velocity potentials form the same pattern as that of a doublet placed in a constant uniform flow as suggested by Acheson [24]. Thus, the complex potential function for a uniform flow past a single cylinder may be expressed as
(1)
W
(
z
)
=
V
(
z
+
1
z
)
,
where V is the flow velocity parallel to the x-axis and z is the complex variable. For flow past a cylinder, the uniform flow field is represented by the flow velocity V. The streamline patterns obtained from the potential theory based on 1
/
V are found to be appropriate. Hence, the uniform flow field is represented by U
=
1
/
V in this study, and the potential function is expressed as ψ
(
U
)
=
ψ
(
1
/
V
).The complex flow past a circular cylinder over thez-plane can be constructed by superposition of potential flow elements. In this case, the authors use a uniform flow and vortex element with certain circulation strength. Here, the vortex element in the potential flow is represented by a rotating body of a certain diameter D and circulation strength Γ. The circulation is either clockwise or counterclockwise depending on the sign assigned to Γ. Consequently, the complex potentials for flow past a single circular cylinder can be visualized as being the resultant of a uniform flow and a rotating vortex element and can be written as
(2)
W
(
z
)
=
U
z
+
U
a
2
z
+
i
Γ
2
π
ln
⁡
z
.Consider the wake flow patterns captured using the particle image velocimetry (PIV) technique and shown in Figures3(a)–3(c); the gyros in the flow are vortices of certain strength in the flow field. These gyros could be represented by placing a circulation element of appropriate strength in the flow, with the corresponding streamlines presented in Figures 3(d)–3(f). The potential equation corresponding to these streamlines is as follows:
(3)
W
(
z
)
=
F
1
(
z
)
+
F
2
(
z
)
+
F
3
(
z
)
+
F
4
(
z
)
+
F
5
(
z
)
=
U
z
+
U
a
2
[
β
1
z
-
L
1
+
β
2
z
-
L
2
+
β
3
z
-
L
3
+
β
4
z
-
L
4
+
β
5
z
-
L
5
]
+
i
Γ
2
π
[
ln
⁡
(
z
-
L
1
)
α
1
+
ln
⁡
(
z
-
L
2
)
α
2
+
ln
⁡
(
z
-
L
3
)
α
3
+
ln
⁡
(
z
-
L
4
)
α
4
+
ln
⁡
(
z
-
L
5
)
α
5
]
,
where F
j
(
z
) represents the complex potential for each circulation element j (j
=
1,2
,
3,4
,
5). β
j and α
j represent the diameter ratio and circulation strength ratio, respectively, for each element. It can be seen from Figure 3 that the principal features of the wake flow patterns could be reproduced by using circulation elements to represent the gyros.Figure 3
Comparison of instantaneous flow patterns around a single cylinder at different Re obtained using PIV (a)–(c), and streamline flow patterns based on potential flow (d)–(f).While the patterns generated reflect physically meaningful flow phenomenon, it is recognized that the above approach requires prior knowledge of the location and strength of the circulations.


### 2.2. Flow Past a Cylinder-Pair in Tandem Configuration

For flow past two tandem cylinders as shown in Figure4, the complex flow can be constructed by superposition of a uniform flow, two vortex elements, and two pertinent circulation strengths with opposite signs. Consider the diameter ratio of the upstream to downstream cylinders is β
=
a
′
/
a (a
′ and a are the radii of the upstream and downstream cylinders, resp.) and circulation at the upstream cylinder is Γ
′
=
β
Γ (where Γ is the circulation of the downstream cylinder). The complex potentials for the upstream and downstream cylinders located, respectively, at z
=
-
l and z
=
l may be written as
(4)
F
1
(
z
)
=
U
z
+
U
(
β
a
)
2
(
z
+
l
)
-
i
β
Γ
2
π
ln
⁡
(
z
+
l
)
,
F
2
(
z
)
=
U
z
+
U
a
2
(
z
-
l
)
+
i
Γ
2
π
ln
⁡
(
z
-
l
)
.Figure 4
A schematic diagram of two cylinders with unequal diameters in tandem arrangement. “O” is the origin. The size of the cylinder is reflected by the strength of circulation. Here, D
=
2
a and d
=
2
a
′, for a
>
a
′.Then, the resultant complex potentialW
(
z
) for two tandem cylinders with unequal diameters in a steady flow (see Figure 4) is
(5)
W
(
z
)
=
F
1
(
z
)
+
F
2
(
z
)
=
U
z
+
U
a
2
(
z
-
l
)
+
U
(
β
a
)
2
(
z
+
l
)
+
i
Γ
2
π
ln
⁡
⁡
(
z
-
l
)
-
i
β
Γ
2
π
ln
⁡
⁡
(
z
+
l
)
=
U
[
z
+
a
2
(
z
-
l
)
+
(
β
a
)
2
(
z
+
l
)
]
+
i
Γ
2
π
ln
⁡
[
(
z
-
l
)
-
β
ln
⁡
(
z
+
l
)
]
.Puttingz
=
x
+
i
y into (5), the complex potential may be expressed as
(6)
W
(
z
)
=
F
1
(
z
)
+
F
2
(
z
)
=
U
x
+
U
a
2
[
(
x
-
l
)
(
x
-
l
)
2
+
y
2
+
β
2
(
x
+
l
)
(
x
+
l
)
2
+
y
2
]
-
Γ
2
π
{
arg
[
(
x
-
l
)
+
y
i
]
-
arg
[
(
x
+
l
)
+
y
i
]
β
}
+
i
{
(
x
-
l
)
2
[
(
x
-
l
)
2
+
y
2
]
β
U
y
-
U
a
2
[
y
(
x
-
l
)
2
+
y
2
+
β
2
y
(
x
+
l
)
2
+
y
2
]
+
Γ
4
π
ln
⁡
(
x
-
l
)
2
+
y
2
[
(
x
-
l
)
2
+
y
2
]
β
}
.
The velocity function ψ and velocity potential ϕ for two tandem cylinders are obtained as
(7)
ψ
=
U
y
-
U
a
2
[
y
(
x
-
l
)
2
+
y
2
+
β
2
y
(
x
+
l
)
2
+
y
2
]
+
Γ
4
π
ln
⁡
(
x
-
l
)
2
+
y
2
[
(
x
+
l
)
2
+
y
2
]
β
,
ϕ
=
U
x
+
U
a
2
[
(
x
-
l
)
(
x
-
l
)
2
+
y
2
+
β
2
(
x
+
l
)
(
x
+
l
)
2
+
y
2
]
-
Γ
2
π
{
arg
[
(
x
-
l
)
+
y
i
]
-
arg
[
(
x
+
l
)
+
y
i
]
β
}
.The equations in (7) are the potential and streamline functions that describe the potential and streamlines of a pair of unequal cylinder in tandem configuration, and they are functions of incident flow velocity U, spacing ratio L
/
D, and diameter ratio d
/
D.Whenβ
=
a
′
/
a
=
1.0, that is, the two tandem cylinders have equal diameter, the velocity function ψ
E and velocity potential ϕ
E for the two cylinders may be simplified as follows:
(8)
ψ
E
=
U
y
-
U
a
2
[
y
(
x
-
l
)
2
+
y
2
+
y
(
x
+
l
)
2
+
y
2
]
+
Γ
4
π
ln
⁡
(
x
-
l
)
2
+
y
2
(
x
+
l
)
2
+
y
2
,
ϕ
E
=
U
x
+
U
a
2
[
(
x
-
l
)
(
x
-
l
)
2
+
y
2
+
(
x
+
l
)
(
x
+
l
)
2
+
y
2
]
-
Γ
2
π
{
arg
[
(
x
-
l
)
+
y
i
]
-
arg
[
(
x
+
l
)
+
y
i
]
}
.


## 2.1. Flow Past a Single Cylinder

When an ideal fluid flows past a single circular cylinder, the streamlines and velocity potentials form the same pattern as that of a doublet placed in a constant uniform flow as suggested by Acheson [24]. Thus, the complex potential function for a uniform flow past a single cylinder may be expressed as
(1)
W
(
z
)
=
V
(
z
+
1
z
)
,
where V is the flow velocity parallel to the x-axis and z is the complex variable. For flow past a cylinder, the uniform flow field is represented by the flow velocity V. The streamline patterns obtained from the potential theory based on 1
/
V are found to be appropriate. Hence, the uniform flow field is represented by U
=
1
/
V in this study, and the potential function is expressed as ψ
(
U
)
=
ψ
(
1
/
V
).The complex flow past a circular cylinder over thez-plane can be constructed by superposition of potential flow elements. In this case, the authors use a uniform flow and vortex element with certain circulation strength. Here, the vortex element in the potential flow is represented by a rotating body of a certain diameter D and circulation strength Γ. The circulation is either clockwise or counterclockwise depending on the sign assigned to Γ. Consequently, the complex potentials for flow past a single circular cylinder can be visualized as being the resultant of a uniform flow and a rotating vortex element and can be written as
(2)
W
(
z
)
=
U
z
+
U
a
2
z
+
i
Γ
2
π
ln
⁡
z
.Consider the wake flow patterns captured using the particle image velocimetry (PIV) technique and shown in Figures3(a)–3(c); the gyros in the flow are vortices of certain strength in the flow field. These gyros could be represented by placing a circulation element of appropriate strength in the flow, with the corresponding streamlines presented in Figures 3(d)–3(f). The potential equation corresponding to these streamlines is as follows:
(3)
W
(
z
)
=
F
1
(
z
)
+
F
2
(
z
)
+
F
3
(
z
)
+
F
4
(
z
)
+
F
5
(
z
)
=
U
z
+
U
a
2
[
β
1
z
-
L
1
+
β
2
z
-
L
2
+
β
3
z
-
L
3
+
β
4
z
-
L
4
+
β
5
z
-
L
5
]
+
i
Γ
2
π
[
ln
⁡
(
z
-
L
1
)
α
1
+
ln
⁡
(
z
-
L
2
)
α
2
+
ln
⁡
(
z
-
L
3
)
α
3
+
ln
⁡
(
z
-
L
4
)
α
4
+
ln
⁡
(
z
-
L
5
)
α
5
]
,
where F
j
(
z
) represents the complex potential for each circulation element j (j
=
1,2
,
3,4
,
5). β
j and α
j represent the diameter ratio and circulation strength ratio, respectively, for each element. It can be seen from Figure 3 that the principal features of the wake flow patterns could be reproduced by using circulation elements to represent the gyros.Figure 3
Comparison of instantaneous flow patterns around a single cylinder at different Re obtained using PIV (a)–(c), and streamline flow patterns based on potential flow (d)–(f).While the patterns generated reflect physically meaningful flow phenomenon, it is recognized that the above approach requires prior knowledge of the location and strength of the circulations.


## 2.2. Flow Past a Cylinder-Pair in Tandem Configuration

For flow past two tandem cylinders as shown in Figure4, the complex flow can be constructed by superposition of a uniform flow, two vortex elements, and two pertinent circulation strengths with opposite signs. Consider the diameter ratio of the upstream to downstream cylinders is β
=
a
′
/
a (a
′ and a are the radii of the upstream and downstream cylinders, resp.) and circulation at the upstream cylinder is Γ
′
=
β
Γ (where Γ is the circulation of the downstream cylinder). The complex potentials for the upstream and downstream cylinders located, respectively, at z
=
-
l and z
=
l may be written as
(4)
F
1
(
z
)
=
U
z
+
U
(
β
a
)
2
(
z
+
l
)
-
i
β
Γ
2
π
ln
⁡
(
z
+
l
)
,
F
2
(
z
)
=
U
z
+
U
a
2
(
z
-
l
)
+
i
Γ
2
π
ln
⁡
(
z
-
l
)
.Figure 4
A schematic diagram of two cylinders with unequal diameters in tandem arrangement. “O” is the origin. The size of the cylinder is reflected by the strength of circulation. Here, D
=
2
a and d
=
2
a
′, for a
>
a
′.Then, the resultant complex potentialW
(
z
) for two tandem cylinders with unequal diameters in a steady flow (see Figure 4) is
(5)
W
(
z
)
=
F
1
(
z
)
+
F
2
(
z
)
=
U
z
+
U
a
2
(
z
-
l
)
+
U
(
β
a
)
2
(
z
+
l
)
+
i
Γ
2
π
ln
⁡
⁡
(
z
-
l
)
-
i
β
Γ
2
π
ln
⁡
⁡
(
z
+
l
)
=
U
[
z
+
a
2
(
z
-
l
)
+
(
β
a
)
2
(
z
+
l
)
]
+
i
Γ
2
π
ln
⁡
[
(
z
-
l
)
-
β
ln
⁡
(
z
+
l
)
]
.Puttingz
=
x
+
i
y into (5), the complex potential may be expressed as
(6)
W
(
z
)
=
F
1
(
z
)
+
F
2
(
z
)
=
U
x
+
U
a
2
[
(
x
-
l
)
(
x
-
l
)
2
+
y
2
+
β
2
(
x
+
l
)
(
x
+
l
)
2
+
y
2
]
-
Γ
2
π
{
arg
[
(
x
-
l
)
+
y
i
]
-
arg
[
(
x
+
l
)
+
y
i
]
β
}
+
i
{
(
x
-
l
)
2
[
(
x
-
l
)
2
+
y
2
]
β
U
y
-
U
a
2
[
y
(
x
-
l
)
2
+
y
2
+
β
2
y
(
x
+
l
)
2
+
y
2
]
+
Γ
4
π
ln
⁡
(
x
-
l
)
2
+
y
2
[
(
x
-
l
)
2
+
y
2
]
β
}
.
The velocity function ψ and velocity potential ϕ for two tandem cylinders are obtained as
(7)
ψ
=
U
y
-
U
a
2
[
y
(
x
-
l
)
2
+
y
2
+
β
2
y
(
x
+
l
)
2
+
y
2
]
+
Γ
4
π
ln
⁡
(
x
-
l
)
2
+
y
2
[
(
x
+
l
)
2
+
y
2
]
β
,
ϕ
=
U
x
+
U
a
2
[
(
x
-
l
)
(
x
-
l
)
2
+
y
2
+
β
2
(
x
+
l
)
(
x
+
l
)
2
+
y
2
]
-
Γ
2
π
{
arg
[
(
x
-
l
)
+
y
i
]
-
arg
[
(
x
+
l
)
+
y
i
]
β
}
.The equations in (7) are the potential and streamline functions that describe the potential and streamlines of a pair of unequal cylinder in tandem configuration, and they are functions of incident flow velocity U, spacing ratio L
/
D, and diameter ratio d
/
D.Whenβ
=
a
′
/
a
=
1.0, that is, the two tandem cylinders have equal diameter, the velocity function ψ
E and velocity potential ϕ
E for the two cylinders may be simplified as follows:
(8)
ψ
E
=
U
y
-
U
a
2
[
y
(
x
-
l
)
2
+
y
2
+
y
(
x
+
l
)
2
+
y
2
]
+
Γ
4
π
ln
⁡
(
x
-
l
)
2
+
y
2
(
x
+
l
)
2
+
y
2
,
ϕ
E
=
U
x
+
U
a
2
[
(
x
-
l
)
(
x
-
l
)
2
+
y
2
+
(
x
+
l
)
(
x
+
l
)
2
+
y
2
]
-
Γ
2
π
{
arg
[
(
x
-
l
)
+
y
i
]
-
arg
[
(
x
+
l
)
+
y
i
]
}
.


## 3. Results and Discussions

### 3.1. Cylinder-Pair with Equal Diameter (Special Case ofV
=
1
/
U
=
1.0)

Based on the potential flow theory, and forV
=
1
/
U
=
1, the authors analyzed the streamlines around two cylinders with equal diameter in tandem at L
/
D
=
1.2 as shown in Figure 5. Figure 5(a) shows the results in the case of an upstream cylinder with a positive vortex and downstream cylinder with a negative vortex, while Figure 5(b) shows the opposite circulation arrangement. It can be seen that, for a small spacing ratio, the two cylinders behave as a single bluff body. The streamline patterns suggest that the flow engulfs the two doublets wholly, and the flow patterns are in good agreement with those observed in experiments (Igarashi [2]; Zdravkovich [11]).Streamlines around two cylinders with equal diameter in tandem atL
/
D
=
1.2. (a) Upstream cylinder with positive vortex; (b) upstream cylinder with negative vortex.
(a)
(b)As can be seen in Figure6, when L
/
D is increased to 1.42, the two cylinders still behave as a single bluff body. In Figure 6, N indicates the innermost point of the upper streamline which is attached to the cylinder with the positive vortex. The counterpart of N, N
′, lies on the streamline attached to the cylinder with the negative vortex. Figure 7 is an enlarged view of the streamlines in the vicinity of N (and N
′), showing the nearest streamlines around the surface of two cylinders for different L
/
D and the pair of tongue-like features near the cylinders. The solid line represents the upper set of streamlines and the dotted line represents the lower set. A sudden change in the flow patterns occurs when L
/
D is in the neighborhood of 1.42. At L
/
D
=
1.43, the tongue disappears and is replaced by streamlines that cut diagonally across the gap between the two cylinders and cross the axis joining the two cylinder centres, indicating the end of type 1 interference and the onset of type 2 interference.Streamlines around two tandem cylinders atL
/
D
=
1.42. (a) Upstream cylinder with positive vortex; (b) upstream cylinder with negative vortex.
(a)
(b)Figure 7
The closest streamlines to the surface of two tandem cylinders at differentL
/
D; the solid line represents the positive value and dotted line represents the negative value. A magnified view is shown at the lower right corner.WhenL
/
D
>
1.42, it can be seen from Figure 8 that the streamline patterns in the gap between the two cylinders are distinctly different from those at smaller L
/
D (see Figures 5 and 6). In Figures 8(a) and 8(b), the streamlines cross diagonally in the region between the two cylinders. This is an interesting observation, suggesting that the alternating reattachment is the result of synchronous alternating circulation about the upstream and downstream cylinders. These intersecting patterns clearly depict the alternating attachment of the shear layers between the cylinders observed at intermediate spacing ratios as reported in the literature (Gao et al. [17]), where the vortices shed from the upstream cylinder reattach on the upper/lower surface of the downstream cylinder in an alternating clockwise and counterclockwise manner. From the potential flow theory formulation, one can see that the phenomenon of alternating reattachment on the downstream cylinder is the result of alternating positive and negative circulations of the doublet-pair.Streamlines around two tandem cylinders atL
/
D
=
2.0. (a) Upstream cylinder with positive vortex; (b) upstream cylinder with negative vortex.
(a)
(b)The slope of the diagonally crossing streamline nearest the surface of the upstream cylinder is considered here as a means of evaluating the range ofL
/
D for which the onset of alternating reattachment flow patterns is observed. As shown in Figure 9, the slope k is defined as the slope of the line joining points P
1 and P
2, where P
1 and P
2 are the horizontal and vertical “intercepts,” respectively, of the above-mentioned streamline.Figure 9
Identification of slope of the streamline nearest to the surface of upstream cylinder.The slopek decreases with L
/
D and depicts an “inflection point” in the vicinity of L
/
D
=
5.5, as shown in Figure 10. The authors hypothesize that there exists a certain slope k
cri at which there is a transition from the pattern of alternating reattachment to the next pattern where the wakes behind the two cylinders are independent of each other. In other words, when k
≥
k
cri (or equivalently L
/
D
≤
(
L
/
D
)
IIcri), the phenomenon of alternating reattachment on the downstream cylinder occurs, and when k
<
k
cri (or equivalently L
/
D
>
(
L
/
D
)
IIcri), the flow patterns correspond to those of two independent cylinder wakes. It can be seen from Figure 10 that, at L
/
D
=
5.3, the curve of k versus L
/
D has a sudden change in slope. One would expect a transition of flow patterns at this critical spacing ratio. This observation is indeed consistent with experimental observations (Igarashi [2]; Zdravkovich [11]).Figure 10
The relationship between slopek and L
/
D at intermediate spacing ratios.It can also be deduced from Figure10 that when L
/
D
>
5.3, the interactions between the two cylinders are rather weak and the two cylinder wakes are independent. Short of any other criterion, L
/
D
=
5.3 may be adopted as the threshold condition for independent wake flow patterns of two equal cylinders in tandem arrangement.The value ofL
/
D for which the two wake flow patterns are independent of each other would be different for unequal cylinder diameters and different incident flow speed. This situation is deliberated later.


### 3.2. Cylinder-Pair with Unequal Diameters (Special Case ofV
=
1
/
U
=
1)

In this section, the authors consider the situation whered
/
D is less than 1 and aim to determine if there is a unique relationship between the diameter ratio and the threshold spacing ratio for the onset of type 2 interference. In order to investigate the relationship between d
/
D and L
/
D, the authors examined the streamline patterns around a pair of cylinders with various d
/
D and L
/
D, drawing on the analogy with the potential flow theory as before.Since the streamlines patterns formed when the positive and negative vortices are reversed are mirror images of each other, only the case with the upstream cylinder having a positive vortex and the downstream cylinder a negative vortex needs to be considered hereafter.Figure11 shows the streamlines around two tandem cylinders with different diameter ratios d
/
D
=
1
/
3, 2
/
3, and 1 at a spacing ratio L
/
D
=
1.0. It can be seen that, at this spacing ratio, the two cylinders behave as a single bluff body in all three cases. Nevertheless, the streamline patterns are different at each diameter ratio, as can be seen from the streamlines in the region between the two cylinders.Streamlines around two tandem cylinders with different diameter ratiosd
/
D at L
/
D
=
1.0.
(a)
(b)
(c)WhenL
/
D is increased to 1.5, as shown in Figure 12, different streamline patterns are again observed at different diameter ratios d
/
D. At d
/
D
=
0.5, the “crossing” pattern between the two cylinders is observed. However, when d
/
D is increased to 1.065, the lower streamline shows the tendency of moving upwards but has not yet extended across the axis joining the two cylinder centres. This is similar to the patterns observed between two cylinders with equal diameter at small spacing ratios, whereupon the two cylinders behave as a single bluff body. In summary, it appears that a certain relationship exists between d
/
D and L
/
D to account for the same flow patterns between the two tandem cylinders of unequal diameters.Streamlines around two tandem cylinders with different diameter ratiosd
/
D at L
/
D
=
1.5.
(a)
(b)
(c)
(d)By examining the variations of this threshold spacing ratio(
L
/
D
)
cri for different diameter ratios d
/
D (see Figure 13), it can be seen that the corresponding (
L
/
D
)
cri is smaller for unequal cylinder-pairs. That the data by other authors scatter widely at d
/
D
=
1 is due to the different flow speed which could not be harmonized in this plot.Figure 13
The variation of critical spacing ratio(
L
/
D
)
cri with diameter ratio d
/
D.Figure13 shows the logarithmic relationship between the threshold spacing ratio (
L
/
D
)
cri and diameter ratio d
/
D. Note that d
/
D is always less than or equal to 1 since D is the larger of the two cylinder diameters. Also included in Figure 13 are the reported values of (
L
/
D
)
cri at which type 1 interference is ended, depending on the range of Re. It can be seen that, for the same U(or Re), the threshold spacing ratio (
L
/
D
)
cri increases with d
/
D. That is to say, the alternating attachment of streamline patterns can be observed at lower centre-to-centre spacing ratio L
/
D for smaller d
/
D, and unequal cylinders should be bundled at closer spacing. Therefore, it can be concluded that in addition to L
/
D, the diameter ratio also has a significant effect on the classification of streamline patterns. The reported data for equal diameter cylinder-pair (
d
/
D
=
1
) are visibly scattered and are logically so because these data are based on different incident flow speed.The relationship between the critical spacing ratio(
L
/
D
)
cri and diameter ratio d
/
D can be fitted based on Figure 13 as follows:
(9)
Log
(
L
D
)
cri
=
0.339
Log
(
d
D
)
+
0.164
.At diameter ratiod
/
D
=
1.0, (
L
/
D
)
cri
=
10
0.164
=
1.458 for the case V
=
1
/
U
=
1, which compares well to the reported results of Carmo et al. [8] or approach the mean value of (
L
/
D
)
cri obtained by Igarashi [1] ((
L
/
D
)
cri
=
1.03–2, depending on the range of Re) and Zdravkovich [11] ((
L
/
D
)
cri
=
1.2–1.8, depending on the range of Re).It is noted that (9) may be expressed as follows, approximating the gradient of 0.339 as 1/3:
(10)
(
L
/
D
)
(
d
/
D
)
1
/
3
=
10
0.164
=
constant
,
which indicates the functional relationship, for the special case of V
=
1
/
U
=
1.0; that is,
(11)
f
[
(
L
D
)
(
D
d
)
1
/
3
]
=
constant
.


### 3.3. Streamline Patterns as a Function of Incident Flow Velocity

In order to investigate the effects of incident flow velocity on the cylinder-pair, the authors examined the streamline patterns about the cylinder-pair at various velocitiesU(=
1
/
V in the mapping function). As an illustration, the various streamline patterns for L
/
D
=
1.42 and d
/
D
=
1 are shown in Figure 14. Different streamline patterns are observed for different U while the spacing and diameter ratios are held constant.Streamlines around two tandem cylinders with different velocityV
=
1
/
U at L
/
D
=
1.42.
(a)
(b)
(c)
(d)At smallU of the order of 0.01–1.0, no “diagonal-crossing” patterns are observed; that is, the two cylinders behave as a single bluff body. However, when U is increased to 2.0, the “crossing” pattern is observed, indicating that the reattachment phenomenon has taken place. The authors repeated the numerical experiment with various diameter ratio (
d
/
D
) and flow velocity U and determined the corresponding threshold spacing ratio (
L
/
D
)
cri. It was observed that (
L
/
D
)
cri decreases with U for a fixed d
/
D, but increases with d
/
D for a fixed velocity. This observation is in good agreement with Igarashi’s [1] suggestion that the critical spacing ratio is sensitive to Reynolds number.Figure15 shows the variation of the dimensional parameter (
L
/
D
)
(
D
/
d
)
1
/
3 as a function of dimensionless velocity in the form of a Reynolds number at different diameter ratio d
/
D, where the dimensional parameter (
L
/
D
)
(
D
/
d
)
1
/
3 is calculated from (10). This Reynolds number is defined as U
(
D
d
)
1
/
2, for which the kinematic viscosity of water is constant and left out of the expression without affecting the quality of the presentation. The root mean square value of D and d is chosen as the characteristic length. It can be seen in Figure 15 that the threshold spacing ratio approaches a constant value (~1.2) when U is large. When U is very small, the reattachment phenomenon is not observed even when L
/
D is large. In terms of physical interpretation, the authors hypothesize that when U is large, the flow is turbulent and reattachment phenomenon would be observed readily. For small U, and in the limit when the flow is laminar, the cylinder-pair would behave like a single bluff body.Figure 15
Threshold spacing ratio as a function ofU
(
D
d
)
1
/
2 at different d
/
D.It can also be seen in Figure15 that the data collapse and fit into a single curve which is fitted as
(12)
(
L
D
)
(
D
d
)
1
/
3
=
0.201
U
(
D
d
)
1
/
2
+
1.163
.Equation (12) is included in Figure 15 as the solid line.


## 3.1. Cylinder-Pair with Equal Diameter (Special Case ofV
=
1
/
U
=
1.0)

Based on the potential flow theory, and forV
=
1
/
U
=
1, the authors analyzed the streamlines around two cylinders with equal diameter in tandem at L
/
D
=
1.2 as shown in Figure 5. Figure 5(a) shows the results in the case of an upstream cylinder with a positive vortex and downstream cylinder with a negative vortex, while Figure 5(b) shows the opposite circulation arrangement. It can be seen that, for a small spacing ratio, the two cylinders behave as a single bluff body. The streamline patterns suggest that the flow engulfs the two doublets wholly, and the flow patterns are in good agreement with those observed in experiments (Igarashi [2]; Zdravkovich [11]).Streamlines around two cylinders with equal diameter in tandem atL
/
D
=
1.2. (a) Upstream cylinder with positive vortex; (b) upstream cylinder with negative vortex.
(a)
(b)As can be seen in Figure6, when L
/
D is increased to 1.42, the two cylinders still behave as a single bluff body. In Figure 6, N indicates the innermost point of the upper streamline which is attached to the cylinder with the positive vortex. The counterpart of N, N
′, lies on the streamline attached to the cylinder with the negative vortex. Figure 7 is an enlarged view of the streamlines in the vicinity of N (and N
′), showing the nearest streamlines around the surface of two cylinders for different L
/
D and the pair of tongue-like features near the cylinders. The solid line represents the upper set of streamlines and the dotted line represents the lower set. A sudden change in the flow patterns occurs when L
/
D is in the neighborhood of 1.42. At L
/
D
=
1.43, the tongue disappears and is replaced by streamlines that cut diagonally across the gap between the two cylinders and cross the axis joining the two cylinder centres, indicating the end of type 1 interference and the onset of type 2 interference.Streamlines around two tandem cylinders atL
/
D
=
1.42. (a) Upstream cylinder with positive vortex; (b) upstream cylinder with negative vortex.
(a)
(b)Figure 7
The closest streamlines to the surface of two tandem cylinders at differentL
/
D; the solid line represents the positive value and dotted line represents the negative value. A magnified view is shown at the lower right corner.WhenL
/
D
>
1.42, it can be seen from Figure 8 that the streamline patterns in the gap between the two cylinders are distinctly different from those at smaller L
/
D (see Figures 5 and 6). In Figures 8(a) and 8(b), the streamlines cross diagonally in the region between the two cylinders. This is an interesting observation, suggesting that the alternating reattachment is the result of synchronous alternating circulation about the upstream and downstream cylinders. These intersecting patterns clearly depict the alternating attachment of the shear layers between the cylinders observed at intermediate spacing ratios as reported in the literature (Gao et al. [17]), where the vortices shed from the upstream cylinder reattach on the upper/lower surface of the downstream cylinder in an alternating clockwise and counterclockwise manner. From the potential flow theory formulation, one can see that the phenomenon of alternating reattachment on the downstream cylinder is the result of alternating positive and negative circulations of the doublet-pair.Streamlines around two tandem cylinders atL
/
D
=
2.0. (a) Upstream cylinder with positive vortex; (b) upstream cylinder with negative vortex.
(a)
(b)The slope of the diagonally crossing streamline nearest the surface of the upstream cylinder is considered here as a means of evaluating the range ofL
/
D for which the onset of alternating reattachment flow patterns is observed. As shown in Figure 9, the slope k is defined as the slope of the line joining points P
1 and P
2, where P
1 and P
2 are the horizontal and vertical “intercepts,” respectively, of the above-mentioned streamline.Figure 9
Identification of slope of the streamline nearest to the surface of upstream cylinder.The slopek decreases with L
/
D and depicts an “inflection point” in the vicinity of L
/
D
=
5.5, as shown in Figure 10. The authors hypothesize that there exists a certain slope k
cri at which there is a transition from the pattern of alternating reattachment to the next pattern where the wakes behind the two cylinders are independent of each other. In other words, when k
≥
k
cri (or equivalently L
/
D
≤
(
L
/
D
)
IIcri), the phenomenon of alternating reattachment on the downstream cylinder occurs, and when k
<
k
cri (or equivalently L
/
D
>
(
L
/
D
)
IIcri), the flow patterns correspond to those of two independent cylinder wakes. It can be seen from Figure 10 that, at L
/
D
=
5.3, the curve of k versus L
/
D has a sudden change in slope. One would expect a transition of flow patterns at this critical spacing ratio. This observation is indeed consistent with experimental observations (Igarashi [2]; Zdravkovich [11]).Figure 10
The relationship between slopek and L
/
D at intermediate spacing ratios.It can also be deduced from Figure10 that when L
/
D
>
5.3, the interactions between the two cylinders are rather weak and the two cylinder wakes are independent. Short of any other criterion, L
/
D
=
5.3 may be adopted as the threshold condition for independent wake flow patterns of two equal cylinders in tandem arrangement.The value ofL
/
D for which the two wake flow patterns are independent of each other would be different for unequal cylinder diameters and different incident flow speed. This situation is deliberated later.


## 3.2. Cylinder-Pair with Unequal Diameters (Special Case ofV
=
1
/
U
=
1)

In this section, the authors consider the situation whered
/
D is less than 1 and aim to determine if there is a unique relationship between the diameter ratio and the threshold spacing ratio for the onset of type 2 interference. In order to investigate the relationship between d
/
D and L
/
D, the authors examined the streamline patterns around a pair of cylinders with various d
/
D and L
/
D, drawing on the analogy with the potential flow theory as before.Since the streamlines patterns formed when the positive and negative vortices are reversed are mirror images of each other, only the case with the upstream cylinder having a positive vortex and the downstream cylinder a negative vortex needs to be considered hereafter.Figure11 shows the streamlines around two tandem cylinders with different diameter ratios d
/
D
=
1
/
3, 2
/
3, and 1 at a spacing ratio L
/
D
=
1.0. It can be seen that, at this spacing ratio, the two cylinders behave as a single bluff body in all three cases. Nevertheless, the streamline patterns are different at each diameter ratio, as can be seen from the streamlines in the region between the two cylinders.Streamlines around two tandem cylinders with different diameter ratiosd
/
D at L
/
D
=
1.0.
(a)
(b)
(c)WhenL
/
D is increased to 1.5, as shown in Figure 12, different streamline patterns are again observed at different diameter ratios d
/
D. At d
/
D
=
0.5, the “crossing” pattern between the two cylinders is observed. However, when d
/
D is increased to 1.065, the lower streamline shows the tendency of moving upwards but has not yet extended across the axis joining the two cylinder centres. This is similar to the patterns observed between two cylinders with equal diameter at small spacing ratios, whereupon the two cylinders behave as a single bluff body. In summary, it appears that a certain relationship exists between d
/
D and L
/
D to account for the same flow patterns between the two tandem cylinders of unequal diameters.Streamlines around two tandem cylinders with different diameter ratiosd
/
D at L
/
D
=
1.5.
(a)
(b)
(c)
(d)By examining the variations of this threshold spacing ratio(
L
/
D
)
cri for different diameter ratios d
/
D (see Figure 13), it can be seen that the corresponding (
L
/
D
)
cri is smaller for unequal cylinder-pairs. That the data by other authors scatter widely at d
/
D
=
1 is due to the different flow speed which could not be harmonized in this plot.Figure 13
The variation of critical spacing ratio(
L
/
D
)
cri with diameter ratio d
/
D.Figure13 shows the logarithmic relationship between the threshold spacing ratio (
L
/
D
)
cri and diameter ratio d
/
D. Note that d
/
D is always less than or equal to 1 since D is the larger of the two cylinder diameters. Also included in Figure 13 are the reported values of (
L
/
D
)
cri at which type 1 interference is ended, depending on the range of Re. It can be seen that, for the same U(or Re), the threshold spacing ratio (
L
/
D
)
cri increases with d
/
D. That is to say, the alternating attachment of streamline patterns can be observed at lower centre-to-centre spacing ratio L
/
D for smaller d
/
D, and unequal cylinders should be bundled at closer spacing. Therefore, it can be concluded that in addition to L
/
D, the diameter ratio also has a significant effect on the classification of streamline patterns. The reported data for equal diameter cylinder-pair (
d
/
D
=
1
) are visibly scattered and are logically so because these data are based on different incident flow speed.The relationship between the critical spacing ratio(
L
/
D
)
cri and diameter ratio d
/
D can be fitted based on Figure 13 as follows:
(9)
Log
(
L
D
)
cri
=
0.339
Log
(
d
D
)
+
0.164
.At diameter ratiod
/
D
=
1.0, (
L
/
D
)
cri
=
10
0.164
=
1.458 for the case V
=
1
/
U
=
1, which compares well to the reported results of Carmo et al. [8] or approach the mean value of (
L
/
D
)
cri obtained by Igarashi [1] ((
L
/
D
)
cri
=
1.03–2, depending on the range of Re) and Zdravkovich [11] ((
L
/
D
)
cri
=
1.2–1.8, depending on the range of Re).It is noted that (9) may be expressed as follows, approximating the gradient of 0.339 as 1/3:
(10)
(
L
/
D
)
(
d
/
D
)
1
/
3
=
10
0.164
=
constant
,
which indicates the functional relationship, for the special case of V
=
1
/
U
=
1.0; that is,
(11)
f
[
(
L
D
)
(
D
d
)
1
/
3
]
=
constant
.


## 3.3. Streamline Patterns as a Function of Incident Flow Velocity

In order to investigate the effects of incident flow velocity on the cylinder-pair, the authors examined the streamline patterns about the cylinder-pair at various velocitiesU(=
1
/
V in the mapping function). As an illustration, the various streamline patterns for L
/
D
=
1.42 and d
/
D
=
1 are shown in Figure 14. Different streamline patterns are observed for different U while the spacing and diameter ratios are held constant.Streamlines around two tandem cylinders with different velocityV
=
1
/
U at L
/
D
=
1.42.
(a)
(b)
(c)
(d)At smallU of the order of 0.01–1.0, no “diagonal-crossing” patterns are observed; that is, the two cylinders behave as a single bluff body. However, when U is increased to 2.0, the “crossing” pattern is observed, indicating that the reattachment phenomenon has taken place. The authors repeated the numerical experiment with various diameter ratio (
d
/
D
) and flow velocity U and determined the corresponding threshold spacing ratio (
L
/
D
)
cri. It was observed that (
L
/
D
)
cri decreases with U for a fixed d
/
D, but increases with d
/
D for a fixed velocity. This observation is in good agreement with Igarashi’s [1] suggestion that the critical spacing ratio is sensitive to Reynolds number.Figure15 shows the variation of the dimensional parameter (
L
/
D
)
(
D
/
d
)
1
/
3 as a function of dimensionless velocity in the form of a Reynolds number at different diameter ratio d
/
D, where the dimensional parameter (
L
/
D
)
(
D
/
d
)
1
/
3 is calculated from (10). This Reynolds number is defined as U
(
D
d
)
1
/
2, for which the kinematic viscosity of water is constant and left out of the expression without affecting the quality of the presentation. The root mean square value of D and d is chosen as the characteristic length. It can be seen in Figure 15 that the threshold spacing ratio approaches a constant value (~1.2) when U is large. When U is very small, the reattachment phenomenon is not observed even when L
/
D is large. In terms of physical interpretation, the authors hypothesize that when U is large, the flow is turbulent and reattachment phenomenon would be observed readily. For small U, and in the limit when the flow is laminar, the cylinder-pair would behave like a single bluff body.Figure 15
Threshold spacing ratio as a function ofU
(
D
d
)
1
/
2 at different d
/
D.It can also be seen in Figure15 that the data collapse and fit into a single curve which is fitted as
(12)
(
L
D
)
(
D
d
)
1
/
3
=
0.201
U
(
D
d
)
1
/
2
+
1.163
.Equation (12) is included in Figure 15 as the solid line.


## 4. Conclusions

Drawing on the analogy of potential flow theory, the flow patterns that are generally seen in experiments for different centre-to-centre spacing ratios and diameter ratios for two cylinders in tandem configuration are reproduced. The authors found that the flow speedU, spacing ratio L
/
D, and diameter ratio d
/
D have significant effects on the flow patterns for flow past a cylinder-pair in tandem arrangement. A cubic relationship f
(
L
/
D
,
d
/
D
,
U
) of the threshold condition for the wake transition from a single bluff body to a reattachment behavior is established. This study has also found that the occurrence of alternating reattachment phenomenon results from the alternating positive and negative circulations of the doublet-pair. Two cylinders of unequal diameters need to be placed closer to one another (compared to cylinders of equal diameter) in order to avoid the occurrence of vortex shedding between the cylinders.



---

*Source: 495179-2014-03-05.xml*