Universidad Tecnológica de Panamá

Facultad de Ingeniería Eléctrica

Licenciatura en Ingeniería Electromecánica

Mecánica de Fluidos II

Velocidad del Sonido y Número de Mach para el Aire como Gas

Ideal
Laboratorio#2

Instructora:
Vielkis Lopéz

Integrantes:
Jorge Molina 20-53-4794
Jameson González 4-816-6
Giordi Salinas 8-917-2420
Einer Mendoza 8-964-1569
Danilo Martínez 2-747-1945

Grupo: 1IE133(A)

Fecha 14 de septiembre de 2022

1. Resumen
Este laboratorio consistió en la utilización de Scilab él cual es un software para análisis
numérico, con un lenguaje de programación de alto nivel para cálculo científico. Cuya
herramienta fue utilizada con el fin de obtener un programa con el fin de determinar diferentes
propiedades termodinámicas del aire como gas ideal al igual que en el laboratorio anterior,
aunado a eso se determinó las velocidades del sonido, del fluido y el número de Mach para
cada estado termodinámico. El programa estará escrito en un lenguaje de alto nivel basado en
el lenguaje C, en el código notará el uso de diversas ecuaciones tales como la ecuación de
estado para gas ideal, la de la entalpía en su forma integral, la entropía en su forma integral,
la relación de calores específicos y la velocidad del sonido en función de la temperatura. En
algunos casos se utilizarán métodos numéricos ya que es necesario iterar para encontrar
soluciones a ciertas ecuaciones no lineales. El método para utilizar será el método de bisección
o método Bolzano.

2. Introducción
La velocidad de propagación de las ondas sonoras es también conocida como la velocidad del
sonido, y en la atmósfera terrestre a una temperatura de 20ºC, ésta es de 343 metros por
segundo. La velocidad con que se transmite el sonido depende, principalmente, de la
elasticidad del medio, es decir, de su capacidad para recuperar su forma inicial. El acero es un
medio muy elástico, en contraste con la plastilina, que no lo es. Otros factores que influyen
son la temperatura y la densidad. En este caso nos interesa determinarla en función de la
temperatura ya que estamos tratando en un gas ideal que es el aire, y en gases ideales solo
depende de la temperatura.
El número de Mach, también representado como Ma, indica la relación entre la velocidad de
un fluido V y la del sonido C en dicho fluido. Esta sencilla relación es de vital importancia,
especialmente en la aeronáutica. La velocidad del sonido en el aire en condiciones ambientales
es de 340 m/s (unos 1238,4 km/h). Sin embargo, para que esto sea totalmente cierto, debemos
puntualizar que esto es así para el aire a una temperatura de 15˚C. Dicho valor depende de
algunas variables, como la temperatura o la presión y la densidad. Por tanto, según
aumentamos la altura respecto al suelo, menor será la velocidad del sonido. De esta forma,
con el número de Mach, sabemos a cuántas veces la velocidad del sonido se mueve un
fluido sin necesidad de conocer dicha velocidad. Por ejemplo, si un fluido viaja con un Ma=1,
sabemos que lo hace a la velocidad del sonido para esas condiciones.

3. Materiales

• Programa Scilab
4. Metodología
Se utilizó el software Scilab y la información suministrada por la profesora con el propósito
de desarrollar un programa que determine diferentes propiedades termodinámicas del aire
como gas ideal.

a) Velocidad del sonido en función de la temperatura (𝑇).

Para determinar la velocidad del sonido en función de la temperatura, se utiliza la
siguiente expresión:

𝐶 = √𝑘𝑅𝑇 (1)
Donde K es la relación de calores específicos, R es la constante específica de gas ideal
𝐽
en donde para el aire tiene un valor de 287(𝑘𝑔∙𝐾).

𝑘 = 𝐶𝑝 /𝐶𝑣 (2)

A su vez la constante de gas ideal es igual a:

𝑅 = 𝐶𝑝 − 𝐶𝑣 → 𝐶𝑣 = 𝐶𝑝 − 𝑅 (3)

Ahora remplazando la ecuación 3 en la 2 tenemos:

𝑘 = 𝐶𝑝 /𝐶𝑝 − 𝑅 (4)

Teníamos del laboratorio pasado que:

5
𝐶𝑝
= ∑ 𝑐𝑖 ∙ 𝑇 𝑖
𝑅
𝑖=0

