Unsteady RANS Modeling of Flow around Two-Dimensional Rectangular Cylinders with Different Side Ratios at Reynolds Number 6.85× 10

In practical engineering, the Reynolds number (Re) of box girder structure is usually very high (Re≥ 10), while most investigations of the flow around bluff bodies are concentrated on relatively lower Reynolds numbers (i.e., Re� 10–10). )is paper presented a numerical study of the unsteady flow around two-dimensional rectangular cylinders under a Reynolds number of 6.85×10 with different side ratios (R� b/h, width to height) ranging from 0.1 to 4.0. )ree unsteady Reynolds-averaged NavierStokes (RANS) two-equation k-ε turbulence models (standard, RNG, and realizable) were adopted in the study. )e realizable k-ε model was chosen because it was found to perform the best among three models in the main aerodynamic integral parameters. According to the distinctions of aerodynamic characteristics with different side ratios, three regimes were divided and discussed in detail. )e distribution of surface pressure over cylinders, the wake parameters, and vorticity contours of the rectangular cylinders with different side ratios were discussed.


Introduction
Bluff structures, which are the most common configurations in numerous practical applications such as tall buildings, bridges, chimneys, and fences, are of significance in both engineering and science when concerning the flow around them. Due to complicated fluid dynamic phenomena, such as separation, reattachment, and vortex shedding, the flow around bluff body has obtained much attractiveness. e fundamental aerodynamic issues of bluff structures like square cylinders have been examined extensively in terms of numerical and experimental studies. Lyn et al. [1] studied the turbulent near-wake flow around a square cylinder at Re � 2.14 × 10 4 with two-component laser Doppler measurements. Shun et al. [2] discussed the flow structure, wake flow characteristics, and drag coefficient of a square cylinder at various Re (ranging from 400 to 3.6 × 10 4 ) and incident angles (adjusting from 0°to 45°) in an open-loop wind tunnel with the particle image velocimetry (PIV) technology. Rodi et al. [3] reviewed the Large Eddy Simulation (LES) and RANS methods for flow over a square cylinder and discussed the factors on the results such as the flow, grids, boundary conditions, and domain size. Sohankar [4] investigated the effects of the Reynolds number on results of the flow over a square cylinder using LES approaches and found that the effect is not much at higher Re, i.e., Re > 2 × 10 4 .
However, bluff bodies with rectangular cross section are employed more universally. Researches on rectangular cylinders' aerodynamic features and flow structures in detail are significant and meaningful for engineering practice. Many researchers have studied the flow around rectangular cylinders in experiments. Nakaguchi et al. [5] studied rectangular cylinders with different side ratios (R � b/h, where b is the chord length along the direction of the flow and h is the depth of the section) experimentally in a wind tunnel and found a peak in drag coefficient (C d ) near R � 0.6 and a discontinuity in Strouhal number (St) at R � 2.8. Okajima et al. [6,7] studied the influence of R on the lift and drag forces, base pressure, and St in the range of Re between 70 and 2 × 10 4 experimentally and discovered a sudden change in St at R � 2.8 and 6.0. Nakamura and Ohya [8] investigated the effect of scale and the intensity of turbulence on the wake of rectangles in a wind tunnel. ey concluded that the small-scale turbulence could raise the growth rate of the separated shear layers and the large-scale turbulence could decline vortex shedding. Both scale and intensity of turbulence affect the mean base pressure (−C pb ) apparently. Ohya et al. [9], Norberg [10], and Nakagawa et al. [11] also obtained similar results at different side ratios and Reynolds numbers.
Besides the experimental studies, numerical studies had also been carried out. Okajima et al. [12,13] investigated flow characteristics of rectangular cylinders with R � 0.6-8 by a finite difference method (FDM) at low Re and by a discrete vortex method (DVM) at high Re. eir numerical results show a good agreement with their experimental ones [6,7]. Taylor et al. [14] studied the two-dimensional (2D) flow around cylinders with R � 0.25-3 at Re � 2 × 10 4 using DVM.
eir numerical methods were verified well with experimental results. Shimada et al. [15] analyzed the aerodynamic characteristics of cylinders with R � 0.6-8 at Re � 2.2 × 10 4 numerically using a modified k-ε model. eir method presents better results than ordinary 2D analysis in vortex shedding, drag coefficients, and distributions of mean surface pressure. ey reported a discontinuity in Strouhal number at both R � 2.8 and R � 6. ey also explained that the reason of the underestimation in pressure and force fluctuations is that the stochastic components are not contained in RANS method. Tian et al. [16] demonstrated some inaccurate results when side ratio is lower (R < 0.6) using k-ω Shear Stress Transport (SST) method. Aerodynamics features of rectangular cylinders (R � 0.4-4) at Re � 10 5 were simulated by Sohankar [17] using two LES approaches. e global integral results and the mean and instantaneous-flow structures were presented in detail in their study. Islam et al. [18] discussed these at Re � 100-250 using incompressible lattice Boltzmann method (ILBM).
Although there have been many researches about the flow over rectangular cylinders before, they are mainly on relatively lower Reynolds numbers, e.g., experimental studies of Nakaguchi et al. [5] (Re � 6 × 10 4 ), Bearman and Trueman [19] (Re � 2.2 × 10 4 ), Mizota et al. [20] (Re � 2.1 × 10 5 ), and Norberg [10] (Re � 1.3 × 10 4 ) and numerical studies of Shimada and Ishihara [15] (Re � 2.2 × 10 4 ), Sohankar [17] (Re � 10 5 ), and Tian et al. [16] (Re � 2.14 × 10 4 ). Most of them were performed at Re ≤ 10 5 . However, the section size of the structure in practical applications (i.e., bridges, buildings, and engineering machineries) is usually larger, so the Reynolds numbers trend to be higher. For example, Re of the beam or leg component of the gantry or port cranes, which can be seen as the typical rectangular column structure, would be more than 5 × 10 5 when wind speed (U) is 10 m/s and when the characteristic height (h) of them is designed to be more than 1 m (e.g., 1.998 m in [21] and 1.82 m in [22]). At some extreme cases (h � 5-10 m, U ≥ 30 m/s), the Reynolds number Re would be up to 10 7 . As a result, it is necessary and meaningful to study the flow around rectangles with various R when Re is much higher. ey are more common in engineering and could give some guidance for industrial design and manufacture.
In this paper, our main purpose is to study the flow around rectangular cylinders with different side ratios when Reynolds number is much higher.
e main global parameters and mean and instantaneous-flow structures would be discussed in this work. e rest of this paper is as follows. Section 2 gives a short presentation of numerical models, governing equations, boundary conditions, and grids. In Section 3, the results and discussion are presented, including the difference in different models, the effect of side ratios for global parameters, analysis of wake, and the vorticity contours. Some concluding remarks are given in Section 4.

