UNIVERSIDAD DE CARABOBO

FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA MECÁNICA
DEPARTAMENTO DE TÉRMICA Y ENERGÉTICA

DINÁMICA DE GASES (Sin editar)

FLUJO ISOTÉRMICO

1.1 Variación de la densidad en la dirección del flujo

Consideremos el flujo compresible de un gas a través de un ducto de sección constante no

aislado y tomemos un volumen de control

Aplicamos las siguientes ecuaciones:

-Ecuación de Cantidad de movimiento

 V  dV   V 
pA  ( p  dp) A  dFw  m
 Adp  dFw  m
 dV

Sabiendo que

dFw  fV 2Ddx

dV
 Adp  fV 2Ddx  AV 2
V

Dividiendo la ecuación por V 2 A

 dp fDdx dV
 
V 2 A V

 dp 4 fdx dV
 
V 2 D V

dp 4 fdx dV
  0 (1)
V 2 D V

-Ecuación de continuidad

dA dV d
  0
A V 
UC. Semestre 1-2011. Dinámicas de Gases. Prof. Gruber A. Caraballo

1
Como A=ctte dA A  0

dV d
 (2)
V 

-Variación de la densidad en la dirección del flujo

De la ecuación de estado p  RT

p  p
 RT  Ctte  d    0
 
1 1
pd    dp  0
 

pd dp
  0
 2

Multiplicando la ecuación por  / p

d dp
  0
 p

d dp
 (3)
 p

Sustituyendo las ecuaciones (2) y (3) en (1) y simplificando se obtiene lo siguiente.

d pd 4 fdx
   0
 V 
2
D

d  p  4 fdx
  1   (4)
  V 2
 D

pV 2
De la ecuación de estado V 2 
RT

Sabiendo que V  Mc  M kRT

V 
p M kRT
2  
2

RT

Multiplicando, dividiendo por k y resolviendo nos queda:

V 2  kpM 2 (5)

Sustituyendo la ecuación (5) en (4) nos queda:

UC. Semestre 1-2011. Dinámicas de Gases. Prof. Gruber A. Caraballo

2
 d  p  4 fdx
  1  
  kpM 2
 D

d   1  kM 2  4 fdx
   
  kM 2
 D

d 4 fkM 2 


dx D kM 2  1  (6)

Tomando la ecuación (6) y analizándola

-Si M  1 / k La densidad disminuye en la dirección del flujo

-Si M  1 / k La densidad aumenta en la dirección del flujo
-Si M  1 / k La densidad tiende a infinito y no puede existir flujo.

1.2 Variación de la entropía en la dirección del flujo

De la segunda ley de la termodinámica

1
Tds  du  pdv  Tds  du  pd  


P
Tds  CvdT  d
2

UC. Semestre 1-2011. Dinámicas de Gases. Prof. Gruber A. Caraballo

3
Como T=ctte entonces dT=0 quedando:

P p
Tds   d  ds   d
 2
 2.T

p d d
ds    R
T  

ds R
 (7)
d 

Si analizamos la ecuación (7) se tiene:

-Si M  1 / k La entropía aumenta en la dirección del flujo

-Si M  1 / k La entropía disminuye en la dirección del flujo
-Si M  1 / k La entropía tiende a infinito y no puede existir flujo.

La ecuación (7) también se puede escribir en función del número de Mach

ds 4 fM  kM 2 
  
 kM  1 
2
dx D

1.3 Variación del número de Mach en la dirección del flujo

Sabiendo que:

V dM 1
M  
c dV c

Como T=ctte entonces c  kRT  ctte

dM dM dV dM dV d
 
dx dx dx dV d dx

dM 1  V 4 f  kM 2 

dx c  D  kM 2  1

dM 4 fM  kM 2 

D  kM 2  1
(8)
dx

Si analizamos la ecuación (8) se tiene:

-Si M  1 / k El Mach aumenta en la dirección del flujo

-Si M  1 / k El Mach disminuye en la dirección del flujo

UC. Semestre 1-2011. Dinámicas de Gases. Prof. Gruber A. Caraballo

4
-Si M  1 / k El Mach tiende a infinito y no puede existir flujo.

Luego de ver como varia las propiedades se puede definir dos tipos de flujo

 

