density-and-air-temperature_COVER

Air density \rho is a typical case of measurement derived from air data. Let’s investigate the relationship between air temperature and air density, which, incidentally, is the base of the external air temperature measurement in moving vehicles.

Commonly the air density can be calculated with an equation of state or with an interpolated best-fit curve. In both cases, we have to deal with an equation of the form \rho=f(T, p, RH, others). Let’s make things a bit simpler and start with a simple relationship between density and temperature:

    \[\rho=f(T, Constant)\]

The temperature in the SI is measured in degrees Kelvin, but many equations of state may also be given with temperature expressed in Celsius=Kelvin - 273.15.

Quoted from www.bipm.org: “The kelvin, unit of thermodynamic temperature, is the fraction 1/273.16 of the thermodynamic temperature of the triple point of water.”.

The thermodynamic temperature is correlated to the state of a thermodynamic system, and according to the third law of thermodynamics, it defines an absolute temperature scale.

For dry air, the ideal gas equation of state is: \rho=\frac{\displaystyle p}{\displaystyle R_{specific} \cdot T} where:

  • \rho is the air density in kg/m^3
  • p is the absolute pressure in Pa
  • T is the absolute temperature in K
  • R_{specific}=287.058\ (J/(kg \cdot K)) is the specific gas constant for air.

The main advantage of the equation of state is that it provides a closed form solution, which in turn can be incorporated in other computations. This is a very convenient property for analysis purposes. Let’s proceed with a density calculation. Refer to figure 1; our idealized test setup is a room with two thermometers and a fan. The first thermometer measures the temperature of the air in the room at a point far away from the fan effects, while the second one is subject to the direct airflow from the fan.

density-and-air-temperature_FIGURE_1

Figure 1 – Example layout

Before any calculation let’s introduce the stagnation temperature, which indicates the temperature at a stagnation point within a flow field. At the stagnation point, the local flow velocity is zero. We find a stagnation point on a Pitot tube tip or in a total air temperature probe.

In our example, the air from the fan which hits the thermometer is brought to rest on the surface of the sensing element. Consider an infinitesimal volume of the overall air stream traveling at a certain speed.

This small volume has a corresponding kinetic and gravitational potential energy. If we assume that this volume does not exchange any form of energy with the surroundings then this energy must be converted when one of its forms is eliminated. In the process where the air slows down, kinetic energy is converted into heat, which causes a temperature rise.

Assuming an adiabatic process we can describe the enthalpy at stagnation:

    \[h_s=h+\frac{V^2}{2}\]

where:

  • the enthalpy term h is that of the nearest point along the streamline
  • V is the speed of the air flow.

For an ideal gas it holds that h=u+pv where:

  • v is the specific volume
  • u is the internal energy of the gas.

By manipulating the enthalpy expression and taking into account the constant pressure heat capacity definition Cp for an ideal gas h=C_p \cdot T, we arrive at:

    \[T_s=T+\frac{V^2}{2C_p}\]

This equation defines how much the temperature rises because is adiabatically brought to rest.

Returning to our example let’s calculate the different values of stagnation temperature in correspondence of different airspeed values. You can find pre-calculated values in the table below.

T [K] Static V [m/s] V [km/h] M T0 [K] Stagnation Deviation [K] Deviation [%] Parameters
288.15
0
0
0
288.15
0
0
Cp
1005
J/Kg/K
288.15
2
7.2
0.01
288.15
0
0
c
340.3
m/s
288.15
3
10.8
0.01
288.15
0
0
288.15
4
14.4
0.01
288.16
-0.01
0
288.15
5
18
0.01
288.16
-0.01
0
288.15
6
21.6
0.02
288.17
-0.02
0.01
288.15
7
25.2
0.02
288.17
-0.02
0.01
Table 1 – Stagnation temperature vs static temperature at different V values. Download spreadsheet

By inspection of the table, the deviation between the two temperature measurements increases with the air speed.

The air density in the room should be calculated using the static temperature, in our case, e.g. p=101325\ Pa the corresponding density is \rho=1.225\ kg/m^3. The same calculation using the total temperature with a fan pushing the air at 55 m/s gives \rho=1.218\ kg/m^3, which is about a 0.5% of error on density.

Inspect the relationship between the static and total temperature for a total air temperature probe; the first formula in this link.

    \[T_{total}/T_s=1+\frac{\gamma-1}{2}M^2\]

Where the “total” subscript indicates the adiabatic temperature and “s” subscript is for the static temperature. To reach this formula from our derived temperature equation, substitute in this formula the definition of the specific heat ratio \gamma=C_p/C_v, Mayer’s relation C_p - C_v = R and the Mach number expression:

    \[M=V/\sqrt{\gamma\cdot R_{specific}\cdot T}\]

The total air temperature computed up to now here is somewhat ideal. All the kinetic energy is converted and contributes to a temperature rise, but it should be accounted with an appropriate recovery factor. The relationship of the recovery factor with airspeed and other environmental factors is not trivial, and for best accuracy, it should be experimentally found. Nonetheless, the tight relationship between airspeed and temperature can be exploited to have a true airspeed indication as per equation 9 in the previous link. Consequently, a total air temperature probe working together with a Pitot-static tube offers a useful redundant airspeed indication.

Even at low airspeeds, to achieve top performance, the temperature measurement rise due to airspeed is existent. However, the impact of this effect on temperature measurements can be neglected in many applications with moderate to no loss of accuracy. Unfortunately, this is not the only source of uncertainty in a total air temperature probe.

Further reading

0 Comments

Leave a Reply

Avatar placeholder