An Extended Methodology for Sizing Solar Unmanned Aerial Vehicles: Theory and Development of a Python Framework for Design Assist

There is a growing interest in using unmanned aerial vehicles (UAVs) in the most diverse application areas from agriculture to remote sensing, that determine the need to project and define mission profiles of the UAVs. In addition, solar photovoltaic energy increases the flight autonomy of this type of aircraft, forming the term Solar UAV. This study proposes an extended methodology for sizing Solar UAVs that take off from a runway. This methodology considers mission parameters such as operating location, altitude, flight speed, flight endurance, and payload to sizing the aircraft parameters, such as wingspan, area of embedded solar cells panels, runway length required for takeoff and landing, battery weight, and the total weight of the aircraft. Using the Python language, we developed a framework to apply the proposed methodology and assist in designing a Solar UAV. With this framework, it was possible to perform a sensitivity analysis of design parameters and constraints. Finally, we performed a simulation of a mission, checking the output parameters.


Introduction
The energy generated by the sun, both as a source of heat and light, is today one of the most promising alternatives to meet the global energy demand. Among the various solar energy applications, we can mention solar power plants, commercial and residential use, transport, rural electrification, solar roads, satellites, and unmanned vehicles.
Unmanned vehicles are motor vehicles that do not have a human operator, are operated autonomously or remotely. These vehicles can carry various loads depending on their type, functionality, characteristics, and mission objectives and profiles [1].
These types of vehicles are classified as unmanned ground vehicles, unmanned underwater vehicles, unmanned surface vehicles (operating on the water's surface), unmanned spacecraft, and unmanned aerial vehicles. Unmanned aerial vehicles (UAVs) present several advantages compared with manned aircraft, mainly in relation to human restrictions. Thus, for missions that can expose the pilot to a high workload or risks, the UAV is highly recommended.
The number of application areas where UAVs can become useful is almost unlimited, and range from military to daily civilian tasks, such as precision agriculture, communications, search and rescue, inspection of facilities and infrastructure, atmospheric satellites, environmental studies, mapping, and continuous border surveillance [1].
These applications have been the subject of interest in several types of research in recent years. In [2][3][4][5], extensive research discusses multiple UAV applications, such as remote sensing (RS), search and rescue (SAR), precision agriculture (PA), infrastructure inspection, road traffic monitoring, and delivery of goods.
Conventional terrestrial communication infrastructure deployment is not workable in some practical scenarios, for example, in dangerous spots, complex terrains, or disaster areas. On the other hand, UAVs equipped with built-in wireless transceivers can be an alternative communication infrastructure because of their high maneuverability and coverage capability [6]. In [7][8][9][10][11], the authors present a detailed survey about UAVs' contribution, benefits, and challenges in wireless communications. In [12,13], the authors present a detailed survey about the use of 5G technology in communication networks based on UAVs.
In communications, satellites are widely used in image capturing missions. However, launching a satellite is not an easy feat, as it requires up-to-date technology and a tremendous investment in time and money. In addition, there is a recent discussion about the growing amount of space debris caused by satellites out of operation in Earth's orbit [14]. For image capturing missions in regions of limited size, the use of a UAV requires less design and implementation time, having a better cost-benefit ratio. The research in [15] compares remote sensing in satellites and UAVs [16].
Using unmanned aerial vehicles in traffic monitoring applications is becoming more and more popular. Compared with traditional monitoring devices, UAVs are more costeffective. In addition, traditional traffic monitoring devices capture traffic conditions at fixed locations; hence, it requires many units to monitor a single road segment. In contrast, a few UAVs can cover a continuous stretch of roadway or even a traffic network.
In natural disasters such as floods, tsunamis, and earthquakes or man-made disasters, such as terrorist attacks, telecommunications systems can be partially or entirely affected. Such situations require rapid solutions to provide communications coverage to support rescue operations [17]. UAVs can provide timely disaster warnings and help speed up rescue and recovery operations when the public communication network is disrupted. In [18], new research was developed using UAVs with pattern recognition and the ability to provide real-time images and advice, for example during forest fires, to the firefighters. Besides that, the UAV can have a router that enables the firefighters to communicate via email, text messaging, and voice mail when the fire destroys the communication infrastructure [18].
There is a growing interest in UAV uses in extensive construction projects monitoring [19] and autonomous navigation for high voltage power lines inspection [20][21][22][23], oil and gas pipelines, wind turbines, and GSM towers infrastructure inspection [17,24].
UAVs can transport food, packages, and other goods [25]. In addition, ambulance drones can deliver medicines, vaccines, and blood samples into and out of unreachable places in the healthcare field. With the rapid demise of snail mail and the massive growth of e-commerce, postal companies have been forced to find new methods to expand beyond their traditional mail delivery business models. Therefore, different postal companies have undertaken various UAV trials to test the feasibility and profitability of UAV delivery services [26].
The mentioned applications require the aircraft to perform long-endurance flights. However, the capabilities of energy sources limit this type of flight. Such constraint has been a subject of interest in various research worldwide, aiming to adopt alternative forms of energy, such as photovoltaic solar energy, giving rise to the term Solar UAV. Some areas of research have been prominent in Solar UAVs, such as improve energy consumption [27] and recharge methods [28], and the development of more efficient components [29]. Solar power as an energy source allows Solar UAVs to perform long-endurance flights, or even perpetual flight, at almost zero emission. Furthermore, Solar UAV is highly efficient because refueling is not required [14].
The operating principle of a Solar UAV is to capture the energy of solar radiation and transform it into electrical energy to supply the propulsive system, actuators, and embedded systems. The surplus electrical energy is stored in a block to supply energy demand when there is no solar radiation, such as on very cloudy days or even in night operations.
The area of the solar cells, from which the energy needed for the aircraft's flight comes, impacts the wing area and the total weight of the aircraft, which in turn affects the amount of energy needed. Understanding and improving the design steps of this type of critical system represents a significant technological advance. Not only for solar-powered UAVs but also for all photovoltaic systems, which can significantly increase overall efficiency and improve its cost-benefit.
Previous studies related to Solar UAVs [30][31][32][33] have performed a proof-of-concept of this type of aircraft. The state-of-the-art solar-powered perpetual flight of low-altitude aircraft was achieved in the UAV developed in [30], based on studies in [31,32]. This UAV achieved a flight of 81 h and 2338 km ground distance. In addition, based on [32,33], the author in [30] developed a Matlab-based conceptual design and analysis framework for solar-powered UAVs.
In [34], the authors proposed a methodology for the conceptual design of Solar UAV, based on [32], and carried out intensive research to verify the project sensitivity for some mission parameters changes from simulations performed on a framework developed in Matlab. The work in [35] applied the methodology proposed in [34], suggesting some simplifications in the input parameters, such as the use of commercially available solar cells and storage blocks. After performing simulations with Solar UAVs, the methodology accuracy was verified using previous authors' aircraft data, with consistent results. A new Matlab framework and a ground model of a Solar UAV were built to validate the methodology proposed in [34].
Recent work in Solar UAVs [36,37] has shown interest in developing conceptual design algorithms, with the possibility of parameter optimization. In [38] the authors proposed an optimization solution using genetic algorithms.
Therefore, this work aims to understand the main parameters in the design of a low altitude long endurance (LALE) Solar UAV, based on [30][31][32], and to extend these methodologies beyond the early contributions proposed in [34,35]. An extended model for sizing Solar UAVs that takeoff and land on a runway, operate in variable endurance, and assemble with commercially available components is proposed. One of the main objectives of the proposed methodology is to simplify the number of parameters provided by the user. Furthermore, this methodology aims to build a Solar UAV with an automatic operating cycle, one that can take off automatically when the batteries reach full charge, necessary for takeoff, cruise flight, and landing. When the battery charge drops to a safety limit, due to the low incidence of solar radiation, the aircraft makes a landing to recharge the batteries. Thus, a new mission cycle begins. Therefore, a detailed study of the energy supply responsible for the aircraft's takeoff and landing is necessary, as well as its impact on the aircraft sizing. The main objectives of this study are summarized as follows: • Develop a complete model for estimating solar radiation on the day of the year with less incidence of sunlight, with no mission site-specific weather data provided by weather stations. The user only needs to inform the operating location, and the algorithm estimates incident solar radiation at that location; • Simplify aerodynamic data input from the model, using an aerodynamic profile database, containing the various aerodynamic parameters necessary for the aircraft sizing; • Develop a model to estimate the power needed for the flight, including the different flight phases (takeoff, climb, cruise descent, and landing). This model must find the optimal relationship between the aircraft's total weight and the power needed to maintain the cruise flight since the highest value of the mass/power ratio is desirable. Such parameters are close to the related mission profile; • Develop a model to estimate the power needed for the flight, including the different flight phases (takeoff, climb, cruise descent, and landing). This model must find the optimal relationship between the total weight of the aircraft and the power needed to carry out the cruise flight since the highest value of the mass/power ratio is desirable.

•
Develop a free open-source framework, with Python language, that is easy to use and that allows the definition of the aircraft route, besides performing aircraft sizing, detailed analysis of design constraints and aircraft optimization.

Solar Unmanned Aerial Vehicles
The operating requirements of the mission determine the characteristics of the UAV, such as shape, size, performance and costs. In addition, other parameters, such as geographic position, aerodynamics, photovoltaic and propulsion characteristics, directly impact the aircraft. From these parameters, the designer creates an initial design of the aircraft as a starting point for a more detailed project or as the first project of iteration [34].

Aerodynamic and Fuselage Settings
The aerodynamic theory used for UAVs is the same as that used for the conceptual design of conventional aircraft [39]. We do not explore the aerodynamic theories deeply, just enough to develop the project focusing on the sizing of the photovoltaic system.
Aircraft wings can take various geometric shapes according to the purpose of the design, but the main types are rectangular, trapezoidal, elliptical, and mixed. Each one has its particular characteristic with advantages and disadvantages when compared with each other.
The rectangular wing has a low aerodynamic efficiency compared with a trapezoidal or elliptical wing. Its advantage is its greater ease of construction and lower manufacturing cost compared with others. In addition, the rectangular shape makes the assembly of solar cells easier.
The planar area of a rectangular wing can be calculated from: where b is the wingspan and c is the profile chord. The aspect ratio (AR) represents the ratio between wingspan and profile chord, as shown in the following equation [40]: (2) Figure 1 shows a rectangular wing and its main geometric features, including its aerodynamic profile, defined by the wing surface in order to obtain an aerodynamic reaction from the flow of the surrounding fluid.
to carry out the cruise flight since the highest value of the mass/power ratio is desirable. • Develop a free open-source framework, with Python language, that is easy to use and that allows the definition of the aircraft route, besides performing aircraft sizing, detailed analysis of design constraints and aircraft optimization.

Solar Unmanned Aerial Vehicles
The operating requirements of the mission determine the characteristics of the UAV, such as shape, size, performance and costs. In addition, other parameters, such as geographic position, aerodynamics, photovoltaic and propulsion characteristics, directly impact the aircraft. From these parameters, the designer creates an initial design of the aircraft as a starting point for a more detailed project or as the first project of iteration [34].

Aerodynamic and Fuselage Settings
The aerodynamic theory used for UAVs is the same as that used for the conceptual design of conventional aircraft [39]. We do not explore the aerodynamic theories deeply, just enough to develop the project focusing on the sizing of the photovoltaic system.
Aircraft wings can take various geometric shapes according to the purpose of the design, but the main types are rectangular, trapezoidal, elliptical, and mixed. Each one has its particular characteristic with advantages and disadvantages when compared with each other.
The rectangular wing has a low aerodynamic efficiency compared with a trapezoidal or elliptical wing. Its advantage is its greater ease of construction and lower manufacturing cost compared with others. In addition, the rectangular shape makes the assembly of solar cells easier.
The planar area of a rectangular wing can be calculated from: where b is the wingspan and is the profile chord. The aspect ratio ( ) represents the ratio between wingspan and profile chord, as shown in the following equation [40]: = .
(2) Figure 1 shows a rectangular wing and its main geometric features, including its aerodynamic profile, defined by the wing surface in order to obtain an aerodynamic reaction from the flow of the surrounding fluid. Each aerodynamic profile has its aerodynamic characteristics, which depend only on its geometric shape, dimensions, bending, thickness, and leading-edge radius. The main aerodynamic characteristics of a profile are the lift, drag, and moment coefficient, the position of the aerodynamic center, and its aerodynamic efficiency [40]. Each aerodynamic profile has its aerodynamic characteristics, which depend only on its geometric shape, dimensions, bending, thickness, and leading-edge radius. The main aerodynamic characteristics of a profile are the lift, drag, and moment coefficient, the position of the aerodynamic center, and its aerodynamic efficiency [40].
The angle of attack α is the term used by aerodynamics to define the angle formed between the profile line chord and the relative wind direction. It represents a parameter that decisively influences the ability to generate profile lift. Usually, the increased angle of attack increases the lift force up to a certain point at which it abruptly decreases. This point is a stall. Increasing the angle of attack also allows the drag to be generated. The dependence of the lift and drag with the angle of attack can be measured by the lift and drag coefficients. Figure 2 shows an aerodynamic profile and its respective angle of attack. between the profile line chord and the relative wind direction. It represents a parameter that decisively influences the ability to generate profile lift. Usually, the increased angle of attack increases the lift force up to a certain point at which it abruptly decreases. This point is a stall. Increasing the angle of attack also allows the drag to be generated. The dependence of the lift and drag with the angle of attack can be measured by the lift and drag coefficients. Figure 2 shows an aerodynamic profile and its respective angle of attack. The lift coefficient represents the efficiency of the profile in generating the lift force, usually determined from a wind tunnel or specific software simulating a wind tunnel. The lift coefficient is a function of the profile model, the Reynolds number, and the angle of attack.
The drag coefficient represents how much drag force the aerodynamic profile produces. The aerodynamic profile is aerodynamically efficient when producing large lift coefficients and small drag coefficients. For an aerodynamic profile, the drag coefficient is also a function of the Reynolds number and the angle of attack.
The forces acting on a wing can be calculated by applying the following equations [40]: where represents the lift force, represents the drag force, represents the air density, is the velocity of the flow, is the wing area, is the drag coefficient and is the lift coefficient.
For a wing of finite dimensions, the drag coefficient is the sum of the profile drag 0 with the induced drag , generated by the wing-tip vortices [40]: The profile drag 0 is experimental and contains all the complex dependencies of the aerodynamic profile.
Mathematically, for a wing with AR ≥ 4, we define the induced drag as [40] where is the wingspan efficiency factor or Oswald efficiency factor. This factor represents a parameter that depends on the geometric model and aspect ratio of the wing. The Oswald efficiency factor can be estimated by [41]: where contains the inviscid part of the induced drag coefficient, and the term is used to express the viscous part of the induced drag coefficient [41]. The lift coefficient represents the efficiency of the profile in generating the lift force, usually determined from a wind tunnel or specific software simulating a wind tunnel. The lift coefficient is a function of the profile model, the Reynolds number, and the angle of attack.
The drag coefficient represents how much drag force the aerodynamic profile produces. The aerodynamic profile is aerodynamically efficient when producing large lift coefficients and small drag coefficients. For an aerodynamic profile, the drag coefficient is also a function of the Reynolds number and the angle of attack.
The forces acting on a wing can be calculated by applying the following equations [40]: where L represents the lift force, D represents the drag force, ρ represents the air density, v is the velocity of the flow, S is the wing area, C D is the drag coefficient and C L is the lift coefficient. For a wing of finite dimensions, the drag coefficient C D is the sum of the profile drag C D0 with the induced drag C Di , generated by the wing-tip vortices [40]: The profile drag C D0 is experimental and contains all the complex dependencies of the aerodynamic profile.
Mathematically, for a wing with AR ≥ 4, we define the induced drag as [40]: where e is the wingspan efficiency factor or Oswald efficiency factor. This factor represents a parameter that depends on the geometric model and aspect ratio of the wing. The Oswald efficiency factor can be estimated by [41]: where Q contains the inviscid part of the induced drag coefficient, and the term P is used to express the viscous part of the induced drag coefficient [41].
Most of the authors provide expressions of the Oswald factor in the form (7) proposed before. Others provide empirical solutions, following wind tunnel data [41].
In [42] the author provides an empirical diagram obtained from flight testing several modern and less modern aircraft. He proposed 1.05 and 0.007 for the values of Q and P respectively, that is:

Energy and Propulsion System
The solar cells, energy storage blocks, propulsion group, payloads, and power electronics define the energy and propulsion systems architecture of a Solar UAV.
The propulsion group comprises the motor, propeller, and motor controller, while the power electronics can use auxiliary systems that optimize power conversions, such as MPPT and DC/DC converters. Finally, we have the payload, actuators, and avionics of the aircraft [43]. Figure 3 shows the basic electrical system of a Solar UAV. posed before. Others provide empirical solutions, following wind tunnel data [41].
In [42] the author provides an empirical diagram obtained from flight testing several modern and less modern aircraft. He proposed 1.05 and 0.007 for the values of and respectively, that is:

Energy and Propulsion System
The solar cells, energy storage blocks, propulsion group, payloads, and power electronics define the energy and propulsion systems architecture of a Solar UAV.
The propulsion group comprises the motor, propeller, and motor controller, while the power electronics can use auxiliary systems that optimize power conversions, such as MPPT and DC/DC converters. Finally, we have the payload, actuators, and avionics of the aircraft [43]. Figure 3 shows the basic electrical system of a Solar UAV.

Sizing of Solar UAV
The sizing process of the aircraft is one of the initial steps from aeronautical design. Normally, it uses historical data and regression analysis from past projects to perform this estimation [44]. Once we create the conceptual design for the aircraft, it must ensure that the estimated weight is sufficiently accurate to meet mission requirements [45].
In designing a Solar UAV, new variables, such as a solar cells panel, make the process more complex and interactive. Furthermore, the addition of a solar cells panel to supply the energy demand directly influences the wing area. Therefore, it directly impacts the total weight of the aircraft and its energy demand.
We should size a Solar UAV according to some design steps: 1. Set mission parameters: payload weight and power demand, flight speed, endurance, mission location, altitude of operation, solar cell and battery model and determine the aerodynamic profile, finding , 0 and e; 2. Adequate survey of the available solar resource at the mission location on the day of the year with less incidence of sunlight; 3. Assume a value for the area of the solar panel, calculate the solar energy collected.
Based on the solar panel area, estimate the wing area and set the wing chord (c) based on the solar cell width. Obtain the wingspan and AR, based on wing chord and wing area. Calculate the lift and drag forces; 4. Assume a value for the total weight of the aircraft; 5. Calculate the required propulsion and total energy demand (power for mission profile) and compare with the value of solar energy collected in 3;

Sizing of Solar UAV
The sizing process of the aircraft is one of the initial steps from aeronautical design. Normally, it uses historical data and regression analysis from past projects to perform this estimation [44]. Once we create the conceptual design for the aircraft, it must ensure that the estimated weight is sufficiently accurate to meet mission requirements [45].
In designing a Solar UAV, new variables, such as a solar cells panel, make the process more complex and interactive. Furthermore, the addition of a solar cells panel to supply the energy demand directly influences the wing area. Therefore, it directly impacts the total weight of the aircraft and its energy demand.
We should size a Solar UAV according to some design steps: 1. Set mission parameters: payload weight and power demand, flight speed, endurance, mission location, altitude of operation, solar cell and battery model and determine the aerodynamic profile, finding C L , C D0 and e; 2.
Adequate survey of the available solar resource at the mission location on the day of the year with less incidence of sunlight; 3.
Assume a value for the area of the solar panel, calculate the solar energy collected. Based on the solar panel area, estimate the wing area and set the wing chord (c) based on the solar cell width. Obtain the wingspan and AR, based on wing chord and wing area. Calculate the lift and drag forces; 4.
Assume a value for the total weight of the aircraft; 5.
Calculate the required propulsion and total energy demand (power for mission profile) and compare with the value of solar energy collected in 3; 6.
Sizing storage block and the total weight of aircraft and compare with the value in 4.
Iteration in steps 3 through 6 may be required several times. For example, in the sizing in step 6, the result for the total weight may require an increase in power for cruise flight since it requires a greater lift to support the heavier weight. Figure 4 shows a flowchart of the UAV sizing. Iteration in steps 3 through 6 may be required several times. For example, in the sizing in step 6, the result for the total weight may require an increase in power for cruise flight since it requires a greater lift to support the heavier weight. Figure 4 shows a flowchart of the UAV sizing.

Mission Parameters (Step 1)
The definition of the operating altitude of the aircraft is one of the fundamental requirements of the mission. From this parameter, we obtain the absolute pressure and air density necessary for the calculation of lift and drag forces on the wing. In addition, the latitude and longitude of the mission are essential to define the solar resource available in the operating region.
The payload is the effective load that an aircraft carrier can perform on a specific mission. The size and mass of the payload and its electric power consumption are the major determinants of the aircraft's layout, size, and overall mass.
The flight endurance of the aircraft may range from, for example, 1 h for a shortrange surveillance system to over 24 h for a long-range surveillance system. The volume and mass of the power storage block determine the endurance. Thus, it matters in the definition of the total mass of the aircraft.
The operational need of the mission determines the flight speed. This speed is a dominant factor in the propulsion configuration of the aircraft.

Mission Parameters (Step 1)
The definition of the operating altitude of the aircraft is one of the fundamental requirements of the mission. From this parameter, we obtain the absolute pressure and air density necessary for the calculation of lift and drag forces on the wing. In addition, the latitude and longitude of the mission are essential to define the solar resource available in the operating region.
The payload is the effective load that an aircraft carrier can perform on a specific mission. The size and mass of the payload and its electric power consumption are the major determinants of the aircraft's layout, size, and overall mass.
The flight endurance of the aircraft may range from, for example, 1 h for a short-range surveillance system to over 24 h for a long-range surveillance system. The volume and mass of the power storage block determine the endurance. Thus, it matters in the definition of the total mass of the aircraft.
The operational need of the mission determines the flight speed. This speed is a dominant factor in the propulsion configuration of the aircraft.
Less weight and greater efficiency of the solar cell is an important aspect of performing aircraft. Furthermore, the addition of rigid solar cells on the aircraft might affect the shape of the UAV wing. Semi-flexible cells, such as the SunPower C60, are ideal, because of their low weight, high efficiency, and ease of assembly.
The specific energy density (Wh/kg) is the main parameter to consider when choosing the power storage block of a UAV, since any increase in weight requires increased propulsion power. Because of its high specific energy density, most of today's solar aircraft, such as Sky-Sailor [32], use lithium-ion batteries, such as the NCR 18650B model.

Solar Resource Assessment (Step 2)
This project stage seeks to quantify the incident global solar irradiation on the photovoltaic panel in a region.
For critical systems, such as those used in UAVs, it is necessary to size the system to operate on the day with the lowest incidence of solar radiation. If the system can operate under these conditions, it can operate on other days of the year.
Some parameters, such as sunshine duration and cloud cover, affect available solar radiation. This radiation is the determining factor in defining the aircraft's wing area. There are several models to estimate the value of solar radiation. Some of them are empirical models which use mathematical formulae. There are also models based on atmospheric characteristics and satellite data such as ozone absorption, and Rayleigh scattering for forecasting solar radiation [46]. It is necessary to point out that the aircraft flight and mission profile will be accomplished according to the design of the aircraft assumptions. To consider general weather conditions for a specific aircraft project, there are necessary tests to feedback the design or match the site conditions to the flight or mission profile. A reasonable solution for general weather is to assume an empirical factor < 1 to multiply the calculated autonomy. The designer or user will obtain more realistic insight after evaluating results.
In this research, we used a model that considers the amount of sunlight available and the possibility of cloud cover. The first parameter found is the declination δ, and the angular position of the sun at the solar noon. To find the declination for a day n of the year, we use [47]: where n is the day of the year, specified as a sequence, ranging from 1 for 1 January to 365 for 31 December. Thus, Equation (8) produces 365 values of declination δ, with a variation in the parameter n. The next parameter is the hourly angle of the sunset ω s , obtained as a function of latitude φ and declination δ and calculated by: From ω s , we obtained the number of hours (N) with the availability of sun, according to [48]: The declination values produce 365 values of ω s , in Equation (9) and N in Equation (10), the smallest of the values of N corresponds to the day of the year with the lowest incidence of solar radiation.
For altitudes less than 2.5 km, we estimated the total average daily extraterrestrial irradiance (H O ), available from the sun using [47]: We estimated the average daily, monthly radiation H on a horizontal surface, considering the effects of the atmosphere, using [48]: where N is the daily monthly average of hours with the availability of sun, i.e., the average duration of the day in the month, calculated as the average of the number of hours N for each month of the year, n is the monthly average daily insolation (hours), available at many hundreds of weather stations in many countries. Figure 5 shows the value of n for the month of June, in Brazil.
We estimated the average daily, monthly radiation � on a horizontal surface, considering the effects of the atmosphere, using [48]: where � is the daily monthly average of hours with the availability of sun, i.e., the average duration of the day in the month, calculated as the average of the number of hours N for each month of the year, � is the monthly average daily insolation (hours), available at many hundreds of weather stations in many countries. Figure 5 shows the value of � for the month of June, in Brazil.

Figure 5.
Monthly average daily insolation for the month of June, in Brazil, source: Adapted from [49].
Here, we propose the use of a data source Global Solar Atlas [50] that considers the daily average global horizontal irradiation to obtain the value of �. We divided the daily average global horizontal irradiation by the nominal solar energy power, i.e., 1000 / 2 , to obtain �. Figure 6 shows the daily average global horizontal irradiation in the city of São Paulo, Brazil. Finally, we estimated the total hourly radiation , per square meter, by [47]: Here, we propose the use of a data source Global Solar Atlas [50] that considers the daily average global horizontal irradiation to obtain the value of n. We divided the daily average global horizontal irradiation by the nominal solar energy power, i.e., 1000 W / m 2 , to obtain n. Figure 6 shows the daily average global horizontal irradiation in the city of São Paulo, Brazil.
We estimated the average daily, monthly radiation � on a horizontal surface, considering the effects of the atmosphere, using [48]: where � is the daily monthly average of hours with the availability of sun, i.e., the average duration of the day in the month, calculated as the average of the number of hours N for each month of the year, � is the monthly average daily insolation (hours), available at many hundreds of weather stations in many countries. Figure 5 shows the value of � for the month of June, in Brazil. Here, we propose the use of a data source Global Solar Atlas [50] that considers the daily average global horizontal irradiation to obtain the value of �. We divided the daily average global horizontal irradiation by the nominal solar energy power, i.e., 1000 / 2 , to obtain �. Figure 6 shows the daily average global horizontal irradiation in the city of São Paulo, Brazil. Finally, we estimated the total hourly radiation , per square meter, by [47]: Finally, we estimated the total hourly radiation I b , per square meter, by [47]: where ω is the hourly angle for each hour of the day (Hour), provided by: Sensors 2021, 21, 7541 10 of 26 for radians, or for degrees using: The constants c and d correspond to [47]: We summed the radiation I b for each hour during the day to find the total solar energy in one day (I T , per square meter).

Solar Energy Collected and Wing Parameters (Step 3)
Initially, we assumed a value for the area of the solar cells panel, and then, we calculated the wing area, considering the total area of the solar cells panel divided by a fill factor, that is: The fill factor is empirical and depends on the solar cell model used. For example, on a solar panel, we consider the effective area of the solar cells panel mounted on that panel.
Solar energy is obtained according to [33,51] by: where I T is the total solar radiation incident on the day of the year with the lowest incidence of solar radiation, A C is the area of the solar cells, η sc and η mppt are the solar cell and maximum power point tracker efficiencies, respectively. The parameter η cbr accounts for solar module level losses mainly caused by the camber of the cell arrangement. These losses are typically in the order of 10% [32].

Aerodynamic Profile
The aerodynamic profile defines lift and drag coefficients. Therefore, a database containing the detailed characteristics of various aerodynamic profiles would be ideal for solar UAV designers. Airfoil tools [52] is an online database that contains detailed information on thousands of aerodynamic profiles. This information can define the values of the wing's lift, drag, and profile drag coefficients.

Power for Mission Profile (Steps 4 and 5)
The first step to determine the total required power for a UAV mission is to define the flight mission segments or mission profile. The mission profile can be divided into takeoff, climb, cruise or level flight, descent, and landing.

Power for Level Flight
In level flight, lift and weight forces are equal, as well as drag and thrust forces. So, we have two equations [53]: The power required for level flight is provided by [30]: where, m tot is the total airplane mass, g is the gravity constant, A wing is the wing area, and ρ(h) is the altitude-dependent air density. Considering the efficiency of the different components of the propulsion system, the total power in level flight is determined by: where η prop includes propeller, gearbox, motor, and motor controller efficiency, P av is the avionics power and P pld is the payload power.

Power for Climb Flight
For the estimation of the power in climb flight, we choose a reasonable climbing angle α cl (between 20 • for hand-launched models and 10 • for aircraft with very large wings) at climbing speed v cl [33].
In cruise flight, the thrust force is the correspondent drag force, but for climb flight, the estimation uses the constant rate of climb. Here, the physical difference concerning level flight is that the thrust not only acts to overcome drag, but on climb flights it also supports a weight component.
Therefore, the required thrust for climbing is provided by [54]: And the net power during the climb flight is: where v cl is climbing speed, provided as a function of stall speed. The stall speed is the minimum speed at which an aircraft can safely lift off from the ground. For safety, it requires the climbing speed to be greater than the stall speed. Considering a value 20% higher, we have [55]: where Cl max is the maximum lift coefficient of the aircraft. It is necessary to multiply the net power during the climb by the time to climb to obtain the power consumed during the climbing flight, that is: We obtained the time to climb from the rate of climb, and the vertical speed of the plane. From [40] comes the expression: where h stands for altitude, T for thrust, D for drag, W for weight, and α cl for the climb angle.
Considering T − D as the excess power, Equation (25) becomes Speed is simply the temporal rate of change in the distance; in this case, the distance is the altitude h, therefore: where dt is a small increment of time required for the small increment dh in altitude. Therefore, the time to climb from an altitude h1 to another h2 is obtained by integrating Equation (27) [54]: For steady climbing flight lift becomes As the altitude increases, L (and therefore C L ) decrease, so the induced drag is also less. As a result, the total drag for the climbing flight is less for the level flight with the same speed [54].
Rearranging Equation (3) to obtain C L From the new value of C L , we calculate C Di from Equation (6), Cd from Equation (5), and D from Equation (4).
With the decrease in the total drag of the aircraft, we have a decrease in excess power, decreasing the R/C according to Equation (26).
We obtained the time to climb from different values of R/C, obtained by the distance traveled between the takeoff altitude and the cruise altitude. To calculate the time of climb graphically, we plotted R/C −1 versus altitude values on a graph, as shown in Figure 7. The area under the curve from takeoff altitude to cruise altitude is the time to climb [54].
where is a small increment of time required for the small increment ℎ in altitude. Therefore, the time to climb from an altitude ℎ1 to another ℎ2 is obtained by integrating Equation (27) [54]: For steady climbing flight lift becomes = cos . (29) As the altitude increases, (and therefore ) decrease, so the induced drag is also less. As a result, the total drag for the climbing flight is less for the level flight with the same speed [54].
From the new value of , we calculate from Equation (6), Cd from Equation (5), and D from Equation (4).
With the decrease in the total drag of the aircraft, we have a decrease in excess power, decreasing the / according to Equation (26).
We obtained the time to climb from different values of R/C, obtained by the distance traveled between the takeoff altitude and the cruise altitude. To calculate the time of climb graphically, we plotted ( / ) −1 versus altitude values on a graph, as shown in Figure 7. The area under the curve from takeoff altitude to cruise altitude is the time to climb [54]. For safety, only the R/C values of the initial and final altitude were plotted, bringing the area of the ( / ) −1 versus altitude graph closer to the area of a trapeze, as shown in the red area of Figure 8. Thus, the area considered in the time to climb estimation is the sum of the blue area and the red area. Finally, the time to climb is multiplied by 2, considering the descending flight as equal to the climbing flight, that is: For safety, only the R/C values of the initial and final altitude were plotted, bringing the area of the (R/C) −1 versus altitude graph closer to the area of a trapeze, as shown in the red area of Figure 8. Thus, the area considered in the time to climb estimation is the sum of the blue area and the red area. Finally, the time to climb is multiplied by 2, considering the descending flight as equal to the climbing flight, that is:  Considering aircraft that take off from a runway, an important parameter is the distance required for takeoff and landing, and the total power required for this stage. The same aircraft may require different runway distances to operate, depending on some parameters, such as air temperature, altitude, wind velocity, and runway declivity. For an aircraft to lift off the ground, the wings must produce a lift force greater than the aircraft's weight. Four variables are essential to produce such force, airspeed, air density, wing area, and lift coefficient [56]. The takeoff and landing distances are provided by [40]: where is the friction coefficient between the aircraft wheels and the ground. This coefficient can vary from 0.02 for paved runways to 0.1 for grass runways. In [40] a list of the friction coefficient with the type of runway material is provided.
As the values of drag and lift forces change as the speed increases, the calculation of Equations (31) and (32) become very complex. As a way of simplifying the solution, the author in [57] suggests using 70% of the takeoff speed for an average flight velocity. Equations (3) and (4) calculate the values of and , considering = 0.7 , therefore: where ϕ is the ground effect, the reduced aerodynamic drag that an aircraft's wings generate when they are close to a fixed surface. Its value is: where ℎ is the height between the wing and the ground, and is the wingspan. The time required to travel the runway during takeoff is and, for landing: Finally, the total time to climb is the sum of the time to climb, time to takeoff, and time to landing, that is: Considering aircraft that take off from a runway, an important parameter is the distance required for takeoff and landing, and the total power required for this stage. The same aircraft may require different runway distances to operate, depending on some parameters, such as air temperature, altitude, wind velocity, and runway declivity. For an aircraft to lift off the ground, the wings must produce a lift force greater than the aircraft's weight. Four variables are essential to produce such force, airspeed, air density, wing area, and lift coefficient [56]. The takeoff and landing distances are provided by [40]: where µ is the friction coefficient between the aircraft wheels and the ground. This coefficient can vary from 0.02 for paved runways to 0.1 for grass runways. In [40] a list of the friction coefficient with the type of runway material is provided.
As the values of drag and lift forces change as the speed increases, the calculation of Equations (31) and (32) become very complex. As a way of simplifying the solution, the author in [57] suggests using 70% of the takeoff speed for an average flight velocity. Equations (3) and (4) calculate the values of L and D, considering v = 0.7v cl , therefore: where φ is the ground effect, the reduced aerodynamic drag that an aircraft's wings generate when they are close to a fixed surface. Its value is: where h is the height between the wing and the ground, and b is the wingspan. The time required to travel the runway during takeoff is and, for landing: Finally, the total time to climb is the sum of the time to climb, time to takeoff, and time to landing, that is: which can now be used in Equation (24).

Weight Prediction (Step 6)
Weight prediction is an essential part of the aircraft design process [44]. It affects the cost and performance characteristics of an aircraft [58]. We can categorize weight prediction into finite element, empirical, and semi-empirical approaches. Finite element weight prediction provides accurate and reliable estimates of airframe mass but they are problem specific; not readily generalized and would increase the computational cost of the process. Empirical methods estimate the mass of the aircraft's main component group, such as a wing or fuselage, using empirical equations that combine geometric parameters, aircraft design speeds, load factor, and statistically derived coefficients. In this approach, we estimated the weight of the individual components of aircraft with the help of curve fitting methods. The sum of these weights is the overall empty weight, or gross takeoff weight. Empirical methods are easier to implement and more efficient than the finite element approach. However, they are of low fidelity. Semi-empirical methods comprise analytically derived equations corrected with the statistical correlation from historical data [59]. The semi-empirical method is the most suitable weight prediction method for the conceptual design synthesis of aircraft. This is because they are easy to implement, as well as sufficiently accurate for the conceptual design phase [60]. Therefore, this is the method chosen for weight prediction in this research.
Weight prediction of Solar UAV starts with the guess of gross takeoff weight, which considers the weights of the structure, propulsion group, solar cells, storage block, and payload, based on the following weight breakdown equation:

Fixed Weights
Some fixed weights do not depend on the sizing of other parts of the aircraft, such as the payload, which is a mission parameter defined at the beginning of the project. Besides the payload, we can consider avionics, such as autopilot, and micro-controlled systems, such as a Raspberry Pi board. Thus, fixed weight is:

Weight Prediction of the Solar Cells
Using commercially available solar cells simplifies the mass prediction of this step, since the manufacturer provides the nominal weight values. The mass of the solar panel is: where Cell area and cell weight are the nominal values of the area and weight of the solar cell used, g is the gravity constant and k is a constant (k > 1) to represent the weight portion of the solar cells panel encapsulation.

Weight Prediction of the Energy Storage Block
We based the weight of the energy storage block on the energy stored for the endurance of the Solar UAV. The energy stored is: where P level and P climb were calculated by Equations (20) and (24), respectively, but now it must consider the efficiency of the energy storage block. T endurance is defined as flight endurance. From Equation (41), we can divide the energy stored in the battery by the specific energy density of the battery (Wh/kg) that is chosen during the definition of the project and multiply by gravity, that is: In [32], a 10% increase in the weight of the energy storage block, propulsion group, and solar cells is proposed due to the weight of the wiring, transformers, controllers, etc.

Weight Prediction of the Aircraft Structure
To estimate the weight of the aircraft structure, the author in [32] proposes a new empirical model, validated using a database containing the parameters of 415 aircraft of various dimensions, divided into 92 radio-controlled unmanned models and 323 manned aircrafts. The following equation is used to estimate the weight of the aircraft structure: where g is the gravity constant, S is the wing area, AR is the aspect ratio. K, x1, and x2 are constants that vary according to the construction quality of the aircraft and where obtained from a table in [32].

Weight Prediction of the Propulsion Group
Aircraft launch conditions have a significant influence on the power to mass ratio of the propulsion group elements. Airplanes that take off from a runway smoothly increase their speed until takeoff, so the difference between the takeoff power and level flight power is low. Aircraft launched manually need to increase their speed and gain altitude quickly. In these cases, the engine must provide a starting power much greater than the power supplied in level flight. This over-sizing is necessary for takeoff, but it also helps to fly in headwinds and turbulence [32].
Therefore, for the contribution of the propulsion group in the total weight of the aircraft, we use the following formula, proposed by [32], based on the over-sizing for handlaunched aircraft, which already considers the contribution of the motor, motor controller, and propeller: W propulsion = 0.008 * P req * g .
The estimated total weight should be compared with the weight assumed initially in the calculation of propulsion power (m tot in Equation (19)), which in this case corresponds to the lift force calculated in Equation (3). If the estimated weight is greater than the assumed weight, a new value is assumed for the initial total weight and the cruise propulsion power is recalculated.

Solar Resource Assessment Model
Many studies related to Solar UAV [32,34,61] use the model proposed by [47] to estimate solar radiation in a given region and period of the year. However, this model uses the following equation, instead of Equation (12), to estimate the monthly average daily solar radiation on a horizontal surface: Equation (45) employs two empirical constants a and b that explain the local climate, such as tropical forest, desert, steppe or semiarid. The work in [47] provides a list of varied climates around the globe, with the values of constants a and b for some regions. In addition, in the equation, the value of n is obtained from meteorological data for the region and is also provided by the user.
As one aim of the proposed methodology is to simplify the number of parameters provided by the user, we employ Equation (12), which eliminates the use of the constants a and b, and we can obtain the value of n at the time of simulation, through a web request from the website global solar atlas. Table 1 shows a comparison between the previous model and the proposed model. We obtained the value of n from [49] and the values of the constants a and b were 0.26 and 0.5, respectively, an approximation of the values provided by [47] for the city of São José do Campos. The day of the year with the lowest incidence of solar irradiation was June 22. As one can see in Table 1 the solar radiation of both models is very close, but the proposed model matches the aim of this extended methodology, to simplify the number of parameters provided by the user.

Parameters Sensitivity and Constraints Analysis
To show the interrelation between the several parameters a multiparameter analysis was performed. Such analysis evidence which designs parameters is more sensitive in the aircraft sizing. The simulations demonstrate the critical parameters of the aircraft mission, identified as the main optimization objectives. To develop the simulations a basic condition with constant parameters was adopted. Table 2 shows the mission parameters that remained fixed during the variation in the other parameters. The proposed algorithm has two main functions: finding the minimum wingspan for constructing the embedded solar panel, producing enough energy to carry out the mission, and optimizing the wingspan to obtain the best relation between aircraft mass and power for cruise flight. Thus, it is desirable that the mass/power ratio is as high as possible. These are the maximum values of the curves, highlighted in these simulations. Figure 9 shows the curves of the mass/power ratio as a function of wingspan for different cruising speeds. Table 3 shows the main parameters of the aircraft for each speed range.  From this analysis, it is possible to verify the great sensitivity of cruise speed in aircraft sizing. This is expected due to the quadratic relationship of speed in the drag force generation (Equation (4)). The higher the speed, the more power the aircraft needs to produce to overcome the drag forces. Thus, this value must be minimized.
The aircraft's mission will usually require a payload with fixed mass; however, it is important to optimize this parameter, since the extra weight in the payload represents greater flexibility in the aircraft's assembly, where additional parts and screws may be needed. Weather characteristics can add weight to the aircraft, such as ice on the wing. Therefore, the payload mass should be maximized. Figure 10 shows the curves of the mass/power ratio as a function of wingspan for different payload mass. Table 4 shows the main parameters of the aircraft for each payload mass range.   From this analysis, it is possible to verify the great sensitivity of cruise speed in aircraft sizing. This is expected due to the quadratic relationship of speed in the drag force generation (Equation (4)). The higher the speed, the more power the aircraft needs to produce to overcome the drag forces. Thus, this value must be minimized.
The aircraft's mission will usually require a payload with fixed mass; however, it is important to optimize this parameter, since the extra weight in the payload represents greater flexibility in the aircraft's assembly, where additional parts and screws may be needed. Weather characteristics can add weight to the aircraft, such as ice on the wing. Therefore, the payload mass should be maximized. Figure 10 shows the curves of the mass/power ratio as a function of wingspan for different payload mass. Table 4 shows the main parameters of the aircraft for each payload mass range.  From this analysis, it is possible to verify the great sensitivity of cruise speed in aircraft sizing. This is expected due to the quadratic relationship of speed in the drag force generation (Equation (4)). The higher the speed, the more power the aircraft needs to produce to overcome the drag forces. Thus, this value must be minimized.
The aircraft's mission will usually require a payload with fixed mass; however, it is important to optimize this parameter, since the extra weight in the payload represents greater flexibility in the aircraft's assembly, where additional parts and screws may be needed. Weather characteristics can add weight to the aircraft, such as ice on the wing. Therefore, the payload mass should be maximized. Figure 10 shows the curves of the mass/power ratio as a function of wingspan for different payload mass. Table 4 shows the main parameters of the aircraft for each payload mass range.   The payload mass values directly impact the total mass of the aircraft, and therefore the cruise power and climb power needed to carry that mass. However, the effect is much smoother than the cruise speed impact.
The longer the aircraft's flight duration, the greater the area it can cover in carrying out the mission before needing to land to recharge its batteries. Therefore, flight endurance must be maximized. Figure 11 shows the curves of the mass/power ratio as a function of wingspan for different flight endurance. Table 5 shows the main parameters of the aircraft for each flight endurance range.  The payload mass values directly impact the total mass of the aircraft, and therefore the cruise power and climb power needed to carry that mass. However, the effect is much smoother than the cruise speed impact.
The longer the aircraft's flight duration, the greater the area it can cover in carrying out the mission before needing to land to recharge its batteries. Therefore, flight endurance must be maximized. Figure 11 shows the curves of the mass/power ratio as a function of wingspan for different flight endurance. Table 5 shows the main parameters of the aircraft for each flight endurance range.  The mass/power ratio curves as a function of wingspan for different flight endurance showed the smallest amplitude compared with the other curves, proving that this factor is the least sensitive in optimization. Figure 12 shows the parameters to be optimized according to design sensitivity.  The mass/power ratio curves as a function of wingspan for different flight endurance showed the smallest amplitude compared with the other curves, proving that this factor is the least sensitive in optimization. Figure 12 shows the parameters to be optimized according to design sensitivity.

Differences between Mission Profiles
In this section, we performed an energy profile and weight comparison between an aircraft with a mission that considers climb power and an aircraft that does not. Figure 13 shows the curves of the mass/power ratio as a function of wingspan for an aircraft with a mission that considers climb power and an aircraft that does not. Table 6 shows the main parameters of each aircraft.  The last parameters analyzed were the total weight and weight of each aircraft component, as a function of the wingspan, for each type of mission profile, with climbing power, and without climbing power, as shown in Figures 14 and 15, respectively.

Differences between Mission Profiles
In this section, we performed an energy profile and weight comparison between an aircraft with a mission that considers climb power and an aircraft that does not. Figure 13 shows the curves of the mass/power ratio as a function of wingspan for an aircraft with a mission that considers climb power and an aircraft that does not. Table 6 shows the main parameters of each aircraft.

Differences between Mission Profiles
In this section, we performed an energy profile and weight comparison between an aircraft with a mission that considers climb power and an aircraft that does not. Figure 13 shows the curves of the mass/power ratio as a function of wingspan for an aircraft with a mission that considers climb power and an aircraft that does not. Table 6 shows the main parameters of each aircraft.  The last parameters analyzed were the total weight and weight of each aircraft component, as a function of the wingspan, for each type of mission profile, with climbing power, and without climbing power, as shown in Figures 14 and 15, respectively.  The last parameters analyzed were the total weight and weight of each aircraft component, as a function of the wingspan, for each type of mission profile, with climbing power, and without climbing power, as shown in Figures 14 and 15, respectively.
From these analyses, it is possible to verify the great impact that the climb power represents on the UAV solar project. For a mission in the same region and with the same characteristics of speed, range, mass, and payload power, the power for the cruise flight represented an increase of about 191% compared with an aircraft that disregards the climb power in the mission profile. This increase is mainly because of the power consumed by the aircraft during the climb and the addition of weight in the battery to store this power.  From these analyses, it is possible to verify the great impact that the climb power represents on the UAV solar project. For a mission in the same region and with the same characteristics of speed, range, mass, and payload power, the power for the cruise flight represented an increase of about 191% compared with an aircraft that disregards the climb power in the mission profile. This increase is mainly because of the power consumed by the aircraft during the climb and the addition of weight in the battery to store this power.

Python Framework
We developed a free opensource framework for sizing Solar UAVs, using the Python language and the open-source libraries Folium [62], Numpy [63], Matplotlib [64], and PyQt [65]. The main aim of this framework is to facilitate the designer's interaction in the main stages of the Solar UAV project, simulating the aircraft's output parameters, such as wingspan, total mass, and propulsion energy, from the mission parameters such as altitude, flying velocity, payload mass, payload power, and latitude. This framework will be available soon, in the following repository [66].
To test the developed Framework, we proposed a surveillance mission along the Paraibuna river, at latitude −23.396648 and longitude −45.598283. Table 7 shows the main parameters of the mission.   From these analyses, it is possible to verify the great impact that the climb power represents on the UAV solar project. For a mission in the same region and with the same characteristics of speed, range, mass, and payload power, the power for the cruise flight represented an increase of about 191% compared with an aircraft that disregards the climb power in the mission profile. This increase is mainly because of the power consumed by the aircraft during the climb and the addition of weight in the battery to store this power.

Python Framework
We developed a free opensource framework for sizing Solar UAVs, using the Python language and the open-source libraries Folium [62], Numpy [63], Matplotlib [64], and PyQt [65]. The main aim of this framework is to facilitate the designer's interaction in the main stages of the Solar UAV project, simulating the aircraft's output parameters, such as wingspan, total mass, and propulsion energy, from the mission parameters such as altitude, flying velocity, payload mass, payload power, and latitude. This framework will be available soon, in the following repository [66].
To test the developed Framework, we proposed a surveillance mission along the Paraibuna river, at latitude −23.396648 and longitude −45.598283. Table 7 shows the main parameters of the mission.

Python Framework
We developed a free opensource framework for sizing Solar UAVs, using the Python language and the open-source libraries Folium [62], Numpy [63], Matplotlib [64], and PyQt [65]. The main aim of this framework is to facilitate the designer's interaction in the main stages of the Solar UAV project, simulating the aircraft's output parameters, such as wingspan, total mass, and propulsion energy, from the mission parameters such as altitude, flying velocity, payload mass, payload power, and latitude. This framework will be available soon, in the following repository [66].
To test the developed Framework, we proposed a surveillance mission along the Paraibuna river, at latitude −23.396648 and longitude −45.598283. Table 7 shows the main parameters of the mission. The first step of the design is the definition of the place of operation of the aircraft, as shown in Figure 16. On the map, the user draws the trajectory of the aircraft. The blue lines represent the route traveled and the red rectangles represent the aircraft's landing and takeoff locations.
The first step of the design is the definition of the place of operation of shown in Figure 16. On the map, the user draws the trajectory of the air lines represent the route traveled and the red rectangles represent the air and takeoff locations. In the next step, the fields for the mission parameters are filled, as sh 17.
Selecting the aerodynamic profile, the solar cell model and the batt framework finds, from a database, the lift and drag coefficient, the densit ciency of solar cell, and density and efficiency of the battery. The user can ponents directly in the database of the framework. The last stage of the framework operation is the sizing itself, as show where it is possible to visualize the aircraft output parameters, as shown in 9 shows the aircraft optimized output parameters. In the next step, the fields for the mission parameters are filled, as shown in Figure 17. Selecting the aerodynamic profile, the solar cell model and the battery model, the framework finds, from a database, the lift and drag coefficient, the density and the efficiency of solar cell, and density and efficiency of the battery. The user can add new components directly in the database of the framework.

Battery Model 18650B
Solar Cell Model Sunpower C60 The first step of the design is the definition of the place of operation of the aircraft, a shown in Figure 16. On the map, the user draws the trajectory of the aircraft. The blu lines represent the route traveled and the red rectangles represent the aircraft's landin and takeoff locations. In the next step, the fields for the mission parameters are filled, as shown in Figur  17.
Selecting the aerodynamic profile, the solar cell model and the battery model, th framework finds, from a database, the lift and drag coefficient, the density and the eff ciency of solar cell, and density and efficiency of the battery. The user can add new com ponents directly in the database of the framework. The last stage of the framework operation is the sizing itself, as shown in Figure 18 where it is possible to visualize the aircraft output parameters, as shown in Table 8. Tabl 9 shows the aircraft optimized output parameters. The last stage of the framework operation is the sizing itself, as shown in Figure 18, where it is possible to visualize the aircraft output parameters, as shown in Table 8. Table 9 shows the aircraft optimized output parameters.

Contributions of the Study
The contributions of this study are summarized as follows: • A complete model for estimating solar radiation on the day of the year with less incidence of sunlight. There is no longer a need for specific data provided by meteorological stations; • The Oswald efficiency factor is now estimated before we should provide this value to estimate the induced drag of the aircraft; • A complete model to estimate the power needed for the flight, including the different flight phases (takeoff, climb, descent and landing), including the type of runway material; • Changes to some details in the weight prediction model to account for commercially available solar cells and batteries; • Regarding the framework, the use of Python language, instead of Matlab, with only open-source libraries, the definition of the aircraft route in an interactive map, detailed analysis of design constraints, and optimal point analysis of sizing measures, varying the mission parameters.
The above contributions are beyond the results achieved early in [34,35]. It is important to highlight that the proposed methodology aiming a mission with limited flight time, instead of a long-range or perpetual flight, as proposed in [30][31][32][33].

Conclusions
After initial considerations and information about Solar UAV designs, an extended methodology for the conceptual design of solar airplanes was presented and compared with original studies.
The proposed methodology enables the evaluation of the role of several parameters to find the optimal relationship between the total weight of the aircraft and the power needed to carry out the cruise and climbing flight. This methodology uses mathematical models that relate the sizing of all aircraft elements, from the concepts of energy and mass balances during the year's day with less incidence of sunlight.
Comparing the solar radiation information data for the proposed methodology with early works of solar models reveals the data are consistent, even with the simplification adopted, which fulfills the aim of the new solar model, to reduce the number of parameters provided by the user.
We performed a multi parameters analysis and showed the interrelation between mass/power ratio, cruise speed, payload mass, and flight endurance. Such analysis evidence, which designs parameters, is more sensitive in the aircraft sizing and obtained the highest sensitivity of the mass/power ratio with the cruise speed.
The methodology enables the authors to develop a framework with a graphical user interface, using Python to simplify user interaction. We show the framework usability and output characteristics after sizing and optimizing a Solar UAV's output parameters from a proposed mission profile and constraints according to the initial discussion and pointed aims.
To demand a small number of entries and data by the user, and taking into account the takeoff, climb, and landing energies, the methodology proves to be adequate for obtaining more realistic results in this type of aircraft, bring new results, and can be considered as a helpful tool for Solar UAV developments.

Conflicts of Interest:
The authors declare no conflict of interest.