Characterization of Flow Separation around Inline Cylinders within Oscillatory Flow

)is paper has examined the effects of Reynolds number (Re), Keulegan–Carpenter number (KC), and gap ratio on flow separation around a cylinder array by PIVmethod in experiment.)e vortex shedding in such situation occurs each half period of the oscillatory motion from the observation. Nomatter howmany cylinders, KC is the key dominant parameter under low Re that has a great impact on the flow regime and flow motion in oscillatory flows. )ere is an influence area of vortex shedding around the cylinder. When two inline cylinders are in the flow, smaller gap ratio may extend the influence area of the vortex shedding. For s/d� 2, the vortex shedding that happens in each cylinder has a simultaneity and independence. For s/d� 1.5, the reduced gap ratio leads the upstream shed vortex to interact with downstream cylinder andmakes the influence area of vortex shedding around downstream cylinder extends further. For s/d� 1, the interference in the central area is significantly obvious and the vortex shedding is suppressed and even follows the cross gap flow.


Introduction
Flow separation around a single cylinder is well understood [1] that hydrodynamics around a single cylinder in both steady current and oscillatory flows from two aspects which are flow regimes and forces. For a cylinder in steady current, the flow regimes are dependent on the cylinder Reynolds number (Re). With increasing Re, vortices would be produced in the wake region of the cylinder.
Vortex shedding occurs at either side at a certain frequency in the wake region when Re is above 40. e turbulent boundary layer separation occurs at one side of cylinder with laminar at the other side in 3 * 10 5 < Re < 3.5 * 10 5 [2], which can cause significant nonzero mean lift on the cylinder. e two-dimensional vortex shedding occurs in the range of 40 < Re < 200 and the model of vortex shedding is laminar [3]. It will remain unchanged in the span direction. e critical Reynolds number for the onset of vortex shedding is determined to be 47 [4]. On the other hand, the same conclusion can be drawn that the vortex shedding turns into three-dimensional, and vortices are shed in cells in the spanwise direction in the range of 200 < Re < 300 [5].
Specifically, the onset of three-dimensionality is already observed at Re � 160 as confirmed numerically and experimentally [6]. e ratio for the friction drag to the total drag force is less than 4% with Re > 10 4 so that the friction drag is unimportant to be ignored in the total mean drag [7]. e aperiodic shedding regime may occur at certain KC and Re within switching between different modes in numerical ways [8].
e Hybrid Lagrangian/Eulerian discrete vortex method was used to simulate oscillatory flow around a cylinder [9], which has a good agreement with the experimental data. e interaction of vortex systems between a pair of sideby-side square cylinders in numerical, which demonstrated that different gap ratio of cylinders, plays an important role in the vortex shedding of side-by-side cylinders for a certain Keulegan-Carpenter number (KC) value [10]. e effects of different KC values have been studied by finite-element simulation, which is based on the oscillatory flows around four cylinders in square arrangement of 2 times gap ratio [11]. A direct-forcing immersed boundary method [12] is developed to simulate the hydrodynamic loading on circular cylinder array in oscillatory flow, from which the results are consistent with a previous study [11]. LBM (Lattice Boltzmann method) in numerical and PIV (particle image velocimetry) is used [13] in experiment to investigate the ow around six inline square cylinders in steady current, which pointed out the gap ratio is important for ow regime [14]. All the research has indicated the importance of gap ratio on the ow regime formation around a cylinder array. However, no one has observed ow separation around a cylinder array in the experiment. From the previous literatures, four leading factors should be vital to the ow separation around a cylinder array in oscillatory ows, which are Re, KC, gap ratio, and the number of cylinders erefore, an experiment was conducted to observe ow separation around a cylinder array, herein limited as two aligned cylinders, in oscillatory ows, and the analysis is presented in this paper.

