FLUJO COMPRESIBLE

Los flujos compresibles se encuentran con frecuencia en dispositivos que

incluyen el flujo de gases a altas velocidades. Los flujos compresibles combinan la
dinámica de fluidos y la termodinámica, ambas son absolutamente necesarias para el
desarrollo de los fundamentos teóricos necesarios.

Para conocer más sobre este tema debemos conocer y tener claro que es el
estado de estancamiento, velocidad del sonido y número de Mach para flujos
compresibles. Las relaciones entre las propiedades estáticas y de estancamiento se
desarrollan para flujos isentrópicos de gases ideales, y se expresan en función de la
razón de calores específicos y el número de Mach. Se tratan los efectos de los
cambios del área en flujos isentrópicos unidimensionales subsónicos y supersónicos.
Estos efectos se ilustran considerando el flujo isentrópico a través de toberas
convergentes y toberas convergentes divergentes. Además debemos conocer las
ondas de choque y la variación de las propiedades de flujo a través de ondas de
choque normales y oblicuas. Y a su vez, consideraremos los efectos de la fricción y
la transferencia de calor en flujos compresibles.

PROPIEDADES DE ESTANCAMIENTO

Cuando se analizan volúmenes de control

es muy conveniente combinar la energía
interna y la energía de flujo de un fluido
en un solo término: entalpía, definida por
unidad de masa como h = u + P/ρ. Cada
vez que la energía cinética y la energía
potencial de un fluido son insignificantes, Figura 1.1. a)
como sucede con frecuencia, la entalpía representa la energía total de un fluido.
Para flujos a altas velocidades, como los que fluyen en motores de propulsión

Rogmer D. Guevara F. Ingeniero Mecánico

 (Figura 1.1 a) b) Los aviones y los motores
de propulsión implican altas velocidades, y
por esta razón, al analizarlos, el término de
la energía cinética siempre debe tomarse en
cuenta.), la energía potencial del fluido es
insignificante, pero la energía cinética no lo
Figura 1.1. b) es. En tal caso, es conveniente combinar la
entalpía y la energía cinética del fluido en un solo término llamado entalpía de
estancamiento (o entalpía total) h0, definido por unidad de masa como:

𝑉2 𝑘𝐽
ℎ𝑜 = ℎ + ( ⁄𝑘𝑔)
2

Cuando la energía potencial de un fluido es insignificante, la entalpía de

estancamiento representa la energía total de un flujo de fluido por unidad de masa.
De esta manera se simplifica el análisis termodinámico de flujos a altas velocidades.

Para este estudio citaremos como entalpía estática la entalpía normal h —

cada vez que sea necesario— para distinguirla de la entalpía de estancamiento.
Observe que la entalpía de estancamiento es una combinación de las propiedades
de un fluido, tal como la entalpía estática. Ambas entalpías son idénticas cuando
la energía cinética de un fluido es despreciable. Considere el flujo estacionario de un
fluido a través de un ducto, tal como una
tobera, un difusor, o cualquier otro conducto
de flujo en el cual el flujo es adiabático y
donde no se realiza el trabajo de flecha o
trabajo eléctrico (Figura 1.2. Flujo
Figura 1.2.
estacionario de un fluido a través de un ducto
adiabático.). Al suponer que el fluido experimenta pequeños cambios o ninguno en
su elevación y en su energía potencial, el balance de energía (𝐸̇𝑒𝑛𝑡 = 𝐸̇𝑠𝑎𝑙 ) para este
dispositivo de flujo estacionario de una entrada y una salida se reduce a:

Rogmer D. Guevara F. Ingeniero Mecánico

 𝑉12 𝑉22
ℎ1 + = ℎ2 + o ℎ01 = ℎ02
2 2

Entonces, en ausencia de cualquier interacción de calor y de trabajo, así como

de cualquier cambio en la energía potencial, la entalpía de estancamiento de un
fluido permanece constante durante un proceso de flujo estacionario. Los flujos a
través de toberas y difusores por lo general satisfacen estas condiciones, y cualquier
aumento en la velocidad del fluido en estos dispositivos origina un decremento
equivalente en la entalpía estática del fluido.

