Numerical Investigation of Flow around Two Tandem Identical Trapezoidal Cylinders

Wang, Yupu; Cheng, Wenming; Du, Run; Wang, Shubiao; Deng, Yong

doi:https://doi.org/10.1155/2020/3759834

Mathematical Problems in Engineering

On this page

Abstract Introduction Results and Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 3759834 | https://doi.org/10.1155/2020/3759834

Numerical Investigation of Flow around Two Tandem Identical Trapezoidal Cylinders

Yupu Wang,^1,2Wenming Cheng,^1,2Run Du ,^1,2Shubiao Wang,^1,2and Yong Deng^1,2

Academic Editor: Sébastien Poncet

Received31 Dec 2019

Revised28 Feb 2020

Accepted11 Mar 2020

Published27 Apr 2020

Abstract

The trapezoidal beam structure is ubiquitous in giant engineering equipment, while their aerodynamic characteristics have not been clearly understood. Numerical simulation method was adopted to investigate the flow around two tandem identical trapezoidal cylinders. The study was conducted using a Reynolds number of 2.2 × 10⁴, and with a spacing ratio varying from 0.5 to 10. The incompressible two-dimensional finite volume method was used for solving Reynolds-Averaged Navier–Stokes (RANS) equations with realizable model. The effects of cylinder geometry and spacing between the cylinders on aerodynamic characteristics, unsteady flow patterns, time-averaged flow characteristics, and flow instability was studied. The results show that the flow around the two tandem trapezoidal cylinders is highly dependent on the spacing ratio. The flow modes can be classified into: extended-body regime (Mode I, ≤ 1), reattachment regime (Mode II, 2 ≤ ≤ 3), and binary regime (Mode III, ≥ 4). We explored their respective flow characteristics and distinctions through the force/pressure coefficients, time-average streamwise velocity, and the flow field evolution.

1. Introduction

The flow over a bluff body is often accompanied by complicated phenomena such as flow separation, reattachment, and unsteady vortex shedding in the wake region. The repeated occurrence of large-scale vortex shedding repetition makes the structure susceptible to vibration, resulting in structural fatigue damage.

For a long time, scholars have focused on the flow problem of square, rectangular, or circular cross-section structures, but less attention has been paid to trapezoidal or other non-rectangular cross-sections. In the experiment, Luo et al. [1] investigated the influence of the afterbody shape and the angle of incidence on the flow structure around the prism by means of wind tunnel measurement and water tunnel flow field visualization. They found the main effect of the fluctuating aerodynamic forces differences is the proximity of the separated shear layers to their corresponding side faces. Venugopal et al. [2] carried out an experimental study on the influence of upstream disturbance and blockage on the flow around the trapezoid cylinder. Blockage refers to the ratio of the characteristic height of the bluff body to the flow field (D/H). The results show that the blockage ratio of 0.3 is found to be the best among all the blockages (0.14, 0.24, and 0.3) under both disturbed and undisturbed conditions.

In numerical simulation, Lee [3] studied the flow around non-rectangular sharp-edged prisms (trapezoidal cylinder) at the Reynolds number from 25 to 1000. The study shows that the structure of bluff bodies and Reynolds number have great influence on the characteristics of flow recirculation, flow separation, and flow interactions. Kahawita and Wang [4] studied the dynamic instability of a trapezoidal cylinder wake through two-dimensional numerical simulation. The spline method of fractional steps was used to study the Bénard-von Kármán hydrodynamic instability behind trapezoidal bluff body. Research found that the effect of trapezoid height seems to be the main effect on the value of Strouhal number, when compared with the effect of the smaller trapezoid bottom width. Cheng and Liu [5] studied the influence of rear-to-front width ratio (λ) from 0 to 1 (section varied from triangle to square) on the flow characteristics. The results show that λ = 0.5 makes the separation state of the shear layer and the wake vortex have the opposite trend. This means, when λ < 0.5, the vortex shedding frequency increases with the increase of λ, but when λ > 0.5, the frequency showed a negative correlation with the width. Study of Chung and Kang [6] showed that the Strouhal number of the trapezoidal cylinder depends on the Reynolds number and height ratio, which is the ratio of the height of the rear face to that of the front face (d/D). Dhiman and Hasan [7] analyzed the flow separation and flow state of the trapezoidal cylinder under low Reynolds number (Re = 1–150). Ding et al. [8] studied the flow induced vibration characteristics of the cylinder at high Reynolds number (Re = 3 × 10⁴∼1.1 × 10⁵) and analyzed the amplitude, frequency, and vortex shape of the cylinder. Miau et al. [9] studied the low-frequency behavior of vortex shedding embedded in the trapezoid cylinder with the Reynolds number ranging from 7.9 × 10³ to 3.24 × 10⁴.