Experimental Setup.
e experiment is conducted in a U-tube with a computer-controlled machine called SmartMotor [15] which can provide the oscillatory ows, as shown in Figure 1. Cylinders are set up inline in the center of bottom of the U-tube. e tracer particles need to be added into the ow before the experiment, which can visualize uid motions under a light source (below the U-tube). During the experiment process, a digital camera is set in front of the U-tube to capture the uid motions, and then data are output to the computer in Digi ow [16]. e uid can have di erent density by dissolving salt into the tap water, and the density can be adjusted for particles to make the uid fully lled with particles so that the best observation performance of uid motion can be achieved. e oscillation was often formed from the oscillation of the cylinder [1,17,18], which has three signi cant shortcomings: (1) the cylinder oscillation can lead to vibration, which may have an e ect on the ow stability and contribute deviation for the experiment results; (2) the movement of the cylinder means that it cannot cling to the wall, so the end e ect of the cylinder cannot be neglected; (3) the cylinder oscillation is more limited to the device than the motor motion used here which causes some inconvenience for oscillation simulation. Hence, the oscillatory ow is realized through the motor motion in the airtight U-tube in this experiment.
Due to vortices formed in the wake region, the ow can be greatly complex and unstable. Particle image velocimetry (PIV) is such a method which can capture instantaneous uid velocity with tracer particles in the whole ow elds. PIV [19] has evolved to be a dominant and standard method for velocimetry in experimental uid mechanics [20][21][22][23], especially in researching turbulent and complex ows. e digital camera can capture the tracer particle motion under a light source with ow to reconstruct the uid motion by interrogation in computer so that the ow velocity can be obtained.

Generation of Oscillatory Flow.
SmartMotor can simulate various oscillate movements easily in many applications by setting the target position, velocity, and acceleration. e oscillatory ow will take a cosinusoid curve for simplicity in this project, where A is the amplitude, T is the period, and t is the time. However, two key parameters of the oscillatory ow can have a great impact on the ow separation around cylinders, Re and KC, as follows: where, Re is Reynolds number, KC is Keulegan-Carpenter number, D is the diameter of the cylinder of 50 mm, U m is the maximum velocity, ] is the kinematic viscosity, and T is the period of the oscillatory ow. In the cosinusoid motion, longer T and smaller A make the oscillatory motion gentler. erefore, the value of Re and KC can be controlled through adjusting frequency and amplitude of motor motion. Meanwhile, the control variate method is applied to gure out how Re and KC a ect ow separation, respectively. U m can be obtained from U m ∼ y' 2π(A/T), so Re can be obtained from Re ∼ A/T and KC can be obtained from KC∼A. As certain amplitude A is decided, KC is de nitely the same for oscillatory ow with di erent periods T, and then the e ect of Re can be obtained. On the other hand, for the same A/T sample group, the e ect of KC can be observed by comparison of di erent groups.
Considering rationality of Re and KC values and the limitation of devices, values of A are set as 0.01, 0.02, 0.05, and 0.1 m and T are set as 3, 6, 9, and 12 s. As a result, the data sets are formed as shown in Tables 1 and 2. However, when SmartMotor performs the designed motions, there will be deviations between the proposed motions and simulation motions. ere is a phase di erence between the designed cosine curves and simulation curves with the time going. Additionally, the phase di erence gets more obvious with a longer period, which indicates that the  simulation frequency may be a little larger than the designed cosine movement. However, the amplitude is exactly the same for both curves. Hence, the Fourier transform is used to figure out the components of simulation motions for different frequencies. FFT (fast Fourier transform) is applied in MATLAB, which is a faster version of the discrete Fourier transform (DFT). e DFT takes a discrete signal in the time domain and transforms that signal into its discrete frequency domain representation. e result from FFT for y � 10 cos(2πt/6) with 20-period samples are demonstrated in Figure 2 e designed frequency against the simulated frequency curve is shown in Figure 2(b). As the order of magnitude of ] is significantly larger than T, the effect of frequency difference can be neglected. Meanwhile, KC is related to A. Hence, the generation of oscillatory flow is appropriate for such experiments.

