PRQ 2202

FLUJO DE FLUIDOS COMPRESIBLES

Número de Mach
En los problemas referentes al flujo de gases es conveniente utilizar un parámetro de
velocidades sin dimensiones denominado número de Mach que se define como la relación entre
la velocidad de flujo y la velocidad sónica a las condiciones locales. Así, para un gas ideal:
𝑣 𝑣
𝑀𝑎 = = (1)
𝑣𝑠 √𝛾𝑃𝑉̂

La velocidad acústica es la velocidad de propagación de una pequeña perturbación bajo

condiciones reversibles ideales, ésta será la velocidad máxima con la que se propaga la
perturbación molecular en un fluido. Sin embargo, grandes perturbaciones de violencia
explosiva, con cambios importantes en presión y densidad pueden viajar a velocidades varias
veces la del sonido. Tales perturbaciones se denominan ondas de choque.

Escurrimiento Adiabático
Cañería aislada térmicamente: Suposiciones

 𝑄̂ = 0, esto indica que no hay transferencia de calor.

 𝑊̂ = 0, esto quiere decir que la cañería no tiene, entre los límites elegidos, maquinas
compresoras, ventiladores, etc.
 𝜌𝑔∆𝑧 = 0, el peso de la columna de gas, aunque la cañería sea vertical, es pequeña.

Así, con las consideraciones anteriores, aplicando la ecuación de energía total a los límites del
sistema.

1ra Ley, sistema abierto:

1
̂ + 𝑃𝑉̂ + 𝑣 2 + Φ) = 𝑄̂ − 𝑊
∆ (𝑈 ̂ (2)
2

̂ 𝑑𝑉
𝑈𝑡𝑜𝑡 = ∫ 𝜌𝑈
𝑉

Φtot = ∫ ρΦdV
V

1 2
𝑘𝑡𝑜𝑡 = ∫ 𝜌𝑣 𝑑𝑉̂
2
𝑉
̂ + 𝑃𝑉̂ = 𝐻
𝑈

𝑑
(𝑈 + 𝑘𝑡𝑜𝑡 + 𝜙𝑡𝑜𝑡 ) = 𝜌1 𝑣1 𝑉̂1 𝑆1 − 𝜌2 𝑣2 𝑉̂2 𝑆2 + ⋯
𝑑𝑡 𝑡𝑜𝑡

1 1
… + 𝜌1 (𝑣1 )3 𝑆1 𝜌2 (𝑣2 )3 𝑆2 + 𝜌1 𝑣1 𝜙1 𝑆1 − 𝜌2 𝑣2 𝜙2 𝑆2 + 𝑄̂ − 𝑊
̂ + 𝑃1 𝑉1 𝑆1 − 𝑃2 𝑉2 𝑆2
2 2

1
 PRQ 2202

𝐸𝑡𝑜𝑡 + 𝑈𝑡𝑜𝑡 + 𝑘𝑡𝑜𝑡 + 𝜙𝑡𝑜𝑡 , con velocidades masicas: 𝑤1 = 𝜌1 𝑣1 𝑆1 𝑦 𝑤2 = 𝜌2 𝑣2 𝑆2 , entonces:

𝑑(𝐸𝑡𝑜𝑡 ) 1
̂ + 𝑃𝑉̂ + 𝑣 2 + Φ) + 𝑄̂ − 𝑊
= −∆ [(𝑈 ̂]
𝑑𝑡 2
0 0 0
1 ∆𝑣 2
̂ + 𝑃𝑉̂ + 𝑣 2 + Φ) = 𝑄̂ − 𝑊
∆ (𝑈 ̂ => ̂ + 𝑃𝑉̂) = −
∆(𝑈 (3)
2 2
2 2
̂ = − ∆𝑣 , ó también 𝑑𝑈
∆𝐻 ̂ + 𝑑(𝑃𝑉̂) = − 𝑑𝑣
2 2

Ecuación de Energía:
̂ + 𝑑(𝑃𝑉̂) + 𝑣𝑑𝑣 = 0
𝑑𝑈 (4)

Ecuación de Energía mecánica: (escrita de manera diferencial)

2𝑓𝑣 2 𝑑𝐿
𝑣𝑑𝑣 + 𝑉̂𝑑𝑃 + =0 (5)
𝐷
𝑅
̂ = 𝑐̂𝑣 𝑑𝑇, para los gases ideales => 𝑑𝑈
𝑑𝑈 ̂= 𝑑𝑇
𝑀(𝛾−1)

𝑅 𝑅 1
̂=
𝑑𝑈 𝑑(𝑃𝑉̂) = 𝑑(𝑃𝑉̂) (6)
𝑀(𝛾 − 1) 𝑀 (𝛾 − 1)

Reemplazando (6) en (4):

1
𝑑(𝑃𝑉̂) + 𝑑(𝑃𝑉̂) + 𝑣𝑑𝑣 = 0 (7)
(𝛾 − 1)

1 𝑣2
𝑑[ (𝑃𝑉̂) + ] = 0 (8)
(𝛾 − 1) 2

Aplicando la ecuación (8) a 2 puntos cualesquiera en la cañería, pero con:

𝐺
𝑣=𝐺= => 𝑑𝑣 = 𝐺𝑑𝑉̂
𝜌
𝐺 = 𝑑𝑒𝑛𝑠𝑖𝑑𝑎𝑑 𝑑𝑒 𝑓𝑙𝑢𝑗𝑜 𝑚𝑎𝑠𝑖𝑐𝑜 = [𝑚𝑎𝑠𝑎/𝑎𝑟𝑒𝑎 ∗ 𝑡𝑖𝑒𝑚𝑝𝑜]

Integrando la ecuación (8):

(𝛾 − 1) 2 2
[𝑃𝑉̂ + 𝐺 𝑉 ] = 𝑐𝑡𝑡𝑒 = 𝐶 (9)
2𝛾

La ecuación (9) puede reajustarse para expresar las relaciones entre la P, la T° y las relaciones
de volumen, dividiendo todo por 𝜌1 𝑉̂1, aplicando a los límites y despejando de manera que nos
convenga:

2
 PRQ 2202
2
𝑃2 𝑉̂1 𝛾 − 1 𝐺 2 𝑉1 𝑉̂2
= {1 + [1 − ( ) ]} (10)
𝑃1 𝑉̂2 2𝛾 𝛾𝑃1 𝑉̂1

Pero:

𝑣1 2 𝐺 2 𝑉̂12 𝐺 2 𝑉̂1
𝑀𝑎 2 = ( ) = = (11)
𝑣2 𝛾𝑃1 𝑉̂1 𝛾𝑃1

Reemplazando ecuación (11) en (10):

𝑃2 𝑉̂1 𝛾−1 𝑉̂2

= {1 + 𝑀𝑎12 [1 − ( )]} (12)
𝑃1 𝑉̂2 2 𝑉̂1

2
𝑃2 𝑉̂2 𝑇2 𝛾−1 2
𝑉̂2
= = {1 + 𝑀𝑎1 [1 − ( ) ]} (13)
𝑃1 𝑉̂1 𝑇1 2 𝑉̂1

En las ecuaciones, las variables son: (𝑃1 𝑉̂1 ), (𝑃2 𝑉̂2 ) y 𝐺 (𝑜 𝑒𝑙 𝑀𝑎1 ), conocidas dos de ellas se
𝛾
puede calcular la tercera. A partir de la ecuación (9), tomando 𝑘𝑡𝑒 = 𝐶 ( )
𝛾−1

𝛾−1 𝐶 𝛾 − 1 𝑘𝑡𝑒 𝐺 2 𝑉̂
𝑃 = −( 𝐺𝑉̂ 2 + ) = − [ − ] (14)
2𝛾 𝑉̂ 𝛾 𝑉̂ 2
𝛾 − 1 𝑘𝑡𝑒 1 2
𝑑𝑃 = ( 2 − 𝐺 ) 𝑑𝑉̂ / ÷ 𝑉̂
2 𝑉̂ 2
𝑑𝑃 𝛾 − 1 𝑘𝑡𝑒 1 2 1
= [ 3 − 𝐺 ] 𝑑𝑉̂
𝑉̂ 2 𝑉̂ 2 𝑉̂
Integrando:
2
𝑑𝑃 𝛾 − 1 𝑘𝑡𝑒 1 1 1 𝑉̂2
∫ = [ ( 2 − 2 ) − 𝐺 2 𝑙𝑛 ( )] (15)
1 𝑉̂ 2 2 𝑉̂2 𝑉̂1 2 𝑉̂1

Sustituyendo el valor de la ctte de la ecuación (9):

2
𝑑𝑃 𝛾 − 1 𝐺 2 𝑉̂12 1 1 𝛾𝑃1 𝑉̂1 1 1 1 𝑉̂2
∫ = [ ( 2 − 2) + ( 2 − 2 ) − 𝐺 2 𝑙𝑛 ( )]
1 𝑉̂ 𝛾 4 𝑉̂2 𝑉̂1 2(𝛾 − 1) 𝑉̂2 𝑉̂1 2 𝑉̂1

2
𝑑𝑃 𝛾 − 1 2 𝑉̂12 𝑉̂2 𝑃1 𝑉̂1 1 1
∫ = 𝐺 ( 2 − 1 − 2𝑙𝑛 ( )) + ( 2 − 2) (16)
1 𝑉̂ 4𝛾 𝑉̂2 𝑉̂1 2 𝑉̂2 𝑉̂1

3
 PRQ 2202

Ahora, combinando las ecuaciones (4) y (5), se tiene:

𝑑𝑃 𝑑𝑉̂ 2𝑓𝐺 2 𝑑𝐿
+ 𝐺2 + =0 (17)
𝑉̂ 𝑉̂ 𝐷

Integrando esta ecuación:

2
𝑑𝑃 𝑉̂2 2𝑓𝐺 2 𝑑𝐿
∫ + 𝐺 2 𝑙𝑛 ( ) + =0 (18)
1 𝑉̂ 𝑉̂1 𝐷

Reemplazando la ecuación (16) en la (18):

𝑉̂2 𝛾 − 1 2 𝑉̂12 𝑉̂2 𝑃1 𝑉̂1 1 1 2𝑓𝐺 2 𝑑𝐿

𝐺 2 𝑙𝑛 ( ) + 𝐺 ( 2 − 1 − 2𝑙𝑛 ( )) + ( 2 − 2) + =0
𝑉̂1 4𝛾 𝑉̂2 𝑉̂1 2 𝑉̂2 𝑉̂1 𝐷

Reordenando:

2
4𝑓𝐿 1 2𝛾𝑃1 + (𝛾 − 1)𝐺 2 𝑉̂1 𝑉̂1 𝛾+1 𝑉̂2
= 𝑘 = 𝑒𝑣 = [ ] [1 − ( ) ] − 𝑙𝑛 ( ) (19)
𝐷 𝛾 2𝐺 2 𝑉̂1 𝑉̂2 2𝛾 𝑉̂1

Esta ecuación (19) expresa completamente las relaciones entre los volúmenes específicos en
función de las condiciones iniciales y el parámetro de fricción. Esta expresión permite calcular
𝑉̂2 , el volumen específico para el extremo de salida de la tubería, del fluido que circula con una
densidad de flujo másico G, desde una presión 𝑃1 , aguas arriba. La G se puede calcular
alternativamente en función del volumen específico del fluido por las 2 presiones, 𝑃1 𝑦 𝑃2 . La
presión 𝑃2 en el extremo de salida de la cañería se obtiene sustituyendo el valor de 𝑉̂2 de la
ecuación (12).

Para las condiciones constantes aguas arriba, el flujo máximo a través de la tubería se obtiene
𝜕𝐺
diferenciando con respecto a 𝑉̂2 y tomando ̂2
𝜕𝑉
= 0. Resulta por tanto, que el flujo máximo
ocurre cuando la velocidad para el extremo de salida es igual a la velocidad del sonido √𝛾𝑃2 𝑉̂2 .

La velocidad de flujo de un gas en condiciones adiabáticas nunca sobrepasa el 20% de la que se

obtiene para la misma diferencia de presiones en condiciones isotérmicas.
̂
𝑉
La solución simultanea de las ecuaciones (12), (13) y (19), eliminando (𝑉̂1), proporciona
2
𝑃2 𝑇2 4𝑓𝐿
relaciones de 𝑃 o también las relaciones de temperatura 𝑇 en función de 𝐷
.
1 1

4
 PRQ 2202

La ecuación (19) se puede expresar en función del número de Mach:

2
4𝑓𝐿 1 2 + (𝛾 − 1)𝑀𝑎12 𝑉̂1 𝛾+1 𝑉̂2
= 𝑘 = 𝑒𝑣 = [ ] [1 − ( ) ]− 𝑙𝑛 ( ) (20)
𝐷 2𝛾 𝑀𝑎1 2
𝑉̂2 2𝛾 𝑉̂1

Bibliografía
1. Coulson J. M., Richardson J. F.; Ingeniería Química, Ed, Reverte S.A., Barcelona, 1984.
2. Hougen O. A., Watson K. M., Ragatz R. A.; Principios de los procesos químicos, Tomo II,
Ed. Reverte S. A., Barcelona, 1974.
3. Bird R. B., Stewart W. E., Lightfoot E. M., Fenómenos de Transporte, Ed. Reverte S. A.,
Barcelona, 1982.