-Si M  1 / k M  Flujo Subcritico
s 


 

-Si M  1 / k M  Flujo Supercrítico
s 


1.4 Longitud Máxima del ducto

Si tomamos la expresión (8) y se escribe de la siguiente manera

dx D  kM 2  1
   (9)
dM 4 fM  kM 2 

Esta ecuación es igual a cero “0” si M  1 / k . Esto significa que la longitud del ducto tiene
un máximo valor posible para unas condiciones bien determinadas a la entrada del mismo
cuando por el pasa un flujo de masa dado. Si integramos la ecuación (9)

4f  kM 2  1 
 dx  
x max 1/ k
 3
dM
x D M
 kM 

4 f ( xmax  x)  1  kM 2 
  2
  Ln kM 2  
D  kM 

4 fLmax  1  kM 2 
  2
  Ln kM 2 (10)  
D  kM 
1.5 Ecuaciones del flujo Isotérmico

Podemos establecer, al igual que en Fanno o en Rayleigh relación entre los parámetro del fluido
en una sección critica de dicho ducto, los cuales se denominan parámetros críticos del fluido

Sabiendo que V  Mc entonces:

V Mc Mc
 
V * * M * *c c / k
UC. Semestre 1-2011. Dinámicas de Gases. Prof. Gruber A. Caraballo

5
 V
M k (11)
V **

-De la ecuación de continuidad

VA   * *V * *A

Como el ducto es de área de sección constante nos queda:

V  **
 M k (12)
V ** 

-De la ecuación de estado p  RT

p RT

p * *  * *RT * *

Por ser el flujo isotérmico eso quiere decir que T=T** entonces la relación anterior queda:

p  1
  (13)
p **  ** M k

-De las relaciones de temperaturas con las de estancamiento:

 k 1 2 
T0  T 1  M  (14)
 2 
Si T  T * * y T0  T0 * * se tiene
 k 1
T0 * *  T * *1  M * *2 
 2 

Sabiendo que M * *  1 / k y sustituyendo en la ecuación anterior nos queda:

 k  1
T0 * *  T * *1 
2k 
(15)


Si dividimos las ecuaciones (14) y (15) nos queda:

 k 1 2 
1 M 
T0 T  2 

T0 * * T * *  k  1
1  2k 

Simplificando:

T0 2k  k  1 2 
 M  1 (16)
T0 * * 3k  1  2 
UC. Semestre 1-2011. Dinámicas de Gases. Prof. Gruber A. Caraballo

6
-De las relaciones de presiones con las de estancamiento:

k
 k  1 2  k 1
p0  p 1  M  (17)
 2 

Si p  p * * y p0  p0 * * se tiene

k
 k 1
M * *2 
k 1
p0 * *  p * *1 
 2 

Sabiendo que M * *  1 / k y sustituyendo en la ecuación anterior nos queda:

k
 k  1 k 1
p0 * *  p * *1 
2k 
(18)


Si dividimos las ecuaciones (14) y (15) nos queda:

k
 k  1 2  k 1
1 M 
p0 p  2 
 k
p0 * * p * *
 k  1 k 1
1  2k 

Sustituyendo la relación (13) y simplificando:

k
p0 1  2k  k  1 2  k 1
  3k  1  2 M  1
p0 * * M k   

UC. Semestre 1-2011. Dinámicas de Gases. Prof. Gruber A. Caraballo

7
APLICACIÓN:

1.-Considere el siguiente arreglo, el flujo es adiabático y con fricción solo en los ductos 1-2, 4-5
de área constante donde f=0,025, En el ducto 5-6 existe transferencia de calor (no adiabático) y
con fricción f=,, hallar:
a. Presión de estancamiento en 2
b. Flujo másico
c. Presión en 1
d. El numero de Mach y la presión en 5
e. El numero de Mach, la presión en 6
f El cambio de entropía de 5-6

Datos:

D1  2,97cm
D3  1,5cm
Dcho  1,7cm
D5  3cm
P3  10kPa
T03  300K

Aire k=1,4
R=287J/(kg.K)

Con los diámetros se calcula el área de cada sección

D 2
A entonces:
2