Preliminary Experiment without Cylinders
e first stage of this experiment is to apply different oscillatory flows without cylinders to verify PIV measurement. At first, the oscillatory flow in the U-tube was considered as oscillatory Poiseuille flow. e analytical solution for oscillatory Poiseuille flow can be derived from the Navier-Stokes equation. e coordinate is set as follows: the x-axis is in the center of the channel parallel to the streamwise direction; the z-axis is perpendicular to the top and bottom walls (the spanwise direction); and the y-axis is perpendicular to the front and rear walls. In the experiment, just consider the two-dimensional flow separations so just consider the x-axis and z-axis. erefore, the Navier-Stokes equation can be expressed as As the velocity should be just horizontal one, only zu/zt is left. Meanwhile, the pressure is up to the different water pressure from the liquid level difference. As a result, is the amplitude of cosine motion, and l is the length of the U-tube and it is 2 m), so use them to replace the Navier-Stokes equation to get en solve this nonhomogeneous differential equation to get the velocity equation for oscillatory Poiseuille flow: where, a � 0.01 m, E � gA/l, ω � 2π/T (T is the period of the oscillatory movement), and ν � 1.5×10 −6 ·m 2 ·s −1 .
As there is a real part and imaginary part in this equation, MATLAB is applied to demonstrate how the velocity of Poiseuille flow changes with z and t. T (period) is set as 3 s and A is 0.1 m (to compare results from experiment). t will change from 0 to 3 s, so the velocity curve at z � 0 in an oscillatory period is demonstrated in Figure 3 z will change from −0.01 to 0.01 m, so the velocity curve at t � 2.3 is demonstrated in Figure 4. erefore, u-t curve has showed that the velocity has the same trend with the cosine oscillatory rate (same period). u-z curve presents the velocities of flows in different areas of U-tube. e flow has the same velocity and is uniform in most areas from bottom to top, but the velocity will become significantly small in the area near the bottom and top. However, the cylinder will be set up in the center of the tube, so the central area can provide the oscillatory flow with a uniform velocity. Similarly, PIV measurement also has the same results.
Oscillatory flows generated by motor with A � 1, 2, 5, and 10 cm and T � 3, 6, 9, and 12 s are simulated, and the PIV measurement is conducted for them, as shown in Figure 5.
e motion with A of 10 cm and T of 3 s is too quick to capture the tracer particles so it is given up. e cross section of the motor is different from the cross section of the flow, so the continuity principle is applied to transform the motor motion to the flow motion. e cross section of the motor is a circular section with diameter D of 0.185 m, and the flow section is a square section of 0.
So, the theoretical ow amplitude A f can be derived. Meanwhile, the uniform part of the velocity is cosine, so the theoretical maximum ow velocity also can be obtained. e whole uid velocity eld can be obtained from PIV measurement, and the average velocity from 80 mm to −80 mm is taken as the measured velocity. So the maximum measured velocity can be read from each curve, as shown in Table 3. Comparing the theoretical maximum ow velocity and the maximum measured velocity, there is a di erence between them, which may be caused by friction of the motor or the interrogation of the PIV method. However, the deviation is considered to be acceptable.

