9. Flujo de Fluidos compresibles. Cuando un fluido es compresible su densidad cambia a lo largo de la dirección del flujo, producto del cambio en la presión y la temperatura. En general, para un sistema determinado se considera que el fluido es compresible, si su densidad varía más de un 10% entre los puntos de entrada y salida. El flujo de gases a presiones atmosféricas o superiores, puede desarrollarse en régimen laminar o turbulento, con transición para Re = 2000 - 3000.
Flujo isotérmico de un gas ideal a través de una cañería horizontal. Se cumple: p
p1
=
ρ
=
ρ 1
p 2
= cte. =
ρ 2
RT
ρ
M
=
ρ 1
p1
⋅ p
Balance de energía diferencial para condiciones de entrada y salida: p 1, ρ1 y p2, ρ2 respectivamente.
Figura 9.1. Flujo isotérmico de un gas ideal en una cañería horizontal.
⎛ v 2 ⎞ dp ⎟⎟ + + dh f = 0 d ⎜⎜ ⎝ 2 g c ⎠ ρ
(9.1)
⎛ v 2 ⎞ dp dx v 2 ⎟⎟ + + f ⋅ ⋅ d ⎜⎜ =0 ρ g D g 2 2 c ⎝ c ⎠ 2
dv v
+
2 g c dp v2
⋅
+ f ⋅
ρ
/⋅
dx D
2gc
(9.2)
v2
=0
(9.3)
Por continuidad. G = v⋅ρ = cte. dG = v⋅dρ + ρ⋅dv = 0 dv v
=−
d ρ
,
ρ
d ρ =
ρ 1
p1
⋅ dp
MECÁNICA DE FLUIDOS
ρ
=
ρ 1
p1
⋅ p
d ρ ρ
=
dp
dv
p
v
=−
dp p
DPQBA, UTFSM
2 g c dp v2
⋅
,
ρ
G
= v ⋅ ρ
G
2
2
2
= v ⋅ ρ
v
2
=
G2
2g c
2 ρ
G
2
2 g c dp 2 g c ρ 1 ⋅ = 2 pdp 2
ρ dp
v
ρ
G p1
Reemplazando en balance de energía:
−2
dp p
+
2 g c ρ 1 G 2 p1
pdp + f ⋅
dx D
=0
(9.4)
Integrando p: entre p 1 y p2 y x: entre x = 0 y x = L: p 2
− 2∫
dp
2 g c ρ 1
p 2
1 L
∫ pdp + f D ∫ dx = 0
(9.5)
⎛ p ⎞ g ρ L − 2 ln⎜⎜ 2 ⎟⎟ + c2 1 [ p 22 − p12 ] + f = 0 D ⎝ p1 ⎠ G p1
(9.6)
p1
p
+
G 2 p1
0
p1
Se define N = número de alturas de velocidad, N = = f
L D
⎛ p ⎞ g M − 2 ln⎜⎜ 2 ⎟⎟ + 2c [ p22 − p12 ] + N = 0 ⎝ p1 ⎠ G RT
(9.7)
2
⎛ p ⎞ g M N − ln ⎜⎜ 2 ⎟⎟ = 2c [ p12 − p 22 ] ⎝ p1 ⎠ G RT
(9.8)
Ecuación general para el flujo másico superficial de un gas ideal isotérmico: G
2
=
g c M RT
[ p12 − p22 ] 2 ⎡ ⎛ p2 ⎞ ⎤ ⎢ N − ln⎜⎜ ⎟⎟ ⎥ ⎢⎣ ⎝ p1 ⎠ ⎥⎦
(9.9)
Ecuación particular si se cumple:
⎛ p ⎞ N >> ln ⎜⎜ 2 ⎟⎟ ⎝ p1 ⎠
2
MECÁNICA DE FLUIDOS
DPQBA, UTFSM
Ecuación de Weymouth: G
2
=
g c M RTN
[ p12 − p22 ]
(9.10)
Según la ecuación general (9.9), en un comienzo si p 2 disminuye G2 aumenta. Se cumple que indistintamente para p 2 →p1, G2→0, y para p 2→0, G2→0. Lo anterior indica la existencia de un valor máximo para G 2 asociado con un valor de p2 crítico p2*.
Figura 9.2. Flujo másico superficial en función de la presión de salida p 2. La región entre 0 < p < p 2* es ficticia, ya que una vez que la presión de salida p 2 alcanza el valor crítico p2* no se obtiene un flujo mayor. La presión p 2* crítica se obtiene derivando la ecuación general con respecto a p 2 para un p1 fijo según: dG
2
dp 2 2
=0
(9.11)
2
⎛ p1 ⎞ ⎛ p ⎞ ⎜⎜ ∗ ⎟⎟ − ln⎜⎜ 1* ⎟⎟ = N + 1 ⎝ p 2 ⎠ ⎝ p 2 ⎠
(9.12)
Resolviendo para p2*. Gmáx2 correspondiente obedece a la ecuación: 2 Gmáx
=
2
Gmáx
MECÁNICA DE FLUIDOS
g c M RT
=
( p2* )2 = g ρ 1 ( p2* )2 c
p1
g c ρ 1 p1 g ρ p 1 + N + ln c 21 1 Gmáx
(9.13)
(9.14)
DPQBA, UTFSM
También se puede demostrar que: G máx
*
=
*
g c p 2 ρ 2
= ρ 2*v 2*
(9.15)
La correspondiente velocidad de salida será: *
v 2*
g c p 2
=
(9.16)
*
ρ 2
La cual se puede interpretar como la velocidad de una hipotética onda de sonido isotérmica a las condiciones de salida, dado que se dispone de las siguientes relaciones para la velocidad del sonido y un gas ideal:
=
c
g c dp
p
d ρ
ρ
=
RT M
Así, la velocidad de una hipotética onda de sonido isotérmica está dada por: c=
g c RT
(9.17)
M
En la práctica, sin embargo, las ondas circulan en forma isentrópica, y la velocidad del sonido es: c
Donde, k =
=
kg c RT
(9.18)
M
c p cv
Flujo adiabático de un gas ideal a través de una cañería horizontal. Se cumple: p k
ρ
=
p1 k
ρ 1
=
p 2 k
ρ 2
= cte.
ρ
=
ρ 1 k p11 /
⋅ p1 / k
k =
c p cv
Balance de energía diferencial para condiciones de entrada y salida: p 1, ρ1, T 1 y p2, ρ2, T2 respectivamente.
⎛ v 2 ⎞ dp ⎟⎟ + + dh f = 0 d ⎜⎜ g 2 ⎝ c ⎠ ρ
MECÁNICA DE FLUIDOS
(9.1)
DPQBA, UTFSM
⎛ v 2 ⎞ dp dx v 2 ⎟⎟ + + f ⋅ ⋅ d ⎜⎜ =0 ρ g D g 2 2 c ⎝ c ⎠ dv
2
v
2 g c dp
+
⋅
v2
+ f ⋅
ρ
/⋅
dx
2gc
(9.2)
v2
=0
D
(9.3)
Por continuidad. G = v⋅ρ = cte. dG = v⋅dρ + ρ⋅dv = 0 dv v
=−
d ρ
,
ρ
d ρ =
2 g c dp v
2
⋅
1
1− k
ρ 1
k
p 1 / k
k p1
dp
ρ
2
2
ρ 1 k p11 /
=
d ρ
⋅ p1 / k
=
ρ
1 dp
dv
k p
v
=−
1 dp k p
,
ρ
G
= v ⋅ ρ
G
2
2
= v ⋅ ρ
v
G
=
2g c
2 2
G
ρ
2
ρ dp
2 g c dp 2 g c ρ 1 1 / k ⋅ = 2 1 / k p dp 2 v
ρ
G p1
Reemplazando en balance de energía:
−
2 dp k p
+
2 g c ρ 1
p 1 / dp + f ⋅ 1 / k k
G 2 p1
dx D
=0
(9.19)
Integrando p: entre p 1 y p2 y x: entre x = 0 y x = L: p
−
2 2 dp
∫
k p1 p
+
2 g c ρ 1 2
p 2
∫ p
1 / k
G p1
1 / k
dp + f
p1
1 L D
∫ dx = 0
(9.20)
0
k +1 ⎤ 2 ⎛ p 2 ⎞ 2 g c ρ 1 k ⎡ k k +1 L k p p − ln⎜⎜ ⎟⎟ + 2 1 / k − ⎢ 2 ⎥ + f = 0 1 k ⎝ p1 ⎠ G p1 k + 1 ⎣ D ⎦
Se define N = número de alturas de velocidad, N = f densidad las condiciones de entrada, ρ 1 = k −1
2 ⎛ p2 ⎞ 2 g c Mp1 k − ln⎜⎜ ⎟⎟ + 2 k
⎝ p1 ⎠
⎛ p ⎞ N − ln⎜⎜ 2 ⎟⎟ ⎝ p1 ⎠ MECÁNICA DE FLUIDOS
2 g c Mp1 k 2
, y se utiliza para la
.
k +1 ⎡ k k +1 ⎤ ⎢ p2 − p1 k ⎥ + N = 0 k + 1 ⎣ ⎦
k −1
=
RT 1
D
k
G RT 1
2 / k
Mp1
L
G RT 1
k +1 ⎡ k k +1 ⎤ k p p − ⎢ 1 2 ⎥ k + 1 ⎣ ⎦
k
(9.21)
(9.22)
(9.23)
DPQBA, UTFSM
Ecuación general para el flujo másico superficial de un gas ideal adiabático equivalente a (9.9): k +1 ⎡ k k +1 ⎤ ⎢ p1 − p 2 k ⎥ k ⎣ ⎦ 2 / k k + 1 ⎡ ⎛ p ⎞ ⎤ ⎢ N − ln⎜⎜ 2 ⎟⎟ ⎥ ⎢⎣ ⎝ p1 ⎠ ⎥⎦
k −1
G
2
=
2 g c Mp1 k RT 1
(9.24)
Ecuación particular si se cumple:
⎛ p ⎞ N >> ln ⎜⎜ 2 ⎟⎟ ⎝ p1 ⎠
2
Ecuación de Weymouth, equivalente a (9.10): k −1
G
2
=
k +1 ⎡ k k +1 ⎤ ⎢ p1 − p 2 k ⎥ k + 1 ⎣ ⎦
2 g c Mp1 k
k
NRT 1
(9.25)
Al igual que con flujo isotérmico existe un G máx asociado con una presión de salida p2 crítica p2*. La presión p 2* crítica se obtiene derivando la ecuación general con respecto a p 2 para un p 1 fijo según: dG
2
dp 2 p ⎞ ⎛ 2 ⎞⎛ ⎜ ⎟⎜⎜ 1∗ ⎟⎟ ⎝ k + 1 ⎠⎝ p 2 ⎠
k +1 k
=0
⎛ p ⎞ − ln⎜⎜ 1* ⎟⎟ ⎝ p2 ⎠
(9.11)
2 / k
2 ⎞ = N + ⎛ ⎜ ⎟ ⎝ k + 1 ⎠
(9.26)
Resolviendo para p2*. Gmáx2 correspondiente obedece a la ecuación: 2
G máx
2 Gmáx
=
=
g c kM RT 1
k −1
p1
k
k +1
( ) p 2*
k
=
g c k ρ 1 k p11 /
k +1
( ) p 2*
k
(9.27)
2 g c k ρ 1 p1 2 ⎡ ⎤ k +1 ⎛ ⎞ ρ g k p 2 ⎞ + N + ln⎜ c 1 1 ⎟ ⎥ (k + 1) ⋅ ⎢⎢⎛ ⎜ ⎟ ⎜ G2 ⎟ ⎥ ⎝ k + 1 ⎠ ⎝ máx ⎠ ⎥ ⎢⎣ ⎦
(9.28)
Todas las ecuaciones para flujo adiabático se transforman en las de flujo isotérmico para k = 1.
MECÁNICA DE FLUIDOS
DPQBA, UTFSM
FLUJO COMPRESIBLE EN TOBERAS. Otro ejemplo de flujo compresible ocurre con la descarga de un gas desde un recipiente a alta presión a través de una tobera, consistente de una sección convergente que conduce a una “garganta”, seguida posiblemente de una sección divergente o “difusor”.
Figura 9.3. Flujo a través de una tobera convergente/divergente. El flujo de alta velocidad descargando a la atmósfera es del tipo adiabático y dado que sólo se desarrolla en una corta distancia se le considera sin fricción. p k
p1
=c=
ρ
⎛ p ⎞ ρ = ρ 1 ⎜ ⎜ p ⎟⎟ ⎝ 1 ⎠
k
ρ 1
1 / k
k =
c p cv
Balance de energía entre el recipiente y alguna posición en la tobera donde la velocidad es v y la presión es p. v
2
2g c p
∫
p1
dp ρ
=
p11 /
k p
ρ 1
dp
∫ p
p1
−
1 / k
v12
p
+∫
dp
=0
(9.29)
⎡⎛ p ⎞ ( k −1) / k ⎤ ⎢⎜ ⎟ = − 1⎥ k − 1 ρ 1 ⎢⎜⎝ p1 ⎠⎟ ⎥⎦ ⎣
(9.30)
2g c
p1
ρ
k p1
Reemplazando en el balance y considerando v 1 = 0. ( k −1) / k ⎡ ⎤ ⎛ ⎞ g k 2 p p c 1 ⎢1 − ⎜⎜ ⎟⎟ ⎥ v2 = k − 1 ρ 1 ⎢ ⎝ p1 ⎠ ⎥⎦ ⎣
MECÁNICA DE FLUIDOS
(9.31)
DPQBA, UTFSM
Figura 9.4. Efecto de la variación de la presión de salida en el flujo de la tobera. Dado que w = ρ⋅v⋅A en cualquier punto, donde w es el flujo másico de gas que descarga y A es el área de sección transversal en la tobera. G2
2 ⎡ ⎛ p ⎞ ( k −1) / k ⎤⎛ p ⎞ 2 / k 2 g c k w ⎞ ⎛ ⎥⎜⎜ ⎟⎟ p1 ρ 1 ⎢1 − ⎜⎜ ⎟⎟ =⎜ ⎟ = A k p − 1 ⎝ ⎠ ⎢⎣ ⎝ 1 ⎠ ⎥⎦⎝ p1 ⎠
(9.32)
El flujo másico superficial w/A (G) es claramente un máximo en la garganta, donde vale w/AT (o GT). Dado que la presión en la garganta es una variable, existe un valor máximo para w/AT (o GT), que se obtiene de: dGT d ( p / p1 )
=
d (w / AT ) d ( p / p1 )
=0
(9.33)
Así la razón de presión crítica es: p cT p1
2 ⎞ = ⎛ ⎜ ⎟ ⎝ k + 1 ⎠
k /( k −1)
(9.34)
La velocidad vcT correspondiente es: 2
vcT
⎛ 2 g k ⎞ p kg p = ⎜ c ⎟ 1 = c cT = ρ cT ⎝ k + 1 ⎠ ρ 1
MECÁNICA DE FLUIDOS
kg c RT cT M
(9.35)
DPQBA, UTFSM
Finalmente el flujo másico crítico wc es:
wc
= AT
⎛ 2 ⎞ kg c p1 ρ 1 ⎜ ⎟ ⎝ k + 1 ⎠
( k +1) /( k −1)
(9.36)
En relación a la figura se tienen varias situaciones según la presión de salida p 2. p2 A B C
D
E
F
Fenómenos. Si la presión de salida es ligeramente inferior a la del recipiente, existe un flujo de descarga pequeño. Se puede evaluar utilizando la ecuación (9.32) sustituyendo A = A 2 (área de salida) y p =p 2. La ecuación (9.32) entrega la variación de presión en la tobera. El flujo es siempre subsónico. Lo mismo que en A, excepto que la velocidad de flujo es más grande. Si la presión de salida se reduce suficientemente, la velocidad en la garganta alcanza el valor dado por la ecuación (9.35). En la sección divergente, la presión crece y el flujo es subsónico. Para una presión de salida entre C y E, no hay solución continua posible. El flujo, que es crítico, es supersónico por una cierta distancia más allá de la garganta, sin embargo existe un incremento repentino de la presión, conocido como “shock”, y a continuación de esto el flujo es subsónico. El “shock” es un fenómeno irreversible, que resulta en bruscos cambios en velocidad, presión y temperatura sobre una extremadamente corta distancia de unas pocas moléculas en espesor. Para el mismo flujo másico crítico como en C, la ecuación (9.32) tiene una segunda raíz, correspondiente a una presión de salida E. En este caso, sin embargo, existe un continuo decrecimiento de la presión en el difusor, donde el flujo es ahora supersónico. Para una presión de salida más baja que E, una posterior expansión irreversible ocurre justo afuera de la tobera.
Si no hay difusor, el flujo es esencialmente aquel que se obtiene a través de un orificio de un estanque a presión. El flujo será subsónico si la presión de salida excede pcT. Si la presión de salida es igual a p cT, entonces se obtiene un flujo crítico con velocidad sónica a través del orificio. Y si cae bajo este valor, aun se obtendrá flujo crítico, pero con una posterior expansión irreversible justo afuera del orificio.
MECÁNICA DE FLUIDOS
DPQBA, UTFSM
SOLUCIÓN GRÁFICA PARA FLUJO COMPRESIBLE Considere la descarga de un estanque de grandes dimensiones, como la situación presentada en la figura 9.5:
Figura 9.5. Descarga de un estanque a través de una tubería Se define el parámetro G ci, como el flujo másico máximo por unidad de área hipotético que se alcanza con una expansión isotérmica del gas a través del sistema mostrado en la Figura 9.5, cuando N = 0, es decir en el punto 1 o de descarga. Gci
=
g C ⋅ p 0
g C ⋅ M ⋅ ρ 0 = p 0 2,718 2,718 ⋅ R ⋅ T 0
⎡ lbm ⎤ ⎢ ft 2 ⋅ s ⎥ ⎣ ⎦
(9.37)
Gráficamente a partir de ecuaciones de balances de energía y de continuidad Lapple determinó para diferentes tipos de flujo (k = 1.0, 1.4 y 1.8), la relación que existe entre la razón G/G ci (flujo másico superficial G en la cañería y el parámetro Gci) con la razón de presiones (p 2/p0) o (p3/p0), para diversos valores de N, número de alturas de velocidad, asociados con la resistencia por fricción en la cañería. • G, Gci, p0, p1, p2 y p3 definidos según esquema de descarga. • Otros flujos según valores de k se pueden obtener por interpolación. Los gráficos de Lapple también identifican un flujo másico superficial máximo Gmáx, obtenido de la razón (G/Gci)c crítica o máxima, de acuerdo a las presiones que se desarrollen a lo largo de la cañería. El valor aludido, para un N cualquiera, se obtiene cuando la curva en el gráfico se hace vertical. Las razones (p2/p0) y (p3/p0) son iguales, para un N cualquiera, si la razón G/Gci es menor que el valor crítico o máximo (G/G ci)c. Sólo cuando se alcanza la razón (G/Gci)c, para un N cualquiera, (p 3/p0) puede ser menor que el valor de (p 2/p0). En ese caso el flujo es el máximo e independiente del valor de (p 3/p0). La razón (p2/p0) se conserva fija, ya que corresponde a la presión p 2* crítica definida antes.
MECÁNICA DE FLUIDOS
DPQBA, UTFSM
Figura 9.6. Gráfico de Lapple para flujo isotérmico k =1.0
Figura 9.7. Gráfico de Lapple para flujo adiabático k =1.4
MECÁNICA DE FLUIDOS
DPQBA, UTFSM
Figura 9.8. Gráfico de Lapple para flujo adiabático k =1.8 Una situación para el cálculo de flujos compresibles se puede plantear conociendo las condiciones al interior de la cañería. Suponga conocidas las condiciones de presión y temperatura en 1 y el flujo másico superficial G que circula por la cañería. 1
2
Para resolver este sistema de modo de conocer la presión y temperatura en un punto 2 cualquiera. Se puede plantear un arreglo ficticio formado por un estanque a presión y temperatura p 0 y T0, para las condiciones hipotéticas de entrada dadas por el punto 1.
Entonces suponiendo un estanque en condiciones “0” desconocidas:
MECÁNICA DE FLUIDOS
DPQBA, UTFSM
Se define G’ci: Gci'
g c · M
= p1 ·
(9.38)
e· R·T 1
Se puede escribir la razón entre las ecuaciones (9.37) y (9.38): Gci Gci'
G GCI
=
p 2 p1
T 2 T 1
p 0
=
G ' GCI
p1 p1
·
p1
T 2
=
T 1
T 1 T 0
p 0
T 1
p 2
=
·
(9.40)
T 0
p 0
(9.41)
p 0
T 0
(9.42)
T 0
Método de cálculo:
Se determina G’ci con (9.38) Se supone G/Gci Se determina para N=0: p1 /p0 y T 1 /T 0 Se calcula G/Gci de (9.40) y se chequea. Si cumple, entonces se determina para N=N : p2 /p0 y T 2 /T 0 Si no cumple se itera con G/Gci
MECÁNICA DE FLUIDOS
DPQBA, UTFSM
PROBLEMAS RESUELTOS
1. Gas natural (metano, se asume como gas ideal) circula estacionariamente a 55 ºF por una cañería horizontal de 12” de diámetro nominal de 20 millas de longitud, con un factor de fricción f = 0.014. Si la presión de entrada es de 100 [psia], que presión de salida correspondería con la velocidad de flujo máximo en la cañería? Si la presión de salida real es de 10 [psia], cual es la velocidad de flujo del gas en [lb m/h]? Repetir para flujo adiabático considerando T1 = 55 ºF, (k = 1.31).
ISOTÉRMICO Datos generales: T [ºR] f
515 0.014
M N
L [millas]
20
p1 [psia]
L [ft] D [ft]
105600 1
16 1478.4 2
p1 [lbf /ft ]
100 14400
Determinación de la presión crítica de salida p 2*: 2
2
⎛ p1 ⎞ ⎛ p ⎞ ⎜⎜ ∗ ⎟⎟ − ln⎜⎜ 1* ⎟⎟ = N + 1 ⎝ p2 ⎠ ⎝ p2 ⎠ Tabla para el cálculo numérico de p2*: p1 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400
2
p2 380.00 379.00 378.00 377.00 376.00 375.00 374.00 373.90 373.80 373.70 373.60 373.50 373.49 373.48
(p1/p2) 1436.01 1443.60 1451.25 1458.96 1466.73 1474.56 1482.46 1483.25 1484.04 1484.84 1485.63 1486.43 1486.51 1486.59 1486.67 1486.75 1486.83
373.47
373.46 373.45
2ln(p1/p2) 7.3 7.3 7.3 7.3 7.3 7.3 7.3 7.3 7.3 7.3 7.3 7.3 7.3 7.3 7.3 7.3 7.3
diferencia 1428.74 1436.32 1443.97 1451.67 1459.44 1467.26 1475.15 1475.95 1476.74 1477.53 1478.33 1479.12 1479.20 1479.28 1479.36 1479.44 1479.52
N+1 1479.40 1479.40 1479.40 1479.40 1479.40 1479.40 1479.40 1479.40 1479.40 1479.40 1479.40 1479.40 1479.40 1479.40 1479.40 1479.40 1479.40
DELTA 50.66 43.08 35.43 27.73 19.96 12.14 4.25 3.45 2.66 1.87 1.07 0.28 0.20 0.12 0.04 -0.04 -0.12
Resultado: *
p2
2
373.47 [lbf /ft ]
2.59 [psia]
MECÁNICA DE FLUIDOS
DPQBA, UTFSM
Determinación del flujo máximo G2máx: 2 Gmáx
=
g c ρ 1 p1 g ρ p 1 + N + ln c 21 1 Gmáx
Tabla para el cálculo numérico de G2máx: 2
Gmáx 90 90.1 90.2 90.3
90.31
90.32 90.33 90.34 90.35 90.36 90.37 90.38 90.39 90.4 90.41 90.42 90.43
gcρ1p1 N+1 134265.714 1479.40 134265.714 1479.40 134265.714 1479.40 134265.714 1479.40 134265.714 1479.40 134265.714 1479.40 134265.714 1479.40 134265.714 1479.40 134265.714 1479.40 134265.714 1479.40 134265.714 1479.40 134265.714 1479.40 134265.714 1479.40 134265.714 1479.40 134265.714 1479.40 134265.714 1479.40 134265.714 1479.40
2
LN(gcρ1p1/Gmáx ) 7.307766386 7.306655891 7.305546629 7.304438595 7.30432786 7.304217136 7.304106425 7.303995726 7.303885039 7.303774364 7.303663702 7.303553052 7.303442414 7.303331789 7.303221175 7.303110574 7.302999985
DELTA -0.31 -0.21 -0.11 -0.01 0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10 0.11 0.12
Gmáx 9.49 9.49 9.50 9.50 9.50
9.50 9.50 9.50 9.51 9.51 9.51 9.51 9.51 9.51 9.51 9.51 9.51
Resultado: 2
Gmáx
9.50 [lbm/ft s]
Determinación de la velocidad de flujo del gas en [lb m/h] para una presión de salida p2 de 10 [psia]: G
2
=
[ p12 − p22 ]
g c M RT
p1 14400
2 ⎡ ⎛ p2 ⎞ ⎤ ⎢ N − ln⎜⎜ ⎟⎟ ⎥ ⎢⎣ ⎝ p1 ⎠ ⎥⎦
p2 1440
gcM/RT 0.000647501
2
2
(p1) - (p2) 2ln(p2/p1) 205286400.0 -4.6
N 1478.4
2
G 89.63
G 9.47 2
[lbm/ft s]
Resultado: 2
Aflujo
0.785
[ft ]
W
7.44
[lbm/s]
W
26768
[lbm/h]
MECÁNICA DE FLUIDOS
DPQBA, UTFSM
ADIABÁTICO Datos generales: T1 [ºR] f
515 0.014
M N
L [millas]
20
p1 [psia]
L [ft] D [ft]
105600 1
16 1478.4 100
2
p1 [lbf /ft ] k
14400 1.31
Determinación de la presión crítica de salida p 2*: p ⎞ ⎛ 2 ⎞⎛ ⎜ ⎟⎜⎜ 1∗ ⎟⎟ ⎝ k + 1 ⎠⎝ p2 ⎠
k +1 k
⎛ p ⎞ − ln⎜⎜ 1* ⎟⎟ ⎝ p2 ⎠
2 / k
2 ⎞ = N + ⎛ ⎜ ⎟ ⎝ k + 1 ⎠
Tabla para el cálculo numérico de p2*: p1 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400 14400
(k+1)/k
p2 220.00 219.00 218.00 217.00 216.00 215.00 214.00 213.00 212.00 211.00
(2/(k+1))(p1/p2) 1379.02 1390.15 1401.41 1412.82 1424.37 1436.08 1447.93 1459.94 1472.10 1484.43 1485.67 1486.91 1488.16 1489.40 1490.65 1491.90 1493.15 1494.40 1495.66
210.90
210.80 210.70 210.60 210.50 210.40 210.30 210.20 210.10
(2/k)ln(p1/p2) 6.38 6.39 6.40 6.40 6.41 6.42 6.43 6.43 6.44 6.45 6.45 6.45 6.45 6.45 6.45 6.45 6.45 6.45 6.45
diferencia 1372.64 1383.76 1395.01 1406.41 1417.96 1429.66 1441.50 1453.50 1465.66 1477.98 1479.22 1480.46 1481.71 1482.95 1484.20 1485.45 1486.70 1487.95 1489.21
N+(2/(k+1)) 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27
DELTA 106.63 95.51 84.25 72.85 61.31 49.61 37.76 25.76 13.60 1.29 0.04 -1.20 -2.44 -3.69 -4.93 -6.18 -7.43 -8.69 -9.94
Resultado: *
p2
2
210.90 [lbf /ft ]
1.46 [psia]
Determinación del flujo máximo G2máx:
2 Gmáx
=
⎛ k ⎞ 2 g c ρ 1 p1 ⎜ ⎟ ⎝ k + 1 ⎠ g k ρ p ⎞ ⎛ 2 ⎞ + N + ln⎛ ⎜⎜ c 2 1 1 ⎟⎟ ⎜ ⎟ ⎝ k + 1 ⎠ ⎝ Gmáx ⎠
MECÁNICA DE FLUIDOS
2 k +1
DPQBA, UTFSM
Tabla para el cálculo numérico de G2máx: 2
Gmáx 100.00 101.00 102.00 102.10 102.20 102.30 102.40 102.41 102.42 102.43 102.44 102.45 102.46 102.47 102.48 102.49 102.50 102.51 102.52
2gcρ1p1(k/(k+1)) 152284.0566 152284.0566 152284.0566 152284.0566 152284.0566 152284.0566 152284.0566 152284.0566 152284.0566 152284.0566 152284.0566 152284.0566 152284.0566 152284.0566 152284.0566 152284.0566 152284.0566 152284.0566 152284.0566
N+(2/(k+1)) 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27 1479.27
2
(2/(k+1))ln(gckρ1p1/Gmáx ) 6.469638967 6.461023962 6.452493835 6.451645427 6.450797849 6.4499511 6.449105178 6.449020631 6.448936093 6.448851562 6.44876704 6.448682527 6.448598021 6.448513524 6.448429035 6.448344554 6.448260082 6.448175618 6.448091161
DELTA -2.49742484 -1.49801917 -0.49860765 -0.39866618 -0.29872466 -0.19878307 -0.09884143 -0.08884727 -0.0788531 -0.06885893 -0.05886476 -0.04887059 -0.03887642 -0.02888225 -0.01888808 -0.00889391 0.00110026 0.01109444 0.02108861
Gmáx 10.00 10.05 10.10 10.10 10.11 10.11 10.12 10.12 10.12 10.12 10.12 10.12 10.12 10.12 10.12 10.12 10.12
10.12 10.13
Resultado: Gmáx
2
10.12 [lbm/ft s]
Determinación de la velocidad de flujo del gas en [lb m/h] para una presión de salida p2 de 10 [psia]: k −1
G
2
=
2 g c Mp1 k RT 1
k +1 ⎡ k k +1 ⎤ k p p − ⎢ 1 2 ⎥ k ⎣ ⎦ 2 / k k + 1 ⎡ ⎛ p ⎞ ⎤ ⎢ N − ln⎜⎜ 2 ⎟⎟ ⎥ ⎢⎣ ⎝ p1 ⎠ ⎥⎦
p1
p2
2gcM/RT
14400
1440
0.001295001
(k-1)/k
(k/(k+1)p1
5.466373405
(k+1)/k
(p1)
(k+1)/k
- (p2)
21141249.6
(2/k)ln(p2/p1) -3.5
N
G
1478.4
2
100.99
G 10.05 2
[lbm/ft s]
Resultado: 2
Aflujo
0.785
[ft ]
W
7.89
[lbm/s]
W
28414
[lbm/h]
MECÁNICA DE FLUIDOS
DPQBA, UTFSM
2. Etileno debe ser bombeado a lo largo de una cañería de 6” de diámetro interno por una distancia de 5 millas a una velocidad de flujo másico de 2.0 [lbm/s]. La presión de suministro al final de la cañería debe ser de 2.0 atmósferas absolutas, y el flujo se puede considerar isotérmico a 60ºF. Si el factor de fricción f = 0.012, calcular la presión de entrada requerida. Asumir comportamiento de gas ideal, y justificar cualquier suposición adicional. Repetir para flujo adiabático considerando T1 = 60ºF, (k = 1.255)
ISOTÉRMICO Datos generales: T [ºR]
520
M
28
W [lbm/s]
f
0.012
N
633.6
Aflujo [ft ]
L [millas]
5
p2 [atm]
2.0
G [lbm/ft s]
L [ft]
26400
p2 [psia]
29.4
D [ft]
2
0.5
p2 [lbf /ft ]
2.0
2
0.196
2
10.19
4233.6
Determinación de la presión de entrada p 1 con ecuación general: G
2
=
[ p12 − p22 ]
g c M RT
2 ⎡ ⎛ p2 ⎞ ⎤ ⎢ N − ln⎜⎜ ⎟⎟ ⎥ ⎢⎣ ⎝ p1 ⎠ ⎥⎦
Tabla para el cálculo numérico de p1: 2
2
G
2 calculado
2
p1
p2
gcM/RT
(p1) - (p2)
2ln(p2/p1)
N
G
G
DELTA
8788.0
4233.6
0.001122231
59305575.0
-1.5
633.6
104.80
103.75
10.19
1.05
8768.0
4233.6
0.001122231
58954455.0
-1.5
633.6
104.18
103.75
10.19
0.43
8758.0
4233.6
0.001122231
58779195.0
-1.5
633.6
103.87
103.75
10.19
0.12
8757.0
4233.6
0.001122231
58761680.0
-1.5
633.6
103.84
103.75
10.19
0.09
8756.0
4233.6
0.001122231
58744167.0
-1.5
633.6
103.81
103.75
10.19
0.06
8755.0
4233.6
0.001122231
58726656.0
-1.5
633.6
103.78
103.75
10.19
0.03
8754.9
4233.6
0.001122231
58724905.1
-1.5
633.6
103.78
103.75
10.19
0.02
8754.8
4233.6
0.001122231
58723154.1
-1.5
633.6
103.77
103.75
10.19
0.02
8754.7
4233.6
0.001122231
58721403.1
-1.5
633.6
103.77
103.75
10.19
0.02
8754.6
4233.6
0.001122231
58719652.2
-1.5
633.6
103.77
103.75
10.19
0.01
8754.5
4233.6
0.001122231
58717901.3
-1.5
633.6
103.76
103.75
10.19
0.01
8754.4
4233.6
0.001122231
58716150.4
-1.5
633.6
103.76
103.75
10.19
0.01
8754.3
4233.6
0.001122231
58714399.5
-1.5
633.6
103.76
103.75
10.19
0.00
8754.2
4233.6
0.001122231
58712648.7
-1.5
633.6
103.75
103.75
10.19
0.00
8754.1
4233.6
0.001122231
58710897.8
-1.5
633.6
103.75
103.75
10.19
0.00
8754.0
4233.6
0.001122231
58709147.0
-1.5
633.6
103.75
103.75
10.19
-0.01
Resultado: p1
2
8754.2 [lbf /ft ]
60.79 [psia]
MECÁNICA DE FLUIDOS
4.14 [atm]
DPQBA, UTFSM
Determinación de la presión de entrada p 1 con ecuación de Weymouth: G2
=
g c M RTN
[ p12 − p22 ] 2
2
gcM/NRT G 1.7712E-06 103.75
2
p2 p1 p1 17923368.96 76501202.3 8746.5
Resultado: 2
p1
8746.5 [lbf /ft ]
60.74 [psia]
4.13 [atm]
ADIABÁTICO Datos generales: T1 [ºR]
520
f
M
0.012
28
N
W [lbm/s]
633.6
2.0
2
0.196
2
Aflujo [ft ]
L [millas]
5
p2 [atm]
2.0
G [lbm/ft s]
10.19
L [ft]
26400
p2 [psia]
29.4
k
1.255
D [ft]
0.5
2
p2 [lbf /ft ]
4233.6
Determinación de la presión de entrada p 1 con ecuación general: k −1
G2
=
2 g c Mp1 k RT 1
k +1 ⎡ k k +1 ⎤ ⎢ p1 − p2 k ⎥ k ⎣ ⎦ 2 / k k + 1 ⎡ ⎛ p2 ⎞ ⎤ ⎢ N − ln⎜⎜ ⎟⎟ ⎥ ⎢⎣ ⎝ p1 ⎠ ⎥⎦
Tabla para el cálculo numérico de p1: p1
p2
2gcM/RT1
(k-1)/k
(k/(k+1)p1
(k+1)/k
(p1)
(k+1)/k
- (p2)
(2/k)ln(p2/p1)
N
G
2 calculado
2
G
G
DELTA
8600.0 4233.6
0.002244461
3.507010805
8452124.6
-1.1
633.6
104.8157244 103.75 10.19
1.06
8590.0 4233.6
0.002244461
3.50618184
8427613.7
-1.1
633.6
104.4873629 103.75 10.19
0.73
8580.0 4233.6
0.002244461
3.505352105
8403125.5
-1.1
633.6
104.1594023 103.75 10.19
0.41
8570.0 4233.6
0.002244461
3.5045216
8378660.0
-1.1
633.6
103.8318427 103.75 10.19
0.08
8569.0 4233.6
0.002244461
3.504438507
8376214.7
-1.1
633.6
103.7991088 103.75 10.19
0.05
8568.0 4233.6
0.002244461
3.504355406
8373769.6
-1.1
633.6
103.7663789 103.75 10.19
0.01
8567.9 4233.6
0.002244461
3.504347095
8373525.1
-1.1
633.6
103.7631061 103.75 10.19
0.01
8567.8 4233.6
0.002244461
3.504338785
8373280.6
-1.1
633.6
103.7598334 103.75 10.19
0.01
8567.7 4233.6
0.002244461
3.504330474
8373036.1
-1.1
633.6
103.7565607 103.75 10.19
0.00
4233.6
0.002244461
3.504322163
8372791.6
-1.1
633.6
103.753288 103.75 10.19
0.00
8567.5 4233.6
0.002244461
3.504313852
8372547.1
-1.1
633.6
103.7500154 103.75 10.19
0.00
8567.4 4233.6 0.002244461
3.504305542
8372302.7
-1.1
633.6
103.7467429 103.75 10.19
-0.01
8567.3 4233.6 0.002244461
3.504297231
8372058.2
-1.1
633.6
103.7434703 103.75 10.19
-0.01
8567.2 4233.6 0.002244461
3.50428892
8371813.7
-1.1
633.6
103.7401979 103.75 10.19
-0.01
8567.1 4233.6 0.002244461
3.504280608
8371569.2
-1.1
633.6
103.7369254 103.75 10.19
-0.02
8567.0 4233.6
3.504272297
8371324.8
-1.1
633.6
103.733653 103.75 10.19
-0.02
8567.6
0.002244461
MECÁNICA DE FLUIDOS
DPQBA, UTFSM
Resultado: 2
p1
8567.6 [lbf /ft ]
59.50 [psia]
4.05 [atm]
Determinación de la presión de entrada p 1 con ecuación de Weymouth: k −1
G
2
=
2 g c Mp1 k NRT 1
k +1 ⎡ k k +1 ⎤ k − p p ⎢ 1 2 ⎥ k + 1 ⎣ ⎦
k
Tabla para el cálculo numérico de p1: p1 8563.0 8562.9 8562.8 8562.7 8562.6 8562.5 8562.4 8562.3 8562.2 8562.1 8562.0
8561.9 8561.8
2
G 103.75 103.75 103.75 103.75 103.75 103.75 103.75 103.75 103.75 103.75 103.75 103.75 103.75
2gcM/NRT1 3.54239E-06 3.54239E-06 3.54239E-06 3.54239E-06 3.54239E-06 3.54239E-06 3.54239E-06 3.54239E-06 3.54239E-06 3.54239E-06 3.54239E-06 3.54239E-06 3.54239E-06
(k-1)/k
(k+1)/k
(k+1)/k
(k/(k+1)p1 (p1) - (p2) 3.503939786 8361547.6 3.503931472 8361303.3 3.503923157 8361058.9 3.503914843 8360814.5 3.503906528 8360570.1 3.503898213 8360325.8 3.503889899 8360081.4 3.503881584 8359837.0 3.503873269 8359592.7 3.503864954 8359348.3 3.503856639 8359103.9 3.503848324 8358859.6 3.503840009 8358615.2
2
G calculado 103.79 103.78 103.78 103.78 103.77 103.77 103.77 103.76 103.76 103.76 103.75 103.75 103.75
DELTA -0.03 -0.03 -0.03 -0.02 -0.02 -0.02 -0.01 -0.01 -0.01 0.00 0.00 0.00 0.01
Resultado: p1
2
8562.0 [lbf /ft ]
59.46 [psia]
4.04 [atm]
3. Calcular la velocidad de descarga de aire a la atmósfera, desde un recipiente a P0=150 psig y T0 = 70ºF a través de tubería de 33 ft de largo, con diámetro de 2”, Sch 40, acero comercial y 3 codos medio.
Figura. Problema ejemplo, flujo de gases Para el análisis, se tiene que determinar primero el valor de N.
MECÁNICA DE FLUIDOS
DPQBA, UTFSM
N = f ⋅
L D
Considerando flujo turbulento para un ε /D de 0.0004, se puede suponer un f = 0.016 Se calcula N:
Resistencia L/D
N=f·L/D
entrada cañería recta 3 codos
0.50 (se asume) 3.06 1.44 5.00
Σ
Haciendo los cambios de unidades adecuados: P0 = (150+14.7) = 164.7 psia = 164.7 · 144 = 23700 lbf/pie2 P3 = 14.7 psia = 14.7 · 144 = 2210 lbf/pie2 T0 = 70 ºF = 70 + 460 = 530 ºR
Se calcula según la ecuación (9.37): Gci = P0
G ci
= 23700 ⋅
g C ⋅ (PM )
2,718 ⋅ R ⋅ T 0
⎛ lb f ⎞ ft 3 ⎜ ft 2 ⎟ ⎝ ⎠ R = 1545 lb-molº R
32 .17 ⋅ 29 2 .718 ⋅ 1545 ⋅ 530
==> Gci = 486 [lbm/ft2·s] D = 2.067/12 = 0.1722 ft A = 0.0233 ft2 Caso
Isotérmico
Adiabático
k (P3/P0) N (G/Gci) Gci G (lb/pie·s) G·A = W T2/T0 T2 (ºF)
1.0 0.0893 5.0 0.545 486 265 6.18 1 79
1.4 0.0893 5.0 0.565 486 275 6.41 0.833 -18
MECÁNICA DE FLUIDOS
DPQBA, UTFSM
Temp. media
µ Re
70
(70-18)/2 = 26
1.21·10-5 lbm/ft·s 3.77·106
1.14·10-5 lbm/ft·s 4.16·106
Presiones P2 y P1
(P1/P0) (P2/P0) G Gci
Isotérmico
Adiabático
0.91 0.32
0.89 0.26
= cte
MECÁNICA DE FLUIDOS
DPQBA, UTFSM
f = 0,016
