Study On Thin Airfoil Theory Performance Test of E
Study On Thin Airfoil Theory Performance Test of E
Study On Thin Airfoil Theory Performance Test of E
net/publication/324967403
CITATIONS READS
0 4,103
3 authors:
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
A Developed Length Based Product Separating Conveyor for Industrial Automation View project
RFID and Finger Print Based Dual Security System: A Robust Secured Control to Access Through Door Lock Operation View project
All content following this page was uploaded by Md. Mostafizur Rahman Komol on 09 September 2018.
I. INTRODUCTION
Aero-dynamics is the perusal of the dynamics associated II. THEORETICAL REDACTION
with the motions of air over a moving or stationary object.
A. Thin Airfoil Theory
When an aircraft or airfoil accelerates down the runway to
take off, it produces a relative motion between the air and Thin airfoil theory is based on the assumption that under
the aircrafts wings or adjacent body of contact. It is this rip certain conditions an airfoil section may be replaced by its
of airflow over and below the wings, originated by the mean camber line (mcl).
aircrafts forward or backward motion, which produces lift Theoretical Assumption
and consequently makes all forms of powered flight 1. Airfoil is thin η << c;
possible. Insect wings annex magnificent stability and high 2. Angles/slopes are small e.g. sinα≈α, cosα≈1,
load-bearing capacity during flapping flight, glide, taxing slope≈angle;
and hover. Scientists have been caballed by them and have 3. Airfoil only slightly disturbs free stream u', v' << V∞.
performed research for biomimetic applications. Because of
[3]
B. Experimental Observation
Published on April 28, 2018.
Md. Nesar Ali was with Chittagong University of Engineering & If airfoil section of the same mcl but different thickness
Technology (e-mail: nesarali455@gmail.com) functions are tested experimentally at the same angle, it is
Md. Mostafizur Rahman Komol was with Khulna University of
Engineering & Technology (e-mail: komolmostafizur@gmail.com)
found that the lift L and the point application of the lift for
Mohammad Takiuddin Saki was with Bangladesh Institute of Marine
Technology (email: bimt.saki@gmail.com)
DOI: http://dx.doi.org/10.24018/ejers.2018.3.4.665 48
EJERS, European Journal of Engineering Research and Science
Vol. 3, No. 4, April 2018
Therefore, in the thin airfoil theory the mcl (mean camber A. Analytical Test of Elliptical Wing and Model Mosquito
line) is replaced by a vortex sheet of varying strength γ(s) Wing
such that the above conditions are satisfied and our aim is to The most important part of analytical studies is to select
determine this ‘γ’ distribution. an airfoil. We select elliptical wing because further more
studies of aerodynamics equation about finite wing theory
Г = ∫ 𝛾𝑑𝑠 (1) we realize that elliptical wing has greater aerodynamics
efficiency than any other wing. On the other hand, as
proceeding our experiment we fabricate model wing of
mosquito wing to analyze the practical result of insect wing
flight in nature. We also fabricate model of NACA 64A012
mod airfoil for comparing those three values. The fabricated
model wing may have a certain or a non-constant cross-
section across the wing.
Fig. 2. Thin airfoil Approximation There are couple of ways to recite wing airfoil section:
1. Airfoil sketch
The fundamental equation for thin airfoil theory is given 2. Airfoil sorting
below: The sketch and fabrication of the model airfoil is an
intricate and lengthy as well as time consuming process
𝑑𝑍 𝐶 𝛾(𝜉)𝑑𝜉 which needs sagacity in fundamentals of aerodynamics
V∞ (α- ) − ∫0 =0 (for 0<x<c)
𝑑𝑥 2𝜋 (𝑥−𝜉) knowledge. Since the fabricated airfoil further need to be
ascertained by testing it in a wind tunnel [4], [12].
After solving equation γ(ξ) to satisfy; Elliptical wing, Mosquito wing & NACA 64A012 airfoil
a) The camber line is a streamline; Fabrication
b) Kutta condition γ(c) = 0 boundary conditions.
We draw the pattern of Elliptical wing, Mosquito wing &
NACA 64A012 mod airfoil in ‘SOLIDWORS’ as we
considered the top view of that wing as the cross-section of
our wing. At the time of wing modelling & fabrication we
don’t follow any conventional parameters subsequent to the
airfoil geometrical shape. We just sketched and fabricated
the Elliptical wing, Mosquito wing & NACA 64A012 mod
airfoil but all these three wing has demonstrated parameter
like length is ‘171 mm’ and width is ’94 mm’ and the
average height of these two in ‘3.2 cm’. Since those wing
has a lot geometrical shape error but we tried to minimize
that problem as our main target is to replicate those wing
approximately look like an Elliptical wing, Mosquito wing
Fig. 3. Position of Maximum Camber
and NACA 64A012 mod airfoil
DOI: http://dx.doi.org/10.24018/ejers.2018.3.4.665 49
EJERS, European Journal of Engineering Research and Science
Vol. 3, No. 4, April 2018
.
(a)
(b)
Fig. 7 (a,b). NACA 64A012 mod airfoil (model)
B. Simulation Setup
The model was made using the ‘Curves using XYZ point’
function. The SOLIDWORKS pattern with proper
dimension of chord length and span of the wing as shown
below figures.
(b)
Fig. 5 (a,b). Elliptical Wing (model)
(a)
(a)
(b)
Fig. 8 (a,b). Schematic drawing Elliptical wing (model)
(b)
Fig. 6 (a,b). Mosquito wing (model)
(a)
(b)
Fig. 9 (a,b). Schematic drawing Mosquito Wing (model)
(b)
DOI: http://dx.doi.org/10.24018/ejers.2018.3.4.665 50
EJERS, European Journal of Engineering Research and Science
Vol. 3, No. 4, April 2018
(b)
Fig. 10 (a, b). Schematic drawing NACA 64A012 mod Airfoil (model)
DOI: http://dx.doi.org/10.24018/ejers.2018.3.4.665 51
EJERS, European Journal of Engineering Research and Science
Vol. 3, No. 4, April 2018
integrate the pressure distribution to find out the net TABLE III: THEORETICAL LIFT OF NACA 64A012 MOD AIRFOIL
theoretical lift force. For theoretical lift this integration is sib Lift
Obs.
Lap no 1 2 3 4 5 force
with the pressure distribution in the vertical direction while No
(N/m2)
for drag the horizontal pressure distribution is important. Curved Pressure
You can see how the surface “pushes back” with a force Surface Head, mm
8 9 8 6 8
against the pressure as interlude with incoming air, which Pressure
Lower
always acts normal to the surface of airfoil. Surface bent is Surface
Head, 9 9 8 8 7
important two mechanisms are used to generate Lift. This is mm
oft used for subsonic flying inflictions. The another step is Curved Pressure
0.09 0.056 0.09 0.068 0.09 0.061
Surface (N/m2)
to incline the airfoil at an angle allied to horizon, which is
Lower Pressure
usually known as the “relative wind angle”. For low values surface (N/m2)
0.056 0.056 0.09 0.09 0.079
of this angle of attack, the fluid flow remains attached on Pressure
0.034 0 0 -.022 0.011
both upper and lower foil surfaces. For higher angles of Difference
attack flow separation occurred that increased drag force
with corresponding reduces lift force. Eventually, the airfoil
B. Comparison with the Analytical Values of those Airfoil
reaches a stall & halt condition occurred where the pressure
distribution over the airfoil on the top and bottom are equal. Our main goal is to compare the analytical or numerical
value of those prototype Elliptical wing, Mosquito wing &
NACA 64A012 mod airfoil. As we proceeding our project
IV. EXPERIMENTAL DATA work firstly we evaluate the experimental value that we get
from wind tunnel manometric height reading at different
A. Theoretical Lift Force angle of attack 00, 50, 100, 120, 140, 150 respectively. For
The theoretical lift force of NACA 64A012 mod is 0.061 furthermore numerical & statistical analysis we tested those
N/m2 whereas the other two prototype airfoil of Elliptical prototype sketched airfoil pattern in simulation software like
wing & Mosquito wing is 0.0165 N/m2 and 0.005 N/m2. ‘SOLIDWORKS’ for actual value which gives also the
This value indicates that the NACA 64A012 mod airfoil has information of velocity contour around those airfoils. These
the greater flying ability than other two prototype airfoil and velocity contour indicates that how air flow interact with
also indicates that the prototype Elliptical wing is more these airfoils and shearing force acting in this airfoil.
flying efficient than prototype model Mosquito wing.
TABLE IV-A: EXPERIMENTAL RESULTS (ELLIPTICAL WING)
TABLE I: THEORETICAL LIFT FORCE OF ELLIPTICAL WING Angle of Lift, Drag,
L/D CL CD CL/CD
Lift Attack L D
Obs.
Lap no 1 2 3 4 5 force 0 0.028 0.06 0.46630 0.137 0.2938 0.46630
no
(N/m2)
Curved Pressure 5 0.045 0.064 0.70288 0.22 0.313 0.70288
12 10 11 10 10
Surface Head, mm 10 0.095 0.075 1.26667 0.46 0.3673 1.26667
Pressure 12 0.11 0.08 1.35550 0.53 0.391 1.35550
Lower
Head, 10 10 10 9 11
Surface 14 0.12 0.085 1.41004 0.587 0.4163 1.41004
mm
Curved Pressure 15 0.145 0.105 1.38152 0.7101 0.514 1.38152
.135 .113 .124 .113 .113 .0165
Surface (N/m2)
Lower Pressure TABLE IV-B: SIMULATION RESULTS (ELLIPTICAL WING)
0.113 .113 .113 0.102 0.124
surface (N/m2)
Pressure Angle of Lift, Drag,
0.022 0 0.011 0.011 -.011 L/D CL CD CL/CD
Difference Attack L D
0 0.028 0.06 0.46630 0.137 0.2938 0.46630
TABLE II: THEORETICAL LIFT FORCE OF MOSQUITO WING 5 0.045 0.064 0.70288 0.22 0.313 0.70288
Obs. Lap no 1 2 3 4 5 Lift 10 0.095 0.075 1.26667 0.46 0.3673 1.26667
no force 12 0.11 0.08 1.35550 0.53 0.391 1.35550
(N/m2)
14 0.12 0.085 1.41004 0.587 0.4163 1.41004
Curved Pressure 8 9 8 7 12
Surface Head, mm 15 0.145 0.105 1.38152 0.7101 0.514 1.38152
Lower Pressure 8 8 9 9 8
Surface Head, TABLE IV-C: EXPERIMENTAL RESULTS (MOSQUITO WING)
mm
Angle
Curved Pressure 0.091 0.102 0.091 0.079 0.135 0.005 Drag,
of Lift, L L/D CL CD CL/CD
Surface (N/m2) D
Attack
Lower Pressure 0.091 0.091 0.102 0.102 0.102
surface (N/m2) 0 0.001 0.004 0.251 0.00489 0.0195 0.251
Pressure 0 0.011 -.011 -.023 0.033 5 0.0038 0.0042 0.9048 0.00186 0.0205 0.905
Difference 10 0.07 0.0043 16.2791 0.342 0.02106 16.28
12 0.088 0.0052 16.9231 0.4312 0.0254 16.92
14 0.0104 0.007 1.4857 0.0509 0.034 1.486
15 0.0138 0.073 0.18904 0.0675 0.357 0.189
DOI: http://dx.doi.org/10.24018/ejers.2018.3.4.665 52
EJERS, European Journal of Engineering Research and Science
Vol. 3, No. 4, April 2018
DOI: http://dx.doi.org/10.24018/ejers.2018.3.4.665 53
EJERS, European Journal of Engineering Research and Science
Vol. 3, No. 4, April 2018
(a)
(b)
(b)
(c)
(c)
(d)
(d)
(e)
(e)
(f)
Fig. 14. Velocity distribution of NACA 64A012 mod airfoil for angle of
attack (a) 00 (b) 50 (c) 100 (d) 120 (e) 140 (f) 150 respectively
(f)
Fig. 13. Velocity distribution of Mosquito wing for angle of attack (a) 0 0 (b)
50 (c) 100 (d) 120 (e) 140 (f) 150 respectively Normally for negative pressure for upper surface
manometer reading is positive and lower surface manometer
reading is negative. But in this experiment for all tapping
points manometric readings were positive. This occurred for
the whole size and shape that gives rise to manometric
pressure. Therefore, the whole curve shifts upwards while
plotting the values. The lift-to-drag ratio L/D is a measure of
the aerodynamic efficiency of an airplane or the aircraft
particle. It only makes sense that maximum aerodynamics
efficiency should be lead to minimum thrust required. So, it
(a) is important to show the lift-to-drag ratio L/D of those three
airfoil.
DOI: http://dx.doi.org/10.24018/ejers.2018.3.4.665 54
EJERS, European Journal of Engineering Research and Science
Vol. 3, No. 4, April 2018
(i)
(i)
(ii)
Fig. 17. Angle of attack vs lift-to-drag ratio L/D for NACA 64A012 mod
(ii) airfoil (i) experimental result (ii) simulation result
Fig. 15. Angle of attack vs lift-to-drag ratio L/D for Elliptical wing (i)
experimental result (ii) simulation result
From graph we can see that the lift-to-drag ratio L/D
which known as aerodynamic efficiency is greater in NACA
64A012 mod airfoil then the other two model like Elliptical
wing & Mosquito wing. So we can assume that NACA
64A012 airfoil has approximately higher fly ability than
prototype Elliptical wing & Mosquito wing. The results
showed the factual influence of the aforecited unsteady
condition parameters corresponding to the maximum lift and
drag co-efficient. The parameters increase or decrease the
aerodynamic forces with respect to the optimum peak values
depending on the surrounding fluid flow composition [13].
During our experiment we use analog or manual wind tunnel
to test pressure distribution across the airfoil, lift and drag
(i) force of it. Using this analog wind tunnel we have some
experimental and analytical procedure error. On the Other
hand, the prototype model wing that we made is not proper
in dimension and corresponding airfoil geometry. This is
also a reason why we have some experimental error over
numerical value or actual value. Besides, the wind tunnel we
used was worn out (Mach no. 0.3) and the data collection
technique has become very ineffective. Today’s wind tunnel
device uses automatic weight balance, digital force
measurement system and pressure sensors for pressure
measurement and that wind tunnel system get value that
auto plot in the graph. In most of the case they are integrated
(ii) with a computer via software with a printer.
Fig. 16. Angle of attack vs lift-to-drag ratio L/D for Mosquito wing (i)
experimental result (ii) simulation result
VI. CONCLUSION
Wing is the part and parcel of aircraft flying. More
accurate shape and size of airfoil more stable in flying and
clearly showing aerodynamic characteristics. For choosing
NACA 64A012 mod airfoil is that L-39 Albatross lite
combat trainer aircraft of Bangladesh Air Force (BAF) using
DOI: http://dx.doi.org/10.24018/ejers.2018.3.4.665 55
EJERS, European Journal of Engineering Research and Science
Vol. 3, No. 4, April 2018
this type of airfoil. On the other hand, choosing of elliptical [12] Rashed Ahammad. Performance Test of NACA 2412 Airfoil. 2015.
ICMERE Conference Proceeding, 2015.
wing because during era of WW II ‘The Supermarine [13] M.R. Amiralai. An investigation into the effects on unsteady
Spitfire’ which has greater aerodynamic performance at that parameters on the aerodynamics of a low Reynolds number pitching
time using this elliptical wing. Lastly, choosing of mosquito airfoil. 2010. Journal of Fluids and Structure 26 (2010), 979-993.
[14] M.A. Ashraf. Reynolds number, thickness and camber effect on
wing because this type insect flying in nature with the help
flapping airfoil propulsion. 2011. Journal of Fluids and Structures 27
of flapping its wing which create thrust to stable during fly. (2011) 145-160.
Since Lift to Drag is known as the aerodynamic efficiency [15] Jiyu Sun. The structure and mechanical properties of dragonfly wings
factor that’s why we are only showing L/D ratio with and their role on fly ability. 2012. Journal of Fluids and Structure 340,
(2012), 3-17
respect to different angle of attack of those three model
wing like Elliptical wing, Mosquito wing and NACA Md. Nesar Ali was born in Bangladesh. He obtained
64A012 mod airfoil. The test corroborated to the high BSc in Mechanical Engineering from Chittagong
University of Engineering & Technology (CUET) in
lift/drag ratio (Aerodynamic Efficiency) and the high 2017. He has published a number of publication on
maximum lift coefficient (CL) but the airfoil was found to be different Peer reviewed article international research
a little sensitive to leading edge surface roughness. journal (AIP) and also published many conference
proceedings in many international conference which
held in Bangladesh. His field of interest includes
ACKNOWLEDGMENT Aerodynamics, Computational Fluid Dynamics,
Renewable energy, sustainable energy,
The authors would like to be obliged to Chittagong Thermodynamics, Heat Transfer and Mechatronics.
University of Engineering and Technology for providing
laboratory facilities and financial assistance. Md. Mostafizur Rahman Komol was born in
Bangladesh. He obtained BSc in Mechanical
Engineering from Khulna University of Engineering
REFERENCES & Technology (KUET) in 2017. He has published a
[1] http://www.codecogs.com/library/engineering/aerodynamics/basic- number of conference proceedings in many
aerodynamics-theory.php international conference held in Bangladesh. His field
[2] http://web.mit.edu/16.00/www/aec/flight.html of interest includes Control engineering,
[3] Fundamental of Aerodynamics/John D Anderson/Chapter3 Aerodynamics, Mechatronics, Computational Fluid
[4] http://faculty.dwc.edu/.../chapter5,’Daniel Webster College’; Wing dynamics, engineering Mechanics, Thermodynamics
Design Robotics and Renewable Energy.
[5] http://www.airfoiltools.com/airfoil/details?airfoil=n64012a-irfoil
[6] http://www.nasa.gov/audience/k-4/stories/nasa-knows/what-are-wind- Mohammad Takiuddin Saki was born in
tunnel-k4.html Bangladesh. He obtained BSc in Mechanical
[7] http://www.education.com/science-fair/article/physics_experiments- Engineering from Chittagong University of
wind-tunnel/ Engineering & Technology (CUET) in 2010. He is
[8] Md. Nesar Ali. Study on Bird’s & Insect’s Wing Aerodynamics and working now as “Senior Instructor” at Bangladesh
Comparison Its Analytical Value with Standard Airfoil.2017. AIP Institute of Marine Technology (BIMT) since 2016
Conference Proceeding: Volume 1851,020056. and Head of the Department, Shipbuilding
[9] Md. Nesar Ali. Study on flow over finite wing with respect to F-22 Technology, BIMT. He has more than five years of
raptor, Supermarine Spitfire, F-7 BG aircraft wing and analyze its working experience in Ship building industry of
stability performance and experimental values.2017. AIP Conference Bangladesh. His field of interest includes Fluid
Proceeding: Volume 1851,020079. Dynamics, Engineering Mechanics, Thermodynamics, Heat Transfer,
[10] Dr. Md. Mahbubul Alam. A Comparative Studies on various model Automation & Automobile Engineering, Engineering Operation
Bird’s wing & standard airfoil & also track down the most effective Management, Ship Building & Ship Propulsion System, Mechanical and
airfoil using CFD solver software. 2017. ICMERE conference Electrical System of water vehicles and Renewable & Sustainable Energy.
proceeding, 2017
[11] Dr. A.B.M Toufique Hasan. Study on Bird’s flapping wing
aerodynamics and compared its analytical values with NACA 64A012
mod Airfoil. 2016. ICMIEE Conference proceeding 2016.
DOI: http://dx.doi.org/10.24018/ejers.2018.3.4.665 56