 Open Access
 Total Downloads : 157
 Authors : Sumit Sharma
 Paper ID : IJERTV5IS110361
 Volume & Issue : Volume 05, Issue 11 (November 2016)
 DOI : http://dx.doi.org/10.17577/IJERTV5IS110361
 Published (First Online): 28112016
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
An Aerodynamic Comparative Analysis of Airfoils for LowSpeed Aircrafts
Sumit Sharma
B.E. (Mechanical) MPCCET, Bhilai, C.G., India
M.E. (Thermal) SSCET, SSGI, Bhilai, C.G., India
Abstract – This paper presents investigation on airfoil S1223, S819, S8037 and S1223 RTL at low speeds to find out the most suitable airfoil design to be used in the low speed aircrafts. The method used for analysis of airfoils is Computational Fluid Dynamics (CFD). The numerical simulation of low speed and highlift airfoil has been done using ANSYS FLUENT (version 16.2) to obtain drag coefficient, lift coefficient, coefficient of moment and LifttoDrag ratio over the airfoils for the comparative analysis of airfoils. The S1223 RTL airfoil has been chosen as the most suitable design for the specified boundary conditions and the Mach number from 0.10 to 0.30.
Key words: CFD, Lift, Drag, Pitching, LifttoDrag ratio.

INTRODUCTION
An aircraft acted upon by four aerodynamic forces; Thrust, Drag, Lift and Weight. Aircrafts are able to fly due to aerodynamic force produced when a fluid passes over the airfoil. An Airfoil is defined as the cross section of a body that is placed in an airstream in order to produce an aerodynamic force in the most efficient manner possible[1]. If the pressure below the wing is higher than the pressure above the wing, there is a net force upwards and this upward force generates lift.
There are two types of airfoils; symmetric and nonsymmetric. The symmetrical airfoil has identical upper and lower surfaces and the mean camber line and the chord line are same in this airfoil. The Non Symmetric airfoil has different upper and lower surfaces, the surface above the chord line has greater curvature than the surface below the chord line. The nonsymmetrical airfoil can produce useful lift at zero angle of attack, that is why nonsymmetrical airfoil has been used for the analysis in this paper.
In this paper comparative analysis has been done on S1223, S819, S8037 and S1223 RTL airfoils at low mach number. The goal of this paper is to find out that, which airfoil is the most suitable to be used in a low speed aircraft on the basis of their LifttoDrag ratio, Lift coefficient, Drag coefficient and Moment coefficient under the specified boundary conditions and the value of mach number varies from 0.10 to 0.30 with an increment of 0.5.