There has been much research on the flow past two tandem cylinders. Zdravkovich [10–12] conducted a series of studies on the effects of interference between cylinders and divided them into three basic flow interference categories based on the position of the downstream cylinder. The three categories they proposed are as follows. (1) Proximity interference: the free shear layers separated from the upstream cylinder do not reattach on the downstream cylinder, but form directly the vortex behind the downstream cylinder. (2) Wake interference: the free shear layer separated from the upstream cylinder reattaches to the downstream cylinder. This mode creates a vortex only behind the downstream cylinder. (3) Binary vortices: two vortex channels are formed behind the two cylinders. The separated shear layers of the upstream cylinder are rolled up alternately in front of the downstream cylinder and form vortices. The vortex street behind the downstream cylinder is called binary flow because each vortex street is actually composed of two vortices formed, respectively, from the upstream cylinder and the downstream cylinder, and the binary vortex causes excessive fluctuations in the downstream cylinder. Meneghini et al. [13] studied numerically the vortex shedding and flow interference between two cylinders in tandem and in side-by-side arrangements. Jiang et al. [14] used the LBM (Lattice Boltzmann Method) to study the flow-induced vibration of one and two tandem circular cylinders and discussed the mechanism of vibration generation. Sumner [15] reviewed the understanding of the flow around two identical cylinders in a stable cross flow. The structure and behavior of near-wake and intermediate-wake of the two cylinders, as well as the general trends of the aerodynamic coefficient and Strouhal number, were summarized. The flow patterns can be roughly divided into (i) extended-body regime; (ii) reattachment regime; (iii) coshedding regime.

Sohankar and Etminan [16] numerically studied the flow characteristics and heat transfer on two equal square cylinders arranged in series at a low Reynolds number (Re = 1∼200). A floating stability analysis by Choi et al. [17] was performed to detect the occurrence of secondary instability in laminar flow around two identical square cylinders. The presence of adjacent cylinders can change the flow topology, resulting in changes in stability characteristics. Parametric studies were conducted with the gap between two square cylinders as the key parameter. Based on the gap, three different basic flow patterns were found and different modes of secondary instability were identified for each basic flow pattern. Yen et al. [18] used the PIV (Particle Image Velocimetry) and flow topology method to investigate the two squares cylinders in tandem experimentally. They divided flow characteristics into three modes: single mode, reattached mode, and binary mode at low Reynolds number (Re = 535). Some researchers carried out studies on the flow around two cylinders, and one of them (upstream/downstream) was forced to oscillate in the lateral direction. Xu and Zhou [19] used a hot-wire anemometer to study a series of square columns with Reynolds number Re = 800∼4.2 × 10⁴ and cylinder center-to-center distance L/d = 1–15. They proposed four distinct St-Re behavior modes. More et al. [20] used PIV, hot-wire anemometers, and flow visualization techniques to study the effect of the spacing between two series-connected cylinders and the vibration of the upstream cylinder on the flow characteristics. The experiments were performed at Reynolds number of 295. The study also captured three different flow patterns similar to previous studies.

To sum up, scholars have made some achievements in the study of the flow around a single trapezoidal cylinder, but the investigation on the flow around two tandem trapezoidal cylinders is still insufficient. In view of this, this paper numerically simulated the flow around two tandem trapezoidal cylinders at Re = 2.2 × 10⁴. We studied the aerodynamic characteristics, flow regime, and development pattern of the tail vortex at different spacing ratios so as to provide reference for the structural design of engineering equipment.

2. Numerical Details and the Governing Equations

2.1. Problem Description

The description of the problem is shown in Figure 1. Two identical trapezoidal cylinders in tandem were exposed to a constant free stream with uniform velocity U (=1.61 m/s). The flow was described in a Cartesian coordinate system (x and y), where the x-axis was aligned with the inlet flow direction and the y-axis was perpendicular to the x-axis. The origin was located at the center point of the upstream cylinder. The height D on the windward side of the trapezoidal cylinder was equal to 200 mm. All dimensions of the computational domain refer to the height of the windward side of trapezoidal cylinder (D). The width B of the trapezoidal cylinder was equal to D, and the height d-the rear face of the trapezoidal cylinder was equal to 0.5D.

In the y-direction, the vertical distance between the upper and lower boundaries of the domain was H (=20D), which determines the blockage rate of the flow field. The distance from the center of the upstream cylinder to the inlet was 10D, and the distance from the center of the downstream cylinder to the outlet was 20D. The spacing ratio was defined as = S/D, where S represents the distance between the rear face of the upstream cylinder and the front face of the downstream cylinder. was selected as 0.5, 1, 2, 3, 4, 6, and 10.

