Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 495179, 13 pages

http://dx.doi.org/10.1155/2014/495179

Research Article

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

Yangyang Gao,1,2 Danielle S. Tan,³ Zhiyong Hao,⁴ Xikun Wang,⁵ and Soon Keat Tan⁶

¹Ocean College, Zhejiang University, Hangzhou 310058, China

²The State Key Laboratory of Fluid Power Transmission and Control, Zhejiang University, Hangzhou 310058, China

³Cornell University, Ithaca, NY 14853, USA

⁴College of Logistics Engineering, Shanghai Maritime University, Shanghai 201306, China

⁵Maritime Research Centre and School of Civil and Environmental Engineering, College of Engineering, Nanyang Technological University, Singapore, Singapore, 639798

⁶Nanyang Environment and Water Research Institute and School of Civil and Environmental Engineering, Nanyang Technological University, Singapore, Singapore, 639798

Received 25 August 2013; Revised 31 December 2013; Accepted 14 January 2014; Published 5 March 2014

Academic Editor: Shihua Li

Copyright © 2014 Yangyang Gao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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, (diameters and , ), spacing ratio, (centre-to-centre distance, , to cylinder diameter, ), and equivalent incident flow speed, A unique functional relationship (, , ) of this threshold condition is established.

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 ( or Re) and depending on the spacing ratio and diameter ratio , these three types of flow interference can also be observed, as illustrated in Figure 1. Here, for a fixed incident flow speed and , the reattachment and alternating vortex shedding patterns in the wake of the upstream cylinder begin to appear at , with the cylinder-pair behaving like a single bluff body for . This observation indicates that the value 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 and , 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.

(a)

(b)

Figure 1: Instantaneous flow patterns around two tandem cylinders with at different , for a fixed incident flow speed (Re). (a) Laboratory observation based on particle image velocimetry (PIV) and (b) schematic flow structure.

(a)

(b)

(c)

Figure 2: Instantaneous flow patterns around two tandem cylinders at for a fixed but different Re.

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 of , or said expression, as a function of and speed of the cross-flow , or Re. Reported spacing ratios 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: ,
Re

Igarashi [1]

1.03–1.91
2.09–3.53

Zdravkovich [11]

Not given
1–1.2~1.8
1.2~1.8–3.4~3.8

Ljungkrona and Sundén [3]

1.25
2.0, 4.0

Xu and Zhou [6]

1.0–2.0
2.0–3.0

Zhou and Yiu [7]

7000
1.3
2.5, 4.0

Gao et al. [16]

1200–4800
1.2
2.4

: centre-to-centre spacing ratio range for single bluff body pattern.

: 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 Figures 1 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

where is the flow velocity parallel to the -axis and is the complex variable. For flow past a cylinder, the uniform flow field is represented by the flow velocity . The streamline patterns obtained from the potential theory based on are found to be appropriate. Hence, the uniform flow field is represented by in this study, and the potential function is expressed as .

The complex flow past a circular cylinder over the -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 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

Consider the wake flow patterns captured using the particle image velocimetry (PIV) technique and shown in Figures 3(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:

where represents the complex potential for each circulation element (). and 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 Figure 4, 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 ( and 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 and may be written as

Figure 4: A schematic diagram of two cylinders with unequal diameters in tandem arrangement. “” is the origin. The size of the cylinder is reflected by the strength of circulation. Here, and , for .

Then, the resultant complex potential for two tandem cylinders with unequal diameters in a steady flow (see Figure 4) is

Putting into (5), the complex potential may be expressed as

The velocity function and velocity potential for two tandem cylinders are obtained as

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 , spacing ratio , and diameter ratio .

When , that is, the two tandem cylinders have equal diameter, the velocity function and velocity potential for the two cylinders may be simplified as follows:

3. Results and Discussions

3.1. Cylinder-Pair with Equal Diameter (Special Case of )

Based on the potential flow theory, and for , the authors analyzed the streamlines around two cylinders with equal diameter in tandem at 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]).

(a)

(b)

Figure 5: Streamlines around two cylinders with equal diameter in tandem at . (a) Upstream cylinder with positive vortex; (b) upstream cylinder with negative vortex.

As can be seen in Figure 6, when is increased to 1.42, the two cylinders still behave as a single bluff body. In Figure 6, indicates the innermost point of the upper streamline which is attached to the cylinder with the positive vortex. The counterpart of , , 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 (and ), showing the nearest streamlines around the surface of two cylinders for different 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 is in the neighborhood of 1.42. At , 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.

(a)

(b)

Figure 6: Streamlines around two tandem cylinders at . (a) Upstream cylinder with positive vortex; (b) upstream cylinder with negative vortex.

Figure 7: The closest streamlines to the surface of two tandem cylinders at different ; the solid line represents the positive value and dotted line represents the negative value. A magnified view is shown at the lower right corner.