LITERATURE SURVEY
S. Kandwal et al [1] (Sept.2012), performs the study of twodimensional NACA 4412 airfoil using ANSYS FLUENT for aerodynamic forces on an airfoil. The flow is considered as inviscid and steady with the inlet temperature 288.18 K and Mach number 0.15. After obtaining different parameters from analysis of
NACA4412 at specified boundary conditions, A comparison is done between results obtained from CFD simulation in ANSYS and results obtained from an experiment done in the wind tunnel. The comparative result shows that, the result from experiment and from CFD simulation showing close agreement between these two methods. So, the CFD simulation using ANSYS FLUENT can be used as an relative alternative to experimental method in determining drag and lift.
Jasminder Singh et al [2] (June 2015), Analysis is done on NACA 4412 and Seling 1223 airfoils using Computational Fluid Dynamics. The working conditions are; Inlet velocity of working fluid is 2.5 m/sec, laminar model and 1e03 residual. The Angle of Attack (AOA) varied from 00 to 150 with an increment of 50 but, the area of airfoil, velocity and density of air were kept constant over various angle of attack. They compared lift and drag coefficients for these two airfoils and suggested that, NACA 4412 is suitable for use in Sports Plane, whereas, S1223 is suitable for Heavy lift Cargo Plane due to its high lift.
Mark D. Maughmer et al [3] (June 2001) , They designed PSU94097 airfoil to use on winglets of high performance Sailplanes. For the best performance of the winglet it must be capable to operate in a wide range of Reynold number. That's why, the analysis of this airfoil has been done on reynold numbers from 2.4×105 to 1.0×106 at different angle of attacks. For the validation of design of PSU94097 airfoil, it was tested in the Penn State University Low Speed, LowTurbulence wind tunnel with the rectangular test section having filleted corners and then the results from the wind tunnel testing compared with the highly regarded twodimensional airfoil codes, PROFILE Code and XFOIL. In all respect this airfoil satisfies all the requirements for the design.
M.M.M. Talukder et al [4] (January 2016), They did analysis on S819 and S821 airfoil using Computational Fluid Dynamics method with the wind speed of 3m/sec, 4m/sec and 5m/sec and the angle of attack varied from 40 to 160 with an increment of 40 for a wind turbine blade profile. The results generated from CFD method by using ANSYSFLUENT software were verified experimentally by testing the wooden model of airfoils in an open circuit suction type wind tunnel with a test section of 380 x 355 x 330 mm. They found that, sliding ratio plays a key role in defining the effectiveness and usability of wind turbine blade profile. The optimum angle of attack should be between 00 to 20 to get the maximum value of sliding ratio.
Rong Ma et al [5], Worked to optimize the design of S1223 airfoil for a high efficiency propeller of low dynamic vehicle in near space. To make a lowspeed, highlift airfoil they used direct search optimization algorithm EXTREM and airfoil flow field solver XFOIL. They generated two profiles of airfoil S1223_OPT1 and S1223_OPT2on the basis of S1223 airfoils profile and compared the results obtained after analysis of these three airfoils; S1223, S1223_OPT1 and S1223_OPT2 at M=0.1 at various angle of attack. They found that, S1223_OPT2 meets all the optimization design requirements among these three airfoils.
Mazharul Islam et al [6] (April 2012), Performed comparative analysis of S8037, SG6040, S1210 and NACA0015 for fixedpitch straightbladed vertical axis wind turbine. For the analysis, the velocities of flowing fluid are 5 m/sec, 10 m/sec and 15 m/sec. During the analysis they found that, NACA0015 produces lowest amount of torque. At the end of their comparative analysis they found S8037 as the best airfoil for the applied boundary conditions among all other candidates, both in terms of starting torque and maximum (Cp)net.
Tarun B Patel et al [7] (April 2015), Performed analysis of Lift and Drag forces of NACA airfoils using PYTHON (x,y) 2.7.9.0 programming.. For a wind turbine blade NACA0012, NACA4412 and NACA6412 Airfoil profile is considered for analysis and with the variation in angle of attack and velocity the comparative analysis has been dome on these three airfoils. They concluded that, Maximum lift forces obtained for the NACA4412 airfoil comparative to all other airfoils under specified boundary conditions.
MD. Safayet Hossain et al [8] (October 2014), They did comparative analysis of NACA 6409 and NACA 4412. The pressure distributions as well as liftto drag ratio of these two airfoils were compared with laminar flow velocity of 1 m/sec. After the comparative analysis, it was found that, NACA 4412 aerofoil is more efficient than NACA 6409 aerofoil.

METHODOOGY

Geometry Modeling and Meshing
The ANSYSFluent (version 16.2) was used to analyze the flow analysis of flowing fluid over the airfoils. Imported coordinates of airfoil in Workbench to generate a two dimensional geometry of airfoil and also generated a two dimensional region around the geometry of airfoil as a flow region for flowing fluid and it is assumed that fluid flow in zdirection is negligible. Structured mesh has been generated over the airfoil for better results. In meshing, element size is 7e005 meter and airfoil model discretized into 27957 nodes and 27600 elements.
Fig. 1: Mesh Model

Boundary Conditions
All the boundary conditions are applied with an option of Double Precision for all the airfoils in 2D dimension. SpalartAllamas (1 equation) flow equation has been chosen for the analysis of flowing fluid. The flowing fluid is considered as Air Ideal gas and the flow of the flowing fluid is steady.
For the analysis of airfoils on above mentioned boundary conditions, The Mach Numbers M= 0.10, M= 0.15, M= 0.20, M= 0.25 and M= 0.30 are applied on each of the airfoil for their comparative analysis for various flow conditions as mentioned in the Table1
Table 1: Boundary Conditions for Airfoil Analysis
Inlet Boundary Conditions
Inlet Type
Pressurefar Field
Temperature
288.17 K
Wall Boundary Conditions
Wall Motion
Stationary Wall
Shear Condition
No Slip
Outlet Boundary Conditions
Outlet Type
PressureOutlet
Gauge Pressure
0 Pa
Solution Method
Scheme
Simple
Gradient
GreenGause Cell Based
Initialization Method
Hybrid Initialization
Properties of Flowing Fluid
Fluid
Air Ideal Gas
Specific Heat
1.006.43 J/kgk
Viscosity
1.7894e05 kg/ms


RESULT AND DISCUSSION

