Abstract

The commercially available unmanned aerial vehicles are not good enough for search and rescue flight at high altitudes. This is because as the altitude increases, the density of air decreases which affects the thrust generation of the UAV. The objective of this research work is to design thrust optimized blade for an altitude range of 3,000–5,000 m with a density of air 0.7364 kg/m3, respectively, and perform thrust analysis. The property of aluminum alloy 1,060 being lightweight is chosen for designing and testing of blade. The blade element theory-based design and analysis code was developed, and user-friendly aerodynamic inputs were used to obtain the desired outputs. The geometry designed for an altitude range of 3,000-5,000 m faced the total stress of 6.0 MPa which was at 70% of the blade span. This stress is within the limit of yield strength of the aluminum alloy, 28 MPa. The modal analysis shows the first natural frequency occurs at around 12,000 RPM which is safe for operating the blade at 0-5,000 RPM. Experimental analysis of the blade gave a thrust of 0.92 N at 2,697 RPM at 1,400 m. The analytical solution for thrust with the same conditions was 1.7 N with 85.6% efficiency. The validation of experimental results has been done by the CFD analysis. The CFD analysis was performed in ANSYS CFX which gave a thrust value of 2.27 N for the same boundary conditions. Thus, the blade designed for high altitude SAR UAV is structurally safe to operate in 0-5,000 RPM range, and its use in search missions could save many lives in the Himalayas.

1. Introduction

Unmanned aerial vehicles (UAVs) have experienced rapid development in recent years. Impressive achievements have been made by the use of the UAV in civilian and military fields, for example, aerial photography, surveillance, remote sensing, agricultural monitoring, and fire control. There are many institutes which have been actively involved in exploring the more potential applications of UAV. These UAVs range from large () [1, 2] to mini (takeoff , ) [3, 4] and even micro-UAVs (takeoff weight 5 kg, altitude 250 m) [5]. High altitude long endurance (HALE) UAV (17-25 km) [6] has been growing interest in recent years. The altitude and payload mass have a huge influence on the size and design of UAV [7]. Helios (wingspan 73 m) developed by NASA and UAS ascended to an altitude of over 29 km [6]. Unfortunately, Helios was destroyed during the flight test on June 26, 2003, because of turbulence and structural failure [8]. Similarly, Zephyr (wingspan 25 m, weight 75 kg) designed and built by British Company QinetiQ holds the official endurance world record of 14 days, 22 min, 8 sec without refueling [9].

The conventional airfoils such as the NACA airfoils were designed to operate at high Reynolds numbers [10]. Thus, for optimum aerodynamics characteristics, the low Reynolds number airfoil should be used for UAVs [11]. A few low Reynolds number (Re < 300,00) airfoil have been designed for the maximum lift to drag () ratio. [12, 13], such as A18 (), S6062 (), SD7032 (), and BW-3 () [14]. For low speed and high lift to drag ratio (, at Re 50,000 and 3.25° AOA), S1223 airfoil is most suitable [15, 16].

Most UAVs have fixed pitch propellers [17, 18]. There are various theories which predict the efficiency of fixed-pitch propellers. Some are blade element theory, vortex theory [19], momentum theory [20, 21], flat plate theory [22], and semiempirical methods [23]. Blade element theory is simple and efficient. It accounts for the Reynolds number effects in low advance ratios [2426]. Thus, the blade element theory is widely used to model propeller blade aerodynamics.

Weerasinghe and Monasor [27] have carried out computational fluid dynamics (CFD) analysis and experimental wind tunnel tests on Syma X5SC quadrotor; even if simulation predictions were away from the experimental results at some flight conditions, CFD simulations were found good in describing complex flows of rotatory movements in overall [28]. Selig and Brandt found out that propeller performances decrease dramatically at lower rotational speeds with the effect of Reynolds number [29]. Past literature has shown that thrust and power coefficient determination methods mainly focused on wind tunnel and PIV measurements, CFD analysis, and analytical models [28].

