Evaluating the Impact of Parasitic Drag on The
Total Page:16
File Type:pdf, Size:1020Kb
1 EVALUATING THE IMPACT OF PARASITIC DRAG ON THE ACCURACY OF ENDURANCE PREDICITIONS Leonid Heide, Chris Regan, Demoz Gebre Egziabher University of Minnesota: Twin Cities August 26, 2018 Abstract The cost-effectiveness and versatility of electric unmanned aerial vehicles (UAVs) has led to an increase in demand for efficient payload-carrying aircraft [4]. Many project criteria require an endurance prediction based on a specific mission profile. Current methods for determining the endurance of fixed-wing UAVs involves estimating the key variable of parasitic drag using empirical equations and approximated drag values. This study seeks to evaluate the impact of estimating parasitic drag using conventional methods on the accuracy of endurance predictions for UAVs. Using glide tests, drag is measured for a small electric aircraft, and compared to predictions made using the component drag build up-method for conventional aircraft. It is demonstrated that the parasitic drag generated by small UAVs is sufficiently large to invalidate endurance predictions based on coefficients of drag calculated with the standard component drag build up method developed for fixed wing aircraft. Background The endurance of electric unmanned aerial vehicle (UAV) can be predicted using the drag and the power system characteristics of an aircraft. Total drag consists of drag caused by the generation of lift 2 (푘퐶퐿 ) and a lift-independent component (퐶퐷0) [6]. This relationship is shown in the expression used to calculate the drag coefficient (퐶퐷): 2 (1) 퐶퐷 = 퐶퐷0 + 푘퐶퐿 . Drag caused by the creation of lift, known as induced drag, can be reliably estimated since the coefficient of lift (퐶퐿) is typically easy to calculate using aircraft characteristics or wind-tunnel data. However the lift-independent component of drag, known as parasitic drag, is difficult to extrapolate analytically as it includes complex factors such as skin friction, pressure drag, and interference drag [2]. 2 In order for an aircraft to maintain steady flight, sufficient thrust must be generated by the power system to overcome the total drag force [1]. In the case of un-accelerated flight thrust must equal drag [6] and the power required is the cross-product of the aircraft’s velocity and drag. In the case of electric aircraft, power is provided by a battery pack. The time over which the battery can reliably output the necessary current can be calculated using the characteristics of the battery, the discharge rate, and the power system’s efficiency (which can be accurately estimated for electric UAVs) [10]. By combining the required power and battery equations, the endurance of an aircraft in hours (E) can be determined by the equation [10]: 푛 1−푛 휂푡표푡푉×퐶 (1) 퐸 = 푅푡 [1 2푊2푘 ] 휌푈3푆퐶 +( ) 2 퐷0 휌푈푆 Equation 1 demonstrates the relationship between the endurance and the coefficient of parasitic drag (퐶퐷0). An accurate drag prediction is therefore essential to predict endurance, and to maximize an aircraft’s flight time and range (since drag must be minimized) [4]. The quickest method to the estimate parasitic drag coefficient is by using aircraft surface area and estimated skin friction coefficients [5]. The skin-friction method is based on the assumption that the parasitic drag for “clean” aircraft is primarily skin-friction drag [7],[3]. However this method only provides a rough estimate of the minimum coefficient of drag. A more accurate drag prediction can be made with the component drag build-up method (CDBM). The CDBM estimates parasitic drag by using modified skin friction coefficients and component “form factors” that estimate the pressure drag due to viscous separation [7], [3]. Effects due to component interference and special features can also be accounted for using the component build-up method. Empirical equations utilized for the CDBM are derived from analyses of full-scale aircraft designs with high Reynolds numbers. The Reynolds number is used as a scaling factor for calculations in aerodynamic, and represents the ratio of the inertial forces to the viscous forces experienced in motion as a body moves through a fluid [2]. The small size of electric UAVs coupled with their lower operational velocities results in them having significantly smaller Reynolds numbers and higher drag relative to their size [9]. The purpose of this study is to demonstrate that estimating drag of a small electric UAV using the component drag build-up method developed for conventional aircraft causes significant inaccuracies in endurance predictions. It is proposed that accurate drag coefficients are necessary to make valid endurance predictions. To demonstrate that parasitic drag on small UAVs is high enough to cause a significant divergence from aircraft endurance predictions, an Eflite Ultra-Stick 25e belonging to the University of Minnesota Uninhabited Aerial Vehicle lab was equipped with an airspeed sensor, GPS, and the Goldy-III 3 flight controller. The aircraft was put through a series of glide tests during which descent rates and true airspeeds were recorded by the Goldy-III flight controller. Using the collected data and standard atmospheric conditions, the glide-angle, lift, and drag equations, the parasitic drag curve was determined for the aircraft. The real drag was then compared to a prediction made using the component drag build up method. It was concluded that the parasitic drag on a small UAV is higher than a prediction made using assumptions that hold for conventional aircraft, thereby invalidating endurance predictions made using coefficients derived from the standard CDBM. Prediction The coefficient of induced drag (퐶퐷푖) was calculated using the equation [6]: 퐶2 (2) 퐶 = 퐿 푑푖 휋퐴푅푒 Where (퐴푅) is the wing’s aspect-ratio. The coefficient of lift (퐶퐿) was calculated using the aircraft’s parameters [6]: 2푊 (3) 퐶퐿 = 2 휌푈 푆푟푒푓 Where (푊) is the aircraft’s weight in, () is density of the free-stream, (푈) is the airspeed, and (푆푟푒푓) is the aircraft’s reference area. Oswald’s efficiency factor (푒) was calculated using the equation provided by Raymer [7]: (4) 푒 = 1.78(1 − 0.0045퐴푅0.68) − 0.64 The coefficient of parasitic drag (퐶퐷0) was calculated using the equation [7]: Σ(푄푐퐶푓퐹퐹푐푆푤푒푡푐) (5) 퐶푑0 = 1.05 [ + 퐶푑0 + 퐶푑0 ] 푆푟푒푓 푝푟표푝푒푙푙푒푟 푙푎푛푑푖푛푔 푔푒푎푟 Where 퐶푓 is the skin friction coefficient, and 퐹퐹푐, 푆푤푒푡푐 represent form factor, and wetted area of a given component respectively. Contributions due to interference drag are represented by 푄푐, which was assumed to be one due to the aircraft being high-wing [7]. The additional 5% accounts for any additional protrusions, such as the pitot-static tube mount. For the wing and the horizontal and vertical stabilizers, the form factor (퐹퐹) was calculated using the following [7]: 4 0.6 푡 푡 0.18 (6) 퐹퐹푤푖푛푔 & 푠푡푎푏푖푙푖푧푒푟푠 = [1 + ( ) + 100 ( ) ] ∗ [1.34푀 ] (푥⁄푐)푚 푐 푐 4 푡 Where (푥⁄푐) is the chord-wise location of the airfoil’s maximum thickness, ( ) is the airfoil thickness 푚 푐 divided by chord, and (푀) is the Mach number. The fuselage’s form factor was calculated using [7]: 60 푓 (7) 퐹퐹 = 1 + + 푓푢푠푒푙푎푔푒 푓3 400 푙 (8) 푓 = 4 √( )퐴 휋 푚푎푥 Where (푙) is the maximum thickness, and (퐴푚푎푥) is the maximum frontal area (assuming a smooth surface). The skin friction coefficient (퐶푓) was calculated using a method from “Boundary Layer Theory” by Hermann Schlichting [8]: 0.455 퐴 (9) 퐶푓 = 2.58 − (푙표푔10푅퐿) 푅퐿 Where the Reynold’s number (푅퐿 ) is based on the wing chord length [2], The constant (A) was determined using a table provided in “Boundary Layer Theory”[8]. A linear fit was assigned to the provided values of (A) and the critical Reynold’s numbers. Variable A was then found for the lower critical Reynold’s number for the NACA 0012 airfoil [11]. To account for a wind-milling prop, an additional factor (퐶푑0푝푟표푝푒푙푙푒푟) was factored on, calculated using: 0.8휎퐴푝푟표푝푒푙푙푒푟 푑푖푠푘 (10) 퐶푑0 = 푝푟표푝푒푙푙푒푟 푆푟푒푓 푛 (11) = 푏푙푎푑푒 퐴푅푏푙푎푑푒휋 Where (푛푏푙푎푑푒) is the number of blades, (퐴푝푟표푝푒푙푙푒푟 푑푖푠푘) is the frontal area of the propeller disk, and (퐴푅푏푙푎푑푒) is the blade aspect ratio. Additionally, the landing gear was accounted for: 0.25퐴푡푖푟푒+0.05퐴푠푡푟푢푡 (12) 퐶푑0 = 1.2 [ ] 푙푎푛푑푖푛푔 푔푒푎푟 푆푟푒푓 Where (퐴푡푖푟푒) and (퐴푠푡푟푢푡) are the frontal cross-sectional areas of the tire and strut respectively. An additional 20% was added to account for interference drag between the components and the fuselage. Using the coefficients, the induced, parasitic, and total drag curves could then be plotted using the drag equation [6]: 1 (13) 퐷 = 푈2푆 퐶 2 푟푒푓 퐷 5 Methodology An Ultra-Stick 25e, a single engine electric aircraft with a 1.25 meter wingspan, was equipped with an airspeed sensor, GPS, and the Goldy-III flight controller and data collector. The aircraft made several controlled glides, where the airspeed was adjusted using through deflection of the elevator. The aircraft altitude and airspeed was recorded by Goldy-III throughout the flight. Using the data, the descent rate (퐻̇ ) and average airspeed (푉퐴푉퐺) was determined over various ranges throughout the glides. Drag was then obtained by first solving for the glide angle (훾) [6]: 퐻̇ (14) 훾 = sin−1 [ ] 푉퐴푉퐺 (15) 퐷 = 푊 sin 훾 Figure 1: “Thor”, the Ultra-Stick 25e used for the glide tests, with the Goldy III flight controller on board 6 Results Figure 2: The bottom graph shows altitude throughout the flight (glides correspond to downward slopes), and the corresponding airspeed on the top graph. Figure 3: A sample range from the third glide test, showing the steady glide descent and the stable velocity 7 Measured Drag 12 11 10 9 8 7 Drag (N) 6 5 4 3 10 12 14 16 18 20 22 24 26 Airspeed (m/s) . Figure 4: The total drag values determined by the glide tests.