MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2015/568176 568176 Research Article Computation of Flow Past an In-Line Oscillating Circular Cylinder and a Stationary Cylinder in Tandem Using a CIP-Based Model Fu Yingnan 1 http://orcid.org/0000-0002-1392-139X Zhao Xizeng 1, 2 Wang Xinggang 2 Cao Feifeng 3 Zhao Ming 1 Ocean College Zhejiang University Hangzhou 310058 China zju.edu.cn 2 State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering Nanjing Hydraulic Research Institute Nanjing 210029 China nhri.cn 3 Institute of Harbor-Channel and Coastal Engineering, Department of Civil Engineering, Zhejiang University of Technology, Hangzhou 310014 China 2015 18102015 2015 13 04 2015 21 05 2015 27 05 2015 18102015 2015 Copyright © 2015 Yingnan Fu 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.

Viscous flow past an upstream in-line forced oscillating circular cylinder with a stationary cylinder downstream at Reynolds number of 100 is investigated using a CIP model. The model is established in a Cartesian coordinate system using a high-order difference method to discretise the Navier-Stokes equations. The fluid-structure interaction is treated as a multiphase flow with fluid and solid phases solved simultaneously. An immersed boundary method is used to deal with the fluid-body coupling. The CFD model is firstly applied to the computation of flow past a fixed circular cylinder for its validation; then flow over two stationary tandem cylinders is investigated and good agreements are obtained comparing with existing ones. Computations are then performed with flow past two tandem cylinders with an upstream in-line oscillating cylinder with a small spacing L=2D. Considerable attention is paid to the spectrum characteristics and vortex modes.

1. Introduction

Flow over multiple bluff bodies is one of the classical issues in fluid mechanics. It is of great importance in engineering applications, such as heat exchangers, adjacent tall buildings, and piles of offshore platforms. The interaction between the multiple bodies and the wake is much more complicated than that of flow past a single body.

The tandem arrangement is one of the simplest configurations of the multiple bluff bodies. Extensive experimental and numerical studies have been reported about flow past two stationary tandem circular cylinders. Zdravkovich (1977)  revealed a discontinuous “jump” at some critical spacing in his measurements and found the discontinuity is caused by the abrupt change from one stable flow pattern to another at the critical spacing. Li et al. (1991)  utilized a finite element program to simulate flow past two tandem circular cylinders at different intervals for Re=100. Discontinuity changes in the flow pattern, the Strouhal number, and the pressure distribution were detected. Mittal et al. (1997)  used a finite element formation to simulate three configurations at Re=100 and 1000. It was observed that the flow characteristics depend strongly on the arrangements of the cylinders and the Reynolds number. Meneghini et al. (2001)  presented the results of flow around tandem circular cylinders for Reynolds number of 200. Vorticity contours and force tine histories were displayed. Jester and Kallinderis (2003)  employed a second-order upwind scheme to investigate this problem. Qualitative and quantitative comparisons were made with published experimental data. The mean and fluctuating lift and drag coefficients were recorded for center-to-center cylinder spacings between 2 and 10 diameters in Sharman et al. (2005) . Results showed that the critical spacing was between 3.75 and 4 diameters. Mussa et al. (2009)  applied the LBE (Lattice Boltzmann Equation) method with MRT (multiple relaxation time) to simulate the flow past two tandem cylinders in an unbounded domain with Re=100 and got perfect results. Numerical studies of Liu et al. (2013)  showed that separated point of the shear layer backwards when L/D increases because of the interference of downstream cylinder.

The study concerning oscillating circular cylinders in tandem arrangement is relatively rare. Direct numerical simulation was carried out by Li et al. (1992)  to study the response of a transversely oscillating cylinder in uniform flow and in the wake of an upstream cylinder. They found that the response of the cylinder wake was either a periodic (lock-in) or a quasi-periodic (nonlock-in) state. Mahir and Rockwell (1996)  investigated the wakes between and downstream of a tandem arrangement of two transversely oscillating cylinders with a phase difference experimentally. Their results showed that the lock-in range is much wider than that of a corresponding single cylinder. Following the experiment of Mahir and Rockwell (1996) , Papaioannou et al. (2006)  conducted a numerical study to investigate the holes in the Arnold tongues. Zou et al. (2004)  simulated the flow past a forced system of a stationary circular cylinder and a forced transversely oscillating cylinder downstream. The vortex switch and the change of vortex modes were observed in the “lock-in” region. Prasanth and Mittal (2009)  investigated the flow-induced oscillation of two circular cylinders in tandem arrangement at low Reynolds number numerically. The phase between the transverse responses of the two cylinders was studied. Nonlinear responses to a transversely oscillating cylinder in the wake of a stationary upstream cylinder were studied by Yang and Zheng (2010) , based on velocity phase portraits. They categorized the nonlinear response states into three states: lock-in, transitional, and quasi-periodic. Huera-Huarte and Bearman (2011)  carried out an experiment to investigate vortex and wake-induced vibrations of two tandem flexible circular cylinders. Different response regimes have been identified based on the type of oscillations exhibited by the cylinders: vortex-induced (VIV), wake-induced (WIV), or combinations of both. Bao et al. (2012)  employed LBM (Lattice Boltzmann Method) to study flow around a transversely oscillating cylinder in tandem with a stationary cylinder downstream. The influences of spacing, oscillation frequency, and amplitude were investigated. Tao et al. (2013)  utilized LBM to study flow past a pair of rotating circular cylinders in tandem arrangement at Re=100. Flow past two cylinders in tandem arrangement under forced vibration has been studied experimentally by Yang et al. (2014)  employing the hydrogen bubble visualization technique at Re=250.

One of the greatest challenges in the simulation of flow around oscillating cylinders is to deal with moving solid boundaries in complex geometries, especially for large amplitude body motions. The ALE (Arbitrary Lagrangian-Eulerian) method has been used in some studies. This means that the domain should be remeshed at every time interval; it is time-consuming. Hence, a numerical method should be proposed without aforementioned restrictions. In this paper, an Immersed Boundary Method  is adopted to treat the movable structure boundaries. Meanwhile, the fluid-structure interaction is treated as a multiphase problem. To do so, the CIP/CCUP (constraint interpolation profile/CIP combined, unified procedure)  is combined with the Cartesian grid system where the multiphase problem is solved in one set of equations. The CIP method has been applied for water-body interaction problems of free surface flow [21, 22]. In the present study, the CIP method is used for the study of flow over a streamwise oscillating circular cylinder with a stationary cylinder downstream.

To date, few studies considering the flow past two tandem circular cylinders with the upstream one forced to oscillate streamwise to the flow could be found. The objective of the present study is to extend the CIP-based method for this gap. In this paper, the Reynolds number remains at 100, oscillating amplitude is A/D=0.3, and the spacing between two tandem cylinders is L/D=2.0. By changing the oscillating frequency (fe) of the upstream cylinder (fe/f0=0.51.7; f0 is the natural shedding frequency of the two stationary tandem cylinders), the response of the wake is investigated. The rest of the paper is organized as follows. Section 2 gives the CIP-based numerical model. The flow solver and the coupling of fluid-body are briefly introduced. In Section 3, the current numerical model is first applied to flow over a fixed circular cylinder and two stationary tandem circular cylinders. Then, flow past an oscillating circular cylinder with a stationary cylinder downstream is computed in Section 4. Summary and conclusions are given in the final section.

2. Mathematical Modeling 2.1. Governing Equations

The basic equations governing the incompressible fluid flow are the mass conservation equation and the Navier-Stokes momentum equations written as(1)uixi=0,uit+ujuixj=-1ρpxi+1ρSijxj+fi,where Cartesian tensor notation is used (i= 1, 2) and t, uj, p, and xj are the time, velocities, hydrodynamic pressure, and spatial coordinates, respectively. fi represents momentum forcing components. Sij is the viscous term given by(2)Sij=μuixj+ujxi,where ρ and τ are the density and viscosity, respectively, appropriate for the phase that is occupying the particular spatial location at a given instant.

