Fast Airfoil Selection Methodology for Small Unmanned Aerial Vehicles

Featured Application: Provided results are intended to be used along with the airfoil selection procedure described in the author M. H. Sadraey’s textbook Design of Unmanned Aerial Systems . For more detailed information regarding the procedure, please refer to textbook Aircraft Design: A Systems Engineering Approach of the same author. Abstract: The purpose of this study is to ﬁll the gap that exists when applying the airfoil selection methodology according to the textbooks that appear in the above featured application section, in the low Reynolds number segment, by providing useful data. Data acquisition software was XFLR5. The major result is the construction of a prototype maximum lift coefﬁcient versus ideal lift coefﬁcient diagram, or (cid:0) C l max − C l i (cid:1) diagram, composed exclusively of low Reynolds number airfoils. In addition, the necessary supplementary airfoil characteristics’ tables are provided, for conducting fast airfoil selection for Small Unmanned Aerial Vehicles (SUAVs). As a conclusion by implementing the proposed methodology, the SUAV designer is disengaged from the time-consuming process of the construction of similar C l max − C l i diagrams and supplementary characteristic tables and the airfoil selection-processing time can be greatly shortened, because the main work of the process is reﬂected by the current ﬁndings. To express the time gain in a percentage manner, authors estimate that 85% of engineering time will be economized in the overall airfoil selection procedure if the current ﬁndings are used, due to the fact that no new airfoil simulations are required. Finally, candidate SUAV designers are encouraged to expand the airfoil database, according to the proposed methodology.


