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

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

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

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 2163928 | https://doi.org/10.1155/2020/2163928

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

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

Academic Editor: Zhengbiao Peng

Received03 Jul 2020

Accepted08 Sept 2020

Published15 Sept 2020

Abstract

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⁴). This 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. Three unsteady Reynolds-averaged Navier-Stokes (RANS) two-equation k-ε turbulence models (standard, RNG, and realizable) were adopted in the study. The 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. The distribution of surface pressure over cylinders, the wake parameters, and vorticity contours of the rectangular cylinders with different side ratios were discussed.

1. 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.

The 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⁴ 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⁴) 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⁴.

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⁴ 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. They 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. Their 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⁴ using DVM. Their 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⁴ numerically using a modified k-ε model. Their method presents better results than ordinary 2D analysis in vortex shedding, drag coefficients, and distributions of mean surface pressure. They reported a discontinuity in Strouhal number at both R = 2.8 and R = 6. They 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⁵ were simulated by Sohankar [17] using two LES approaches. The 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⁴), Bearman and Trueman [19] (Re = 2.2 × 10⁴), Mizota et al. [20] (Re = 2.1 × 10⁵), and Norberg [10] (Re = 1.3 × 10⁴) and numerical studies of Shimada and Ishihara [15] (Re = 2.2 × 10⁴), Sohankar [17] (Re = 10⁵), and Tian et al. [16] (Re = 2.14 × 10⁴). Most of them were performed at Re ≤ 10⁵. 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⁵ 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⁷. As a result, it is necessary and meaningful to study the flow around rectangles with various R when Re is much higher. They 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. The main global parameters and mean and instantaneous-flow structures would be discussed in this work.

The 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.

2. 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. They 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. The 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]. That 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.1. RANS k-ε Turbulence Model

The Reynolds-averaged incompressible Navier-Stokes equations [29] are expressed aswhere u_i is the velocity component in x_i coordinate directions and, for the present 2D simulation, i, j = 1, 2, x₁ and x₂ are the streamwise and cross-stream directions, respectively; p is the pressure; ρ is the density of the flow; t is the time; 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. The 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.

2.1.1. Standard k-ε Model

The standard k-ε model proposed by Launder and Spalding [30] contains two transport equations of k and ε. The 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.

2.1.2. RNG k-ε Model

The 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. The RNG k-ε is given bywhere the additional item R_ε is given bywhere , = 4.38, β = 0.012, and the other constants are given by

2.1.3. 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

The constants in (11) and (12) are given by C₂ = 1.9, σ_k = 1.0, and σ_ε = 1.2.

2.2. Numerical Schemes and Settings

The case is in a Cartesian coordinate (x, y) (Figure 1). The 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. The height of the box girder is 1 m; thus Re is 6.85 × 10⁵.

An incompressible flow with constant fluid properties is assumed. All geometrical lengths are scaled with h. The normalized dimensions of the domain (L_U, L_D, h, and H, seen in Figure 1) are selected as 12, 30, 1, and 20, respectively. The R is selected as 0.1, 0.2, 0.3, 0.4, 0.5, 0.62, 0.7, 0.8, 0.9, 1.0, 1.5, 2.0, 2.8, 3.0, and 4.0.

The 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 () are used at the upper and lower boundaries.

Mesh treatments are shown in Figure 2. The O-type grid was adopted for the simulations. The 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. The Courant number (Co) was used to adjust the stability and convergence of the calculation. It is defined as , where y₀ 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.

(a)

(b)

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. The RANS-based modeling technique is well embedded in FLUENT. The scheme applied in the present study was the pressure implicit with splitting of operators (PISO). The 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. The first order implicit was used for time integration.

2.3. Definitions of Characteristic Parameters

The formula of Reynolds number and the Strouhal number is defined as follows:where f is the frequency of the vortex shedding.

The 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) arewhere F_d and F_l are the streamwise and cross-stream time-averaged 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.

3. Results and Discussion

3.1. 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⁵ 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₀ were discussed. The details of grids and computational results are shown in Table 1. The 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]. The 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.

The cases in Table 1 with different grids were signed with G1∼G4, and different turbulence models were signed with M1∼M3. The minimum relative deviation of drag coefficient, fluctuation of lift coefficient, and Strouhal number were the cases M3G4, M3G4, and M1G1, respectively. The maximum relative deviation of them is M1G1 (or M1G3 and M1G4), M2G1, and M1G3. The detailed comparative analysis of the methods and grids is shown more intuitively in Figure 3.

