Velocidades verdaderas, equivalentes y calibradas

  • La velocidad verdadera (TAS) True Air Speed es la velocidad real a la que el avión se desplaza con respecto al aire.
  • Por el hecho de volar más alto, la densidad es menor y para una misma velocidad, la diferencia entre presión estática y dinámica es menor que volando más bajo. La velocidad equivalente (EAS) Equivalent Air Speed indica qué velocidad provoca ese salto de presiones a nivel del mar.
  • La velocidad calibrada (CAS) Calibrated Air Speed trata de corregir, la EAS. Es la que un instrumento de vuelo debería mostrar si no tuviese errores.
  • La velocidad indicada (IAS) Indicated Air Speed es la que finalmente muestra el intrumento.

Realción entre EAS y TAS.

$$\frac{1}{2}\rho V_{TAS}^{2}=\frac{1}{2}\rho V_{EAS}^{2}\rightarrow V_{EAS} = \sqrt{\sigma}\cdot V_{TAS}$$

Velocidad CAS.

Se parte de la ecuación de Saint-Venant:

$$p_{T}=p_{S}\left(1+\frac{\gamma-1}{2}M^{2}\right)^{\frac{\gamma}{\gamma-1}}$$

Sustituyendo el número de Mach por $M = V_{TAS} \cdot \sqrt{\gamma R_g T} = V_{TAS} \cdot \sqrt{\gamma p / \rho}$ y despejando $V_{TAS}$ se obtiene:

$$V_{TAS}^{2}=\frac{2\gamma}{\gamma-1}\frac{p_{s}}{\rho}\left[\left(\frac{\bigtriangleup p}{p_{s}}+1\right)^{\frac{\gamma-1}{\gamma}}-1\right]$$

La CAS es la misma que la TAS sería a nivel del mar:

$$V_{CAS}^{2}=\frac{2\gamma}{\gamma-1}\frac{p_{0}}{\rho_0}\left[\left(\frac{\bigtriangleup p}{p_{0}}+1\right)^{\frac{\gamma -1}{\gamma}}-1\right]$$

In [1]:
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline
In [2]:
#Parámetros problema

gamma = 1.4

p_0 = 101325 #Pa
T_0 = 288.15 #K
rho_0 = 1.225 #Kg/m^3
a_0 = np.sqrt(gamma * p_0 / rho_0)
In [3]:
#ISA
# h siempre en metros h<11000 m

def p_ISA(h):
    return p_0 * (1 - 22.558e-6 * h) ** 5.2559

def T_ISA(h):
    return T_0 * (1 - 22.558e-6 * h)

def rho_ISA(h):
    return rho_0 * (1 - 22.558e-6 * h) ** 4.2559

def a_ISA(h):
    return a_0 * (1 - 22.558e-6 * h) ** 0.5
In [17]:
def dif_presion(V_TAS, h):
    ''' Devuelve la diferencia de presión en Pa que
    mediría el instrumento si fuese perfecto'''
    rho = rho_ISA(h)
    p = p_ISA(h)
    V = V_TAS
    
    y = (gamma - 1) / (2 * gamma) * rho / p * V ** 2 + 1
    y = (y ** ( gamma / (gamma -1)) - 1) * p
    return y

def V_TAS(dif_presion, h):
    p = p_ISA(h)
    rho = rho_ISA(h)
    dp = dif_presion
    
    y =  2 * gamma / (gamma - 1.) * p / rho
    y = y *( (dp/p + 1) ** ((gamma - 1) / gamma) - 1)
    return np.sqrt(y)
    
def V_CAS(h, Vel_TAS):
    dp = dif_presion(Vel_TAS, h)
    return V_TAS(dp, 0)

def V_EAS(h, V_TAS):
    return rho_ISA(h) / rho_0 * V_TAS

Diferencia de presión que genere una velocidad según la altura

In [13]:
h = np.linspace(0, 11000, 1000)
V = np.arange(50, 350, 50)

fig, ax = plt.subplots()

for ii in V:
    ax.plot(h, dif_presion(ii, h)/p_0, label=u'$V_{TAS} = %i m/s$' %ii )

ax.set_ylabel(u'$\Delta p /p_0$')
ax.set_xlabel(u'$h(m)$')
ax.legend()
ax.grid()
In [16]:
h = np.linspace(0, 11000, 1000)
V = np.linspace(0, 340, 1000)

VV, hh = np.meshgrid(V,h)

plt.contourf(VV, hh, dif_presion(VV,hh)/p_0,200)
plt.colorbar()
plt.contour(VV, hh, dif_presion(VV,hh)/p_0, 10, colors = 'black')

plt.title(u'$\Delta p /p_0$')
plt.ylabel(u'$h(m)$')
plt.xlabel(u'$V_{TAS}(m/s)$')
Out[16]:
<matplotlib.text.Text at 0x131bc400>

Diferencias entre los valores de las velocidades

In [19]:
V = 150
h = np.linspace(0, 11000, 1000)

plt.plot(h, V*np.ones_like(h))
plt.plot(h, V_EAS(h, V))
plt.plot(h, V_CAS(h, V))
Out[19]:
[<matplotlib.lines.Line2D at 0xa375898>]
In [49]:
V = np.arange(50, 400, 50)
h = np.linspace(0, 11000, 1000)

fig, ax = plt.subplots(nrows= np.size(V))
fig.set_size_inches(10,15)

counter = 0

for ii in V:
    axn = ax[counter]
    
    axn.plot(h, ii*np.ones_like(h), label = u'$V_{TAS}$', lw = 2)
    axn.plot(h, V_EAS(h, ii), label = u'$V_{EAS}$')
    axn.plot(h, V_CAS(h, ii), label = u'$V_{CAS}$')
    
    axn.set_ylim(0,ii*1.10)
    axn.set_xlim(0, 11000)
    axn.set_xlabel(u'$h(m)$')
    axn.set_ylabel(u'$V(m/s)$')
    
    axn.legend(loc=3)
    axn.grid()
    
    counter += 1