Numerical Methods
Turbulence is a time-dependent three-dimensional process whose simulation by means of computational techniques is particularly challenging. Several turbulence models, comprising modified versions, have been proposed.
ey are different in their physical foundations, degree of complexity, and computer power demands. And the most commonly used models at present are the Reynolds-averaged Navier-Stokes (RANS) and Large Eddy Simulation (LES).
Rodi et al. [3], Murakami [23], Hanjalic [24], Menter and Florian [25], and Blocken [26] have done lots of investigations about the comparisons and discussions of these two methods. Although there is some difference between them, their basic views are relatively consistent. e LES approach applied in a 3D domain might be more accurate and reliable, but it is also the most computationally burdensome because of its high cost in wall boundary layers. It could not be widely used in industrial applications especially when Reynolds number is higher. And, as for the RANS models, due to the use of wall functions, the cost is less, and the accuracy of calculation is relatively acceptable; hence, it is still the most widely used turbulence model in engineering. Concerning a 3D flow domain, RANS models cannot resolve any other turbulence details except for single large-scale structures [25]; their results are practically similar to 2D simulations and show weak three dimensionalities [27,28]. at is, 3D simulations cost more but do not receive significantly better results than 2D ones for global parameters; thus 2D RANS is a better choice.
Taking into account the fact that the primary goal of this piece of research is the CFD-based computation of flow around rectangles with various side ratios in much higher Re, 3D simulations have been discarded, and the computer cost linked with LES approaches could not be afforded by the authors; the 2D RANS k-ε turbulence models based on the Boussinesq assumptions [29] are chosen.