Si el fluido se lleva al reposo, entonces la velocidad del estado 2 es 0 y la

𝑉12
ecuación es: ℎ1 + = ℎ2 = ℎ02
2

Así, la entalpía de estancamiento representa la entalpía de un fluido cuando

se lleva al reposo adiabáticamente.

Durante el proceso de
estancamiento, la energía cinética
de un fluido se convierte en entalpía
(energía interna energía de flujo), la
cual da como resultado un aumento
en la temperatura y la presión del
fluido (Fig. 1.3. La energía cinética
Figura 1.3. se convierte en entalpía durante un
proceso de estancamiento). Las propiedades de un fluido en estado de estancamiento
se llaman propiedades de estancamiento (temperatura de estancamiento, presión
de estancamiento, densidad de estancamiento, etc.). El estado de estancamiento y las
propiedades de estancamiento se indicarán con el subíndice 0. El estado de
estancamiento se llama estado de estancamiento isentrópico cuando el proceso de
estancamiento es reversible y adiabático (es decir, isentrópico).

Rogmer D. Guevara F. Ingeniero Mecánico

 La entropía de un fluido permanece
constante durante el proceso isentrópico de
llevar el fluido al estado de estancamiento. El
proceso real (irreversible) y el proceso
isentrópico de llevar al reposo un flujo de
fluido se ilustra en el diagrama h-s (Figura
1.4. El estado real, el estado de
estancamiento real y el estado de
estancamiento isentrópico de un fluido en un
diagrama h-s) Observe que la entalpía de
Figura 1.4. estancamiento del fluido (y la temperatura de
estancamiento si el fluido es un gas ideal) es la misma para ambos casos. Sin
embargo, la presión de estancamiento real es menor que la presión de estancamiento
isentrópica porque la entropía aumenta durante el proceso real de estancamiento
como resultado de la fricción del fluido. Con frecuencia, los procesos de
estancamiento se aproximan a isentrópicos y a las propiedades de estancamiento
isentrópico se les llama simplemente propiedades de estancamiento.

Cuando el fluido se aproxima como un gas ideal con calores específicos

constantes, su entalpía puede reemplazarse por cpT y puede expresarse como:

𝑉2 𝑉2
𝐶𝑝 𝑇0 = 𝑐𝑃 𝑇 + o 𝑇0 = 𝑇 +
2 2𝐶𝑝

En ésta, T0 se llama temperatura de estancamiento (o temperatura total), y

representa la temperatura que alcanza un gas ideal cuando se lleva al reposo
adiabáticamente. El término V2/2cp corresponde al incremento de la temperatura
alcanzado durante tal proceso y se llama temperatura dinámica. Por ejemplo,
la temperatura dinámica del aire que fluye a 100 m/s es (100 m/s)2/(2 x 1.005
kJ/kg · K) = 5.0 K.

Rogmer D. Guevara F. Ingeniero Mecánico

 Por lo tanto, cuando el aire a 300 K y 100 m/s
se lleva al reposo adiabáticamente (en la punta
de la sonda de medición de temperatura, por
ejemplo), su temperatura alcanza el valor de
estancamiento de 305 K (Figura 1.5). Observe
que para los flujos a bajas velocidades, las
temperaturas de estancamiento y estática (o
normal) son prácticamente iguales. Pero, para
flujos a altas velocidades, la temperatura medida por una sonda en reposo colocada
en el fluido (temperatura de estancamiento) puede ser considerablemente mayor que
la temperatura estática del fluido.

La presión que alcanza un fluido cuando se lleva al reposo isentrópicamente

se llama presión de estancamiento P0. Para un gas ideal con calores específicos
constantes, P0 está relacionado con la presión estática del fluido mediante:

𝑃0 𝑇0 𝑘 /(𝑘−1)
=( )
𝑃 𝑇

Se observa que ρ = 1/v y al utilizar la relación isentrópica Pvk = P0vk0, el cociente

entre la densidad de estancamiento y la densidad estática pueden expresarse como:

𝜌0 𝑇0 𝑘 /(𝑘−1)
=( )
𝜌 𝑇

