Numerical and experimental study of the impact on aerodynamic characteristics of the NACA0012 airfoil

: Using computational models and low-speed wind tunnel tests, the aerodynamic characteristics of the NACA 0012 airfoil with low Re numbers of (8 × 10 4 , 2 × 10 5 , 3 × 10 5 , and 4 × 10 5 ) and angle of attack (AOA) ranging from 0° to 18° by two steps are examined. Using the same 3-D wind tunnel dimensions, numerical simulations were run. The software program ANSYS FLUENT was used to solve the mathematical model using the continuity equation, the Navier ‒ Stokes equations, and the k – ω shear ‒ stress transport turbulence model. Findings demonstrate that at all AOAs, there is a direct relationship between Reynolds numbers (Re), lift and drag coe ﬃ cients, kinetic energy, and stall angle. The lift coe ﬃ cient rises linearly as the AOA increases, peaking at 14°, the stall angle at higher Reynolds number. The lift coe ﬃ cient was found to decline when the AOA was increased further, reaching its minimal value at an AOA of 18°. With a greater AOA, the airfoil ’ s drag coe ﬃ cient rises, creating turbulent ﬂ ow. The eddies produced by the turbulence cause the ﬂ ow to start separ-ating from the airfoil surface as turbulence increases. As a result, the airfoil lift coe ﬃ cient drops, and its drag coe ﬃ - cient rises at the same time, leading to poor performance. The validation of the numerical results through wind tunnel experiments provided con ﬁ dence in the ﬁ ndings of the study.


