UNIVERSIDAD TECNOLÓGICA DE PANAMÁ

Facultad de Ingeniería Eléctrica

Licenciatura en Ingeniería Electromecánica

Laboratorio de Mecánica de Fluidos 2

Instructor Gabriel Ayú Prado

Informe de Laboratorio N°1

Determinación de las Propiedades del Aire como Gas

Ideal

BURGOS, KEILYN 6-722-1538

CASTILLO, JESÚS 4-808-2303

JUSTAVINO, ANGEL 8-954-2425

Código de Grupo (1IE-133)

Fecha (10/09/2020)
Marco teórico
Un gas ideal es un gas teórico que esta conformado por un conjunto de partículas
que tienen un desplazamiento de manera aleatoria sin interactuar la una con la otra.
Para obtener una descripción sobre el comportamiento de los gases en función de
su capacidad calorífica se encuentran distintos modelos que presentamos a
continuación:
• Gas perfecto: 𝑐𝑃 y 𝑐𝑉 son constantes.
• Gas ideal: 𝑐𝑃 y 𝑐𝑉 son dependientes de la temperatura 𝑐𝑃(𝑇) y 𝑐𝑉(𝑇).
• Gas real: 𝑐𝑃 y 𝑐𝑉 son dependientes de la temperatura y presión 𝑐𝑃(𝑇, 𝑃) y
𝑐𝑉(𝑇, 𝑃).

Comúnmente, los modelos se distinguen en función de si la presión del calor

especifico es constante o si el calor especifico es a volumen constante, pero
son dependientes de parámetros como la temperatura T o la presión P.
Para este laboratorio nos enfocaremos en el modelo de gases ideales por
como mencionamos en el punto anterior, el calor especifico de estos es
únicamente en función de la temperatura.
El 𝑐𝑃 para el aire ha sido aproximado a un polinomio de quinto grado por medio de
una serie de Taylor.

5
𝐶𝑝
= ∑ 𝑎𝑖 𝑇 𝑖
𝑅
𝑖=0
𝑎𝑖 son los coeficientes del polinomio y la temperatura 𝑇 debe estar expresada
en Kelvin. Los coeficientes 𝑎𝑖 están dados en la Tabla 1. Debido a que el
polinomio ha sido aproximado por medio de una serie de Taylor, el rango de
resultados aceptables solo es válido para un rango de 200𝐾 < 𝑇 < 2500𝐾.

Tabla 1. Coeficientes del polinomio de quinto grado usado para aproximar el

𝑐𝑃 del aire como gas ideal.

𝑎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 −3.02383 × 10−16 𝐾−5
Resultados
Dentro de la experiencia de laboratorio pudimos programar las diferentes
ecuaciones que nos brindarían estos resultados dependiendo del caso especificado.
T[𝐾] P[𝑃𝑎] 𝐽 𝐽 𝑚3
h[ ] s[ ] v[ ]
𝑘𝑔 𝑘𝑔∗𝐾 𝑘𝑔
300 101325 300426.94 5701.6312 0.8497409
459.1 103426 461469.34 6126.1257 1.2739708
269.4 99268 269735.20 5599.6112 0.7788794
522.61405 102456 526845 6262.1865 1.4639478
1382.9520 100741 1495267 7348.9718 3.9398777
625.4 186201.14 634368.77 6278.5 0.9639565
1843.5 108446.43 2057556.6 7678.5 4.8787636
345.17907 102526.12 345856 5839.3 0.9662551
648.7780 104674.15 659146 6482.7 1.7788469
1205.6762 101017.92 1284597 7185.2 3.4254228
Tabla 1. Resultados de los datos recolectados con el programa Scilab.
Análisis de Resultados
Durante la experiencia se pudo apreciar como los resultados podrían variar
dependiendo de las ecuaciones dependiendo de cual estuviésemos trabajando por
ejemplo en el caso de:
𝑇
ℎ(𝑇) = ℎ𝑟𝑒𝑓 + ∫ 𝑐𝑝(𝑇)𝑑𝑇
𝑇𝑟𝑒𝑓

podíamos ver como la entalpia dependía totalmente de la temperatura, entre más

alta la temperatura también se elevaba el valor de la entalpia.
En el caso de la entropía:

