The Seventh Asia-Pacific Conference on Wind Engineering, November 8-12, 2009, Taipei, Taiwan EFFECTS OF SECTION SIZE ON AERODYNAMIC FORCES ON AN ELLIPTIC CYLINDER UNDER SHORT-RISE-TIME GUSTS Takashi Takeuchi1, Junji Maeda2, Tomohiko Hayata3 and Hiromasa Kawashita4 JSPS Research Fellow (Graduate), 3,4 Graduate, Graduate School of Human-Environment Studies, Kyushu University, 6-10-1 Hakozaki, Higashi-ku, Fukuoka 812-8581, Japan, [email protected], [email protected], [email protected] 2 Professor, Faculty of Human-Environment Studies, Kyushu University, 6-10-1 Hakozaki, Higashi-ku, Fukuoka 812-8581, Japan, [email protected] 1 ABSTRACT The characteristics of aerodynamic forces acting on an elliptic cylinder with an angle of attack subjected to short-rise-time gusts were investigated using a specially-equipped gust wind tunnel. Remarkable overshoot phenomena of aerodynamic forces were indicated, and we analyzed their patterns as they changed with the section-size of the cylinder. The overshoot coefficients tended to increase with a decrease in the non-dimensional rise time regardless of the section-size of cylinder. The overshoot coefficients were explainable as one pattern by the non-dimensional rise time regardless of the section-size of cylinder. Finally, we showed that a computational fluid-dynamics-analyzed relationship between the overshoot coefficient and the non-dimensional rise time agreed well with the experimental results. KEYWORDS: AERODYNAMIC FORCE, SHORT-RISE-TIME GUST, OVERSHOOT OF WIND FORCE, SECTION-SIZE Introduction Many reports have clarified detailed characteristics of flow over a bluff body and wind force on such a body under a stationary or a quasi-steady turbulent flow that has an intensity of turbulence of about 10%. However, there is a possibility that wind gusts with a very short rise time occur in tornadoes or high winds that are associated with typhoons or frontal passages. There have been several reports on the characteristics of unsteady aerodynamic forces on structures under such gusts. [Taneda (1972)] investigated the unsteady lift acting on an elliptic cylinder rapidly started at an angle of attack using a water tank test and reported that a remarkably big lift appeared just after starting. [Sarpkaya (1966)] showed that the drag of a body grew up to about 25 percent higher during the growth of the first pair of vortices than in a steady flow, using an impulsive flow test over circular cylinders in a vertical water tunnel, adding some results of potential flow analyses around the circular cylinders. [Morison and O’Brien et al. (1950) and Keulegan and Carpenter (1958)] evaluated the unsteady aerodynamic characteristics using a semi-empirical equation added to an inertia term proportional to flow acceleration. [Nomura and Kitamura et al. (2008)] computed the unsteady drag acting on a square cylinder under a sudden change of flow speed and reported that the drag component proportional to flow acceleration played quite an important role in the total unsteady drag when the flow speed was relatively low. [Shiraishi and Matsumoto et al. (1982)] evaluated unsteady aerodynamic characteristics using aerodynamic indicial functions. [Matsumoto and Shimamura et al. (2007)] measured a transient drag on a twodimensional cylindrical model under a suddenly-changing wind speed using a wind tunnel test The Seventh Asia-Pacific Conference on Wind Engineering, November 8-12, 2009, Taipei, Taiwan with a working section of 200mm by 200mm square and reported an overshoot phenomenon in which the drag increased by about 20 % compared to the force in a steady flow. [Yamaguchi and Kashiwagi et al. (2008)] investigated unsteady aerodynamic characteristics of rectangular prisms and a truck model subjected to quickly-accelerated flow using a multifan type blow-down wind tunnel with a working section of 2540mm by 1800mm which was composed of 99 small fans, each of which was independently controlled by a computer. [Takeuchi and Maeda et al. (2008)] investigated the effects of the rise time of a step-functionlike gust on the overshoot phenomenon of wind forces acting on a railcar-like body using a specially-equipped wind tunnel, which could generate gusts with a rise time of 0.2 to 5 seconds from a flat calm by controlling the rotation speed of the blade rows, and reported that the overshoot phenomenon was more remarkable in the case of gusts with a shorter rise time. [Nakamura and Maeda et al. (2008)] investigated effects of a step-function-like gust on local wind pressures of a gable-roof body using the same specially-equipped wind tunnel and reported that the step-function-like gust produced an overshoot phenomenon of wind pressure bringing much larger local wind pressures than a steady flow on the surface of the body, and that the overshoots of the surface pressure related to those of the drag acting on the body. [Kawashita and Maeda et al. (2008)] confirmed some overshoot phenomena of wind forces acting on an elliptic cylinder with some angles of attack using the same specially-equipped wind tunnel, and reported that a non-dimensional parameter of a rise time of step-functionlike gust considerably influenced the occurrence of the overshoot phenomenon. As mentioned, it was apparent that an overshoot phenomenon bringing a much bigger wind force than in a steady flow occurs in the case of a structure subjected to a short-rise-time gust. To apply the phenomenon to a gust resistant measure for full-scale structures, it is necessary to clarify the effects of the scale of a structure on the overshoot phenomenon. In this paper, with the purpose of verifying the accuracy of the evaluation of the overshoot phenomenon by the non-dimensional parameter of rise time of a gust, as an example, we investigated and here discuss the effects of elliptic cylinder section-size on overshoot phenomena of wind forces under short-time-rise gusts using a specially equipped wind tunnel. And we weigh the results of our experiment against the results from simulation using CFD software. General Specifications of the Wind Tunnel Test We used a wind tunnel of the Eiffel type at Kyushu University. The site plan of our testing system is illustrated in Figure 1. The section area of the working space was 1.5m by 1.5m and the available length was 3m. A step-function-like gust was generated by pulsating rapid rotation of flat blade-rows. The drag and lateral force acting on the specimen were measured using an aerodynamic balance, and the wind speed in the working space was confirmed by a hot-wire anemometer and an ultrasonic anemometer. An example of the time 200 1,300 1,500 Dummy model Hot-wire anemometer Supersonic anemometer Wind Specimen Aerodynamic balance 700 900 400 a) Site plan Plane board Blade row (top) Blade row (front) Blade row (bottom) Dummy Specimen (mm) b) Installation of specimen and dummy Figure 1: Site plan of wind tunnel The Seventh Asia-Pacific Conference on Wind Engineering, November 8-12, 2009, Taipei, Taiwan U Target wind velocity, Ut (mm) 0.5D D Wind Small-cylinder D = 70 mm Rise time, tr t Specimen α Large-cylinder D = 140 mm Wind Aerodynamic balance Figure 2: Example of time evolution of b) Bird’s-eye view wind velocity measured under step- a) Section view of specimen Figure 3: Specimens function-like gust evolution of the wind velocity generated by control of the blade rows is shown in Figure 2. The rise time of gust, tr, is defined as the time required for an approaching wind to reach a target wind, and the target wind velocity, Ut, is referred as the reached wind. The measured data was scanned at a frequency of 1000Hz. The specimens were two elliptic cylinders of axial ratio 2:1, 500mm high, and 70mm and 140mm in major axis, as shown in Figure 3, and referred to as “Small-cylinder” and “Large-cylinder”, respectively. They were fixed to the aerodynamic balance with an angle of attack, α, from 0 degrees to 90 degrees and the drag and lateral forces acting on them were measured. The angle of attack was defined as 0 degrees when the elliptic cylinder was fixed to the aerodynamic balance in such a way that its major axis was parallel to the wind direction. Additionally, a dummy 350mm in height was set to realize a two-dimensional flow around the specimens, as shown in Figure 1. The dummy was placed at a distance of several millimeters from the specimens so as not to measure the wind force on the dummy. They were subjected to step-function-like gusts with a target wind velocity, Ut, between 2.0m/s and 7.0m/s and rise time, tr, between 0.2sec and 1.4sec. Experimental Results and Discussions Overshoot of wind forces Figures 4 and 5 show the time evolutions of the wind velocity, drag and lateral force on Small-cylinder and Large-cylinder at a 45-degree angle of attack respectively, when the wind velocity reached 2.0m/s in the rise times of 0.2sec and 1.4sec. In the case of the rise time of 0.2sec, we confirmed the overshoot phenomena of aerodynamic forces, which reached much larger values than the steady values of the drag and lateral forces. The overshoot phenomena achieved peaks when wind velocity reached the target wind velocity, and they decreased immediately and approached a steady state gradually. The time evolution of the wind velocity for Large-cylinder was almost the same as that for Small-cylinder but the time evolutions of the wind forces on Large-cylinder were different from those on Small-cylinder in that the drag on Large-cylinder had two peaks in 0.3 seconds from starting, and the time the lateral force on Large-cylinder achieved a peak was later than that on Small-cylinder. Figure 6 shows effects of the dummy on the overshoot of wind force acting on Largecylinder at a 45-degree angle of attack in the case of a target wind velocity of 2.0m/s and a rising time of 0.2sec. The time evolution of the drag was not affected by the installation of the dummy. However, the peak value of the lateral force with the dummy was larger than that without the dummy. And after the peak, the lateral force which vibrated in around 3.5Hz with the dummy was bigger than that without the dummy. The lateral force acting on the specimen vibrated under steady flow too, and the frequency of the vibration was proportional to the gust speed and inversely proportional to the section-size of specimen. Thus, it would appear that The Seventh Asia-Pacific Conference on Wind Engineering, November 8-12, 2009, Taipei, Taiwan the vibration was caused by vortex shedding. In this case, the Strouhal number of Largecylinder at a 45-degree angle of attack under steady flow was estimated to be about 0.2. The reference length for the Strouhal number was defined as the length of the flow direction of the elliptic cylinder. Effects of the angle of attack of elliptic cylinder on the overshoot phenomenon Figure 7 shows the angle of attack plotted against the wind force coefficient of Largecylinder, when the target wind velocity reached 5.0m/s in a rising time of 0.2sec. The wind force coefficient is defined by the reference area of the model and the steady velocity pressure, as shown in the following equations: Steady wind force (1) 0.5ρU 2 A Peak wind force Peak wind force coefficient = (2) 0.5 ρU 2 A where ρ is air density, U is the steady wind velocity and A is the projection area of the Steady wind force coefficient = elliptic cylinder. The steady drag coefficient increased with an increase in the angle of attack. The peak drag coefficient was the same as the steady drag coefficient at angles of attack less than 30 degrees, and was larger at angles of attack from 50 degrees to 80 degrees. The steady lateral force coefficient was largest at a 10-degree angle of attack. The peak lateral force coefficient was larger at angles of attack from 25 degrees to 50 degrees. Figure 8 shows the -0.4 0.0 0.4 0.8 1.2 Time (sec) tr = 0.2 sec tr = 1.4 sec 1.6 -0.4 0.0 a) Wind velocity 0.4 0.8 1.2 Time (sec) 0.5 0.4 0.3 0.2 0.1 0.0 -0.1 tr = 0.2 sec tr = 1.4 sec Lateral force (N) tr = 0.2 sec tr = 1.4 sec 0.5 0.4 0.3 0.2 0.1 0.0 -0.1 Drag (N) Wind velocity (m/s) 2.5 2.0 1.5 1.0 0.5 0.0 -0.5 -0.4 0.0 1.6 b) Drag 0.4 0.8 1.2 Time (sec) 1.6 c) Lateral force Figure 4: Time history data of Small-cylinder (D =70mm, α =45deg., Ut =2.0m/s, tr =0.2sec, 1.4sec, with dummy) -0.4 0.0 0.4 0.8 1.2 Time (sec) tr = 0.2 sec tr = 1.4 sec 1.6 -0.4 0.0 a) Wind velocity 0.4 0.8 1.2 Time (sec) 1.6 0.5 0.4 0.3 0.2 0.1 0.0 -0.1 Lateral force (N) tr = 0.2 sec tr = 1.4 sec 0.5 0.4 0.3 0.2 0.1 0.0 -0.1 Drag (N) Wind velocity (m/s) 2.5 2.0 1.5 1.0 0.5 0.0 -0.5 tr = 0.2 sec tr = 1.4 sec -0.4 b) Drag 0.0 0.4 0.8 1.2 Time (sec) c) Lateral force Figure 5: Time history data of Large-cylinder (D =140mm, α =45deg., Ut =2.0m/s, tr =0.2sec, 1.4sec, with dummy) Drag (N) With dummy model Without dummy model -0.4 0.0 0.4 0.8 1.2 Time (sec) a) Drag 1.6 0.5 0.4 0.3 0.2 0.1 0.0 -0.1 With dummy model Without dummy model Lateral force (N) 0.5 0.4 0.3 0.2 0.1 0.0 -0.1 -0.4 0.0 0.4 0.8 1.2 Time (sec) 1.6 b) Lateral force Figure 6: Comparison between time history of wind forces with dummy and without (D =140mm, α =45deg., Ut =2.0m/s, tr =0.2sec) 1.6 The Seventh Asia-Pacific Conference on Wind Engineering, November 8-12, 2009, Taipei, Taiwan relationship between the overshoot coefficient and the angle of attack and clarifies the effect of the angle of attack on the overshoot phenomenon. The overshoot coefficient is defined as the ratio of the peak value to the steady value, as shown in the following equation: Overshoot coefficient = Peak wind force Steady wind force (3) where the peak wind force is defined as the actual maximum force, and the steady wind force is defined as a mean of 20 seconds in a steady flow. The overshoot of drag was remarkable at angles of attack larger than 50 degrees. The overshoot of lateral force was remarkable at angles of attack from 30 degrees to 80 degrees. Effects of the section-size of elliptic cylinder on the overshoot phenomenon Figure 9 shows the relationship between the target wind velocity and the overshoot coefficients of wind force on Small-cylinder and Large-cylinder at a 45-degree angle of attack in the case of the rise time of 0.2sec. The overshoot coefficients of Small-cylinder and Largecylinder increased with a decrease in the target wind velocity. The overshoot coefficient of Large-cylinder was larger than that of Small-cylinder for all target wind velocities. We reported that the overshoot phenomenon depended on the non-dimensional rise time, tr', defined by Eq. (4) in [Kawashita and Maeda et al. (2008) and Takeuchi and Maeda et al. (2008)]: t r′ = U t ⋅ t r / d (4) in which Ut is the target wind velocity, tr is the rise time and d is a reference length (here, the length of the flow direction of the elliptic cylinder). Figure 10 shows the relationship between the non-dimensional rise time and the overshoot coefficient. The overshoot coefficients ordered according to the non-dimensional rise time are plotted on one line regardless of the Steady drag coefficient Peak drag coefficient 2.5 2.0 Peak not found 1.5 Steady lateral force coefficient Peak lateral force coefficient 3.0 Wind force coefficient Wind force coefficient 3.0 1.0 0.5 0.0 -0.5 2.5 Peak not found 2.0 1.5 1.0 0.5 0.0 -0.5 0 10 20 30 40 50 60 70 Angle of attack, α (degree) 80 90 0 10 20 30 40 50 60 70 Angle of attack, α (degree) 80 90 80 90 b) Lateral force a) Drag Figure 7: Wind force coefficient with angle of attack (D =140mm, Ut =5.0m/s, tr =0.2sec, with dummy) 5 Overshoot coefficient Overshoot coefficient 5 4 3 Peak not found 2 1 0 4 3 2 1 Peak not found 0 0 10 20 30 40 50 60 70 Angle of attack, α (degree) a) Drag 80 90 0 10 20 30 40 50 60 70 Angle of attack, α (degree) b) Lateral force Figure 8: Overshoot coefficient with angle of attack (D =140mm, α =45deg., Ut =5.0m/s, tr =0.2sec, with dummy) The Seventh Asia-Pacific Conference on Wind Engineering, November 8-12, 2009, Taipei, Taiwan target wind velocity. The figure indicates a strong relationship between the overshoot phenomenon and the non-dimensional rise time. The overshoot coefficient tended to increase with a decrease in the non-dimensional rise time regardless of the section-size of cylinder. The overshoot coefficients are plotted on one line regardless of the section-size of cylinder and are explained by the non-dimensional rise time. The overshoot phenomenon may be strongly affected by the non-dimensional rise time. Numerical Investigation The flow around Small-cylinder subjected to step-function-like gusts was simulated using CFD software (CDAJ products) using a k-ε model, as shown in Figure 11a. A twodimensional incompressible flow was assumed and the differential equations were discretized by the finite volume method. The MARS scheme was applied to the convective term. The 8 Small-cylinder Large-cylinder 6 Overshoot coefficient Overshoot coefficient 8 4 2 Small-cylinder Large-cylinder 6 4 2 0 0 2.0 3.5 5.0 Target wind velocity (m/s) 2.0 7.0 3.5 5.0 Target wind velocity (m/s) 7.0 b) Lateral force a) Drag Figure 9: Overshoot coefficient with target wind velocity (D =140mm, α =45deg., tr =0.2sec, with dummy) 10 U t = 2.0m/s Overshoot coefficient Overshoot coefficient 10 1 Small-cylinder experiment (D = 70 mm) Large-cylinder experiment (D = 140 mm) 1 1 10 100 Non-dimensional rise time U t = 5.0m/s U t = 7.0m/s U t = 2.0m/s U t = 3.5m/s U t = 5.0m/s U t = 7.0m/s 1 10 100 Non-dimensional rise time b) Lateral force a) Drag U t = 3.5m/s c) Graph legend Figure 10: Relationships between overshoot coefficient and non-dimensional rise time (α =45deg., with dummy) y (mm) 1500 U (m/s) Nesting1 Nesting3 Ut Nesting2 Nesting4 1050 U = −2 Ut 3 U t + 3 2t t 2 t r3 tr 450 0 Inflow boundary x (mm) 0 700 13000 2000 a) Analytic region tr b) Inflow wind velocity Figure 11: Analytic region and inflow wind velocity t (sec) The Seventh Asia-Pacific Conference on Wind Engineering, November 8-12, 2009, Taipei, Taiwan -0.4 0.0 0.4 0.8 1.2 Time (sec) Experimental Calculated 1.6 -0.4 a) Wind velocity 0.0 0.4 0.8 1.2 Time (sec) 0.5 0.4 0.3 0.2 0.1 0.0 -0.1 Lateral force (N) Experimental Calculated 0.5 0.4 0.3 0.2 0.1 0.0 -0.1 Drag (N) Wind velocity (m/s) 2.5 2.0 1.5 1.0 0.5 0.0 -0.5 1.6 b) Drag Experimental Calculated -0.4 0.0 0.4 0.8 1.2 Time (sec) 1.6 c) Lateral force Figure 12: Comparisons between calculated and experimental results with dummy (D =70mm, α =45deg. , Ut =2.0m/s, tr =0.2sec) 10 Small-cylinder with dummy experiment (D = 70 mm) Large-cylinder with dummy experiment (D = 140 mm) Overshoot coefficient Overshoot coefficient 10 1 1 1 10 100 Non-dimensional rise time a) Drag 1 10 100 Non-dimensional rise time b) Lateral force CFD Simulation (D = 70 mm) c) Graph legend Figure 13: Relationships between overshoot coefficient and non-dimensional rise time (α =45deg.) pressure and velocities were coupled by means of the PISO method. The standard k-ε model was used for the turbulence model. The short-rise-time wind velocity measured in the wind tunnel test was modeled on an approaching flow in the CFD simulation, as shown in Figure 11b. Figure 12 shows the comparisons between the CFD calculation and experimental results in the case of a target wind velocity of 2.0m/s and a rising time of 0.2 seconds, as an example. In this case, we confirmed the overshoot phenomenon in the simulation as well as in the experimental results, but in some cases, we were not able to simulate the overshoot phenomenon in the experimental results. The steady wind forces of the simulation corresponded to those of the experimental results in every case. The simulated relationship between the non-dimensional rise time and the overshoot coefficient in Figure 13 agreed well with the experimental relationship in that the overshoot phenomenon of lateral force was more remarkable than of drag. Conclusions The characteristics of aerodynamic force on an elliptic cylinder subjected to short-risetime gusts were investigated by using a specially equipped wind gust tunnel. We investigated mainly the effects of section-size of an elliptic cylinder on overshoot phenomena of wind forces using two elliptic cylinders of axial ratio 2:1, 70mm and 140mm in major axis. This revealed the features of the transient wind forces acting on an elliptic cylinder subjected to a short-time-rising gust as follows: • We confirmed a remarkable overshoot phenomenon generating a much larger wind force than in a steady flow. • The overshoot of drag was little affected by the installation of a dummy to realize a two-dimensional flow around the specimens. On the other hand, the overshoot of lateral force with the dummy was larger than that without the dummy. • For the elliptic cylinder subjected to the gust of the target wind velocity of 5.0m/s and the rise time of 0.2sec, the overshoot of drag was remarkable at angles of attack larger The Seventh Asia-Pacific Conference on Wind Engineering, November 8-12, 2009, Taipei, Taiwan • • • than 50 degrees, and the overshoot of lateral force was remarkable at angles of attack from 30 degrees to 80 degrees. The overshoot coefficient, which is defined as the ratio of the maximum to the steadystate value of wind force, was determined by a non-dimensional rise time, which is defined by rise time, gust speed and section-size, regardless of the rise time and speed of gust. The overshoot coefficient increased with a decrease in the non-dimensional rise time regardless of section-size of cylinder. And the overshoot coefficients followed one pattern regardless of the section-size of cylinder and were explained by the nondimensional rise time. Using CFD we estimated the flow around the cylinder subjected to step-function-like gusts and confirmed that the overshoot phenomenon was the same as in the experimental results. Finally, it was indicated that the CFD-analyzed relationship between the overshoot coefficient and the non-dimensional rise time agreed well with the experimental results. Acknowledgment This study was supported by Grant-in-Aid for JSPS Fellows 21-4043, Japan Society for the Promotion of Science. References Kawashita, H., Maeda, J., Hayata, T., Takeuchi, T., Nakamura, S. and Tsuru, N. (2008), “The Characteristics of Wind Forces on an Elliptic Cylinder subjected to Step-function-like Gusts”, Proceedings of the 20th National Symposium on Wind Engineering, Tokyo, Japan, December 3~5, 325-330 (in Japanese) . Keulegan, G. H. and Carpenter, L. H. (1958), “Forces on Cylinders and Plates in an Oscillating Fluid”, Journal of Research of the National Bureau of Standards, Vol. 60, No. 5, 423-440. Matsumoto, M., Shimamura, M., Maeda, T., Shirato, H., Yagi, T., Hori, K., Kawashima, Y. and Hashimoto, M. (2007), “Drag forces on 2-D cylinders due to sudden increase of wind velocity”, Proc. of 12th International Conference on Wind Engineering, 1727-1734. Morison, J. P., O’Brien, M. P., Johnson, J. W. and Schaaf, S. A. (1950), “The Force Exerted by Surface Waves on Piles”, Petroleum Trans. AIME, 189, 149-154. Nakamura, S., Maeda, J., Kawashita, H., Takeuchi, T., Hayata, T. and Tsuru, N. (2008), “Unsteady Pressure on a Gable Roof Body Subjected to a Short Time Rise Gust”, Proceedings of the 20th National Symposium on Wind Engineering, Tokyo, Japan, December 3~5, 337-342 (in Japanese). Nomura, T., Kitamura, N. and Kitagawa, T. (2000), “Characteristics of unsteady drag on a square cylinder under sudden change of wind speed”, Proceedings of Computational Wind Engineering, 3-6. Sarpkaya, T. (1966), “Separated Flow about Lifting Bodies and Impulsive Flow about Cylinders”, AIAA Journal, Vol. 4, No. 3, 414-420. Shiraishi, N., Matsumoto, M. and Shirato, H. (1982), “A Fundamental Study about Unsteady Aerodynamic Characteristics of Structures due to Fluctuating Wind”, Proceedings of JSCE, No.328, 19-30 (in Japanese). Takeuchi, T., Maeda, J., and Kawashita, H. (2008), “The Overshoot of Aerodynamic Forces on a Railcar-like Body under Step-function-like Gusty Winds”, Sixth International Colloquium on Bluff Body Aerodynamics & Applications, Milano, July 20~24, 193-196. Taneda, S. (1972), “The Development of the Lift of an Impulsively Started Elliptic Cylinder at Incidence”, Journal of the Physical Society of Japan, Vol. 33, No. 6, 1706-1711. Yamaguchi, Y., Kashiwagi, Y. and Ozono, S. (2008), “Aerodynamic characteristics of vehicles under quicklyaccelerated winds”, Proceedings of the 20th National Symposium on Wind Engineering, Tokyo, Japan, December 3~5, 349-354 (in Japanese).
© Copyright 2026 Paperzz