Symbols and acronyms
Aerodynamics is the branch of aeronautical mechanics concerned with the study of the forces exerted by air and other gases affecting object motion.In military and civil fields, structures and equipment such as gliders, wind turbines, small aerial vehicles, and unmanned aerial vehicles are subject to low Reynolds (Re) numbers.In aerodynamics, an airfoil is a cross-sectional structure with a rounded leading edge and a curved surface tapering in the flow direction.The airfoil generates a lifting force at an angle perpendicular to the fluid flow and a drag force along the fluid flow direction.These dynamic properties are important in flight and power generation; hence, the airfoil gains vital scientific importance in fluid dynamics.Examples of airfoils are cross-sections of wings, fan blades, windmill blades, compressors, turbine blades in jet engines, water cylinders, aircraft vertical stabilizers, submarine fins, rotary wings, and certain fixed wings.Thin and streamlined airfoils, owing to their low lift and drag, are generally preferred in high-velocity aircraft, while thick airfoils with high lift and drag are often used in aircraft for carrying heavy weights.Airfoils can be divided into symmetrical and asymmetrical types in terms of their surface shape.Symmetrical airfoils have an upper surface and a lower surface of similar camber, whereas asymmetric airfoils have an upper surface that differs from the lower surface.Lift and drag have been widely studied to obtain accurate results for various engineering applications.Sadikin et al. [1] numerically studied the dynamic properties of a symmetric NACA0012 airfoil at angles of attack (AOAs) ranging between −10°and 15°and Re of 3 × 10 6 using three perturbation models, namely, the Spalart- Allmaras, k-Realizable, and k-ω shear-stress transport (SST) models.They concluded that the Realizable k-ε turbulence model results in smooth airflow over the upper surface of the airfoil.This delayed flow separation contributes to an increase in the lift force.In contrast, the Spalart Allmaras and k-ω SST turbulence models exhibit earlier flow separation due to the presence of a large separation bubble.Furthermore, the study found that increasing the AOA leads to a linear increase in the lift coefficient until the stall angle is reached.However, as the AOA continues to increase, the flow velocity decreases at the trailing edge of the upper surface of the airfoil.This decrease in flow velocity causes a reduction in lift force.Shabur et al. [2] numerically examined the turbulent kinetic energy (TKE) distribution of a symmetrical NACA0012 airfoil with a chord length of 1 m using different AOAs and Re numbers and compared the specimen with another symmetrical but thicker NACA0018 airfoil.They found a linear relationship between TKE distribution and AOA, indicating that higher speeds result in increased disturbance.The NACA 0018 airfoil was observed to generate more turbulence compared to the NACA 0012 airfoil due to increased friction with air particles resulting from its thicker design.In aircraft applications, minimizing turbulence is crucial, leading to the recommendation of using the NACA 0012 winglet.Conversely, the NACA 0018 winglet is recommended for wind turbines, where turbulence is not a significant concern.
Martínez-Aranda et al. [3] experimented on a closed wind tunnel device to examine the effect of AOA and Re numbers (3.33 × 10 4 , 6.67 × 10 4 , 10 5 , and 1.33 × 10 5 ) on the NACA0012 airfoil.As the Re numbers and AOA increased, the lift and drag coefficients increased until stall angles ranging between 12°and 14°were reached.
Singh [4] experimentally studied an open-circuit wind tunnel device with a symmetric NACA0012 airfoil composed of acrylic materials and 29 ports.The parameters included AOAs from −14°to +15°increased in 2°increments and air velocities of 6.94, 7.76, 8.5, and 9.82 m/s.The pressure outlets were connected to a pressure gauge that contained gasoline liquid whose density is less than that of water.As air was passed through the test chamber, the pressure on the airfoil surfaces varied, which led to different pressure gauge readings.The values were positive for the airfoil's lower surface but negative for the upper surface, and the generated lifting force increased linearly with increasing AOA until the stall angle of 12°was reached.Furthermore, the lift coefficient presented inverse and direct relationships with increasing Re values at negative and positive AOA values, respectively.Koshy and Jacob [5] numerically studied a symmetric NACA0012 airfoil and compared its pressure and velocity distributions and flow characteristics to those of an asymmetric NACA 2424 airfoil by varying the AOAs at 0°, 6°, and 9°with the Re number set to 3 × 10 5 .The symmetric airfoil did not produce a lift at an AOA of 0°in contrast to the asymmetric airfoil.The difference can be explained by the specific designs of the symmetric and asymmetric airfoils (i.e., with and without a camber, respectively); nonetheless, the pressure and velocity distributions surrounding the airfoil surfaces were similar.The asymmetric airfoil, which produced more lifting force, was more efficient than the symmetric airfoil.
Raval et al. [6] numerically analyzed the flow characteristics of a symmetric NACA0012 airfoil with a chord length of 1 m and selected V = 51.45m/s and AOAs ranging between 10°and 22°for the k-ε turbulence model.They concluded that as the AOA increases, the lift coefficient shows a linear increase until it reaches its maximum value at an angle of approximately 18.4°, which is known as the stall angle.Beyond the stall angle, the lift coefficient rapidly decreases.Additionally, it was observed that as the AOA increases, the drag coefficient gradually increases until reaching the stall angle, after which it increases even more rapidly.This behavior is attributed to the adverse pressure gradient that continuously increases, ultimately leading to flow separation and a loss of lift.
Kumar [7] investigated the flow behavior of a symmetric NACA0012 airfoil with a chord length of 1 m using a k-ω SST turbulence model set to an air velocity of 30 m/s and AOAs of 0°, 3°, 5°, 8°, 11°, 13°, 16°, and 18°.At an AOA of 0°, the pressure and velocity distribution around the airfoil were similar, resulting in no lift generation.As the AOA increased, both the lift and drag coefficients gradually increased until reaching the stall angle of 16°, which represents the optimal performance of the airfoil.However, further increasing the AOA caused the separation of the hydrodynamic boundary layer, creating a turbulence zone filled with vortices.This turbulence negatively impacted the performance of the airfoil, leading to decreased efficiency and potential suspension of its operation.
Mallela et al. [8] numerically analyzed the flow behavior of the NACA0012 airfoil in terms of lift, drag, and pressure and velocity distributions.The Mach number of the k-ω SST turbulence model was 0.44, and the AOAs varied from −6°to 10°in 2°increments.It was observed that at negative AOAs, the pressure increases and the velocity decreases on the upper surface of the airfoil.Conversely, at positive AOAs, the pressure decreases and the velocity increases on the upper surface of the airfoil.It was concluded that the lift coefficient initially has negative values for negative AOAs, but these values increase and become positive as the AOA increases.The lift coefficient reaches its optimal value at an angle of 10°.On the other hand, the drag coefficient decreases at negative AOAs but gradually increases as the AOA becomes positive.
Paper and Muramatsu [9] experimentally studied the effect of AOAs varying from 0°to 15°on the pressure distribution and bubble separation behavior of the NACA0012 airfoil at low Re numbers ranging from 1 × 10 4 to 5 × 10 4 .For the closed wind tunnel device, the airfoil had a chord length of 75 mm, and its 70 pressure ports were distributed on the two surfaces.Based on the pressure distributions, the study concluded that before reaching the stall angle, a short separation bubble appeared at Re numbers of 3 × 10 4 and 5 × 10 4 .However, after the stall angle, there was a transition from a short bubble to a long bubble at Re numbers of 3 × 10 4 and 5 × 10 4 .At Re = 1 × 10 4 , it was observed that the separation bubble moved toward the leading edge with an increasing AOA, indicating the formation of a long bubble on the airfoil surface.In addition, it was concluded that the lift coefficient curve exhibits different behaviors within specific AOA ranges.From 0°to 2°, the lift coefficient curve is zero.From 3°to 4°, it increases rapidly.From 5°to 9°, it increases slightly.From 9°to 12°, it becomes negative.Finally, from 12°to 15°, it remains constant.
Kabir et al. [10] numerically studied the aerodynamic properties of a symmetric NACA0012 airfoil.The k-ω and turbulence models used parameters such as Mach number and AOAs ranging between −15°and 15°in 5°increments.The increasing AOA caused a pressure difference above and below the airfoil's surface, and a vortex was observed at the trailing edge of the airfoil.Additionally, the drag coefficient decreased before the AOA reached 0°, but it gradually increased with increasing AOAs.
Shahariar [11] numerically studied and compared the flow behaviors of two symmetric airfoils, namely, NACA0012 and asymmetric NACA4412, with the same chord lengths of 1 m.In the k-ε turbulence model, the AOAs were varied from 0°to 25°in 5°increments at Re = 3 × 10 5 .A linear relationship was observed between the lift and drag coefficients and the AOAs until the stall angle of 16°was reached.Then, the lift coefficient decreased, while the drag coefficient increased rapidly for both airfoils.The findings were more prominent for the upper and lower surfaces of the NACA4412 airfoil.
Sahoo and Maity [12] numerically analyzed three airfoils with the same chord lengths of 1 m: one of them was symmetrical (NACA0012), whereas the other two were asymmetrical (NACA4412 and S809).Six turbulence models were used, the AOAs varied between 0°and 20°, and the air velocity was set to 51 m/s.They concluded that the (k-ω SST) turbulence model provides results that closely match experimental data and is particularly effective in capturing shear transport.Additionally, it was observed that the NACA 4412 airfoil, with its thin profile and camber height, generates a higher lift coefficient compared to the NACA 0012 and S809 airfoils.
Jha et al. [13] numerically investigated a symmetric NACA0012 airfoil at Re numbers of 2.21 × 10 5 and 2.81 × 10 5 and AOAs of 0°, 5°, 10°, and 20°and validated the find- ings by comparing them with the experimental results.An SAS turbulence model on a simulated adaptive scale was used for the numerical study, while an open-circuit subsonic wind tunnel apparatus was used for the experimental study.The relationship between the lift coefficient and AOAs was linear until the peak of 13°was reached; this finding was consistent with the experimental results.Moreover, the lift coefficient decreased, while the drag coefficient increased with increasing Re numbers.
Yousefi and Razeghi [14] examined the varying effects of Re numbers and AOAs on the turbulent laminar transitions of three symmetric airfoils, namely, NACA0012, NACA0015, and NACA0018.They found that the location of laminar-turbulent transition on an airfoil is significantly influenced by the Reynolds number.As the Reynolds number decreases, the transition point moves closer to the downstream region of the airfoil.For the NACA 0012, 0015, and 0018 airfoils, the lower critical Reynolds numbers for flow were determined to be 1 × 10 5 , 5 × 10 4 , and 3.5 × 10 4 , respectively.Furthermore, it was found that the transition location varies with the AOA.As the AOA increases, the transition point shifts toward the leading edge of the airfoil.
Eftekhari and Al-Obaidi [15] studied the effect of AOAs (0°and 90°) on the aerodynamic properties of an airfoil at low Re numbers (1 × 10 5 , 2 × 10 5 , and × 10 5 ).The findings were validated by comparing them with the results of an experimental study involving a wind tunnel device.The flow separation started upon reaching 5°; at this time, vortices were formed and the lift coefficient decreased.Furthermore, the drag coefficient gradually increased up to an angle of 74°and decreased slightly between 74°and 82°, after which it settled at 90°.The rapidly increasing drag coefficient resulted in a decreasing lift coefficient, which agreed well with the experimental results.
The aerodynamics of airfoil surfaces, particularly the NACA0012 airfoil, have been the subject of numerous earlier investigations and studies.However, other important factors have not been fully examined, including pressure distribution, low Re numbers, and AOA.Many studies concentrated on examining the effects of certain factors, including low Reynolds number or AOA.These previously unconsidered factors are now given importance in the current study, which covers a wider range of low Re numbers (8 × 10 4 , 2 × 10 5 , 3 × 10 5 , and 4 × 10 5 ) and various AOAs (0°-18°with 2°increments).The goal of this study is to examine how they affect the pressure and velocity distribution on the airfoil as well as the lift and drag coefficients.This study offers designers a more complete understanding of airfoil aerodynamics by combining data on lift and drag coefficients as well as pressure distribution.It promotes the design of high-performance, more effective airfoil shapes.

