alexa Optimal Design and CFD Analysis of Wing of a Small-Scale UAV to Obtain Maximum Efficiency

ISSN: 2168-9792

Journal of Aeronautics & Aerospace Engineering

Reach Us +44-1904-929220

Optimal Design and CFD Analysis of Wing of a Small-Scale UAV to Obtain Maximum Efficiency

Khuntia SK* and Ahuja AS
Department of Aeronautics, Naval Institute of Aeronautical Technology Naval Base, Kochi, India
*Corresponding Author: Khuntia SK, Head, Department of Aeronautics, Naval Institute of Aeronautical Technology Naval Base, Kochi-682004, India, Tel: +919969613377, Email: [email protected]

Received Date: Jan 20, 2018 / Accepted Date: Jan 31, 2018 / Published Date: Feb 04, 2018

Abstract

In the recent past small-scale Unmanned Aerial Vehicles (or Unmanned Aerial System as it is called presently) have evolved as an important tool in non-conventional fields like agriculture, e-Commerce, policing, medical logistics in addition to military applications. This paper presents the complete methodology (very limited article exists) applied to optimally design the wing of a small scale Unmanned Aerial Vehicle with help of widely used CFD software, ANSYS to maximize its efficiency. In this study, the application of computational methods in the iterative design process is successfully explored. The various design parameters and features of wing also have been explored with help of CFD analysis to derive the advantages. The strength and stiffness analysis of the UAV wing has also been carried out using ANSYS. The modelling of wing and associated parts were carried out on CATIA. The final design of the wing has two spars, 20 ribs a side with mid wing root strengthening, near elliptical planform and has an all composite structure. The wing is lighter in weight as compared to a similar wing made from Aluminum view weight optimization, and strong enough to meet all in-flight load conditions with safety margins. This paper demonstrates the design process/methodology to optimally design an efficient small-scale UAV.

Keywords: UAV wing; CATIA; CFD analysis; CFX; ANSYS; Contour of turbulence; Aspect ratio; Static pressure distribution; Flow analysis; Directional deformation

Introduction

Unmanned Aerial Vehicles (UAV) innovations started in the early 1900s which originally focused on providing practice targets for training military personnel. UAV development continued during World War I which led to the rise of early powered unmanned aerial combat vehicles (UCAV) in the year 1915 by Nikola Tesla. It was during World War II during which Nazi Germany used UAV for fly attack missions. Over last decade the impetuous on UAV or RPA (Remote Pilotless Aircraft) has increased manifold.

The age of the small UAV started in 1987, when Dr Paul MacCready's company AeroVironment developed the Pointer, the first hand-launched, backpack-carried UAV. Pointer combined the airframe technology of a high-performance model sailplane with an electric motor and propeller, a consumer video camera and a radio datalink. In 1999, the US Army bought four AeroVironment Pointer, small UAVs for testing in the service's "Military Operations in Urban Terrain" and was enthusiastic about the usefulness of the Pointer. In recent years small scale Drones have found their utility in a variety of important fields like Agriculture, e-Commerce, Logistics, urban warfare, Policing, Traffic management, climate monitoring, Land survey, Media coverage etc. The list is endless. Be it a large-scale aircraft or a small-scale UAV, the design process is iterative and follows a certain methodology. The process is iterative because the design parameters have to be evaluated over and over again so as to reach an optimal design. It is one of the most challenging fields in engineering and requires highly creative thinking. It requires a compromise between the various design parameters and eventually optimization [1-5].

Phases of aircraft design

The most important part of any aircraft design is the wing that provides lift to make it airborne. Hence Wing design has been analysed critically in this paper. However, this section gives a glimpse of complete design process. Aircraft design can be broken into three major phases namely Conceptual Design, Preliminary and Detailed design.

Conceptual design: This phase consists of selecting a basic design depending upon the requirements. It includes steps like initial requirements, approximate weight, draft design and the type of propulsion.