LifttoDrag Ratio
The higher value of this ratio is the most desirable factor for any aircraft design. The value of dragtolift ratio can be increased either by increasing the value of lift coefficient or by decreasing the value of drag coefficient but in case of an aircraft design the lift directly depends on the weight of an aircraft and the drag depends upon the aerodynamic design of aircraft and its wings.
Mach Number
Cl/Cd at S1223
Cl/Cd at S819
Cl/Cd at S8037
Cl/Cd at S1223 RTL
0.10
0.5020
0.1632
0.1252
0.8215
0.15
0.5989
0.1745
0.1737
0.9239
0.20
0.7801
0.1697
0.2122
1.0911
0.25
0.9323
0.1803
0.2608
1.3437
0.30
1.0485
0.2254
0.3005
1.4143
Mach Number
Cl/Cd at S1223
Cl/Cd at S819
Cl/Cd at S8037
Cl/Cd at S1223 RTL
0.10
0.5020
0.1632
0.1252
0.8215
0.15
0.5989
0.1745
0.1737
0.9239
0.20
0.7801
0.1697
0.2122
1.0911
0.25
0.9323
0.1803
0.2608
1.3437
0.30
1.0485
0.2254
0.3005
1.4143
Table 2: Lift to Drag Ratio of airfoils at different Mach number
1.5
Cl/Cd at
S1223
1.5
Cl/Cd at
S1223
1
1
0.5
Cl/Cd at
S819
0.5
Cl/Cd at
S819
0
M= M= M= M= M=
0.10 0.15 0.20 0.25 0.30
Cl/Cd at
S8037
0
M= M= M= M= M=
0.10 0.15 0.20 0.25 0.30
Cl/Cd at
S8037
Fig. 2: LifttoDrag ratio verses Mach Number
1.40E03
1.20E03
1.00E03
8.00E04
6.00E04
4.00E04
2.00E04
0.00E+00
M= 0.10
M= 0.15
M= 0.20
M= 0.25
M= 0.30
Cd at S1223
Cd at S819
Cd at S8037
Cd at S1223 RTL
It can be seen from the table 2, that the value of lift to drag ratio for all the airfoils increases from mach number 0.10 to
0.30 except S819 airfoil. From the table 2 and figure 2, It can also be observed that, At the mach number of 0.10 the lift to drag ratio for S8037 is lesser than that of S819 but as the value of mach number increases, the value of lift to drag ratio of S8037 gets higher than that of S819 as the value of mach number crosses mach number 0.15.
At all the mach number the highest value of lift to drag ratio is shown by S1223 RTL and the S1223 came on the second place in order to show the highest value of lift to drag ratio.

Drag Coefficient
The drag force acts in the opposite direction of the moving object in a medium of a fluid. It not only opposes the motion of an object in a medium of fluid but also reduces its lift. The drag depends on the density of the fluid, velocity of flowing fluid or an object, compressibility and viscosity of flowing fluid or a fluid around a moving object and the size and shape of the object.
The Coefficient of Drag is a dimensionless quantity, used to evaluate resistance of an moving object in a fluid.
Fig. 3: Coefficient of Drag verses Mach Number
From the table 3, it can be observed that, for all the airfoils value of drag coefficient continuously decreases from 0.10 to 0.30 mach number. From figure 3, it is clear that, S1223 produces the highest amount of drag coefficient at various mach number whereas, S8037 produces lowest amount of drag coefficient at various mach numbers. S1223 RTL airfoil is at the second lowest position in the graph of drag coefficient at different Mach numbert.