2
Mathematical Problems in Engineering

RANS k-ε Turbulence Model.
e Reynolds-averaged incompressible Navier-Stokes equations [29] are expressed as where u i is the velocity component in x i coordinate directions and, for the present 2D simulation, i, j � 1, 2, x 1 and x 2 are the streamwise and cross-stream directions, respectively; p is the pressure; ρ is the density of the flow; t is the time; u i ′ u j ′ is the Reynolds stress component, where the superscript of u i is defined as the velocity fluctuation.
In the present study, the two-equation k-ε turbulence models were adopted as our main methods. e k-ε models which determine a turbulent length and time scale by solving two separate transport equations k (the turbulence kinetic energy) and ε (its rate of dissipation) are generally applied with three types: the standard, RNG, and realizable k-ε models.

Standard k-ε Model.
e standard k-ε model proposed by Launder and Spalding [30] contains two transport equations of k and ε. e transport equation k comes from the exact equation, but the transport equation ε was gotten from physical reasoning and has little matter to its mathematically exact counterpart. And they are obtained from the following transport equations: where G k is the generation of turbulence kinetic energy of the mean velocity gradients and, in a manner accordance with the Boussinesq hypothesis [29], G k is evaluated by the following first equation; G b is generated by turbulence kinetic energy due to buoyancy, and Y M is the parameter of the fluctuating dilatation and overall dissipation rate in compressible turbulence; those two parameters are usually set as zero when the fluid considered is incompressible; C 1ε , C 2ε , and C 3ε are constants, given usually by the following second equation (6); σ k and σ ε are the turbulence Prandtl numbers for k and ε, given by default as (6); S k and S ε are user-defined source terms.

RNG k-ε Model.
e RNG k-ε model, which was proposed by Yakhot and Orszag [31], was obtained through renormalization group theory. And the biggest change among the equations is the additional item R ε in ε-equation to be more responsive to the effects of rapid strain and streamline curvature than the standard k-ε model. e RNG k-ε is given by where the additional item R ε is given by where η � Sk/ε, η 0 � 4.38, β � 0.012, and the other constants are given by

Realizable k-ε Model.
Compared with standard k-ε turbulence model, the realizable k-ε model [32] has made some improvements and adjustments, containing an alternative formulation in the viscosity and modifying the ε-equation with deriving from an exact equation. And the modeled transport equations for k and ε are given as follows: where Here, C μ is no longer constant in comparison with the standard and RNG k-ε models; it is calculated as e constants in (11) and (12) are given by C 2 � 1.9, σ k � 1.0, and σ ε � 1.2.