Cuando se usan entalpías de estancamiento no es necesario referirse a la energía

cinética de manera explícita. Entonces el balance de energía 𝐸̇𝑒𝑛𝑡 = 𝐸̇𝑠𝑎𝑙 para el
dispositivo de flujo estacionario con una entrada y una salida puede expresarse
como:

𝑞𝑒𝑛𝑡 + 𝑊𝑒𝑛𝑡 + (ℎ01 + 𝑔𝑧1 ) = 𝑞𝑠𝑎𝑙 + 𝑊𝑠𝑎𝑙 + (ℎ01 + 𝑔𝑧2 )

Rogmer D. Guevara F. Ingeniero Mecánico

donde h01 y h02 son las entalpías de estancamiento en los estados 1 y 2,
respectivamente. Cuando el fluido es un gas ideal con calores específicos constantes
se convierte en:

(𝑞𝑒𝑛𝑡 − 𝑞𝑠𝑎𝑙 )(𝑊𝑒𝑛𝑡 − 𝑊𝑠𝑎𝑙 ) = 𝐶𝑝 (𝑇02 − 𝑇01 ) + 𝑔(𝑧2 − 𝑧1 )

donde T01 y T02 son las temperaturas de estancamiento.

VELOCIDAD DEL SONIDO Y NÚMERO DE MACH

Un parámetro importante en el estudio de flujo compresible es la velocidad

de sonido (o velocidad sónica), que es la velocidad a la cual una onda de presión
infinitesimalmente pequeña viaja a través de un medio. La onda de presión puede ser
originada por un pequeño disturbio, el cual crea un ligero aumento en la presión
local.

Para obtener la expresión para la velocidad de

sonido en un medio, considere un conducto lleno
de fluido en reposo, como se muestra en la figura
1.6 (Propagación de una pequeña onda de presión
a lo largo de un ducto.). Un émbolo colocado en el
conducto se mueve a la derecha con una velocidad
infinitesimal constante dV y crea una onda sónica.
El frente de onda se mueve a la derecha a través del
fluido a la velocidad del sonido c y separa el fluido
adyacente al pistón que ya está en movimiento del
fluido que aún está en reposo. El fluido a la
Figura 1.6. izquierda del frente de onda experimenta un
cambio infinitesimal de sus propiedades termodinámicas, mientras que el fluido a la
derecha del frente de onda mantiene sus propiedades termodinámicas originales,
como se muestra en la figura 1.6.

Rogmer D. Guevara F. Ingeniero Mecánico

Para simplificar el análisis, considere un
volumen de control que encierra al frente de
onda y que se mueve con él, como se muestra
en la figura 1.7 (Volumen de control que se
mueve con la pequeña onda de presión en un
ducto.) Para un observador que viaja con el
frente de onda, el fluido a la derecha parecería Figura 1.7.
moverse hacia el frente de onda con una velocidad c y el fluido a la izquierda
parecería alejarse del frente de onda con una velocidad c – dV. Por supuesto, el
observador pensará que el volumen de control que encierra el frente de onda (y a él
mismo o ella misma) está en reposo, y el observador es testigo de un proceso
estacionario. El balance de masa para este proceso de flujo estacionario de una
entrada y una salida puede expresarse como:

𝑚̇𝑑𝑒𝑟 = 𝑚̇𝑖𝑧𝑞 o 𝜌𝐴𝑐 = (𝜌 + 𝑑𝜌)𝐴(𝑐 − 𝑑𝑉)

Al cancelar el área de la sección transversal A y despreciar los términos de

orden superior; esta ecuación se reduce a: 𝑐 𝑑𝜌 − 𝜌𝑑𝑉 = 0 (𝒂)

Si no hay transferencia de energía en forma de calor o trabajo a través de las

fronteras del volumen de control durante este proceso estacionario, y el cambio en la
energía potencial puede despreciarse, entonces el balance de energía para el proceso
estacionario eent – esal se convierte en:

𝑐2 (𝑐 − 𝑑𝑉 )2
ℎ+ = ℎ + 𝑑ℎ +
2 2

que da como resultado:

𝑑ℎ − 𝑐 𝑑𝑉 = 0 (𝒃)

donde se ha despreciado el término de segundo orden dV2. La amplitud de la onda

sónica que se presenta normalmente es muy pequeña y no causa ningún cambio

Rogmer D. Guevara F. Ingeniero Mecánico

apreciable en la presión y la temperatura del fluido. Por eso, la propagación de una
onda sónica no es solamente adiabática, también es casi isentrópica. Entonces, la
relación termodinámica T ds = dh – dP/ρ, se reduce a:

𝑑𝑃
𝑑ℎ =
𝜌

Con la combinación de las ecuaciones a, b y c se obtiene la expresión deseada para

la velocidad del sonido como:

𝑑𝑃 𝜕𝑃
𝑐2 = 𝑎 𝑠 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 o 𝑐2 = ( )
𝑑𝜌 𝜕𝜌 𝑥

Se deja como ejercicio al lector demostrar que al usar las relaciones de las
propiedades termodinámicas, que la ecuación antes planteada también puede
𝜕𝑃
escribirse como: 𝑐2 = 𝑘 ( )
𝜕𝜌 𝑇

donde k es la razón de calores específicos del fluido. Observe que la velocidad del
sonido en un fluido es una función de las propiedades termodinámicas de éste.

Cuando el fluido es un gas ideal (P = ρRT), la diferenciación con la última ecuación

planteada puede hacerse con facilidad para obtener:

𝜕𝑃 𝜕(𝜌𝑅𝑇)
𝑐2 = 𝑘 ( ) = 𝑘 [ ] = 𝑘𝑅𝑇 o 𝑐 = √𝑘𝑅𝑇
𝜕𝜌 𝑇 𝜕𝜌 𝑇

Observe que la constante del gas R tiene valor

fijo para un gas ideal en particular y que la razón de
calores específicos k de un gas ideal es, a lo mucho,
una función de la temperatura, se observa que la
velocidad del sonido en un gas ideal dado es función
de la temperatura solamente (Fig. 1.8. La velocidad
del sonido varía con la temperatura y con el fluido.).

Figura 1.8.

Rogmer D. Guevara F. Ingeniero Mecánico

 Un segundo parámetro importante en el análisis de flujo compresible es el
número de Mach Ma, llamado así en honor al físico austriaco Ernst Mach (1838-
1916). Es el cociente de la velocidad real del fluido (o de un objeto que se mueve en
el fluido en reposo) entre la velocidad del sonido en el mismo fluido, en el mismo
𝑉
estado: 𝑀𝑎 =
𝑐

Observe que el número de Mach depende de la velocidad del sonido, la cual

depende del estado del fluido. Por lo tanto, el número de Mach de un avión que
vuela a velocidad constante a través de aire quieto puede ser diferente en regiones
distintas (Fig.1.9. El número de Mach puede ser diferente a temperaturas distintas,
inclusive si la velocidad de vuelo es la misma.). Con frecuencia, el régimen de flujos
se describe en términos del número
de Mach. El flujo se llama sónico
cuando Ma = 1, subsónico cuando
Ma < 1, supersónico cuando Ma > 1,
hipersónico cuando Ma >> 1 y
transónico cuando Ma ≅ 1.
Figura 1.9.

ONDAS DE CHOQUES Y ONDAS DE EXPANSIÓN

Se ha visto que las ondas de propagación de sonido se deben a cambios de

presión infinitesimalmente pequeños, y a que éstas viajan a través del medio a la
velocidad del sonido. También se ha visto que para algunos valores de presión del
receptor ocurren cambios abruptos en las propiedades de los fluidos en una sección
muy delgada de una tobera convergente-divergente en condiciones de flujo
supersónico y crean una onda de choque. Es interesante estudiar las condiciones en
las cuales se forman las ondas de choque y cómo afectan al flujo.

Rogmer D. Guevara F. Ingeniero Mecánico

CHOQUES NORMALES

Se considera primero a las ondas de choque que ocurren en un plano normal a

la dirección de flujo, llamadas ondas de choque normales. El proceso de flujo a
través de la onda de choque es muy irreversible y no puede aproximarse como un
proceso isentrópico.

