FLUJOS COMPRESIBLES ESTACIONARIOS

OBJETIVOS
- Presentar algunas caractersticas de los flujos compresibles.

Introduccin:
*Repaso de conceptos de termodinmica.
*Caractersticas de la dinmica de gases compresibles.

- Analizar flujos compresibles en conductos de seccin variable.

Desarrollo:
*Velocidad del sonido y velocidad de propagacin de
perturbaciones..
*Flujo estacionario en un gas compresible.
*Parmetros crticos.
*Efectos de la variacin del rea en flujos unidimensionales
compresibles.
*Flujo en una tobera convergente.
*Flujo en una tobera Laval.

Aspectos salientes de flujos Compresibles

Aspectos salientes de flujos Compresibles

Flujos Incompresibles >>> constante

Flujos Incompresibles >>> T constante

Flujos Compresibles >>> vara

Flujos Compresibles >>> T vara

Cambios en
la densidad
en funcin
del nmero
de Mach

La energa interna y la energa cintica

alcanzan similar importancia
Ma=2
En cintica= 2.3 105Joules
En Interna= 2.0 105 Joules

Es necesario incorporar en el anlisis la

ecuacin de conservacin de la energa

Descripcin de un GAS PERFECTO

h = ei +

ei = ei (T )
Esta ecuacin es vlida para todo proceso termodinmico ( isbaro,
isotermo, etc. )
Considerando esta definicin, el cambio en la energa interna en un
proceso cualquiera:
T2

ei = Cv dT
T1

de
calor especfico a volumen constante
Cv = i
dT V =cte

Tomando como hiptesis

gas perfecto:

Se define entalpa :

Se puede definir como gas perfecto a aquella sustancia cuya energa

interna es funcin solo de la temperatura.

= ei + RT
T2

El cambio de entalpa
para un gas perfecto:

h = C p dT
T1

Donde

h
Cp =

T p =cte

Existe una relacin entre

es el calor especfico a presin constante, y esta

expresin es vlida aun para procesos no isobaros.

Cp

Cv dada por:

Cv = Cv ( p, T )

C p Cv = R
Cvo =

Esta ecuacin es vlida para cualquier proceso, a volumen cte. o no.

ei = ei (T )

pV = NRT
p
R=
= cte
T

En general se pueden definir valores medios:

C po =

1
T

T2

1
T

T2

C dT
v

T1

C dT
p

T1

Surge la relacin entre

los calores especficos
que es independiente
de la temperatura:

k=

k = 1,67

Gases monoatmicos

k = 1,4

Mayora gases poli atmicos

Distintos tipos de procesos:

C po
Cvo

Cuando consideramos los cambios que ocurren durante un proceso,

para gases ideales se puede aplicar la ecuacin politrpica:

p p
n = n = cte
1 2

La entropa de un gas perfecto se puede expresar como :

p
S = ln k

p1/ k

= Cp ln

Se reconocen los siguientes casos:

Los gases perfectos tienen como gran ventaja que se conocen todas las relaciones
entre variables termodinmicas, y ellas adoptan una forma simple.

En virtud de lo anterior es posible la resolucin completa de las ecuaciones generales

de la dinmica de gases perfectos.

n=0
n =1

Procesos isotrmicos

n=k

Procesos isoentrpicos

Tambin se puede rescribir la ecuacin como:

p2

p1

Caractersticas de la dinmica de gases compresibles:

Cuando la velocidad de movimiento de un fluido se hace comparable a la velocidad
del sonido, los efectos inherentes a la compresibilidad del fluido avanzan al primer
plano.
Remarquemos que los flujos compresibles toman lugar a elevados nmeros de
Reynolds.
= viscocidad
La teora cintica de gases nos indica que:

cl

c = vel. sonido
l = libre camino medio

Procesos isobricos

n 1
n

= 2
1

n 1

T2
T1

Velocidad del sonido y propagacin del perturbaciones:

Analizamos la propagacin de una perturbacin en una geometra unidimensional:
Supongamos que el fluido esta
en reposo en el conducto y que
el mbolo se pone en
movimiento de repente con una
velocidad pequea (V)

r
V

r
c

0 +

0
p0

p0 + p

a) - OBSERVADOR ESTATICO

Si la velocidad caracterstica del sistema es del orden de la velocidad del sonido:

Lc

Re =

r
c
r
r
c V

valor que resulta muy alto casi siempre