Numerical Schemes and Settings.
e case is in a Cartesian coordinate (x, y) ( Figure 1). e x-axis is aligned with the inlet flow direction, and the y-axis is perpendicular to x and is aligned with the height direction. Rectangular cylinders with various side ratios (R � b/h) are exposed to a uniform stream velocity U � 10 m/s. e height of the box girder is 1 m; thus Re is 6.85 × 10 5 .
An incompressible flow with constant fluid properties is assumed. All geometrical lengths are scaled with h. e boundary conditions used for the computations are applied as follows (seen in Figure 1): a uniform flow (u � U, v � 0) is prescribed at the inlet boundary, the turbulence intensity is set as I � 5% (medium turbulence intensity), and the pressure is set as zero normal gradients. At the outlet, the velocity and turbulence fields (k and ε) are specified as zero normal gradients, and the zero-pressure outlet boundary condition is used. No-slip conditions are prescribed at the surface of the cylinder; the pressure and turbulence fields are set as zero. Symmetry conditions with a frictionless wall (zu/zy � 0, v � 0) are used at the upper and lower boundaries.
Mesh treatments are shown in Figure 2. e O-type grid was adopted for the simulations.
e O-type grids were distributed around the cylinder with a distance of 5 h, according to the hyperbolic tangent laws. Outside of the O-type grids, they were extended far away from the cylinder at a constant growth ratio. e Courant number (Co) was used to adjust the stability and convergence of the calculation. It is defined as Co � UΔt/y 0 , where y 0 is the distance of the first gird to the wall of cylinder and Δt is the time step. In this paper, through adjusting the time step, the Courant number was set as to be around 0.5 (Co < 1 is required according to the criterion), and the calculation of each step could be completed in 20 iterations.
Numerical simulations had been performed by using the commercial software Fluent [33], which allows the analysis of flow around rectangles, by employing an ALE (Arbitrary Lagrangian-Eulerian) formulation based on a colocated finite-volume approach.
e RANS-based modeling technique is well embedded in FLUENT. e scheme applied in the present study was the pressure implicit with splitting of operators (PISO). e solver was pressure based and transient. Gradient is discretized spatially by the least square cell based scheme. Pressure is dealt with a second order spatial discretization scheme and momentum is a second order upwind scheme. e first order implicit was used for time integration.

Definitions of Characteristic Parameters.
e formula of Reynolds number and the Strouhal number is defined as follows:  Mathematical Problems in Engineering where f is the frequency of the vortex shedding. e aerodynamic parameters of time-averaged drag and lift coefficient (C d and C l ), the distribution of mean pressure coefficient (C p ) around the cylinders, and the RMS value or fluctuation of lift coefficient (C lrms ) are where F d and F l are the streamwise and cross-stream timeaveraged force components on the cylinders per unit length, respectively; p and p ∞ are the certain point of the cylinder surface and the reference pressure components of location, respectively; C l (t) is the instantaneous value of lift coefficient; T is the fluctuation period.