2.2. Numerical Calculation Governing Equations

Based on continuum theory and Newton’s second law, the Navier–Stokes equation of fluid motion can be derived. Due to the presence of turbulent flow, the Navier–Stokes equation is highly nonlinear and it is difficult to obtain a theoretical solution. In the RANS (Reynolds-Averaged Navier–Stokes) equations, the variables are decomposed into the mean and fluctuating components. For the incompressible flow in this paper, the RANS equations can be written as follows:where is strain rate tensor, .

Therefore, to make the equation closure, an effective method which employs Boussinesq hypothesis [21] to relate the Reynolds stresses to the mean velocity gradients is introduced:

Then, how to determine the turbulent viscosity becomes a key to the solution of problems in calculation. Shih et al. [22] proposed that the calculation of the turbulent viscosity of turbulent flow is directly related to the strain rate, and the turbulent viscosity can be expressed aswhere can be defined as the standard value of 0.09 in the equilibrium boundary layer, and 0.05 in a strong homogeneous shear flow.

Based on the Boussinesq hypothesis, is an isotropic scalar quantity. Although this hypothesis is not strictly validated. On many occasions, models based on the Boussinesq hypothesis perform very well. The realizable model has shown substantial improvements over the standard model. Because of the flow in this paper is thermostatic and incompressible, the transport equations for and in the realizable model arewhere

2.3. Calculation Area and Boundary Conditions

The calculation of flow area and boundary conditions are shown in Figure 1. The tandem trapezoidal cylinders in this study were placed in a two-dimensional, unsteady, viscous, and incompressible flow. The boundary of the flow field should be sufficiently far from the cylinders to achieve a lower blockage rate to meet the requirement of free-stream conditions, and the blockage of the simulation examples in this paper was 0.05. The calculation boundary conditions for the flow field were set as follows: Reynolds number: 2.2 × 10⁴ Uniform velocity inlet: = U, = 0 Pressure outlet: 0 Pa No-slip wall on the surface of the trapezoidal cylinder Top and bottom boundaries: symmetry Turbulence intensity of the flow: 5%

2.4. Definitions of Characteristic Quantities

The non-dimensional drag, lift, and pressure coefficient of the cylinder are defined as

The formula for the Reynolds and Strouhal number is

The time average and root mean square of the force coefficient fluctuation are calculated as follows (take as an example):where T is the fluctuation period.

2.5. Code Validation and Grid Independence Study

Figure 2 shows the non-uniform grid for the tandem trapezoid. In order to ensure the accuracy of the numerical simulation, the flow around a square and around a trapezoidal cylinder at Re = 2.2 × 10⁴ was calculated first. The calculation area was divided based on orthogonal grids, and the thickness of the first grid layer near the cylinder surface () was 0.0213D and 0.0355D, respectively. The Courant number (Co) is used to adjust the stability and convergence of the calculation. It is defined as follows:

According to the criterion that the Courant number is less than 1, the time steps were 1.33 × 10⁻³ s and 2.12 × 10⁻³ s, respectively. The calculation used the PISO algorithm, and the convergence value of the continuity equation was set to 10⁻⁵. Numerical results are shown in Table 1.

The drag coefficient for the calculation of the flow around the square-cylinder measured in this study was consistent with the experimental results of Lyn et al. [23] and Norberg [24]. At the same time, the calculation results matched well with the simulation results of Sohankar [25], Srinivas et al. [26], and Zhou et al. [27]. The Strouhal number of the square cylinder measured in the present study accord with the research findings of Atsushi and Okajima [28]. The time average drag coefficients and St of trapezoidal cylinders were in good agreement with the experimental values of Luo et al. [1], and the difference was less than 10%. It can be seen that the numerical simulation performed in this paper is accurate and reliable.

In both cases, grid thickness of the first layer near the cylinder surface was less than 3.6%, which was consistent with the literature value. Considering the computational cost of the calculation example, grid thickness of the first layer was determined as = 0.0355D in this paper for subsequent calculations.

3. Results and Discussion

3.1. Drag and Lift Coefficients

For the case of blunt body flow in engineering, the average drag coefficient () can be used to measure the wind load on bridges, buildings, or mechanical structures. The lift oscillation also influences safety in engineering. Figure 3 shows the variation of the average drag coefficient () and root-mean-square lift coefficient () for the two tandem trapezoidal cylinders at different spacing ratios .

(a)

(b)