Experimental study
Laboratory wind tunnels are frequently used for experimental investigations on aerodynamic qualities.Laboratory wind tunnels offer an efficient and practical alternative to airborne laboratory wind tunnel equipment that is limited in accessibility while producing reliable data.There are two different types of wind tunnels that are described: closed wind tunnels, where air circulates inside the tunnel, and open-circuit wind tunnels, where air enters and leaves the tunnel from the atmosphere.
The experimental study focuses on investigating the low-speed aerodynamic properties of a NACA0012 airfoil at Re numbers of 8 × 10 4 and 2 × 10 5 using an open-circuit wind tunnel apparatus.The experiment was carried out in the Al-Furat Al-Awsat Technical University/Engineering Technical College's Fluid Mechanics Laboratory in Najaf.The main objective of the experimental research is to validate the numerical results.Wind tunnel parts are depicted in Figure 1.

Stages of the experimental work
As shown in Figure 1, a separate control and instrumentation unit control the wind tunnel apparatus.The dynamic conical form of the structure causes the air to flow and accelerate linearly as it enters via the effuser.The air then goes to the work section, a square area with dimensions of 305, 305, 600 mm and was constructed of acrylic material.The sides of the work section are moveable, while the upper and lower portions are fixed.Pitot tubes attached to two manometers, which are used to detect static pressure, can be inserted through two holes in the upper portion.A special holder for installing the airfoil is located on one side, and a special holder for connecting the other portion of the airfoil to a device that measures angles so that the aileron angle may be adjusted is located on the other side.The air then travels via a grid to the axial fan (the grid safeguards the fan from harm), then it travels through the silencer unit before reaching the atmosphere.
The pressure display unit, which has 32 pressure ports, displays the pressure results for each point of the airfoil.A single-component lift and drag balance unit determines the lift and drag forces for each angle.Bernoulli's equation ( 1) is used to compute the changing air velocity, which is tracked by a separate control and instrumentation unit, while equation ( 2) is used to determine the Re numbers [6]: A symmetrical NACA0012 airfoil was used as the study's model.The airfoil was 300 mm wide and 150 mm in chord.Twenty pressure tappings were installed along its chord as part of its equipment.The tiny diameter tubes that protruded from the airfoil's end were attached to these tappings.The tubes might connect to a bigger diameter pipe by being further attached to labeled flexible pipes using adaptors.By adjusting the opposite end of the airfoil, the AOA of the airfoil may be changed.
The flow was set in static and incompressible states, and the airfoil was in 3D.The turbulent viscosity coefficient was calculated using the k-ω SST turbulence model that consisted of two equations, namely, the conservation equation and two transport equations.The two input variables were the TKE (k) for determining the energy in turbulence and the specific turbulent dissipation rate (ω) for determining the dissipation rate per unit of perturbed kinetic energy.ω is also referred to as the turbulence measure, which is suitable for modeling thick boundary layers for flows with low Re numbers.The k-ω SST turbulence model is one of the most widely used models for determining the effect of turbulent flow conditions, and it is best used for near-wall processing.The k-ω SST turbulence model can also predict excessive and early separation, has a better convergence rate and behavior than other turbulence models for opposite pressure gradients, and requires low-memory computation.
For the four digits of the NACA0012 non-cambered airfoil, 00 indicates the non-cambered symmetric biplane airfoil and 12 indicates the airfoil thickness of 12% chord length.