Introduction
Unmanned Aerial Systems (UAS), are being used increasingly over the years, compared to manned air vehicles, in order to accomplish certain aerial missions.The main advantage of the unmanned type of flight is the reduced cost of usage [1].Specifically, the Small Unmanned Aerial Vehicles (SUAV) category is of special interest, because it can be used in a plethora of short-range aerial missions, such as in meteorological applications [2], forestry applications [3], agricultural applications [4], or in general purpose security surveillance and reconnaissance applications.The current airfoil selection methodology is a result of the «Euclid» SUAV development, which has a primary role of executing security surveillance and reconnaissance missions.
The UAV's classification can be completed considering various criteria, such as (i) the presence of a main wing or not, (ii) the mass of the vehicle, (iii) the vehicle's maximum altitude, (iv) the vehicle's range and endurance, and (v) the nature of the mission that the vehicle is designed to execute [5].Qualitative and quantitative widths, and also the naming of the UAV's categories, are not strictly defined bibliographically, resulting in a mixed classification between various authors and organizations worldwide [6].For example, the «Small» prefix in front of the UAV acronym, used in the title of the present article, is interpreted as a UAV of mass between 0.1-150 kg according to Australia's Civil Aviation Safety Authority (CASA), and a UAV of mass ≤ 20 kg according to the Civil Aviation Authority (CAA) of Great Britain [7].
In order to be more precise in the current article, the «Small» prefix relates to the Reynolds numbers that the main wing of the vehicle faces.These numbers, of course, are, respectively, in analogy with the physical dimensions and overall performance characteristics of the UAV.
The airfoil, or the different type of airfoils, that will be placed in the main wing of an aircraft, can be either selected among a plethora of well-known geometries, or they can be designed from scratch.Usually, for the purpose of the design of a small aircraft such as a small UAV, the airfoil selection procedure is preferred, in order to preserve resources, as well as to reduce the complexity and development times [8].For general information about the airfoil design procedure, the reader can refer to the textbook [9].For the recent advances in the field, the reader can refer to scientific journal articles [10][11][12].
This study focuses on the airfoil selection means of obtaining an airfoil, which is an absolutely critical procedure for the final performance of the final product [13,14].Different airfoil selection procedures are available when selecting an airfoil for other purposes, such as wind turbines [15][16][17], various propellers [18], special purpose airfoils such as supersonic airfoils [19], unique flying machines such as flapping wings [20], etc.
Figure 1 depicts the contribution of the current study within the broadest context, following the path with the green marking.Before the discussion of the details of the novelty of the current study, all of the possible internal processes of Figure 1 are discussed in detail, along with the necessary recent literature review.
Appl.Sci.2022, 12, 9328 2 of 19 naming of the UAV's categories, are not strictly defined bibliographically, resulting in a mixed classification between various authors and organizations worldwide [6].For example, the «Small» prefix in front of the UAV acronym, used in the title of the present article, is interpreted as a UAV of mass between 0.1-150 kg according to Australia's Civil Aviation Safety Authority (CASA), and a UAV of mass ≤ 20 kg according to the Civil Aviation Authority (CAA) of Great Britain [7].
In order to be more precise in the current article, the «Small» prefix relates to the Reynolds numbers that the main wing of the vehicle faces.These numbers, of course, are, respectively, in analogy with the physical dimensions and overall performance characteristics of the UAV.
The airfoil, or the different type of airfoils, that will be placed in the main wing of an aircraft, can be either selected among a plethora of well-known geometries, or they can be designed from scratch.Usually, for the purpose of the design of a small aircraft such as a small UAV, the airfoil selection procedure is preferred, in order to preserve resources, as well as to reduce the complexity and development times [8].For general information about the airfoil design procedure, the reader can refer to the textbook [9].For the recent advances in the field, the reader can refer to scientific journal articles [10][11][12].
This study focuses on the airfoil selection means of obtaining an airfoil, which is an absolutely critical procedure for the final performance of the final product [13,14].Different airfoil selection procedures are available when selecting an airfoil for other purposes, such as wind turbines [15][16][17], various propellers [18], special purpose airfoils such as supersonic airfoils [19], unique flying machines such as flapping wings [20], etc.
Figure 1 depicts the contribution of the current study within the broadest context, following the path with the green marking.Before the discussion of the details of the novelty of the current study, all of the possible internal processes of Figure 1 are discussed in detail, along with the necessary recent literature review.Possible routes to result with an aircraft airfoil and the exact spot of current's study contribution.Main airfoil selection procedure for this study is based on [8,21].Figure 1.Possible routes to result with an aircraft airfoil and the exact spot of current's study contribution.Main airfoil selection procedure for this study is based on [8,21].
When following the airfoil selection route, the designers have to decide if they are willing to apply a strict airfoil selection methodology which can be obtained from the international aviation literature textbooks, or to attempt applying a custom method.Various Appl.Sci.2022, 12, 9328 3 of 18 well known authors in the aviation design field, describe analytical methodologies in their textbooks [21][22][23].
In the following paragraphs, a literature review is undertaken among recent scientific journal articles where an airfoil selection procedure is taking place.The purpose of this literature review is to collect all of the common attributes from the examined studies, in order to express some useful conclusions for the current study's significance.It should be noted that not under any circumstances are the authors of this study disputing the methodologies used or the results provided in these studies.
In study [24], an airfoil selection procedure for a "Flying Wing Micro Aerial Vehicle" was conducted, where 10 low Reynolds number airfoils were examined.A series of airfoil's aerodynamic measured values, were presented in the tables.In study [25], the airfoil selection for a "Morphing and growing micro Unmanned Air Vehicle" was conducted, without providing enough information for the airfoil measurements or the process itself that could be reused by others.Study [26] discusses a "mini Unmanned Aerial Vehicle" design where three airfoils were subject to selection.Little reusable information is given.Reference [27] examines six NACA airfoils for a "UAV with m tow = 200 kg".No reusable information is given regarding the airfoils.In study [28], six airfoils participated in the selection for a "long-range Vertical Take-Off and Landing (VTOL)" aircraft.An interesting point in this study was the proposal of the "bucket width" characteristic to participate as a variable in the airfoil selection.Our study also gives a great importance to this specific attribute and also displays it in a graphical manner in the C l max − C l i diagram.The importance of the "bucked width" attribute, is discussed later in our study.
Reference [29] provided a satisfactory amount of reusable information for 15 under selection airfoils, for a "Solar-Powered UAV".It also provided a diagram which can be useful for any designer willing to replicate the study's airfoil selection procedure, in other similar projects.This diagram displays C l cruise in the x axis and C l 3/2 /C d in the y axis.This diagram, unfortunately, was not useful for someone applying the base methodology of [8,21].Finally, there were some published articles such as [30,31], where no useful information was given at all regarding the airfoil selection procedure.These studies are written in a manner of an overall vehicle design demonstration, where such details were skipped for practical publishing purposes.
All of the above studies surely accomplished their main target, which was to select the best airfoil for the given circumstances.Some common attributes of the examined studies are the following: 1.
For the majority of the studies, the simulation software used was either Xfoil or XFLR5; 2.
There are studies where useful information cannot be obtained for reuse, due to the small amount of the participating airfoils; 3.
In most cases, studies are not offering well organized and sufficient information regarding airfoil aerodynamic parameters which can be reused in similar new applications.As a result, there are a lot of overlapping simulations for the same airfoils; 4.
None of the above reviewed, indicative studies, but also none in a more extensive literature review that took place before taking the decision of conducting the current study, offers a C l max − C l i diagram along with the necessary supplementary tables for conducting a low Reynolds number airfoil selection, that strictly follows the rules in the literature [8,21].
Continuing with the literature review of the various processes within Figure 1, it can be noticed that after an airfoil is selected or designed, using any of the possible methods, then this airfoil may be a subject of optimization in order to better suit the original requirements.In [32], the authors propose the use of a Genetic Algorithm for the purpose of optimizing UAV airfoils.An application of similar Genetic Algorithms is carried out in [33], where the authors are optimizing the airfoil of an agricultural aircraft.Optimizing processes can also be applied for airfoils of any type, such as wind turbine airfoils [34].
As long as all of the internal details of Figure 1 are now discussed, it is comprehensible that the aircraft's design will result in an airfoil, more or less efficient for the given needs, whatever possible route is followed by the designers.
This study focuses on a certain path, which is marked in green color within Figure 1.Specifically, it concerns selecting an airfoil with the strict application of the methodology described in Dr. Mohammad H. Sadraey's textbook Design of Unmanned Aerial Systems [21].For the readers willing to acquire more information regarding this method, the same author's textbook Aircraft Design: A Systems Engineering Approach [8], is the appropriate source.
This study offers the solution to the following problem: Any high Reynolds number-operating aircraft designer will not face any difficulty applying the aforementioned airfoil selection method, due to the fact that the literature already offers all of the required data to do so.
In reverse, a low Reynolds number aircraft designer, who is willing to select an airfoil for a small UAV, will discover that the data available in the literature are not sufficient.This is one of the reasons that the overlapping simulations of the same airfoils are present in the indicative literature review previously conducted.The designer will eventually be forced to execute a new cycle of airfoil simulations, which will cost in terms of development time.
The novelty of this study consists in the following facts: 1.
It fills the gap that exists in the low Reynolds section when applying the airfoil selection methodology according to [8,21], as it provides the only bibliographically available C l max − C l i diagram until now that presents the following unique characteristics in comparison with the similar diagrams that exists, of course, only for the high Reynolds section.
a.The diagram additionally displays the "bucket width" each studied airfoil has in its C l max − C l i plot.The usefulness of this characteristic is discussed later in this study.b.
The diagram's x and y axes are expressing different Reynolds numbers comparing to each other for increased accuracy, while the x axis corresponds to the cruise phase and the y axis corresponds to the take-off phase, where the vehicle's velocities are different by nature; 2.
It provides supplementary airfoil tables containing the full set of the necessary variables that procedure [8,21] needs, for later use on airfoil comparisons.
The numerical performance gain that a small UAV designer will obtain by using the current study's findings is not related to the performance of the produced UAV, however it is exclusively related to the reduction in the development times because: 1.
If the designer is satisfied with current study's participating airfoils, no additional time-consuming simulations are needed.