In Figure 3(a), compared with the average drag coefficient of a single trapezoidal cylinder [1], the of the upstream cylinder was slightly reduced when = 0.5, 1. As the increased, the of the upstream trapezoidal cylinder remained basically unchanged and was similar to that of a single one. The of the downstream trapezoidal cylinder was smaller than that of the single one due to the influence by the wake of the upstream cylinder. The value of were negative when = 0.5, 1. When the spacing ratio increased to = 2, the changed from negative to positive and maintained at a relatively high value ( = 2, 3). As the spacing ratio continued to increase ( ≥ 4), the of the downstream trapezoidal cylinder first decreased and then gradually remained unchanged. Referring to study of Yen et al. [18] on tandem square cylinders, we divided the flow modes around two identical tandem trapezoidal cylinders into three categories: extended-body regime (Mode I, ≤ 1), reattachment regime (Mode II, 2 ≤ ≤ 3), and binary regime (Mode II ≥ 4). The distinction between mode I and mode II was a positive-negative conversion of the average drag coefficient of the downstream trapezoidal cylinder.

From Figure 3(b), the of the upstream trapezoidal cylinder maintained a low value in mode I, while it remained basically unchanged in mode II and mode III. The value of was equal to that of a single cylinder [1] in mode II and mode III. This shows that when ≥ 2, the upstream trapezoidal cylinder was less affected by the downstream trapezoidal cylinder. In mode I, the of the upstream and downstream trapezoidal cylinders were close. The of the downstream trapezoidal cylinder had a large rise from = 1 to = 2, and then reduced gradually with the increase of the spacing ratio ( ≥ 2). This meant that the effect of the upstream trapezoidal cylinder on the downstream trapezoidal cylinder gradually weakened as the spacing ratio increased. The of the downstream trapezoidal cylinder was larger than that of the upstream trapezoidal cylinder in mode II and reversed in mode III. Compared with the of a single trapezoidal cylinder, the results of mode I and mode III were smaller, while the results of mode II were larger.

3.2. Vortex Frequency and Strouhal Number

The lift force spectrogram of the two tandem trapezoidal cylinders with different spacing ratios are shown in Figure 4. The lift force spectrogram was related to the vortex frequency and Strouhal number of the flow field. For the upstream cylinder, the lift force oscillation was relatively stable at = 0.5, which means the interference between the two trapezoidal cylinders was weak. When 2 ≤ ≤ 4, the spectrogram changed drastically. As the spacing ratio increased, the interference between the two cylinders gradually disappeared, and the oscillation became stable again.

(a)

(b)

For the downstream cylinder, its lift force spectrogram had the same trend as the upstream cylinder. The difference was that when the pitch ratio was 2 ≤ ≤ 4, the vibration of the downstream cylinder was more severe. This phenomenon occurs because the vortex generated by the upstream cylinder directly acts on the downstream cylinder.

Vortex frequency and Strouhal number (St) for different spacing ratios are shown in Figure 5, as referring to the study by Wang et al. [29] and Kim et al. [30] on the tandem square cylinders. According to the research of the variation of St with on the two tandem square cylinders, there is a dip of St with a minimum value of ≈ 3. As the increases, the interaction between the two cylinders weakens, and the value of St is close to that of the single cylinder [23, 24].

When the spacing was small, the St of the two tandem trapezoidal cylinders reduced first and then enlarged with the increase of the space ratio, which was similar to the finding of the study on the tandem square cylinders. When ≥ 8, the St was equal to the finding, in the case of single trapezoidal cylinder by Luo et al. [1], studied. This shows that, in the case of large spacing, two cylinders exhibited an independent mode like a single cylinder (Mode III).

3.3. Aerodynamic and Flow Characteristics

Figure 6(a) shows the time-averaged surface pressure coefficient () around the upstream trapezoidal cylinder for different spacing ratios and the comparison with the single trapezoid. The on the windward side AB of the upstream trapezoidal cylinder was basically not affected by the spacing ratio, and the of the upper face BC and the rear face CD were greatly affected by the spacing ratio. The time-averaged surface pressure coefficient distribution trend was similar to the study results of the trapezoidal cylinder by Venugopal et al. [2] and Kahawita and Wang [4].

(a)

(b)

When the spacing ratio was relatively small ( ≤ 1), the pressure coefficient of the CD surface was affected by the reaction of the downstream cylinder, which resulted in a higher value. This result also explained why the of the upstream cylinder was slightly reduced in model I, as shown in Figure 3. With the increase of the spacing ratio, the circumferential distribution of the upstream trapezoidal cylinder approached to that of the single one. When = 0.5, 1, 3, 4, 6, and 10, the (pressure coefficient at point D) was −1.02, −1.22, −1.25, −1.50, −1.61, and −1.64, respectively. This meant that the effect of the downstream trapezoidal cylinder on the upstream one gradually decreased as the spacing ratio increased, as shown in Figure 3.