The rescue operation of missing tourists in the Himalayas of Nepal is challenging. The authorities have to travel either by helicopter or on foot to the remote areas. The air being thin at higher altitudes makes it difficult to perform the rescue operation via UAVs. Liu et al. [30] found that with the increase in altitude, the UAV climb rate gradually decreases. It can be seen from Figure 1 that at 5,000 m altitude, thrust generation decreases by about 40%. This research work will help to improve propulsive efficiency for future SAR UAVs. The use of SAR will not only be limited in the field of search and rescue but also in high-altitude surveillance and mapping. The search and rescue UAV used by Mario and the team operated its flight from 300 masl to 2,600 masl. [31]. There are various UAVs which have been developed and tested with the purpose of surveillance and remote sensing for very high altitude [3235]. Regarding the development of propeller for high altitude SAR UAVs, Korea Aerospace Research Institute have built the propeller with a diameter of 0.61 m to fly at 10 km ASL altitude with propeller efficiency of 68% (at 3,120 RPM) [36]. The other 3-bladed propeller with 72% efficiency climbed at 22 km altitude with 0.5588 m diameter propeller generating 7 N thrust at 50 m/s free stream velocity [37].

In Nepal, UAVs have played a major role during an earthquake in 2015. But lack of rescue operations in the Himalayan has taken many innocent lives [38]. On April 30, 2018, the nonprofitable organization National Innovation Centre tested its medical drone in remote village Rangi of Maygdi district, Nepal (altitude 2,300 masl) [39]. Although the flight was successfully tested, but the quadcopter performance was not satisfactory. There are various factors that affect the working of the propeller. The major ones are air density, size of the propeller blade, and pitch of the blade. We can see from Figure 1 that propeller thrust is directly proportional to the air density. In the atmosphere, air density decreases as the altitude increases. The relation between attitude and density is fairly exponential. This loss of thrust is due to the change in density of air which can be optimized by a differently designed propeller. Thus, the main objective here is to design and test the high-altitude fixed-pitch propeller of search and rescue unmanned aerial vehicles for maximum thrust.

The assumptions are as follows: (i)Fixed-pitch blade(ii)Chord (b) is given by equation (4)(iii)Blade span radius is 0.15 m(iv)Two number of blades(v)Variations of the angle of attack do not affect the efficiency(vi)The chosen airfoil S1223 is an efficient design for higher altitude

2. Methodology

2.1. Simple Blade Element Theory

Simple blade element theory is the most used propeller theory. In this theory, the blade is divided into different elements. Lift along with torque is calculated by each of these elements which is then summed to find total thrust and torque [4042]. The total thrust for the propeller having number of blades is given by equation (1). Similarly, the torque is given by equation (2), where is the number of blades, is the air density, is the velocity, is the chord, is the radius of blade, and is coefficient of lift. The efficiency is given by equation (3). Figure 2 [43] shows the velocity vector diagram with reactions in the blade.

It is assumed that the propeller tip which has the highest relative velocity does not exceed the speed of sound; thus, all the location does not exceed the speed of sound. It can be clearly seen from equation (3) that increasing the free stream velocity aids in increased efficiency.

2.2. Geometric Modelling

The design model is a linear and explicit type that the inputs for calculating thrust and other forces are solved by known input parameters such as density and earth’s altitude. It is time-dependent and the approach is deterministic. Thrust calculation was based on equation (1); other assumed parameters are the density of air at 4,000 m is 0.8194 kg/m3, the number of blades is 2, and for AOA 4° is 1.35 and 0.05, is 33,000, the radius of the blade is 0.15 m, and the chord is given by equation (4) [44]. The propeller will be used with 0-4000 RPM limit, whereas for the factor of safety, certain analysis has been performed till 5,000 RPM.

In Table 1, aluminum alloy properties give the aluminum alloy 1,060 physical properties. The compositions include , , and . This value is used for calculations and analysis here. Figure 3 shows the maximum thrust received by the blade is 3.35 N at 4,000 RPM. The figure also shows thrust generation for 500-3,500 RPM. The increasing trend suggests that the thrust generation will increase rapidly if the rotational velocity is increased. But vibration analysis should be carried out for increasing the RPM.