𝑇
𝑐𝑝(𝑇)𝑑𝑇 𝑃
𝑠(𝑇, 𝑃) = 𝑠𝑟𝑒𝑓 + ∫ − 𝑅𝑙𝑛( )
𝑇𝑟𝑒𝑓 𝑇 𝑃𝑟𝑒𝑓

Podemos observar en los resultados que la entropía depende totalmente de la

presión y la temperatura. En este caso también si la temperatura y la presión se
elevaban la entropía también tendía a elevarse.
Y por último el caso del volumen especifico, el cual la ecuación está dada por:
𝑃𝑣 = 𝑅𝑇

Con esta ecuación tan conocida dada cuando trabajamos con gases ideales
podemos observar que depende de tres factores, los cuales son: La Presión, La
Temperatura, Y la contante R del Gas.
Mientras que la Temperatura aumenta y la presión disminuye el volumen especifico
tiende a ser menor en valor y en el caso contrario donde la presión es mayor que la
temperatura el volumen especifico es menor, esto se puede comprobar mediante
un ejemplo sencillo donde encerramos un fluido dentro de una cámara incluyendo
un pistón para se pueda someter a más presión se puede observar como el fluido
tiende a bajar su temperatura y por ende también su volumen especifico.
Cálculos experimentales
Cálculo de Presión
𝐽 𝐽
Para T=625.4K y s=6278.5𝑘𝑔∗𝐾, 𝑠𝑟𝑒𝑓 = 6599.2 𝑘𝑔∗𝐾 , 𝑃𝑟𝑒𝑓 = 100000 𝑃𝑎, 𝑇𝑟𝑒𝑓 =
𝐽
298.15𝐾,R=287 𝑘𝑔∗𝐾

𝑇 𝑐𝑝(𝑇)𝑑𝑇 𝑃
𝑠(𝑇, 𝑃) = 𝑠𝑟𝑒𝑓 + ∫𝑇 − 𝑅𝑙𝑛(𝑃 )
𝑟𝑒𝑓 𝑇 𝑟𝑒𝑓

Dado la ecuación de la entropía despejamos para P

𝑇 𝑐𝑝(𝑇)𝑑𝑇 𝑃
𝑠 − 𝑠𝑟𝑒𝑓 − ∫𝑇 = −𝑅𝑙𝑛( )
𝑟𝑒𝑓 𝑇 𝑃𝑟𝑒𝑓

𝑇
𝑐𝑝(𝑇)𝑑𝑇 𝑃
−𝑠 + 𝑠𝑟𝑒𝑓 + ∫ = 𝑅𝑙𝑛( )
𝑇𝑟𝑒𝑓 𝑇 𝑃𝑟𝑒𝑓

Desarrollando la parte izquierda de la ecuación

𝐽 𝐽
−6278.5 + 5699.2
𝑘𝑔 ∗ 𝐾 𝑘𝑔 ∗ 𝐾
𝐽
625.4K 287 (𝑎0 + 𝑎1𝑇 + 𝑎2𝑇 2 + 𝑎3𝑇 3 + 𝑎4𝑇 4 + 𝑎5𝑇 5 )𝑑𝑇
𝑘𝑔 ∗ 𝐾
+∫
298.15𝐾 𝑇
𝑃
= 𝑅𝑙𝑛( )
𝑃𝑟𝑒𝑓
𝐽 𝐽 𝑃
178.4 = 287 𝑙𝑛( )
𝑘𝑔 ∗ 𝐾 𝑘𝑔 ∗ 𝐾 100000 𝑃𝑎

Despejando para P
𝐽
178.4 𝑃
𝑘𝑔 ∗ 𝐾
= 𝑙𝑛( )
𝐽 100000 𝑃𝑎
287
𝑘𝑔 ∗ 𝐾
𝑃
e0.6216 =
100000𝑃𝑎
(100000e0.6216 )𝑃𝑎 = 𝑃
Finalmente
𝑃 =186201.14 Pa
Cálculo de Entalpía
𝑇
ℎ(𝑇) = ℎ𝑟𝑒𝑓 + ∫ 𝐶𝑝(𝑇)𝑑𝑡
𝑇𝑟𝑒𝑓
 𝐽
