JOURNAL OF AIRCRAFT
Solar and Gas Powered Long-Endurance Unmanned Aircraft Sizing via Geometric Programming
Michael Burton∗ and Warren Hoburg† Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 DOI: 10.2514/1.C034405 Fueled by telecommunication needs and opportunities, there has been a recent push to develop aircraft that can provide long-endurance (days to weeks) persistent aerial coverage. These aircraft present a complicated systems engineering problem because of the multifaceted interaction between aerodynamics, structures, environmental effects, and engine, battery, and other component performance. Using geometric programming, models capturing the interaction between disciplines are used to analyze the feasible limits of solar-electric and gas-powered, long- endurance aircraft in seconds to a level of detail and speed not previously achieved in initial aircraft sizing and design. The results show that long-endurance, gas-powered aircraft are generally more robust to higher wind speeds than solar-powered aircraft but are limited in their endurance by the amount of fuel that they can carry. Although solar- electric-powered aircraft can theoretically fly for months, they are operationally limited by reduced solar flux during the winter and wind speeds at higher latitudes. A detailed trade study between gas-powered and solar-powered aircraft is performed to discover which architecture is best suited to meet a given set of requirements and what is the optimum size and endurance of that platform.
Nomenclature I0 = tail boom root moment of inertia A = wing aspect ratio Kq = wing loading constant A = nondimensional cross-sectional area k = tail boom taper index BSFC = brake specific fuel consumption Lh = horizontal tail lift lfuse = fuselage length BSFC100% = brake specific fuel consumption at full throttle b = wing span lh = horizontal tail moment arm C = aircraft drag coefficient lv = vertical tail moment arm D M = wing bending moment Cd0 = nonwing drag coefficient C = skin friction coefficient m = tail boom mass f m 0 = local wing mass CL = lift coefficient mfac = weight margin factor CLmax = maximum lift coefficient c = wing chord mmotor = electric motor mass N = number of flight segments cdp = wing profile drag coefficient c = lift slope coefficient Nmax = safety load factor lα n = number of wing segments Dboom = tail boom drag DOY = day of the year Pavionics = avionics power d = tail boom diameter Poper = aircraft operating power e = span efficiency factor Pshaft = engine/motor shaft power ’ Psun = useable solar power per unit area E = Young s modulus ’ Psun surface = power emitted at the sun s surface Ebatt = energy stored in battery P0 = magnitude of available solar power Esun = total solar energy available ∕ ∕ P S min = minimum operational solar power E S day = daily required solar energy ∕ pwind = percentile wind speed E S sun = available solar energy ∕ R = aircraft range E S twilight = twilight solar energy Rearth orbit = distance from Earth to sun fsolar = planform area fraction covered in solar cells Rfuse = fuselage radius fstructural = fractional structural weight g = gravitational constant Rspec = specific gas constant h = flight altitude Rsun = radius of the sun h = battery specific energy r0 = average distance from Earth to sun batt q = distributed wing loading hcap = spar cap separation _ q = normalized distributed wing loading hmin = minimum climb rate I = cap spar moment of inertia Re = aircraft Reynolds number Reboom = tail boom Reynolds number S = wing planform area S Presented as Paper 2017-4147 at the 18th AIAA/ISSMO Multidisciplinary = shear force Analysis and Optimization Conference, Denver, CO, 5–9 June 2017; received Sfuse = fuselage wetted surface area 2 February 2017; revision received 16 May 2017; accepted for publication 1 Sh = horizontal tail planform area July 2017; published online 25 August 2017. Copyright © 2017 by the Ssolar = solar cell area American Institute of Aeronautics and Astronautics, Inc. All rights reserved. Sv = vertical tail planform area All requests for copying and permission to reprint should be submitted to CCC T = aircraft thrust at www.copyright.com; employ the ISSN 0021-8669 (print) or 1533-3868 t = flight time (online) to initiate your request. See also AIAA Rights and Permissions t = spar cap thickness www.aiaa.org/randp. cap *Ph.D. Candidate, Department of Aeronautics and Astronautics tday = hours of daylight Engineering; [email protected]. Student AIAA. tnight = hours of darkness † Assistant Professor, Department of Aeronautics and Astronautics tsunrise = time of sunrise Engineering; [email protected]. Member AIAA. tsunset = time of sunset Article in Advance / 1 2 Article in Advance / BURTON AND HOBURG
t0 = tail boom root wall thickness wind speed, and other disciplines and environmental models. The V = true airspeed complexity of sizing these kinds of aircraft makes the tradeoffs Vfuse = fuselage volume between gas-powered and solar-powered architectures potentially Vgust = gust velocity nonintuitive. Vh = horizontal tail lift coefficient This paper presents a physics-based optimization model that uses Vref = reference wind airspeed geometric programming as a systematic and rapid approach to Vwind = wind speed evaluate tradeoffs between the two architectures for a given mission. Vv = vertical tail lift coefficient Geometric programming, a form of convex optimization, is chosen as W = aircraft weight a means of evaluating the design space because of its rapid solve time Wave = average weight over a flight segment and guaranteed convergence to a global optimum [5]. Models and Wbatt = battery weight equations representing the interaction of the various disciplines are Wcent = center of aircraft weight expressed in a geometric programming form and are then combined Wemp = empennage weight to form an optimization model. The optimization can be solved in a Wengine = engine weight fraction of a second and is used quantify to the feasible limits of the Wfadd = additional wing structural weight solar-electric and gas-powered long-endurance aircraft. Wfinal = end of flight segment aircraft weight The driving requirement that sizes solar-electric-powered aircraft Wfuel = fuel weight is the ability to operate at multiple locations and during all seasons. To Wfuselage = fuselage weight fly multiple days, solar-electric-powered aircraft must carry enough W h = horizontal tail weight solar cells and batteries to fly during the day while storing enough W = start of flight segment aircraft weight initial energy to fly through the night [6]. If this condition can be achieved W = max takeoff weight MTO during the winter solstice, when the solar flux is at a minimum, then a W = payload weight payload solar-electric-powered aircraft can theoretically fly during any time W = wing skin weight skin of the year for long durations [6]. This becomes more difficult to Wsolar = solar cell weight W = wing spar weight achieve at higher latitudes as solar flux during the winter solstice spar decreases. Additionally, to station keep, the aircraft must fly faster Wstructural = structural weight W = vertical tail weight than the local wind speeds. Wind speeds are a function of latitude, v altitude, and season and tend to increase at higher latitudes and during Wwing = wing weight w = wing deflection winter months. Therefore, a key sizing study of solar-electric- powered aircraft is the effect of latitude on aircraft size. wcap = spar cap width The key driving requirement for a gas-powered aircraft is wmax = maximum deflection limit y = distance from wing root along span endurance. Because gas-powered aircraft are endurance-limited by the amount of fuel that they can carry, longer endurance requires more zbre = Breguet range helper variable α = angle of attack fuel and hence a larger aircraft. Because gas-powered aircraft are not α gust = local angle of attack from gust velocities affected by the solar flux, their station-keeping ability at different Δ = solar declination angle latitudes only depends on the local wind speed. If a gas-powered ΔW = wing section weight aircraft can fly at the latitude with the worst wind speed, it can Δy = wing section length theoretically fly anywhere else in the world. Θ = wing bending deflection angle Using this geometric programming methodology, tradeoffs θ = angle normal to surface between different aircraft configurations, power sources, and θ boom = tail boom deflection angle requirements can be calculated in seconds. The results show that η charge = battery charging efficiency gas-powered aircraft can generally be built lighter and can fly faster ηdischarge = battery discharging efficiency and can therefore fly at higher latitudes and in higher percentile ηmotor = electric motor efficiency wind speeds. Solar-powered aircraft have greater endurance but are η prop = propulsive efficiency limited operationally by their ability to reach higher speeds. The η solar = solar cell efficiency presented optimization methodology can quantify the difference in ρ = air density weight and endurance between and gas and solar-electric-powered ρ A = area density of carbon fiber aircraft for the same set of requirements. Varying mission ρ cfrp cfrp = density of carbon fiber requirements reveals how the best architecture changes depending ρ foam = density of foam on the requirements. ρsolar = solar cell density σ cfrp = maximum carbon-fiber stress τ h = horizontal tail thickness-to-chord ratio τt = wing thickness-to-chord ratio II. Mission Requirements τ v = vertical tail thickness-to-chord ratio A comparison of the gas and solar-electric-powered aircraft is τ w = cap spar width-to-chord ratio achieved by comparing their respective capabilities at satisfying a ϕ = latitude common set of requirements shown in Table 1. These requirements can and will be changed in Sec. VI to observe which architecture is best suited to meet various sets of requirements. A discussion of each I. Introduction requirement is included in this section. ONG-ENDURANCE station-keeping aircraft are potential L solutions to providing internet, communication, and persistent aerial coverage. Recently, technology companies such as Google [1] Table 1 Mission requirements and Facebook [2] have explored the possibility of using solar- electric-powered aircraft to provide internet to parts of the world Requirement Value where 4G or 3G network is not available. Engineering firms like Payload 10 lb Aurora Flight Sciences [3] and Vanilla Aircraft [4] have developed Station keeping 90% winds gas-powered long-endurance platforms as a means of providing Endurance >5 days Season all seasons continual surveillance for days at a time. Sizing of such long- Altitude >4600 m endurance aircraft is complicated because of the multifaceted Latitude 30° interaction between aerodynamics, structural weight, solar energy, Article in Advance / BURTON AND HOBURG 3
A. Station Keeping To station keep, an aircraft must fly at least as fast as the wind speed:
≥ V Vwind (1)
Distributions of wind-speed data [7] indicate that it is impractical for long-endurance unmanned aerial vehicles to station keep in 100 percentile wind speeds. Therefore, the station keeping requirement is parameterized by a percentile wind speed pwind, above which the aircraft is allowed to drift off station. The wind speed at a station is also a function of latitude, altitude, season, or day of the year: