24 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES EXPERIMENTAL INVESTIGATION OF THE EFFECT OF VARIOUS WINGLET SHAPES ON THE TOTAL PRESSURE DISTRIBUTION BEHIND A WING Mohammad Reza Soltani, Kaveh Ghorbanian, Mehdi Nazarinia Department of Aerospace Engineering, P. O. Box: 11365-8639 Sharif University of Technology, Tehran, 14566, Iran Keywords: Winglets, Pressure Distribution, Tip Vortex, Winggrid, Turbulator Abstract A series of low speed wind tunnel tests to investigate the effect of different winglet shapes and orientations on the flowfield of a swept wing were conducted. Both surface pressure and wake data at various angles of attack, α = -1 to 15 degrees were measured. Further, the flowfield around the wing is visualized using tufts. Four types of winglet shapes, spiroid (forward and aft), blended and winggrid are used in this investigation. The wake profile behind the wing at various angles of attack with different winglets is measured, too. Furthermore, the effect of Turbulator strip on the wing upper surface is examined. Results are compared with conditions of no winglet. The results indicate that the presence of winglets change the flowfield over and around the wing significantly. The total pressure in the wake of the model varies drastically when the wing is equipped with winglets. 1 Introduction From the beginning of aviation, designers are searching for methods and technologies to reduce the required fuel consumption of commercial aircraft. Wingtip devices aim the reduction of induced drag, which are responsible for 3-4% of the total drag of a transport-aircraft at longrange cruise condition and for considerably downgrading the climb performance of an aircraft. Winglets along tip tanks, raked wingtips, and aligned fans belong to this class of devices. During 198s, NASA Dryden Flight Research Center was involved with general aviation research with the KC-135 aircraft. A winglet, developed by Richard Whitcomb of the Langley Research Center, was tested on a KC-135A tanker loaned to NASA by the Air Force. The flight test showed that the winglets could increase the aircraft s range by as much as 7 percent at cruise speeds. The first application of NASA s winglet technology in industry was on General Aviation business jets, but winglets are now incorporated into most commercial and military transport jets, i.e., the Gulfstream III, IV, and V business jets, the Boeing 747-4, McDonnell Douglas MD-11 and Airbus A34, A33, A31-3, A321, airliners, the McDonnell Douglas C-17 military transport, and Embraer aircraft. In recent years, many modification kits are offered for installing winglets to aircrafts, which originally did not have them. A close look at the drag breakdown of a typical civil transport aircraft (Fig.1), reveals that the skin friction and lift induced drag together represent more than 8% of the total drag and may offer the highest potential for drag reduction. The remaining components, while of no less importances represent only about 2%, of the total drag. Induced drag is one of the major contributors (about 35%) of the total aircraft drag, Figure 1. Various research programs are concentrated on using wing tip devices to increase the effective span of the wing. Integrated winglet concept, with continuous evolution of the wing tip shape has been studied by several research centers worldwide. Among the various types studied, new concepts such as spiroid winglets seems to be the most promising ones [1], which is the focus of the present experimental work as well. 1
Mohammad Reza Soltani, Kaveh Ghorbanian, Mehdi Nazarinia Fig. 1. Drag Breakdown of a Typical Transport Aircraft [1] Finally, while wingspan extensions can reduce the drag of a given airplane, an increase in weight of the original wing is unavoidable. Additional weight penalties can also occur if the extension aggravates flutter for instance. The potential for drag savings has led to several new tip designs appearing on commercial aircraft [2]. 1.1 Technical Discussion Producing induced drag on a finite wing is accompanied by spanwise flow over the wing surface. In particular, the pressure gradients caused by the lower pressures on the upper surface relative to the higher pressures on the lower surface lead to inward spanwise flow on the upper surface and outward spanwise flow on the lower. At the trailing edge, the merging of these two flows with different directions generates a vorticity that is shed from the finite wing and is the origin of induced drag. It has been known for over a century that an endplate at the tip of a finite wing can reduce the spanwise flow and thereby reducing the induced drag. Unfortunately, to be effective, the endplate must be so large that the increase in wetted area drag outweighs any drag reduction. A winglet carries an aerodynamic load which produces a flowfield interacting with that of the main wing thereby reducing the amount of spanwise flow rather than being a simple fence which limits the spanwise flow. In essence, the winglet diffuses the influence of the tip vortex so that the downwash and the induced drags are reduced. In this way, the winglet acts like an endplate in reducing the spanwise flow; however by carrying the proper aerodynamic loading, this is accomplished with much less wetted area. A different way of explanation is that the winglet produces a vertical diffusion of the vorticity in tip region. This diffusion process is also realized as an expansion of the wake in the far field due to the induced velocities from the non-planar components of the winglet. The out of plane bound vortex on the upward winglet induces horizontal velocities on the free wake that cause a spanwise spreading of the wake field [3]. 2 Experimental Apparatus Experiments were conducted in the low speed wind tunnel at Sharif University of Technology (SUT), department of Mechanical Engineering. The tunnel is an open circuit indraft with a rectangular test section of 45 45 16 cm 3. The maximum obtainable speed in the test section is approximately 45 m/sec. A half-span wing with sweep angle of 2 degrees is mounted in the test section. The angle of attack is varied manually from zero to 9 degrees. Figures 2 and 3 show the tunnel and the model as well as different winglet shapes used in this investigation. Fig. 2. SUT wind tunnel 2
EXPERIMENTAL INVESTIGATION OF THE EFFECT OF VARIOUS WINGLET SHAPES ON THE TOTAL PRESSURE DISTRIBUTION BEHIND A WING Fig. 3. The wing model including winglet shapes Four types of winglet shapes, spiroid [4,5] (forward and aft), blended winglets [6], and winggrid [7,8] are used in this investigation. The wing model used in this study has an average chord of 15.7 cm made of fiberglass and manufactured using a numerically controlled milling machine. The Spiroid-typed winglets are made of aluminum and the remaining winglets are made of plexi-glass. Two types of spiroid winglets are used as shown in Figures 4 and 5, Forward (FWD) spiroid winglet and Aft (AFT) spiroid winglet. A 32-tube rake probe is mounted horizontally from the back of the model. A traversing mechanism positions incrementally the probe across the wake. The probe is positioned at the tunnel centerline and the wake survey data are measured at a single traverse plane (at a distance of 2. times average chord) downstream of the model trailing edge. Both surface pressures and wake data are obtained using HCXPM1D6V and HCXPM5D6V pressure transducers with a quoted accuracy of.2% of fullscale pressure range. Each transducer data is collected via a multiplexer and transferred to the computer through a 16-bit Analog-to-Digital (A/D) board. Various sampling rates are performed and finally the best one is selected. 3 Results and Discussions All measurements are conducted at a Reynolds number of 2, using four different types of winglet at various angles of attack. The wake data for all winglets are measured. Each flowfield data set acquired in the tunnel (the traverse plane) comprised of approximately 25 data points. Reduction of the raw pressure data yielded axial velocity component and total as well as static pressure coefficients. The repeatability of the pressure data from the winglet models is also analyzed. For further information the reader is referred to reference [9]. 3.1 Flow Visualization Flow Direction Fig. 4. FWD spiroid winglet model Flow Direction Fig. 5. AFT spiroid winglet model In order to recognize the wing flow pattern, tufts are used. Tufts are small, visible, and light threads attaching in well-positioned places over the wing so that the flow direction is clearly seen. Figures 6a e show one of the positions of flow pattern being photographed. The angle of attack of the wing, for this case is 1 degrees. The flow direction over a typical swept and tapered wing on the upper surface is from tip to root, which is exactly the same as its stall pattern (Fig. 7.). The flow pattern over the model used in this investigation follows the same manner. The stall characteristic is supposed to be similar to those shown in Fig. 7., for the wing without winglet. However, for any kind of installed winglets, the 3
Mohammad Reza Soltani, Kaveh Ghorbanian, Mehdi Nazarinia flow field will change drastically, thus changing the stall characteristics. It should be noted that not all winglets changed the flow field in the same way (Fig. 6.). From this figure it is clearly seen that the flow direction near the tip of wing, as marked in the figure, differs for each winglet, while the angle of attack for all cases are the same. Inspections show that the tip vortex for the cases in Figures 6b and 6c are stronger than other cases and the flow is more complicated. As seen in Fig. 6b and c, the flow direction on the trailing edge (TE) of the wingtip is contrary to what is expected in Fig. 7. The flow direction normally on the upper surface of a swept wing is from tip to root; however, for the wing equipped with Spiroid winglet the opposite is true - that is, from root to tip. However, for the wing grid case, Fig. 6e, the tufts are not affected by the winglet significantly. For the bare wing one, Fig. 6a, it is seen that the flow is almost separated in the tip region. (c) AFT Spiroid (d) Blended (e) Wing Grid Fig. 6. Flow Visualization with Tufts, α= 1 (b) FWD Spiroid 4
EXPERIMENTAL INVESTIGATION OF THE EFFECT OF VARIOUS WINGLET SHAPES ON THE TOTAL PRESSURE DISTRIBUTION BEHIND A WING figure by inspection it is seen that for this angle of attack, the tip vortex for the blended winglet, Fig.9.b, is slightly weaker than other vortices. For other winglets, FWD spiroid, AFT spiroid, and winggrid one the shape and width of the tip vortex differs significantly from those of the bare and blended ones..15.1 Fig. 7. Stall Patterns 3.2 Total Pressure Distribution A 32-tube rake probe mounted horizontally from the back of the model measured the total pressure distribution behind the wing. The measurement nodes comprised of approximately 25 data points, which form a grid shown in Fig. 8. Note that the total pressure distribution is measured in the wake of the model at a distance of x/c=2..5 -.5 5.65.7.75.8.85.9.95.7 86921.5 86917.4 86911.2 8695 86898.8 86892.6 86886.3 8688.1 86873.9 86867.7 86861.5.6.5.4.3.2.1 -.2 P T (Pa) 86919.6 86917.5 86912.5 86911.4 8697.5 8696.3 8692.5 8691.6 86897.5 86892.5 86887.5 Fig. 8. Domain of total pressure measurements behind the wing Figure 9.a-e. shows total pressure distribution behind the wing with four different winglets at zero angle of attack. The pressure distribution for the bare wing (wing without winglet) is also shown as comparison (Fig. 9.a). It can be seen by inspection that at zero angle of attack, the tip vortex for the bare wing is quite clear. The airfoil used in this wing is not symmetric, thus at zero angle of attack it produces some lift which is the result of different pressure distribution over its upper and lower surfaces. As a result of this lift, a tip vortex will be generated as seen by Fig.9.a. The tip vortices for the cases with winglet are different from that of the bare wing, Fig.9.b-e. From this.3.2.1 -.2 1 1.5 (b) Blended Winglet.8 1 1.2 (c) FWD Spiroid 87339.6 87331.3 87323.7 87316 8738.4 873.7 87293.1 87285.4 87277.8 8727.1 87262.5 5
Mohammad Reza Soltani, Kaveh Ghorbanian, Mehdi Nazarinia.2.1.3.2.1 -.2 -.2.8 1 1.2 (d) AFT Spiroid.7.8.9 1 1.1 1.2 87339.6 87333.3 87327 8732.7 87314.4 8738.1 8731.8 87295.5 87289.2 87282.9 87722.1 87714.8 8777.6 877.4 87693.1 87685.9 87678.7 87671.4 87664.2 87657 87649.8 (e) Wing Grid Fig. 9. Effect of different winglet shapes on total pressure distribution behind the wing, α= Effects of angle of attack on the tip vortices formed by those winglets are shown in figures 1 and 11. Again the data for the bare wing are shown for both angle of attack, α=1, 15 for comparison. As the angle of attack increases, the height, width, and strength of the tip vortices varies significantly, figures 1.a and 11.a. As seen from this figure, the tip vortex for the clean model moves inward (toward the root) and the vortex core widens, figures 1.a and 11.a. Figure 1.a shows that the maximum height of the tip vortex core for the clean model at α=1 is about, while for the other two cases shown in Figures 1.c and 1.d, and for the same angle of attack, = and =.12, respectively. Further comparing Figures 1.c and 1.d, it is clear that the FWD spiroid winglet displaces the wake upward, i.e. above the primary tip vortex, while the AFT spiroid case moves the wake below the primary tip vortex. The wake formation behind these two winglets, FWD spiroid and the AFT one as seen from Figure 1, differs from each other significantly. More experiments are needed to better understand these phenomena. The effect of Blended winglet on the total pressure distribution in the wake of the model measured at x/c=2. and α=1 is shown in Figure 1.b. From this figure two relatively weak tip vortices are seen. The total pressure loss caused by these vortices in the wake of the model is again less than that of the bare wing, shown in Figure 1.a. Also, comparing Figures 1.a and 1.b, it could be seen that the height of the tip vortex of the bare wing is lower than that of the blended one. The blended winglet creates two tip vortices, one reaches maximum height of about =.5 and the second one reaches =.45. While the maximum height of the core of the tip vortex for the bare wing is near =, Figure 1.a. Further, it should be mentioned that the first tip vortex of the blended winglet is much weaker than that of the bare wing, Figure 1.a. Comparing Figure 1.a and Figure 1.b, it seems likely that for both cases, the primary tip vortices have almost the same height,, although the blended winglet has displaced the primary tip vortex slightly, Figure 1.b. However, the blended winglet has created another vortex that is located at a height of about =.45. This vortex is probably due to the blended winglet itself, which acts as a half model wing attached at one end to the primary model tip instead of the tunnel wall. Thus it seems likely that this winglet will generate a tip vortex itself when the flow passes over it. However, for the spiroid cases shown in Figure 1, the situation is different. Also note that both and for the tip vortex of the bare wing are much larger than those of the blended winglet ones, i.e. the bare wing has a much larger tip vortex. Further comparing Figures 1.c and 1.d with Figure 1.b, it is clearly seen that the width of the wake generated by the blended winglet is much smaller than those of spiroid ones. Note that all cases shown here are for the same angle 6
EXPERIMENTAL INVESTIGATION OF THE EFFECT OF VARIOUS WINGLET SHAPES ON THE TOTAL PRESSURE DISTRIBUTION BEHIND A WING of attack and the measured total pressure are at the same distance behind the model, x/c=2.. Figure 1.e shows total wake pressure distribution at x/c=2. behind the model equipped with the winggrid. Although the tip vortex formed by this winglet is similar to that of the bare wing one, Figure 1.a, the wake behind the model is different. Further it seems that this type of winglet has displaced the tip vortex too, =.5. Note for the bare wing case, when at the same condition, the height of the vortex core was about. Other data at different angles of attack show similar trends [1]. For higher angle of attack, α=15, Fig.11, the shape of vortices formed behind the wing equipped with winglets vary significantly with the corresponding ones but at lower alpha, Fig.1. It is seen by inspection that the main vortex formed by the blended winglet is widened and has moved inward, Fig.11.b. However, the second vortex formed by this winglet, due to the winglet itself, is much weaker than that of figure.1.b. Furthermore, comparing figure.11.c and 1.c, it is seen that the changes in the tip vortex of the FWD spiroid is much less than that of the AFT one, Fig.11.d and 1.d. The flowfield behind the wing equipped with winggrid also shows some variations with angle of attack. As alpha increases, the tip vortex moves toward the wing root and shifts downward slightly, Fig.11.e and 1.e..1 87898.5 87869 87839.5 8781.1 8778.6 87751.1 87721.7 87692.2 87662.7 87633.3 8763.8 87574.3 87544.9 87515.4.7.6.5.4.3.2.1 -.2 -.3.4.3.2.1 -.2 -.3.4.3.2.1 -.2 -.3 1 1.5 (b) Blended Winglet 87539.2 8753.2 87521.1 87512.1 8753 87494 87484.9 87475.9 87466.9 87457.8 87448.8 8726.9 87195.8 87184.7 87173.6 87162.6 87151.5 8714.4 87129.3 87118.3 8717.2 8796.1 8785 8774.8 1 1.2 1.4 (c) FWD Spiroid.8 1 1.2 1.4 (d) AFT Spiroid 87133.8 87125.4 87117 8718.6 871.3 8791.9 8783.5 8775.1 8766.7 8758.3 875 8741.6 -.2.7.8.9 1 7
Mohammad Reza Soltani, Kaveh Ghorbanian, Mehdi Nazarinia.2.4.1 -.2 87744.2 87731 87717.8 8774.6 87691.4 87678.2 87665 87651.9 87638.7 87625.5 87612.3 87599.1.3.2.1 -.2 87394.4 87388.4 87381.4 87367.9 87354.4 8734.9 87327.4 87313.9 873.4 87286.9 87273.4 87259.9 87246.4 87233 -.3 -.3.8 1 1.2 (e) Wing Grid Fig. 1. Effect of different winglet shapes on total pressure distribution behind the wing, α=1.4.75 1 1.25 1.5 (c) FWD Spiroid.2 8775.9 87746.8 87739.6 87727.7.1 87697.7 87667.7 87637.7 8767.7 87577.7 87547.7 87517.7 87487.7 87457.7 87427.7 87397.7 87367.7 -.2.7.8.9 1 1.1 1.2 1.3.3.2.1 -.2 -.3.3.2.75 1 1.25 1.5 (d) AFT Spiroid 87214 8723.3 87192.6 87181.9 87171.1 8716.4 87149.7 87138.9 87128.2 87117.5 8716.8 8796.6.5.4.3.2.1 -.2 -.3 1 1.5 (b) Blended Winglet 86923.2 8699.1 86895.1 86881 86867 86853 86838.9 86824.9 8681.8 86796.8 86782.8 86768.7 86754.7.1 -.2 -.3 -.4 87481.5 87468 87454.6 87441.1 87427.7 87414.2 874.8 87387.4 87373.9 8736.5.75 1 1.25 1.5 (e) Wing Grid Fig. 11. Effect of different winglet shapes on total pressure distribution behind the wing, α=15 8
EXPERIMENTAL INVESTIGATION OF THE EFFECT OF VARIOUS WINGLET SHAPES ON THE TOTAL PRESSURE DISTRIBUTION BEHIND A WING.3.4.2.1 87454.7 87453.6 87452.4 8745.1 87449 87444.4 87439.8 87435.2 8743.6 87426 87421.4.3.2.1 8733.9 87323.4 87315.8 8738.3 873.8 87293.3 87285.8 87278.3 8727.8 87263.2 87255.7.7.8.9 1 1.1 1.2 -.2.8 1 1.2 1.4 (d) AFT Spiroid.6.2.5.4.3.2.1 86973.6 86967.8 86961.9 86956 8695.2 86944.3 86938.5 86932.6 86926.8.1 87832.2 87825.9 87819.7 87813.4 8787.2 878.9 87794.7 87788.5 -.2.4.75 1 1.25 1.5 (b) Blended Winglet.7.8.9 1 1.1 1.2 1.3 (e) Wing Grid Fig. 12. Effect of different winglet shapes on total pressure distribution behind the wing, α=-5.3.2.1 -.2 8736 87299.9 87293.9 87287.9 87281.9 87275.8 87269.8 87263.8 87257.8 87251.7.8 1 1.2 1.4 (c) FWD Spiroid 3.3 Turbulator Strip The high Reynolds number test corresponds to a typical flight condition. To anchor our data obtained in the wind tunnel, the wing with winglets are tested at a matching low Reynolds number condition with the boundary-layer tripping (forced transition) strategy used in a conventional wind tunnel. Boundary-layer transition strip (Turbulator tape) is placed on the upper surface of the wing. The turbulator strip consists of a.5 cm wide band of sandpaper grains set in a plastic adhesive. The turbulator tape is applied at 1% average chord from the wing leading edge [11]. In order to demonstrate that the flow over 9% of the aft chord of the wing model is changed to turbulent, the Eppler Analysis Code [12] is used to calculate the transition point. The 9
Mohammad Reza Soltani, Kaveh Ghorbanian, Mehdi Nazarinia transition point was calculated to be located at about 1% of the average chord from LE. Figure 13. shows variation of total pressure behind the bare wing and wing with winggrid at α=1. For both cases shown in this figure the model is equipped with turbulator strip installed at x/c=.1, hence creating early flow transition. Comparing the total pressure distribution over the wing surface for this case, figure 13, with the corresponding ones but without turbulator, Fig. 1.a and 1.e, it is seen that the turbulator has not changed the shape and orientation of the vortices significantly. However, the wake behind the model with turbulator is slightly thickened and the total pressure has been reduce, figure 13. Surface pressure data [1] shows similar trend and indicates that the transition has occurred at a place where the turbulator is installed. The total pressure reduction is due to the change of the laminar flow to turbulent, while the difference in flow regime has only quantitative effect. Boundary Layer Turbulator Strip.1 -.2.7.8.9 1 87425.1 87417.6 87414.1 87411.4 8743.3 8741.9 87398.5 87391.6 87377.3 87352.7 87328 8733.4 87278.7 87254.1 87229.4 8724.8 8718.2 87155.5 8713.9 8716.2 8781.6 Boundary Layer Turbulator Strip.2.1 -.2 -.3 P T (Pa) 87277.6 87265.6 87253.7 87241.8 87229.9 87218 8726.1 87194.1 87182.2 8717.3 87158.4 87146.5 87134.6.7.8.9 1 1.1 1.2 1.3 (b) Wing Grid Fig. 13. Effect of Turbulator strip on the wing total pressure distribution, α=1 4 Conclusions The effect of different winglet shapes on the flowfield behind a tapered wing has been experimentally investigated. The measurement included total wake survey behind the wing at various angles of attack with different winglet shapes, flow visualization over the wing surface using tufts, etc. The effect of Reynolds number was investigated by forcing flow transition over the model surface at about x/c.1. The total pressure data in the wake show significant changes for different winglets, but at the same angle of attack. The tip vortices formed by different winglet shapes are totally different from each other and seen that they may have significant effect on the flowfield over the wing surface. These vortices vary differently while the model angle of attack changes. The value of the total pressure measured behind these winglets differs from each other too, indicating their effect on the induced drag. Furthermore, winglets displaced the primary tip vortex and changed its sharp in comparison with that of the bare model. However, these variations differ for each winglet; i.e. different winglets changed the abovementioned characteristics differently. More detailed experimental as well as computational investigations are needed to better understand these phenomena and their effect on the 1
EXPERIMENTAL INVESTIGATION OF THE EFFECT OF VARIOUS WINGLET SHAPES ON THE TOTAL PRESSURE DISTRIBUTION BEHIND A WING overall force and moments generated by the wing when equipped with different winglets. 5 Acknowledgments The present work has been funded in part by the Institute for Aeronautical Research, Bureau of Aircraft Design, HESA. The authors would like to thank the Director of the Design Bureau for his support. References [1] Thiede, P., (Editor), Aerodynamic Drag Reduction Technologies. Proceedings of the CEAS/DragNet European Drag Reduction Conference, 19-21 June 2, Potsdam, Germany, 1 st Ed., Springer, 21. [2] Gold, N. and Visser, K.D. Aerodynamic Effects of Local Dihedral on a Raked Wingtip. 4 th AIAA Aerospace Sciences Meeting & Exhibit, Reno, NV, AIAA-22-831, January 14-17, 22. [3] Maughmer Mark D. The Design of Winglets for High- Performance Sailplanes. AIAA-21-246, 21. [4] Gratzer, Louis, B. Spiroid-Tipped Wing. United States Patent # 5,12,68, April 7, 1992. [5] The Spiroid-Tipped Wing. 21, available on-line; URL: http://www.aviationpartners.com/company/concepts.html [6] Gratzer, Louis, B. Blended Winglet. United States Patent # 5,348,253, September 2, 1994. [7] La Roche. Wing With a Wing Grid As The End Section. United States Patent # 5,823,48, October 2, 1998. [8] La Roche et al. Wing Comprising a Distal Wing Grid. United States Patent # 6,431,499, B1, August 13, 22. [9] Nazarinia, M. Flow Analysis Over and Behind a Wing with Different Winglets. M.Sc. Thesis, Sharif University of Technology, Department of Aerospace Engineering, 23. [1] Soltani M.R., Ghorbanian K. and Nazarinia M. Flow Analysis Over and Behind a Wing with Different Winglet Shapes. 42 nd Aerospace Sciences Meeting & Exhibit, Reno, NV, AIAA 24-723, January 4-8 24. [11] Albert E. von Doenhoff; Elmer A. Horton. A Low- Speed Experimental Investigation of the Effect of a Sandpaper Type of Roughness on Boundary-Layer Transition. NACA Report 1349, 1956, (Supersedes NACA Technical Note 3858). [12] Eppler; Richard. Airfoil Design and Data, Springer- Verlag, 199. 11