Mathematical model
The first stage of the numerical study was the mathematical model, which included integral equations and boundary conditions.The 3D continuity equation and momentum equations are given by [16]: Momentum equations: In order to compare the results with the experimental test inside the wind tunnel, the same dimensions of the test tunnel and the available wing model in the laboratory were used.Table 1 shows the input parameters of the airfoil.For the modeling, the symmetric airfoil data of NACA0012 were imported from the website [17] into the ANSYS FLUENT platform.Figure 2 shows the boundary conditions.The boundary conditions are defined as follows: • Inlet section: the velocity of the inflow is set at 8.08, 20.21, 30.31, and 40.42 m/s.• Outlet section is the outlet pressure (P = 0).
• The no slip wall condition is selected.

Mesh configuration
The mesh was built using the ANSYS FLUENT platform to simulate the NACA0012 airfoil.Figure 3 shows the construction and intensification of the mesh around the airfoil, allowing for the visualization of the flow separation while obtaining the most accurate results (Table 2).

Setting up FLUENT
The geometry and mesh data were imported into the ANSYS FLUENT platform for initialization and solver operations.Increasing the kinetic energy would affect the flow of lowviscosity fluids and form vortices of various sizes and thus can be regarded as unstable.Table 3 shows the simulation parameters.Here, "Coupled" was used as a method for pres-sure_velocity coupling, "Least squares cell based" was used for the spatial discretization section, and "First-order Upwind" and "Second-order upwind" were used to solve the momentum and obtain acceptable and good solutions for the turbulent dissipation rate and TKE.The program was initialized before calculation, and then, the calculation was run.