ℎ = 298570 ⁄𝑘𝑔
𝐽 3.635888502𝐾(𝑇) −0.001360125𝐾(𝑇)
+ 287 ⁄𝑘𝑔 ∗ 𝐾 [( )+( )
1 2
3.81443𝑥10−6 𝐾(𝑇) −3.02383𝑥10−9 𝐾(𝑇)
+( )+( )
3 4
300
−1.03878𝑥10−12 𝐾(𝑇) −1.32878𝑥10−16 (𝑇)
+( )+( )]
5 5 298.15
 𝐽
ℎ = 300500.115 ⁄𝑘𝑔

Cálculo de volumen específico

𝑃𝑣 = 𝑅𝑇
𝑅𝑇
𝑣=
𝑃

(287 𝐽⁄𝑘𝑔 ∗ 𝐾 )(300𝐾)

𝑣=
101325𝑃𝑎
3
𝑣 = 0.849740933 𝑚 ⁄𝑘𝑔

Cálculo de entropía
𝑇
𝐶𝑝(𝑇)𝑑𝑡 𝑃
𝑠(𝑇, 𝑃) = 𝑠𝑟𝑒𝑓 + ∫ − 𝑅𝑙𝑛 ( )
𝑇𝑟𝑒𝑓 𝑇 𝑃𝑟𝑒𝑓

𝐽 101325𝑃𝑎
𝑠(𝑇, 𝑃) = 5699.2 ⁄𝑘𝑔 ∗ 𝐾 + 0.022417127 − 287𝑙𝑛 ( )
100000𝑃𝑎
𝐽
𝑠(𝑇, 𝑃) = 5695.44464 ⁄𝑘𝑔 ∗ 𝐾

Cálculo de Temperatura
𝐽 𝐽 𝐽
Para h=345856𝑘𝑔, 𝑇𝑟𝑒𝑓 = 298.15𝐾,R=287 𝑘𝑔∗𝐾, ℎ𝑟𝑒𝑓 = 298570 𝑘𝑔
𝑇
ℎ(𝑇) = ℎ𝑟𝑒𝑓 + ∫ 𝑐𝑝(𝑇)𝑑𝑇
𝑇𝑟𝑒𝑓

despejando para T
𝑇
ℎ− ℎ𝑟𝑒𝑓 = ∫ 𝑐𝑝(𝑇)𝑑𝑇
𝑇𝑟𝑒𝑓

𝐽 𝑎1𝑇 2 𝑎2𝑇 3 𝑎1𝑇 4 𝑎1𝑇 5 𝑎1𝑇 6 𝑎1𝑇𝑟𝑒𝑓 2

ℎ− ℎ𝑟𝑒𝑓 = 287 ((𝑎0𝑇 + + + + + ) − (𝑎0𝑇𝑟𝑒𝑓 +
𝑘𝑔 ∗ 𝐾 2 3 4 5 6 2

𝑎2𝑇𝑟𝑒𝑓 3 𝑎1𝑇𝑟𝑒𝑓 4 𝑎1𝑇𝑟𝑒𝑓 5 𝑎1𝑇𝑟𝑒𝑓 6

+ + + + ))
3 4 5 6
 𝐽 𝑎1𝑇 2 𝑎2𝑇 3 𝑎1𝑇 4 𝑎1𝑇 5 𝑎1𝑇 6
0 = −ℎ+ ℎ𝑟𝑒𝑓 + 287 ((𝑎0𝑇 + + + + + ) − (𝑎0𝑇𝑟𝑒𝑓
𝑘𝑔 ∗ 𝐾 2 3 4 5 6

𝑎1𝑇𝑟𝑒𝑓 2 𝑎2𝑇𝑟𝑒𝑓 3 𝑎1𝑇𝑟𝑒𝑓 4 𝑎1𝑇𝑟𝑒𝑓 5 𝑎1𝑇𝑟𝑒𝑓 6

+ + + + + ))
2 3 4 5 6

Ya que tenemos exponentes hasta 6 no podemos despejar de manera

convencional por lo cual utilizaremos el método de newton Raphson
𝑓(𝑥𝑖 )
𝑥𝑖+1 = 𝑥𝑖 −
𝑓´(𝑥𝑖 )

Comenzando con un valor de 𝑇𝑖 =300K