Preliminary design: This phase involves freezing of the configuration although it is open to minor revisions. Major components are designed in this stage followed by “lofting”. Lofting is modelling of outer skin of the aircraft with sufficient accuracy to ensure proper fit between different parts [6-8].

Detail design: This is the final phase in which the actual pieces to be built are designed which included detailed parts like individual ribs, spars etc. Another important part of detail design is called production design which determines how the airplane will be fabricated, starting with the smallest and simplest sub-assemblies and building up to the final assembly process.

The conceptual design process starts with the selection of suitable configuration, which determines the results of entire conceptual design process. Therefore, special care must be taken to choose a suitable but yet reasonable configuration in this stage. The entire conceptual design process is dependent on configuration. Hence, while choosing the configuration we must take into account the following subjects:

(a) Design requirements.

(b) Historical trends.

(c) Operational ease

Selecting the correct configuration is very difficult in initial phases because nothing is fixed at this stage. If the chosen design does not provide required performance criterion then the entire conceptual design has to be discarded. Therefore, care must be taken in selecting the right configuration. Once the theoretical design is complete on paper the concepts of Computer Aided modeling of chosen aircraft wing were understood with help of CATIA V5R21 and the flow analysis was carried out using CFD Solver Fluent in ANSYS. Flow Analysis was carried out on the wing section iteratively to offer an efficient design for the same. During the design process, various high lift devices and aspects to maximize the efficiency of the wing were also analyzed viz wing taper, dihedrals etc. [9-11].

CFD Analysis

Governing equations

The set of equations solved by ANSYS CFX are the unsteady Navier-Stokes equations in their conservation form. For all the following equations, static (thermodynamic) quantities are given unless otherwise stated. Transport Equations, the instantaneous equation of mass, momentum, and energy conservation majorly contribute to Computational Fluid Dynamics. For turbulent flows, the instantaneous equations are averaged leading to additional terms. The instantaneous equations of mass, momentum and energy conservation can be written as follows in a stationary frame:

Continuity equation equation(1)

Momentum equation

equation(2)

Where equation(3)

I: Local energy change with time

II: Convective term

III: Pressure work

IV: Heat flux (diffusion)

V: Irreversible transfer of mechanical energy into heat

Generalized Navier stroke equation

To simplify the Navier-Stokes equations, we can rewrite them as the general form.

equation(4)

When equation, we can respectively get continuity equation, momentum equation and energy equation. Another general form of NSE is

equation (5)

The Navier-Stokes equations are analytical equations. Human can understand and solve them, but if we want to solve them by computer, we have to transfer them into discretized form. This process is discretization. The typical discretization methods are finite difference, finite element and finite volume methods.

Finite volume approach: Integrating the general form of Navier- Stokes equation over a control volume and applying Gauss Theory:

equation(6)

We obtain integral form of Navier-Stoke equation:

equation(7)

Once purpose and CFD codes are chosen, “CFD process” is to assemble the system of algebraic equations and solve the system to get approximate solutions for all flow properties over wing. Figure 1 explains a step by step CFD process.

aeronautics-aerospace-engineering-flow-chart

Figure 1: Flow chart of CFD process.

The approximate solution in an CFD process is obtained by Integration of the governing equations of fluid flow over all the (finite) control volumes of the solution domain. This is equivalent to applying a basic conservation law (e.g. mass or momentum) to each control volume and time marching. Typically, we have two types of interpolations, one is upwind interpolation, and the other one is central interpolation (Figures 2-10).

aeronautics-aerospace-engineering-air-foil

Figure 2: Air-foil coordinates.

aeronautics-aerospace-engineering-ribs

Figure 3: Internal design of wing with ribs and spar.

aeronautics-aerospace-engineering-span

Figure 4: Full wing span.

aeronautics-aerospace-engineering-uav

Figure 5: Modeling of a small scale UAV.

aeronautics-aerospace-engineering-nodes

Figure 6: After applying nodes.

aeronautics-aerospace-engineering-geometry

Figure 7: Full modelled geometry of wing.

aeronautics-aerospace-engineering-pressure