The three different blades that have been designed for the same thrust output (3.35 N) are shown in Figure 4. If the blade with the same thrust generation of 3.35 N is used at sea level altitude, then the required radius reduces to 134 mm. Similarly, the blade is used at an altitude of 8,000 m for the same thrust output that the required radius is increased to 167 mm. The blade used for altitude 4,000 m was further analyzed at different free stream velocity, i.e., 10 m/s, 50 m/s, 100 m/s (shown in Figure 5). With a higher velocity of air, a higher twist angle was required as shown in Figure 6. For the analysis purpose, the blade designed for 4,000 m altitude and 10 m/s velocity is chosen here. The designed blade has 7.1 mm thick root and 1.5 mm thin tip.

3. Structural Design and Analysis

3.1. Shear Force and Bending Moment

The shear force and bending moment diagram of the blade can be seen in Figures 7 and 8. The bending moment is based on resultant thrust loading. Since the thrust force acts upwards, the bending moment acts downward or anticlockwise. Furthermore, as the radius of the blade increases, the shear force and bending moment increases as well. All this analysis is limited to 5,000 RPM, because this is the maximum speed the blade may face. The geometry of the blade is designed for 4,000 RPM only. So, for the factor of safety, structural analysis was performed till 5,000 RPM. Increasing the speed higher than 5,000 RPM will cause the blade root to face more stress and may fail.

3.2. Stress Concentration

The load produced by lift force causes the blade to withstand bending stress and the rotational speed causes centrifugal stress. The stress concentration is calculated using equation (5). where is the stress produced due to bending moment, is the stress produced due to centrifugal force, is the calculated bending moment, is the vertical distance, is the moment of inertia, is mass of blade, and is the area of the airfoil at a particular section.

Figure 9 shows total stress distribution and the combination of centrifugal and bending stress along the span of the blade. The maximum stress (6.5 MPa) faced by the blade is at 70% of the blade span towards the root and 2.7 MPa at the root. This stress compared to aluminum yield strength 28 MPa is in the limit; thus, the blade is safe to operate with a rotational velocity of 0-5,000 RPM. In Figure 10, we can see the results of total stress distribution which is obtained from ANSYS 18. It can be seen that with the same boundary conditions and at 4,000 RPM, the maximum stress is 9.35 MPa. Thus, the blade will operate with out failure within the rotational speed of 0-4,000 RPM.

3.3. Modal Analysis

Figure 11 shows the Campbell diagram for the blade with propeller rotation speed in -axis and system frequency in -axis. It can be seen that the first natural frequency of the blade occurs at 200 Hz (), operating the blade at 12,000 RPM will excite the critical frequency at 4,500 RPM which may lead to failure. But the blade operating range is 0-4,000 RPM; thus, it is safe from resonance and other modes of vibrations.

4. CFD Analysis

4.1. Computational Model

The computational domain is divided into 2 parts (Figure 12). The stationary domain takes 10 times the diameter of the propeller in upstream and downstream. The span-wise diameter is 10 times. Inlet, outlet, and propeller outer region are stationary. The rotating domain with cylindrical shape has diameter 400 mm and height 400 mm. Figure 13 shows individual mesh for stationary and rotating domains along with the propeller blade. Figure 14 shows the mesh independence test. Three different mesh structure was used to guarantee the results’ independency at the propeller rotation speed of 2,679 rpm. The number of course mesh is 0.39 million; the used mesh (fine mesh) is made by 0.42 million cells. The maximum discrepancy between the fine mesh and the course mesh is only 0.26% for the generated thrust. Thus, the results are not affected by the used mesh. Hence, the fine mesh is used in the CFD simulation of the designed propeller.

The mesh quality metrics are to assess the suitability of a computational domain for simulations. There are various parameters that define the mesh quality. Such as most of the element’s aspect ratio is 7-9. Furthermore, Figure 15 shows the element quality and Figure 16 shows the orthogonal quality. It can be noticed from Figures 15 and 16 that not all the 0.4 million elements are shown. This is because the bad quality elements are very small in number, and it is hard or near impossible to spot in a full-scale figure.

4.2. Computational Method

The flow field of the propeller is stationary under the rotating coordinate system. The ANSYS CFX turbo mode is used to solve the problem. Reynolds-Averaged Navier-Stokes equations are regarded as the governing equations. The finite volume method with the pressure-based solver is used to discretize the governing equations. Second-order upwind (high resolution) is used for the advection scheme. Pérez et al. carried out the CFD of the quadrotor’s propeller and found out that k-w and Spalart Allmaras turbulence models both overestimated thrust and torque compared with experimental results taken from flight tests [45]. In order to simulate the turbulence, the standard k-e model [46] which requires flow to be fully turbulent is used [47]. The convergence is also guaranteed by monitoring the residual value which drops below .