(a)

(b)

(c)

The 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₀/h ranging from 0.0005 to 0.002, while the distinction between the three models at a certain y₀/h is significantly different. The 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₀/h is 0.002. The 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. Their minimum relative deviation is 19.5% and 10.5%, respectively. Thus, 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. The values of C_lrms for each model are almost unchanged when y₀/h increases, too. For every y₀/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. The realizable k-ε model is slightly less than the reference values [4, 15] and nearly the same as Rodi’s [34] result. The results of realizable k-ε model are the best.

Figure 3(c) shows St of three models with [1, 4, 15]. The results change irregularly with the increase of y₀/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.

The distributions of C_p on the surfaces of three models when y₀/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 numerical [15, 36] results for all of models. On the side-surface (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. The 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. The 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–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.

Overall, the time-averaged integral components change a little when y₀/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.

3.2. Effect of Various R

3.2.1. 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. The 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. They 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.

(a)

(b)

(c)

(d)

Figure 5(b) shows the base pressure coefficient (−C_pb) at various R. The 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. The variation is also in a good agreement with the available results [5, 10, 17].

The 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. The 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. The 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. The 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. The 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. The 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].

3.2.2. Flow Field Characteristics with Various R

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). The 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 defined above, the difference of C_p varying along the back side mainly happens in Regimes I and II. These 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). The 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. The 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.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

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. These stochastic perturbations, corresponding to the whole averaged flow fields, are small-scale eddies and are treated as a part of eddy viscosity [15]. Therefore, such small-scale eddies were not visible in the present flow patterns in Figure 8. The 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. Thus, −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. Therefore, 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.

3.2.3. Discussion of the Three Regimes

Much attention has been paid to the effect of side ratios of rectangular cylinders on the aerodynamic characteristics. There are mainly two critical side ratios for aerodynamic parameters. The 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–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.

Thus, three regimes at various side ratios were divided artificially in the present study which has been mentioned above. The main distinctions of three regimes are shown in Figure 9. Regimes I and II are divided on the basis of C_d and −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. The 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. The wake characteristics of l_R and U_BF decrease with the increasing of R. The 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. The 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. The 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. The recirculation length l_R and the recirculation velocity U_BF grow up with R but change a little. The vorticity contours show a reattachment at R ≥ 2.8, and the vortex moves further away from cylinder.

The characteristic of every regime is obvious and presents a difference with each other, and their boundaries represent a critical position. The 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. The drag force, base pressure, and the lift fluctuation appear to be at peak when side ratio equals 0.62. Thus, this side ratio should be avoided to decline the structural drag force and to increase its stability in the wind fields. The 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.

4. Conclusions

In this paper, the flow around the rectangular cylinders with various side ratios for much higher Reynolds numbers (Re = 68.5 × 10⁵) was presented. Considering the cost and applying for the practical engineering, the two-dimensional RANS turbulence models were used for the computations. Three typical k-ε turbulence models were compared and discussed firstly. Then, 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. The main conclusions in the research are summarized as follows:(1)Three typical RANS k-ε turbulence models were compared for the prediction of aerodynamic integral parameters. Through 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)The flow field characteristics around the rectangular cylinders with various side ratios could also be distinguished with the three regimes. The contribution of C_p along the back side presents a similar tendency as C_d and–C_pb in the three regimes. The wake fields parameters l_R and U_BF present an opposite tendency of C_d. The distinction of vortex in three regimes is presented with the visual vorticity contours. The closest vortex at R = 0.62 and the reattachment of vortex at R = 2.8 are shown more clearly and intuitively through the vorticity.

The 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. Therefore, 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

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

Conflicts of Interest

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

Acknowledgments

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

