## Aerodynamics Example #1: Finding the Critical Mach Number of an Airfoil |

**Some Background**

In the field of aerodynamics, airspeeds are usually stated in terms of *Mach number*, which is the local airspeed divided by the local speed of sound.

For example, assume we are considering airflow at a point, say, "A" on the surface of an airfoil.

Assume we know the local airspeed is **v _{A}**.

The local speed of sound, **a**, depends upon the local temperature according to the relation

**a** = √(** γ R T **)

**γ** is defined as

**γ ≡ C _{p}/C_{v}**

where **C _{p}** and

**C**are the specific heats of a gas at constant pressure and volume, respectively. For most calculations, standard air conditions are assumed, and a value of

_{v}**γ**= 1.4 is used.

**R** is the specific gas constant. For air at standard conditions, **R** = 287 J/(kg K).

**T** is the temperature, in degrees Kelvin.

Assume we know the speed of sound at "A" is **a _{A}**.

The Mach number at "A" would be calculated as

**M _{A}** =

**v**/

_{A}**a**.

_{A}Further analysis would continue using **M _{A}** rather than

**v**or

_{A}**a**.

_{A}
**Why is the critical Mach number important?**

For values of **M _{∞}** slightly greater than

**M**, the airfoil experiences a dramatic increase in drag coefficient. As a result, for more detailed analysis, it is important to know when this condition occurs.

_{cr}
**On With the Present Example**

Consider an airfoil in sub-sonic airflow. The undisturbed airflow far ahead of the airfoil is called the *freestream* airflow, whose airspeed is usually denoted **M _{∞}**.

As the airflow passes over the airfoil, its speed increases and, somewhere on the surface of the airfoil, its speed reaches a maximum value before completely passing over the airfoil and returning to its freestream airspeed. Now increase the the speed of the freestream airflow until the maximum airspeed achieved over the airfoil is M = 1. This condition defines the *critical Mach number*. By definition, the critical Mach number, **M _{cr}**, is the

**freestream**Mach number at which sonic flow is first achieved on the airfoil.

**Calculating the critical Mach number**

The critical Mach number is determined by equating two pressure coefficients: (**C _{p}**)

_{min}and

**C**.

_{p,cr}
(**C _{p}**)

_{min}is given by

(**C _{p}**)

_{min}= (

**C**)

_{p,0}_{min}/√(1 - (

**M**)

_{∞}^{2})

where (**C _{p,0}**)

_{min}is the minimum pressure coefficient for the airfoil at a 0 degree angle of attack. It is usually supplied for a particular airfoil, determined experimentally (in a wind tunnel).

In turn, **C _{p,cr}** is calculated from an equation:

**C _{p,cr}** = 2/(

**γ M**

_{cr}^{2}) {[(1 + ((γ - 1)/2)

**M**

_{cr}^{2})/(1 + (

**γ**- 1)/2)]

^{γ/(γ - 1)}- 1}

The critical Mach number is then calculated by setting

(**C _{p}**)

_{min}=

**C**

_{p,cr}which can be rearranged as

(**C _{p}**)

_{min}-

**C**= 0

_{p,cr}Expanding the terms yields:

(**C _{p,0}**)

_{min}/√(1 - (

**M**)

_{∞}^{2}) - 2/(

**γ M**

_{∞}^{2}) {[(1 + ((γ - 1)/2)

**M**

_{∞}^{2})/(1 + (

**γ**- 1)/2)]

^{γ/(γ - 1)}- 1} = 0

Note that all instances of **M _{cr}** were replaced by

**M**; for the desired condition, they are the same.

_{∞}
**A Numerical Example**

Assume you are trying to determine the critical Mach number of the NACA 0012 airfoil.

For the given conditions, you are given that (**C _{p,0}**)

_{min}= -0.43.

You plan to use the equation just above to determine the critical Mach number:

This equation is solved by using Root-Finder Utility #2; however, the data must be altered a little to get it in the correct form for use with the utility.

For Root-Finder Utility #2,

constant **a _{1}** is entered as

**a**= (

_{1}**C**)

_{p,0}_{min}(γ/2) = (-0.43)(1.4/2) = -0.301.

constant

**a**is entered as

_{2}**a**= (γ - 1)/2 = (1.4 - 1)/2 = 0.2.

_{2}constant

**a**is entered as

_{3}**a**= γ/(γ - 1) = 1.4/(1.4 - 1) = 3.5.

_{3}The endpoints

**b**and

**c**can be suggested by graphing the above equation, or simply trial-and-error. The upper endpoint can be set to 1, because we know

**M**MUST BE LESS THAN 1. Let us set

_{∞}**b**= 0.5 and

**c**= 1 for the present example.

Plugging these numbers into Root-Finder Utility #2 yields M

_{cr}= 0.7371059