Lift Coefficient
Lift is a mechanical aerodynamic force that is generated by a solid object passing through a fluid and this force opposes the weight of flying object and hold it in the air. It is a vector quantity and it acts through the center of pressure of the flying object. The lift is generated by the difference in velocity of flying object and fluid, around that flying object. It takes no difference whether the object is passing through the fluid or the fluid is flowing over an object.
Lift is generally measured as a non dimensional coefficient, Coefficient of Lift (CL).
1 2
1 2
CL = L
( )V A
Cd = D
(12)V2S
where, D is the Drag Force, is the density of fluid, V is the flow speed and S is the reference area.
M
Cd at S1223
Cd at S819
Cd at S8037
Cd at S1223 RTL
0.10
1.3081e03
9.8369e04
8.0983e04
8.9127e04
0.15
1.1230e03
9.0849e04
6.9873e04
8.0573e04
0.20
9.0151e04
7.6648e04
6.2320e04
7.0369e04
0.25
7.6365e04
6.4104e04
5.0989e04
5.6043e04
0.0
6.9307e04
5.3447e04
4.5262e04
5.3113e04
M
Cd at S1223
Cd at S819
Cd at S8037
Cd at S1223 RTL
0.10
1.3081e03
9.8369e04
8.0983e04
8.9127e04
0.15
1.1230e03
9.0849e04
6.9873e04
8.0573e04
0.20
9.0151e04
7.6648e04
6.2320e04
7.0369e04
0.25
7.6365e04
6.4104e04
5.0989e04
5.6043e04
0.30
6.9307e04
5.3447e04
4.5262e04
5.3113e04
Table 3: Coefficient of Drag at different Mach Number
2
where, L is the Lift Force, is the density of fluid, V is the flow speed and A is the relative plan area.
Here, Coefficient of Lift is generated over the four airfoils at different Mach Number to perform a comparative analysis of these airfoils on the basis of their predicted lift coefficient at mach number from 0.10 to 0.30. Table 4 and Figure 4 represents the variation of lift coefficient at different mach number for these four airfoils.
Fig. 4: Coefficient of Lift at different Mach Number
M
Cl at S1223
Cl at S819
Cl at S8037
Cl at S1223 RTL
0.10
6.5673e04
1.6060e04
1.0143e04
7.3219e04
0.15
6.7267e04
1.5862e04
1.2138e04
7.4447e04
0.20
7.0332e04
1.3010e04
1.3227e04
7.6784e04
0.25
7.1202e04
1.1564e04
1.3298e04
7.5282e04
0.30
7.2674e04
1.2047e04
1.3604e04
7.5120e04
9.00E04
8.00E04
7.00E04
6.00E04
5.00E04
4.00E04
3.00E04
2.00E04
1.00E04
0.00E+00
M= 0.10
M= 0.15
M= 0.20
M= 0.25
M= 0.30
Cl at S1223
Cl at S819
Cl at S8037
Cl at S1223 RTL
The airfoil moment (M) is dependent on air speed, wing area and the chord length. The aircraft becomes unstable in pitch if the tail is required to produce positive or upward lift in order to balance the sum of the airfoil and lift moments. This would be indicated by a negative or nose down tail moment, and is most likely to occur at low airspeeds.
M
Cm at S1223
Cm at S819
Cm at S8037
Cm at S1223 RTL
0.10
4.9707e07
1.0022e07
9.8015e08
4.8075e07
0.15
5.1371e07
9.8571e08
9.2766e08
4.8140e07
0.20
4.8763e07
8.6479e08
9.1921e08
4.7896e07
0.25
4.6713e07
7.8680e08
6.7716e08
4.7199e07
0.30
4.4377e07
7.8757e08
6.9833e08
4.6982e07
M
Cm at S1223
Cm at S819
Cm at S8037
Cm at S1223 RTL
0.10
4.9707e07
1.0022e07
9.8015e08
4.8075e07
0.15
5.1371e07
9.8571e08
9.2766e08
4.8140e07
0.20
4.8763e07
8.6479e08
9.1921e08
4.7896e07
0.25
4.6713e07
7.8680e08
6.7716e08
4.7199e07
0.30
4.4377e07
7.8757e08
6.9833e08
4.6982e07
Table 5: Coefficient of Moment of airfoils at different Mach number
Fir. 4: Coefficient of Lift verses Mach Number
As Shown in the table 4, For S8037, S1223 and S1223 RTL the value of coefficient of lift keep increasing from mach number 0.10 to 0.30, but in the case of S819 airfoil, the value of coefficient of lift keep decreasing from mach number 0.10 to 0.25 but there is a slight increment in this value at 0.30 mach number
In the figure 4, The value of lift coefficient for S8037 airfoil is lowest at mach number 0.10 but it keep increasing with the increment of mach number and crosses the curve of S819 airfoil at mach number 0.20. The maximum lift coefficient is generated by S1223 RTL among all these four airfoils under the selected range of mach number. The value of lift coefficient for S1223 RTL
airfoil keep increasing from mach number of 0.10 to 0.20
6.00E07
5.00E07
4.00E07
3.00E07
2.00E07
1.00E07
0.00E+00
M= 0.10
M= 0.15
M= 0.20
M= 0.25
M= 0.30
Cm at S1223
Cm at S819
Cm at S8037
Cm at S1223 RTL
but for further increment of mach number lift coefficient keep decreasing up to the value of mach number 0.30. S1223 airfoil, produces the second highest amount of lift coefficient that can be seen in the table 4 and figure 4. For S1223 airfoil, the value of coefficient of lift keep increasing from mach number 0.10 to 0.30.