As Figure 6(b) shows, the circumferential pressure distribution of the downstream trapezoidal cylinder varied greatly under different spacing ratios. When the spacing ratio was very small ( ≤ 1), the on the windward side (EF) of downstream trapezoidal cylinder was negative, which lead to a negative of the downstream cylinder. The upper face (FG) of the downstream trapezoidal cylinder was little affected by the change of spacing ratio, and it was basically around −0.7. However, when = 3, the negative pressure of the FG face increased significantly, and the increased sharply around point F. This was due to the fact that the boundary layer separated from upstream cylinder reattached around point F. When = 3, 4, 6, and 10, the (pressure coefficient at point H) was −1.45, −1.52, −0.73, and −0.50, respectively, which meant the vortex in the wake region gradually moved away from the cylinder.

Figure 7(a) presents the time-averaged streamwise velocity distribution along the wake centerline, as shown in Figure 6(a), between the two trapezoidal cylinders. Velocity stagnation at points D and E was 0. For smaller spacing ratios ( = 0.5, 1), the time-average streamwise velocity was U-shaped, which meant that there was no vortex forming between the two cylinders. As the spacing ratio increases, the line shape changed to “S.” The shape of time-average streamwise velocity distribution is similar to the result in the case of tandem circular and square cylinders studied by Zhang et al. [31] There was a stationary point (X) where the velocity changed from negative to positive. The stationary point in the fluid can indicate the position of vortex. At the spacing ratios = 2, 3, the point X appeared near . As the spacing ratio increased, the point X moved forward and stabilized around . In summary, as the spacing ratio increased, the wake flow field of the upstream cylinder changed from no vortices, to distant vortices () and closer vortices ().

(a)

(b)

Figure 7(b) shows the streamwise velocity distribution at the rear of the downstream trapezoidal cylinder. When the spacing ratio = 0.5, 1, the time-average streamwise velocity of the downstream had a similar shape. As the spacing ratio increased, the downstream trapezoidal cylinder entered into mode II and mode III, and the wake speed decreased due to the influence of the upstream cylinder.

It can be seen from Figure 8 that the pulsating velocity () in the wake direction and vertical pulsating velocity () gradually decreased as the interval ratio increased. The number of wave crests was related to the generation of vortices. At the spacing ratio of = 4, the RMS velocity of the downstream appeared to be an asymmetric peak (y < 0). It can be inferred that there were asymmetric vortices. When = 10, the double peaks show that there was a double-row tail vortex, as shown in Figure 9.

(a)

(b)

(c)

(a)

(b)

(c)

3.4. Flow Field Evolution

Figure 9 gives the instantaneous flow visualization images obtained during the evolution of flow field around two tandem trapezoidal cylinders. The three spacing ratios in Figure 9 represented the flow field characteristics in mode ( = 0.5), mode ( = 3), and mode ( = 6), respectively. When = 0.5, the downstream trapezoidal cylinder was located in the trailing vortex formed by the boundary layer separation from the upstream cylinder (Figure 9(a)). The two trapezoidal cylinders behaved as a whole cylinder in the flow field (Figure 9(b)). Hence, model was called an extended-body regime. In this mode, the wake was a single-row vortex (Figure 9(c)). The phenomenon was similar to the single trapezoidal cylinder.

When the spacing ratio increases to = 3, the boundary layer separated from the upstream trapezoidal cylinder reattach on the front or upper face of the downstream trapezoidal cylinder (Figure 9(a)). The tail of the upstream trapezoidal cylinder formed a slightly backward vortex (compare to the tail image of = 6 in Figure 9(a)). The boundary layer flow, attaching to the downstream trapezoidal cylinder, made the wake disorder (Figure 9(b)). This also explained the reason for the high of the downstream trapezoidal cylinder in mode in Figure 3. In this mode, a double-row vortex street appeared in the wake of the downstream trapezoidal cylinder (Figure 9(c)).

When = 6, the wake vortex of the upstream trapezoidal cylinder was formed. Under these circumstances, the flow around the upstream trapezoidal cylinder was similar to that around single one, and it was hardly affected by the downstream trapezoidal cylinder (Figure 9(a)). The downstream trapezoidal cylinder was still influenced by the upstream cylinder, appearing to be a double-row tail vortex.

3.5. Characteristics and Distinction of the Three Modes