2.
If the designer is willing to expand the airfoil database or include other Reynolds numbers, reduced times will also be needed due to the fact that the current study offers its own expansion, step by step, methodology.
To express the time gain in a percentage manner, the authors, after implementing the overall procedures into Euclid Small UAV, estimate that 85% of engineering time will be economized in the above circumstance outlined in 1.The time gain for the above circumstance 2 varies, as it is related to the expansion depth that the designer is willing to apply.

The Euclid Small UAV
«Euclid» SUAV, in which the proposed airfoil selection methodology was implemented for the first time, is the aerial subsystem of the homonymous Unmanned Aerial System (UAS) Euclid.It is a rapid deployment, fixed wing vehicle, for security surveillance and reconnaissance applications, which is being developed by the General (Core) Department of the National and Kapodistrian University of Athens.
Until now, the conceptual design phase and preliminary design phase of the system are completed, and research is entering into the detail design phase of the aerial vehicle.
After stability and control simulations are implemented to the vehicle, valuable stability and control derivatives, along with a plethora of other useful parameters, are imported to the Flight Dynamics Model (FDM) of a flight simulator, in order for the vehicle to be evaluated regarding its handling qualities according to the Cooper and Harper scale.The evaluation results indicated that all of the scores fluctuated in the 1-3 score region.This score region is translated as satisfactory handling qualities, without improvement needed to the system, according to [35].The authors have the ambition to announce the detailed results of this research in the near future and to move on to the construction phase of the vehicle.
Table 1 presents some useful information regarding Euclid Small UAV.In Figure 2, a 3-view of Euclid SUAV is demonstrated.«Euclid» SUAV, in which the proposed airfoil selection methodology was implemented for the first time, is the aerial subsystem of the homonymous Unmanned Aerial System (UAS) Euclid.It is a rapid deployment, fixed wing vehicle, for security surveillance and reconnaissance applications, which is being developed by the General (Core) Department of the National and Kapodistrian University of Athens.
Until now, the conceptual design phase and preliminary design phase of the system are completed, and research is entering into the detail design phase of the aerial vehicle.
After stability and control simulations are implemented to the vehicle, valuable stability and control derivatives, along with a plethora of other useful parameters, are imported to the Flight Dynamics Model (FDM) of a flight simulator, in order for the vehicle to be evaluated regarding its handling qualities according to the Cooper and Harper scale.The evaluation results indicated that all of the scores fluctuated in the 1-3 score region.This score region is translated as satisfactory handling qualities, without improvement needed to the system, according to [35].The authors have the ambition to announce the detailed results of this research in the near future and to move on to the construction phase of the vehicle.
Table 1 presents some useful information regarding Euclid Small UAV.In Figure 2, a 3-view of Euclid SUAV is demonstrated.

Main Airfoil Selection Procedure Description
The author, Sadraey, describes in his book entitled Design of Unmanned Aerial Systems [21], the steps that the UAV designer has to execute, in order to select the optimum airfoil for each case.Detailed information on the procedure can be found at reference [8] of the same author.
In this session, the internal procedures of the main airfoil selection procedure are discussed.The internal steps of the author Sadraey's main selection procedure can be divided into three main groups, as Figure 3 depicts.
The author, Sadraey, describes in his book entitled Design of Unmanned Aerial Systems [21], the steps that the UAV designer has to execute, in order to select the optimum airfoil for each case.Detailed information on the procedure can be found at reference [8] of the same author.
In this session, the internal procedures of the main airfoil selection procedure are discussed.The internal steps of the author Sadraey's main selection procedure can be divided into three main groups, as Figure 3 depicts.