Coefficient of Moment
Coefficient of Moment is generally known as the Pitching Moment. Pitch Moment is a moment which acts on the pitching axis of the object moving or traveling through a fluid medium. It acts on the aerodynamic center of the airfoil rather than the center of pressure of the airfoil. The pitching moment coefficient does not vary with the lift coefficient on the Aerodynamic Center. The aerodynamic center of an airfoil is generally close to 0.25 times of chord behind the leading edge of the airfoil. The aerodynamic center is a point on the chord line of the airfoil at which the pitching moment coefficient does not vary with the change in angle of attack.
The Pitching Moment (Cm) is generally
defined as;
Cm = P S c
where, M is the Pitching Moment, P is the dynamic pressure, S is the relative plan area and c is the chord length of the airfoil.
Fig. 5: Coefficient of Moment verses Mach Number
From the table 5 and figure 5 it can be observed that, from mach number 0.10 to 0.20 S1223 airfoil showing the maximum moment coefficient than other airfoils, because the moment coefficient of S1223 airfoil keep decreasing from mach number 0.15 to 0.30, the values of moment coefficient of S1223 airfoil gets lower to that of S1223 RTL airfoil at mach number 0.25 and 0.30. The lowest value of moment coefficient has been shown by the S8037 airfoil among all these airfoils.


CONCLUSION
In this paper, analysis of four airfoils has been carried out at low mach number from 0.10 to 0.30 by using SIMPLE scheme, GreenGause Cell Based gradient and Spalart Allamas Model. For the design of an aircraft it is most desirable to choose a suitable airfoil to design its wings. so that, desirable amount of lift can be generated under the specified boundary conditions for that aircraft. For this purpose, LifttoDrag ratio and coefficient of moment are the most important valuesto be taken care for proper aerodynamic design of the wing
From the analysis, It can be concluded that. Airfoil S1223 RTL showing the maximum dragtolift ratio, maximum value of coefficient of lift and the second lowest value of drag coefficient at mach numbers from
0.10 to 0.30. The coefficient of moment or pitching moment is another factor which needs to be considered
while selecting airfoils for the wing of an aircraft. So on the basis of these values at different mach numbers, the selection must be S1223 RTL airfoil among these four airfoils for the specified conditions.

REFERENCES

S.Kandwal and Dr. S. Singh, on "Computational Fluid Dynamics Study Of Fluid Flow And Aerodynamic Forces On An Airfoil", International Journal of Engineering Research & Technology (IJERT), ISSN: 2278 0181, Vol. 1 Issue 7, September – 2012.

Jasminder Singh, Dr. Jaswinder Singh, Ampritpal Singh, Abhishek Rana, Ajay Dahiya, on " Study of NACA 4412 and Selig 1223 airfoils through computational fluid dynamics", SSRG International Journal of Mechanical Engineering (SSRGIJME) volume 2 Issue 6June 2015.

Mark D. Maughmer, Timothy S. Swan, and Steven M. Willits, on " The design and testing of a winglet airfoil for lowspeed aircraft", Journal of Aircraft, Vol. 39, No. 4 (2002), pp. 654661.

M. M. M. Talukder*, M. K. Islam and M. R. Rukan, on " Comparative Aerodynamic Analysis of Wind Turbine Blade Profiles", International Journal of Engineering Research & Technology (IJERT), ISSN: 22780181, IJERTV5IS010031, Vol. 5 Issue 01, January2016.

Rong Ma, Bowen Zhong, Peiqing Liu, Dimitris Drikakis, on " MultiObjective optimization design of low Reynolds Number airfoils S1223", 27th International congress of the aeronautical sciences.

Mazharul Islam, Rupp Carriveau and Amir Fartaj, on " Performance analyses of a fixedpitch straight bladed VAWT with selected low Reynolds number airfoils", International Journal of Environmental Studies, Vol. 69, No. 2, April 2012, 289298.

Tarun B Patel, Sandip T Patel, Divyesh T Patel and Maulik Bhensdadiya, on " An Analysis of Lift and Drag Forces of NACA Airfoils Using Python", International Journal of Application or Innovation in Engineering & Management (IJAIEM), ISSN 2319 – 4847, Volume 4, Issue 4, April 2015.

MD. Safayet Hossain, Muhammad Ferdous Raiyan, Mohammed Nasir Uddin Akanda and Nahed Hassan Jony, on " A comparative flow analysis of NACA 6409 and NACA 4412 Aerofoil", International Journal of Research in Engineering and Technology, eISSN: 23191163, pISSN: 23217308, Volume: 03, Issue: 10, Oct2014.