4.3. Boundary Conditions

All the boundary conditions are applied based on the experiment carried out at 1,400 m. The flowing fluid is considered as air at 20°C and the flow is steady. The boundary conditions are shown in Table 2 which are the assumptions based on the past authors’ work [15, 4850].

5. Manufacturing and Testing

The designed blade was 3D printed first, and sand casting was performed. The final results can be seen in Figure 17. The aluminum alloy 1,060 was used for casting. A static balance check was performed before the final thrust test. The thrust testbed can be seen in Figure 18. The test rig was chosen as length 1 m and height 1 m. The test condition includes the free stream velocity of air as negligible as the testbed was fixed. The thrust test results are shown in Table 3. The blade developed the thrust of 0.932 N at 2,679 RPM. Due to the testbed limitations, higher RPM could not be tested. Although the blade was designed for higher altitude, the test was carried out at 1,400 m. The propeller was made up of aluminum, and to test it, we needed a testbed enclosed inside a box with one side open to atmosphere. Carrying the testbed to a higher altitude was empirically challenging from both financial and technical perspectives. Thus, the experiment was carried out at 1,400 m altitude.

6. Results and Discussions

The results obtained from analytical values and experimental and numerical values have been compared. The comparison of all values can be seen in Table 4. The experimental value gave 0.92 N thrust at 2,679 RPM. The corresponding value from analytical calculations were 1.7 N. The CFD analysis gave the corresponding value of 2.27 N at 2,679 RPM. The experimental results and analytical value gave 45% test results. This could be due to various factors involved in the experiment. Such as for safety purposes, an experiment was conducted inside a closed room so the free stream air velocity was ignored; aerodynamics profile of the fabricated propeller may not be accurately the same to the designed propeller.

The pressure contours as seen in Figures 19 and 20 show the pressure distribution profile. The blade faced the highest pressure of 651 Pa at the pressure side and the lowest pressure of -700 Pa at the suction side.

7. Conclusion

The blade was designed for an altitude range of 3,000-5,000 m which is suitable for search and rescue operations in Nepal. Although the blade could not be tested at high altitude, but the test results performed at the IOE laboratory and CFD analysis show that it can be used for SAR missions in Nepal. Here are the conclusions that can be drawn from this research work. (i)The CFD thrust analysis gave thrust generation of 2.27 N for an altitude of 1,400 m with a rotational velocity of 2,679 RPM. Experimental thrust 0.92 N was less accurately predicted by CFD. Past literature has shown the error between CFD, and the experiment varied from 14 to 24% [49, 51]. However, here, numerical analysis results differ largely from experimental results. From Figure 14, we can see that with more element size mesh, we get more accurate results. Thus, further investigations are suggested in order to reduce the differences in the results(ii)The total stress faced by the blade is 6 MPa which is around 70% span of the blade which is safe for the blade compared to aluminum yield strength of 28 MPa(iii)The thrust rating test showed the thrust generation of 0.92 N at 2,679 RPM. The analytical value based on test conditions for 2,679 RPM would be 1.7 N. This gives 45% of the test results

Abbreviations

AOA:Angle of attack
:Chord length
FEA:Finite element analysis
IOE:Institute of Engineering
masl:Meter above sea level
PIV:Particle image velocimetry
:Revolution per sec
Re:Reynolds number
SAR:Search and rescue
UAS:Aeronautical environment corporation
:Stress produced due to bending moment
:Stress produced due to centrifugal force
:Moment of inertia
:Efficiency
:Blade/helix angle
:Density of air
:Radius of blade
:Angle of attack
:Flow angle
:Angle between reaction force (lift and drag components) and lift
:Relative velocity
:Free stream velocity
:Rotational velocity.

Data Availability

The data used to support the observations of this study are available from the corresponding author upon request.

Disclosure

The funding institute had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Conflicts of Interest

The authors declare that there is no conflict of interest with regard to the publication of this paper.

Acknowledgments

The authors would like to thank, Nepal Academy of Science and Technology (NAST) for providing financial support to complete this work.