𝑓(𝑇𝑖 )
𝑇𝑖+1 = 𝑇𝑖 −
𝑓´(𝑇𝑖 )
𝑓(𝑇𝑖 ) = −ℎ+ ℎ𝑟𝑒𝑓
𝐽 𝑎1𝑇 2 𝑎2𝑇 3 𝑎1𝑇 4 𝑎1𝑇 5 𝑎1𝑇 6
+ 287 ((𝑎0𝑇 + + + + + ) − (𝑎0𝑇𝑟𝑒𝑓
𝑘𝑔 ∗ 𝐾 2 3 4 5 6

𝑎1𝑇𝑟𝑒𝑓 2 𝑎2𝑇𝑟𝑒𝑓 3 𝑎1𝑇𝑟𝑒𝑓 4 𝑎1𝑇𝑟𝑒𝑓 5 𝑎1𝑇𝑟𝑒𝑓 6

+ + + + + ))
2 3 4 5 6

𝐽 2 ∗ 𝑎1𝑇 1 3 ∗ 𝑎2𝑇 2 4 ∗ 𝑎1𝑇 3 5 ∗ 𝑎1𝑇 4 6 ∗ 𝑎1𝑇 5

𝑓´(𝑇𝑖 ) = 287 ((𝑎0 + + + + + ) − (𝑎0
𝑘𝑔 ∗ 𝐾 2 3 4 5 6

2 ∗ 𝑎1𝑇𝑟𝑒𝑓 1 3 ∗ 𝑎2𝑇𝑟𝑒𝑓 2 4 ∗ 𝑎1𝑇𝑟𝑒𝑓 3 5 ∗ 𝑎1𝑇𝑟𝑒𝑓 4 6 ∗ 𝑎1𝑇𝑟𝑒𝑓 5

+ + + + + ))
2 3 4 5 6

Haciendo iteraciones hasta tener un porcentaje de error de aproximadamente

0.00005
T=345.17907K
Conclusiones
- En este laboratorio se pudo determinar cómo las ecuaciones de los gases
ideales cumplen para todas sus propiedades, quedo demostrado como las
ecuaciones por ejemplo de entalpia dependen únicamente de la temperatura
y del caso de la entropía dependerá de la temperatura y de la presión, en mi
caso la preferida fue la ecuación típica de los gases ideales que una vez más
pude comprobar que el volumen especifico cumple para sus efectos de
presión y temperatura el cual con el ejemplo de la cámara de gas y el pistón
se puede comprobar como la ecuación cumple. Además, pude aprender a
programar de manera básica en el programa Scilab donde pude codificar
integrales y derivadas donde comprobé los resultados de los problemas a
mano con un simple clic.
Jesús Castillo

- En el presente informe logramos adentrarnos un poco más sobre el modelo

de gases ideales que son en función únicamente de la temperatura.
Utilizamos la ecuación del estado, la entalpia y la entropía en un gas ideal
donde las utilizamos dentro de una aplicación, en este caso utilizamos scilab
(programa que se especializa en cálculos matemáticos), donde codificamos
un programa para resolver los distintos datos solicitados en la tabla de la guía
del laboratorio. La parte más interesante del informe fue realizar el método
de newton-raphson, un método visto en materias anteriores pero difícil de
realizar en el programa para que cumpliera la función deseada.
Ángel Justavino

- En este informe podemos observar cómo un gas realmente cumple con sus
propiedades, convirtiéndolo en un gas ideal. Conocemos como gas ideal, a
un gas que se comporta de la misma manera en diferentes situaciones
porque sus propiedades y condiciones son estables. Es interesante cómo
con conocimientos de asignaturas anteriores podemos relacionarlo con
Mecánica de Fluidos, y así ahorrarnos muchísimo tiempo en cálculos,
haciendo uso de la herramienta Scilab, con conocimientos de programación
y métodos numéricos podemos crear un programa que nos dará un mínimo
margen de error en la obtención de resultados.
Keilyn Burgos
Referencias
• Mecánica de Fluidos; Yunus A. Cengel, John M. Cimbala; McGraw Hill,
Tercera Edición, 2014.
• Termodinámica; Yunus A. Cengel, Michael A. Boles; McGraw Hill, Octava
Edición, 2015.