A continuación seguimos los pasos de Pierre Laplace (1749-1827), G. F.

Bernhard Riemman (1826-1866), Willian Rankine (1820-1872), Pierre Henry
Hugoniot (1851-1887), Lord Rayleigh (1842-1919) y G. I. Taylor (1886-1975) y
desarrollaremos las relaciones para las
propiedades de flujo antes y después Figura 1.10.

del choque. Para ello aplicaremos las

ecuaciones de la conservación de masa,
cantidad de movimiento y energía, así
como algunas relaciones entre las
propiedades del flujo en el volumen de
control que contiene la onda de choque,
como se muestra en la figura 1.10.
(Volumen de control para el flujo a través de una onda de choque normal.)

Las ondas de choque normales son extremadamente delgadas, de tal manera

que las áreas de flujo entrante al volumen de control y saliente de éste son casi
iguales (Figura 1.11. Imagen de Schlieren (estriograma)
de un choque normal en una tobera de Laval. El número
de Mach en la tobera corriente arriba (a la izquierda) de
la onda de choque es casi 1.3. Las capas límite
distorsionan la forma de la onda de choque normal en
cercanía a las paredes y llevan el flujo a separarse por

Figura 1.11. abajo de la onda de choque).

Rogmer D. Guevara F. Ingeniero Mecánico

 Se considera flujo estacionario sin transferencia de calor o interacciones de
trabajo y sin cambios de energía potencial. Se denotan las propiedades del flujo
inmediatamente corriente arriba de la onda de choque con el subíndice 1 y las
propiedades del flujo inmediatamente corriente abajo, con el subíndice 2, y se tiene
lo siguiente:

Conservación de masa:

ρ1 A𝑉1 = ρ2 A𝑉2 o ρ1 𝑉1 = ρ2 𝑉2

Conservación de energía:

𝑉12 𝑉22
h1 + = h2 + o h01 = h02
2 2

Conservación de cantidad de movimiento:

A(𝑃1 − 𝑃2 ) = 𝑚̇ (𝑉2 − 𝑉1 )

Incremento de la entropía: 𝑠2 − 𝑠1 ≥ 0

Las relaciones de conservación de masa y energía pueden combinarse en una

sola ecuación y graficar en un diagrama h-s si se utilizan las relaciones entre las
propiedades. La curva resultante se llama línea de Fanno, y a lo largo de esta curva
se localizan los estados que tienen el mismo valor de entalpía de estancamiento y
flujo de masa por unidad de área. De la misma manera, al combinar las ecuaciones
de conservación de masa y de cantidad de movimiento en una sola ecuación y
graficarla en un diagrama h-s se obtiene la curva llamada línea de Rayleigh. Ambas
líneas se muestran en el diagrama h-s de la figura 1.12 (Variación de las
propiedades normalizadas de un fluido y del área de sección transversal en un
ducto con la caída de presión de 1 400 a 200 kPa.)

Rogmer D. Guevara F. Ingeniero Mecánico

 Figura 1.12.

Las líneas de Fanno y de Rayleigh se intersecan en dos puntos (puntos 1 y 2),

que representan los dos estados donde las tres ecuaciones de conservación se
satisfacen. Uno de éstos (estado 1) corresponde al estado antes del choque, y el otro
(estado 2)
corresponde al estado
después del choque.
Se observa que el
flujo es supersónico
antes del choque y
subsónico después.
Por lo tanto, el flujo
cambiará de
supersónico a
subsónico si ocurre
una onda de choque
normal. Cuanto Figura 1.13.

Rogmer D. Guevara F. Ingeniero Mecánico

mayor sea el número de Mach antes del choque, más fuerte será el choque. En el
caso límite de Ma = 1, la onda de choque simplemente se convierte en una onda de
propagación de sonido. Se muestra en la figura 1.13. (Diagrama h-s para el flujo a
través de un choque normal) que la entropía aumenta: S2 > S1. Esto se espera porque
el flujo a través del choque es adiabático, pero irreversible.

El principio de la conservación de energía exige que la entalpía de

