Challenge

Flow over an airfoil

Numerical and Experimental Investigation by Flow over NACA0017 airfoil.

Aim:

To determine model flow over a NACA0017 airfoil. For the angle of attacks 0,2,4,6,8,10 compare the lift and drag forces.

Abstract:

Airfoil is the most basic unit of an aircraft wing, tail, or propeller blade the main function of the airfoil is the creation of a major part of the lift and drag of the aircraft as a whole. This airfoil generates a lifting force that acts at right angles to the airstream and a dragging force that acts in the same direction as the airstream.

There are many aero-foil profiles available on the NACA database, these are unique in characteristics. These characteristics are obtained by varying one or more parameters of the aerofoil shown in the below figure.

Airfoil Terminologies:

  • Leading-edge the front edge of an airfoil.
  • Trailing edge—the rearmost edge of an airfoil.
  • The chord line is the line joining the leading and trailing edge of the airfoil profile
  • Blade span—the length of the rotor blade from the center of rotation to the tip of the blade.
  • The Camber line is the curve formed by the locus of the points halfway between the upper and the lower surface.
  • Mean camber line—a line drawn halfway between the upper and lower surfaces of the airfoil. The chord line connects the ends of the mean camber line. Camber refers to the curvature of the airfoil and may be considered curvature of the mean camber line. The shape of the mean camber is important for determining the aerodynamic characteristics of an airfoil section. Maximum camber (displacement of the mean camber line from the chord line) and its location help to define the shape of the mean camber line. The location of the maximum camber and its displacement from the chord line is expressed as fractions or percentages of the basic chord length. By varying the point of maximum camber, the manufacturer can tailor an airfoil for a specific purpose. The profile thickness and thickness distribution are important properties of an airfoil section.
  • Induced flow—the downward flow of air through the rotor disk.
  • Resultant relative wind—relative wind modified by the induced flow.
  • The angle of attack (AOA)—the angle measured between the resultant relative wind and chord line.
  • The angle of incidence (AOI)—the angle between the chord line of a blade and the rotor hub. It is usually referred to as the blade pitch angle. For fixed airfoils, such as vertical fins or elevators, angle of incidence is the angle between the chord line of the airfoil and a selected reference plane of the helicopter.
  • Resultant relative wind—relative wind modified by the induced flow.
  • Center of pressure—the point along the chord line of an airfoil through which all aerodynamic forces are considered to act. Since pressures vary on the surface of an airfoil, an average location of pressure variation is needed. As the AOA changes, these pressures change, and the center of pressure moves along the chord line.
  • Total Aerodynamic force - It is the resultant force of the force components Lift and Drag.
  • Point of Impact - The point at which the air separates to flow about the airfoil. A high-pressure area or stagnation point is formed at the point of impact.

Theory:

The lift force of the airfoil depends on the angle of attack (AOA). Generally, Airfoil requires a positive angle of attack for lift (But a cambered airfoil can generate lift at zero angles of attack as well). Particularly subsonic airfoils are more round leading edge making the airfoil very insensitive to the angle of attack. Supersonic airfoils on the other hand have a sharper leading edge which is very sensitive to the angle of attack. So, in modern-day aviation, airfoils are selected on various parameters each of them optimized for their area of application.

  • Drag coefficient: Resistive force against the body caused by the body moving through the fluid. It is in the opposite direction to the velocity of the flow. The relative velocity between body and flow is called the drag coefficient. `Cd = frac{2Fd}{rho*A*v^2}` `

, where Fd = Drag force

ρ = mass density of fluid,

v = flow speed of the object relative to the fluid.

A = Reference Area

  • Lift Coefficient: Lift Coefficient is a dimensionless coefficient that relates the lift generated by a lifting body to the fluid density around the body, the fluid velocity, and an associated reference area A lifting body is a foil. `Cl = frac{2Fl}{rho*A*v^2}`

    , where FL = Lifting force

    ρ = mass density of fluid,

    v = flow speed of the object relative to the fluid.

    A = Reference Area

Objective: In this report, a detailed study is done on varying angles of attack (0,2,4,6,8,10) and its effect on lift forces (FL) and drag forces (Fd).

Simulation Setup:

    • Initially download the sample NACA file. Then create a new part file and insert the curve by clicking at insert command Curve Curve through XYZ points.
    • A dialog box opens browse the .txt file of the downloaded curve of NACA0017.
    • After inserting the curve, use the convert entity option to convert the curve into the sketch.
    • Use a line command at the trailing edge to make the sketch a close contour.
    • Now use the command of “Move entities” to rotate entities for varying the Angle of Attack (AOA). [So, given in this report the curve(spline) is rotated at angles 0,2,4,6,8,10] w.r.t origin.
    • After setting the angle of attack (AoA)the sketch of this curve is extruded by 0.5meter.
    • Then by setting the Flow simulation wizard we set up an external flow and the fluid that is selected is Air which flows across X-direction over the airfoil with 600 m/s.
    • Automatically a computational domain is created. Now set the computational domain in 2-dimensional in XY plane and adjust (Xmin, Xmax) such that it is 5-10 times of chord length and 10-15 times chord length respectively, and (Ymin, Ymax) with at least 6-8 times of chord length.
    • Now add Global Goals
  • Force in X and Y direction.
  • Frictional in X and Y direction.

Also, add Surface Goals by selecting the face of the airfoil

  • Surface Goals Force in X (drag force)
  • Surface Goals Force in Y (Lift force)
    • Run the simulation
    • After the solver comes to a stalling situation add a cut-plot velocity in XY direction to the airfoil.

Simulation Results:

Case – 1: When the angle of attack (AOA) of Airfoil is 0 degrees.

The Velocity Cut-plot is the lowest at the front side and is maximum at the tail side. This is due to the boundary layer formation across the airfoil. As the (AoA = 0) velocity profiles will be the same on both the sides of wing.

Pressure plot across the airfoil when it is in (AoA = 0) shows as the airfoil is symmetric and the pressure pattern in the above and below the airfoil is the same.

SG Force(X) – It is the drag force that is faced by the airfoil and the average value is 796.15N

SG Force(Y) – It is the lift force that is felt by the airfoil and the average value is -2.9N.

This means at 0 degrees the total Lift force is 0(zero) N. Because of the same pressure found at both the regions of the wing creating no lift force.

  • Case – 2: When the angle of attack (AOA) of Airfoil is 2 degrees.

The Velocity Cut-plot shows that velocity of airflow on the upper side of wing flows faster than the lower side

Pressure plot across the airfoil when it is in (AoA = 2) shows as the airfoil is symmetric and the pressure pattern is seen that there will be high pressure in the windward side of the wing and low pressure in the leeward side about (Pmin, Pmax) = (63661 Pa, 479869 Pa)

SG Force(X) – It is the drag force that is faced by the airfoil and the average value is 834.6N

SG Force(Y) – It is the lift force that is felt by the airfoil and the average value is 2060N.

  • Case – 3: When the angle of attack (AOA) of Airfoil is 4 degrees.

The Velocity Cut-plot shows that velocity of airflow on the upper side of wing flows faster than the lower side

Pressure plot across the airfoil when it is in (AoA = 4) shows as the airfoil is symmetric and the pressure pattern is seen that there will be high pressure in the windward side of the wing and low pressure in the leeward side about (Pmin, Pmax) = (63290 Pa, 473931 Pa)

SG Force(X) – It is the drag force that is faced by the airfoil and the average value is 1036N

SG Force(Y) – It is the lift force that is felt by the airfoil and the average value is 4097.924N.

  • Case – 4: When the angle of attack (AOA) of Airfoil is 6 degrees.

The Velocity Cut-plot shows that velocity of airflow on the upper side of wing flows faster than the lower side

Pressure plot across the airfoil when it is in (AoA = 6) shows as the airfoil is symmetric and the pressure pattern is seen that there will be high pressure in the windward side of the wing and low pressure in the leeward side about (Pmin, Pmax) = (44199 Pa, 480505 Pa)

SG Force(X) – It is the drag force that is faced by the airfoil and the average value is 1401.93N

SG Force(Y) – It is the lift force that is felt by the airfoil and the average value is 6110.92N.

  • Case – 5: When the angle of attack (AOA) of Airfoil is 8 degrees.

The Velocity Cut-plot shows that velocity of airflow on the upper side of wing flows faster than the lower side

Pressure plot across the airfoil when it is in (AoA = 8) shows as the airfoil is symmetric and the pressure pattern is seen that there will be high pressure in the windward side of the wing and low pressure in the leeward side about (Pmin, Pmax) = (52683 Pa, 310903 Pa)



SG Force(X) – It is the drag force that is faced by the airfoil and the average value is 1772.31N

SG Force(Y) – It is the lift force that is felt by the airfoil and the average value is 7825.07N.

  • Case – 6: When the angle of attack (AOA) of Airfoil is 10 degrees.

The Velocity Cut-plot shows that velocity of airflow on the upper side of wing flows faster than the lower side

Pressure plot across the airfoil when it is in (AoA = 10) shows as the airfoil is symmetric and the pressure pattern is seen that there will be high pressure in the windward side of the wing and low pressure in the leeward side about (Pmin, Pmax) = (35088 Pa, 505352 Pa)

SG Force(X) – It is the drag force that is faced by the airfoil and the average value is 2532.83N

SG Force(Y) – It is the lift force that is felt by the airfoil and the average value is 10159.9N.

Relationship between (AoA) vs Drag Force vs Lift Force

The Angle of Attack (AoA)

Drag force (SG Force X) [Newton]

Lift force (SG Force Y) [Newton]

0

796.15N

0N

2

834.6N

2060N

4

1036N

4097.924N

6

1401.93N

6110.92N

8

1772.31N

7825.07N

10

2532.83N

10159.9N

Calculation of Drag Coefficient and Lift Coefficient

Given Formulas:

`Cd = frac{2Fd}{rho*A*v^2}` & `Cl = frac{2Fl}{rho*A*v^2}`

These are formulas used for calculating the Drag coefficient and lift coefficient.

Given Data: = mass density of fluid(air) = 1.2754 kg/m³ [IUPAC standards]

A = Reference area = 1.01m2

v = Relative velocity of body to fluid = 600m/s (X-direction)

The Angle of Attack (AoA)

Drag coefficient

Lift coefficient

0

0.0034

0

2

0.0036

0.0089

4

0.0045

0.0177

6

0.0060

0.0264

8

0.0076

0.0337

10

0.0109

0.0438

Summary:

As the angle of attack (AoA) increases so does the lift force & drag forces and correspondingly lift coefficients and drag coefficients. Thus, this continues to increase until the critical or stalling angle of attack is achieved which is mostly around 15° - 20° for airfoils generally. Although the critical angle depends on various factors such as aircraft weight, load factor, C.G of the aircraft, etc.,

So, therefore, it is very necessary to study these parameters before selecting the correct airfoil for application with a best-case optimized scenario to make the overall system efficiency increased.

Your portfolio, next

Ship work like this
in your own course.

Talk to a career counsellor

Free 20-minute call. We'll map the right course to where you want to work.

By submitting you agree to our Terms and Privacy Policy.

Talk to a career counsellor

Free 20-minute call. We'll map the right course to where you want to work.

By submitting you agree to our Terms and Privacy Policy.