Great attention has been paid to the flow around two tandem cylinders in the literature. The studies of the flow around two tandem circular cylinders at different Reynolds numbers and with different spacing ratios were conducted by Zdravkovich [10–12], Sumner [15], and Xu and Zhou [19]. Zdravkovich [10] classified the flow around two circular cylinders into three regimes based on the spacing ratio. Regime 1, the free shear layers separating from the upstream cylinder overshoot the downstream cylinder, forming one vortex street. Regime 2, the free shear layers separated from the upstream cylinder is reattached on the downstream cylinder. A single vortex street is formed behind the downstream cylinder. Regime 3, the free shear layers roll up alternately to form vortices between the cylinders. Vortices generate at both cylinders. Sumner [15] summarized the research studies on the flow around two tandem circular cylinders and divided the flow patterns into (i) extended-body regime; (ii) reattachment regime; (iii) co-shedding regime. The study of Xu and Zhou [19] found that the flow patterns around the two tandem circular cylinders is not only related to the spacing ratio but also to the Reynolds number. The change of the Reynolds number will lead to the change of the critical spacing ratio under different modes.

Similar phenomena were observed in the studies of the flow around two square cylinders in tandem by Sohankar and Etminan [16], Choi et al. [17], More et al. [20], and Yen et al. [18]. The flow pattern changes with the spacing ratio. Yen et al. [18] also divided the flow characteristics into three modes: single mode, reattached mode, and binary mode. The aerodynamic parameters of cylinders and critical spacing ratio of each mode are different from those of the circular cylinders. The most notable to the flow around two circular cylinders is that the wake of the downstream square cylinder forms a double vortex street in reattached and binary modes.

The flow around two tandem trapezoidal cylinders was similar to that around square cylinders. We divided the flow modes around two identical tandem cylinders into three categories: extended-body regime (Mode , ≤ 1), reattachment regime (Mode , 2 ≤ ≤ 3), and binary regime (Mode , ≥ 4), as above. The characteristics and distinction of the three modes are shown in Table 2.

In mode , the of the windward side (EF) of downstream trapezoidal cylinder was negative. The of the upstream one was slightly smaller than that of the single trapezoidal cylinder, and the of the downstream trapezoidal cylinder was negative. As the spacing ratio increased, there was a critical spacing ratio that changed the of the downstream trapezoidal cylinder from negative to positive (also the of the EF surface). This transition could be a demarcation between mode and mode . The time-average streamwise velocity between the two trapezoidal cylinders was U-shaped. There was no vortex forming between two cylinders; at the same time, there was only one vortex in the wake of the downstream trapezoidal cylinder. In mode , the two trapezoidal cylinders behave like a single bluff-body.

In mode , the of the downstream trapezoidal cylinder was larger than that of the upstream trapezoidal cylinder. The of the EF surface reduces significantly. The time-average streamwise presents a horizontal S-shape (Figure 7 (a)) and the stationary point was around . The boundary layer separated from the upstream cylinder was reattached to the downstream cylinder.

In mode , the distribution of the aerodynamic coefficients (, , and ) were similar to that of the single one. The time-average streamwise velocity between the two trapezoidal cylinders was S-shaped, where the stationary point was around . In this mode, the separated shear layers of the upstream cylinder are rolled up alternately in front of the downstream cylinder and form vortices. The vortex street behind the downstream cylinder is called binary flow because each vortex street is actually composed of two vortices formed behind the upstream cylinder and the downstream cylinder.

To sum up, the flow pattern around two tandem cylinders could divide into three modes according to the differences of the aerodynamic characteristics and flow field structure under different spacing ratios. Under different cross sections (circular, square, and trapezoidal), the aerodynamic parameters and critical spacing ratio of each mode were different.

4. Conclusions

In this paper, the numerical simulation of the flow was carried out around two identical trapezoidal cylinders in a tandem arrangement at Re = 2.2 × 10⁴. This study investigated the interference of flow past two tandem trapezoidal cylinders at different spacing ratios of = 0.5–10. The aerodynamic characteristics, wake flow, and vortex shedding patterns of the two cylinders were considered. The spacing ratio had a significant effect on the flow field around two tandem identical trapezoidal cylinders. According to the differences of the aerodynamic characteristics and flow field structure under different spacing ratios, we divided the flow modes into three categories: extended-body regime (Mode , ≤ 1), reattachment regime (Mode , 2 ≤ ≤ 3), and binary regime (Mode , ≥ 4).

At the spacing ratio of = 0.5, 1, the downstream cylinder was located in the trailing vortex formed by the boundary layer separated from the upstream cylinder, which led to a relatively low amplitude of the lift oscillations. The average drag coefficient of the downstream trapezoidal cylinder was negative, and the sum of drag coefficient of the two cylinders was minimum. The separated boundary layer of the upstream cylinder only formed a single-row vortex at the tail of the downstream cylinder. There was no vortex between the two cylinders. In mode , the downstream trapezoidal cylinder acts as an extended-body of the upstream one.