The numerical model considers the fluid-body interaction as a multiphase problem that includes fluid and structure. A fixed Cartesian grid that covers the whole computation domain is used. A volume function (or color function) ϕm (m=1 and 2 indicate fluid and solid, resp.) is defined to represent and track the interface. The total volume function for fluid and body is solved by using the following advection equation:(3)ϕ12t+uiϕ12xi=0.Here, ϕ12=ϕ1+ϕ2. The density and viscosity of the solid phase are assumed to be the same as those of a liquid phase to ensure stability. The volume function for solid body ϕ2 is determined by a Lagrangian method in which a rigid body is assumed. The volume function for fluid ϕ1 is then determined by ϕ1=1.0-ϕ2. After all volume functions have been calculated, the physical property λ, such as the density and viscosity, is calculated by the following formula: (4)λ=m=12ϕmλm.

2.2. Numerical Methods

Following Zhao and Hu (2012) , the governing equations are discretized using a high-order finite difference method on a Cartesian grid system. A staggered grid configuration is used to discretize the dependent variables. The governing equations (1) are solved using a fractional step scheme:(5)uit+ujuixj=0,ui-uiΔt=2μρSijxj+fi,xi1ρpxi=1Δtuixi,uin+1=ui-Δtρpxi.In the advection phase the spatial derivatives of the velocities are needed to apply the CIP method  for this phase. In the first nonadvection phase, terms of diffusion and external forces are solved by explicit method and renewed intermediate velocities are obtained. In the final phase, the velocities are updated coupling with pressure field.

To model the body motions, the fluid-structure interaction is coupled by using the Fractional Area Volume Obstacle Representation (FAVOR) method. The FAVOR was shown to be one of the most efficient methods to treat the immersed solid bodies. The effect of a moving solid body on the flow is considered by imposing the velocity field of the solid body into the flow at the solid edge. The following equation is introduced to update the local information of the fluid domain covered by the body:(6)U=ϕ2Ub+1-ϕ2u.Here, Ub is the local velocity of the solid body and u denotes the flow velocity obtained from the fluid flow solver. The solid phase ϕ2 is obtained using a virtual particle method . Ub, the local velocity of the solid body, is tracked by a Lagrange method. By integrating the pressure on the body surface, the hydrodynamic forces acting on the body are first calculated. With Newton’s Law, the body motion accelerations and velocities are calculated. More details can be found in previous references [21, 24].

3. Numerical Results 3.1. Flow Past a Stationary Circular Cylinder

Convergence study has first been performed for flow past a stationary circular cylinder at low Reynolds number (Re=100) as a reference for further investigation. The Reynolds number was defined on the diameter D of the circular cylinder and the dynamic viscosity coefficient μ. Three different nonuniform meshes are used with a minimum grid size of Δx=Δy=0.01D, 0.02D, and 0.04D, respectively. Results of the mean value of drag coefficient (Cd) and the Strouhal number (St) are shown in Table 1 and compared with previous experimental and numerical results. It can be seen that the computation converges fast and the results of mesh 1 and mesh 2 show excellent agreements with both the experimental and numerical data in existing literatures. Thus mesh 2 is adopted in the following calculations.

Comparison of mean drag coefficient (Cd) and Strouhal number (St) with those of other authors.

Data Min size of mesh Number of meshes C d S t
Mesh 1 0.01D 307800 1.38 0.166
Mesh 2 0.02D 92800 1.34 0.166
Mesh 3 0.04D 32175 1.28 0.166
Exp. (1959)  1.26
Exp. (1988)  0.164
Herfjord (1996)  1.36 0.168
Calhoun (2002)  1.33 0.175
Linnick and Fasel (2005)  1.34 0.166
Berthelsen and Faltinsen (2008)  1.38 0.169
Rajani et al. (2009)  1.34 0.157
3.2. Flow Past Two Stationary Cylinders in Tandem

Generally, the flow patterns of flow past two stationary tandem cylinders can be divided into two regimes by the critical spacing Lc. When L<Lc, there is no clear vortex shedding from the upstream cylinder (vortex suppression regime, VS). When L>Lc, the vortex sheds from both the upstream and downstream cylinders (vortex formation regime, VF).

