A E R O D Y N A M I C S A N D D Y N A M I C S O F B L U F F B O D I E S IN P R E S E N C E OF T H E M O V I N G SURFACE B O U N D A R Y - L A Y E R C O N T R O L S A N D E E P R. M U N S H I B.E., Maharaja Sayajirao University, India, 1987 M.Sc.(Engg.), Indian Institute of Science, Bangalore, India, 1990 A THESIS S U B M I T T E D IN P A R T I A L F U L F I L M E N T O F T H E REQUIREMENTS FOR T H E DEGREE OF D O C T O R OF PHILOSOPHY in The Faculty of Graduate Studies Department of Mechanical Engineering We accept this thesis as conforming to the required standard T H E U N I V E R S I T Y O F BRITISH C O L U M B I A January 1996 © Sandeep R. Munshi, 1996 in presenting degree this thesis in at the University of partial fulfilment of the requirements of department this or thesis by for scholarly his publication of this thesis or her purposes may be representatives. It is M££riAiilCRL The University of British Columbia Vancouver, Canada Date DE-6 (2/88) 3ANUAR 7 ^1 advanced permission for extensive granted by the head of understood that copying my or for financial gain shall not be allowed without my written permission. Department of an British Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that copying for ^1^6, (yXhlfrgRXN<k ABSTRACT Moving Surface Bpundaryvlayer Control (MSBC) was applied to several two dimensional bluff bodies using a high speed rotating cylinder as a momentum injecting device. Flow past a symmetric airfoil; a D-section; as well as square and rectangular prisms, representing a family of shapes with progressively increasing bluffness were studied in presence of the M S B C . In the case of the airfoil, the leading edge was replaced by a rotating cylinder; while the cylindrical element formed the top and bottom upstream corners of the D-section, square and rectangular prisms. Extensive wind tunnel investigation gave data about the effect of system parameters like rate of the momentum injection, angle of attack and the surface condition of the cylinder on steady and fluctuating components of the pressure distribution around'the body, vortex shedding frequency (Strouhal number), and the lift and drag coefficients. A gain in the Strouhal number with increasing momentum injection suggest a decrease in the effective bluffness of the body. A significant reduction in the drag (up to 80%) was observed for the prisms at a maximum rate of momentum injection, UcjU = 4 (Uc = cylinder surface speed, U — freestream wind speed). In the case of the airfoil, the lift coefficient increased by 160% and the stall angle was delayed from 10° to 48°. A rough criterion in terms of the location of the stagna1 tion point was established to help decide the reversal in the direction of momentum injection as a function of angle of attack to ensure continued benefit. Effect of momentum injection in suppressing the vortex resonance and galloping type of instabilities were studied by mounting the bluff.prism models on a dynamic-test rig inside the wind tunnel test section. The measurement of amplitude and frequency of the transverse oscillations over a range of wind speeds showed complete vibration suppression for momentum injection rates Uc/U < 1. Asymmetric momentum injection (e.g. top cylinder rotating, bottom cylinder stationary), was also found to be effective in disrupting the vortex shedding process and thereby inhibiting vibrations. The suppression of galloping instability in presence of the M S B C was also predicted by the quasisteady analysis. A numerical panel method was developed to simulate bluff body fluid dynam- n ics in presence of the M S B C . The body is descretized into a large number of panels (100 - 150) with each panel comprising of a continuous distribution of linearly varying vorticity and a constant source strength. A set of linear algebraic equations approximates the Fredholm type integral equation derived from ideal fluid flow assumption. The wake is modelled by upper and lower 'free vortex layers' emanating from the separation points on the body. Vorticity is allowed to be shed and dissipated as it is convected downstream along the panels on the 'free vortex layers'. A n analytical expression relates the point vortex modelling a.rotating cylinder to the rate of momentum injection. The panel method is capable of treating multielement configurations (e.g. a rotating cylinder and the truncated airfoil). A n iterative scheme based on the convergence of the wake shape is used to obtain the final solution. The numerically obtained pressure distribution, the lift and drag coefficients agree well with the experimental results. Flow visualization studies in a water channel were performed to obtain better physical insight into the M S B C process. Plexiglas models with rotating cylinders in conjunction with a fine suspension of polyvinyl chloride particles and slit lighting were used to visualize the streaklines. The still photographs and video movie recorded, rather dramatically, the effectiveness of the M S B C in suppressing separation and vortex shedding, making the flow approach the potential character. Overall, the present research firmly establishes potential of the M S B C as a versatile tool for lift augmentation, drag reduction and vibration suppression of several bluff bodies encountered in industrial engineering practice. 111 T A B L E OF C O N T E N T S ABSTRACT . . ii T A B L E OF C O N T E N T S . iv LIST O F S Y M B O L S viii . LIST O F F I G U R E S . . . . . . . xiii LIST O F T A B L E S xxi ACKNOWLEDGEMENT . . xxii DEDICATION 1. INTRODUCTION . . . . . . . . . . . . . 1.1 Preliminary. Remarks 1.2 Physics of Moving Surface Boundary-layer Control 1.3 A Brief Review of the Relevant Literature 1'.4 2. xxiii 1 1 . . . . . . 16 . 1.3.1 Dynamics and fluid dynamics of bluff bodies 1.3.2 Vibration and flow control devices 1.3.3 Numerical simulation of bluff body 17 17 . . :. . flows 24 26 Scope of the Present Investigation 41 E X P E R I M E N T A LINVESTIGATION AND N U M E R I C A L MODELLING. . . / . . 44 2.1 Preliminary Remarks 44 2.2. Experimental Investigations 44 2.2.1 Bluff body shapes 44 2.2.2 Wind tunnel static experiments iv . 53 2.3 3. 2.2.3 Static test parameters . 53, 2.2.4 Dynamical experiments in the wind tunnel 59 2.2.5 Dynamic test parameters 64 2.2.6 Flow visualization 64 Numerical Analysis . 66 2.3.1 Formulation of the problem . 66 2.3.2 Single element configuration . . . : . . 70 2.3.3 Multi-element configuration . . . . . ;. . . . . . . . 78 2.3.4 Computer implementation . 80 A E R O D Y N A M I C S O F A N AIRFOIL W I T H M O M E N T U M INJECTION . 82 3.1 Preliminary Remarks 3.2 Wind Tunnel Investigation . . . \ 3.3 4. . . . . . . . . . . . . . . . . . . . . . . . 82 . . . . . . . . . . . . . 82 3.2.1 Lift and drag characteristics 82 3.2.2 Pressure distribution 92 Numerical Simulation '. 3.3.1 Pressure distribution 3.3.2 Comparison of the various numerical schemes . . . . . 3.3.3 Prediction of the lift coefficient . . . . . . . . . . . . 3.3.4 Iteration characteristic and the wake shape 3.4 Flow Visualization 3.5 Summary 94 94 97 ... ,97 . . . . . . 103 . . 107-.' 110 F L O W PAST A D-SECTION W I T H M O M E N T U M INJECTION 4.1 . . Preliminary Remarks ... . . . . . . . . . . . . . .; . . . . . . v . ...... . . . . Ill Ill 4.2 5. Wind Tunnel Investigation . . 113 4.2.1 Pressure distribution 113 4.2.2 Wake pressure and vortex shedding 4.2.3 Aerodynamic coefficients . 118 123 4.3 Numerical Simulation 124 4.4 Effect of Surface Condition of a Rotating Cylinder 4.5 Flow Visualization 136 4.6 Summary 139 . ...... . . F L O W P A S T R E C T A N G U L A R PRISMS W I T H M O M E N T U M INJECTION .141 5.1 Preliminary Remarks 5.2 Wind Tunnel Investigation . . . . . 5.3 5.4 o.o 6. 131 141 5.2.1 Pressure distribution 5.2.2 Aerodynamic coefficients 5.2.3 Wake pressure and Strouhal number . . . . . 142 142 Universal Strouhal Number . . . . . . . . . . . . . . . . . 152 156 .• 160 5.3.1 Correlation between some nondimensional numbers . . . 160 5.3.2 Experimental confirmation 164 Flow Visualization . bummary 164 . , . . . . . . . . . • • 165 B L U F F B O D Y D Y N A M I C S IN P R E S E N C E O F M O M E N T U M INJECTION . . . . . . . 171 6.1 Preliminary Remarks 171 . 6.2 Control of Vortex Resonance Type Vibrations 172 6.2.1 Mathematical description vi . 172 6.2.2 6.3 Suppression of the vortex shedding . . . Control of Galloping Type of Instability 6.4 178 6.3.1 Criterion for determining the galloping instability . . . . 179 6.3.2 Instability prediction based on static aerodynamic tests Confirmation Through Dynamic Experiments . . . 7. ' 184 188 6.4.1 Dynamical respose of the square prism 189 6.4.2 Dynamic response of the D-section and the rectangular prism 6.5 176 195 Summary 199 CLOSING C O M M E N T S : 200 7.1 Summary of Conclusions . . . . . . . . . . . . . . . . . . . 200 7.2 Suggestions for Future Research 202 REFERENCES . . . . . . . . . 205 APPENDICES A . I N S T R U M E N T A T I O N USED IN T H E S T U D Y B ..226 SINGULARITY E L E M E N T S AND INFLUENCE COEFFICIENTS C POWER CONSUMPTION FOR T H E MSBC • D 228 • i • DYNAMIC SYSTEM PARAMETERS vii . . . . . . . . . 235 • . . . . . . . . . ". . 238- LIST O F S Y M B O L S AR aspect ratio of the bluff body, Aji, Bji,Fji h/hi influence coefficients Cp mean pressure coefficient, (p — p )/(l/2)pV CPB mean base pressure coefficient Cpj mean pressure coefficient at the mid-point of the jth panel C D drag coefficient, C L lift coefficient, Cy lateral force coefficient, Fy / (l / 2) pDVi CJJCD lift to drag ratio \Cjy\ absolute value of the drag coefficient \C-L\ absolute value of the lift coefficient C airfoil chord D characteristic dimension of the bluff body . F force Fr> drag force FL lift force Fy lateral force, (Ft cos a + FD sin a) Hinner total pressure inside the wake Houter total pressure outside the wake K' proportionality constant M total number of panels descretizing each 'free vortex line' 2 QO F /(l/2)pDU 2 D F /(l/2)pDU 2 L 2 , . , , viii • N total number of panels descretizing the body for single element airfoil N R . cylinder rotational speed in rpm Pi coefficients of the polynomial expansion;"i=l,3,5,7 R characteristic radius of the bluff body, D/2 Re Reynolds number, UD/u S circumferential length along the contour of the bluff body St Strouhal number, St,p Strouhal number based on the projected height normal to the flow, fyh/U Sr wake Strouhal number, ffi'/Uj, Su, SL lengths along the upper and lower 'free vortex layers', respectively U freestream fluid speed fyD/U U mean velocity at the edge of the boundary-layer at separation U rotating cylinder surface speed Uc/U momentum injection parameter V nondimensional wind speed, U/u) h B c n —#• V velocity at any point * freestream fluid velocity Vi relative wind velocity V critical wind speed necessary to initiate galloping X x-coordinate direction x/c. nondimensional length along the chord of the airfoil Y nondimensional transverse amplitude, y/oj h\ a longitudinal spacing between two successive vortices in a row c n IX ' '"" ' : b lateral spacing between the two rows of vortices b/a vortex spacing ratio d diameter of the rotating cylinder / frequency f frequency of vortex shedding f" structure natural frequency h projected height of the bluff body normal to the flow h' distnace between two shear layers enclosing the wake hi height of the bluff body k stiffness ki base pressure parameter, \ / l — Cps 7i. length of the bluff body m structure mass per unit length n mass parameter, ft vector denoting normal direction p pressure at any point in the flow Poo freestream pressure v n qj fluid ph l\/2m 2 speed at the jth panel dinner fluid speed inside the wake pouter fluid speed outside the wake r radius of the rotating cylinder, d/2 t time ut translational speed of the vortex street x u'i, u'j fluctuating velocity in the i u, w x, y components of the velocity due to a line vortex c th c < u'jU - > 1 and ] i h directions respectively; i,j = 1,2,3 Reynolds stress u surface velocity of the rotating cylinder Uj,Wj perturbation velocity components at the j to circular natural frequency, 2irf t h u) structure circular natural frequency, 2 x / UJ . circular frequency of vortex shedding, 2irf x x-coordinate direction y y-coordinate direction yi initial vibration amplitude n V panel = y/k/.m n v ; vibration amplitude after fi cycles T circulation strength of the point vortex Tc circulation strength of thepoint vortex AH decrease in the total pressure in the wake with respect to that at infinity %ACD % change in. the drag coefficient %ACx % change in the lift coefficient . angular velocity of the rotating cylinder a angle of attack /3 nondimensional damping coefficient, r i / 2 m w „ 7 vortex strength 7i vortex strength at the i 6 fraction of the orginal circulation surviving the vortex formation ih , panel xi 6± logarithmic decrement for the decaying vibrations, l/[(n' — C system damping ratio C.f viscous damping ratio Ct total effective damping ratio 77 effciency of momentum injection A vortex dissipation factor v fluid viscosity p • fluid density a source strength T nondimensionalized time, u) t 4> potential function for the ideal, incompressible fluid <$>i generalized potential function n Xll l)ln(yi/y )]. n LIST O F F I G U R E S 1-1 • Several examples of bluff body flows: (a) normal flat plate with a wide wake; (b) square section showing' separating shear layers, reattachment bubble, wake and vortex shedding; (c) shear-layer reattachment on a long rectangular section and effect of aspect ratio on the drag. 1-2 3 Schematic diagrams showing principles of: (a) vortex resonance; and (b) galloping instabilities. Here: St, Strouhal number; V', freestream velocity; V , resultant velocity; fv, frequency of vortex shedding; fn, natural frer quency of the system; F(a), force on the prism due to motion; y, transverse displacement during vibration. 1-3 ...... . . . . 5 Some examples of the devices used for the suppression of wind-induced oscillations 1-4 9 Several applications of a rotating cylinder to an aircraft wing: (a) boundarylayer control function of a rotating cylinder flap; (b) high-lift leading-edge rotating cylinder in combination with a droop nose; (c)'high-lift rotating cylinder in combination with a slat 1-5 12 The North American Rockwell OV-10A aircraft in flight demonstrating a successful application of the rotating cylinder as a high-lift device. 1-6 13 Schematic flow illustrating development of circulation around a rotating cylinder and the associated vortex shedding into the wake. The circulation around the circuit Cw embracing the wake but not the body must have a positive value r e - Therefore, the circulation around a'circuit CB embracing the body but not the wake has a negative value -Vc, as required by the fact that circulation around the closed curve C must be zero.. 1-7 Velocity profiles in the neighbourhood of separation: (a) waif moving downstream; (b) wall moving upstream [65] 1-8 19 A bluff body with a rotating cylinder: (a) streamline patterns, (b) pressure distribution; i) 1-9 18 , no rotation, ii) , with rotation. : 20 Applications of the Moving Surface Boundary-Layer Control (MSBC).pro-, xiii - cedure:, (a) airfoil type geometries; (b) bluff bodies such as a flat plate at a large angle of attack; rectangular prisms; and tractor-trailer truck configurations. Boundary-layer control and the associated drag reduction through fence type tripping device is also indicated 27 The M S B C concept presents an exciting possibility of applications to next generation of civil engineering structures such as bridge-towers and. supertall buildings. . . . . . . . . . . . . . , . . . . : . . ; . A schematic diagram showing the scope of the investigation. . . . 29 43 Bluff body geometries, without and with momentum injection, considered for investigation: (a) schematic diagrams; (bi) photograph of the Joukowski airfoil model; (b2) photograph showing details of the pressure taps and pressure conducting tubings; (c) D-section and rectangular prism; (d) square prism and the dynamical test model. \ . . . . . . . . 46 Detailed schematic of the rotating cylinder and drive mechanism. . 51 Schematic diagram showing details of the pressure taps near the rotating element. . . . . . . . . . . . . . . . . 52 A schematic diagram of the low, speed, low turbulence, closed circuit wind tunnel used in the test-program. 54 Experimental set-up for static aerodynamic investigation. . . . . 55 Constructional details.and calibration plots of the disc probe. . . . 57 Instrumentation layout for vortex shedding frequency and fluctuating pressure measurement. . . . . . . . . . . . . . . . 58 Rotating cylinders with three different types of surface condition: (a) schematic diagram; (b) photographs showing details 60 A schematic diagram, of the dynamic experiment set-up involving the wind tunnel, test-rig and data acquisition system. . . . . . . . . . . 62 A three dimensional view of the light weight balsa wood model utilised for dynamic experiments. 63 A schematic diagram of the closed circuit water channel facility used in the flow visualization study. Slit lighting was used to minimize distortion due to three dimensional character of the flow. Long exposure provided pathxiv lines with ploy vinyl chloride particles serving as tracers. The dimensions are in mm • 65 An arbirtrary two dimensional body. . . . . . . . . 69 Numerical panel formulation describing, (a) single element configuration and, (b) multi-element configuration. 72 Numerical modelling of a rotating cylinder using a line vortex. Numerical algorithm for the panel method. . . . . . . . . . . . . . . . 77 81 Experimentally obtained lift coefficient for the Joukowski airfoil in presence of momentum injection: (a) smooth leading-edge rotating cylinder; (b) rough cylinder; (c) splined cylinder. . . . . . . . . . . . . 84 Experimentally obtained drag coefficient for the Joukowski airfoil in presence of momentum injection: (a) smooth leading-edge rotating cylinder; (b) rough cylinder; (c) splined cylinder. 87 Experimentally obtained lift/drag ratio of the Joukowski. airfoil in presence of momentum injection: (a) smooth cylinder; (b) rough cylinder; (c) splined cylinder. . . . . 88 Variation of CJJ with CD for the Joukowski airfoil in presence of momentum injection: (a) smooth cylinder; (b) rough cylinder; (c) splined cylinder. 8 A comparative study showing the effect of cylinder surface during momentum injection on the variation of: (a). C£ with a; (b) CL/CD with a 90 A comparative study showing the effect of surface roughness and momentum injection on the lift coefficient at high angles of attack. . . . 91 Effect of momentum injection with the splined surface cylinder on the variation of CD with a. The smooth cylinder results are also presented to serve as reference ". . . . . . . . . . 93 Experimentally obtained pressure distribution for the Joukowski airfoil with the M S B C provided by a smooth surface cylinder: (a) a = 5° ; (b) a = 10° ; (c) a = 15°; (d) a = 20° ; (e) a = 30° ; (f) a = 4 0 ° , . . 95 Comparison of the numerically obtained pressure distribution with the xv experimental data for a Joukowski airfoil with the M S B C : (a) a = 5° (b) a = 10°; (c) a = 15° ;'(d) a = 20°. ... . . ' 3-10 98 Comparison between the results obtained through the single and two element numerical panel formulations in presence of the M S B C : (a) a = •: 10° U /U c = 0; (b) a = 10°, U /U c = 1. The experimental data are also included. 3-11 .'. : . . . . . 100 Comparative studies of pressure plots in'presence of the M S B C : (a) attached versus separated flow solutions for a single element formulation; "(b) corresponding results for the two element airfoil! The experimental data are also included. . . . . . . . . . 3-12 101 Relative performance of the attached and separated flow models as applied to the two element airfoil:(a) a = 20°, U /U c = 2; (b) a = 2 0 ° , U /U c 4. 3-13 102 Numerically obtained lift coefficient as compared to the experimental data for a Joukowski airfoil in presence of the M S B C . 3-14 = . 104 Iteration characteristic and convergence of the wake shape for a flow past the Joukowski airfoil in presence of the M S B C . The results were obtained using the two element separated flow numerical model. 3-15 105 Effect of number of panels on the convergence accuracy of the numerical scheme. 3- 16 . . . . . . . . . . . . . . . . . . •: . . . 106 Representative flow visualization pictures showing, rather dramatically, effectiveness of the boundary-layer control through momentum injection. The observed results compare well with the delay in separation predicted by the panel code: (a) a = 15°; (b) a = 45°. '. 108 4- 1 The 38 storey Prince Hotel in Otsu, Japan; height 137 m. . . . . . . 4-2 Surface pressure plots for the D-section showing effects of the angle of attack and momentum injection a = 0 - 3 0 ° , 4 0 ° - 7 5 ° , 9 0 ° - 1 8 0 ° . 4-3 . 114 Effect, of three different modes of momentum injection on the surface pressure distribution at a = 90° 4-4 . 112 • • • 119 Wake pressure and the Strouhal number for the D-section as affected by the momentum injection 122 xvi Drag and lift characteristics of the D-section as affected by the momentum injection . . . . . . . . . . . 125 The effect of three different modes of momentum injection on the drag coefficient of the D-section as a function of angle of attack. . . . . 126 Comparison between the numerically and experimentally obtained pressure plots for the D-sectiori. Note, the numerical scheme accurately predicts presence of large suction peaks, as well as rise in the wake pressure in presence of momentum injection', a •= 0, 10°, 20°, 30°. . . . . .. 128 Numerical simulation data showing the effect of momentum injection: (a) drag coefficient; (b) wake geometry as represented by the 'free vortex lines'. . . . . . . . . . . . . 130 Effect of cylinder surface condition on: (a) drag coefficient; (b) wake pressure 132 Effect of cylinder surface roughness on:(a) lift coefficient; (b) lift/drag ratio. . 133 Effect of cylinder surface characteristic on pressure distribution over the D-section:(a) a = 0; (b)a = 45°. 134 Flow visualization photographs showing effectiveness of the Moving Surface Boundary-layer Control (MSBC) as applied to the D-section. Note, at Uc/U = 4, the fluid motion tends to approach the potential character 137 Effect of momentum injection on the characteristic features of the wake and the shed vortex system, as inferred from the flow visualization study 138 Surface pressure plots for the rectangular prism showing effects of the angle of attack and momentum injection: a = 0, 15°, .30°, 45,° 7 5 ° , 90°, 105°, 135°, 180°. 143 Surface pressure plots for the square prism as affected by the angle of attack and momentum injection: a = 0, 20°, 45°, 75°, 90°, 105°. . 147 Effect of reversing the direction of upstream cylinder on the surface pressure distribution at a — 90°: (a) rectangular prism; (b) square xvn prism 5-4 . 149 Drag and lift characteristics of the rectangular prism as affected'by the momentum injection. 5-5 154 The drag and lift characteristics of the square prism as affected by momentum injection . . 5-6 155 Effect of momentum injection on the wake pressure and Strouhal number of the rectangular prism. . . . . ' 5-7 . . . Variation in the wake pressure and Strouhal number for the square prism as affected by the momentum injection. 5-8 158 . . 159 Effect of the angle of attack and bluff body geometries on the universal Strouhal number Su: (a) no momentum injection; (b) with momentum injection 5-9 163 Effect of momentum injection on the characteristic features of the wake and the shed vortex system associated with the rectangular prism ( A R = 0.5) as observed during the flow visualization study 5-10 166 Flow visualization pictures showing effect of the momentum injection on the wake of the square prism at a = 15°. . . . . . . . . . . . . . . . 167 5-11 Flow visualization pictures for a square prism at a = 90°. 5- 12 • A schematic diagram showing possible strategies of momentum injection for a rectangular prism undergoing a change in the angle of attack. 6- 1 168 169 Schematic diagram showing a bluff structure experiencing vortex excitation: (a) physical model; (b) variation of vortex shedding frequency and resonance; (c) phase space representation of the dynamical model. . 6-2 175 Power spectral representation of the wake pressure showing reduction in the strength of the shed vortex in presence of the momentum injection: (a) Uc/U = 0; (b) U /U c 6-3 . . . . 177 A bluff structure experiencing galloping and the criterion for instability. 6-4 = 1, 2. . . . . . . . . . . 180 (a)Influence of the momentum injection on the lateral force coefficient for the D-section. The criterion for galloping instability is dC /da y > 0. Note a marked improvement in the stability with the momentum injection xvm (Uc/U = 4); (b) The rectangular prism is neutrally, stable (dC /da y for a < 20° in absence of momentum injection (Uc/U is significantly improved (dC /da y = 0). The stability < 0) with the cylinder rotation (Uc/U = 2, 4); (c) The square prism is unstable (dC /da > 0) for a < 5° in y absence of momentum injection (Uc/U (Uc/U = 0) = 0). With the cylinder rotation = 2, 4) the stability is restored completely (dC /da y < 0). . 185 Dynamic response of the square prism without and with symmetric momentum injection. Note the presence of both vortex-resonance and galloping in absence of the momentum injection (Uc/U = 0). ; The square prism becomes completely stable in presence of the momentum injection as small as (7c/(7 = 1. . . . 190 Effect of an increase in damping on the dynamic response of the square prism without and with symmetric momentum injection: Note, a significant reduction in the amplitude, particularly in galloping, in presence of the momentum injection. 'As before, the stability is restored for Uc jU — 1 . . . . . . . . . . 192 Effect of asymmetric momentum injection on the dynamic response of the square prism at various wind speeds. Note, it becomes necessary to increase the asymmetric momentum injection to at least Uc/U restore complete stability =1.5 to • 193 Schematic diagrams showing dynamics of bluff bodies in the extended parameter space (Y, Y, V, Uc/U): (a) an.unstable bluff structure displays both the vortex-resonance and galloping in absence of the momentum injection (Uc/U = 0). Stability is completely restored with small amount of momentum injection (Uc/U = 1); (b) increase in the system damping affects the critical wind speed and system amplitude. . . . . . . 194 Effect of direction of the momentum injection on. the galloping instability: (a) for V > Vc, an inherently unstable system can be stabilized by momentum injection in the right direction (b) galloping instability at subcritical wind velocity (V < Vc) under adverse direction of the momentum injection (U JU C < 0).. . . . . . . . . . . xix .196 6-10 Dynamic response of the D-section as affected by the momentum injection. The D-section experiences vortex-resonance, but is stable in galloping.197 6-11 Effect of momentum injection on the galloping instability of a rectangular prism (AR=0.5) at a = 0. Note, the prism experiences vortex resonance, but. is stable in galloping. . . . . . . . . . . . . . . . . .• . ; 198 B-l A n infinitely long vortex filament 233 B-2 i-th panel with constant vortex distribution. B-3 Reference coordinate system.. 234 B-4 i-th panel with linear vortex distribution. 234 C-l Experimentally measured power consumption of a momentum injection . . . . . , . . . rotating cylinder D-l 237 Amplitude vs. time history of the bluff body model (square-section) undergoing free vibrations. D-2 : .233 : 240 Natural frequency of the vibrating bluff body (square-section) as obtained through the F F T analysis of the amplitude vs. time diagram. D-3 Calibration of the displacement transducer. D-4 Calibration of the spring used for the dynamic test set-up. xx . . , . . • 241 , 242 . . . . 243 LIST O F T A B L E S Table 3.1 Comparative performance of the three rotating cylinders in terms of the maximum lift coefficient and the stall angle for Uc/U —3 and 4 85 Table 3.2 Average values of dCi,/dct for three rotating cylinders. Table 3.3 Average values of dCj^/da as affected by the cylinder surface condition and rate of. momentum injection. . . . . . . . . . . . . . . . . . . . 85 86 Table 4.1 Drag reduction with momentum injection Table 5.1 Drag reduction for the rectangular prism with the momentum injection at various a Table 5.2 152 Change in the lift of the rectangular prism with the momentum injection at various a Table 5.3 153 Drag reduction for the square prism with the momentum injection at various a Table 5.4 123 153 Change in the lift of the square prism with the momentum injection at various a. . . . v 153 xxi ACKNOWLEDGEMENT Here Iwish to express sincere appreciation to those who directly or indirectly helped me, while I was involved with the present research. First and foremost, I would like to thank my supervisor, Prof. Vinod J . Modi, for all that he has done for me in last five years. His invaluable insight and technical advise has .been the guiding spirit behind this research. I have also imbibed from him : the art of better organization and presentation. Those friendly discussions oh world affairs provided useful balance to the academic pursuit as a graduate student. I wish to thank Dr. Gautam Bandyopadhyay (IIT Kharagpur, India) for helpful discussions during development of the numerical code when he spent one year as a research associate in our group in 1991-92. I am also grateful to Prof. Toshio Yokomizo for making available the excellent flow visualization facility at the Kanto Gakuin University, Japan. Special thanks go to Mrs. Kyoko Yokomizo for arranging several sightseeing trips which made my stay in Japan an unforgettable experience. Comments of the committee members, Prof. Ian Gartshore, Prof. Sander Calisal and Prof. Michael Quick, are gratefully acknowledged. I wish to extend my appreciation to Mr. Tony Basic, Mr. E d Abel and Mr. John Richards at the department workshop for fabricating the bluff body models as well as the experimental set-up. Of course, I cannot overlook my past and present colleagues and friends, Satyabrata Pradhan, Mae Seto, Anant Grewal, Spiros Kalantzis, Mathieu Caron, Gary Lim, Mark Chu, Shinji Hokamoto, Itzhak Marom, Yuan Chen, Khan Muneer, Zhiyong Ma, Harry Mah, Afzal Suleman, Alfred Ng, Oliver Ying, Simon St.Hill, Andrew Dobric, Faqurrudin Karray, Nagaraja Rudrapatna, and Dipak Ghosh who made my stay at the U B C a memorable one. The investigation reported here was supported by the Natural Sciences and Engineering Research Council of Canada, Grant No. A-2181. xxn Dedicated to my parentsj whose affection and guidance have always motivated rhe xxin 1. I N T R O D U C T I O N 1.1 Preliminary Remarks Fluid dynamics of bluff bodies, of varying shapes, has been a subject of con- siderable study ever since the pioneering contributions by Strouhal [1], Prandtl [2], Karman [3], and others. The concept of the boundary-layer was first introduced by Prandtl in 1904. Ever since, engineers and scientists have tried to minimize its adverse effects and control it to advantage. Delay of the wing stall to high angles of attack, as well as increase in lift and reduction in drag of aircraft have been of considerable importance to the aeronautical industry. Rise in the fuel cost and growing environmental awareness have led to design of energy efficient road vehicles such as cars, buses, trucks, etc. Furthermore, wind effects on bluff structures may lead to vortex resonance and galloping type of instabilities. Flow-induced vibrations of tall buildings, bridges, smoke-stacks, air traffic control towers, transmission lines, etc. have been reported frequently by the engineering industry. Collapse of the Tacoma Narrows Bridge in U.S.A. in 1940, due to torsional instability, is a well known example. Tall buildings have been known to oscillate at low frequency (< 1 Hz) and large amplitude. Recent advances in material science and computer aided design have led to structures with reduced stiffness making them prone to wind, earthquake, as well as ocean waves and current excited oscillations. The twenty first century would witness super tall buildings (> 1000m) and extra long span bridges (> 2000m), as pointed out by Kubo et al. [4,5]. A vast body of literature accumulated over the years has been reviewed by several authors including Wille [6], Marris [7], Morkovin [8], Parkinson [9], Cermak [10], Welt [11,12], and Modi et al. [13]. Streamlined objects such as a thin airfoil aligned with the flow suffer little flow 1 1 separation, have negligible wake and the drag is primarily due to viscous skin friction. In contrast, a bluff body immersed in a fluid stream exhibits flow separation resulting in a wide wake downstream. A large contribution to the drag can be ascribed to the presence of a significant difference in pressure between the freestream and the wake. The wake in general carries vorticity, is turbulent and persists for many diameters downstream of the body. At low Reynolds numbers and depending on the shape as well as orientation'of the body, the shear layers separating from the top and bottom surface roll-up to form alternately shedding vortices resulting in the wellknown Karman vortex street. In the case of a bluff body with sharp corners (e.g. a rectangular prism) location of the boundary-layer separation is fixed, but for a body with curved surface, as in the case of a circular cylinder, the point of separation is sensitive to changes in the Reynolds-number. Figure 1-1 illustrates some typical examples of bluff body-flows. Flow past a circular cylinder [14,15] is a good example illustrating the fluid dynamical phenomena associated with bluff bodies. At a very low Reynolds number (Re < 1), the flow is purely viscous and the wake is absent. For 5 < Re < 15, an attached symmetric pair of Foppl vortices appear in the immediate vicinity of the cylinder. As the Reynolds number is further increased to around 40, there is a laminar separation of the shear layers and periodic vortex shedding initiates. The shed vortices remain laminar up to Re « 150. A transition to turbulence in vortices occurs for 150 < Re < 300. The vortex street becomes fully turbulent as the Reynolds number increases further (300 < Re < 3 x l 0 ) . Around 3 x l 0 5 5 the boundary-layer undergoes laminar to turbulent transition and becomes fully turbulent at Re « 3.5 x 10 . Now 6 the wake narrows significantly, is disorganized, and there is a sudden drop in the total drag force. For Re > 3.5 x 10 the turbulent vortex shedding is re-established and 6 the wake widens with a gradual increase in drag, which reaches a stationary value at some higher Reynolds number. 2 4 Figure 1-1 Several examples of bluff body flows: (a) normal flat plate with a wide wake; (b) square section showing separating shear layers, reattachment bubble, wake and vortex shedding; (c) shear-layer reattachment on a long rectangular section and effect of aspect ratio on the drag. The portion of the bluff body lying inside the wake is referred to as the afterbody. In the normal range of the Reynolds number of engineering interest (10 3 < Re < 10 ), a bluff body experiences significant vortex shedding. This cre8 ates time-dependent pressure distribution on the top and bottom faces of the body resulting in a fluctuating transverse force. If the frequency of the vortex shedding coincides with the natural frequency of the structure it experiences resonance, a forced vibration phenomenon (Figure l-2a). The vortex shedding from a bluff body is mainly a two-dimensional phenomena but the shed vortices are not tw6^in^n^i0nal-.[16]. Hence spanwise effects play an important role in real-life slender structures which may be considered essentially twodimensional. As noted by Bearman [16], oscillation of a bluff body dramatically alters the vortex shedding process. The oscillations provide a means for coupling the flow along the span of the body, which is reflected in a higher spanwise correlation of various physical quantities (e.g. pressure fluctuations). The body motion is also responsible for the capture of the vortex shedding frequency by the body frequency over a certain range of fluid velocity, known as frequency 'lock-in'. In fact, in case of an oscillating square prism, the vortex formation process is comprised of the Karman vortex street as well as a second wake vortex system created by the transverse motion of the body [17]. Thus, the term 'lock-in' is usually applied to the the region where there is a strong magnification of the body-frequency signal at the expense of the Strouhal signal. In-line oscillations have also been observed for structures immersed in a dense fluid (e.g. water). Here the frequency of vibration is twice that of the transverse case. For a circular cylinder, the vortices are first shed in symmetric pairs and the familiar Karman vortex street develops further downstream [16]. Ocean engineering structures (e.g. off-shore oil exploration rigs) are more susceptible to in-line oscillations. As against the vortex resonance, galloping is a different form of instability 4 t STRUCTURE NATURAL FREQUENCY, f 1 Figure 1-2 n Schematic diagrams showing principles of: (a) vortex resonance; and (b) galloping instabilities. Here: St, Strouhal number; V , freestream velocity; V,., resultant velocity; fv, frequency of vortex shedding; fn, natural frequency of the system; F(a), force on the prism due to motion; y, transverse displacement during vibration. 5 referred to as the self-excited oscillations. It is exhibited by only a certain family of bluff body geometries (e.g. a square prism). Here the excitation force depends on the motion itself. The amplitude of oscillation grows until the rate at which the energy is extracted from the fluid stream balances the rate at which it is dissipated by the system (Figure l-2b). In the quasi-steady theory proposed to study the galloping motion [18], the excitation force is approximated through a nonlinear negative damping term in the equation of motion. The theory assumes that the bluff body, at any instant during the galloping cycle, experiences the same lateral force as that when held stationary at the same angle of attack. In galloping, the plunging amplitude could be very large under certain circumstances. A classical example is the galloping of ice-covered transmission cables in an extremely cold weather. A bluff body could also undergo galloping in the torsional mode. In general, experimental results have compared well with prediction of the quasi-steady approach. Another interesting form of self-excited oscillation is flutter, a multidegrees of freedom phenomenon which requires phase difference between the system's generalized coordinates for energy dissipation. One of the important shape parameters of a bluff geometry is the afterbody length. A section with a relatively short afterbody (e.g. a D-section) is generally stable in galloping. Those with somewhat larger afterbodies (e.g. a square and rectangular prisms) would gallop; but the amplitude decreases as the afterbody becomes longer. Stability is restored for bluff bodies with an aspect ratio (length/height ratio, h/h>i) ° f 3 or more [19]. Irrespective of the afterbody shape, all bluff bodies experience vortex-excited oscillations. Furthermore, the aeroelastic behaviour of bluff bodies is significantly affected by the geometry of the corners. In a study by Bokaian and Geeola [20], it was found that the galloping vibrations of a square section were sensitive to the ratio of corner radius to body dimension. A n increase in the rounding radius made the square section less susceptible to galloping instability; and beyond a critical value (0.318) no vibrations could be detected. From the structural dynamics point of view, response of a slender elastic body is determined by the various vibration modes involved in the motion. Furthermore, besides Reynolds number and body shape, the bluff body flows are also influenced by factors such as surface roughness, presence of turbulence in the approaching fluid stream and interference effects due to neighbouring structures. A bluff body which*, behaves as a soft oscillator in a smooth flow becomes less susceptible to galloping instability as the intensity of turbulence is increased [21]. In the same vein, a hard oscillator tends to become relatively soft. Widespread occurrence of wind induced instabilities of structures with the associated problems of comfort and structural damage have forced engineers to predict and suppress these vibrations in advance. A common approach to control wind induced instabilities is to modify the flow field. This has led to the design of helical strakes, perforated shrouds, slats arid similar devices as shown in Figure'1-3. They contribute to aerodynamic damping. Strakes have been frequently installed on smokestacks and ocean based structures, although the associated increase in the drag is a limitation. The effectiveness of helical strakes have been known to deteriorate with an increase in intensity of turbulence [22]. Structural response to wind excitation is also quite sensitive to the damping. This.has led to the design of a variety of energy dissipation devices such as tuned mass dampers, hydraulic dashpots, impact dampers, etc. (Figure 1-3). Approximately 3-4% of the structure weight is necessary for such devices to be effective. Of course, there is a limitation as. to how large the weight could be as the structure grows taller. For example [23], the 60 storey tall John Hancock Tower Building in Boston, U.S.A., requires two 300 tonnes concrete blocks (one each for two orthogonal directions) for effective, vibration control at a cost of US$ 3 million (in 1975). Similar dampers have also been installed on the 110,storey World Trade Center Towers, in New York, U.S.A. In the same category 7 of passive devices belongs a relatively simple concept involving sloshing of a liquid within a closed container. The dissipation of energy takes place through the action of viscous and turbulent stresses associated with wave breaking. The presence of a free surface permits significant wave motion of the sloshing liquid. A n axisymmetric torus shaped container filled with liquid, known as nutation damper, has been studied by Welt and Modi [24,25]'. As the frequency encountered in wind, induced oscillations of a large structure is usually less than 1 Hz, the nutation damper is. ideally suited to this class of problems. Generally a large array of nutation dampers (where each damper is a fraction of a meter in dimension), stacked in a modular fashion, are used in' real-life structures. Nutation dampers have been successfully installed in Japan [26].' ; Apart from the passive control, there is yet another class of active devices requiring a supply of energy for their operation.' Boundary-layer suction or blowing fall in this category. These concepts are more suitable where an appreciable boundarylayer growth occurs. In case of bluff bodies, like a rectangular prism, there is no significant growth of the boundary-layer and the flow separates at the top and bottom . edges of the front face. In this situation, application of the Moving Surface Boundarylayer Control (MSBC) appears quite promising [5]. As the name implies, a portion of the body-surface exposed to the flow is moving and thus imparts tangential velocity to the approaching, fluid. This injection of momentum helps in keeping the flow attached to the surface in the region of adverse pressure gradient. It could result in ; a substantial increase in lift, reduction in drag and assist in suppression of the flowinduced vibrations. Extensive studies of the concepts like boundary-layer suction, blowing, vortex generators,.turbulence promoters, etc. have been carried out over the years. A n indepth review of the subject has been presented by Goldstein [27], Lachmanri [28], Rosenhead [29], Schlichting [30], Chang [31] and others. Despite the specific recommendations by the Associate Committee on Aerodynamics appointed 8 Helical Strakes 2 optimal pitch «5 d p r o t u s i o n = 0.10 - 0.12ds l o c a t i o n : top 33% of height d r a g p e n a l t y : C » 1.3 f o r 8 x 1 0 ' £ Re £ 2 x 1 0 (i.e. 3 times higher i n s u p e r c r i t i c a l range) s d s Slats Shrouds ODD gap : shroud-chimney = 0 . I 2 d open a r e a r a t i o = 20 - 3 6 % l o c a t i o n = t o p 25% o f h e i g h t d r a g p e n a l t y : C * 0.9 f o r 8 x 1 0 * £ Re £ 2 x 1 0 ' (i.e. 2 times higher i n s u p e r c r i t i c a l range) d Internal damping s s l a t width = d /l1.5 gap : s l a t — c h i m n e y =' d / 7 open a r e a r a t i o = -40% d r a g p e n a l t y : C * 1.05 t o I . I f o r 1.5x10* £ Re i 1.5x10 (i.e. e x p e c t e d 2.5 t i m e s h i g h e r i n s u p e r c r i t i c a l range) s s d s ; out of phase v i b r a t i o n s Hydraulic Dashpot Nutation Damper Figure 1-3 Some examples of the devices used for the suppression of windinduced oscillations. 9 by the National Research Council [32], the subject of moving surface boundary-layer control has received relatively little attention. The astute mind of Newton [33] was probably the first to notice the effect of spin on the trajectory of a ball without any inkling to its cause. At least two centuries elapsed before Magnus [34] applied the concept of a rotating cylinder to generate lift, replacing the sails of a ship. Further details on the 'Magnus Lift' concept are available in the reviews by Swanson [35] and Iverson [36]. Prandtl [37] demonstrated his "ship of zero resistance" in 1910, which had two counter-rotating cylinders. Flettner [38] in 1924 applied the principle to the ship 'Bachau' for propulsion. Goldstein [27] demonstrated the use of a rotating cylinder at the leading-edge of a flat plate as a boundary-layer control device. Favre [39] applied the concept to an airfoil by replacing the upper surface with a belt moving over two rollers. The airfoil achieved the maximum lift coefficient of 3.5 with the stall delayed to 55°! No significant research on the concept of the M S B C was reported between 1938-'60. A renewed interest in the subject is indicated through contributions by Alvarez-Calderon [40,41] who reports a successful application of leading-edge rotating cylinder flap to generate high lift for STOL-type aircraft. Flight tests were performed on a single engine high wing research aircraft designed by the Aeronautics Division of the Universidad Nacional de Ingenieria in Lima, Peru. He also suggested possible application of the rotating cylinder for an aircraft wing (Figure 1-4). During the same period Brooks [42] presented his tests on a hydrofoil with a leading or trailing edge rotating cylinders. Between the two configurations, the latter was found to provide a substantial gain in the lift. The primary motivation behind this research was the improvement in fin performance for torpedo control. Steele and Harding [43] applied a rotating cylinder to a ship rudder to increase its manueverability. Force measurements and flow visualization studies were carried out using a water tunnel and a large circulating water channel. Three different configurations of 10 rudder were used. The rotating cylinder : (i) in isolation; (ii) at the leading edge of a rudder; (iii) combined with a flap-rudder, the cylinder being at the leading edge of the flap. From the overall consideration of hydrodynamic performance, mechanical complexity and power consumption, the configuration in (ii) was preferred. A n application to a massive tanker (250,000 tonnes) showed the power requirement for 1 m diameter cylinder rotating at 350 rpm to be around 400 kW^. Flight tests were also performed on the North American Rockwell's OV-10A by NASA's Ames Research Center. The leading-edge of the flap was replaced by a high speed rotating cylinder developed by Alvarez-Calderon [40,41]. The flight test program was designed to assess handling qualities of the propeller-driven S T O L type aircraft at higher lift coefficients. The aircraft was flown at a landing speed of 29-31 m/s, with an approach angle of -8°, corresponding to a lift coefficient of 4.3. The pilot's report suggested that any further decrease in the approach speed were limited due to the lateral-directional stability and control characteristics. Excellent photographs (Figure 1-5) of airplane on ground and in flight have been published in the Aviation Week and Space Technology [44,45]. As evident from the discussion above, the use of a rotating cylinder was so far targeted at specific applications and few attempts were made to study the underlying fluid mechanics involved. A n effort was made by Tennant in 1971 to assess the effectiveness of a rotating cylinder as a boundary-layer control device when applied to a diffuser with a step change in area [46,47]. The diffuser incorporated rotating cylinders to form a part of its wall at the station of the area change. Preliminary experiments were conducted for the area ratio up to 1:2.5, which showed no separation 11 HIGH SPEED POSITION Figure 1-4 AIR FLOW Several applications of a rotating cylinder to an aircraft wing: (a) boundary-layer control function of a rotating cylinder flap; (b) highlift leading-edge rotating cylinder in combination with a droop nose; (c) high-lift rotating cylinder in combination with a slat. 12 CO Figure 1-5 The North American Rockwell OV-10A aircraft in flight demonstrating a successful application of the rotating cylinder as a high-lift device. for appropriate moving surface to diffuser inlet velocity ratio. Tennant et al. [48] have also conducted tests with a wedge shaped flap having a rotating cylinder as the leading edge. Flap deflection was. limited to 15° and the critical cylinder velocity necessary to suppress separation was determined. Effects of increase in the gap-size (between the cylinder and the flap surface) were also assessed. No effort was made to observe the influence of an increase in the cylinder surface velocity beyond Uc/U = cylinder surface velocity, U — free stream velocity). =1.2 (Uc Subsequently, Tennant et al. [49,50] have reported circulation control for a symmetrical airfoil with a rotating cylinder forming its trailing edge. For zero angle of attack, the lift coefficient (Cx) of 1.2 was attained with Uc/U=Z. Of interest is their study concerning boundary- layer growth on moving surfaces accounting for gap effects [51]. Ericsson [52] has also reported use of a rotating cylinder for boundary-layer control on an aircraft wing. Modi et al. [53-60] were the first to generalize the concept of a rotating cylinder as a versatile boundary-layer control device for diverse applications such as high-lift airfoils, drag reduction of static as well as moving bluff bodies, control of wind induced instabilities of civil engineering structures, reduction in snow deposition over a roof, suppression of wave and current induced oscillations of offshore platforms and marine risers, etc. Mokhtarian and Modi [55] have studied, in a comprehensive manner, application of the Moving Surface Boundary-layer Control (MSBC) with reference to two-dimensional Joukowsky airfoils having one or two cylinders acting as momentum injecting elements. A Joukowski airfoil fitted with a rotating cylinder on its upper surface near the leading edge resulted in an increase in the lift coefficient by 200 % with the stall angle delayed to 48° . As a further application of the M S B C concept at a more fundamental level, a flat plate and several rectangular prisms with momentum injection through rotating cylinders were also studied. A drag reduction of about 75 % was achieved for a normal flat plate with twin cylinders, one at each edge (Uc/U = 3). Wind tunnel tests on a tractor-trailer truck model with a rotating cylinder at the top 14 leading edge of the trailer resulted in a significant drag reduction of about 26 %. Flow visualization pictures for the flat plate with the M S B C revealed considerable narrowing of the wake at a cylinder speed of about Uc/U = 3. Even for a Joukowski airfoil at a high angle of attack ( 2 5 ° ) , the momentum injection corresponding to Uc/U — 4 was successful in making the flow essentially potential (no separation). A preliminary analysis of the power requirement for a rotating cylinder suggested that the benefits of drag reduction would easily outweigh the energy demand. This fact is also confirmed in the report by Brown [61] on the application of a leading-edge rotating cylinder flap on an aircraft. The flight tests revealed that the total power required to rotate the cylinder was less than 2% of the propeller power. Increasing air traffic is responsible for the growth in size and speed of airplanes. This, in turn, can only be accomplished by generating higher lift force. Powerful engines producing considerable noise would be required for such applications. Reduction in jet engine noise is difficult to accomplish without paying a price in terms of performance. For communities surrounding an airport, noise-level can be kept to a minimum by employing high-lift take-off and landing at reduced speeds on shorter runways. Moving surface boundary- layer control represents an attractive solution in 1 such a situation, not to mention highly maneuverable military airplanes demanding exacting and superior performance. Every year large amount of goods are transported by the trucking industry using the network of highways.- Even a small reduction in the drag of a tractor-trailer truck could result in a significant saving of the fuel cost. Application of rotating elements for boundary-layer control presents one possibility to that end. In the ever competitive construction industry, the engineers are looking for effective, economical and implementable solutions to fluid-structure interaction problems. Here again, application of M S B C represents a potential candidate. The discussion above touched upon usefulness of a rotating cylinder as an effective boundary-layer control device in the field of applied aerodynamics. Further 15 advances in this technology can only emerge from an understanding of the basic fluid dynamics involved. The' present research aims at exploring M S B C at a fundamental level with particular emphasis on lift augmentation, drag reduction and vibration control of bluff bodies through a comprehensive program involving wind tunnel tests, numerical investigation, and flow visualization. 1.2 Physics of Moving Surface Boundary-layer Control Since, this thesis deals with boundary-layer control using a rotating cylinder, it would be appropriate to outline the basic flow physics involved. A n elegant explanation on the generation of circulation is given by Lighthill [62] as follows (Figure 1-6). Consider, a cylinder, rotating in a clockwise sense, immersed in a fluid flowing from left to right with velocity U. The tangential surface, velocity of the rotating cylinder is Uc- It is known that the strength of a vortex sheet is equal to the jump in the tangential velocity across the sheet. On the upper surface of the cylinder the vortex sheet strength is U — Uc, whereas on the bottom surface it is U + Uc- Thus a more of positive vorticity is shed from the lower surface as compared to the top. The wake receives a net positive amount of vorticity. The Kelvin's theorem on the conservation of the to;tal circulation dictates that the body (i.e the rotating cylinder) in turn receive a net negative vorticity. This imparts a lift to the rotating cylinder in accordance with the Magnus effect. Another effect of rotation is to move the point of separation on the upper surface downstream. The external flow velocity at separation is thus further reduced and so contributes additionally to the reduction in the rate of shedding of vorticity from the upper surface. It is also important to note the presence of a thin boundary-layer on the surface of a rotating cylinder. Vorticity is generated in the boundary-layer around the rotating cylinder. As shown in Figure 1-7, if the wall is moving in the downstream direction, then separation occurs at some height above the wall, and there would be boundary16 layer flow beneath the separated region. The velocity profiles are also shown in the inset. If the wall is moving upstream, separation occurs at the wall and there is a small sublayer of recirculating flow ahead of the separation. The wake behind a rotating cylinder is significantly smaller, the streamlines curve more sharply around the cylinder and the flow approaches an ideal (potential) one. If the rotating cylinder is integral to a bluff body (Figure 1-8), the streamlines bunch together near the rotating element generating a suction pressure field. The separation is delayed and the wake becomes smaller. The pressure in the wake increases. Overall, there is a significant reduction in the drag of a bluff body. The strength of the shed vorticity is reduced considerably and could disappear altogether if the momentum injection is sufficiently strong. The changes produced in the flowfield are proportional to the size and total surface area of the exposed cylinder. The total circulation generated around a rotating cylinder in a real fluid is always less than that computed using the ideal flow considerations. A n experimental plot of circulation against rate of cylinder rotation is given by Rauscher [63]. Associated lift and drag coefficients are given by Goldstein [64]. A rotating cylinder can provide extremely high lift but the drag penalty could become severe as the rate of rotation increases. An airfoil with a leading-edge rotating cylinder could provide a good compromise with a higher lift and lower drag. 1.3 1.3.1 A Brief Review of the Relevant Literature Dynamics and fluid dynamics of bluff bodies A large number of objects encountered in the industrial aerodynamics practice are not streamlined, e.g. buildings, smokestacks, road vehicles, etc. Even a nominally streamlined airfoil at a high angle of attack presents a bluff geometry. Strouhal [1] was the first to correlate the frequency of vortex shedding to the free stream velocity of fluid and the diameter of a circular cylinder. The ratio, known as the 'Strouhal 17 M c • c According to Kelvin's Theorem, D(J u • dx)/Dt = 0 C < •'• Sc w Figure 1-6 u-dx = - f u.dx, C=C + C CB B W Schematic flow illustrating development of circulation around a rotating cylinder and the associated vortex shedding into the wake. The circulation around the circuit Cw embracing the wake but not the body must have a positive value Tc- Therefore, the circulation around a circuit CB embracing the body but not the wake has a negative value -Tc, as required by the fact that circulation around the closed curve C must be zero. 18 ENT EMENT i F i g u r e 1-7 Velocity profiles in the neighbourhood of separation: (a) wall moving downstream; (b) wall moving upstream [65]. 19 Figure 1-8 A bluff body with a rotating cylinder: (a) streamline patterns, (b) pressure distribution; i) , no rotation, ii) , with rotation. number', is independent of the Reynolds number in a certain range. The Strouhal number varies with shape and orientations of the object. There have been attempts to define a universal Strouhal number which is independent of the angle of attack for a given bluff body. Modi and Wiland [66,67] have used the projected diameter in case of a two dimensional ellipse and the results indicate that the variation of the Strouhal number with angle of attack is significantly reduced. Pioneering contribution in understanding of the bluff body flows was due to Karman [3]. He presented simple but accurate analysis of the stable vortex wake. It showed that, for stability, the alternately shedding vortices should arrange themselves in a ratio of transverse to longitudinal spacing, b/a, equal to 0.281. Kirchoff [68] was the first to use the concept of free streamlines in a hodograph model of the potential flow past a normal flat plate. But he assumed incorrectly that the base pressure in the wake is equal to the freestream static pressure and obtained a drag coefficient of 0.88 (the correct value is 1.98). Roshko [69] proposed a better|^d^r^h^model for a normal flat plate and obtained a value of 2.0 for the drag coefficient. The hodograph models are quite useful but difficult to construct for complicated bluff bodies. A much simpler and more versatile wake source potential flow model was presented by Parkinson [70]. It uses the concept of wake vortices located on a mapping circle to force separation of the incoming flow. The method is applicable to an airfoil, ellipse, circle, normal flat plate, etc. The wake base pressure downstream of the body is required as an input in the analysis. A number of reviews relating to the bluff body aerodynamics have appeared in past. The book by Blevins [71] on flow-induced vibrations provides useful information to the practicing engineers. Other books related to bluff body aerodynamics, flowinduced vibrations and other topics of interest to wind engineers include those by Simiu and Scanlan [72], Houghton and Carruthers [73], Sachs [74], and Walshe [75]. Some of the recent reviews on the subject are by Bearman [16], and Parkinson [19]. 21 Bearman covers mostly experimentally observed vortex shedding phenomenon and its effect on vibrations of a circular cylinder. A paper by Bearman and Davies [76] reports that, for a D-section and a triangular geometry with vertex pointing downstream, there is a decrease in base pressure at lock-in of vortex shedding frequency with the system natural frequency. The review by Parkinson [19] covers flow-induced vibrations of circular, D-section and rectangular cylinders together with a discussion on the^quasisteady theory modelling galloping instability. The department of mechanical engineering at the University of British Columbia has been quite active in studying the flow-induced vibrations of bluff bodies. Several studies have been carried out by Parkinson and Modi et al. over the years. The first series of experiments were carried out by Brooks [77] on the galloping characteristics of rectangular prisms with the aspect ratio (AR) ranging from 0.38 to 2.5. The wake survey and Strouhal number measurements were also conducted. Smith undertook additional experiments up to an aspect ratio of 6 and also studied the effect of varying damping on the vibrations [78]. These results indicate that bluff sections with short afterbodies (AR < 0.5), like a D-section, do not gallop from rest but are susceptible to vortex resonance type vibrations. D-section has been found to behave as a hard oscillator requiring a threshold amplitude to initiate galloping over a certain range of wind speeds. Slightly larger afterbodies ( A R between 0.75 to 1.5) like a square section gallop from rest as the wind speed is increased from zero (soft oscillator). Long rectangular prisms (AR > 3) are stable in the plunging mode at all wind speeds. This is probably due to the reattachment of the separated shear layers on the top and bottom sides of the rectangle. Heine [79] carried out measurements of fluctuating pressure on both stationary and oscillating models of rectangular sections with A R of 0.5 and 2. His results indicate that the self-excited oscillations of a bluff body do not affect the frequency of vortex shedding and pressure fluctuations. Measurements on a circular cylinder by Ferguson [80] suggest presence of amplitude modulations in 22 the fluctuating pressure signal which are in phase with the vortex-excited response. The results also indicate that the modulations are critically dependent on the wind speed. Santosham [81] conducted force measurements on rectangular and D-section cylinders and compared the dynamical results with predictions of the quasi-steady theory as proposed by Parkinson [19]. It was assumed that the bluff body at any instant during the galloping cycle experiences the same lateral force as that when held stationary at the same angle of attack. The experimental results compared well with the theory. Feng [82] carried out measurements of span-wise correlation of wake velocity for stationary and oscillating circular as well as D-section cylinders. The results indicate that the vortex wake is highly three dimensional. The response measurements showed hysterisis type behaviour, i.e. higher amplitudes were obtained when the wind velocity was increased over a certain range than when it was decreased back over the same range. The presence of hysterisis is not fully explained yet, but it is believed to be due to nonlinear interaction between the wake and the afterbody, and is not associated with the elastic behaviour of the system. Wind-induced vibrations of a structural angle section were studied in detail by Slater and Modi [83]. The frequent use of structural angle section in construction of towers, antenna masts and bridges, etc. warrants that their galloping and vortex resonance type response to wind forces be controlled. The wind tunnel tests showed that the angle section is susceptible to galloping in either plunging or torsional degree of freedom. The results correlated well with the predictions of the quasi-steady theory [19]. The quasi-steady theory was applied by Novak [84] to obtain a universal response curve which can be calculated from static force measurements on a rigid section or derived from response of an arbitrary model having the same type of cross-section. From these studies several important points, as listed by Parkinson [19], emerge: a) Conditions in the wake are primarily responsible for all forms of bluff body 23 vibrations and needs further investigation. b) The Karman vortices are the dominant source of excitation. During galloping, their effect is to cause an asymmetric mean pressure on the bluff body to produce transverse vibrations at small angles of attack. c) The most important shape parameter for a bluff body is the afterbody length. d) The large amplitude oscillations of a bluff body are essentially nonlinear in nature. 1.3.2 Vibration and flow control devices The suppression of wind induced vibrations has received considerable attention over the years. A n excellent review of various techniques developed to this end has ' been presented by Zdravkovich [85,86]. He has classified and compared a wide variety of aerodynamic and hydrodynamic means for suppression of vortex shedding and associated vibrations of bluff bodies. The three general categories described are: (i) surface protrusions (helical strakes, wires, fins); (ii) shrouds (perforated shrouds, gauze shrouds, axial slats); (iii) Near wake stabilizers (splitter plates, guiding vanes, base-bleed). In general these vibration control devices are passive in character. The surface protrusions like helical strakes, which are probably the most effective, affect the separation lines and the shear layers. condition in the wake. The shrouds influence the entrainment In contrast, the nearwake stabilizers prevent interactions between the entrainment layers. Generally, surface protrusions and shrouds are omnidirectional while the near-wake stabilizers are uni-directional. It is of interest to note that most of these devices need to be applied only to the top one-third portion of the structure (e.g. a chimney). Although some of the devices mentioned above are quite successful in suppressing the vortex shedding and the associated vibrations, their application may result in a significant increase in the drag. It has been observed that the effectiveness of strakes decreases with an increase in the intensity of turbulence as well as free stream velocity. 24 Another novel approach to control oscillations of structures has been studied by Modi et al. [13]. It uses the concept of energy dissipation due to sloshing of liquid in torus shaped nutation dampers. Since the civil engineering structures like buildings, bridges, smokestacks, airport control towers etc. have a natural frequency of about 1 Hz or less, the nutation dampers are especially suitable. Originally they were designed to control the very low frequency librational oscillations of orbiting satellites. A systematic parametric study involving overall mass of the damper as a percentage of the structure mass, damper geometry, liquid height and excitation amplitude, forced and natural frequencies, etc. was carried out by Welt [11-13]. Several designs of nutation dampers have proved to be quite successful in controlling the vibrations of structures and have been implemented in practice. As mentioned earlier, the concept of moving surface boundary layer control (MSBC) has been studied in detail by Modi et al. [53-60]. Several applications of the M S B C are shown in Figure 1-9. A rotating cylinder is used to inject momentum into the boundary-layer in order to keep it attached to the surface in the region of adverse pressure gradients. A systematic wind tunnel study of the M S B C as applied to a flat plate, an airfoil, rectangular prisms and tractor-trailer truck models has been carried out by Modi et al. [56,58-60]. A rotating cylinder located at the leading edge of a Joukowski airfoil increased the lift by about 200 percent and delayed the stall angle to 48° . Further lift enhancement could be achieved by locating one extra cylinder on the top surface of the airfoil [53-55,57]. Two counter-rotating cylinders located at the top and bottom leading edges of a normal flat plate could reduce the drag coefficient by 75% at Uc/U = 3, where Uc/U is the ratio of the cylinder surface speed to the free stream speed. Similar results were achieved for rectangular prisms of varying aspect ratios. The flow visualization pictures taken during the experiments dramatically revealed effectiveness of the M S B C in promoting the attached flow and narrowing the wake. In some case (e.g. a Joukowski airfoil at a high angle of attack of 2 5 ° ) , the 25 flow field became essentially potential, i.e.the boundary-layer separation was virtually eliminated. A 1:12 model of a tractor-trailer truck fitted with a rotating cylinder at the top leading edge of the trailer resulted in a spectacular drag reduction of 26 % [56,58-60]. The flow-visualization pictures clearly showed that the turbulent vortex wake adjacent to the bluff body is almost eliminated at a sufficiently high speed of rotation of the cylinders. These encouraging results, though preliminary in character, suggest a need for more serious closer look at the concept from both aerodynamic and dynamic considerations. Figure 1-10 shows possible future applications of the M S B C to a bridge-tower and a supertall building. 1.3.3 Numerical simulation of bluff body flows The past four decades have seen revolutionary progress in the field of com- puter technology. Today, the subject of computational methods in fluid mechanics and wind engineering has evolved into a major field of research. Historically, Computational Fluid Dynamics (CFD) received serious attention since introduction of the panel method simulating the potential flow around a physical object. At present, supercomputing capability has made it possible to simulate the complete Navier-Stokes equations in three dimensions, accounting for the turbulent flow field. Today, significant amount of design and development activities in the aerospace industry are based on computer simulations. Over the years, several journals and conferences dedicated to C F D have been established to meet the growing demands of the research community. Unprecedented growth in computer technology and sophisticated numerical methods promise a future where considerable amount of research and technological developments in fluid mechanics would be routinely performed using numerical simulation tools. This is especially true for the bluff body flows. Hence in any serious discussion on bluff body fluid dynamics, it is appropriate to provide at least some appreciation as to the state of the art in this important area. Of course, to review the vast body of literature on the C F D , accumulated over the years, is an impossible 26 •-. (»•) Leading-Edge The Basic Configuration Cylinder 0.05 C Figure 1-9 Applications of the Moving Surface Boundary-Layer Control ( M S B C ) procedure: (a) airfoil type geometries. 27 Two dimensional flat plate Figure l-9(cont.) Applications of the Moving Surface Boundary-Layer Control (MSBC) procedure: (b) bluff bodies such as a flat plate at a large angle of attack; rectangular prisms; and tractor-trailer truck configurations. Boundary-layer control and the associated drag reduction through fence type tripping device is also indicated. Figure 1-10 The MSBC concept presents an exciting possibility of applications to next generation of civil engineering structures such as bridge-towers and supertall buildings. task! Here an attempt is made to touch upon only those aspects of the C F D directly relevant to the problem of bluff body fluid dynamics and boundary-layer control. The first ever symposium dedicated to Computational Wind Engineering (CWE) was held in 1992, in Tokyo, Japan. The proceedings of the symposium represent a window to appreciate advances in the field [87]. The Navier-Stokes equations provide a complete mathematical description of the physics involved in a fluid dynamic phenomenon. A set of initial and boundary conditions complement the equations which are unique for every fluid dynamic process and thus impart distinct character to the problem under consideration. A solution would predict fluid flow behaviour such as boundary-layer formation and separation, vortex shedding, momentum transfer, in the wake, etc. The phenomenon of turbulence, suspected to be a chaotic behaviour of a deterministic dynamical system, is also governed by the Navier-Stokes equations in conjunction with an appropriate turbulence model. Low speed bluff body flows are governed by a mixed set of elliptic-parabolic form of the Navier-Stokes equations. A number of books and review articles [87-90] provide a general introduction to the subject of numerical simulation of the NavierStokes equations. There are several popular approaches among the researchers to numerically model the Navier-Stokes equations. One of the earliest and still the most widely used approach is the Finite Difference Method (FDM). The governing nonlinear partial differential equations are descretized using finite difference approximations. A n iterative scheme is generally necessary to solve the resulting, nonlinear set of algebraic equations. In the last two decades many algorithms have been proposed to incorporate some kind of finite difference scheme. Also for bodies with complicated shapes various mesh generating schemes have been developed. A good introduction to the subject has been given by Roache [88], Patankar [89], Anderson et al [90], and others. 30 A n alternative approach is to use some kind of finite element formulation procedure [91-94] popularly known as the Finite Element Method ( F E M ) . The F E M technique originated in 1950s to analyze complex structural problems of the aircraft industry. Due to nonlinear and dissipative (viscous effects) nature of the Navier- Stokes equations, use of the F E M in fluid mechanics was not preferred for a long time. But the situation is changing rapidly and most of the commercially available C F D softwares have been developed based on the finite element method. The F E M has turned out to be more versatile, since one can readily solve complex problems involving dynamic interactions between the fluid and a structure. This suggests its potential in simulation of bluff body fluid mechanics and fluid-structure interaction problems encountered in wind engineering. As compared to the F D M , the F E M does not require mesh generation programs to descretize a complicated flow domain involving arbitrarily shaped geometries. In addition to the F D M and the F E M , a totally different class of numerical approaches have been developed which are commonly referred to as the vortex methods. The vorticity transport equation is simulated by a collection of distinct vortices. These methods range from simple to highly sophisticated, requiring the use of a supercomputer. One of the methods (or its variation) known as the Discrete Vortex Method ( D V M ) has been routinely employed to study vortex shedding past a flat plate, rectangular and circular cylinders, and other bluff bodies. A numerical study of the stability of rows of vortices individually free to move at induced velocities was made by Abernathy and Kronauer [95]. The results clearly show the pattern of vortices observed in the wake of bluff bodies. Similar methods have been used, notably by Gerrard [96] and Sarpkaya [97], to model the unsteady flow around a circular cylinder. A n analogous model was developed by Laird [98] in which the attached boundary-layers were also reproduced and new vortices were added at the separation points. Similar to Laird, Chaplin [99] modelled the flow around a 31 cylinder free to execute transverse oscillations. The frequency of oscillations obtained numerically agreed with the experiments quite well, but the amplitude of the lift and drag forces exceeded that of the experiments, due to the absence of turbulent entrainment in the shear layers. Sarpkaya [100] has also modelled vortex shedding from an inclined flat plate using the discrete vortex method. The vortex shedding and Strouhal number re- sults agree favourably with the experimental data, but the calculated normal force coefficient is 20-25% larger than that measured by Fage and Johansen [101]. A similar study of vortex shedding past an inclined plate by Kiya and Arie [102] shows good agreement with the flow visualization results. Their results suggest that the vortex street behind the plate inclines as a whole towards the direction of the timeaveraged lift force exerted on the plate. It is also indicated that the vortex shedding from one edge of the plate will not occur at the mid-interval of the successive vortex shedding at the other edge. In another paper,.Kiya and Arie [103] have performed discrete-vortex calculation of the near wake velocity profiles, turbulence intensities and Reynolds stresses. They have also introduced a simple constant to achieve the vorticity transfer from large scales to small scales of turbulence. Several other bluff geometries have been studied using the D V M , such as the flow past a square section [104], a rectangular prism [105], a normal flat plate [106], two parallel plates [107], and two circular cylinders arranged side by side [108]. Chein and Chung [109] have performed discrete vortex calculations on an inclined and vertical plates using a different approach to determine the strength and location of the vortices just shed from the edges of the plate. In their model, the vortices are shed in such a way that their presence in the flow field offsets the singularities that exist at the plate edges as predicted by the potential flow solution. A review of various vortex methods is presented by Leonard [110] and Maull [111]. In a recent study by Leonard and Koumoutsakos [112], new advances in vortex methods for the simulation of un- 32 steady incompressible flow past a bluff body, the fast-vortex algorithm for convection, and the method of particle-exchange for viscous diffusion are discussed. Results of a transient flow past a circular cylinder for Re = 40-9500 have been reported [112] using a vector supercomputer. This method promises to become a powerful tool for accurate numerical simulation of bluff body flows in three dimensions by utilizing massively parallel computer architecture. Pan, Chew and Lee, [113] have proposed a new vortex calculation method using "generalized conservation of circulation model". By introducing a "generalized velocity field", the generalized circulation for a viscous flow is obtained. This simple algorithm is aimed at reducing the computational effort. In an interesting paper, Bienkiewicz and Kutz [114] discuss the simulation of unsteady separated flow past a bridge-deck section employing the discrete vortex method. Their numerical results show only a fair agreement with experimental data due to difficulties encountered in modelling flow reattachment and identification of the dominant separation nodes. A paper, by Fernando and Modi [115] describes in detail a mathematical model based on the D V M for predicting the performance of a Savonius wind turbine. Results show good agreement with flow visualization and wind tunnel data for both stationary and rotating Savonius configurations. In general, fluid flows in nature exhibit turbulent behaviour to some extent. Turbulence is strongly affected by "extra strains" including stratification, buoyancy, rotation, chemical reaction and compression [116]. The structure of the turbulence varies considerably depending on the flow regime. A unified model of turbulence versatile enough to describe the multitude of phenomena associated with turbulence is still not available. Bluff body flows encountered in wind engineering could be quite complex due to one or more structures immersed in a highly turbulent flow, without any preferred direction of the wind, and may require a full three dimensional treatment. 33 Although the steady force on an object such as a building can be computed reasonably well by a Reynolds-averaged Navier-Stokes (RANS) turbulence model, e.g. the k-e model (k=turbulent kinetic energy, €=dissipation rate), the intensity and frequency of the unsteady forces are difficult to predict by such a method. Since the geometry demands nothing less than a three dimensional treatment, the computational effort becomes too expensive [116]. Generally, simulation of turbulent flows . through the solution of the unsteady Navier-Stokes equations in three dimensions can be accomplished in two different ways. In the Direct Numerical Simulation (DNS) of a turbulent flow, the NavierStokes equations are solved for the largest scale of motion to the smallest scales in the dissipation range (Kolmogorov range). If the boundary and initial conditions are applied accurately and the numerical errors are kept to a minimum, the DNS is capable of reproducing a physical experiment in all of its details. The limitation of the DNS arises from its inability to deal with complex geometries and high Reynolds number flows due to the skyrocketing computational effort. In terms of accurate modelling of turbulence, potential of the DNS is unsurpassed, but it is unlikely to become a practical engineering tool [116]. In general, at low Reynolds numbers, the flow is a laminar, two dimensional periodic motion which can be studied fairly well with the present-day numerical methods. At higher Reynolds numbers, which usually occur in practice, stochastic three dimensional turbulent fluctuations may be superimposed on the two dimensional periodic vortex shedding motions [117]. Although, in principle, this can be accomplished by resorting to the DNS, presently it is feasible only at relatively low Reynolds numbers. The dissipative part of the turbulent motion has such small scales compared to the bluff body dimensions that it cannot be resolved in a numerical calculation. The number of grid points required to resolve this motion increases approximately as R e [117]- 34 3 A method that can be applied to situations at high Reynolds numbers is the Large-Eddy Simulation (LES) which consider only the larger-scale motions. However, the effect of the smaller-scale motions, which cannot be resolved on a given grid, needs to be accounted for. This effect is mainly dissipative. It can be achieved in two ways : (i) through a sub-grid scale model for determining turbulent stresses introduced by the smaller-scale fluctuations; (ii) the energy dissipation through damping introduced by the numerical scheme [117]. LES is a very powerful tool but it requires enormous amount of computing power and hence very expensive to implement. Attempts have been made to simulate turbulence in vortex shedding flows with statistical models which do not resolve any of the stochastic turbulent motion but average it out altogether [117]. Several statistical turbulence models are available to compute the Reynolds stresses < u^u'j > appearing in the ensemble averaged Navier-Stokes equations simulating the vortex, shedding flow past a bluff body. One approach is known as the eddy-viscosity model [117]. Alternatively, the Reynolds-Stress-Equation (RSE) solves the transport equations for the individual Reynolds stresses < u'^u'j > [118]. With the various turbulence models, different approaches were tested for handling the nearwall region. Several strategies are discussed by Deng [119], Nagano and Tagawa [120], Franke and Rodi [118], Jansson [121], Norris and Reynolds [122], and others. A good introduction to the LES method can be found in the review paper by Ferziger [123]. Until now, the numerical methods most frequently used for descretization of spatial derivatives in the DNS and the LES have been spectral in character. For geometrically simpler domains with simple boundary conditions, they offer maximum accuracy at a given cost. This has been reviewed at length in the book by Canuto [124]. It is expected that in the near future the vector supercomputer would reach its ultimate performance limit at the level of 100 G F L O P S (Giga Floating Point Opera- 35 tions per second) [125]. Any further increases in computational speed is likely to come from parallel processing technology. The C F D applications represent probably the most computational intensive tasks ever presented to a computer to perform. Hence, it is natural that research efforts are being directed towards accelerated implementation of the C F D to other disciplines (e.g structural analysis) on massively parallel computer systems. As noted by Murakami [126], in future, large scale computation of turbulent flows past bluff structures will of necessity rely on the parallel computation methodology. It would be appropriate to have at least some appreciation of the present and future developments in parallel processing, which is soon expected to become a norm. In U S A , the Numerical Aerodynamic Simulation (NAS) Systems Division at N A S A Ames Research Center has embarked upon an ambitious program of development and implementation of parallel processing (both hardware and software) technology [127]. The 'High Performance Computing and Communication Program ( H P C C P ) ' at N A S A focusses on the use of parallel processing for extremely computation intensive, multidisciplinary applications involving aerodynamic, propulsion, controls, structural and acoustical analysis necessary in aerospace missions of the future. These grand challenge problems may involve: 5,000,000 grid points; 50,000 iterations; 5,000 operations per point per second; and 10 15 operations per problem! For an automated design process requiring 0.1-0.01 hour turn-around time, it would involve 3-30 T F L O P S (Tera Floating Point Operations per second) capability with 4-40 gigawords of memory storage per problem [127]. A n example of the current state of the art in parallel processing as applied to the C F D can be found in the paper by Kato et al. [125]. Computations were carried out for an airflow around a two dimensional square rib placed within a boundarylayer by means of a massive parallel computer consisting of 1024 processors arranged in a M I M D (Multi-Instructions Multi-Data) network. Since shared memory paral- 36 lelism is inapplicable due to a conflict in memory access, a network of loosely-coupled multiprocessors with each having its own memory was preferred. The result of their investigation confirm the fact that the local memory parallelism requires a great deal of data exchange between the processors. The efficiency of parallel computa- tion is significantly affected by the amount of inter-processor communication. It is also necessary to take into account effects of granularity (ratio of calculation time to communication time for each processor). For a network of 512 processors the computational speed was boosted by a factor of 420 compared to that of a single processor, i.e. an efficiency of 82% (i.e 420/512 x 100) of parallel computation. The future of parallel processing critically depends on resolving several key issues regarding the implementation on parallel architecture machines, development of numerical algorithm for the parallel C F D , performance evaluation and the necessary computer science involved in achieving Teraflops capability. Such a technology would be extremely useful in carrying out the DNS or LES simulation of turbulent flows past complex bluff body configurations. As pointed out before, the precise numerical solution of a given problem can be obtained by solving the general time-dependent Navier-Stokes equations incorporating a suitable turbulence model. However, for realistic values of the Reynolds number, this demands enormous computational effort and cost. On the other hand, judicious modelling of the flow character can provide information of sufficient accuracy for all engineering design purposes with nominal computational tools and insignificant cost. To that end extension of the well developed panel code to multielement bluff bodies with moving surface boundary-layer control appears quite attractive. Numerical modelling of inviscid flows can be accomplished using the panel method. Over the past several decades, this approach has evolved to a sophisticated level where it can successfully tackle two and three dimensional geometries with flow separation, shock wave formation and unsteady aerodynamic phenomena. Even though the 37 panel method uses potential flow assumptions, various schemes have been devised to account for viscous dissipation, flow separation and existence of a finite wake. In the classical paper by Hess and Smith [128], the authors describe the original panel method, which replaces the body surface by source type singularities to obtain a numerical solution of the governing Fredholm integral equation of the second kind. The method is quite general in the sense that it can be easily extended to realistic flow regimes past ensemble of bodies, with nonrigid surfaces, moving with respect to each other; internal flows; suction for boundary-layer control; etc.. Nevertheless, certain classes of potential flow problems like partially unknown location of boundaries (e.g. trailing-vortex wake of unknown position) are excluded [128]. Recently, Farassat and Brown [129] have applied a particular form of the integral representation of the velocity potential to study helicopter rotor and propeller induced noise during flight. Using a similar approach, Lee and Yang [130] have developed a panel formulation to handle complex geometries (e.g. rotors, propellers etc.) undergoing complicated motions. Yon, Katz and Plotkin [131] have investigated the effect of trailing-edge thickness on the solution accuracy of a given panel method. The panel method under consideration uses a combination of source/doublet distribution with Dirichlet type boundary conditions. The results suggest that both the use of larger panels making a bigger angular separation at the trailing-edge or an alternate form of Kutta condition requiring the upper and lower trailing-edge velocities to be equal give much better results. The information obtained suggests the sensitivity of the solution to the size and distribution of panels. Kida et al. [132] discuss the accuracy of the panel method with distributed sources applied to two dimensional bluff bodies. Their investigation shows that the panel method with point sources is mathematically superior compared to that with linear sources. A panel method based on linear sources would diverge if the number of panels is even, which results in a singular coefficient matrix. Katz 38 and Plotkin [133] provide a good introduction to various panel methods and relevant computational details. Mokhtarian [134] has used a type of panel method with boundary-layer correction technique to predict the separation point on an airfoil with moving surface boundary-layer control. It accounts for the wall confinement and involves replacement of the airfoil and wind tunnel walls with vorticity distribution in conjunction with appropriate constraint relations. Inclusion of a source within the contour of the airfoil models the wake when there is flow separation from the surface. A finite difference boundary-layer scheme is used to introduce viscous corrections. The scheme employs potential flow pressure distribution results to calculate the boundary-layer characteristics at the top and bottom surfaces starting from the stagnation point until the point of separation. The procedure uses the displacement thickness to construct an equivalent airfoil and then iterates between the potential flow and boundary-layer scheme to converge to the final pressure distribution. Thus the objective is to match the outer potential flow solution with the inner boundary-layer prediction. The thin shear layer approximations of the Navier-Stokes equations for steady, two dimensional, incompressible flow are used. The finite difference method employed for viscous corrections is due to Keller and Cebeci [135,136]. The eddy viscosity term is expressed as suggested by Cebeci and Smith [137] which treats the turbulent boundary-layer as a cdmposite layer consisting of inner and outer regions with separate expressions for eddy viscosity in each region. The details of the formulation and the finite difference procedure followed are those given by Cebeci and Bradshaw [138]. Celik et al. [139] have reported a viscous-inviscid interaction method to calculate the mean flow past a smooth circular cylinder. The method is applicable to flows with a Reynolds number as high as 10 . 8 To handle separated flow problems as well as other complications arising due 39 to geometry, flow characteristics and boundary conditions, many extensions of the basic panel method of Hess and Smith have been proposed. In the method developed by Maskew and Dvorak [140], the interface between the potential flow region and the trailing-vortex wake have been modelled,as a simple shear layer. The enclosed wake region between the upper and lower free-vortex lines is interpreted as a potential flow region with a reduced total pressure. The model is essentially formulated as mixed boundary condition problem. On the active part of the solid surface the position of the boundary is known but the strength of the panel vorticity is unknown. On the 'free-vortex lines' the inverse is true: the vortex strength is a known constant but the position is an unknown. A n iteration scheme is required to establish the wake shape, i.e. the position of the 'free-vortex lines' by integrating the induced velocity vectors. The particular advantage of the 'free vortex lines' model is that the low pressure in the separated zone can be calculated directly. An insight to the problem of viscous separated flow involving dispersion of vorticity field in the wake was given by Pearcy et al. [141] based on experimental observation of various cylindrical bodies having different surface roughness. In an attempt to explain the movement of separation point from bodies having smooth contours, it was argued by the authors that a balance must be maintained between the rate at which the vorticity field is shed into the wake and its subsequent dispersion by diffusion, dissipation and cancellation among mutually opposite shear layers. In the improved model of Ribaut [142], the vorticity field is allowed to dissipate along the 'free vortex lines'. The vorticity dissipation model thus indirectly accounts for the presence of viscosity. Applying this model to a"flatplate and blunt trailing edge section, Ribaut demonstrated that extent of, and pressure in, the wake are essentially determined by the amount of viscous diffusion and dissipation of the vorticity. As a result a finite closed wake is obtained. A numerical approach, based on the first-order panel method employing lin- 40 early varying vorticity plus a constant source strength distribution, has been developed by Mukherjea and Bandyopadhyay [143,144] to model the separated flow past a wedge configuration. The solution obtained by this method, using a linear model of vorticity dispersion in the wake, shows close agreement with the experimental results over a range of the wedge angle. 1.4 Scope of the Present Investigation A rather brief review of the literature, more relevant to the present investi- gation, provides some understanding of the present state in the area of industrial aerodynamics. With that as background, a comprehensive program of study was formulated to explore the effect of Moving Surface Boundary-layer Control (MSBC) as applied to a family of slender as well as bluff bodies. The two-dimensional models used in the investigation include a symmetrical airfoil, a D-section and a set of rectangular prisms representing geometries of increasing bluffness. The focus is on the effect of the M S B C on the aerodynamics of the models when they are stationary, as well as their dynamical response during vortex resonance and galloping type of instabilities. The problem is approached using wind tunnel tests, numerical analysis and flow visualization thus providing better appreciation of the subject. To begin with, details of the models used, wind tunnel test methodology, modified panel method approach and flow visualization procedure are described in Chapter 2. This is followed by the results and their discussion for the airfoil, D-section and rectangular prisms in Chapter 3, 4 and 5, respectively. As all the results (wind tunnel test, flow visualization, etc.) for the bluff body under study are presented together, the individual chapter tends to be essentially self-contained. Finally, Chapter 6 looks at the dynamical response of the D-section and rectangular prisms in vortex resonance and galloping with the M S B C . Miscellaneous supporting information is contained in Appendices. The thesis ends with a summary of important results, original contribu41 tions, and suggestions for future investigation which are likely to be informative as well as satisfying. Figure 1-11 presents outline of the project. 42 Aerodynamics and Dynamics of Bluff Bodies in Presence of Moving Surface Boundary-layer Control ( M S B C ) Experimental Investigation T Flow Visualization Without and With MSBC Numerical Approach Airfoil, D-section, Rectangular and Square Prisms Dynamic Tests in the Wind Tunnel Static Measurements in the Wind Tunnel Steady-state Pressure Distribution Spectral Analysis of Unsteady Pressure Signals through FFT Aerodynamic forces Vibration Response during Vortex Resonance and Galloping Determination of the Vortex Shedding Frequency and Strouhal Number Figure 1-11 Airfoil, D-section, and Rectangular Prisms Vibration Suppression through MSBC Panel Method and Numerical Code One-Element Model Multi-Element Model Pressure Distribution, Aerodynamic Forces A schematic diagram showing the scope of the investigation. 2. E X P E R I M E N T A L I N V E S T I G A T I O N A N D N U M E R I C A L M O D E L L I N G 2.1 Preliminary Remarks As pointed out before, the investigation has three distinct phases: wind tun- nel tests, flow visualization and numerical simulation. A set of four two-dimensional bodies, representing increasing bluffness, was used during the study: Joukowski airfoil; D-section; rectangular prism; and square prism. Wind tunnel tests to assess aerodynamics and dynamics of the models are described first. This is followed by an explanation of the.flow visualization study using a closed circuit water channel. Finally, numerical analysis of the problem using the simple panel method is touched upon. 2.2 , Experimental Investigations 2.2,1 Bluff body shapes The set of bluff body shapes investigated is shown in Figure 2-1. Joukowski Airfoil , The Joukowski airfoil (Figure 2-1) is 16% thick with a chord length of 370 mm and is symmetric about the chord. The rotating cylinder diameter is 10% of the chord length and is located at the leading edge of the airfoil. It is driven by a 1/4 H.P., variable speed, A . C . motor. 42 pressure taps, 0.5 mm in diameter, are located at the top and bottom surfaces of the airfoil. The bluffness of the airfoil is a function of the angle of attack (a). At zero angle of attack the airfoil behaves as a.streamlined body. As the angle of attack is increased the bluffness increases. Variation of,the lift, drag and lift/drag ratio with the angle of attack.in presence of the M S B C are of primary interest in this study. The details of the motor connection to the rotating cylinder are shown in Figure 2-2. 44 The measurement of pressure on the surface of a rotating cylinder posed a problem. It was overcome in the following manner. Five pressure transmission tubes were located inside circumferential grooves as shown in Figure 2-3. The tubes remained stationary when the cylinder rotated. Thus the pressure taps remained immersed.in the boundary-layer and provided accurate estimate of the static pressure. D-Section The front face of the D-section is flat and normal to the free stream (at a . = 0) whereas the afterbody is semi-circular. Rotating cylinders are located at the two corners of the D-section. The presence of the cylinders causes rounding of the corners, i.e. departure from a sharp edged bluff body. This in turn, affects the vortex shedding process and drag of the body. The diameter of each rotating cylinder is 21% of the D-section width. Pressure taps are provided around the contour of the D-section, including' taps on the rotating portion. As before, the cylinders are driven by a 1/4 H.P., A . C . Motors, capable of very high rates of rotation (upto 22000 rpm). The aspect ratio (as defined in the list of symbols) of the D-section is 0.5. The D-section is likely to exhibit vortex resonance, galloping as well as autoration [145]: The present study is aimed at understanding of the D-section fluid dynamics and flow-induced instability in presence of the M S B C . Rectangular Prisms Two rectangular prisms selected for the study are also shown in Figure 2-1. The aspect ratios of the prisms are 1/2 (vertical rectangle) and 1 (square). The rotating cylinders are located at the two corners of the prism. The rotating cylinder diameter is 27% of the prism width. Pressure taps are provided around each model. The dynamics of a rectangular prism depends on its aspect ratio [19]. It is susceptible to both vortex resonance and galloping instabilities. The square prism, being a soft oscillator at all angles of attack, is particularly suitable for assessing the effectiveness of the M S B C in suppressing the vortex resonance and galloping vibrations. Note, 45 AR = 0.5 2R R R 2R AR = 1.0 AR = 0.5 (AR = Aspect Ratio = Width / Height) Figure 2-1 Bluff body geometries, without and with momentum injection, considered for investigation: (a) schematic diagrams. 46 Figure 2-l(cont.) Bluff body geometries, without and with momentum injection, considered for investigation: (bi) photograph of the Joukowski airfoil model. 47 Figure 2-l(cont.) Bluff body geometries, without and with momentum injection, considered for investigation: (b2) photograph showing details of the pressure taps and pressure conducting tubings. 48 Figure 2-1 (cent.) Bluff body geometries, without and with momentum injection, considered for investigation: (c) D-section and rectangular prism. 49 2-l(cont.) Bluff body geometries, without and with momentum injection, considered for investigation: (d) square prism and the dynamical test model. 50 A.C. MOTOR 1/4 H.P. , 3.8 AMP (VARIAC CONTROL) NO LOAD RPM : 22,000 FENNER COUPLING HIGH SPEED BEARING CLEARANCE, 1MM ROTATING LEADING EDGE CYLINDER WIND ROTATING CYLINDER (DIA., 38 MM) PRESSURE TAPS (TOTAL, 36) LOWER SUPPORT BRACKET HIGHSPEED BEARING Figure 2-2 Detailed schematic of the rotating cylinder and drive mechanism. 51 Figure 2-3 Schematic diagram showing details of the pressure taps near the rotating element. 52 the vertical rectangular prism represents a transition between the D-section and the square prism, thus providing an important insight into the bluff body fluid dynamics. 2.2.2 W i n d tunnel static experiments During static wind tunnel tests the bluff body model is held stationary. The wind tunnel used in the test-program is shown in Figure 2-4. It is a closed circuit tunnel with a test-section of 0.914m x 0.686m x 2.6m. The test-section is provided with 4 5 ° corner fillets varying from 0.152m x 0.152m to 0.121m x 0.121m" which compensate for the boundary-layer growth. The tunnel contraction ratio is 7:1. The wind speed can be varied form 1 - 4 5 m/s. A wire screen at the inlet to the test-section ensures "uniform velocity with a spatial variation of less than 0.25%. The intensity of turbulence in the test-section is less than 0.1%. The wind speed is calibrated using a Betz manometer which can read a variation of 0.2 mm of water. The wind tunnel is powered by a 230V, 57.2A, 15 H.P., D.C. motor driving a commercial axial-flow fan. A Ward-Leonard system is used for the speed control of the electric motor. A six component pyramidal force balance sits below the test-section and supports the model under test. The output of the force balance load cells are six electrical signals (one each for the lift, drag, and side forces; and the pitch, yaw and roll moments) which are transmitted to the computer based data acquisition system. Details of the instruments used during the experiments can be found in Appendix A. 2.2.3 Static test parameters (a) The mean pressure distribution around the body was measured using pressure taps provided on the circumference of the model. The pressure taps were connected to a scanivalve pressure transducer system. The scanivalve has 48 pressure ports and is driven by a solenoid controller. The pressure conveying tubes had little effect on attenuation and phase response of the signal [79]. The output signal is directed to the signal conditioner, amplifier and finally to the computer through an analog to digital ( A / D ) interface card. The experimental set-up is shown in Figure 2-5. The pressure signals were sampled at a frequency of 100 Hz 53 Figure 2-4 A schematic diagram of the low speed, low turbulence, closed circuit wind tunnel used in the test-program. Signal Conditioner and Amplifier A/D Interface and Data Bus Data Acquisition, Processing and Analysis I . I Scanivalve Pressure Transducer F i g u r e 2-5 Experimental set-up for static aerodynamic investigation. ing post processing, each signal is averaged to obtain the mean and rms (root mean square) components of the pressure. The mean pressure field around the surface of the body was integrated to obtain the aerodynamic force coefficients (the lift and drag). (b) The fluctuating component of the pressure is obtained by calculating the rms value from the time dependent signal. The rms value of the pressure signal gives estimation of the unsteadiness in the flow field in presence of a bluff body. (c) As the flow separates around the corners of a bluff body, a vortex is formed and is subsequently shed in the wake. This vortex creates suction pressure field which can be detected by pressure taps lying in its vicinity. By performing Fast Fourier Transform ( F F T ) analysis on the pressure against time signals, a spectrum of frequencies is obtained, with a dominant peak at the vortex shedding frequency. At times, due to weak strength of the vortex, none of the pressure taps could capture the vortex shedding frequency. In that disk probe was inserted in the wake close to the body to determine the vortex shedding frequency. The geometry and calibration charts for the disk probe as well as the test set-up are shown in Figures 2-6 and 2-7, respectively. The Strouhal number (St) was calculated from the vortex shedding frequency. (d) The experiments were performed in the angle of attack range 0-180°. Important changes in the flow character occur as the angle of attack is increased from zero'. It was anticipated that as the angle of attack approaches 45° - 60° range, reversal in the direction of momentum injection of the upstream cylinder would become necessary. Similarly, beyond a = 90° the afterbody faces upstream and the forebody turns towards the downstream flow direction. Again this necessitates changes in the direction of momentum injection. (e) The momentum injection parameter (U /U) was varied from 0 to 4 in steps of c one unit. 56 DISK P R O B E P O L Y E T H Y L E N E TUBES L=1.5M, D=1.68mm Q BAROCEL PRESSURE TRANSDUCER DAMPING BOTTLE SIGNAL o POWER SUPPLY CONDITIONER O FILTER OSCILLOSCOPE I Figure 2-7 1 Instrumentation layout for vortex shedding frequency and fluctuating pressure measurement. 58 (f) For the Joukowski airfoil and the D-section prism, the effect of roughness of the cylinder surface on efficiency of the momentum injection process was also investigated. Rotating cylinders with three different surface roughness were fabricated. The baseline case is the cylinder with smooth surface. The splined cylinder has 10 axially parallel surface grooves. The depth of these grooves is about 5% of the cylinder diameter. The third cylinder is characterized by roughness elements (5% thick, 18% wide and 55% long - based on cylinder diameter). The three types of rotating cylinders are shown in Figure 2-8. The longitudinal and circumferential grooves act as tiny vortex generators and it was anticipated that this would increase effectiveness of the momentum injection process. 2.2.4 Dynamical experiments in the wind tunnel The study of vortex resonance and galloping instability of a bluff body involved wind tunnel experiments simulating the dynamic response of a body. The wind tunnel utilized for static tests was also used for the dynamic experiments. The model was suspended on a dynamic test-rig with an arrangement allowing transverse oscillations (Figure 2-9). The rig is constructed of thick steel frame, on which two pairs of air bearings are fixed on the top and bottom support members as shown in the figure. Each air bearing was pressurized to around 40 N / c m . 2 The support frame is positioned in such a way that the the model is located inside the test-section of the wind tunneL The model can oscillate transverse to the freestream direction. The desired stiffness is provided by four springs attached between the model and the supporting structure. The controlled damping is provided by four eddy current dampers. The damping is proportional to the current level in the damper. The vibration amplitude of the model is measured by a displacement transducer. The time dependent signal also provides frequency of the response through spectrum analysis. The bluff body model design was similar to that in the static experiments, except 59 teal SMOOTH CYLINDER Figure 2-8 R O U G H CYLINDER SPLINED C Y L I N D E R Rotating cylinders, with three different types of surface conditions, used in the test program: (a) schematic diagram. c 60 COMPUTER DATA ACQUISITION SYSTEM Figure 2-9 A schematic diagram of the dynamic experiment set-up involving the wind tunnel, test-rig and data acquisition system. 63 for the pressure taps. The model was constructed from balsa wood to conserve weight (Figure 2-10). The photograph of the balsa wood model is given in Figure 2-1. The electric motors were mounted on top of the model and are coupled directly to the rotating cylinders. The rotating cylinders were constructed from hollow aluminium to keep the system light. 2.2.5 Dynamic test parameters (a) As mentioned above, the amplitude of vibration was measured by a displacement transducer. To identify regions of vortex resonance and galloping it was necessary to plot the vibration amplitude as a function of wind speed. Effect of momentum injection on the dynamical characteristics of the system was assessed by analyzing the amplitude versus wind speed plots at various speeds of the rotating cylinders. (b) Frequency of the vibration was obtained directly through the spectrum analyzer and also through F F T analysis of the amplitude versus time signal. (c) A disk probe, inserted in the wake downstream of the body, was. used to measure the vortex shedding frequency. This is further described in Section 2.3.3. (d) The experiments were performed at zero angle of attack. (e) The rate of momentum injection was varied from 0-4 in steps of one unit. (f) Rotating cylinders were located at. the leading (upstream) edges of each bluff body model (Figure 2-ld). Instead of rotating both the cylinders, it was decided to activate only one of them in order to investigate the effect of asymmetric momentum injection. 2.2.6 Flow visualization Flow visualization represents a powerful tool in the study of fluid dynamical problems as it provides physical appreciation of the complex character of the flow. Flow visualization pictures and video movie of the fluid dynamics of bluff bodies in presence of the M S B C were taken. The tests were conducted by the author using the facility at the Kanto Gakuin University, Japan. 64 ONLINE TV-VCR UNIT VIDEO CAMERA SOURCE OF L I G H T SLIT MIRROR 01 Cn Figure 2-11 A schematic diagram of the closed circuit water channel facility used in the flow visualization study. Slit lighting was used to minimize distortion due to three dimensional character of the flow. Long exposure provided path-lines with ployvinyl chloride particles serving as tracers. The dimensions are in mm. The flow visualization setup used is shown in Figure. 2-11. The experiments were carried out in a closed circuit water channel. Bluff body plexiglass model fitted with rotating cylinders was placed in the water channel. The water was seeded with fine polyvinyl chloride (PVC) powder. The P V C particles are almost neutrally buoyant in the flow and hence do not rise or fall as they move in the water channel. As shown in the figure, a slit lighting arrangement was used to illuminate a very thin plane of the test-section. The steady flow of P V C particles resulted in streaklines which were viewed from the front of the test-section. The flow velocity in most cases was fixed at 0.3 m/s. The experiments were carried out in a dark room resulting in a sharp contrast between the bright streaklines and a dark background. A series of long exposure photographs dramatically captured the streaklines, showing formation and subsequent dissipation of alternate vortex shedding process. A video was also taken to assess time evolution of the process. A sensitive C C D (charged coupled device) video camera was utilized for the purpose. 2.3 2.3.1 Numerical Analysis Formulation of the problem The numerical approach is based on the well-established vortex panel method. The surface of a given body is descretized into a large number of panels with vortex and source type singularities distributed on each panel. The overall flow at any point in the domain is a result of elemental contributions from each panel. Simulation of the rotating cylinder integral to a bluff body poses a problem. As a first approximation the rotating cylinder and the rest of the body can be considered a single object and is descretized into panels accordingly. This formulation is referred to as the "single element configuration". In reality, the rotating cylinder is separate from the rest of the bluff body, with the presence of a gap. In the multi-element panel method, the rotating cylinder and the rest of the body are treated as two distinct physical objects. 66 The details of both the formulations are described below. The general two dimensional Navier-Stokes equation is given by [146] ^ + (V • S7)V = y V P + "V V, (2.1) 2 where: V = fluid velocity at any given point; p = pressure; v = kinematic viscosity of the fluid; p = density of the fluid. The continuity equation for an incompressible fluid is y ^ = 0. (2.2) All body forces (e.g. gravity) have been assumed conservative and their potentials have been absorbed in the pressure term. For the potential flow situation the viscous terms vanish and the above set of equations reduces to Euler's equation dV ¥ -1 + (F-V)^=yVP. (2-3) The general fluid dynamic problem under consideration is characterized by a known boundary surface and the normal components of fluid velocity are prescribed on these boundaries. The boundary condition is V-n = 0, (2.4) where n is the outward normal from the surface 'S'. The set of equations (2.2) and (2.3) do not define the potential flow. To obtain potential form of the equation —• —• the condition of irrotationality is invoked. Let V = V q o + v, where: —• = free- stream velocity, satisfies continuity equation; v = perturbation velocity field due to the presence of the boundary, satisfies irrotationality condition, v — — \j4>\ and <f> = 67 disturbance velocity potential. Recognizing that y • v = 0 gives V V = 0- (2.5) This is a linear, elliptic second order partial differential equation. Now V — — y<^>, and the boundary condition takes the form (Voo - V0) • "I* = 0, d(b ^ with | \7 = V r o 6 -n|„ (2.6) —> 0 at infinity. The set of equations (2.5), (2.6) represents a well-posed elliptic boundary value problem with Neumann-type velocity boundary condition. Once the velocity field is obtained, the pressure is calculated by intergrating'equation (2.3) to give Bernoulli's equation where p and are the free-stream pressure and velocity, respectively. Reduction of the problem into an integral equation with singularity distribution on the body surface is carried out as follows. Consider an arbitrary body as sketched in Figure 2-12. It shows a combination of unit source and a unit vortex singularity at a point q(x ,z ). The potential at any point P(x,z) is given.by q q <b = In r(p q) + 0(p q) ) where: r(p,q) = x ) + (z - z ) ; 2 q t (2.8) 1 6(p,q) = tan~ [(z - z )/(x - x )]. 2 1 q q q In case of a continuous distribution of the singularities, the total potential at any point P(x,z) is a linear combination of such contributions. For a constant source strength (a) and linearly varying vorticity strength ( 7 ) distributed on the contour of 68 69 the body surface, the perturbation potential can be written as <t>=<f [(T(q)-lnr(p,q)+j(q)-e(p,q)^dl. (2.9) c Applying the boundary conditions and rearranging the terms gives [128], 2TT [a(p) + (p)] -<f>j7 \<r{q) • In r(p, q) + (?) • 0(p, q)] dl = -n(p) 7 • V^. (2.10) («) (0 This is a Fredholm integral equation of the second kind where:(i) represents contribution of the singularity obtained through the limiting process' [147]; (ii) is the kernel of the integral equation. The solution of the integral equation is obtained by approximating it as a set of linear algebraic equations. The body contour is descretized into a large number of line elements (panels), whose characteristic dimensions are small compared to those of the body. Over each panel a distribution of unknown singularities is prescribed. Thus the problem reduces to that of obtaining a finite number of singularity strengths, one for each panel. The form of the singularity distribution over each panel is a known geometric function and hence integration can be performed. Application of normal velocity boundary condition on each panel results in a set of linear algebraic equations relating unknown singularity strengths. The coefficient matrix consists of normal velocities induced by the panels at each collocation point for unit value of the singularity strength. Once the linear equations are solved, the flow velocities and pressures are evaluated readily. 2.3.2 Single element configuration In the single element panel formulation, the rotating cylinder and the rest of the bluff body are treated as a single physical object. The numerical formulation is described here using the airfoil as a bluff body model (Figure 2-13a). The body surface is descretized into a large number (N) of panel elements. Each element is a 70 combination of linearly varying vortex strength ( 7 ) and a constant source strength (a). The end points of any panel are known as nodes, whereas the mid-point of a panel is called the collocation point. The collocation point on a given panel serves as a convenient reference to evaluate velocity and pressure after implementing the boundary conditions. On an i-th panel, the nodal vortex strength varies from ( 7 ; ) to ( 7 i 4 _ i ) in a linear fashion. The source strength (cr) remains constant over all the panels. The source distribution on each panel prevents leakage of the fluid across the panel. For a total of N panels, there are (N+l) unknown vortex strengths and an unknown source strength, making a total of (N+2) unknowns. • Modelling of the 'Free Vortex Sheets' As shown in Figure 2-13(a), the separated flow is represented by the upper and lower 'free vortex sheets' issuing from the separation points. They are also descretized in to a large number (M) of panels. However, the vorticity strength in the wake does not constitute an additional unknown. At the beginning of the 'free vortex sheet', the vorticity strength equals that at the point of separation. If the shed vortex strength is taken to remain constant as it is convected along the 'free vortex sheet', an infinitely long wake would result with a constant wake pressure. Realistically, the wake is finite and closed, with a gradual rise in pressure. This can be attributed to the presence of viscosity in the flow. To account for this effect, the shed vortex is allowed to dissipate as it is convected downstream along the 'free vortex sheet'. The rate of dissipation is an empirical input to the program. Originally this idea was proposed by Ribaut [142]. The Kutta condition requires zero net sheddingof vorticity in the wake. Thus as we move along the panels on the 'free vortex sheets', the upper and lower shed vorticities remain equal and opposite in strength at all time. If a linear variation of vorticity dissipation is assumed, then the vorticity strength on the 'free vortex sheets' can be obtaining in terms of the vortex strength at the upper and lower separation 71 Y-, i 1 + 1 PANEL i y VARIES LINEARLY f j CONSTANT (b) F i g u r e 2-13 ' Numerical panel formulation'describing, (a) single element configuration and, (b) multi-element configuration. from the experimental results. 72 Note, 'S' is obtained points. Mathematically it can be expressed as 7 c = T i ( l - A - ) + W l - A - ) = 0, c c . (2.11) where A is the vorticity dissipation ratio, and Su and SL are the distances along the upper and lower 'free vortex sheets' from the separation points. The perturbation velocity components can be calculated by taking the effects of surface singularities together with contribution from panels on the 'free vortex sheets' as: JV+l Uj N =E ^ + M i=l JV+l C i=l JV 'i = Z) i=l M ' w E j'+E B M E^ +E ^ - +E ^ ^ ; \ CT +IR I) (- ) 2 12 i=l M a<i + J2 ^<i- F F (- ) 2 13 For N number of panels on the airfoil contour, the (N+2) unknowns are obtained by satisfying the boundary condition of zero normal velocity in eq.(2.14) at N collocation points together with two auxiliary conditions. The condition of zero normal velocity at the j-th collocation point on the airfoil surface can be written as dz *- idx W U dz Ucosa— dx - Usina (2.14) 3 Satisfying the boundary condition at the mid-points of N panels gives N linear algebraic equations involving (N+2) unknown singularity strengths. The other two equations are obtained as follows. The Kutta condition of zero net shedding of vorticity in the wake is satisfied by setting the vorticity values at the upper and lower surfaces the same (and equal to the value of the 'free vortex sheet'). Mathematically, this can be expressed as: 7 S = -7i J (2-15) 7JV+I = 0; (2.16) lL=7s; (2.17) 73 JLI = 7i; (2.18) where 7 5 is the vorticity strength at the upper separation point and 7 ^ u and 7 ^ are the vorticity strengths on the upper and lower wake panels, respectively. Once the solution is obtained, the perturbation velocity components may be reobtained from equations (2.12) and (2.13). In the case of the D-section the auxiliary conditions are Ti = 7JV+1 and j' wu = -y' . wl The total velocity can be obtained from qj = yj (UJ + Ucosa) + (WJ + Usina) , 2 (2.19) 2 and the pressure coefficient is given by C. v = 1 - ( /U) . (2.20) 2 qj Geometry of the Wake Since both the shape and the strength of the 'free vortex sheet' are unknown, the solution of the problem formulated above is obtained by an iterative procedure. The wake shape is initially assumed to be a straight line inclined at an angle the mean angle at which the separation vortices are convected. a/2, Once the solution is obtained based on this assumed wake shape, velocities are recalculated at the mid-points of the panels on the wake. The new wake shape is then determined by piece-wise integration starting at the upper and lower separation points. At each iteration, the wake influence coefficients only need to be calculated. The computation is terminated when a convergence accuracy of 1 0 - 6 is achieved. This is usually accomplished within 50-60 iterations. Wake Pressure The pressure in the wake was calculated using the procedure outlined by Mukherjea and Bandyopadhyay [143]. The base pressure coefficients in the separated 74 region can be calculated directly from the Bernoulli equation as _ r (P-Poo) 2' 1_ where qinner 9 t n n e r _ A (o 21) H is the velocity inside the separated zone and total pressure from that at infinity. Since qinner is small, C, = l - i ^ . AH C p is the decrease in the can be written as (2-22) The jump in the total pressure across the free shear layer can be calculated by assuming zero static pressure drop across it. If the average velocity along the free shear layer is denoted by Qav = -(? outer Hmner ), (2.23) then, for upper shear layer: Since the vorticity y' , wu qouter = qav + " ^ - J (2.24) qinner — qav (2.25) on the upper shear layer is Ifwu ~ Qouter AH T • Qinneri (2.26) can now be calculated as AH — Houter ~~ H^ inner j Pouter + ^P(qav + -^I'vjuf 1 * Pinner + ^Pylav = PlavlLu- 1 i o ~ ~^lwu) ) (2-27) Using eqs. (2.25) and (2.27), the coefficient of pressure in the separated region can 75 be calculated as inner = i- i 2^tuu u KJu) , 2 =• 1 - (<ys/U) . 2 ( (2.28) Modelling of the Vorticity Dissipation The wake in the present model is finite and its extension is primarily dependent upon viscous dissipation A. From eq. (2.27), AH = pqavl' u W = Planter ^)(1 ~ - A^) 7 5 . (2.29) At the closure of the wake, its width becomes zero, i.e. s u = S u - c/X. (2.30) With decreased dissipation, i.e. with reduction in A, the wake length increases. The wake becomes infinite in the case of ideal fluid when A = 0. Numerical Simulation of the Rotating Cylinder The effect of the rotating cylinder can be mathematically simulated either by a line vortex (Figure 2-14) at the center of the cylinder or by line vortices distributed over the discretized cylinder surface. The strength of the line vortex can be related, by a simple analogy, to the surface velocity and diameter of the cylinder. Mathematically the total circulation around an arbitrary body is defined as = judl, where: \u\ — \Uc\ =constant; \dl\ — (2.31) \r-d&\. r 76 Thus: u -r§;' de, z c A bluff body with a rotating cylinder Numerical Idealization U = cylinder surface speed C f = cylinder radius P = equivalent point vortex 0 = exposed portion of the cylinder surface £1 = angular velocity of the cylinder F i g u r e 2-14 Numerical modelling of a rotating cylinder using a line vortex. 77 = *Z-r-U Ci = %d-U . (2.32) c wherex d=2r. This circulation can be mathematically simulated by a line vortex of strength T. The effect of the moving wall can now be incorporated in a simple manner through the boundary condition by modifying the right hand side of eq. (2.11) for upper surface collocation points as dz (u + Ucosa)^dx c - (w + Usina), c (2.33) where u , w are the velocity components due to the line vortex of strength T situated c c at the centre of the cylinder. In this way it is possible to directly relate the circulation strength T to the speed ratio (U /U) through eq. (2.32). c 2.3.3 Multi-element configuration ' The next logical step is to develop a model for the moving wall problem based on the multi-element configuration by incorporating the effects of separation from both the cylinder and rear airfoil surface. This model, as before, involves distribution of linearly varying vorticity and sources on (N1-I-N2) number of panels descretizing the multi-element configuration (Figure 2-13b). Flow is assumed to separate from the cylinder surface at top (T) and bottom (B), and pass over the airfoil surface. Flow within the gap is assumed to be local with no effect on the external flow field. On the airfoil, the flow separates at point S on the upper surface and at the trailing edge (point 1) on the lower surface. The vortex sheets separating from the cylinder are assumed to undergo dissipation so that the vorticity strengths (jwu^wi) become zero at points S and 1. This results in different dissipation ratios for upper and lower separation panels given by: K = c/S ; (2.34) Xi = c/Si; (2.35) u 78 where S , 5; are the upper and lower wake lengths from the cylinder top and bottom u separation points to points S and 1 on the airfoil surface, respectively. For.a flow separating from the airfoil surface, 'free vortex'sheets' emanate from separation points S and 1. The vorticity strengths of the upper and lower separation panels can be obtained, as before, in terms of js and 71 through eqs. (2.15)-(2.18). However, since the shape of the separating panels are not known 'a priori', an iterative procedure, similar to the one described in Section 2.3.2, is adopted. The perturbation velocity components, in this case, may be expressed as: JVi+JVa+3 Nj_ i=l 2 i=l wl + ^ E i j i=l w l + Y^Eji7 (2-36) i=l N N i=l i=l ± 2 i=l i=l M wu M + Y^ ^u Y^ ^- (- ) F F wl i=l + Y^ ii<i-> N 2 + ^Fja u E wu i=l Ni+N +2 Hyji w M , 7 N i=l M N 1=1 N + 1=1 2 37 1=1 Using eqs. (2.36) and (2.37), boundary condition of zero normal flow can be satisfied at (N1+N2) collocation points on the multi-element configuration surface giving (N1+N2) equations. To obtain (N1+N2+4) unknowns, auxiliary relations, pertaining to flow separation and Kutta condition, are used. For the cylinder surface,, these conditions require the vorticity strengths just upstream and downstream of the top and bottom separation points be zero. Furthermore, 7(^+2) = 7(JV-Nv"2+2) 1 = ^' ^ o r the airfoil surface, these conditions require that vorticity strengths at the separation points be equal and opposite to the strength of the separating vortex sheets. Also, vorticity strengths just downstream of top separation point and upstream of lower separation point are taken to be zero. Numerical solution can be obtained, as before, by an iterative scheme starting 79 ' ; with an assumed wake shape. Once the convergence is achieved, velocity and pressure can be calculated using eqs. (2.19), (2.20), (2.36), and (2.37). 2.3.4 Computer implementation Computer codes for single and multi-element configurations were developed using F O R T R A N 77. The programs can allow for zero and non-zero dissipation rate, desired cylinder rotation speed, angle of attack, etc. The total number of panels varies between 100-200 depending upon the configuration used. The infinite wake in case of zero vorticity dissipation is fixed at 150 times the airfoil chord. The programs require about 50-60 iterations to achieve a convergence accuracy of about 1 0 , and - 6 take about 15 minutes on a SUN SPARC 2 work station. A flowchart outlining the panel method is described in Figure 2-15. Details of the mathematical expressions for calculating the influence coefficients are given in the Appendix B. 80 ^START ) [ R E A D INPUT PARAMETERS > t DESCRETIZE THE BODY SURFACE INTO LARGE NUMBER OF PANELS. REPEAT THE PROCEDURE FOR THE FREE-VORTEX LINES MODELLING THE WAKE > A S S U M E INITIAL WAKE SHAPE CALCULATE GEOMETRICAL PARAMETERS FOR PANELS ON THE BODY AND WAKE APPLY BOUNDARY CONDITIONS, OBTAIN THE R.H.S. OF THE EQUATION — , V SET-UP THE INFLUENCE COEFFICIENT MATRIX i INVERT THE MATRIX TO OBTAIN THE SINGULARITY STRENGTHS * CALCULATE P R E S S U R E S VELOCITY ON THE BODY AT EACH COLLOCATION POINT * CALCULATE OVER ALL FORCE AND MOMENT COEFFICIENTS : : *k OBTAIN THE NEW WAKE SHAPE BY INTEGRATION OF THE VELOCITY ON UPPER & LOWER FREE-VORTEX LINES P R E P A R E OUTPUT DATA, PRINT & PLOT ( Figure 2-15 STOP ) Numerical algorithm for the panel method. 81 3. A E R O D Y N A M I C S O F A N A I R F O I L W I T H M O M E N T U M 3.1 INJECTION Preliminary Remarks The Moving Surface Boundary-layer Control (MSBC) was applied to a two- dimensional symmetric airfoil. The leading-edge of the Joukowski airfoil was replaced by a rotating cylinder, which injects momentum into the boundary-layer on the upper surface of the airfoil. Cylinders with smooth, splined (axial grooves) or rough (axial and circumferential grooves) surface were tested to assess their relative performance. To begin with wind tunnel results are presented. Effect of various flow parameters on the aerodynamic coefficients is discussed to emphasize physical changes in the flow brought about by application of the M S B C . Next, the numerical simulation results are presented and compared with the experimental data. Results of single and two element panel configurations with attached or separated flow model are presented and their relative accuracy is compared. The wake shape iterations and the effect of number of panels on the convergence accuracy of the numerical scheme are also discussed. Finally, the flow visualization photographs are presented which clearly demonstrate effectiveness of the M S B C in delaying the flow separation. 3.2 3.2.1 W i n d Tunnel Investigation Lift and drag characteristics Variations in the lift coefficient with the angle of attack for smooth, rough and splined rotating cylinders are presented in Figure 3-1. The airfoil with a smooth leading edge cylinder is the benchmark case (Figure 3-la) and all other results are compared with it. In absence of cylinder rotation (Uc/U = 0), the airfoil stalls at around a = 10° regardless of the type of rotating cylinder employed. In case of the smooth cylinder (Figure 3-la), as the cylinder rotation speed is increased from Uc/U = 0 to 4, the stall is delayed progressively to higher angles of attack: 15° at 82 Uc/U = 1; 25° at U /U c = 2; and 3 0 ° - 3 5 ° at U /U c increase in the lift coefficient (CL) with a. range Uc/U = 3-4. There is a corresponding For a momentum injection rate in the = 3-4, the CL,max of 1.84 for the smooth cylinder represents an increase of around 160%! The corresponding gain with the splined cylinder surface was equally impressive at 210% (CL,max = 2.2). As a is increased beyond 35°, CL decreases, but the stall is very gradual (i.e. it does not fit the nominal definition of a stall, since a rapid loss in lift is avoided). Even for a as high as 5 0 ° , CL is about 1.48 (smooth cylinder), which is quite significant. The rate of change in CL with a can be estimated dC^/da. There is no significant change in the slope of the lift curve from the slope with the smooth rotating cylinder case. However, with the rough and splined rotating cylinders, there is a small rise in dC^/da with an increase in the momentum injection. Note, there is no appreciable further benefit in terms of an increase in lift when the momentum injection is raised from Uc/U = 3 to 4. This indicates a practical limit to the beneficial influence of the momentum injection. Compared to the smooth surface rotating cylinder, the rough cylinder has inferior performance (Figure 3-lb). The maximum CL is about 1.52 at a = 30° which is lower than that obtained with a smooth cylinder (CL,max = 1.84). The axial grooves on the splined cylinder scoop the incoming fluid, and thereby improve momentum injection in the boundary-layer. In the case of the rough surface cylinder, the axial splines are supplemented by circumferential grooves (Figure 2-8), which seem to disrupt the effective scooping mechanism. The axial splines extend along the axis of the rotating cylinder from end to end, maintaining a two-dimensional flow. On the otherhand, the circumferential grooves introduce three-dimensional disturbances upstream of the two-dimensional airfoil, resulting in a deteriorated performance. Table 3.1 lists Cz, ,a. >rn x and stall values for the three different cylinder surface conditions. In Table 3.2, average values of the slope dC^/da are listed. It is apparent that the performance of the spline surface cylinder is relatively superior. 83 Figure 3-1 Experimentally obtained lift coefficient for the Joukowski airfoil in presence of momentum injection: (a) smooth leading-edge rotating cylinder; (b) rough cylinder; (c) splined cylinder. 84 Table 3.1 Comparative performance of the three rotating cylinders in terms of the maximum lift coefficient and the stall angle for( U /U c Uc/U = 0 Cylinder Surface Uc/U = 3,4 dsta.ll smooth rough splined Table 3.2 Olstall 10° 15° 10° Average values of Range of a, deg. 0 - 20° 0 - 30° 30 - 50° = 3 and 4. 0.70 0.80 0.75 35° 30° 35° 1.75 1.50 2.15 dC^/da for three rotating cylinders (da in radian). Smooth Cylinder 4.30 3.44 -0.86 Rough Cylinder 3.87 2.87 -1.15 Splined Cylinder 5.01 4.11 -1.29 The corresponding drag characteristics for the airfoil are shown in Figure 3-2. There is an almost linear rise in the drag coefficient (CD) with a for'all the three cases. For the Joukowski airfoil with splined cylinder (Figure 3-2c), CD remains quite small for a < 15°. It attains a maximum value of around 2.5 at a = 50° for the airfoil with smooth or rough cylinder. On the other hand, maximum CD is about 1.95 in the case of the splined cylinder. The lower drag coupled with a higher lift in case of the airfoil with a splined cylinder appears quite promising for practical aerodynamic applications. In Table 3.3, the slope dCo/da is listed for Uc/U — 0 and 4. Note, by application of the momentum injection, flow separation near the leading edge is avoided and the local pressure on the upper surface remains much lower compared to that at the lower surface of the airfoil. This results in a large normal force acting perpendicular to the airfoil. Since the drag of an airfoil is proportional to the sine component (i.e. CD oc sina) of this normal force, it contributes more to CD as a is increased. The rate of momentum injection increases the slope dCo/da in all three cases (Table 3.3). 85 Table 3.3 Average values of dCn/da as affected by the cylinder surface condition and rate of momentum injection. Type of Cylinder Smooth Rough Splined Uc/U The lift to drag ratio (CL/CQ) = 0 (0 < a < 10°) 0.997 0.229 0.516 Uc/U = 4 (0 < a < 5 0 ° ) 2.865 2.865 2.292 for the airfoil with three different cylinder geometries is plotted in Figure 3-3. The CL/CD has a peak in the range 0 < a < 10°, before stall. In-the case of splined cylinder, large values of CL/CD are achieved, resulting from a combination of high lift and low drag coefficients. Beyond stall, the CL/CD ratio drops off rapidly and there is only little effect of momentum injection. At high angles of attack (a > 2 0 ° ) , both the lift and the drag coefficients show an increase keeping the CL/CD ratio low. Further appreciation in this regard can be obtained by examining CL VS. Cp plots for the three cases (Figure 3-4). For the splined cylinder, the steep gradient of the plot suggests a favourable performance. In Figure 3-5, variation of both CL as well as CL/CD- with a are plotted for the three cases. These results clearly establish significant changes brought about by the surface condition of a momentum injecting rotating cylinder. Results in Figure 3-6 focus on the variation of CL with a beyond the nominal stall angle of 10° (reference case), as affected by the surface roughness of the rotating element and Uc/U. For the smooth cylinder expected, the lift coefficient first dips at stall followed by a monotonic rise with an increase in the angle of attack in the range tested (a < 5 0 ° ) , in absence of the momentum injection. The effect of Uc/U « is to delay the stall to around 35° for Uc/U 160% as mentioned before. = 4 with an increase in CL by Now the discontinuity in the CL VS. a plot at the stall is eliminated. Note, the flat peak suggests onset of a gradual stall, a desirable feature. The presence of the splined surface accentuates this behaviour resulting in 86 Figure 3-2 Experimentally obtained drag coefficient for the Joukowski airfoil in presence of momentum injection: (a) smooth leading-edge rotating cylinder; (b) rough cylinder; (c) splined cylinder. 87 / 15.0 C /C L -i 60.0 L r D J C /C 1—i T 1 I i I i T L 1 •^•-"••O-.--.-030 40 i — i 1 r —rOr-r-6 50 T Re=143000 — • — u /u = o c - _ A - - U /U = 2 - _ _ . _ U /U = 4 D C 0 C 30.0 - (c) A_ 20 Figure 3-3 a o 30 40 50 Experimentally obtained lift/drag ratio of the Joukowski airfoil in presence of momentum injection: (a) smooth cylinder; (b) rough cylinder; (c) splined cylinder. 88 ~i 2.0 L 1 r (a) — O - A o-._. •o -o I A 1.0 _|_ O.Q. -I I l_ 2.0 1.0 0.0 2.0 o- -r~ • (b) i 3.0 -i I 1 r o~._ "O- _L 2.0 ~i „ ,.-o „o^ -' -A ..-O' ^A- ' o A 1 i r ~—I _i i 3.0 L. -I 2.0 — U /U = 0 A -- U /U = 2 O _ U /U = 4 C •O 1.0 i_ Re = 143000 — I i 1 I ^ —o-.^ A" ' _l ——o C - C I L. 3.0 D Figure 3-4 Variation of CL with C p for the Joukowski airfoil in presence of momentum injection: (a) smooth cylinder; (b) rough cylinder; (c) splined cylinder. 89 T — i — i — r i—r i . O ,o-_ J O ' ' _ -o- --•- -A- -A A •• ~ -A _ ~ 0 - . - -o- - -A - . (a): -I I I I I I ~ i — i — i — I — i — i — i — r l_ cc° 30 60.0 I -I I l_ 40 50 i — r SYMBOL CYLINDER SQUARE SMOOTH TRIANGLE ROUGH CIRCLE SPLINED Re = 143000 =A=J u c /u = o - A- {U/U = 4 -o-> F i g u r e 3-5 0 A comparative study showing the effect of cylinder surface during momentum injection on the variation of: (a) with a. 90 CL with a; (b) CL/CD r Figure 3-6 A comparative study showing the effect of surface roughness and momentum injection on the lift coefficient at high angles of attack. higher slope of the lift curve, increased CL trnax and slightly further delay in the onset of stall (w 32°), which continues to be gradual. The corresponding plots for the CD variation are presented in Figure 3-7. The relatively small increase in the drag for the splined cylinder case suggests effectiveness of the M S B C even at moderately high angles of attack. At an angle of attack of 30°, a reduction in the drag coefficient from CD ~ 1-6 for the smooth cylinder case in absence of momentum injection to CD ~ 0.7 for the airfoil with a splined cylinder at Uc/U = 2 is indeed impressive. It amounts to 56%. 3.2.2 Pressure distribution Extensive wind tunnel tests provided a wealth of information in terms of the pressure distribution with respect to different flow parameters. Some typical results for the smooth and splined surface cylinders are presented in Figure 3-8 to explain the physical effects of momentum injection on the flow past the airfoil. As mentioned earlier, the airfoil is symmetric about its chord and hence has similar pressure distribution on the top and bottom surfaces resulting in a zero lift at a = ,0. At a = 5° (Figure 3-8a), the pressure on the top surface of the airfoil is lower compared to that on the bottom surface. As a is increased further, the pressure on the top surface becomes even lower (Figure 3-8b). On the other hand, there is no change in the pressure distribution on the bottom surface. Also note the presence of a near stagnation value (Cp « 0.9) at the bottom of the leading edge. The airfoil has not yet stalled at a — 10°. Introduction of the momentum injection does not bring about any noticeable changes in the pressure distribution, as the flow is not separated. With a further increase in the angle of attack (a = 15°), the airfoil stalls completely in absence of the momentum injection (Uc/U increased from Uc/U — 0). As the momentum injection is — 0 to 4, the flow separation on the top surface is eliminated completely. The pressure coefficient (Cp) reaches a minimum value of -4.0 near the upper leading edge. Similar trend persists when a is further increased upto 40° (Fig92 3.0 i i i T Re =143000 — • u/u = o — c 0 ^ } U c / U = /U = 2 - - D —o I ' . " } c U 2.0 so co D 50 EfFect of momentum injection with the splined surface cylinder on the variation of Co with a. The smooth cylinder results are also presented to serve as reference. ure 3-8d to 3-8f). The peak negative Cp is about -6.0 at the top face of the rotating cylinder {Uc/U = 4), which significantly contributes to a large increase in the lift. As expected, in presence of the momentum injection, the point of separation on the top surface moves toward the trailing edge with an increase in the cylinder rotation speed. 3.3 3.3.1 Numerical Simulation Pressure distribution Figure 3-9 compares pressure distribution plots on the surface of the airfoil obtained through the numerical scheme and wind tunnel tests. The experimental data pertain to the case of the Joukowski airfoil with smooth cylinder. The wind tunnel test Reynolds number based on the airfoil chord was 143,000. The numerical results were obtained using the two element panel configuration as described earlier. At a = 5° and in absence of momentum injection {Uc/U = 0), there is a good agreement between the theory and experimental data (Figure 3-9a). However, at a higher rate of momentum injection {Uc/U = 4), the numerical scheme somewhat overpredicts the negative pressure on the top surface of the airfoil. On the other hand, presence of- a large suction peak at the rotating cylinder is successfully captured by the numerical scheme. Due to limited number of pressure taps on the surface of the rotating cylinder, the suction peak could be captured only partially during the experiment, though a tendency towards large negative C p is evident. Overall, there is a satisfactory agreement between the numerical and experimental pressure distribution plots over the complete airfoil. A n increase in a to 10° (Figure 3-9b), makes the suction peak even more pronounced. Close to the trailing edge the flow is almost separated at Uc/U = 0, indicating a tendency to stall. The complete stall at still higher angles of attack (a = 15°) is accurately captured by the numerical scheme. In presence of momentum injection (Uc/U 94 = 4), the point of separation is pushed a = 5° 0.2 0.0 _J I 0.5 0.8 1.0 -i i r a = 10° -i i i_ _i 0.0 L 0.2 _i i i i_ 0.8 0.5 i i i i- Re = 143000 • _ u /u = o _ _A- - U /U = 2 —O U /U = 4 — c i I I i i I 1.0 a = 15 c c c -o- -1 Figure 3-8 J I L_ 0.5 x/c _J I I 0.8 1.0 Experimentally obtained pressure distribution for the Joukowski airfoil with the M S B C provided by a smooth surface cylinder: (a) a = 5° ; (bya = 10° ; (c) a = 15°. 95 T 6.0 - ( d ) • P C 1 a = 20 - c £ 3.0 f - ^ ^ o . — i — i — i I I .00 6.0 O - a 3.0 4: i i 0.2 o.o i i i i • . . 1 1.0 ^-tt 1 r Re=143000 __• r°v -A- - —O— % . u /u = o U /U =2 U /U =4 c -9= 0.2 0.5 x/c a = 40 c c c -O 0.0 . 0.8 a = 30° 6 . 0 p T 3.0|\ i I -i P i 0.5 Q. -U •G 1 3 5 ^ ' 0.8 1.0 F i g u r e 3-8(cont.) Experimentally obtained pressure distribution for the Joukowski airfoil with the M S B C : (d) a = 20° ; (e) a = 30° ; (f) a = 4 0 ° . 96 further downstream ( X / C « 0.6). The effect of the momentum injection diminishes as X / C increases and becomes negligible beyond X / C « 0.6, which is expected due to the vorticity dispersion and dissipation effects. 3.3.2 C o m p a r i s o n of the various n u m e r i c a l schemes Figures 3-10 to 3-12 present comparison of results obtained using the various numerical schemes described in Section 2.3. Numerically obtained pressure plots using a single and two element panel formulations with the attached flow model are shown in Figure 3-10. The results are essentially similar except for the suction peak, which is larger for the two element case. Figure 3-ll(a) shows comparison between attached and separated flow solutions for the single element configuration. As expected, the separated flow solution agrees better with the experimental data since it is more realistic. The suction peak near the leading edge diminishes due to flow separation. The discrepancy near the trailing edge (Figure 3-10) is also eliminated. In Figure 3-ll(b), the results for the single element airfoil are compared with those obtained using the two element model. Note, the prediction by the two element case shows a slight improvement. The spike in the suction peak (which is physically unlikely) is replaced by a smooth curve indicating a gradual change in the pressure. Figure 3-12 presents results at a = 20° for the two element airfoil with attached and separated flow models. As can be expected, there is a significant improvement in prediction with the separated flow model. The unrealistically high suction peak predicted by the attached flow model is now substantially, reduced. The modelling of flow separation near the leading edge agrees more closely with the actual physical process. The discrepancies between the numerical and experimental pressure values near the leading and the trailing edges are also eliminated. 3.3.3 P r e d i c t i o n of the lift coefficient In Figure 3-13, the numerically obtained lift coefficient is compared with that 97 F i g u r e 3-9 Comparison of the numerically obtained pressure distribution with the experimental data for a Joukowski airfoil with the M S B C : (a) a = 5° ; (b) a = 10°. 98 Figure 3-9(cont.) Comparison of the numerically obtained pressure distribution with the experimental data for a Joukowski airfoil with the MSBC: (c) a = 15° ; (d) a = 20°. 99 Figure 3-10 Comparison between the results obtained through the single and two element numerical panel formulations in presence of the M S B C : (a) a = 10°, Uc/U = 0; (b) a = 10°, U /U c are also included. 100 = 1. The experimental data Figure 3-11 Comparative studies of pressure plots in presence of the M S B C : (a) attached versus separated flow solutions for a single element formulation; (b) corresponding results for the two element airfoil. The experimental data are also included. 101 15.0 T 1 r T 1 1 r T r experimental data (Re = 143000) attached flow solution (two element aerofoil) separated flow solution (two element aerofoil) X/C T T T experimental data (Re = 143000) - attached flow solution (two element aerofoil) • separated flow solution (two element aerofoil) a = 20°, U / U = 4 r -i X/C F i g u r e 3-12 0.8 i i_ 1.0 Relative performance of the attached and separated flow models as applied to the two element airfoil:(a) a = 20°, Uc/U 20°, Uc/U = 4. 102 — 2; (b) a = from the wind tunnel tests. The Cx against a curve agrees closely with the experimental results, especially at low and moderate angles of attack (ct < 20°). The airfoil stall around a = 10° is also captured by the numerical scheme accurately. However, for a > 30° the numerical results start to diverge from the experimental data. This can be expected due to relatively complex character of the flow. It indicates a need for improvement in the 'free vortex layer' model of the wake at very high angles of attack. 3.3.4 ' Iteration characteristic and the wake shape The iteration character of the numerical scheme is indicated in Figure 3-14. The results are for the two element airfoil with separated flow model in presence of momentum injection (a = 1 5 ° , Uc/U = 1). The lift coefficient converges quite rapidly to its final value in less than 15 iterations, which is considered very good. Each iteration cycle affects the shape of the wake. Initially, the flow is represented by parallel upper and lower 'free vortex layers' emanating from the separation points. This is shown in the inset box. If the wake vorticity dissipation rate (A) is considered zero then an infinite wake is obtained, which is unrealistic. In case of the nonzero dissipation rate (A = 0.62), the 'free vortex layers' converge to a closed wake in about 15 iterations. This takes about 2 minutes on a SUN SPARC 2 workstation. The effect of number of panels on the convergence of the numerical scheme was also investigated (Figure 3-15). With just 57 panels (airfoil + wake) the residual (i.e. the maximum difference in wake shape between the two consecutive iterations) is about 1 0 - 2 when a cyclic variation sets in. If the number of panels is increased to 104, the accuracy is improved to 1 0 . - 5 Finally, with 204 panels, an accuracy of the order of 10~ is achieved in 100 iterations which is sufficient for all engineering 6 purposes. This suggests sensitivity of the solution to the distribution of panels. 103 2.5 I I I I I I I I I I I I I I I I I I ] Re = 143000 - I I I SYMBOLS - EXPERIMENT LINES -THEORY 2.0 • • A }u /u = o }U /U =2 c o C 1.5 - O 0 } U /U=4 c A o 20 o 25 30 a Figure 3-13 Numerically obtained lift coefficient as compared to the experimental data for a Joukowski airfoil in presence of the M S B C . 35 1.06 U /U = 1 c a=15° 1.04 o wake shape {X = 0.0) 1.02 1.00 wake shape [X = 0.62) NUMBER OF ITERATION 5 F i g u r e 3-14 10 15 Iteration characteristic and convergence of the wake shape for a flow past the Joukowski aerofoil in presence of the M S B C . The results were obtained using the two element separated flow numerical model. 10° Figure 3-15 10 1 NUMBER OF ITERATION 10 2 Effect of number of panels on the convergence accuracy of the numerical scheme. 3.4 Flow Visualization To get better appreciation as to the physical character of the complex flow field as affected by the angle of attack and momentum injection parameters, extensive flow visualization study was undertaken. It also gave useful information about relative importance of the various system parameters and hence assisted in planning of the experiments as well as the numerical analysis. The study showed, rather dramatically, effectiveness of this form of boundarylayer control (Figure 3-16). For the leading edge smooth cylinder at a = 15° (Re = 29,400) and in absence of the cylinder rotation, a well-defined early separation, almost at the leading edge, results in a wide wake with large-scale vortices sweeping away downstream (Figure 3-16a). As the momentum injection is increased to Uc/U — 1, there is flow reattachment of the separated shear-layer at about X / C ~ 0.6 forming a large bubble with recirculating flow. However, the momentum injection is not sufficient to keep the flow attached and the flow separates again rather quickly, forming a well organized wake. Compared to the case of Uc/U = 0, the streamlines move closer on the upper surface of the airfoil, diminishing the wake width considerably ( « 0.48C). At Uc/U = 2, the streamlines are very close to the upper surface and the wake is quite small (~ 0.3C). The separation bubble has disappeared completely. Note the presence of a stagnation point near the leading edge on the bottom surface of the airfoil. The separation point is further pushed towards the trailing edge. At Uc/U = 4, an essentially attached flow is established over most of the upper surface of the airfoil, even at such a high angle of attack, considerably beyond the nominal stall angle of around 10°. Similar trends persists even at a = 45° (Figure 3-16b). The flow pattern in the wake of the airfoil was found to be quite unsteady with the vortex layer separating and forming a bubble on reattachment, the whole structure drifting downstream, diffusing, and regrouping at different scales of vortices. Ultimately the flow sheds large as well as small vortices. This unsteady character of the 107 Figure 3-16 Representative flow visualization pictures showing, rather dramatically, effectiveness of the boundary-layer control through momentum injection. The observed results compare well with the delay in separation predicted by the panel code: (a) a = 15°. 108 Figure 3-16(cont.) Representative flow visualization pictures showing, rather dramatically, effectiveness of the boundary-layer control through momentum injection. The observed results compare well with the delay in separation predicted by the panel code: (b) a = 45°. 109 separating shear layer and the wake is clearly evident in the flow visualization video. Thus the flow character indicated by the experimentally obtained time-averaged plots appear to be a fair description of the process. 3.5 Summary Based on the wind tunnel, numerical and flow visualization results it can be concluded that significant improvement in the lift coefficient as well as delay in the stall angle could be achieved by judiciously selecting the rate of the momentum injection, and the cylinder surface geometry. Among the three different surface roughness conditions studied, the cylinder with axial splines was found to be the most effective. In presence of the M S B C the stall becomes gradual and could be delayed to 50°. The numerical scheme provides good approximation of the actual fluid dynamics of the M S B C as applied to an airfoil at high angles of attack. The flow visualization results dramatically confirm the effectiveness of the M S B C in delaying separation. 110 4. F L O W P A S T A D - S E C T I O N W I T H M O M E N T U M I N J E C T I O N 4.1 Preliminary Remarks Over the years, considerable attention has been directed towards the aerodynamics and dynamics of circular cylinders and rectangular prisms. On the other hand, the D-section has received, relatively, less attention. This is surprising because it does exhibit behaviour of importance to design engineers. The D-section has a frontal flat face and a semicircular afterbody. A n interesting toy known as "Lanchesters's aerial tourbillion" [148] has a D-shaped cross-section and exhibits autorotation [145]. The D-section is susceptible to both the vortex resonance and galloping instabilities. Being a hard oscillator due to nonlinear character of the system, the D-section would not gallop from rest, but requires an initial threshold disturbance at a given freestream speed, to execute transverse vibrations. Dynamics of ice-covered transmission lines exposed to wind has been studied by employing the D-section as the likely model [149, 150]. In fact, Prince hotel in Otsu, Japan, around 135 m tall, has a cross-section of a D-geometry (Figure 4-1). With this as background, the present study focuses on: (i) Wind tunnel tests with a two-dimensional model of the D-section to assess its aerodynamics as affected by the angle of attack and momentum injection. Results show the M S B C to be effective in delaying the boundary-layer separation resulting in a significant reduction in the drag, a desirable feature from the structural design consideration. It also proves to be quite successful in suppressing the vortex shedding even with a small amount of momentum injection. (ii) Numerical simulation of the fluid dynamics of the D-section, in presence of momentum injection, using an extension of the vortex panel method. Results 111 OTSU i ® PKP^OE HOTEL Mttt3NR«> /t»-^»i ^ * * » , « » t t . *«t-r»> X * * * Hi Sir J It I w U U 4 n t O t t t < c t m : | t . R f frill'*. Figure 4-1 • %lTtfriT1r'WT~fl'»'-f^^*f-ft v The 38 storey Prince Hotel in Otsu, Japan; height 137 m. 112 compare well with the experimental data promising considerable saving in time, effort and cost associated with wind tunnel tests. (iii) Extensive flow visualization study in a closed circuit water channel using slit lighting in conjunction with polyvinyl chloride tracer particles. It shows, rather dramatically, effectiveness of the M S B C as applied to the D-section. 4.2 W i n d Tunnel Investigation Here results are presented for a D-section with smooth rotating cylinders, which also serve as the benchmark case for comparison with the numerical data. 4.2.1 Pressure distribution Experimentally obtained pressure distribution around the D-section is pre- sented for a complete range of angle of attack (0-180°) in Figure 4-2. Experiments were carried out for the rate of momentum injection ranging from Uc/U steps of one unit, however, for clarity, results are presented for Uc/U = 0 to 4 in — 0, 2 and 4. At a = 0, the pressure distribution is symmetrical about the stagnation point as expected. The stagnation point is located midway on the front flat face of the Dsection. The D-section has rounded corners due to the presence of rotating cylinders. The pressure starts decreasing gradually as the flow approaches the corners and finally it separates near the rotating cylinders in absence of momentum injection (Uc/U = 0). The semicircular afterbody lies completely in the wake with a negative pressure coefficient, which is constant everywhere on the afterbody. This indicates presence of a large wake, which is responsible for a higher drag (CD = 1.3) experienced by the D-section. As the momentum injection is increased (Uc/U = 0 —> 4) the pressure distribution is altered substantially. The rotating cylinders have a small effect on the upstream flow, except that the fluid is accelerated towards the rotating cylinder resulting in a slightly lower pressure. As seen from the pressure plot for a = 0, large 113 Q--*-Q Re = 66000 • a - = 30° _1 I L. -2.0 Figure 4-2 -1.0 j 11(711=0 - A - - —o I 0.0 1 S/R . 1.0 . . Uc/U=2 U(AI=4 . 1 2.0 Surface pressure plots for the D-section showing effects of the angle of attack and momentum injection (a = 0 - 3 0 ° ) . •114 0.0 - 0 Q—-6=- P -4.0 h a = 60 c -2.0 P— —r 1 0.0 -1.0 i I i 0.0 1.0 i i 2.0. r 'p -6.0 i Re = 66000 a = 75 -2.0 O! \I c \l i -1.0 i i -O-L i 0.0 S/R 1.0 D - -A- --O— • • - 11,711=0 11,711=2' Uc/U=4 • T - 2.0 Figure 4-2(cont.) Surface pressure plots for the D-section showing effects of the angle of attack and momentum injection (a = 4 0 ° - 7 5 ° ) . 115 i I i i r i I o.o A- -6.0 h i o.ok.t. 'P &---0- i r .p.--=3E^3£^y • -6.0 - oc = 105 c -2.0 _i i i • ' -1.0 t_ -1 1 I I I 1 0.0 T 1.0 2.0 1 1 1— I I 1 1 1 1 1 1 0.0 - -3.0 - a =135° Re = 66000 -A- •-0I -2.0 Figure 4-2(cont.) -1.0 I U(/U=2 U^U=4 I 0.0 S / R 1.0 2.0 Surface pressure plots for the D-section showing effects of the angle of attack and momentum injection (a = 9 0 ° - 1 8 0 ° ) . 116 suction peaks appear on the surface of the rotating cylinders. The strength of the suction peak grows with an increase in the rate of momentum injection. A part of the suction peak is located towards the front of the D-section, which also contributes to the drag reduction. The injection of momentum also results in a higher wake, pressure, i.e. the pressure at the base of the afterbody lying in the wake. In absence of momentum injection, the wake pressure coefficient Cps ~ -0.90, whereas at Uc/U =4, CPB ~ -0.74, i.e. a rise of about 22%, which is significant. As the angle of attack is increased from zero, the symmetry in the pressure distribution is disturbed. The stagnation point no longer remains at the center of the front face. It gradually shifts towards the rotating cylinder facing upstream. The asymmetry in the flow field results in an adjustment of pressure distribution over the whole body including the wake. As a result, the suction peaks on the rotating cylinders are no longer of equal strength. The upstream suction peak is diminished while the downstream suction peak remains almost unaffected at low angles of attack (a < 30°). Thus the movement of the stagnation point has important consequences on the overall pressure distribution. At a = 40°, the stagnation point has moved close to the upstream corner, merging with the original location of the suction peak. Note, now the suction peak is no longer visible on the upstream rotating cylinder; even at the highest.rate of momentum injection used during the experiment (Uc/U = 4). But the wake pressure still shows a positive increase in presence of the momentum injection. The suction peak on the downstream rotating cylinder progressively diminishes. A n increase in a upto 60° produces no further changes in the overall pressure distribution. For 75° < a < 105°, a large suction peak reappears at the upstream rotating cylinder due to local acceleration of the flow. The rest of the body lies completely in the wake on either side of the stagnation point. At a = 90°, there is no rise in the wake pressure with momentum injection. In fact, at a = 105°, the wake pressure becomes more negative due to adverse effect of upstream cylinder rotation, and the resulting 117 drag is higher. This indicates that for a > 9 0 ° , the direction of upstream cylinder rotation needs to be reversed. At a = 180° the semicircular afterbody of the D-sectioh faces upstream and the flat face turns towards the downstream direction. This situation is of some interest, since at a = 180°, the D-section has no afterbody lying in the wake. The rotating cylinders are now located at the trailing edges as against at the leading edge. Note, the front portion of the D-section behaves like a circular cylinder. The original scheme of cylinder rotation is obviously detrimental and a large decrease in the wake pressure coupled with a very high drag (160% increase) is observed. As seen earlier, a = 90° represents a special case. The flat face is now aligned with the flow. The upstream rotating cylinder encounters the flow first. By examining the schematic of the D-section (Figure 4-3), it is apparent that the upstream cylinder should be rotated clockwise so as not to oppose the incoming flow. It was decided to assess the effect of both clockwise and counter-clockwise rotation of the upstream cylinder. In addition, a third case involving reversal in the direction of the downstream cylinder was also studied. The results are presented in Figure 4-3. The reversal in the downstream cylinder rotation has virtually no effect since it already operates in a separated flow environment (upstream cylinder rotation counter-clockwise). The other two cases had little effect on the overall pressure distribution, though a small drag reduction was achieved for both the cases. Thus, for a = 90°, the results suggest that the D-section has essentially reached the minimum drag value and the injection of momentum has little effect on further reduction of Cp. 4.2.2 Wake pressure and vortex shedding The fluid in the wake is separated from the main flow by a thin layer of vorticity which originate on the body and detaches at the separation points. These surfaces of discontinuity represent a jump in the velocity and could be idealized as "free vortex lines" in two-dimensions. In reality, the demarcation between the wake and the outer 118 —i— — — — —r 1 0.0 1 -O O 1 1 -i 1 r- ^r- 'p -4.0 ot = 90 -8.0 c -2.0 ' ' -1.0 0.0 ' • • ^ 0.0 r -4.0 a = 90 -8.0 c _l -2.0 ROUGH CYLINDER I I I L_ -1.0 - r — i — i — i — i — i — i o.ok 0.0 i 1— I -I I I 1_ 1.0 J 2.0 L. r *0" 'p Re = 66000 -4.0 - -A- —o- a = 90° -8.0 Figure 4-3 ' -2.0 -1.0 0.0 ' • J S/R 1—1 1.0 I I I -u /u = 0 u /u = 2 u /u = 4 c c c I I Xm 2.0 Effect of three different modes of momentum injection on the surface pressure distribution at a = 90°. 119 potential flow is not so sharp as the viscosity and turbulence smear this boundary. The mean freestream velocity at separation (UB) can be approximated as [16] U = U (1-C ), 2 where Cps (4.1) 2 B PB — base pressure in the wake and U = free stream velocity. Thus the two shear layers, carrying opposite sign of vorticity, are.responsible for vortex shedding. The presence of a bluff body merely modifies the process by allowing a feedback mechanism between the wake and the shedding of circulation at the separation points • [16]. The rate at which the circulation (T) is shedded at the separation points can be approximated as [16] * = f- ( 4 ' 2 ) It is also known that vortex formation involves mixing of flows of oppositely signed vorticity [96]. The strengths of individual wake vortices will be less than the total circulation shed from one side of a bluff body during a shedding cycle. Hence the fraction of the original circulation that survives vortex formation is [16] , 6 = 2 - ^ ( l - C P B ) - \ • (4.3) where St = Strouhal number, I V = strength of a wake vortex, and D = characteristic dimension of the body. • The value of 6 has been determined mostly through experiment and range from 0.2 to 0.6 for different bluff body shapes [6, 16, 151-153].' Even for a given bluff shape, value of 8 is not consistent [16]. The set of equations (4.1) to (4.3) emphasize the strong link that exists between the vortex shedding and base pressure in the wake. A higher base pressure, in general, would correspond to a lower rate of shedding. This fact could be inferred qualitatively from the present experiments with momentum injection. Though no quantitative estimates were sought, it was obvious 120 that momentum injection disrupts the vortex shedding process and, at a sufficiently high rate, it completely suppresses it. Variation of the base pressure as affected by the angle of attack and momentum injection is shown in Figure 4-4. Upto a = 105°, pressure tap A (as shown in the schematic of the D-section) was found to lie in the wake. For a = 135° and 180°, pressure tap B was used for the purpose. The results show that the base pressure is much lower than the freestream value (Cpg < 0), and increases significantly due to momentum injection. On average, almost 30% rise in the base pressure was achieved for 0 < a < 9 0 ° . At a = 90°, the rise in base pressure was about 8%. For a > 90°, the wake pressure reduces, indicating a need for the reversal in the direction of rotation of the upstream cylinder. For a > 135°, it becomes necessary to reverse the direction of both the rotating cylinders (i.e. inward mode of rotation). This raises the wake pressure by 10 - 25 %. In Figure 4-4(b), variation in the Strouhal number is plotted against the angle of attack. The Strouhal number remains essentially unchanged except in the range 60° < a < 105°, with a peak at a = 90°. The flat face of the D-section is now aligned with the flow thus making the body relatively streamlined and causing a jump in the vortex shedding frequency. As the momentum injection is increased from Uc/U = 0 - 4, a dramatic change occurs in the wake. At a low rate of momentum injection (Uc/U « 1), the vortex shedding process was extremely weak and hence no clear frequency could be identified through spectral analysis of pressure signal from the surface of the body. Even a disk probe (Section 2.2.3) in conjunction with an extremely sensitive Barocel pressure transducer (resolution 1 0 - 6 l b / i n , 6.9 x 2 10~ N / M ) as well as a hot wire anemometer system (DISA Elektronik, Denmark), 5 2 both traversing the wake to identify discernible signals, failed to detect a coherent vortex shedding frequency. Perhaps the small amount of vorticity shed was quickly dissipated in the turbulent wake (Re « 66,000). 121 Of course, with the momentum 0.5 = Strouhal Number (St). r T Re •= 66000 r .(b) u /u = o c 0.2 ~ JO -o- "a— "-o _ _ } St (Based on diameter) _0—St,p (Based on projected height) "T • A • A 0.0 A ^ W * L & - A - L ^ _ J L A - A - L . U /U = 2 A 120 60 Figure 4-4 c 180 Wake pressure and the Strouhal number for the D-section as affected by the momentum injection. 122 injection, the shed vorticity diminishes in strength, and for Uc/U > 2 it was virtually eliminated with the Strouhal number designated as zero (Figure 4-4). The changes in the wake brought about by momentum injection can be better appreciated through flow visualization photographs, which are presented later in this chapter. The elimination of vortex shedding is also indicated by a higher base pressure in the wake. A n increase in the base pressure can be attributed to a smaller width of the wake and a lower drag coefficient. The rise in the base pressure is similar to that obtained when a splitter plate of appropriate length is placed behind a circular cylinder at a suitable location [151]. The momentum injection is responsible for delaying and even suppressing the boundary-layer separation. 4.2.3 A e r o d y n a m i c coefficients Variations in the drag and lift coefficients with the angle of attack are shown in Figure 4-5. The net drag reduction is summarized in Table 4.1. The drag reduction (%ACx>) is defined as %AC = ( C/V=°) Cd u D ~ °Vcl™) C (u /u=o) C D x 1 0 0 . ( 4 . 4 ) c The direction of rotation is kept fixed as at a = 0 throughout. Table 4.1 Drag reduction with momentum injection for various a a° 0 30 45 60 90 105 180 C (U /U = 0) 1.293 0.997 0.806 0.583 0.209 0.306 0.130 D c %AC D (Uc/U 42 19 22 29 =4) 75 -13 -158 The drag coefficient reduces with the angle of attack and reaches a minimum at a = 9 0 ° . In presence of momentum injection the drag coefficient is lowered considerably. The drag reduction continues for 0 < a < 105°. Beyond a > 105°, the 123 need for reversal in the direction of cylinder rotation is clearly indicated. Note, the minimum drag coefficient at a = 90° is also coupled with the highest vortex shedding frequency (Figure 4-4). The decrease in drag by as much as 75% (a = 9 0 ° ) is indeed impressive. The absolute change in the lift coefficient by itself is of little significance for a bluff body like a D-section, but the rate of change of lift coefficient, dC^/da, is quite important in determining the susceptibility to flow-induced vibrations. A negative dCjj/da indicates a likelyhood for galloping instability at low angles of attack [154]. Figure 4-5 clearly shows dC^/da < 0 at Uc/U = 0, but becomes positive for Uc/U > 2 indicating that the stability is restored. This results are discussed in more detail in Chapter 6 in connection with dynamic experiments. At a = 0 and 180°, the D-section is symmetric with respect to the flow direction and the lift coefficient is zero. In the range a = 40° - 45°, the lift coefficient is negative but large in value (CL ~ -1.0) in absence of momentum injection. However, with Uc/U = 4, the lift coefficient is increased substantially (CL « - 0.55, i.e. the absolute value is reduced). Thus in general, momentum injection reduces transverse load on the structure. This could be used to advantage in their design. Figure 4-6 shows effect of various strategies of cylinder rotation on the drag coefficient. As pointed out earlier, for a < 105°, the outward rotation mode is most effective in reducing the drag. For a = 135°, the direction of rotation of the upstream cylinder needs to be reversed (i:e. the downstream rotation mode). For a > 135°, both the cylinders are affected (i.e. inward rotation mode) to obtain drag reduction. In practice, for a tall building , this can be accomplished quite readily through a microprocessor based switching actuator operating in conjunction with a wind direction transducer. 4.3 N u m e r i c a l Simulation A comparison between numerically and experimentally obtained pressure plots, 124 90 45 1.0 T 1 1r T 1 1- 180 135 i 1 ~i r 1—:—r Re .= 66000 —a— u /u = o • " _ c _ _A- -2.0 Figure 4-5 _i i i _i i_ 45 i i_ _l 90 I I a - U ^ U = 1 — V — • Uc/U = 2 ....<>... Uc/U = 3 -O-- Uc/U = 4 L. 135 180 Drag and lift characteristics of the D-section as affected by the momentum injection 125 DOWNSTREAM ROTATION OUTWARD ROTATION Re = 66000 - • — u /u = o c 180 INWARD ROTATION Figure 4-6 The effect of three different modes of momentum injection on the drag coefficient of the D-sectioh as a function of angle of attack. 126 as presented in Figure 4-7, shows good agreement. At a — 0, the numerical scheme predicts separation at the corners and a constant pressure wake. It successfully captures the existence of suction peaks, and predicts a rise in the base pressure in presence of momentum injection. For a > 0, an asymmetric pressure distribution coupled with the movement of the stagnation point matches with that obtained from the wind tunnel tests. At higher a, the pressure distribution on the upstream face of the D-section agrees more closely with the wind tunnel results, whereas in the wake, the difference between the theory and experiments is relatively higher, though the trend is quite similar. This, of course, is understandable due to potential character of the panel method. Figure 4-8(a) compares the corresponding drag data. The numerical scheme accurately predicts drag reduction in presence of momentum injection. The comparison was extended upto a — 40° to assess the limit of applicability of the panel method used here. For a > 30°, discrepancy between the theory and the experiment is detected. This suggests that, at high angles of attack, the panel method needs an improved model for the wake. Figure 4-8(b) shows numerically obtained wake geometry at a = 0. In absence of momentum injection, a wide wake results. It should be pointed out that only the near wake region (upto X / R = 10) is shown here. The actual wake continues for many diameters downstream, with the 'free vortex lines' drawing near each other and eventually closing. A finite wake is obtained due to the dissipation of the vorticity modelled in the panel code. The vorticity dissipation factor determines the wake length and is an empirical input to the program. The value of the base pressure in the wake was found to be a function of the vorticity dissipation factor. As the momentum injection is increased from zero, the "free vortex lines" converge toward each other resulting in a smaller wake width. This feature of the numerically obtained wake shape is in good qualitative agreement with the flow visualization results. 127 i i a =0 0.0 C'P WftLi-frfr p -2.0 -4.0 h . I -2.0 2.0 T o.o ffe a=10° . ^-6^ A ^ ^ „ „ • „ A..A, A,A.A.-A | 'A / i Ai. i -.3.0 A / A / — / 1 \ • Re = 66000 • }u /u = o c -6.0 LINES-THEORY - SYMBOLS - EXPERIMENT -I -2.0 Figure 4-7 A I I _l u 1 . l_ -1.0 0.0 S/R 1.0 }U /U = 4 C I • 2.0 Comparison between the numerically and experimentally obtained pressure plots for the D-section. Note, the numerical scheme accurately predicts presence of large suction peaks, as well as rise in the wake pressure in presence of momentum injection (a = 0, 10°). 128 T — i — i — I — r r a = 20 o.o P . i -2.0 - -4.0 - c ' -2.0 2.0 Re = 66000 -6.0 ;/U = 0 LINES-THEORY • - SYMBOLS - EXPERIMENT -2.0 A ju /y =4 r 0.0 -1.0 S/R 1.0 2.0 Figure 4-7(cont.) Comparison between the numerically and experimentally obtained pressure plots for the D-section. Note, the numerical scheme accurately predicts presence of large suction peaks, as well-as rise in wake pressure in presence of momentum injection ..(<* = 20°,.30°). 129 2.0 T 1 1 T" LINES-THEORY SYMBOLS -EXPERIMENT D }u /u = o c _ _ 0 _ }U /U=4 C (b) 2.5 L L / U - 0 I W U = 1 Y/R 0.0 -2.5. 0.0 Figure 4-8 2.5 5.0 X/R 7.5 10.0 Numerical simulation data showing the effect of momentum injection: • (a) drag coefficient;'(b) wake geometry as represented by the 'free vortex lines'. , 130- Thus even a simple panel code in conjunction with vorticity dissipation gives surprisingly good results, even for such a complex flow phenomena as represented by a bluff body with momentum injection. 4.4 Effect of Surface Condition of a Rotating Cylinder A change in the surface conditions of the rotating cylinder could result in important benefits in terms of more efficient momentum injection. As described in Chapter 3, a splined rotating cylinder-was found most effective in increasing the lift of an airfoil. Hence it was decided to assess the relative performance of cylinders with different surface conditions (smooth, rough and splined). Figures 4-9 and 4-10 present relative performance of the three different rotating cylinders in terms of the drag and lift coefficients, lift/drag ratio and the base pressure coefficient. In absence of momentum injection (Uc/U = 0), the splined and rough rotating cylinder have a higher drag due to wider wake caused by early separation. However, the surface roughness does improve efficiency of the momentum injection process resulting in further decrease in the drag and increase in the base pressure (compared to the smooth cylinder case), particularly for a < 20° (Figure 4-9). On the other hand, so far as the lift is concerned, the cylinder surface roughness has little effect even upto a as large as 30° (Figure 4-10). However, beyond that there are some discernible differences which can be better understood by referring to the pressure distribution plots (Figure 4-11). At a = 0, in absence of momentum injection, the base pressure is lower for the rough and splined rotating cylinders, which explains the measured higher drag coefficient. However, with Uc/U = 4, the wake pressure rise is more for splined and rough cylinders, and strength of the suction peak is slightly diminished ( S / R = l ) . The resultant effect is the decrease in drag. At a = 4 5 ° , the strength of the suction peak at S/R = -1 is significantly increased with the momentum injection leading to a higher C^. 131 0 o — I — — ' — — ' — — ' — ' ' '—I 1 0 1 1 10 1 I L — 20 i I 30 I I I I I 40 I I • I I 50 I • . . I 60 •J 0 Figure 4-9 10 20 30 40 50 60 Effect of cylinder surface condition on: (a) drag coefficient; (b) wake pressure. 132 134 Re = 66000 - SMOOTH - ROUGH _ SPLINED o / - -A - —•o-8.0 -2.0 a = 45° 0' A I -1.0 • 0.0 g • / R • 1.0 2.0 Figure 4 - l l ( c o n t . ) Effect of cylinder surface characteristic on pressure distribution over the D-section: (b)a = 45°. 135 4.5 Flow Visualization Flow visualization has always been a powerful tool in studying complex fluid dynamics phenomena. A set of revealing photographs can convey important information about the physical character of the flow, at least qualitatively. Thus it helps in understanding of possible mechanisms at work that characterize the flow field. Depending on the accuracy of the flow visualization process and parameters of interest, even quantitative information can be obtained through analysis of photographs and video movies. The details of flow visualization procedure adopted was described in Section 2.2.6. Here, results for the D-section are presented (Figures 4-12 and 4-13). At a = 0 and in absence of momentum injection, the first thing that captures the attention is the existence of a wide wake behind the D-section (Figure 4-12). Note, the wake is significantly wider than the diameter of the D-section. The Reynolds number was 20,000 and the boundary-layer at separation is turbulent. The wake vortex structure displays considerable turbulence, with associated mixing and dissipation. Immediately behind the D-section, there is a small region of essentially stagnant fluid. The separating shear layer from the top has a sinusoidal character and divides the wake into well-defined pockets of diffused vorticity, generated from the viscous dissipation of the shed vortices. The overall configuration of the wake t reverses itself in an alternate manner in phase with the vortex shedding process. The video revealed considerable unsteadiness of the wake, oscillating from side to side due to the low pressure region associated with the nascent vortex. At a given instant, the wake is slightly asymmetric (Figures 4-12 and 4-13). The shed vortices grow rapidly as they are convected downstream, diffusing vorticity in the surrounding fluid at a faster rate. With an increase in the momentum injection from Uc/U = 0 —• 4, the wake becomes progressively narrower, turning into a small bubble at Uc/U = 4. This suggests considerable influence of the momentum injection on the inclination of the 136 a =0 Uc/U = 0 Uc/U = 2 CO Figure 4-12 Flow visualization photographs showing effectiveness of the Moving Surface Boundary-layer Cor (MSBC) as applied to the D-section. potential character. Note, at Ug/U = 4, the fluid motion tends to approach Uc/U Uc/U =0 D-SECTION = 0 S E P A R A T E D SHEAR L A Y E R RAPIDLY DIFFUSING VORTEX A CLOSED WAKE GROWING VORTEX RELATIVELY SMALL M O M E N T U M INJECTION ELEMENT WAKE BUBBLE T U R B U L E N T WAKE WITHOUT V O R T E X SHEDDING Uc/U = 2 Figure 4-13 Uc/U = 4 Effect of momentum injection on the characteristic features of the wake and the shed vortex system, as inferred from the flow visualization study. separating shear layers. For example, the separation angle at the top corner of the D-section is around 10° at Uc/U — 0 (measured counter-clockwise with respect to the flow direction), changes to « - 1 0 ° at Uc/U = 2 and reaches a value of approximately - 3 5 ° at Uc/U = 4 thus resulting in a progressively smaller cavity in the wake. The oscillatory behaviour of the wake manifests itself in a much weaker form further downstream. At low rates of momentum injection {Uc/U < 1), the video clearly showed a significant increase in the vortex shedding frequency. However, at sufficiently high rates of momentum injection (Uc/U > 2), vortex shedding is suppressed completely and the fluid motion approaches the potential flow character. Since the "bluffness" of a body and the width of the wake are related, a reduction in the wake-width implies, effectively, a more streamlined body with an associated reduction in the drag. Number of researchers [153, 155, 156] have pointed out the discrepancy between the experimentally observed vortex spacing ratio b/a (i.e. transverse to longitudinal spacing of a staggered row of vortices) and that predicted by Karman's stability analysis (b/a fa 0.281). A close scrutiny of a large number of flow visualization photographs revealed considerable scatter (b/a = 0.23-0.58 for Uc/U = 0 and depending on a) away from the value predicted by Karman. The Reynolds number being sufficiently high ( « 20,000), the separation is turbulent. This could introduce intermittency and randomness in the vortex formation. Three dimensional effects are also known to alter the vortex spacing ratio. At a = 45° (Figure 4-12) the essential features of the wake remain similar. Even at Uc/U = 0, the wake is narrower compared to that for the a = 0 case, and this explains the lower drag coefficient. In presence of momentum injection the bubble containing separated recirculating fluid is nearly eliminated. Note, the freestream leaves the body at an angle due to the momentum injection. 4.6 Summary The D-section in presence of momentum injection shows significant decrease 139 in drag as revealed from the wind tunnel, numerical and flow visualization results. The angle of attack beyond a certain range necessitates reversal in the direction of the momentum injection for one or both the cylinders. The reversal in the direction of momentum injection is governed primarily by the location of the upstream stagnation point. A small increase in the drag reduction can be achieved by changing the cylinder surface condition. Vortex shedding is completely eliminated even at a moderate rate of momentum injection (Uc/U > 2). A reduction in the drag as well as the lift coefficient in presence of momentum injection coupled with the dynamically more stable configuration, free from vortex resonance and galloping, are features quite attractive to design engineers. 140 FLOW PAST R E C T A N G U L A R PRISMS W I T H M O M E N T U M INJECTION 5.1 Preliminary Remarks Rectangular prisms represent a family of two-dimensional bluff bodies defined by a non-dimensional parameter called the aspect ratio (length to height ratio abbreviated as A R ) . This is not the case with a circular cylinder or a D-section, the geometries which remain unaffected by changes in their dimensions. By varying the aspect ratio of a rectangular prism, a variety of aerodynamic phenomena, like complete flow separation (AR < 3), reattachment ( A R > 3), galloping instability (0.8 < A R < 2.5) and vortex resonance (AR > 0) can be studied. A rectangular prism with A R « 0 is the familiar normal flat plate, while A R —• oo (i.e. for very large values of aspect ratio, say 20) results in a flat plate.aligned with the flow. Rectangular shapes are encountered frequently in actual engineering practice, e.g. cross-sectional shapes of majority of tall buildings and ground vehicles (tractor-trailer trucks; buses and trains). From the structural dynamics point of view, a varying aspect ratio implies changes in the stiffness (longitudinal, transverse and torsional); natural vibrational frequency; load bearing capacity and internal damping. Though much attention has been directed towards the basic fluid mechanics of rectangular prisms, their performance in presence of the M S B C has received little attention. Important questions regarding the effectiveness of momentum injection in terms of the overall drag reduction, suppression of vortex shedding, susceptibility to galloping type of wind-induced instability as well as effect on the wake structure, remain unanswered. Two representative rectangular bluff body models, a vertical rectangle ( A R = 0.5) and a square prism (AR = 1.0), were subjected to extensive wind tunnel investigation. The rectangular prism at a = 90° is aligned with the flow and has A R = 2.0. The wind tunnel results in terms of pressure distribution, drag and lift 141 coefficients, wake pressure and the Strouhal number are discussed as affected by the angle of attack and the rate of momentum injection. The aerodynamic coefficients are also utilized for quasi-steady analysis to assess susceptibility of the prisms to galloping type of instability. A n attempt is made to arrive at a universal Strouhal number. The chapter ends with flow visualization results which reveal intricate character of the flow field and the effect of momentum injection on it. 5.2 5.2.1 W i n d Tunnel Investigation Pressure distribution The surface pressure distribution results as affected by the angle of attack and rate of momentum injection are presented in Figures 5-1 and 5-2 for the rectangular and the square prisms, respectively. The angle of attack (a) was varied from 0 to 180° for the rectangular prism and from 0 to 105° for the square prism. Experiments were carried out for the rate of momentum injection ranging from Uc/U = 0 to 4 in steps of one unit as before. Only some typical results useful in establishing trends are presented here for conciseness and clarity. At a = 0, the pressure distribution is symmetric about the stagnation point. The entire^ front faces of the prisms are stagnation zones. Note, both the prisms have rounded upstream corners due to presence of the rotating cylinders. The pressure starts decreasing, at first gradually, as one moves away from the mid-point, followed by a rapid reduction on the surface of the rotating cylinder resulting in large' suction peaks. In absence of momentum injection, the flow separates near the corners and •the rest of the afterbody lies in the wake. The pressure coefficient (Cp) in the wake is negative and remains nearly constant on the surface of the afterbody. The presence of a large wake is primarily responsible for the high drag coefficient (CD « 1.2) experienced by both the rectangular and the square prisms. Note that the pressure plots for the rectangular and the square prisms at a = 0 are almost identical indicating 142 - i — i — i — r 0.0 - •o- — • — -a— 4.0 - CX — 0 • . . . i -2.0 . . . . i •1.0 0.0 1.0 2.0 0.0 -6.o -ra = 15 c -2.0 0.0 Figure 5-1 S / R 1.0 Surface pressure plots for the rectangular prism showing effects of the angle of attack and momentum injection: a = 0, 15°, 30°. 143 -4.0 cc = 7 5 l o -2.0 L I -1.0 0.0 -i—i—i—|— 1 -I o.od-—<>--• v . X I I I I I I I I L 2.0 —I—i—i—i—i |—i—i—i- I , &^3^=^=? \A- '"Y '&-~ ® _~ C ^ 7Z r < J l A" ^ Re =52000 p I 1 -4.5 h -J I 1.0 11 I I cc = 90 • — u /u =0 - - u /u =2 -,-o- — u /u =4 c A l c c _. -2.0 Figure 5-l(cont.) . I I l _ -1.0 o.o 5/R 1.0 2.0 Surface pressure plots for the rectangular prism showing effects of the angle of attack and momentum injection: a = 45,° 75°, 90°. 144 I I I I 1- 0.0 ^ N / 'p -6.0 if V a = 105° Re = 52000 — u /u =0 - - u /u =2 — o - — u /u =4 - - A c c c -1.0 -2.0 F i g u r e 5-l(cont.) 0.0 S / R 1.0 1 I I l_2.0 Surface pressure plots for the rectangular prism showing effects of the angle of attack and momentum injection: a = 105°, 135°, 180°. 145 that in absence of reattachment, it is frontal geometry that primarily determines the aerodynamics of a bluff body. As the momentum injection is increased from Uc/U = 0 —> 4, the pressure distribution is altered substantially. The rotating cylinder has a negligible effect in the upstream direction (i.e. on the stagnation zone), except for the acceleration experienced by the fluid in its vicinity. The strength of the suction peak grows with an increase in the rate of momentum injection. A part of the suction peak located towards the front of the prism may also contribute to the drag reduction. Of course, as can be expected, injection of the momentum by the rotating cylinders results in a higher near wake pressure, i.e. the pressure on the base of the afterbody in the wake. At a — 0, the rise in the wake pressure is about 33% and 20% (Uc/U = 4) for the rectangular and square prism, respectively. This rise in the wake pressure was reflected in the lower drag coefficients for both the prisms. As the angle of attack is increased from zero, the symmetry in the pressure distribution is no longer present. The stagnation point starts moving towards the rotating cylinder facing upstream. The asymmetry in the flow demands an adjustment in the pressure distribution over the whole body. The suction peaks on the rotating cylinders are no longer of equal strength: the one on the upstream rotating cylinder is weaker while the downstream suction peak becomes more pronounced (0 < a < 30°). Thus the movement of the stagnation point has important consequences on the overall pressure distribution. For a w 15° - 20°, there are two small stagnation regions, on either side of the rotating cylinder, for both the rectangular and the square prisms. In case of the rectangular prism this phenomena persists upto a = 45°. For a > 20°, the stagnation point moves very close to the upstream rotating cylinder and merges with the suction peak, causing the upstream suction peak to disappear even at the highest rate of momentum injection (Uc/U = 4). For 20° < a < 40°, a significant difference develops between the surface pressure character of the rectangular and square prisms. 146 Figure 5-2 Surface pressure plots for the square prism as affected by the angle of attack and momentum injection: a = 0, 20°, 4 5 ° . 147 I I r 0.0 a = 75° -3.0 -I -2.0 I •1.0 l_ • ' ' 0.0 1.0 2.0 0.0 1.0 2.0 I II I 1— — 1 0.0 - 30 "a = 9 0 ° I I I I -2.0 I • • •1.0 _l I I L 0.0 - *. o Re = 52000 -a— - u / u = o c " 3 0 Ha = 1 0 5 ° c -i -2.0 - U /U=2 U /U = 4 i i I -1.0 i i i -•-« i_ o.o 5/R i.o • • r • i 2.0 i . . Figure 5-2(cont.) Surface pressure plots for the square prism as affected by the angle of attack and momentum injection: a = 7 5 ° , 9 0 ° , 105°. 148 u r § ) + S/R ' ' - S / R ( g R OUTWARD ROTATION T1 1 ~ DOWNSTREAM ROTATION 11 -i r- 'P -4.0 <c> -1.0-I -2.0 1 0.0 1 I I I 1 I0.0 _J • ' , 1.0 2.0 I f— -1— — 1 0Re = 52000 -4.0 h i —•—u /u = o c a = 90° i i i -2.0 U /U = 2 - - -A i i c U /U = 4 Q i •1.0 0.0 S / R 1.0 c 2.0 (a) Figure 5-3 Effect of reversing the direction of upstream cylinder on the surface pressure distribution at a = 90°: (a) rectangular prism. 149 u a u a a R R OUTWARD ROTATION -i 1 D O W N S T R E A M ROTATION 1 r 1 1 1 1 . • r- 1 — —'— —I— — 1 1 1 0.0 -2.0 _i_ -2.0 -1.0 0.0 1.0 2.0 - i — i — | — i 0.0 'P -4.0 O. i o 3v o o a =90° • -2.0 Re = 52000 o— - U / U = 0 A- - - u / u = 2 0—• - u / u = 4 .' • • • I i o.o •1.0 S i i_ i / R j 1.0 i_ c c c 1 ' 1 2.0 ' 1 (b) Figure 5-3(cont.) Effect of reversing the direction of upstream cylinder on the surface pressure distribution at a = 90°: (b) square prism. 150 As against the rectangular prism where the wake pressure rose uniformly for Uc/U > 0, with the square'prism, the wake pressure rises only on the upstream side (S/R > 0) and becomes lower on the downstream side (S/R < 0). This is reflected in the minimum reduction in the drag of ~ 21% compared to the average drag reduction of « 36% (0 < a < 75°) for the square prism. For a > 4 5 ° , the situation near the upstream rotating cylinder becomes quite complicated. The stagnation point has clearly moved to the side face for both the rectangular and the square prisms. Also, multiple regions of primary and secondary separations as well as partial reattachment are present, as indicated by the irregular pressure plots. Similar trend persists in presence of the momentum injection, except for variations in the magnitude. In case of the rectangular prism, the suction peak on the downstream rotating cylinder maintains itself even at a as high as 90°, however, this is not the case for the square prism. For a = 90° and 105°, almost the entire body lies in the wake and rotating the cylinders in the original direction does not make any difference. Note, the pressure distribution on the square prism at a — 90° is quite similar to that at a = 0 in absence of momentum injection, except for a shift of the stagnation point to the upstream facing side. It is of interest to recognize that, at a = 9 0 ° , the upstream cylinder is rotating in the opposite sense to the direction of the flow. This suggests a need for reversal in the direction of rotation. Figure 5-1 also shows pressure results for the rectangular prism in the range a > 135°. Clearly by not reversing the direction of momentum inection there is an adverse effect on the overall pressure distribution. The pressure in the wake becomes even more negative in presence of the momentum injection. This causes a large increase in the drag coefficient. Direction of rotation of both the cylinders should now be reversed in order for the momentum injection to be effective. The a = 90° is a special case. The face carrying the rotating elements is now 151 aligned with the flow. Note, the aspect ratio of the rectangular prism is now 2. It is apparent that the upstream cylinder should be rotated clockwise in order not to oppose the incoming flow. It was decided to test both the clockwise and counterclockwise rotation of the upstream cylinder. The results are shown in Figure 5-3. It is of interest to point out that reversal in the downstream cylinder rotation had virtually no effect, since it is submerged in the separated flow field. For both the prisms, the pressure distribution is altered in such a way that a net drag reduction is achieved by reversing the upstream cylinder rotation. 5.2.2 A e r o d y n a m i c coefficients Variation in the drag and the lift coefficients with the angle of attack and rate of momentum injection are shown in Figures 5-4 and 5-5 for the rectangular and the square prisms, respectively. The net changes in the drag ( A C p ) and the lift ( A C x ) are given in Tables A C J D and % ACL 5.1-5.4 for both the rectangular and square prisms. Note, the % are defined as: %AC = Vclv=») ~ W=*) CD D C x 1 Q 0 . ( 5 1 ) C {U IU=0) D %AC L = — v C L C ' CL(U /U=O) '- x 100. (5.2) C Table 5.1 Drag reduction for the rectangular prism with the momentum injection at various a. a 0 C {U /U D = C 0) %AC (U /U D = C 0 1.234 61 30 1.528 46 45 1.436 27 60 1.143 25 75 0.527 18 90 0.334 -15 135 0.540 -24 180 1.551 -17 152 4) T a b l e 5.2 Change in the lift of the rectangular prism with the momentum injection at various a. a° L C (U /U L 0 30 45 60 75 90 135 180 T a b l e 5.3 = 0) C %AC (U /U L 0.0 -0.310 -0.415 -0.539 -0.221 -0.113 -0.751 0.0 C = 4) 0 542 260 232 270 143 -1 0 Drag reduction for the square prism with the momentum injection at various a. a° C {U /U D C 0 30 45 60 75 . 90 105 T a b l e 5.4 = 0) %AC (U /U D 1.201 1.755 1.846 1.568 1.272 1.231 1.107 C = 4) 61 24 26 40 26 -10 -1 change in the lift of the square prism with the momentum injection at various a. ct° 0 30 45 60 • 75 90 105 C (U /U L C 0.0 0.304 0.101 0.031 0.093 0.0 -0.709 = 0) %AC (U /U L C = 4) 0 210 377 1213 1159 0 -5 For 0 < a < 75°, a significant decrease in drag is achieved for both the. rectangular and the square prisms. The drag reduction was found to increase directly with the rate of momentum injection. As explained earlier, at a = 9 0 ° , it necessary to reverse the direction of the upstream rotating cylinder. The projected height of the rectangular prism is minimum at a = 90° and the drag coefficient attains a minimum value (Figure 5-4a). In case of the square prism the projected height is maximum 153 154 a u a u u U u ^r-r^ c R R OUTWARD ROTATION 2.0 -i a 11 R DOWNSTREAM ROTATION INWARD ROTATION r -~0 -J I J 1 I I I l_ 0.0 20 - 40 60 - DOWNSTREAM. ROTATION u /u = o c - -o- ) — J U /U = 3 c 20 F i g u r e 5-5 100 Re = 52000 A — — 80 -I I 40 i i I I I I I I i_ 60 a ° 80 100 The drag and lift characteristics of the square prism as affected by momentum injection 155 at a. — 4 0 ° - 4 5 ° , and the drag coefficient is the highest (Figure 5-5a). Also, the drag and the lift coefficients at a = 0 and 90° remain essentially the same for the square prism as expected. For 0 < a < 75°, an average reduction in drag was around 36% for both the rectangular and the square prisms, which is indeed remarkable. For 90° < a < 135°, introduction of the downstream rotation is necessary to establish a favourable-condition for drag reduction. For a > 135°, the inward mode of rotation (i.e. reversal in the direction of rotation of both the upstream and downstream cylinders) is required to assure a decrease in drag. For a — 0, 90° and 180°, the rectangular and the square prism are symmetric with respect to the flow direction and the lift coefficient is zero. But for 0 < a < 75°, a large increase in the lift is observed as the momentum injection is increased from Uc/U — 0 to 4 (Figure 5-4b and 5-5b). This is also evident from Tables 5.2 and 5.4. This means a substantial increase in the load acting transverse to the wind.direction and should be taken into consideration while designing a real life structure. But the maximum value of \ C L \ at Uc/U = 4 is not greater than \CD\ at Uc/U = 0 and the structure should be able to withstand the load. A negative dCL/da indicates a likelihood for galloping instability at small angles of attack. A n explanation of the criterion for stability is described in Section 6.3.1. At small angles of attack (a < 20°), dCi/da (Uc/U 5.2.3 < 0 (Uc/U = 0) and becomes positive in presence of momentum injection > = 2) indicating that the stability is restored. W a k e pressure and Strouhal number A strong link between the base pressure and the vortex shedding was emphasized in Section 4.2.2. Based on the experience gained through the study of D-section, it was expected that, for sufficiently high rates of momentum injection (Uc/U > 2), vortex shedding would be eliminated completely. Dependence of the base pressure on the angle of attack and rate of momentum injection is shown in Figure 5-6 (rectangular prism) and Figure 5-7 (square prism). 156 Pressure tap - A was found to lie in the wake for a < 105° as shown in the schematic diagrams at the top of the figures. For a > 135°, pressure tap - B was used for the purpose. In general, the wake pressure coefficient is much lower than the freestream value (CPB < 0), and increases significantly in presence of the momentum injection. In case of the rectangular prism, an average increase of 24% in the wake pressure was observed for 0 < a < 75° and Uc/U = 4. Similarly, for the square prism, the average rise in the wake pressure was around 23% under similar conditions. If one continues to rotate the upstream cylinder in the counter-clockwise direction, the wake pressure becomes lower in presence of momentum injection at a = 90°. Obviously, this has an adversed effect on the drag of the body. By adopting the downstream rotation mode, the situation is improved significantly in the range of 75° < a < 105° (Figures 5-6a and 5-7a). Figures 5-6(b) and 5-7(b) show variation in the Strouhal number (St) for the rectangular and the square prisms, respectively. In case of the rectangular prism, an interesting feature is the sudden jump in the vortex shedding frequency at a = 90°. In fact this is very similar to that observed for the D-section (Figure 4-3b). The D-section and the vertical rectangular prism have a measure of similarity in this configuration resulting in similar Strouhal number characteristics. At a = 9 0 ° , the rectangular prism is aligned with the flow direction and behaves as a relatively streamlined body. Hence a jump in the vortex shedding frequency can be expected. The behaviour of the square prism is somewhat different. It shows several small peaks over a range of the angle of attack investigated (0 < a < 105°). Since the rectangular and the square prisms have well defined side faces (unlike the D-section), it was easier to capture the vortex shedding frequency using the pressure taps.located on the sides. As evident from the results, a small but consistent increase in the Strouhal number is detected in presence of the momentum injection (Uc/U = 1). But the amplitude of the dominant vortex shedding frequency is much smaller as obtained from the 157 Figure 5-6 Effect of momentum injection on the wake pressure and Strouhal number of the rectangular prism. 158 Figure 5-7 Variation in the wake pressure and Strouhal number for the square prism as affected by the momentum injection. 159 power spectral analysis. At Uc/U = 2, no vortex shedding could be detected. As confirmed from the flow visualization, at Uc/U > 2, the vortex shedding is completely eliminated. 5.3 5.3.1 Universal Strouhal Number Correlation between some nondimensional numbers The problem of flow past a bluff body has eluded complete theoretical treat- ment so far. Notable advances in this direction were made by Strouhal [1], Karman [3] and Roshko [69,151,158-159]. In recent times, the advent of numerical computation procedures have, to some extent, overshadowed the value of simple, elegant, carefully designed experiments. It is well known that all bluff bodies are characterized by wide wakes, vortex shedding and lower than freestream pressure in the wake. As remarked by Roshko [159], "the wake structure of different bluff bodies is similar. That is, the shape, or bluffness, of the body has no characteristic effect on the wake other than to determine its geometrical and velocity scales". As seen from the results of bluff bodies under study, the vortex shedding frequency and the wake pressure are functions of the angle of attack. Hence it is logical to look for a parameter which can be used to compare the wakes of different bluff bodies, or a single bluff body at different orientations. Furthermore, it would be of interest to assess the effect of momentum injection on such a parameter. Attempts to formulate a universal Strouhal number are not new. The present approach is based on the contributions by Roshko [159] and Bearman [155]. The wake of a given bluff body can be modelled as two shear layers, a distance h! apart with the velocity outside the layers equal to Z7j, the velocity at the edge of the boundary layer at the point of separation. The frequency with which the vortices are formed can be considered proportional to U\,/b! and the 'Wake Strouhal Number', 160 Sr, can be denned as Sr = — (5.1) , where / y = frequency of vortex shedding. By applying Bernoulli's equation to the flow at the separation point, just outside the shear layer, U = Uy/[1-C ],. b (5.2) PB where: U = free stream velocity; and Cps = base pressure. Putting k\ = — CPB], the base pressure parameter [158] U = 17*1. (5.3) b Now a new Strouhal number St is defined as . (5.4) where D= characteristic diameter of the body. Then - ^ (-) ' The value of h' was originally obtained by Roshko [159] using the notched hodograph theory. It gives the spacing of the shear layers when they become parallel. In general, it is difficult to apply the notched hodograph theory to complex bluff geometries. One alternative would be to solve, numerically, the Navier-Stokes equation, but this would require considerable effort and time. As a first approximation, one may replace the distance h' by h (the height of the bluff body projected normal to the flow). This is reasonable because the separation of the shear layers is directly related to the projected height of the bluff body. Also, h is a function of the angle of attack and can be evaluated quite readily. The Strouhal number based on the projected height (St ) p is defined as _f h_(St)h v ~~U~~ P 161 ~D~- ( 5 - 6 ) By introdcucing '5t ' in equation (5.6), one can define a universal nondimensional p parameter Su as St fh P For a given bluff body the quantities fy, v h and k± are all functions of the angle of attack (a). Hence for a given wind speed (U), The next step would be to assess sensitivity of the parameter Su to the angle of attack and geometry of the bluff body. 5.3.2 E x p e r i m e n t a l confirmation Figure 5-8(a) shows variation of the universal Strouhal number, Su, with a for two distinct bluff bodies (rectangular and square prisms) in absence of momentum injection. It is apparent from the figure that the points show scatter around a representative Su of around 0.155. This suggests relative insensitivity of Su to changes in the geometry and angle of attack. For a circular cylinder, a flat plate and a 90° wedge Roshko [159] obtains a mean value of about 0.163. Similarly, Bearman [155] gave a value of 0.181 for a blunt trailing-edge body with either a splitter plate or base bleed. It should be pointed out that their definition of the universal parameter (Su = fvh'/Ub) is somewhat different than the one used here. Figure 5-8(b) show effect of the momentum injection. A small increase in the mean value of Su is observed for the rectangular prism at Uc/U = 1. In case of the square prism, the universal Strouhal number is shifted to a higher value with a mean of 0.219. For a bluff body at a given angle of attack, St p and Su can be obtained from Figure 5-8. is known from equation (5.6) Now ki, and hence the base pressure coefficient, can be calculated from equation (5.7). Similarly if the wake pressure is known then the Strouhal number can be obtained at a given angle of attack. Thus correlation between St p and k± confirms the observation made by Roshko [159] that 162 0.6 -i 1 1 r- (a) I - 0.2 - 1 1- Re = 52 - 6 6 X 1 0 U /U = 0 c 0.4 -I • A 3 VERTICAL RECTANGLE SQUARE Su 6 •66 0 A A £ 0.0 0.6 • A 60 (b) 120 a 180 Re = 5 2 - 6 6 X 1 0 U /U = 4 c • A 0.4 3 VERTICAL RECTANGLE SQUARE Su A 0.2 . A ga A A • D o.o Figure 5-8 _l A • n I • • I A . A • A i I _ 60 a -J i 120 i i_ 180 Effect of the angle of attack and bluff body geometries on the universal Strouhal number Su: (a) no momentum injection; (b) with momentum injection. 163 the wakes of various bluff bodies are similar except the changes in the geometric and velocity scales. 5.4 Flow Visualization Figure 5-9 shows the effect of momentum injection on the wake of the rectangu- lar prism (Re = 16,200). For Uc/U = 0, the wake is wider than the projected height (h) of the rectangular prism. In fact, measurements from the photographs show that the maximum wake width is.around 2h, and occurs, approximately, the distance ' / J ' downstream from the base of the body. The region immediately downstream of the body is also the scene of intense periodic pressure and velocity fluctuations caused by the alternately shed vortices. The size of the shed vortex (Figure 5-9a) is nearly comparable to the size of the body. The shed vortices undergo rapid dissipation due to the viscous action. As explained in Section 4.5, the overall wake displays considerable unsteadiness. Introduction of the momentum injection (Uc/U > 0) causes the wake to become progressively narrower until it becomes the same as the vertical face of the rectangular prism. Any further increase in the momentum injection beyond Uc/U = 4 does not have significant effect on the wake. The extent of the wake along the flow direction is also progressively reduced however, even at Uc/U = 4, the wake does not close completely. For sufficiently high rates of the momentum injection (Uc/U > 2) vortex shedding is completely eliminated, the wake becomes quite small and the overall fluid motion approaches the potential flow character. The reduction in the wake width implies drag reduction, which was confirmed by the wind tunnel tests results presented earlier. The results for the square prism at a = 15° and 90° are shown in Figures 5-10 and 5-11, respectively (Re = 23,400). As seen in Figure 5-10(a), in absence of any momentum injection the top and bottom faces are immersed in a separated flow field. A large vortex can be seen in the process of formation: As Uc/U is increased, the separation at the top disappears completely and the vortex formation is suppressed. 164 There is an associated rise in the wake pressure which causes the drag coefficient to be lowered as seen before. Note, the stagnation point on the front face has moved away from the mid-point towards the bottom corner. In Figures 5-11 (a = 90°), both the cylinders are rotating in the same direction (i.e. clockwise) to achieve a maximum benefit. For Uc/U = 4, the wake is nearly eliminated. Figure 5-ll(c) shows that the fluid at the top is accelerated, where as at the bottom the speed remains essentially unchanged. A region of low pressure is created on the top of the square prism which explains a large increase in the lift coefficient (Figure 5-5b). The flow visualization study brings to attention the fact that to obtain a further control of the flow separation, it is necessary to provide rotating elements at all four corners of the body. Figure 5-12 attempts to show, schematically, anticipated flow patterns for a rectangular prism as affected by the momentum injection and angle of attack. At a — 0 it is necessary rotate at least cylinders A and B , but the wake does not close completely. To achieve the wake closure, it may be necessary to inject a small amount of momentum (Uc/U = 1 or 2) through rotation of cylinders C and D. For 0 < a < 90°, cylinders B, C and D are the most important ones to close the wake. Similarly at a = 90°, cylinders A , B would have to rotate clockwise and the cylinders C , D in the counter-clockwise direction to eliminate the wake completely. Such a scheme can be easily implemented using a microprocessor controlled switching actuator in conjunction with a wind direction sensor. 5.5 Summary The results from the wind tunnel tests and the flow visualization study confirm that application of the momentum injection to the rectangular and the square prisms lead to significant drag reduction. In addition, the vortex shedding is completely suppressed, thereby eliminating the problem of vortex resonance. The width and the length of the wake is reduced considerably, thereby reducing the interaction between the wake and the afterbody, which is primarily responsible for the galloping instability. 165 Figure 5-9 Effect of momentum injection on the characteristic features of the wake and the shed vortex system associated with the rectangular prism (AR = 0.5) as observed during the flow visualization study. 166 167 Figure 5-11 Flow visualization pictures for a square prism at a — 90°. 168 a = 0 U /U c a = 45° a = B C A D = 0 Uc/U = 0 = 0 90° Figure 5-12 =4 Uc/U A schematic diagram showing possible strategies of momentum injection for a rectangular prism undergoing a change in the angle of attack. 169 A universal Strouhal number (Su), based on the St and the base pressure parameter k, is proposed. The results firmly establish effectiveness of the momentum injection to achieve boundary-layer control of rectangular prisms. 170 6. B L U F F B O D Y D Y N A M I C S I N P R E S E N C E O F M O M E N T U M I N J E C T I O N 6.1 Preliminary Remarks Fluid induced vibrations of bluff structures encountered in practice may be classified into three distinct categories: uncoupled vortex resonance and galloping instabilities, and coupled flutter vibrations. At times, such oscillations could have catastrophic consequences. One of the celebrated and disastrous case on record [160,161] is that of the Tacoma suspension bridge built across the Narrows at Puget Sound in the State of Washington, U.S.A. It had a central span of about 853 meters and a width of 12 meters. Based on several studies [162-164], it is believed that the bridge was destroyed by a combination of more than one type of vibratory phenomena. At the outset, the bridge was found to be prone to vertical oscillations (amplitude « 1.27m, frequency ~ 0.62 Hz). In less than five months after opening to the public (on July 1, 1940), the bridge collapsed on 7 November, 1940, in a strong gale of about 67 km/hr. Few hours before failure, torsional oscillations set-in, which eventually led to the collapse. The edges of the bridge had a maximum vertical displacement of about 8.5m with a ± 4 5 ° twist. The frequency of torsional vibrations was about 0.23 Hz. A rough estimate of the power supplied by the wind is about 4000 kW! A steady wind blowing around a structure can pump energy continuously into the system. Such a fluid loading is nonconservative due to the presence of an energy source. Furthermore, the system dissipates energy through structural and external damping. For a nonconservative system the fluid force can not be derived from a potential function and the Hamiltonian is not a constant of the motion [165]. In the equations of motion describing a structure undergoing flow-induced vibrations, the energy source is normally modelled as a nonlinear function of the structural velocity in the form of negative damping. Such self-excited systems begin to vibrate spontaneously when subjected to a small disturbance, the amplitude increasing, until 171 a balance is struck between the rate at which the energy is introduced into system and the rate at which it is dissipated. In self-excited vibrations the force sustaining the motion exists because of the motion itself. The governing equations of motion describing flow-induced vibrations are generally nonlinear and their solutions display a rich variety of physical phenomena generally not found in a linear system. Examples include frequency entrainment or lock-in during vortex resonance, hysterisis and jump in amplitude for galloping instability, as well as relaxation type oscillations due to velocity dependent friction, and others. The investigation reported in this chapter is mainly experimental. Accompanying discussion attempts to clarify the physical mechanism of the process and emphasize the challenge involved in modelling and controlling flow-induced vibrations in practice. The phenomena of vortex-resonance is briefly touched upon to begin with. Next, the effect of momentum injection in suppressing vortex-resonance is discussed based on the vortex shedding frequency and vibration amplitude measurements. The basic equation of motion governing the galloping instability in plunging is outlined and a simple criterion useful in identifying regions of instability is explained. The criterion is applied using the wind tunnel static test data to predict the dynamic behaviour of bluff bodies (D-section, square and rectangular prisms). Finally, results are presented to demonstrate the effectiveness of momentum injection in suppressing wind induced instabilities, both vortex resonance and galloping. Effect of asymmetric rotation is also touched upon. The chapter concludes with a short summary. 6.2 6.2.1 Control of Vortex Resonance Type Vibrations Mathematical description In order to better appreciate how momentum injection affects the vortex- excitation response of a structure, a short description of the model describing vortex172 resonance would be appropriate. The vibrations are in the plane perpendicular to the incident flow under the action of oscillatory pressure loading created by the alternately shed vortices. Ideally, the governing equations for this problem are the dynamic equation of motion for the bluff body and the time-dependent Navier-Stokes equations for the fluid flow. Flow separation and vortex shedding would be the natural outcome of-their solution, which would also provide the forcing function acting on the bluff body. This could be accomplished through an appropriate numerical scheme implemented on a computer. It, of course, would entail enormous amount of time and effort. As a closed-form solution of the Navier-Stokes equations is not available, considerable efforts have been expended in developing simpler models based on physical considerations. Among others, the wake-oscillator model by Hartlen and Currie [166] and the Iwan-Blevins model [167,168] are quite well known. Blevins- model in brief. Consider the Iwan- It utilizes a van der Pol type equation [169-171] with a 'hidden' flow variable to simulate the effect of the vortex shedding (Figure 6-la). The governing equations are: where: a" = pD Oi/(m 2 + a pD ), 2 3 i = 3,4; u> = v/fc7W[l + (a D )/m\; ( = 2 n f 3P a[ = Q i / ( o o + 0,3), i = 1,2, 3, 4; ( t = [((y/k/m)/u n + ( ]/(l f + a pD /m); 2 3 aipDU/(2mw ); n ( = system damping ratio; (t — total effective damping ratio; (f = viscous fluid damping ratio; k = stiffness; m = structure mass per unit length; p = density of fluid; D = body characteristic dimension; LJ V = circular frequency of vortex shedding; u) = n 173 system circular natural frequency; U = freestream velocity; u = translational velocity t of the vortex street; i = 'hidden' fluid variable (dz/dt) which can be considered as the weighted average of the transverse component of the flow velocity within the control volume; y = coordinate describing transverse displacement, t = time; a-i = dimensionless constants; K ' = K / (ao + as), and K = proportionality constant. Equation (6.1) describes an autonomous nonlinear fluid oscillator coupled with equation (6.2), which is a linear differential equation describing motion of the bluff body. The first term on the right-hand-side of equation (6.1) represents negative damping, meaning extraction of energy from the freestream fluid and delivering it to the structure in motion, which causes build-up of the vibrations. The second term on the right is cubic (in velocity) which eventually limits the amplitude of the vibrations. The terms on the left-hand-side of equation (6.1) represent the feedback between the wake and the boundary-layer on the bluff body. The terms on the right-hand side of equation (6.2) represent the excitation which couple the two equations. The frequency of vortex shedding (fv) varies proportionally as the wind speed is increased from zero. At some critical value of the wind speed, it coincides with the natural frequency (/ ) of the bluff structure resulting in resonance (Figure 6-la). n This can also be explained by employing the concepts of bifurcation and catastrophe as applicable to nonlinear dynamical systems in the phase space [172-179]. The change in behaviour as described by transition from stability to instability and vice versa due to variation in a given parameter is known as bifurcation. Mathematically, "Catastrophe" describes a point in some parameter space at which a dynamical system loses structural stability, so that its solutions change their topological type [174]. Nonelementary catastrophes have a more dynamic nature, one example being the Hopf bifurcation [179] to a limit cycle. For a real life bluff structure, generally all parameters (e.g. stiffness, mass, internal damping etc.) remain fixed, except the fluid velocity. If the fluid velocity is increased form zero, the solution in the phase space (y, y, t) in 174 H I D D E N FLUID V A R I A B L E (z) J ] CONTROL V O L U M E ! SEPARATING BOUNDARY-LAYER R O L L I N G U P INTO D I S C R E T E V O R T I C E S FLUID DENSITY VERTICAL VELOCITY OF T H E F L U I D IN W A K E CONSTANT OF PROPORTIONALITY DETERMINED FROM THE EXPERIMENT V O R T E X SHEDDING F R E Q U E N C Y oTURE N A T U R A L F R E Q U E N C Y = in at resonance Figure 6-1 Schematic diagram showing a bluff structure experiencing vortex excitation: (a) physical model; (b) variation of vortex shedding frequency and resonance; (c) phase space representation of the dynamical model. 175 the form of a singularity at the origin (i.e. a stable focus representing the equilibrium position) can transform itself into a growing limit cycle (enclosing an unstable focus) at some critical velocity. A n unstable focus under the action of small perturbations (e.g. a wind gust) gives rise to a limit cycle whose amplitude is a function of the fluid velocity. Figure 6-l(b) shows schematically the topological behaviour of a self-excited system. A stable focus/limit cycle acts as an attractor since all perturbations away from it are damped out. Similarly, an unstable focus/limit cycle acts as a repeller as all perturbations away from it tend to grow. The characteristic equation derived from the equation of motion has complex roots (eigenvalues) with positive real part in the region of instability [173]. This situation is the exact opposite of a damped linear oscillator (which has negative real roots). The model also successfully explains the entrainment of the vortex shedding frequency by the natural frequency of the system. It is apparent that the vortex resonance vibrations arise out of the nonlinear interactions between the vortex wake and the bluff body. Unfortunately, the fluid dynamics during vortex shedding, especially in the near-wake region, is rather complex and no exact solution to the fluid-structure interaction problem has yet been found. Any attempt to control the vortex resonance could greatly benefit by preventing the very process of vortex formation. This is where momentum injection could prove useful. The control of boundary-layer separation and resulting suppression of alternate vortex shedding would mathematically imply vanishing of the hidden fluid variable dz/dt and hence an elimination of the forcing terms in equation (6.2). 6.2.2 Suppression of the vortex shedding Figure 6-2 presents typical set of results showing the effect of the momentum injection on the vortex shedding for a D-section. Similar results were obtained for all the bluff bodies (D-section, square and rectangular prisms) at various angles of attack, without and with momentum injection. As described in Section 2.2.3, the vortex shedding frequency and the signal amplitude were obtained by performing 176 5 1 4.5 1 1 r— r- 1 1 — 1 1 • 4 - S/R\ 3.5 - U + S/R . ., 3 Normalized / A m p l i t u d e 2.5 2 1.5 1I 1 0 0 10 12 5 ~i 1 14 16 18 20 f, Hz ( ) a — r ~ -I r - 4.5 4 3.5 Uc/U = 1 + + + Uc/U = 2 u 3 Normalized A m p l i t u d e 2.5 12 14 20 f, Hz Figure 6-2 Power spectral representation of the wake pressure showing reduction in the strength of the shed vortex in presence of the momentum injection: (a) U /U c = 0; (b) U /U c 177 = 1, 2. the Fast Fourier Transform ( F F T ) analysis on the time dependent pressure signals, obtained through the data acquisition system. The trailing vortices cause pressure fluctuations as they are convected downstream in the wake. The dominant frequency , in the pressure signal is the vortex shedding frequency. The amplitude of the signal is a measure of the strength of the shed vortex. In absence of momentum injection (Uc/U = 0), Figure 6-2(a) clearly shows the vortex shedding at f v ~ 7.0 Hz, with a large signal amplitude. For a small amount of momentum injection (Uc/U = 1), the vortex shedding becomes much weaker as indicated by the drop in the signal strength (Figure 6-2b). Also, the frequency of vortex shedding has increased to about f v 8.0 Hz. The corresponding Strouhal number (St) increases from 0.208 (Uc/U to 0.237 (Uc/U = = 0) = 1) suggesting a reduction in the effective bluffness. At Uc/U = 2, the vortex shedding is virtually eliminated since no pressure fluctuations could be detected. This fact is also confirmed by the flow visualization photographs and a video movie. The complete suppression of the vortex shedding in presence of momentum injection can be interpreted in terms of the interaction between the wake and the bluff body. The introduction of the momentum injection reduces the extent of the wake and delays the separation of the boundary-layer. Mathematically this amounts to setting z and dz/dt equal to zero in equations (6.1) and (6.2). The fluid in the wake no longer rolls-up into vortices. The momentum injection prevents transfer of the convective energy in the wake to the structure thus preventing growth of structural vibrations.' 6.3 Control of Galloping Type of Instability All bluff bodies with noncircular cross-sections are susceptible to self-excited oscillations referred to as galloping. Pure galloping is a single degree of freedom phenomena in which the body experiences plunging vibrations in a plane perpendicular to the direction of the flow. Because of the inherent difficulties in modelling unsteady 178 aerodynamics around a vibrating structure, it becomes necessary to invoke the concept of quasi-steady flow, i.e. the aerodynamic force on a bluff body in motion, having a certain instantaneous orientation, can be approximated by a force acting on a fixed body having the same angle of attack. Any changes in the aerodynamic force arising out of the unsteadiness in the flow are neglected. This assumption has proved generally valid as demonstrated by Parkinson and Brooks [180], Novak [84], and others. The quasi-steady analysis works well at higher fluid velocities in the range of U / a i D n > 3 where U = freestream velocity; u} = circular natural frequency of the. structure; n and D = characteristic dimension of the structure. For U / w D < 3, quasi-steady n analysis is questionable due to the likelihood of vortex-resonance. The problem of galloping instability has received considerable attention in the past with important 7 contributions from Parkinson [19], Novak [84,181], Richardson [182], Blevins [71], Lawson [183], and others. In the next section mathematical description of galloping instability is briefly reviewed to emphasize key parameters governing the vibrations [19,71, 84, 157, 181]. Further details pertaining to nonlinear dynamics of self-excited vibration can be found in references [172-179]. 6.3.1 Criterion for determining the galloping instability Consider a tall bluff structure with square cross-section undergoing galloping type of oscillations (Figure 6-3). The stiffness (k) and the internal damping (ri) of the structure are assumed to be linear. By convention, the structural motion is considered positive (+y) downwards. The resultant fluid velocity (-y) acts upwards. Please note that this transverse motion in reality lies in a horizontal plane. The instantaneous angle of attack (a) is given by a = tan- ^, x (6.3) where y — dy/dt. The lateral aerodynamic force ( - F ) is in the upward direction and opposes the y motion. The absolute value of Fj, is not of concern here, but its variation with respect 179 my + r y + ky = ± Figure 6-3 -F {a) Y A bluff structure experiencing galloping and the criterion for instability. 180 to. a determines the stability of the structure. If F increases with a, then a sudden y transverse motion (+y) of the structure (possibly caused by a wind gust) would increase a and hence decrease the aerodynamic force (~F ) opposing the motion. y This would leave the structure more vulnerable to further displacement. situation is clearly unstable. structure remains stable. For most aerodynamic shapes dF /doc y Such a < 0, and the But, for some rectangular geometries (0.75 < A R < 3) there is a likelihood of dF /da y to be negative in a certain range of a. The D-section is neutrally stable in a smooth flow (dF /dct y — 0) for 0 < a < 3 0 ° , but becomes unstable at a = 11° in a turbulent flow [184]. The criterion for instability is dF , v . The aerodynamic forces are given as: F =C ^phlV ; (6.5) 2 D D F x = C ^phlV ; (6.6) 2 L L 1 F =C lphlU ; (6.7) 2 y y 2 with cosa — — . (6-8) The lateral aerodynamic force can be expressed as F Y = -(F cosa + F sina). L (6.9) D Substituting the expressions for FD and FL in equation (6.9) and simplifying gives C Y = —(CL + C£>tana)seca. (6.10) Equation (6.11) can be used to evaluate the lateral force coefficient (C ) y from the drag and the lift coefficients obtained through static wind tunnel tests. At a = 0, the 181 instability criterion becomes dC L da + C D . da > 0, (6.11) i.e. dC L + C D da < 0. (6.12) Equation (6.12) represents the Den Hartog criterion [154] for instability. The equation of motion governing the vibrations is my + r y + ky = F . x (6.13) y The lateral force coefficient (C ) is a function of a, and therefore of tana — y/U. y C y can be approximated by a polynomial in y/U over the pertinent range of a. It is an odd function by the reason of symmetry, and requires at least a seventh degree polynomial for a close fit [18], *<•> = * ( £ ) - * ( £ ) + * $ - * ( (6.14) The equation of motion now becomes my + r y + ky = x -phlU' *(£)-«(f)'+»(£)"-*(£ (6.15) The P{ (i—1,3,5,7) are the coefficients of the expansion and can be obtained from the curve of C against a using the least square technique. Defining nondimensional y variables as: Y = y/h; r = ut; oo = k/m; V = U/(uh); n = ph l/(2m); (3 = r/(2mw); 2 2 i.e. Y + 20Y + Y = i.e. 182 nV C , 2 Y (6.16) Y + Y = nPi where Y = dY/dr. This is a nonlinear ordinary differential equation. As a first approximation, neglecting the third and higher degree terns in Y gives the total damping as 2(3 nPAV Here the term (3 represents the internal damping of the structure, which remains constant. As long as the fluid velocity V < 2(3/(nP\), there is a net positive damping, i.e. dissipation of energy. At some critical velocity, Vc — 2(3/(nP\), the total damping is zero and beyond that the damping becomes negative leading to the galloping instability. In general, for a given structure, the system parameters (e.g. stiffness, mass, internal damping, etc.) remain fixed, and the only control variable is the fluid velocity. The fluid in motion can pump energy into the system causing build-up of large, low frequency vibrations. The higher order terms in equation (6.17) eventually limit the amplitude of motion. The resultant oscillatory motion has a character of a limit cycle in the phase plane (Y, Y). As V increases from zero and crosses Vc, a stable focus (attractor) bifurcates into an unstable focus (repeller) with a growing limit cycle. In fact, as demonstrated by Parkinson and Smith [18], the characteristic equation representing the galloping motion is a cubic and its roots describe a family of stable and unstable limit cycles in the phase plane. In equation (6.17), the highest derivative being of the second order, it is not possible to observe a class of chaotic bifurcations known as strange attractors. To put it differently, the system loses structural stability (in some parameter space) so that the solution changes its topological character. With momentum injection the system has one additional control parameter, Uc/U. Section 6.4 discusses how an originally unstable.bluff body (V > Vc) can be stabilized by a systematic variation of 183 Uc/U. 6.3.2 Instability prediction based on static aerodynamic tests The plots of C vs. a for the D-section as well as the rectangular and square y prisms are presented in Figure 6-4. As explained in Section 6.3.1, these results were obtained from the lift and drag data obtained through the wind tunnel tests. The test Reynolds number ranged form 52 x 10 - 66 x 10 . The D-section is neutrally 3 stable (Figure 6-4a) as 3 dC /dct = 0 in the range 0 < a < 2 0 ° . The presence of y rotating cylinders has the effect of rounding the corners of the D-section. Parkinson and Brooks [180] have reported stability study of a D-section with sharp corners. The results for the two cases (model with sharp corners and round corners) were found to be identical upto a = 20°. As obvious from Figure 6-4, the stability of the D-section significantly improves (dC /dct < 0) in presence of the momentum injection (Uc/U y = 2, 4). In the range 20° < a < 3 0 ° , there is a mild tendency towards instability. The region of concern is 30° < a < 4 0 ° , where dC /da = +2.1 at U /U = 0, and the y c D-section would exhibit strong instability. The unstable behaviour of the D-section at around 40° has been reported by Parkinson and Brooks [180]. However, in presence of the momentum injection (Uc/U = 4), the stability is restored completely. Overall, there is a marked improvement in the stability of the D-section with momentum injection. Results for the rectangular prism (Figure 6-4b) show similar trends in the range a < 2 0 ° , with the system maintaining neutral stability. This is expected due to similarity in their shapes. This is in accord with the results obtained by Parkinson and Brooks [180] for a rectangular prism with sharp edges. In presence of momentum injection (Uc/U = 4), there is a marked improvement in the stability as dC /da fa y -3.7. Even in absence of momentum injection (Uc/U stable in the range 20° < injection, dC /da y = 0), the rectangular prism is a < 4 5 ° (dC /da fa -2.7). In presence of the momentum y becomes more negative (fa -3.0) suggesting further increase in stability. 184 Figure 6-4 (a) Influence of the momentum injection on the lateral force coefficient for the D-section. The criterion for galloping instability is dC /dct y > 0. Note a marked improvement in the stability with the momentum injection (Uc/U = 4). 185 00 Figure 6-4(b) The rectangular prism is neutrally stable (dC /da 20° in absence of momentum injection (Uc/U is significantly improved (dC /da y (Uc/U = 2, 4). 186 y — 0) for a < = 0). The stability < 0) with the cylinder rotation Figure 6-4(c) The square prism is unstable (dC /da > 0) for a < 5° in absence of momentum injection (Uc/U y — 0). With the cylinder rotation (Uc/U = 2, 4) the stability is restored completely (dC /da < 0). y 187 ) The results for the square prism are shown in Figure 6-4(c). For 0< a < 5 ° , the system shows strong instability injection. For 5° < (dC /da « +2.6) in absence of the momentum y a < 4 5 ° , dC /da < 0, and the square prism is stable. In presence of the momentum injection y (Uc/U = 4), dC /da ; becomes -4.6 y (a < 5 ° ) and the stability is restored. It may be pointed out that C vs. a curve for the square prism y with rounded corners is somewhat different from that obtained by Parkinson and Brooks [180], though the general trend is similar. In the present case, the region of instability extends from 0 to 5° for the prism with round corners, whereas it is 0 to 15° for the prism with sharp corners. Thus, interestingly, mere presence of cylinders at the corners of the upstream face (Uc/U = 0) significantly reduces the range of galloping instability of a square prism. 6.4 Confirmation Through Dynamic Experiments Though the wind tunnel data provide sufficient information to predict possible susceptibility of a structure to galloping instability, it was considered appropriate to conform the predictions through experiments. There are several reasons for this. In actual practice, a structure (e.g. a tall building with complicated geometry) is exposed to a a fluid stream that is normally turbulent. The structure generally has nonlinear elastic and damping characteristics. It is also known that vortex resonance could set in for any structure if the vortex shedding and structural frequencies are close enough. Thus it is necessary to confirm the effect of momentum injection in controlling both the vortex resonance and the galloping type of instabilities through actual dynamic experiments. Such experiments were performed on elastically mounted bluff body models (D-section, rectangular and square prisms), free to oscillate transverse to the flow in the wind tunnel. The description of the experimental set-up was explained in Section 2.2.4. The determination of the system parameters (total damping, natural frequency, mass, stiffness etc.) and the calibration of the displacement transducer 188 are given in Appendix D. The damped natural frequency of the system was found to be very close to the undamped frequency. The experiments were carried out at zero angle of attack. 6.4.1 - - Dynamical response of the square prism Figures 6-5 to 6-7 show the amplitude response of the square prism at a = 0 as affected by the momentum injection and damping. In absence of the momentum injection, as the wind speed (V) is increased form zero, the square prism experiences vortex- resonance (Figure 6-5) in the range of 0.9 < V < 1.25. This finite band of velocities indicate the entrainment of the vortex shedding frequency by the system natural frequency. For V > 1.25, the square prism becomes stable as the vortex shedding frequency is now higher than the system natural frequency. With a further increase in the wind velocity, there is an onset of galloping instability at Vc ~ 2.0. Beyond this the amplitude of vibration grows quickly to its limiting value (Y — 0.26). Any further increase in the wind velocity does not change the overall dynamical response. In contrast to this, with the momentum injection corresponding to Uc/U = 1, the square prism does not exhibit either vortex resonance or galloping. The stability of the prism is restored completely. This is in accord with the prediction based on the simple stability criterion mentioned earlier. It is remarkable that despite completely different mechanisms responsible for vortex resonance and galloping type of instabilities, the momentum injection is effective in suppressing both of them. With an increase in the damping coefficient (/3) from 0.2693% to 0.9311% (Figure 6-6), the maximum amplitude of the vibrations is considerably reduced, both for the vortex resonance and galloping. As before, a small amount of momentum injection again restores the stability. In the previous two cases, both the cylinders were made to rotate so as to inject momentum downstream. One can disrupt formation of the vortex street by rendering the flow asymmetric. This can be achieved by rotating only one of the 189 III III I III u "i 0.4 1 1 1 1 1 1 I 11111 U 1 1 1 11II r—| 1 1 1 r 1 1 1 1 r (3=0.002693 a = 0 -•— U /U = 0 c U /U= 1 A c S Y S T E M LIMITING A M P L I T U D E 0.3 Y a • a •- VORTEX RESONANCE •A Figure 6-5 1 £,—i L-A-J 1—A< L. 10 Dynamic response of the square prism without and with symmetric momentum injection. Note the presence of both vortex-resonance and galloping in absence of the momentum injection (Uc/U = 0). The square prism becomes completely stable in presence of the momentum injection as small as Uc/U 190 = 1. cylinders. Results for this case are presented in Figure 6-7. Note, the vortex excited vibrations disappear completely for Uc/U =1 and the galloping is delayed to V = 6.4. With a small further increase in the injected momentum (Uc/U =1.5), the galloping oscillation are suppressed for V as large as 10. Recently, Kubo et al. [185] have also carried out similar experiments on a two- dimensional square prism. Their results also demonstrate that symmetric momentum injection is effective in suppressing galloping response of the square prism. But the vortex-resonance could not be eliminated completely in presence of symmetric momentum injection, though the peak amplitude showed significant reduction. To eliminate vortex-resonance completely, it was necessary to introduce asymmetric momentum injection (i.e. only one cylinder rotating). This difference in response may be attributed to the difference in the relative size of the rotating element with respect to the prism. In the present case, ratio of the rotating cylinder diameter to the side of the square prism is about 0.3, whereas for Kubo's model this ratio is 0.1 (i.e. the momentum injection elements are much smaller). This would mean that at a given rate of rotation, the cylinders in the present model can deliver larger amount of energy to the boundary-layer, and hence are more effective. The oscillatory behaviour of the square prism without and with momentum injection is represented schematically in the phase space in Figure 6-8. mentum injection, there are two control parameters, V and Uc/U. With mo- The oscillations are represented by trajectories in the phase space (Y, Y, V for a given Uc/U). Uc/U With = 0, the dynamic behaviour of the square prism can be observed qualitatively as the wind speed is increased from zero (Figure 6-8a). The system demonstrates bifurcation twice, once during the vortex-resonance and later during galloping. The bifurcation represents transformation of a stable focus (attractor) into an unstable focus (repeller) with the change in the wind velocity. A stable limit cycle (attractor) encompasses an unstable focus (repeller). For V > Vc (i.e. in the region of galloping), 191 F i g u r e 6-6 Effect of an increase in damping on the dynamic response of the square prism without and with symmetric momentum injection. Note, a significant reduction in the amplitude, particularly in galloping, in presence of the momentum injection. As before, the stability is restored for Uc/U = 1. 192 F i g u r e 6-7 Effect of asymmetric momentum injection on the dynamic response of the square prism at various wind speeds. Note, it becomes necessary to increase the asymmetric .momentum injection to at least Uc/U =1.5 to restore complete stability. 193 (a) <y = 1 Uc/U A v U /U c ATTRACTOR COMPLETE ^STABILITY BIL Y = 0 V ATTRACTING LIMIT C Y C L E V ATTRACTOR REPELLER ATTRACTOR VORTEX RESONANCE GALLOPING (b) Y . U /U c = = 1 Uc/U 0 COMPLETE STABILITY Y o Figure 6-8 Schematic diagrams showing dynamics of bluff bodies in the extended parameter space (Y, Y, V, Uc/U): (a) an unstable bluff structure displays both the vortex-resonance and galloping in absence of the momentum injection (Uc/U = 0). Stability is completely restored with small amount of momentum injection (Uc/U = 1); (b) increase in the system damping affects the critical wind speed and system amplitude. 194 with the momentum injection at Uc/U = 1, the galloping is completely eliminated. This is shown by an inwardly spiralling (stable) trajectory. With an increase in damping (Figure 6-8b), limit cycles representing vortex-resonance and galloping are separated further, and galloping occurs at a higher critical wind velocity. Also, the peak response is reduced. Figure 6-9 is based on experimental observations and clearly emphasizes importance of maintaining the correct direction of rotation of the momentum injecting cylinders. For V > Vc, the square prism would gallop in absence of the momentum injection {Uc/U = 0). However, at Uc/U > 0, the vibrations become progressively smaller and vanish eventually when Uc/U > 1.5 (Figure 6-9a). Reversing the di- rection of momentum injection of both the cylinders, galloping can be induced even if V < Vc! Such behaviour was indeed confirmed during the dynamic experiments when the direction of cylinder rotation was reversed. 6.4.2 Dynamic response of the D-section and the rectangular prism Dynamical responses for the D-section and the rectangular prism are shown in the Figures 6-10 and 6-11, respectively. Both the bluff bodies display vortex-resonance in absence of the momentum injection. For the D-section it is necessary to increase the momentum injection to Uc/U = 1.5 to completely suppress the motion. On the other hand, the rectangular prism becomes stable at Uc/U = 1 as the top and bottom faces promote reattachment. Increasing the wind speed to a value as high as 10, does not make either the D-section or the rectangular prism susceptible to galloping. This is in agreement with the theoretical predictions based on the static wind tunnel test data. Both the D-section and the rectangular prism are neutrally stable to galloping (dCy/da = 0) for a — 20°, in absence of the momentum injection. With an appropriate level of momentum injection the stability is maintained as suggested by the results in Figures 6-10 and 6-11. 195 v > v c ATTRACTING, LIMIT C Y C L E (b) •Figure 6-9 Effect of direction of the momentum injection on the galloping instability: (a) for V > Vc, an inherently unstable system can be stabilized by momentum injection in the right directionj(b) galloping instability at subcritical wind velocity (V < Vc) under adverse direction of the momentum injection (Uc/U 196 < 0). F i g u r e 6-10 Dynamic response of the D-section as affected by the momentum injection. The D-section experiences vortex-resonance, but is stable in galloping. 197 F i g u r e 6-11 Effect of momentum injection on the galloping instability of a rectangular prism (AR=0.5) at a = 0. Note, the prism experiences vortex resonance, but is stable in galloping. 198 6.5 Summary The results pertaining to the dynamic response of three bluff bodies under consideration confirm effectiveness of the momentum injection in suppressing both the vortex-resonance and galloping type of instabilities. Even a small amount of the momentum injection (Uc/U « 1-1.5) is sufficient to restore stability. The qualitative phase space representation provide insight into the complex character of the flowinduced vibrations and their dependence on various system parameters. 199 7. C L O S I N G C O M M E N T S 7.1 Summary of Conclusions The thesis presents results of a rather fundamental study aimed at the Moving Surface Boundary-layer Control (MSBC) as applied to a family of two-dimensional bluff bodies represented by a symmetrical airfoil at high angles of attack, a D-section, and a set of rectangular prisms. The comprehensive character of the investigation is revealed by the four distinct approaches involved: static model wind tunnel tests, numerical simulations, flow visualization, and dynamical response of the system in vortex resonance and galloping, without as well as with the M S B C . The study shows original contributions in every phase. It reveals rather remarkable effectiveness of the M S B C in increasing lift, reducing drag, and suppression of both vortex resonance and galloping. Successful elimination of wind-induced instabilities establishes the M S B C as a vibration suppression method based on the control of excitation, and hence energy input to the system, rather than dissipation. More important conclusions based on the study may be summarized as follows: (a) One of the major parameters affecting the M S B C is the ratio Uc/U. most applications the maximum value of Uc/U In required to obtain desirable effect of the M S B C is limited to 4. In fact, the wind induced instabilities can be completely suppressed for Uc/U as small as 1.5, which agrees with the predictions based on the quasi-steady analysis. (b) The direction of the momentum injection depends on the location of the upstream stagnation point and requires reversal at critical orientations to promote reattachment of the separated boundary-layer. (c) <- Surface roughness of the rotating elements has significant effect on the efficiency of the momentum injection. performed the best. 200 In the present study, splined cylinders Effectiveness of the M S B C suffers when the rotating element is located in the separated shear layer. Measurement of the vortex shedding frequency and signal amplitude shows that the shed vortices become weaker and eventually disappear as the momentum injection is increased from Uc/U = 0 to 2. This is accompanied by reduction in the size of the wake, rise in wake pressure (average 23% over 0 < a < 90° at Uc/U =4), delay of flow separation, and reattachment of the separated shear layer. In presence of the M S B C (splined cylinder), the peak lift coefficient for a symmetric Joukowsky airfoil increased by 210% with the stall delayed from 10° to 50°! Furthermore, the stall sets in rather gradually, i.e. dCz,/da is relatively small, which is a desirable feature. The presence of suction peaks at the rotating elements contribute significantly to the drag reduction of bluff bodies. The peak drag reduction was found to be as high as 75% with the average over 0 < a < 90° of around 36%. Results of dynamic studies also reveal effectiveness of asymmetric momentum injection in eliminating wind-induced instabilities, which is based on the disruption of organized vortex shedding. Previous as well as present estimates suggest that the power consumption in implementing the M S B C is quite small. This makes the M S B C essentially a semi-passive flow control procedure. The analysis based on the universal Strouhal number (Su) confirms the observation that the wakes of different bluff bodies are similar except the changes in the geometric and velocity scales. The simple modified panel method successfully predicts correct trends even for such complex flow field associated with multielement geometries in presence 201 of the momentum injection. This suggests that judiciously applied numerical approach, capturing physics of the system under consideration, can provide results of sufficient accuracy saving considerable computational time and effort. (1) The flow visualization study substantiates, rather spectacularly, effectiveness of the M S B C . (m) Such a comprehensive investigation of the bluff body aerodynamics and dynamics in presence of the momentum injection, based on extensive wind tunnel, numerical and flow visualization experiments, has never been reported before. It presents an exciting possibility of application to a wide variety of problems of contemporary interest. 7.2 Suggestions for Future Research There are several avenues for further exploration of the fluid dynamics and dynamics of bluff bodies in presence of the M S B C . Some of the more important aspects, which are likely to advance the field further, are briefly touched upon below. (i) The present study focussed on the M S B C as applied to two-dimensional flows. The next logical step would be to consider its application to three-dimensional objects which would be closer to the real life situation. This suggests performance evaluation of the M S B C in presence of end effects. (ii) Application of the M S B C to supertall buildings of the future (> 1,000 m), which are in the planning stage, should prove to be challenging both aerodynamically as well as structurally. Immersed in the atmospheric boundary-layer, they will be exposed to random gust-type loading as observed in the case of the Empire State Building. The reduced stiffness of such gigantic structures (< 1 Hz.) will only accentuate the problem. This may require hybrid strategy involving the M S B C for reducing force excitation as well as passive nutation dampers for dissipation of energy. It is also imperative to obtain quantitative 202 estimate of the optimum size of a rotating cylinder as compared to the overall dimensions of a tall building. The next generation of bridge-towers for long suspension bridges have shown them to be prone to wind induced instabilities during wind tunnel tests. Application of the M S B C (without and with nutation dampers) may prove to be one of the viable solutions. The present study focussed on application of the M S B C to straight, rectangular wings. There is an exciting possibility of its implementation on a delta wing and its control surfaces. Injection of momentum may stabilize the twin helical vortices, delay or even suppress the vortex burst phenomenon and improve the overall aerodynamic performance. When applied to control surfaces (aileron, elevator, rudder), it can lead to highly maneuverable airplanes. A . real wing is a complex system. Incorporation of a rotating cylinder on such a wing would require a careful analysis of its effect on the aerodynamics and dynamic stability of the aircraft. A numerical study of a rotating cylinder having time-dependent rate of rotation has been reported by Ou [186]. For an aircraft undergoing unsteady manoeuvres, such a time-dependent momentum injection (Uc(t)/U) may be- come necessary. Further information regarding this concept as applied, to an actual aircraft is necessary. Application of the concept to control wing tip vortices and promote their dissipation may prove to be interesting. This may ease the problem faced by airtraffic controllers, particularly at large busy airports, who are forced to increase the time interval between successive landings to minimize the effect of wing tip vortices. Some wind tunnel model test results aimed at reduction of the aerodynamic drag experienced by highway trucks, with the M S B C applied at the top lead- 203 ing edge of the trailers, have been reported in the literature [59]. However, systematic road tests are necessary to confirm the results. Even 1% reduction in the drag can literally translate into' millions of dollars of saving in the fuel cost. (viii) One aspect that is relevant in all the above mentioned studies concerns the effect of surface roughness of the rotating cylinder on the efficiency of the momentum injection process.. In the present study, the cylinder with splined surface performed the best. There is a scope to arrive at an optimum configuration of the splined surface geometry. (ix) Recent advances in digital image processing could be useful in analyzing the flow visualization data to obtain quantitative estimates of the vortex shedding process. This would provide further information as to the effect of the M S B C on the structure of the wake. (x) There is a scope for numerical simulation of such multielement systems using the complete Navier-Stokes equations in conjunction with appropriate turbulence model. It can provide a rather powerful and versatile tool, applicable to a large class of systems, when combined with the dynamics of the system. 204 REFERENCES Rosenhead, L . , "Vortex Systems in Wakes," Advances in Applied Mechanics, Vol. 3, 1966, Academic Press, New York, p. 189-198. Prandtl, L . , and Tietjens, 0. G . , Applied Hydro and Aeromechanics, Dover Publications, New York, 1934, p. 58. Thwaites, B., Incompressible Aerodynamics, Dover Publications, New York, 1960, p. 587. Kubo, Y . , Yasuda, H . , and Kotsubo, C . , Hirata, K . , "Active Control of Super-tall Structure Vibrations under Wind Action by a Boundary-layer Control Method," Journal of Wind Engineering and Industrial Aerodynamics, Vol. 50, 1995, pp. 361-372. Kubo, Y . , Hayashida, K . , Kotsubo, C , and Kato, K . , "On the Suppression of Aerodynamic Vibrations of a Square Prism in High Angles of Attack by boundary-Layer Control Method," First World Conference on Structural Control, Los Angeles, U.S.A., August 1994. Wille, R., "Karman Vortex Streets," Advances in Applied Mechanics, Academic Press, New York, 1960, Vol. 6, pp. 273-287. Marris, A . W . , "A Review on Vortex Streets, Periodic Wakes, and Induced Vibration Phenomena," Transactions of the ASME, Journal of Basic Engineering, Series D., Vol. 86, 1964, pp. 185-196. Morkovin, M . V . , "Flow Around a Circular Cylinder - A Kaleidoscope of Challenging Fluid Phenomena," ASME Symposium on Fully Separated Flows, Philadelphia, 1964, pp. 102-118. 205 Parkinson, G . V . , "Wind Induced Instability of Structures," Philosophical Transactions of the Royal Society of London, Series A, Vol. 129, 1971, pp. 115-135. Cermak, J . E . , "Application of Fluid Mechanics to Wind Engineering," Freeman Scholar Lecture, Transactions of the ASME, Journal of Fluids Engineering,, Vol. 97, No.l, 1975, pp. 9-38. Welt, F . , "A Parametric Study of Nutation Dampers," M.A.Sc. Thesis, University of British Columbia, Vancouver, Canada, 1983. Welt, F . , "A Study of Nutation Dampers with Application to Wind Induced Oscillations," Ph.D. Thesis, University of British Columbia, Vancouver, Canada, 1988. Modi, V . J . , Welt, F . , and Seto, M . L., "Control of Wind Induced Instabilities through Application of Nutation Dampers: A Brief Review," Engineering Structures: Auxiliary Damping Systems for Vibration Mitigation, Special Issue Guest Editor: J . D. Holms, in press. Tri.tton, D.. J . , Physical Fluid Dynamics, Clarendon Press, New York, 1988, 2 nd Edi- tion, pp. 18-29. Fluid Vortices, Editor: Green, S. I., Kluwer Academic Publishers, Boston, 1995, p. • 537. Bearman, P. W . , "Vortex Shedding from Oscillating Bluff Bodies," Annual Review of Fluid Mechanics, Vol. 16, 1984, pp. 195-222. , Bearman, P. W . , Gartshore, I. G . , Maull, D . J . , and Parkinson, G . V . , "Experiments on Flow-Induced Vibrations of a Square-Section Cylinder," Journal of Fluids and Structures, Vol. 1, 1987, pp. 19-34. Parkinson, G . V . , and Smith, J . D., "The Square Prism as an Aeroelastic Non-Linear 206 Oscillator," Quarterly Journal of Mechanics and Applied Mathematics, Vol. 17, Part 2, May 1964, pp. 225-239. Parkinson, G . V . , "Phenomenon and Modelling of Flow-Induced Vibrations of Bluff Bodies,"Progress in Aerospace Sciences, Vol. 26, 1989, pp. 169-224. Bokaian, A . R., and Geoola, F . , "Hydroelastic Instabilities of Square Cylinders," Journal of Sound and Vibrations, Vol. 92, No. 1, 1984, pp. 117-141. Laneville, A . , and Parkinson, G . V . , "Effects of Turbulence on Galloping of Bluff Cylinders," Proceedings of the 3rd International Conference on Wind Effects on Build-, ings and Structures, Editors: Japan Organizing Comittee, Tokyo, 1971, pp. 787-797. Gartshore, I. S., Khanna, J . , and Laccinole, S., "The Effectiveness of Vortex Spoilers on Circular Cylinders in Smooth and Turbulent Flow," Proceedings of the 5th International Conference on Wind Engineering, Fort Collins, Colorado, U S A , 1978. Engineering News - Record, McGraw-Hill's Construction Weekly, October 30, 1975, pp. 11. Welt, F . , and Modi, V . J . , "Vibration Damping Through Liquid Sloshing: Part 1 A Nonlinear Analysis," Transactions of ASME, Journal of Vibration and Acoustics, Vol. 114, N o . l , 1992, pp. 10-16. Welt, F . , and Modi, V . J . , "Vibration Damping Through Liquid Sloshing: Part 2 Experimental Results," Transactions of ASME, Journal of Vibration and Acoustics, Vol. 114, N o . l , 1992, pp. 17-23. Tamura, Y . , Kousaka, R., and Modi, V . J . , "Practical Application of Nutation Damper for Suppressing Wind-Induced Vibrations of Airport Towers," Journal of Wind Engineering and Industrial Aerodynamics, Vol. 41-44, 1992, pp. 1919-1930. 207 Goldstein, S., Modern Developments in Fluid Mechanics, Vol. I and II , Oxford University Press, Oxford, 1938. Lachmann, G . V . , Boundary-Layer and Flow Control, Its Principles and Application, Vol. I, Pergamon Press, New York, 1961, pp. 1-2. Rosenhead, L . Laminar Boundary Layers, Clarendon Press, Oxford, 1966, pp. I l l , 332-335, 547. Schlichting, H . , Boundary Layer Theory, McGraw-Hill Book Company, New York, 1968, pp. 362-390. Chang, P.K., Control of Flow Separation, Hemisphere Publishing Corporation, McGrawHill Book Company, Washington, 1970, pp. 154-177. Report of the Research Co-ordination Group on Boundary Layer Control to Suppress Separation, Associate Committee on Aerodynamics, National Research Council, September 1966. Thwaites, B., Incompressible Aerodynamics, Dover Publications, New York, 1960, pp 105-106, 215. Massey, B. S., Mechanics of Fluids, Van Nostrand Reinhold, Berkshire, U . K . , 5th Edition, 1983, p. 346. Swanson, W . M . , "The Magnus Effect: A summary of Investigation to Date," Transactions of the ASME, Journal of Basic Engineering, Vol. 83, September 1961, pp. 461-470. Iverson, J . D., "Correlation of Magnus Force Data for Slender Spinning Cylinders," AIAA 2nd Atmospheric Flight Mechanics Conference, Palo Alto, California, September 1972, Paper No. 72-966. 208 Betz, A . , "History of Boundary Layer Control in Germany," Boundary Layer and Flow Control, Editor: G . V . Lachmann, Pergamon Press, New York, Vol. I, pp. 1-20. Flettner, A . , "The Flettner Rotor Ship," Engineering, Vol. 19, January 1925, pp. 117-120. Favre, A . , "Contribution a l'Etude Experimentale des Mouvements Hydrodynamiques a Deux Dimensions," Thesis presented to the University of Paris, 1938. Alvarez-Calderon, A . , and Arnold, F. R., "A study of the Aerodynamic Characteristics of a High Lift Device Based on Rotating Cylinder Flap," Stanford University Technical Report RCF-1, 1961. Alvarez-Calderon, A . , "Rotating Cylinder Flaps for V / S T O L Aircraft," Aircraft Engineering, October 1964, pp. 304-309. Brooks, D. A . , "Effect of a Rotating Cylinder at the Leading and Trailing Edges of a Hydrofoil," U.S. Naval Ordanance Test Station, Department of the Navy, N A V W E P S Report 8042, April 1963. Steele, B. N . , and Harding, M . H . , "The Application of Rotating Cylinder to Ship Maneuvering," National Physical Laboratory, Ship Division, U . K . , Report No. 148, December 1970. Aviation Week and Space Technology, Vol. 95, No. 16, October 18, 1971, p. 19. Aviation Week and Space Technology, Vol. 95, No. 24, December 13, 1971, front cover and p. 7. Tennant, J . S., "The Theory of Moving Wall Boundary Layer Control and its Experimental Applications to Subsonic Diffusers," Ph.D. Dissertation, Clemson University, May 1971. 209 Tennant, J . S., "A Subsonic Diffuser with Moving Walls for Boundary-Layer Control," AIAA Journal, Vol. 11, No. 2, February 1973, pp. 240-242. Johnson, W . S., Tennant, J. S., and Stamps, R. E . , "Leading- Edge Rotating Cylinder for Boundary-Layer Control on Lifting Surfaces," Journal of Hydronautics, Vol. 9, No. 2, April 1975, pp. 76-78. Tennant, J . S., Johnson, W. S., and Krothapalli, A . , "Rotating Cylinder for Circulation Control on an Airfoil," Journal of Hydronautics, Vol. 10, No. 3, July 1976, pp. 102-105. Tennant, J . S., Johnson, W. S., and Keaton, D. D., "On the Calculation of Boundary Layers Along Rotating Cylinders," Journal of Hydronautics, Vol. 11, No. 2, April 1977, pp. 61-63. Tennant, J . S., Johnson, W. S., and Keaton, D. D., "Boundary Layer Flows from Fixed to Moving Surfaces including Gap Effects", Journal of Hydronautics, Vol. 12, No. 2, April 1978, pp. 81-84. Ericsson, L. E . , "Moving Wall Effects in Unsteady Flow," Journal of Aircraft, Vol. 25, No. 11, November 1988, pp. 977-990. Modi, V . J . , Sun, J . . L . C , Akutsu, T . , Lake, P., McMillan, K . , Swinton, P. G. and Mullins, D., "Moving Surface Boundary Layer Control for Aircraft Operation at High Incidence," Journal of Aircraft, AIAA, Vol. 18, November 1981, pp. 963-968. Modi, V . J . , Mokhtarian, F . , Yokomizo, T . , Ohta, G . , and Oinuma, T . , "Bound Vortex Boundary-Layer Control with Application to V / S T O L Airplanes," Proceedings of the IUTAM Symposium on Fundamental Aspects of Vortex Motion, Tokyo, Japan, 1987, pp. 196-199; also Vortex Motion, Editors: H . Hashimoto and T . Kambe, North- Holland Publishing Company, Amsterdam, 1987, pp. 225-230. 210 Mokhtarian, F . , and Modi, V . J . , "Fluid Dynamics of Airfoils with Moving Surface Boundary-Layer Control," Journal of Aircraft, Vol. 25, No. 2, February 1988, pp. 163-169. Modi, V . J . , Ying, B., and Yokomizo, T . , "Boundary-Layer Control of Bluff Bodies Through Momentum Injection," SAE Transactions, Journal of Commercial Vehicles, Section 2, Vol. 99, 1990, pp. 778-794. Modi, V . J . , Mokhtarian, F . , and Fernando, M.S.U.K., "Moving Surface Boundary Layer Control as applied to Two-dimensional Airfoils," Journal of Aircraft, Vol. 28, No. 2, 1991, pp. 104-112. Modi, V . J . , Fernando, M.S.U.K., and Yokomizo, T . , "Moving Surface BoundaryLayer Control : Studies with Bluff Bodies and Application," AIAA Journal, Vol. 29, No..9, 1991, pp. 1400-1406. Modi, V . J . , Shih, E . , Ying, B. and Yokomizo, T . , "Drag Reduction of Bluff Bodies Through Momentum Injection," Journal of Aircraft, Vol. 29, No. 3, May-June 1992, pp. 429-436. Modi, V . J . , Ying, B., and Yokomizo, T.,-"Effect of Momentum Injection and Fences on the Drag of a Tractor-Trailer Truck Configurations," AIAA 10th Applied Aerodynamics Conference, Palo Alto, C A , U S A , June 1992, Paper No. AIAA-92-2640. Brown, D. A . , "Peruvians Study Rotating-Cylinder Flap," Aviation Week and Space Technology, Vol. 88, No. 23, December 1964, pp. 70-76. Lighthill, M . J . , An Informal Introduction to Theoretical Fluid Dynamics, Clarendon Press, Oxford, 1986, pp. 174-182. Rauscher, M . , Introduction to Aeronautical Dynamics, John Wiley & Sons, New York, 211 1953, Chapter VI, pp. 244-252. Goldstein, S., Modern Developments in Fluid Mechanics, Vol. II, Oxford UniversityPress, Oxford, 1938, pp. 545-549. Chang, P. K . , Control of Flow Separation, Hemisphere Publishing Corporation, McGrawHill Book Company, Washington, 1976, Chapter 5, pp. 286-298. Modi, V . J . , and Wiland, E . , "Unsteady Aerodynamics of Stationary Elliptic Cylinders in Subcritical Flow," AIAA Journal, Vol. 8, October 1970, pp. 1814-1821. Wiland, E . , "Unsteady Aerodynamics of Stationary Elliptic Cylinders in Subcritical Flow," M.A.Sc. Thesis, Department of Mechanical Engineering, University of British Columbia, April 1968. Lamb, H . , Hydrodynamics, Cambridge University Press, Cambridge, 6 th Edition, 1932, pp. 99-102. Roshko, A . , "A New Hodograph for Free-Streamline theory," NACA TN 3168, July 1954. Parkinson, G . V . , and Jandali, T . , "A Wake Source Model for Bluff Body Potential Flow," Journal of Fluid Mechanics, Vol. 40, Part 3, pp. 577-594. Blevins, R. D., Flow-Induced Vibration, McGraw-Hill Book Company, New York, 1990. Simiu, E . , and Scanlan, R. H . , Wind Effects on Structures: An Introduction to Wind Engineering, John Wiley & Sons., New York, 1978. Houghton, E . L . , and Carruthers, N . B., Wind Forces on Buildings and Structures An Introduction, Edward Arnolds, London, U . K . , 1976. 212 Sachs, P., Wind Forces in Engineering, Pergamon Press, New York, 1972. Walshe, D. E . J . , Wind-Excited Oscillations of Structures, National Physical Laboratory, Department of Trade and Industry, Her Majesty's Stationary Office, London, U . K . , 1972. Bearman, P. W . , and Davies, M . E . , "The Flow about Oscillating Bluff Structures," Proceedings of the J^th International Conference on Wind Effects on Buildings Structures, London, U.K., Editor: K . J . Eaton, Cambridge University Press, pp. 195-222. Brooks, N. P. H . , "Experimental Investigation of the Aeroelastic Instability of Bluff Cylinders," M.A.Sc. Thesis, Department of Mechanical Engineering, University of British Columbia, 1960. Smith, J . D., "An Experimental Study of the Aeroelastic Instability of Rectangular Cylinders,"M.A.Sc. Thesis, Department of Mechanical Engineering, University of British Columbia, 1962. Heine, W . , "On the Experimental Investigation of Vortex Excited Pressure Fluctuation," M.A.Sc. Thesis, Department of Mechanical Engineering, University of British Columbia, August 1964. Ferguson, N . , "The Measurement of Vortex Induced Effects in Flow Past Stationary and Oscillating Circular and D-section Cylinders," M. A. Sc. Thesis, Department of Mechanical Engineering, University of British Columbia, 1965. Santosham, T . V . , "Force Measurements on Bluff Cylinders and Aeroelastic Galloping of Rectangular Cylinder," M.A.Sc. Thesis, Department of Mechanical Engineering, University of British Columbia, 1966. Feng, C . C , "The Measurement of Vortex Induced Effects in Flow Past Stationary and Oscillating Circular and D-section Cylinders," M.A.Sc. 213 Thesis, Department of Mechanical Engineering, University of British Columbia, October 1968. Slater, J . E . , and Modi, V . J . , "On the Wind Induced Vibrations of Structural Angle Sections," Proceedings of the 3rd International Conference on Wind Effects on Buildings and Structures, Tokyo, Japan, 1971, Editors: Japan Organizing Comittee, pp. 1195-1205. Novak, M . , "Aeroelastic Galloping of Prismatic Bodies," Proceedings of the ASCE, E M I , 1969, pp. 115-142. Zdravkovich, M . M . , "Review and Classification of various Aerodynamic and Hydrodynamic Means for Suppressing Vortex Shedding," Journal of' Wind Engineering and Industrial Aerodynamics, Vol. 7, No. 2, March 1981, pp. 145-189. Zdravkovich, M . M . , "Reduction of Effectiveness of Means For Suppressing WindInduced Oscillatio'n," Journal of Engineering Structures, Vol. 6, October 1984, pp. 344-349. Leschziner, M . A . , "Computational Modelling of Complex Turbulent Flow - Reality and Prospects", Proceedings of the First International Symposium on Computational Wind Engineering (CWE 92), Institute of Industrial Science, University of Tokyo, Tokyo, 1992, Editors: Japan Association for Wind Engineering, pp. 180-194. I Roache, P. J . , Computational Fluid Dynamics, Hermosa Publishers, Albuquerque, • U S A , 1972. Patankar, S. V . , Numerical Heat Transfer and Fluid Flow, Hemisphere Publishing Corporation, McGraw-Hill, New York, 1980. Anderson, D. A . , Tannehill, J . C , and Pletcher, R. H . , Computational Fluid Mechanics and Heat Transfer, Hemisphere Publishing Corporation, McGraw-Hill, Washington, 1984. 214 Zienkiewicz, 0. C , and Taylor, R. L . , The Finite Element Method, Vol. I & II, McGraw-Hill, New York, 1989. Chung, T . J . , Finite Element Computational Fluid Dynamics, McGraw-Hill, New York, 1978, pp. 170-265. Baker, A . J . , Finite Element Computational Fluid Mechanics, McGraw-Hill, New York, 1983, pp. 21-83. Reddy, J . N . , An Introduction to Finite Element Method, McGraw-Hill, New York, 1993. Abernathy, F . H . , and Kronauer, R. E . , "The Formation of Vortex Streets," Journal of Fluid Mechanics, Vol. 13, Part 1, 1962, pp. 1-20. Gerrard, J . H . , "The Mechanics of the Formation Region of Vortices Behind Bluff Bodies," Journal of Fluid Mechanics, Vol. 25, Part 2, 1966, pp. 401-413. Sarpakaya, T . , "An Analytical Study of Separated Flow about Circular Cylinders," Transactions of the ASME, Journal of Basic Engineering, Vol. 90, No. 4, 1968, pp. 511-520. Laird, A. D. K . , "Eddy Formation Behind Circular Cylinders," Journal of the Hydraulics Division, ASCE, Vol. 97, No. HY6, Proceedings Paper 8170, June 1971, pp. 763-775. Chaplin, J . R., "Computer Model of Vortex Shedding From a Cylinder," Journal of the Hydraulics Division, Proceedings of the American Society of Civil Engineers, January 1973, HY1, pp. 155-165. Sarpakaya, T . , "An Inviscid Model of Two-Dimensional Vortex Shedding for Transient and Asymptotically Steady Separated Flow Over an Inclined Plate," Journal of Fluid 215 Mechanics, Vol. 68, 1975, pp.. 109-128. [101] Fage, A . , and Johansen, F . C , "On the Flow of Air Behind an Inclined Flat Plate of Infinite Span," Proceedings of the Royal Society, A116, 1927, pp. 170-197. [102] Kiya, M . , and Arie, M . , "Discrete-Vortex Simulation of Unsteady Separated Flow Behind a Nearlly Normal Plate," Bulletin of the JSME, Vol. 23, No. 183, 1980, pp. 1451-1458. [103] Kiya, M . , Sasaki, K . , and Arie, M . , "Discrete-Vortex Simulation of a Turbulent Separation Bubble," Journal of Fluid Mechanics, Vol. 120, 1982, pp. 219-244. [104] Clements, R. R., "An Inviscid Model of Two-Dimensional Vortex Shedding," Journal of Fluid Mechanics, Vol. 57, 1973, pp. 321-336. [105] Nagano,.S., Naito, M . , and Takata, H . , "A Numerical Analysis of Two-Dimensional Flow Past a Rectangular Prism by a Discrete Vortex Model," Computers and Fluids, Vol. 10, 1983, No. 4, pp. 243-259. [106] Richards, W . B., "Discrete-Vortex Calculation of Flow Past a Plate," Report, NASA Lewis Research Center, 1984. [107] Kamemoto, K . , Naohiro, O., and Aizawa, M . , "Numerical Experiments on the Unsteady Separated Flow Behind Two Flat Plates in a Side by Side Arrangement Using the Discrete Vortex Method," Bulletin of JSME, Vol. 27, 1984, pp. 2118. [108] Stanby, P. K . , "A Numerical Study of Vortex Shedding Form One and Two Circular Cylinders," Aeronautical Quarterly, Vol. 32, Part 1, 1981, pp. 48-68. [109] Chein, R., and Chung, J . N . , "Discreet-Vortex Simulation of Flow Over Inclined and Normal Plates," Computers and Fluids, Vol. 16, No. 4, 1988, pp. 405-427. [110] Leonard, A . , "Vortex Methods for Flow Simulation," Journal of Computational Physics, 216 Vol. 37, 1980, pp. 289-335. [111] - Maull, D. J . , An Introduction to the Discrete Vortex Method, Cambridge University, Cambridge, 1981. [112] Leonard, A . , and Koumoutsakos, P., "High Resolution Vortex Simulation of Bluff Body Flows," First International Symposium on Computational Wind Engineering, Tokyo, 1992, Journal of Wind Engineering and .Industrial Aerodynamics, No. 52, August 1992, pp. 345-355. [113] Pan, L . S., Chew, Y . T . , and Lee, T . S., "A Proposed New Vortex Calculation Method Using Generalized Conservation of Circulation Method," Proceedings of the Sixth Asian Congress of Fluid Mechanics, Singapore, May 1995, Editors: Y . T . Chew and C. P. Tso, Vol. 1, pp. 445-448. [114] Bienkiewicz, B., and Kutz, R. F . , "Applying the Discrete Vortex Method to Flow about Bluff Bodies," Journal of Wind Engineering and Industrial Aerodynamics, Vol. 36, 1990, pp. 1011-1020. [115] Fernando, M . S. U . K . , and Modi, V . J . , "A Numerical Analysis of the Unsteady Flow Past a Savonius Wind Turbine," Journal of Wind Engineering and Industrial Aerodynamics, Vol. 32, 1989, pp. 303-327. [116] Ferziger, J . , "Simulating of Complex Turbulent Flows: Recent Advances and Prospects in Wind Engineering," Journal of Wind Engineering and Industrial Aerodynamics, Vol. 46 - 47, 1993, pp. 195-212. [117] Rodi, W . , "On the Simulation of Turbulent Flow Past Bluff Bodies," Journal of Wind Engineering and Industrial Aerodynamics, Vol. 46 - 47, 1993, pp. 3-19. [118] Franke, R., Rodi, W . , and Schonung, B., "Analysis of experimental Vortex-Shedding Data with Respect to Turbulence Modelling," Proceedings of the 7th Symposium on 217 Turbulent Shear Flows, Stanford University, 1989, pp. 24.4.1 - 24.4.5. [119] Deng, D. W . , Picquet, J . , Queutey, P., and Visonneau, M . , "Vortex-Shedding Flow Predictions with Eddy-Viscosity Models," Engineering Turbulence Modelling and Measurements 2, Elsevier, New York, 1993. [120] Nagano, Y . , and Tagawa, M . , "An Improved k-e model for Boundary-Layer Flows," Journal of Fluids Engineering, Vol. 100, 1990, pp. 33-39. [121] Jansson, S;, "Numerical Investigation of Steady and Unsteady Flows Comparing Turbulence Models and Different Near-Wall Models," Re-port ISRN CTH-TFD-PB- 92/1/SE, Chalmers University of Technology, Goteberg, Sweden, 1992. [122] ' Norris, H . L . , and Reynolds, W . C , "Turbulent Channel Flow with a Moving Wavy Boundary," Department of Mechanical Engineering Report, FM-10, Stanford University, 1975. [123] Ferziger, J . H . , "Higher Level Simulations of Turbulent Flow," Computational Methods for Turbulent, Transonic, and Viscous Flows, Editors: J . - A . Essers, Hemisphere, New York, 1983. [124] Canuto, C , Hussaini, A . , Quarteroni, A . , and Zang, T . A . , "Spectral Methods in Fluid Mechanics," Springer Verlag, New York, 1987. [125] Kato, S., Murakami, S., Utsumi, Y . , and Mizutani, K . , "Application of Massive Parallel Computer to Computational Wind Engineering," Journal of Wind Engineering and Industrial Aerodynamics, Vol. 46 - 47, 1993, pp. 393-400. [126] Murakami, S., "The Future of C F D in Civil Engineering: Large-Scale Computation with Vector and Massive Parallel Computers," U.S.-Korea-Japan Trilateral Seminar on Frontier R & D for Constructed Facilities, Honolulu, Hawai, 1991. 218 [127] Simon, H . D., van Dalsem, W. R., and Dagum, L . , "Parallel C F D : Current Status and Future Requirements," Parallel Computational Fluid Dynamics, Editor: H . D. Simon, M I T Press, Massachusetts, 1992, pp.1-29. [128] Hess, J . L . , and Smith, A . M . O., "Calculation of Potential Flow About Arbitrary Bodies," Progress in Aeronautical Sciences, Vol. 8, Pergamon Press, New York, 1966, pp. 1-138. [129] Farassat, F . , and Brown, T . J . , "A New Capability for Predicting Helicopter Rotor and Propeller Noise Including the Effect of Forward Motion," N A S A TM-74037, June 1977. [130] Lee, Y . J . , and Yang, J . Y . , "A Panel Method for Arbitrary Moving Boundaries Problems," AIAA Journal, Vol. 28, No. 3, March 1970, pp. 432-438. [131] Yon, S., Katz, J . , and Plotkin, A . , "Effect of Airfoil (Trailing-Edge) Thickness on the Numerical Solution of Panel Methods Based on the Dirichlet Boundary Condition," AIAA Journal, Vol. 30, No. 3, March 1992, pp. 697-702. [132] Kida, T . , Nagata, T . , and Nakajima, T . , "Accuracy of the Panel Method with Distributed Sources Applied to Two-Dimensional Bluff Bodies," Computational Fluid Dynamics Journal, Vol. 2, No. 1, April 1993, pp. 73-90. [133] Katz, J . , and Plotkin, A . , Low-Speed Aerodynamics - From Wing Theory to Panel Methods, McGraw-Hill Book Company, New York, 1991, 237-264. [134] Mokhtarian, F . , "Fluid Dynamics of Airfoils with Moving Surface Boundary-Layer Control," Ph.D. Thesis, Department of Mechanical Engineering, University of British Columbia, February 1988. [135] Keller, H . B., and Cebeci, T . , "Accurate Numerical Methods for Boundary-Layer Flows, Part 1, Two-dimensional Laminar Flows," Lecture Notes in Physics, 8, Pro- 219 ceedings of the Second International Conference on Numerical Methods in Fluid Dynamics, Springer-Verlag, New York, 1971, p. 92. [136] Keller, H . B., and Cebeci, T . , "Accurate Numerical Methods for Boundary- layer Flows, Part 2, Two-Dimensional Turbulent Flows," AIAA Journal, Vol. 10, 1972, p. 1193. [137] Cebeci, T . , and Smith, A . M . 0., Analysis of Turbulent Boundary Layers, Academic Press, New York, 1974, p. 170. [138] Cebeci, T . , and Bradshaw, P., Momentum Transfer in Boundary Layers, HemisphereMcGraw Hill, Washington, 1977, pp. 213-234. [139] Celik, I., Patel, V . C , and Landweber, L . , "Calculation of the Mean Flow Past Circular Cylinders by Viscous-Inviscid Interaction," Transactions of the ASME, Journal of Fluids Engineering, Vol. 107, June 1985, pp. 218-223. [140] Maskew, B., and Dvorak, F . A . , "The Prediction of C i m a : c Using a Separated Flow Model," Journal of American Helicopter Society, Vol. 23, , April 1978, pp. 2-8. [141] Pearcy, M . M . , Singh, S., Cash, R. F . , and Matten, R.B., "Fluid Loading on Roughened Cylindrical Members of Circular Cross Section," NMI Report 191, January, 1985. [142] Ribaut, M . , "A Vortex Sheet Method for Calculating Separated Two-dimensional Flows," AIAA Journal, Vol. 21, August 1983, pp. 1079-1084. [143] Mukherjea, S., and Bandyopadhyay, G . , "Separated Flow about a Wedge," The Aeronautical Journal of the Royal Aeronautical Society, June-July 1990, pp. 196-202. [144] Mukherjea, S., "Study of Incompressible Separated Flow Problems in Fluid Mechanics," Ph.D. Thesis, Department of Aerospace Engineering, I.I.T. Kharagpur, India, 1990. 220 [145] Lugt, H . J . , "Autorotation", Annual Review of Fluid Mechanics, Vol. 15, 1983, pp. 123-147. [146] Currie, I. G . , Fundamental Mechanics of Fluids, McGraw-Hill Book Company, New York, pp. 220-221. [147] Kellog, 0. D., Foundations of Potential Theory, Frederick Ungar Publishing Company, New York, 1929, pp. 277-338. o [148] Lanchester, F . W . , Aerodynamics, Constable, London, 1907, pp. 44. [149] Cheers, F . A . , "A Note on Galloping Conductors," National Research Council of Canada Report No. [150] MT-14, June 30, 1950. Edwards, A . T . , and Madeyski, A . , "A Progress Report on the Investigation of Galloping of Transmission Line Conductors," A.I.E.E. Transactions, Vol. 75, Part III, 1956, pp. 666-686. [151] Roshko, A . , "On the Drag and Shedding Frequency of Two-dimensional Bluff Bodies", NACA [152] TN 3169, 1954. Mair, W . A . , and Maull, D. J . , "Bluff Bodies and Vortex Shedding - A Report on Euromech 17", Journal of Fluid Mechanics, Vol. 45, Part 2, 1971, pp. 209-224. [153] Berger, E . , and Wille, R., "Periodic Flow Phenomena", Annual Review of Fluid Mechanics, Vol. 4, 1972, pp. 313-340. [154] Den Hartog, J . P., Mechanical Vibrations, McGraw-Hill, New York, Fourth Edition, 1956, p. 299. [155] Bearman, P. W . , "On Vortex Street Wakes", Journal of Fluid Mechanics, Vol. 28, Part 4, 1967, pp. 625-641. 221 [156] Sirovich, L., and Lim, C , "Comparison of Experiment with the Dynamics of the Von Karman Vortex Trail", Studies of Vortex Dominated Flows, Editors: M . Y . Hussaini and M . D. Salas, Springer-Verlag, New York, 1987, pp. 44-60. [157] Every, M . J . , King, R. and Weaver, D. S., "Vortex Excited Vibrations of Cylinders and Cables and Their Suppression," Ocean Engineering, Vol. 9, 1982, pp. 135-157. [158] Roshko, A . , "On the Development of Turbulent Wakes from Vortex Streets", NACA Re-port No. 1191, 1954. [159] Roshko, A . , "On the Wake and Drag of Bluff Bodies", Journal of Aeronautical Sciences, February, 1955, pp. 124-132. [160] Mclachlan, N. W . , "Ordinary Nonlinear Differential Equations in Engineering and Physical Sciences," Clarendon Press, Oxford, 2nd Edition, 1955, pp. 140-141. [161] "The Failure of Tacoma Narrows Bridge,"' Federal Works Agency, Washington, D . C . , U.S.A., March 1941. [162] Farquharson, F . B., "Aerodynamic Stability of Suspension Bridges with Special Reference to the Tacoma Narrows Bridge," The University of Washington Press, Seattle, 1949. [163] Steinman, D. B., Famous Bridges of the World, Random House, New York, 1953, pp. 69, 90-91. [164] Parkinson, G . V . , "Wind-Induced Instability of Structures," Philosophical Transactions of the Royal Society of London, A269, pp. 395-409. [165] Synge, J . L . , and Griffith, B. A . , Principles of Mechanics, McGraw-Hill Book Company, New York, 3rd edition, 1959, p. 435. [166] Hartlen, R. T . , and Currie, I. G . , "Lift-Oscillator Model of Vortex-Induced Vibra- 222 tion," Journal of the Engineering Mechanics, Proceedings of the ASCE, October 1970, E M 5 , pp. 577-591. [167] Iwan, W . D., and Blevins, R. D., "A Model for Vortex Induced Oscillation of Structures," Transactions of the ASME, Journal of Applied Mechanics, September 1974, pp. 581-586. [168] Blevins, R. D., Flow-Induced Vibrations, McGraw-Hill Book Company, New York, 2nd edition, 1990, pp. 67-74. [169] Verhulst, F . , Nonlinear Differential Equations and Dynamical Systems, SpringerVerlag, Berlin, 1990, pp. 15-16. [170] Minorsky, N . , Nonlinear Oscillations, D. Van Nostrand Company, New York, 1962, p. 219. [171] Nayfeh, A . H . , and Mook, D. T . , Nonlinear Oscillations, John Wiley & Sons, New York, 1979, pp. 3-6. [172] Thompson, J . M . T . , Bifurcatibnal Aspects of Catastrophe Theory, Annals New York Academy of Sciences, Vol. 316, 1979, pp. 553-571. [173] Thompson, J . M . T . , Instabilities and Catastrophes in Science and Engineering, John Wiley & Sons, New York, 1982, pp. 1-34, 155-182. [174] Stewart, I., "Applications of Catastrophe Theory to the Physical Sciences," Physica 2D, 1981, pp. 245-305. [175] Poore, A . B., and Al-Rawi, A . , "Some Applicable Hopf Bifurcation Formulas and an Application in Wind Engineering," Annals New York Academy of Sciences, Vol. 316, 1979, pp. 590-604. [176] Holmes, P., and Marsden, J . E . , "Qualitative Techniques for Bifurcation Analysis 223 of Complex Systems," Annals New York Academy of Sciences, Vol. 316, 1979, pp. 608-622. [177] Poore, A . B., Doedl, E . J . , and Cermak, J . E . , "Dynamics of the Iwan- Blevins Wake Oscillator Model," International Journal of Nonlinear Mechanics, Vol. 21, No. 4, 1986, pp. 291-302. [178] Hall, S. A . , and Iwan, W . D., "Oscillations of a Self-excited, Nonlinear System," Transactions of the ASME, Journal of Applied Mechanics, Vol. 51, December 1984, pp. 892-898. [179] Thompson, J . M . T . , and Stewart, H . B., Nonlinear Dynamics and Chaos - Geometrical Methods for Engineers and Scientists, John Wiley & Sons, New York, 1986, pp. 59-62. [180] Parkinson, G . V . , and Brooks, N. P. H . , "On the Aeroelastic Instability of Bluff Cylinders," .Transactions of the ASME, Journal of Applied Mechanics, June 1961, pp. 252-258. [181] Novak, M . , "Galloping and Vortex Induced Oscillations of Structures," Proceedings of the 3rd International Conference on Wind Effects on Buildings and Structures, Tokyo, 1971, Editors: Japan Organizing Comittee, pp. 799-809. [182] Richardson, A . S., "Predicting Galloping Amplitudes," A S C E Journal of Engineering Mechanics, Vol. 14, 1988, pp. 716-723. [183] Lawson. T . V . , Wind Effects on Buildings - Volume 1, Design Applications, Applied Science Publishers, London, 1980, Chapter 6, pp. 97-167. [184] Novak, M . , and Tanaka, H . , "Effect of Turbulence on Galloping Instability," ASCE Journal of the Engineering Mechanics Division, Vol. 100, pp. 27-47. 224 Kubo, Y . , Modi, V . J . , Yasuda, H . , and Kato, K . , "On the Suppression of Aerodynamic Instabilities Through the Moving Surface Boundary-layer Control," Journal of Wind Engineering and Industrial Aerodynamics, Vol. 41-44, 1992, pp. 205-216. Ou, Y . R., "Control of Oscillatory forces on a Circular Cylinder by Rotation," Proceedings of the Fourth International Symposium on Computational Fluid Dynamics, Davis, California, U.S.A., September 1991, Vol. 2, pp. 897-902. 225 A P P E N D I X A : I N S T R U M E N T A T I O N USED IN T H E S T U D Y (1) Pressure Measurement System ( S C A N I V A L V E C O R P O R A T I O N , USA) - Model J (500 psi), scanivalve transducer - CTLR2/S2-S6, solenoid controller - SCSG2, signal conditioner - accuracy, 2.5381 x 1 0 - 3 mV/N/m 2 (2) Linear Variable Differential Transformer ( C H A E V I T Z E N G I N E E R I N G , USA) - model 3000 HR-1722 - range = ± 3.0 inch (7.5cm), excitation = 3 volts @ 2500 Hz. - linearity = ± 0.14% of fullrange output (3) Real Time Analyzer ( S P E C T R A S C O P E , USA) - model SD335 - range = 20 Hz. to 50 K H z . - bandwidth = 0.06 Hz.(@ 20 Hz.) to 150 Hz.(@ 50 KHz.) (4) Variable Transformer ( O H M I T E M A N U F A C T U R I N G C O M P A N Y , USA) - catalogue no. V T 8 - F - input = 120 volts, 50-60 Hz. - output = 0-140 volts, 7.5A, D . C . (5) D . C . Power Supply (EPSCO INC., USA) - input = 115 volts, 50-60 Hz. - output = 0-16 volts, 8A maximum or 0-32 volts, 4A @ 0.1% ripple (6) High Speed A . C . Motors ( D U M O R E C O R P O R A T I O N , U S A ) - catalogue no. 6-021, Model 8551, Superflex-623G - Specifications: 115 Volts, 3.8A, 0-60Hz., 1/4 h.p., noload 22000 rpm 226 (7) Digital Tachometer (SHIMPO, JAPAN) - Model DT-205B - range = 6-30,000 rpm - accuracy = 2 rpm (8)" A / D Input^Output System (DATA T R A N S L A T I O N , INC., U.S.A.) - D T 2801 Series - Number of Channels = 8 or 16 - A / D Resolution = 12 bits - Maximum Gain = 8 - A / D input = 13.7 K H z . (9) I B M PC/486 Compatible Computer (10) Barocel Pressure Sensor ( D A T A M E T R I C S , U.S.A.) - Type 550-5 - Range 10 PSI (690 N / m ) 2 - Resolution, 10- 6 PSI ( 6 . 9 x l 0 - N / M ) 5 2 (11) Hotwire System ( D A N T E C ; Denmark, DISA Elektronik; Denmark) - D A N T E C Probe Type 55 P l l - Sensor Resistance at 20°C, R20 — 3.4fi - Leads Resistance, R i = 0.5Q - Sensor T C R , a 2 0 = 0.36% - Total Resistance, R = RTOTa2oR2o(T aor To) _ sen - measure KTOT at ambient conditions - DISA Type 55D 25 Auxilliary Unit - DISA Type 55D 31 Digital Voltmeter - DISA Type 55D 01 Constant Temperature Anemometer - DISA Type 55D 35 RMS Unit 227 APPENDIX B : SINGULARITY ELEMENTS AND INFLUENCE COEFFICIENTS B.l Velocity Induced by a Two-dimensional Straight Vortex Filament Figure B - l shows a straight vortex filament with a constant vortex strength T ' . According to the Biot-Savart Law, induced velocity due to an element dl is r di x f 6V = S 4^r Irl (B.l) 3 For an infinitely long (two dimensional) vortex filament the Biot-Savart law gives dl — rid/3/cos 8; and for —oo < I < + 0 0 we Substituting, r\ = rcos8, I = r\tanB, 2 get 0 < )3 < 7T, reulting in r 6V = d r / 4irr J 1 B.2 r (B.3) sinBd/3 = Q Panel with Uniform Vortex Distribution Figure B-2 shows a two-dimensional panel element composed of infinitely long vortex filaments distributed along the length of the panel. Each filament is perpendicular to the plane of the page. In this case one can apply equation (B.3). The strength of the vortex filament, ji per unit length, remains constant from one end of the panel to the other. The velocity components du e and dw e at the point (x,z) due to the small element of vortex distribution on d£ are: du (x,z) e = (BA) (*-0 dw (x,z) = e (x-iY 228 +* + 2 z\ 2 di. (5.5) Therefore: U {x,z) = e -iJ-Si/2 2TT Ii I tan * u (Is , z )\ = — (5.6) -1 + X e f z 27T u (x,z) i tan z =jiF(x,z,6i); e J-Si/2 (5.7) (5.8) (5-9) 27T .21 io (z:,2:) = - J - Z n e Si-, - - 1 ^ + 2 (5.10) 27T ( z + ^ +Z iu (a5,-z) = 2 (5.11) 'YiF(x,z,6 ). e i Since the panels are at different orientations, the expressions of velocity components must be defined with respect to the coordinate system (X,Z) of the body (Figure B-3). For the vorticity distribution along a panel i, the velocity components u and w, taken relative to the reference coordinate system (X,Z), are: u(X, Z) — yi^F(x,z,Si)cosdi w(X, Z) — ji^F(x,z,6i) — G(x, z, Si) sin$i|; sinBi - G(x, z, <$;) cos 0,-|; (5.12) (5.13) where x — (X — Xmi) cos.Qi + (Z — Zmi) sinOi] z = -{X-X )8in6i + (Z-Z dco8 mi Here X i m and Z i m th (5.15) are the mid-points of panel i in the (X,Z) system. The velocity components on the } due to the i 0i. m (-B-14) tfl panel (at the center of the panel, called the collocation point) panel are: Uji = ji(Fji cos 0i — Gji sin 6i)\ (5.16) Wji = ji(Fjisin (5.17) Bi + Gji cos 9i). For a collocation of N panels descretizing a body surface, this becomes: N N (5.18) »=i t=i 229 N W 3 N H X)(*y< ° i + GjiCosOi) = ^ s i n t=l Bjiji] (5.19) i=l where: Aji = A(x,z,6i,6i); (5.20 Bji = B(x,z,6i,6i)i (5.21 Xji = ( X y — Xmi) cos 9i + (Zmj — Zmi) sin 6n (5.22 m Xji = —(Xmi ~ Xrni) sin 6i + (Z j — Z i) COS 6i] m Si = [(X - Xi) + (Z 2 i+1 cos6i = ( X - Zi) ]; (5.24 2 i+1 i + 1 (5.23 m - Xi)/6i\ (5.25 sin8i = (Z -Zi)/6i. (5.26 i+1 B.3 Panel with Linear Vortex Distribution Consider a linearly varying vortex distribution across the panel i between - 6i/2 < x < Si/2 with strengths 7^ and 7i+i at x = ±Si/2, respectively, as shown in Figure B-4. The variation of the vortex strength across the panel can be described as 7 i = 7i±Lp._ii±ifyi t 2. ( B . 2 7 ) Oi Using equations (B.4) and (B.5), (B.14) and (B.15), and (B.27), the velocity components can be written as: u(X, Z)= l [ {(F(x, 7i+1 . -(G(x,z,6i) z, 6i) - zG(x, z, 6i)/( j) 6 - xG(x,z,6i)/( -±) 6 - xF(x, z, + zF(x,z,6i)/( -i) 6 +7i{(F(x, z, 6i) + zG(x, z, - cos 9i l/*)sin8i} + xF(x, z, Si) / (^))cos9i -(G(x, z, Si) + xG(x, 2, S i ) / ( | ) - zF(x, z, & ) / ( § ) + l / ^ m * ? , } ] ; 230 (5.28) w(X, Z)= l [n+i{(F(x, z, 6i) - zG(x, z, 6i)/( j) - xF(x, z, sin 6 6 +(G(x, z,6i)+ji{(F(x, xG(x, z, £)/(|) + zF(x, z, z, 6i) + zG(x, z, 6i)/(j) £)/(|) - + xF(x, z, -(G(x, z, 6i) + xG(x, z, 6 ; ) / ( | ) - zF(x, z, { l/ir)cos9i} Si)/(^))sin9i (5.29) + l/*)cos9i} The velocity components at the j-th collocation point due to N panels on the body surface can be written in the form: N i = u JV+l N ^Z im+-i. + J^&;»7» = a (5.30) 'ii > A 1=1 i=l i=l N JV JV+l i=l i=l 7i (5.31) jilfi'i where: (F(x, z, Si) -zG(x, z, -(G(x,z,Si) - xG(x,z„Si)/( j) fc)/(f) hi = 2 (F(x, z, Si) + zG(x, z, -(G(x, z, Si) + xG{x, z, &)/(|) - zF(x, z, {F(x, z, Si) - zG(x, z, +(G(x,z, Si) - xG{x, z, - + - xF(x, z, + zF(x, z, - 1 r Si. (F(x, z, Si) + zG{x, z, Si)f(^) + xF(x, z, -(G(x, z, Si) + xG(x, z, - zF(x, z, &)/(|) + Here: A'ji = ji-i + bji', a Bji = cos 9 { (5.32) l/ir)sin9i + xF(x, z, Si)/(^))cos6i ~ 2L da = xF(x, z, £)/(§)) + zF(x,z,Si)/(^) 6 C j i - ji-l + c 231 (5.33) l/ir)sin9i} sin 0 t lMcosBi (5.34) ,Si. 6i)/(^))sin9, l/^)cos9i} (5.35) - 'ji = B 'jN+l B fyu ji A = jN- c The quantities Xji, Zji, 6i, cosOi, and sinOi are given by equations (B.20) to (B.26). 232 233 234 APPENDIX C : POWER CONSUMPTION FORT H E MSBC The M S B C being an active control method a remark concerning the amount of energy consumption would be appropriate. A n aspect of importance is the relative gain in terms of power when the M S B C is used as a drag reduction device. Similarly it would be necessary to know the power required for lift augmentation on a wing. To answer such questions it was decided to obtain an estimate of the power requirements for real-life applications. The rotating cylinder used for the bluff body model had a diameter of 38mm, a length of 733mm and weighed 2043 gm. It was fabricated from solid aluminium. The power consumption was calculated from the measured voltage and current required to rotate the cylinder. The cylinder rotation rate was fixed at various values to obtain appropriate rates of the momentum injection (Uc/U). The results of the experiment are shown in Figure C - l . At the maximum rate of the momentum injection (Uc/U = 4), the power consumption is about 0.22 h.p. Most of the energy is expended in overcoming the bearing friction, the contribution due to aerodynamic resistance being marginal. For a full-scale truck, it is estimated that the rotating cylinder would be ten times as heavy (assuming a hollow cylinder made of thin sheet of aluminium). The power consumption also would be ten times as much, « 2.5 h.p., The truck engine power is around 400 - 500 h.p. Similar estimates for a rotating cylinder applied to a 60 storey (240m) tall building for vibration control would require a cylinder about 80m long (l/3rd the height of the building). Such a cylinder would weigh 100 times (sa 200 kg) and would need about 25 h.p.. For a building one requires cylinders at all the four corners, and thus the total power requirement would be about 100 h.p. Note, all the estimates are highly conservative. A rotating cylinder applied to the leading-edge of the wing of a Cessna type aircraft has been flight-tested to evaluate its performance. The results suggest that the power consumption is about 2% of the engine power. N A S A has flown the OV-10A aircraft built by the R O C K W E L L international, which was fitted with rotating cylinders at the leading-edge of the flap. The OV-10A has two engines, each 235 about 730 h.p., so' a 2% power requirement translates into about 15 h.p. per each rotating cylinder. 236 0.3 -i 1 r power consumption (h.p.) Cylinder Data diameter = 38 mm length = 733 mm mass = 2043 grams 2 3 rate of momentum injection ( U / U) c Figure C - l Experimentally measured power consumption of a momentum injecting rotating cylinder. 237 APPENDIX D : DYNAMIC SYSTEM PARAMETERS In the dynamic experiments to study the flow-induced vibrations it is necessary to determine the system natural frequency and damping factor. The set-up for the dynamic experiments is shown in Figure 2-9. The system natural frequency and damping factor were determined from a typical amplitude vs. time signal of the bluff body model undergoing free vibrations (Figure D-l). The signal sampling period (T) was 100 seconds and the sampling rate was 100 times per second. The F F T analysis showed that the natural frequency of the vibration (/ ) is 1.39 Hz (Figure D-2). The n amplitude of vibrations in Figure D - l is shown in terms of the computer output which can be calibrated in terms of the displacement. The calibration of the displacement transducer is shown in Figure D-3. The logarithmic decrement for the decaying vibrations is defined as, $ = _L-/ _d2_ p—1 n / i 0,1,2,..., = (D.l) where: AQ — initial amplitude of the vibration; Ay, — amplitude of vibration after /x cycles; and p — number of cycles. From Figure D - l , the logarithmic decrement 6 is determined to be 0.16923. The circular natural frequency w = 27r/ = 8.7336 rad/s. n n The stiffness (k ) of each spring (Figure 2-9) is 40,712 gm/s . There are four springs 2 x in parallel, giving a total stiffness (k) of 4x40712 gm/s . Figure D-4 shows the force 2 vs. displacement characteristics of the spring, as well as the operating range. The circular natural frequency is given by (D.2) The system effective mass (m ) as determined from equation (D.2) is 2,135 gm. e The damping factor (3) is given by P=\l 2- V 1 + (2*/6) 2 The value of 8 was found to be 0.002693. The system parameters may be summarized as follows: 238 (• model length; / i = 671 mm; model height, hi = 127 mm; system mass, m e =. 2135 grams; natural frequency, f n = 1.39 Hz; damping factor, 8 = 0.002693; total stiffness, k = 4 x 40712 gm/s ; 2 air density, p = 1.2 k g / m ; 3 mass parameter, m = ph\ lx/2m 2 e = 0.003042; Strouhal number of the square section model, St = 0.161; critical wind velocity for resonance, V = f h\/St c 239 n — 1.1 m/s. sampling rate = 100 Hz sampling time = 100 sec 500 1000 1500 2000 2500 3000 3500 4000 Non-dimensional time Figure D - l Amplitude vs. time history of the bluff body model (square-section) undergoing free vibrations. 240 180 sampling rate = 100 Hz 160 sampling time = 100 sec 1) -a 3 Oh S a c 140 120 100 .2 a 80 C3 o 60 40 20 0^= 0 0.5 1.5 2.5 Frequency (Hz) F i g u r e D-2 Natural frequency of the vibrating bluff body (square-section) as obtained through the F F T analysis of the amplitude vs. time history. 241 1500.0 h 1250.0 COMPUTER ~ OUTPUT 1000.0 - 750.0 h 500.0 -35.0 1 1 1 1 1 DISPLACEMENT, MM Figure D-3 Calibration of the displacement transducer 242 0.30 SPRING CONSTANT = 40.712 N/m 0.25 DISPLACEMENT (m) v ' 0.20 0.15 0.10 4.0 6.0 FORCE (N) Figure D - 4 Calibration of the spring used for the dynamic test set-up. ) 243 10.0
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Aerodynamics and dynamics of bluff bodies in presence...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Aerodynamics and dynamics of bluff bodies in presence of the moving surface boundary-layer control Munshi, Sandeep R. 1996
pdf
Page Metadata
Item Metadata
Title | Aerodynamics and dynamics of bluff bodies in presence of the moving surface boundary-layer control |
Creator |
Munshi, Sandeep R. |
Date | 1996 |
Date Issued | 2009-02-19T20:15:44Z |
Description | Moving Surface Boundary-layer Control (MSBC) was applied to several two dimensional bluff bodies using a high speed rotating cylinder as a momentum injecting device. Flow past a symmetric airfoil; a D-section; as well as square and rectangular prisms, representing a family of shapes with progressively increasing bluffness were studied in presence of the MSBC. In the case of the airfoil, the leading edge was replaced by a rotating cylinder; while the cylindrical element formed the top and bottom upstream corners of the D-section, square and rectangular prisms. Extensive wind tunnel investigation gave data about the effect of system parameters like rate of the momentum injection, angle of attack and the surface condition of the cylinder on steady and fluctuating components of the pressure distribution around the body, vortex shedding frequency (Strouhal number), and the lift and drag coefficients. A gain in the Strouhal number with increasing momentum injection suggest a decrease in the effective bluffness of the body. A significant reduction in the drag (up to 80%) was observed for the prisms at a maximum rate of momentum injection, Uc/U = 4 (Uc = cylinder surface speed, U = freestream wind speed). In the case of the airfoil, the lift coefficient increased by 160% and the stall angle was delayed from 110° to 48°. A rough criterion in terms of the location of the stagnation point was established to help decide the reversal in the direction of momentum injection as a function of angle of attack to ensure continued benefit. Effect of momentum injection in suppressing the vortex resonance and galloping type of instabilities were studied by mounting the bluff prism models on a dynamic-test rig inside the wind tunnel test section. The measurement of amplitude and frequency of the transverse oscillations over a range of wind speeds showed complete vibration suppression for momentum injection rates Uc/U < 2. Asymmetric momentum injection (e.g. top cylinder rotating, bottom cylinder stationary), was also found to be effective in disrupting the vortex shedding process and thereby inhibiting vibrations. The suppression of galloping instability in presence of the MSBC was also predicted by the quasisteady analysis. A numerical panel method was developed to simulate bluff body fluid dynamnics in presence of the MSBC. The body is descretized into a large number of panels (100 - 150) with each panel comprising of a continuous distribution of linearly varying vorticity and a constant source strength. A set of linear algebraic equations approximates the Fredholm type integral equation derived from ideal fluid flow assumption. The wake is modelled by upper and lower 'free vortex layers' emanating from the separation points on the body. Vorticity is allowed to be shed and dissipated as it is convected downstream along the panels on the 'free vortex layers'. An analytical expression relates the point vortex modelling a.rotating cylinder to the rate of momentum injection. The panel method is capable of treating multielement configurations (e.g. a rotating cylinder and the truncated airfoil). An iterative scheme based on the convergence of the wake shape is used to obtain the final solution. The numerically obtained pressure distribution, the lift and drag coefficients agree well with the experimental results. Flow visualization studies in a water channel were performed to obtain better physical insight into the MSBC process. Plexiglas models with rotating cylinders in conjunction with a fine suspension of polyvinyl chloride particles and slit lighting were used to visualize the streaklines. The still photographs and video movie recorded, rather dramatically, the effectiveness of the MSBC in suppressing separation and vortex shedding, making the flow approach the potential character. Overall, the present research firmly establishes potential of the MSBC as a versatile tool for lift augmentation, drag reduction and vibration suppression of several bluff bodies encountered in industrial engineering practice. |
Extent | 16506412 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
File Format | application/pdf |
Language | eng |
Collection |
Retrospective Theses and Dissertations, 1919-2007 |
Series | UBC Retrospective Theses Digitization Project |
Date Available | 2009-02-19 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0080849 |
Degree |
Doctor of Philosophy - PhD |
Program |
Mechanical Engineering |
Affiliation |
Applied Science, Faculty of Mechanical Engineering, Department of |
Degree Grantor | University of British Columbia |
Graduation Date | 1996-05 |
Campus |
UBCV |
Scholarly Level | Graduate |
URI | http://hdl.handle.net/2429/4818 |
Aggregated Source Repository | DSpace |
Download
- Media
- [if-you-see-this-DO-NOT-CLICK]
- ubc_1996-091392.pdf [ 15.74MB ]
- Metadata
- JSON: 1.0080849.json
- JSON-LD: 1.0080849+ld.json
- RDF/XML (Pretty): 1.0080849.xml
- RDF/JSON: 1.0080849+rdf.json
- Turtle: 1.0080849+rdf-turtle.txt
- N-Triples: 1.0080849+rdf-ntriples.txt
- Original Record: 1.0080849 +original-record.json
- Full Text
- 1.0080849.txt
- Citation
- 1.0080849.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Country | Views | Downloads |
---|---|---|
United States | 64 | 4 |
China | 17 | 36 |
Canada | 15 | 3 |
United Kingdom | 7 | 0 |
France | 5 | 0 |
India | 4 | 1 |
Japan | 4 | 0 |
Hong Kong | 2 | 0 |
Unknown | 2 | 0 |
Iran | 2 | 0 |
Sweden | 2 | 0 |
Tunisia | 1 | 0 |
Australia | 1 | 0 |
City | Views | Downloads |
---|---|---|
Washington | 40 | 0 |
Unknown | 20 | 13 |
Beijing | 7 | 0 |
Chengdu | 4 | 2 |
Boulder | 4 | 0 |
Edmonton | 4 | 0 |
Sacramento | 4 | 0 |
Mountain View | 4 | 3 |
Shenzhen | 3 | 34 |
Montreal | 3 | 0 |
Tokyo | 3 | 0 |
Air Force Academy | 3 | 1 |
Ottawa | 2 | 0 |
{[{ mDataHeader[type] }]} | {[{ month[type] }]} | {[{ tData[type] }]} |
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0080849/manifest