A1  6,928 *10 4 m 2
A3  1,767 *10 4 m 2
Ach  2,27 *10 4 m 2
A5  7,069 *10 4 m 2

Como existe choque en la parte divergente de la tobera quiere decir que se alcanza Mach igual a
1 en la garganta del mismo (sección 3)

En TFI con M 3  1

UC. Semestre 1-2011. Dinámicas de Gases. Prof. Gruber A. Caraballo

8
 p3
 0,528..........(a)
p03
T3
 0,833...........(b)
T03

De la relación (a) p03  18,939kPa

Como el flujo es isoentrópico antes del choque p03  p02  18,939kPa

Como el flujo alcanza mach igual 1 en la garganta este se encuentra estrangulado por lo tanto el
flujo másico que pasa por la tobera es el flujo máximo

A3 * p03
m max  0,6847
RT03

(1,767 *10 4 m 2 )(18939Pa)

m max  0,6847
(287J / kgK )(300K )
 max  7,809 *103 kg / s
m

Ax Ach
En TFI con   1,285
A * A3

M x  1,64
px
 0,221...........(c)
p0 x
Tx
 0,65.............(d )
T0 X

En TCN con M x  1,64

M y  0,65
py
 2,979...........(e)
px
Ty
 1,417............( f )
TX
p0 y
 0,87..............(g )
p0 x

p03  p0 X  18,939kPa
De la relación (c) p X  4,1855kPa
De la relación (d) TX  195K
UC. Semestre 1-2011. Dinámicas de Gases. Prof. Gruber A. Caraballo

9
De la relación (e) pY  12,47kPa
De la relación (f) TY  276,32K
De la relación (g) p0Y  16,652kPa

De la relación de choque p0Y AY*  p0 X AX*

p0 X AX*
A 
*
Y  2 *10 4 m 2
p0Y

Como A4  A5  A6  7,69 *10 4 m 2

A4
En TFI con  3,98
AY*

M 4  0,166
p4
 0,981.................(h)
p04

En TFF con M 4  0,16

fL max 4
 22,28.............(i)
D
p4
 6,8291......................( j )
p F*

fL max 5 fL max 4 fL45

 
D D D

fL45 4 * 0,025 * 0,5

  1,666
D 3 *102

fL max 5 fL max 4 fL45

   20,614
D D D

fL max5
En TFF con  20,614
D

M 5  0,172

p5
 6,35.................(k )
p *f

De la relación (h) p4  16,33kPa

De la relación (j) p F*  2,39kPa
UC. Semestre 1-2011. Dinámicas de Gases. Prof. Gruber A. Caraballo

10
De la relación (k) p5  15,18kPa

A2 A
En TFI con  1  3,92
A * A3
M 2  0,12
p2
 0,989..........(l )
p02

En TFF con M 2  0,12

fL max 2
 45,408.............(m)
D
p2
 9,1156......................(n)
p F*

fL max1 fL max 2 fL12

 
D D D

fL12 4 * 0,025 *10

  33,67
D 2,97 *10 2

fL max1 fL max 2 fL12

   79,08
D D D

fL max1
En TFF con  79,08
D

M 1  0,09

p1
 12,162.................(o)
p *f

De la relación (l) p2  18,73kPa

De la relación (j) p F*  2,05kPa
De la relación (k) p1  24,99kPa

En TFI M 5  0,172

4 fL max 5
 19,96
D
p5
 4,914...........( p)
p **

UC. Semestre 1-2011. Dinámicas de Gases. Prof. Gruber A. Caraballo

11
De la relación (p) p * *  3,308kPa

fL max 6 fL max 5 fL56

 
D D D
fL56 4 * 0,005 * 5
  1,667
D 3 *10 2

fL max 6 fL max5 fL65

   18,293
D D D

fL max 6
En TFIso  18,293
D

M 6  0,179

p6
 4,722...........(q)
p **

De la relación (q) p6  15,62kPa

Se observa que p disminuye p6  p5 en el sentido de flujo subcrítico

p 
s   RLn 7 
 p6 

 15,62kPa 
s  (287J / kgK ) Ln 
 16,256kPa 

s  11,454J / kgK

UC. Semestre 1-2011. Dinámicas de Gases. Prof. Gruber A. Caraballo

12