Results of Two Inline Cylinders in
Oscillatory Flow 4.1. Low KC. Two inline cylinders are set up along the center line with di erent distances between them. e gap ratio s/d (s is the distance between two cylinders' surfaces) is 1, 1.5, and 2. Similarly, KC is the key parameter for ow separation around two inline cylinders in low Re. When KC 0.8 and 1.7, the same ow pattern for every single cylinder occurs as just one cylinder set in [1]. e ow around each cylinder is almost attached to the cylinder, and ow separation cannot happen with low KC. So, the interference from the upstream cylinder cannot reach to the downstream one within such spacing s/d of 1, 1.5, and 2.

KC 4.2.
When KC reaches to 4.2, a pair of vortices will form in the wake region of each cylinder with respect to streamwise. e vortex pair is still attached to each cylinder as what happens in single cylinder. e ow pattern does not change with such KC and spacing. Figure 6(a) is the typical ow pattern in such condition. However, the spacing between cylinders can have an e ect on the velocity between them. According to the horizontal velocity contours for di erent gap ratios, the velocity in the gap is clearly decreased by 10 mm/s with the gap ratio reduced by 0.5, as shown in Figures 6(b)-6(d). Generally, low KC cannot change the ow pattern around the cylinder, but the smaller gap can decelerate the gap ow between cylinders. 8.4. In this situation, vortex shedding is supposed to occur around single cylinder, and the same is true of each one of two inline cylinders from observation, depicted in Figure 7. All experiments are conducted in Re 3515.1 and KC 8.4. Although the vortex shedding forms as before, the ow pattern will change in some respects with the di erent gap ratios. For s/d 2, the spacing between two cylinders is far enough to form an unperturbed vortex shedding around each cylinder. Additionally, the vortex shed from left cylinder moves the same distance with that shed from right  cylinder. As a result, the second downstream cylinder can be considered to have no e ect on upstream vortex shedding formation within s/d 2. Meanwhile, the gap ow depends on the ow motion in the wake region of each cylinder. Similarly, the vortex shedding occurs at each half period of the oscillatory ow motion as single cylinder. For s/d 1.5, the spacing is no longer long enough for each cylinder not to interfere with each other. e shed vortex from upstream cylinder reaches to the second downstream cylinder, so the vortex can interact with the downstream cylinder. It can be depicted in Figure 7(b) in which the downstream vortex shedding has propagated out of the camera range, but the downstream vortex just moves close to another cylinder. Furthermore, gap ow is a ected by the downstream cylinder. For s/d 1, such gap ratio makes a great di erence on the vortex shedding. As the blockage of the second cylinder, the ow will cross the gap, which can lead to the direction change of the vortex shedding. e vortex was washed towards the top boundary by the ow crossing the gap in Figure 7(c).

KC
However, a strange phenomenon is observed occasionally when s/d is 1.5, which is that the vortex ought to shed away from cylinder just to propagate a tiny distance, and the vortex shedding even disappears sometimes. e reason for this can be described as follows: rstly, the downstream cylinder blocks the path of vortex shedding and re ect the vortex back to upstream; secondly, the side ows put some pressure to vortex and make vortices attach to the   cylinder; and thirdly, the newly formed vortex may conflict with the old vortex formed in last half period, so the old vortex can push new one to attach the cylinder. e phenomenon will not happen within s/d � 2 because the gap is big enough while it also does not happen within s/d � 1 because the cross flow will make the vortex shed in the vertical direction.

Analysis of Velocity
Profile. e coordinate is set with a null point at the midpoint in the center line between two cylinders. One second is 30 frames. e fluid velocity field is provided by PIV measurement, as well. As above, the average velocity of the edge area is regarded as the average flow velocity without cylinder or the average flow velocity not affected by the cylinder (U 0 ) at each frame. All the velocity here is horizontal velocity. e case in such a condition with Re � 3515.1 and KC � 8.4 is focused on to figure out the influence of gap ratio on the flow separation. e variation of the normalized velocity U/U 0 along three lateral locations is investigated. e three locations are the left and right sides of two cylinders (which are the front and wake region) and the center lateral liner between these two cylinders. In order to compare the flow separation in two wake regions of these two cylinders, the distance from the two side locations to each cylinder should be equal to distance from center to each cylinder. So, x takes ±3d for s/d � 2, x takes ±2.5d for s/d � 1.5, and x takes ±2d for s/d � 1. Figures 8-10 depict the normalized velocity and velocity contour in a half period for s/d � 2, 1.5, and 1, respectively.
For s/d � 2, these two cylinders are considered not to have an impact on the flow separation around each other, which means that variation of velocity field in the wake region of each cylinder should show a similarity, as shown in Figure 8. x � 0 and x � −3d is the wake region as the flow comes from right side, while x � 0 and x � 3d is the wake region as the flow comes from left side. ese two wake regions are the same with respect to each cylinder. For both flow directions, the normalized velocity in wake region has a similar variation along y-axis. Moreover, the variation is the same with that of single cylinder. In addition, the difference of normalized velocity of two wake regions is relatively small, especially compared to that of s/d � 1 case. In other words, the same shape of these two wake-region-normalized velocity variations means that the vortex shedding happening in each cylinder has a simultaneity and independence. In the front area, the velocity at this lateral location is beyond the area of flow separation and vortex influence and keeps the same with the whole flow motion.
For s/d � 1.5, these two cylinders can affect each other, which can result in a chaotic state. However, in this situation, something strange occurred. From the first figure in Figure 9, the velocity near the cylinder in front area significantly decreased. e explanation is provided here: the distance is close enough to make the interaction of upstream vortex shedding and downstream cylinder possible. When the shed vortex is propagated to hit the downstream cylinder, it can make a contribution to cut off the vorticity source in downstream vortex shedding. So, the downstream vortex shedding is formed easily and quickly, which means that the vortex can move further. When the flow direction changes, the vortex will affect this area that converts from wake region to front region where the velocity will be decelerated. Certainly, the cylinder can block the flow, which also can cause the decelerated flow motion in front of the cylinder. Hence, the velocity decrease in front area can be observed. It can be seen in the velocity contour that a narrow band in the front area of the upstream cylinder is the slow-down area. On the other hand, the location along x � −2.5d is where vortex shedding is supposed to occur and vortex shedding is happening in the center area. However, the difference

Shock and Vibration
between them is considerably great compared with the case of s/d 2, which means that smaller gap leads to the interference on the vortex shedding in the central area. In the next half period, the velocity in front area also behaves as before but it generally keeps pace with the whole ow motion. Meanwhile, the velocities in both wake regions change in the same trend and there is still a di erence between them. In other words, the smaller gap ratio has an obvious e ect on the vortex shedding in the gap but the e ect is still weak. erefore, the gap ratio change can not only have an in uence on the center area between two cylinders but also make a great contribution to the ow regime of front and wake region. For s/d 1, velocity variation in the front area still exists as above described. Furthermore, the interference in the central area is signi cantly obvious. At rst, the velocity is reduced more in the downstream wake region. It can be derived that vortex shedding has already occurred in the downstream wake region while it is suppressed in the central area. Nevertheless, the velocity is reduced more in central  Figure 11. Figure 12 demonstrates the process of normalized velocity with different gap ratio s/d. According to the process of normalized velocity with time, the flow regime also shows a periodicity, and the period is still the half period of the whole flow motion. Point D and E have the same trend to some degree, so the effect of the gap ratio is relatively lower on the area of both sides. On the contrary, the velocity in the central area becomes completely different with a variation of gap ratio. As point B is located in the center of two cylinders, the velocity becomes more and more chaotic with two cylinders moving closer due to the interference of gap flow. Point A and C are located at the top and bottom area, and they have an opposite process compared with each other, which may be resulted from the vortex shedding formed at one side. However, with the gap smaller, the velocity on the top and bottom is no longer just periodically decreased or increased monotonously, which can be attributed to the gap cross flow formed. e interaction of gap flow and vortex shedding makes the motion more complex.

Conclusion
e particle image velocimetry (PIV) method is applied to measure the flow velocity field in different configurations to explore the effects of Re, KC, and gap ratio on flow separation around a cylinder array in oscillatory flows in this project. e first configuration is oscillatory flow without cylinder in which the whole flow in the tube acts as oscillatory Poiseuille flow. Oscillatory flow without cylinder is uniform in the most spanwise area except for a considerably small boundary layer near top and bottom wall where the velocity is significantly small. So, the whole flow without cylinder can be regarded as a uniform flow, and it is also a cosine oscillatory flow which has a same period with the oscillatory motor motion. e second configuration is two inline cylinders in oscillatory flow in which the effect of the gap ratio (s/d) is explored. Low KC cannot change the flow pattern around the cylinder, but smaller gap can decelerate the gap flow between cylinders. When KC � 8.4, vortex shedding will occur around each cylinder, and flow regime can be influenced significantly by the gap ratio. x =-2d x =2d x =0 In fact, the flow separation pattern under such a configuration is mainly determined by the cylinder-to-cylinder gap ratio according to the above investigation. More specifically, the gap flow between two incline cylinders would shape the track of vortex shedding.
(a) At s/d � 2, gap flow follows the general cosine oscillatory Poiseuille flow, and gap flow velocity keeps pace with velocity in wake region. Consequently, unperturbed independent vortex shedding around each cylinder simultaneously is produced. (b) At s/d � 1.5, the gap flow leads upstream vortex to attaching to the downstream cylinder so that upstream vortex interacts with downstream one and downstream vortex shedding goes further. Moreover, front velocity is significantly decreased due to the vortex shedding formed in last half period. (c) As s/d decreases to 1, the gap flow transforms from horizontally to vertically. Velocity pattern at top of the gap is opposite to that at bottom of the gap, while velocity in wake region exhibits resemblances with that of s/d � 2. And the cross gap flow would orientate the vortex shedding and push vortices generated from upstream cylinder away from the downstream cylinder.
Generally speaking, this paper provides a new way to analyze flow separation around a cylinder array with PIV measurement. Observation of vortex shedding is transferred to observation of velocity field variation, and it can provide a quantitative analysis to some extent.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.