Mesh independence
High-accuracy results of the numerical solution were obtained using several mesh elements, taking into account computer memory and time limitations.In general, large mesh elements require more computer memory and time.Here, the mesh size and its effect on the numerical simulation were verified by testing five mesh element sizes (1,034,566, 1,343,782, 2,246,431, 2,676,481, and 3,476,882) and subsequently obtaining the most accurate results within reasonable time and suitable computer memory.The mesh independence test was performed at Re = 8 × 10 4 , and the ratio (C C / L D ) with AOAs ranging between 0°and 18°w as analyzed for the symmetric NACA0012 airfoil.The test findings presented 2,246,431, 2,676,481, and 3,476,882 as the  three mesh element sizes.These sizes attained almost the same results, but with varying times and computer memory requirements.Finally, 2,676,481 for the mesh element size with the most accurate results and reasonable time and computer memory requirements was selected (Figure 4).

Validation
The numerical results of the published research of Patel and Thakor [7] were examined on the ANSYS FLUENT platform with the objective of ensuring the accuracy and reliability of the current research's results.In particular, the test was conducted by comparing the results of the lift and drag coefficients on the symmetric NACA0012 airfoil with a chord length of 1 m at different AOAs (0°, 3°, 5°, 8°, 11°, 13°, 16°, and 18°) and u s = 30 m/s.The trends in Figures 5  and 6 indicate the high likelihood of adopting the code for the study.The results of the current research are highly compatible with those obtained by the previous work [7], indicating high reliability.The TKE behavior, pressure, and velocity distributions on the surfaces of the NACA0012 airfoil are affected by the changes in both AOAs and Re numbers.Figures 7-9 show the TKE behavior, pressure, and velocity distributions surrounding the upper and lower surfaces of the airfoil.At an AOA of 0°, the pressure and velocity are distributed similarly around the airfoil.When increasing AOA, the pressure increases more prominently at the lower surface than at the upper surface; this trend is inconsistent with the changing velocity distribution, which more prominently increases at the upper surface than at the lower surface.The finding is in line with Bernoulli's principle, which states that fluid pressure rises as velocity decreases, and vice versa.The TKE starts to increase from the trailing edge until the stall angle is reached.However, further increasing the AOA causes the flow to separate and move toward the front edge of the airfoil, which can be explained by the changing pressure distribution.As the pressure increases at the upper surface and decreases at the lower surface, the velocity distribution changes (i.e., it increases at the lower surface and decreases at the upper surface), and the TKE increases along the upper surface of the airfoil.This situation results in the formation of vortices and the occurrence of chaos and turbulence, eventually causing airfoil failure; this phenomenon is called stalling.The pressure and velocity distributions and the TKE values are directly proportional to the increasing Re numbers.