Como para Re elevados los

efectos viscosos pueden ser
despreciados

El comportamiento del gas puede ser

descripto con las ecuaciones de un flujo
perfecto en la mayora de los casos.

r
c

0 +

p 0 + p

p0

b) Observador que se mueve con la onda de presin

A pesar de compartir las ecuaciones, los flujos compresibles son muy distintos de los
incompresibles.

Analicemos el movimiento desde un referencial que se mueve con

la onda generada. En este referencial el flujo es estacionario y
unidimensional. Planteando la ecuacin de continuidad:

V (t )

dV + (v n )dS = 0
S (t )
t

De lo anterior se obtiene:

V = c
(1)

SC c = ( + )SC (c V )

1
=c

+
+1

Surge entonces que si

los cambios en V son
pequeos, tambin los
cambios en la
densidad deben ser
pequeos.

Despreciando las fuerzas volumtricas y las fuerzas

viscosas se obtiene:

Observar entonces que a

menor V, menor p

SC c((c V ) c ) = pSC ( p + p )SC

Combinando (1) y (2):

c2 =

p = cV

(2)

p
1 +

Para movimientos que impliquen un / ms elevados tendremos que la velocidad

Consideremos ahora la ecuacin de conservacin de la cantidad de movimiento:

( )

( )(

de propagacin ser ms grande que la que se corresponde con / moderados .

( )

a
V (t ) t dV + S (t ) a w n dS + V (t ) adV = V (t ) f v dV + S (t ) n dS
Donde:

a=

velocidad absoluta = velocidad desde terna fija

w=

velocidad relativa = velocidad desde terna mvil

velocidad de giro de terna = 0

En el lmite de 0, tendremos lo que

hemos definido como velocidad del sonido:

c2 =

Analizamos a continuacin la propagacin de una perturbacin puntual en un gas

en movimiento.

o v
CASO N1

Entonces, en un flujo supersnico la perturbacin originada en un punto

cualquiera se propaga segn un cono cuyo ngulo de apertura decrece si
aumenta v.

Existen 2 configuraciones posibles:

Este ngulo se lo llama ngulo de Mach y la superficie del cono superficie

de Mach o superficie caracterstica.

cn

v<c

cn

v>c

CASO N2

En ambos casos la perturbacin es transportada por el fluido pero adems

es propagada por el mismo a una velocidad c en todas las direcciones.
En el CASO N1 la velocidad de propagacin de la perturbacin producida
en O (v + c.n) tiene distintos valores conforme a la direccin n
considerada. Los vectores v y c.n pueden aqu tener cualquier direccin en
el espacio. Esto significa que en un flujo subsnico la perturbacin se
propaga en todas direcciones.

Caso en que
la fuente se
mueve en
un fluido en
reposo

En el CASO N2 la direccin de los vectores v + c.n se circunscriben a un

cono de vrtice O cuyo ngulo de apertura mximo es: sin( ) = c
v

Flujo estacionario de un gas compresible:

Esta propiedad del flujo supersnico le confiere un carcter totalmente

distinto al de un flujo subsnico.

Cuando vimos la ecuacin de la conservacin de la energa para un flujo perfecto

obtuvimos:

En un flujo subsnico, la presencia de un objeto se hace sentir aguas

abajo y agua arriba del objeto.

Dh Dp
=
Dt Dt

Como el flujo es estacionario

uh = up

Considerando la ecuacin de EULER:

Du
1
= f v p
Dt

Por ser el flujo estacionario y

despreciando fuerzas en volumen:

En un flujo supersnico, el fluido ataca al objeto en forma ciega, la

influencia del objeto se produce solo aguas abajo, en tanto que en la zona
aguas arriba el fluido no se entera de la presencia del objeto.

Para un flujo irrotacional tenemos:

uu =

h =

p = h

u2

+ h = 0
2

Considerando:

uu =

u2
+ u
2

u2
+ h = cte
2

cte es la misma
para todos los
fluidos si es flujo
irrotacional.

Si el flujo es no irrotacional, esta expresin es vlida a lo largo de una lnea

de corriente donde u=0 se verifique siempre.

Si aplicamos esta expresin a un flujo compresible, tendremos entre 2 secciones

sobre una misma lnea de corriente:
2

h1 +
Considerando ahora que:

u1
u
= h2 + 2
2
2

u2 u1 =
2

h = C pT =

Si cambia v cambian los

parmetros termodinmicos y
viceversa.

Volviendo a la ecuacin (4) tenemos:

2

(5)

u1
u
+ C pT1 = 2 + C pT2 = C pT0
2
2

C p p0
R

J = u

(4)

Cp p
Cp
p
=
R C p Cv

T
2k p1 p2 2k
=
RT1 1 2
k 1 1 2 k 1
T1

u1max = 2C pT0 = 2

Consideremos a continuacin como es la variacin de densidad de flujo msico a lo

largo de una lnea de corriente:

(3)

u2 u1 = 2(h1 h2 ) = 2C p (T 1T2 )

h = C p T

Si es gas es perfecto, entonces

De la ecuacin (5):

T0 es la temperatura del punto

de estancamiento o de
reposo.

La ecuacin (3) muestra que la velocidad es mayor donde la entalpa es menor, para una lnea de
corriente dada o para el todo el flujo en caso irrotacional.

La ecuacin de Euler nos permite escribir a lo largo de una lnea de corriente que la
relacin entre du y dp es:

udu =

d
u
= 2
du
c

dp

Recordando que:

Considerando que el diferencial

u2
dJ d (u )
=
= 1 2
du
du
c

dp
= c2
d

d (u ) = du + ud

Lo que muestra que en flujo subsnico (Ma<1) si nos

desplazamos a lo largo de una lnea de corriente la
densidad del flujo crece si aumenta la velocidad.

En el dominio de flujos supersnicos, la densidad de flujo decrece al aumentar la

velocidad.
Esta diferencia de comportamiento puede interpretarse como que en flujos
subsnicos las lneas de corriente se acercan al aumentar la velocidad y divergen
para flujos supersnicos.

Parmetros crticos:
J

p = cte. k

Si el proceso es isoentrpico:

J*
1.0
Para el caso de gases perfectos:

p
p
= cte.k . k 1 = k .

p
p
= k . = kRT

c = kRT = k

Se ve que la velocidad de propagacin del sonido es funcin de la Temperatura

absoluta.

u
1.0

2.5

La velocidad del sonido en un flujo supersnico puede diferir de un punto a otro en

funcin de los valores de estas magnitudes termodinmicas locales.

c*

El flujo J adopta su valor mximo cuando la velocidad del gas iguala la velocidad
del sonido local.
En ese caso:

J = J * = *c*

Considerando:

El asterisco indica magnitudes correspondientes a este punto.

T
=
0 T

k 1

p k 1
=
p0 0

Se obtiene las siguientes expresiones

T* =

Considerando que:

como
2

c =
2

dh = C p dT
2

u1
= h0
2
y que sobre una lnea de corriente:
2
u2
2
2
h2 +
= h0
k p
c
2
h = C pT =
=
h1 +
2

k 1

k 1

c*
c
c
+ * = 0
k 1 2 k 1

0 pto de reposo
o estacionario

c* = c0

de lo anterior

= 0 1

Independiente
de la geometra
del tubo.

2
k +1

k 1 u 2
k 1 u2
= T0 1
T = T0 1
k + 1 c 2
2 c0 2

k 1 u2
k 1 u 2 k 1
= 0 1
k + 1 c 2
2 c0 2

k 1 u 2

p = p0 1
2 c0 2

k
k 1

k 1 u 2

= p0 1
2
k +1 c

No ocurre lo mismo para un flujo compresible donde:

.u. A = cte

Vamos a ver a continuacin que el comportamiento de un fluido

compresible depende del valor que adopta el valor del nmero de Mach.
Considerando la ecuacin de continuidad,
tomando logaritmos y diferenciando tendramos:
Para flujos isoentrpicos:

c2 =

dp
d

1 du 1 dA 1 d
+
+
=0
u dx A dx dx

X es la direccin de la lnea de corriente

1 du 1 dA
1 dp
+
+
=0
u dx A dx c 2 dx
Considerando la ecuacin de Euler sobre la lnea de corriente:

1 du 1 dA u du
+
=
u dx A dx c 2 dx

3 FLUJO SNICO ( Ma = 1 )

u
dp
=
dx
x

1 du
1 dA
1 Ma 2 =
u dx
A dx

Es necesario que

dA
=0
dx

Esta condicin se produce en la garganta de un conducto conv div o div conv.

Sin embargo el flujo ser snico solo si el p entre aguas arriba y aguas debajo de
la garganta es suficientemente grande como para acelerar el fluido, caso contrario,
en la garganta se tiene una situacin de Ma < 1.

2 k 1
p* = p0
= p0 .0,528
k +1
1

2 k 1
= 0 .0,634
k +1

* = 0

k
k 1

Si u = c* se
obtienen los
parmetros
crticos:

Valores para el AIRE

Los parmetros crticos son

independientes de la geometra.

1 du
1 dA
1 Ma 2 =
u dx
A dx

Efectos de la variacin del rea en flujo unid. Compresible:

En flujo unidimensional incompresible la velocidad de un fluido varia
inversamente proporcional al rea transversal de la corriente:
u. A = cte

1
k 1

2T0
= T0 .0,83
k +1

1 FLUJO SUBSNICO ( Ma < 1 )

dA
<0
dx

du
>0
dx

dA
>0
dx

du
<0
dx

u
p

2 FLUJO SUPERSNICO ( Ma > 1 )

dA
<0
dx

du
<0
dx

dA
>0
dx

du
>0
dx

u
p

Flujo compresible a travs de una tobera convergente:

Supongamos el flujo compresible a travs de la tobera de la figura
conectada en su entrada a un recipiente muy grande donde se dan las
condiciones de reposo:

Los posibles flujos en canales o

conductos convergentes
divergentes se indican en la figura.
Estos flujos se logran si las
presiones de ensayo se ajustan
correctamente.
Para flujos que tengan lugar en toberas
en turbomquinas, existen 2 problemas:
1- Se conoce A(x) y se desconocen las
magnitudes del flujo en funcin de x
(problema directo).
2- Se especifica u(x) y se desea
determinar A(x) (problema indirecto).

Si la velocidad inicial es despreciable, la presin en un punto cualquiera de

la seccin (1) va a prcticamente coincidir con la presin de reposo (o
estacionaria) p0 .

Entonces tendramos:

u2 =
2

2 k p0 p 2

k 1 1 2

Si el proceso es isoentrpico:

p0

1k

p2

2k

Supongamos ahora que Ma2=1,0 sustituyendo para la relacin de presiones

entrada - salida nos queda:

c2 = k

p2

k 1

u2
2 p0 k
2

c = Ma2 = k 1 p 1
2
2

p2 2 k +1
=

p0 k + 1

p
1 = 0 2
p2

Y entonces se obtiene la formula de Saint Venant

Wantzel:
Teniendo en cuenta que:

1
k

u2 =
2

2 k p2
k 1 2

p0
p2

k 1
k

Como

Independiente de la geometra del conducto.

J . A = m& = cte J max . Amin = cte

Por lo tanto Jmax solo puede

ocurrir donde la garganta es
mnima, en Amin.

Entonces el nico punto donde la velocidad puede ser snica es

en la garganta ya que
J max J *

Se obtiene:

Existen 2 configuraciones posibles para toberas:

Si la diferencia de presiones
entre p0 y la presin de salida es
suficientemente grande, se
alcanza la velocidad snica en la
seccin 2.

Si:

patm > p*
p0 = patm

Si el flujo en la garganta es subsnico la presin en la garganta coincide

con la presin aguas debajo de la tobera.
Si el flujo en la garganta es snico la presin en la garganta puede ser =
generalmente mayor que la presin aguas abajo.

El flujo msico en la tobera es

Si:

patm < p*
p0 = patm

m& subs = A2
Expansin,
ondas despus
de la tobera.

k 1
2

p k p k
2k
p0 0 2 2
p0
k 1
p0

Expresin vlida
mientras la velocidad
en la garganta sea
menor que la snica.

Si la velocidad en la garganta es snica entonces:

k +1

El valor mximo de p
Para el aire

m& = uA = J . A

Reemplazando en S. V. W. tendremos si la seccin de la garganta es Ai ( p2>p* ):

p = p0 pc = p0 p* = p0 1

k +1

p* = 0,53. p0

2 k 1
m& son = A2 kp0 0
= J * . A2
k +1

k
k 1

Expresin que corresponde

al caudal mximo y puede
rescribirse como

k +1

m& son =

A2 p0
T0

k 2 k 1

R 0 k +1

Estas expresiones son independientes de la presin en la atmsfera a

la que sale el flujo.

pmax = 0,47. p0

El caudal depende del p entre el recipiente y el ambiente hasta que

se alcanza el valor mximo.
Mas all de este valor, el flujo msico no depende de p, por lo que
disminuir la presin ambiente por debajo de p* no implica mayores
caudales. Se dice que la tobera se atora.

Caso en que la presin

ambiente coincide con la
presin en el extremo de la
tobera.

Toberas convergentes divergentes (toberas Laval):

Las toberas convergentes divergentes surgen como consecuencia de la
inaccesibilidad de velocidades supersnicas al utilizar toberas del tipo convergentes
(donde la velocidad del sonido solo se alcanza en la salida).

Si la presin ambiente es aumentada, surge

entonces lo que se suele llamar chorro
sobreexpandido. El gas en la tobera se
expande desde una presin pe<patm

J*
1.0

Inicialmente el flujo en la tobera no cambia. Afuera

de la tobera el flujo no es mas unidimensional. La
descripcin de este flujo requiere del concepto de
onda de choque (que trataremos la clase que viene).

x
Si en la garganta se llega a flujo snico en el difusor se produce disminucin de
presin y aumento de velocidad.

k 1
1

p k
2 k p0 p k
J = u = 0
1
k 1 0 p0
p0

1/ 2

Por ahora nos contentamos diciendo que la

onda de choque es una superficie donde la
presin y la temperatura son discontinuas.
Esta onda de choque emana del borde de la
tobera, permitiendo aumentar la presin de
descarga a la presin ambiente. Las ondas
de choque en este caso no son simples, ya
que hay reflexiones y formacin de ondas
secundarias, llamadas ondas de expansin.

Si la presin ambiente se aumenta aun ms, la onda se mueve hacia el

interior de la tobera y forma la onda de choque normal (curva 2).

Aguas abajo de esta onda de choque normal, el flujo es subsnico.

Esta parte de la tobera funciona entonces como un difusor subsnico, el cual
tericamente aumenta la presin y disminuye la velocidad.
En realidad existe separacin de la capa lmite y la presin despus de la onda de
choque permanece casi igual a la presin ambiente.
Si la presin ambiente aumenta aun ms, la onda de choque se introduce aun ms en
la tobera, pero se debilita ya que el nmero de Ma delante de la onda es menor.
Si la presin aumenta aun ms, la onda de choque alcanza finalmente la garganta, la
intensidad de la onda o del salto disminuye a cero y el conjunto del flujo se vuelve
subsnico. (curva 5 de la pgina anterior).
Si aun incrementamos la presin ambiente, el nmero de Ma es mximo en la
garganta pero nunca llega a 1 (curva 4 de la pgina anterior).
El chorro se lo llama subexpandido si la presin de salida es mayor que la presin
ambiente:

CONCLUSIONES
La reduccin a la presin ambiente se logra a partir de ondas de
expansin que al reflejarse provocan ondas de compresin. El
flujo en la tobera permanece sin embargo inalterado.
En las toberas convergentes divergentes tambin aparece el
fenmeno de atoramiento.

Como en la garganta

Hemos visto que la velocidad de propagacin de perturbaciones depende de la

magnitud de la misma, y que si las perturbaciones son pequeas el fenmeno
difiere segn que la propagacin pueda ser en una o en todas las direcciones. En
este caso la velocidad del fluido limita las direcciones en las cuales se puede
propagar la perturbacin.
Despreciando los efectos viscosos, hallamos una relacin muy til que vincula la
velocidad con la entalpa en distintos puntos del fluido.

u c*

J max = J * = .c*

En esta clase nos concentramos en analizar algunas caractersticas salientes de

los flujos compresibles.

Los flujos compresibles con seccin del conducto variable presentan

comportamiento muy distinto al de los flujos incompresibles.

m& max = J * S g arg = J * .A2

Entonces en tanto que el flujo sea supersnico, el caudal queda

limitado por el atoramiento. Las expresiones de los caudales
son entonces las mismas que para los flujos en toberas
convergentes.

Los flujos snicos presentan comportamiento inverso en lo que respecta a presin

y velocidad, que los subsnicos.
La velocidad mxima en toberas convergentes est limitada a la velocidad snica
local funcin de las condiciones de entrada en la tobera.
En las toberas Laval el flujo puede ser supersnico, pero un fenmeno adicional de
saltos bruscos de presin puede tomar lugar.
En ambas toberas puede ocurrir el fenmeno de atoramiento donde el caudal
msico adopta un valor que solo depende de las condiciones de entrada.