estancamiento permanezca constante durante el choque; h01 = h02. Para gases ideales
h = h(T), y así: T01 = T02

Esto significa que la

temperatura de estancamiento de
un gas ideal también permanece
constante durante el choque. Sin
embargo, se nota que la presión de
estancamiento disminuye durante
el choque debido a las
irreversibilidades, mientras que la
temperatura normal (estática) Figura 1.14.

aumenta drásticamente debido a la conversión de energía cinética en entalpía y

causa un gran descenso en la velocidad del fluido (Figura 1.14.).

Ahora se desarrollarán las relaciones entre varias propiedades antes y después

del choque para un gas ideal con calores específicos constantes. Una relación para la
razón de las temperaturas estáticas T2/T1 se obtiene:

𝑇01 𝑘−1 𝑇02 𝑘−1

=1+( ) 𝑀𝑎12 y =1+( ) 𝑀𝑎22
𝑇1 2 𝑇2 2

Se divide la primera ecuación entre la segunda y al notar que T01 = T02, se

tiene:

Rogmer D. Guevara F. Ingeniero Mecánico

 𝑇1 1 + 𝑀𝑎12 (𝑘 − 1)/2
=
𝑇2 1 + 𝑀𝑎22 (𝑘 − 1)/2

A partir de la ecuación de estado del gas ideal:

𝑃1 𝑃2
𝜌1 = y 𝜌2 =
𝑅𝑇1 𝑅𝑇2

Se sustituyen éstas en la ley de conservación de masa ρ1V1 = ρ2V2 y al notar que

Ma V/c y c = √𝑘𝑅𝑇, se tiene que:

La razón de presiones a través del choque:

Esta es una combinación de las ecuaciones de conservación de la masa y la

energía; por lo tanto, es también la ecuación de la línea de Fanno para un gas ideal
con calores específicos constantes. Una relación similar para la línea de Rayleigh
puede obtenerse si se combinan las ecuaciones de conservación de la masa y la
ecuación de cantidad de movimiento.

Combinando ecuaciones obtenemos:

Rogmer D. Guevara F. Ingeniero Mecánico

 Esta ecuación representa la intersección de las líneas de Fanno y Rayleigh y
relaciona el número de Mach corriente arriba del choque al número de Mach
corriente abajo del choque.

El choque de ondas no se limita solamente a toberas supersónicas. Este

fenómeno también se observa en la entrada del motor de un avión supersónico,
donde el aire pasa por un choque y desacelera a velocidades subsónicas antes de
entrar al difusor del motor. Las explosiones también producen la propagación de
ondas de choque normales muy poderosas que pueden ser muy destructivas.

Existe una tabla donde se indican razones de varias propiedades de flujo

después y antes del choque para un gas ideal con k = 1.4. Al revisar esta tabla se
encuentra que Ma2 (el número de Mach después del choque) es siempre menor que 1
y cuanto mayor es el número de Mach supersónico antes del choque, menor es el
número de Mach subsónico después del choque. También se observa que la presión
estática, la temperatura y la densidad aumentan después del choque, mientras que la
presión de estancamiento disminuye.

El cambio de entropía en el choque se obtiene cuando se aplica éste a la

ecuación del cambio de entropía de un gas ideal:

la cual puede expresarse en términos de k, R y Ma1 al incluir las relaciones

desarrolladas antes en esta sección. La gráfica del cambio de entropía adimensional
en el choque normal (s2 – s1)/R contra Ma1, se muestra en la figura 1.15 (Cambio de
entropía a través del choque normal). Puesto que el flujo a través de la onda de
choque es adiabático e irreversible, la segunda ley de termodinámica exige que el
flujo aumente su entropía al atravesar la onda de choque. Por lo tanto, una onda de
choque no puede existir para valores de Ma1 menores que la unidad, donde el

Rogmer D. Guevara F. Ingeniero Mecánico

cambio de entropía sería negativo. Para flujos adiabáticos, las ondas de choque
existen solamente para flujos supersónicos, Ma1 > 1.

Figura 1.15

Rogmer D. Guevara F. Ingeniero Mecánico