Effect of the AOA on C L and C D
The variations in AOA affect the rise and fall of the airfoil as well as its stopping state.Tables 4 and 5 show the numerical and experimental results for the Cl and Cd for AOAs ranging from 0°to 18°, with a 2°increase at various Re numbers.The NACA 0012 airfoil is surrounded on both of its surfaces by the same pressure and velocity distributions due to its non-cambered form, so the pressure difference is equal to zero.According to equation (7), which represents the calculation formula for the lift force with pressure differences, the airfoil does not produce lift at an AOA of 0°.With increasing AOA, the lift coefficient increases and is directly proportional to the AOA, as illustrated in Figure 10, while the drag coefficient gradually increases until the stall angle is reached (Figure 11).This finding can be explained by the difference in pressure and velocity distributions between the lower and upper surfaces.However, further increasing the AOA causes a sudden and rapid decrease of the lift coefficient, and the drag coefficient exhibits dramatic behavior.This finding can be explained by the formation of vortices that extend from the trailing edge to the leading edge of the airfoil, further increasing the turbulence intensity and causing complete flow separation.Then, the lift coefficient is decreased at high AOAs, whereas the drag coefficient is increased, finally causing airfoil failure: The solution was validated by comparing the numerical results with the experimental results.Figures 12 and 13 show the behavior of C L and C D with AOA at Re = 8 × 10 4 .The C L and C D curves for AOAs of 0°-8°agree well between

Effect of Re on C L and C D
The Re number, a nondimensional number with great importance in fluid mechanics applications, is defined as the ratio of the inertial force to the viscous force.Thus, the Re number can help to determine the relative importance of these forces under certain flow conditions.The effect of Re on the stall angle and the lift and drag coefficients is illustrated in Figures 14 and 15.In Figure 15, the lift coefficient reaches its maximum of 0.896, 0.979, and 1.0188 at the stall angle of 12°and Re numbers of 8 × 10 4 , 2 × 10 5 , and 3 × 10 5 ; it also reaches the maximum at 1.1235 at the stall angle of 14°at Re = 4 × 10 5 , which can be attributed to    The findings indicate that changing the Re number in small increments has a negligible effect on the aerodynamic properties of airfoils because it entails only slight changes in the locations of the perturbed laminar transition.