When the spacing ratio was 2 ≤ ≤ 3, the flow field around the structure was in reattachment regime. The downstream cylinder was affected by the reattachment of the wake flow of the upstream cylinder. There were many small-scale vortices around the downstream cylinder, and the lift amplitude spectrum became messy and multiple low peaks appear. The maximum value of the downstream trapezoidal cylinder appeared in this mode. In mode , the boundary layer separated from the upstream cylinder and reattached on the downstream cylinder. The large negative pressure on the upper surface of the downstream cylinder was one of the important characteristics that distinguishes mode from the other two modes.

As the spacing ratio increases to ≥ 4, the influence of the upstream trapezoidal cylinder on the downstream cylinder was further weakened. The aerodynamic parameters of the upstream cylinder were similar to that of the flow around a single one. Compared to mode , the stationary point between the two trapezoidal cylinders changed from the distant vortices () to the closer vortices (). In this mode, the separated shear layers of the upstream cylinder were rolled up alternately in front of the downstream cylinder and form vortices. Because each vortex street was actually composed of two vortices formed behind the upstream cylinder and the downstream cylinder, this mode was called binary regime.

Nomenclature

U:	Upstream velocity (m/s)
D:	Height of the windward side of the trapezoidal cylinder
B:	Width of the trapezoidal cylinder
d:	Height of the rear end of the trapezoidal cylinder
S:	Spacing length
:	Spacing ratio
k:	Kinetic energy of turbulence
:	Pressure and average pressure
:	Drag force
:	Lift force
:	Drag coefficient
:	Lift coefficient
:	Pressure coefficient
:	Mean back pressure coefficient
L _r:	Swing length
f:	Frequency of vortex shedding
Re:	Reynolds number
St:	Strouhal number
:	Fluid velocity
:	Average and fluctuation component of velocity
x, y:	Cartesian coordinates
:	Fluid kinematic viscosity

Greek letters

:	Fluid density
:	Rate of dissipation of turbulent kinetic energy
:	Fluid dynamic viscosity

Subscripts

avg:	Time averaged
rms:	Root mean square
i, j:	Component in x and y directions.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (NSFC), project no. 51675450.

References