Lift Coefficients Calculation
The two main computational parameters the designer has to know at this point are the average weight ( ), of the aircraft and the main wing's area ().These values are assumed to be already computed during preliminary design.The average weight can be determined using the simple Equation (1): where  represents the initial weight of the aircraft before cruise and  represents the final weight of the aircraft after the cruise.If the UAV is electric powered, which is the most possible situation depending on the Reynolds numbers the article refers to, then  =  =  .
All of the required parameters regarding the various lift coefficients, which are used in the procedure, can be calculated by applying the equations of the following Table 2.

Lift Coefficients Calculation
The two main computational parameters the designer has to know at this point are the average weight (W avg ), of the aircraft and the main wing's area (S).These values are assumed to be already computed during preliminary design.The average weight can be determined using the simple Equation (1): where W i represents the initial weight of the aircraft before cruise and W f represents the final weight of the aircraft after the cruise.If the UAV is electric powered, which is the most possible situation depending on the Reynolds numbers the article refers to, then All of the required parameters regarding the various lift coefficients, which are used in the procedure, can be calculated by applying the equations of the following Table 2.

Use of Airfoil C l max − C l i Diagrams and Supplementary Tables
The next step, after the designer finishes the lift coefficients' calculations, is to use corresponding airfoil diagrams in order to focus on a certain region of the airfoils that satisfy these values.These diagrams are referred to as the maximum lift coefficient versus ideal lift coefficient diagrams or C l max − C l i diagrams.Such diagrams and relative information regarding the high Reynolds number airfoil section are given in [8,21], but no such organized information are available for the low Reynolds number section, either in [8,21], or in the broadest aircraft design literature.In the Materials and Methods section, we discuss the steps that were taken to collect this information.In the Results section, we present these reusable findings.

Airfoil Comparison
Having the valuable C lmax − C li diagram in hand, the designer can focus on an area of interest, according to the calculated C lmax and C li values.If there are a lot of airfoils found to be suitable, then the designer has to organize a «weight table», a type of comparison table, for further airfoil processing based on supplementary airfoil characteristics.
The airfoil characteristics to be evaluated, and also the weight of each characteristic, are subject to the discretion of the designer.An indicative/suggestive characteristics' weight composition for the low Reynolds number section is presented in the following Table 3.The explanation of each characteristic variable can be found in Section 5.1.Total: 100 In the case of large aircraft design, other airfoil characteristics may contribute to the selection process, such as manufacturability, room adequacy for the fuel tanks, etc.In this case, such characteristics are not participating because all of the twenty examined airfoils are easily manufacturable and there is no need for a large space inside the wing of a small UAV.The airfoil to be finally chosen is the one which will concentrate the best overall score.

Materials and Methods
In this section we discuss the steps that were taken in order to collect all of the critical values for the variables needed and for the construction of the C lmax − C li diagram for the low Reynolds number airfoils.Figure 4 visualizes these steps.
UAV.The airfoil to be finally chosen is the one which will concentrate the best overall score.

Materials and Methods
In this section we discuss the steps that were taken in order to collect all of the critical values for the variables needed and for the construction of the Clmax − Cli diagram for the low Reynolds number airfoils.Figure 4 visualizes these steps.A SUAV designer who is not satisfied with the airfoils participating in this study, or a designer who wants to contribute and expand the current airfoil database, has to repeat Figure 4 steps 1-6 for each new airfoil.
Data acquisition for the proposed methodology is based on XFLR5 v6.48 software.XFLR5 is a free wing and aircraft computational software operating at low Reynolds numbers.Design and analysis are based on Lifting Line Theory, Vortex Lattice Method, and 3D Panel Method [36].XFLR5 is based on Xfoil, a software inspired by Prof. Mark Drela.The validity in the XFLR5 results is scientifically documented since the Xfoil era [37], and continues to be verified by the scientific community until today, by comparing the software results to those acquired by computational fluid dynamics (CFD) simulations and wind tunnel tests [24,38].
The starting point of an airfoil selection is the Reynolds number determination.This number can be determined using the following Equation ( 2): describes the fluid velocity,  the characteristic length of the studied arrangement which in this case is the length of the airfoil's chord, and  is the kinematic viscosity of the fluid.In this article, a typical chord length of  = 0.2 m is considered, along with the three typical SUAV speeds appearing in Table 4.The kinematic viscosity of the fluid (air), is set at 1.511 ⋅ 10 m s , which corresponds to an ambient temperature of 20 °C.A SUAV designer who is not satisfied with the airfoils participating in this study, or a designer who wants to contribute and expand the current airfoil database, has to repeat Figure 4 steps 1-6 for each new airfoil.
Data acquisition for the proposed methodology is based on XFLR5 v6.48 software.XFLR5 is a free wing and aircraft computational software operating at low Reynolds numbers.Design and analysis are based on Lifting Line Theory, Vortex Lattice Method, and 3D Panel Method [36].XFLR5 is based on Xfoil, a software inspired by Prof. Mark Drela.The validity in the XFLR5 results is scientifically documented since the Xfoil era [37], and continues to be verified by the scientific community until today, by comparing the software results to those acquired by computational fluid dynamics (CFD) simulations and wind tunnel tests [24,38].
The starting point of an airfoil selection is the Reynolds number determination.This number can be determined using the following Equation ( 2): V describes the fluid velocity, l the characteristic length of the studied arrangement which in this case is the length of the airfoil's chord, and ν is the kinematic viscosity of the fluid.In this article, a typical chord length of l = 0.2 m is considered, along with the three typical SUAV speeds appearing in Table 4.The kinematic viscosity of the fluid (air), is set at 1.511 • 10 −5 m 2 /s, which corresponds to an ambient temperature of 20 • C. Setting the precise calculated Reynolds number in the simulation software is not critical due to the fact that this number can change slightly between different missions.This is because, for example, the UAV operator may choose to fly at a slightly different speed from the recommended, or because of a different mass payload is attached.For these reasons, Table 4 contains the Reynolds numbers in increments of 50,000.

