Fancher & Brown PDF
Fancher & Brown PDF
Fancher & Brown PDF
Caracas, 2007
TRABAJO ESPECIAL DE GRADO
Caracas, 2007
ii
iii
iv
DEDICATORIA
v
AGRADECIMIENTOS
vi
Cravino C., Annabella y Pérez B., Al D.
vii
ÍNDICE DE CONTENIDO
DEDICATORIA....................................................................................................................................V
AGRADECIMIENTOS ...................................................................................................................... VI
ÍNDICE DE CONTENIDO ............................................................................................................. VIII
ÍNDICE DE FIGURAS.........................................................................................................................X
ÍNDICE DE TABLAS....................................................................................................................... XII
NOMENCLATURA, SÍMBOLOS Y ABREVIATURAS............................................................. XIV
INTRODUCCIÓN................................................................................................................................. 1
CAPÍTULO I: PROBLEMA EN ESTUDIO ...................................................................................... 4
PLANTEAMIENTO DEL PROBLEMA ........................................................................................... 4
OBJETIVOS Y ALCANCES ............................................................................................................. 6
METODOLOGÍA .............................................................................................................................. 7
CAPÍTULO II: REVISIÓN BIBLIOGRÁFICA .............................................................................. 12
ASPECTOS FUNDAMENTALES .................................................................................................. 12
ANTECEDENTES........................................................................................................................... 27
CAPÍTULO III: DESCRIPCIÓN DETALLADA DE MODELOS ................................................ 28
MODELOS HOMOGÉNEOS .......................................................................................................... 28
Modelo Homogéneo de Mc Adams et al. (1942).......................................................................... 28
Correlación de Poettman y Carpenter (1952) ............................................................................. 29
Modelo Homogéneo de Cicchitti et al. (1960)............................................................................. 32
Correlación de Baxendell y Thomas (1961) ................................................................................ 33
Correlación de Fancher y Brown (1963)..................................................................................... 35
Modelo Homogéneo de Wallis (1969) ......................................................................................... 37
Correlación de Oliemans (1976) ................................................................................................. 38
Correlación de Beattie y Whalley (1982) .................................................................................... 41
Modelo Homogéneo de Ouyang (1997)....................................................................................... 43
MODELOS DE CORRELACIÓN ................................................................................................... 46
Correlación de Lockhart y Martinelli (1949) .............................................................................. 47
Correlación de Reid et al. (1957) ................................................................................................ 52
Correlación de Hoogendorn (1959) ............................................................................................ 55
Correlación de Griffith y Wallis (1961) ...................................................................................... 58
Correlación de Duns y Ros (1963) .............................................................................................. 64
Correlación de Dukler et al. (1964) ............................................................................................ 74
Correlación de Hagedorn y Brown (1965).................................................................................. 78
Correlación de Orkiszewski (1967) ............................................................................................. 84
Correlación Eaton et al. (1967)................................................................................................... 91
Correlación de Aziz et al. (1972)................................................................................................. 95
Correlación de Beggs y Brill (1973).......................................................................................... 100
Correlación de Griffith et al. (1973) ......................................................................................... 107
Correlación de Chierici et al. (1974) ........................................................................................ 112
Correlación de Kadambi (1981)................................................................................................ 115
Correlación de Müller y Heck (1986)........................................................................................ 118
viii
Correlación de Chen (2002)...................................................................................................... 119
Correlación Universal de García (2006) .................................................................................. 122
MODELOS MECANICISTAS ...................................................................................................... 127
Modelo de Taitel y Barnea (1990)............................................................................................. 127
Modelo de Ansari et al (1994) ................................................................................................... 132
Modelo Mecanicista de Ouyang (1999)..................................................................................... 141
Modelo de Gómez et al (2000)................................................................................................... 150
CAPÍTULO IV: PRESENTACIÓN Y ANÁLISIS DE RESULTADOS ...................................... 159
RESULTADOS OBTENIDOS....................................................................................................... 159
Base de Datos Experimentales .................................................................................................. 159
Evaluación del Desempeño de Modelos y Correlaciones en la Predicción del Gradiente de
Presión....................................................................................................................................... 160
SÍNTESIS Y CONCLUSIONES ...................................................................................................... 177
RECOMENDACIONES................................................................................................................... 179
REFERENCIAS BIBLIOGRÁFICAS ............................................................................................ 180
APÉNDICES...................................................................................................................................... 189
A-1: ECUACIONES DE AJUSTE DE DATOS EXPERIMENTALES ......................................... 189
Ecuaciones correspondientes a la correlación de Fancher y Brown (1963): ........................... 192
Ecuaciones correspondientes a la correlación de Lockhart y Martinelli (1949): ..................... 193
Ecuaciones correspondientes a la correlación de Hagedorn y Brown (1965):......................... 194
Ecuaciones correspondientes a la correlación de Griffith y Wallis (1961): ............................. 196
Ecuaciones correspondientes a la correlación de Duns y Ros (1963): ..................................... 197
Ecuaciones correspondientes a la correlación de Eaton et al.(1967):...................................... 200
Ecuaciones correspondientes a la correlación de Griffith et al. (1973): .................................. 201
A-2: ECUACIONES COMPLEMENTARIAS DEL MODELO DE ANSARI ET AL. (1994)......... 205
A-3: ECUACIONES COMPLEMENTARIAS DEL MODELO DE TAITEL Y BARNEA (1990):
........................................................................................................................................................ 214
A-4: ECUACIONES COMPLEMENTARIAS DEL MODELO MECANICISTA DE OUYANG
(1999): ............................................................................................................................................ 218
A-5: ECUACIONES COMPLEMENTARIAS DEL MODELO DE GÓMEZ ET AL. (2000): ........ 222
ix
ÍNDICE DE FIGURAS
x
Fig. 33: Gradiente de presión calculado vs. experimental de los 10 modelos o
correlaciones con mejor desempeño para toda la base de datos experimentales. ..... 165
Fig. 34: Gradiente de presión calculado vs. experimental de los 10 modelos o
correlaciones con mejor desempeño en flujo tipo tapón........................................... 167
Fig. 35: Gradiente de presión calculado vs. experimental de los 10 modelos o
correlaciones con mejor desempeño en flujo tipo espuma........................................ 169
Fig. 36: Gradiente de presión calculado vs. experimental de los 10 modelos o
correlaciones con mejor desempeño en flujo tipo ondulado..................................... 171
Fig. 37: Porcentaje de datos experimentales sobrepredichos y sobrepredicción
porcentual promedio de los 10 modelos o correlaciones con mejor desempeño en
flujo tipo tapón. ......................................................................................................... 173
Fig. 38: Porcentaje de datos experimentales sobrepredichos y sobrepredicción
porcentual promedio de los 10 modelos o correlaciones con mejor desempeño en
flujo tipo espuma....................................................................................................... 174
Fig. 39: Porcentaje de datos experimentales sobrepredichos y sobrepredicción
porcentual promedio de los 10 modelos o correlaciones con mejor desempeño en
flujo tipo ondulado. ................................................................................................... 175
Fig. 40: Comparación entre la curva original del factor K de Griffith et al. y la
ecuación de ajuste, Griffith (1973)............................................................................ 190
Fig. 41: Comparación entre la curva original del factor L2 de Duns y Ros y la
ecuación de ajuste, Beggs y Brill (1988). ................................................................. 191
xi
ÍNDICE DE TABLAS
xii
Tabla 24: Valores de los coeficientes de las funciones de ajuste de la correlación de
Hagedorn y Brown (1965). ....................................................................................... 195
Tabla 25: Intervalos de aplicación de las ecuaciones de ajuste de la correlación de
Hagedorn y Brown (1965). ....................................................................................... 195
Tabla 26: Funciones de ajuste de datos experimentales de la correlación de Griffith y
Wallis (1961)............................................................................................................. 196
Tabla 27: Valores de los coeficientes de las funciones de ajuste de la correlación de
Griffith y Wallis (1961). ........................................................................................... 196
Tabla 28: Intervalos de aplicación de las ecuaciones de ajuste de la correlación de
Griffith y Wallis (1961). ........................................................................................... 196
Tabla 29: Funciones de ajuste de datos experimentales de correlación de Duns y Ros
(1963). ....................................................................................................................... 197
Tabla 30: Valores de los coeficientes de las funciones de ajuste de la correlación de
Duns y Ros (1963). ................................................................................................... 198
Tabla 31: Intervalos de aplicación de las ecuaciones de ajuste de la correlación de
Duns y Ros (1963). ................................................................................................... 199
Tabla 32: Funciones de ajuste de datos experimentales de la correlación de Eaton et
al.(1967).................................................................................................................... 200
Tabla 33: Valores de los coeficientes de las funciones de ajuste de la correlación de
Eaton et al.(1967)...................................................................................................... 200
Tabla 34: Intervalos de aplicación de las ecuaciones de ajuste de la correlación de
Eaton et al.(1967)...................................................................................................... 200
Tabla 35: Funciones de ajuste de datos experimentales de la correlación de Griffith et
al.(1973).................................................................................................................... 201
Tabla 36: Valores de los coeficientes de las funciones de ajuste de la correlación de
Griffith et al.(1973)................................................................................................... 202
Tabla 37: Intervalos de aplicación de las ecuaciones de ajuste de la correlación de
Griffith et al.(1973)................................................................................................... 204
xiii
NOMENCLATURA, SÍMBOLOS Y ABREVIATURAS
Nomenclatura:
xiv
f f: Factor de fricción de Fanning.
f i: Factor de fricción de Fanning en la interfase.
f1 : Factor de fricción de Moody para flujo monofásico de líquido.
f2 : Correlación de fricción para burbuja.
f3 : Factor de corrección.
Fe: Fracción de líquido presente en el núcleo
Fr: Número de Froude.
F1,2,3,4: Números de la velocidad del deslizamiento para burbuja.
F5,6,7: Números de la velocidad del deslizamiento para tapón.
′ Factores de corrección.
F3,6 :
xv
La : Frontera de patrón de flujo.
Lb : Frontera entre flujo burbuja y tapón, longitud para la burbuja.
Lm : Frontera entre flujo transitorio y anular.
Ls : Frontera entre flujo tapón y transitorio, longitud del tapón.
M: Peso molecular.
m: Factor adimensional.
m& : Flujo másico.
n: Número de experimentos.
Nd: Número del diámetro.
Ne: Factor adimensional.
Ngv: Número de la velocidad del gas.
Ngv’: Número de la velocidad del gas de Griffith et al.
Nl : Número de la viscosidad del líquido.
Nlv: Número de la velocidad del líquido.
Nlv’: Número de la velocidad del líquido de Griffith et al.
Nv: Factor adimensional.
Nx,y: Factor de régimen de flujo.
Nμ : Número adimensional dependiente de la viscosidad del líquido.
N1,2,3: Fronteras de flujo.
P: Presión.
Q: Flujo volumétrico.
qI : Tasa volumétrica radial por longitud de tubería.
xvi
Rl: Relación de los factores de densidad por el multiplicador de fricción.
Rv: Factor de deslizamiento.
S: Perímetro.
s: Parámetro usado para el cálculo de la relación de factores de fricción.
Sa: Velocidad del deslizamiento adimensional.
T: Temperatura.
t: Parámetro de ajuste.
T1 : Factor dependiente del influjo o del reflujo.
U: Velocidad.
Ubf: Velocidad de ascenso de la burbuja alongada.
Ubs: Velocidad de ascenso de la burbuja alongada en una columna
estacionaria de líquido.
Uo: Velocidad de deslizamiento de las burbujas.
Us: Velocidad del deslizamiento.
Usc: Velocidad superficial en el núcleo de gas.
Ut: Velocidad traslacional de las burbujas alongadas.
v: Velocidad relativa.
Vb : Volumen de la burbuja.
w: Tasa de flujo de masa total.
We: Número de Weber.
x: Fracción de flujo de masa de gas.
XM: Factor de Lockhart y Martinelli.
XMll: Factor de Lockhart y Martinelli para flujo laminar-laminar.
XMtt: Factor de Lockhart y Martinelli para flujo turbulento-turbulento.
YM: Factor de Lockhart y Martinelli.
Y: Factor de proporcionalidad.
y: Factor de cálculo.
Z: Factor de correlación.
xvii
Símbolos griegos:
α: Fracción de gas.
β: Ángulo de inclinación de la tubería con respecto a la horizontal, diferencia
entre las fracciones volumétricas de líquido con y sin deslizamiento.
βaE : Coeficiente de expansión.
δ: Coeficiente de distribución del líquido, espesor de la película de líquido,
coeficiente de correlación.
Δ: Variación o diferencia.
Δt: Período de una burbuja.
ε: Rugosidad total, tasa de disipación de energía por unidad de masa.
ξ: Variación entre el valor calculado y el experimental.
φ: Error relativo porcentual, factor adimensional.
φ: Parámetro adimensional de presión.
μ: Viscosidad dinámica.
λ: Fracción volumétrica de gas o líquido sin deslizamiento.
π: Número pi.
ρ: Densidad.
ρ ’: Densidad modificada.
ρa: Densidad promedio.
Ω: Factor de corrección.
σ: Tensión superficial.
τ: Esfuerzo de corte.
τ ’: Esfuerzo de corte que existiría si la tubería estuviera llena de líquido.
υ: Volumen específico.
ψ: Factor de corrección.
xviii
Subíndices y superíndices:
a: Debido a la aceleración.
c: En el núcleo.
calc: Calculado.
CB: Criterio de “creaming” de la burbuja.
CD: Crítico para evitar aglomeracion y coalescencia.
b: Burbuja.
exp: Experimental.
f: Película de líquido, debido a la fricción.
g: Gas, debido al cambio de energía potencial.
gls: De gas en el tapón de líquido.
gdb: De las burbujas dispersas en el tapón de líquido.
h: Hidráulico.
hom: Homogéneo.
i: Contador, interno, interfase.
l: Líquido.
lls: De líquido en el tapón.
ltb: Del líquido en la burbuja de Taylor.
lu: De lìquido en la unidad de tapón.
m: De la mezcla.
máx: Máximo.
min: Mínimo.
n: Sin deslizamiento.
o: De ascenso libre.
o∞ : De ascenso de la burbuja.
s: Relativo al tapón.
sg: Superficial del gas.
sc: Superficial en el núcleo.
xix
sl: Superficial del líquido.
t: Traslacional.
tot: Total.
tb: En la burbuja de Taylor.
tp: Bifásico, de las dos fases.
u: Unidad de tapón.
w: En la pared.
Abreviaturas:
xx
INTRODUCCIÓN
1
flujo bifásico, se han desarrollado modelos simplificados, tales como los modelos
homogéneos, las correlaciones empíricas y los modelos mecanicistas, entre otros.
2
En el capítulo II, se presentan los aspectos fundamentales del flujo bifásico y
se muestra una amplia revisión bibliográfica, que incluye a los principales modelos
para evaluar el gradiente de presión de flujo bifásico, presentes en la literatura
especializada.
3
Capítulo I: Problema en Estudio
Los primeros modelos para el cálculo de flujo vertical de dos fases surgieron
alrededor de 1950 y se distinguen por haber sido desarrollados de manera muy
limitada, es decir, se restringían a ciertos diámetros de tubería y a ciertos intervalos
de las propiedades de los fluidos, tales como la densidad o la viscosidad. En general,
estos modelos se caracterizan por asumir que el flujo es casi homogéneo y, por tanto,
se le asignan propiedades promedio a la mezcla.
4
Capítulo I: Problema en Estudio
5
Capítulo I: Problema en Estudio
OBJETIVOS Y ALCANCES
1. GENERAL
2. ESPECÍFICOS
6
Capítulo I: Problema en Estudio
METODOLOGÍA
7
Capítulo I: Problema en Estudio
1 n
E1 = ∑ ϕi
n i =1
Ec. 1
8
Capítulo I: Problema en Estudio
1 n
E2 = ∑ ϕi
n i =1
Ec. 2
1 1
E3 = ∑ (ϕ i − E1 )2 Ec. 3
n − 1 i =1
1 1
E4 = ∑ (ϕi )2 Ec. 4
n − 1 i =1
El error promedio:
1 n
E5 = ∑ ξi
n i =1
Ec. 5
1 n
E6 = ∑ ξi
n i =1
Ec. 6
9
Capítulo I: Problema en Estudio
1 n
E7 = ∑ (ξ i − E5 )2 Ec. 7
n − 1 i =1
1 n
E8 = ∑ (ξ i )2 Ec. 8
n − 1 i =1
⎛ ΔPcalc − ΔPexp ⎞
Siendo: ϕ i = ⎜ ⎟ ⋅ 100 , ξ i = ΔPcalc − ΔPexp y n el número de datos
⎜ Δ P ⎟
⎝ exp ⎠
experimentales.
10
Capítulo I: Problema en Estudio
11
Capítulo II: Revisión Bibliográfica.
ASPECTOS FUNDAMENTALES
Donde:
m& : tasa de flujo másico total
12
Capítulo II: Revisión Bibliográfica.
Q = Ql + Qg = U sl A + U sg A Ec. 10
Siendo:
Q: tasa de flujo volumétrico total
Ql: tasa de flujo volumétrico de líquido
Qg: tasa de flujo volumétrico de gas
De la fase líquida:
Ql
Ul = Ec. 11
Al
De la fase gaseosa:
Qg
Ug = Ec. 12
Ag
13
Capítulo II: Revisión Bibliográfica.
De la fase líquida:
Ql
U sl = Ec. 13
A
De la fase gaseosa:
Qg
U sg = Ec. 14
A
U m = U sl + U sg Ec. 15
Al
Hl = Ec. 16
A
14
Capítulo II: Revisión Bibliográfica.
Ql
Al Al Ql
Ul U sl U
λl = = = = = = sl Ec. 17
A Al + Ag Ql Q g Ql + Q g U sl + U sg U m
+
Ul U g
Ag
αg = = 1 − Hl Ec. 18
A
Ag
αg = = 1 − λl Ec. 19
A
De líquido:
m& l
Gl = = ρ lU sl Ec. 20
A
De gas:
m& g
Gg = = ρ gU sg Ec. 21
A
15
Capítulo II: Revisión Bibliográfica.
Total:
m& l + m& g
Gtot = Ec. 22
A
64
f = Ec. 23
Re
⎡ε ⎡ ⎡ ε ⎤ ⎤⎤
1,11
1 ⎢ Di 4,518 ⎢ 6,9 ⎢ Di ⎥ ⎥ ⎥
= −2 log ⎢ − log ⎢ +
⎢ ⎥ ⎥⎥ Ec. 24
f ⎢ 3,7 Re ⎢ Re ⎢⎣ 3,7 ⎥⎦ ⎥ ⎥
⎣ ⎣ ⎦⎦
ρ l DiU sl
Re l = Ec. 25
μl
ρ g DiU sg
Re g = Ec. 26
μg
16
Capítulo II: Revisión Bibliográfica.
ρ m DiU m
Re m = Ec. 27
μm
17
Capítulo II: Revisión Bibliográfica.
• Flujo Burbuja: A bajas tasas de líquido la fase gaseosa puede encontrarse como
pequeñas burbujas en el líquido, distribuidas uniformemente en la sección
transversal. Existe deslizamiento entre las fases.
• Flujo Tapón: En este tipo de flujo el gas se concentra en una gran burbuja
conocida como Burbuja de Taylor, la cual tiene casi el mismo diámetro de la
tubería. Entre una burbuja y la siguiente se forma un tapón de líquido. La película
de líquido que se encuentra entre la burbuja y la pared de la tubería penetra el
tapón formando burbujas en el mismo.
• Flujo Anular: A altas tasas de gas se puede formar una película de líquido
ascendente en la pared de la tubería con gas fluyendo en la parte interior,
arrastrando el líquido. Los esfuerzos que se producen en la interfase son muy
elevados.
18
Capítulo II: Revisión Bibliográfica.
Modelos Homogéneos
19
Capítulo II: Revisión Bibliográfica.
⎛ dP ⎞ ⎛ dP ⎞ ⎛ dP ⎞ ⎛ dP ⎞
⎜ ⎟ =⎜ ⎟ +⎜ ⎟ +⎜ ⎟ Ec. 28
⎝ dz ⎠ tot ⎝ dz ⎠ f ⎝ dz ⎠ a ⎝ dz ⎠ g
⎛ dP ⎞ f ⋅ ρ m ⋅U m2 dU m
⎜ ⎟ = + ρ m ⋅U m + gρ m ⋅ sin β Ec. 29
⎝ dz ⎠ tot 2 ⋅ Di dL
Con deslizamiento:
ρ m = ρ l ⋅ H l + ρ g ⋅ (1 − H l ) Ec. 30
μ m = μ l ⋅ H l + μ g ⋅ (1 − H l ) Ec. 31
20
Capítulo II: Revisión Bibliográfica.
Sin deslizamiento:
ρ m = ρ l ⋅ λl + ρ g ⋅ (1 − λl ) Ec. 32
μ m = μ l λl + μ g (1 − λl ) Ec. 33
Modelos de Correlación
Los modelos de correlación son aquellos en los que sus autores proponen una
serie de ajustes de datos experimentales para correlacionar una variable determinada.
Estas correlaciones pueden considerar tanto el deslizamiento como la existencia de
patrones de flujo; por tanto, requieren de métodos para determinar el patrón de flujo
presente. Una vez que se ha determinado el patrón de flujo correspondiente a unas
condiciones dadas, se determina la correlación apropiada para el cálculo del factor de
fricción como el holdup, las cuales, generalmente, son distintas dependiendo del
patrón.
Modelos Mecanicistas
Los modelos mecanicistas son aquellos en los que sus autores presentan un
análisis bastante minucioso de las variables que supuestamente influyen en el flujo
bifásico y simulan los mecanismos de flujo en función del patrón existente. Estos
modelos están basados en leyes fundamentales y generalmente buscan ser aplicables
en cualquier patrón de flujo, diámetro de tubería y condiciones de la mezcla bifásica.
21
Capítulo II: Revisión Bibliográfica.
Algunos modelos mecanicistas son: Taitel et al. (1980), Barnea (1987), Hasan
y Kabir (1988), Ansari et al. (1994), Fernandes et al. (1983), entre otros.
dP d
= ρ mU m (U m ) Ec. 34
dz a dz
Considerando que:
Gtot
Um = Ec. 35
ρm
dP G d ⎛ Gtot ⎞
= ρ m tot ⎜ ⎟ Ec. 36
dz a ρ m dz ⎜⎝ ρ m ⎟⎠
22
Capítulo II: Revisión Bibliográfica.
dP 2 d ⎛ 1 ⎞
= (Gtot ) ⎜ ⎟ Ec. 37
dz a dz ⎜⎝ ρ m ⎟⎠
d ⎛ 1 ⎞ 1 d
⎜ ⎟= (ρ m ) Ec. 38
dz ⎜⎝ ρ m ⎟⎠ ρ m2 dz
d ⎛ 1 ⎞ 1 d
⎜ ⎟= (ρ l ⋅ H l + ρ g ⋅ (1 − H l )) Ec. 39
dz ⎜⎝ ρ m ⎟⎠ ρ m2 dz
23
Capítulo II: Revisión Bibliográfica.
d ⎛ 1 ⎞ (1 − H l ) d
⎜ ⎟= (ρ ) Ec. 40
dz ⎜⎝ ρ m ⎟⎠ ρ m2 dz g
PM RT
ρg = óυg = Ec. 41
RT PM
dP 2 (1 − H l ) M dP
= −(Gtot ) Ec. 42
dz a ρ m2 R T dz tot
Gg
x= Ec. 43
Gtot
1 x 1− x
= + Ec. 44
ρm ρg ρl
1
= υ g x + υ l (1 − x ) Ec. 45
ρm
24
Capítulo II: Revisión Bibliográfica.
d ⎛ 1 ⎞ d
⎜⎜ ⎟⎟ = (υ g x ) = x d (υ g ) Ec. 46
dz ⎝ ρ m ⎠ dz dz
d
(υ g ) = d ⎛⎜⎜ RT ⎞⎟⎟ = − RT 2 dP Ec. 47
dz dz ⎝ MP ⎠ MP dz tot
d
(υ g ) = − M 2 dP Ec. 48
dz RTρ g dz tot
dP 2 x M dP
= −(Gtot ) 2 Ec. 49
dz a ρ g R T dz tot
ρm = ρ g
(1 − λl )2 + ρl
λl 2
Ec. 50
1 − Hl Hl
25
Capítulo II: Revisión Bibliográfica.
dP (
= −Gt2 ⎜
l l m g l sl ) (
l l m g g l l ) (
⎛ 2 H (1 − H )U G ρ U 2 + 2 H (1 − H )U G 3υ − H 2 (1 − H )U 2 G 2
m g ) ⎞⎟ ⋅ ⎛⎜
M ⎞ dP
⎟
dz a
⎜
⎝ (
(1 − H l )ρ l U l + H l G g υ g
2 2 2
) ⎟ ⎜ R Tρ 2
⎠ ⎝ g
⎟ dz
⎠ Ec. 51
ρ l λl + ρ g (1 − H l )
ρ tp = Ec. 52
(1 − β l )
Se aplica un proceso similar a los anteriores para obtener una ecuación que
contiene las variables a derivar en forma más directa, finalmente se obtiene el término
de aceleración del gradiente de presión para la ecuación de Oliemans, expresado
como:
⎛
⎜
(
⎜ (1 − H )G G + (1 − H )2 G ρ U + (1 − H )2 G 2
l l g l g l sl l g )⎛⎜⎜ R−TMρ 2
⎞
⎟
⎟
dP ⎝ g ⎠+
= −Gt2 ⎜
dz a ⎜ ((1 − H l )ρ gU sl + Gl + (1 − H l )G g )2
⎜
⎜
⎝ Ec. 53
( )
(1 − H l ) U mU sl + (1 − H l )U sl ⎛⎜ M ⎞⎟ ⎞⎟
2 2
⎝ R T ⎠ ⎟ dP
− ⋅
((1 − H l )ρ g U sl + Gl + (1 − H l )G g )2 ⎟⎟ dz tot
26
Capítulo II: Revisión Bibliográfica.
ANTECEDENTES
27
Capítulo III: Descripción Detallada de Modelos.
MODELOS HOMOGÉNEOS
Para los autores, la caída de presión se calcula aplicando la Ec. 28, la cual se
define para el modelo de Wallis (1969). La principal diferencia entre este último
modelo y el de Mc Adams radica en la definición de las propiedades del fluido. Dicho
esto, la viscosidad de mezcla μ m se calcula de la forma:
1 x 1− x
= + Ec. 54
μm μg μl
PROCEDIMIENTO DE CÁLCULO:
28
Capítulo III: Descripción Detallada de Modelos.
(b) Obtener el factor de fricción para la mezcla con la Ec. 24, para el número de
Reynolds de mezcla Ec. 27.
29
Capítulo III: Descripción Detallada de Modelos.
5. La pérdida de energía varía a lo largo del tubo pero se tomó como un valor
promedio constante a lo largo del mismo.
dP 1 ⎡g w2 f ⎤
= ⎢ ρm + 5⎥ Ec. 55
dz tot 144 ⎣ g c 7,413 × 1010 ρ m Di ⎦
Donde:
(dP dz )tot : caída de presión [psi/ft].
30
Capítulo III: Descripción Detallada de Modelos.
Di ρ mU m
Donde:
31
Capítulo III: Descripción Detallada de Modelos.
PROCEDIMIENTO DE CÁLCULO:
Para los autores la caída de presión se calcula aplicando la Ec. 29, la cual se
define para el modelo de Wallis (1969). La principal diferencia entre este último
modelo y el de Cicchitti radica en la definición de las propiedades del fluido. Dicho
esto, la viscosidad de mezcla se calcula de la forma:
μ m = μ g x + μ l (1 − x ) Ec. 59
PROCEDIMIENTO DE CÁLCULO:
32
Capítulo III: Descripción Detallada de Modelos.
(f) Obtener el factor de fricción para la mezcla con la Ec. 24, para el número de
Reynolds de mezcla Ec. 27.
El valor del factor de fricción de Moody fue graficado versus el numerador del
número Reynolds ( Di ρ mU m , expresado en [ lbm ft ⋅ s ]) como se muestra en la Fig. 4:
33
Capítulo III: Descripción Detallada de Modelos.
Di ρ mU m
Fig. 4: Correlación para el factor de fricción propuesta por Baxendell y Thomas (1961).
PROCEDIMIENTO DE CÁLCULO:
34
Capítulo III: Descripción Detallada de Modelos.
35
Capítulo III: Descripción Detallada de Modelos.
Di ρ mU m
Fig. 5: Correlación para el factor de fricción propuesta por Fancher y Brown (1963).
36
Capítulo III: Descripción Detallada de Modelos.
PROCEDIMIENTO DE CÁLCULO:
2
fGtot
+ ρ m g sin( β )
dP 2 Di ρ m
− = Ec. 60
dz tot 2 d ⎛ 1 ⎞
1 + Gtot ⎜ ⎟
dz ⎜⎝ ρ m ⎟⎠
37
Capítulo III: Descripción Detallada de Modelos.
aunque hace referencia a otras ecuaciones, como las propuestas por McAdams y
Cicchitti.
PROCEDIMIENTO DE CÁLCULO:
(a) Con los valores anteriores determinar el número de Reynolds para la mezcla.
(b) Calcular el factor de fricción de Moody para la mezcla con la Ec. 24.
38
Capítulo III: Descripción Detallada de Modelos.
como también afecta la tasa de flujo másico total por unidad de área, la viscosidad y
la densidad bifásicas.
m&
Gtot = Ec. 62
(1 − β l )A
μ l λl + μ g (1 − H l )
μm = Ec. 63
(1 − β l )
39
Capítulo III: Descripción Detallada de Modelos.
Como se puede notar, las cuatro ecuaciones anteriores dependen del factor βl,
el cual es la fracción volumétrica ocupada por el líquido estático, y viene dado por la
ecuación:
β l = H l − λl Ec. 64
Di m&
Re tp = Ec. 65
A (1 − β l )μ m
2
dP Gtot
= f Ec. 66
dz f 2 Def ρ m
40
Capítulo III: Descripción Detallada de Modelos.
PROCEDIMIENTO DE CÁLCULO:
(f) Calcular la caída de presión por elevación como se muestra en la Ec. 29.
Estos autores desarrollaron una correlación en la cual los patrones de flujo son
parcialmente considerados, es decir, solo se toman en cuenta en la definición de la
viscosidad de la mezcla. Beattie y Whalley reconocen los patrones anular y burbuja,
para los cuales determinaron una ecuación de viscosidad correspondiente. Sin
embargo, en este mismo modelo, sus autores proponen que las definiciones
particulares de las viscosidades dependientes del patrón de flujo sean reemplazadas
por una definición híbrida, que es la siguiente:
41
Capítulo III: Descripción Detallada de Modelos.
ρl x
α= Ec. 68
ρ l x + ρ g (1 − x )
2
Gtot xM
Ek = − 2 Ec. 70
ρ g RT
PROCEDIMIENTO DE CÁLCULO:
(a) Calcular el flujo másico total por unidad de área y la fracción de gas con la Ec. 22
y la Ec. 43 respectivamente.
42
Capítulo III: Descripción Detallada de Modelos.
(c) Determinar la fracción de vacío homogénea usando la Ec. 68 y luego sustituir este
resultado en la Ec. 67 para obtener la viscosidad de la mezcla bifásica.
dP τ wS
=− Ec. 71
dz f A
Donde:
fρ mU tp
2
τw = Ec. 72
8
43
Capítulo III: Descripción Detallada de Modelos.
dP β aE ⎡ dP dP ⎤ βaE ⎡ dP ⎤
= ⎢ + ⎥= ⎢ ⎥ Ec. 73
dz a 1 − β aE ⎣⎢ dz f dz g ⎦⎥ 1 − βaE ⎣ dz tot ⎦
ρ mU mU sg
βaE = Ec. 74
P
U sg
Hl = 1− Ec. 75
C oU m + U o
gσ (ρ l − ρ g )
U o = 1,54 ⋅ 4 ⋅ sin β Ec. 76
ρ l2
ρg
Co = C w − 0,2 Ec. 77
ρl
44
Capítulo III: Descripción Detallada de Modelos.
PROCEDIMIENTO DE CÁLCULO:
⎛ρ ⎞ ⎛ρ ⎞
U tp = ⎜ l ⎟U sl + ⎜ g ⎟U sg Ec. 78
⎜ρ ⎟ ⎜ρ ⎟
⎝ tp ⎠ ⎝ tp ⎠
(h) Calcular la caída de presión por elevación como se muestra en la Ec. 29.
45
Capítulo III: Descripción Detallada de Modelos.
MODELOS DE CORRELACIÓN
ρl
N gv = U sg ⋅ 4 Ec. 79
g ⋅σ
ρl
N lv = U sl ⋅ 4 Ec. 80
g ⋅σ
g
N l = μl ⋅ 4 Ec. 81
ρl ⋅ σ 3
ρl ⋅ g
N d = Di ⋅ Ec. 82
σ
46
Capítulo III: Descripción Detallada de Modelos.
(a) Factor adimensional que determina la transición entre burbuja y tapón (Lb):
Lb = 1,071 −
(0,2218 ⋅ U ) 2
m
Ec. 83
Di
Donde:
U m : Velocidad de la mezcla [ft/s].
(b) Factor adimensional que establece la frontera entre tapón y flujo transitorio (Ls):
Ls = 50 + 36 ⋅ N lv Ec. 84
(c) Factor adimensional que determina el cambio de flujo transitorio a anular (Lm):
Lm = 75 + 84 ⋅ N lv0,75 Ec. 85
47
Capítulo III: Descripción Detallada de Modelos.
f l ρ lU sl
2
dP
2 dz sl 2 Di
XM = = Ec. 86
f g ρ gU sg
2
dP
dz sg 2 Di
48
Capítulo III: Descripción Detallada de Modelos.
Fig. 6: Curvas que relacionan φ g , φl , Rl , Rg y X para todos los mecanismos de flujo, Lockhart
Martinelli (1949)
dP dP 2 dP
= φl = φg
2
Ec. 87
dz f dz sl dz sg
49
Capítulo III: Descripción Detallada de Modelos.
2. El volumen ocupado por el líquido más el ocupado por el gas debe ser igual a
volumen total de la tubería.
En tal sentido los autores graficaron la data experimental para los cuatro tipos
de flujo, obteniendo así una relación gráfica de las variables adimensionales, las
cuales fueron aproximadas por Chisholm (1967), de la forma:
Cc 1
φl 2 = 1 + +
XM XM2 Ec. 88
φ g = X M + Cc X M + 1
2 2
Donde el valor de Cc depende del numero de Reynolds para el líquido y del gas.
Re l Re g Cc
<2300 <2300 5
>2300 <2300 10
<2300 >2300 12
>2300 >2300 20
50
Capítulo III: Descripción Detallada de Modelos.
PROCEDIMIENTO DE CÁLCULO:
51
Capítulo III: Descripción Detallada de Modelos.
(j) Calcular la caída de presión por elevación como se muestra en la Ec. 29.
Los autores observaron que había muy poca data sobre caída de presión para
tuberías de diámetro mayor a 3,5 pulgadas. El propósito de su investigación fue
extender los resultados de Baker (1954) y validar las correlaciones para grandes
diámetros de Chenoweth-Martin (1955) y Lockhart-Martinelli.
52
Capítulo III: Descripción Detallada de Modelos.
53
Capítulo III: Descripción Detallada de Modelos.
dP dP −0,86
= λl Ec. 89
dz f dz l
dP
Siendo la caída de presión calculada como si el líquido fluyera solo en la tubería
dz l
a la misma tasa que las dos fases:
2
dP fG
= l l Ec. 90
dz l 2 Di ρ l
54
Capítulo III: Descripción Detallada de Modelos.
PROCEDIMIENTO DE CÁLCULO:
(d) Con la densidad de la mezcla Ec. 32, calcular la caída de presión por elevación
como se muestra en la Ec. 29.
55
Capítulo III: Descripción Detallada de Modelos.
⎤⎡ ⎤
Hg
dP ⎡
0 ,84
dP ⎛ Gg ⎞ ⎛ρ ⎞
= ⋅ ⎢1 + 230⎜⎜ ⎟⎟ ⎥ ⎢0,00138⎜ l ⎟ ⎥ Ec. 91
dz dz l ⎢ ⎥⎦ ⎢ ⎜ρ ⎟ ⎥
f
⎣ ⎝ Gtot ⎠ ⎣ ⎝ g ⎠ ⎦
0 ,85
Hg ⎡ ⎛ U sg U sl ⎞⎤
= A⎢ H g ⎜ − ⎟⎥ Ec. 92
1− Hg ⎜ ⎟
⎢⎣ ⎝ H g 1 − H g ⎠⎥⎦
Donde A = 0,60(m s )
−0 ,85
. Para obtener el valor de la fracción volumétrica de gas se
diseñó un proceso iterativo basado en el método de bisección.
La ecuación propuesta por Hoogendorn para el caso del flujo estratificado liso
u ondulado no es objeto de este trabajo, sin embargo, en el caso de flujo, anular el
mencionado autor recomienda se utilice la siguiente fórmula:
⎛ Gg 2 ⎞
= 0,12(G g )
dP −0 , 25 ⎜ ⎟ Ec. 93
dz ⎜ 2 Di ρ g ⎟
f ⎝ ⎠
56
Capítulo III: Descripción Detallada de Modelos.
Donde:
Gg : tasa de flujo de flujo másico por unidad de área [ kg m 2 ].
ρ g : densidad de la mezcla [ kg m 3 ].
PROCEDIMIENTO DE CÁLCULO:
a) Determinar los flujos másicos por unidad de área de líquido, de gas y total con las
ecuaciones Ec. 20, Ec. 21 y Ec. 22 respectivamente.
57
Capítulo III: Descripción Detallada de Modelos.
El flujo tapón se caracteriza por largas burbujas, las cuales ocupan casi toda la
tubería, y están separadas entre ellas por tapones de líquido. La parte superior de la
burbuja es redondeada mientras que la inferior es plana. Es importante señalar que en
el tapón se puede conseguir pequeñas burbujas. La forma típica de una burbuja en el
flujo tipo tapón se muestra en la Fig. 10:
58
Capítulo III: Descripción Detallada de Modelos.
(Ls + Lb )A
Δt = Ec. 94
Ql + Q g + U b A
Vb
Qg = Ec. 95
Δt
59
Capítulo III: Descripción Detallada de Modelos.
⎡ Vb ⎤ ⎡ Vb ⎤
ρ a = ρ l ⎢1 − ⎥ + ρg ⎢ ⎥ Ec. 96
⎣ (Ls + Lb )A ⎦ ⎣ (Ls + Lb )A ⎦
⎡ Ql + Vb A ⎤ ⎡ Qg ⎤
ρ a = ρl ⎢ ⎥ + ρg ⎢ ⎥ Ec. 97
⎣⎢ Ql + Q g + Vb A ⎦⎥ ⎣⎢ Ql + Q g + Vb A ⎦⎥
U b = C1C 2 ( gDi )
0,5
Ec. 98
60
Capítulo III: Descripción Detallada de Modelos.
Fig. 11: Constante adimensional C1 vs. Número de Reynolds de la burbuja, Griffith y Wallis
(1961).
ρ lU b Di
Re b = Ec. 99
μl
ρ lU m Di
Re l = Ec. 100
μl
61
Capítulo III: Descripción Detallada de Modelos.
Fig. 12: Constante adimensional C2 vs. Número de Reynolds de líquido para varios Reynolds de
la burbuja, Griffith y Wallis (1961).
62
Capítulo III: Descripción Detallada de Modelos.
dP 4τ
= ρ a g sin (β ) + w Ec. 101
dz tot Di
Donde:
⎡ Ql ⎤ '
τw = ⎢ ⎥τ w Ec. 102
⎣⎢ Ql + Q g ⎦⎥
fρ U
2
τ = l m
'
w
Ec. 103
8
PROCEDIMIENTO DE CÁLCULO:
63
Capítulo III: Descripción Detallada de Modelos.
incluyendo una de transición. Se desarrolló una correlación para obtener, para cada
tipo de flujo, el coeficiente de fricción y la velocidad del deslizamiento, la cual
permite calcular la fracción volumétrica. Es importante destacar que, para asegurar la
correcta aplicación elas correlaciones, los datos fueron convertidos al sistema ingles
de unidades.
64
Capítulo III: Descripción Detallada de Modelos.
Donde:
μ : viscosidad [cp].
σ : tensión superficial [ dynes cm ].
Di : diámetro de la tubería [ft].
ρ m : densidad [ lbm ft 3 ].
U : velocidad [ ft s ].
P: presión [psf].
65
Capítulo III: Descripción Detallada de Modelos.
Fig. 14: Números de la velocidad de deslizamiento para burbuja, Duns y Ros (1963).
66
Capítulo III: Descripción Detallada de Modelos.
F5
Fig. 15: Números de la velocidad de deslizamiento para tapón, Duns y Ros (1963).
67
Capítulo III: Descripción Detallada de Modelos.
0 , 25
⎛ρ ⎞
S a = 1,938U s ⎜⎜ l ⎟⎟ Ec. 104
⎝σl ⎠
Donde:
U sg U sl
Us =U g −Ul = − Ec. 105
1− Hl Hl
Hl = s
(
U − U m + (U s − U m ) + 4U sU sl
2
)0,5
Ec. 106
2U s
68
Capítulo III: Descripción Detallada de Modelos.
Flujo Burbuja
2
′ ⎛ N gv ⎞
S a = F1 + F2 N lv + F3 ⎜⎜ ⎟⎟ Ec. 107
⎝ 1 + N lv ⎠
dP f ρU U g
= m l sl m + ρ m sin (β ) Ec. 108
dz tot 2 g c Di gc
f2
f m = f1 Ec. 109
f3
Donde f1 es el factor de fricción de Moody para flujo monofásico como función del
número de Reynolds para el líquido Ec. 25, f 2 factor de corrección obtenido
correlacionando la relación de velocidades de las fases y un grupo adimensional,
como se muestra en la Fig. 16. Y el factor de corrección:
0,5
⎛ U sg ⎞
f 3 = 1 + f1 ⎜⎜ ⎟⎟ Ec. 110
⎝ 50U sl ⎠
69
Capítulo III: Descripción Detallada de Modelos.
Flujo Tapón
⎛ N 0,982 + F ′ ⎞
S a = (1 + F5 )⎜⎜
gv 6 ⎟
Ec. 111
⎜ 1 + 7 N lv ⎟⎟
2
⎝ ⎠
Para este tipo de flujo, los autores concluyeron que no hay deslizamiento entre
las fases. Es por esto que la caída de presión se calcula:
fρ gU sg
2
+ ρm g
dP 2 g c Di Ec. 112
=
dz t 1 − Ek
U mU sg ρ m
Ek = Ec. 113
gc P
70
Capítulo III: Descripción Detallada de Modelos.
εU sg 2 ρ g
We = Ec. 114
σ
Y:
μl 2
Nμ = Ec. 115
ρ lσε
ε 0,0749σ
WeN μ ≤ 0,005 ⇒ = Ec. 116
Di ρ gU sg 2 Di
Y:
ε 0,3713σ
WeN μ > 0,005 ⇒ = (WeN μ )0,302 Ec. 117
Di ρ gU sg Di
2
En el caso que ε
Di > 0,05 el factor de fricción de Moody se calcula:
⎡ 1 1, 73 ⎤
f = 4⎢ + 0,067( ε Di ) ⎥
⎢⎣ (4 log(0,27 ε Di ))
Ec. 118
⎥⎦
Flujo Transición
En este caso los autores propusieron hacer una interpolación lineal entre las
fronteras de flujo Ls y Lm , es decir entre flujo tapón y tipo “mist”. Se utiliza la
ecuación:
71
Capítulo III: Descripción Detallada de Modelos.
dP Lm − N gv ⎛ dP ⎞ N gv − Lm ⎛ dP ⎞
= ⎜ ⎟+ ⎜⎜ ⎟⎟ Ec. 119
dz tot Lm − Ls ⎜⎝ dz ⎟ L −L
tapon ⎠ m s ⎝ dz mist ⎠
ρ g N gv
ρ g' = Ec. 120
Lm
PROCEDIMIENTO DE CÁLCULO:
0 ≤ N gv ≤ L1 + L2 N lv Ec. 121
L1 + L2 N lv ≤ N gv ≤ Ls Ec. 122
72
Capítulo III: Descripción Detallada de Modelos.
de F3' se calcula:
F4
F3' = F3 − Ec. 125
Nd
2) Usar las Ec. 107, Ec. 104 y Ec. 106 para determinar Sa, U s y H l .
F6' se calcula:
73
Capítulo III: Descripción Detallada de Modelos.
2) Usar las Ec. 111, Ec. 104 y Ec. 106 para determinar Sa, U s y H l .
3) Repetir el paso 3), 4) y 5) para flujo burbuja.
1) Calcular la caída de presión para flujo tapón y tipo mist siguiendo los
pasos c) y d) utilizando como densidad del gas la calculada en la Ec. 120.
2) Obtener la caída de presión de la Ec. 119.
Este modelo se caracteriza por estar dividido en dos partes: en la primera, sus
autores comparan ciertas correlaciones para el cálculo de la presión y del holdup
frente a una data experimental; y en la segunda, realizan un análisis de similitud con
74
Capítulo III: Descripción Detallada de Modelos.
2
dP 2Gtot ff
= α (λ )ϕ Ec. 127
dz f Di ρ m
ρ l λl 2 ρ g (1 − λl )2
ϕ= + Ec. 128
ρ m H l ρ tm (1 − H l )
En la cual, la densidad bifásica debe ser calculada según la Ec. 30, utilizando
la fracción de líquido sin deslizamiento. El holdup debe ser calculado a través de un
proceso iterativo en el cual se utilizan las ecuaciones de Bankoff (1960) y Hughmark
(1962). Bankoff propuso la siguiente ecuación para correlacionar la fracción de
líquido con deslizamiento:
1− Hl
K= Ec. 129
1 − λl
1 1
⎛ ρ mU m Di ⎞ ⎛ Um ⎞
6 8
⎜ ⎟ ⎜⎜ ⎟⎟
⎜ μ H + μ (1 − H ) ⎟
⎝ gDi Ec. 130
1 1
Re 6 Fr ⎠
=⎝
l ⎠
8
l l g
Z=
λl λl
1 1
4 4
75
Capítulo III: Descripción Detallada de Modelos.
− ln λl
α (λ ) = 1 +
1,281− 0,478(− ln λl ) + 0,444(− ln λl ) − 0,094(− ln λl ) + 0,00843(− ln λl ) Ec. 132
2 3
Otro aspecto importante de este modelo tiene que ver con la definición del
gradiente de presión debido a los cambios de energía cinética. Haoulo (2004)
desarrolló detalladamente las ecuaciones propuestas por Dukler y compañía, hasta
llegar a la relación mostrada en la Ec. 51.
PROCEDIMIENTO DE CÁLCULO:
(a) Determinar el flujo másico total por unidad de área con las Ec. 20, Ec. 21 y la Ec.
22.
76
Capítulo III: Descripción Detallada de Modelos.
(d) Determinar la fracción de líquido con deslizamiento, para lo cual debe realizarse
un proceso iterativo que involucre las ecuaciones Ec. 129, Ec. 130 y Ec. 131.
(e) Calcular el factor adimensional φ según la Ec. 128, el flujo másico total m& con
las Ec. 9 y luego sustituir en la siguiente fórmula para determinar el número de
Reynolds de la mezcla homogénea.
4m&
Re tp = ϕ Ec. 133
πDi μ m
0,125
f f = 0,014 + 0 , 32 Ec. 134
Re tp
(g) Determinar el factor α(λ) según la Ec. 132 y luego calcular el volúmen específico
del gas de esta forma:
RT
υg = Ec. 135
PM
(h) Determinar el gradiente de presión por fricción con la Ec. 127 y el gradiente de
presión por aceleración con la Ec. 51.
(i) Con la densidad de la mezcla, calcular la caída de presión por elevación como se
muestra en la Ec. 29.
77
Capítulo III: Descripción Detallada de Modelos.
Las correlaciones propuestas por estos dos autores fueron desarrolladas a partir de
datos experimentales provenientes de pozos de prueba cuyas tuberías tenían a lo
sumo 1 ½ pulgadas de diámetro nominal. No obstante, Hagedorn y Brown aseguran
que es posible utilizar sus correlaciones en el cálculo de gradientes de presión en
tuberías de mayor diámetro que aquellas donde fueron desarrolladas dichas
correlaciones. Según los autores esta extrapolación es posible con cierto grado de
precisión. Es importante destacar que, para asegurar la correcta aplicación elas
correlaciones, los datos fueron convertidos al sistema ingles de unidades.
Donde:
μ : viscosidad [cp].
σ : tensión superficial [ dynes cm ].
Di : diámetro de la tubería [ft].
ρ m : densidad [ lbm ft 3 ].
U : velocidad [ ft s ].
P, dP/dz: presión [psi].
viscosidad del líquido N l . Para ello en este trabajo se hizo un ajuste de los datos
reales, obtenidos de la siguiente figura, a una función polinómica de décimo orden, la
cual puede verse en el apéndice A-1, en la Tabla 23.
78
Capítulo III: Descripción Detallada de Modelos.
Fig. 17: Coeficiente de la viscosidad del líquido versus Nl, según Hagendorn y Brown(1965).
⎛ N lv ⎞⎛ P ⎞ 0,1 ⎛ C nl ⎞
C2 = ⎜ ⎟⎜ ⎟ ⎜ ⎟⎟ Ec. 136
⎜ N 0,575 ⎟⎝ 14,7 ⎠ ⎜⎝ N d ⎠
⎝ gv ⎠
79
Capítulo III: Descripción Detallada de Modelos.
Fig. 18: Factor de Holdup vs. C2, por Hagendorn y Brown (1965).
⎛ N gv ⋅ N lv0,38 ⎞
C3 = ⎜ ⎟ Ec. 137
⎜ N 2,14 ⎟
⎝ d ⎠
80
Capítulo III: Descripción Detallada de Modelos.
⎛H ⎞
H f = ⎜ ff ⎟ ⋅ψ Ec. 138
⎝ ψ ⎠
PROCEDIMIENTO DE CÁLCULO:
N lv (Ec. 80).
N l (Ec. 81).
N d (Ec. 82).
(b) Calcular el valor de la velocidad de la mezcla (Ec. 15) y de la frontera entre flujo
burbuja y tapón (Ec. 83).
81
Capítulo III: Descripción Detallada de Modelos.
(d) Calcular el número de Reynolds del líquido (Ec. 25) y el factor de fricción de
Moody (Ec. 24).
N gv ≤ 0 Ec. 139
λ g ≤ Lb Ec. 140
dP ρ U 2 g ρl
= 1,0792 ⋅ 10 −4 ⋅ f ⋅ l sl + sin( β ) Ec. 141
dz tot Di g c 144
U s = 0,7874016 [ ]
ft
s Ec. 142
⎛ ⎛ Um ⎞
2
U ⎞
⎜ Um ⎟
H l = 1 − 0,5 ⋅ ⎜1 + − ⎜⎜1 + ⎟⎟ − 4 ⋅ sg ⎟⎟ Ec. 143
⎜ Us ⎝ Us ⎠ Us
⎝ ⎠
U sl
Ul = Ec. 144
Hl
Calcular el factor de fricción de Moody según (Ec. 24) usando el resultado anterior.
82
Capítulo III: Descripción Detallada de Modelos.
dP 1 ⎛ g ρ lU sl2 ⎞
= ⎜⎜ ρ m sin β + f ⋅ ⎟ Ec. 145
dz tot 144 ⎝ gc 2 g c Di H l2 ⎟⎠
C2 (Ec. 136).
H ff
ψ (Ec. 328).
C3 (Ec. 137).
ψ (Ec. 330).
H f (Ec. 138).
λl (Ec. 17).
ρm = ρl ⋅ H l + ρ g ⋅ H g Ec. 146
ρ n = ρ l ⋅ λl + ρ g ⋅ λ g Ec. 147
83
Capítulo III: Descripción Detallada de Modelos.
ρn2
ρ tp = si ρ m > ρ n
ρm Ec. 148
ρ tp = ρ n si ρ m < ρ n
λl :
dP 1 ⎡⎛ g ⎞ g fρ tpU m2 ⎤
= ⎢ ⎜ ρ sin β ⎟ + ⎥ Ec. 150
dz tot 144 ⎢⎣⎜⎝ g c ⎟g
tp
⎠ c 2 g c Di ⎥⎦
84
Capítulo III: Descripción Detallada de Modelos.
patrones de flujo que Orkiszewski reconoce son: burbuja, tapón, transitorio y anular;
y para cada uno de ellos, o bien el autor recomienda el uso de una determinada
correlación que según su criterio es la más adecuada, o bien el autor desarrolla una
serie de correlaciones que conducen al cálculo del gradiente de presión. Es importante
destacar que, para asegurar la correcta aplicación de las correlaciones, los datos
fueron convertidos al sistema ingles de unidades.
Donde:
μ : viscosidad [cp].
σ : tensión superficial [ dynes cm ].
Di : diámetro de la tubería [ft].
ρ m : densidad [ lbm ft 3 ].
U : velocidad [ ft s ].
dP/dz: presión [psf].
85
Capítulo III: Descripción Detallada de Modelos.
dP ρ lU m2 ⎛ (U sl + U b ) ⎞
= f⋅ ⋅ ⎜⎜ + δ ⎟⎟ Ec. 151
dz f 2 g c Di ⎝ (U m + U b ) ⎠
86
Capítulo III: Descripción Detallada de Modelos.
Ec. 98. En dicha ecuación, C1 se expresa en la Fig. 11 como función del número de
Reynolds de flujo burbuja y C2 se extrae de la Fig. 12 como función tanto del número
de Reynolds de flujo burbuja como del número de Reynolds del líquido. Como se dijo
previamente, se hallaron funciones de ajuste de los datos de ambas gráficas para ser
utilizados en las correspondientes subrutinas.
U b ⋅ Di ⋅ ρ l
Re b = Ec. 152
μl
Cuando Re b ≤ 3.000 :
Cuando Re b ≥ 8.000 :
87
Capítulo III: Descripción Detallada de Modelos.
ψ = (0,251 + 8,74 × 10 −6 ⋅ Re )⋅ g ⋅ Di
1 13,59 ⋅ μ l Ec. 155
U b = ⋅ψ + ψ 2 +
2 ρ l ⋅ Di
FASE ECUACIÓN
U m [ ft s ]
CONTINUA PROPUESTA
⎛ 0,013 log μ l ⎞
δ = ⎜⎜ ⎟⎟ − 0,681 + 0,232 logU m − 0,428 log Di Ec. 156
⎝ Di7 ,99 ⎠
⎛ 0,045 log μ l ⎞
δ = ⎜⎜ ⎟⎟ − 0,709 − 0,162 logU m − 0,888 log Di Ec. 157
⎝ Di7 ,99 ⎠
⎛ 0,0127 log(μ l + 1) ⎞
δ = ⎜⎜ 1, 415
⎟⎟ − 0,284 + 0,167 logU m + 0,113 log Di Ec. 158
⎝ D i ⎠
88
Capítulo III: Descripción Detallada de Modelos.
⎛ 0,0274 log(μ l + 1) ⎞
δ = ⎜⎜ ⎟⎟ + 0,161 + 0,569 log Di
⎝ Di1,371 ⎠
Ec. 159
⎡⎛ 0,01 log(μ l + 1) ⎞ ⎤
− log U m ⎢⎜⎜ 1, 571
⎟⎟ + 0,397 + 0,63 log Di ⎥
⎣⎝ Di ⎠ ⎦
Ub ⎛ ρ tp ⎞
δ ≥− ⋅ ⎜⎜1 − ⎟⎟ Ec. 161
Um + Ub ⎝ ρ l ⎠
Donde:
ρ l (U sl + U b ) + ρ gU sg
ρ tp = + ρ lδ Ec. 162
Um +Ub
PROCEDIMIENTO DE CÁLCULO:
89
Capítulo III: Descripción Detallada de Modelos.
b) Si el flujo es tipo tapón, calcular el número de Reynolds del líquido (Ec. 25),
utilizando como velocidad la de la mezcla, el factor de fricción de Moody (Ec.
24) y la rugosidad relativa.
c) Asignar los valores iniciales del proceso iterativo para calcular la velocidad de la
burbuja. Asumir una primera aproximación a la velocidad de la burbuja y calcular
con ella el número de Reynolds de la burbuja (Ec. 152).
e) Comparar esta nueva velocidad con la supuesta en el paso (c) y seguir iterando
hasta que la diferencia entre el valor supuesto y el calculado sea menor que una
tolerancia especificada o se alcance el número máximo de iteraciones.
90
Capítulo III: Descripción Detallada de Modelos.
Donde:
μ : viscosidad [cp].
σ : tensión superficial [ dynes cm ].
Di : diámetro de la tubería [ft].
ρ m : densidad [ lbm ft 3 ].
U : velocidad [ ft s ].
dP/dz: presión [psf].
m& : flujo másico [lb/s].
91
Capítulo III: Descripción Detallada de Modelos.
Fig. 20: Correlación para el factor de pérdida de energía de Eaton et al. (1967).
92
Capítulo III: Descripción Detallada de Modelos.
0 , 05
1,84 N lv0,575 ⎛ P ⎞
⎜⎜ ⎟⎟ N l0,1 Ec. 164
N gv N d0,0277 ⎝ 14,65 psi ⎠
Fig. 21: Correlación para el Holdup propuesta por Eaton et al. (1967).
93
Capítulo III: Descripción Detallada de Modelos.
PROCEDIMIENTO DE CÁLCULO:
(a) Calcular los flujos másicos de gas, líquido y total. Determinar la abscisa de la Fig.
20.
(c) Calcular los números adimensionales necesarios para utilizar la Ec. 164.
Determinar la abscisa de la Fig. 21.
2
dP Gtot
= f Ec. 165
dz f 2 g c Di ρ m
(g) Con el holdup calculado en (d) determinar la densidad bifásica, Ec. 30, y con ésta
el gradiente de presión debido al cambio de energía potencial. Aplicar el modelo
de Wallis (1969) para determinar la caída de presión debida al cambio de energía
cinética.
94
Capítulo III: Descripción Detallada de Modelos.
Donde:
μ : viscosidad [cp].
σ : tensión superficial [ dynes cm ].
Di : diámetro de la tubería [ft].
ρ m : densidad [ lbm ft 3 ].
U : velocidad [ ft s ].
dP/dz: presión [psf].
m& : flujo másico [lb/s].
Para definir los patrones de flujo, los autores definen las siguientes variables,
(no se debe hacer cambio de unidades):
1 1
⎛ ρg ⎞ 3
⎛ 72 ρ l ⎞ 4
N x = U sg ⎜⎜ ⎟ ⎜⎜ ⎟⎟ Ec. 166
⎟
⎝ 0, 0764 ⎠ ⎝ 62,4σ l ⎠
1
⎛ 72 ρ l ⎞ 4
N y = U sl ⎜⎜ ⎟⎟ Ec. 167
⎝ 62,4σ l ⎠
N 1 = 0,51(100 N y )
0 ,172
Ec. 168
95
Capítulo III: Descripción Detallada de Modelos.
N 3 = 70(100 N y )
−0 ,152
Ec. 170
Flujo Burbuja:
U sg
Hl =1− Ec. 171
U bf
Siendo:
⎡ σg (ρ l − ρ g )⎤
1
4
U bs = 1,41⎢ ⎥ Ec. 173
⎣ ρl 2 ⎦
fρ U
2
dP
= m m Ec. 174
dz f 2 g c Di
96
Capítulo III: Descripción Detallada de Modelos.
Flujo Tapón:
Para este caso se utilizan las Ec. 171 y Ec. 172, y la siguiente:
⎡σ g (ρ l − ρ g )⎤
1
4
U bs = C ⎢ ⎥ Ec. 175
⎣ ρl 2 ⎦
Donde:
⎡ 3, 37 − N e
⎤
[
C = 0,345 1 − e − 0,029 N v ⎢1 − e m ] ⎥ Ec. 176
⎢⎣ ⎥⎦
Siendo:
gDi (ρ l − ρ g )
2
Ne = Ec. 177
σ
Nv =
[gD ρ (ρ − ρ )]
i
3
l l g
1
2
Ec. 178
μl
Nv m
≥ 250 10
−0 , 35
250 > N v > 18 69 N v
≤ 18 25
97
Capítulo III: Descripción Detallada de Modelos.
fρ H U
2
dP
= l l m Ec. 179
dz f 2 g c Di
Donde f es el factor de fricción de Moody correspondiente al número de Reynolds de
líquido calculado para una velocidad igual a la de la mezcla.
Flujo Mist:
Flujo Transición:
Para este tipo de patrón de flujo para el cálculo de la caída de presión se debe
usar las ecuaciones de flujo tapón y “mist”. Luego se aplica un método de
interpolación similar al propuesto por Duns y Ros (1963):
dP N 3 − N x dP ⎛ N − Nx ⎞ dP
= + ⎜⎜1 − 3 ⎟⎟ Ec. 180
dz tot N 3 − N 2 dz slug ⎝ N3 − N2 ⎠ dz mist
98
Capítulo III: Descripción Detallada de Modelos.
PROCEDIMIENTO DE CÁLCULO:
99
Capítulo III: Descripción Detallada de Modelos.
(d) Caída de presión para flujo “mist”: en este caso se utiliza el m´todo descrito para
Duns y Ros (1963).
100
Capítulo III: Descripción Detallada de Modelos.
f tp ρ mU m
2
dP
[ ]
g sin θ ρ l H l + ρ g (1 − H l ) +
2 Di
− =
dz
1−
[ ]
ρ l H l + ρ g (1 − H l ) U mU sg
Ec. 185
101
Capítulo III: Descripción Detallada de Modelos.
Fig. 22: Fracción Volumétrica de Líquido vs. Ángulo de inclinación, Beggs y Brill (1973).
Los patrones de flujo considerados por Beggs y Brill son los siguientes:
Fig. 23: Patrones de flujo según Beggs y Brill, Beggs y Brill (1973).
102
Capítulo III: Descripción Detallada de Modelos.
PROCEDIMIENTO DE CÁLCULO:
(a) Determinar la velocidad media (Ec. 15) y la fracción de líquido sin deslizamiento
Ec. 17.
(b) Calcular:
N lv (Ec. 80).
2
U
Fr = m Ec. 186
gDi
L1 = 316λl
0 , 302
Ec. 187
−2 , 4684
L2 = 0,0009252λl Ec. 188
−1, 4516
L3 = 0,1λl Ec. 189
−6 , 738
L4 = 0,5λl Ec. 190
SEGREGADO (segregated):
λl < 0,01 Y Fr < L1
O
λl ≥ 0,01 Y Fr < L2
TRANSICIÓN (transition):
λl ≥ 0,01 Y L2 ≤ Fr ≤ L3
103
Capítulo III: Descripción Detallada de Modelos.
INTERMITENTE (intermittent):
0,01 ≤ λl < 0,4 Y L3 < Fr ≤ L1
O
λl ≥ 0,4 Y L3 < Fr ≤ L4
DISTRIBUIDO (distributed):
λl < 0,4 Y Fr ≥ L1
O
λl ≥ 0,4 Y Fr > L4
aλl
b
Donde:
PATRÓN DE FLUJO a b c
(e) Cálculo de C:
(
C = (1 − λl ) ln dλl N lv Fr g
e f
) Ec. 192
104
Capítulo III: Descripción Detallada de Modelos.
Donde:
PATRÓN DE FLUJO d e f g
(f) Cálculo de Ψ:
[
Ψ = 1 + C sin (1,8β ) − 0,333 sin 3 (1,8β ) ] Ec. 193
TRANSICIÓN
L3 − Fr
A= Ec. 195
L3 − L2
105
Capítulo III: Descripción Detallada de Modelos.
1
fn = 2
⎡ ⎛ Re n ⎞⎤ Ec. 198
⎢2 log⎜⎜ − 3,8215 ⎟⎟⎥
⎣ ⎝ 4,5223 log(Re n ) ⎠⎦
(j) Clcular el factor de fricción de Moody bifásico, el cual depende del parámetro s:
Donde:
λl
y= 2 Ec. 200
Hl
R f = es Ec. 201
f tp = f n R f Ec. 202
106
Capítulo III: Descripción Detallada de Modelos.
Con respecto a los patrones de flujo, Griffith y compañía consideran solo flujo
tipo tapón y anular, y para discriminar entre uno y otro redefinen los siguientes
parámetros adimensionales:
ρl
N lv′ = U sl
gDi (ρ l − ρ g )
Ec. 203
ρg
′ = U sg
N gv
gDi (ρ l − ρ g )
Ec. 204
⎛ ρ ⎞
La = ⎜ 7 + 0,06 l ⎟ N lv′ Ec. 205
⎜ ρg ⎟
⎝ ⎠
107
Capítulo III: Descripción Detallada de Modelos.
Y:
U sg U sg
Hl = 1− = 1− Ec. 207
1,15U m + U o 1,15U m + K 2 gDi
108
Capítulo III: Descripción Detallada de Modelos.
1 1
Hl = 1− = 1−
U g ρ g (1 − x ) ρ g (1 − x ) Ec. 208
+1 Rv +1
Ul ρl x ρl x
dP ρ mU m2 Rl
= f
dz f 2 Di
Ec. 209
ρ
Donde : Rl = m S
ρl
109
Capítulo III: Descripción Detallada de Modelos.
Fig. 26: Factor multiplicador de la fricción propuesto por Griffith et al. (1973).
110
Capítulo III: Descripción Detallada de Modelos.
ρ mU m Di
Re = Ec. 211
μl
PROCEDIMIENTO DE CÁLCULO:
(b) Hallar el factor de fricción bifásico de Moody utilizando la Ec. 210, para el
número de Reynolds calculado en la Ec. 211.
(d) Calcular los flujos másicos del gas y del líquido según la Ec. 21 y la Ec. 22.
Determinar la calidad con la Ec. 43.
111
Capítulo III: Descripción Detallada de Modelos.
(e) Con la calidad y la relación entre las densidades del líquido y del gas, utilizar la
respectiva ecuación de ajuste de las curvas que se muestran en la Fig. 26 para
hallar el factor multiplicador de la fricción.
112
Capítulo III: Descripción Detallada de Modelos.
Donde:
μ : viscosidad [cp].
σ : tensión superficial [ dynes cm ].
Di : diámetro de la tubería [ft].
ρ m : densidad [ lbm ft 3 ].
U : velocidad [ ft s ].
dP/dz: presión [psf].
113
Capítulo III: Descripción Detallada de Modelos.
la densidad del mismo, cuando se requiera calcular los gradientes por elevación y
fricción en flujo anular.
PROCEDIMIENTO DE CÁLCULO:
Flujo Tapón
(a) Calcular el número de Reynolds del líquido (Ec. 25) sustituyendo U m en lugar de
(b) Asignar los valores iniciales del proceso iterativo para calcular la velocidad de la
burbuja. Asumir una primera aproximación a la velocidad de la burbuja y calcular
con ella el número de Reynolds de la burbuja (Ec. 152).
114
Capítulo III: Descripción Detallada de Modelos.
de ajuste de las Fig. 11 y Fig. 12 para hallar los factores necesarios para el cálculo
de la nueva velocidad de la burbuja.
(d) Comparar esta nueva velocidad con la supuesta en el paso (c) y seguir iterando
hasta que la diferencia entre el valor supuesto y el calculado sea menor que una
tolerancia especificada o se alcance el número máximo de iteraciones.
U sg
Hg = Ec. 215
Um +Ub
fρ U 2
= gρ m sin (β ) + (1 − H g ) l m
dP
Ec. 216
dz tot 2 Di
115
Capítulo III: Descripción Detallada de Modelos.
1 + 0,6684 X Mll
Hg = 2 Ec. 217
1 + 1,2424 X Mll + 0,201X Mll
2
1 + 1,719 X Mll + 0,3502 X Mll
φ g2 = Ec. 218
1 + 0,01436 X Mll
116
Capítulo III: Descripción Detallada de Modelos.
2
1 + 8,4493 X Mtt + 4,6431X Mtt
φ g2 = Ec. 220
1 + 5,7798 X Mtt
PROCEDIMIENTO DE CÁLCULO:
(a) Calcular el flujo másico total por unidad de área y la fracción másica de gas, con
las ecuaciones Ec. 21 y Ec. 22.
(b) Determinar el número de Reynolds para el gas, Ec. 26 y, con este valor,
determinar el factor de fricción de Moody para el gas f g con la Ec. 24
f g ρ gU sg
2
dP Ec.
=
dz sg 2 Di 221
(f) Calcular el gradiente de presión por fricción como en la Ec. 87 y el valor obtenido
en el paso (c).
117
Capítulo III: Descripción Detallada de Modelos.
(h) Con la fracción de vacío obtenida, determinar la densidad de la mezcla con la Ec.
30, recordando que H l = 1 − H g , y con ella determinar el gradiente de presión por
Este modelo esta basado en dos características importantes: la calidad del flujo
y la tasa de flujo másico total por unidad de área. Sus autores proponen inicialmente
calcular los números de Reynolds del gas y del líquido. Luego, si dicho número (del
gas o del líquido) es menor o igual a 1.187, se determina el factor de fricción laminar
respectivo según la Ec. 23, en caso contrario, el factor de fricción correspondiente se
determina con:
0,3164
f g ,l = Ec. 222
4 Re
g ,l
2
dP Gtot
= f g ,l ⋅ Ec. 223
dZ g ,l 2 Di ρ g ,l
118
Capítulo III: Descripción Detallada de Modelos.
dP ⎡ dP ⎛ dP dP ⎞⎟⎤ 3 3 dP dP
=⎢ + 2 x⎜
⎜
−
⎟ ⎥ ⋅ (1 − x ) + x + ρ m g sin (β ) + Ec. 224
dz tot ⎣⎢ dz l ⎝ dz g dz l ⎠⎦⎥ dz g dz a
Donde x es la calidad del flujo, es decir, la relación entre el flujo másico por
unidad de área del gas y del líquido.
PROCEDIMIENTO DE CÁLCULO:
(a) Calcular el flujo másico total por unidad de área y la fracción másica de gas, con
las ecuaciones Ec. 21 y Ec. 22.
(b) Determinar el número de Reynolds del líquido, Ec. 25, y del gas, Ec. 26. Luego,
calcular los respectivos factores de fricción de Moody, ya sea en régimen laminar
con la Ec. 23 o turbulento con la Ec. 222.
(e) Aplicar el modelo de Wallis para calcular el gradiente de presión por aceleración,
usando la suma a la izquierda de la componente de caída de presión por
aceleración, Ec. 224, y el holdup sin deslizamiento.
119
Capítulo III: Descripción Detallada de Modelos.
2
⎛ Di ⎞
⎜ 2⎟
Bo = g (ρ l − ρ g ) ⎝ ⎠ Ec. 225
σ
2
Gtot Di
We = Ec. 226
σρ m
2
G
Fr = tot Ec. 227
Di gρ m
120
Capítulo III: Descripción Detallada de Modelos.
dP dP
= Ω Ec. 228
dz f dz hom
−0 , 25
f m = 0,0791 Re m Ec. 231
PROCEDIMIENTO DE CÁLCULO:
121
Capítulo III: Descripción Detallada de Modelos.
(d) Determinar la caída de presión por elevación como se vio para otros modelos, y la
caída de presión por caeleración aplicando la Ec. 49.
122
Capítulo III: Descripción Detallada de Modelos.
f f = F2 +
(F1 − F2 )
d
⎛ ⎛ Re m ⎞ c ⎞ Ec. 233
⎜1 + ⎜ ⎟
⎜ ⎝ t ⎟⎠ ⎟
⎝ ⎠
Donde:
b1
F1 = a1 Re m
b2
Ec. 234
F2 = a2 Re m
Intervalo a1 b1 a2 b2 c d t
0,5 > HL ≥ 0,4 16,0077 -0,9815 0,3082 -0,2790 2,4364 0,0445 295
0,4 > HL ≥ 0,3 16,0136 -0,9198 0,0399 -0,1611 5,2260 0,1142 295
0,3 > HL ≥ 0,2 15,9999 -0,8619 0,1786 -0,2649 3,2130 0,1364 10.000
0,2 > HL ≥ 0,1 16,0110 -0,8385 0,0158 -0,0949 4,8460 0,1079 10.000
0,1 > HL ≥ 0,05 16,2165 -0,7890 0,1421 -0,2986 2,7274 0,2588 10.000
123
Capítulo III: Descripción Detallada de Modelos.
Y=F+
(1 − F )
d
⎛ ⎛ ⎛ 1 ⎞ Qg ⎞
c
⎞ Ec. 236
⎜1 + ⎜ ⎜ ⎟ ⎟⎟ ⎟
⎜ ⎜⎝ ⎝ t ⎠ Ql ⎠ ⎟
⎝ ⎠
b
F = a Re m Ec. 237
124
Capítulo III: Descripción Detallada de Modelos.
Donde:
Intervalo a b c d t
2 f f ρ mU m
2
dP
= Ec. 238
dz tot Di
PROCEDIMIENTO DE CÁLCULO:
lugar de U sl .
125
Capítulo III: Descripción Detallada de Modelos.
(b) Calcular la fracción volumétrica de líquido utilizando las Ec. 235 y Ec. 236,
utilizando la fracción volumétrica de líquido sin deslizamiento, Ec. 17. Para esto,
antes es necesario determinar F de la Ec. 237, utilizando la Tabla 8.
(d) Determinar la caída de presión por fricción aplicando la aplicando la Ec. 238 y la
de elevación como se muestra en la Ec. 29.
126
Capítulo III: Descripción Detallada de Modelos.
MODELOS MECANICISTAS
⎛ ⎛ hf ⎞⎞
f i = 0,005 ⋅ ⎜⎜1 − 300 ⋅ ⎜⎜ ⎟⎟ ⎟
⎟ Ec. 239
⎝ ⎝ Di ⎠⎠
127
Capítulo III: Descripción Detallada de Modelos.
τ fSf τ g Sg ⎛ 1 1 ⎞⎟
− − τ i Si ⎜ + + (ρ l − ρ g )g sin β
dh f Af Ag ⎜A A ⎟ Ec.
= ⎝ f g ⎠
dz
(ρl − ρ g )g cos β − ρ l vl (U t − U l )H lls dH ltb
− ρ g vg t
(U − U b )(1 − H lls ) dH ltb 240
2
H ltb dh f (1 − H ltb )2 dh f
La Ec. 240 puede resolverse de diversas maneras, sin embargo, los autores se
enfocan en 3 casos particulares: el primero, asumiendo la formulación general con
una aproximación a un canal de flujo unidimensional; el segundo, suponiendo que la
película de líquido es un canal de superficie libre; y finalmente, asumiendo que en la
128
Capítulo III: Descripción Detallada de Modelos.
zona de la burbuja existe una película de espesor constante. Los dos primeros
métodos requieren de integración numérica y son considerablemente complejos; no
obstante, a pesar de que el último método de resolución es mucho más simple, según
los autores los resultados en el cálculo de la caída de presión obtenidos por esta vía
están bastante cercanos a los obtenidos con los otros dos métodos. Por consiguiente,
ellos recomiendan el uso de ese enfoque y así se asumió en este trabajo.
τ fSf τ g Sg ⎛ 1 1 ⎞⎟
− − τ i Si ⎜ + + (ρ l − ρ g )g sin β = 0 Ec. 241
Af Ag ⎜A A ⎟
⎝ f g ⎠
dP 1 ⎛ τ πD τ S τ S ⎞
= ⎜⎜ ρ s g sin (β )Ls + s i Ls + ρ f g sin (β )L f + f f L f + g g L f ⎟⎟ Ec. 242
dz tot Lu ⎝ A A A ⎠
PROCEDIMIENTO DE CÁLCULO:
(a) Determinar la velocidad media (Ec. 15), el número de Reynolds del líquido (Ec.
25), el factor de fricción (Ec. 24) y la rugosidad relativa.
129
Capítulo III: Descripción Detallada de Modelos.
(c) Obtener la velocidad traslacional con la siguiente relación, en la cual los autores
indican que el valor de C es igual a 2 si el flujo es laminar e igual a 1,2 si el flujo
es turbulento:
Ut = C ⋅U m + U d Ec. 244
(e) Hallar la velocidad de las burbujas dispersas con las siguiente relación, en la cual
el valor de Bd es 1,2 sugerido por Wallis (1969):
U b = Bd ⋅ U m + U o Ec. 245
130
Capítulo III: Descripción Detallada de Modelos.
2
⎡ ⎛ 0,4σ ⎛ 2 fU m3 ⎞ ⎛ ρ l ⎞ 5 ⎞⎟
2
5 3
⎤
H lls = 1 − 0,058 ⋅ ⎢2⎜ ⋅ ⎜⎜ ⎟⎟ ⋅ ⎜ ⎟ − 0,725⎥
(ρ l − ρ g )g
Ec. 247
⎢ ⎜
⎣ ⎝ ⎝ Di ⎠ ⎝ σ ⎠ ⎟⎠ ⎥
⎦
(i) Hallar las áreas total At , ocupada por el líquido Al y ocupada por el gas Ag .
U m − U b ⋅ (1 − H lls )
Ul = Ec. 248
H lls
(k) Una vez calculados todos estos valores, asumir un espesor de película inicial y
luego utilizar el método de Aproximaciones Sucesivas para hallar la raíz de la Ec.
240 y de esta manera obtener el espesor de película correspondiente.
131
Capítulo III: Descripción Detallada de Modelos.
Fig. 28: Mapa de patrones de flujo utilizado por Ansari et al. (1994).
132
Capítulo III: Descripción Detallada de Modelos.
Flujo Burbuja
Flujo Tapón
Este modelo trata el flujo tapón totalmente desarrollado por separado al flujo
en desarrollo. Fernandes et al. (1986) el primer modelo completo para flujo tipo
tapón, el cual fue simplificado por Sylvester (1987), quien presentó una correlación
para calcula la fracción volumétrica de gas en el tapón. Ambos modelos son aplicados
para el caso de patrón desarrollado. McQuillan y Whalley (1985) introdujeron el
concepto de tapón en desarrollo.
133
Capítulo III: Descripción Detallada de Modelos.
Fig. 29: Esquema del flujo tapón desarrollado y en desarrollo, Ansari (1994).
cuando sale de la película con una velocidad U ltb hacia un tapón de líquido que se
mueve a una velocidad U lls . Si el cambio de velocidad sucede dentro del tapón
entonces no ocurrirá una caída neta de presión en el mismo, mientras el flujo tapón
sea estable. Es por esto que Ansari no toma en cuenta la componente de caída de
presión debido a la aceleración.
Flujo Anular
Wallis presentó una ecuación para el cálculo del arrastre de líquido y el factor
de fricción en función del espesor de película. Luego Hewitt y Hall-Taylor (1970)
ofrecieron un análisis exhaustivo de los mecanismos involucrados en el flujo tipo
anular. El valor de gradiente de presión total puede ser calculado por dos ecuaciones
134
Capítulo III: Descripción Detallada de Modelos.
⎛ dP ⎞ ⎛ dP ⎞ 2 ⎛ dP ⎞
⎜ ⎟ =⎜ ⎟ = φc ⎜ ⎟ + gρ c sin (β ) Ec. 249
⎝ dz ⎠ tot ⎝ dz ⎠ c ⎝ dz ⎠ sc
⎛ dP ⎞ ⎛ dP ⎞ 2 ⎛ dP ⎞
⎜ ⎟ = ⎜ ⎟ = φ f ⎜ ⎟ + gρ l sin (β ) Ec. 250
⎝ dz ⎠ tot ⎝ dz ⎠ film ⎝ dz ⎠ sl
PROCEDIMIENTO DE CÁLCULO:
(b) Para definir el patrón de flujo se deben calcular los puntos de corte de las curvas
de las fronteras de transición.
⎡ gσ (ρ l − ρ g )⎤
1
4
U sg = 0,25 ⋅ 1,53⎢ ⎥ + 0,333U sl Ec. 251
⎣ ρl 2 ⎦
135
Capítulo III: Descripción Detallada de Modelos.
⎡ (ρ l − ρ g )σ ⎤
1
2
Dmin = 19,01⎢ ⎥ Ec. 252
⎣ ρl g ⎦
2
1 3 2
⎡ 0,4σ ⎤ 2 ⎛ ρ l ⎞ 5 ⎛ f ⎞ ⎛ U sg ⎞
(U + U sg )
5
= 0,725 + 4,15⎜ ⎟
1, 2
2⎢ ⎥ ⎜ ⎟ ⎜ ⎟ Ec. 253
⎣⎢ (ρ l − ρ g )g ⎦⎥ ⎝ σ ⎠ ⎝ 2d ⎠
sl ⎜U +U ⎟
⎝ sg sl ⎠
U sg
U sl = − U sg Ec. 254
0,76
⎡ g (ρ l − ρ g )σ ⎤
1
4
136
Capítulo III: Descripción Detallada de Modelos.
sin (β )g (ρ l − ρ sc )
YM =
dP Ec. 256
dz sc
dP dP
Donde y corresponden a la caída de presión por fricción en el
dz sl dz sc
líquido y en el núcleo respectivamente, y se calculan:
f sl ρ lU sl
2
dP
= Ec. 257
dz sl 2 Di
f ρ U
2
dP
= sc sc sc Ec. 258
dz sc 2 Di
Donde los factores de fricción de Moody anteriores se calculan con la Ec. 24 para
número de Reynolds de líquido, Ec. 25, y el número de Reynolds para el núcleo:
ρ scU sc Di
Re sc = Ec. 259
μ sc
F (δ ) = YM −
(2 − 1.5 A1 )X M
2
Ec. 260
A1 (1 − 1.5 A1 )
3
Siendo A1:
A1 = 4δ (1 − δ ) Ec. 261
137
Capítulo III: Descripción Detallada de Modelos.
(c) Caída de presión para flujo burbuja y burbuja dispersa. Se itera, como se muestra
en el apéndice A-2, para un valor inicial de flujo volumétrico de líquido H l = 0,2
la ecuación:
⎛ gσ (ρ l − ρ g ) ⎞
0 , 25
U sg
1,53⎜⎜ ⎟ 0,5
⎟ Hl = + 1,2U m Ec. 262
⎝ ρ l
2
⎠ 1 − Hl
138
Capítulo III: Descripción Detallada de Modelos.
dP
= ρ m g sin (β ) Ec. 263
dz e
f tp ρ mU m
2
dP
= Ec. 264
dz f 2 Di
27.
(d) Caída de presión para flujo tapón. Se consideran dos casos, flujo tapón en
desarrollo y desarrollado, en ambos casos se calcula la relación entre la longitud
de la burbuja de Taylor y la longitud total, es decir la longitud de la burbuja más
la del tapón. Para calcular estos valores es necesario determinar las fracciones
volumétricas de líquido y gas de la burbuja y del tapón. Las ecuaciones que se
utilizaron se muestran en el apéndice A-2.
− U sl + U lls H lls
Re l = Ec. 265
U lls H lls + U ltb H ltb
Ls
Lu = Ec. 266
1 − Re l
139
Capítulo III: Descripción Detallada de Modelos.
dP
= ((1 − Re l )ρ lls + Re lρ )g sin (β ) Ec. 268
dz e
dP ρ lls f lsU m 2
= (1 − Re l ) Ec. 269
dz f 2 Di
ρ ls DiU m Ec.
Re ls =
μ ls 270
(e) Caída de presión para flujo anular. Se itera para valores de espesor de película, la
función:
140
Capítulo III: Descripción Detallada de Modelos.
2
Z XM
F (δ ) = YM − + Ec. 271
4δ 1 (1 − δ 1 )(1 − 4δ 1 (1 − δ 1 )) (4δ 1 (1 − δ 1 ))3
2,5
Luego se obtienen:
⎛ Z ⎞
φc 2 = ⎜⎜ ⎟
5 ⎟ Ec. 272
⎝ (1 − 2δ ) ⎠
φc 2 − YM Ec. 273
φf 2 = 2
XM
Los patrones de flujo considerados por este autor son: anular, burbuja, burbuja
dispersa, intermitente y estratificado. Este último no se toma en cuenta para este
estudio ya que el flujo vertical no asume esta configuración.
141
Capítulo III: Descripción Detallada de Modelos.
radial serán nulos, ya que se considera que no existe flujo (de gas o de líquido) en
esta dirección.
Flujo Anular:
Para este patrón de flujo, el esfuerzo por fricción con las paredes ocurre en la
película, es por esto que el balance de momento para el núcleo y para el líquido en la
película es:
dP 2ρ c
Ac = −τ i S i − ρ c Ac g sin (β ) − U c q Ig Ec. 274
dz 1 − Fe
dP
Af = τ i S i − ρ l A f g sin (β ) − τ wl S l − 2 ρ lU f T1 Ec. 275
dz
142
Capítulo III: Descripción Detallada de Modelos.
τ w Sl ⎛ 1 1 ⎞
− + τ i S i ⎜⎜ + ⎟⎟ − (ρ l − ρ c )g sin (β )
Al ⎝ Ac Al ⎠
Ec. 276
2ρ ⎛ Fe ⎞ 2 ρ cU c q Ig
− l U l ⎜ q Il − q Ig ⎟ + =0
Al ⎝ 1 − Fe ⎠ Ac 1 − Fe
La cual es equivalente a:
X 2 F2 − F1 − 4Y − 4 I = 0 Ec. 277
Donde:
⎧ dP
⎪
⎪ X 2 = dz sl = f slU sl ρ l
2
⎪ dP f sg U sg ρ g
2
⎪ dz sg
⎨ Ec. 278
⎪ (ρ − ρ c )g sin (β )
⎪Y = l
⎪ dP
⎪ dz sg
⎩
143
Capítulo III: Descripción Detallada de Modelos.
⎧ f i ρ c ⎛⎜ U c ⎞⎟
2
⎛ 1 1 ⎞
⎪ F1 = Di S i ⎜⎜ + ⎟⎟
⎪ f sg ρ g ⎜⎝ U sg ⎟⎠ A
⎝ l A c ⎠
⎪ 2
⎪⎪ f wl ⎛ U f ⎞ Di S l
⎨ F2 = ⎜⎜ ⎟⎟ Ec. 279
⎪ f sl ⎝ U sl ⎠ Al
⎪ −1
⎪ ⎧ 2ρ l ⎡ Fe ⎤ 2 ρ c U c q Ig ⎫⎛⎜ dP ⎞
⎟
⎪ I = ⎨ A U l ⎢⎣q Il − 1 − Fe q Ig ⎥⎦ − A 1 − Fe ⎬⎜ dz ⎟
⎪⎩ ⎩ l c ⎭⎝ sg ⎠
0 , 085
⎛ σ ⎞
f i = 0,24 f c Re f
0 , 305
⎜ ⎟ Ec. 280
⎜ ρ U 2D ⎟
⎝ c c c ⎠
El flujo tipo burbuja dispersa ocurre para altas tasas de líquido. En este caso se
estudia la mezcla como homogénea, en la cual puede o no existir deslizamiento entre
las fases, lo cual afecta únicamente al valor de la fracción volumétrica de líquido. En
la Fig. 31 se puede ver de manera esquemática todas las variables que el autor de este
modelo considera.
144
Capítulo III: Descripción Detallada de Modelos.
Fig. 31: Esquema del Flujo Burbuja Dispersa según Ouyang (1999).
ρl ρg
U tp = U sl + U sg Ec. 281
ρm ρm
⎛ dP ⎞ τ S 2
−⎜ ⎟ = w + ρ mU tp q I tp + ρ m g sin (β ) Ec. 282
⎝ dz ⎠ tot A A
145
Capítulo III: Descripción Detallada de Modelos.
⎡ 6,8 0 , 278 ⎛ D ⎞
0 ,112
⎤
U sl ≥ ⎢ 0, 444 (gσ (ρ l − ρ g )) ⎜⎜ i ⎟⎟ ⎥ Ec. 283
⎢⎣ ρ ⎝ μl ⎠ ⎥⎦
U sg
Hl = 1− Ec. 284
CoU m + U o
Flujo Intermitente:
Este patrón se caracteriza por la presencia de burbujas alargadas junto con una
película de líquido, seguido por un tapón de líquido que contiene burbujas dispersas.
La Fig. 32 presenta claramente las variables y consideraciones hechas por el autor:
146
Capítulo III: Descripción Detallada de Modelos.
⎛ dP ⎞ 2
fρ mU tp + ρ m g sin (β )
2
−⎜ ⎟ = Ec. 285
⎝ dz ⎠ tot Di
ρ m DiU tp Ec.
Re m =
μm 286
PROCEDIMIENTO DE CÁLCULO:
0,6
⎛σ ⎞
d b = 1,14⎜⎜ ⎟⎟ ⋅ ε −0 , 4 Ec. 287
⎝ ρl ⎠
(b) Asumir un valor para el espesor de película, iterar hasta hallar el correcto y
calcular el holdup según alguna correlación para flujo anular (p. e. Ansari et al.).
147
Capítulo III: Descripción Detallada de Modelos.
(
d CB ± = 0,5 A ± A 2 + 4 B ) Ec. 288
d CB − < d b < d CB +
0,4σ
d CD = 2
(ρ l − ρ g )g Ec. 289
d b < d CD
ρl ρg
U tp = U sl + U sg Ec. 292
ρ tp ρ tp
148
Capítulo III: Descripción Detallada de Modelos.
dP S
− = τ w + gρ tp sin β Ec. 293
dz tot A
149
Capítulo III: Descripción Detallada de Modelos.
Una característica importante de este modelo es que sus autores afirman que el
régimen es laminar si el número de Reynolds es menor o igual a 1.500 y no 2.300
como la mayoría de los investigadores.
150
Capítulo III: Descripción Detallada de Modelos.
⎡ ⎛ U sg ⎞
0,5
⎤⎛ σ ⎞
0, 6
⎛2f U3 ⎞
−0 , 4
3 ρl f U3
DCB = ⋅ ⋅ m m
8 (ρ l − ρ g ) g cos β
Ec. 298
Luego, Gómez et al. (2000) sugieren que la transición entre flujo burbuja y
tapón viene dada por varias condiciones. La primera de ellas es que el ángulo de
inclinación de la tubería tiene que estar entre 60º y 90º; en segundo término, el caudal
de líquido tiene que ser bajo, por tanto, la velocidad superficial máxima del gas está
dada por:
⎡ g (ρ l − ρ g )σ ⎤
1
1−α
4
U sg − 1,53(1 − α ) ⎥ ⋅ sin β
0,5
U sl = ⎢ Ec. 299
α ⎣ ρ l2 ⎦
151
Capítulo III: Descripción Detallada de Modelos.
Dmin = 19
(ρ l − ρ g )σ
Ec. 300
gρ l2
Flujo Tapón
(c) Uso de una nueva correlación (Gómez et al., 1999) en el cálculo de la fracción
de líquido con deslizamiento en el tapón de líquido.
152
Capítulo III: Descripción Detallada de Modelos.
dP πDi Ls τ wf S f + τ wg S g L f
= ρ u g sin β + τ s + Ec. 302
dZ tot A Lu A Lu
Siendo:
ρ u = ρ g (1 − H lu ) + ρ l H lu Ec. 303
Como se puede ver, la Ec. 242 y la Ec. 302 son muy parecidas, salvo por las
diferentes densidades a las que están referidas. El resto del modelo para el cálculo del
gradiente de presión en flujo tipo tapón es prácticamente idéntico al expuesto
anteriormente en la descripción del modelo de Taitel y Barnea (1990). Las ecuaciones
que los diferencian se pueden ver en el apéndice A-5.
Flujo Anular:
Sf ⎛ 1 1 ⎞⎟
τ wf − τ i Si ⎜ + + (ρ l − ρ c )g sin β = 0 Ec. 304
Af ⎜A ⎟
⎝ f Ac ⎠
153
Capítulo III: Descripción Detallada de Modelos.
f i = f sc ⋅ I Ec. 305
Flujo Burbuja:
Los autores extendieron el trabajo de Hasan y Kabir (1988) para poder aplicar
este modelo a todo el intervalo de inclinaciones. Por consiguiente, modifican la
ecuación de la velocidad del gas, haciéndola depender tanto del ángulo de inclinación
como del holdup. De manera que, la fórmula propuesta para el cálculo de la velocidad
del gas es:
dP τ
= ρ m g sin β + 4 m Ec. 307
dz tot Di
154
Capítulo III: Descripción Detallada de Modelos.
En el caso de que el patrón de flujo presente sea burbuja dispersa, los autores
recomiendan el uso del modelo homogéneo de Wallis (1969) para el cálculo del
gradiente de presión.
PROCEDIMIENTO DE CÁLCULO:
Di ≥ Dmin
Ec. 309
(1 − H l ) < 0,25
155
Capítulo III: Descripción Detallada de Modelos.
⎛ 1 1 ⎞⎟ τ wf S f τ wg S g
(ρ − ρ g )g sin β − τ s S i ⎜
⎜A
+ + − =0 Ec. 310
Ag ⎟⎠
l
⎝ f Af Ag
H l ≥ 0,24
Ec. 311
STAB > 0
fm ρm U m U m
τm = Ec. 312
2
dP τm4
= Ec. 313
dz f Di
dP
= ρ m g sin (β ) Ec. 314
dz e
156
Capítulo III: Descripción Detallada de Modelos.
L f = Lu − L s Ec. 318
ρ lU ltb U ltb
τ f = fl Ec. 319
2
ρ gU gtb U gtb
τ g = fg Ec. 320
2
ρ sU m 2
τs = f Ec. 321
2
157
Capítulo III: Descripción Detallada de Modelos.
dP − τ i Si
= − ρ c g sin (β ) Ec. 322
dz t Ag
ρ cU i 2
τ i = fi Ec. 323
2
158
Capítulo IV: Presentación y Análisis de Resultados
RESULTADOS OBTENIDOS
*
Valores aproximadamente constantes, según los respectivos autores.
159
Capítulo IV: Presentación y Análisis de Resultados
160
Capítulo IV: Presentación y Análisis de Resultados
161
Capítulo IV: Presentación y Análisis de Resultados
162
Capítulo IV: Presentación y Análisis de Resultados
Tabla 12: Comparación de la precisión de los 30 modelos o correlaciones de diferentes autores contra
los 74 datos experimentales.
Parámetros Estadísticos
Modelo o
E1 E2 E3 E4 E5 E6 E7 E8
Correlación
[%] [%] [%] [%] [Pa/m] [Pa/m] [Pa/m] [Pa/m]
CGR 4,44 13,95 17,96 18,50 120,76 435,24 558,25 571,33
MGO 9,05 16,49 20,67 22,59 268,53 514,65 632,68 688,03
CLM 16,19 17,64 17,81 24,15 483,32 538,26 531,75 720,81
MMO 3,48 18,08 24,37 24,62 108,11 583,86 775,80 783,40
MTB 15,24 18,57 21,34 23,84 483,24 592,15 653,53 739,04
CGA 11,76 21,24 25,23 27,87 342,55 678,05 837,92 906,13
CBB 11,36 22,09 27,10 29,42 378,45 707,52 846,02 927,87
CHB -0,08 22,48 27,89 27,89 -42,47 760,90 964,75 965,70
CAZ -1,39 24,58 29,81 29,84 7,68 740,61 875,07 875,10
CEA 15,99 25,14 30,80 34,75 547,75 823,84 994,74 1.137,38
CDR -10,14 25,62 28,31 30,10 -290,30 783,32 872,57 920,22
CDU -25,78 27,01 18,24 31,72 -858,19 893,53 622,49 1.064,93
MAN 1,82 28,16 32,45 32,51 99,68 895,29 1.042,95 1.047,77
CGW -25,82 29,35 25,92 33,30 -810,92 927,92 807,05 1.041,55
COR 6,98 36,41 58,36 58,79 284,87 1.251,15 2.151,12 2.170,16
CCHi 39,12 39,20 36,99 54,03 1.161,84 1.165,22 1.045,53 1.568,92
COL -37,85 39,36 21,82 43,91 -1.178,82 1.226,50 658,09 1.357,10
MHC -43,13 43,71 15,68 46,17 -1.376,21 1.394,27 508,34 1.475,91
CCH -43,08 44,78 21,49 48,40 -1.378,54 1.424,60 655,96 1.535,15
CMH -47,14 47,71 16,06 50,11 -1.499,56 1.517,30 510,27 1.593,70
MHMc -47,75 48,33 16,50 50,83 -1.516,00 1.534,01 516,25 1.611,29
MHO 48,61 48,61 23,38 54,24 1.473,15 1.473,15 574,24 1.590,49
CBW -48,37 48,95 16,87 51,54 -1.534,60 1.552,66 526,01 1.632,16
CRE -49,29 49,85 17,12 52,49 -1.560,95 1.578,68 525,59 1.657,17
MHW -50,82 51,38 17,66 54,13 -1.609,86 1.627,44 548,10 1.711,02
CHO 55,64 56,34 38,39 67,92 1.799,63 1.825,71 1.211,88 2.179,83
CKA 67,14 75,78 65,47 90,85 1.801,05 2.193,25 1.781,37 2.453,97
CFB 315,47 332,34 694,15 763,37 10.815,86 11.443,23 25.742,71 27.951,25
CPC 1.088,72 1.097,56 2.263,18 2.514,67 37.420,05 37.741,75 84.815,86 92.807,17
CBT 2.226,07 2.226,07 3.164,91 3.878,13 75.757,64 75.757,64 119.698,1 141.934,76
163
Capítulo IV: Presentación y Análisis de Resultados
(1949) con 17,64%. Igualmente se puede notar que el peor desempeño lo obtuvo la
correlación de Baxendell y Thomas (1961) con 2.226,07% de error absoluto
promedio. Es notable que, cuando se comparan contra toda la base de datos, 7 de los
10 modelos con el mejor desempeño son correlaciones empíricas. Asimismo, resalta
el hecho de que 3 de los 4 modelos mecanicistas estudiados se encuentren entre los 10
con mejor comportamiento frente a toda la base de datos. Finalmente, vale la pena
destacar que más de la mitad de las correlaciones o modelos estudiados obtuvieron
errores absolutos promedio por debajo de 40%.
164
Capítulo IV: Presentación y Análisis de Resultados
7.000
CGR
MGO
CLM
6.000
MMO
CGA
5.000 CGA
CHB
Presión calculada [Pa\m]
CAZ
4.000 CEA
CDR
-25%
3.000 +25%
Identidad
2.000
1.000
0
0 1.000 2.000 3.000 4.000 5.000 6.000 7.000
Fig. 33: Gradiente de presión calculado vs. experimental de los 10 modelos o correlaciones con mejor
desempeño para toda la base de datos experimentales.
165
Capítulo IV: Presentación y Análisis de Resultados
Tabla 13: Comparación de la precisión de los modelos o correlaciones contra los 33 datos experimentales
que presentan patrón de flujo tapón (SL).
Parámetros Estadísticos
Modelo o
E1 E2 E3 E4 E5 E6 E7 E8
Correlación
[%] [%] [%] [%] [Pa/m] [Pa/m] [Pa/m] [Pa/m]
CGW -9,11 11,78 13,51 16,38 -309,49 404,78 447,02 546,45
CGR 13,94 15,46 13,55 19,59 435,08 494,40 413,47 605,12
CDR 12,90 15,79 14,97 19,89 408,22 509,32 467,96 625,17
CGA 8,41 16,41 18,00 19,92 197,70 537,49 589,42 622,67
COR -7,43 17,47 19,30 20,72 -343,12 658,79 724,58 804,01
MMO 20,84 21,01 14,25 25,52 671,54 679,44 409,07 795,23
COL -20,86 21,09 12,68 24,69 -730,15 737,98 457,06 871,03
MGO 21,47 21,60 14,77 26,34 688,10 694,47 423,28 816,97
CAZ 22,84 22,91 14,09 27,14 739,29 743,05 403,47 852,30
CHB 0,06 23,58 28,58 28,58 -148,23 840,85 1.013,40 1.024,51
CDU -24,22 24,25 13,24 27,94 -895,96 896,69 563,32 1.070,12
CLM 24,53 25,41 17,83 30,64 751,14 796,83 509,16 917,11
MTB 26,27 26,35 16,83 31,54 845,65 850,33 467,69 977,86
CCHi 27,23 27,29 15,88 31,89 865,00 868,49 414,50 971,30
MAN 29,82 30,87 21,19 36,96 997,88 1.029,68 690,70 1.226,35
CBB 32,71 32,71 17,59 37,59 1.071,32 1.071,32 477,03 1.187,92
CCH -33,08 33,08 12,44 35,82 -1.171,03 1.171,03 494,72 1.287,99
MHC -35,06 35,06 10,74 37,19 -1.240,48 1.240,48 460,55 1.341,26
MHMc -35,95 35,95 10,74 38,05 -1.266,99 1.266,99 449,81 1.362,99
CBW -36,01 36,01 11,00 38,19 -1.268,29 1.268,29 453,45 1.365,45
CMH -36,37 36,37 10,80 38,48 -1.280,76 1.280,76 448,75 1.375,86
CRE -36,56 36,56 10,97 38,71 -1.286,08 1.286,08 449,72 1.381,28
MHW -37,24 37,24 11,06 39,40 -1.306,95 1.306,95 444,15 1.399,56
CKA 16,74 37,30 41,19 41,69 280,57 1.213,74 1.401,89 1.338,27
CEA 39,86 39,86 17,90 44,26 1.347,52 1.347,52 557,51 1.477,62
MHO 54,58 54,58 22,50 59,82 1.774,80 1.774,80 452,14 1.858,16
CHO 85,64 85,64 27,12 91,10 2.863,86 2.863,86 524,03 2.955,10
CFB 107,90 130,91 321,36 339,53 3.215,89 4.168,76 9.814,28 10.343,37
CPC 397,91 409,98 955,43 1.037,37 12.181,51 12.664,22 28.929,84 31.463,66
CBT 1.411,68 1.411,68 1.350,21 1.969,31 47.610,44 47.610,44 40.400,61 63.006,34
En la Tabla 13, donde se evalúa la precisión de todos los modelos frente a los
datos que presentan un patrón de flujo tapón, se observa que el mejor desempeño lo
tiene la correlación de Griffith y Wallis (1961) con un error absoluto promedio de
166
Capítulo IV: Presentación y Análisis de Resultados
7.000
CGW
CGR
CDR
6.000
CGA
COR
5.000 MMO
COL
MGO
4.000 CAZ
Presión calculada [Pa\m]
CHB
-25%
3.000 +25%
Identidad
2.000
1.000
0
0 1.000 2.000 3.000 4.000 5.000 6.000 7.000
Fig. 34: Gradiente de presión calculado vs. experimental de los 10 modelos o correlaciones con
mejor desempeño en flujo tipo tapón.
167
Capítulo IV: Presentación y Análisis de Resultados
Tabla 14: Comparación de la precisión de los modelos o correlaciones contra los 28 datos experimentales
que presentan patrón de flujo espuma (FR).
Parámetros Estadísticos
Modelo o
E1 E2 E3 E4 E5 E6 E7 E8
Correlación
[%] [%] [%] [%] [Pa/m] [Pa/m] [Pa/m] [Pa/m]
CLM 8,37 10,87 17,00 19,01 247,83 331,03 530,23 587,24
MGO -3,17 11,12 18,81 19,09 -108,61 352,18 593,48 603,69
CBB -0,92 11,74 22,12 22,14 -33,72 366,85 696,49 697,34
CGR -5,96 13,57 18,61 19,58 -199,29 428,16 587,05 621,14
MMO -6,24 13,77 20,42 21,38 -209,17 442,88 657,04 690,71
CEA -0,04 14,12 29,23 29,23 0,61 433,83 912,54 912,54
CAZ -11,55 18,92 22,25 25,16 -335,86 566,44 668,21 750,65
CGA 7,52 19,93 25,01 26,16 231,94 628,95 805,55 839,46
CHB -2,83 22,90 29,28 29,42 -55,81 740,18 979,85 981,50
MAN -17,59 25,29 23,26 29,36 -541,31 781,25 732,09 916,42
CCHi 26,87 27,01 26,22 37,89 797,73 802,54 759,55 1.112,14
CDR -22,88 29,97 24,54 33,84 -699,91 921,52 764,88 1.045,50
CDU -35,05 36,98 18,53 40,22 -1.085,35 1.145,85 594,38 1.254,95
MHO 41,03 41,03 19,46 46,09 1.225,62 1.225,62 527,30 1.354,92
CHO 41,08 41,08 26,84 49,70 1.266,92 1.266,92 847,77 1.543,77
COR 9,09 41,42 58,38 59,11 376,67 1.346,85 1.989,26 2.025,91
COL -50,42 54,14 21,00 55,48 -1.551,26 1.668,04 663,88 1.713,56
MHC -52,91 54,43 16,79 56,43 -1.630,81 1.678,56 543,21 1.747,32
CCH -53,12 54,75 20,39 57,81 -1.634,93 1.685,98 641,81 1.784,35
CMH -57,26 58,75 16,43 60,58 -1.766,27 1.813,13 539,13 1.877,74
MHMc -57,68 59,20 16,47 61,01 -1.778,05 1.825,67 537,22 1.888,69
CBW -58,37 59,89 16,57 61,71 -1.799,48 1.847,21 540,91 1.910,67
CRE -59,21 60,70 16,65 62,55 -1.824,36 1.871,20 540,82 1.934,95
MHW -61,06 62,54 16,81 64,41 -1.883,21 1.929,66 554,19 1.996,24
CKA 81,63 81,63 38,16 89,81 2.411,25 2.411,25 798,84 2.535,48
CFB 341,90 356,84 719,38 799,21 11.703,71 12.168,86 25.533,74 28.178,39
CPC 1.175,06 1.183,36 2.247,57 2.546,27 39.937,70 40.196,20 79.943,89 89.694,60
CBT 2.200,85 2.200,85 3.122,94 3.843,94 73.073,80 73.073,80 111.936,58 134.414,83
168
Capítulo IV: Presentación y Análisis de Resultados
empíricas y que todas ellas presentan errores absolutos promedio por debajo de 23%.
Asimismo, también es resaltante que todos los modelos mecanicistas que fueron
evaluados contra esta serie de puntos se encuentran entre los 10 con mejor desempeño
y que sus errores absolutos promedio están por debajo de 26%. Igualmente, en la Fig.
35 se observa gráficamente la comparación de los resultados obtenidos contra la data
experimental de los 10 primeros modelos en este patrón de flujo.
7.000
CLM
MGO
CBB
6.000
CGR
MMO
5.000 CEA
CAZ
CGA
4.000 CHB
Presión calculada [Pa\m]
MAN
-25%
3.000 +25%
Identidad
2.000
1.000
0
0 1.000 2.000 3.000 4.000 5.000 6.000 7.000
Fig. 35: Gradiente de presión calculado vs. experimental de los 10 modelos o correlaciones con
mejor desempeño en flujo tipo espuma.
169
Capítulo IV: Presentación y Análisis de Resultados
Tabla 15: Comparación de la precisión de los modelos o correlaciones contra los 11 datos experimentales
que presentan patrón de flujo ondulado (RI).
Parámetros estadísticos
Modelo o
E1 E2 E3 E4 E5 E6 E7 E8
Correlación
[%] [%] [%] [%] [Pa/m] [Pa/m] [Pa/m] [Pa/m]
CGR 5,48 10,68 13,12 14,32 103,47 265,85 330,46 347,82
CEA -9,48 11,16 10,30 14,31 -261,28 297,37 270,47 385,03
CLM 12,72 13,47 9,38 16,31 314,96 335,73 210,42 391,65
CDU -11,30 14,08 18,38 21,87 -314,68 375,30 497,59 597,09
MGO 7,46 14,91 20,17 21,63 135,66 365,32 480,50 501,12
CHO 14,57 15,79 13,66 20,50 361,22 394,46 343,70 511,52
MMO -16,06 16,37 22,57 28,17 -468,19 475,47 648,39 813,35
CBB -15,21 16,43 10,63 19,17 -416,43 446,17 314,79 538,38
CHB 0,78 16,48 22,70 22,71 77,40 488,45 721,62 726,17
MAN -26,79 26,79 5,40 28,61 -718,49 718,49 232,62 788,65
CGA 22,89 31,04 33,36 41,10 645,75 872,77 993,09 1.202,05
CDR -40,49 40,49 6,53 42,97 -1.068,44 1.068,44 257,98 1.149,90
CAZ -42,74 42,74 14,10 46,99 -1.078,71 1.078,71 288,12 1.167,47
MHC -45,30 45,30 11,23 48,82 -1.187,42 1.187,42 333,83 1.289,35
COR 6,41 47,36 68,74 69,07 395,45 1.420,58 2.327,57 2.364,24
CCH -45,89 53,21 33,42 58,59 -1.223,09 1.402,98 892,54 1.562,74
CMH -54,11 54,11 7,52 57,25 -1.424,91 1.424,91 302,91 1.524,85
COL -54,99 54,99 5,57 57,94 -1.444,59 1.444,59 256,34 1.536,63
MHMc -57,61 57,61 6,56 60,78 -1.514,57 1.514,57 285,31 1.613,91
MHO 57,83 57,83 24,49 65,42 1.428,93 1.428,93 450,30 1.564,86
CBW -59,29 59,29 5,53 62,43 -1.559,40 1.559,40 281,90 1.659,63
CRE -61,54 61,54 5,39 64,77 -1.615,45 1.615,45 268,05 1.715,37
MHW -63,81 63,81 4,14 67,05 -1.680,60 1.680,60 301,77 1.788,27
CCHi 92,73 92,73 46,28 107,71 2.372,57 2.372,57 1.087,89 2.715,79
CKA 163,64 163,64 49,27 178,56 4.122,39 4.122,39 749,36 4.388,05
CFB 443,04 449,42 721,20 857,93 13.979,40 14.157,27 27.090,27 30.803,39
CPC 1.451,01 1.453,09 2.366,81 2.813,85 46.062,56 46.120,58 88.561,65 100.881,62
CBT 2.611,55 2.611,55 3.347,65 4.325,39 79.793,50 79.793,50 126.939,33 152.043,73
170
Capítulo IV: Presentación y Análisis de Resultados
empíricas, todas ellas con errores absolutos promedio por debajo de 17%. Igualmente,
es resaltante que todos los modelos mecanicistas que fueron evaluados contra esta
serie de datos se encuentran entre los 10 con mejor desempeño y que sus errores
absolutos promedio están por debajo de 27%. En la Fig. 36 puede verse la
distribución de las predicciones de los 10 primeros modelos en este patrón de flujo.
7.000
MGO
CGR
CEA
6.000
CLM
CDU
5.000 CHO
Presión calculada [Pa\m]
MMO
CBB
4.000 CHB
MAN
-25%
3.000 25%
Identidad
2.000
1.000
0
0 1.000 2.000 3.000 4.000 5.000 6.000 7.000
Fig. 36: Gradiente de presión calculado vs. experimental de los 10 modelos o correlaciones con
mejor desempeño en flujo tipo ondulado.
171
Capítulo IV: Presentación y Análisis de Resultados
172
Capítulo IV: Presentación y Análisis de Resultados
CHB
CAZ
MGO
COL
MMO
COR
CGA
CDR
CGR
CGW
0 10 20 30 40 50 60 70 80 90 100
% de datos experimentales sobrepredichos % de Sobreprediccion
En la Fig. 37 se observa que para los datos que presentaron un patrón de flujo
tipo tapón, la correlación de Aziz et al. (1972), el modelo mecanicista de Gómez et
al. (2000) y el modelo mecanicista de Ouyang (1999) sobrepredijeron el gradiente de
presión en más del 90% de los datos disponibles, aunque esa sobrepredicción se
encuentra en promedio por 23%. En cambio, si bien la correlación de Hagedorn y
Brown (1965) sólo sobrepredijo los resultados en el 42,4% de los datos, su
sobrepredicción promedio es de 27,9%, que es la mayor de este grupo de modelos.
173
Capítulo IV: Presentación y Análisis de Resultados
MAN
CHB
CGA
CAZ
CEA
MMO
CGR
CBB
MGO
CLM
0 10 20 30 40 50 60 70 80 90 100
% de datos sobrepredichos % Sobreprediccion
174
Capítulo IV: Presentación y Análisis de Resultados
MAN
CHB
CBB
MMO
CHO
MGO
CDU
CLM
CEA
CGR
0 10 20 30 40 50 60 70 80 90 100
% de datos sobrepredichos % de Sobreprediccion
175
Capítulo IV: Presentación y Análisis de Resultados
Tabla 16: Comparación porcentual de la caída de presión por fricción, elevación y aceleración para cada
modelo.
176
Capítulo IV: Presentación y Análisis de Resultados
SÍNTESIS Y CONCLUSIONES
• De los datos extraídos de Govier y Leigh (1958) y Brown et al. (1960) se observó
que la mayoría eran para flujo tipo tapón (44,6 % de la base de datos), luego para
el patrón tipo espuma (37,84% de la base de datos) y, en tercer lugar, para el
régimen ondulado (14,86% de la base de datos). Los treinta modelos y
correlaciones fueron evaluados para la base de datos completa y por patrón de
flujo, obteniendo así cuatro casos de estudio.
177
Capítulo IV: Presentación y Análisis de Resultados
178
Capítulo IV: Presentación y Análisis de Resultados
RECOMENDACIONES
• En este trabajo se utilizó una base de datos compuesta por sólo 74 puntos, de los
cuales la mayoría presentan flujo tapón. Se recomienda el uso de una base de
datos más amplia y con mayor variedad de patrones de flujo.
• También se recomienda hacer un estudio similar a éste, pero utilizando una base
de datos con diferentes condiciones de operación como podrían ser mayores
diámetros de tubería, altas viscosidades o flujo inclinado. Otra variante de este
estudio que se podría desarrollar sería la comparación de modelos y correlaciones
de flujo bifásico, pero en espacios anulares.
179
Referencias Bibliográficas
REFERENCIAS BIBLIOGRÁFICAS
Alves, I. N. et al. (1988, Nov. 27 – Dec. 2). Modeling Annular Flow Behavior for
Gas Wells. Paper presentado en el Winter Annual Meeting of ASME,
Chicago, Illinois, EUA
Ansari, A., Sylvester, N., Sarica, C., Shoham, O., Brill, J. (1994). A Comprehensive
Mechanistic Model for Upward Two-Phase Flow in Wellbores. SPE
Production & Facilities J., (Vol. 9, No. 2), 142-152.
Aziz, K., Govier, G. W., Fogarasi, M. (1972). Pressure Drop in Wells Producing Oil
and Gas. Journal of Canadian Pet. Tech., (S.V./S.N.), 38-48
Baker, O. (1954). Design of Pipelines for the Simultaneous Flow of Oil and Gas. Oil
and Gas Journal, (Vol. 53, S.N.), 185
Bankoff, S. (1960). Variable Density Single Fluid Model for Two-phase Flow with
Particular Reference to Steam-Water Flow. Trans. ASME (Vol. 82, S.N.), 265
Barnea, D. (1987). A Unified Model for Predicting Flow Pattern Transition for the
Whole Range of Pipe Inclinations. Intl. Journal Multiphase Flow (Vol. 13, No.
1), 1-12
Barnea, D., Shoham, O., Taitel, Y. (1982). Flow Pattern Transition for Vertical
Downward Two-Phase Flow. Chem. Eng. Sci. (Vol. 37, No. 5), 741-744
180
Referencias Bibliográficas
Barnea, D., Tailtel, Y., Dukler, A.E. (1980). Modelling Flow Pattern Transitions for
Steady Upward Gas-Liquid Flow in Vertical Tubes, AIChE Journal (Vol. 26,
No. 3), 345-354.
Beattie, D., Whalley, P.D. (1982). A Simple Two-Phase Frictional Pressure Drop
Calculation Method, Int. J. Multiphase Flow (Vol. 8, No. 1), 83-87
Beggs, H., Brill, J. (1973). A Study of Two-Phase Flow in Inclined Pipes, J. Pet.
Tech. (Vol. 25, No. 5), 607-617.
Brill, J., Beggs, H. (1988). Two Phase Flow in Pipes. 6th edition. Copyright 1978 by
J. Brill and H. Beggs.
Brown, R. A. S., Sullivan, G. A., Govier, G. W. (1960). The Upward Vertical Flow of
Air-Water Mixtures: III. Effect of Gas Phase Density on Flow Pattern, Holdup
and Pressure Drop. The Canadian Journal of Chemical Engineering,
(S.V./S.N.), 62-66
Chen, I., Yang, K., Wang, C. (2002). An Empirical Correlation for Two-phase
Frictional Performance in Small Diameter Tubes. Int. Journal of Heat Mass
Transfer (Vol. 45, No. 17), 3667-3671
Chenoweth, J. M., Martin, N. W. (1955). Petroleum Refiner (Vol. 34, S.N.), 151
181
Referencias Bibliográficas
Chisholm, D. (1967). A Theoretical Basis for the Lockhart Martinelli Correlation for
Two-phase Flow. Int. Journal of Heat and Mass Transfer (Vol. 10, No. 12),
1767-1778
Cicchitti, A, Lombardi, C., Silvestri, M., Soldaini, G., Zavattareui, R., (1960). Two-
Phase Cooling Experiments-Pressure Drop, Heat Transfer and Burnout
Measurements, Energi Nucl. (Vol. 7, No. 6), 407-425.
CTES, L.C. (1998). Multiphase Flow Models Range of Applicability. Conroe, Texas,
EUA: Bharath Rao.
Davis, R. M., Taylor, G. I. (1950). The Mechanics of Large Bubbles Rising through
Extended Liquids and through Liquids in Tubes. Proc. Of the Royal Society,
London, England (Vol. 200, A) 375-390
Duns, H. Jr., Ros, N. C. J. (1963). Vertical Flow of Gas and Liquid Mixtures in
Wells. Proc., 6th. Worl Pet. Congress, 451
182
Referencias Bibliográficas
Dukler, A. E., Wicks III, M., Cleveland, R. (1964). Frictional Pressure Drop in Two-
Phase Flow: A. A Comparision of Existing Correlations for Pressure Loss
and Holdup, AIChE Journal (Vol. 10, No. 1), 38-43.
Dukler, A. E., Wicks III, M., Cleveland, R. (1964). Frictional Pressure Drop in Two-
Phase Flow: B. An Approach through Similarity Analysis, AIChE Journal
(Vol. 10, No. 1), 44-51.
Eaton, B., Andrews, D., Knowles, C., Silberberg, I., Brown, K. (1967). The
Prediction of Flow Patterns, Liquid Holdup and Pressure Losses Occurring
During Continuous Two-phase Flow in Horizontal Pipelines. Trans. AIME,
(S.V./S.N.), 815
Fernandes, R. C., Semait, T., Dukler, A. E. (1986). Hydrodynamic Model for Gas-
Liquid Slug Flow in Vertical Tubes. AIChE Journal, (Vol. 29, S.N.) 981-989
García, F., García R., Joshep, D. (2005). Composite Power Law Holdup Correlations
in Horizontal Pipes. Int. Journal of Multiphase Flow (Vol. 31, No. 12), 1276-
1303
García, F., García, R., Padrino, J. C., Mata, C., Trallero, J., Joseph, D.(2003). Power
Law and Composite Power Law Friction Factor Correlations for Laminar and
183
Referencias Bibliográficas
García, F., García, J. M., García, R., Joseph, D. (2006). Friction Factor Improved
Correlations for Laminar and Turbulent Gas-Liquid Flow in Horizontal
Pipelines. In press.
Griffith, P., Lau, C. W., Hon, P. C., Pearson, J. F. (1973). Two-phase Pressure Drop
in Inclined and Vertical Pipes. Tech. Report No. 80063-81, Heat Transfer
Laboratory, Dept. of M.E., MIT, Cambridge, Massachusetts, EUA
Griffith, P., Wallis, G. B. (1961). Two-phase Slug Flow. Journal of Heat Transfer,
Trans. ASME, (Vol.83, S.N.), 307-320
Gómez, L., Shoham, O., Schimdt, Z., Chokshi, R., Northug, T. (2000a). Unified
Mechanistic Model for Steady-State Two-Phase Flow: Horizontal to Vertical
Upward Flow. SPE 65705, SPE Journal (Vol. 5, No. 3) 339-350
Gómez, L., Shoham, O., Taitel, Y. (2000b). Prediction of Slug Liquid Holdup:
Horizontal to Upward Vertical Flow. Int. Journal of Multiphase flow (Vol. 26
No. 3), 517-521
Govier, G. W., Leigh, W. (1958). The Upward Vertical Flow of Air-Water Mixtures:
II. Effect of Tubing Diameter on Flow Pattern, Holdup and Pressure Drop.
The Canadian Journal of Chemical Engineering, (S.V./S.N.), 195-202
184
Referencias Bibliográficas
Hoogendorn, C. (1959). Gas-Liquid Flow in Horizontal Pipes, Chem. Eng. Sci. (Vol.
9, No. 4), 205-217
Kadambi, V. (1981). Void Fraction and Pressure Drop in Two-phase Stratified Flow,
The Canadian Journal of Chemical Engineering (Vol. 59, S. N.), 584-589.
Lockhart, R., Martinelli, R. (1949). Proposed Correlation of Data for Isothermal Two-
Phase Two Component Flow in Pipes, Chem. Eng. Prog. (Vol. 45, No. 1), 39-
48.
185
Referencias Bibliográficas
Mishima, K., Ishii, M. (1984). Flow Regime Transition Criteria for Upward Two-
phase Flow in Vertical Tubes. Int. Journal of Heat and Mass Transfer (Vol.
27, No 5), 723-737
Ouyang, L. (1998). Single Phase and Multiphase Fluid Flow in Horizontal Wells. Phd
Dissertation Thesis. Department of Petroleum Engineering. School of Earth
Sciences. Standford University. Standford, California, EUA, 248
186
Referencias Bibliográficas
Ouyang, L., Aziz, K. (1999, Oct. 3-6). A Mechanistic Model for Gas-Liquid Flow in
Pipes with Radial Influx or Outflux. Paper SPE 56525 presentado en el
Annual Technical Conference and Exhibition, Houston, Texas, EUA
Poettman, F. H., Carpenter, P. G. (1952). The Multiphase Flow of Gas, Oil and Water
Through Vertical Flow Strings with Application to the Design of Gas-Lift
Installations, Drilling and Production Practice API, 257-317
Reid, R., Reynolds, A., Diglio, A., Spiewak, I., Klipstein, D. (1957). Two-Phase
Pressure Drops in Large-Diameter Pipes, AIChE J. (Vol. 3, No. 3), 321-324.
Shoham, O., (2001). Two-Phase Flow Modeling, The University of Tulsa, Oklahoma,
EUA
Taitel, Y., Barnea, D. (1990). A Consistent Approach for Calculating Pressure Drop
in Inclined Slug Flow, Chemical Engineering Science (Vol. 45, No. 5), 1199-
1206.
Taitel, Y., Barnea, D., Dukler, A. E. (1980). Modelling Flow Pattern Transitions for
Steady Upward Gas-Liquid Flow in Vertical Tubes. AIChE Journal (Vol. 26,
S. N.), 345-354
187
Referencias Bibliográficas
Taitel, Y., Barnea, D., (1990). Two-phase Slug Flow, Advances in Heat Transfer
(Vol. 20, S. N.), 83-132.
Wallis, G. B. (1969). One-Dimensional Two-Phase Flow. New York City, New York,
EUA. McGraw-Hill Book Co. Inc., 409
Xiao, J., Shoham, O., Brill, J. (1990, Sept. 23-26). A Comprehensive Mechanistic
Model for Two-phase Flow in Pipelines. The 65th SPE Annual Technical
Conference and Exhibition, New Orleans, Louisiana, EUA. Paper SPE 20631,
167-180
188
Apéndices
APÉNDICES
Muchos de los modelos incluidos en este trabajo requerían del uso de gráficas
o figuras para la determinación de parámetros necesarios para calcular la caída de
presión. Es por ello que, con la finalidad de hacer automáticos estos procedimientos,
se hizo una digitalización de las imágenes necesarias para poder luego aproximarlas
con una ecuación, de tal manera que dicha fórmula de ajuste se pudiera incluir en los
códigos de los programas.
189
Apéndices
Fig. 40: Comparación entre la curva original del factor K de Griffith et al. y la ecuación de ajuste,
Griffith (1973).
190
Apéndices
Fig. 41: Comparación entre la curva original del factor L2 de Duns y Ros y la ecuación de ajuste,
Beggs y Brill (1988).
191
Apéndices
Tabla 17: Funciones de ajuste de datos experimentales de la correlación de Fancher y Brown (1963).
Modelo o
ECUACIÓN DE AJUSTE E2 E2min E2max Ec.
Correlación
⎛ b ⎞
CFB ⎜ a + +c ln( n ) ⎟ -0,350 -1,740 1,700 Ec. 325
f = e⎝ n ⎠
, donde : n = ρUDi
n
CFB f = , donde : n = ρUDi 0,903 -2,018 5,718 Ec. 326
a + bn − cn 2
Tabla 18: Valores de los coeficientes de las funciones de ajuste de la correlación de Fancher y Brown (1963).
Ec. a b c d e f g h i j k
Ec. 324 2,78E-11 -3,19E-9 1,61E-7 -4,74E-6 9,00E-5 -1,16E-3 1,03E-2 -0,06 0,27 -0,73 1,05
192
Apéndices
Tabla 20: Funciones de ajuste de datos experimentales de la correlación de Lockhart y Martinelli (1949).
x
CLM Hl = 2,952 -7,313 16,275 Ec. 327
a + bx + c x
Tabla 21: Valores de los coeficientes de las funciones de ajuste de la correlación de Lockhart y Martinelli (1949).
Ec. a b c d e f g h i j k
Ec. 327 7,19E-01 8,18E-01 2,93 - - - - - - - -
193
Apéndices
Tabla 23: Funciones de ajuste de datos experimentales de la correlación de Hagedorn y Brown (1965).
Modelo o
ECUACIÓN DE AJUSTE E2 E2min E2max Ec.
Correlación
H ff = C 2 ⋅ e b−c2c
a
CHB 1,208 -11,601 37,610 Ec. 328
CHB ψ = aC310 + bC39 + cC38 + dC37 + eC36 + fC35 + gC34 + hC33 + iC32 + jC3 + k 0,127 -1,402 2,197 Ec. 330
194
Apéndices
Tabla 24: Valores de los coeficientes de las funciones de ajuste de la correlación de Hagedorn y Brown (1965).
Ec. a b c d e f g h i j k
Ec. 328 4,18E-1 2,77 1,22E02 - - - - - - - -
Ec. 329 -4,87E04 1,16E05 -1,19E05 6,82E04 -2,39E04 5,23E03 -7,01E02 5,34E01 -2,03 7,74E-2 1,74E-3
Ec. 330 -5,18E13 3,50E13 -9,99E12 1,59E12 -1,57E11 1,00E10 -4,14E08 1,08E07 -1,67E05 1,40E03 -3,77
195
Apéndices
Tabla 26: Funciones de ajuste de datos experimentales de la correlación de Griffith y Wallis (1961).
Modelo o
ECUACIÓN DE AJUSTE E2 E2min E2max Ec.
Correlación
CCHi
CGW C = a Re10 + b Re9 + c Re8 + d Re7 + e Re6 + f Re5 + g Re 4 + h Re3 + i Re 2 + j Re + k -0,590 -21,582 11,306 Ec. 331
1 b b b b b b b b b b
COR
Tabla 27: Valores de los coeficientes de las funciones de ajuste de la correlación de Griffith y Wallis (1961).
Ec. a b c d e f g h i j k
Ec. 331 -16,14 14,66E1 -56,60E1 1,22E3 -1,60E3 1,34E3 -71,68E1 24,40E1 -51,33 6,30 -0,01
196
Apéndices
Tabla 29: Funciones de ajuste de datos experimentales de correlación de Duns y Ros (1963).
Modelo o
ECUACIÓN DE AJUSTE E2 E2min E2max Ec.
Correlación
10 9 8 7 6 5 4 3 2
CDR L1 = aN d + bN d + cN d + dN d + eN d + fN d + gN d + hN d + iN d + jN d + k -4,901 -8,085 -3,323 Ec. 332
10 9 8 7 6 5 4 3 2
CDR L 2 = aN d + bN d + cN d + dN d + eN d + fN d + gN d + hN d + iN d + jN d + k 0,002 -0,710 1,234 Ec. 333
10 9 8 7 6 5 4 3 2 Ec. 334
CDR F1 = aN l + bN l + cN l + dN l + eN l + fN l + gN l + hN l + iN l + jN l + k 0,025 -0,726 0,961
Ec. 335
10 9 8 7 6 5 4 3 2 Ec. 336
CDR F2 = aN l + bN l + cN l + dN l + eN l + fN l + gN l + hN l + iN l + jN l + k -0,011 -1,109 1,909
Ec. 337
10 9 8 7 6 5 4 3 2
CDR F3 = aN l + bN l + cN l + dN l + eN l + fN l + gN l + hN l + iN l + jN l + k 0,026 -0,712 0,541 Ec. 338
Nl
CDR F3 = 2,030 -6,298 7,902 Ec. 339
a + bN l + c N l
197
Apéndices
10 9 8 7 6 5 4 3 2
CDR F5 = aN l + bN l + cN l + dN l + eN l + fN l + gN l + hN l + iN l + jN l + k -0,219 -11,084 14,593 Ec. 343
F7 = aN l e − cN l
b
CDR -0,112 -5,706 5,345 Ec. 344
Tabla 30: Valores de los coeficientes de las funciones de ajuste de la correlación de Duns y Ros (1963).
Ec. a b c d e f g h i j k
Ec. 332 -5,34E-16 2,22E-13 -3,61E-11 2,53E-09 6,01E-09 -1,47E-05 1,19E-03 -4,82E-02 1,10 -13,51 70,57
Ec. 333 -2,20E-16 1,00E-13 -1,96E-11 2,17E-09 -1,48E-07 6,39E-06 -1,75E-04 2,87E-03 -2,51E-02 9,43E-02 0,30
Ec. 334 -8,54E16 2,62E16 -3,54E15 2,75E14 -1,37E13 4,52E11 -1,00E10 1,47E08 -1,37E06 7,24E03 -15,39
Ec. 335 -11,84 121,21 -535,91 134,15E01 -209,41E01 211,92E01 -140,28E01 600,04 -159,32 22,35 0,91
Ec. 336 -13,68E12 77,78E11 -19,24E11 27,17E10 -24,12E09 13,67E08 -5,27E07 12,55E05 -17,20E03 120,28 -8,41E-02
Ec. 337 -2,19 22,50 -100,97 259,44 -421,30 450,54 -320,79 150,18 -44,16 6,98 6,03E-1
Ec. 338 -1,18E24 6,85E22 -1,74E21 2,54E19 -2,36E17 1,46E15 -6,10E12 1,69E10 -2,96E07 2,97E04 -12,10
Ec. 340 6,57E-14 -1,72E-11 1,93E-09 -1,22E-07 4,80E-06 -1,21E-04 1,98E-03 -2,06E-02 0,13 -0,47 1,18
Ec. 343 -0,76 7,07 -26,73 51,02 -46,09 2,04 36,01 -34,31 14,31 -2,68 0,22
198
Apéndices
199
Apéndices
Tabla 33: Valores de los coeficientes de las funciones de ajuste de la correlación de Eaton et al.(1967).
Ec. a b c d e f g h i j k
Ec. 345 -4,87E4 1,16E5 -1,19E5 6,82E4 -2,39E4 5,23E3 -7,01E2 5,34E1 -2,03 7,74E-2 1,74E-3
200
Apéndices
10 9 8 7 6 5 4 3 2
CGR S = aRl + bRl + cRl + dRl + eRl + fRl + gRl + hRl + iRl + jRl + k , para : x = 1,5% 0,011 -3,805 2,297 Ec. 347
10 9 8 7 6 5 4 3 2
CGR S = aRl + bRl + cRl + dRl + eRl + fRl + gRl + hRl + iRl + jRl + k , para : x = 2% 0,009 -3,957 2,838 Ec. 348
10 9 8 7 6 5 4 3 2
CGR S = aRl + bRl + cRl + dRl + eRl + fRl + gRl + hRl + iRl + jRl + k , para : x = 3% -0,028 -3,582 2,315 Ec. 349
10 9 8 7 6 5 4 3 2
CGR S = aRl + bRl + cRl + dRl + eRl + fRl + gRl + hRl + iRl + jRl + k , para : x = 4% -0,020 -4,940 3,637 Ec. 350
10 9 8 7 6 5 4 3 2
CGR S = aRl + bRl + cRl + dRl + eRl + fRl + gRl + hRl + iRl + jRl + k , para : x = 6% -0,005 -4,369 2,009 Ec. 351
10 9 8 7 6 5 4 3 2
CGR S = aRl + bRl + cRl + dRl + eRl + fRl + gRl + hRl + iRl + jRl + k , para : x = 8% 0,008 -5,603 1,777 Ec. 352
10 9 8 7 6 5 4 3 2
CGR S = aRl + bRl + cRl + dRl + eRl + fRl + gRl + hRl + iRl + jRl + k , para : x = 10% -0,089 -6,614 3,023 Ec. 353
9 8 7 6 5 4 3 2
CGR S = aRl + bRl + cRl + dRl + eRl + fRl + gRl + hRl + iRl + j , para : x = 15% 0,271 -2,656 6,468 Ec. 354
10 9 8 7 6 5 4 3 2
CGR S = aRl + bRl + cRl + dRl + eRl + fRl + gRl + hRl + iRl + jRl + k , para : x = 20% 0,095 -2,078 2,487 Ec. 355
201
Apéndices
10 9 8 7 6 5 4 3 2
CGR S = aRl + bRl + cRl + dRl + eRl + fRl + gRl + hRl + iRl + jRl + k , para : x = 30% 0,209 -2,804 4,365 Ec. 356
10 9 8 7 6 5 4 3 2
CGR S = aRl + bRl + cRl + dRl + eRl + fRl + gRl + hRl + iRl + jRl + k , para : x = 40% -0,087 -6,330 2,175 Ec. 357
7 6 5 4 3 2
CGR S = aRl + bRl + cRl + dRl + eRl + fRl + gRl + h, para : x = 50% 0,697 -2,856 6,924 Ec. 358
10 9 8 7 6 5 4 3 2
CGR S = aRl + bRl + cRl + dRl + eRl + fRl + gRl + hRl + iRl + jRl + k , para : x = 60% -0,279 -8,576 4,376 Ec. 359
10 9 8 7 6 5 4 3 2
CGR S = aRl + bRl + cRl + dRl + eRl + fRl + gRl + hRl + iRl + jRl + k , para : x = 80% -0,583 -22,538 4,242 Ec. 360
6 5 4 3 2
CGR S = aRl + bRl + cRl + dRl + eRl + fRl + g , para : x = 100% -0,216 -15,282 4,610 Ec. 361
c
CGR Rv = a ⋅ b Rl ⋅ Rl 0,024 -4,090 5,615 Ec. 363
Tabla 36: Valores de los coeficientes de las funciones de ajuste de la correlación de Griffith et al.(1973).
Ec. a b c d e f g h i j K
Ec. 346 -5,38E-26 2,08E-22 -3,43E-19 3,16E-16 -1,79E-13 6,48E-11 -1,53E-8 2,40E-6 -2,61E-4 2,82E-2 8,70E-1
Ec. 347 -9,71E-26 3,47E-22 -5,23E-19 4,30E-16 -2,08E-13 5,97E-11 -9,52E-9 7,72E-7 -7,11E-5 2,84E-2 8,65E-1
Ec. 348 1,37E-25 -5,35E-22 8,86E-19 -8,04E-16 4,35E-13 -1,43E-10 2,77E-8 -2,89E-6 8,41E-5 3,30E-2 8,61E-1
Ec. 349 6,15E-26 -2,22E-22 3,31E-19 -2,65E-16 1,24E-13 -3,47E-11 6,23E-9 -7,87E-7 5,40E-6 4,69E-2 7,48E-1
Ec. 350 -3,81E-25 1,47E-21 -2,41E-18 2,19E-15 -1,19E-12 4,00E-10 -8,03E-8 8,84E-6 -5,20E-4 7,09E-2 6,95E-1
Ec. 351 -2,79E-25 1,15E-21 -2,01E-18 1,95E-15 -1,14E-12 4,08E-10 -8,65E-8 1,01E-5 -6,76E-4 1,04E-1 5,25E-1
202
Apéndices
Ec. 352 -1,73E-25 7,38E-22 -1,34E-18 1,35E-15 -8,14E-13 3,01E-10 -6,65E-8 8,35E-6 -7,23E-4 1,40E-1 3,48E-1
Ec. 353 -1,28E-24 4,76E-21 -7,51E-18 6,52E-15 -3,40E-12 1,09E-9 -2,08E-7 2,25E-5 -1,47E-3 1,85E-1 1,19E-1
Ec. 354 2,16E-22 -7,87E-19 1,19E-15 -9,73E-13 4,55E-10 -1,21E-7 1,69E-5 -1,25E-3 2,32E-1 5,65E-2 -
Ec. 355 6,55E-26 1,00E-22 -8,27E-19 1,44E-15 -1,21E-12 5,72E-10 -1,54E-7 2,32E-5 -2,06E-3 3,26E-1 -4,59E-1
Ec. 356 -4,37E-25 1,91E-21 -3,49E-18 3,47E-15 -2,06E-12 7,41E-10 -1,58E-7 1,89E-5 -1,51E-3 4,36E-1 -8,67E-1
Ec. 357 2,36E-24 -8,46E-21 1,26E-17 -1,00E-14 4,57E-12 -1,19E-9 1,60E-7 -6,78E-6 -1,12E-3 5,92E-1 -1,69
Ec. 358 8,02E-18 -2,31E-14 2,91E-11 -2,02E-8 8,24E-6 -2,16E-3 7,52E-1 -2,17 - - -
Ec. 359 2,88E-25 -1,38E-21 2,58E-18 -2,51E-15 1,38E-12 -4,19E-10 5,72E-8 2,59E-6 -2,43E-3 9,23E-1 -3,04
Ec. 360 -3,18E-24 1,16E-20 -1,77E-17 1,49E-14 -7,52E-12 2,39E-9 -4,96E-7 7,02E-5 -7,43E-3 1,35 -5,14
Ec. 362 -7,66E-17 3,61E-14 -7,33E-12 8,37E-10 -5,92E-8 2,68E-6 -7,79E-5 1,43E-3 -1,62E-2 1,25E-1 7,83E-3
203
Apéndices
204
Apéndices
⎛ gDi (ρ l − ρ g ) ⎞
0,5
U tb = 1,2U m + 0,35⎜⎜ ⎟⎟
⎝ ρ l ⎠
U sg
H gls =
0,425 + 2,65U m
Donde:
⎛ ⎛ σ g (ρ − ρ ) ⎞
0 , 25
⎞
Altb = H glsU tb + (1 − H gls )⎜U m − H gls 1,53⎜⎜ ⎟ (1 − H gls )0,5 ⎟
l l g
⎜ ρl 2 ⎟ ⎟
⎝ ⎝ ⎠ ⎠
(
U ltb = 9,916 gDi (1 − H gtb ) )
0 , 5 0.5
205
Apéndices
ρ ls = ρ l H lls + ρ g (1 − H lls )
μ ls = μ l H lls + μ g (1 − H lls )
H ltb
U lls = U tb − (U tb + U ltb )
H lls
⎛ gσ l (ρ l − ρ g ) ⎞
0 , 25
ρ 2 ⎟
⎝ l ⎠
(1 − H lls )
U gtb = U tb − (U tb + U gls )
(1 − H ltb )
• Espesor de la película:
δ 3 g (ρ l − ρ g ) ⎛
2
2δ ⎞
F (δ ) = − ⎜⎜1 − ⎟⎟ U tb + (U tb − U gls )(1 − H lls )
3 μ ⎝ Di ⎠
4 l
206
Apéndices
2
1 ⎛ U U ⎞
Lc = ⎜⎜U tb + gNtb (1 − H lNtb ) − m ⎟⎟
2g ⎝ H lNtb H lNtb ⎠
Donde:
2
⎛ 2δ ⎞
H l Ntb = 1 − ⎜⎜1 − ⎟⎟
⎝ Di ⎠
(1 − H lls )
U l Ntb = U tb − (U tb − U gls )
(1 − H l Ntb )
U sg
Aa = 1 −
U tb
U sg − U gls (1 − H lls )
Bb = Lls
U tb
U sg − U gls (1 − H lls )
Bb = Lls
U tb
U tb − U lls
Cc = H lls
(2 g )0,5
2 2
⎛ − 2 AaBb − 4Cc 2 ⎞ ⎛ Bb ⎞
LtbA = ⎜⎜ ⎟⎟ − 4⎜ ⎟
⎝ Aa 2 ⎠ ⎝ Aa ⎠
207
Apéndices
• Si LtbA > 0 se calculan los valores de longitud LtbB y LtbC para asignarle uno
de estos dos valores a la longitud de la burbuja de Taylor en desarrollo:
2
⎛ 2 AaBb + 4Cc 2 ⎞
⎜⎜ ⎟⎟ + LtbA 0,5
⎝ Aa 2 ⎠
LtbB =
2
2
⎛ 2 AaBb + 4Cc 2 ⎞
⎜⎜ ⎟⎟ − LtbA 0,5
Aa 2 ⎠
LtbC =⎝
2
• Relación de longitudes:
Ltbo
Re l o =
Lsuo
208
Apéndices
Flujo anular
• Velocidad critica:
10000U sg μ g
U cri = 1
⎛ρ ⎞ 2
σ⎜ g ρ ⎟
⎝ l ⎠
• Arrastre de líquido:
Fe = 1 − e −0,125(U cri −1,5 )
• Densidad en el núcleo:
ρ sc = ρ l λlc + ρ g (1 − λlc )
• Viscosidad en el núcleo:
μ sc = μ l λlc + μ g (1 − λlc )
209
Apéndices
1 3 2
⎡ 0,4σ l ⎤ 2
⎛ ρl ⎞ 5⎛ f ⎞ ⎛ U sg ⎞
(U + U sg ) = 0,725 + 4,15⎜
5
2⎢ ⎥ ⎜⎜ ⎟⎟ ⎜ ⎟ 1, 2
⎟
⎣⎢ (ρ l − ρ g )g ⎦⎥
⎜U +U ⎟
⎝σl
sl
⎠ ⎝ 2d ⎠ ⎝ sg sl ⎠
210
Apéndices
5,0E-01
4,5E-01
4,0E-01
3,5E-01
3,0E-01
FF 2,5E-01
2,0E-01
1,5E-01
1,0E-01
5,0E-02
0,0E+00
1,3 1,8 2,3 2,8 3,3 3,8
Vmezcla
211
Apéndices
⎛ gσ l (ρ l − ρ g ) ⎞
0 , 25
U sg
1,53⎜⎜ ⎟ 0,5
⎟ Hl = + 1,2U m
⎝ ρ l
2
⎠ 1− Hl
212
Apéndices
6,0
5,0
4,0
3,0
F(Hl) 2,0
1,0
0,0
0 0,2 0,4 0,6 0,8 1 1,2 1,4 1,6 1,8 2
-1,0
-2,0
Hl
F(Hl) vs Hl Asintota
Al graficar esta función se observó que la misma tenía raíz entre 0 y 1, como
se esperaba, siendo la abscisa correspondiente al valor de la fracción volumétrica de
líquido. Por esto se propuso en este trabajo un lazo para reubicar la semilla de la
iteración. Este aplica cuando el valor de fracción calculado por Newton-Raphson es
mayor a 1 o menor que 0. Cuando se cumplen estas condiciones se define un
diferencial, el cual va a ir disminuyendo a medida que aumenten las iteraciones. El
valor que se asignará a la fracción no será el calculado ya que se sale del intervalo, si
no el anterior incrementado por el diferencial. Este método asegura que el valor de la
fracción nunca va a ser mayor que la unidad además de una rápida convergencia.
213
Apéndices
H gls = 1 − H lls
Ag = A ⋅ H gls
Al = A ⋅ H lls
H ltb =
π
⋅
(D − (D − 2h ) )
i
2
i f
2
4 A
S f = π ⋅ Di
S g = π ⋅ (Di − 2h f )
Si = S f − S g
214
Apéndices
⎛ 2 ⎞
1 ⎜ ⎛ 2h f ⎞ ⎛ 2h f ⎞ ⎛ ⎛ 2h f ⎞ ⎞⎟ ⎟
H ltb = ⋅ ⎜ π − arccos⎜⎜ − 1⎟⎟ + ⎜⎜ − 1⎟⎟ ⋅ ⎜1 − ⎜ − 1⎟⎟ ⎟
π ⎜ ⎝ Di ⎠ ⎝ Di ⎠ ⎜ ⎜⎝ Di ⎠ ⎟⎠ ⎟
⎝ ⎝ ⎠
⎛ ⎛ 2h f ⎞⎞
S f = Di ⋅ ⎜⎜ π − arccos⎜⎜ − 1⎟⎟ ⎟⎟
⎝ ⎝ Di ⎠⎠
2
⎛ 2h f ⎞
S g = Di ⋅ 1 − ⎜⎜ − 1⎟⎟
⎝ Di ⎠
S i = π ⋅ Di − S f
U f = Ut −
(U t − U l )H lls
H ltb
U m − U f ⋅ H ltb
Ug =
1 − H ltb
4 Af
Dhf =
Sf
4 Af
Dhg =
S g + Si
215
Apéndices
⎛ 0,0002 ⋅ ε 10 6 ⎞⎟
f wf = 0,001375 ⋅ ⎜1 + 3 +
⎜ Dhf Re l ⎟⎠
⎝
⎛ 0,0002 ⋅ ε 10 6 ⎞⎟
f wg ⎜
= 0,001375 ⋅ 1 + 3 +
⎜ Dhg Re g ⎟⎠
⎝
• En la interfase gas-líquido:
⎛ 300 ⋅ h f ⎞
f i = 0,005 ⋅ ⎜⎜1 − ⎟⎟
⎝ Di ⎠
ρl U f U f
τ f = f wf ⋅
2
ρg U g U g
τ g = f wg ⋅
2
• En la interfase gas-líquido:
ρ g U g − U f (U g − U f )
τ i = fi ⋅
2
216
Apéndices
• En el tapón de líquido:
ρ sU m2
τs = f ⋅
2
ρ f = ρ g ⋅ (1 − H ltb ) + ρ l ⋅ H ltb
Ls ⋅ (U l ⋅ H lls − U f ⋅ H ltb )
Lu =
U sl − U f ⋅ H ltb
L f = Lu − Ls
217
Apéndices
Flujo Burbuja:
3ρ l f mU m2
A=
8(ρ l − ρ g )g cos β
2 f mU m3
ε=
Di
Flujo Intermitente:
Número de Bond:
Bo =
(ρ l − ρ g )gDi2
σ
218
Apéndices
gDi (ρ l − ρ g )
U dh∞ = 0,54 − 1,76 Bo −0,56 ⋅
ρl
gDi (ρ l − ρ g )
(
U dv∞ = 0,345 ⋅ 1 − e 0,337−0,1Bo ⋅ ) ρl
1
H lls = 1, 39
⎛ Um ⎞
1 + ⎜⎜ ⎟⎟
⎝ 8,66 m s ⎠
Velocidad de la burbuja:
⎛ g (ρ l − ρ g )σ
0 , 25
⎞
U b = 1,53 sin (β )⎜⎜ ⎟
⎟
⎝ ρl 2 ⎠
U gdb = 1,2U m + U b
219
Apéndices
Densidad de la mezcla:
ρ m = ρ l H l + ρ g (1 − H l )
Viscosidad de la mezcla:
μ m = μ l H l + μ g (1 − H l )
Flujo Anular:
1
τi = f i ρ cU i U i
2
1
τ wl = f wl ρ lU f
2
Densidad en el núcleo
FeU sl
λlc =
U sg + FeU sl
ρ c = ρ l λlc + ρ g (1 − λlc )
μ c = μ l λlc + μ g (1 − λlc )
220
Apéndices
0 , 074
Fe ⎛ μ l 2U sg 2 ρ g ⎞ ⎛ U sg ⎞
0, 2
= 0,735⎜ ⎟ ⎜⎜ ⎟⎟
1 − Fe ⎜ σ 2 ρl ⎟ ⎝ U sl ⎠
⎝ ⎠
U c = (U sg + FeU sl )
1
(1 − 2δ )2
1
U f = U sl (1 − Fe )
δ (1 − δ )
Ui = Uc − U f
221
Apéndices
Flujo Anular:
Velocidad Uo :
⎛ g (ρ l − ρ g )σ
0 , 25
⎞
U o = 1,4⎜⎜ ⎟
⎟ sin (β )
⎝ ρl 2 ⎠
Si: U o ≥ 0,001 ⇒ H g =
(
U o + U sl + U sg − (U o + U sl + U sg ) − 4U oU sg
2
)0,5
2U o
U sg
Para otros valores: H g =
Um
Condición de estabilidad:
−0 , 2
U sl
2−0, 2 ((H f −Hf) + (1 − 2H ) )
2
f
2
(H − H )
f f
2 3
U sc = U sg + Fe ⋅ U sl
222
Apéndices
Fe = 1 − e − (0,125( f −1,5 ))
Donde f es igual a:
U sg μ g ρg
f = 10 4
σ ρl
I = I h ⋅ cos 2 β + I v ⋅ sin 2 β
Donde:
El parámetro de corrección interfacial horizontal se calcula:
I h = 1 + 850 f A
( ) ( )
0, 4
⎡ 0,707 Re sl 2,5
+ 0,0379 Re sl
2,5
⎤
⎢ ⎥⎦ ⎛ Ul ⎞ ρl
fA = ⎣ ⎜ ⎟
Re 0,9 ⎜U ⎟ ρ
sg ⎝ g ⎠ g
δ
I v = 1 + 300
Di
Flujo Burbuja:
gσ (ρ l − ρ g )
U o∞ = 1,53 ⋅ 4
ρ l2
223