Figure 8: (a) Static pressure and (b) Force report of SA 7035.

aeronautics-aerospace-engineering-senig

Figure 9: (a) Pressure coefficient (b) Force report of Senig A 7035 at 4° angle of attack.

aeronautics-aerospace-engineering-coefficient

Figure 10: (a) Pressure coefficient (b) Force report of Senig A 7035 at 8° angle of attack.

Weight estimation

Design take off gross weight” is the total weight of the aircraft as it begins the mission for which it is designed. This is not necessarily the same as the “maximum take-off weight”. Design take off gross weight can be broken into crew weight, payload weight, fuel weight, and empty weight which include weight of structure, engine, landing gear, avionics etc.

equation(8)

The crew and payload weights are both known since they are given in the design requirements. The unknowns are the fuel weight and empty weight but they are dependent on the total aircraft weight. Thus, an iterative process must be used for aircraft sizing. Thus following equation is obtained.

equation (9)

This can be solved for W0 as follows:

equation (10)

equation (11)

Now W0 can be determined if (We/W0) and (Wf/W0) can be estimated called empty weight estimation and fuel fraction estimation respectively.

Weight of fuel and crew is zero in our case as we are considering a battery operated remote controlled aircraft. Also the weight of the motor is considered separately. So,

equation(12)

Equation 1 becomes,

equation(13)

In the conceptual design phase of the airplane whether it is manned or unmanned, first step is to estimate the take-off gross weight (Figures 11-15). It is an initial rough estimate that gives the “ball-park” figure. (Eq. 6).

aeronautics-aerospace-engineering-report

Figure 11: (a) Pressure coefficient (b) Force report of Senig A 7035 at 12° angle of attack.

aeronautics-aerospace-engineering-pressure

Figure 12: (a) Pressure coefficient (b) Force report of Senig A 7035 at 16° angle of attack.

aeronautics-aerospace-engineering-density

Figure 13: Flow density over wing.

aeronautics-aerospace-engineering-pressure

Figure 14: Total pressure coefficient.

aeronautics-aerospace-engineering-turbulent

Figure 15: Turbulent velocity over wing.

Since, payload weight is known from requirements, fuel and empty weights are both dependent on total take-off weight.

Therefore the equation becomes:

equation(14)

Considering the following weights as rough estimates based on the various small scale UAVs available:

Rough weights of components: From the above Table 1,

equation

S. No Description Weights (Grams)
1 Motor 350
2 Payload 150
3 Propeller 30
4 Servos 120
5 First person view camera 50
6 Telemetry 30
7 Structure 2000
8 Battery 370

Table 1: Rough estimate of UAV parts.

The WGuess can be any weight for the calculation purpose.

Substituting all the values in equation 6.6, we get:

(a) First estimate: Consider WGuess = 3010 grams, so W0 becomes

equation

(b) Second estimate: Consider WGuess = 3030 grams, so W0 becomes

equation

(c) Third estimate: Consider WGuess = 3050 grams, so W0 becomes

equation

(d) Fourth estimate: Consider WGuess = 3080 grams, so W0 becomes

equation

(e) Fifth estimate: Consider WGuess = 3095 grams, so W0 becomes

equation

(f) Sixth estimate: Consider WGuess = 3100 grams, so W0 becomes

equation

So it is estimated that the gross take-off weight of the aircraft would be about 3100 grams.

Air-foil selection

Initially three airfoils viz., Clark Y, Senig-Donovan 7032 and Senig-Ashok Gopalaratnam (SA) 7035 were considered for flow analysis using ANSYS software to find out lift and drag produced by them at 20 m/s and SA 7035 was the best choice due to its high lift and low drag as compared to other options. This airfoil was finalized for the UAV. A rectangular mono wing with high wing mounts of a span 04 metres and chord length of 0.3 metres was selected. The other characteristics of SD 7035 are as follows:-

(a) Aspect ratio: equation

Cl max of airfoil is 1.4

(b) Stall speed: equation

(c) Take off speed: equation

(d) Load factor: equation

equation