In this section, flow past two stationary cylinders in tandem at Re=100 is computed with different spacing ratios L/D. The computational domain is scheduled in Figure 1. The inlet is located 10D upstream of the front cylinder, and the outlet is located 25.5D downstream of the front cylinder. The upper and lower boundaries are placed 10D away from the centerline of the domain. On the inlet, the streamwise velocity is set to make Re=100 for the fixed tandem circular cylinders and the transverse velocity is zero. On the outlet, a zero gradient boundary condition is specified. A free-slip boundary condition is imposed on the lateral boundaries and no-slip condition is applied on the surface of the two cylinders. Figure 2 shows the computational mesh and the detailed view of the mesh with L/D=5.0.

Sketch of the computational domain for flow past two tandem circular cylinders.

The computational mesh for flow past two tandem circular cylinders and local mesh near the cylinders with L/D=5.0.

Figure 3 illustrated the variation of coefficients at different spacing ratios L/D, ranging from 1.5 to 7.0. The computed St is shown in Figure 3(a), where the results of Li et al. (1991) , Sharman et al. (2005) , and Zheng et al. (2008)  are also presented for comparison. It shows a good agreement with existing results except for the results of Li et al. (2008)  a little larger. Figure 3(b) displays the mean value of the drag coefficient Cd and shows a good consistency with that of Sharman et al. (2005) . A sudden jump can be observed at L/D=3.54.0, which means that the critical spacing Lc is between 3.5 and 4.0. Figure 3(c) shows the same property as Figures 3(a) and 3(b). Figure 4 plots the vorticity contours at L/D=3.5 and 4.0. Positive contours are plotted in slid lines and negative contours in dashed lines. There is no clear vortex shedding from the upstream cylinder at L/D=3.5. After the critical spacing, the vortex shedding from both the cylinders can be observed at L/D=4.0.

Variation of coefficients with spacing ratio L/D for flow past two stationary circular cylinders: (a) Strouhal number St; (b) mean value of the drag coefficient Cd; (c) rms value of the lift force coefficient Clrms.

Vorticity contours of flow past two stationary tandem circular cylinders: (a) L/D=3.5; (b) L/D=4.0.

3.3. Flow Past a Streamwise Oscillating Circular Cylinder with a Stationary Cylinder Downstream in Tandem

In this section, flow past a streamwise oscillating circular cylinder with a stationary cylinder downstream in tandem with the spacing ratio L/D=2.0 at Re=100 is investigated. L/D=2.0 is in the regime of VS, where the wake interaction of upstream and downstream is much more complicated than that of VF because of the close distance. The computational domain and mesh are similar to that of flow past two stationary cylinders, as shown in Figures 1 and 2. The upstream cylinder oscillates streamwise with an imposed motion X(t)=Asin(2πfet). The oscillating amplitude A is set to be 0.3D, relatively large amplitude for vortex mode analysis . The natural shedding frequency of flow past two stationary cylinders f0 is 0.122 from the results of Section 2.2. The oscillating frequency fe/f0 varies from 0.5 to 1.7. In the present study, the computational time is set as 300 s and time interval is 0.002 s. The actual consuming time is no more than 83 min per condition using PC Intel i7-3770.

Figure 5 plots the spectrum of lift force for both the cylinders. It can be noticed that the upstream and downstream wakes have peaks at f0, fe, or their harmonics. The corresponding vorticity contours are presented in Figure 6. When fe/f0=0.5, the natural shedding frequency f0 plays a dominant role in the wake of both the cylinders. The vortex shedding is familiar with that of flow past two stationary tandem cylinders as shown in Figure 6(a). As the oscillating frequency fe/f0 increases to 0.7, the natural shedding frequency f0 dominates the upstream wake, while the forced frequency fe plays an important role in the downstream wake, and its harmonics becomes stronger. Corresponding vorticity contour displays a symmetric S pattern the same as flow past an oscillating cylinder . The vortex shed simultaneously from both the upper and lower sides of the cylinders. At fe/f0=0.91.1, the significant peaks appear at the excitation frequency, and 2fe exists in a much lower energy. The formation of vortex remains the same as that of fe/f0=0.7. There is a sudden change at fe/f0=1.3, 0.5fe emerges and plays the significant role in both the upstream and downstream wake, and its harmonics 1.5fe and 2.5fe also behave at fe/f0=1.31.7. The vortex sheds in an irregular 2P structure at fe/f0=1.3 and turns to a regular 2P at fe/f0=1.5 as displayed in Figures 6(f)-6(g). Figure 7 shows the vorticity contours in one period. The vorticity contours at 0T and T are nearly identical in the near-wake region at fe/f0=0.5 and 1.0, but it plots a good coincidence of T and 2T at fe/f0=1.5. It can be seen from Figures 8(c)-8(d) that a limit cycle is temporally repeated in the lift portrait plots for both the upstream and downstream cylinder at fe/f0=1.5, which shows a lock-in phenomenon. However, in Figures 8(a)-8(b), there are no single orbits but multiple loops, which indicate a quasi-periodic characteristic for both the upstream and downstream cylinder.