Comparison with Different k-ε Models.
In order to study the effects of the three k-ε turbulence models on the results, the square cylinder at Re � 6.85 × 10 5 was taken as an example for comparison and analysis. Considering the fact that the accuracy and the efficiency of the results have a lot to do with the resolution of the grids, four grids with various distances y 0 were discussed. e details of grids and computational results are shown in Table 1. e main aerodynamic integral parameters such as C d , C lrms , and St were presented and compared with available experimental [1,10] and numerical results [4,15,16,34]. e results of C d by Lyn et al. [1], C lrms by Rodi [34], and St by Lyn et al. [1] were chosen as the reference values; the relative deviations of all cases were shown in Table 1.
e cases in Table 1 with different grids were signed with G1∼G4, and different turbulence models were signed with M1∼M3. e minimum relative deviation of drag coefficient, fluctuation of lift coefficient, and Strouhal number were the cases M3G4, M3G4, and M1G1, respectively. e maximum relative deviation of them is M1G1 (or M1G3 and M1G4), M2G1, and M1G3. e detailed comparative analysis of the methods and grids is shown more intuitively in Figure 3.
e comparison of C d with different models is presented in Figure 3(a). For each turbulence model, the value of C d changes a little as y 0 /h ranging from 0.0005 to 0.002, while the distinction between the three models at a certain y 0 /h is significantly different. e difference of the realizable k-ε model with the reference values [1,4,15] is the smallest; the minimum relative deviation is only 1% when y 0 /h is 0.002. e result of RNG k-ε model is slightly larger than the reference values, and the standard k-ε model is smaller; the relative deviation of the standard k-ε is the largest.
eir minimum relative deviation is 19.5% and 10.5%, respectively. us, for the prediction of C d among the three models, the realizable k-ε is better than the others. Figure 3(b) presents the comparison of C lrms . e values of C lrms for each model are almost unchanged when y 0 /h increases, too. For every y 0 /h, the RNG k-ε model is always larger than the reference value [4,15,34]; the standard k-ε model is smaller than the reference values. e realizable k-ε model is slightly less than the reference values [4,15] and nearly the same as Rodi's [34] result. e results of realizable k-ε model are the best. Figure 3(c) shows St of three models with [1,4,15]. e results change irregularly with the increase of y 0 /h in all three models, but it still could be observed that the realizable k-ε model performs the best and the standard and RNG k-ε perform the worst.
e distributions of C p on the surfaces of three models when y 0 /h � 0.002 are plotted in Figure 4. As shown in Figure 4, C p on the front side of the surface (A-B) is positive and is in a good agreement with the experimental [35] and  Mathematical Problems in Engineering numerical [15,36] results for all of models. On the sidesurface (B-C) and back surface (C-D), C p is negative, which indicates that the direction and components of the flow velocity have changed when flow passes by the cylinder and its kinetic energy has been consumed because of the viscous sublayer near the surfaces. So, the relative pressure declines.
e predictions of C p with the three models are different: the standard k-ε model underestimates the pressure (absolute value) on B-C-D sides, and the RNG and realizable k-ε models show a good agreement with the reference studies [15,35,36].
It can be seen that C d , C lrms , and C p of realizable k-ε model show a good agreement with the references compared to the other two models, while the relative deviation of St is still large. e difference is related to the fact that the RANS k-ε models are two-dimensional simplified methods and the three-dimensional effect is ignored which has been discussed in [24][25][26]. For large structures in engineering practice, the drag force and the distribution of surface wind pressure are more important and need to be considered, so that the deviation is acceptable.

Mathematical Problems in Engineering
Overall, the time-averaged integral components change a little when y 0 /h increases except for St; thus the three turbulence models are not sensitive in grids. Among the three turbulence models, for the predicted values of C d , C lrms , St, and C p on the surfaces, the realizable k-ε model's performance is always the best and the standard k-ε model is the worst. As a result, the realizable k-ε is the final model that we decided to choose, and the case M3G4 was taken as the reference for the rectangular cylinders' simulations in the following computations.