When , 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 (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.

(a)

(b)

Figure 8: Streamlines around two tandem cylinders at . (a) Upstream cylinder with positive vortex; (b) upstream cylinder with negative vortex.

The slope of the diagonally crossing streamline nearest the surface of the upstream cylinder is considered here as a means of evaluating the range of for which the onset of alternating reattachment flow patterns is observed. As shown in Figure 9, the slope is defined as the slope of the line joining points and , where and 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 slope decreases with and depicts an “inflection point” in the vicinity of , as shown in Figure 10. The authors hypothesize that there exists a certain slope 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 (or equivalently ), the phenomenon of alternating reattachment on the downstream cylinder occurs, and when (or equivalently ), the flow patterns correspond to those of two independent cylinder wakes. It can be seen from Figure 10 that, at , the curve of versus 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 slope and at intermediate spacing ratios.

It can also be deduced from Figure 10 that when , the interactions between the two cylinders are rather weak and the two cylinder wakes are independent. Short of any other criterion, may be adopted as the threshold condition for independent wake flow patterns of two equal cylinders in tandem arrangement.

The value of 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 of )

In this section, the authors consider the situation where 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 and , the authors examined the streamline patterns around a pair of cylinders with various and , 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.

Figure 11 shows the streamlines around two tandem cylinders with different diameter ratios , , and 1 at a spacing ratio . 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.

(a)

(b)

(c)

Figure 11: Streamlines around two tandem cylinders with different diameter ratios at .

When is increased to 1.5, as shown in Figure 12, different streamline patterns are again observed at different diameter ratios . At , the “crossing” pattern between the two cylinders is observed. However, when 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 and to account for the same flow patterns between the two tandem cylinders of unequal diameters.

(a)

(b)

(c)

(d)

Figure 12: Streamlines around two tandem cylinders with different diameter ratios at .

By examining the variations of this threshold spacing ratio for different diameter ratios (see Figure 13), it can be seen that the corresponding is smaller for unequal cylinder-pairs. That the data by other authors scatter widely at is due to the different flow speed which could not be harmonized in this plot.

Figure 13: The variation of critical spacing ratio with diameter ratio .

Figure 13 shows the logarithmic relationship between the threshold spacing ratio and diameter ratio . Note that is always less than or equal to 1 since is the larger of the two cylinder diameters. Also included in Figure 13 are the reported values of at which type 1 interference is ended, depending on the range of Re. It can be seen that, for the same (or Re), the threshold spacing ratio increases with . That is to say, the alternating attachment of streamline patterns can be observed at lower centre-to-centre spacing ratio for smaller , and unequal cylinders should be bundled at closer spacing. Therefore, it can be concluded that in addition to , the diameter ratio also has a significant effect on the classification of streamline patterns. The reported data for equal diameter cylinder-pair are visibly scattered and are logically so because these data are based on different incident flow speed.

The relationship between the critical spacing ratio and diameter ratio can be fitted based on Figure 13 as follows:

At diameter ratio , for the case , which compares well to the reported results of Carmo et al. [8] or approach the mean value of obtained by Igarashi [1] (–2, depending on the range of Re) and Zdravkovich [11] (–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:

which indicates the functional relationship, for the special case of ; that is,

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 velocities ( in the mapping function). As an illustration, the various streamline patterns for and are shown in Figure 14. Different streamline patterns are observed for different while the spacing and diameter ratios are held constant.

(a)

(b)

(c)

(d)

Figure 14: Streamlines around two tandem cylinders with different velocity at .

At small 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 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 and flow velocity and determined the corresponding threshold spacing ratio . It was observed that decreases with for a fixed , but increases with 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.

Figure 15 shows the variation of the dimensional parameter as a function of dimensionless velocity in the form of a Reynolds number at different diameter ratio , where the dimensional parameter is calculated from (10). This Reynolds number is defined as , 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 and 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 is large. When is very small, the reattachment phenomenon is not observed even when is large. In terms of physical interpretation, the authors hypothesize that when is large, the flow is turbulent and reattachment phenomenon would be observed readily. For small , 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 of at different .

It can also be seen in Figure 15 that the data collapse and fit into a single curve which is fitted as

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 speed , spacing ratio , and diameter ratio have significant effects on the flow patterns for flow past a cylinder-pair in tandem arrangement. A cubic relationship 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.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research was supported by the Open Foundation of the State Key Laboratory of Fluid Power Transmission and Control (GZKF-201310), the State Key Laboratory of Satellite Ocean Environment Dynamics (Second Institute of Oceanography, SOA), China, and the National Research Foundation, Singapore (CRP-5-2009-01); funding support from the Maritime and Port Authority of Singapore (MPA) is gratefully acknowledged. Zhiyong Hao also would like to acknowledge the support from Science and Technology Commission of Shanghai Municipality (Pujiang Program, Grant no. 10PJ1404700), Shanghai Education Committee (Program Grant no. 12ZZ149 and J50604), and Shanghai Maritime University (Science & Technology Program, Grant no. 20100089).

References

T. Igarashi, “Characteristics of the flow around two circular cylinders arranged in tandem, 1st report,” Bulletin of the JSME B, vol. 24, no. 188, pp. 323–331, 1981.

T. Igarashi, “Characteristics of the flow around two circular cylinders arranged in tandem, 2nd report,” Bulletin of the JSME B, vol. 27, no. 233, pp. 2380–2387, 1984.

L. Ljungkrona and B. Sundén, “Flow visualization and surface pressure measurement on two tubes in an inline arrangement,” Experimental Thermal and Fluid Science, vol. 6, no. 1, pp. 15–27, 1993.

D. Sumner, S. J. Price, and M. P. Païdoussis, “Tandem cylinders in impulsively started flow,” Journal of Fluids and Structures, vol. 13, no. 7-8, pp. 955–965, 1999.

J.-C. Lin, Y. Yang, and D. Rockwell, “Flow past two cylinders in tandem: instantaneous and averaged flow structure,” Journal of Fluids and Structures, vol. 16, no. 8, pp. 1059–1071, 2002.

G. Xu and Y. Zhou, “Strouhal numbers in the wake of two inline cylinders,” Experiments in Fluids, vol. 37, no. 2, pp. 248–256, 2004.

Y. Zhou and M. W. Yiu, “Flow structure, momentum and heat transport in a two-tandem-cylinder wake,” Journal of Fluid Mechanics, vol. 548, pp. 17–48, 2006.

B. S. Carmo, J. R. Meneghini, and S. J. Sherwin, “Possible states in the flow around two circular cylinders in tandem with separations in the vicinity of the drag inversion spacing,” Physics of Fluids, vol. 22, no. 5, Article ID 054101, 2010.

B. S. Carmo, J. R. Meneghini, and S. J. Sherwin, “Secondary instabilities in the flow around two circular cylinders in tandem,” Journal of Fluid Mechanics, vol. 644, pp. 395–431, 2010.

D. Sumner, “Two circular cylinders in cross-flow: a review,” Journal of Fluids and Structures, vol. 26, no. 6, pp. 849–899, 2010.

M. M. Zdravkovich, “The effects of interference between circular cylinders in cross flow,” Journal of Fluids and Structures, vol. 1, no. 2, pp. 239–261, 1987.

T. Igarashi, “Characteristics of the flow around two circular cylinders of different diameters arranged in tandem,” Bulletin of the JSME B, vol. 25, no. 201, pp. 349–357, 1982.

M. Mahbub Alam and Y. Zhou, “Strouhal numbers, forces and flow structures around two tandem cylinders of different diameters,” Journal of Fluids and Structures, vol. 24, no. 4, pp. 505–526, 2008.

M. Zhao, L. Cheng, B. Teng, and D. F. Liang, “Numerical simulation of viscous flow past two circular cylinders of different diameters,” Applied Ocean Research, vol. 27, no. 1, pp. 39–55, 2005.

M. Zhao, L. Cheng, B. Teng, and G. Dong, “Hydrodynamic forces on dual cylinders of different diameters in steady currents,” Journal of Fluids and Structures, vol. 23, no. 1, pp. 59–83, 2007.

Y. Gao, S. Etienne, D. Yu, X. Wang, and S. Tan, “Bi-stable flow around tandem cylinders of different diameters at low Reynolds number,” Fluid Dynamics Research, vol. 43, no. 5, Article ID 055506, 2011.

Y. Gao, S. Tan, E. Stephane, and X. Wang, “Experimental study on the flow around two tandem cylinders of different diameters,” Internal Report of MRC, 2012.

A. Millward, D. Nicolaou, and S. G. Rigby, “Numerical modeling of the water flow around a fast ship with a transom stern,” International Journal of Maritime Engineering, The Royal Institute of Naval Architects, pp. 181–194, 2003.

M. Shahjada Tarafder and K. Suzuki, “Numerical calculation of free-surface potential flow around a ship using the modified Rankine source panel method,” Ocean Engineering, vol. 35, no. 5-6, pp. 536–544, 2008.

Y. Ahmed and C. Guedes Soares, “Simulation of free surface flow around a VLCC hull using viscous and potential flow methods,” Ocean Engineering, vol. 36, no. 9-10, pp. 691–696, 2009.

E. R. Johnson and N. R. McDonald, “The motion of a vortex near two circular cylinders,” Proceedings of The Royal Society of London A, vol. 460, no. 2044, pp. 939–954, 2004.

D. A. Burton, J. Gratus, and R. W. Tucker, “Hydrodynamic forces on two moving discs,” Theoretical and Applied Mechanics, vol. 31, no. 2, pp. 153–188, 2004.

D. G. Crowdy, “Analytical solutions for uniform potential flow past multiple cylinders,” European Journal of Mechanics B, vol. 25, no. 4, pp. 459–470, 2006.

D. J. Acheson, Elementary Fluid Dynamics, Oxford University Press, Oxford, UK, 1990.