Spectrum analysis of the lift force for flow around upstream oscillating cylinder and downstream fixed cylinder: (a) upstream cylinder; (b) downstream cylinder.

Vortex contours of flow around upstream oscillating cylinder and downstream fixed cylinder: (a) fe/f0=0.5; (b) fe/f0=0.7; (c) fe/f0=0.9; (d) fe/f0=1.0; (e) fe/f0=1.1; (f) fe/f0=1.3; (g) fe/f0=1.5; (h) fe/f0=1.7.

Vortex contours in one period: (a) fe/f0=0.5; (b) fe/f0=1.0; (c) fe/f0=1.5.

Phase portrait plot at fe/f0=0.5 and 1.5: (a, c) upstream cylinder; (b, d) downstream cylinder.

4. Conclusions

A CIP-ZJU model is utilized to simulate the viscous flow past one and two circular cylinders in this paper. The model is established in a Cartesian coordinate system, using a CIP method as the base flow solver to discretise the N-S equations. The fluid-structure interaction is treated as a multiphase flow with fluid and solid phases solved simultaneously. An Immersed Boundary Method was used to deal with the fluid-body coupling.

Flow past a stationary circular cylinder at Re=100 is computed firstly to validate the method. Good agreements are obtained comparing with existing literatures. Then flow around two stationary tandem circular cylinders is computed at Re=100 of different spacing ratios L/D=1.57.0. The critical spacing Lc is between 3.5 and 4.0 as the results shown, which agrees well with previous data. Benchmark test cases showed that the CIP-based numerical discretization is accurate enough for low Reynolds number flows.

Flow around an upstream in-line oscillating cylinder and a stationary cylinder in tandem with spacing ratio L/D=2.0 at Re=100 is investigated. The oscillating amplitude is set as 0.3D, and the oscillating frequency fe/f0 is varied from 0.5 to 1.7. The spectrum analysis of lift force for both cylinders is discussed. When fe/f0=1.1, the natural shedding frequency f0 and the imposed frequency fe appear in the peaks. As the excitation frequency increases, 0.5fe plays a dominate role in both the upstream and downstream of the tandem system and 1.5fe and 2.5fe appear in a very low power. At the low oscillating frequency fe/f0=0.5, the shape of vortex is similar to that of flow past two stationary cylinders and the phase portrait of lift force shows a quasi-periodic characteristic. With the frequency increasing, S mode is observed at fe/f0=0.71.2 and the vorticity contours show a good accordance of 0T and T. A sudden change of the vortex formation is displayed at fe/f0=1.3, when vortex of irregular 2P mode shed from the cylinders. Regular 2P mode is formed at fe/f0=1.5, and the vorticity contours show a good identity every 2T. The portrait lift force displays a lock-in phenomenon with a single limit cycle.

Conflict of Interests

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

Acknowledgments

This work was supported by National Natural Science Foundation of China (51209184 and 51479175), Open Research Foundation of State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering (2013490211), Specialized Research Fund for the Doctoral Program of Higher Education (20120101120150), and Key Laboratory of Water-Sediment Sciences and Water Disaster Prevention of Hunan Province (Grant no. 2013SS03).