Aerodynamic Characteristics with Various R.
Drag coefficient (C d ) varying with different side ratios (R) is shown in Figure 5(a). Compared with the related results [5,10,12,17,20], the present predicted C d shows a good agreement with them. As side ratio (R) grows up from 0.1, the drag coefficient (C d ) increases and reaches its peak at R ≈ 0.62 and then decreases with R which continues to increase; finally C d changes a little after R ≥ 2.8. e changes of C d at various R can be explained as in Bearman and Trueman's [19] studies. When R reaches its critical value (about 0.62), the base cavity gradually decreases due to the interaction of the separation shear layer near the surface of the cylinder, resulting in the maximum of C d . When R exceeds the critical value, with the side-surface increasing longer, the trailing vortex generated by the back corner of the cylinder makes C d smaller. But Nakamura and Ohya [8] explained it in another way. ey claimed that the turbulence scale of the wake, the extra small-scale vortex, and the transverse vibration (sometimes the resonance) might be the main factor when R reaches the critical value of 0.62.
According to the trend of C d , we divided the side ratios into three parts artificially: Regime I (R ≤ 0.62), Regime II (0.62R ≤ 2.8), and Regime III (R ≥ 2.8). Regimes I and II are distinguished by the local maximum of the drag coefficient. Figure 5(b) shows the base pressure coefficient (−C pb ) at various R. e tendency of −C pb is the same as C d . It goes up in Regime I, then goes down in Regime II, and finally becomes nearly stable in Regime III. e variation is also in a good agreement with the available results [5,10,17]. e fluctuation of the lift coefficient (C lrms ) varying with R is shown in Figure 5(c). Lack of enough experimental data and some typical published results with RANS and LES models [15,17] were compared. e fluctuation of lift coefficient C lrms increases with the increasing of R in Regime I, then it decreases in Regime II, and next it increases suddenly at R � 2.8 and finally declines again in Regime III. e predicted C lrms is almost the same as the results of Shimada and Ishihara [15] and is also close to Sohankar's [17] results. However, the values of C lrms in the present study are relatively smaller compared to the LES results [17] in Regimes II and III. According to Reynolds and Hussain's [37] studies, the fluctuation value consisted of periodic and stochastic fluctuation component, but RANS models ignore the stochastic turbulent fluctuation, and the proportion of stochastic fluctuation component depends a lot on R, so that the present results of C lrms are always underestimated and show a difference at different R [15]. As seen in Figure 5(c), the proportion of stochastic fluctuation in Regime II is larger than the other two regimes so that the present result of C lrms is much lower than the LES results [17] in Regime II. Figure 5(d) shows that the Strouhal number (St) changes with various R. e predicted St changes a little in Regime I and decreases with R increases in Regimes II and III. But St presents a discontinuity at R � 2.8; this is due to the reattachment of the vortex when R is around 2.8. e result is just the same as the results of Nakaguchi et al. [5], Norberg [10], Okajima et al. [6,12], and Sohankar [17].
Generally, the aerodynamic characteristics obtained from present RANS realizable k-ε model are able to give reasonable predictions at R � 0.1-4.0. e present results also indicate that the Reynolds dependence of the results has become much smaller for higher Reynolds numbers, which has been mentioned by Sohankar et al. [4,17] and Nakamura and Ohya [8]. Figure 6 shows the distribution of the surface pressure coefficient (C p ) along the front and back side of cylinders with different side ratios (R � 0.1, 0.2, 0.4, 0.62, 0.8, 1.0, 1.5, 2.0, 2.8, and 3.0).

Flow Field Characteristics with Various R.
e pressure coefficient C p on the front side for the side ratios collapses together, while the tendency of C p on the back side with various R varies apparently, which suggests that the differences in drag are mainly due to the differences in the pressure distribution along the back side.
Meanwhile, the pressure distribution on the back side varies with different side ratios: the pressure changes apparently when R < 1.0, especially when R � 0.62, and it is almost constant when R ≥ 1.0. As the regimes have been Mathematical Problems in Engineering defined above, the difference of C p varying along the back side mainly happens in Regimes I and II. ese might come from the differences of vortex distributions of the wake with various R, concerning the fact that C p is actually generated by vortex shedding of the cylinder. Figure 7 presents two parameters of the wake with different R: recirculation length (l R ) and the negative maximum velocity of back-flow (U BF , in absolute value). e recirculation length l R can be seen as the index of the wake distance and U BF indicates the strength of the wake. As seen in Figure 7, the values of l R and U BF varying with different R can also be divided into three regimes as mentioned above. In Regime I, l R and U BF decrease when R increases. In Regimes II and III, l R and U BF also increase with R but grow up slowly in Regime III. e tendencies of l R and U BF varying with side ratios are just the opposite of C d and −C pb . Figure 8 shows the vorticity contours in the wake zone for various cylinders at R � 0.1, 0.2, 0.4, 0.62, 0.8, 1.0, 1.5, 2.0, 2.8, and 3.0 as the lift coefficient reaches its maximum value. Because the size and position of the vortex in the wake zone   would affect the pressure around the cylinders and the drag forces acting on the cylinder, it can be clearer to understand the aerodynamic characteristics mentioned above through the instantaneous vorticity contours. As can be seen in Figure 8, the fluid separates from the surface and forms the shear layer at the leading edge of the cylinders and then produces a base cavity near the back sides.
It is worth mentioning that the two-dimensional vorticity in Figure 8 is quite different from the three-dimensional vortex structure which was presented by Norberg [10], Sohankar et al. [4,17], and Cao and Tamura [38]. As mentioned above, the exact value of turbulent fluctuation consisted of periodic and stochastic component, but the stochastic turbulent fluctuation is not modeled in RANS models. ese stochastic perturbations, corresponding to the whole averaged flow fields, are small-scale eddies and are treated as a part of eddy viscosity [15]. erefore, such smallscale eddies were not visible in the present flow patterns in Figure 8. e present vorticity contours were more of a qualitative analysis of the flow field.
In Regime I, as seen in Figures 8(a)-8(d), the location of the generation of the vortex moves closer to back side when R increases. Almost at R � 0.62, the generation of the vortex becomes the strongest. us, −Cpb and C d reach their local maximum. When R continues to grow up to Regimes II and III, as seen in Figures 8(e)-8(j), the location of the vortex begins to move further away from the cylinder, corresponding to the decline of -C pb and C d . On the other hand, it could be observed that the flow separates from the front corners and is not reattached to the side faces in Regimes I and II (Figures 8(a)-8(i)). As seen in Figures 8(i) and 8(j), when side ratios are in Regime III, there is a flow reattachment phenomenon of wake vortex around the cylinder, which gives a good interpretation for the discontinuity of St and C lrms in Figures 5(c) and 5(d).
In addition, the wake vortex has an important impact on the distributions of the pressure in the back side: the closer the vortex to the back side, the larger the pressure (in absolute) in the back side. erefore, the distribution of the pressure in back side is changing when the vortex is closer to the cylinder at R < 1.0, shown in Figure 6, and pressure maintains almost constant when the vortex moves further away from the cylinder at R ≥ 1.0. Particularly, when R � 0.62, seen in Figure 8(d), the vortex attaches to the back side, and the pressure of the central position in the back side is the largest (in absolute value), because it just is located at the center of the vortex.

Discussion of the ree Regimes.
Much attention has been paid to the effect of side ratios of rectangular cylinders on the aerodynamic characteristics. ere are mainly two critical side ratios for aerodynamic parameters. e first one is located at R � 0.62 for the maximum of mean drag and base pressure coefficient, and the other is for the discontinuity of Strouhal number when the side ratio equals 2.8 [5,6,10,12,[15][16][17]19].
However, there are little clear subregional discussion on the side ratios, and the characteristics of each region are only simply described. Nakaguchi et al. [5] listed numbers of experimental results at side ratios ranging from 0.1 to 4.0, but correlation analysis and discussion were little. Okajima et al. [6,7,12] discussed C d , St, and flow patterns at R � 0.6-8.0 containing wider range and got similar results. But some other characteristics such as-C pb and C lrms were not analyzed, as well as the range of R < 0.6. Norberg [10] conducted C d , C p , and St at R � 1, 1.67, 2.5, and 3 at angles of attack 0°∼90°; the series of side ratios are not enough to distinguish them. Shimada and Ishihara [15] researched the cylinders at R � 0.6∼8.0. Tian et al. [16] did not contain side ratios which are larger than 1.0. us, three regimes at various side ratios were divided artificially in the present study which has been mentioned above. e main distinctions of three regimes are shown in Figure 9. Regimes I and II are divided on the basis of C d and