S. C. Luo, M. G. Yazdani, Y. T. Chew, and T. S. Lee, “Effects of incidence and afterbody shape on flow past bluff cylinders,” Journal of Wind Engineering & Industrial Aerodynamics, vol. 53, no. 7, pp. 375–399, 1994.
View at: Publisher Site | Google Scholar
A. Venugopal, A. Agrawal, and S. V. Prabhu, “Influence of blockage and upstream disturbances on the performance of a vortex flowmeter with a trapezoidal bluff body,” Measurement, vol. 43, no. 4, pp. 603–616, 2010.
View at: Publisher Site | Google Scholar
T. S. Lee, “Early stages of an impulsively started unsteady flow past non-rectangular prisms,” Computers & Fluids, vol. 27, no. 4, pp. 435–453, 1998.
View at: Publisher Site | Google Scholar
R. Kahawita and P. Wang, “Numerical simulation of the wake flow behind trapezoidal bluff bodies,” Computers & Fluids, vol. 31, no. 1, pp. 99–112, 2002.
View at: Publisher Site | Google Scholar
M. Cheng and G. R. Liu, “Effects of afterbody shape on flow around prismatic cylinders,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 84, no. 2, pp. 181–196, 2000.
View at: Publisher Site | Google Scholar
Y. J. Chung and S.-H. Kang, “Laminar vortex shedding from a trapezoidal cylinder with different height ratios,” Physics of Fluids, vol. 12, no. 5, pp. 1251–1254, 2000.
View at: Publisher Site | Google Scholar
A. Dhiman and M. Hasan, “Flow and heat transfer over a trapezoidal cylinder: steady and unsteady regimes,” Asia-Pacific Journal of Chemical Engineering, vol. 8, no. 3, pp. 433–446, 2013.
View at: Publisher Site | Google Scholar
L. Ding, L. Zhang, and D. Y. Jiang, “Flow induced motion of bluff bodies with different cross sections in flow field with high Reynolds number,” Journal of Vibration and Shock, vol. 34, no. 12, pp. 176–181, 2015.
View at: Google Scholar
J. J. Miau, J. T. Wang, J. H. Chou, and C. Y. Wei, “Low-frequency fluctuations in the near-wake region of a trapezoidal cylinder with low aspect ratio,” Journal of Fluids and Structures, vol. 17, no. 5, pp. 701–715, 2003.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
M. M. Zdravkovich, “Review-review of flow interference between two circular cylinders in various arrangements,” Journal of Fluids Engineering, vol. 99, no. 4, pp. 618–633, 1977.
View at: Publisher Site | Google Scholar
M. M. Zdravkovich, “Smoke observations of wakes of tandem cylinders at low Reynolds numbers,” Aeronautical Journal, vol. 76, no. 734, pp. 108–114, 1972.
View at: Google Scholar
J. R. Meneghini, F. Saltara, C. L. R. Siqueira, and J. A. Ferrari, “Numerical simulation of flow interference between two circular cylinders in tandem and side-by-side arrangements,” Journal of Fluids and Structures, vol. 15, no. 2, pp. 327–350, 2001.
View at: Publisher Site | Google Scholar
R. J. Jiang, J. Z. Lin, and X. K. Ku, “Flow-induced vibrations of two tandem circular cylinders in a parallel-wall channel,” Physics of Fluids, vol. 26, no. 10, p. 413, 2014.
View at: Publisher Site | Google Scholar
D. Sumner, “Two circular cylinders in cross-flow: a review,” Journal of Fluids and Structures, vol. 26, no. 6, pp. 849–899, 2010.
View at: Publisher Site | Google Scholar
A. Sohankar and A. Etminan, “Forced-convection heat transfer from tandem square cylinders in cross flow at low Reynolds numbers,” International Journal for Numerical Methods in Fluids, vol. 60, no. 7, pp. 733–751, 2009.
View at: Publisher Site | Google Scholar
C. B. Choi, Y. J. Jang, and K. S. Yang, “Secondary instability in the near-wake past two tandem square cylinders,” Physics of Fluids, vol. 24, no. 2, pp. 1–22, 2012.
View at: Publisher Site | Google Scholar
S. C. Yen, K. C. San, and T. H. Chuang, “Interactions of tandem square cylinders at low Reynolds numbers,” Experimental Thermal and Fluid Science, vol. 32, no. 4, pp. 927–938, 2008.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
B. S. More, S. Dutta, M. K. Chauhan, and B. K. Gandhi, “Experimental investigation of flow field behind two tandem square cylinders with oscillating upstream cylinder,” Experimental Thermal and Fluid Science, vol. 68, pp. 339–358, 2015.
View at: Publisher Site | Google Scholar
J. O. Hinze, Turbulence, McGraw-Hill Publishing Co., New York, USA, 1975.
T.-H. Shih, W. W. Liou, A. Shabbir, Z. Yang, and J. Zhu, “A new k–ε eddy viscosity model for high Reynolds number turbulent flows,” Computers & Fluids, vol. 24, no. 3, pp. 227–238, 1995.
View at: Publisher Site | Google Scholar
D. A. Lyn, S. Einav, W. Rodi, and J.-H. Park, “A laser-Doppler velocimetry study of ensemble-averaged characteristics of the turbulent near wake of a square cylinder,” Journal of Fluid Mechanics, vol. 304, pp. 285–319, 1995.
View at: Publisher Site | Google Scholar
C. Norberg, “Flow around rectangular cylinders: pressure forces and wake frequencies,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 49, no. 1–3, pp. 187–196, 1993.
View at: Publisher Site | Google Scholar
A. Sohankar, “A numerical investigation of the flow over a pair of identical square cylinders in a tandem arrangement,” International Journal for Numerical Methods in Fluids, vol. 70, no. 10, pp. 1244–1257, 2012.
View at: Publisher Site | Google Scholar
Y. Srinivas, G. Biswas, A. S. Parihar, and R. Ranjan, “Large-eddy simulation of high Reynolds number turbulent flow past a square cylinder,” Journal of Engineering Mechanics, vol. 132, no. 3, pp. 327–335, 2006.
View at: Publisher Site | Google Scholar
Q. Zhou, H. L. Liao, and S. Y. Cao, “Numerical study of flow characteristic around square cylinder at high Reynolds number,” Journal of Southwest Jiaotong University, vol. 53, no. 3, pp. 533–539, 2018.
View at: Google Scholar
Atsushi and Okajima, “Strouhal numbers of rectangular cylinders,” Journal of Fluid Mechanics, vol. 123, pp. 379–398, 1982.
View at: Publisher Site | Google Scholar
X. K. Wang, Z. Hao, J. X. Zhang, and S. K. Tan, “Flow around two tandem square cylinders near a plane wall,” Experiments in Fluids, vol. 55, no. 10, pp. 1–14, 2014.
View at: Publisher Site | Google Scholar
M. K. Kim, D. K. Kim, S. H. Yoon, and D. H. Lee, “Measurements of the flow fields around two square cylinders in a tandem arrangement,” Journal of Mechanical Science and Technology, vol. 22, no. 2, pp. 397–407, 2008.
View at: Publisher Site | Google Scholar
W. Zhang, H.-S. Dou, Z. Zhu, and Y. Li, “Unsteady characteristics of low-Re flow past two tandem cylinders,” Theoretical and Computational Fluid Dynamics, vol. 32, no. 4, pp. 475–493, 2018.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Yupu Wang 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.

PDF Download Citation

Download other formats

Order printed copies