Literature Review in Low Reynolds Number Airfoils for Small UAV's
An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.

.dat File Construction with the Polars of Each Airfoil
Certain airfoil families, such as NACA, have the advantage that some of their properties are known only by reading the numbers of their name.XFLR5 can also take advantage of this convention and display an airfoil's coordinates only by entering their name.For the special purpose airfoils studied in this article, such an automation is not present because there is not a correlation between the names and physical properties of the airfoils.In order to display them correctly in the software, the user must create a .datfile containing the numerical polars of each one.Some popular airfoil's .datfiles can be found on the Internet.In the case of absence, the .datfiles can be created manually using the values of the polar tables in the corresponding source reference (all twenty participating airfoil's .datfiles can be found in .datformat under Supplementary Materials File S1, for methodology repeatability and verification purposes).

Construction of a Clone File of Each Airfoil and of Trailing Edge Flap
Wishing to study the effect of the deflection of a trailing edge flap for each airfoil, a clone of each airfoil was created and a trailing edge flap was integrated at 80% of the cord, with a 30 • down deflection angle.

Defining XFLR5 Simulation Parameters and Creation of .xfl Files
Multithreaded Batch Analysis was the selected type of simulation throughout the process.Three discrete Reynolds numbers were used (100,000, 200,000, 350,000) for the airfoils without flap deflection.The airfoils with a deflected flap were simulated at Re = 100,000.The angle of attack was studied within −6 • + 25 • range, with a 0.5 • step.

XFLR5 Simulations
Twenty simulations were carried out, examining a corresponding number of airfoil pairs (with and without flaps).Generated diagrams were C d -C l , C l -α, C m -α and (C d /C l )-α (all twenty participating airfoil's XFLR5 simulation files can be found in .xflformat under Supplementary Materials File S2, for methodology repeatability and verification purposes).An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.from the recommended, or because of a different mass payload is attached.For these reasons, Table 4 contains the Reynolds numbers in increments of 50,000.

Literature Review in Low Reynolds Number Airfoils for Small UAV's
An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.from the recommended, or because of a different mass payload is attached.For these reasons, Table 4 contains the Reynolds numbers in increments of 50,000.

Literature Review in Low Reynolds Number Airfoils for Small UAV's
An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.from the recommended, or because of a different mass payload is attached.For these reasons, Table 4 contains the Reynolds numbers in increments of 50,000.

Literature Review in Low Reynolds Number Airfoils for Small UAV's
An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.Low Reynolds is because, for example, the UAV operator may choose to fly at a slightly different speed from the recommended, or because of a different mass payload is attached.For these reasons, Table 4 contains the Reynolds numbers in increments of 50,000.

Literature Review in Low Reynolds Number Airfoils for Small UAV's
An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.

E205
Richard Eppler [33] Low Reynolds is because, for example, the UAV operator may choose to fly at a slightly different speed from the recommended, or because of a different mass payload is attached.For these reasons, Table 4 contains the Reynolds numbers in increments of 50,000.

Literature Review in Low Reynolds Number Airfoils for Small UAV's
An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.

MH114
Martin Hepperle [41] Ultralight Airfoil is because, for example, the UAV operator may choose to fly at a slightly different speed from the recommended, or because of a different mass payload is attached.For these reasons, Table 4 contains the Reynolds numbers in increments of 50,000.

Literature Review in Low Reynolds Number Airfoils for Small UAV's
An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.RC Sailplane Airfoil ical due to the fact that this number can change slightly between different missions.This is because, for example, the UAV operator may choose to fly at a slightly different speed from the recommended, or because of a different mass payload is attached.For these reasons, Table 4 contains the Reynolds numbers in increments of 50,000.

Literature Review in Low Reynolds Number Airfoils for Small UAV's
An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.

E387
Richard Eppler [38] Low Reynolds ical due to the fact that this number can change slightly between different missions.This is because, for example, the UAV operator may choose to fly at a slightly different speed from the recommended, or because of a different mass payload is attached.For these reasons, Table 4 contains the Reynolds numbers in increments of 50,000.

Literature Review in Low Reynolds Number Airfoils for Small UAV's
An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.Hand Launch RC Airfoil ical due to the fact that this number can change slightly between different missions.This is because, for example, the UAV operator may choose to fly at a slightly different speed from the recommended, or because of a different mass payload is attached.For these reasons, Table 4 contains the Reynolds numbers in increments of 50,000.

Literature Review in Low Reynolds Number Airfoils for Small UAV's
An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.Setting the precise calculated Reynolds number in the simulation software is not crit-ical due to the fact that this number can change slightly between different missions.This is because, for example, the UAV operator may choose to fly at a slightly different speed from the recommended, or because of a different mass payload is attached.For these reasons, Table 4 contains the Reynolds numbers in increments of 50,000.

Literature Review in Low Reynolds Number Airfoils for Small UAV's
An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.

S4083
Michael Selig [45] Hand Launch RC Airfoil Setting the precise calculated Reynolds number in the simulation software is not critical due to the fact that this number can change slightly between different missions.This is because, for example, the UAV operator may choose to fly at a slightly different speed from the recommended, or because of a different mass payload is attached.For these reasons, Table 4 contains the Reynolds numbers in increments of 50,000.

Literature Review in Low Reynolds Number Airfoils for Small UAV's
An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.Setting the precise calculated Reynolds number in the simulation software is not critical due to the fact that this number can change slightly between different missions.This is because, for example, the UAV operator may choose to fly at a slightly different speed from the recommended, or because of a different mass payload is attached.For these reasons, Table 4 contains the Reynolds numbers in increments of 50,000.

Literature Review in Low Reynolds Number Airfoils for Small UAV's
An extensive literature review took place, aiming to locate airfoil families or discrete airfoils, whose effectiveness is proven in the SUAV field.From the literature review, twelve references qualified for the next step.

Selection of Adequate Number of Airfoils for Later Processing
Based on the qualified literature, twenty discrete airfoils were chosen to participate.Table 5 presents those airfoils.

.dat File Construction with the Polars of Each Airfoil
Certain airfoil families, such as NACA, have the advantage that some of their properties are known only by reading the numbers of their name.XFLR5 can also take advantage of this convention and display an airfoil's coordinates only by entering their name.For the special purpose airfoils studied in this article, such an automation is not present because there is not a correlation between the names and physical properties of

Results
Figure 5 contains the three steps that were taken for the results to be prepared for presentation.As stated at the beginning of this work, these steps represent the main results from the current study.These findings are sufficient to fast conduct an airfoil selection among the participating airfoils.
airfoils without flap deflection.The airfoils with a deflected flap were simulated at Re = 100,000.The angle of attack was studied within −6° + 25° range, with a 0.5° step.

XFLR5 Simulations
Twenty simulations were carried out, examining a corresponding number of airfoil pairs (with and without flaps).Generated diagrams were Cd-Cl, Cl-α, Cm-α and (Cd/Cl)-α (all twenty participating airfoil's XFLR5 simulation files can be found in .xflformat under Supplementary Materials File S2, for methodology repeatability and verification purposes).

Results
Figure 5 contains the three steps that were taken for the results to be prepared for presentation.As stated at the beginning of this work, these steps represent the main results from the current study.These findings are sufficient to fast conduct an airfoil selection among the participating airfoils.

Results Export from XFLR5 and Import to Spreadsheet Software
Unfortunately, XFLR5 does not provide any automated routine for user-defined result exporting.Despite this fact, the user can acquire the precise numerical pairs of any given point on a plot curve.Given this, all of the measurements of interest were transferred manually to spreadsheets.The following Tables 6-9, were constructed in order to be used by the SUAV designer at the airfoil comparison procedure, if needed (all twenty participating airfoil's aerodynamic measurements can be found in .xlsxformat under Supplementary Materials Table S2, for methodology repeatability and verification purposes).The measured attributes of the above Tables 6-9 are explained in detail in the following numbered paragraphs, according to [8,21]  «bucket width» C lbucket_width : the authors thought that it could be useful to conduct a series of measurements under the C lbucket_width name.It is observed that regarding the classic airfoils, such as the NACA series, the bucket always has a smaller width compared to the low Reynolds number airfoils, where the bucket is more spread out in the x axis.An example is demonstrated in Figure 6, where the SG6043 low Reynolds number airfoil is simulated at Re = 350,000, contrary to a classic NACA 66 2 -415 at Re = 3,000,000.The low Reynolds number-widened bucket, is an excellent advantage both for the designer and the operator of the SUAV, as the possible suitable airfoils are increased for the designer and the deviation margin in the flight envelope is also bigger for the operator, without a big penalty in the range and the endurance of the vehicle.(Units: dimensionless);

Airfoil Stall Quality Judgement
All of the above variables are quantitative except for the stall quality variable.Quantitative values can be inserted with absolute accuracy to Tables 6-9.A way to quantify the stall quality is described in Figure 7.

Airfoil Stall Quality Judgement
All of the above variables are quantitative except for the stall quality variable.Quantitative values can be inserted with absolute accuracy to Tables 6-9.A way to quantify the stall quality is described in Figure 7.

Airfoil Stall Quality Judgement
All of the above variables are quantitative except for the stall quality variable.Quantitative values can be inserted with absolute accuracy to Tables 6-9.A way to quantify the stall quality is described in Figure 7.    Obviously, the preference is green > orange > red.We consider a vertical line at distance ∆x from the α s point that intersects with the three stall curves.Afterwards, in the magnified right part of the diagram, we record ∆ y1 , ∆ y2 , ∆ y3 points.The curve which records the smallest ∆y is the one to receive the best score to contribute to its selection process.
Unfortunately, things are not quite so easy, because with the aforementioned method we are studying ideal curves that are unlikely to exist in real simulation results because: 1.
Real under study points are close to each other so inaccuracy will appear in the measurements; 2.
Real simulation curves sometimes exhibit secondary or/and tertiary fluctuations (for example, while a curve seems to fade steeply, then it is being regularized, or vice versa. For these reasons, a completely different method was chosen for the stall quality judgement.The proposed method faces the stall curves in a holistic way and is based on the blind judgement of three individual reviewers.A simple Likert scale was formed, where score 1 corresponded to the worst stall quality, while score 5 corresponded to excellent stall quality behavior.For each Reynolds number, an average value was calculated from the three reviewers' opinion.Eventually, a final stall quality score was calculated for each airfoil, using Table 10 weights.The final stall quality scores that each airfoil received are presented in Table 11.The detailed judgement of each reviewer along with the stall curves diagrams, can be found in .xlsxformat in Table S1.The average score value is 2.25.The main differentiation that the proposed C lmax − C li diagram (Figure 8; Table S2 in .xlsxformat) introduces over against the corresponding high Reynolds number diagrams that can be found in the literature, consists of the following statements: 1.
It introduces two different Reynolds numbers across the two axes, instead of one.C lmax is calculated at Re = 100,000, while C li at Re = 200,000; 2.
Aside of the points that represent each airfoil, a faded line appears to the left and/or right of each point.This line represents the bucket width variable that we discussed earlier.S2).
1.The big bucket widths are now visually verified; 2. The Cli points' position is not mandatory to be in the middle of the bucket.The exact position has to do with the physical shape that the bucket has each time; 3. The four airfoils (S1123, FX63-137, MH114 and SG6043) that appear at the upper part  S2).

1.
The big bucket widths are now visually verified; 2.
The C li points' position is not mandatory to be in the middle of the bucket.The exact position has to do with the physical shape that the bucket has each time; 3.
The four airfoils (S1123, FX63-137, MH114 and SG6043) that appear at the upper part of the diagram (those with high C lmax ), are indeed characterized by their designers as High Lift, Human power aircraft airfoil, Ultralight Airfoil, and High L/D, respectively (see Table 5).All of these airfoils are of high camber.In combination with Table 9, ∆C Lf measurement, it can be stated that flap deflection does not have a great effect on the lift coefficient improvement of those high camber airfoils.The S8064 airfoil which has the lowest C lmax is a counterexample, as it is characterized by its designer as a fast RC aircraft; 4.
C li dispersion is located around the 0.5-0.6 region and C lmax dispersion is located around the 0.9-1 region.

Conclusions
This article's findings disengage the SUAV designer from the time-consuming process of the construction of C l max − C l i diagrams and supplementary characteristic tables when applying the main airfoil selection procedure [8,21] regarding a low Reynolds number aerial vehicle.To express the time gain in a percentage manner, the authors, after implementing the overall procedures into Euclid Small UAV, estimate that 85% of engineering time will be economized using the proposed C l max − C l i diagram and the proposed supplementary characteristic tables.The work left for the SUAV designer is to focus on a certain region of the map and to compose an airfoil characteristics' comparison table, in order to select the best airfoil for the given application.
The validity of the proposed methodology is reinforced with the successful use of it in the Euclid SUAV project.Finally, candidate SUAV designers are encouraged to expand the airfoil database, according to the proposed expansion methodology discussed in Section 4.

Figure 1 .
Figure 1.Possible routes to result with an aircraft airfoil and the exact spot of current's study contribution.Main airfoil selection procedure for this study is based on [8,21].Figure 1. Possible routes to result with an aircraft airfoil and the exact spot of current's study contribution.Main airfoil selection procedure for this study is based on [8,21].

Figure 3 .
Figure 3.The internal procedures of the author Sadraey's main airfoil selection procedure, divided in three main groups.

Figure 3 .
Figure 3.The internal procedures of the author Sadraey's main airfoil selection procedure, divided in three main groups.

Figure 4 .
Figure 4.The steps that were taken for the airfoil data collection.

Figure 4 .
Figure 4.The steps that were taken for the airfoil data collection.

Figure 5 .
Figure 5.The steps were taken for results' organization.

Figure 5 .
Figure 5.The steps were taken for results' organization. :

1 .
Minimum drag coefficient C dmin : Is the minimum C d value in the C d -C l diagram (Units: dimensionless); 2.

3 .
Ideal lift coefficient C li : It corresponds to the minimum drag coefficient.When studying an ideal airfoil with a symmetric bucket, C li should be in the middle of the bucket.In the real world, C li can be anywhere close to the curve that we consider as the bucket.(Units: dimensionless); 4.Ideal angle of attack αi : Is the angle of attack that corresponds to C li .(Units: degrees); 5. Pitching moment coefficient C m : The Variations of Pitching Moment Coefficient versus Angle of Attack.The design objective is to have C m close to zero as far as possible.(Units: 1/rad); 6. (C l /C d ) max : Is the maximum value of the C d /C l -α diagram.Large values are desirable.(Units: dimensionless); 7. Stall angle α s : Is the critical angle of attack beyond which the airfoil will suffer a decrease in lift force.Large values are desirable, especially at lower Reynolds numbers and when using flaps.(Units: degrees); 8. Airfoil gross maximum lift coefficient C lmaxg : It corresponds to the maximum C l value on a C l -α diagram.(Units: dimensionless); 9. Stall quality: It is a prediction of what is going to happen in the uncomfortable situation where the airfoil will enter the stall region.More details on the stall quality judgment procedure can be found in the next Section 5.2.(Units: Likert 1-5 scale, where: 1 corresponds to the worst stall quality, while score 5 corresponds to excellent stall quality behavior); 10.Flap deflection lift coefficient difference ∆C Lf : It exclusively relates to the flap deflection data (Table 9).It is synonymous with ∆C lHLD.It measures the difference in lift coefficients between the flap and absence of flap deflection.(Units: dimensionless); 11.C lmax : Is the maximum lift coefficient, calculated by C lmax = (C lmaxg − ∆ Lf ). (Units: dimensionless).Appl.Sci.2022, 12, 9328 14 of 19 11.Clmax: Is the maximum lift coefficient, calculated by Clmax = (Clmaxg − ΔLf).(Units: dimensionless).

Figure 6 .
Figure 6.An example of the outspreaded bucket of a low Reynolds number airfoil (SG6043) operating at Re = 350,000, in contrast with the narrow bucket of a NACA 662-415 operating at Re = 3,000,000.

Figure 6 .
Figure 6.An example of the outspreaded bucket of a low Reynolds number airfoil (SG6043) operating at Re = 350,000, in contrast with the narrow bucket of a NACA 66 2 -415 operating at Re = 3,000,000.

Figure 6 .
Figure 6.An example of the outspreaded bucket of a low Reynolds number airfoil (SG6043) operating at Re = 350,000, in contrast with the narrow bucket of a NACA 662-415 operating at Re = 3,000,000.

Figure 7 .
Figure 7.A way to quantify the stall quality.

Figure 7
Figure 7 represents an ideal Cl-α diagram.Left part contains αs and Clmax points.Further to the right than the critical as point, three hypothetical stall curves displayed.Obviously, the preference is green > orange > red.We consider a vertical line at distance Δx from the αs

Figure 7 .
Figure 7.A way to quantify the stall quality.

Figure 7
Figure 7 represents an ideal C l -α diagram.Left part contains α s and C lmax points.Further to the right than the critical a s point, three hypothetical stall curves are displayed.Obviously, the preference is green > orange > red.We consider a vertical line at distance ∆x from the α s point that intersects with the three stall curves.Afterwards, in the magnified right part of the diagram, we record ∆ y1 , ∆ y2 , ∆ y3 points.The curve which records the smallest ∆y is the one to receive the best score to contribute to its selection process.Unfortunately, things are not quite so easy, because with the aforementioned method we are studying ideal curves that are unlikely to exist in real simulation results because:

5. 3 .
Construction of the Airfoil «Map-Diagram» The necessary visual aid to help the designer in the airfoil selection process is the maximum lift coefficient versus ideal lift coefficient diagram or C lmax − C li diagram.This diagram is a form of an airfoil map.C li values lies at x axis and C lmax values at y axis.Within the diagram, points exist, each one representing a certain airfoil.The designer can focus on a certain region of the diagram to further analyze the airfoils by comparison.

Figure 8 .
Figure 8.The constructed C lmax − C li diagram (the diagram along with raw data can be found in .xlsxformat in TableS2).
li Preferred are those airfoils whose Cli is nearest to the calculated Cli.

Table 5 .
The twenty airfoils analyzed.Table 4 contains the Reynolds numbers in increments of 50,000.4.1.Literature Review in Low Reynolds Number Airfoils for Small UAV's

Table 5 .
The twenty airfoils analyzed.Table4contains the Reynolds numbers in increments of 50,000.
4.1.Literature Review in Low Reynolds Number Airfoils for Small UAV's

Table 5 .
The twenty airfoils analyzed.Table4contains the Reynolds numbers in increments of 50,000.
4.1.Literature Review in Low Reynolds Number Airfoils for Small UAV's

Table 6 .
Measurements of interest at Re = 100,000 and no flap usage.All variables' units are dimensionless, except C m units which are 1/rad and α i and α s which are degrees.

Table 7 .
Measurements of interest at Re = 200,000 and no flap usage.
All variables' units are dimensionless, except C m units which are 1/rad and α i and α s which are degrees.

Table 8 .
Measurements of interest at Re = 350,000 and no flap usage.All variables' units are dimensionless, except C m units which are 1/rad and α i and α s which are degrees.

Table 9 .
Measurements of interest at Re = 100,000 and 30 • flap deflection at 80% of chord.All variables' units are dimensionless, except α s which are degrees.

Table 10 .
Weights of each Re, towards final airfoil stall quality calculation.