Mathematical Problems in Engineering
−C pb at R � 0.62. Regimes II and III are distinguished by the discontinuity of St at R � 2.8. As shown in Figure 9, in Regime I, the aerodynamic parameters of C d , −C pb , and C lrms increase when R grows up and reach their local maximum at R � 0.62. St keeps almost constant in Regime I. e pressure coefficient C p (in absolute value) on back side of different cylinders also shows similar increasing tendency, and the distribution of C p along the back side changes to be larger when R increases. e wake characteristics of l R and U BF decrease with the increasing of R. e vortex in the wake (seen as the flow visualization in Figure 9) moves closer to the cylinder and becomes the closest at R � 0.62.
In Regime II, C d , −C pb , C lrms , St, and C p present a decreasing tendency with increase of R. e distribution of C p along the back side changes to be smaller when R grows up and becomes almost constant at R > 1.0. When R increases, l R and U BF increase. e vortex in the wake begins to move away from the back side. In Regime III, C d , −C pb , and C p continue to decrease but change a little when R increases. C lrms and St present a sudden jump at R � 2.8 and still show the same decreasing tendency at R > 2.8. e recirculation length l R and the recirculation velocity U BF grow up with R but change a little. e vorticity contours show a reattachment at R ≥ 2.8, and the vortex moves further away from cylinder. e characteristic of every regime is obvious and presents a difference with each other, and their boundaries represent a critical position. e division of the three regimes could give a clearer and more systematic cognition of the flow around rectangular cylinders. Moreover, the tendency of aerodynamic parameters varying with side ratios in high Reynolds number could give some guidance for practical engineering in design and manufacture. e drag force, base pressure, and the lift fluctuation appear to be at peak when side ratio equals 0.62. us, this side ratio should be avoided to decline the structural drag force and to increase its stability in the wind fields. e discontinuity of Strouhal number and change of the lift coefficient fluctuation at side ratio of around 2.8 should also be noticed when referring to the vibrations of the structures. In other words, the side ratios near the boundaries of the three regimes should be taken into account in practical engineering.

Conclusions
In this paper, the flow around the rectangular cylinders with various side ratios for much higher Reynolds numbers (Re � 68.5 × 10 5 ) was presented. Considering the cost and applying for the practical engineering, the two-dimensional RANS turbulence models were used for the computations. ree typical k-ε turbulence models were compared and discussed firstly. en, the aerodynamic characteristics such as the drag coefficient, the fluctuation of lift coefficient, the back pressure coefficient, and Strouhal number at various side ratios were presented and compared with the wind tunnel experimental results [5,6,10,20] and numerical results [15,17]. Finally, the flow field characteristics around the cylinders were presented and illustrated in detail. e main conclusions in the research are summarized as follows: (1) ree typical RANS k-ε turbulence models were compared for the prediction of aerodynamic integral parameters. rough the comparisons of C d , C lrms , St, and C p over the cylinder, the realizable k-ε model was found to perform the best and was finally chosen as the main method for the rectangular computations.
(2) According to the tendency of aerodynamic parameters varying with side ratios, three regimes (Regime I for R ≤ 0.62; Regime II for 0.62R ≤ 2.8; Regime III for R ≥ 2.8) were divided. In Regime I, C d , -C pb , and C lrms increase with the increasing of R and reach their local maximum at R � 0.62, and St changes a little; in Regime II, C d , −C pb , C lrms , and St decrease with R; and, in Regime III, C d , −C pb , C lrms , and St decrease less with R, but C lrms and St show a discontinuity at R � 2.8. (3) e flow field characteristics around the rectangular cylinders with various side ratios could also be distinguished with the three regimes. e contribution of C p along the back side presents a similar tendency as C d and-C pb in the three regimes. e wake fields parameters l R and U BF present an opposite tendency of C d . e distinction of vortex in three regimes is presented with the visual vorticity contours. e closest vortex at R � 0.62 and the reattachment of vortex at R � 2.8 are shown more clearly and intuitively through the vorticity. e two-dimensional RANS realizable k-ε model adopted in this paper could give a reasonable prediction for the basic aerodynamic parameters and flow field characteristics with the comparison with published literatures. However, the derivative is large for the lift fluctuation and vortex shedding frequency when referring to the structural vibration and stability. erefore, in the following study, the practicability and reliability of the model still need to be verified when it involves vortex-induced vibration, flutter, and other vibration problems.

Data Availability
e data used to support the findings of the study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest as regards the publication of this paper. Mathematical Problems in Engineering 11