Effect of AOA and Re on C C / L D
Simply relying on aircraft design is insufficient in generating lift, as a high C C / L D ratio must be gained to increase the flying distance of aircraft, especially when heavy weights are involved.The C C / L D ratio, generally regarded as the efficiency of the airfoil, is the ratio between the lift resulting from the airfoil to the drag resulting from its movement in air.The C C / L D ratio increases with increasing lift coefficient, i.e., it is directly proportional to each of the Re numbers and the AOA until the stall angle is reached.ratio decreases due to the sudden decrease in the lift coefficient and the rapid increase in the drag coefficient.

Conclusion
In the current research, experimental and numerical study was conducted for the symmetrical NACA 0012 airfoil at various AOAs and low Re numbers.The computational analysis was carried out on the ANSYS FLUENT R2 platform using the k-ω SST turbulence model, while the experimental investigation was carried out in a wind tunnel device with dimensions of 305 × 305 × 600 mm.The results of this study may be summed up as follows: • The NACA 0012 airfoil has a non-cambered airfoil shape, which prevents lift at an AOA equal to 0°.Also, between airfoil surfaces, there is no difference in pressure.• The lift coefficient rises linearly with the AOA increased until it reaches its maximum value (C L max) at a certain angle known as the stall angle.However, as the AOA increases further, the flow encounters high oppositepressure gradients it is hard to beat, leading to flow separation and the formation of eddies.This causes the airfoil to experience lower flow velocity and higher turbulence intensity on its upper surface.Consequently, the pressure on the top of the airfoil increases, resulting in a continuous decrease in the lift coefficient.• The drag coefficient on the airfoil increases with the AOA increased, where this increase leads to transitioning the flow from laminar to turbulent.As turbulence intensifies, the flow separates from the airfoil surface due to eddies created by the turbulence.This leads to a decrease in the lift coefficient and an increase in resistance, leading to poor airfoil performance.• As the Reynolds number rises, the free current's disturbance also rises, giving the boundary layer more TKE and causing it to stay connected to the surface for a longer period compared with lower Re numbers.The lift coefficient is somewhat lower at the lesser Reynolds number at higher AOAs due to the trailing edge flow separation.Additionally, increasing the Reynolds number results in an increase in the stall angle (C L max).• Based on the aforementioned findings, additional thorough investigation into the AOA or investigating the connection between asymmetric airfoils and the current study is required, clarifying their applicability.Enhancing airfoil performance with the inclusion of lifting devices, such as flaps and slats, is another area that merits investigation.Researchers in the discipline may want to consider these suggestions as potential future study topics.
Acknowledgements: The authors express gratitude to the Department of Aeronautical Engineering, Technical Al-Furat Al-Awsat Technical University (ATU)/Republic of Iraq, for providing wind tunnel and technical support that enabled successful study completion and the achievement of desired results.They express gratitude for their cooperation and distinction. Abbreviation

Figure 3 :
Figure 3: Mesh for the NACA0012 airfoil and the airfoil domain.

Table 3 : 7 6. 1
Input data for FLUENT simulation Viscous model k-ω SST (two equations) Operating temperature 293 K Operating pressure 0.0 Pa Density of fluid 1.204 kg/m 3 Kinematic viscosity 1.825 kg/m s AOAs 0°-18°steps by 2°N umerical and experimental study of NACA0012 airfoil  Effects of pressure, velocity, and TKE on the NACA0012 airfoil

Figure 5 :
Figure 5: Effect of AOA on the lift coefficient at u s = 30 m/s.Figure 6: Effect of AOA on the drag coefficient at u s = 30 m/s.

Figure 6 :
Figure 5: Effect of AOA on the lift coefficient at u s = 30 m/s.Figure 6: Effect of AOA on the drag coefficient at u s = 30 m/s.

Figure 10 :
Figure 10: Lift coefficient versus AOA at different Re numbers.

Figure 11 :
Figure 11: Drag coefficient against AOA at different Re numbers.

Figure 16 :
Figure 16: C C / L D against AOA at different Re numbers.

Figure 17 :
Figure 17: C C / L D against Re value at AOAs of 10°and 12°.

Table 1 :
Input parameters for the NACA0012 airfoil

Table 2 :
Mesh input parameters