𝐶𝑝 = 𝑅(𝑐0 + 𝑐1 ∙ 𝑇 + 𝑐2 ∙ 𝑇 2 + 𝑐3 ∙ 𝑇 3 + 𝑐4 ∙ 𝑇 4 + 𝑐5 ∙ 𝑇 5 )

Sustituyendo ambas ecuaciones en la 4:

𝑅 (𝑐0 + 𝑐1 ∙ 𝑇 + 𝑐2 ∙ 𝑇 2 + 𝑐3 ∙ 𝑇 3 + 𝑐4 ∙ 𝑇 4 + 𝑐5 ∙ 𝑇 5 )
𝑘= (5)
𝑅(𝑐0 + 𝑐1 ∙ 𝑇 + 𝑐2 ∙ 𝑇 2 + 𝑐3 ∙ 𝑇 3 + 𝑐4 ∙ 𝑇 4 + 𝑐5 ∙ 𝑇 5 ) − 𝑅

Por último, remplazamos en la primera ecuación 5 en la 1 y obtenemos la ecuación definida

para determinar la velocidad del sonido a diferentes temperaturas en el aire como medio.
 (𝑐0 + 𝑐1 ∙ 𝑇 + 𝑐2 ∙ 𝑇 2 + 𝑐3 ∙ 𝑇 3 + 𝑐4 ∙ 𝑇 4 + 𝑐5 ∙ 𝑇 5 )
𝐶 = √( ) 𝑅𝑇 (6)
(𝑐0 + 𝑐1 ∙ 𝑇 + 𝑐2 ∙ 𝑇 2 + 𝑐3 ∙ 𝑇 3 + 𝑐4 ∙ 𝑇 4 + 𝑐5 ∙ 𝑇 5 ) − 1

Hay que recordar que los valores que pueden tomar los coeficientes del polinomio de
Taylor son los siguientes:

Tabla#1 valores de cada coeficiente

𝑐0 3.635888502
𝑐1 −0.001360125 𝐾 −1
𝑐2 3.81443 × 10−6 𝐾 −2
𝑐3 −3.02383 × 10−9 𝐾 −3
𝑐4 1.03878 × 10−12 𝐾 −4
𝑐5 −1.32878 × 10−16 𝐾 −5

b) Temperatura como función de la velocidad del sonido 𝑇(c).

Para determinar la temperatura en función de la velocidad del sonido se necesita hacer uso
de un método numéricos que nos permite determinar las raíces de la ecuación 6, que es la
temperatura que necesitamos determinar. Se utilizó el método de bisección para realizar
este cálculo.
En este método se eligen 2 valores de temperatura que dentro de ese rango se encuentre
una raíz, dichos valores son 200 y 2500 grados kelvin. El método consiste en calcular el
valor medio de las temperaturas mencionadas y las 3 evaluarlas en la función igualada a la
velocidad del sonido. En las funciones evaluadas, donde haya un cambio de signo se
encuentra una raíz y el valor independiente que resultó en la función con el cambio de signo
con respecto a la evaluada en el valor medio, se mantiene constante para la siguiente
iteración, y el otro valor se tomará del valor medio de la iteración anterior.

c) Número de Mach como función de la temperatura y velocidad local (T, V).

El número de Mach es la relación de la velocidad local del fluido y la velocidad del sonido,
la cual está en función de la temperatura.
𝑀𝑎 = 𝑉/𝐶(𝑇) (7)
Para determinar la velocidad del sonido en función de la temperatura se utiliza la
ecuación 6.

d) Velocidad como función de la temperatura y número de Mach V (𝑇, Ma).

Simplemente se despeja la ecuación 7 para la velocidad del fluido y se realiza el mismo
procedimiento que en el caso anterior para determinar la velocidad del sonido en función
de la temperatura.
𝑉 = 𝐶(𝑇) ∙ 𝑀𝑎
e) Determinación de otras propiedades.
Se utilizaron las mismas ecuaciones y métodos del laboratorio anterior para determinar ya
sea temperatura, entalpía, presión, entropía, en donde estén en función de alguna otra
propiedad de estas mismas mencionadas.

Programa en Scilab para determinar las propiedades, velocidad del aire y del
sonido en diferentes estados como gas ideal:

function lab2()
R=287 //J/Kg-K
a0=3.635888502 //1/K
a1=-0.001360125
a2=0.00000381443
a3=-0.00000000302383
a4=0.00000000000103878
a5=-0.000000000000000132878
Pref=100000 //Pa
href=298570 //a T=298.15K en J/Kg
tref=298.15 //K
sref=5699.2 // aT=298.15K y P=100000 en J/Kg-K
printf('El siguiente programa calcula las propiedades del aire como gas ideal ')
printf('\nconocidas dos propiedades intensivas e independientes, mas una propiedad extra ')
printf('\ny calcula las propiedades desconocidas')
printf('\n1. Si conoce la Temperatura (K), la presion (Pa) y la velocidad promedio (m/s)')
printf('\n2. Si conoce la temperatura (K), la presion (Pa) y el numero mach')
printf('\n3. Si conoce la presion (Pa), la velocidad promedio (m/s) y la velocidad del sonido (m/s)')
printf('\n4. Si conoce la Presion (Pa), la Velocidad promedio (m/s) y si conoce la entalpia (J/kg)')
printf('\n5. Si conoce la presion (Pa), el numero mach y si conoce la entalpia (J/kg)')
printf('\n6. Si conoce la velocidad promedio(m/s), si conoce la velocdiad del sonido (m/s) y si conoce la
entropia (J/kg*k)')
printf('\n7. Si conoce la velocidad promedio en (m/s), si conoce la entalpia (J/kg) y si conoce la entropia
(J/kg*K)')
printf('\n8. Si conoce Mach, si conoce la entalpia (J/kg) y si conoce la entropia (J/kg*K)')

/*Si conoce la entalpia (J/Kg) y la entropia (J/Kg*k)

Si conoce la presion (Pa), la velocidad promedio (m/s) y la velocidad del sonido (m/s)*/

printf('\n')
printf('\n')
x=input (' Ingrese el caso requerido: ')
if x==1 then
t=input('Ingrese la temperatura en K: ')
p=input('Ingrese la presion en Pa: ')
vel=input('Ingrese la velocidad promedio en m/s: ')

deff ('cp=f(T)','cp=R*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
hb=href+intg(tref,t,f)
deff ('cp=g(T)','cp=(R/T)*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
s=sref+((intg(tref,t,g)))-(R*log(p/Pref))
v=(R*t)/p
cp=R*(a0+(a1*t)+(a2*t^2)+(a3*t^3)+(a4*t^4)+(a5*t ^5))
k=(cp)/(cp-R)
c=sqrt(k*R*t)
ma=vel/c
printf('La h es= %f\n',hb)
printf('La s es= %f\n',s)
printf('La v es= %f\n',v)
 printf('c es= %f\n',c)
printf('Ma es= %f\n',ma)
end
if x==2 then
p=input('Ingrese el valor de la presion en Pa: ')
t=input('Ingrese el valor de la temperatura en K ')
ma=input('Ingrese el valor de Mach: ')
deff ('cp=f(T)','cp=R*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
hb=href+intg(tref,t,f)
deff ('cp=g(T)','cp=(R/T)*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
s=sref+((intg(tref,t,g)))-(R*log(p/Pref))
cp=R*(a0+(a1*t)+(a2*t^2)+(a3*t^3)+(a4*t^4)+(a5*t ^5))
k=cp/(cp-R)
c=sqrt(k*R*t)
vel=ma*c
printf('La h es= %f\n',hb)
printf('La s es= %f\n',s)
printf('La velocidad promedio es= %f\n',vel)
printf('La velocidad del sonido es c= %f\n',c)

end
if x==3 then
p=input('Ingrese la presion en Pa: ')
vel=input('Ingrese la velocidad promedio en (m/s): ')
c=input('Ingrese la velocidad del sonido (m/s): ')
a=200
b=2500
iter=1
imax=input('Ingrese el numero de iteraciones que desea realizar para T:')
while iter<=imax
deff ('cp=f(T)','cp=R*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
T=(a+b)/2
cp=R*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))
ca=c
cb=sqrt((cp/(cp-R))*R*T)
z=cb-ca
iter=iter+1
if z>0 then
b=T
else
a=T
end
end

Ma=vel/c
deff('cp=e(T)','cp=R*((a0)+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
hb=href+intg(tref,T,e)
deff('cp=g(T)','cp=(R/T)*((a0)+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
s=sref+intg(tref,T,g)-(R*log(p/Pref))
printf('La Temperatura= %f\n',T)
printf('El numero mach es= %f\n',Ma)
printf('La s es= %f\n',s)
printf('La h es= %f\n',hb)

end
if x==4 then
p=input('Ingrese el valor de la presion en Pa: ')
vel=input('Ingrese el valor de la velocidad promedio en (m/s): ')
h=input('Ingrese el valor de la entalpia en J/Kg: ')
a=200
b=2500
iter=1
imax=input('Ingrese el numero de iteraciones que desea realizar para T:')
while iter<=imax
deff ('cp=f(T)','cp=R*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
T=(a+b)/2
ha=h
hb=href+intg(tref,T,f)
z=hb-ha
iter=iter+1
if z>0 then
b=T
else
a=T
end
end

printf('La T es = %f\n',T)
cp=R*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))
k=cp/(cp-R)
c=sqrt(k*R*T)
Mach=vel/c
deff('cp=g(T)','cp=(R/T)*((a0)+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
s=sref+intg(tref,T,g)-R*log(p/Pref)
v=(R*T)/p
printf('El volumen es= %f\n',v)
printf('La velocidad del sonido es= %f\n',c)
printf('El numero Mach es= %f\n',Mach)
printf('La s es= %f\n',s)
end
if x==5 then
p=input('Ingrese la presión en Pa= ')
Ma=input('Ingresse el numero Mach= ')
h=input('Ingrese la entalpia en (J/kg)= ')
a=200
b=2500
iter=1
imax=input('Ingrese el numero de iteraciones que desea realizar para T:')
while iter<=imax
deff ('cp=f(T)','cp=R*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
T=(a+b)/2
ha=h
hb=href+intg(tref,T,f)
z=hb-ha
iter=iter+1
if z>0 then
b=T
else
a=T
end
end
cp=R*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))
k=cp/(cp-R)
c=sqrt(k*R*T)
vel=Ma*c
deff('cp=g(T)','cp=(R/T)*((a0)+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
s=sref+intg(tref,T,g)-R*log(p/Pref)
printf('La temperatura en K es= %f\n',T)
printf('La Velocidad promedio es (m/s)= %f\n',vel)
printf('La velocidad del sonido es c (m/s)= %f\n',c)
printf('La entropia s (J/kg*K) es= %f\n',s)
end
 if x==6 then
vel=input('Ingrese la velocidad promedio (m/s)= ')
c=input('Ingresse la velocidad del sonido (m/s)= ')
s=input('Ingrese la entropia en (J/kg*K)= ')
a=200
b=2500
iter=1
imax=input('Ingrese el numero de iteraciones que desea realizar para T:')
while iter<=imax
T=(a+b)/2
ca=c
cp=R*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))
k=cp/(cp-R)
cb=sqrt(k*R*T)
z=cb-ca
iter=iter+1
if z>0 then
b=T
else
a=T
end
end
Ma=vel/c
deff ('cp=g(T)','cp=(R/T)*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
p=(%e^((sref+(intg(tref,T,g))-s)/R))*Pref
deff ('cp=h(T)','cp=(R)*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
h=href+((intg(tref,T,h)))
printf('La temperatura en K es= %f\n',T)
printf('La presion en Pa es= %f\n',p)
printf('El numero mach es= %f\n',Ma)
printf('La entalpia en J/kg es= %f\n',h)
end
if x==7 then
v=input('Ingrese la velocidad en m/s= ')
h=input('Ingrese la entalpia en J/kg= ')
s=input('Ingrese la entropia en J/kg*K= ')
a=200
b=2500
iter=1
imax=input('Ingrese el numero de iteraciones que desea realizar para T:')
while iter<=imax
deff ('cp=f(T)','cp=R*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
T=(a+b)/2
ha=h
hb=href+intg(tref,T,f)
z=hb-ha
iter=iter+1
if z>0 then
b=T
else
a=T
end
end
deff ('cp=g(T)','cp=(R/T)*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
p=(%e^((sref+(intg(tref,T,g))-s)/R))*Pref
cp=R*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))
k=cp/(cp-R)
c=sqrt(k*R*T)
Ma=v/c
printf('La temperatura en K es= %f\n',T)
printf('La presion en Pa es= %f\n',p)
 printf('La velocidad del sonido en m/s es= %f\n',c)
printf('El numero Mach es= %f\n',Ma)

end
if x==8 then
Ma=input('Ingrese el numero Mach= ')
h=input('Ingrese la entalpia en J/kg= ')
s=input('Ingrese la entropia en J/kg*K= ')
a=200
b=2500
iter=1
imax=input('Ingrese el numero de iteraciones que desea realizar para T:')
while iter<=imax
deff ('cp=f(T)','cp=R*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
T=(a+b)/2
ha=h
hb=href+intg(tref,T,f)
z=hb-ha
iter=iter+1
if z>0 then
b=T
else
a=T
end
end
deff ('cp=g(T)','cp=(R/T)*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))')
p=(%e^((sref+(intg(tref,T,g))-s)/R))*Pref
cp=R*(a0+(a1*T)+(a2*T^2)+(a3*T^3)+(a4*T^4)+(a5*T^5))
K=cp/(cp-R)
c=sqrt(K*R*T)
v=Ma*c
printf('La temperatura en K= %f\n',T)
printf('La presion en Pa= %f\n',p)
printf('La velocidad en m/s= %f\n',v)
printf('La velocidad del sonido en m/s= %f\n',c)

end

endfunction
 Consola del programa

5. Análisis y Resultados
Tabla#2 propiedades termodinámicas del aire como gas ideal en diferentes estados.

𝑻(𝑲) 𝑷(𝑷𝒂) V (𝒎/𝒔) C (𝒎/𝒔) Ma 𝑱 𝑱

𝒉( ) 𝒔( )
𝒌𝒈 𝒌𝒈 ∙ 𝑲

456.5 105512 213.6 426.731204 0.500549 458808.772973 6114.583107

1534.5 156183 121.493935 759.337094 0.16 1678283.179130 7348.691448

851.421653 321451 98.6 574.3 0.171717 879149.829300 6455.505801

725.457357 105316 321.6 532.5 0.603944 741287.404127 6600.596451

1672.167078 99513 361.6 791.532422 0.456835 1846365 7582.944066

1014.540532 113175 386.598891 623.546599 0.62 1062851 6952.421141

726.323543 245233.414646 453.3 532.8 0.850788 742223.002001 6359.3

1542.734607 127914.806685 943.6 761.3 1.239459 1688291.870788 7412.5

2087.382130 121175.425232 895.6 881.380727 1.016133 2361558 7801.5

315.447138 96150.913550 523.301414 355.987356 1.47 315941 5767.1

6. Conclusiones
Durante el proceso de la realización del experimento, hemos tenido la oportunidad de conocer
otras alternativas para conocer las diferentes propiedades del aire, con la ayuda del programa
Scilab, que es un programa de lenguaje c destinado a la programación en general. Hemos
utilizado para obtener ecuación que nos ayudara a encontrar diferentes propiedades, la
velocidad del sonido, del flujo y su relación todo a través del aire como medio. Para esto se
utilizaron las ecuaciones de la velocidad del sonido en función de la temperatura y la de calores
específicos ya que como ambos dependen de la temperatura, podemos determinar la propiedad
o propiedades desconocidas a través de ambas, con la ayuda de ecuaciones del laboratorio
anterior como la determinación de la temperatura en función de la entalpía, esto con la ayuda
del método de bisección. Al final del laboratorio se utilizó la ecuación general del gas ideal
ya que en este caso asumimos que el aire es un gas ideal y tuvimos que manipular la ecuación
para obtener otras propiedades de la ecuación.
Como resultado del laboratorio se logró visualizar como la velocidad del sonido aumenta
cuando la temperatura es mayor, y a su vez la velocidad local disminuye, esto es importante
ya que en ocasiones donde la temperatura era bastante grande, el Número de Mach tendía a
ser menor que uno, siempre y cuando la presión no fuera tan pequeña, si esta decrece, el fluido
tiende a ser transónico.

7. Referencias bibliográficas

• Cengel, Y. (2019). Termodinámica . México D. F.: McGRAW-HILL.

• Çengel, Y. A., & Cimbala, J. M. (2018). Mecánica de fluidos: Fundamentos Y

aplicaciones.
• Chapra, S. C., & Canale, R. P. (2015). Métodos

• ¿Cuál es la velocidad del sonido en el Aire? (2021, December 23). Página

principal.

• ¿Qué es el número de Mach? (2021, June 20). Aerodinámica F1.