Zdravkovich M. M. Review of flow interference between two circular cylinders in various arrangements Journal of Fluids Engineering 1977 99 4 618 633 Li J. Chambarel A. Donneaud M. Martin R. Numerical study of laminar flow past one and two circular cylinders Computers and Fluids 1991 19 2 155 170 10.1016/0045-7930(91)90031-C ZBL0721.76049 2-s2.0-0025721618 Mittal S. Kumar V. Raghuvanshi A. Unsteady incompressible flows past two cylinders in tandem and staggered arrangements International Journal for Numerical Methods in Fluids 1997 25 11 1315 1344 10.1002/(SICI)1097-0363(19971215)25:11<1315::AID-FLD617>3.0.CO;2-P ZBL0909.76050 2-s2.0-0031413304 Meneghini J. R. Saltara F. Siqueira C. L. R. Ferrari J. A. Jr. Numerical simulation of flow interference between two circular cylinders in tandem and side-by-side arrangements Journal of Fluids and Structures 2001 15 2 327 350 10.1006/jfls.2000.0343 2-s2.0-0035527643 Jester W. Kallinderis Y. Numerical study of incompressible flow about fixed cylinder pairs Journal of Fluids and Structures 2003 17 4 561 577 10.1016/s0889-9746(02)00149-4 2-s2.0-0037870552 Sharman B. Lien F. S. Davidson L. Norberg C. Numerical predictions of low Reynolds number flows over two tandem circular cylinders International Journal for Numerical Methods in Fluids 2005 47 5 423 447 10.1002/fld.812 ZBL1085.76044 2-s2.0-13444250888 Mussa A. Asinari P. Luo L.-S. Lattice Boltzmann simulations of 2D laminar flows past two tandem cylinders Journal of Computational Physics 2009 228 4 983 999 10.1016/j.jcp.2008.10.010 MR2489939 2-s2.0-58149171905 Liu J. Guo H. Zhao J. Numerical simulation of flow around two tandem circular cylinders of equal diameter Periodical of Ocean University of China 2013 43 12 92 97 Li J. Sun J. Roux B. Numerical study of an oscillating cylinder in uniform flow and in the wake of an upstream cylinder Journal of Fluid Mechanics 1992 237 457 478 10.1017/s0022112092003495 2-s2.0-0026792607 Mahir N. Rockwell D. Vortex formation from a forced system of two cylinders. Part I: tandem arrangement Journal of Fluids and Structures 1996 10 5 473 489 10.1006/jfls.1996.0032 2-s2.0-0030480977 Papaioannou G. V. Yue D. K. P. Triantafyllou M. S. Karniadakis G. E. Evidence of holes in the arnold tongues of flow past two oscillating cylinders Physical Review Letters 2006 96 1 014501 10.1103/physrevlett.96.014501 2-s2.0-32644454741 Zou J.-F. Ren A.-L. Deng J. A. Forced system of two cylinders with various spacings China Ocean Engineering 2004 18 3 391 402 2-s2.0-8844228330 Prasanth T. K. Mittal S. Flow-induced oscillation of two circular cylinders in tandem arrangement at low Re Journal of Fluids and Structures 2009 25 6 1029 1048 10.1016/j.jfluidstructs.2009.04.001 2-s2.0-68849096866 Yang X. Zheng Z. C. Nonlinear spacing and frequency effects of an oscillating cylinder in the wake of a stationary cylinder Physics of Fluids 2010 22 4 10.1063/1.3372169 2-s2.0-77953337467 Huera-Huarte F. J. Bearman P. W. Vortex and wake-induced vibrations of a tandem arrangement of two flexible circular cylinders with near wake interference Journal of Fluids and Structures 2011 27 2 193 211 10.1016/j.jfluidstructs.2010.11.004 2-s2.0-79151482672 Bao S. Chen S. Liu Z. H. Li J. Wang H. F. Zheng C. G. Simulation of the flow around an upstream transversely oscillating cylinder and a stationary cylinder in tandem Physics of Fluids 2012 24 2 023603 10.1063/1.3683565 2-s2.0-84863235017 Tao S. Zhou L. Zong Z. Multi-block lattice Boltzmann method for flows around two tandem rotating circular cylinders Chinese Journal of Computational Physics 2013 30 2 159 168 2-s2.0-84877139797 Yang Y. C. Aydin T. B. Ekmekci A. Flow past tandem cylinders under forced vibration Journal of Fluids and Structures 2014 44 292 309 10.1016/j.jfluidstructs.2013.11.001 2-s2.0-84891634362 Peskin C. S. Flow patterns around heart valves: a numerical method Journal of Computational Physics 1972 10 2 252 271 10.1016/0021-9991(72)90065-4 2-s2.0-49649145203 Yabe T. Xiao F. Utsumi T. The constrained interpolation profile method for multiphase analysis Journal of Computational Physics 2004 169 2 556 593 10.1006/jcph.2000.6625 Hu C. H. Kashiwagi M. A CIP-based method for numerical simulations of violent free-surface flows Journal of Marine Science and Technology 2004 9 4 143 157 10.1007/s00773-004-0180-z 2-s2.0-12844258898 He G. A new adaptive Cartesian-grid CIP method for computation of violent free-surface flows Applied Ocean Research 2013 43 234 243 10.1016/j.apor.2013.10.004 2-s2.0-84887214119 Zhao X. Hu C. Numerical and experimental study on a 2-D floating body under extreme wave conditions Applied Ocean Research 2012 35 1 13 10.1016/j.apor.2012.01.001 2-s2.0-84863161281 Zhao X. Z. Ye Z. T. Fu Y. N. Cao F. F. A CIP-based numerical simulation of freak wave impact on a floating body Ocean Engineering 2014 87 50 63 10.1016/j.oceaneng.2014.05.009 Zheng T. Fei B. Yang J. Tang S. K. Numerical analysis of flow around two tandem circular cylinders Journal of Southwest Jiaotong University 2008 43 6 747 750, 798 2-s2.0-58149388942 Al-Mdallal Q. M. Lawrence K. P. Kocabiyik S. Forced streamwise oscillations of a circular cylinder: locked-on modes and resulting fluid forces Journal of Fluids and Structures 2007 23 5 681 701 10.1016/j.jfluidstructs.2006.11.001 2-s2.0-34248192834 Ongoren A. Rockwell D. Flow structure from an oscillating cylinder. Part 2. Mode competition in the near wake Journal of Fluid Mechanics 1988 191 225 245 10.1017/s0022112088001570 2-s2.0-0023738505 Tritton D. J. Experiments on the flow past a circular cylinder at low Reynolds numbers Journal of Fluid Mechanics 1959 6 4 547 567 10.1017/s0022112059000829 Williamson C. H. K. Defining a universal and continuous Strouhal-Reynolds number relationship for the laminar vortex shedding of a circular cylinder Physics of Fluids 1988 31 10 2742 2744 10.1063/1.866978 Herfjord K. A study of two-dimensional separated flow by a combination of the finite element method and Navier-Stokes equations [Dr.-Ing. thesis] 1996 Trondheim, Norway Norwegian Institute of Technology, Department of Marine Hydrodynamics Calhoun D. A Cartesian grid method for solving the two-dimensional streamfunction-vorticity equations in irregular regions Journal of Computational Physics 2002 176 2 231 275 10.1006/jcph.2001.6970 MR1894768 ZBL1130.76371 2-s2.0-0036501622 Linnick M. N. Fasel H. F. A high-order immersed interface method for simulating unsteady incompressible flows on irregular domains Journal of Computational Physics 2005 204 1 157 192 10.1016/j.jcp.2004.09.017 MR2121908 2-s2.0-14544280218 Berthelsen P. A. Faltinsen O. M. A local directional ghost cell approach for incompressible viscous flow problems with irregular boundaries Journal of Computational Physics 2008 227 9 4354 4397 10.1016/j.jcp.2007.12.022 MR2406541 ZBL05276027 2-s2.0-40849114423 Rajani B. N. Kandasamy A. Majumdar S. Numerical simulation of laminar flow past a circular cylinder Applied Mathematical Modelling 2009 33 3 1228 1247 10.1016/j.apm.2008.01.017 MR2478559 2-s2.0-55549119297