CFD Analysis NACA0012
CFD Analysis NACA0012
Abstract
The airfoils are an important component of many modern-day machines that aerodynamic services
like airplanes, engines, helicopters, windmills, etc. The airfoils are analyzed for aerodynamic
forces like lift, drag, and thrust using experimental and CFD tools. The comparison of the CFD
work of airfoil with the provided experimental work is part of the current study. The CFD analysis
of airfoil NACA0012 was carried out by creating the geometry in the Design Modeler of ANSYS
Fluent, followed by meshing, and preprocessing steps like turbulence models, solver methods,
discretization schemes, and others. The addition of y+ to capture the viscous effects of the flow
near the walls of the airfoil was carried out by calculating the y+ through the given procedure. The
mesh convergence study was also carried out using h-refinement to obtain the mesh-insensitive
analysis. The converged mesh was exported to analysis for all the cases of the angle of attack and
the values of lift coefficient, drag coefficient, and pressure coefficient at given locations on the
airfoil were noted to carry out a comparative study with the experimental results. The airfoil for
the current study was analyzed for the angle of attack ranging from 0º to 20º to sort out the optimum
angle and stall angle. The optimum angle for the airfoil came out to be 7º for the experimental
results where the value of the lift-drag ratio is 15.33, whereas the numerical results came out with
an optimum angle of 7º with a ratio of 13.833. The stall angle for the airfoil resulted in 14º for the
experimental results whereas the lift at this angle is 1.364. For the numerical results, the stall angle
came out to be 15º with the value of the lift as 1.2037. The difference between the experimental
and numerical results is not very much prominent. The pressure coefficient was noted at the
provided pressure-tapping positions and plotted along the walls of the airfoil. The results of
pressure coefficients for different angle of the attack shows that the increment in the angle of attack
of the flow increases the pressure coefficient near the leading edge which predicts that the flow
separation keeps shifting from the trailing edge towards the leading edge when the angle of attack
is increased.
Chapter 1
Introduction
1.1 Airfoil
An airfoil is an object or a structure with a specific cross-section designed to generate the desired
reaction force when it is moved through a gas. The desired reaction forces are lift and drag, while
the gas is air because the application of airfoils exists in the objects that move through the air. The
airfoils are used in designing the cross-sections of airplane wings and stabilizers, the rotor blades
of the helicopter, compressors, propellers, and fans. When an object like a wing or a blade having
an airfoil cross-section is moved in a gas medium like air, it generates an aerodynamic force on
the airfoil (Abzalilov, Il’inskii and Stepanov, 1996). The component of the wing that is
perpendicular to the free stream velocity of the air is called the lift force while the component of
the force parallel to the free stream velocity is called the drag force. An airfoil is designed in such
a way that when it is moved through the air at a high speed, the airfoil shape deflects the air in
such a way the air gets maximum velocity on the upper surface and lower velocity on the bottom
surface. Consequently, the value of the pressure on the upper surface becomes very less as
compared to the pressure on the bottom surface. That high pressure on the bottom surface generates
the lift force that is utilized in aircraft, helicopters, and other applications (Bartaria and Sharma,
2015).
Figure 1 Side view mirror (left) old design, (right) optimized design (Olosson, 2011)
The effects of this drag force directly entangle the efficiency of the car. Therefore, the utilization
of CFD to reduce the drag on the mirrors has optimized the mirror designs, and therefore, modern
manufacturers have developed many efficient models of the car’s side view mirror. Figure 2 shows
the comparative performance of a car’s side-view mirror
Figure 2 Results of total pressure on old and modern side view mirrors (Olosson, 2011)
Chapter 2
CFD Modeling
The CFD modeling of airfoil NACA0012 was carried out by creating the geometry in the Design
Modeler of ANSYS Fluent, followed by meshing, and preprocessing steps like turbulence models,
solver methods, discretization schemes, and others. The addition of y+ to capture the viscous
effects of the flow near the walls of the airfoil was carried out by calculating the y+ through the
given procedure. The CFD modeling of the airfoil for the current study is explained next.
The domain created for the flow of air around the airfoil and the face splitting to obtain a structured
mesh is shown in figure 2.
Figure 4 Domain around the airfoil and face splitting
Therefore,
𝐶𝑓 = 0.058𝑅𝑒𝑙 −0.2
𝐶𝑓 = 0.0046089
Now,
1
𝜏𝑤 = 𝐶 𝜌𝑈 2
2 𝑓 ∞
1
𝜏𝑤 = × (0.0046089) × (1.2) × (32)2
2
𝐾𝑔
𝜏𝑤 = 2.83
𝑚𝑠 2
Also,
𝜏𝑤 2.83
𝑢𝜏 = √ = √
𝜌 1.2
𝑚
𝑢𝜏 = 1.535
𝑠
Hence
𝑦+𝜇 1 × 1.825 × 10−5
𝑦= =
𝑢𝜏 𝜌 1.536 × 1.2
𝒚 = 𝟗. 𝟗 × 𝟏𝟎−𝟔 𝒎
𝒚 = 𝟗. 𝟗 × 𝟏𝟎−𝟑 𝒎𝒎
The height of the first element created in the mesh on the airfoil is also shown in figure 6.
The wall y+ along the walls of the airfoil was also verified in the post-processing and is shown in
figure 7.
Figure 7 Wall y+ along the airfoil
0.023
0.022
0.021
Cd
0.02
0.019
0.018
0 20000 40000 60000 80000 100000
Number of elements
The value of Cd kept on decreasing and became constant at about 60000 elements and kept constant
for further increase in the number of elements. Therefore the mesh created on the airfoil and
domain consists of elements and nodes given in table 3.
Table 3 Configuration of converged mesh
Constraint Setting
Solver type Pressure-based
Time setting Steady-state
Turbulence model k-w SST
Method Simple
Discretization 1st order upwind
Reference values Area = 0.045𝑚2 , Length = 0.15 m,
Depth = 0.3 m
Boundary conditions Inlet: Velocity(magnitude and
direction)
Outlet: Pressure outlet
Residuals 10−5(All)
Report definitions 𝐶𝑑 , 𝐶𝑙
Iterations 20000
1.400
1.200
1.000
0.800
Cl
0.600
0.400
0.200
0.000
0 5 10 15 20 25
AOA
experimental numerical
The prediction of actual results based on CFD results for the current study is shown in the figure
above which shows that the error between the experimental and numerical values is comparatively
less for the initial AOAs while the difference increases in the values of upper AOAs. The
comparison of the drag coefficient for the experimental and numerical study is also given in figure
15.
0.600
0.500
Cd 0.400
0.300
0.200
0.100
0.000
0 5 10 15 20 25
AOA
experimental numerical
The difference between the values of Cd for the experimental and numerical values is very small
in the initial values of AOAs and not significantly increased in higher AOAs which means the
prediction of Cd is more comparative w.r.t Cl.
Experimental results
1.600
1.400
1.200
1.000
Cd,Cl
0.800
0.600
0.400
0.200
0.000
0 5 10 15 20 25
AOA
Numerical results
1.4
1.2
0.8
Cd,Cl
0.6
0.4
0.2
0
0 5 10 15 20 25
AOA
Table 5 shows the value of the stall angle for both experimental and numerical results.
Table 5 Comparison of stall angle
The optimum angle is the angle at which an airfoil gives maximum values of the Lift-Drag ratio.
Therefore, the experimental and numerical results are plotted in figure 18 for the Lift-Drag ratio
to identify the maximum value of the ratio that identifies the value of the optimum angle.
Cl/Cd
18
16
14
12
10
Cl/Cd
0
0 2 4 6 8 10 12 14 16 18 20
AOA
Experimental Numerical
Figure 18 Optimum angle identification for the experimental and numerical results
Table 6 shows the value of the optimum angle for both experimental and numerical results.
Table 6 Comparison of stall angle
Experimental 7º 15.3384
Numerical 7º 13.833
Lower surface
1.20
1.00
0.80
0.60
0.40
Cp
0.20
0.00
-0.20 0.00 20.00 40.00 60.00 80.00 100.00 120.00 140.00 160.00
-0.40
-0.60
Position along chord length (mm)
Experimental Numerical
Upper surface
1.20
1.00
0.80
0.60
0.40
Cp
0.20
0.00
-0.20 0.00 20.00 40.00 60.00 80.00 100.00 120.00 140.00 160.00
-0.40
-0.60
-0.80
Position along chord length (mm)
Experimental Numerical
Upper surface
0.50
0.00
0.00 20.00 40.00 60.00 80.00 100.00 120.00 140.00 160.00
-0.50
Cp
-1.00
-1.50
-2.00
Position along chord length (mm)
Experimental Numerical
Lower Surface
1.20
1.00
0.80
0.60
Cp
0.40
0.20
0.00
0.00 20.00 40.00 60.00 80.00 100.00 120.00 140.00 160.00
-0.20
-0.40
Position along chord length (mm)
Experimental Numerical
Lower surface
1.20
1.00
0.80
0.60
Cp
0.40
0.20
0.00
0.00 20.00 40.00 60.00 80.00 100.00 120.00 140.00 160.00
-0.20
-0.40
Position along chord length (mm)
Experimental Numerical
Upper Surface
0.50
0.00
0.00 20.00 40.00 60.00 80.00 100.00 120.00 140.00 160.00
-0.50
-1.00
Cp
-1.50
-2.00
-2.50
-3.00
-3.50
Length across airfoil (mm)
Experimental Numerical
Upper Surface
0.00
0.00 20.00 40.00 60.00 80.00 100.00 120.00 140.00 160.00
-1.00
-2.00
-3.00
Cp
-4.00
-5.00
-6.00
-7.00
-8.00
Length across airfoil (mm)
Experimental Numerical
Lower Surface
1.50
1.00
0.50
0.00
Cp
-1.00
-1.50
-2.00
Length across airfoil (mm)
Experimental Numerical
Lower Surface
2.00
1.00
0.00
-1.00 0.00 20.00 40.00 60.00 80.00 100.00 120.00 140.00 160.00
-2.00
Cp
-3.00
-4.00
-5.00
-6.00
-7.00
Length across airfoil (mm)
Experimental Numerical
Upper Surface
0.00
-1.00 0.00 20.00 40.00 60.00 80.00 100.00 120.00 140.00 160.00
-2.00
-3.00
-4.00
Cp
-5.00
-6.00
-7.00
-8.00
-9.00
Length across airfoil (mm)
Experimental Numerical
Upper surface
2
0
0 20 40 60 80 100 120 140 160
-2
Cp
-4
-6
-8
-10
Length across airfoil
0 4 7 8 12 14 15
Figure 29 Pressure coefficient on the upper wall of the airfoil for numerical results
Lower surface
2
1
0
-1 0 20 40 60 80 100 120 140 160
-2
Cp
-3
-4
-5
-6
-7
Length across airfoil
0 4 7 8 12 14 15
Figure 30 Pressure coefficient on the lower wall of the airfoil for numerical results
4.4.3 Discussion
The visual distribution of pressure and velocity is shown in the form of contours in figures 31-40
for the important angles of attack. The results show that the flow separation is minimum for the
AOA = 0º and it keeps on increasing and moving from the trailing edge toward the leading edge
when the angle of attack is increased. The pressure coefficient near the leading edge increases
when the angle of attack is increased.
Chapter 5
Conclusions
The CFD analysis of airfoil NACA0012 was carried out by creating the geometry in the Design
Modeler of ANSYS Fluent, followed by meshing, and preprocessing steps like turbulence models,
solver methods, discretization schemes, and others. The addition of y+ to capture the viscous
effects of the flow near the walls of the airfoil was carried out by calculating the y+ through the
given procedure. The mesh convergence study was also carried out using h-refinement to obtain
the mesh-insensitive analysis. The converged mesh was exported to analysis for all the cases of
the angle of attack and the values of lift coefficient, drag coefficient, and pressure coefficient at
given locations on the airfoil were noted to carry out a comparative study with the experimental
results.
The airfoil for the current study was analyzed for the angle of attack ranging from 0º to 20º to sort
out the optimum angle and stall angle. The optimum angle for the airfoil came out to be 7º for the
experimental results where the value of the lift-drag ratio is 15.33, whereas the numerical results
came out with an optimum angle of 7º with a ratio of 13.833. The stall angle for the airfoil resulted
in 14º for the experimental results whereas the lift at this angle is 1.364. For the numerical results,
the stall angle came out to be 15º with the value of the lift as 1.2037. The difference between the
experimental and numerical results is not very much prominent. The pressure coefficient was noted
at the provided pressure-tapping positions and plotted along the walls of the airfoil. The results of
pressure coefficients for different angle of the attack shows that the increment in the angle of attack
of the flow increases the pressure coefficient near the leading edge which predicts that the flow
separation keeps shifting from the trailing edge towards the leading edge when the angle of attack
is increased.
References
Abzalilov, D. F., Il’inskii, N. B. and Stepanov, G. Y. (1996) ‘Airfoil design with external flow
suction’, Fluid Dynamics, 31(6), pp. 814–818. doi: 10.1007/bf02030097.
Al-Obaidi, A. S. M. and Wei, T. C. (2018) ‘Lift and Drag of Non-conventional Wings at
Subsonic Speeds and Zero Angle of Attack - An Experimental Investigation’, MATEC Web of
Conferences, 152, pp. 1–15. doi: 10.1051/matecconf/201815202017.
Bartaria, V. N. and Sharma, S. (2015) ‘Aifoil Terminology, Its Theory and Variations as Well as
Relations with Its Operational Lift Force and Drag Force in Ambient Conditions’, International
Journal of Recent Research in Civil and Mechanical Engineering (IJRRCME), 2(1), pp. 268–
277.
Hossain, M. A., Mashud, M. and Asaduzzaman, M. (2013) ‘FLOW SEPARATION CONTROL
OVER AN AIRFOIL BY USING CO-FLOW JET’, in.
Jameson, A. and Fatica, M. (2005) ‘Using Computational Fluid Dynamics for Aerodynamics’, p.
10. Available at: http://aero-comlab.stanford.edu/fatica/papers/jameson_fatica_hpc.pdf.
Olosson, M. (2011) ‘Designing and Optimizing Side-View Mirrors’, pp. 1–97.
Post, S. and Morris, M. (2010) ‘Force Balance Design For Educational Wind Tunnels’, in. doi:
10.18260/1-2--15891.
Rauf, W. et al. (2020) ‘Effect of Flow Separation Control with Suction Velocity Variation: Study
of Flow Characteristics, Pressure Coefficient, and Drag Coefficient’, Universal Journal of
Mechanical Engineering, 8(3), pp. 142–151. doi: 10.13189/ujme.2020.080302.
Roychowdhury, D. (2020) Computational Fluid Dynamics for Incompressible Flows. doi:
10.1201/9780367809171.