ANALYTICAL AND EXPERIMENTAL STUDIES OF WING TIP VORTICES By Shizhong Duan B. Sc. Tsinghua University, Beijing, 1984 M. Sc. Beijing University of Aeronautics and Astronautics, Beijing, 1987 A THESIS SUBMIITED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY in THE FACULTY OF GRADUATE STUDIES MECHANICAL ENGINEERING We accept this thesis as conforming to the required standard THE UNIVERSITY OF BRITISH COLUMBIA September 1995 ©Shizhong Duan, 1995 In presenting this thesis in partial fulfillment of the requirements for an advanced degree at the University of British Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that permission for extensive copying of this thesis for scholarly purposes may be granted by the head of my department or by his or her representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission. Mechanical Engineering The University of British Columbia 2324 Main Mall Vancouver, Canada V6T 1Z4 Date: I M5 Abstract Wing tip vortices, and their relationship to wing tip geometry and wing total drag, are investigated here both analytically and experimentally. The purpose of the analysis is to answer the basic question — what are the effects of wing drag on tip vortex structure? The purpose of the experiments is to determine the effect of wing tip geometry on wing tip vortices (tip vortex cavitation) and wing drag (vortex induced drag). In the analytical work, a new method has been developed for wing tip vortices. This approach is referred to as the “quasi-similarity” method. This new method combines a polynomial solution with a similarity variable technique. A non-linear analytical tip vortex model is achieved for the first time. The first order (linear) velocity components and pressure of the polynomial solution for the new model can be expressed in complete function form. Higher order (non-linear) velocity components and pressure can be obtained analytically or numerically. The most important feature of this new tip vortex model is that the predicted wing drag due to such a single vortex is finite. To the author’s knowledge, no other vortex model has this property. At distances fairly far downstream of the wing, tip vortex structure is found to be an explicit function of wing total drag, vortex circulation, freestream velocity, downstream distance and fluid properties. It is verified that the first order tangential velocity component and pressure of the tip vortex model are exactly the same as that of the linear wing tip vortex model proposed by Batchelor (1964). The decay of the axial and tangential velocity components predicted by this theory compare well with experimental measurements. In the experimental work, a novel ducted tip device was tested in a wind tunnel and a water tunnel. The ducted tip consists of a hollowed duct attached to the tip of a rectangular untwisted wing. This novel tip device was found to improve the Lift/Drag ratio by up to 6% at elevated angles of attack compared with a conventional round tip configuration with the same span. The wing tip vortex cavitation was substantially delayed by the ducted tip device. In view of its superior cavitation characteristics and aerodynamic performance, the ducted tip has potential application to marine propellers. H The ducted tip is effective because it redistributes the shed vorticity in the transverse plane (Trefftz plane) behind the wing. The shed vorticity is distributed along the wing and duct trailing edge for the ducted tip configuration, rather than solely along the wing trailing edge, which would be the case for a conventional wing tip configuration. ill Table of Contents Abstract ii List of Tables vii List of Figures viii Nomenclature xi Acknowledgements Xiv Chapter 1 Introduction 1 1.1 Motivation for tip vortex research 2 1.2 Previous measurements on tip vortex structure 4 1.3 Previously developed mathematical models of vortices 5 1.4 Scope of the present work 11 Chapter 2 Quasi-Similarity Method for Vortices in a Freestream 14 2.1 Overview of quasi-similarity method 14 2.2 The governing equations 15 2.2.1 General governing equations 15 2.2.2 Non-dimensional governing equations 17 2.2.3 Structure of the solutions 19 2.2.4 Governing equations of the first order 20 2.2.5 Governing equations of the second order 21 2.2.6 Governing equations of the arbitrary i-th order 21 2.3 The first order problem 22 2.3.1 Derivation of the governing quasi-similarity equations 22 2.3.2 Boundary conditions on the quasi-similarity equations 25 2.4 The second order problem 27 2.5 The i-th order problem 30 iv 2.6 The solutions in complete function form 2.7 Discussions of assumptions 36 2.8 Conclusions 38 Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 3.1 The first order solutions .35 39 39 3.1.1 General first order solution in complete function form 39 3.1.2 Tangential velocity and pressure 42 3.1.3 Effects of the constant 43 13 on the axial similarity functions 3.2 Drag due to a single quasi-similar vortex in viscous flow 44 3.3 The second order solutions 48 3.3.1 Second order tangential velocity and pressure 49 3.3.2 Second order axial and radial velocity components 51 3.4 The procedure for i-th order solution (i >2) 53 3.4.1 Conservation of the vortex circulation F 55 3.4.2 The i-th order drag 58 3.5 The axial velocity on the vortex centerline 60 3.6 Comparison with Batchelor’s vortex 60 3.7 Comparison with experimental data 64 3.8 3.7.1 Axial velocity decay on the vortex centerline 64 3.7.2 Maximum tangential velocity decay 65 Conclusions 66 Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 68 4.1 Previous work on tip modification 69 4.2 The ducted tip wing geometry 70 4.3 Experimental apparatus and procedures 72 V 4.4 .78 Results and discussion 4.4.1 Flow visualization 78 4.4.2 Cavitation characteristics 85 4.4.3 Aerodynamic performance 92 4.4.3.1 The partial chord ducted tip 4.4.3.2 The hi-wing tip 92 108 4.5 Experimental errors and applicability of the tip device 109 4.6 Conclusions 110 112 Chapter 5 Summary and Conclusions 5.1 5.2 Summary and conclusions 112 5.1.1 Analytical method 112 5.1.2 Tip vortex structure 113 5.1.3 Wing tip modification 114 115 Suggestions for future work 117 Bibliography vi List of Tables 4.1 Lift coefficients of the ducted-tip wing and the conventional wing 92 4.2 Drag coefficients of the ducted-tip wing and the conventional wing 93 4.3 Comparison of the hi-wing tip to the conventional tip at 8° angle of attack 108 4.4 Comparison of the bi-wing tip to the conventional tip at ct=14° 109 vii Lists of Figures 1.1 Helmholtz vortex laws interpretation of tip vortices 1 1.2 Shed vorticity from a finite wing 1 1.3 Rolling-up of a trailing vortex sheet 2 1.4 Induced drag acting on a finite wing 3 1.5 Burgers vortex model 6 2.1 Vortex in a freestream 19 3.1 The first order quasi-similarity functions. 13=0.361 42 3.2 Axial velocity excess in wing tip vortices 43 3.3 Axial velocity deficit in wing tip vortices 43 3.4 Control volume around the wing tip 44 3.5 The effects of wing total drag on the axial velocity 47 3.6 The effects of wing total drag on the radial velocity 47 3.7 The second order axial quasi-similarity function. 13<0 52 3.8 The second order axial quasi-similarity function. 13>0 52 3.9 The second order radial quasi-similarity function. 13<0 53 3.10 The second order radial quasi-similarity function. 3>0 53 3.11 Axial velocity decay with downstream distance in a wing tip vortex 64 3.12 Maximum tangential velocity decay with downstream distance in a wing tip vortex 4.1 65 Shed vorticity in the Trefftz plane from (a) conventional and (b) ducted tip wings 4.2 71 Schematics of wing tip devices. Flow is right to left. (a) Ring wing tip. (b) Bi wing tip 72 4.3 The Low Speed Wind Tunnel at the University of British Columbia 76 4.4 The data acquisition system for lift and drag measurement in the LSWT 76 4.5 The data acquisition system for pressure measurements in the LSWT 77 4.6 Tip vortex rollup in the Trefftz plane 78 4.7 Surface flow visualization on the round tip geometry VIII 79, 80 4.8 Surface flow visualization on the ducted tip wing 4.9 Separated flow near the entrance of the ducted tip. a=l2° 4.10 Developed trailing vortex cavitation behind the rounded tip geometry. a=7°, 82, 83 84 . a=l.8, DAC=5.4ppm 6 Re=l.3x10 4.11 88 Cavitation inception index versus angle of attack for the rounded tip geometry. , DAC=5.4 ppm. (Uncertainty in a=±0.2°, and in =±lO% at the 6 Rel.4xl0 95% confidence level) 4.12 89 Ducted tip vortex made visible by migration of cavitation nuclei, generated at the duct leading edge, into the vortex. . a=1.5, 6 a=7°, Re=l.2x10 DAC=7.Oppm 4.13 90 Cavitation inception index versus angle of attack for the ducted tip geometry. , DAC=7.0 ppm 6 Re1.4xl0 91 4.14 Drag coefficients of rounded and ducted tip wings. Re=7.lxlO 5 95 4.15 A simple model for the reduction in induced drag of the ducted tip geometry 95 4.16 Lift/Drag ratio improvement of the ducted wing relative to the conventional tip 96 4.17 Pressure taps in the surfaces of the basic wing 4.18 Pressure distribution on the pressure side of the wing with the conventional tip. a=8°. (a) the 3-D view, (b) the tip view 4.19 102 The difference between the ducted and conventional tip pressure coefficients on the pressure side of the basic wing. a=8° 4.21 101 Pressure distribution on the pressure side of the wing with the ducted tip. x=8°. (a) the 3-D view, (b) the tip view 4.20 100 103 Pressure distribution on the suction side of the wing with the conventional tip, c=8°. (a) the 3-D view, (b) the leading edge view, (c) view from the wing root 4.22 104 Pressure distribution on the suction side of the wing with the ducted tip, a=8°. (a) the 3-D view, (b) the leading edge view, (c) view from the wing root ix 105 4.23 The difference between the ducted and conventional tip pressure coefficients on the suction side of the basic wing, c8° 4.24 106 Pressure distribution on the suction side of the basic wing without a tip device, 107 x Nomenclature a A constant for the i-th order problem b Wing span 1 b A constant for the i-th order problem c Wing chord CD Drag coefficient Cf, A constant for the i-th order problem 1 C/i A constant for the i-th order problem c, A constant for the i-th order problem c, Pressure coefficient CL Lift coefficient D Drag due to a single tip vortex DAC Dissolved air content [ppm] D Wing drag f(q) An arbitrary function f(z) An arbitrary function (i) 1 F i-th order similarity function of axial velocity Oi) An arbitrary function (i1) i-th order similarity function of pressure 1 G (ri) i-th order similarity function of function 1 H ji K(x) Vortex circulation L A constant in Batchelor’s vortex 0 L Operator for ordinary differential equations L Operator for partial differential equations p Pressure Pv Water vapour pressure pr,,. Pressure at infinity P(z) i-th order amplitude function of pressure P() Similarity function for pressure in Batchelor’s vortex xi Pf Reference pressure Q A similarity function in Batchelor’s vortex Q2 A similarity function in Batchelor’s vortex R Radius of the control volume around a tip vortex r The radial distance from the vortex centerline 1 R The radius of a tip vortex Re Reynolds number R•(i) i-th order similarity function of radial velocity s Distance between two trailing vortices (i) 1 T i-th order similarity function of tangential velocity U Freestream velocity U(z) i-th order amplitude function of axial velocity (z) 1 V i-th order amplitude function of radial velocity z The axial distance along the vortex Greek Symbols A constant related to wing drag /r 9 A parameter defined as tangential velocity over the radius, 1 A positive number not far from unity Angle of attack p Fluid density v Kinematic viscosity Similarity variable TI Similarity variable Similarity variable Stream function 8 The azimuthal coordinate F Vortex circulation xii e 1 Tangential velocity p(z,r) An arbitrary function To Bound circulation at mid-span of a wing A constant for the i-th order problem a A constant for the i-th order problem 1 a Cavitation inception index, a 1 P(z) i-th order amplitude function of parameter N hr Radial velocity i Axial velocity co The axial vorticity (pt,. — XIII , )/o. 5p u 1 p Acknowledgements I would like to express my gratitude to my supervisor, Professor S. I. Green. His guidance throughout this program has been invaluable in finishing the work successfully. I wish to express my sincere appreciation to Professor I. S. Gartshore, Professor P. G. Hill, Professor S. M. Causal and Professor G. W. Bluman for their insightful discussions and suggestions. Thanks are also due to A. Schreinders, T. Besic, D. Camp, D. Bysouth and other technicians in the Mechanical Engineering Workshop for their skillful technical assistance in design and construction of the instruments and test models. I wish also to thank Professor G. S. Schajer for proof-reading my thesis. Financial support for this work by the Natural Science and Engineering Research Council of Canada and the University of British Columbia is gratefully acknowledged. I would especially like to thank my parents, and my wife Belinda for their patience, love, faith and encouragement. It is to them that I dedicate this work. xiv Chapter 1 Introduction Chapter 1 Introduction Trailing vortices occur whenever a lifting surface terminates in a fluid. A simple explanation of the existence of trailing vortices derives from the application of Helmholtz vortex laws. Consider a finite length wing impulsively started from rest. For - Flow the simplest case, a lift line (with a constant circulation along the wing span) must exist in order for the wing to generate lift force. Kelvin’s Theorem demands that this circulation Figure 1.1 Helmholtz vortex laws interpretation of tip vortices. be matched by an equal and opposite shed circulation (the “starting vortex”). Because vortex lines can never end in a fluid, these two vortices must be connected by tip vortices, as illustrated in Figure 1.1. A more realistic description involves a number of lifting lines distributed along the wing chord and span, as depicted in Figure 1.2. As a result of this distribution of lifting lines, a sheet of vorticity is shed by the wing. Such a vortex sheet is not stable; roll-up of the sheet occurs. The free edge of the sheet curls over, under the influence of the induced velocity field of the vortex sheet, and takes the form of a spiral with a continually increasing number of turns in the downstream direction (Figure 1.3). It has now been well established generated that by trailing wings vortices roll up exceptionally quickly. In general, roll up is complete, i.e. vortex circulation is virtually independent of downstream distance, 1-2 chords downstream of the trailing edge of an airfoil (Shekarriz et al. 1992, Stinebring et a!. 1991, Green and Acosta 1991). Contraction of vorticity in a plane transverse to the freestream during roll-up results in an axial velocity in the vortex core. After the shed vorticity is fully 2 Chapter 1 Introduction rolled up, there are only two concentrated trailing vortices behind a finite wing — consistent, coincidentally, with the conclusion that resulted from the simplest interpretation of Helmholtz vortex laws. 1.1 Motivation for Tip Vortex Research The wing tip vortex flow is one of great importance because it occurs often in practical problems. These problems include: landing separation distances for aircraft at runways — one airplane may experience a dangerous rolling moment if it accidentally flies into the tip wake of another (Donaldson 1971, Kantha et Figure 1.3 Rolling-up of a trailing vortex sheet. al. 1971, Snedeker 1972, Barber et al. 1975); blade/vortex interactions on helicopter blades, which is an undesirable source of noise and vibration (Martin et al. 1984, Summa 1982, Mosher and Peterson 1983, Lewy and Caplot 1982, Widnall and Wolf 1980); and propeller cavitation on ships due to the relatively low pressure in the vortex core (Ligneul and Latorre 1993, Arndt and Keller 1992, Chahine et al. 1993, Arndt et al. 1991). In the following airplane problem, the tip vortices generated by one aircraft may be sufficiently strong that a following plane accidentally entering into one of the vortices experiences a loss of control. This accidental interaction is most likely to occur during aerial refueling and over airport runways. Critchley (1991) has documented that serious incidents involving such interactions occur at an average rate of 9 per year at Heathrow airport! The cost of maintaining sufficient plane separation distances on landing to limit the frequency of such hazardous interactions is about $10 million annually for every major airport (Winter 1991). Helicopter rotor blades generate strong tip vortices. If the helicopter is hovering, following rotor blades will pass through the tip vortex generated by leading blades. In passing through the tip vortex, the following blade experiences large fluctuating forces Chapter 1 Introduction 3 (due to the unsteady aerodynamic loading) that can cause premature rotor blade fatigue and excessive blade noise (Poling et al. 1989). Finally, there is a tendency for tip vortices generated in marine applications to cavitate even at elevated inception indices (refer to section 1.4 for the definition of cavitation inception). The reason for the high values of cavitation inception index, , may be traced to the fact that the vortical motion is particularly strong (large 0 e / U) in tip vortices. In order to sustain these large tangential velocities (and hence extremely high centripetal accelerations), there exists a large radial pressure gradient from the vortex core to the surrounding fluid. If the freestream pressure is low enough, or the tangential velocity is high, the pressure in the vortex core will locally fall below the water vapour pressure and cavitation may occur. In every application, tip vortices act Induced drag as lifting surface inefficiencies. Flow downwash caused by the vortices - decreases the effective incidence angle of the wing. This tilting of the local incoming flow causes the lift force, which is perpendicular to the incoming Downwash Figure 1.4 Induced drag acting on a finite wing. flow direction, also to be tilted. The component of the lift force in the undisturbed freestream direction is the wing induced drag (Figure 1.4). Induced drag on wings is a particularly salient issue. Approximately 35% of the drag on a typical aircraft is lift-induced drag (Webber and Dansby 1983). Due to the small aspect ratio of marine propellers, induced drag on propellers is likely to be a greater percentage of the overall drag than in aeronautical applications. The potential savings to be reaped by reducing airplane induced drag, even fractionally, are staggering — roughly $100 million/year worldwide for each 1% induced drag reduction. Tip clearance flow in axial flow fans (Ruden 1974) and compressors (Raines 1954) decreases their efficiency. An understanding of the tip vortex flow is also pertinent to the design of propfan blades (Vaczy and McCormick 1987). Tip flow even plays an important role in the design of America’s Cup yacht keels (Devoss 1986). Chapter 1 Introduction 4 1.2 Previous Measurements on Tip Vortex Structure Attempts have been made to measure the tangential and axial velocity distributions around trailing vortices using Laser Doppler Velocimetry (LDV) (Orloff and Grant 1973, Baker et al. 1974, Higuchi et al. 1986a), hot wire anemometers (Corsiglia et al. 1973, Chigier and Corsiglia 1972, Zalay 1976), five-hole or seven-hole probes (McCormick 1968, Logan 1971, Chow et al., 1993), and Particle Image Velocimetry (PIV) (Green 1987, Shekarriz et al. 1992, 1993). Any intrusive method used to measure tip vortices is subject to large experimental errors because tip vortices are extremely sensitive. LDV is accurate locally, but there is a risk, due to its small measurement volume and the tendency of trailing vortices to meander, that the measurement volume may miss the target. This unsteady movement of trailing vortices may lead to a large overall measurement error when LDV is employed. Particle Image Velocimetry (PIV) has the distinct advantage that particle tracers move with the vortex. Although the local accuracy of Ply is not as high as LDV, the overall error of PIV is possibly the lowest among all the methods discussed above. Jn view of the sensitivity of vortices to disturbances and their tendency to meander, one must view with suspicion many of the early experimental studies of tip vortices. More recent studies, though, have produced some reliable and interesting results. For example, Pauchet et al. (1993) have measured vortex core axial velocities as big as 2.5 times the freestream velocity within two chord-lengths downstream of a particular hydrofoil. This axial velocity increases with the angle of attack, cx. Typical vortex core radii (the radius at which the tangential velocity is a maximum) are approximately 0.01 where , is the wing average chord (Green and Acosta 1991). This radius increases only slowly with downstream distance. Further downstream, viscosity leads to a slow diffusive increase of the core diameter. The gradual slowing-down of the azimuthal motion by viscous action leads to an increased pressure at the axis (there being less centrifugal force then), and so to an axial deceleration of the core fluid (Batchelor 1964). It is found experimentally that the Chapter 1 Introduction 5 tangential velocity decays slowly with downstream distance, whereas the axial velocity decays much faster. (Green 1991, Pauchet et al. 1993, Shekarriz et al. 1993). Maximum tangential velocities around the core in excess of (Arndt et al. 1991) have been measured. Much further downstream, trailing vortices typically “demise” after undergoing the Crow instability (Crow 1970). This instability takes the form of long wavelength sinusoidal disturbances on the vortex in fixed planes 48° from the wing planform. Eliason, Gartshore and Parkinson (1975) have verified experimentally the Crow instability. The amplitude of the disturbance grows in time. Following growth of the sinusoidal instability, the trailing vortices “link” in the regions where the two vortices approach most closely. Linking of vortices is followed closely by the formation of vortex rings. Vortex rings usually form hundreds or thousands of chord-lengths downstream of a wing; once vortex ring formation has occurred, the trailing vortex dissipates rapidly (Sarpkaya and Daly 1987). The wing tip vortices are also known to be highly unsteady. The axial velocity on the centerline can fluctuate with an r.m.s amplitude up to 0.2U, 0 (Green and Acosta 1991). Wing tip or strake vortices generated by delta wings may undergo an unusual phenomenon known as vortex “breakdown” or “bursting” (Peckham and Atkinson 1957). A stagnation point appears suddenly in the vortex core as vortex breakdown happens (normally at high angles of attack) and the flow downstream immediately becomes unsteady and irregular. The lift drops drastically and an asymmetrical rolling moment appears which may result in a loss of control. 1.3 Previously Developed Mathematical Models of Vortices Mathematical models of vortex behaviour have become increasingly more sophisticated and realistic with the passage of time. The earliest and simplest model of an isolated vortex is the “Rankine vortex.” The Rankine vortex model is a two region model of a vortex. The vortex core moves in solid body rotation (r t), and the region outside Chapter 1 Introduction 6 the core is irrotational. The vortex is purely two-dimensional and axisymmetric, with only a tangential velocity: ii.,(z,8,r)=O 0 2 TrI2itr rr F/2tr r>r “-ye (z,O ,r) = p—p(z,8,r)= 2 pF .r 2 pf’ r 2 4t 4 r 2 8it rr CpF2 r 2 8it r>rC The Rankine vortex does not satisfy the Navier-Stokes equations; at r = i the velocity gradient is discontinuous. Nonetheless, the pressure field result for a Rankine vortex is instructive — the pressure at the centerline falls with 1/r . The pressure (and thus 2 cavitation characteristics) of a Rankine vortex is very sensitive to r. The diffusion of a two-dimensional axisymmetric vortex was analyzed by Oseen(1912). At time t = 0 an irrotational line vortex (‘us = F / 2itr) is placed in a viscous fluid. Viscosity causes the vorticity to diffuse radially in time: =0 r(Z,0,r,t) = 0 (z,O ,r,t) = (F122tr) [1_exp(_r2 /4vt)] The “Lamb-Oseen” vortex is an exact solution to the Navier-Stokes equation. The solution shows that viscous growth of a vortex is very slow; the core size varies as The Rankine and Lamb-Oseen models are both vortex models for which the axial velocity is zero. More complex vortex models that incorporate axial flow in the vortex have ‘liz r 0 ‘ li been developed. The “Burgers Vortex” (Burgers 1948) models the flow around a vortex that is I0 z stretched in the axial direction. The three velocity components in the Burgers vortex are Figure 1.5 Burgers vortex model. Chapter 1 Introduction 7 (refer to Figure 1.5): u(z,r,t) = 2Cz i(z,r,t) = —Cr i. (z,r,t) f F =___.[l_exP_i-Jj This velocity field has associated vorticity distribution: r where 2 +62 I ex j-r _-JexP(_2Ct). The Burgers vortex solution reveals that axial stretching of vortex lines (e.g. as occurs immediately downstream of a wing) decreases the growth in vortex core size that arises from viscosity. Long (1961) extended the work of Burgers by considering a viscous vortex in an infinite liquid. Similarity arguments lead to a solution for this viscous vortex, r = = where = F r__, and f( 0 is the stream function. The similarity functions f, k, and s satisfy a set of ordinary differential equations. Long integrated these ODEs numerically, and found that the self-similar profiles have different characteristic shapes, depending on the value of the non-dimensional momentum transfer, M, given by: M= J (1 p 0 T --. + . 2itr dr j Batchelor (1964) has deduced the dynamical necessity of axial flow in a tip vortex core. In order to maintain the centripetal acceleration around the vortex, the pressure in the core of a trailing vortex is lower than that upstream of the wing. Thus, in the absence of viscous effects, fluid particles on a streamline originating ahead of the wing and entering the core are accelerated in the downstream direction (based on Bernoulli’s Chapter 1 Introduction 8 equation, flow velocity increases on a favorable-pressure-gradient streamline). Relative to the fluid at rest at infinity, the fluid in the core moves in a direction opposite to the direction of flight. In a coordinate system fixed to the wing, the component of velocity in the direction of the flow at infinity is increased in the core. Batchelor showed that a Rankine vortex with maximum tangential velocity u in the vortex core (1), )mac = 2 -.Ji + 2k U,, . = ku , 0 , has maximum axial velocity In Batchelor’s paper, a partly-linear theory was developed to predict the axial flow in a wing tip vortex core. The axial velocity deficit or excess was considered to be much smaller than the freestream velocity. Then, the tangential momentum equation could be linearized, and the following asymptotic form of the tangential velocity was obtained: 0 c (u_.r21 ;vz ] — where C 0 is a constant, U,,, is the freestream velocity, and v is the kinematic viscosity. By integrating all terms of the linearized axial momentum equation, Batchelor suggested the following asymptotic solution for the axial velocity in a wing tip vortex. where = 20 C 00 ZU 8vz v 2 8vz 172 U,., 8vz 2 U r and L is a constant related to any initial velocity defect which may be 4vz independent of the circulation. Batchelor studied the vortex wake far downstream, where it decays under the action of viscosity. The vortex radius is O(vz I U,., Batchelor found that viscous effects produce an axial velocity deficit in the far field. The principal limitations of Batchelor’s theory are that it is only valid at large distances downstream of the vortex-generating wing, and that the theory requires that the axial velocity deficit or excess be much smaller than the freestream velocity. As described previously, the latter assumption is often violated by the flow around real wing tip vortices. Moore and Saffman (1973) considered the structure of laminar trailing vortices rolling-up behind a wing. Their vortex structure consists of an inviscid roll-up region and Chapter 1 Introduction 9 a viscous vortex core. For an elliptical loading case, the wing is replaced by a bound vortex along the wing span, 0 <x < b, of strength K(x) = 2.. [x(b — x)] Kaden (1931) had remarked that roll up starts at the wing tips so that in the initial stages of the process, the rolled up portions at each end cannot significantly interact. Thus one can usefully consider a semi-infinite lifting line of circulation K(x) = 2y x, where = 0 b. The vortex sheet rolls up into a spiral with an infinite number of turns and the F velocity field near the centre of this spiral can be obtained by a simple argument — the velocity field, relative to the centre of the spiral, can be characterized by a distribution of circulation F(r) = 27[. r] where 2 is a constant. Thus if ij is the transverse velocity relative to the centre of the spiral, t r4 The pressure in the roll-up region is found by the radial equilibrium equation and is used to calculate the axial velocity, which is produced by streamwise pressure gradients. where n is a positive constant and Moore assumed further that the tip loading is 0 <n 1. The case n < = 1/2 represents elliptic loading. The axial velocity excess is then given by w=(1 —U provided w << U,. ) ic .(±_l.r2n 2Un atr0 ) is a constant. Note that w general argument. However w —* oo as r — > 0, in agreement with Batchelor’s 0. Therefore, the effect of viscosity has to be considered. Upon consideration of the viscosity, Moore and Saffman deduced that roll-up is 90% complete only tens of chords downstream of a wing, a result that differs from experiments (Shekarriz et al. 1992, Green and Acosta 1991). The discrepancy between theory and experiments may be due to their mathematical requirement that U, which virtually never occurs in practice. I U I Chapter 1 Introduction 10 A class of self-similar solutions of the steady, axisymmetric Navier-Stokes equations, representing the flows in slender vortices was reported by Mayer and Powell (1992a). Such self-similar vortices will have constant axial velocity in the vortex core if the external axial flow is a constant (i.e. In a wing tip vortex, the axial velocity in the vortex core decays along the downstream distance (based on experimental observations) while the freestream velocity is a constant, and hence wing tip vortices are not part of the family of self-similar vortices. The inviscid and viscous stability of a vortex, which has the asymptotic structure obtained by Batchelor (1964), was first studied in the papers of Lessen et al. (l974a) and Lessen & Paillet (1974b). They used a Runge-Kutta-type scheme along with the asymptotic behavior of the solutions at large and small radial distance to integrate the disturbance equations. A spectral collocation and matrix eigenvalue method was used by Mayer and Powell (1992b) to study the linear stability of Batchelor’s trailing line vortex. Their paper yielded the first global pictures of inviscid instabilities of the trailing vortex. Gartshore (1963) studied the quasi-cylindrical approximation of a trailing vortex using an integral approach. He found that a singularity may appear in the course of integration of the equations and the solution cannot be continued beyond the singular point. He suggested that these singularities may be mathematical indications of vortex breakdown. This work was continued by Mager (1972) who also used the integral approach to solve the quasi-cylindrical approximation equation. He made an attempt to connect the solutions he obtained with the concept of the breakdown as a supercritical-subcritical jump (suggested by Benjamin (1962)). Trigub et al. (1994) performed a consistent asymptotic study of steady axisymmetrical trailing vortices by using the quasi-cylindrical approximation. They presented numerical solutions for unbounded vortex breakdown parabolically expanding far downstream. As discussed previously, pairs of trailing vortices do not decay by simple diffusion. Usually they undergo a symmetric and nearly sinusoidal instability, until eventually they join at intervals to form a train of vortex rings. Crow (1970) developed a theory that describes the early stage of this instability. The equation relating induced velocity to vortex displacement gives rise to an eigenvalue problem for the growth rate of sinusoidal Chapter 1 Introduction 11 perturbations. The perturbation that grows most rapidly is a sinusoid of wavelength approximately 6 wing spans, which makes an angle of 48° to the wing planform. Most practical wing tip vortex flows are turbulent. The mechanisms involved in turbulent flows may be entirely different from those of laminar flows. Hoffmann et a!. (1963) predicted by theory, and confirmed by experiment, that the circulation in a turbulent vortex is proportional to the logarithm of the radius under certain conditions. They described a universal distribution of circulation in the inner region of the vortex, analogous to the turbulent boundary layer. Govindaraju and Saffman (1971) showed that a turbulent vortex must develop an overshoot of circulation to F,. as r —4 °°. — the circulation rises above F,. for a finite r and then falls back They stressed this prediction is independent of any hypothesis about the Reynolds stress. Any numerical calculation involving any closure approximation should demonstrate an overshoot of circulation, provided the closure approximation is selfconsistent with the conservation of angular momentum. The overshoot was found independently by Donaldson and Sullivan (1970), who studied the turbulent vortex numerically with a particular closure approximation. However, experiments do not indicate the presence of any significant overshoot in the circulation profile. 1.4 Scope of the Present Work The decay of wing tip vortices due to viscosity is investigated analytically in this thesis. As previously mentioned, Batchelor’s linear vortex model is the most advanced wing tip vortex model so far, and is valid far downstream of the wing where the axial velocity excess or deficit is much less than the freestream velocity. This downstream distance is estimated to be thousands of wing chords for a practical wing tip vortex (Moore and Saffman 1973). To understand the vortex structure fairly close to the vortex generating wing, non-linear effects must be considered. An analytical quasi-similarity method is developed for wing tip vortices. This new approach incorporates non-linear effects and velocity decay in a 3-D tip vortex. The central idea of the quasi-similarity method is that the solution consists of a polynomial in Chapter 1 Introduction which each term looks like a self-similar solution in complete function form 12 — an amplitude function times the similarity function. The amplitude function is a function of downstream distance only, and the similarity function depends on a similarity variable. The similarity variable in turn is a non-dimensional combination of the downstream distance and radius coordinates. The limitation of this new method is that the axial velocity deficit or excess must not be significantly greater than the freestream velocity. This is much less restrictive than the limitation on Batchelor’s linear theory. The geometric details of the wing that generates a tip vortex has a substantial influence on the vortex flow. Changing simply the airfoil section shape (e.g. from an NACA 0020 to an NACA 16020), while keeping the wing planform constant, can dramatically change the tip vortex flow (Pauchet et al. 1993). Similarly, changing the wing planform has a marked influence on the tip vortex (Fruman et al. 1995). For example, elliptic wing planforms tend to behave quite differently from rectangular planforms (e.g. compare Figure 4 of Arndt et al. 1991 with Figure 1 of Green and Acosta 1991). Changing the wing tip shape and even roughness (Green et al. 1988, Stinebring et al. 1991) can also dramatically alter the tip vortex. The search for a new wing geometry that avoids cavitation and/or reduces aerodynamic drag has been attractive to researchers. A novel ducted tip device, consisting of a small diameter flow-through duct, was attached to a rectangular untwisted wing and tested experimentally. Because the wing tip vortex generated by a ducted tip configuration would be expected to have a larger core than the tip vortex of a conventional tip, the pressure in the core would be higher than that of a conventional tip configuration. Therefore one might expect, and the expectation has been confirmed experimentally (Green et al. 1988), that the cavitation inception could be delayed by addition of a ducted tip device. A ducted tip is believed to reduce the induced drag of a wing as well, although it may cause a parasitic drag increase due to the extra wetted area at the tip. Experiments on two kinds of ducted tip devices mounted on a NACA 66-209 rectangular planform wing have been done in a wind tunnel. The two tip devices studied were a ring-wing tip and a ducted bi-wing tip. Several configurations were considered for Chapter 1 Introduction 13 each specific device. The lift and drag dependence on angle of attack, cc and pressure distribution on the wing surface were measured. All the results were compared with the baseline case of a rounded tip wing of identical span. The cavitation characteristics of the ring-wing tip device were investigated in a water tunnel. With the hydrofoil and tip device in the tunnel, gradually U was increased and decreased. At the moment of cavitation inception (as evidenced by the appearance of at least one macroscopic bubble per second in the trailing vortex core under stroboscopic illumination) IL. and p. were measured. The cavitation performance is characterized by a parameter called the cavitation inception index, which is defined as (pa, — p)/O.5pU, where p, is the water vapor pressure. The smaller the index, the better the cavitation performance. In the next chapter, the quasi-similarity method for vortices in a freestream (e.g. wing tip vortices) is discussed. Chapter 3 contains quasi-similarity solutions for wing tip vortex structure. The cavitation and aerodynamic performance of the novel ducted tip device is documented in chapter 4. Chapter 5 presents a summary and the main conclusions of the present work. Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream 14 Chapter 2 Quasi-Similarity Method for Vortices in a Freestream In this chapter, a quasi-similarity method for 3-D viscous vortices in a freestream (e.g. wing tip vortices) is discussed. Polynomial solutions for the velocity components and pressure distribution in such a vortex are obtained. The analysis shows that the tangential velocity in a tip vortex decays much more slowly along the downstream direction than the axial velocity, an observation that is consistent with experimental findings. 2.1 Overview of Quasi-Similarity Method The main idea of the quasi-similarity method is described briefly below. We start with the Navier-Stokes (N-S) equations in cylindrical coordinates. Some assumptions are made to simplify the N-S equations. The simplified governing equations, including the continuity equation, are a nonlinear system, which can be represented mathematically as (2.1.1) L is an operator for partial differential equations, and (2.1.1) can be either the continuity or momentum equations. The following polynomials are found to be the solution of (2.1.1), provided the summations are convergent. )=[1+.())] z (2.1.2a) (2.1.2b) 1=1 r(Zñ 1 . 5 (— ( 11) R p(z,r)-p =pU.-L.G(1l)J (2.l.2c) (2.1.2d) 1=I Here 1 = r is a “similarity variable” and the coefficients denoted by capital C are constants (we will verify that only one constant is independent in the first order of the polynomial, and two constants are independent in each higher order of the polynomials). F is the circulation of the vortex, v is the kinematic viscosity, U , is the freestream 0 Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream 15 velocity, z is the axial distance along the vortex, r is the radial distance from the vortex centerline, and p,,, is the pressure at infinity. Note that this is not a true similarity solution because the velocity components and pressure cannot be expressed solely in terms of i. Substitution of (2.1.2) into (2.1.1) yields, 1=1 Z [F(r), 1 0 L (1),G 7(ri),R ( r)] = 0 (2.1.3) where L 0 is the operator for an ordinary differential equation. Because z is an arbitrary positive variable, for each of i = 1,2,3,..., the following ordinary differential equation (ODE) must hold: [F(ri), 7(rI), J? (i), G, (i)] = 0 0 L (2.1.4) The ODE (2.1.4) can be solved with proper boundary conditions. Only one constant is involved in the first order problem. It will be shown in the next chapter that this constant in the first order problem (i = 1) is related to the total drag of a wing. The constants in each of the higher order problems (i 2) will be discussed later. It will be shown in Chapter 3 that the quasi-similarity solutions meet the following two requirements. One requirement is the conservation of the vortex circulation. The other requirement is the independence of wing total drag on the downstream distance. With this overview of the quasi-similarity method completed, we proceed in the next sections to a detailed discussion thereof. 2.2 The Governing Equations 2.2.1 General governing equations The incompressible Navier-Stokes equations in a cylindrical coordinate system are: Continuity equation: r r ao az (2.2.1) Chapter 2 Quasi-Similarity Method for Vortices in a Freestream 16 Tangential momentum equation: - + 9 -+V.V1) ( 1 ap 1.) rr. 2 h “1 (2.2.2) Radial momentum equation: ia -. +V•V1)r V=1)r 1 ‘ 1 2” (2.2.3) Axial momentum equation: at (2.2.4) paz - where, a ua - V’V=u 2 V —+-—-—+1) ar r O 1 a(r— a =—.—I r ar ar 1 a za 1 a 2 ra8 — 2 a az In the above equations p is the density,1.L the axial velocity, l)r the tangential velocity, the radial velocity, p the pressure, z the axial distance along the vortex, r the radial distance from the vortex centerline, and 9 the azimuthal coordinate. We can simplify the above equations by making the following assumptions for a vortex in a freestream: (1) The flow is incompressible, steady, and laminar (2) The flow is axisymmetric (3) r <<‘09 (4) az 2 << u9 2 a u 2 a ar- az- << ar Then, the governing equations can be simplified and rewritten as: !. a(-r) r ar az =0 a 9 a ar (2.2.5) r c’z 9 2 (a cir 2 9 ia rcir r (2.2.6) Chapter 2 Quasi-Similarity Method for Vortices in a Freestream 17 (227) par • “ar àz ar = paz ar + (2.2.8) rar It is convenient to define a new variable N = i I r. With this definition of N, (2.2.6) and (2.2.7) become: (2.2.9) (2.2.10) NJ2r=P The function i has a definite physical meaning. For an axisymmetric vortex flow, or r=O r=O where c is the axial vorticity. ji(r = 0) is thus half the value of the vorticity at r = 0. Because the vorticity is non-zero in a real vortex at r =0, so too is N• (In fact, because the central core of a vortex is forced by the action of viscosity into solid body rotation, the vorticity on the axis of a vortex is 22, where 2 is the angular frequency of the rotation.) This form of i is useful because it makes the tangential momentum and axial momentum equations take a similar form. It is therefore convenient in numerical calculations. 2.2.2 Non-dimensional governing equations The dimensional information in the governing equations (2.2.5)-(2.2.8) can be expressed as (2.2.11) z r rO ccv.- z r (2.2.12) r (2.2.13) r pr Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream lii) ccZZ r z 18 (2.2.14) locV- pz r The non-repeated information is listed below r (2.2.15) z (2.2.16) p’u (2.2.17) pocp1) (2.2.18) p cc Now we assume that the circulation of the vortex, f, is a constant along the axis of the vortex. Thus (2.2.19) .rccF 0 u Furthermore, it is reasonable to nondimensionize the axial velocity by the freestream velocity U. (2.2.20) cc Therefore, the variables describing a viscous vortex in a freestream can be non dimensionized using the following characteristic physical parameters — the freestream velocity U, the circulation of the vortex F, the fluid density p. and the kinematic viscosity v: F 2 U•v z= r= (2.2.21) (2.2.22) —. Uo (2.2.23) = = U•v p=pU = (2.2.24) r (2.2.25) • or (2.2.26) = where the variables with a symbol “—“ are non-dimensional parameters. The dimensionless governing equations are therefore: Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream 0 1a(Ir)z r ar 19 (2.2.27) = aN! r (2.2.28) = (2.2.29) or a- a— a 2 ciz ciz (2.2.30) m-2-or rcir The boundary conditions for the above equations will be discussed later. 2.2.3 Structure of the solutions In the derivation that follows, we will assume that a “quasi-similarity” solution UOO describes the vortex in a freestream. The phrase “quasi-similarity” is here used to describe a polynomial solution in which each term of the polynomial is a power product of Figure 2.1 Vortex in afreestream some function of the axial distance (the amplitude function), and a different function of the similarity variable (the similarity function). If a quasi-similarity solution exists for the flow in such a vortex in a freestream, one might expect the similarity variable to take much the same form as that of a free jet or boundary layer, i.e.: (2.2.31) With this choice of similarity variable, the quasi-similarity solution for the three velocity components and pressure is then posited to take the form: = l+U ().FO1) 1 r( ) = V (i). R (ii) = l+i (2.2.32) 1)r(j) (2.2.33) Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream = 20 ‘P)• H ) 1 (il) = j’J( 1 = = (2.2.34) (2.2.35) P(j) where the amplitude functions (those to be functions of only ) will be verified in section 2.5 to take the form: V()= C• , F()=2& P ( 1 )=-’, z (2.2.36) The coefficients denoted by capital C are constants. Substitution of (2.2.32) and (2.2.33) into (2.2.27), yields 1 [ 0 L 0 F(ii),R 1)]=0 ji+L where L 0 is an ODE operator. L[F(r),R (ri)j=0 must hold because 1 is arbitrary. Similarly, substituting (2.2.33)-(2.2.36) into (2.2.28)-(2.2.30) results in three other ordinary differential equations for each order (i.e. --, i = 1,2,3,.”) of the polynomial solution. 2.2.4 Governing equations of the first order The first order partial differential equations consist of all the terms inversely proportional to the distance that occur when (2.2.32)-(2.2.36) are substituted into (2.2.27)-(2.2.30): r ar a az =0 (2.2.37a) a 2 2 a (2.2.37b) •F= (2.2.37c) = _L a ( 2 i) +!. a(i3), F F Note that F / can be replaced by the similarity variable 2 (2.2.37d) 12 in the above equations. Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream 21 2.2.5 Governing equations of the second order Like the first order differential equations, the second order partial differential equations consist of all the terms inversely proportional to the square of the distance : 2 aQ) 2 a(i5r) O (2.2.38a) } (22.38b) (2.2.38c) P = Similarly, a2() 2/ } OZ)2 (2.2.38d) can be replaced by the similarity variable r 2 in the above equations. 2.2.6 Governing equations of the arbitrary i-th order Substitution of (2.2.32)-(2.2.36) into (2.2.27)-(2.2.30), and retention of all the terms proportional to --, yields the governing equations for the ith order problem: + + = r 2 a + — 0 J. (2.2.39a) z )j (2.2.39b) .[T’-i + 2]} +‘{(i3r)i (fr = . (2.2.39c) Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream 2() 1 +La) a(i), _ 1 a() _L+’f( Z 2 a 22 J — i-Il (2.2.39d) a(i):. ZJ +1(’Or).• In like fashion to the first order problem, the combination j2/ can be replaced by the sinlilarity variable r 2 in the above equations. 2.3 The First Order Problem 2.3.1 Derivation of the Governing Quasi-Similarity Equations To simplify subsequent mathematical manipulations, we state below some simple identities. First, we recognize that the following fundamental relations hold for the similarity variable defined by r = 1 ai 1 i ar z r 3 2 (2.3.2) Next, considering an arbitrary function p(, ) f() g(r), the following differential . relations are found to be valid: = 2z z = = r { 2 f’() g(i) -1. (2.3.3) g’(i) (2.3.4) . f(z) !Lf. g”(r) (2.3.5) Making use of the above relations, the continuity equation (2.2.37a) for the first order problem can be rewritten as 2Ji (z) 1 U Ro)+ RiO1)1 1 j 2 U) (z) 1 11 (-l’(i) = 0 (2.3.6) Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream 23 A “quasi-similarity” solution can exist for the continuity equation (Equation 2.3.6) if the equation can be made to be solely a function of the independent variable, Equation 2.3.6 is a function of 17 alone provided: z 1 (2.3.7) -.Iz (z) 1 U = Z_al 1 C (2.3.8) where C 1 and K are constants. The constant K does not affect the computed velocities, , a 1 and is set to 1 for simplicity. Using equations (2.3.7) and (2.3.8), the continuity equation can be rewritten as (1)+ 2R(11) 1 211R = 12. 1 Th (1) (11)— 2a (2.3.9) The tangential momentum equation (2.2.37b) can similarly be simplified to read 2 ‘“ H (ii) -iH(1) = 2H1)$ H(1) (z) 1 ‘P ii (2.3.10) A quasi-similarity solution exists for the tangential momentum equation if and only if (2.3.11) where C 2 and b are constants. The simplified ODE of the tangential momentum equation is thus, 1 (1)+[12+6].Hj(fl)_2bl.1.Hi(1)=0 21H (2.3.12) In similar fashion to the derivation of (2.3.6) above, the ODE form of the radial momentum equation (2.2.37c) can be shown to be — 1 G 11 ( )=Z lZ.l.H2(l) F(z) (2.3.13) for which a quasi-similarity solution exists only if (2.3.14) Finally, the axial momentum equation (2.2.27d) can be rewritten as 21.F(1)+[12 = P() i - (z) 1 U +2].1 (11)—2a 1 i(i) (2.3.15) (i)i .G’(l)j 1 .[2(1+2b)i.G 2 The following relation must hold for a quasi-similarity solution of (2.3.15) to exist, 24 Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream F(z) ) U ( 1 (2.3.16) = where 3 is a constant. The ODE form of the axial momentum equation is then, 2i I (11)+[12 1 •11I(1i) +2j.I (11)_2a (2.3.17) =[2(1+2b)1.G(1)_.T2 .G(r)] Combining (2.3.8) with (2.3.11), (2.3.14) and (2.3.16), one obtains F(!) 2 )_ (l_aj _i + 1 h Because (2.3.18) 1 C (z) 1 U (z) 1 U is a constant, we have (2.3.19) and =2b+1 1 a The constants b 2 can also be determined by applying the condition that the 1 and C circulation around the vortex far from the vortex axis is independent of downstream distance (a consequence of Kelvin’s Theorem): = (2.3.20) constant The dimensionless form of (2.3.20) is (F•). ( = ( =——, thus =1 = It can be verified that for the second or higher order, (2.3.2 1) (2 . , ) [12 1 (1])] H 0 (refer to section 4.1 in Chapter 3). Thus, (2.3.21) can be rewritten as 2 C .[12 1 and therefore, b thus = (ii)j 1 .H —1 1 =2b a +l=—1, 1 (2.3.22) =- and 1 C 2 C 21t[1]2 1 1 (ii)] H 47t213[12 (1)] 1 .H = (2.3.23a) =---= f3 2 (2.3.23b) To summarize, the following amplitude functions have been obtained for the first order problem: Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream )= U ( 1 1 42[2 H (n)] (2.3.24a) 2 I TJI—’_ Zj— 1 v 25 Jz ,— (2.3.24c) = 24112Hj = 21t[ where ()] 1 H (n)] (2.3.24d) 1 () U . f3 is a constant. 2.3.2 Boundary Conditions on the Quasi-Similarity Equations With all the unknown constants (except 6) evaluated, we may rewrite the governing ODEs that determine the first order similarity functions: (1) 1 ii)+ 2R 2iiR ( 1 = 2 ,, (2.3.25a) IOi)÷ 2i1• F(11) I 1 (rj)+21•H 1 (‘ri)+0i 2rj•H (rj)=O 1 2 +6)•H (2.3.25b) i) = i•H(11) G ( 1 (2.3.25c) 2 +2) 21•()+(fl H(i)] (2.3.25d) Equations 2.3.25a-d are four coupled ordinary differential equations in the four unknown , and second order 1 1 and G 1 ,G ,F 1 functions R 1 and H . The equations are first order in R 1 , , and two 1 , and consequently require one boundary condition on R and G 1 in F 1 and H 1 and H . The boundary conditions that apply to these equations are: 1 each on the F (2.3.26a) (O)=O 1 R 1 (0) = 1, H H(0) = 0 (2.3.26b) (2.3.26c) (oo)=0 1 G (2.3.26d) (0) is 1 (2.3.26a) must hold to ensure the flow is symmetrical about the vortex centerline. H set equal to one for simplicity (2.3.26b); a constant different from one will not affect the calculated tangential velocity. The second equation of (2.3.26b) ensures the second Chapter 2 Quasi-Similarity Method for Vortices in a Freestream 26 derivative at the origin has a finite value for symmetrical flow = 0). The first equation of (2.3.26c) holds because f3 in the amplitude function 1 U ( z) is a free constant. A non-unity value of the constant F (O) will not affect the final results. The second 1 condition of (2.3.26c) ensures the second derivative exists at the origin for a symmetrical flow (-ar = 0). (2.3.26d) arises from the freestream velocity condition. r=O The ODEs given by Equations (2.3.25a-d) form a decoupled system. We can solve (2.3.25b) for H (rj), then (2.3.25d) for 1 (rI) first, then (2.3.25c) for G 1 fj(11), and finally (2.3.25a) for R (rI), as will be discussed in Chapter 3. 1 The limit on 112 . 1 (i) in the amplitude functions (2.3.24a-d) can be evaluated. The H non-singular solution of (2.3.25b) with the boundary conditions (2.3.26b) is 4( (H 1 rl)——--I 1—e - (2.3.27) J Thus, [12 . =4 1 (TI)] H (2.3.28) The amplitude functions (2.3.24a-d) for the first order problem can therefore be simplified and rewritten as: (2.3.32a) 64it z V (z) P() JJ() lz = (2.3.32b) ..Jz .U ( 1 ) (2.3.32c) 1 () 1 =8i3.U 8it z (2.3.32d) — 64ic . . where 3 is a constant. There is only one constant, 13, which is undetermined in the first order problem. We will show in the next chapter that this constant is closely related to the total drag of the wing that generates the wing tip vortices. Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream 27 2.4 The Second Order Problem The ODE form of the second order continuity equation (2.2.38a) is similar to the form of the first order continuity equation: 2iR(i)+2R ( 2 fl)= 12. F’(1)—2a 2 11F(1i) where a., = z - (2.4.1) ,(z) () 7 U (2.4.2) A quasi-similarity solution exists for (2.4.1) if and only if (z)= 2 V (2.4.3) ,— qz The tangential momentum equation for the second order problem (2.2.38b) can be simplified to: 2P) H,(1)-1H(11) = 2H’(1)$H(1) 1 - 1 ((i) U (244 {[ii + 4H 1 (11)]R 1 (1)- [21H, (1) + 12H(1)]1 (1)} The following relation has to hold for a quasi-similarity solution of (2.4.4) to exist. ()P 1 U ( ) ‘P,(z) where thus 2 =a 2 (2.4.5) is a constant. •W() (z) 2 ‘P •U ’ 1 () (z) 1 I’ U(z) b =—2 (2.4.6) The tangential momentum equation thus can be rewritten as 2 +6)•H(1)+41•H 21•H’(1)+(1 (1) 2 = c f{2iiH(1) + 4IJ (n)f R1 (1) _[12 (2.4.7) H() + 21H 1 (1)] F (ii)} The radial momentum equation (2.2.38c) can be rewritten as the following ODE: (2.4.8) G(i) = ( 2 ( 1 21H 1)•H 11) Quasi-similarity solutions exist for (2.4.8) as long as 2 ( 1 )=P ( )P ) ±.p().U() Finally, the axial momentum equation (2.2.38d) can be rewritten as (2.4.9) Chapter 2 Quasi-Similarity Method for Vortices in a Freestream 2 211F ()÷O.2 ‘ — (z) 2 U + U, (z) 2 U + 2)F’(q)— 2 fri) 2a ( lF 28 = (2.4.10) 1 (11) •112 (i)] H 3 [4o) + 2ii (1)1’(ii) iiI (ia) [2I (11) + ii1’(ii)j} - A quasi-similarity solution exists only if the following relations hold (2.4.11) (z) 2 U U) U-,(z) 2 =c (2.4.12) where J 2 are constants. 6 and c thus — 2 a ZUC) () 2 U = 2 ’ 1 2.U () () 1 U (2.4.13) By applying the relations (2.4.9), (2.4.11), (2.4.12) and (2.3.16), one finds F) 1 F 2 c ( (z)-U z) C213 (2.4.14) Comparing (2.4.11) with (2.4.14), one obtains C (2.4.15) = f2 To summarize, the following relations must be satisfied in order for quasi-similarity solutions to exist: ()=±.U=( U ) 2 i 6 .1 (2.4. 16a) 2 C ( 1 1 (2.4.16b) 2 C 2 2 13 1 .P().U()=) 2 C2 12 2 •c uJ 1 (2.4.16c) 2.C 2 I(U)=(M2) 2 8ic (2.4.16d) •2c Notice that there are only two independent constants in the second order problem (except fi. It is permissible to select c 2 = 1 and a = 1, and let 112(0) = 2 and F Ch (0) 2 = , 2 Cf 29 Chapter 2 Quasi-Similarity Method for Vortices in a Freestream without loss of generality. From equation (2.4.15), /3 = /3. Therefore, the amplitude functions can be expressed as () =U() 2 U (2.4.17a) () =— 2 V =•U) 1 (2.4.l7b) 41Z P,(z) = 2 (z) P = 3•U(z) (2.4.l7c) 1 (z) 8t13 U (2.4. l7d) The second order quasi-similarity solutions of the vortex flow in a freestream can be obtained by solving the following ODEs with their boundary conditions: (1)= 2 2r1R(1)+2R 2 .F’(r)+4rj. F (ri) 2 (2.4.18a) +6)•H,’(1)+4T1.H 2r1•H,”(1)+(T1 ( 2 1) = 1 (‘rj)j. R 1 (ii) _{12 1 1 (ii)]. F (rj) {2riH(i) +4H H ( 1) + 2rH G(i) 21 F, (2.4.18c) (1)H 1 211H ( 2 1) (i) + (2 + 2) ‘(ii) + 4 . F () (rj)+2i,3 2 2 =-_.[4r.G (ri).H 1 .H ( i)] + 2rR (r) j7’() 1 . where — (2.4.1 8b) (2.4.18d) rF (ri).[2F (i) + /3 is the constant that appeared in the first order problem. The boundary conditions of this second order problem are (O)=O 2 R (O)= C/i 2 H 1 1(O) = Cf, (2.4.19a) = = constant, constant, H(O)=O (2.4.19b) F’(O) = 0 (2.4.19c) (oo)=0 2 G (2.4.19d) Analogously to the first order boundary conditions, (2.4. 19a) must hold to ensure the flow is symmetrical about the vortex centerline. The constant in the first equation of (2.4. 19b) is related to the second order vorticity at the origin r = 0. The second equation of (2.4. 19b) ensures the second derivative at the origin has a finite value for symmetrical flow = ar 0). The constant in the first equation of (2.4.19c) represents the amplitude r=O of the second order axial velocity on the vortex centerline. The second condition of Chapter 2 Quasi-Similarity Method for Vortices in a Freestream 30 (2.4.1 9c) ensures the second derivative exists at the origin for a symmetrical flow (- 0). (2.4. 19d) arises from the freestream velocity condition. = r=O There are two constants, Ch 2 and Cf 2 (in addition to /), involved in the second order problem. They will be discussed in the next chapter. 2.5 The i-th Order Problem (i = 2,3,4...) We will prove by mathematical induction the proposition that the i-th order quasisimilarity solution of the tip vortex can be obtained by solving the following ODEs with their boundary conditions: 1 2riR’+2R 12 (2.5.la) •F’+2i•rj•1 1 21H[’+ (112 + 6)II’+ 2111 Fl ‘— 1 + 4fl )R 2w 1 -( + 2(i - j)11n )i 1 (2.5. ib) = (2.5.lc) .H+ G1=i(H ) 1 2 +2)’+2i1 2iF/(i = _.[2iflGi +11(H .H )] + 1 j=1 +{2iR where /3 is -l . (2.5.ld) +‘j} the constant that was established in the first order problem solution. Like the second order problem, the boundary conditions of this i-th order problem are: R(0) = (2.5.2a) 0 (O)=ch H , 1 H,’(O)=O (2.5.2b) , 1 cf F’(O) =0 (2.5.2c) I (0) = G,(oo)=O (2.5.2d) The following relations have to be satisfied if the quasi-similarity solutions exist: = U,’() (2.5.3a) 31 Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream l’ () = =. (2.5.3b) U () (2.5.3c) ‘{ (2.5.3d) () = 8it• U () We begin our proof by first substituting the quasi-similarity parameters into the governing equations of the i-th order problem (2.2.39), yielding: ).2z (2R[+2R.)+U ( 1 (z) 1 U • 2.’I’() -‘ I 2P.() 0 (2.5.4a) .FJ1 (2.5.4b) ] 1 .11H Ti H — 1 •1 • 2 ].R 1 +4H } (2.5.4c) U().[22 +6)•’- 2•U) [2’()G2G]P() (2.5 .4d) F 2 .R]} Now, we proceed with the proof by induction by checking that the proposition is correct for i = 2. For i =2, we have shown in section 2.4 that ) U ( 2 = ) U ( 2 1 V,() =—=•U() (2.5.5a) (2.5.5b) Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream P,(!) 2 = 32 f3.U) (2.5.5c) () = 8itf3. U () (2.5.5d) while the quasi-similarity functions satisfy the following ODEs, 2 2rjR+2R 2ri• H’+ (2 F’+4iF 2 r +6) H + 4ri F! 2 ’+4H 1 =[211H j .R, G = (2.5.6a) _[12 (2.5.6b) .H[+21,H j 1 .F (2.5.6c) 1 H 2 21 H ” 2 211F + 2)•F’+4ri.F +(’ri j +21i3.H 2 =—13.[411•G . 1 H (2.5.6d) 1 •F’—llF .{2F +2rR +iiF’j The proposition is obviously correct for i =2. Next, we assume the proposition is correct for i = 1, 1 k. This means that for the 1- th order problem: U, () = U () (2.5.7a) () =_J=.U(C) 1 V (2.5.7b) p () f3. U ’ () 1 (2.5.7c) 8ic3• U’ () (2.5.7d) P, () = = and the quasi-similarity functions satisfy 2TjR[+2R, =12 (2.5.8a) F’-f-21•1F 2 + 6)H[+ 2liJl, 21H[’+ (rj = {(2TIH + 4H,)R — (112H + 2(1 — j)i,H,)F} G=n.(H .H,÷ ) 1 21F’- (12 + 2)J’+ 211F (2.5.8b) (2.5.8c) = [3 .[2liGi + ii(H ÷i)] 1 H (2.5.8d) - . .{2(l_j) +TI’j} Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream 33 Now, we prove the proposition is correct for i = k + 1. The general governing equations for (k÷ 1)-st order problem (2.5.4), replacing i with k+ 1, can be written as • k+1 [2 Uk’ (z) () (2R÷ + 2Rk÷ ) + Uk+l () 1 () [21H k÷1—1 = j=I { +(12 ‘ + 6) k+I — 2 . :+ () w .lHk+l] ‘k÷I() 12. H+ii]. .P().[ U ‘Pk+1J(z) k÷1—l ( + (2.5.9a) _12. Uk÷l() k+lJ(Zi.[2h1’k÷lJ } (2.5.9b) } 1 J]• R 1 +4Hk÷ j=l k÷1 1 G÷ () {pj () = k+l-j+1 () Hk+l . j=l . [21 Uk÷l + (12 + — Uk+1 L •11•Gk+j (2.5.9c) (z) 12 +() k÷1—l + () 1 E2 .U÷ UJ().Uk+lf().[ Uk+IJ() L. I—I + j} 2 U’÷ 1 () ‘ 6) k÷l k+l-j+l 11 {. ]} (2.5 .9d) F ).uk÷lJ().[21F;+lJ .R]} By using (2.5.7), it is apparent the following equations must be satisfied for quasi similarity solutions to exist, (2.5.lOa) Uk+l() Ck +1 (2.5. lOb) = C k+I Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream () = ‘k+l(Z) f3 c • 34 () U () (!)•U)= 1 =—---Y (2.5. lOc) ‘ (2.5.lOd) k+I where Ck+l , ,6+ and a’k÷l are constants. ,8 is the constant appearing in the first order problem. Like the second order problem, there are only two independent constants (except /5j. It is permissible to select Ck+I=l and ak+=l, and at the same time let 1 (O)=Ch H and F (O)=Cf without loss of generality. Thus fk÷ 1 1 = ---- = J3. Therefore the amplitude k÷1 functions (2.5. lOa-d) can be rewritten as (2.5.lla) Uk+l() =U’() +I 1 () u () = (2.5.1 ib) (2.5.llc) () 1 =U (2.5.lld) Thus, (2.5.9) can be rewritten as: (2.5.12a) 1+ 2TjR÷ 1 +2(k+1)1F+ 1 Rk+l 1 •F’÷ 2 1 2i1H + [2 + 6]H÷ 1 + 2(k + 1)ilHk÷l k+1—1 (2iI-I+_ 1 +4Hk+l_J)RJ [12Hz (2.5.12b) +2(k+1—j)’flHk+{JIFJI = 1 G÷ = ‘ri (H Hk÷I_J÷,) . (2.5. 12c) 35 Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream 1 1 +(r2 +2)Fk’+I +2(k+1)1F,+ 2i1F = ( 1 + + l)1G÷ . ] i+ 1 H÷ ) (2.5. 12d) k+1—I + - iF .[2(k + 1- j)F + where j3 is the constant of the first order problem. The proposition is correct for i = k + 1, provided that it is correct for i =1, 1 k. Since the proposition is correct for i=2, this concludes the proof by mathematical induction. 2.6 The Solutions in Complete Function Form Combining all the results obtained previously, the quasi-similarity solution of vortex flow in a freestream is given by: (z,r)=U .[1+U(z).(i)] .U((z).R,(i) 1)r(z,r) = i1) U(z)H ( 8 1 W(z,r)= or (z,r) 8 = 8•U VZ U(z)() (2.6.1) (2.6.2) (2.6.3) (2.6.4) where 7(ii)=iiH (ri) 1 p(z, r) where r = — p. r!—-, V Vz (z)= 1 U = 1 (11) pU13 U (z) G . (2.6.5) and v•z 64 U 2 (2.6.6) Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream 36 where 3 is a constant. The solutions (2.6. 1)-(2.6.5) are convergent provided that the similarity functions are of the order of unity and 1 U (z) <1 (2.6.7) The quasi-similarity functions of arbitrary i-th order H(rI), G ) and R 1 (), F(r 1 (1) 1 can be solved one by one from lower order to higher order by solving the ODEs (2.5.1) with the boundary conditions (2.5.2). Details will be discussed in the next chapter. 2.7 Discussion Of Assumptions Five assumptions have been made to derive the governing equations (2.2.5) - (2.2.9) used in this analysis. Those assumptions are: (1) incompressible, steady, laminar flow (2) axisymmetric flow (3) 1•) (4) (5) v U U 2 az v U 2 z the circulation of the vortex is independent of the downstream distance The first two assumptions are straightforward, although they will be violated by any mechanism that causes vortex instability. Consider now the third and fourth assumptions. Using (2.6.2) and (2.6.4), one finds IU —°i V I . z v 1oc—.U(z) r . 8iqU.vzJ (2.7.1) F Also, be:ause T = rifZ, (2.3.3) together with (2.3.5) yields (2 a2 ar 2 j-’j 2 a2 j’} —‘ ar- 2 (2 j’} i(2 4• r’} U (z) (2.7.2) r—2 One can verify that the terms in the governing equations neglected due to assumptions (3) and (4) are proportional to 3.) 1 (z). U In contrast, other terms in the governing Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream 37 equations are of order unity. Assumptions (3) and (4) are justified provided that the vortex circulation is much greater than the kinematic viscosity of the fluid. Because the kinematic viscosity is commonly very small (V v lO m 2 / s for the laminar flow of water, 2 Is for the laminar flow of air), these assumptions will almost always be m 5 10 justified for real wing tip vortex flows. The fifth assumption downstream distance — that the circulation of the tip vortex is independent of — is a consequence of the Helmholtz vortex law requirement that vortex lines not terminate in a fluid. The assumption could be violated if vortex lines doubled back on themselves to form vortex rings, but this behaviour does not occur in the quasi-similarity solution. Alternatively, it can be viewed as a consequence of Kelvin’s Theorem (DF I Dt = 0). Kelvin’s Theorem applies here because on a material loop far from the vortex axis there is negligible effect of viscosity, and this flow is incompressible and without body forces. The assumptions made in the derivation of the quasi-similarity solutions are thus not particularly restrictive. Rather, convergence problems of the quasi-similarity solutions more significantly limit their application. In particular, these limitations can be represented mathematically by, IUt(z)I< 1 (2.7.3) The first order axial velocity excess or deficit in the vortex in a freestream is thus constrained to being less than the magnitude of the freestream velocity. Finally, let us discuss the quasi-similarity solutions. The quasi-similarity expressions (2.6.l)-(2.6.5) can be verified to be the solution of the governing system (2.2.5)-(2.2.8) consisting of the continuity equation and simplified N-S equation. Substitution of (2.6.l)-(2.6.5) into (2.2.5)-(2.2.8) yields: (2.7.4a) 2z {u; ()• M (I1 ,f3)} 0 (2.7.4b) (2.7.4c) 38 Chapter 2 Quasi-Similarity Methodfor Vortices in a Freestream (2.7.4d) H,)}O M,G , ±.{u(). 1 provided U ( < 1, because• 1 H,f3)= MS,(G,I1) M ( 11 F)= (R 0 C , 1 where C , 0 = M(F,G,H3) O, i=l,2,3,... represents the continuity similarity equation (2.5.1 a), M the tangential momentum similarity equation (2.5. ib), M the radial momentum similarity equation 1 the axial momentum similarity equation (2.5. id). (2.5.lc), and M Therefore, the first order terms in the polynomial are the approximate solution of the governing equation with accuracy of O( U (i)) (order U ()). The residue for the second order solution is O(U()). To summarize, the quasi-similarity solution exists provided 1( U < 1. The more terms in the polynomial we obtain, the more accurate is the solution. 2.8 Conclusions A quasi-similarity method for 3-D vortices in a freestream is discussed in this chapter. The general solution in polynomial form has been found. The quasi-similarity solutions are divergent when z < 2 F 1 (or U 1 (z) 2 64it f U, v > 1). Thus, these solutions are most accurate far downstream. The quasi-similarity solutions have some interesting features. All three velocity components decrease with increasing downstream distance. The radial velocity diminishes mostly rapidly, followed by the axial velocity. The tangential velocity changes most slowly. For example, for the first order problem, the tangential velocity decays as one over the square root of the downstream distance, while the axial velocity is proportional to one over the downstream distance, and the radial velocity is proportional to one over the downstream distance to the power of 1.5. 39 Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices In the previous chapter, the quasi-similarity method was discussed. The solutions in polynomial form of the nonlinear governing equations have been presented. Each term in the polynomial consists of the product of two functions. One function is called the amplitude function, and is a function of downstream distance only. The other function, called the similarity function, is a function of a similarity variable that is a combination of the radial coordinate and the axial coordinate. The amplitude functions of each polynomial term have been obtained in the previous chapter. The ODEs and boundary conditions for the similarity functions in each order of the polynomial have also been listed. There is one constant, fi, involved in the first order problem, whereas two 1 and Cf. are involved in the higher order problems (i constants, C’h 2). In this chapter, we solve the first and second order problems numerically and analytically. It is found that the constant ,6 is closely related to the total wing drag. The two constants Ch and C’f (i 2) can be evaluated as well. The quasi-similarity solutions meet the requirements that the vortex circulation and wing total drag are independent of the downstream distance. Fortunately, first order velocity components and pressure in completed functional form can be obtained. The tangential velocity component and pressure in the second order problem can be presented in completed functional form as well. Second order axial velocity and radial velocity similarity functions can only be calculated numerically. The numerical procedure to be used for solving higher order problems will also be mentioned briefly. 3.1 The First Order Solutions 3.1.1 General first order solution in complete function form The first order problem consists of four ODEs for the similarity functions, along with the proper boundary conditions: (1)+2R 1 2iiR ( i1)= 12 .I (1)+21 1(11) (3.1.1) Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 1 211• H ()= 0 1 1 (1)+21 H +6) H (11)+(12 40 (3.1.2) (i1)=11•H(i1) 1 G (3.1.3) 211.J(1)+(1,2 (3.1.4) where +2).) 3 Oi)+ii .HO1)] 1 (ii)+2m1(ii)=—f3.{2ii.G is a constant. The boundary conditions of this first order problem are (O)=O 1 R (3.1.5) H(0)=LO, H(O)=O (3.1.6) F(O)=1.0, F’(O)=O (3.1.7) (oo)0 1 G (3.1.8) (3.1.5) must hold to ensure the flow is symmetrical about the vortex centerline. We apply (O) does not affect the final 1 the first equation of (3.1.6) because the choice of value of H results (one of the characteristics of the similarity method). The second equation of (3.1.6) ensures the second derivative at the origin has a finite value for an axisymmetric flow. The first equation of (3.1.7) holds because J3is a free constant, the second condition ensures the second derivative exists at the origin for a symmetrical flow. (3.1.8) arises from the freestream velocity condition. The following relations must hold if quasi-similarity solutions exist: ()= 1 U (3.1.9) 64it f z 3 (3.1.10) (3.1.11) 1 () 1 =83U 8itz (3.1.12) The ODEs (3.1.1 )-(3. 1.4) are decoupled. The non-singular solution of (3.1.2) with the boundary conditions (3.1.6) is 4( (rI)=—- 1—e 1 H .ii } Thus, the solution of (3.1.3) with boundary condition (3.1.8) is (3.1.13) Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 1 (i) G where Ei(n,x) 2 —8+ 16e —8e = 41 (3.1.14) $----.dt is the exponential integral function. Substituting (3.1.13) and (3.1.14) into (3.1.4), one finds that the solution of (3.1.4) with the boundary conditions (3.1.7) is (i) (3.1.15) =4e .ç{e The solution of (3.1.1) with the boundary condition (3.1.5) is (3.1.16) where F (i) is given by (3.1.15). 1 By combining the amplitude functions and similarity functions obtained above, one can get the first order solutions in completed function form. = z() = F vz 2 l6it •Ie4.J!e4EEi12,_2Ei12,Y_+1l.dy+±.e4 (3.1.17) 4) J 2 ) y 413 [ j i._L_.11e4__.11e4 2tr vz 27tT ) r(1) = U 2 32 (3.1.18) ) •s .{ie . (3.1.19) P(l) =16’ -e + Ei[l]_ Ei[1} (3.1.20) Equations (3.1.1 7)-(3. 1.19) show that the radial velocity in a tip vortex falls most rapidly (as z’ ) with downstream distance. The axial velocity changes more slowly with z 5 (roughly as z ), and the tangential velocity varies most slowly with downstream distance 1 (as z° ). 5 Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 42 It is interesting to note that this first order tangential velocity distribution is identical to that of the Lamb-Oseen vortex provided one replaces the time variable in the Lamb in the tip vortex problem. Oseen problem with its equivalent, The Runge-Kutta method was also used to solve the first order problem. The Runge Kutta results agree to within numerical error with the exact analytical results given by (3.l.13)-(3.1.16). The results corresponding to 3 = 0.361 (this is the case when the second derivative of the axial velocity component in the radial direction at the origin is zero) are shown in Figure 3.1. 7 (rj) i = 1 (Tj) is the similarity function for the tangential velocity H component in the tip vortex. 3.1.2 Tangential velocity and pressure Notice that the first order tangential velocity component and distribution are independent of the constant 1 First Order Similarity Functions pressure ST 1 (q) 1 T We note in addition that (3.1.21) (0)—2.773 1 G and this result is also independent of Furthermore, 7 (‘ii) = . 1 (ri), ‘qH the function describing the radial distribution of tangential velocity, is a maximum when = 2 Figure 3.1 The first order quasi-similarity functions. I = 0.361 2.25, and its maximum value is = 2.25) = (3.1.22) 1. 2763 Thus, to first order, the radius of a tip vortex (R ) is given by 1 1 =2.25W R I—, vUo or R. z 2.25 (3.1.23) —-= The similarity between the growth of a tip vortex and the growth of a laminar, zero pressure-gradient boundary layer (for which /x 5/.J) is apparent. Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 43 The first order vortex radius is independent of the constant 1. This independence implies that the radius of the wing tip vortex far downstream (where the first order solution is a good approximation) increases as the square root of the downstream distance. 3.1.3 Effects of the constant [3 on the axial similarity functions The first order axial velocity consists Effect of BETA on Axial Similarity Functions of an amplitude function multiplied by a similarity function. The amplitude 1.4 1.2 13= function is always positive if the constant 0.5 0.8 13 is positive. The similarity functions 0.7 0.36 1 0.1 0.6 0.4 corresponding to several positive values of the constant j3 are plotted in Figure 3.2. When [3 0.2 0 o equals 0.361, the second order derivative at the origin is zero. For 8 >0.361 the peak in the axial velocity ‘f so ‘ c ‘t- C’ - t - 00 V — TI Figure 3.2 Axial velocity excess in wing tip vortices. (The curves from top to bottom correspond with the values of J3 listed in the column to the right) distribution occurs away from the vortex axis. This finding is consistent with observations by Mason and Marchman (1972) of a real wing tip vortex. For 0<fl <0.36 1 there is still an axial velocity excess in the tip vortex, although the axial Effect of BETA on Axial Similarity Functions F() 13= -0.1 -0.361 0.8 0.6 0.4 0.2 0 -0.5 -0.7 00 — velocity decreases monotonically from the -0.4 -0.6 centerline. Green and Acosta (1991) have -0.8 — so — observed similar behaviour.* When * [3 is negative, the first order Figure 3.3 Axial velocity deficit in wing tip vortices. (The curves from top to bottom correspond with the values of/i listed in the column to the right) One must be careful with these comparisons. As we shall see, 3>O corresponds with negative wing drag, which was not the case for either of the two experiments mentioned. In addition, both experiments were carried out for small z, where this theory likely does not apply. Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 44 amplitude function is always negative. The first order similarity function is shown in Figure 3.3. Because the axial velocity is the product of these two functions, the behaviour shown in Figure 3.3 implies a centerline axial velocity deficit combined with a small core edge velocity excess. This kind of axial velocity distribution was reported by Logan (1971) and Thompson(1975). Baker (1974) and Green and Acosta (1991) also observed an axial velocity deficit by using LDV and PIV respectively. However, they did not report a core edge velocity excess, possibly because of the small downstream distance in their experiments, or their comparatively high (and therefore turbulent) vortex Reynolds number. 3.2 Drag due to a single quasi-similar vortex in viscous flow A single trailing vortex is considered in p this section. The axial momentum equation z will be applied to a control surface enclosing the wing tip region. The control surface has the form of a right cylinder with DRAG generator parallel to the z-axis, which marks the centerline of the wing tip vortex, and the area A at each end face (Figure Figure 3.4 Control volume around the wing tip. 3.4). The upstream end face and the curved surface of the cylinder are both at a large distance from the wing tip (compared with the vortex radius ), so that conditions there are approximately as in the freestream. Then, in the usual way, we find z-momentum flux outwards across curved surface z-momentum flux outwards across end face resultant normal force on end faces * = (p = — = pJ p) pU., — J(u u) . — . dA dA Tip vortex radius is just a few percent of the chord (Green and Acosta 1991, Fruman et a!. 1995), so it is possible to select a radius that is large relative to the vortex yet small relative to the wing. Batchelor (1964) used the same control volume. Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices where i 45 and p are, respectively, the z-component of the velocity and the pressure at the downstream end face, a distance z from the wing tip. One can verify that the ratio of the shearing force on the curved surface of the control volume over the force due to pressure gradient and momentum exchange is proportional to /R2. control surface is therefore small and can be neglected, when The viscous stress at the .& approaches infinity. The drag D on the wing due to this single vortex is thus given by D =J{p. —p+p .(u —i)}.dA (3.2.1) The integral of (3.2.1) can be represented as a polynomial by using the previously derived formulae for pressure and axial velocity components in the vortex. —p+pu, •(U—1) r 1 i-I (3.2.2) Substituting (3.2.2) into (3.2.1), the drag of this single vortex can be represented as D = -2p• u 1 (1) + 0i) + F 1 G (i) . F1iOi)] 1• di} or 2 D=-2pF .J[G.i+F+F1].d} 1 .4 {2] (3.2.3) Because the wing drag is independent of downstream distance D= urn f[Gi(ii)+ 32Jt1R- and At first glance, 1 G (rj) oc i •dii, , (3.2.3) implies for i =1 (3.2.4) fori2 (3.2.5) -4- at large radius, and therefore the integration in (3.2.4) 11 seems to be logarithmically divergent in the outer field (which is why the limit has been 46 Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices taken in that integration). However, we will verify that the integration is finite because ‘ —- 13 for large radii too, and fortuitously the 1/112 terms cancel. T Now let us consider the integration in the first order drag. Because urn J’ 13 ii’di,= urn o = = I [3 [ 8 limf{J i] .dY}.d[e4J+ 2,+ 1 1 2Ei R. di 413 $ .EEjI2,— 4) jj 2) urn j (3.2.6) ROOO11L Thus, urn [G (i)÷ 1 0 flR UR = urn $ I [3 l)]lldll r 2 2N / l—e I +4rI Eu l,— 1—Eu —i r [ t 2) ) 7 2\1 F 1 I+--I Eu 4,)J ru[ i,— / / 2 2,— 1—2EiI 2,— i+i j!>.dTj+ 4) Jj 2) 13 2 (3.2.7) 11 The integration on the right side of the above equation was proved to be zero using Maple V. a symbol manipulation software package. The drag due to a single quasi-similar wing tip vortex is therefore inversely proportional to the constant D=— lim J[G (ii)+ 1 32it L 0 l1idi,=— l6itf3 j 13 [3: (3.2.8) It has been assumed in the analysis of Chapter 2 that the distance between the two wing tip vortices is much larger than the vortex radius, and therefore that interference between the two vortices (e.g. the Crow instability) is negligible. The total drag on the wing is thus two times the drag due to a single vortex: (3.2.9a) 8t[3 Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices — Thus, = 47 2 pF (3.2.9b) 8tD In Figure 3.5 the non-dimensional Effects of Dntg on Axial Velocity axial velocity (defined as the axial 7 velocity divided by the term = l6it 2 V•Z is plotted against the non-dimensional drag (defined as 2it•D pF 2 ). The larger the • 0 ‘> -I 0 z Non-dm2nsional Drag= -‘-5 -0.69 -0.5 -0.36 0.36 0.5 0.69 2.5 --3 drag, the greater the deficit of the axial velocity on the vortex centerline. In a real trailing vortex an axial velocity is excess often observed just Figure 3.5 The effects of wing total drag on the axial velocity. (The different curves from top to bottom correspond with the different non-dimensional drags tabulated to the right) downstream of the trailing edge (Green and Acosta 1991). This velocity excess occurs as a result of the wing tip roll up. After the shed vorticity has been fully rolled-up into the wing tip vortices, usually far downstream of the wing, only an axial velocity deficit can exist for a wing experiencing a positive drag force. An axial velocity excess can exist far downstream of a wing generating thrust (e.g. with tip mounted engines). Figure the shows 3.6 nonEffects of Drag on Radial Velocity dimensional radial velocity (defined as Non-rjinnsionaI Drag— the radial velocity divided by the term F 2 .JiJc.z’ 2 16it 5 ) as a function of the non-dimensional drag. The amplitude of velocity radial rc ‘Uz(I) -2.5 -4169 -0.5 -0.36 1 () —), —..Ju F F oc is small and it decays —t 0.36 0.5 (169 2.5 Figure 3.6 The effects of wing total drag on the radial velocity. (The different curves from top to bottom correspond with the different non-dimensional drags tabulated to the right) more rapidly along z than does the axial velocity. Despite its small value, the radial velocity has important effect on the second order terms. As will be mentioned in section 3.3, the existence of a radial velocity leads Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices to one important conclusion of this chapter — 48 the second order tangential velocity must be zero. 3.3 The second order solutions Substitution of the first order solution (3.1.16) into ODEs (2.4.18a-d), which specify the second order problem, yields 2 F,’(rI)+4r. F,(rj) 2rR(ri)+2R,(rj) = rj (3.3.la) 1 6).H(rI)+4rI.H,(1)=O 2 21,.H’(11)+( + (3.3.lb) rI) 1)H 2 ( 1 G(T1)=2T.H ( (3.3.lc) l(1) )+(fl 2).F,’(11)+4f 2 21•(11 + where fi is (33 id) the constant that was established in the first order problem. The boundary conditions of this second order problem are (3.3.2a) O)=O (R 2 O) H ( 2 F, (0) = Ch., Cf, constant, H(O) constant, F,’(O) = = co)=0 G ( 2 = = 0 (3.3.2b) 0 (3.3.2c) (3.3.2d) (3.3.2a) must hold to ensure the flow is symmetrical about the vortex centerline. The constant in the first equation of (3.3.2b) is related to the value of second order vorticity at r=0. The second equation of (3.3.2b) ensures the second derivative at the origin has a finite value for symmetrical flow. The constant in the first equation of (3.3.2c) is related to the second order axial velocity on the vortex centerline, the second condition ensures the second derivative exists at the origin for a symmetrical flow. (3.3.2d) arises from the freestream velocity condition. The following relations have to be satisfied in order for the quasi-similarity solutions to exist, (3.3 .3a) Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 49 (3.3.3b) (3.3.3c) = 8iU() (3.3.3d) There are two constants, Ch 2 and Cf . in this second order problem. Two conditions 2 must be satisfied by the second order solution — the downstream independence of both the vortex circulation and the total drag of the wing. These conditions are represented mathematically by lim[ri2 (3.3.4a) .H(i)J=O (3.3.4b) $[2 3.3.1 Second order tangential velocity and pressure The non-singular solution of ODE (3.3. ib) with boundary condition (3.3.2b) is = Cli, e (3.3.5) The solution of (3.3. ic) therefore is G ( 2 i) = Ch .[4Ei(1,1)_4Ei(1,)j (3.3.6) The tangential velocity component and pressure distribution are therefore U, •Ffr( )(z,T1) = 2 = = thus, 1)e(,)(z,T1) (ij) 2 8iti3U(z) H pU .f3.U(z)G,(ri) (3.3.8) 3 • 2 Ch F = = 5 J 3 512ic . (vz) 3U Ch p 2 T 1024ic (vz) 2 (3.3.7) 11e .[Ei(1,c)... Ei(1,!)] (3.3.9) (3.3.10) Now let us consider the constraints (3.3.4a). Since lim[i2 . 112(11)] = Ck li(112 .e-5) = 0 (3.3.11) Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 50 . (3.3.4b) can be 2 (3.3.4a) is satisfied without placing any constraint on the constant Ch evaluated by integrating all the terms in (3.3. id), yielding F (1)+(12 (3.3.12) =J{.[41.G,(1)+21 .Hj(11).H,(1)J_2n12(1)}.d11 Because $o + 2). F,’(l) d = (12 + 2). F, (1: - 2$ F, (1) ° (3.3.13) ri).ridri =—2Cf,—2$F ( 2 $ 2iF, and (1) d = 2{1F’(1) . = -j F’(1) di} (3.3.14) 0 0 2 -2{F,O1)1} = 2Cf thus, (3.3.12) can be rewritten as (l)+[l2 2 j{1.F,(1+121}.d1 =_J{2l.G o o = .G(i)]}.d (3.3.15) —13$ G,(1)1•dfl That is 1)+ F(1)+I2(1)}.1.d1 = 0 ${13.G ( 2 (3.3.16) (3.3.16) is exactly (3.3.4b). The equivalence implies that the second order drag equals . 2 2 and Cf zero, independent of the constants Ch To determine these two constants, let us consider the physics behind the quasisimilarity method. The first order problem Consists of the most important terms. Because the first order solutions include information on the vortex circulation and drag due to the vortex, they construct the first order approximation to the solution. The second order problem can be considered to be an ODE system with the lower order results as inputs. Higher order solutions are used to compensate for the non-linear terms (or higher order Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 51 terms with respect to the downstream distance) in the governing equation. For the second order ODE system, complementary functions, which are the solutions of the homogeneous equation, exist. The non-singular complementary functions of (3.3. ib) and (3.3.ld) are Cli, e and Cf, .(rj2 —4).e respectively. These solutions are unrelated to the circulation of the vortex and the wing total drag, and are also not functions of the 2 and lower order solutions as well. Physically, it is reasonable to select the constants Ch 2 as zero so that only the part of the second order solution beyond the first order results Cf remains: 2 Ch = (3.3.17) Cf, =0 The second order tangential velocity and pressure are thus, = (3.3.18) 0 (3.3.19) Therefore the first order tangential velocity and pressure are actually the second order approximation of the solution. In other words, the first order solution are accurate up to the third order, O(U(z)). 3.3.2 Second order axial and radial velocity components The governing equations (3.3.1 a) and (3.3.1 d) can be rewritten as (3.3.20) () = i 2 ’(l)+ 411. F(1) 2 2flR(1)+ 2R 2F 2Th F (ii)+0i 2 +2). F’(1)+41 F,(i) = _2112(11) (3.3.2 1) The boundary conditions of this second order problem are (3.3.22) (0)=0 2 R (0)=Cf F = 2 0, F’(0)=0 (3.3.23) The axial and radial velocity components can be obtained numerically by the Runge 2 must be a 1 is a function of the constant J3 F Kutta method. Notice that because F function offlas well. Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 52 By using the Runge-Kutta method and the first order functions, one can obtain the second order functions for any value of ,& The numerical results for the axial similarity function are shown in Figures 3.7 and 3.8. For all values of /3, the second order axial similarity function is negative near 11=0, whereas the second order axial amplitude function is always positive. Therefore, the first order solution underpredicts the axial velocity defect in the vortex core when /3<0 (i.e. when the wing experiences a positive drag). For the case /3 >0 (a wing generating thrust), the first order solution overestimates the axial velocity excess in the vortex core. Second Order Similarity Function 0.05 Second Order Similarity Function 0.10 0.00 “ -0.20 0 . ‘.0 t 03 .—05 . -0.50 -0.60 -0.30 11 Figure 3.7 The second order axial quasisimilarity function. fi <0 11 Figure 3.8 The second order axial quasisimilarity function. /3 >0 The second order radial similarity function can be expressed in terms of the axial similarity function, 11 ( R ( 2 )+.JF ii)=flF,( x).xdx 2 110 (3.3.24) Numerical values of the radial similarity function for different values of 3 are plotted in Figures 3.9 and 3.10. Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 53 Second Order Similarity Function Second Order Similarity Function 0.00 0.00 -0.30 : -0.40 -1.50 -0.50 -2.00 -0.fO ‘ -0.20 11 T Figure 3.9 The second order radial quasisimilarity function. JJ <0 Figure 3.10 The second order radial quasisimilarity function. fi >0 3.4 The procedure for i-th order solution ( i> 2) The i-th order solution can be obtained by solving the following ODEs with the boundary conditions and constraints: (3.4.1) •’+2i•TI•F 2 1 =1 21R[+2R ’4-(i +6)H,’+2ilH, 1 2iiH 2 = {(2iH G= ( +4H)R _(2HF )} 1 +2(i—j)i (3.4.2) (3.4.3) . 2iF (12 + 2)’+ 2i1 ’ 1 = — .[2i1G + 1=’ - (3.4.4) .[2(i_j)F . where fi is the constant that was established in the first order problem solution. Like the second order problem, the boundary conditions of this i-th order problem are R(O)=O (3.4.5) 1 =0, H,(0)=ch (3.4.6) F(0) = = 0, F’(0) =0 (3.4.7) Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 54 (oo)=0 1 G (3.4.8) There are two conditions the i-th order solutions must satisfy — conservation of the vortex circulation and the independence of the wing total drag with the downstream distance. Respectively, and (3.4.9) 0 lim[12 H (ii)] (ii)+ (+(i)• F (i)}.di = 0 1 G (3.4.10) Meanwhile, the following relations have to be satisfied if the quasi-similarity solutions exist, (3.4.11) ()=U U ( 1 ) (3.4.12) I(z) =f3U() (3.4.13) ‘() =8icI3•U(() (3.4.14) The following procedure could be used to solve these coupled ordinary differential equations. Suppose that the lower order (less than i-th order) solutions have been obtained. Then the i-tb order problem can be calculated numerically by solving ODEs ,H 1 , and G 1 (3.4.1)-(3.4.4) with the boundary conditions (3.4.5)-(3.4.8) to find R,, F . 1 In summary, the i-th order results can be expressed as = = r(i) (z,r) . U = U (z) F (Ti) = (3.4.15) 8U VZ )(z,r) 1 rf( V 1 (i) U( (z) R p(l)(z,r)—pU,l3Ul(z)GI(T1) where (z)= 1 U 1 1-2 13 U 2 64t U(z)Ti(1) (3.4.16) (3.4.17) (3.4.18) 55 Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 3.4.1 Conservation of the vortex circulation F We have assumed in section 3.1 of Chapter 2 that, in order to conserve the vortex circulation, the following equation must hold, {2HO)] =0, (i2) (3.4.19) Now we will verify that this relation does hold for the higher order problem. According to the results of the first order problem, the following asymptotic characteristics can be found, (r) 1 H —-- (3.4.20a) 11 (rj) 1 G )—- (3.4.20b) 11 (il) 1 R ‘° -11 (3.4.20c) -11 (3.4.20d) Based on the ODEs governing the high order problems, one can prove by mathematical induction the proposition that the asymptotic results for the i-th order problem are 411 R(ii) (3.4.21b) (ii) 1 G I(i) (1) 1 R (3.4.21a) 11 °° —- (3.4.21c) 11 (3.4.21d) ‘ 11 We begin our proof by checking the first order problem (3.4.20). Obviously the proposition (3.4.2 1) is correct when i=1. Next, we assume the proposition is correct for i=l, lk. This assumption means that for the 1-th order problem, Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices H ( 1 ) 56 (3.4.22a) >—- TI G,(rl) >—TI (3.4.22b) —- (3.4.22c) °° F(ii) TI R ( 1 1) (3.4.22d) 1-I- Now, we prove the proposition is correct for i=k+1. We are going to use L’Hopital’s Rule in our proof. For an arbitrary function f(TI) if limf(TI) =0 and limf’(TI) = 0, then one has, lim ‘TI urn = — = °TV -‘-. urn f”i) 2m- TI provided these limits exist or are +00 or -oo (3.4.23) (Barnett and Ziegler 1991). The above relation can be rearranged as 2 f’(TI)=—IimTIf(TI) limTI (3.4.24) lim (3.4.25) f”() = 2 urn TI f(TI) Replacing i with k+1 in equation (3.4.2) and taking the limit on both sides, one can obtain the following equation using (3.4.22a-d): ] ] 1 + 2(k +1). lim[2TIH, + lim[(ri + 6)H÷ The value 6 is small compared with 1 is much smaller than TI H m Hk+I (3.4.26) n{TIHk÷l] 1 TI 772 when 77 approaches infinity, and based on (3.4.25) at large i Therefore, (3.4.26) can be rewritten as urn TIHk÷I (TI) or (3.4.27) Hk÷l(TI) Similarly, replacing i in (3.4.3) with (k+l), multiplying by taking the limit yields lirn[TI2 . 1 (TI)] G’÷ ii TI’ 172 on both sides, and 57 Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices or (‘n)] lim[i — Therefore, Gk+j (11) (3.4.28) > We can do the same thing to equation (3.4.4) as we did to equation (3.4.2): ] 1 n{iiF+ ijm[21,Fij+iim[(i12 +2)F]+2(k+ 1) ‘(k+l-j) (3.4.29) The first term on the left side of the above equation is much less than the third term. The value 2 in the second term can be neglected compared with 772 at large 17. Consequently, (3.4.28) can be rewritten as 1 1 lim rF÷ or 2k÷1 (3.4.30) ÷(11) ) (3.4.1), along with boundary condition (3.4.5), can be solved. The solution is (3.4.3 1) (11) = 12. ()+2(i— 1)$x(x).dx 1 211R Replacing i with k+ 1 in the above equation, and dividing by 277 on both sides, yields Rkl (1) = () — 2 +J x(x) dx (3.4.32) . 11Q Taking the limit on both sides of the above equation, one has lim Rk+I (i) = 1 (ii)] + k 1im[. F÷ 1im[’rF (3.4.33) (n)j Therefore, Rk÷l (1) (3.4.34) > 1 2(k÷I)-1 We have therefore proved that the proposition is correct for k+ 1, provided that it is correct for i—i, lk. Since the proposition is correct for = 1, this concludes the proof by mathematical induction. Hence* The second order tangential equation happened to be a homogeneous ODE. Therefore the second order similarity functions described by (3.3.5) and (3.3.6) are different from that described by (3.4.2 la-b). The Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices for 0, [2.H(1)} 58 i2 Thus, we have verified that (i2) =0 [2H(1)] 3.4.2 The i-th order drag In this sub-section, we will prove that the i-th order solutions satisfy the condition that the higher order drag is zero. Integrating the axial momentum equation from zero to infinity, yields j{2o2 +2)’+2i1} •dii .[2iiG +ii(H = . 0 0 F’ + . - - j) (3.4.35) d + j=I Because 52 () di = 2{’()I: = - S () (3.4.36) -2{(1)I} =2CJ and 5(2 + 2)’Oi).dii = (2 +2). - 2J(i).idii (3.4.37) =—2Cf —2$I(1).r1dT1 the left side of the equation (3.4.35) is thus, + 2)I’+2iiI} .di = 2(i — 1).f I(ii)• i,dii (3.4.38) Furthermore difference does not affect the validity of the assumption, it only hastens the approach of the similarity functions to zero as the similarity variable approaches infinity. 59 Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices = _{2(i_ 1)G, j=I 0 = - {s 0 For the higher order problem (i — 2), } + [2G 2(i -1). G 1 () [112G (11)] ].d + ) .J[2iiGi +ii(H 1 (3.4.39) = 0. Therefore —2(i — 1).JG(i).idii (3.4.40) Now, consider the last term in equation (3.4.35). We can replace R with F based on the continuity equation (3.4.1). (3.4.1) can be rewritten as 2[iR (i)] 1 or 1) 2iiR ( 1 = = [12F (11)] + 2(i - 1)i (11) (3.4.41) )+2(i— l)JxFj(x)dx 2 ( Thus, . 1=’ 0 = 0 { - 1)• f - 2(i i)1F . 0 } di (3.4.42) ={2i1.s .d_2(ij)fflFF .d} =—2(i— 1)Ji •d j=l 0 Substitution of (3.4.38), (3.4.40) and (3.4.42) into (3.4.35), yields 2(i_1){J[. 1 1). F()].d} = 0 ()÷(1)+F G ( For the higher order problem (i (3.4.43) 2) 0, i 2 (3.4.44) 60 Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices This result means that the higher order solutions satisfy the requirement that the wing total drag is independent of the downstream distance because the higher order drags are all zero. 3.5 The axial velocity on the vortex centerline The axial velocity on the vortex centerline is (z)]=U i_ 1 1)(z,O)=U .[l+U [ 1 D 8itpvzUJ (3.5.1) (O)=O for i2. The result given by (3.5.1) is valid provided 1 since F 1 (z)I 1U = D 8itpvzU,, <1 (3.5.2) where D is the wing total drag, U, is the freestream velocity, z is the downstream distance from the wing tip, p is the density, and v is the laminar kinematic viscosity. The axial velocity component on the tip vortex centerline can be evaluated by using (3.5.1) if the wing total drag is known. The axial velocity deficit on the vortex centerline is proportional to the wing total drag, and inversely proportional to the downstream distance. At a cursory glance it would appear that the axial velocity is not related to the circulation of the wing tip vortices. This impression is not correct. The total drag of the wing is the sum of the friction drag and the tip vortex induced drag, and because the induced drag is a function of the wing circulation, so too is the vortex axial veloctiy. 3.6 Comparison with Batchelor’s vortex Batchelor (1964) made a first effort to solve the wing tip vortex flow analytically. Batchelor’s analysis involves a linearization of the governing equations. Is there any connection between Batchelor’s vortex and the first order quasi-similarity solutions? As we shall see, the tangential velocity component and pressure distribution in Batchelor’s vortex are exactly the same as our first order tangential velocity and pressure, because the Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 61 ODEs and boundary conditions for the tangential velocity and pressure are the same for both models. In contrast, the axial velocity is different between the two models. The axial velocity field of Batchelor’s vortex results in a single vortex with an infinite drag. Batchelor resolves this problem by invoking the (somewhat artificial) concept of a “drag associated with the vortex core of a trailing vortex.” There is no equivalent limitation in the quasi-similar solutions; the axial velocity distribution does not lead to infinite values of drag. The limitation on our quasi-similarity method is U 1 (z)I < 1, which is a significant improvement over Batchelor’s vortex, for which the constraint is equivalent to (z)I <<1. 1 1U A summary of the procedure Batchelor used to obtain his similarity solution, along with commentary, follows. Far downstream, where the boundary-layer-type approximations (assumptions (3) and (4) in section 2 of Chapter 2) are supplemented by the approximation U, — (3.6.1) <<U the governing equations (2.2.5)-(2.2.8) reduce to 0 ia(!1)r)al)z r 3r az r (3.6.2a) 3z ar (3.6.2b) r3r (3.6.2c) par az az ar r ar) (3.6.2d) Batchelor solved the tangential and radial momentum equations by similarity methods. The tangential velocity and pressure obtained by him are found to be exactly the same as our first order solutions. Batchelor found 1)9 =_—(i_-e) (3.6.3) Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices - = = 2 pUT where the similarity variable = Ur 2 4vz {(i_e 62 2 . . = —, 4 (3.6.4) +2ei()_2ei(2)} and the function ei () = J dy. It is yj not difficult to show that (3.6.3) and (3.6.4) are the same as the first order quasi-similar solutions. The axial momentum equation (3.6.2d) now becomes u ._v.1à2z 3r iz E()+] r ar ) 32itvz [ d j (3.6.5) In order to obtain a clearer view of the asymptotic dependence of v on z, all the terms of above equation were integrated over a cross-sectional plane: d i (u — Z — ) rdr = F 2 U 2 l6it — since F 2 , l6ivU 7d(P) ‘ d d (366) z P— as—>co.Thus j(U, —ii).rdr= o F zU , log +const. l6rU,. v (3.6.7) where v/U, has been used as a convenient unit of length in the logarithm. The relation (3.6.7) suggests the following form for an asymptotic solution of the axial momentum equation: F gZUo 2 10 . 32it VZ V F .Q ( 2 )—Le 32it v 8vz (3.6.8) where L is a constant with the dimensions of area, and the last term, the complementary function, accounts for any initial velocity defect that may be independent of the circulation. The similarity functions Q 1 and Q2 can be obtained and expressed as ()=e 1 Q ()=e1[ 2 Q (3.6.9) +Pp().e].d (3.6.10) Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 63 The axial velocity given by (3.6.8) is obviously different from our first order axial velocity. Now, let us consider the drag on the wing: J - p = P — P p + 2icrdr — (3.6.11) J Batchelor found that the integral of (3.6.11) is logarithmically divergent as r—oo, and therefore that the drag associated with his isolated trailing vortex is infinite. Batchelor argued that the divergence does not dependent on the structure of the core of the trailing vortex, and on this basis he defined the “drag associated with the core of a trailing vortex” as = p limr2j R[ +i(U — p 1 u).rdrL.log RU 4it vj j (3.6.12) D differs from half the total drag on the wing by an amount that depends only on f and the distance between the two trailing vortices far downstream. In the case of two trailing vortices a distance s apart far downstream, it may be shown quite readily, from an evaluation of the kinetic energy of the motion in the lateral plane, that the total drag D is D W =2DC +2—log_- (3.6.13) v Substituting the pressure and axial velocity solutions ((3.6.3) and (3.6.4)) into (3.6.12), Batchelor found 8ic p (3.6.14) 2 where 2 is a positive number not far from unity. The solution (3.6.8) can now be rewritten as = U 8vz ir p 4it v *, 32itvz () 2 Q (3.6.15) If the trailing vortex system from a wing on which the total drag D consists of two vortices with centres a distance s apart, the axial velocity can be written in the further alternative form Chapter 3 = ut,, 64 Quasi-Similarity Solutions for Wing Tip Vortices _ -- 8v j,ic p v 4it v j F 32it vz (3.6.16) In summary, Batchelor’s linearized analysis yields the same tangential velocity and pressure distribution as our first order quasi-similarity solution. Batchelor’s vortex has an axial velocity distribution that results in a single vortex with an infinite drag, a problem that Batchelor resolves by defining the “drag associated with the vortex core.” Our axial velocity distribution is different from Batchelor’s and does not have infinite drag. Batchelor’s analysis is limited to restrictive limitation - UJtL <<1; by way of contrast, our analysis has the less (U 1 z)I <1. 3.7 Comparison with experimental data 3.7.1 Axial velocity decay on the vortex centerline The axial velocity defect in the Axial Velocity Decay in Tip Vortex cores of trailing vortices behind a lifting airfoil of rectangular planform was measured using a scanning laser 0.4 _. 0.12 _ ° Doppler •‘ I ‘ — velocimeter (Ciffone and Orloff, 20) 40) 600 .— 8(4) 10(4) 120) -0.04 1974). Data were obtained at several different angles of attack and downstream distances ranging from 30 to 1000 chord lengths. -o. / Lie 111 Figure 3.11 Axial velocity decay with downstream distance in a wing tip vortex. (The square symbols represent the experimental data, the solid line is the quasi-similar result and the dashed line is Batchelor’s result) The experiments were performed in a water tow-tank facility 61 m long, 2.44m wide, and 1.7m deep. The test Reynolds number based on the wing chord was nominally 2.43x 1O. The wing used has an NACA 0015 airfoil section, an aspect ratio of 5.33, and a span of 0.61m. The test was designed to obtain continuous data from the near field to the far field by the hydrogen bubble technique. Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 65 The axial velocity defect data at 4° angIe of attack are selected for comparison. Good agreement between the quasi-similarity theory and the experimental results occurs for zic > 300 (Figure 3.1l)* This agreement was achieved by using a laminar viscosity v and a . drag coefficient, CD = 0.0005. This drag coefficient is about twenty times smaller than a typical aircraft drag coefficient of CD = 0.01. There is a plausible explanation for the unreasonably low drag coefficient. The flow in a real wing tip vortex is probably turbulent. If one can model the effects of turbulence by means of an eddy viscosity VT that is constant spatially, then the quasi-similarity solution is still valid provided V is replaced by this way, but enforcing the constraint CD = VT. Modeling the flow of Figure 3.11 in 0.01, one determines VT 20v. Although no measurements of turbulent Reynolds stresses have ever been made in the far field of a tip vortex, a turbulent eddy viscosity twenty times the laminar viscosity is typical of many flows. Thus, although the quasi-similarity theory has effectively a free constant, which has been used here to match the theory with the experimental results, the constant (related to the drag coefficient) selected to give good agreement with experiment is reasonable. Batchelor’s theory gives much poorer agreement with the axial velocity results (Figure 3.11) at these small downstream distances owing to the logarithmic terms in (3.6.16). At yçjy far downstream distances Batchelor’s theory and the quasi-similar theory show the same dependence of axial velocity deficit on the downstream distance z because iim(logz)/z = liml/z. 3.7.2 Maximum tangential velocity decay A study of aircraft trailing vortex system, involving actual flight testing, was reported by McCormick et al.(1968). Detailed velocity measurements were made through the vortex immediately behind a test aircraft up to a distance of approximately 1000 chord lengths downstream of the aircraft. Flight testing was performed using a U.S. Army 0-1 aircraft and a Piper Cherokee. Measurement of the trailing vortices of each aircraft was * The disagreement for zJc < 300 may be caused by varying eddy viscosity in the tip vortex near field and the limitations of the theory itself. Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices 66 made using a vortimeter and a tuft grid. The vortimeter consisted of a vertical array of 30, 1-ft-long, l/4-in.-diam cylinders. Each cylinder was supported on strain-gaged flexures so that the aerodynamic drag on each cylinder could be measured. Hence, the velocity distribution through the vortex could be determined as the vortex moves through the array, by measuring the forces acting on each cylinder as a function of time. Velocity distributions through the Maximum Tangential Velocity Decay trailing vortices shed by the 0-1 and Cherokee aircraft at several flight speeds 2 1.5 were measured. The non-dimensional l . . data are compared with our analytical 0.5 results in Figure 3.12. For i c < 200 the 1 0 0 200 agreement is poor, likely because of turbulent flow in the tip vortex. j contrast, the agreement far downstream, zJc >200, is satisfactory. This agreement 400 600 800 1000 200 1400 71C(Ij Figure 3.12 Maximum tangential velocity decay with downstream distance in a wing tip vortex. (The round symbols represent the experimental data, the solid line is the analytical result) confirms that the maximum tangential velocity decays as one over the square root of the downstream distance — one of the results predicted by the quasi-similarity method. 3.8 Conclusions There are two major practical limitations on the quasi-similarity solution. One is the assumption that the axial velocity deficit or excess must not be significantly greater than the freestream velocity. Very near the wing this condition may not apply. The quasisimilarity method is limited as well by the assumption that the flow must be laminar. The limitation on the quasi-similarity solution from the convergence requirement (2.7.3) can be rewritten (assuming a rectangular wing) as Z>DC l6t (3.8.1) where CD is the drag coefficient, Re is the Reynolds number based on the chord, and S is the span of the wing. This distance is about a hundred meters downstream of the wing Chapter 3 Quasi-Similarity Solutions for Wing Tip Vortices trailing edge for a typical drag coefficient CD 67 0.01, Reynolds number Re 5x and a one meter span wing (for which the flow should be laminar). We can also make the following specific conclusions: 1. The quasi-similarity method represents an improvement over the solution of Batchelor. The limitation on the quasi-similarity method is — 1U(z)I <1, instead of u /U.. <<1 as required by Batchelor. Thus, the quasi-similarity method can be used in trailing vortices comparatively close to the wing. 2. The tangential velocity component and pressure in the Batchelor’s vortex are the same as the first order quasi-similarity solution. The axial velocity field of Batchelor’s vortex differs from the quasi-similarity vortex. Batchelor’s vortex has an infinite drag due to a single vortex whereas a quasi-similar vortex has a finite drag due to a single vortex. 3. The first order tangential velocity and pressure distribution have the third order accuracy, O(U(z)). That might explain why Batchelor’s theory is a good model of the tangential velocity even at fairly short distances downstream of a wing. 4. The maximum tangential velocity results have been compared with experimental data. The agreement is satisfactory far downstream. The agreement confirms that the tangential velocity decays as one over the square root of the downstream distance, as predicted by the theory. 5. The axial velocity predictions on the vortex centerline has been compared with measurements. The agreement is good far downstream, which confirms that the axial velocity deficit on the vortex centerline decays as one over the downstream distance. 6. There are two main restrictions on the quasi-similarity method. One is the laminar flow assumption. The other is that quasi-similar solutions are only valid at fairly large downstream distances. Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 68 Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing A novel hydrofoil design, consisting of a small diameter flow-through duct affixed to the tip, has been studied. The tip vortex cavitation inception index, cy., of this hydrofoil geometry is about a factor of 2 lower than that of a conventional rounded hydrofoil tip. This inception improvement comes with little associated performance penalty. For angles of attack greater than 8° the non-cavitating lift-drag ratio is actually superior to that of an unducted hydrofoil of equal span, although with lower wing loadings the hydrofoil performance is diminished by application of the ducted-tip. The ducted tip is effective at reducing the tip vortex inception index because, in contrast with the rounded tip for which shed vorticity in the transverse plane behind the wing (Trefftz plane) is confined to a line, the ducted tip shed vorticity at the trailing edge is distributed over a line and circle. Distributing the vorticity in this fashion causes the trailing vortex to roll up less tightly, and hence have a higher core pressure and lower a , 1 than a conventional hydrofoil tip. It is also suspected that the interaction at the microscale level between the flow through the duct, and the flow around it, makes the vortex core size larger, and therefore a 1 smaller. The superior lift performance of the ducted-tip wing at elevated angles of attack results from the redistribution of the shed vorticity as well. The downwash caused by the vorticity on one portion of the ring tip in the Trefftz plane is smaller than that by the same amount of vorticity positioned in the spanwise direction. Less downwash produces less induced drag. At small angles of attack the benefit from the ducted tip device of reducing the induced drag is more than offset by the increased parasite drag due its larger wetted area. The ducted tip design has many potential marine applications, including to ship and submarine propellers, submarine control fins, and ship rudders. The following section reviews the previous work on tip modifications. The second section explains the rationale for the novel ducted tip geometry. The third section describes the experimental apparatus and procedures. Flow visualization, cavitation characteristics and aerodynamic performance of the ducted tip configuration are presented Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 69 and discussed in section four. The experimental uncertainties are estimated in section five. Finally, conclusions are presented in section six. 4.1 Previous Work on Tip Modifications The purposes of tip modification can be grouped into three categories: attenuating vortex strength to alleviate the rolling moment on a following aircraft, reducing the vortex induced drag, and increasing the vortex core pressure to avoid cavitation. Various kinds of tip devices have been tested to reduce the rolling moment. A leading edge disk flow spoiler, a trailing edge disk flow spoiler and a porous wing-span extension were investigated by Scheiman and Shivers (1971). They achieved vortex cross-sectional redistribution at the expense of wing lift and/or drag characteristics. “Vortex-attenuating splines” were tested by Hastings et al. (1975) and found to cause a considerable noise increase and a lift/drag ratio reduction. Snedeker (1972) measured the rolling moment on a simulated following aircraft caused by a tip vortex with and without axial injection (injection of fluid in the chordwise direction from the wing tip). He found that the injection reduced the rolling moment by only 13% (for vortex generating wings at a 6°). The applicability of such a scheme to commercial aircraft is dubious. The Whitcomb wing tip (Whitcomb 1976 and Flechner et al. 1976) is the only wing tip modification device that has found commercial application. This wing tip consists of a short (about 1/2 chord high) lifting surface mounted almost normal to the wing at the tip. The wing is connected to the suction surface of the wing with a smooth fillet. Some of the aircraft that have flown with such a winglet are the MD-il, the Gulfstream III, the DC-1O and the B 767 (Webber and Dansby 1983 and Devoss 1986). The Whitcomb wing tip reduces marginally, by -1% to 2%, the drag on some aircraft (Shevell 1989). Closely related in concept to the Whitcomb winglet is a series of winglets set at varying angles to the planform plane, which are referred to as “wing tip sails.” Spiliman (1978) has reported up to a 29% reduction in lift dependent drag in flight tests of (small aspect ratio) wings fitted with sails. Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 70 Platzer and Souders (1979) is a good, though now somewhat dated, review of tip vortex cavitation knowledge. Cavitation inception is known to occur in the region immediately downstream, to as far as two chords downstream, of the generating wing (Arndt et al. 1991, Maines and Arndt 1993, Stinebring et al. 1991, Green 1991). Inception is generally highly unsteady, occurring first at one location and subsequently up or downstream from it. Upon further reducing the flow cavitation number below g, a long portion of the tip vortex cavitates. These observations are consistent with the known minimum of vortex core pressure near the hydrofoil trailing edge (owing to rapid rollup), and with the small axial pressure gradient along a tip vortex centerline (resulting from the small variation in i with downstream distance). Tip cavitation inception is highly dependent on small details of the flow near the tip. Stinebring et al. (1991) and McCormick (1962) have shown that merely roughening the pressure surface tip region of a hydrofoil substantially reduces c. Drilling small holes in a hydrofoil from the pressure to the suction surfaces at the tip also reduces a (Arakeri et al. 1985, Sharma et al. 1990). Recently, Chahine et al. (1993) have shown that injection of modest amounts of Polyox solution (a viscoelastic polymer) into a propeller tip vortex can, under optimum conditions, reduce the cavitation inception index by up to 35%. They attribute this reduction to a significant thickening of the vortex core caused by the viscoelasticity of the solution. The dependence of tip vortex a on small details of the tip flow is congruent with the known sensitivity of the single phase tip flow to the same details. The size and quantity of freestream nuclei also has a substantial impact on o. Arndt and Keller (1992) have measured a doubling of c 1 when the freestream fluid is changed from ‘strong’ water (few and small nuclei) to ‘weak’ water (many, large nuclei). 4.2 The ducted tip wing geometry This chapter documents one successful tip cavitation inception delay device — the ducted tip. Before discussing our experimental studies of the ducted tip geometry, we begin with a brief explanation of the rationale for this novel geometry. Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 71 As a consequence of Helmholtz vortex laws, the total shed circulation from a wing is essentially established once the wing lift is specified. Therefore, the only way to modify the tip vortex core pressure (and thus the cavitation inception index, c ) is to redistribute this shed circulation. In particular, the less concentrated the shed circulation from a wing, the larger the tip vortex core will be, and hence the lower its inception index. A conventional planar wing sheds (in theory, see e.g. Mime-Thomson (1958)) a line of vorticity in the Trefftz plane (Figure 4.1(a)). Any wing with an elliptical or near elliptical loading (a desirable characteristic for maximizing L/D) has this shed Wing Root Wing Tip vorticity + concentrated at the wing tips, (a) which accounts in part for the small vortex core radius and hence high y, Wing Root + + of elliptically loaded wings*. Non-planar wings (Cone 1963) are not subject to the same near- (b) + + Figure 4.1 Shed vorticily in the Trefftz plane from (a) conventional and (b) ducted tip wings elliptical loading constraint as conventional wings, and hence do not necessarily shed concentrated vorticity at the wing tip. They therefore potentially have much lower tip vortex cavitation inception indices and, according to Cone, less induced drag, than conventional wings. The ducted-tip wing is one example of a non-planar wing. It consists of a conventional planar wing to which is mounted a hollow duct, approximately aligned with the wing chord (Figure 4.2). Figure 4.1 (b) is a Trefftz plane schematic of the configuration. Provided that the flow through the duct has little rotation (a reasonable assumption if the entrance to the duct is near the wing leading edge), and the flow around the exterior of the duct has a significant swirl (again, a reasonable assumption because the * Real fluid effects that cause the formation of secondary vortices have only a slight impact (addition of a small amount of countersign vorticity near the tips, and smearing of the distribution) on the shed vorticity distribution. Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 72 large pressure gradient that exists near the wing tip, from the pressure surface to the suction surface, will cause a rapid tangential flow around the tip), there will be significant vorticity shed from the duct, as illustrated in Figure 4.1(b). Because this shed vorticity is spread over a region of radius comparable to the (large) duct radius, it seems plausible that the resulting tip vortex core (a) (b) radius Figure 4.2 Schematics of wing tip devices. Flow is right to left. (a) Ring wing tip. (b) Bi-wing tip. the distribution of shed vorticity encouraged us to explore experimentally the ducted tip will also be large. This insight into geometry. 4.3 Experimental Apparatus and Procedures The line-of-reasoning described in section 4.2 leads one to hope that the ducted tip wing geometry has superior tip cavitation characteristics. From a practical standpoint, one is not merely interested in improvements in tip a., one must also know that these improvements are not coincident with a large lift reduction or drag increase. Consequently, the experimental program described here consisted of two discrete parts — measurement of the inception characteristics of a ducted tip wing, and measurement of the lift and drag behaviour of such a wing. The tip vortex inception measurements were carried out by my supervisor, S. I. Green, in the Low Turbulence Water Tunnel (LTWT) at Caltech. The water tunnel facility has a test section of 30.5 cm x 30.5 cm x 2.5 m long, with a freestream velocity adjustable up to 10 mIs. The freestream pressure, p, can be varied from slightly above atmospheric pressure to 25 kPa by means of a vacuum pump. By using a deaeration system and diatomaceous earth filtration, the freestream nuclei content of the water was controlled. A van Slyke device was used to measure the water’s total dissolved air content. Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 73 A rectangular planform, untwisted, constant airfoil section (NACA 64-309 modified), hydrofoil was reflection-plane mounted in the LTWT. The chord of this hydrofoil was 15.2 cm, and the semispan was 17.8 cm, resulting in an effective aspect ratio of 2.3. The Reynolds number based on chord length for all the cavitation tests was in the range 1.1 x <Re< 1.6x 106. 6 10 This basic wing was fitted with two different tips. The first tip was a duct comprised of a 2.9 cm outside diameter (with a wall thickness of 0.2 cm) cambered brass pipe, 10.2 cm long, attached flush with the hydrofoil trailing edge, with axis aligned with the camberline (Figure 4.2(a)). The size of the duct selected for these studies is probably not optimum. A duct was designed based on intuition guided by the following considerations: too large a duct will cause excessive parasite drag, while a duct that is too small will have reduced flow through the duct and will cause little change to the shed vorticity distribution. The second tip, which was used for comparison purposes, was a rounded tip geometry with approximately semicircular cross-sections perpendicular to the chordwise direction, fitted onto the end of the basic wing (Figure 4.7). Two different procedures were used to study the flow around these wing tip geometries. Surface flow visualization was accomplished using the paint drop technique (Green et al. 1988). With the hydrofoil out of the water tunnel, drops of oil-based paint were dotted on its surface. The hydrofoil was then returned to the tunnel, and the LTWT was quickly accelerated up to a set velocity, causing the paint drops to be smeared in the direction of the local shearing stress. The second procedure involved cavitation inception measurement and cavitation photography. The large range of inception numbers encountered precluded us from maintaining a constant Reynolds number for all the tests. Instead, cavitation inception testing consisted of setting the hydrofoil angle of attack and gradually increasing U and decreasing p, until cavitation inception was observed under stroboscopic illumination. Hydrofoil leading edge inception was defined to occur when cavitation was observed on the foil surface half way between the wing root and tip; trailing vortex inception when at least one cavitation event per second was seen. When the leading edge inception index, 1 )ie was larger than the tip vortex inception index, (a, ),. it was necessary to achieve (a Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 74 1 ), quickly after (a, )k to avoid erroneous measurements due to the recirculation of (a cavitation nuclei in the LTWT. Photographs of cavitation were taken by illuminating the hydrofoil and tip vortex in the spanwise direction, and recording on film the light reflected normal to the hydrofoil planform. The cavitation behaviour of the rounded tip was studied first, and several weeks later, following attempts with different ducted tip geometries, testing of the final ducted tip geometry began. The aerodynamic performance of the two tip geometries was evaluated in the Low Speed Wind Tunnel (LSWT) at the University of British Columbia (Figure 4.3). The truncated rectangular cross-section test section of this facility has dimensions 69 cm x 91 cm and is 4 m long. The maximum velocity attainable in the test section is 40 m/s, although the velocity was maintained at 30 m/s for these experiments. A six-axis straingauged force balance, combined with filters, amplifiers and a data acquisition sub-system (Figure 4.4), measures lift and drag forces on models in the test section. The sub-system consists of a Computerscope model ISC-16 for the IBM-PC family of computers, provided on a single plug in card. The ISC-16 model has a 12 bit 1 llsec A/D conversion, accepts analog input signals between -10 Volt and +10 Volt, and has a 1 to 16 channel programmable input multiplexer. About 10,000 samples were taken and averaged for one particular lift or drag measurement. This data acquisition system was calibrated and the results are summarized by the following relations: Drag Force [N] = 0.579472 x [Volt] + 0.029056 Lift Force [N] = 10.07694 x [Volt] + 0.26435 where [Volt] is the corresponding digital reading in voltage from the A/D converter. The data acquisition sub-system, together with a scanivalve/pressure transducer sub system, which consists of four pressure transducers, a solenoid driver and a four-cap scanivalve (Figure 4.5) was used for pressure distribution measurements. Each pressure sensor converts the pressure difference (in the range -5 inches to +5 inches 2 H 0 ) between its two input connectors to an analog signal of between 1 and 6 Volts. The calibration results of the four pressure sensors are respectively, Pt — J= 491. 6743 x[Volt]— 1730.78 2 Pref[N / m Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing P2 — 75 ] = 500.3630 x[Volt]— 1738.41 2 Pref[N / m 2 P 3 P 1 =484.0616><[’101tb1730.47 rej[N/m 4 P ]= 495.5708 x[VoltJ— 1745.80 2 Pref[N / m where [Volt] is again the digital output reading in voltage from the A/D convertor. Aerodynamic testing in the LSWT was done on a rectangular planform, untwisted, constant airfoil section (NACA 66-209) basic wing. This acrylic wing is 30.5 cm in chord and 35.6 cm in semi-span, and is fitted with 92 surface pressure taps. Several different tips were attached to the basic wing, including a series of cambered circular aluminum ducts 3.8 cm in outside diameter (with a wall thickness of 0.2 cm) with lengths ranging from 19.8 cm to 30.5 cm (i.e. 65% of wing chord to 100% of wing chord). These ducts were aligned with the wing trailing edge and were affixed approximately parallel to the wing camberline (Figure 4.2(a)). An elliptical wing tip duct, with major axis 10.5 cm (aligned normal to the wing planform), minor axis 3.9 cm, and length 19.8 cm was also tested. The last type of ducted tip tested was a Thi-wing” configuration. The “bi-wing” tip (Figure 4.2(b)) was constructed of two parallel NACA 0006 short airfoils of 15 cm chord and 3.5 cm span joined together at both span ends with thin aluminum plates. The chordline-to-chordline separation of these airfoils was 3.5 cm. This bi-wing tip was attached to the tip of the basic wing with the short airfoils aligned parallel, and in subsequent tests ±10° from parallel, with the basic wing. All plate/wing junctions were then faired with putty before aerodynamic testing. For comparison purposes, two other tip geometries were also wind tunnel-tested. A semicircular tip constructed of acrylic and body filler 3.8 cm in span was one conventional geometry tested, and a square cut tip (i.e. constant airfoil section to the tip, where the wing chord falls abruptly to zero) 3.8 cm in span was also tested. In summary, a number of different wing tip goemetries were tested in the wind tunnel, all with semi-spans of 39.4 cm, though with aspect ratios varying from 2.58 for the square cut tip, to 2.60 for the rounded tip, to 2.67 for the 65% partial chord ducted tip. The aerodynamic test configurations were not geometrically similar with the cavitation test models, but they were nearly so. Chapter 4 Cavitation and Aerodynamic Peiformance of the Ducted-Tip Wing 76 Three aerodynamic measurements were made in the LSWT: the lift coefficient, CL, versus angle of attack a, the drag coefficient, C , versus a, and the pressure distribution 0 on both sides of the basic wing. Flow visualization of the flow field around the ducted tip in the LSWT was carried out by using tufts on the model. Figure 4.3 The Low Speed Wind Tunnel at the University of British Columbia. E Cl) Cl) Amplifier 60,000 Gain CtJ 1) 0 I. Cl) Amplifier 3,000 Gain —. —‘-——-r - . Drag Force Signal Lift Force Signal Figure 4.4 The data acquisition system for lift and drag measurement in the LSWT. Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing Data Acquisition Sub-system Scanivalve/Pressure Transducer Sub-system Figure 4.5 The data acquisition system for pressure measurements in the LSWT 77 Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 78 4.4 Results And Discussion This section includes the cavitation results of the ducted tip, which were part of the work done by my supervisor at Caltech. A completed picture on the novel ducted tip will be presented in this manner. 4.4.1 Flow Visualization Examining the surface flow over an airfoil is perhaps the easiest way to assess qualitatively the performance of different wing tip geometries. Figure 4.7(a) is a view of the surface flow over the pressure side of the rounded tip geometry. The flow direction is nearly streamwise near the reflection-plane mount (near the bottom of Tip Vortex the photograph) and acquires an increasing spanwise component as the tip is approached. Over the majority of the suction side of the wing (Figure 4.7(b)) there is also an increasing spanwise velocity component (for this side, directed Figure 4.6 Tip vortex rollup in the Trefftz plane. from the tip to the root) as the tip is approached. Streaklines directed toward the tip on the suction side result from, respectively, the tip vortex rollup in the region near the tip (refer to Figure 4.6), and a secondary vortex near the wing root. Note that the reflection plane mount seems to affect the wing only over the bottom 20% of the semi-span. The tip flow in particular is basically unaffected by the presence of the water tunnel floor. Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing (a) (b) Figure 4.7 Surface flow visualization on the rounded tip geometry. (a) Pressure side. Flow is right to left. (b) Suction side. Flow is left to right. (Photograph source: Green 1988) 79 Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing (c) Figure 4.7 Surface flow visualization on the rounded tip geometry. (c) Inboard view of tip. Flow is left to right. The wing pressure surface is at the top. (Photograph source: Green 1988) 80 Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 81 The surface flow over the ducted tip geometry is quite different from that over the conventional tip geometry. The spanwise velocity component at the trailing edge, on both the suction and pressure surfaces (Figure 4.8), is substantially less than that of the rounded tip wing. At one third of the spanwise distance from the reflection plane mount, for example, the angle the rounded tip streaklines make with the freestream direction is about 50% greater than the strealdine angles of the ducted tip. This difference in spanwise velocity component implies that the ducted tip wing sheds less bound circulation over the majority of the wing span than does the rounded wing tip. Smear lines in the duct (Figure 4.8(d)) have little swirl, whereas those on the duct exterior indicate a large tangential velocity. This observation confirms that substantial vorticity is shed from the duct, as indicated in Figure 4.1(b). One may be concerned that the flow will separate from the ducted tip. Vibration of tufts located on the outside suction side close to the entrance of the duct implies that local flow separation occurs there at elevated angles of the attack (Figure 4.9). This flow separation will cause elevated parasite drag. Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing (a) (b) Figure 4.8 (a) Surface flow visualization on the pressure side of the ducted tip wing. The flow is right to left. oL=7°, Re=1.2x10 . 6 (b) Suction side view of the wing in Figure 4.8(a). Flow is left to right. (Photograph source: Green 1988) 82 Chapter 4 Cavitation and Aerodynamic Peiformance of the Ducted-Tip Wing (c) (d) Figure 4.8 (c) Inboard view towards tip of the flow in Figure 4.8(a). The flow is left to right. The pressure surface is at the top of the photograph. (d) Upstream view through the duct of the flow in Figure 4.8(a). The pressure surface is to the right. (Photograph source: Green 1988) 83 Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing Figure 4.9 Separatedflow near the entrance of the ducted tip. c=]2°. The flow is right to left. 84 Chapter 4 Cavitation and Aerodynamic Peiformance of the Ducted-Tip Wing 85 4.4.2 Cavitation Characteristics Figure 4.10 is a photograph of developed tip vortex cavitation behind the rounded tip geometry. The cavity cross-section is oval immediately downstream of the tip, and becomes nearly circular within one chord of the trailing edge. This shape is consistent with the known transition of the tip vortex velocity field from asymmetric at the hydrofoil trailing edge to nearly axisymmetric a short distance downstream (Fruman et al. 1991, Green and Acosta 1991). Development of this continuous vapour-filled tip vortex cavity occurs when the cavitation number is decreased by just 10-20% below y,. The rapidity of the cavitating vortex development implies that the axial pressure gradient along the core of a tip vortex is small, in agreement with recent measurements by Green (1991). Unlike a developed trailing vortex cavity, which is a steady phenomenon, cavitation inception is highly unsteady. Cavitation inception in the rounded tip vortex is observed 0.3 to 1.5 chords downstream of the hydrofoil trailing edge; the exact location of trailing vortex inception fluctuates in an apparently random manner within this range. It is not known if the fluctuation of the inception location is due to variability in the location at which freestream nuclei are captured by the vortex (Ligneul and Latorre 1993), or due to the fluctuating pressure field within the tip vortex core (Green 1991). Trailing vortex inception occurs for o of attack, a = = 1.8 ± 0.2 near the hydrofoil operating angle 7° (Figure 4.11). As a is varied about the operating angle of attack, the inception index varies much less rapidly than the (ct—c 2 dependence predicted by a ) 0 simple Rankine vortex model. Arndt and Keller (1992) have observed the same effect, which they attribute to the significant tension that can be sustained by low-nuclei-content water prior to cavitation inception. Hydrofoil leading edge surface cavitation inception occurs for cy = 0.9 ± 0.1, substantially lower than the value for tip vortex inception. This observation implies that tip vortex cavitation occurs in many situations for which surface cavitation on the hydrofoil is not present. The inception indices plotted in Figure 4.11 were determined in ‘strong’ (i.e. low nuclei content) water. Water with more freestream nuclei exhibited much higher Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 86 inception indices (greater by a factor of 2 when the dissolved air content of the water was doubled from 4 to 8 ppm). Arndt and Keller (1992) and Arndt et al. (1991) observed a similarly strong dependence àf a on nuclei content. Cavitation inception from the hydrofoil ducted tip was first observed on the outside edge of the suction surface portion of the duct. This surface cavitation is caused by the separation (Figure 4.9) that occurs as the flow, which has a high tangential velocity around the tip, encounters the duct leading edge. Figure 4.12 is a photograph of cavitation nuclei generated by the duct surface that have migrated into the ducted tip vortex. Note the diffuse appearance of the vortex. True trailing vortex cavitation occurs only when o is reduced well below the inception index of both the ducted tip surface cavitation and the hydrofoil leading edge cavitation. As a consequence, cavitation nuclei generated by both the wing leading edge and duct surface cavitation are swept into the tip vortex prior to inception there. Because the tip vortex is exposed to water of much higher nuclei content than the freestream water, the tip vortex will cavitate at relatively higher . The ducted tip vortex inception measurements are thus biased to high values. Despite this unfavourableB bias, the ducted tip vortex data (Figure 4.13) show a remarkable decrease in c5 relative to the rounded tip geometry. For an angle of attack of 7°, the inception index of the ducted tip vortex is, angle of attack, the rounded tip vortex has, = = 0.9 ± 0.1. In contrast, at the same 1.8 ± 0.2 (i.e. a of the ducted tip is 50% less). For angles of attack of 5° and 10°, the ducted tip vortex a. is 0.5 ± 0.05 and 1.8 ± 0.2. These values are respectively 61% and 28% less than the rounded tip a values (1.3 ± 0.1 and 2.5 ± 0.1) at the same a. As explained previously, the inception index of tip vortex cavitation is a strong function of the nuclei content of the water. One might therefore posit that the large reduction of a,, when the ducted tip is installed, results from a decrease in the nuclei content of the water between the time of the rounded tip tests and the ducted tip tests. In fact, two observations point to an increase in freestream nuclei between the tests. One such observation is the increase in the dissolved air content (from 5 ppm to 7 ppm with an Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 87 error of ± 0.5 ppm) between the rounded and ducted tip tests. The observed increase in hydrofoil leading edge a 1 (from 0.3 to 0.8 at a 3.0 to 3.5 at a = 5°, from 1.0 to 1.5 at a 7°, and from 10°) between the rounded and ducted tip tests is another strong indicator of an increase in freestream nuclei concentrations. In summary, two factors tend to bias the ducted tip a. to elevated values — the presence of more naturally-occurring freestream nuclei during testing of that geometry, and the generation of nuclei by cavitation from the hydrofoil leading edge. Despite these two influences, cavitation inception from the ducted tip geometry occurred at much lower inception indices than from the rounded tip geometry. One may be concerned that the gains in cavitation inception improvement resulting from the ducted tip modification are negated by the increase in other forms of cavitation. For example, one might postulate that the ducted tip geometry would redistribute the hydrofoil loading, causing the wing leading edge a. increase referred to above. This is not the case. Except for a small region near the wing tip, the measured pressure distributions on the basic wing were virtually unchanged (refer to pressure distribution measurements shown in section 4.4.3; c the same to ±7%) with the addition of the ducted tip. The unchanged pressure distribution implies an unchanged wing loading. Therefore, leading edge cavitation is not significantly affected by the addition of the ducted tip, whereas tip vortex cavitation is significantly reduced by the ducted tip. Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing Figure 4.10 Developed trailing vortex cavitation behind the rounded tip geometry. , a=1.8, DAC=5.4 ppm. 6 a=7°, Re=1.3x10 (Photograph source: Green 1988) 88 Chapter 4 Cavitation and Aerodynamic Peiformance of the Ducted-Tip Wing 89 CAVITATION INCEPTION INDEX (ROUND TIP) -10 0 10 20 ATTACK ANGLES [DEGREES] • Leading Edge • Tip Vortex Figure 4.11 Cavitation inception index versus angle of attackfor the rounded tip geometry. Re 1.4x10 , DAC 6 = 5.4 ppm (uncertainty in a the 95% confidence level) = .i 0.2° and in a, = 10% at Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing Figure 4.12 Ducted tip vortex made visible by migration of cavitation nuclei, generated at the duct leading edge, into the vortex. , y=].5, DAC=7.O ppm. 6 cx=7°, Re=1.2x10 (Photograph source: Green 1988) 90 Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 91 CAVITATION INCEPTION INDEX (RING WING TIP) -10 -5 0 5 10 15 ATTACK ANGLES [DEGREES] Tip Voex A Lding Edge ading Edge of Wing Figure 4.13 Cavitation inception index versus angle of attackfor the ducted tip geometry. Re 1.4 x ]Q6, DAC = 7.0 ppm. (Uncertainty in a at the 95% confidence level). = .t 0.2° and in a. = 10% Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 92 4.4.3 Aerodynamic Performance These encouraging cavitation inception results prompted a study of the aerodynamic performance of the ducted tip and other tip geometries. Wind tunnel testing demonstrated that the aerodynamic performance of the square cut and rounded tip wings was identical to within the experimental error (±1% drag and ± 0.5% in lift). The partial chord circular duct tip had the best performance of any of the ducted tip geometries. Consequently, the aerodynamic characteristics of the partial chord ducted tip and the rounded tip was studied in section 4.4.2 — — the same tip geometries whose cavitation behaviour will be emphasized here. 4.4.3.1 The partial chord ducted tip For all elevated angles of attack (a > 100) the lift coefficient of the ducted tip geometry (where CL is based on the planform area of the rounded tip) is about 4% less than that of the rounded tip (Table 4.1). However, the partial chord ducted tip geometry has a 3.4% smaller planform area than the rounded tip, and thus the actual lift coefficients for the two geometries are identical to within experimental error. Surface flow visualizations described in section 4.4.1 suggest that the ducted tip generates greater lift near the wing root. The implication of the congruity of the two lift coefficients is that the additional lift gained from the root of the ducted tip wing is almost precisely offset by the reduced lift attained in the vicinity of the tip. Table 4.1: Lift Coefficients of the Ducted-tip Wing and the Conventional Wing. a=6° 8° 10° 11° 12° 14° 15° Conventional Tip Ducted Tip 0.3797 0.4864 0.5946 0.6272 0.6614 0.7158 0.7484 0.7864 0.3730 0.4828 0.5703 0.6014 0.6309 0.6832 0.718 1 0.757 1 Improvement -2% -0.7% -4% -4% -4.6% -4.6% -4% -4% 17° Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 93 Unlike the virtually unchanged wing lift, a significant difference in wing drag is observed when the ducted tip is installed (Figure 4.14). For a < 8° the ducted tip geometry has greater drag than the rounded tip, but for larger a the ducted tip geometry has less drag: 6% less at 100, 8% less at 12°, and fully 10% less at 14° (Table 4.2). Even when allowance is made for the slightly reduced planform area of the ducted tip geometry, there is nonetheless an up to 6% drag coefficient benefit arising from the tip modification. This reduction in drag coefficient is somewhat surprising in view of the obviously increased parasite drag of the ducted wing (arising from both its additional wetted area and the likely existence of interference vortices at the duct-wing junction); the only explanation is that the duct attachment to the wing must reduce substantially its induced drag. This conclusion is borne out by the observation that the drag reduction of the ducted tip is most pronounced at elevated angles of attack, for which the induced drag is a larger fraction of the overall drag. Table 4.2: Drag Coefficients of the Ducted-tip Wing and the Conventional Wing. c=6° 8° 10° 110 12° Conventional Tip DuctedTip 0.0107 0.0183 0.0326 0.0488 0.0680 0.0110 0.0183 0.0307 0.0449 Improvement -3% +0.0% +6% +8% 14° 15° 170 0.1135 0.1418 0.1932 0.0623 0.1018 0.1279 0.1770 +8.3% +10.3% +10% +8.4% A simple model to explain the reduced induced drag of the ducted tip is shown in Figure 4.15. The downwash caused by shed circulation from the duct is less (by the factor 2 0, where 0 is the angle between a portion of the duct and a portion of the planar cos wing) than the downwash would be if the same shed circulation were redistributed along the duct diameter, in the spanwise direction. The reduction in downwash causes a consequent reduction in induced drag, as suggested by Cone (1963). Owing to the large parasite drag of the ducted tip geometry, the lift-drag ratio of this geometry is inferior to the conventional tip geometry at low angles of attack. However, for a> 8° the aerodynamic performance of the ducted tip is superior. As shown in Figure Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 94 4.16, the improvement in liftldrag ratio with application of the ducted tip varies from 0% to 6% for a between 8° and 15°. Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 95 DRAG COEFFICIENT 0.25 0.2 CD -10 -5 0.15 0 10 5 20 15 INCIDENT ANGLE IDEGREES1 * Rounded tip • Ducted tip Figure 4.14 Drag coefficients of rounded and ducted tip wings. Re uncertainty in CD is = 7.1 x 10 (The 5 0.0003 and in a is t 0.1°, both at the 95% confidence level). + + + + Figure 4.15 A simple model for the reduction in induced drag of the ducted tip geometry Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 96 Lift-Drag Ratio Improvement of Ducted Tip 7 —6 a 1 —1 9 10 11 12 13 14 15 16 17 INCDENT ANGLE [DEGREES] Figure 4.16 Lift/Drag ratio improvement of the ducted wing relative to the conventional tip (The uncertainty in L/D is 1%, and in a is .i O.1°at the 95% confidence level). Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 97 Figure 4.17 shows the locations of the 92 pressure taps fitted on both the pressure and suction sides of the basic wing. Note that there are more taps on the suction side than the pressure side, and also that the area near the tip has a dense population of taps. From these measured point values of surface pressure, one can obtain the pressure distribution on both sides of the wing by interpolation and extrapolation. Because no pressure taps are installed on tip devices, only the pressure distributions on the basic wing can be compared for different tip configurations. As can be seen in Figure 4.17, the pressure taps are quite sparse in some areas of the wing surface. There were only five pressure taps at the leading edge on the suction side and the pressure tap spacing was also very large in the centre of the wing. We had difficulty getting reasonable interpolated/extrapolated pressure data in those regions. For example, even though the experimental results from the five pressure taps at the leading edge decreased monotonically along the spanwise direction, the 2-D inter/extrapolation method we employed yielded a wavy pressure distribution. This difficulty was overcomed by combining l-D and 2-D inter/extrapolation techniques in this region. To give some indication of the accuracy of the interpolation procedure, the surface pressure distribution was numerically integrated to yield the lift force. The lift coefficients evaluated in this way at 5°, 8° and 12° angles of attack were compared with the force balance measurements. The maximum error was 7%. The error was likely due to poor resolution of the suction surface pressure peak owing to the large pressure tap spacing near the wing leading edge. The pressure distribution on the pressure side of the basic wing with a conventional wing tip configuration is shown in Figure 4.18. Figure 4.19 shows the pressure distribution for the ducted tip configuration. Both of the two figures show the local reduction in pressure caused by the outboard crossflow near the tip region (fluid flow from high pressure to low pressure). It is interesting to observe in Figure 4.19 a pressure peak near the 35% chord position at the tip. This location marks the entrance to the duct, and the pressure peak is likely due to the existence of a stagnation point at the entrance of the ducted tip. The pressure behaviour in the leading edge-tip region is obviously different (Figure 4.20) for the two configurations because the conventional tip had a tip Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 98 extension all the way to the leading edge, whereas the ducted tip did not extend to the leading edge and therefore had more three-dimensional flow near the leading edge. The average pressure gradient from the wing root to the tip on the pressure side is related to the magnitude of the outboard crossflow and is characteristic of the strength of the tip vortex. We can evaluate the average pressure gradient qualitatively by looking at the tip views (Figures 4.18 (b) and 4.19 (b)). There is no significant difference between this pressure gradient on the conventional wing and on the ducted tip wing except near the wing tip. Figure 4.20 verifies this conclusion. The three dimensional flow in the tip region is also evident on examination of the pressure distribution on the wing’s suction side (Figure 4.21 and 4.22). There is no significant difference, on the majority of the wing surface, between the pressure distribution of the conventional configuration and the ducted tip configuration (Figure 4.23). Close to the midspan of the leading edge, the difference in between the two configurations is about 10%. This difference is within the error caused by the extrapolation method we used. Consider pressure taps numbered 43 and 44 (Figure 4.17(b)) for example. A 7% error in the pressure measurement at those two taps will cause a maximum 10% error in the estimated c,, at the leading edge, if a linear extrapolation method is used. There is possibly a local suction peak in the trailing edgetip region, as shown more clearly in the leading edge view (Figures 4.21(b) and 4.22(b)). The measurements do not capture the whole local peak because, due to the add-on tip, no pressure measurements were made within 10% of the semi-span from the tip. Figure 4.24 shows the pressure distribution on the suction side of the basic wing without any tip device installed. There is obviously a suction peak near the wing tip, near the locus of 3D stagnation points (Point B in Figure 4.6), caused by the roll-up of the primary tip vortex. Such a suction peak has also been observed by Chow et al. (1993) in their pressure distribution measurements on a wing with a rounded tip. In summary, no pressure distribution difference between tip geometries was observed over the majority of the wing surface. However, near the wing tips the pressure distribution was significantly changed by the addition of tip devices. The pressure measurements are not accurate enough (about ±7%) to confirm the LiftlDrag ratio Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing improvement, but they are supportive to the cavitation results — 99 tip vortex cavitation is significantly modified by the ducted tip, whereas leading edge cavitation is not significantly affected by addition of the tip device. o. . 2 f ccE 1 -1.8 -1.6 -1.2 -8.8 -8,6 -8.2 .ESSL -0.1-0.00.1 0.2 0.3 0.4 0.5 Chapter 4 Cajtatjon andAerodyflam. Peopma,lce ofthe Ducted Tip Wing 101 I (a) L Co c•’J Leading Edge 0.99 0.71 B4 0 . TIp 0.i (b) Figure 4.18 Pressure distribution on the pressure side of the Wing with the conventional tz, _—8° (a) the 3-D view (b) the ti view. Chapter 4 Cavitation and Aerodynamic Peiformance of the Ducted-Tip Wing 102 (a) d C-, c\j d Leading Edge I I 0.89 0.71 0.54 0.36 0.19 0.00 TIP (b) Figure 4.19 Pressure distribution on the pressure side of the wing with the ducted tip, a=8°. (a) the 3-D view, (b) the tip view. Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 103 Uj (1) (1) 0. ‘0 Figure 4.20 The dzfference between the ducted and conventional tip pressure coefficients on the pressure side of the basic wing, a = 8 . . ( . — % ‘— m E & & (_) & cD & - - - o II -1.8 I / Co (D 0.’ rn-i —i__i__i_I. 1 • / CQE. -0.8 -0.5 -1.1 -0.8 COE. 5/ “ -0.2 “\ “5” \“5’ -0.5 -0.2 1:’ D -tic’) JTh IM” 7 • J / /// II PRESSURE /7 Ill,. -1.1 PRESSURE -1.4 -1.4 /1 -1.8 / 0 III -2.1 - // & 0’ & -2.1 1% ‘1 ‘I’ Iy 0 p’ 4 9 0 h f 1 -1.8 -1.5 -1.2 -9.8 -8.6 -9.2 . 2 PRESSIJF CcE. II _4 . o 5-- —S 5-- r I-’ . -4 o Co fr 5— —S —I — - — H So —, -2.0 6)- 6) & c_I 0’ 6) 6) 6) 6) —4 & & & —4 & & I / C C ‘-3 (D I / )I/I/ (ii III I/I ) 1 U)) ( ( / / //// / I/I — I PRESSURE COE. .—_..——7////I /1 I )) -0.2 ‘. -0.2 I 1.4-1.1-0.8-0.5 / ///// — -1.1 -0.6 -0.6 I ,,.,. ii,.. — /1 /1 — / I,, ,‘ / I/il / / ii,,, -1.4 I, —.- / I I -1.7 PRESSURE COE. D)C) UCI) (DC 0 *7 0 leo ‘I’ 0 V’ / ,i,,,, i,,,i,,,, i,,,, ,j-2.9-f.7-I.4-f.f -8.8-0.6-8.2 PRESSLRE CLE . o. 1% 0 -e31 -0.17 -0.04 0.10 0.23 PRESSURE COE. 0.37 C Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing ,— %‘ 107 Suction Peak Lu Figure 4.24 Pressure distribution on the suction side of the basic wing without tip device, a 60 Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 108 4.4.3.2 The bi-wing tip In order to improve the poor lift efficiency of the ducted wing surface, a hi-wing configuration was constructed, which was an improvement over the ducted tip in two respects. One respect was the ease with which the angle of attack of the hi-wing could be varied with respect to the wing chord. Another advantage of the hi-wing tip was that it consists of two smooth lifting surfaces. We chose NACA 64-0006 airfoil sections to build the hi-wing because that airfoil section has good stall characteristics. Three variants of the hi-wing tip were studied — “underloaded”, “meanloaded” and “overloaded”. The “underloaded” bi-wing has a local angle of attack that is less than that of the basic wing, and vice versa for the “overloaded” attachment. Experimental results for two angles of attack show that the underloaded arrangement of the hi-wing tip is the best of the three geometries. A comparison between hi-wing tip configurations and the simple extended wing tip is shown in Table 4.3. The overload attachment produces a 1.4% increase in lift coefficient due to both the relatively higher angle of attack at the tip and the increased area of the hiwing structure. The total drag coefficients of the three hi-wing tip configurations were bigger than that of a conventional tip at this angle of attack. At the same angle of attack, as mentioned before, the drag coefficients were virtually identical for the ring wing tip and conventional tip configurations. This difference in performance is due to the existence of sharp edges and corners in the hi-wing tip configuration. Those sharp edges and corners will generate interference vortices and possibly cause local flow separation. Table 4.3: Comparison of Bi-Wing Tip to Conventional Tip at 8° AngIe of Attack. Conventional Tip CD CL CL/CD 0.0188 0.4699 25.00 Overload Bi-Wing Tip 0.0234 +24% 0.4767 +1.4% 20.4 -18% MeanloadBi-WingTip 0.0211 +12% 0.4699 0% 22.3 -11% +7% 0.464 1 -1.2% 23.0 -8% Underload Bi-Wing Tip 0.0202 Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 109 The total drag coefficient at a=14° was successfully reduced by 4% by using an underloaded bi-wing tip configuration (Table 4.4). This reduction implies that there is an identical drag reduction mechanism for both ring wing tip and bi-wing tip configurations at elevated angles of attack. The higher lift coefficient attained by the bi-wing tips relative to the ducted tip (at a=14°) implies that the rationale for the bi-wing tip design diminished separation at high c — — was valid. Table 4.4: Comparison of the Bi-Wing Tip to the Conventional Tip at x=14°. Conventional Tip CD CL CL/CD 0.1140 0.7093 6.22 Overload Bi-Wing Tip 0.1142 +0.2% 0.7 141 +0.7% 6.25 +0.5% Meanload Bi-Wing Tip 0.1141 +0.1% 0.7021 -1% 6.15 -1% Underload Bi-Wing Tip 0.1091 -4.4% 0.6949 -2% 6.37 +2.4% 4.5 Experimental Errors and Applicability of the Tip Device There are several factors that affect the accuracy of lift and drag force measurements. These factors are: (1). angles of attack; (2). wind speed in the test section of the wind tunnel; (3). the electronics in data acquisition system. The effect of the incident angle can be eliminated by testing different tip devices with the angle of attack fixed. The wind speed in the wind tunnel is measured by a manometer, that reads the pressure difference between the test section and the settling chamber (refer to Figure 4.3) accurately to 0.05mm alcohol. The wind speed can thus be controlled at the testing speed, 3OmIs, with an accuracy of better than 0.1%. The data acquisition system collects 10,000 samples within one second. The datum for lift or drag force is the average of those samples. In several minutes, a number of data for each set-up can be obtained. It Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 110 was found that lift force can be measured to ±0.5% and drag force to ±1% at elevated angle of attack. All experiments were redone several weeks later to confirm the repeatability of data. The lift and drag measurements described in this chapter are of high precision (lift typically measured to ±0.5% and drag to ±1% at elevated angle of attack), but fairly low accuracy because of the interaction of the wind tunnel boundary layer with the wing. However, it is known (Green and Acosta 1991) that in the experimental configuration the domain of influence of the wing/wall interaction is confined to a region within 0.2 chords of the wing root. Consequently, it is highly unlikely that modifying the tip geometry would have a significant impact on this interference drag. It is therefore reasonable to argue that although the lift and drag measurements described here are likely to suffer from large errors for any particular wing geometry owing to wind tunnel floor interference, comparison between different wing geometries is valid because this interference effect is nearly constant. The ring-wing tip modification improves the Lift/Drag performance of a small aspect ratio rectangular planform wing. Therefore, it seems worthwhile to speculate on the potential applications of this device. The most apparent immediate application of this technology is to marine propellers. Marine propellers, which are of small aspect ratio and often operate at high lift coefficients, would benefit from both the superior cavitation behaviour and the improved Lift/Drag ratio of the ring-wing tip. The potential for applying this technology to aircraft wings is much less evident. Aircraft wings commonly operate at a moderate lift coefficient, and have much larger aspect ratios than studied here. Both these factors would mitigate against the installation of ring-wing tips on normal aircraft wings. 4.6 CONCLUSIONS A novel hydrofoil tip geometry, consisting of a flow-through duct attached to the hydrofoil tip, has been tested for both cavitation inception and lift/drag performance. The principle behind the ducted tip design is to cause the Trefftz-plane shed circulation of the Chapter 4 Cavitation and Aerodynamic Performance of the Ducted-Tip Wing 111 hydrofoil to take the form of a line with an attached circle, rather than the simple line of a conventional hydrofoil. The tip vortex inception index of the ducted tip geometry is 50 ± 15% less than that of a conventional tip at normal operating angles (e.g. a = 7°), and is at least 30% less for all positive angles of attack. Because the ducted tip hydrofoil shed vorticity has been redistributed relative to the baseline hydrofoil, the induced drag on the hydrofoil should also be reduced. Such a reduction in the induced drag has been observed. For a > 8° the total drag of the ducted tip hydrofoil is less than that of a conventional hydrofoil, despite the fact that the parasite drag of the conventional hydrofoil is less, owing to its much reduced wetted area. The overall lift-drag ratio of the ducted tip hydrofoil, for a > °, is up to 6 ± 1% greater than that of a conventional tip. Thus, the much improved cavitation behaviour of the ducted tip comes at no cost (for elevated a), and even some benefit, in terms of the hydrodynamic (non-cavitating) performance of the hydrofoil. The ducted tip geometry is most advantageous when the lifting surface is highly loaded. In practice, propellers on tugs and fishing trawlers are often highly loaded. Studies are presently underway to explore the effectiveness of the ducted tip geometry on such marine propellers. Chapters Summary and Conclusions 112 Chapter 5 Summary and Conclusions This final chapter is divided into two parts. The first section summarizes the “quasi similarity” methodology and conclusions. The second section makes some suggestions for future work. 5.1 Summary and Conclusions A wing tip vortex was studied both analytically and experimentally. In the analytical work, a new method was developed that describes the tip vortex structure from far downstream to a downstream distance fairly close to the vortex-generating wing. Wing drag is the key parameter affecting the tip vortex structure. In the experimental work, a hydrofoil with a novel ducted tip device was found to have superior cavitational and aerodynamic performance over a conventional rounded tip configuration. This section is divided into three sub-sections: analytical method, tip vortex structure and wing tip modification. 5.1.1 Analytical method A new analytical method, referred to as the “quasi-similarity” method, has been developed for modelling vortices in a freestream. The new approach combines a polynomial solution with the similarity variable technique. By using this method, the non linear partial differential governing equations are reduced to sets of ordinary differential equations. The fundamental idea of the “quasi-similarity” method can be applied to a broad range of non-linear problems that cannot be described by self-similar solutions. The “quasi-similarity” method is used here to model wing tip vortex flow. It provides the first non-linear analytical wing tip vortex model. In this new method, each resulting set of ordinary differential equations describes the flow with varying degrees of faithfulness. The first order set of ODEs gives solutions that are second order accurate in the axial variable z; higher order sets of ODEs provide higher order accuracy. Fortunately, ChapterS Summary and Conclusions 113 the first order terms of the solution polynomial can be obtained analytically in complete function form. The second order tangential velocity component and pressure distribution are also obtained in complete function form. Other higher order terms must be calculated numerically by solving a set of ODEs with appropriate boundary conditions. There are three major assumptions of the quasi-similarity method. These assumptions are: (1). the axial velocity deficit or excess in tip vortices is not significantly greater than the freestream velocity; (2). the vortex circulation remains constant with downstream distance; (3). the flow is laminar. The first limitation can be expressed mathematically (assuming a rectangular wing) as z> CRe.S/l67u where z is the downstream distance, CD is the wing total drag coefficient, Re is the Reynolds number based on the wing chord, and S is the span of the wing. This distance is about a hundred meters downstream of the wing trailing edge for a typical drag coefficient CD O.O1, Reynolds number Re 5xlO , and a one meter span wing (for 5 which the flow should be laminar). The second assumption is a reasonable one after the vortex roll-up is complete. Therefore, very near the wing the first two conditions may not apply. The flow in a practical trailing vortex is usually turbulent, at least near the wing, due to large Reynolds number. This theory is still valid for a turbulent flow only if an eddy viscosity model is used and the eddy viscosity is a constant everywhere in the flow. 5.1.2 Tip vortex structure All three velocity components in tip vortices decrease with increasing downstream distance at fairly large downstream distance of the vortex-generating wing. The radial velocity diminishes most rapidly, followed by the axial velocity. The tangential velocity decays most slowly. This prediction is consistent with many experimental observations. Chapters Summary and Conclusions 114 At a fairly large distance downstream of the wing, only an axial velocity deficit can exist for a wing experiencing a positive drag force. An axial velocity excess can exist far downstream of a wing generating thrust (e.g. with wing tip-mounted engines). A unique and important feature of the quasi-similar wing tip vortex model is that the predicted drag on a vortex-generating wing due to a single quasi-similar tip vortex is finite. To the author’s knowledge, no other vortex model has this property. The tangential velocity component and pressure of the first order quasi-similarity solution are exactly the same as those in an advanced linear tip vortex model (Batchelor 1964). The axial velocity field of the linear vortex model differs from that of the quasisimilarity model. This discrepancy in the axial velocity causes a single vortex in the linear model to have an infinite drag whereas the quasi-similar model indicates a finite drag. The axial velocity and tangential velocity predicted by quasi-similarity theory have been compared with experimental measurements. The agreement is good at large downstream distances from a wing. The experimental results confirm that the axial velocity deficit on the vortex centerline decays as the reciprocal of the downstream distance; the tangential velocity in a tip vortex is proportional to the inverse square root of the downstream distance. 5.1.3 Wing tip modification A novel ducted tip device has been tested in a wind tunnel and a water tunnel. Aerodynamic and cavitation performance of the ducted tip has been compared with that of a conventional rounded tip. The ducted tip has up to 6% better Lift/Drag performance than a square cut wing tip* of equal span at higher angles of attack, roughly 8° or more. The principal cause of these performance improvements is the greatly reduced induced drag afforded by the ducted tip at high lift coefficients, which more than offsets the increase in parasitic drag resulting from the increased wetted area. Experimental data show that a square cut wing has similar aerodynamic performance to a rounded tip wing. The difference between those two conventional tip configurations is within experimental error. ChapterS Summary and Conclusions 115 The tip vortex cavitation inception index of the ducted tip geometry is at least 33% lower than that of a conventional rounded hydrofoil tip at all positive angles of attack. No pressure distribution difference between tip geometries was observed over the majority of the wing surface (differences from one tip geometry to another were less than the experimental error, which is about 7%). However, near the wing tips the pressure distribution was significantly modified by the addition of tip devices. The pressure measurements are not accurate enough to confirm the Lift/Drag ratio improvement, but they are supportive of the cavitation results — tip vortex cavitation is significantly modified by the ducted tip, whereas leading edge cavitation is not significantly affected by addition of the tip device. The ducted tip is effective at reducing the tip vortex inception index and the vortex induced drag. This tip design spreads the shed vorticity in the transverse plane behind the wing (Trefftz plane) over a line and circle, while a conventional wing tip spreads the vorticity only over a line. The vastly reduced tip cavitation index, together with the improved Lift/Drag performance, suggests there is a significant role for the ducted tip device in both military and civilian marine propeller design. 5.2 Suggestions for future work A careful experiment should be designed in order to further verify the quasi-similarity vortex model. The wing total drag and vortex circulation should be measured accurately in a direct or indirect way with the Reynolds number maintained low enough that the vortex flow is laminar. Application of the “quasi-similarity” method to other non-self-similar problems could be possible. In fact, jets with a tangential velocity component (twisted jets) were studied in detail by Loitsyanski (1966) and Goldshtik et al. (1979 and 1986). They employed a non-similar power series to construct their analytical solutions. However, twisted jets are different from wing tip vortices. The swirling jets are characterized by the following conservation integrals: flowrate, axial momentum and angular momentum. In contrast, Chapter 5 Summary and Conclusions 116 the tip vortices are characterized by a freestream velocity, and the conservation of vortex circulation and axial momentum. Related problems, for instance, jet flows (with or without swirl) in a freestreàm seem amenable to treatment using the quasi-similarity technique. A different set of polynomials, with an amplitude function proportional to z 1 instead of the liz’ used in the theory developed in this thesis, is a possible solution of the simplified governing equation (2.1.1). This amplitude function has the general feature of vortex rollup — the magnitude of velocity components increases with increasing downstream distance. Is this a solution pertinent to the roll up of the wing tip vortex? Are there any physical problems related to the new polynomial solution? Further investigations are needed. Instability of a trailing line vortex has attracted great interest in recent years. Some of the stability studies based on Batchelor’s partly-linearized theory include works by Stewartson and Capell (1985), Stewartson and Brown (1985), Mayer and Powell (1992b), and Duck et al. (1992). The results in this thesis may prompt studies of the effect of wing drag on the stability of a single tip vortex, and the non-linear effect on tip vortex instability. The potential for applying the novel ducted tip to aircraft wings is quite small, because aircraft wings commonly operate at moderate lift coefficients and have much larger aspect ratios than studied here. The most apparent immediate application of ducted tips is to marine propellers. Marine propellers, which have small aspect ratio and large operating angles of attack, would benefit from both the superior cavitation performance and the improved Lift/Drag ratio of a ducted tip. Future work should concentrate on further improving the tip geometry and establishing whether the ducted tip has a practical marine application. Bibliography 117 Bibliography Arakeri, V.H., Sharma, S.D., and Mani, K., 1985, “A Technique to Delay the Inception of Tip Vortex Cavitation from Marine Propellers,” ASME Cavitation and Multiphase Flow Forum. Arndt R. E. A. and Keller A. P., 1992, “Water Quality Effects on Cavitation Inception in a Trailing Vortex,” J. Fluid Eng. Vol. 114, pp. 430-438. Arndt R. E. A., Arakeri V. H. and Higuchi H., 1991, “Some Observations of Tip-Vortex Cavitation,” J. Fluid Mech. vol.229, pp. 269-289. Baber, M.R. et al., 1975, “Flight Test Investigation of the Vortex Wake Characteristics Behind a Boeing 727 During Two Segment and Normal 1LS Approaches”, NASA TM X-62,398. — Baker, G.R., Barker, S.J., Bofah, K.K., and Saffman, P.G., 1974, “Laser Anemometer Measurements of Trailing Vortices in Water”, J. Fluid Mech., Vol. 65, Part 2, pp. 325336. Bandyodpadhyay, P.R., Stead, D.F., and Ash, R.L., 1990, “The Organized Nature of a Turbulent Trailing Vortex”, 21st Fluid Dynamics, Plasma Dynamics And Lasers Conference, Seattle,Wa. Barnett R. A. and Ziegler M. R., 1991, Applied Mathematics, Dellen Macmillan. Batchelor, G.K., 1964, “Axial Flow in Trailing Line Vortices”, J. Fluid Mech., Vol. 20, Part 4, pp. 645-658. Benjamin, T. B., 1962, “Theory of the Vortex Breakdown Phenomenon.” I. Fluid Mech. 14, . 629 593 pp. Betz, A., 1933, “Behavior of Vortex Systems”, NACA TM 713. Burgers 3. M., 1948, “A Mathematical Model Illustrating the Theory of Turbulence”, Advances in Applied Mechanics, pp.171-196. Carlin, G., Dadone, L. and Spencer, R., 1989, “Results of an Experimental Investigation of Blade Tip Vortex Modifcation Devices”, NASA CR- 181853. Chahine G. L. Frederick G. F. and Bateman R. D., 1993, “Propeller Tip Vortex Cavitation Suppression Using Selective Polymer Injection,” 3. Fluids Eng. vol. 115, pp. 497-503. Bibliography 118 Chigier, N.A. and Corsiglia, V. R., 1972, “Wind Tunnel Studies of Wing Wake Turbulence”, AJAA #72-41. Chow, J. S., Zilliac, G. G. and Bradshaw, P., 1993, “Measurements in the Near-Field of a Turbulent Wing Tip Vortex”, AIAA-93-055 1. Ciffone, D. L., and Orloff, K. L., 1974, “Axial Flow Measurement in Trailing Vortices,” AIAA Journal, Vol.12, No.8, . 1155 1154 pp. Cone, C.D., 1963, “The Theory of Induced Lift and Minimum Induced Drag of Non Planar Lifting Systems,” NASA Technical Report, R-139, N63 16288. Corsiglia, V. R., Schwind, R. G. and Chigier, N. A., 1973, “Rapid Scanning, 3-D HotWire Anemometer Surveys of Wing Tip Vortices,” Journal of Aircraft, Vol.10, No.12, pp.752-757. — Critchley, J.B., 1991, “Wake Vortex The Program in the United Kingdom,” Proceedings of the FAA Wake Vortices Conference, Washington, D.C. — Crow, S. C., 1970, “Stability Theory for a Pair of Trailing Vortices”, AIAA Journal, Vol. 8, No. 12, pp.2172-2179. Devoss, D., 1986, “The Race To Recover The Cup,” Los Angeles Times Magazine, Aug. 31, pp.9+. Donaldson, C. and Bilanin, A. J., 1975, “Vortex Wakes of Conventional Aircraft”, AGARD-AG-204. Donaldson, C. du P. and Sullivan, R. D., 1970, “In Aircraft Wake and Its Detection,” Plenum, New York, p389. Donaldson, C., 1971, “A Brief Review of the Aircraft Trailing Vortex Problem”, AFOSR TR-71-1910. Duck, P.W. and Khorrami M. R., 1992, “A Note on the Effects of Viscosity on the Stability of a Trailing-Line Vortex,” JFM, vol. 245, . 189 175 pp. El-Ramly, Z., and Rainbird, W. J., 1977, “Flow Survey of the Vortex Wake Behind Wings.” Journal of Aircraft, Vol.14, No.11, pp. 1 102-1108. Eliason, B. G., Gartshore, I.S. and Parkinson, G. V., 1975, “Wind Tunnel Investigation of Crow Instability.” J. Aircraft 12, pp.985-988. Falcao de Campos, J.A.C., George, M.F., and MacKay, M., 1989, “Experimental Investigation of Tip Vortex Cavitation for Elliptical and Rectangular Wings,” ASME Cavitation and Multiphase Flow Forum, La Jolla, California. Bibliography 119 Flechner, S.G., Jacobs, P.F. and Whitcomb, R.T., 1976, “A High Subsonic Speed Wind Tunnel Investigation of Winglets on a Representative Second Generation Jet Transport Wing” NASA TN D-8264. , Fruman D. H., Cerrutti P., Pichon T. and Dupont P., 1995, “Effect of Hydrofoil Planform on Tip Vortex Roll-up and Cavitation,” J. Fluids Eng., vol. 115, pp. 162-169. Fruman, D.H., Dugue, C. and Cerruti, P., 1991, “Tip Vortex Roll-up and Cavitation,” ASME Cavitation and Multiphase Flow Forum, Portland, Oregon. Gartshore, I. S., 1963, “Some Numerical Solutions for the Viscous Core of an Irrotational Vortex,” NRC Can. Aero. Rep. LR-378. Goldshtik M. A., 1979, “On Twisted Jets”, Izv. Akad. Nauk. SSSR, MzhG, 1. Goldshtik, M. A. and Yavorskii, N. I., 1986, “On Submerged Jets.” PMM USSR, Vol.50, No. 4, . 445 438 pp. Gorsigtlia, V.R., Schwing, R.G. and Chigier, N.A., 1973, “Rapid Scanning, Three Dimensional Hot-Wire Anemometer Surveys of Wing Tip Vortices”, J. Aircraft, Vol. 10, No. 12, pp. 752-757. Govindaru, S.P. and Saffman, P.G., 1971, “Flow in a Turbulent Trailling Vortex”, Physics Of Fluid, Vol.14, Number 10, Oct. 1971, pp. 2074-2080. Green, S. I., 1988, “Tip Vortices Signie Phase And Cavitating Flow Phenomena”, Ph.D. Thesis, California Institute Of Technology, Report No. Eng. 183.17. — Green, S. I. and Duan, S. Z., 1995, “The Ducted Tip A Hydrofoil Tip Geometry with Superior Cavitation Performance,” accepted by the ASME Journal of Fluids Engineering. — Green, S. I. and Acosta, A.J., 1991, “Unsteady Flow in Trailing Vortices,” Journal of Fluid Mechanics, Vol. 227, pp. 107-134. Green, S. I., 1991, “Correlating Single Phase Flow Measurements with Observations of Trailing Vortex Cavitation,” Journal of Fluids Engineering, Vol. 113, No. 1, pp. 125-129. Green, S. I., 1994, “Wing Tip Vortices,” to appear in Fluid Vortices, S. Green (ed.), Kiuwer, Dordrecht. Green, S. I., Acosta, A.J. and Akbar, R., 1988, “The Influence of Tip Geometry on Trailing Vortex Rollup and Cavitation,” ASME Cavitation and Multiphase Flow Forum, Cincinnati, Ohio. Green. S. I., 1987, “Three Tip Vortex Examination Techniques”, International Symposium on Cavitation Research and Techniques, pp. 155-160. Bibliography 120 Hall, H.W., Arndt, R.E.A. and Billit, M.L., 1972, “Limited Cavitation and the Related Scale Effects Problem”, JSME Proc. 2nd mt. Symp. Fluid Mechanics And Fluidics, Tokyo, pp. 303-314. Hastings, E.C. Patterson, J.C., Shanks, R.E., Champine, R.A., Copeland, W.L., and Young, D.C., 1975, “Development and Flight Tests of Vortex Attenuating Splines”, NASA TN D- 8083. Higuchi, H, Quadrelli, J.C. and Farell C., 1986, “Vortex Rollup From an Elliptical Wing at Moderately Low Reynolds Numbers”, AIAA #86-0562. Hoeijmakers, H.W.M., 1983, “Computational Vortex Flow Aerodynamics”, AGARD Conference Proceedings #342, 1 pp. 8.1-18.35. Hoffmann, E.R., and Joubert, P. N., 1963, “Turbulent Line Vortices,” Journal of Fluid Mechanics, Vol.16, . 411 395 pp. Huang, T., 1987, David Taylor Naval Ship Research And Development Center, Personal Communication. Ikohagi, T., Higuchi, H. and Avndt, R.E.A., 1986, “The Structure of Trailing Vortices”, In Advancement In Aerodynamics, Fluid Mechanics And Hydraulics, Joint AIAA, ASME And ASCE Symposium, Univ. of Minnesota. Chow, J. S., Zilliac, G. G. and Bradshaw, P., 1991, “Initial Roll Up of a Wing Tip Vortex”, Personal Communication. Kaden, H., 1931, “Aufwickling einer unstabilen unstetigkeitsflache”, Ing. Arch. 2, 140. (English translation, R. A., Library translation no. 403) Kantha, H.L., Lewellen, W.S. and Durgin,F.H., 1971, “Qualitative Responses of a Vortex Core to Tip Blowing and Intersecting Airfoils”, ASRL TR 153-4. Katz, I. and Galdo, J. B., 1989, “The Effect of Roughness on Rollup of Tip Vortices on a Rectangular Hydrofoil.” J. Aircraft 26, . 253 247 pp. Kuiper, G., 1981, Cavitation Inception on Ship Propellor Models, M. Veenamnen Zonen B.V. Wageningen, pp. 130-159. Lessen, M. and Paillet, F., 1974b, “The Stability of a Trailing Line Vortex, Part 2, Viscous Theory”, J. Fluid Mech., 65, . 779 769 pp. Lessen, M., Singh, P. J. and Paillet, F., 1974a, “The Stability of a Trailing Line Vortex, Part 1, Inviscid Theory”, J. Fluid Mech., 63, pp.753-763. Bibliography 121 Lewy, S. and Caplot, M., 1982, “Theoretical and Experimental Study of Helicopter Rotor Noise”, ONERA TP 1982-74. Lignueul P. and Latorre R., 1993, “Study of Nuclei Distribution and Vortex Diffusion Influence on Nuclei Capture by a Tip Vortex and Nuclei Capture Noise,” J. Fluid Eng, vol. 115, pp. 504-507. Logan A. H., 1971, “Vortex Velocity Distributions at Large Downstream Distance”, J. Aircraft, vol.8, No. 11. pp.930-932. Loitsyanski, L. G., 1966, “Mechanics of Fluids and Gases.” Oxford University Press. Long, R. R., 1961, “A Vortex in an Infinite Viscous Fluid”, J. Fluid Mech., 11, pp.611624. Mager, A., 1972, “Dissipation and Breakdown of a Wing-Tip Vortex”, I. Fluid Mech. 55, pp.609-628. Maines, B.H. and Arndt, R.E.A., 1993, “Bubble Dynamics of Cavitation Inception in a Wing Tip Vortex,” ASME Cavitation and Multiphase Flow Forum, Washington, D.C. Martin, R.M., Elliott, J.W., and Head, D.R., 1984, “Comparison of Experimental and Analytical Predictations of Rotor Blade Vortex Interactions Using Model Scale Acoustic Data”, AIAA #84-2269. Mason, W.H. and Marchman J. F., 1972, “Far Field Structure of an Aircraft Trailing Vortex, Including Effects of Mass Injection”, NASA CR-62078. Mayer, E. W. and Powell, K. G., 1992a, “Similarity Solutions for Viscous Vortex Cores,” JFM, vol.238, . 507 487 pp. Mayer, E. W. and Powell, K. G., 1992b, “Viscous and Inviscid Instabilities of a Trailing Vortex”, JFM, vol. 245, . 114 91 pp. McCormick, B.W., 1962, “On Cavitation Produced by a Vortex Trailing from a Lifting Surface,” ASME Journal of Basic Engineering, pp. 369-379. McCormick, B.W., Tangier, J. L., and Sherrieb, H. E., 1968, “Structure of Trailing Vortices,” Journal of Aircraft, Vol.5, No.3, . 267 260 pp. Milne—Thomson, L.M., 1958, Theoretical Aerodynamics, Dover, New York. Moordzij, L., 1977, “A Note on the Scaling of Tip Vortex Cavitation Inception”, International Shipbuilding Progress, No. 24, pp. 233-236. Moore, D. W. and Saffman, P. G., 1973, “Axial Flow in Laminar Trailing Vortices,” Proceedings of the Royal Society of London, Series A. Vol.333, . 508 491 pp. Bibliography 122 Moore, D.W., 1974, “A Numerical Study of the Rolling Up of a Finite Vortex Sheet”, J. Fluid Mech., Vol. 63., Part 2, pp. 225-235. Moore, D.W. and Saffman, P.G., 1973, “Axial Flow in Laminar Trailing Vortixes”, Proc. Rey. Soc. Lond. A. Vol. 333, pp. 491-508. Mosher, M. and Peterson, R.L., 1983, “Acoustic Measurements of a Full-Scale Coaxial Helicopter”. AIAA #83 -0722. Olsen, J.H., 1971, “Results of Trailing Vortex Studies in a Towing Tank”, Publ. In “Aircraft Wake Turbulence and Its Detection”, Ed. By J.H. Olsen et al. Plenum Press, New York, pp. 455. Ooi, K.K., 1981, “Scale Effects on Cavitation Inception in Submerged Jets”, Ph.D. Thesis, California Institute Of Technology, Report No. 183-6. Orloff, K.L. and Grand, G.R., 1973, “The Application of a Scannig Laser Doppler Velocimeter to Trailing Vortex Definition and Alleviation”, AIAA #73-680. Oseen, C. W., 1912, Uber wirbelbewegung ineiner reibenden flussigkeit, Art. f. Mat. Astron. Fys. 7.14. Pauchet, A., Briancon-Marjollet, L. and Fruman, D.H., 1993, “Recent Results on the Effect of Cross-Section on Hydrofoil Tip Vortex Occurrence at High Reynolds Numbers,” ASME Cavitation and Multiphase Flow Forum, Washington, D.C. Peckham, D. H. and Atkinson, S. A. 1957, “Preliminary Results of Low Speed Wind Tunnel Test on a Ghotic Wing of Aspect Ratio 1.0”, British Aeronaut. Res. Council, CP508. Phillips W., 1981, “Turbulent Trailing Mechanics, Vol. 105, . 467 451 pp. Vortex During Roll Up”, Journal of Fluid Platzer, G.P. and Souders, W.G., 1979, “Tip Vortex Cavitation Delay with Application to Marine Lifting Surfaces A Literature Survey,” DTNSRDC Report 79/051. — Poling, D.R., Dadone, L. and Telionis, D.P., 1989, “Blade-Vortex Interaction,” AIAA Journal, Vol. 27, No. 6, pp. 694-699. Prandtl, L., 1920a, “Theory of Lifting Surface Part I” NACA TN #9. , Prandtl, L., 1920b, “Theory of Lifting Surface Part II”, NACA TN #10. Raines, D.A., 1954, “Tip Clearance Flows in Axial Compressors and Pumps”, Report M.E.-5, Jet Propulsion Center, California Institute Of Technology. Bibliography 123 Ruden, p., 1974, “Investigation of a Single Stage Axial Fan” NACA TM-1062. , Sarpkaya, T., 1983, “Trailing Vortices in Homogeneous and Density Stratified Media”, J. Fluid Mechs., Vol. 136, pp. 85-109. Sarpkaya, T., 1992, “Interaction of a Turbulent Vortex With a Free Surface”, Nineteenth Symposium On Naval Hydrodynamics. Seoul, Korea. Sarpkaya, T. and Daly, J. J., 1987, “Effect of Ambient Turbulence on Trailing Vortices.” J. Aircraft 24, . 404 399 pp. Scheiman, J. and Shivers, J., 1971, “Exploratory Investigation of the Structure of the Tip Vortex of a Semispan Wing for Several Wing Tip Modifcation”, NASA TN D-6 101. Sharma, S. D., Mani, K and Arakeri, V. H., 1990, “Cavitation Noise Studies on Marine Propellers,” Journal of Sound and Vibration, Vol .138, No. 2, pp. 255-283. Shekarriz, A., Fu, T.C., Katz, J., Liu, H.L. and Huang, T.T., 1992, “Study of Junction and Tip Vortices Using Particle Displacement Velocimetry,” AIAA Journal, Vol. 30, No. 1, pp. 145-152. Shekarriz, A., Fu, T. C. and Katz, J., 1993, “Near-Field Behavior of a Tip Vortex”, AIAA J., vol. 31, No. 1, . 118 112 pp. Shevell, R. S., 1989, Fundamentals of Flight, Prentice-Hall, p154. Snedeker, R.S., 1972, “Effect of Air Injection on the Torque Produced by a Trailing Vortex”, J. Aircraft. Vol. 9, No. 9, pp. 682-686. Sounders, W.G. and Platzer, G.P., 1981, “Tip Vortex Cavitation Characteristics and Delay of Inception on a Three Dimensional Hydrofoil,” David Taylor Maval Ship Research And Development Center, Rep. No. 8 1/007. Spiliman, 3., 1978, “The Use of Wing Tip Sails to Reduce Vortex Darg”, Aeronautical J., pp. 387-395. Spivey, W.A., Morehouse, G.G., 1970, “New Insights into the Design of Swept Tip Rotor Blades”, 26th Annual National Forum Of The American Helicopter Society. Sponagle, N.C. and Leggat, L. J., 1984, “The Spatial Distribution of Acoustic Sources at Propellor Vortex Cavitation Inception”, ASME International Symposium On Cavitation Inception, pp.191-198. Srinvivasan, G.R. McCrosky, W.J., Baeder, F.D. and Edwards, T.A., 1980, “Numerical Simulation of Tip Vortices of Wings in Subsonic and Transonic Flows”, AIAA Journal, Vol. 26, Number 10, pp. 1153-1162. Bibliography 124 Stewartson K. and Brown S. N., 1985, “Near-Neutral Centre-Modes as Inviscid Perturbations to a Trailing Line Vortex,” JFM, vol.156, pp.387-399. Stewartson, K. and Capell, K., 1985, “On the stability of Ring Modes in a Trailing Line Vortex: the Upper Neutral Points,” JFM, vol. 156, pp.369-386. Stinebring, D.R., Farrell, K.J. and Billet, M.L., 1989, “The Structure of a ThreeDimensional Tip Vortex at High Reynolds Number,” 22nd American Towing Tank Conference, St. John’s, Newfoundland. Stinebring, D.R., Farrell, K.J. and Billet, M.L., 1991, “The Structure of a ThreeDimensional Tip Vortex at High Reynolds Number,” Journal of Fluids Engineering, Vol. 113, No. 3, pp. 496-503. Summa, J.M., 1982, “Advanced Rotor Analysis Methods for the Aerodynamics of Vorex/Blade Interactions in Hover”, Association Aeronautique et Astro-nautiqe de France, 8th European Rotorcraft Forum. Takami, H., 1964, “A Numerical Experiment with Discrete Vortex Approximation with Reference to the Rolling Up of a Vortex Sheet”, Stanford University, SUDAER 202. Thompson, D. H., 1975, “Experimental Study of Axial Flow in Wing Tip Vortices”, J. Aircraft, vol.12, No. 11, pp. 910-911. Trigub, V. N., Blokhin, A. B. and Simakin, I. N., 1994, “The Asymptotic Study of Dissipation and Breakdown of a Wing-Tip Vortex”, J. Fluid Mech., Vol.274, pp.293-337. Vaczy, C.M. and McCormick, D.C., 1987, “A Study of the Leading Edge Vortex and Tip Vortex on Prop-Fan Blades”, J. Turbomachinery, Vol. 109, pp. 325-331. Vortex Wake Alleviation Workshop, 1980, NASA Langley Research Center, Hampton, VA. Webber, G.W. and Dansby, T., 1983, “Wing Tip Devices for Energy Conservation and Other Purposes Experimental and Analytical Work in Progress at the LockheedGeorgia Company,” Canadian Aeronautics and Space Journal., Vol. 29, No. 2, pp. 105120. — Westwater, F.F., 1935, “Rolling Up of the Surface Discontinuity behind an Aerofoil of Finite Space” Aero. Res. Counc. R. And M. # 1962. , Whitcomb, R.T., 1976, “A Design Approach and Selected Wind Tunnel Results at High Subsonic Speed for Wingtip Mounted Winglets”, NASA TN D-8260. Widnall, S.E. and Wolf, T.L., 1980, “Effect of Tip Vortex Structure on Helicopter Noise due to Blade Vortex Interaction”, J. Aircraft, Vol. 17, pp. 705-7 11. Bibliography 125 William, G. M., 1974, “Viscous Modeling of Wing-generated Trailing Vortices,” Aeronautical Quarterly, Vol.25, . 154 143 pp. Winter, H., 1991, “R and D for a Wake Vortex Warning System in Germany,” Proceedings of the FAA Wake Vortices Conference, Washington, D.C. Zalay, A.D., 1976, “Hot Wire and Velocity Meter Wake Vortex Surveys”, AIAA 3., Vol. 15, No.5, pp. 694-696. Zhao, Z. J., 1993, “A System for the Efficient Automated Analysis of Reconstructed Double-Pulsed Holograms.” M.Sc. Thesis, Dept. of Mech. Eng., University of British Columbia.
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Analytical and experimental studies of wing tip vortices
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Analytical and experimental studies of wing tip vortices Duan, Shizhong 1995-12-31
pdf
Page Metadata
Item Metadata
Title | Analytical and experimental studies of wing tip vortices |
Creator |
Duan, Shizhong |
Date | 1995 |
Date Issued | 2009-04-16T23:39:01Z |
Description | Wing tip vortices, and their relationship to wing tip geometry and wing total drag, are investigated here both analytically and experimentally. The purpose of the analysis is to answer the basic question — what are the effects of wing drag on tip vortex structure? The purpose of the experiments is to determine the effect of wing tip geometry on wing tip vortices (tip vortex cavitation) and wing drag (vortex induced drag). In the analytical work, a new method has been developed for wing tip vortices. This approach is referred to as the “quasi-similarity” method. This new method combines a polynomial solution with a similarity variable technique. A non-linear analytical tip vortex model is achieved for the first time. The first order (linear) velocity components and pressure of the polynomial solution for the new model can be expressed in complete function form. Higher order (non-linear) velocity components and pressure can be obtained analytically or numerically. The most important feature of this new tip vortex model is that the predicted wing drag due to such a single vortex is finite. To the author’s knowledge, no other vortex model has this property. At distances fairly far downstream of the wing, tip vortex structure is found to be an explicit function of wing total drag, vortex circulation, freestream velocity, downstream distance and fluid properties. It is verified that the first order tangential velocity component and pressure of the tip vortex model are exactly the same as that of the linear wing tip vortex model proposed by Batchelor (1964). The decay of the axial and tangential velocity components predicted by this theory compare well with experimental measurements. In the experimental work, a novel ducted tip device was tested in a wind tunnel and a water tunnel. The ducted tip consists of a hollowed duct attached to the tip of a rectangular untwisted wing. This novel tip device was found to improve the Lift/Drag ratio by up to 6% at elevated angles of attack compared with a conventional round tip configuration with the same span. The wing tip vortex cavitation was substantially delayed by the ducted tip device. In view of its superior cavitation characteristics and aerodynamic performance, the ducted tip has potential application to marine propellers. The ducted tip is effective because it redistributes the shed vorticity in the transverse plane (Trefftz plane) behind the wing. The shed vorticity is distributed along the wing and duct trailing edge for the ducted tip configuration, rather than solely along the wing trailing edge, which would be the case for a conventional wing tip configuration. |
Extent | 3094251 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-04-16 |
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.0080887 |
URI | http://hdl.handle.net/2429/7256 |
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 | 1995-11 |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-ubc_1995-059480.pdf [ 2.95MB ]
- Metadata
- JSON: 831-1.0080887.json
- JSON-LD: 831-1.0080887-ld.json
- RDF/XML (Pretty): 831-1.0080887-rdf.xml
- RDF/JSON: 831-1.0080887-rdf.json
- Turtle: 831-1.0080887-turtle.txt
- N-Triples: 831-1.0080887-rdf-ntriples.txt
- Original Record: 831-1.0080887-source.json
- Full Text
- 831-1.0080887-fulltext.txt
- Citation
- 831-1.0080887.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
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-0080887/manifest