(e) Velocity during flare: 1.15*Vstall = 12.054 ft/sec

(f) Touch down velocity: 1.23*Vstall = 11.27 ft/sec.

Note: Dry air has a density of 0.076 lb/ft3, S = Wing Area = 13.13 ft2

CFD analysis

CFD is the art of replacing the differential equation governing the Fluid Flow, with a set of algebraic equations (the process is called discretization), which in turn can be solved with the aid of a digital computer to get an approximate solution. After importing UAV components into Ansys we define the analysis types by applying loads and initial conditions for the finite element solution. During mesh generation loading boundary conditions of inlet outlet wall and symmetry conditions are then applied to this elements and nodes of the UAV wing structure. For simulation of fluid flow analysis the UAV model was transferred to CFX/ FLUENT preprocessor. In this process details for fluid and solid domains were assigned. Material used for fluid type domain in our case is air at 250 C. Inlet and outlet details in fluid and Boundary Conditions were applied prior flow analysis we obtain the pressure counter and velocity streamline of fluid flow (Figures 16-20).

aeronautics-aerospace-engineering-dihedral

Figure 16: Force report of SA 7035 with a dihedral angle 3.50 and 0° angle of attack.

aeronautics-aerospace-engineering-taper

Figure 17: Force report of SA 7035 with a taper ratio of 0.5 and 0° angle of attack.

aeronautics-aerospace-engineering-intensity

Figure 18: Stress intensity over wing.

aeronautics-aerospace-engineering-directional

Figure 19: Directional deformation of wing under stress.

aeronautics-aerospace-engineering-strain

Figure 20: Strain energy.

Grid generation

The sub- division of the flow domain into a number of smaller, non-overlapping sub-domains is the concept behind Grid generation. Accuracy of a solution, calculation time and cost in terms of necessary computer hardware are dependent on the fineness of the grid. Over 50% of time spent in industry on a CFD project is devoted to the definition of domain geometry and grid generation.

Flow analysis of wing

One of the most important criteria for any aircraft design is airfoil selection for its major lift producing surface called main plane or wing. The three airfoils viz., Clark Y, SD7032 and SA7035 were studied in ANSYS for the lift and drag at 20 m/s. The best results for Lift and Drag were that of SA7035. The Static pressure and the force reports are as follows

SA7035 was found to have the highest lift and lowest drag among the three, so this airfoil was finalized for the UAV. The wing was again studied for various flow parameters at different angles of attack and structural loads. The critical values were studied for setting up a safe flying envelope.

As per analysis, the stall angle is between 8 degrees and 12 degrees when the drag increases drastically. At 12 degrees there is a clear indication of flow separation near quarter chord region which can be considered as the main cause of Stall. It was also found that at 16 degrees, drag is more than lift which clearly indicates stalling of the wing. After simulating the stall angle with help of ANSYS Flow density, Contour of Total Pressure and effect of turbulence was also observed for any unusual behavior. With the chosen airfoil, chord, span and wind loading it was found that the aerodynamic design of wing is optimal and efficient without any unusual flow separation at higher angle of attacks.

Wing designed with SA7035 airfoil was also considered at a dihedral of 3.5 degrees for analysis and the lift generated (net force of 1.08) was found to be lesser than the wings with nil dihedral angle (1.5). Hence a better design for low speed small scale UAV would be a wing without any dihedral. The force report with dihedral is as follows:

The wing was subsequently considered for a taper ratio (λ) of 0.5 for analysis to find out if the design was efficient but the lift produced was very less when compared to a taper ratio of 1 i.e., a rectangular wing. The modeling and force report is as follows

Structural analysis of wing

In order to determine the effects of inertial loads during manoeuvres and natural frequencies, a Dynamic analysis was also performed. Dynamic analysis comprised of deformation calculations due to inertial loads at higher load factors and estimation of natural frequencies and mode shapes.

Discussion

With an aim to obtain an efficient wing for a long range and long endurance small scale UAV an iterative design an selection process was followed. The salient features of the wing are as follows:

(i) Static pressure distribution

It was observed that the static pressure rapidly increases as the angle of attack increased from 8 degrees to 12 degrees indicating that the stall angle lies in that region. Also at 16 degrees, the drag is more than lift which clearly indicates stalling of the wing. Stall should develop progressively from root to tip which is ideal Stall Progression (Figures 8-12 and Table 2).

Angle of Attack (in Deg) Lift (in N) Drag (in N)
15.44 3.37
22.56 13.56
12° 30.08 29.08
16° 33.08 36.36

Table 2: The salient features of the wing.

(ii) Aspect ratio

To obtain long range, therefore, the designer is driven to design an aircraft which will cruise at high altitude and have a long wing in order to reduce the induced drag at high altitude. This results in a very slender wing of aspect ratio perhaps in the range 20–25 which then gives a structural design challenge to achieve it without incurring excess weight. The Aspect Ratio in our case (small scale UAV) was optimally obtained (13.34) keeping in mind the structural strength through stress analysis in an iterative process.

(iii) Contour of turbulence

It is required to study the contour of turbulence on wing to obtain the advantage of lower response to air turbulence and the versatility of hover around a selected way-point and low-speed flight (for better endurance and low noise). An indicative measure of the response to air turbulence of an aircraft, and to some degree its relative aerodynamic efficiency, may be given by the ratio ˄ of its surface area to its mass. The larger the surface area, the more it may be disturbed by aerodynamic forces. The greater its mass, so greater will be its inertia (resistance) to the imposed forces. Using the scaling laws, it may be seen that the area/ mass ratio ˄ will vary as n/ρD, the linear dimension ratio n divided by the packaging density ρD (Figures 13 and 15).

(iv) Dihedral lateral stability

Dihedral provides lateral stability to aircraft during its rolls. However, it comes with a cost that is decreased roll rate and increased drag which was quite evident from the Force Report of our analysis. Hence for mall scale UAV there is nil requirement of anhedral or dihedral (Figure 16).

(v) Taper ratio

The selection of optimum wing span thus requires an analysis of much more than just cruise drag and structural weight. Here, again, a diverse set of considerations are important. The major design goal is to keep the taper ratio as small as possible (to keep the wing weight down) without excessive Cl variation or unacceptable stalling characteristics. Since the lift distribution is nearly elliptical, the chord distribution should be nearly elliptical for uniform Cl's. Reduced lift or t/c outboard would permit lower taper ratios (Figure 17).

(vi) Structural strength

Reported studies on aircraft failures indicate that the main cause of wing failure is due to fatigue cracks which nucleate from the wing root region. The stress intensity shows a red zone at the root and half chord region. Hence it is apt to strengthen the mid spar of the wing (Figures 18-20).

Conclusion

An old military adage which is quite apt for civilian usage too links the use of UAVs to roles which are dull, dirty or dangerous (DDD). There may be much truth in that but if you add surveillance, covert, diplomatic, research, environmental like critical roles then the dynamism of present day UAV is complete. In addition, the economics of operation are often to the advantage of the UAV [11]. This study of ours proves that though the design of any aircraft is a compromise between several aims the various important parameters and features can be tweaked to achieve desired result in terms of maximum efficiency, better L/D ratio and stall characters.

References

Citation: Khuntia SK, Ahuja AS (2018) Optimal Design and CFD Analysis of Wing of a Small-Scale UAV to Obtain Maximum Efficiency. J Aeronaut Aerospace Eng 7: 207. DOI: 10.4172/2168-9792.1000207

Copyright: ©2018 Khuntia SK, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Select your language of interest to view the total content in your interested language

Post Your Comment Citation
Share This Article
Relevant Topics
Recommended Conferences
Viewmore
Article Usage
  • Total views: 1338
  • [From(publication date): 0-2018 - Oct 18, 2018]
  • Breakdown by view type
  • HTML page views: 1175
  • PDF downloads: 163

Post your comment

captcha   Reload  Can't read the image? click here to refresh
Leave Your Message 24x7