v1 Covered
v1 Covered
v1 Covered
Research Article
Keywords: Computation uid dynamics, Aerodynamic characteristics, Wing planforms, Taper ratio, Aspect
ratio, Sweptback wing, Angle of Attack
DOI: https://doi.org/10.21203/rs.3.rs-1806010/v1
License: This work is licensed under a Creative Commons Attribution 4.0 International License.
Read Full License
To investigate and compare the wing planform’s effect on the
aerodynamic parameters of aircraft wings using
computational fluid dynamics (CFD)
Abstract
A high lift wing’s design is a critical aerodynamic parameter of aircraft. Research has been
performed to enhance the wing performance by altering the wing planform, out-of-plan
transformation, and airfoil modification. Herein, the wing planforms, airflow velocity (5 – 35 m/s),
and their effect on aerodynamic parameters as a function of an angle of attack (-20° - +50°) are
investigated numerically. The NACA 0012 profile 3D wing with varying aspects, taper ratios, and
sweptback angles is investigated using computational fluid dynamics (CFD). The study exposed
that, a rectangular wing has a high lift and drag coefficient among all the wing planforms. It also
reveals that there will be no lift generation for symmetric airfoil at an angle of attack (AoA) of 0°.
Furthermore, the effect of changing the leading-edge sweptback angle (LESA) on the lift and drag
is prominent at higher AoA. It changes airflow from wing root to wing tip however, it also
increases the tip vortices. The lift and drag increase with the increase of flow velocity and taper
ratio (TR) and becomes more prominent after an AoA of 9°. Similarly, with the increase of aspect
ratio (AR), lift and drag decrease to a certain limit (in this case AR 10), beyond this limit, it starts
increasing with an increase in AR and its effect is more visible at an AoA of 15°. The TR and LESA
are the best parameters for the aircraft wing to get better flight conditions depending upon the
applications.
1. Introduction
The wing is made up of an airfoil that generates lift. The design of low-weight and high-lift
generating wings is one of the most important needs of the modern world. An efficient wing can
be generated by optimizing its geometrical (wing planforms) and aerodynamic parameters (lift
and drag). It has different configurations and geometries depending upon the application. With
the increase of aircraft speed, the air friction around the wing also increases. Due to this, the
aerodynamic parameters like turbulence and drag increase, and hence cause an imbalance of the
2. Methodology
2.1. Planforms of wing
A planform is the top view of a wing that interacts with three-dimensional airflow. The wing
planforms that are critical to optimizing the overall aerodynamics parameters of the wing are
LESA, TR, AS, AoA, the geometry of the airfoil, and its configuration and airflow velocity. The
ratio of wingspan to its chord length is expressed as an aspect ratio. Although increasing the aspect
ratio helps reduce drag, it also increases the weight of the wing. The ratio of the airfoil tip chord
(Ct) to its root chord (Cr) is known as the taper ratio. Tapering the wing can also help to reduce
drag as well as the wing’s weight. Sweep back is the wing's rearward inclination from root to tip
of the wing [15]. The LESA of the wing can also enhance the aerodynamic parameter by shifting
the flow pattern from root to tip. The various parameters of the airfoil are given in figure 2.
Sweep
Planform Span
Morphing wing types
chord
Twist
Camber
Airfoil
Thickness (b)
(a)
Figure 1. (a) The various types of morphing in the aircraft wing (b) The various parameters of the airfoil.
(a) (b)
Figure 2. The various types of wing planforms (a) Sweptback of the rectangular wing (b) Sweptback of
the wing with a taper ratio of 0.5 (c) The taper ratios of the wing (d) The aspect ratios of the wing.
The main aerodynamic parameters of the wing are the lift (Cl) drag (Cd) coefficients, lift
(Fl) and drag (Fd) forces and the airfoil efficiency (Cl/Cd). A fluid flowing around an airfoil
can exerts a force on it. This force can be resolved into two components, lift force and drag
force as shown in figure 4. Fl is the component that is normal to the fluid flow. The airfoil
shape is design so that it generates a pressure difference between up and lower surface
flow, due to which lift is generated. Fl is acting upward to counter the effect gravity,
although it can be acting in any direction normal to the fluid flow. In contrast, the Fd is the
component that is acting parallel to the fluid flow. Drag force is acting in opposite direction
to the fluid velocity. The Cd is the sum of zero-lift drag (Cdo) and the induced drag (Cdi).
The ratio of lift (Cl) to drag coefficient (Cd) determines the aerodynamic performance of the
wing. The Cl /Cd shows the amount of lift generated by an airfoil to its drag. The Cl/Cd ratio
represents the efficiency or performance of the airfoil. The Cl/Cd determine the stall angle
where the wing efficiency is maximum.
2𝐹𝐹𝑙𝑙
𝐶𝐶𝑙𝑙 = (3)
𝜌𝜌𝜌𝜌𝑉𝑉 2
3. Numerical analysis
The main aim of the investigation is to design wing which generates high lift and low drag to
reduced fuel consumptions. For this purpose, the effect of wing planforms mostly TR, AR and
LESA edge are investigated numerically. For the outer profile of the wing, NACA 0012 airfoil was
considered in this study with chord of 0.3 m and span of 1.1 m. However, these investigations are
carried out at a lower Re of 20.54 × 103. The air density (ρ) is 1.225 kg/m3 and its dynamic viscosity
(μ) of 1.7894 × 10−5 kg/m-s. The fluent turbulence model of k-ω SST with two-equation was used.
The k-ε model is the turbulent flow model and assumes that the molecular viscosity effect is
minimal. The transport equation was derived by using the assumption. This model is valid for
fully developed turbulent flow only. While the SST K-ω assumes that the turbulent flow viscosity
is changing. The turbulent viscosity of fluid is varying to account for the transfer of the primary
turbulent shear stress. By comparing the SST K-ω model to the K-ω model, this changing flow
viscosity gives an advantage to SST K-ω over both the standard k-ω model and the standard k-ε
model. The transport equation of SST k − ω model is given by [10]:
𝜕𝜕 𝜕𝜕 𝜕𝜕 𝜕𝜕𝜕𝜕
(𝜌𝜌𝜌𝜌) + (𝜌𝜌𝜌𝜌𝑢𝑢𝑑𝑑 ) = (𝛤𝛤𝑘𝑘 ) + 𝐺𝐺𝜕𝜕 − 𝑌𝑌𝜕𝜕 + 𝑆𝑆𝜕𝜕 (4)
𝜕𝜕𝜕𝜕 𝜕𝜕𝑥𝑥𝑖𝑖 𝜕𝜕𝑥𝑥𝑗𝑗 𝜕𝜕𝑥𝑥𝑗𝑗
𝜕𝜕 𝜕𝜕 𝜕𝜕 𝜕𝜕𝜔𝜔
(𝜌𝜌𝜌𝜌) + (𝜌𝜌𝜌𝜌𝑢𝑢𝑑𝑑 ) = (𝛤𝛤𝜔𝜔 ) + 𝐺𝐺𝜔𝜔 − 𝑌𝑌𝜔𝜔 + 𝐷𝐷𝜔𝜔 + 𝑆𝑆𝜔𝜔 (5)
𝜕𝜕𝜕𝜕 𝜕𝜕𝑥𝑥𝑖𝑖 𝜕𝜕𝑥𝑥𝑗𝑗 𝜕𝜕𝑥𝑥𝑗𝑗
Where effective diffusivity of K and ω are ΓK and Γω the respectively. Yk and Yω represent the
dissipation of K and ω due to turbulence and Dω represents the cross-diffusion term. Sk and Sω are
user defined source terms. GK represents the generation of turbulence kinetic energy due to mean
velocity gradients and Gω represents the generation of ω. The boundary condition for the
investigation is:
S. No Chord length (m) Taper ratio Span length Sweptback angle (ᴧ)
Root Tip (TR) (m)
chord chord
1 0.3 0.3 1 1.1 0°
2 0.3 0.3 1 1.1 5.2°
3 0.3 0.3 1 1.1 10.3°
4 0.3 0.3 1 1.1 15.3°
5 0.3 0.3 1 1.1 20°
𝐶𝐶𝜕𝜕
𝜆𝜆 = (6)
𝐶𝐶𝑟𝑟
Results
Lift and drag Forces are dependent upon the density (ρ) of the fluid in which the airfoil is moving,
area (A) of the airfoil, the square of the fluid velocity (V2) and Lift (Cl), drag (Cd) coefficients. Lift
(Cl), drag (Cd) coefficients and the airfoil performance (Cl/Cd), for the wing of 1.1 m span and 0.3
m chord length are computed using k-ω SST turbulence model of Fluent at various wing
planforms. The results are plotted as function of AOA. These investigations are carried out at a
(a) (b)
(c)
(d)
Figure 5. The drag coefficient (Cd)with angle of attack (a) The Cd of the rectangular wing (b) The Cd of
rectangular wing with taper ratio of 0.5 (c) The Cd of taper wing (d) The Cd of aspect ratio (e) The Cd
with varying velocity.
The effect of various wing planforms and flow velocity on the aerodynamic lift coefficient (Cl) as
function of AoA are shown in figure 3. The Cl increases with increase in AoA up to stall angle of
airfoil. Beyond stall angle, the Cl decreases sharply due to the flow separation form the airfoil
separation [16].The results show that, at 0° AoA, there is no lift generation for symmetric airfoil as
NACA 0012 used in the study. There is no pressure difference generated between the upper and
lower surfaces of the wing. The Cl for rectangular wing versus AoA with LESA is presented in
Figure 7(a). The Cl increase with increasing of the sweptback angle of the wing. The sweptback
angle of 20°, gives the maximum lift
coefficient but is more effective at higher AOA of 22°. The rectangular wing has no LESA with TR
0.5 gives the maximum Cl. Figure 7(b) shows the combine effect of wing with LESA and TR of 0.5
on Cl. The sweptback wing with TR has no more significant effect on Cl of wing. The TR has
significant effect on the Cl of wing and give maximum efficiency at an AoA of 6° as presented in
figure 8(c). With the increase of the TR, the Cl increases. However, another advantage increase in
TR is reduction of wing weight. The thin larger wing has higher Cl. The Cl decreases with increase
of AS up to certain limit (AR 10). After that, further increase in aspect ratio can increase the Cl as
given in figure 7(d). But there is also an increase in weight of the wing occur. The Cl also increases
with the increase of flow velocity as given in figure 7(e).
(d)
(c)
(e)
Figure 6. The drag coefficient (Cl)with angle of attack (a) The Cl of the rectangular wing (b) The Cl of
rectangular wing with taper ratio of 0.5 (c) The Cl of taper wing (d) The Cl of aspect ratio (e) The Cl with
varying velocity.
(a) (b)
(c) (d)
2. Conclusions
An investigation on the effects of wing planforms on the aerodynamic characteristics of
NACA0012 airfoil wing with a chord length of 0.3 m and span length of 1.1m has been carried out.
ANSYS FLUENT with two-equation k-ω SST turbulence model is used for the RE of 20.54 × 103.
The major conclusions of the investigation are:
1. The rectangular wing has high lift and drag coefficient with no or little vertex indices.
2. As NACA 0012 airfoil is symmetric. So, for the symmetric airfoil, there is no lift generation at
an
3. angle of 0°. This is due to the same flow on the upper and lower surface of wing. Hence no
pressure difference is generated due to no lift is produced.
4. The sweptback wing has nearly the same or a little effect on the values of lift and drag
coefficient. However, its effect is more visible at higher angle of attack (30°).
5. Taper ratio is most effective planform parameter which nor only improve the aerodynamic
characteristics to a greater extend but there is also weight reduction. With the increase of taper
ratio, the aerodynamic characteristics increases up to stall and vice versa. However, it is effect
on drag coefficient is more visible an AOA of 6.2° while for lift coefficient is at an AOA of 2°.
6. Aspect ratio has also tremendous effect on the aerodynamic characteristics of wing. With the
increase of aspect ratio, the aerodynamic characteristics decreases up to an AR of 10. After
that, an increase in aspect ratio, increase the aerodynamic characteristics. But there is also an
increase in weight of the wing occurs with the increase of aspect ratio.
7. The analyzed wing has a stall angle of 6. The efficiency of the wing is increasing and becomes
maximum at this angle. Prior to stall angle, the value of Cl is greater than the Cd, however its
value is decreases after the stall angle.
Supplementary Materials: All the figures, graphs and image are the copyright of this paper
except figure 1 (b) and 2. It has been taken from the paper which is cited there.
3. References
[2] H. Zaini and N. I. Ismail, “A Review of Morphing Wing(综述),” Mech. Eng. Colloq. (MEC
2016), no. August, pp. 110–115, 2016.
[5] S. P. Setyo Hariyadi, Sutardi, Sukahir, and Jamaludin, “Characteristics of Separation and
Induced Drag in the Use of Swept-back Wing Unmanned Aerial Vehicle,” J. Phys. Conf.
Ser., vol. 2117, no. 1, 2021, doi: 10.1088/1742-6596/2117/1/012013.
[6] İ. H. Güzelbey, Y. Eraslan, and M. H. Doğru, “Effects of Taper Ratio on Aircraft Wing
Aerodynamic Parameters: A Comperative Study,” Eur. Mech. Sci., vol. 3, no. 1, pp. 18–23,
2019, doi: 10.26701/ems.487516.
[7] Tejas Dhekane and Dr. Nitin Sherje, “CFD Simulation of Different Taper Ratio Wings,
Performing Trade-off Assessment and Development of A New Methodology to Plot Lift
Distribution Curve and 3D Local Coefficient of Lift Distribution Graph,” Int. J. Eng. Res.,
vol. V9, no. 04, pp. 517–535, 2020, doi: 10.17577/ijertv9is040471.
[10] Z. Siddiqi and J. W. Lee, “A computational fluid dynamics investigation of subsonic wing
[13] P. Hospodář, A. Drábek, and A. Prachař, “Aerodynamic Design and Strength Analysis of
the Wing for the Purpose of Assessing the Influence of the Bell-Shaped Lift Distribution,”
Aerospace, vol. 9, no. 1, 2022, doi: 10.3390/aerospace9010013.
[16] L. D. Ratio, “4 Lift / Drag,” Textbook, pp. 4–10, 2008, [Online]. Available:
https://www.google.com/url?sa=t&rct=j&q=&esrc=s&source=web&cd=&ved=2ahUKEwjvl
pHQxuHqAhUS4zgGHT5jBZsQFjABegQIMxAF&url=https%3A%2F%2Ffanyv88.com%3A443%2Fhttps%2Fwww.nordian.n
et%2FREPOSITORY%2F117_jaa_principles_of_flight_demo.pdf&usg=AOvVaw1SMFIvOc
TFxl3R1GxdyQf2.