References

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, no. 1, pp. 285–319, 1995.
View at: Publisher Site | Google Scholar
S. C. Yen and C. W. Yang, “Flow patterns and vortex shedding behavior behind a square cylinder,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 99, no. 8, pp. 868–878, 2011.
View at: Publisher Site | Google Scholar
W. Rodi, J. H. Ferziger, M. Breuer, and M. Pourquie, “Status of large eddy simulation; results of a workshop,” ASME Journal of Fluids Engineering, vol. 2, no. 119, pp. 248–262, 1997.
View at: Publisher Site | Google Scholar
A. Sohankar, “Flow over a bluff body from moderate to high reynolds numbers using large eddy simulation,” Computers & Fluids, vol. 35, no. 10, pp. 1154–1168, 2006.
View at: Publisher Site | Google Scholar
H. Nakaguchi, K. Hashimoto, and S. Muto, “An experimental study on aerodynamic drag of rectangular cylinders,” Journal of the Japan Society for Aeronautical and Space Sciences, vol. 168, no. 16, pp. 1–5, 1968.
View at: Publisher Site | Google Scholar
A. Okajima, “Strouhal numbers of rectangular cylinders,” Journal of Fluid Mechanics, vol. 123, 1982.
View at: Publisher Site | Google Scholar
A. Okajima, “Flow around a rectangular cylinder with a section of various width/height ratios,” Wind Engineers, JAWE, vol. 1983, no. 17, pp. 1–19, 1983.
View at: Publisher Site | Google Scholar
Y. Nakamura and Y. Ohya, “The effects of turbulence on the mean flow past two-dimensional rectangular cylinders,” Journal of Fluid Mechanics, vol. 149, no. 149, pp. 255–273, 1984.
View at: Publisher Site | Google Scholar
Y. Ohya, “Note on a discontinuous change in wake pattern for a rectangular cylinder,” Journal of Fluids and Structures, vol. 8, no. 3, pp. 325–330, 1994.
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. 1-3, no. 49, pp. 187–196, 1993.
View at: Publisher Site | Google Scholar
S. Nakagawa, K. Nitta, and M. Senda, “An experimental study on unsteady turbulent near wake of a rectangular cylinder in channel flow,” Experiments in Fluids, vol. 27, no. 3, pp. 284–294, 1999.
View at: Publisher Site | Google Scholar
A. Okajima, “Numerical simulation of flow around rectangular cylinders,” Journal of Wind Engineering & Industrial Aerodynamics, vol. 33, no. 1-2, pp. 171–180, 1990.
View at: Publisher Site | Google Scholar
A. Okajima, T. Nagahisa, and A. Rokugoh, “A numerical analysis of flow around rectangular cylinders,” JSME International Journal. Ser. 2, Fluids Engineering, Heat Transfer, Power, Combustion, Thermophysical Properties, vol. 33, no. 4, pp. 702–711, 1990.
View at: Publisher Site | Google Scholar
M. V. Ian Taylor, “Prediction of unsteady flow around square and rectangular section cylinders using a discrete vortex method,” Journal of Wind Engineering & Industrial Aerodynamics, vol. 1-3, no. 82, pp. 247–269, 1999.
View at: Publisher Site | Google Scholar
K. Shimada and T. Ishihara, “Application of a modified k - ε model to the prediction of aerodynamic characteristics of rectangular cross-section cylinders,” Journal of Fluids and Structures, vol. 1, no. 16, pp. 465–485, 2002.
View at: Publisher Site | Google Scholar
X. Tian, M. C. Ong, J. Yang, and D. Myrhaug, “Unsteady RANS simulations of flow around rectangular cylinders with different aspect ratios,” Ocean Engineering, vol. 58, no. 58, pp. 208–216, 2013.
View at: Publisher Site | Google Scholar
A. Sohankar, “Large eddy simulation of flow past rectangular-section cylinders: side ratio effects,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 96, no. 5, pp. 640–655, 2008.
View at: Publisher Site | Google Scholar
S. U. Islam, C. Y. Zhou, A. Shah, and P. Xie, “Numerical simulation of flow past rectangular cylinders with different aspect ratios using the incompressible lattice boltzmann method,” Journal of Mechanical Science and Technology, vol. 26, no. 4, pp. 1027–1041, 2012.
View at: Publisher Site | Google Scholar
P. W. Bearman and D. M. Trueman, “An investigation of the flow around rectangular cylinders,” Aeronautical Quarterly, vol. 23, no. 3, pp. 229–237, 1972.
View at: Publisher Site | Google Scholar
A. Okajima and T. Mizota, “Aerodynamic characteristics of fundamental structures (part 1),” Wind Engineers JAWE, vol. 1988, no. 36, 1988.
View at: Publisher Site | Google Scholar
W. Yupu, C. Wenming, R. DU, W. Shubiao, and D. Yong, “Bionic design method of crane box girder wind load reduction based on ostracion-cubicus,” Journal of Southwest Jiaotong University, vol. 55, no. 3, pp. 664–671, 2020.
View at: Google Scholar
A. Scarabino, J. Marañón Di Leo, J. S. Delnero, and F. Bacchi, “Drag coefficients and strouhal numbers of a port crane boom girder section,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 93, no. 6, pp. 451–460, 2005.
View at: Publisher Site | Google Scholar
S. Murakami, “Overview of turbulence models applied in CWE-1997,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 74-76, no. 98, pp. 1–24, 1998.
View at: Publisher Site | Google Scholar
K. Hanjalic, “Will RANS survive LES? A view of perspectives,” Journal of Fluids Engineering, vol. 127, no. 127, p. 831, 2005.
View at: Publisher Site | Google Scholar
F. R. Menter and R. Florian, “Review of the shear-stress transport turbulence model experience from an industrial perspective,” International Journal of Computational Fluid Dynamics, vol. 23, no. 4, pp. 305–316, 2009.
View at: Publisher Site | Google Scholar
B. Blocken, “LES over RANS in building simulation for outdoor and indoor applications: a foregone conclusion?,” Building Simulation, vol. 11, no. 5, 2018.
View at: Publisher Site | Google Scholar
D. Sun, J. S. Owen, and N. G. Wright, “Application of the k-ω turbulence model for a wind-induced vibration study of 2D bluff bodies,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 97, no. 2, pp. 77–87, 2009.
View at: Publisher Site | Google Scholar
C. Mannini, A. Šoda, and G. Schewe, “Unsteady RANS modelling of flow past a rectangular cylinder: investigation of Reynolds number effects,” Computers & Fluids, vol. 39, no. 9, pp. 1609–1624, 2010.
View at: Publisher Site | Google Scholar
J. O. Hinze, Turbulence, McGraw-Hill Publishing Co, New York, NY, USA, 1975.
B. E. Launder and D. B. Spalding, Lectures in Mathematical Model of Turbulence, Academic Press, London, England, UK, 1972.
V. Yakhot and O. A. Orszag, “Renormalization group analysis of turbulence, I. Basic theory,” Journal of Scientific Computing, vol. 1, no. 1, pp. 3–51, 1986.
View at: Publisher Site | Google Scholar
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
F. Inc, Fluent 6.3 User’s Guide, Fluent Inc, Lebanon, NH, USA, 2006.
W. Rodi, “Comparison of LES and RANS calculations of the flow around bluff bodies,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 69-71, pp. 55–75, 1997.
View at: Publisher Site | Google Scholar
B. E. Lee, “The effect of turbulence on the surface pressure field of a square prism,” Journal of Fluid Mechanics, vol. 69, no. 2, pp. 263–282, 1975.
View at: Publisher Site | Google Scholar
S. Oka and T. Ishihara, “Numerical study of aerodynamic characteristics of a square prism in a uniform flow,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 97, no. 11-12, pp. 548–559, 2009.
View at: Publisher Site | Google Scholar
W. C. Reynolds and A. K. M. F. Hussain, “The mechanics of an organized wave in turbulent shear flow. Part 3. theoretical models and comparisons with experiments,” Journal of Fluid Mechanics, vol. 54, no. 2, pp. 263–288, 1972.
View at: Publisher Site | Google Scholar
Y. Cao and T. Tamura, “Large-eddy simulations of flow past a square cylinder using structured and unstructured grids,” Computers & Fluids, vol. 137, pp. 36–54, 2016.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Shubiao 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

Views

1141

Downloads

792

Citations

Mathematical Problems in Engineering

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

Abstract

1. Introduction

2. Numerical Methods

2.1. RANS k-ε Turbulence Model

2.1.1. Standard k-ε Model

2.1.2. RNG k-ε Model

2.1.3. Realizable k-ε Model

2.2. Numerical Schemes and Settings

2.3. Definitions of Characteristic Parameters

3. Results and Discussion

3.1. Comparison with Different k-ε Models

3.2. Effect of Various R

3.2.1. Aerodynamic Characteristics with Various R

3.2.2. Flow Field Characteristics with Various R

3.2.3. Discussion of the Three Regimes

4. Conclusions

Data Availability

Conflicts of Interest

Acknowledgments

References

Copyright

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