Tesis Sobre Flujo Slug
Tesis Sobre Flujo Slug
Tesis Sobre Flujo Slug
Trabajo de ascenso presentado ante la Universidad Simn Bolvar como requisito para
ascender a la categora de Profesor Asociado
RESUMEN i
CONTENIDO ii
CAPITULO 1: INTRODUCCIN 1-1
CAPITULO 2: DISEO EXPERIMENTAL PARA LA MEDICIN DE LA 2-1
TENSIN SUPERFICIAL E INTERFACIAL A ALTAS PRESIONES
INTRODUCCIN 2-1
MARCO TEORICO 2-3
Resultados experimentales de la tensin interfacial de 2-3
estudios previos
Esquemas experimentales propuestos en la literatura 2-6
DISEO EXPERIMENTAL PROPUESTO 2-10
RESULTADOS Y DISCUSIONES 2-15
CONCLUSIONES 2-28
RECOMENDACIONES 2-29
REFERENCIAS BIBLIOGRAFICAS 2-30
CAPITULO 3: DETERMINACIN DEL COEFICIENTE DE DIFUSION 3-1
POR MEDIO DEL MTODO DE LA GOTA COLGANTE
INTRODUCCIN 3-1
MARCO TEORICO 3-3
Modelos de clculo de la difusividad 3-3
Modelo de transferencia de masa Yang y Gu (2005) 3-3
Volumen de la gota colgante dinmica 3-6
Funcin objetivo y de minimizacin 3-8
Mtodo numrico 3-12
Resultados experimentales de la difusividad 3-14
Modelo propuesto para el clculo de la difusividad 3-14
Modelo de difusin radial 3-15
Modelo de difusin radial y angular 3-22
Ecuaciones empricas 3-27
iii
INTRODUCCIN
comunes la inundacin por gases miscibles (por ejemplo con CO2) y la recuperacin trmica
de crudo pesado por inyeccin de gas-soluto.
El estudio de flujos multifsicos es muy importante a nivel industrial, especialmente en
la industria petrolera donde toma gran importancia, pues cuando se extrae el petrleo del
yacimiento junto a l viene asociada una mezcla de gas natural, agua y arena. Esta mezcla es
llevada a alguna planta donde es separada y finalmente refinada. En estos tramos, se presenta
flujo multifsico de slidos, lquido y gases. Por ejemplo en la industria petrolera, la
prediccin de un patrn de flujo puede ser muy importante, ya que dependiendo del patrn
agua-crudo que se presente, un oleoducto puede presentar mayor corrosin causada por el
agua, o puede haber cadas de presin elevadas, lo que aumenta los costos de bombeo y
mantenimiento.
En este trabajo se realizaron modificaciones al montaje experimental de tensin
superficial e interfacial ya existente en el laboratorio. Se realizaron mediciones de tensin
superficial empleando el mtodo de la gota colgante o emergente para el sistema CO2-Agua a
presiones manomtricas entre 0 y 5 MPa y a temperaturas entre 20 y 50 C, para compararlo
con datos reportados en la literatura y verificar la reproducibilidad de las mediciones con el
equipo experimental. Tambin se presentan los antecedentes de la tcnica experimental para la
medicin y calculo del coeficiente de difusin mediante el mtodo de la gota colgante. El
trabajo experimental se realiz a presin y temperatura constante, se proponen modelos para
el clculo de la difusividad y se presentan los resultados experimentales obtenidos para el
clculo del coeficiente de difusin. Los resultados obtenidos son posteriormente comparados
con los trabajos reportados en la literatura y a travs de correlaciones empricas desarrolladas
por otros autores para el clculo de la difusividad. Con respecto al estudio de flujo multifsico
para el sistema gas-lquido se emple la tcnica de visualizacin de imgenes para evaluar el
comportamiento de los diferentes parmetros del rgimen Slug (espesor de la pelcula a ambos
lados de la burbuja, la frecuencia de aparicin de las Burbujas de Taylor y la longitud y
velocidad de ascenso de las Burbujas de Taylor y del Slug de lquido), con el fin de determinar
la influencia de los caudales de aire y agua, y del ngulo de inclinacin sobre el
comportamiento de estos parmetros. Para el caso de flujo bifsico lquido-lquido se estudio y
midi la retencin de fases o holdup en tuberas a 90 y 45 de inclinacin, as como tambin
1-3
inclinadas, para cinco caudales de flujo de una mezcla aire-agua que ascienden por una tubera
a diferentes ngulos de inclinacin (entre 0-90 medidos desde la horizontal), se midi el
espesor de la pelcula a ambos lados de la burbuja, la frecuencia de aparicin de las Burbujas
de Taylor y la longitud y velocidad de ascenso de las Burbujas de Taylor y del Slug de lquido;
para as determinar la influencia de los caudales de aire y agua y del ngulo de inclinacin
sobre el comportamiento de estos parmetros. En el ultimo captulo del libro el cual
corresponde al Captulo 5 se muestran los resultados obtenidos del estudio en sistemas
multifsicos lquido- lquido, y en el cual se reporta la retencin de fases o holdup en flujo
bifsico lquido-lquido en tuberas a 90 y 45 de inclinacin, as como tambin se reportan
los diferentes patrones de flujo caractersticos los cuales fueron registrados mediante tcnicas
de visualizacin. Al final de los captulos se presentan los Apndices, los cuales contienen los
proceedings arbitrados. El Apndice A de este captulo contiene el trabajo arbitrado titulado:
Determination of the interfacial dynamic tension of the systems CO2-Decane, CO2-Dodecane
and H2O-Dodecane to high pressures with T = 20, 30, 40 and 50 C. El cual fue presentado en
el 16th Symposium on Thermophysical Properties en Boulder, Colorado (USA), en Julio 2006
y Determination of superficial dynamic tension of CO2+water from 293 K to 323 K with
pressures of 0-2 MPa, presentado en el VII Iberoamerican Conference on Phase Equilibria
and Fluid Properties for Process Desing. Morelia, Mxico, Octubre 2006. En el Apndice B de
este capitulo se anexa el trabajo arbitrado titulado: Modelo simplificado para la estimacin
del coeficiente de difusin mediante la tcnica de la gota colgante, y el cual fue presentado en
el IX Congreso Internacional de Mtodos Numricos y Ciencias Aplicadas. Margarita,
Venezuela, Febrero 2008. En el Apndice C de este capitulo se anexa el trabajo arbitrado
titulado: Estudio experimental de parmetros caractersticos del flujo slug para sistemas
lquido-gas en tuberas inclinadas, presentado en el XXII Congreso Interamericano de
Ingeniera Qumica. Buenos Aires, Argentina, Octubre 2006. En el Apndice D de este
capitulo se anexa el trabajo arbitrado titulado: Estudio experimental sobre patrones de flujo
lquido-lquido en tuberas verticales, presentado en el 8vo Congreso Iberoamericano de
Ingeniera Mecnica. Cusco, Per, Octubre 2007. En el Apndice E se presentan las memorias
arbitradas y no arbitradas de Congresos que se participaron durante el periodo de ascenso, y
los cuales estn relacionados con el rea de fenmenos interfaciales y las tcnicas de
visualizacin para el estudio interfaces y de sistemas multifsicos.
CAPITULO 2
INTRODUCCIN
inyectado en el crudo bajo las condiciones del yacimiento (Sun y Chen, 2004, 2005; Boustani
y Maini, 2001). Es por esto que se considera de gran importancia el estudio de la variacin de
la tensin interfacial con un cambio de presin y temperatura en sistemas lquido-gas y
lquido-lquido.
En esta seccin se har una revisin de los diferentes montajes experimentales
propuestos en la literatura para sistemas a altas presiones empleando la tcnica de la gota
colgante o emergente. Se mostrara de igual forma en este capitulo las modificaciones
realizadas al equipo de medicin de tensin interfacial existente en el laboratorio de
Visualizacin, ya que con el equipo solo se poda trabajar a presin atmosfrica y a
temperaturas no mayores a 40 C. Esto es debido a que a temperaturas por encima de est
ocurra empaamiento de los visores de la celda experimental, lo que originaba la distorsin
del contorno de la gota, obtenindose una imagen difusa, la cual era imposible emplear para el
calculo de la tensin por medio de la medicin de los parmetros geomtricos de la gota. El
objetivo en este capitulo fue realizar la modificaciones del montaje experimental de tensin
superficial e interfacial ya existente en el laboratorio. Se realizaron mediciones de tensin
superficial para el sistema CO2-Agua a presiones manomtrica entre 0 y 5 MPa y a
temperaturas entre 20 y 50 C, para compararlo con datos reportados en la literatura y verificar
la reproducibilidad de las mediciones con el equipo experimental. En el Apndice A se
encuentran los trabajos arbitrados como memorias de congresos, relacionados con el tema
desarrollado en este captulo los cuales se titulan:
Determination of superficial dynamic tension of CO2 + water from 293 K to 313 K
with pressures of 3, 4 and 5 MPa. El cual fue presentado en el 16th Symposium on
Thermophysical Properties en Boulder, Colorado (USA), en Julio 2006. (Melndez J.,
Zeppieri S. y Lpez A., 2006).
Determination of superficial dynamic tension of CO2+water from 293 K to 323 K
with pressures of 0-2 MPa. VII Iberoamerican Conference on Phase Equilibria and
Fluid Properties for Process Desing. Morelia, Mexico, Octubre 2006. (Melndez J.,
Zeppieri S. y Lpez A., 2006).
2-3
MARCO TEORICO
Tewes y Boury (2004) presentaron como parte de sus resultados el cambio de la tensin
interfacial para el sistema agua-dixido de carbono en un intervalo de 60.000 s, bajo las
condiciones de presin de 30 a 80 bar y de temperaturas de 20, 30 y 40 C tambin
presentaron la variacin de dicha tensin interfacial tanto inicial como de equilibrio con la
presin del sistema. Tewes y Boury emplearon la tcnica de la gota colgante con temperatura
y volumen constantes tomando cambios de presin. Los resultados obtenidos se presentan en
las Figs. 2.5 a 2.7
Figura 2.5 Tensin interfacial vs. tiempo del sistema agua-CO2 para
la temperatura de 40 C (Tewes y Boury, 2004).
2-6
MPa, cambiando las composiciones de cada gas. Dichos autores desarrollaron el experimento
manteniendo la presin y la temperatura constantes y modificando el volumen de la gota
colgante durante un intervalo de tiempo para cada mezcla de gases. El esquema experimental
desarrollado por dichos autores se muestra en la Fig. 2.8:
Figura 2.8 Diagrama esquemtico de la unidad experimental: (1) Celda de alta presin, (2)
Bao trmico, (3) Cilindros de muestra, (4) Bomba JEFRI 100-1-1 HB, (5) Bomba JEFRI 10-
1-12-NA, (6) Lente
Yan et al. (2001) realizaron estudios de la tensin interfacial para los sistemas (metano
+ nitrgeno) + agua y (dixido de carbono + nitrgeno) + agua aplicando la teora del
gradiente lineal (LGT) (Jha, 1986; Zuo Stenby, 1996). Para ello modificaron el esquema
experimental propuesto por Ren et al. 2000 para trabajar con rangos de presin de 1 a 30 MPa
y temperatura de 298-373 K. El esquema experimental utilizado se muestra en la Fig. 2.9.
2-8
Figura 2.9 Diagrama esquemtico de la unidad experimental: (1) Celda de alta presin,
(2) Bao trmico, (3) Cilindros muestra, (4) Bomba JEFRI 100-1-1 HB, (5) Tanque de Gas,
(6) Bomba JEFRI 10-1-12-NA
Hebach et al. (2002) realizaron sus estudios tensin empleando el sistema CO2-H2O,
trabajaron en un rango de presin 0,1 a 20 MPa y temperatura de 278 a 335 K, con el diseo
experimental propuestos por estos autores logrando disminuir el error experimental al 2%. El
diseo experimental empleado se muestra en la Fig.. 2.10
Figura 2.10 Celda de alta presin (A), Ventanas de inspeccin (B), Lmpara (C), Cmara
CCD (D), Bomba de engranajes (E), Filtro (F), Separador de Fase (G), Termostato (J), Bomba
de pistn (K), Bomba de jeringa (L), Compresor (M) (Hebach et al., 2002).
2-9
Tewes y Boury (2004) realizaron estudios sobre las propiedades termodinmicas del
sistema H2O-CO2 a altas presiones en un rango de temperatura de 20-40 C entre dichas
propiedades termodinmicas incluyen la tensin interfacial. Tewes y Boury midieron la
tensin interfacial registrando un cambio de presin con el volumen de gota y temperatura
constante, en un intervalo de tiempo de 60.000 s, suficiente para que el sistema alcanzara su
estado estacionario. Los materiales que emplearon fueron: un cilindro de dixido de carbono
de una pureza de 99 % Vol. y agua tri-destilada; como equipos utilizaron una celda de alta
presin para formar la gota de agua, una lmpara, una cmara CCD, una computadora y una
jeringa con un capilar de 1mm de dimetro interno. El esquema experimental empleado se
muestra en la Fig. 2.11.
Figura 2.11 Montaje experimental de la gota colgante con temperatura y volumen constante
con la presin variable (Tewes y Boury, 2004).
2-10
Figura 2.12 Diagrama esquemtico de la unidad experimental: (1) Celda de alta presin, (2-a,
2-b) Baos trmicos, (3) Cilindros muestra, (3-1, 3-2 y 3-3) Cuerpo, pistn y tapa del cilindro
muestra, (4) Cilindro de dixido de carbono, (5) Cilindro de aire comprimido, (6) Lente
microscpico, (7) Video cmara, (8) Mesa antivibracin, (9) Computadora; (10) Lmpara
de fibras ptica, (11) Bomba de vaco, (P) Manmetro digital, (T) Termmetro digital,
(V-M) Vlvula micro-mtrica, (V-2, V-3, V-4, V-6) Vlvulas de aguja,
(V-1, V-5) Reguladores de presin.
2-11
El equipo experimental consta de una celda de acero inoxidable, diseada para soportar
presiones de hasta 35 MPa, en la parte frontal y posterior de la celda se encuentra un orificio
de cada lado sellado con cristales de cuarzo de 6 mm de espesor, dichos orificios permiten
visualizar la gota por un lado y por el otra la iluminacin de la gota, esto ultimo se realiza con
el fin de lograr una distribucin uniforme de la imagen de la gota, Fig. 2.13. El diseo de dicha
celda fue elaborado en la Universidad de Tulsa (Lopz de Ramos, 1993) para realizar
mediciones de la tensin superficial del sistema aire-agua, dicha celda tambin fue usada para
medir tensiones interfaciales del sistema alcanos-agua en la Universidad Simn Bolvar,
Lpez de Ramos et al. (1993, 1995) y Zeppieri et al., (2001, 2002). El cuerpo principal de la
celda posee cuatro canales que la atraviesan hasta la cmara principal, los canales presentan
conexiones roscadas que permiten la colocacin del capilar en el tope, el drenaje del lquido
por el fondo de la celda y la circulacin de agua transversalmente para mantener la
temperatura del sistema. La celda est colocada sobre una mesa antivibraciones (TMC 78A-
455-12) para evitar perturbaciones al formar la gota.
Figura 2.14 Cilindro de alta presin que contienen los fluidos del sistema en estudio.
(a) (b)
Figura 2.18 Vlvula de aguja (a). Vlvula micro-mtrica de aguja (b). Vlvula regulador de
presin.
El sistema de visualizacin, Fig. 2.19, consta de una cmara de video (SONY CCD
XC-ST70), una lente microscopica de 60mm (Kikkon) con dos anillos expansores (Nikkon
PK-13), una lmpara de fibra ptica (Dolan-Jenner MI-150), una tarjeta capturadota de video
(Pinnacle System DC-1000) y una computadora para la visualizacin y captura de la imagen
de la gota en el tiempo, as como para la adquisicin de los datos de presin y temperatura del
sistema.
RESULTADOS Y DISCUSIONES
observar que son muy prximos, presentando una diferencia de 0,67 mN/m. Para la condicin
de P = 3 MPa con T = 50 C se realizaron las mediciones 6 veces con el fin de verificar la
reproducibilidad de los experimentos obtenindose una desviacin estndar de 0,25 mN/m. En
la Fig. 2.21 se muestran los resultados obtenidos de tensin en tiempo para diferentes
presiones y temperaturas, se observar que el comportamiento de esta isoterma vari
asintticamente hasta alcanzar un valor de equilibrio.
Figura 2.21 Calibracin del equipo y mtodo experimental de la tensin interfacial para el
sistema CO2-agua. P =1 MPa y T = 25,2 C, P =2 MPa y T = 25,2 C,
P =3 MPa y T = 50 C
Figura 2.22 Variacin de la tensin interfacial del sistema CO2-agua a 20 C con respecto al
tiempo para algunas presiones del CO2. (+) 0 MPa, ( ) 1 MPa,
() 2 MPa, () 3 MPa, () 4 MPa, () 5 MPa.
Figura 2.23 Variacin de la tensin interfacial del sistema CO2-agua a 30 C con respecto al
tiempo para algunas presiones del CO2. (+) 0 MPa, ( ) 1 MPa,
() 2 MPa, () 3 MPa, () 4 MPa, () 5 MPa.
2-19
Figura 2.24 Variacin de la tensin interfacial del sistema CO2-agua a 40 C con respecto al
tiempo para algunas presiones del CO2. (+) 0 MPa, ( ) 1 MPa,
() 2 MPa, () 3 MPa, () 4 MPa, () 5 MPa.
Figura 2.25 Variacin de la tensin interfacial del sistema CO2-Agua a 50 C con respecto al
tiempo para algunas presiones del CO2. (+) 1 MPa, ( ) 2 MPa,
() 3 MPa, () 4 MPa, () 5 MPa.
Se puede observar que las isotermas para cada presin se aproximan ms al aumentar
la temperatura del sistema; esto es causado por la disminucin de la tensin interfacial al
2-20
mezclarse por difusin el CO2 con las molculas de agua en la gota, y por el equilibrio msico
que produce el efecto de evaporacin al aumentar la temperatura, ya que las molculas de agua
en la interfase se transfieren con mayor facilidad a la fase gaseosa reduciendo el tamao de la
gota. Esto se puede observar con claridad para las isotermas de presiones de 0 y 1 MPa, donde
el cambio de dicha tensin con el tiempo es mayor que para las isotermas de altas presiones.
Analizando los cambios de la tensin interfacial mostrados en las Figs. 2.26 a la 2.31
que corresponden a las condiciones de temperatura variable y presin constante, se observa un
comportamiento similar para las condiciones de temperatura constante y presin variable
reportados en la literatura, tanto en la tensin interfacial inicial como la de equilibrio, en
donde la tensin alcanza un valor asinttico en tiempos superiores a los 30.000 s.
Figura 2.26 Variacin de la tensin interfacial del sistema CO2-agua a 0 MPa con respecto al
tiempo para algunas temperaturas. () 20 C, () 30 C, () 40 C, ( ) 50 C.
2-21
Figura 2.27 Variacin de la tensin interfacial del sistema CO2-agua a 1 MPa con respecto al
tiempo para algunas temperaturas. () 20 C, () 30 C, () 40 C, ( ) 50 C.
Figura 2.28 Variacin de la tensin interfacial del sistema CO2-agua a 2 MPa con respecto al
tiempo para algunas temperaturas. () 20 C, () 30 C, () 40 C, ( ) 50 C.
2-22
Figura 2.29 Variacin de la tensin interfacial del sistema CO2-agua a 3 MPa con respecto al
tiempo para algunas temperaturas. () 20 C, () 30 C, () 40 C, ( ) 50 C.
Figura 2.30 Variacin de la tensin interfacial del sistema CO2-agua a 4 MPa con respecto al
tiempo para algunas temperaturas. () 20 C, () 30 C, () 40 C, ( ) 50 C.
2-23
Figura 2.31 Variacin de la tensin interfacial del sistema CO2-agua a 5 MPa con respecto al
tiempo para algunas temperaturas. () 20 C, () 30 C, () 40 C, ( ) 50 C.
En las Figs. 2.26 a 2.31 se muestran la variacin de eq durante los 30.000 s para cada
una de las temperaturas en estudio, en estas figuras se puede observar la dispersin de algunos
puntos en las isotermas, que es causada por el constante movimiento de la interfase, producto
de la transferencia del CO2 hacia el bulto de la gota (Fig. 2.20), pero a pesar de dicho
movimiento, las isotermas describen un claro comportamiento al disminuir con respecto al
tiempo a medida que el sistema alcanza el equilibrio.
En las Fig. 2.32 y 2.33 se reporta la tensin interfacial de equilibrio (la cual
corresponde al valor constante de tensin interfacial alcanzado en un determinado tiempo) y la
tensin inicial (la cual es el valor de tensin reportada al inicio de cada experimento) en
funcin de la presin, para diferentes temperaturas. Se puede observar que a medida que
aumenta la presin la tensin interfacial del sistema disminuye, lo mismo ocurre para la
temperatura, a medida que esta aumenta la tensin disminuye; sin embargo podemos observar
que el efecto de la temperatura en la tensin se invierte en un punto de presin en el cual a
partir de este, a medida que aumenta la temperatura del sistema la tensin interfacial tambin
aumenta. En el punto de intercepcin los efectos causados por la temperatura se igualan a los
efectos producidos por la presin del CO2 sobre las molculas de la interfase, originndose
2-24
Figura 2.32 Variacin de la tensin interfacial en equilibrio del sistema CO2-agua con
respecto a la presin del CO2. ( ) 20 C, ( ) 30 C, ( ) 40 C, ( ) 50 C.
Figura 2.33 Variacin de la tensin interfacial inicial del sistema CO2-agua con respecto a la
presin del CO2. ( ) 20 C, ( ) 30 C, ( ) 40 C, ( ) 50 C.
2-25
Comparando los resultados obtenidos en las Figs. 2.32 y 2.33 con los publicados por
Tewes y Boury (2004) (Figs 2.6 y 2.7) para las temperaturas de 20, 30 y 40 C, se observa que
los valores de la tensin interfacial en la presin de intercepcin de las isotermas son los
mismos, a diferencia del comportamiento de la isoterma a 40 C, que a bajas presiones reporta
un descenso brusco causado por el efecto de evaporacin; tambin se puede observar dicho
comportamiento para la isoterma a 50 C. Es importante acotar que a pesar de que Tewes y
Boury (2004) trabajaron a temperatura y volumen de la gota constante con presin del sistema
variable, en este trabajo se pudo reproducir la tendencia y comportamiento observado por
estos autores, trabajando a presin y temperatura del sistema constante con el volumen de la
gota variable, en un intervalo de tiempo menor. Las desviaciones estndar de los resultados de
este trabajo y los publicados por Tewes y Boury (2004) se presentan en la tabla 2.1.
3,0 MPa
20 4,87 % 4,61 %
30 7,59 % 2,60 %
40 8,42 % 3,15 %
4,0 MPa
20 0,88 % 1,33 %
30 1,14 % 3,82 %
40 0,71 % 3,52 %
5,0 Mpa
20 0,15 % 0,31 %
30 2,28 % 0,31 %
40 0,92 % 0,89 %
Haciendo la comparacin con los datos publicados por Chun y Wilkinson (1995),
quienes emplearon la tcnica de capilaridad para la medicin de la tensin interfacial a
diferentes presiones y temperaturas tal como se muestra en la Fig. 2.34, se observa que las
isotermas a 30 C en equilibrio presentan una gran diferencia, pero describe un
comportamiento similar. Esto se debe a que las tensiones interfaciales reportadas por dichos
autores no son para un estado de equilibrio con un intervalo de tiempo tan amplio como se
desarroll en este trabajo, sino para valores iniciales, tal como lo muestra la isoterma a 30 C
inicial; tambin se observa el punto de intercepcin de las isotermas descrito anteriormente en
la Fig. 2.34, pero a una presin ligeramente mayor, posiblemente por causa de los errores que
contiene la tcnica o por la poca exactitud de los equipos usados para la poca.
Figura 2.34 Comparacin de la variacin de la tensin interfacial de este trabajo con los datos
publicada por B. Chun y G. Wilkinson para el sistema CO2-gua con respecto a la presin del
CO2. ( ) Este trabajo a 30 C en Equilibrio, ( ) Este trabajo a 30 C Inicial, ( ) Chun
a 25 C, ( ) Chun a 30 C, ( ) Chun a 38 C, ( ) Chun a 53 C.
Al comparar estos resultados con los obtenidos por Da Rocha et al., (1998) y Hebach
et al. (2002) , Figs. 2.2 y 2.4, se puede observar que el comportamiento de las isotermas en
2-27
dichas figuras es semejante al obtenido en este trabajo de tesis reportado en la Fig. 2.32. En los
resultados de Hebach et al. (2002) se distingue claramente que la tensin interfacial inicial de
la isoterma de 30 C (303 K) es mayor que para la de 45 C (318 K), luego se describe un
solapamiento entre las mismas hasta llegar a una presin de 2,5 MPa, donde se interceptan y
cambian de valores, siendo la de mayor tensin la de 45 C y la de menor para la de 30 C a la
presin de 5 MPa.
2-28
CONCLUSIONES
La gota colgante presenta una interfase mvil que vara para todas las presiones y
temperaturas durante el tiempo requerido para alcanzar las condiciones de equilibrio.
RECOMENDACIONES
Colocar mantas de calentamiento en las lneas de distribucin del CO2 con el fin de
garantizar que el gas se encuentre a la misma temperatura del sistema.
2-30
REFERENCIAS BIBLIOGRAFICAS
Adamson A. W.; Gast A. P., Physical Chemistry of Surfaces, 6th Ed., John Wiley & Sons
Inc., New York 1967.
Boustani A.; Maini B., The Role of Diffusion and Convective Dispersion in Vapour
Extraction Process, J. Can. Pet. Technol., 40, 68-77, 2001.
Chun B.; Wilkinson G., Interfacial Tension in High-pressure Carbon Dioxide Mixtures,
Ind. Eng. Chem. Res., 34, 4371-4377, 1995.
Chun B.; Wilkinson G., Interfacial Tension in High-Pressure Carbon Dioxide Mixtures, Ind.
Eng. Chem. Res., 34, 4371-4377, 1995.
Cussler E.L., Diffusion, Mass Transfer in Fluid System, 1a Ed., The Press Syndicate of the
University of Cambridge, New York 1984.
Da Rocha S.; Harrison K. and Johnston K., Effect of Surfactants on the Interfacial Tension
and Emulsion Formation between Water and Carbon Dioxide, Langmuir, 15, 419-428, 1998.
Hebach A.; Oberhof A.; Dahmen N.; Kgel A.; Ederer H. and Dinjus E., Interfacial Tension
at Elevated Pressures Measurements and Correlations in the Water + Carbon Dioxide
System, J. Chem. Eng. Data, 47, 1540-1546, 2002.
Jennings H. y Harley Y., The effect of temperature and pressure on the interfacial tension of
benzene-water and normal decane-water, J. of Colloid and Interface Science, 24, 323-329,
1967.
Jennings J. W.; Pallas N. R., An Efficient Method for the Determination of Interfacial
Tensions from Drop Profiles, Langmuir, 4, 4, 960-967, 1988.
Jha K. N., A Laboratory Study of Heavy Oil Recovery with Carbon Dioxide, J. Can. Pet.
Technol., 25, 54-63, 1986.
2-31
Lpez de Ramos A., Capillary Enhanced Diffusion of CO2 in Porous media, PhD.
Dissertation, The university of Tulsa,1995.
Lpez de Ramos A.; Redner R. A. y Cerro R. L., Surface Tension from Pendant Drop
Curvature, Langmuir, 9, 3691-3694, 1993.
Moritis G., CO2 and HC Injection Lead EOR Production Increase, Oil & Gas, Journal, Apr
23, 49-82, 1990.
Orr F. M. Jr.; Taber J. J., Use of Carbon Dioxide in Enhanced Oil Recovery, Science, 224,
563-569, 1984.
Ren Q.; Chen G.; Yan W. and Guo T., Interfacial Tension of (CO2 + CH4) + Water from 298
K to 373 K and Pressures up to 30 MPa, J. Chem. Eng. Data, 45, 610-612, 2000.
Smith J.; Van Ness H., Termodinmica en Ingeniera Qumica, sexta ed., Mc Graw Hill,
Mxico 2003.
Sun C. Y.; G. J. Chen, Measurement of Interfacial Tension for the CO2 Injected Crude Oil +
Reservoir Water System, J. Chem. Eng. Data, 50, 936-938, 2005.
Sun C.; Chen G. and Yang L., Interfacial Tension of Methane + Water with Surfactant near
the Hydrate Formation Conditions, J. Chem. Eng. Data, 49, 1023-1025, 2004.
Sun C.; Chen G., Measurement of Interfacial Tension for the CO2 Injected Crude Oil+
Reservoir Water System, J. Chem. Eng. Data, 50, 936-938, 2005.
2-32
Tewes F.; Boury F., Thermodynamic and Dynamic Interfacial Properties of Binary Carbon
Dioxide-Water Systems, J. Phys. Chem. B, 108, 2405-2412, 2004.
Yan W.; Zhao G.; Chen G. and Guo T., Interfacial Tension of (Methane + Nitrogen) + Water
and (Carbon Dioxide + Nitrogen) + Water Systems, J. Chem. Eng. Data, 46, 1544-1548,
2001.
Yang C.; Gu Y., A New Method for Measuring Solvent Diffusivity in Heavy Oil by Dynamic
Pendant Drop Shape Analysis (DPDSA), Society of petroleum Engineers (SPE) Annual
Technical Conference and Exhibition, Denver, Co, October 5-8,; Papers SPE 84202, 2003.
Yang C.; Gu Y., Modeling of the Adsorption Kinetics of Surfactants at the Liquid-Fluid
Interface of a Pendant Drop, Langmuir, 20, 2503-2511, 2004.
Yang C.; Gu Y., New Experimental Method for Measuring Gas Diffusivity in Heavy Oil by
the Dynamic Pendant Drop Volume Analysis (DPDVA), Ind. Eng. Chem. Res., 44, 4474
4483, 2005.
Yang D.; Tontiwachwuthikul P. and Gu Y., Interfacial Tensions of the Crude Oil + Reservoir
Brine + CO2 Systems at Pressures up to 31 MPa and Temperatures of 27 C and 58 C, J. of
Chem. and Eng. Data, 50, 1242-1249, 2005.
Zuo X.; Stenby H., A Linear Gradient Theory Model for Calculating Interfacial Tensions of
Mixtures, Journal of Colloid and Interface Science, 182 (1), 126-132, 1996.
CAPITULO 3
INTRODUCCIN
MARCO TEORICO
Figura 3.1 (a) Diagrama esquemtico de la gota colgante simtrica rodeada por gas descrita
en coordenadas cilndricas (r,z). (b) Mallado triangular de la gota
colgante usado en el mtodo computacional (Yang y Gu, 2005).
Donde rn, hn y n son el radio, la altura y el espesor del capilar; es la fase que
constituye a la gota de petrleo pesado, junto con la fase que est dentro del capilar, n la capa
lmite formada entre el capilar y la fase de la gota, int entre dicha fase y el CO2.
Al no ser voltil la gota, se asegura que existe slo transferencia de masa desde el gas
hasta el bulto de dicha gota. Durante el proceso de transferencia de masa, la tensin interfacial
disminuye gradualmente y el volumen de la gota aumenta lentamente por el efecto de
hinchamiento. La ecuacin diferencial de transferencia de masa por difusin en coordenadas
cilndricas que describe la distribucin de la concentracin del CO2 dentro de la gota ()
puede ser expresada como:
c
r 2
c r c
= D
1
+ 2 , (r,z) , t > 0 (3.1)
t r r z
donde Csat es la concentracin molar de saturacin del gas en la gota de petrleo pesado a
presin y temperatura experimental. Csat est definida como el nmero de moles del gas
disuelto por unidad de volumen de petrleo pesado a presin y temperatura de saturacin.
El capilar es impermeable al CO2, por lo tanto, el plano cortante tambin puede ser
considerado como un lmite impermeable, as que la consideracin de la no penetracin puede
ser aplicada a la frontera n:
c c
D r + z = 0, (r,z) n, t >0 (3.4)
r z
donde r y z son los cosenos direccionales en las componentes r y z externas del vector
normal unitario al contorno, esta condicin de borde es de segundo tipo, la cual es conocida
como condicin de Neumann.
Para determinar la difusividad (D) con la ec. 3.1, la condicin inicial de borde ec. 3.2,
de equilibrio ec. 3.3 y la condicin de lmite de la no penetracin ec. 3.4, es necesario
adimensionar las variables como se muestra a continuacin:
3-6
c r z Dt
C= , R= , Z = , = 2 (3.5)
Csat rn rn r n
C 1 C 2C
= R + , (R,Z), >0 (3.6)
R R R 2
C ( R, Z , ) =0 = 0, (R,Z) (3.7)
R y Z son los cosenos direccionales en las componentes R y Z externas del vector normal
unitario al contorno.
La Ec. (3.6), junto con las Ecs. (3.7) y (3.8), se resuelven numricamente aplicando el
mtodo semidiscreto de elemento infinito de Galerkin. En este mtodo numrico la
discretizacin consiste en dos pasos, uno la discretizacin espacial y el otro la temporal. En la
discretizacin de las variables espaciales se usa el mtodo antes mencionado, donde la
ecuacin parcial diferencial es transformada a un sistema de ecuacin diferencial ordinaria
dependiente del tiempo (sEDO), mientras que en la discretizacion temporal el sistema
resultante de la sEDO es resuelto usando el mtodo (Burnett, 1987).
En la Fig. 3.1b se muestra el mallado triangular de la simulacin numrica, usando un
mallado fino cerca de la interfase CO2-petrleo pesado. Los valores usados en el clculo
numrico son los siguientes: radio del capilar rn= 0,79 mm, espesor de la pared n=0,59 mm,
altura del capilar hn=2,4 mm.
inicial de dicha gota a un tiempo t=0, gas representa el cambio del volumen por unidad de
k sw 1
gas = (3.12)
M gas C sat
donde C sat es el nmero de moles del CO2 disuelto por unidad de volumen de la gota a
presin y temperatura de saturacin. Sustituyendo las Ecs. 3.10 y 3.12 en la Ec. 3.11 se tiene:
3-8
Dt
Vc (t ) = V0 + (k sw 1)rn3 C R, Z , 3 RdRdZ (3.13)
( R , Z ) rn
t=0, Vm (t ) es medido con la primera imagen digital tomada, una vez formada completamente
la gota de petrleo, el cual indica que Vm (0) = V0 . Los moles de gas disuelto a cualquier
tiempo t despus de tiempo inicial, pueden ser determinados por la siguiente ecuacin:
D(t + t f ) Dt
n(t + t f ) n(t ) = C sat rn3 C R, Z , C
R , Z , RdRdZ (3.14)
( R , Z )
rn2 rn2
n(t f ) representa el nmero de moles del CO2 disuelto en la gota durante el tiempo t f . Al
emplear la Ec. 3.14 para determinar el volumen de la gota colgante desde un tiempo t f hasta
D (t + t f ) C R, Z , Dt RdRdZ
Vc (t ) = V0 + (k sw 1)rn3 C R, Z , (3.15)
( R , Z )
rn2 rn2
tiempo de 0 t t m :
2
1 tm Vm (t ) Vc (t )
E=
tm 0
Vm (t )
dt .100%
(3.16)
D (t + t f ) C R, Z , Dt RdRdZ
Q( D, t ) = rn3 C R, Z , (3.17)
rn2 rn2
( R , Z )
Una vez encontrado el mnimo de la funcin objetivo Emin se pueden determinar los
valores correspondientes del D y k sw . Para esto se deriva la funcin objetivo con respecto a
k sw :
3-10
1 tm [V (t ) V ] Q ( D, t )
2
tm Q ( D , t )
( k sw 1) 0
m 0
dt dt
E ( D, k sw ) tm
2 2
Vm (t ) Vm (t )
0
= (3.19)
k sw 2
1 tm Vm (t ) V0 (ksw 1)Q( D, t )
tm 0 Vm (t ) dt
Cuando la derivada parcial de la ec. 3.19 se iguala a cero, el factor de hinchamiento del
crudo pesado k sw queda determinado por:
tm [Vm (t ) V0 ] Q( D, t ) dt
0 Vm2 (t )
k sw = 1 + 2 (3.20)
tm Q ( D , t )
0 Vm2 (t ) dt
1. Se toma la secuencia de fotos de la gota colgante rodeada por CO2 para determinar Vm (t )
3. La Ec. (3.6) de transferencia de masa conjuntamente con las Ecs. (3.7) a la (3.9) se
resuelven aplicando elemento finito semidiscreto de Galerkin, para obtener el perfil de
concentraciones C ( R, Z , ) del CO2 dentro del dominio de la gota de crudo pesado.
4. El intervalo supuesto inicial de coeficiente de difusin D del CO2 en la gota de crudo
pesado se selecciona como [ D0 , D10 ].
ptimo de k sw,i , i = 0,1, 2,...,10 que se determina por la Ec. (3.20), y con ambos parmetros
se calcula E ( Di , k sw,i ) .
7. Para cada par de ( Di , k sw,i ), i = 0,1, 2,...,10 se calcula el volumen terico Vc (t ) de la gota
difusin.
9. Con el nuevo intervalo supuesto [Dn-1,Dn+1] se realiza el procedimiento hasta que la
diferencia entre el intervalo sea lo suficientemente pequeo para satisfacer la exactitud
especificada.
Se puede ver que para la determinacin de la difusividad (coeficiente de difusin) Dn y
k sw , son los valores globales ptimos de la difusividad del gas y el factor de hinchamiento, los
contrario se seguir el procedimiento desde el paso 5 con el nuevo intervalo [Dn-1,Dn+1] hasta
llegar al nivel de tolerancia buscada para la difusividad.
3-12
Mtodo numrico
Para resolver la Ec. (3.1) con la condicin inicial de la Ec. (3.2) y de frontera, de las
Ecs. (3.3) y (3.4) es conveniente expresar todas las ecuaciones de transferencia de masa en
ecuaciones adimensionales como lo describe la Ec. (3.5).
Las Ecs (3.6) (3.9) presentan un problema tpico de valor inicial de frontera el cual
describe una distribucin transiente en dos dimensiones de la concentracin del CO2 en la gota
colgante de crudo pesado. Debido a que el mtodo de elemento finito es especialmente
apropiado para trabajar con lmites curvos, Yang y Gu (2005) emplearon el mtodo
semidiscreto de elemento de Galerkin para resolver la ec. (3.1), con la condicin inicial de la
Ec. (3.2) y las de frontera con las Ecs. (3.3) y (3.4); obteniendo ms precisin en la
distribucin de la concentracin del CO2 cerca de la interfase usando un mallado fino en sta,
y un mallado ms amplio en el centro de la gota colgante de crudo pesado. La solucin de la
Ec. (3.1) (ecuacin de difusin) por elemento finito, se obtiene usando aproximacin
semidiscreta de Galerkin.
~ 3
C ( e ) ( R, Z , ) = C (j e ) ( ) ( e ) ( R, Z ) (3.21)
j =1
~
donde C ( e ) ( R, Z , ) es la aproximacin de la concentracin adimensional del gas en un
en los tres nodos del elemento triangular; y (j e ) ( R, Z ), j = 1, 2 y 3 son las llamadas funciones
de forma. Yang y Gu (2005) aplicaron una interpolacin lineal para determinar la solucin
aproximada en cada elemento, y la correspondiente funcin que puede ser encontrada en Wait
y Mitchell (1985). Las formulaciones elementales estn deducidas para la adaptacin del
mtodo residual de Galerkin (The Galerkin weighted residual method) (Burnett, 1987). La
ecuacin residual de Galerkin correspondiente a la ecuacin de difusin est dada por:
~
~ (e) ~
C (e)
1 C C (e)
2
RSD ( e ) ( R, Z , ) = R , ( R, Z ) , > 0 (3.22)
R R R Z 2
3-13
dC ( e ) ( )
[ H ]( e ) {(e) (e)
} {
+ [ K ] C ( ) = F ( )
(e)
} (3.23)
d
Donde:
i( e ) j( e ) i( e ) j( e ) i( e ) (j e )
K (e)
= + dRdZ , i, j = 1, 2 y3 (3.25)
( e ) R
ij
R Z Z R R
~ (e) ~
C C (e)
(e) (e)
Fi ( e ) = i nR +
(e) (e)
i n Z dS , i = 1, 2 y3 (3.26)
R Z
(e)
Donde n Z(e ) y n R(e ) son los cosenos direccionales, R y Z del vector exterior unitario
normal para el elemento de borde; y S = s / rn es la coordenada adimensional a lo largo de
dicho elemento en la direccin en sentido contrario a las aguja del reloj. Cuando la ecuacin
de los elementos es ensamblada en un elemento triangular, el sistema de ecuaciones resultante
puede ser escrito de la siguiente forma matricial:
dC ( )
[ H ] + [ K ]{C ( )} = {F ( )} (3.27)
d
1
[H ] + [K ]{C}n = (1 ){F }n1 + {F }n + 1 [ H ] (1 )[ K ]{C}n1 (3.28)
n n
Figura 3.2 Gota esfrica con transferencia de masa por difusin radial.
c m c 1 c 1 c
+ U r + Um + Um
t r r r sen( )
(3.29)
1 c 1 c 1 2c
= DAB 2 r 2 + 2 sen( ) + 2 2
+ Rv
r r r r sen( ) r sen2
( )
c 1 c 1 c 1 2c
= DAB 2 r 2 + 2 sen( ) + (3.30)
t r r r r sen( ) r 2 sen2 ( ) 2
La esfera se puede obtener al rotar una circunferencia en el eje Z (por ser una figura
simtrica) como se puede observar en la Fig. 3.6. Si se desarrollan todas las ecuaciones para
esta circunferencia definida con las coordenadas r, y la giramos en el eje Z, para este caso,
los ngulos y en la Ec (3.30) de transferencia de masa son iguales; pero adems si
consideramos que slo hay cambios de la concentracin del CO2 en sentido radial y no angular
(como lo muestra la Fig. 3.5) la Ec. (3.30) queda de la siguiente manera:
c 1 c
= DAB 2 r 2 (3.31)
t r r r
3-17
c r t
C= R= =
c sat rn r
n
2 (3.32)
D
C 1 2 C
= R
R2 R R
(3.33)
C 1 2 C
= R
R2 R R
(3.34)
C (R, ) = 0 = 0 (3.35)
R
U= (3.37)
S ( )
rad
Donde S () es para cada tiempo y R es el radio adimensional definido en la Ec.( 3.32).
rn
Esto hace que para todos los tiempos la nueva variable U abarque un intervalo de 0
hasta 1. Las ecuaciones finales quedan de la siguiente forma:
C 1 2 C
= U
U2 U U
(3.38)
C (U , ) =0 = 0 (3.39)
C (U , ) U =1 = 1 (3.40)
Esta ecuacin diferencial con sus condiciones (inicial y de borde) se resuelve por el
mtodo de Crank-Nicholson de diferencias finitas, de esta forma se podr calcular el perfil de
concentraciones de CO2 en la gota colgante de agua. El mtodo de Crank-Nicholson es un
mtodo de resolucin de ecuaciones diferenciales de una variable espacial con la variable
temporal usando la metodologa de diferencias finitas. Para esto se generan o se crean nodos
3-19
en ciertas regiones del dominio, luego se discretiza la ecuacin diferencial usando derivadas
centradas de segundo orden para la variable espacial y derivadas adelantadas para la temporal.
Al discretizar las derivadas parciales stas quedan de la siguiente forma:
C Cin+1 Cin
= (3.41)
C C in++11 C in+11 C in+1 C in1
= + (3.42)
U 4 U 4 U
2 C C in++11 2 C in +1 + C in+11 C in+1 2 C in + C in11
= + (3.43)
U 2 2 (U ) 2 (U )
2 2
Donde i es el nodo radial en el cual se est evaluando la derivada. Se observa que las
derivadas se evalan para un tiempo n y para un tiempo n+1, el promedio de estos dos tiempos
ser la interpolacin lineal entre ellos.
Al realizar la sustitucin de las Ecs. (3.41), (3.42) y (3.43) en la Ec. (3.38) se obtendr
la ecuacin de derivadas parciales discretizadas.
Cin+1 Cin 1 Cin++11 Cin+11 Cin+1 Cin1 Cin++11 2 Cin+1 + Cin+11 Cin+1 2 Cin + Cin11
= + + + (3.44)
U 2 U 2 U 2 (U ) 2 (U )
2 2
Se agrupan los trminos del tiempo n+1 (tiempo desconocido) al lado izquierdo de la ec.3.44
y el tiempo n (tiempo conocido) al lado derecho.
Donde:
1 1
A= (3.46)
2 U U 2 (U )2
1 1
B= + (3.47)
(U )2
3-20
1 1
C= (3.48)
2 U U 2 (U )2
1 1
D= + (3.49)
2 U U 2 (U)2
1 1
E= (3.50)
(U )2
1 1
F= + (3.51)
2 U U 2 (U )2
r 3 (rn R )
4 4
V =
3
(3.52)
3 3
Por otra parte, en los estudios de comportamiento de fase se usa un factor llamado
factor de hinchamiento ( k sw ). Este factor es calculado mediante la Ec. (3.12), y a partir de ste
se determina el volumen terico mediante la Ec. (3.11) quedando de la siguiente manera:
K 1
Vc = V0 + sw 4 rn3 c sat C R 2 dR (3.57)
c sat
la gota colgante dinmica a cualquier tiempo t , al conocer los valores del D y k sw empleando
la Ec. (3.45). Sin embargo, estas ecuaciones poseen un error, ya que a medida que se forma la
gota colgante hay difusin; por lo tanto, la cantidad de materia calculada no es la real. Para
disminuir este error se toma un tiempo llamado tf, que se define como el tiempo que se tarda
en formar la gota colgante. Entonces la verdadera cantidad de masa que se difundi dentro de
3-22
la gota desde que se toma la primera imagen est representada por n (t+ tf) n (tf), lo que
genera un ligero cambio en las ecuaciones anteriores.
k 1 D(t + t f ) Dt f 2
Vc (t + t f ) = V0 + sw 4 rn3 c sat C R,
C R, 2
R dR
2
(3.58)
sat
c rn rn
D(t + t f ) D t 2
Q(t , D) = 4 rn3 C R, C R, 2 f R dR (3.59)
2
rn rn
c 1 c 1 c 1 2c
= DAB 2 r 2 + 2 + 2 2 (3.61)
t r r r r sen( ) r r
c r t
C= R= =
c sat rmax 2
r
max
(3.62)
D
C 2 C 2C 1 C 1 2C
= + 2+ 2 + 2 2 (3.63)
R R R R sen( ) R R
n+1
C Ci, j Ci, j
n
= (3.64)
n +1 n +1
C Ci +1, j Ci 1. j
= (3.65)
R 2 R
n +1 n +1 n +1
2C Ci +1, j 2 Ci , j + Ci 1, j
= (3.66)
R 2 (R )2
C C i , j +1 C i. j 1
n n
= (3.67)
2
2 C C i , j +1 2 C i , j + C i , j 1
n n n
= (3.68)
2 ( )2
n+1
C Ci, j Ci, j
n
= (3.69)
C Ci +1, j Ci 1. j
n n
= (3.70)
R 2 R
2C Ci +1, j 2 Ci , j + Ci 1, j
n n n
= (3.71)
R 2 (R )2
n +1 n +1
C C i , j +1 C i. j 1
= (3.72)
2
n +1 n +1 n +1
2 C C i , j +1 2 C i , j + C i , j 1
= (3.73)
2 ( )2
Donde
3-25
1 1
A= 2 (3.75)
R R R
1 2
B= + (3.76)
(R)2
1 1
C= 2 (3.77)
R R R
1 1
D= (3.78)
R 2 R tan( )
2 2 2
1 2
E= 2 (3.79)
R 2
1 1
F= + (3.80)
R 2 R tan()
2 2 2
Para el caso de la discretizacin del radio en un tiempo conocido tenemos los siguientes
coeficientes:
1 1
A= + (3.81)
R 2 R tan( )
2 2 2
1 2
B= + 2 (3.82)
R 2
1 1
C= (3.83)
R 2 R tan( )
2 2 2
1 1
D= + 2 (3.84)
R R R
1 2
E= (3.85)
(R)2
1 1
F= + 2 (3.86)
R R R
Para simplificar la resolucin de este modelo se puede observar que tanto del lado
derecho como del lado izquierdo de la gota el proceso que ocurre es exactamente el mismo
(figura simtrica), por lo que se puede aprovechar esta simetra y resolver el modelo slo para
3-26
un lado (esto es resolverlo de un ngulo desde 0 radianes hasta radianes), ya que el otro lado
ser exactamente igual.
Con el perfil de concentraciones calculado se debe determinar el volumen terico de la
gota, el cual se compara con el volumen experimental (con la funcin objetivo), hasta lograr
minimizar la funcin objetivo dependiendo del grado de tolerancia que se le defina a dichos
volmenes.
Para calcular el volumen terico se necesita conocer la cantidad de moles que se
encuentran difundidos dentro de la gota en un intervalo de tiempo dado. Para esto se conoce
que:
dV = 2 r 2 sen( ) dr d (3.88)
Por otro lado, se conoce que la concentracin de una sustancia se define como la
relacin entre el nmero de moles y el volumen que estos ocupan; por lo tanto, al plantearla en
forma diferencial se tiene:
dn(t )
c= dn(t ) = c dV (3.89)
dV
n(t ) = 2 rmax
3
c sat C R 2 sen( ) dR d (3.92)
k 1
Vc = V0 + sw 2 rmax
3
c sat C R 2 sen( ) dR d (3.93)
c sat
Vc = V0 + (k sw 1) 2 rmax
3
C R 2 sen( ) dR d (3.94)
D(t + t f ) D t f 2
Vc (t + t f ) = V0 + (k sw 1) 2 rmax
3
C R, C R, 2
R .sen( ).dR.d
(3.95)
2
rn rn
Ecuaciones empricas
A continuacin se presentan algunas ecuaciones empricas para el clculo del
coeficiente de difusin:
k BT
D= (3.96)
6 R0
g.cm
Donde k B es la constante de Boltzmann (1,38x10-16 ), es la viscosidad a las
s 2 .K
condiciones de trabajo y R0 el radio molecular del soluto.
3-28
La viscosidad para altas presiones se determina por la ecuacin emprica de Lucas (Raid et al.,
1977) descrita a continuacin:
Pr
A
1+ D
2.118 (3.97)
=
sl 1 + CPr
g
Donde es la viscosidad del solvente a la presin P en , sl es la viscosidad a la
cm.s
g P Pvp
presin de vapor (a 25 C sl = 0, 0127 ), Pr es la relacin de presin de , Pc
cm.s Pc
Pvp =0,0317 bar), es el factor acntrico (para el agua = 0,344 ), Pr es la presin reducida
4, 674 x104
A = 0,9991
(1, 0523Tr 1, 0513)
0,03877
(3.98)
0,3257
D= 0, 2086 (3.99)
(1, 0039 T 2,573 )0,2906
r
C = 0, 07921 + 2,1616Tr 13, 4040Tr 2 + 44,1706Tr 3 84,8291Tr4 + 96,1209Tr5
(3.100)
59,8127Tr6 + 15, 6719Tr7
8, 2 x108 T 3vlq
2/3
D= 1+ (3.101)
lq v1/3 vgas
gas
3-29
cm3
Donde T es la temperatura del sistema en Kelvin, v es el volumen molar en
g .mol
(para el agua es de 18 y para el dixido de carbono es de 34).
Donde M es el peso molecular del lquido, vgas el volumen molar del gas, T
RESULTADOS Y DISCUSIN
Los dos modelos numricos para la transferencia de masa antes desarrollados fueron
evaluados por dos programas independientes, utilizando Matlab como herramienta de
cdigo. Cada modelo numrico fue evaluado con el mismo nmero de nodos (120 nodos),
valor de = 500 y =0.001, durante un tiempo mximo de 3600 s.
La gota inicialmente se forma de una sustancia pura; por lo tanto, se toma como
condicin inicial que la concentracin de dixido de carbono en el bulto de la gota en el
tiempo inicial es nula (igual a cero), para la condicin de borde se asume que toda la gota se
encuentra saturada de dixido de carbono en la interfase, por ser la superficie donde ocurre
primeramente la transferencia de masa. Tambin se supone que por simetra de la esfera, la
difusin slo se ve afectada en las coordenadas radiales y angulares, dependiendo del modelo
numrico de transferencia de masa. Adems que el cambio del volumen de la gota, no afecta al
perfil de concentracin del CO2 en sta (condicin asumida por Yang y Gu, 2005). Al analizar
la gota colgante se particulariza la condicin de borde a que no existe difusin en las paredes
del capilar al realizarse los clculos numricos en esta zona de la imagen de la gota. Por otro
lado, se asume que la sustancia que est dentro del capilar no difunde hacia la gota esfrica.
Para comparar los resultados de los modelos numricos de transferencia de masa
desarrollados en este trabajo, se emplearon los datos de difusividad y factor de hinchamiento
obtenidos por Yang y Gu (2005).
poder disminuir an ms el error, la gota se debe formar con geometra esfrica, como lo
muestra la Fig. 3.6.
Los datos experimentales del sistema CO2-agua utilizados para la validacin del mtodo
numrico radial y angular empleado en el clculo del coeficiente de difusin desarrollado en
este trabajo fueron de Hebach et al. (2002, 2005), a las condiciones de temperatura de 25 C y
presin de 6,6 MPa, los cuales son: cambio del volumen de la gota de agua con el tiempo (Fig.
3.7), dimetro interno del capilar (0,25 mm) y dimetro externo (1,60 mm), volumen inicial de
la gota (0,04015 cm3).
Figura 3.8 Cambio del volumen con el tiempo para el sistema CO2-agua a condiciones
atmosfricas.
Tabla 3.3 Resultados del la difusividad para el sistema CO2-agua a las
condiciones de P=6,6 MPa y T=25 C.
Correlaciones empricas
Modelo propuesto en Wilke y
Stokes-Einstein Sutherland
este trabajo Chang
2
D (m ) 1,016x10-9 1,094x10-9 1,762x10-9 1,674x10-9
s
% de
- 7,102 42,323 39,315
Desviacin
3-35
Comparando el valor obtenido del coeficiente de difusin de este trabajo con los
reportados en las bibliografas (Cussler, 1984 y Raid et al., 1977) se puede observar que
presenta un error del 9,89 % respecto al valor reportado por Cussler, (1984) y un error de 14
% respecto al valor reportado por Reid et al. (1977). Los porcentajes de error que se presentan
con respecto a la bibliografa se deben al movimiento constante de la interfase gota de agua,
adems que no se lleg al estado de saturacin de CO2, debido a que el tiempo que transcurri
el experimento es menor al tiempo requerido para que el sistema llegue al estado de equilibrio.
Los valores experimentales de la difusividad reportados en la tabla 3.5 para el sistema
3-36
CONCLUSIONES
El modelo numrico de transferencia de masa radial puede ser til para obtener un valor
de la difusividad como una aproximacin inicial, mientras que el modelo numrico de
transferencia de masa radial y angular es un modelo ms exacto y preciso en el clculo del
coeficiente de difusin (difusividad), a pesar de la aproximacin de geometra de la gota
colgante a una esfera simtrica, y de la suposicin del perfil de concentracin constante para
cualquier cambio de volumen en dicha esfera con el tiempo.
Al formar la gota ms esfrica que alargada la desviacin del segundo mtodo numrico
disminuye considerablemente, ya que el modelo se ajusta ms a los cambios experimentales.
3-39
RECOMENDACIONES
Para evitar dispersin en la obtencin de los valores experimental del cambio del
volumen en el tiempo se recomienda cambiar la vlvula reguladora de presin por un de
mayor precisin, con el fin de evitar fluctuaciones en la presin las cuales se traducen en una
variacin del volumen de la gota.
3-40
REFERENCIAS BIBLIOGRAFICAS
Aguilera M. E.; Cerro R. L. and A. Lpez de Ramos, Enhanced CO2 Diffusion in Wedges.
Chem. Eng. J., 31-40, 2002.
Boustani A.; Maini B., The Role of Diffusion and Convective Dispersion in Vapour
Extraction Process, J. Can. Pet. Technol., 40, 68-77, 2001.
Civan, F. y Rasmussen, M.L., "Improved Measurement of Gas Diffusivity for Miscible Gas
Flooding Under Nonequilibrium vs. Equilibrium Conditions," SPE 75135, SPE/DOE
Thirteenth Symposium on Improved Oil Recovery held in Tulsa, Oklahoma,13-17 April 2002.
Cussler E.L., Diffusion, Mass Transfer in Fluid System, First Ed., The Press Syndicate of
the University of Cambridge, New York 1984.
Das S. K.; Butler R. M, Mechanism of the Vapor Extraction Process for Heavy Oil and
Bitumen, J. Pet. Sci. Eng., 21, 43-59, 1998.
Grogan A.; Pinczewski V., Diffusion of Carbon Dioxide at Reservoir Conditions: Models
and Measurements, Society of Petroleum Engineers, SPE-14897, 1986.
3-41
Hebach A.; Oberhof A.; Dahmen N.; Kgel A.; Ederer H. and Dinjus E., Interfacial Tension
at Elevated Pressures Measurements and Correlations in the Water + Carbon Dioxide
System, J. Chem. Eng. Data, 47, 1540-1546, 2002.
Hebach A.; Martin G.; Kgel A., and Dahmen N.,Interfacial Tension during Mass Transfer of
CO2 into Water in a Water-Saturated CO2 Atmosphere at 298 K and 6.6 MPa , J. Chem.
Eng. Data, 50, 403-411, 2005.
Jimnez A., Melndez J., Zeppieri S. y Lpez de Ramos A., Modelo simplificado para la
estimacin del coeficiente de difusin mediante la tcnica de la gota colgante. IX Congreso
Internacional de Mtodos Numricos y Ciencias Aplicadas. Margarita, Venezuela, Febrero
2008.
Reid R.; Sherwood T; Prausnitz J., Properties of Gases and Liquids, 3ra ed. New York,
McGraw-Hill, 1977.
Tharanivasan, A.T, Yang, C., y Gu, Y.. Comparison of three different interface mass transfer
models used in the experimental measurement of solvent diffusivity in heavy oil. University
of Regina. J. Pet. Sci. Eng., 44 (3-4 y 15), Pages 269-282, 2004.
Tharanivasan, A.T, Yang, C., y Gu, Y.. Measurements of Molecular Diffusion Coefficients
of Carbon Dioxide, Methane, and Propane in Heavy Oil under Reservoir Conditions. Energy
Fuels, 20, (6), 25092517, 2006.
Wait, R.; Mitchell R., Finite Element Analisys and Aplication, Jhon Wiley & Sons,
Chichester, UK., 2000.
Yang C.; Gu Y., Modeling of the Adsorption Kinetics of Surfactants at the Liquid-Fluid
Interface of a Pendant Drop, Langmuir, 20, 2503-2511, 2004.
3-42
Yang C.; Gu Y., New Experimental Method for Measuring Gas Diffusivity in Heavy Oil by
the Dynamic Pendant Drop Volume Analysis (DPDVA), Ind. Eng. Chem. Res., 44, 4474
4483, 2005.
Zhang Y. P.; Hyndman C.L.; Maini B. B., Measurement of gas Diffusivity in Heavy Oils, J.
Pet. Sci. & Eng., 25, 37-47, 2000.
CAPITULO 4
INTRODUCCIN
MARCO TEORICO
Figura 4.1 Patrn de Flujo Burbuja, observado por Gayn et al. (2003), para una tubera de
4.445 cm de dimetro y QG=0,25 L/s y QL=1,88 L/s.
Flujo churn: La transicin de flujo Slug a Churn tiene lugar cuando aumenta la
cantidad de gas inyectado a la tubera (Fig. 4.4), lo que produce la aceleracin de la nariz del
Slug de gas, la cual se comienza a deformar debido al choque que se produce con las burbujas
no asociadas a la burbuja de Taylor.
Figura 4.3 Deformacin de la cabeza de la Burbuja de Taylor observada por Gayn et al.
(2003), para una tubera de 1,905 cm de dimetro.
4-5
Flujo anular: se caracteriza por la abundancia de la fase gas. El gas se desplaza como
una gran columna que ocupa casi toda la seccin de tubera rodeado por una fina pelcula de
lquido, la cual asciende adyacente a la pared.
Figura 4.5 Patrn de Flujo Anular observado por Gayn et al. (2003), para una tubera
de 3,175 cm de dimetro y QG=2,9 L/s y QL=0,46 L/s.
Debido a que las fuerzas de flotacin no actan en la direccin del flujo, en las
primeras investigaciones de flujo slug, se supona que la velocidad drift era nula para el caso
de flujo horizontal, sin embargo posteriormente se demostr que s esta presente, e inclusive
en algunas ocasiones puede ser mayor que para el flujo vertical (Anselmi, 2002). Se cree que
esta velocidad en el flujo horizontal es inducida por la diferencia de elevacin en la nariz de la
burbuja (Taitel y Barnea, 1990).
Hasan y Kabir (1992), propusieron una correlacin para la estimacin de la velocidad de
ascenso de la Burbuja de Taylor sobre lquido estancado para tuberas inclinadas, que
bsicamente es la misma propuesta anteriormente para tubera vertical, pero con la adicin de
un factor de correccin que es funcin del ngulo de inclinacin respecto de la horizontal. Esta
expresin muestra buenos resultados en el intervalo de inclinacin de 30 a 90.
4-7
0.25
D g
VBR = 1, 4 1 B (4.6)
L
2
D
3
DB = (4.7)
g
VM VGLS (1 H LLS )
VLLS = (4.9)
H LLS
0.1
V
VSG T
H GLS = 10VSG 1
(4.10)
1, 2VM + VBR (1 + cos )0.35
donde VT es la velocidad superficial de transicin entre flujo Burbuja y flujo Slug, definida
por Barnea (1987) como funcin de la velocidad superficial de lquido (VSL) y la velocidad de
ascenso de las burbujas sobre el lquido estancado(VBR) en la siguiente correlacin:
Combinando las ecuaciones (4.12), (4.13) y (4.14), se obtiene la siguiente expresin para la
estimacin del Holdup de lquido en la Burbuja de Taylor:
Y el Holdup de lquido promedio en toda la unidad del Slug puede ser estimado por la
expresin propuesta por Barnea (1990).
4-10
afectado bsicamente por la combinacin de estas variables. A medida que los Slug son ms
cortos o las velocidades de ascenso son mayores, aumenta la frecuencia de aparicin de stos.
La zona de inclinacin (50-70) reportado por Van Hout et al. (2003) como la de mayor
frecuencia de Slug, est relacionada a su vez, con la mayor velocidad y menor longitud de las
Burbujas de Taylor.
Para la estimacin de la frecuencia del Slug, Fernndez et al. (1983) propusieron una
correlacin que se basa en la longitud de la unidad del Slug y de la velocidad de ascenso de la
Burbuja de Taylor, de la siguiente forma:
VTB VTB (1 )
f = = (4.17)
L LS
En el modelo original de Taitel et al. (1980), los autores trabajaron con longitudes
promedios de Slug lquido de 20D. De esta manera, insertando la relacin de la velocidad de la
Burbuja de Taylor suministrada por Hasan y Kabir (1992), para el caso de tubera vertical la
ecuacin se expresa de la siguiente manera:
donde se obtiene del balance de masa para la unidad del Slug, Ec. (4.12), obtenindose la
siguiente expresin:
VLLS H LLS VSL
= (4.19)
VLLS H LLS + VLTB H LTB
S S
h A A A L G (
L L - G G - S 1 + 1 + - gsen
i iA )
L = L G G L (4.20)
x u 2 u 2
F
( )
- gcos-S L L + G G
L G i A
G
A
L
h h h
2
dp dp dp dp
= + + (4.22)
dz TP dz Gr dz Fr dz AC
dp
= [ H G .G + (1 H G ). L ] .g (4.23)
dz Gr
y el componente friccional generalmente se expresa como (considerando una pseudofase):
dp .VM 2
= fD . (4.24)
dz Fr 2.Dt
Por ltimo, el componente aceleracional se desprecia si y slo si el clculo de la cada
de presin se realiza en una seccin axial de tubera pequea, en la cual las propiedades de los
fluidos se consideren prcticamente constantes. No obstante, el efecto de las propiedades de
los fluidos se toma en cuenta para secciones de tuberas mayores.
Para este modelaje se consider el flujo slug como una celda que contiene la burbuja
de Taylor, la pelcula descendente de lquido y el slug de lquido (Fig. 4.6).
La burbuja de Taylor se mueve a una velocidad VTB, la longitud de stas LTB permanece
prcticamente constante si se considera el efecto expansivo despreciable (para columnas
relativamente pequeas o sistemas de alta presin). El gas dentro de la burbuja de Taylor se
supone que est a presin constante, por esta razn no se considera esfuerzo interfacial entre la
burbuja y la pelcula descendente que fluye en cada libre. El slug de lquido fluye detrs de la
burbuja de Taylor conteniendo pequeas burbujas de gas dispersas, las cuales se distribuyen
uniformemente, excepto en la zona adyacente a la burbuja de Taylor donde el holdup de gas
HGH es considerablemente ms alto que en el resto HGLS. Esta alta concentracin de burbujas
es producto del fenmeno de entrainment desde el fondo de la burbuja de Taylor debido al
efecto de la pelcula descendente de lquido. La longitud del slug de lquido LLS permanece
constante.
4-15
L QLTB
PA = (VLTB + VLLS ) = LVLTB (1 H GTB ) . (VLTB + VLLS ) (4.25)
A
2 f LS LVLLS
2
Pf = LLS = 40 f LS LSVLLS
2
(4.26)
Dt
donde LLS se supone igual a 20D.
La primera ecuacin que se plantea representa el holdup en una unidad slug HGSU:
La velocidad del gas en el slug de lquido VGLS se relaciona con la velocidad del lquido
en el slug de lquido suponiendo un comportamiento de flujo burbuja:
g ( L G ) 1/4 1/ 2
VGLS = VLLS + 1,53( ) H GLS (4.34)
L2
La pelcula de lquido descendente, como se mencion antes, no presenta esfuerzo
interfacial. Brotz (1954), citado por Fulford (1964), present una correlacin emprica para el
clculo del espesor de la pelcula para sistemas con agua, pentadecano y aceite refrigerante y
nmero de Reynolds entre 100 y 4300:
1/3 1/3
g 3Re 2F
L 2 = (4.35)
L 590
El nmero de Reynolds de la pelcula de lquido se determina de la siguiente manera:
LVLTB
Re F = (4.36)
L
La velocidad promedio de la pelcula de lquido descendente VLTB se calcula a partir de
las dos ecuaciones anteriores:
L =
Dt
2
(1 H GTB
0,5
) (4.38)
QA + QB = QC (4.39)
donde QA es el flujo de gas desde el slug superior de lquido que no puede entrar a la pelcula
de lquido descendente y coalesce con la burbuja de Taylor, QB es el flujo de gas que se
introduce en la burbuja de Taylor desde el fondo y QC es el flujo de gas desde la burbuja de
Taylor desprendido por el fenmeno de entrainment (Fig. 4.7).
La correlacin tiene los valores de C2 = 0,425 y C3 = 2,65. Una vez que se resuelve el
sistema de ecuaciones no lineales, es posible calcular los componentes de la cada de presin.
El componente aceleracional de la cada de presin es debido al cambio de direccin y
posterior aceleracin de la pelcula de lquido descendente. Sylvester (1987) expresa esta cada
de presin de la siguiente manera:
( P ) A = L (VLTB + VTB )(1 H GTB )(VLTB + VTB + VLLS ) (4.41)
GU TB 2 fTB .LS
( P ) FTB = (4.43)
2 DH
donde,
DH = Dt 2 L (4.44)
L = ( Dt / 2 ) (1 H GTB
1/ 2
) (4.45)
1
fTB = (4.46)
(1 H GTB )
2
1/2
2, 0 log
7, 4
( P ) FLS = (4.47)
2D
El factor de friccin se determina por la correlacin de Zigrang-Sylvester (1982):
1
f LS = 2
(4.48)
/ Dt 5, 02 / Dt 13
2, 0 log log +
3, 7 Re LS 3, 7 Re LS
donde:
4-19
L (1 H GLS ) .VLLS Dt
Re LS = (4.49)
LS
LS L (1 H GLS ) (4.50)
Con las ecuaciones anteriores se puede estimar la cada de presin. Sin embargo,
primero es necesario conocer el valor de LS, el cual depende de las propiedades de los fluidos,
del dimetro de la tubera y de las velocidades superficiales. Sylvester (1987) propuso el valor
de 40D, ya que con ste se ajusta mejor el modelo a los datos experimentales de cada de
presin.
4-20
MARCO EXPERIMENTAL
Para llevar a cabo el desarrollo de los experimentos se cont con un sistema bifsico de
agua-aire, un equipo por donde circula el sistema y por ltimo una seccin de visualizacin
que permite estudiar el comportamiento y las caractersticas del fluido.
El banco consta de tres tuberas de 0,75; 1,25 y 1,75 pulgadas de dimetro interno.
Cada una formada por segmentos de acrlico de 1,83 m de longitud, que se unen mediante
juntas de polietileno para tuberas lisas. La armadura est soportada en dos puntos por una
base anclada en el suelo, hecha de tubos estructurales. Uno de los puntos es fijo y constituye el
pivote, el cual permite la modificacin del ngulo de inclinacin en un intervalo de 0 a 90.
4-21
La mezcla bifsica que estudiada estuvo formada por agua y aire. Este tipo de mezcla puede
ser observada claramente con el uso de la cmara de visualizacin, lo que constituye una
ventaja al momento de realizar los clculos.
En la Tabla 4.1 se muestran la viscosidad y la densidad a condiciones estndares (Presin
atmosfrica y 25C).
Las imgenes fueron filmadas con la cmara a una velocidad de 2250 cuadros por
segundo (fps) para la medicin de velocidad, longitud de la unidad de Slug y espesor de la
pelcula; y se fij a una velocidad de 60 fps para la medicin de frecuencia. Se realiza de
forma simultnea la medicin de cada de presin.
4-22
RESULTADOS Y DISCUSIN
En los resultados se observa como las velocidades de ascenso ms elevadas para cada
ngulo corresponden a las pruebas de mayor relacin QG/QL dicha relacin pertenece al
sistema VSG=91,32 cm/s y VSL=48,50 cm/s. En trminos de la expresin propuesta por Nickiln
et al. (1962), VTB = CVM (Q G , Q L , D) + VD (D, ) , el aumento de los caudales de las fases, se
traduce directamente en el incremento de la velocidad de ascenso de la Burbuja de Taylor. Se
observa adems una desviacin porcentual entre las velocidades de ascenso de la TB, para
cada prueba y ngulo, que oscila entre 44% y el 8%. Esto se debe a que los slug tomados no
fueron lo suficientemente largos como para haber desarrollado un perfil de velocidades que
permitiese un comportamiento ms estable. (Taitel y Barnea, 1993). A su vez el tamao de la
muestra (nmero de slug) debe ser ms grande a fin de obtener valores ms confiables de la
velocidad de ascenso de la TB.
200
y = 1,1223x + 68,949
R2 = 0,987
VTB (cm/s)
150
100
40 60 80 100 120
VM (cm/s)
coeficiente para el caso de flujo turbulento (C=1,2), ya que las capas de lquido se desplazan
de manera uniforme y su empuje es ms efectivo.
En la Figura 4.10, se presenta la velocidad de la Burbuja de Taylor y la velocidad
superficial de la mezcla para las pruebas correspondientes a los tres caudales de flujo. En esta
se confirma la relacin lineal entre estas dos variables, y la pendiente obtenida del ajuste para
el caso de flujo horizontal es 1,12 en comparacin con un C=1,2 reportado por Nickiln et al.
(1962) para flujo turbulento. Calculando el nmero de Reynolds para la fase gas se obtienen
los siguientes valores:
Tabla 4.4.- Nmero de Reynolds para la fase gas
VSG (cm/s) VSL (cm/s) ReL ReG
15,92 56,08 1,53E+04 3,22E+04
47,87 54,82 1,52E+04 9,65E+04
91,32 48,50 1,38E+04 1,85E+05
De la Tabla 4.4 se puede observar que los valores del ReG 8000 y ReL 8000
indican que la fase gaseosa se encuentra en rgimen turbulento y que los valores reportados en
la literatura para la constante C corresponden a los obtenidos experimentalmente.
Para el caso de flujo vertical e inclinado se obtuvo valores para dicha constante C entre
1,17-1,3. Estos resultados fueros ordenados, tabulados y a continuacin se presentan en la
Tabla 4.5.
Tabla 4.5 Valores de la constante C para diferentes ngulos de
inclinacin de la tubera
() C
0 1,12
10 1,44
20 1,14
30 1,13
40 1,06
50 0,89
60 1,17
70 1,21
80 1,21
85 1,29
90 1,35
Con los datos reportados en la Tabla 4.5 se obtiene que el valor promedio de la
constante C, para los distintos ngulos de inclinacin de la tubera, es de 1,19 con una
desviacin porcentual entre los valores de 12%.
4-26
350
300
250
VTB (cm/s)
200
Caudal 1
150 Caudal 2
Caudal 3
100
50
0
0 10 20 30 40 50 60 70 80 90 100
()
En la Figura 4.11 se puede observar que la VTB aumenta linealmente a medida que
aumenta el ngulo de inclinacin entre 0 y 60 y luego permanece constante para ngulos
entre 70-90. Se encuentra un mximo entre 45-60, esto se debe a que la velocidad de
flotacin (otra componente de la velocidad de ascenso) tiene un mximo entre estos ngulos.
De la Ec. (4.3) se sabe que la Velocidad drift (VD) es directamente proporcional al ngulo, de
all el comportamiento descrito anteriormente. No se toma en consideracin la velocidad
superficial de la mezcla debido a que esta componente se mantiene constante.
Los valores predichos por el modelo de Hasan y Kabir se calculan utilizando las Ecs.
(4.3) (para estimar la velocidad drift) y (4.4) (para estimar la velocidad de ascenso de la
burbuja de Taylor). Los valores calculados y obtenidos experimentalmente se comparan a
continuacin en un grfico de dispersin XY para el caudal 6 de estudio, Fig. 4.12. La
desviacin porcentual est alrededor del 10%-15%, encontrndose el mayor error entre los
4-27
ngulos de 0 y 30, esto se debe a que el modelo predice con mayor exactitud para flujos
verticales
300
250
VTB (cm/s)
200
150
0 10 20 30 40 50 60 70 80 90
() Valores Experimentales
Modelo de Hasan y Kabir
Los valores experimentales son mayores a los valores estimados por la correlacin. Esto
puede ser consecuencia de que las burbujas analizadas no fueron lo suficientemente largas
como para que hayan alcanzado un perfil de velocidades uniforme. Las velocidades calculadas
empricamente cuando no se ha estabilizado el flujo son mayores a las obtenidas si el perfil de
velocidades fuese uniforme a las mismas condiciones de operacin-
80
60
LTB (cm)
40
20
0
200 220 240 260 280 300
VTB (cm/s)
verific al trabajar con los diferentes caudales (1, 2 y 6), ya que las longitudes de las burbujas
de Taylor que se obtuvieron, presentan una desviacin promedio del 55% para cada ngulo de
estudio. Sin embargo cuando el flujo de aire es aproximadamente igual al caudal de lquido
(QL/QG 1) las longitudes tienden a estabilizarse y la desviacin porcentual entre los valores
de cada prueba para un ngulo disminuye notablemente. De esta manera para una
VSG=47,87cm/s y VSL =54,82cm/s la desviacin porcentual es del 35%.
Al estudiar la longitud promedio de las burbujas de Taylor de los caudales 1, 2 y 6
variando el ngulo de inclinacin (Fig. 4.16) se puede apreciar que las menores longitudes se
consiguen cuando la VSG=15,92cm/s y VSL =56,08cm/s. Esto puede ser atribuido a la fuerte
turbulencia creada por la gran cantidad de lquido presente lo cual ocasiona un continuo
rompimiento de las burbujas que no permite la coalescencia de las mismas y por ello se
dificulta la formacin de largos tapones de gas. A continuacin se presenta el comportamiento
de las LTB a diferentes inclinaciones de la tubera.
20
15
LTB (cm)
10
0
0 20 40 60 80 100
() VSG=15,92cm/s VSL=56,08cm/s
VSG=47,87cm/s VSL=54,82cm/s
VSG=91,32cm/s VSL=48,50cm/s
30
25
20
LF (cm)
15
10
0
0 20 40 60 80 100
() VSG=15,92cm/s VSL=56,08cm/s
VSG=47,87cm/s VSL=54,82cm/s
VSG=91,32cm/s VSL=48,50cm/s
Figura 4.15. Comportamiento de la Longitud del Slug de Lquido a medida que vara
el ngulo de Inclinacin de la Tubera
25
20
15
LF (cm)
10
0
0 20 40 60 80 100
() VSG=47,87cm/s VSL=54,82cm/s
la horizontal disminua la longitud de la burbuja de Taylor y por ende disminua el valor del
factor (Figura 4.19).
60
50
40
LU (cm)
30
20
10
0
0 20 40 60 80 100
VSG=15,92cm/s VSL=56,08cm/s
() VSG=47,87cm/s VSL=54,82cm/s
VSG=)91,32 cm/s VSL=48,50cm/S
1,0
0,8
= L T B /L u
0,6
0,4
0,2
0,0
0 20 40 60 80
()
Figura 4.19 Comportamiento del factor para el Caudal 6 a medida que vara el ngulo de
Inclinacin de la Tubera
Comparando las Figs. 4.14 y 4.20 para los valores de LTB y F correspondientes al
caudal 6 (VSG=91,32cm/s y VSL=48,50cm/s), se observ que a medida que la longitud de la
burbuja de Taylor iba disminuyendo aumentaba la frecuencia de aparicin. A su vez los
valores de LTB ms pequeos corresponden a los valores de Frecuencia ms altos para este
caudal (6). No obstante este comportamiento no fue el mismo para el caudal 1 y 2, debido a
4-34
que la irregularidad en los valores de LTB no permite dar una tendencia para el
comportamiento de este parmetro
4,0000
3,0000
Frecuencia (slugs/s)
Caudal 1
Caudal 2
2,0000
Caudal 3
Caudal 4
Caudal 6
1,0000
0,0000
0 10 20 30 40 50 60 70 80 90 100
()
De la Fig. 4.20 se puede observar que a medida que aumenta el ngulo de inclinacin
con respecto a la horizontal aumenta el valor de la frecuencia, obtenindose valores mnimos
entre 70-85. Al comparar las Figs. 4.20 y 4.11 se puede afirmar que al aumentar la velocidad
de ascenso de la burbuja de Taylor aumenta su frecuencia de aparicin. Dicha tendencia
corresponde a lo predicho por la literatura, en donde velocidad de ascenso es una variable
directamente proporcionales con la frecuencia de aparicin del bolsillo de gas.
1,8478
y = 0,0004x
3
2
R = 0,9391
F (slug/s)
2
0
60 70 80 90 100 110 120 130 140 150
VM (cm/s)
() k p
0 2e-5 2,4832
20 5e-6 2,666
40 2e-5 2,3767
45 9e-5 2,1088
60 0,0006 1,707
80 0,0015 1,4995
90 0,0004 1,8478
Con los valores obtenidos y reportados en la Tabla 4.8 la ecuacin 4.51 tiene la
siguiente forma:
2,1
V 1,91 VM2
Fs = 0, 00045* LS + (4.52)
VM D gD
Los resultados obtenidos dentro de este marco experimental permitieron estimar los
parmetros de la ecuacin 4.51, lo cual significa que con la ecuacin 4.52 se pueden estimar
valores de frecuencia ms cercanos que los obtenidos con ecuaciones reportadas en la
literatura
Los resultados experimentales a ngulo de inclinacin 0, se compararon con las
correlaciones de frecuencia propuestas por Greskovich y Shrier (1972), y Heywood y
Richardson (1979). En la Fig. 4.22 se muestra la comparacin realizada:
(a) (b)
Figura 4.22 Comparacin de resultados experimentales de frecuencia con estimados por
correlaciones. (a) Modelo de Greskovich y Shrier (1972); (b) Heywood y Richardson
(1979)
Para conocer cuantos nmeros de Slug era necesario tomar para que la muestra fuese
4-37
6,00 6,00
5,00 5,00
4,00 4,00
3,00
3,00
2,00
1,00 2,00
0 20 40 60 80 100 120 140 160 0 20 40 60 80 100 120 140
Ns Ns
(a) (b)
Figura 4.23. Estabilizacin de frecuencia en funcin del nmero de slug para el caudal 6.
(a) Tubera horizontal; (b) Tubera Vertical
6 Angulo 10 / 1a
Angulo 10 / 4a
5
Angulo 10 /3a
4
Fs
(slug/s) 3
2
0
0 50 100 150 200
Ns
Figura 4.24 Estabilizacin de frecuencia en funcin del nmero de slug para el distintos
caudales y un mismo ngulo de inclinacin de la tubera (=10).
4-38
4,5
4 Angulo 10 / 3a
3,5 Angulo 80 / 3a
3 Angulo 45 / 3a
Fs
2,5
(slug/s) 2
1,5
1
0,5
0
0 50 100 150 200
Ns
Figura 4.25 Estabilizacin de frecuencia en funcin del nmero de slug para distintos
ngulos y un mismo caudal (3)
espesor de la pelcula de lquido que rodea a la burbuja, se va haciendo cada vez mas delgado
a lo largo de esta (Figs. 4.26 y 4.27).
Altura de la pelcula
1,4000
1,2000
1,0000
pelicula /R
0,8000 dsup
0,6000 dinf
0,4000
0,2000
0,0000
0,0000 1,0000 2,0000 3,0000 4,0000
LB/D
90
70
45
30
Figura 4.30 Evolucin del cuerpo de la Burbuja de Taylor a medida que se inclina la tubera
para el sistema VSG=91,32cm/s y VSL=48,50cm/s
1,8000
1,6000
1,4000
1,2000
inf/R
1,0000 0
45
0,8000
90
0,6000 30
70
0,4000
0,2000
0,0000
0,0100 0,1000 1,0000 10,0000 100,0000
LB/D
(a)
1,4000
1,2000
1,0000
0
0,8000
sup/R
30
0,6000 45
70
0,4000
90
0,2000
0,0000
0,0100 0,1000 1,0000 10,0000 100,0000
LB/D
(b)
Figura 4.31. Variacin del espesor de la pelcula a ambos lados de la burbuja a medida que se
vara el ngulo de inclinacin de la tubera. (a) inf; (b) sup
2,5 2,0
2,0
1,5
1,5
inf/R
inf/R
1,0
1,0
0,5 0,5
0,0 0,0
0 1 2 3 4 5 6 0 5 10 15
VSG=47,87 cm/s VSL=54,82 cm/s VSG=15,92 cm/s VSL=56,08
LB/D VSG=15,92 cm/s VSL=56,08 cm/s LB /D VSG=47,87 cm/s VSL=54,82 cm/s
VSG=91,32 cm/s VSL=48,50 VSG=91,32 cm/s VSL=48,50 cm/s
(a) (b)
Figura 4.32 Espesor de la pelcula inferior a lo largo de la Burbuja de Taylor (a) =30; =0
Comparando los valores experimentales con los estimados por el modelo de Dukler y
Hubbard (1975) se obtiene la Fig. 4.33.
4-44
0,030 0,030
0,025 0,025
0,020 0,020
pelicula (m)
pelicula (m)
0,015 0,015
0,010 0,010
0,005 0,005
0,000 0,000
0,0 0,1 0,2 0,3 0,4 0,5 0,0 0,1 0,1 0,2 0,2 0,3
Valores Experimentales Valores Experimentales
Xf (m) Xf (m)
Modelo Cook Behnia (1997) Modelo Cook Behnia (1997)
(a) (b)
0,030 0,025
0,025
0,020
0,020
pelicu la (m)
pelicula (m)
0,015
0,015
0,010
0,010
0,005 0,005
0,000 0,000
0,0 0,0 0,0 0,1 0,1 0,1 0,1 0,1 0,0 0,1 0,2 0,3 0,4 0,5
Valores Experimentales Valores Experimentales
Xf (m) Xf (m)
Modelo Cook Behnia (1997) Modelo Cook Behnia (1997)
(c) (d)
Figura 4.33. Valores experimentales del espesor de la pelcula inferior y los estimados por el
modelo de Cook y Behnia (1997) (a)0; (b)10; (c)20; (d)90
En la Fig. 4.33 se observa que el modelo de Cook y Behnia (1997) predice valores para
el espesor de la pelcula aceptables (la desviacin porcentual entre los valores estimados y los
obtenidos experimentalmente es del 9%) siempre y cuando el ngulo de inclinacin se
encuentre en un intervalo de (0, 20). Esto se debe a que la velocidad de ascenso de la burbuja
de Taylor se incrementa a medida que aumenta el ngulo de inclinacin de la tubera (Taitel y
Barnea, 1997) y el modelo es sensible a este parmetro.
La coincidencia entre la teora y las mediciones es sensible a la escogencia de la
velocidad de ascenso de la burbuja de Taylor usada en los clculos. Una de las
recomendaciones que sugieren los investigadores (Cook y Behnia, 1997) es que la VTB sea
cercana a 1,2VM (la Vd es despreciable) para que los valores estimados sean cercanos a los
obtenidos experimentalmente. Esta condicin se presenta para flujos horizontales o
4-45
ligeramente inclinados De all que para 90 la desviacin porcentual entre los resultados
tericos y prcticos es alta.
El espesor de la pelcula calculado se comparan con espesores experimentales
promedios, por tanto es de esperarse que los resultados no estn en total coincidencia con los
calculados; y que la desviacin porcentual entre los valores estimados y los calculados
experimentalmente sea mayores a los que podran obtenerse si se tomara una unidad slug
independiente.
4-46
CONCLUSIONES
El comportamiento de los parmetros del Flujo Slug est determinado por las
caractersticas del sistema empleado para la realizacin de las pruebas experimentales
(dimetro e inclinacin de tuberas).
Para caudales altos de lquido la Longitud de la Burbuja de Taylor (LTB) y la VTB son
inversamente proporcionales (con una tendencia lineal).
Si la Burbuja de Taylor era el doble de larga que las restantes para la misma prueba, su
velocidad siempre fue mayor con respecto a las otras.
Mientras ms alto sea el caudal de lquido las burbujas de Taylor poseen menor longitud.
A medida que se aumenta el caudal de lquido las burbujas son ms delgadas, la altura del
espesor de la pelcula ms grande y el decrecimiento del mismo se lleva a cabo ms
lentamente.
RECOMENDACIONES
Para obtener valores confiables de la FTB se recomienda tomar entre 70-100 Slug por cada
prueba
Utilizar una velocidad para la cmara ms lenta y as obtener un mayor nmero de Slug e
incrementar la muestra para que sea significativa y el valor reportado de la VTB sea confiable.
4-49
REFERENCIAS BIBLIOGRAFICAS
Abdul-Majeed G., A comprehensive mechanistic model for vertical and inclined two flow.
PhD Disertation. The University of Baghdad ,1997.
Alves, I. N., Slug flow phenomena in inclined pipes, Tesis Doctoral, Universidad de Tulsa,
1991.
Ansari M. R., Dynamical behavior of slug initiation generated by short waves in two-phase
airwater stratified flow, ASME HTD 361, 289-295, 1998.
Akagawa K. y Sakaguchi T., Fluctuation of void ratio in two phase JSME 9, 104-110, 1966.
Barnea D., A unified model for predicting flor pattern trnsition for the whole range of pipe
inclinations, Internacional Journal of Multiphase Flow, 1, 1-12, 1987.
Barnea D. y Sller L., Void-fraction measurements in vertical slug flow: applications to slug
caracteristics an transition, Internacional Journal of Multiphase Flow, 15 (5), 495-504, 1989.
Barnea D., Effect of bubble shape on pressure drop calculations in vertical slug flow,
Internacional Journal of Multiphase Flow, 16 (1), 79-89, 1990.
Barnea D. y Taitel Y., A model for slug lenght distribution in gas liquid slug flow,
Internacional Journal of Multiphase Flow, 19 (5), 829-838, 1993.
Caetano E. F., Shoham O. y Brill J. P., Upward vertical two-phase flow through an annulus.
Part I. Single-phase friction factor, Taylor bubble rise velocity, and flow pattern prediction, J.
Energy Resour. Technol. 114, 1-13, 1992.
Chokshi R., Schmith Z. y Doty D., Experimental studies and the development of a
mechanistic model for two-phase flow through vertical tubing, SPE 35676. Western Regional
Meeting, Alaska, 255267, 1996.
4-50
Cook M. y Behnia M. Investigation of Film Profiles Behind Liquid Slugs in Horizontal Gas-
Liquid Pipe Flow, AIChE Journal, 43, 2180-2186, 1997.
Davies R. y taylor G., The mechanics of large bubbles rising through extended liquid and
throgh liquids in tubes, Proc. Roy Soc., A200, 375-347, 1975.
Dukler A. y Hubbard M., A model for gas-liquid slug flor in horizontal and near horizontal
tubes, Ind. Eng. Chem. Fundam. 14, 337-347, 1975.
Dumitrescu D., Stromung an Einer Luftblase im senkrechten rohr, Math. Mech, 23, 139,
1943.
Felizola H. y Shoham O., A unified model for slug flow in upward inclined pipes, J. Energy
Res. Technol., 117, 7-11, 1995.
Fernandez R., Experimental and theorical studies of isothermal upward gas-liquid flows
vertical tubes, Tesis Doctoral, Univ. of Houston, 1981.
Fernandez R., Semial R., y Dukler A., Hydrodinamic model for gas-liquid slug flor in
vertical tubes, AICHE Journal, 29 (6), 981-989, 1983.
Fulford G.D., The flow of liquids in thin films, Advances in Chemical Engineering, editado
por Drew, T.B., Hoopes, J.W., y Vermeulen, T., 5, Academic Press, 151-235, 1964.
Gayon J., Gonzalez A. y Vargas P., Uso de imgenes de video digitales para estimar el
holdup de lquido en tuberas verticales y reconocer los patrones de flujo, Miniproyecto,
Universidad Simn Bolvar, 2003.
Ghassan H., Abdul-Majeed G. y Ali M. AL Mashat, A mechanistic model for vertical and
inclined two-phase slug flow, Journal of Petroleum Science Engineering, 27, 59-67, 2000.
Gmez L. E., Shoham O y Taitel Y., Prediction of slug liquid holdup: horizontal to upward
vertical flor, Internacional Journal of Multiphase Flow, 26, 1999.
Gopal, M. y Jepson, W.P., Development of digital image analysis techniques for the study of
velocity and void profiles in slug flow, Int. J. Multiphase Flow, 23, 945965, 1997.
Gregory G. t Scott D., Correlation of liquid slug velocity and frequency in horintal cocurrent
gas-liquid flow, AICHE Journal, 15 (6), 47-55, 1969.
4-51
Greskovich E. J. y Shrier A. L., Slug Frequency in Horizontal Gas-Liquid Slug Flow, Ind.
Eng. Chem. Proc. Des. Dev., 11 (2), 317 - 318, 1972.
Hasan A. y kabir C., Two phase flow in vertical and inclined annuli, Internacional Journal
of Multiphase Flow, 18, 279-293, 1992.
Hewitt G. y Hall-Taylor N., Anular two phase flow, Pergamon Press, Oxford, 1970.
Heywood N.I. y Richardson J.F, Slug flow of air-water mixtures in a horizontal pipe
determination of liquid holdup by -ray absorption, Chemical Engineering Science, 34, 17
30, 1979.
Liu T. J., Bubble size and entrante Light effects on void development in a vertical channel
Internacional Journal of Multiphase Flow, 19, 99-113, 1993.
Moissis R. y Grifith P., Entrance effects in a two phase slug flow, Journal of Heat Transfer,
84, 29-39, 1962.
Moissis R., The transition from slug to homogeneous two phase slug flow, Journal of Heat
Transfer, 85, 366-370, 1963.
Moreno J., Gonzalez D. y Zeppieri S., Estudio experimental de parmetros caractersticos del
flujo slug para sistemas lquido-gas en tuberas inclinadas. XXII Congreso Interamericano de
Ingeniera Qumica. Buenos Aires, Argentina, Octubre 2006.
Nicklin D., Wilkes D. y Davidson J., Two phase flow in vertical tubes, Trans. Instrum.
Chem, 40, 60-68, 1962.
Orell A. y Rembrand R., A Model for Gas-Liquid Slug Flow in a Vertical Tube. Ind.
Eng. Chem. Fundam., 25 (2), 196-206, 1986.
Polonsky S., Barnea D. y Shemer L., Averaged and time-dependent characteristics of the
motion of an elongated bubble in a vertical pipe, 25, 795-812, 1999.
S
4-52
Shemer L. y Barnea D., Visualization of the instantaneous velocity profiles in gasliquid slug
flow, Physico Chem. Hydrodynam, 8, 243253, 1987.
Shosho C. y Ryan M., An experimental for prediction flow regime transition in horizontal
and near horizontal gas-liquid flow, AICHE Journal, 22, 47-55, 1976.
Shosho C. y Ryan M., An experimental study on the motion of the long bubbles in inclined
tubes, Chemical Engineering Science, 56 (6), 2191-2204, 2001.
Sylvester N. D., A Mechanistic Model for Two-Phase Vertical Slug Flow in Pipes, ASME J.
Energy Resources Tech., 109, 206-213, 1987.
Taitel Y. y Dukler A., A model for prediction flow regime transition in horizontal and near
horizontal gas-liquid flow, AICHE Journal, 22, 47-55, 1976.
Taitel Y., Barnea D. y Dukler A., Modeling flow pattern trnasitions for steady upward gas-
liquid flor in vertical tubes, AICHE Journal, 26, 245-354, 1980.
Taitel Y. y Barnea D., Two phase slug flow, Academia Press Inc., 20, 83-132, 1990.
Tony R., Chen X. y Brill J., Ratio and bridge capacitance transductor for two phase flow
measurements, 4th International Instrumentation Symposium Heldin Aurora, 3, 7-11, 1995.
Van Hout R., Barnea D. y Shemer L., Evolution of statistical of gas-liquid slug flow along
vertical pipes, International Journal of Multiphase Flow, 27, 1579-1602, 2001.
Van Hout R., Shemer L y Barnea D., Evolution of hidrodinamics and statistical parameters of
gas-liquid slug flow along vertical pipes, Chemical Engineering Science, 58, 115-133, 2003.
INTRODUCCIN
que aumenta los costos de bombeo y mantenimiento. El flujo anular, el cual es un tipo de
patrn de flujo, es conveniente para el transporte de aceites o crudos pesados, puesto que el
agua acta como lubricante, y no permite que el crudo el cual posee mayor viscosidad toque
las paredes, lo cual ocasionara una cada de presin mucho mayor, y el riesgo de corrosin
por otro.
Como todas estas variables caractersticas dependen de las propiedades de los fluidos
que se empleen, la informacin disponible es bastante escasa. Para tuberas horizontales y
verticales se han realizado diversos estudios, pero para tuberas inclinadas los estudios
realizados son sustancialmente menores. Generalmente la informacin sobre patrones de flujo
se organiza en mapas de flujo, donde cada patrn es representado en un sistema de ejes
coordenados y se emplean las velocidades superficiales de los fluidos como coordenadas.
Tambin se han implementado modelos matemticos, para predecir las transiciones
entre patrones de flujo, y se aproximan a data experimental, pero el error de estos modelos es
elevado (mayor a 30%), por lo cual se debe trabajar en este tema para mejorar las
predicciones, y esto se logra precisamente generando ms data experimental.
En este capitulo se muestra los resultados obtenidos del estudio en sistemas multifsicos
lquido- lquido, y en el cual se reporta la retencin de fases o holdup en flujo bifsico lquido-
lquido en tuberas a 90 y 45 de inclinacin, as como tambin se reportan los diferentes
patrones de flujo caractersticos los cuales fueron registrados mediante tcnicas de
visualizacin. El Apndice D contiene la memoria en congreso arbitrado titulado:
MARCO TEORICO
Patrones de Flujo
Un patrn de flujo es la distribucin que presenta un par de fluidos no miscibles entre
s cuando son transportados a lo largo de una tubera. Esta distribucin va a depender de las
propiedades que posea cada fluido, de su velocidad y la inclinacin que tenga la tubera.
Flujo Horizontal
En esta geometra se presenta la mayora de los estudios realizados en el rea de flujo
bifsico agua-aceite, motivo por el cual se tiene mucha informacin, por lo que es muy
importante emplear un criterio unificado, con el fin de facilitar la identificacin y
denominacin de patrones de flujo observados experimentalmente. Russel et al. (1959)
encontraron en su trabajo los siguientes patrones de flujo para el flujo de agua y aceite en
tuberas, que son: burbuja (B), estratificado (ST) y mezcla (M). Estos patrones se observan
para velocidades superficiales de agua de 1,08 m/s (3,55 ft/s) y a diferentes valores de Ro/w
(relacin de entrada aceite-agua).
Los estudios realizados por Charles et al. (1961) se caracterizan por las propiedades de
sus dos fluidos de trabajo, pues tienen aproximadamente la misma densidad (998 kg/m3) y
viscosidades diferentes que son de 0,00629 Pas; 0,0168 Pas y 0,065 Pas para los aceites y
0,000894 Pas para el agua. Esto se hace con la intencin de reducir al mximo el efecto de la
gravedad, para solo tomar en cuenta el efecto de la viscosidad del aceite en la formacin de
patrones de flujo, cada de presin y retencin de fase (holdup). Los patrones de flujo
observados por estos autores fueron: gotas de agua en aceite, flujo concntrico de aceite en
agua, flujo tapn de aceite en agua, burbujas de aceite en agua y gotas de aceite en agua.
El trabajo de Guzhov et al. (1973) contiene informacin bastante detallada sobre los
patrones de flujo existentes en flujo horizontal. En su trabajo la inversin de fases ocurre
cuando la fraccin de entrada del aceite (Co) es de 60%. El aceite utilizado en su trabajo posee
una densidad de 998 kg/m3 y viscosidad de 0,0218 Pas y el agua posee densidad de 998 kg/m3
y viscosidad de 0,001 Pas. Se reporta la tensin interfacial que es de 0,0448 N/m y la tubera
utilizada tiene una longitud de 18 m y dimetro de 0,0381 m. Trabajaron con intervalos de Um
entre 0,201 m/s y 1,798 m/s y Co comprendidos entre 0,1 y 0,9. Los patrones de flujo
5-4
Como se pudo observar cada autor posee una clasificacin de los patrones de flujo
observados. Para facilitar el estudio en este trabajo se necesita unificar los criterios en los
patrones de flujo observados. La clasificacin de Trallero (1995) es la ms general y completa,
por consiguiente es la que se va a utilizar. Para flujo anular, se va a utilizar la clasificacin
dada por Bannwart et al. (2004), para el caso principal, que es (A). En la tabla 5.3 se
comparan los patrones de flujo propuestos por los distintos autores.
De esta manera quedan establecidos los patrones de Trallero (1995) como
predeterminados, y su correspondencia con los patrones observados por el resto de los autores,
excepto el flujo anular donde se utiliza la nomenclatura de Bannwart et al. (2004) de flujo
anular, subpatrones de flujo anular y patrones mezcla con flujo anular.
Tabla 5.3 Comparacin entre Patrones de Flujo Horizontales Observados por Distintos
Autores
Trallero Russel et al. Charles et al. Guzhov et al. Bannwart
(1995) (1959) (1961) (1973) et al. (2004)
Estratificado
(ST) No se observa Estratificado (E)
(ST)
Estratificado
con mezcla en
(ST & MI) Mezcla (M) No se observa la interfase y (BD)
capa inferior de
agua
Emulsin de
Burbujas de aceite en agua (BE) y
(Do/w & w) Burbuja (B)
aceite en agua y capa inferior (BDE)
de agua
5-7
Gotas de agua
en aceite y Emulsin de No se
(o/w) No se observa
flujo tapn de agua en aceite observa
agua en aceite
Emulsin de
agua en aceite
y de aceite en
agua y
(Dw/o & estratificado No se
No se observa No se observa
Do/w) con mezcla en observa
la interfase y
capa inferior de
dispersin de
aceite y agua
Emulsin de
aceite en agua
y capa inferior
Gotas de de una No se
(w/o) No se observa
aceite en agua dispersin de observa
aceite en agua
y emulsin de
aceite en agua
Flujo
concntrico
No se observa No se observa No se observa Anular (A)
de aceite en
agua
Flujo Vertical
La informacin que se dispone para flujo vertical es ms escasa que para el flujo
horizontal. El trabajo de Govier et al. (1961, 1972) Corresponde a uno de los trabajos ms
antiguos, por lo cual se sospecha que sea el primer estudio realizado sobre este tema. Los
patrones de flujo observados fueron: burbuja, tapn, espuma, rizado o semianular, pelcula o
anular y niebla. La nica propiedad reportada en su trabajo es la viscosidad del aceite
utilizado, la cual es de 0,0201 Pas, en un intervalo de Usw y Uso de 0,003 m/s hasta 3,048 m/s
en ambos casos.
Zavareh et al. (1988) presenta una clasificacin de 4 patrones de flujo observados, los
cuales son: aquellos dominados por agua (burbuja y burbuja disperso) y los dominados por
aceita (burbuja invertido y burbuja disperso invertido). La longitud de la tubera utilizada es de
5-8
12,8 m, y dimetro de 1,842 m., para un intervalo de Uso entre 0,05 m/s y 0,5 m/s y Usw entre
0,02 m/s y 0,6 m/s.
El trabajo de Flores et al. (1997) es uno de los ms completos que se han realizado en
flujo vertical bifsico agua-aceite, pues se observaron diversos patrones de flujo que de alguna
manera pueden compararse con otros trabajos en esta rea. El patrn de flujo que no fue
observado es el anular. Las propiedades de los fluidos utilizados son de o/w igual a 0,85,
o/w igual a 20, tensin interfacial de 0,029 N/m y tubera de 15,5 m de longitud y 0,0508 m
de dimetro.
Dividen los patrones de flujo en dos grupos, flujo segregado y flujo disperso. El flujo
segregado conserva la continuidad de las dos fases en direccin vertical mientras que el
disperso implica que la continuidad se pierde para alguna fase. Los patrones de flujo que
observaron estos autores son de flujo disperso, algunos dominados por agua y otros por aceite,
y son seis en total, en un intervalo de Uso comprendido entre 0,2 m/s y 8 m/s al igual que Usw.
En la tabla 5.4 y en la figura 5.2 se pueden ver los patrones de flujo observados por estos
autores y sus caractersticas.
Figura 5.2 Patrones de Flujo Verticales Propuestos por Flores et al. (1997)
5-9
Tabla 5.4 Patrones de Flujo Vertical Observados por Flores et al. (1997)
Los patrones de flujo verticales observados por Bannwart et al. (2004) fueron: burbuja
B, burbuja dispersa BD, intermitente I, anular A, anular ondulado AO, anular ondulado
distorsionado AOD, anular ondulado perturbado AOP, anular perfecto AP, burbuja + burbuja
dispersa BBD, intermitente + burbuja dispersa IBD, anular ondulado distorsionado + burbuja
dispersa AODBD, anular ondulado + burbuja dispersa AOBD, anular ondulado perturbado +
burbuja dispersa AOPBD y anular perfecto + burbuja dispersa APBD. Este trabajo forma parte
del trabajo sealado en la seccin de flujo horizontal. La tubera utilizada tiene 5,43m de
longitud y 0,0284m de dimetro. Los intervalos de velocidades superficiales son los mismos
que para el caso horizontal. En la figura 5.3 se puede apreciar los patrones de flujo estudiados
por los presentes autores.
5-10
Figura 5.3 Patrones de Flujo en Tuberas Verticales Observados por Bannwart et al. (2004)
En la tabla 5.5 se presenta una comparacin entre los diferentes patrones de flujo
observados por autor, con la finalidad de establecer una nomenclatura comn entre ellos para
lograr unificar un solo criterio con el fin de facilitar el trabajo experimental.
Observando la tabla 5.5 se establecen los 6 patrones de flujo observados por Flores et
al. (1997) como predeterminados para los flujos dispersos y segregados. Para los flujos
anulares e intermitentes es conveniente agruparlos bajo un solo trmino, que son el flujo
tapn, rizado y pelcula o anular que propone Govier et al. (1961). Esto se dispone as, ya que
los subpatrones y patrones combinados de Bannwart et al. (2004) son muy especficos, y bajo
la nomenclatura de Govier et al. (1961) se hace un poco ms genrico y fcil de manejar.
Flujo Inclinado
Sin haber realizado estudios de flujo bifsico agua-aceite en tuberas inclinadas,
algunos autores pensaban que los patrones de flujo inclinados podan predecirse al combinar el
caso vertical y el horizontal. Esta hiptesis se descart por completo cuando se obtuvieron los
primeros resultados de experimentos en los que se observaron los patrones de flujo para el
sistema agua-aceite en tuberas inclinadas. Se encontraron patrones de flujo que no
correspondan con ningn caso de flujo horizontal o vertical. Tambin se pudo concluir que
los patrones de flujo variaban con el ngulo de inclinacin con respecto a la horizontal, y
dependiendo si el flujo es ascendente o descendente.
Para este trabajo interesa solamente el flujo ascendente, y la informacin que se tiene
para este caso es mucho menor que para el flujo vertical y horizontal, pues se tienen a
disposicin dos autores que tratan el tema en profundidad. Esto se debe a que el flujo
inclinado es mucho ms especfico, pues no solo va a depender de las propiedades de los
fluidos con los que se trabaje, sino tambin del ngulo de inclinacin, ya que este influye en
los patrones que se puedan observar.
Zavareh et al. (1988) trabajaron en flujo bifsico de agua-aceite tanto en tuberas
verticales como en inclinadas. Para tuberas inclinadas trabajaron con ngulos de 5 y 15
respecto al eje horizontal. En la siguiente tabla se describen los diferentes patrones de flujo
obtenidos por estos autores. El intervalo de velocidades superficiales y la tubera utilizada por
estos autores se coment en la seccin e flujo vertical. Para una tubera inclinada 5 se observ
burbuja dispersa y burbuja en contracorriente, para el caso de 15 se observ burbuja en
contracorriente y burbuja dispersa en contracorriente. En las figuras 5.4 y 5.5 se visualizan los
patrones de flujo reportados por estos autores.
5-12
Figura 5.4 Patrones de Flujo en Tuberas de 5 de Inclinacin Observados por Zavareh et al.
(1988)
Figura 5.5 Patrones de Flujo en Tuberas de 15 de Inclinacin Observados por Zavareh et al.
(1988)
Flores et al. (1997) trabajaron con ngulos de inclinacin de 75, 60 y 45. Observaron
algunos patrones de flujo ya descritos por Zavareh et al. (1988) y otros patrones de flujo
diferentes. En la seccin de flujo vertical se indicaron las propiedades de los fluidos que
utilizaron, los intervalos de velocidades superficiales y las dimensiones de la tubera.
En total son siete patrones de flujo los observados, donde cuatro son dominados por
agua, dos dominados por aceite y uno de transicin. En la tabla 5.6 se describen los patrones
de flujo observado por estos autores.
La diferencia entre los patrones de flujo que se observaron en el caso vertical con
respecto al caso inclinado, se debe a la componente gravitacional que se encuentra normal a la
direccin del flujo, sumado a la presin y fuerzas viscosas que son importantes gracias a las
altas viscosidades del aceite.
5-13
Tabla 5.6 Patrones de Flujo Inclinado Observados por Flores et al. (1997)
Figura 5.6 Patrones de Flujo Inclinado Observados por Flores et al. (1997)
5-14
La clasificacin dada por Flores et al. (1997) para patrones de flujo agua-aceite en
tuberas inclinadas es ms completo y simple en comparacin con la clasificacin que presenta
Zavareh et al. (1988), pues Flores et al. (1997) observa los mismos patrones de flujo que
Zavareh et al. (1988) y tres patrones ms. En la tabla 5.7 se establece la comparacin entre los
patrones de flujo observados por ambos grupos de autores.
Tabla 5.7 Cuadro Comparativo entre Patrones de Flujo Observados por Zavareh et al. (1988)
y Flores et al. (1997)
pues emplea los mismos patrones de flujo propuestos por Trallero (1995). A continuacin en
la tabla 5.8 se muestra un resumen de los mapas de patrones de flujo por autor, que contiene
propiedades de fluidos, dimensiones de tuberas y coordenadas empleadas segn cada caso.
(5.1)
(5.2)
Tabla 5.11 Rangos de Aplicacin para las Correlaciones de Mukherjee et al. (1985)
Rango
Variable Mnimo Mximo
Fr 0,95 1,4
Re 6495 16724
Cw 0,2 0,75
30 90
(5.3)
Tabla 5.12 Propiedades de los Fluidos Utilizados por Hasan y Kabir (1988)
MARCO EXPERIMENTAL
El equipo utilizado es el mismo empleado para el estudio de flujo bifsico gas lquido
desarrollado en el capitulo 4. Como ya se mencion antes el equipo se encuentra ubicado en el
Laboratorio de Operaciones Unitarias en la Universidad Simn Bolvar. Consiste en un banco
de tuberas soportado por una estructura metlica de 8,3 m de altura. Dicha estructura tiene la
capacidad de girar de 0 a 90, lo cual permite realizar estudios a cualquier ngulo dentro de
este intervalo en el sentido que se desee (Figura 5.7).
Los fluidos utilizados para los experimentos fueron agua y aceite. Las propiedades del
agua se obtuvieron de tablas. El aceite utilizado es Hidrovenoco S-68. El rango de temperatura
en el laboratorio varia entre 20C y 24C, y las propiedades del agua se calcularon a 23C, al
igual que la del aceite excepto la densidad que est a 21C. La tensin interfacial se obtuvo del
miniproyecto realizado por Semprn et al. La tabla 5.13 resume las propiedades.
RESULTADOS Y DISCUSIN
Patrones de Flujo
En esta seccin se presenta los patrones de flujo observados. La seccin de divide en
dos partes, una dedicada al flujo vertical y la otra al flujo inclinado. Es importante destacar
que los patrones de flujo que se mencionan a continuacin son resultado tanto de observacin
directa durante los experimentos ms la ayuda de las fotografas obtenidas de la cmara de alta
velocidad.
Flujo Vertical
Los patrones de flujo vertical observados en este trabajo son cuatro de los cuales tres
son dominados por la fase acuosa y el restante posee dos fases continuas. Los patrones
dominados por la fase acuosa son: Tapn, dispersin de aceite en agua (Do/w) y dispersin muy
fina de aceite en agua (VFDo/w), mientras que el patrn restante observado es el anular.
5-21
Para poder apreciar como van evolucionando los patrones de flujo a medida que
aumentan las velocidades superficiales de aceite y de agua, se obtuvieron fotos con la cmara
de alta velocidad fijando una velocidad de alguno de los dos fluidos y variando la otra. Este
barrido da como resultado cinco grficos en los que se puede apreciar tanto los patrones de
flujo observados como su evolucin. En la primera transicin, apreciable en las figuras 5.8,
5.9 y 5.10 se eligi una Uso fija de 0,08m/s, y Usw vara desde 0,10m/s hasta 0,64m/s
aproximadamente. En la segunda transicin se deja el agua fija a Usw de 0,24m/s y con valores
de Uso comprendidos entre 0,02 y 0,39m/s, apreciables en las figuras 5.11 y 5.12.
Figura 5.8 Patrn de Flujo Tapn a Uso =0,02m/s para valores de Usw : (A) 0,10m/s, (B)
0,13m/s, (C) 0,16m/s, (D) 0,21m/s y (E) 0,24m/s.
un tapn y otro suele estar ocupado por burbujas de aceite pequeas. Es posible ver dos
tapones seguidos sin la presencia de burbujas de aceite entre uno y otro.
Figura 5.9 Transicin Entre Patrones de Flujo Tapn y Dispersin de Aceite en Agua a
Uso=0,02m/s para valores de Usw : (A) 0,29m/s, (B) 0,32m/s, (C) 0,35m/s, (D) 0,38m/s y (E)
0,41m/s.
Se aprecia que los tapones de aceite de la figura 5.8 presentan mucha similitud en
cuanto a longitud. A medida que aumenta Usw, sera lgico esperar la disminucin del tamao
de los tapones, lo cual no es apreciable en este intervalo de velocidades que es muy pequeo, a
penas de 0,14 m/s aproximadamente.
En la figura 5.9 se puede apreciar los patrones de flujo tapn, dispersin de agua en
aceite y patrones de transicin. En los apartados (A) y (B) de la figura 5.9 se aprecian
burbujas de aceite grandes, con forma de tapn que poseen aproximadamente la mitad de la
longitud de los tapones observados en la figura 4.1, antecedidas y precedidas por burbujas de
aceite de menor tamao, esto ocurre a Usw de 0,29 m/s en (A) y 0,32 m/s en (B); de la figura
5.9. En los apartados (C), (D) y (E) de la figura 4.2 con Usw de 0,35 m/s; 0,38 m/s y 0,41m/s
respectivamente, se aprecian burbujas de aceite de gran tamao entre burbujas de aceite de
menor tamao. La transicin ocurre gracias al aumento de Usw, que ocasiona que la frecuencia
de los tapones aumente, choquen unos con otros y colapsen dando paso a burbujas de aceite
grandes y pequeas pues el aumento de Usw aumenta las fuerzas turbulentas en la tuberas,
capaz de separar los tapones.
5-23
Finalmente lo que se aprecia en la figura 5.10 es la variacin entre los patrones de flujo
dispersin de aceite en agua y dispersin muy fina de aceite en agua. Para los apartados de la
figura 5.10 (A), (B), (C), (D) y (E) con Usw de 0,48 m/s; 0,53 m/s; 0,56 m/s; 0,61 m/s y 0,64
m/s respectivamente se puede apreciar desde el primer patrn de flujo que se presentan
burbujas de aceite de tamao mediano acompaado de otras burbujas de aceite de menor
tamao. A medida que se aumenta Usw se puede ver como el tamao de las burbujas va
disminuyendo y se va haciendo ms uniforme. A partir del patrn de flujo ilustrado en el
apartado (C) se puede apreciar que ya el patrn de flujo corresponde a una dispersin muy fina
de aceite en agua.
De manera general se observa que al aumentar Usw y conservar Uso, el tamao de los
tapones va disminuyendo hasta formar burbujas de aceite de gran tamao y progresivamente
estas burbujas se van haciendo ms pequeas conforme aumenta la velocidad hasta que el
tamao de las burbujas es prcticamente constante. Tanto en el patrn dispersin de aceite en
agua y dispersin muy fina de aceite en agua la distribucin de las burbujas en el rea
transversal y longitudinal de la tubera es prcticamente constante.
Figura 5.10 Variacin del Patrn Dispersin de Aceite en Agua hacia Dispersin muy Fina de
Aceite en Agua a Uso=0,02m/s para valores de Usw : (A) 0,48m/s, (B) 0,53m/s, (C) 0,56m/s,
(D) 0,61m/s y (E) 0,64m/s.
5-24
finalmente forman una fase continua de aceite en el centro de la tubera, rodeado por la fase
acuosa.
Figura 5.12 Transicin de Patrones de Flujo en Tubera Vertical de Tapn hacia Anular a
Usw=0,24m/s para valores de Uso de: (A) 0,23m/s; (B) 0,29m/s; (C) 0,34m/s y (D) 0,39m/s.
5-26
Flujo Inclinado
Para flujo en 45 de inclinacin se observaron 4 patrones de flujo, donde la fase que
domina en todos es la acuosa. Estos patrones son: dispersin de aceite en agua-contracorriente
(CT), dispersin de aceite en agua-pseudopatrn (PS), dispersin de aceite en agua-cocorriente
(CC) y dispersin muy fina de aceite en agua (VDFo/w).
El patrn de flujo (PS) observado consiste en un cmulo de burbujas que se ubican en
la regin superior de la tubera formando una burbuja de mayor tamao, y se observa que el
agua circula en direccin contracorriente al flujo de ascenso del aceite, gracias a la
componente gravitacional, como se aprecia en la figura 5.13. Ocurre a Uso bajas y Usw de bajas
a moderadas. Al aumentar Usw a 0,56 m/s se observa el patrn de flujo (CC), que consiste en
burbujas de dimetros pequeos que viajan en cocorriente con el flujo de agua en la regin
superior de la tubera. Los cmulos observados en el patrn (PS) desaparecen, gracias al
aumento de fuerzas turbulentas causadas por el aumento de Usw. Para Usw superiores a 0,56
m/s se observa el patrn (VDFo/w), que es de apariencia similar al (CC), la nica diferencia es
que la accin de las fuerzas turbulentas causan que la dispersin de las burbujas de aceite sea
ms uniforme en la regin transversal de la tubera.
Si se aumenta Uso a 0,18 m/s manteniendo valores de Usw bajos, las burbujas de aceite
que forman el patrn (PS) coalescen y forman en la regin superior una fase continua delgada,
y en la interfase entre el aceite y el agua se encuentra una capa de burbujas muy pequeas de
aceite y el agua viaja en sentido contracorriente. Al ir aumentando Usw, manteniendo Uso a
velocidades de 0,18 m/s aproximadamente, la fase continua de aceite va desapareciendo hasta
que se puede apreciar un patrn transicin entre (CC) y (PS). Este patrn presenta tapones de
aceite de vez en cuando, y alterna con burbujas de aceite entre tapn y tapn como se explica
en el patrn (CC). Si se aumenta Usw se obtiene el patrn (VDFo/w).
5-27
Flujo Vertical
El mapa de flujo vertical obtenido esta formado por veintisis puntos experimentales,
en donde todos los patrones de flujo reportados fueron observados a simple vista. El rango de
velocidades superficiales del mapa es de 0,02 < Uso(m/s) < 0,39 y de 0,16 < Usw(m/s) <0,72.
Como se aprecia en la figura 5.14, hay zonas bien diferenciadas donde se ubican las
regiones de patrones de flujo, en la parte inferior el flujo tapn, en la regin media la
dispersin de aceite en agua, en la regin superior la dispersin muy fina de aceite en agua, y
en la regin inferior derecha el flujo anular. El patrn de flujo ms observado fue el de
dispersin muy fina de aceite en agua, seguido por el flujo tapn y la dispersin de aceite en
agua. El que menos se presenta fue el patrn de flujo anular, por las condiciones de operacin
del equipo. Tambin es importante resaltar que entre el patrn de flujo tapn y dispersin de
5-28
aceite en agua estn dos puntos del mapa que corresponden a la transicin entre estos dos
patrones de flujo. El patrn anular combinado con dispersin muy fina de aceite en agua no se
tiene registro visualizado, pues las condiciones de operacin a las que se obtiene causan
mucha inestabilidad en el equipo de flujo multifsico.
Figura 5.15 Comparacin entre el Mapa de Patrones de Govier et al. (1961) y el Obtenido
Experimentalmente a 90
5-29
Figura 5.18 Comparacin entre el Mapa de Patrones de Flujo de Bannwart et al. (2004) y el
Obtenido Experimentalmente a 90
En las figuras 5.14, 5.15, 5.16 y 5.17 se ilustra cada mapa de patrones de flujo de 90
encontrados en la bibliografa comparado con el obtenido en este trabajo. Para hacer esto, se
llev el mapa obtenido en este trabajo a la escala de cada mapa.
En la figura 5.14 se coteja el mapa de Govier et al. (1961) con el obtenido en este
trabajo. Los patrones que observan Govier et al. (1961) corresponden prcticamente a los
observados en este trabajo. Incluso la distribucin de los patrones de flujo es muy parecida,
pero los intervalos de Usw y Uso son menores. Por dar un ejemplo de esto, el flujo anular
aparece para Uso mayores a 0,30m/s, mientras que en el mapa de Govier et al. (1961), se
presenta para Uso mayores a 3m/s. Las propiedades de los fluidos y tuberas utilizados, son
diferentes a las que se emplea en este trabajo, esto justifica la diferencia en cuanto a
resultados. Si se conociera el dimetro de la tubera utilizada por Govier et al. (1961), se
pudiera concluir la razn de la discrepancia. Por poseer mayor viscosidad el aceite que se
emplea en este trabajo, es complicado separar las burbujas grandes en otras de menor tamao,
pero si el dimetro de la tubera utilizada es ms pequea que la utilizada por estos autores, las
velocidades superficiales aumentan por lo cual se pudieran vencer las fuerzas viscosas y
facilita la ruptura de burbujas grandes de aceite.
5-31
Flujo Inclinado
El mapa de flujo para 45 obtenido en este trabajo consta de once puntos
experimentales. A pesar de no tener tanta informacin como el mapa de flujo de 90, se puede
observar las regiones de los diferentes patrones de flujo de manera muy general. El rango de
velocidades superficiales estudiadas experimentalmente para generar los puntos del mapa son
de 0,02 < Uso(m/s) < 0,39 y de 0,16 < Usw(m/s) <0,72, iguales que para el caso de 90.
En la figura 5.18 se puede observar que a baja Uso y de Usw, se presenta el patrn de
flujo (PS), que es el que ms se aprecia junto con el patrn de flujo (CT) y (VFDo/w). Este se
observa hasta Usw de 0,40m/s. El patrn (VFDo/w), para Usw ms elevados, mientras que (CT)
aparece para Uso mayores a 0,1m/s, y a Usw moderadas. El patrn de flujo menos observado es
el (CC), que se aprecia entre el patrn (PS) y (VDFo/w) a Uso de 0,02m/s. El ltimo punto
corresponde a una transicin entre el patrn (CC) y (PS), ubicado entre el patrn (CT) y
(VDFo/w).
Figura 5.20 Variacin del Holdup de Agua con respecto a Usw a Uso constante, Tendencia
Terica y Experimental a 90
Anlogamente que para el caso de 90, se tiene una grfica con la tendencia que
presenta el holdup de agua, a medida que aumenta Usw, como se observa a continuacin:
5-36
Figura 5.22 Variacin del Holdup de Agua con respecto a Usw a Uso constante, Tendencia
Terica y Experimental a 45
Claramente se puede apreciar en la figura 5.22 que para Uso constante, Hw tiende a la
unidad cuando aumenta Usw. La diferencia ms notable se encuentra para Uso de 0,02m/s,
donde para valores de Usw mayores a 0,16m/s y 0,24m/s, la tendencia experimental se aleja un
poco de la terica. Lo importante es que la tendencia experimental es igual a la esperada,
ilustrada muy bien por la tendencia que predice la ecuacin de Hasan y Kabir (1988). La
ecuacin predice valores de Hw ligeramente mayores a los experimentales.
Las correlaciones obtenidas por Mukherjee et al. (1985) fueron evaluadas al igual que
en el caso de 90 para as poder conocer la desviacin que presenta con los valores
experimentales. Se repite la tendencia observada en el caso de 90, no se presentan resultados
sobre estos clculos pues las condiciones de este trabajo se encuentran fuera del rango de
aplicacin de dichas ecuaciones
La correlacin de Hasan y Kabir (1988) tambin fue evaluada para el caso de flujo
inclinado. El error relativo promedio para flujo inclinado es de 2%, menor que en el caso
vertical, por lo cual se recomienda utilizar esta ecuacin para predecir holdup en flujo aceite-
agua, bajo las condiciones de operacin del laboratorio. A continuacin en la figura 5.23 se
aprecia el diagrama de paridad que corresponde. Una vez ms se verifica que la metodologa
para medir holdup fue la adecuada.
5-37
CONCLUSIONES
Para flujo vertical se observaron los siguientes patrones de flujo: Tapn, dispersin de
aceite en agua, dispersin muy fina de aceite en agua y anular.
Para flujo vertical se observ que a medida que aumenta Usw a Uso constante, se va
desde patrones de flujo tapn hasta dispersin muy fina de agua en aceite. Cuando se aumenta
Uso a Usw constante, se va desde un patrn de flujo dispersin de aceite en agua hasta anular.
La tendencia que se observ en flujo inclinado es que al aumentar Usw a Uso constante,
los patrones que presentan dispersiones de aceite en agua con burbujas de tamao promedio o
fase continua de aceite en la regin superior de la tubera, tiende a separar las burbujas
pequeas y a interrumpir la fase continua en burbujas grandes, hasta llegar a un patrn de flujo
de dispersin muy fina de aceite en agua y al aumentar Uso a Usw fija la tendencia es ir desde
burbujas dispersas del patrn (PS) hacia (CT).
El mapa de patrones de flujo obtenido por Bannwart et al. (2004), es el que presenta
mayor similitud con el obtenido en este trabajo para flujo vertical, pese a que en el obtenido
experimentalmente, los intervalos de aparicin de los patrones de flujo fue menor. Esto se
debe a que las dimensiones de las tuberas empleadas en ambos casos son similares, pero la
viscosidad es menor en el caso de este trabajo.
El mapa de flujo obtenido para flujo inclinado se pudo comparar slo con el elaborado
por Flores et al. (1997) y se observaron todos los patrones de flujo excepto los que son
dominados por la fase de aceite y la transicin. La ubicacin en ambos mapas de los patrones
de flujo comunes son similares, pero no en los mismos rangos. Se puede decir que el mapa
5-39
obtenido en este trabajo con respecto al elaborado por Flores et al. (1997) se encuentra
desplazado hacia la izquierda y hacia abajo.
Se observ que el holdup de agua aumenta a medida que aumenta Usw a Uso constante y
tiende hacia uno y coincide con la tendencia esperada por la definicin de holdup tanto para
flujo vertical como para flujo inclinado.
La correlacin obtenida por Hasan y Kabir (1988) es aceptable para predecir holdup
tanto a 90 como a 45, pues los errores relativos promedio de los valores experimentales con
respecto a los ofrecidos por la correlacin son de 6% y 2% respectivamente. Esto implica que
la metodologa que se utiliz para medir holdup es adecuada.
5-40
RECOMENDACIONES
En lo posible reforzar la seccin de prueba para poder trabajar con mayores caudales
de aceite y as poder cubrir ms puntos experimentales, tanto para holdup como para patrones
de flujo.
REFERENCIAS BIBLIOGRAFICAS
Bannwart A., Rodrguez O. y Carvalho C. Flow Patterns in Heavy Crude Oil-Water Flow,
Journal of Energy Resources Technology, Vol. 126, 184-189, 2004.
Charles, M.E., G. W Govier. y G. W. Hodson, The horizontal pipeline flow of equal density
oil-water mixture, The Canadian Journal of Chemical Engineering, 39, 27-36 ,1961.
Govier, G. W., G. A Sullivan y R.K. Wood, The Upward Vertical Flow of Oil-Water
Mixtures, The Canadian Journal of Chemical Engineering, 39, 67-72, 1961.
Govier, G. W. y Aziz. K., The Flow of Complex Mixtures in Pipes, Van Nostrand Reinhold
Co, New York, 563, 1972.
Hasan, A. R y Kabir. C. S., A new Model for Two-Phase Oil-Water Flow: Production Log
Interpretation and Tubular calculations, Proceeding Annual Technical Conference and
Exhibition, 3, 369-376, 1988.
Matamoros L., Estudio Experimental del Flujo Bifsico gas-lquido ascendente y descendente
en tuberas verticales, Proyecto de Postgrado, Universidad Simn Bolvar, 2003.
Mukherjee, H., J.P Brill y H.D. Beggs, An Experimental Study of Oil-Water Flow in Inclined
Pipes, Journal of Energy Resources Technology, 103, 56-66, 1981.
Mukherjee H., J. Brill y H. Beggs, Experimental Study of Oil-Water Flow in Inclined Pipes,
Journal of Energy Resources Technology, ASME, Vol. 103, pp 229-215, June 1985.
5-42
Romero C., Sanchez E., Gonzlez Mendizbal D. y Zeppieri S., Estudio experimental sobre
patrones de flujo lquido-lquido en tuberas verticales. 8VO Congreso Iberoamericano de
Ingeniera Mecnica. Argentina, 2007.
Trallero J.L., Oilwater flow patterns in horizontal pipes. Tesis Doctoral, University of
Tulsa, 1995.
Transport Phenomenal Laboratory, University Simon Bolivar, Post Office Box 89.000
Superficial tension for CO2 + Water system is reported for temperature from 293 to 313 K
and pressure range of (3-5) MPa. We have measured the change of superficial tension ()
between H2O and CO2 with a pedant drop method, the experiment was realized three times
to analyze the deviation associated to the method and the experimental system, and to
every experimental data the obtained results were very next, presenting a small standard
minor diversion of 2%. We observed that () decreased with time and with the CO2
pressure. But when increase temperature to a same high pressure the superficial tension in
Introduction
Superficial tension plays an important role in the modeling and design of extraction,
cleanliness, the main factors to determine the efficiency of the mass transference are the
number of interactions between the molecules in the interface and the drop size, which is
strongly dependent on superficial tension, between two phases reflects that their forces
1
Corresponding autor. Tel.: +58 212 9064110, Fax: +58 212 9064137, E-mail: [email protected]
cohesion and interaction [1]. These forces meet affected by the pressure and temperature of
the system, diminishing this way the value of the above mentioned tension. High-pressure
interfacial and surface-tension phenomena govern the migration and recovery of oil and gas
separation and mass transfer in light hydrocarbon fractionation plants [2]. Carbon dioxide
injection is a technique widely used in enhanced oil recovery processes, and the reduction
of superficial tension plays an important role in improving the mobility of reservoir fluids
The objective of this work is to measure the change superficial tension of CO2 + Water
system with the time for different temperature T= 293, 303 and 313 K to high-pressure, in a
experimental unit of low cost in the system pressure, registered the variation the drop
volume in the time with constant pressure and temperature. And to compare the results with
Experimental Section
Materials. The vessel was filled with carbon dioxide (purity 99.9% Vol.) Which were
supplied by AGA GAS, C.A. Distilled water was distilled prior to using was produced a by
our Laboratory.
Apparatus. The apparatus shown in Figure 1 is similar to that used in the reference [3]
and the designs cell similar to the reference [4], except the pressure system has been
largely modified. The modified pressure system consists of a cylinder with a piston
designed under the norms ASMT [5] to high-pressure, where this contained phase disperses
pressurize with compressed air to calculate superficial tension inside of the high-pressure
cell [6]. The optical system consists in a high-resolution Sony photographic camera
Drop 98 [7] developed by our laboratory. This allows convenient observation and recording
of the drop forming process and the effect of drop-age on interfacial tension.
The apparatus consisted in a high-pressure cell, a light source, a CCD camera, a computer.
Pendant drops were formed at the end of a stainless steel tube 1,3mm inside diameter,
connected to a sample cylinder. The pressure inside the cell was controlled during all
the pressure in the drop during the time of experience, so that the surface tension variation
was for the effect of the pressure in the whole interface of the drop. This experimental
Needle Valve.
solution, places the water sample cylinder of the dispersed phase (3-a); to displease all the
air of the pipeline, to open the valves purge and later the Valve V-1 allowing to pass
slowly the water until to form a drop in the capillary. To check in the regulator of pressure
(P1) in the cylinder of compressed air that over the pressure (P2) of the cylinder carbon
dioxide. With the cell closed (1), to open the valve (V-2) that allows the step of the CO2 to
the cell until the wanted pressure, after to open the purge valve (V-3) to displace all the air
contained in the cell. To fix a low flow in the micrometric valve (V-1) to form the drop,
while the system arrives to stationary state to the pressure and wanted temperature,
Certainly, it is necessary to open the purge valve (V-3) slowly to prevent the cell from
filling with the dispersed phase, after reaching the above mentioned stationary condition, a
drop is formed in the cell, the contour being captured of this one in an definite time
interval, and this way to calculate the superficial tension of the system applying the Laplace
equation.
Fig. 2 Profile of drop analyzed with Laplace equation in the cylindrical coordinate system.
Calculation of superficial tension. The interfacial tension was determined the profile of the
droplet (Figure 2) using a CCD camera coupled to a video image profile digitizer board
connect to a computer. The drop profile was processed according to the fundamental
As illustrated in Fig. 2 , a pendent drop of water is formed inside a high-pressure cell that is
filled with the test solvent. The syringe needle has a radius of rn . The profile of the pendent
1 1
Po Ps = + (1)
R1 R2
Where Po is the pressure in the pendent water drop; Ps is the pressure in the solvent phase;
is the Superficial tension between the water and the solvent; R1 and R2 are the principal
elevation. For an arbitrarily chosen point P ( r , z ) on the interface, the pressure is Eq. (1)
z
Po = Po 0 g o dz (2a)
0
and
Ps = Ps 0 s gz (2b)
Where o is the density of the water solvent mixture; s is the density of the solvent [8] at
the experimental pressure and the temperature; g is the gravitational acceleration; Po 0 and
Ps 0 are the respective reference pressures of water and solvent at the origin O of the r-z
cylindrical coordinate system. Substituting Eqs. (2a 2b) into Eq. (1) and rearranging it
gives:
1 1
( Po0 Ps 0 ) gz = + (3)
R1 R2
z
Where, = o dz / z s oave s is the density difference between the Water [9]
0
and solvent; oave is the average density of the water-solvent mixture inside the pendent
drop. During the measurement, oave is a function of time as the solvent gradually dissolves
into the water. The difference between the two reference pressures can be readily
determined by applying the Laplace equation of capillarity at the origin O of the r-z
coordinate system:
2
Po 0 Ps 0 = (4)
R0
Where R0 is the principal radius of curvature of the profile at the apex of the pendant drop.
Recalling the differential geometry, the principal radii of curvature of the interfacial profile
d 2z
1 dr 2
= 3
(5a)
R1
dz
2 2
1 +
dr
dz
1 dr
= 1
(5b)
R2 dz 2 2
r 1 +
dr
d 2z dz
2 dr dr
gz = 3
+ 1
(6)
R0 dz 2 2 dz 2 2
1 + r 1 +
dr dr
yields:
d 2Z dZ
2 dr 2 dR
Bo Z = 3/ 2
1/ 2
(7)
R0 dZ 2 dZ 2
1 + R 1 +
dR dR
Where R0 = R0 / rn is the dimensionless principal radius of curvature of the profile at the
apex of the pendant drop; Bo = grn2 / is defined as the Bond number [10], which is the
ratio of the gravity to the capillary force. It is seen from the above equation that the shape
Equation (7) is a second-order ordinary differential equation and thus requires two
boundary conditions. These two boundary conditions can be easily found at the apex the
pendant drop:
Z R =0
=0 (8a)
dZ
=0 (8b)
dR R =0
The dimensionless principal radius R0 of curvature of the profile at the apex of the pendant
drop in the Eq. (7) is unknown and needs to be determined from the following volume
constraint:
z Vm
0
R 2 dZ =
rn3
= Vm (8c)
Here, Z d is the dimensionless height of the pendant drop; Vm and Vm are the measured
Equation (7) with its boundary conditions expressed by Eqs. (8a-8b) can be further
transformed into two first-order ordinary differential equations with their corresponding
initial conditions by using the so-called variable substitution method. Furthermore, these
two first-order ordinary differential equations can be solved numerically by using the
The numerical procedure for solving Eq. (7) subject to Eqs. (8a-8c) is briefly summarized
as follows.
where R = 0 and Z = 0 .
(2) Equation (7) together with its boundary conditions expressed by Eqs. (8a) and
Z
volume of the pendant drop is calculated by 0
R 2 dZ .
Z
(4) The calculated
0
R 2 dZ is then compared with the measured dimensionless
(5) If the dimensionless volume constraint expressed by Eq. (8c) is satisfied within
interfacial profile of the pendant drop can be readily obtained by inversely transforming
These results are accordance with previous studies [1,2]. The water-CO2 superficial
tensions () versus time have been measured for different pressures (3.0-5.0 MPa) and
pressures are shown in Figure 3, 4 and 5. If we compare these results with those obtained
with the previous studies, we can observe that they are very next when the data obtained for
a same temperature are compared, of same form that the previous studies it is possible to
identify two processes. First, at time 0 of the kinetics (0) decreased with CO2 pressure;
second, between pure water and CO2 decreased with time for the different pressures. It is
good to highlight, that every experiment was realized three times to analyze the deviation
associated to the method and the experimental equipment, and to every experimental data
the obtained results were very next, presenting a small standard minor diversion of 2%, for
pressure variable.
55
55
50
50
45
45
40
(mN/m)
40
(mN/m)
35 35
30 30
25 25 3.0 MPa 4.0 MPa 5.0 MPa
3.0 MPa 4.0 MPa 5.0 MPa
20 20
15 15
0:00 2:24 4:48 7:12 9:36 0:00 2:24 4:48 7:12 9:36
t (h) t (h)
Fig. 4 Superficial Tension CO2+Water to
Fig. 3 Superficial Tension CO2+Water
30 C
to 20 C
During the development of the experiment we observe the important in the precision and
exactitude of supporting the initial volume for every drop in the same conditions of
pressure and temperature. In the Fig.3, 4 and 5 we can observe that decreased more
intensively and rapidly when CO2 pressure increased until to arrive equilibrium value (eq).
Comparison of values obtained of the superficial tension in equilibrium value (eq), for a
same CO2 pressure with temperature variable (Fig. 6, 7 and 8), indicated that in 3.0 MPa,
the value eq is major for the temperature of 20 C followed by the of 30 C and finally the
of 40 C (Fig. 6); and to the 4.0 MPa eq are very next to the three temperature (Fig. 7) and
finally to the pressure of 5.0 MPa the value eq is major for the temperature of 40 C
followed by the of 30 C and finally the of 20 C (Fig. 8). If we compare these results with
those obtained with the previous studies, we can observe that they are very next in values
55 50
50
45
45
40
(mN/m)
40
(mN/m)
35
35
30
30
25 20 C 30 C 40 C
25 3.0 MPa 4.0 MPa 5.0 MPa
20 20
15 15
0:00 2:24 4:48 7:12 9:36 2:18 4:42 7:06 9:30 11:54
t (h)
t (h)
Fig. 5 Superficial Tension CO2+Water to Fig. 6 Superficial Tension CO2+Water to
40 C 3.0 MPa
Comparing the fig 6 with the Fig 8, we can observe that the changes of the superficial
tension are more rapid in the Fig. 8, therefore arrive to the equilibrium in minor time that
the Fig. 6, in consequence the mass transference is more quick to high pressures specially
55 55
50 50
45 45
40 40
(mN/m)
(mN/m)
35 35
30 30
25 25
20 C 30 C 40 C 20 C 30 C 40 C
20 20
15 15
0:00 2:24 4:48 7:12 9:36 12:00 0:00 2:24 4:48 7:12 9:36 12:00
t (h) t (h)
Fig. 7 Superficial Tension CO2+Water to
4.0 MPa Fig. 8 Superficial Tension CO2+Water to
5.0 MPa
With these results we can see the precision and exactitude of the experimental System, only
manipulating the pressure of the cylinder compressed air, without need to acquire
equipments of high costs for pressurize the phases. The disadvantage is it of not being able
Conclusion
In the system CO2 + Water the changes superficial tension at 20, 30 and 40 C to
different Pressures (Fig.3, 4 and 5) in a time 0 hour the initial superficial tension (0)
decreased with CO2 pressure, and superficial tension () water and CO2 system decreased
with the time for the different pressures. The mass transfer is quicker at 40 C for this
In the CO2-Water system the change of the superficial tension in equilibrium value (eq),
for a same CO2 pressure with temperature variable, at the low-pressure the decrease to 40
40 C than to 20 and 30 C, this allows to say that the superficial tension in equilibrium
References
[1] Frdric Tewes and Frank Boury, Thermodynamic and Dynamic Interfacial
2405-2412.
Surfactant near the Hydrate Formation Conditions. J. Chem. Eng. Data 2004, 49,
1023-1025.
[4] A.S. Michaels and E.A. Hauser, Interfacial Tension at Elevated Pressure and
[7] J. Fuentes and A.L. Lpez de Ramos, A New Method for the Determination of
Suface Tension Using the Pendant Drop Technique. Technologic information 2001,
Vol. 2, N 2.
[8] R. Span and W. Wagner, A New Equation of State for Carbon Dioxide Covering the
[9] W. Wagner and A. Prub, The IAPWS Formulation 1995 for the Thermodynamic
Properties of Ordinary Water substance for General and Scientific Use. J. Phys.
[10] Adamson, A.W. and Gast, A.P., Physical Chemistry of Surfaces, 6th Ed., John
1997, 112,135-149.
[12] Gu, Y. and Li D., Equilibrium Shape of Axisymmetric Liquid Solder Menisci around
a Vertical Cylinder on a Horizontal Substrate. J. Chem. Eng. Jpn. 1997, 30, 302-
309.
VII Iberoamerican Conference on Phase Equilibria and Fluid roperties for Process Design
. October 21st 25th , 2006 Morelia, Mich. (Mxico)
Transport Phenomenal Laboratory, University Simon Bolivar, Post Office Box 89.000
Caracas 1080-A Venezuela.
Corresponding autor. Tel.: +58 212 9064110, Fax: +58 212 9064137, E-mail:
[email protected]
Superficial tension for CO2 + Water system is reported for temperature from 293 to 323 K
and pressure range of (0-2) MPa. We have measured the change of superficial tension ()
between H2O and CO2 with a pedant drop method, the experiment was realized three times
to analyze the deviation associated to the method and the experimental system, and to every
experimental data the obtained results were very next, presenting a small standard minor
diversion of 2%. We observed that () decreased with time and with the CO2 pressure. But
when increase temperature to a same high pressure the superficial tension in equilibrium
value (eq) increase.
966
VII Iberoamerican Conference on Phase Equilibria and Fluid roperties for Process Design
. October 21st 25th , 2006 Morelia, Mich. (Mxico)
1. Introduction
Superficial tension plays an important role in the modeling and design of extraction,
cleaning, or emulsification processes. For instance, in the processes of extraction and
cleanliness, the main factors to determine the efficiency of the mass transference are the
number of interactions between the molecules in the interface and the drop size, which is
strongly dependent on superficial tension, between two phases reflects that their forces
cohesion and interaction [1]. These forces meet affected by the pressure and temperature of
the system, diminishing this way the value of the above mentioned tension. High-pressure
interfacial and surface-tension phenomena govern the migration and recovery of oil and gas
from hydrocarbon reservoirs. The phenomena are of particular relevance to phase
separation and mass transfer in light hydrocarbon fractionation plants [2]. Carbon dioxide
injection is a technique widely used in enhanced oil recovery processes, and the reduction
of superficial tension plays an important role in improving the mobility of reservoir fluids
in the reservoir [3]. The objective of this work is to measure the change superficial tension
of CO2 + Water system with the time for different temperature T= 293, 303, 313 and 323 K
to high-pressure, in a experimental unit of low cost in the system pressure, registered the
variation the drop volume in the time with constant pressure and temperature. And to
compare the results with Tewes and Boury work [1].
2. Experimental Section
Materials. The vessel was filled with carbon dioxide (purity 99.9% Vol.) Which were
supplied by AGA GAS, C.A. Distilled water was distilled prior to using was produced a by
our Laboratory.
Apparatus. The apparatus shown in figure 1 is similar to that used in the reference [3]
and the designs cell similar to the reference [4], except the pressure system has been
largely modified. The modified pressure system consists of a cylinder with a piston
designed under the norms ASMT [5] to high-pressure, where this contained phase disperses
pressurize with compressed air to calculate superficial tension inside of the high-pressure
cell [6]. The optical system consists in a high-resolution Sony photographic camera
connected to a computer is capable of processing the photographic data by using a program
Drop 98 [7] developed by our laboratory. This allows convenient observation and recording
of the drop forming process and the effect of drop-age on interfacial tension.
The apparatus consisted in a high-pressure cell, a light source, a CCD camera, a computer.
Pendant drops were formed at the end of a stainless steel tube 1,3mm inside diameter,
connected to a sample cylinder. The pressure inside the cell was controlled during all
experiment by regulatory valve of pressure. Therefore, it was possible to maintain constant
the pressure in the drop during the time of experience, so that the surface tension variation
was for the effect of the pressure in the whole interface of the drop. This experimental
scheme allows to realize investigations in superficial or interfacial tension to a low cost of
system pressure equipments.
967
VII Iberoamerican Conference on Phase Equilibria and Fluid roperties for Process Design
. October 21st 25th , 2006 Morelia, Mich. (Mxico)
968
VII Iberoamerican Conference on Phase Equilibria and Fluid roperties for Process Design
. October 21st 25th , 2006 Morelia, Mich. (Mxico)
Figure 2. Profile of drop analyzed with Laplace equation in the cylindrical coordinate
system.
Calculation of superficial tension. The interfacial tension was determined the profile of the
droplet (figure 2) using a CCD camera coupled to a video image profile digitizer board
connect to a computer. The drop profile was processed according to the fundamental
Laplace equation of capillarity (Eq. 1)
As illustrated in figure 2 , a pendent drop of water is formed inside a high-pressure cell that
is filled with the test solvent. The syringe needle has a radius of rn . The profile of the
pendent drop is governed by the famous Laplace equation:
1 1
Po Ps = + (1)
R1 R2
Where Po is the pressure in the pendent water drop; Ps is the pressure in the solvent phase;
is the Superficial tension between the water and the solvent; R1 and R2 are the principal
radii of the solvent-water interface.
In the absence of external forces other the gravity, the pressure is a linear function of the
elevation. For an arbitrarily chosen point P ( r , z ) on the interface, the pressure is Eq. (1)
can be expressed as:
z
Po = Po 0 g o dz (2a)
0
and
Ps = Ps 0 s gz (2b)
Where o is the density of the water solvent mixture; s is the density of the solvent [8] at
the experimental pressure and the temperature; g is the gravitational acceleration; Po 0 and
Ps 0 are the respective reference pressures of water and solvent at the origin O of the r-z
cylindrical coordinate system. Substituting Eqs. (2a 2b) into Eq. (1) and rearranging it
gives:
1 1
( Po 0 Ps 0 ) gz = + (3)
R1 R2
969
VII Iberoamerican Conference on Phase Equilibria and Fluid roperties for Process Design
. October 21st 25th , 2006 Morelia, Mich. (Mxico)
z
Where, = o dz / z s oave s is the density difference between the Water [9] and
0
solvent; oave is the average density of the water-solvent mixture inside the pendent drop.
During the measurement, oave is a function of time as the solvent gradually dissolves into
the water. The difference between the two reference pressures can be readily determined by
applying the Laplace equation of capillarity at the origin O of the r-z coordinate system:
2
Po 0 Ps 0 = (4)
R0
Where R0 is the principal radius of curvature of the profile at the apex of the pendant drop.
Recalling the differential geometry, the principal radii of curvature of the interfacial profile
at any point P(r,z) can be expressed as:
d 2z
1 dr 2
= 3 (5a)
R1
dz
2 2
1 +
dr
dz
1 dr
=
R2 1 (5b)
dz 2 2
r 1 +
dr
970
VII Iberoamerican Conference on Phase Equilibria and Fluid roperties for Process Design
. October 21st 25th , 2006 Morelia, Mich. (Mxico)
Where R0 = R0 / rn is the dimensionless principal radius of curvature of the profile at the
apex of the pendant drop; Bo = grn2 / is defined as the Bond number [10], which is the
ratio of the gravity to the capillary force. It is seen from the above equation that the shape
of the interfacial profile is dependent on the Bond number Bo and to-be-determined
dimensionless principal radius R0 .
Equation (7) is a second-order ordinary differential equation and thus requires two
boundary conditions. These two boundary conditions can be easily found at the apex the
pendant drop:
Z R =0 = 0 (8a)
dZ
=0 (8b)
dR R =0
The dimensionless principal radius R0 of curvature of the profile at the apex of the pendant
drop in the Eq. (7) is unknown and needs to be determined from the following volume
constraint:
V
0 R dZ = rnm3 = Vm
z
2
(8c)
Here, Z d is the dimensionless height of the pendant drop; Vm and Vm are the measured
dimensional and dimensionless volumes of the pendant drop, respectively.
Equation (7) with its boundary conditions expressed by Eqs. (8a-8b) can be further
transformed into two first-order ordinary differential equations with their corresponding
initial conditions by using the so-called variable substitution method. Furthermore, these
two first-order ordinary differential equations can be solved numerically by using the
fourth-order Runge-Kutta-fehberg method. The variable substitution approach and
complete mathematical transformation procedure for a similar problem can be found
elsewhere [11, 12]:
The numerical procedure for solving Eq. (7) subject to Eqs. (8a-8c) is briefly summarized
as follows.
(1) For a given Bond number Bo , a guessing value R0 ( )
is assumed for the
guess
dimensionless principal radius R0 of curvatura at the apex of the pendant drop, where
R = 0 and Z = 0 .
(2) Equation (7) together with its boundary conditions expressed by Eqs. (8a) and (8b) is
numerically solved using the fourth-order Runge-Kutta-fehberg method.
(3) After a complete interfacial profile Z = Z ( R) is determined, the dimensionless
Z
volume of the pendant drop is calculated by
0
R 2 dZ .
971
VII Iberoamerican Conference on Phase Equilibria and Fluid roperties for Process Design
. October 21st 25th , 2006 Morelia, Mich. (Mxico)
Z
(4) The calculated
0
R 2 dZ is then compared with the measured dimensionless volume
Vm of the pendant drop, as shown in Eq. (8c)
(5) If the dimensionless volume constraint expressed by Eq. (8c) is satisfied within a
prescribed accuracy, an accurate interfacial profile is obtained and the numerical
simulation is terminated. Otherwise, a new guessing value R0 ( )is generated
guess
972
VII Iberoamerican Conference on Phase Equilibria and Fluid roperties for Process Design
. October 21st 25th , 2006 Morelia, Mich. (Mxico)
75
70
65
60
/ (mN/m)
55
50
45
40
35
30
00:00 01:12 02:24 03:36 04:48 06:00 07:12 08:24
t /h
75
70
65
60
/ (mN/m)
55
50
45
40
35
30
00:00 01:12 02:24 03:36 04:48 06:00 07:12 08:24
t /h
the figures 3 and 4 present a rectilinear behavior with little variation in the tension
superficial initial, while that the figures 5 and 6 they show a differ for effect of the
temperature and CO2 pressure; especially the pressure of 0.0 MPa describes an abrupt
decrease, for the evaporation effect, because the transfer of the water molecules increases in
the interface to the gaseous phase (CO2).
973
VII Iberoamerican Conference on Phase Equilibria and Fluid roperties for Process Design
. October 21st 25th , 2006 Morelia, Mich. (Mxico)
75
70
65
/ (mN/m) 60
55
50
45
40
35
30
00:00 01:12 02:24 03:36 04:48 06:00 07:12 08:24
t /h
75
70
65
60
/ (mN/m)
55
50
45
40
35
30
00:00 01:12 02:24 03:36 04:48 06:00 07:12 08:24
t /h
During the development of the experiment we observe the important in the precision and
exactitude of supporting the initial volume for every drop in the same conditions of
pressure and temperature. Comparison of values obtained of the superficial tension in
equilibrium value (eq), for a same CO2 pressure with temperature variable (figures 7, 8 and
974
VII Iberoamerican Conference on Phase Equilibria and Fluid roperties for Process Design
. October 21st 25th , 2006 Morelia, Mich. (Mxico)
9), indicated that in 0.0 MPa, the value eq going to decrease with temperature like show
the figure 7 in the intervals time; the difference between eq is major to 1.0 MPa that to 2.0
MPa
70
65
60
55
/ (mN/m)
50
45
40
35
30
00:00 01:12 02:24 03:36 04:48 06:00 07:12 08:24
t /h
Figure 7. Variation superficial tension of CO2-Water to 0.0 MPa with the time to several
temperatures. , 20 C; , 30 C; , 40 C; , 50 C.
70
65
60
55
/ (mN/m)
50
45
40
35
30
00:00 01:12 02:24 03:36 04:48 06:00 07:12 08:24
t /h
Figure 8. Variation superficial tension of CO2-Water to 1.0 MPa with the time to several
temperatures. , 20 C; , 30 C; , 40 C; , 50 C.
975
VII Iberoamerican Conference on Phase Equilibria and Fluid roperties for Process Design
. October 21st 25th , 2006 Morelia, Mich. (Mxico)
70
65
60
55
/ (mN/m)
50
45
40
35
30
00:00 01:12 02:24 03:36 04:48 06:00 07:12 08:24
t /h
Figure 9. Variation superficial tension of CO2-Water to 2.0 MPa with the time to several
temperatures. , 20 C; , 30 C; , 40 C; , 50 C.
With these results we can see the precision and exactitude of the experimental System, only
manipulating the pressure of the cylinder compressed air, without need to acquire
equipments of high costs for pressurize the phases. The disadvantage is it of not being able
to work to pressures bigger than 7.0 MPa.
4. Conclusion
In the system CO2 + Water the changes superficial tension at 20, 30, 40 and 50 C to
different Pressures (figures 3, 4, 5 and 6) in a time 0 hour the initial superficial tension (0)
decreased with CO2 pressure, and superficial tension () water and CO2 system decreased
with the time for the different pressures. The mass transfer is quicker at 50 C for this
reason the system arrive to equilibrium in minor time.
In the CO2-Water system the change of the superficial tension in equilibrium value
(eq), for a same CO2 pressure with temperature variable, at the low-pressure the decrease
to 50 C is major that to 20, 30 and 40 C, but to high-pressure the decrease of the eq is
minor to 50 C than to 20, 30 and 40 C, this allows to say that the superficial tension in
equilibrium value decrease when decrease the temperature.
References
[1] F. Tewes, F. Boury, J. Phys. Chem. B 108 (2004) 2405.
[2] B. Chun, G. Wilkinson, Ind. Eng. Chem. Res. 34 (1995) 4371.
[3] C. Yu Sun, G. Chen, L. Yang, J. Chem. Eng. Data 49 (2004) 1023.
[4] A.S. Michaels, E.A. Hauser, J. Phys. & Colloid. Chem. 52 (1948)1157.
[5] American Society for Testing and Materials (ASTM), Standard Specification for
Pressure Vessel Plates, A 736/A 736M 88.
[6] A. Lpez de Ramos, R. A. Redner, R. L. Cerro, Langmuir 9 (1993) 3691.
976
VII Iberoamerican Conference on Phase Equilibria and Fluid roperties for Process Design
. October 21st 25th , 2006 Morelia, Mich. (Mxico)
[7] J. Fuentes, A.L. Lpez de Ramos, Technologic information, 12 (2001) ISSN 0716.
[8] R. Span, W. Wagner, J. Phys. Chem. Ref. Data 25 (1996) 1509.
[9] W. Wagner and A. Prub, J. Phys. Chem. Ref. Data 31 (2002) 387.
[10] A.W. Adamson and A.P Gast, Physical Chemistry of Surfaces, 6th Ed., John Wiley &
Sons Inc., New York 1997.
[11] Y. Gu, D. Li, P. Cheng, Colloids Surf. 112 (1997) 135.
[12] Y. Gu, D. Li, J. Chem. Eng. Jpn. 30 (1997) 302.
977
Apndice B:
Alejandro Jimnez
Juan Melndez
Susana Zeppieri
Aura Lpez de Ramos
[email protected]
[email protected]
[email protected]
[email protected]
Dpto. de Termodinmica y Fenmenos de Transferencia. Universidad Simn Bolvar
Valle de Sartenejas. Caracas-VENEZUELA.
Pedro Vargas
[email protected]
Dpto. de Energtica. Universidad Francisco de Miranda
Complejo Acadmico El Sabino, Av. Tchira, Punto Fijo-VENEZUELA
1. INTRODUCCIN
2. MODELO MATEMTICO
Las suposiciones del modelo de transferencia de masa para el interior de la gota son: la
transferencia de masa ocurre solo por difusin, la variacin del volumen de la gota no afecta
el perfil de concentraciones en el interior de la gota, la superficie externa de la gota se
encuentra saturada por el gas y la concentracin del gas en el interior de la gota en el inicio
del experimento es cero.
Figura 1- Contorno de la gota Figura 2- Esfera de rea de Figura 3- Mallado de la
colgante. transferencia de masa equivalente. gota.
c 1 c 1 c 1 2 c
= DAB 2 r 2 + 2 + 2 2 (1)
t r r r r tan( ) r
c(r , , t ) t =0 = 0 0 r rg 0 Iniciales
c(r , , t ) r = r = c sat 1 t>0 Borde
g
c(r , , t )
=0 0<r<rg t>0 Borde
=
c(r , , t )
=0 0 1 t>0 Borde
r r = rg
c r z t
C= R= Z= = 2
c sat rn rn rn
D AB
csat: Concentracin de saturacin
Donde
1 1 1 2 1 1
a= 2 b= + c= 2
R R R (R)2 R R R
1 1 1 2 1 1
d= 2 e= 2 f= 2 +
R 2 R tan( )
2 2
R 2 R 2 R tan( )
2 2
Iteracin implcita en
Donde
1 1 1 2 1 1
a= 2 + b= + 2 c=
R 2 R tan( )
2 2
R 2 R 2 R tan( )
2 2 2
1 1 1 2 1 1
d= + 2 e= f= + 2
R R R (R)2 R R R
2
1 m Vm (t ) Vc (t )
t
t m 0 Vm (t )
E= dt 100% (6)
3. Resultados y discusin
Como observa en la tabla 1, para los valores de los ngulos que definen la pared interna
y externa del capilar, calculados a partir del volumen inicial experimental de la gota colgante
asumindola como una esfera perfecta (Fig. 2), se obtuvo un porcentaje de error absoluto
menor a 2% para el coeficiente de difusin y 1% para el factor de hinchamiento. Los ngulos
se determinaron por simples relaciones y clculos trigonomtricos, empleando como datos el
radios externo y el espesor del capilar as como el radio de la esfera equivalente, calculado
con el volumen inicial experimental de la gota colgante de crudo (t=0 s con un V0=5,764
mm3).
El modelo numrico planteado en este trabajo tambin fue evaluado para el sistema
Agua-CO2 a las condiciones de temperatura de 25 C y de presin de 6,6 MPa con los datos
experimentales obtenidos por Hebach [7,8] de donde se obtuvieron los valores de cambio de
volumen de la gota de agua con respecto al tiempo, dimetro interno (0,25 mm), externo (1,60
mm) del capilar y volumen inicial de la gota colgante de agua (0,04015 cm3).
El resultado para el coeficiente de difusin estimado mediante el modelo simplificado de
este trabajo es mostrado en la tabla 2. Los resultados obtenidos mediante el modelo se
compararon con un grupo de correlaciones empricas [9,10]. Los resultados muestran que
existe una desviacin del 7% respecto al valor determinado por la ecuacin emprica de
Stokes-Einstein [9], mientras que las correlaciones empricas de Sutherland y Wilke-Chang
[9] presentan una desviacin mayor de 39 %, esto ocurre debido a que los factores de
correccin que se emplean en ambas ecuaciones no se ajustan a la condicin de alta presin.
Tabla 2. Resultados del la difusividad para el sistema CO2-agua (P=6,6 MPa y T=25 C).
NOMENCLATURA
REFERENCIAS
La zona de la pelcula est constituida por una pelcula de lquido y una burbuja
alongada de gas denominada burbuja de Taylor. La velocidad del lquido en la pelcula se
designa como vf y la del gas como vG. Aunque la regin lquida del slug se considera
axialmente homognea, las velocidades del lquido y del gas en la zona de la pelcula
cambian a lo largo de la tubera debido a la variacin del grosor de la misma, (z), detrs
del tapn de lquido.
Los anlisis hidrodinmicos del flujo slug se han basado generalmente en el modelo
unidimensional desarrollado originalmente por Dukler y Hubbard (1975). En este modelo
las causas de prdida de presin fueron restringidas al cuerpo del tapn y las prdidas en
la pelcula de lquido se consideraron despreciables. Taitel y Barnea (1990) demostraron
que la cada de presin para este patrn, puede ser calculada a partir de un balance de
fuerzas global en el que los trminos de diferencia de presin por friccin e hidrosttica
son considerados tanto en la regin de la burbuja como en la del tapn. Sin embargo, los
modelos antes mencionados tienden a subestimar los valores experimentales de cada de
presin en aproximadamente un 30%. Cook y Benhia (2000) presentaron un modelo
donde consideran cuatro fuentes de cada de presin: la prdida hidrosttica, la prdida
por friccin causada por los esfuerzos cortantes en la pared, las prdidas viscosas
asociadas con el movimiento en la regin de mezcla, y la prdida por aceleracin debido
a la expansin del gas.
Una gran cantidad de modelos para flujo slug ascendente estn disponibles en la
literatura. Estos pueden ser de dos tipos: El primer tipo consiste en ecuaciones empricas
para la prediccin del patrn de flujo, holdup y cada de presin las cuales han
contribuido al diseo de sistemas de flujo de dos fases por ms de 30 aos. Estos modelos
nunca han sido apropiadamente manejados para el fenmeno de flujo de dos fases, pero
han sido ampliamente utilizados en la industria [Hagedorn y Brown (1965), Brill y Beggs
(1991). La segunda categora est basada en el modelado de los mecanismos bsicos del
flujo. Estos pueden ser aplicados a condiciones de flujo diferentes a las que se utilizaron
para su desarrollo [Orell y Rembrand (1986), Nicklin et al. (1962), Hasan y Kabir
(1992)] sin embargo, requieren del conocimiento de variables como la longitud de la
pelcula de lquido, su forma, el perfil de velocidades a lo largo de la misma y,
especialmente, su altura y velocidad al final de esta zona.
La forma de la pelcula de lquido es muy compleja, particularmente la zona cercana a
la nariz del tapn. Es un problema tridimensional, con una superficie libre de flujo
turbulento. El modelo propuesto por Taitel y Barnea (1990), al igual que el de Dukler y
Hubbard (1975), emplea una aproximacin unidimensional llamada teora de flujo de
canal, en donde se desprecia la cada de presin en la pelcula. No obstante esta
aproximacin no es vlida cuando la burbuja de Taylor es muy larga generando una
desviacin del 30-40% respecto del valor experimental. Otra simplificacin comn, para
evitar una integracin numrica, es suponer que la pelcula de lquido tiene una altura
constante en equilibrio; sin embargo, esto slo debe hacerse para burbujas de gas
relativamente pequeas
En el trabajo de van Hout et al. (2003) se estudiaron las longitudes del slug de lquido
y los bolsillos de gas para diferentes posiciones a lo largo de la tubera y a diferentes
ngulos de inclinacin. Observaron que dichos valores varan ampliamente con los
ngulos de inclinacin, encontrando que para los slug de lquido, los mayores se observan
a 90 y decrecen cuando el ngulo se acerca a la horizontal. La tendencia de la longitud
de las burbujas de Taylor es inversa, ya que las burbujas mas largas se observan para el
caso de tubera horizontal y disminuyen a medida que aumenta la inclinacin. Estos
autores tambin concluyeron que las longitudes de los slug lquidos no dependen
significativamente de los caudales sino bsicamente del dimetro de la tubera; sin
embargo, la longitud de la burbuja alongada depende de los caudales de las fases y el
dimetro de tubera.
Xiao et al. (1990) y Gmez et al. (1999) estudiaron la frecuencia de aparicin del slug
a lo largo de tuberas de 0.024 y 0.054 m de dimetro interno para un sistema aire-agua
en el intervalo de inclinaciones de 2 a 90 de la horizontal. Encontraron que la frecuencia
alcanza un mximo cercano a 60 y que guarda una fuerte dependencia con la longitud y
velocidad de ascenso de las burbujas alongadas.
En el estudio del flujo slug, las tcnicas de visualizacin avanzada se han convertido
en una herramienta muy importante, ya que permiten obtener una secuencia de imgenes
de alta calidad, a partir de las cuales se pueden determinar con gran precisin el
comportamiento de los principales parmetros de este rgimen. Gopal y Jepson (1987) y
van Hout et al. (1998) utilizaron la visualizacin de imgenes para determinar el perfil de
velocidades del slug y la distribucin del gas, y a partir de esos datos lograron plantear
modelos para la prediccin de holdup de lquido.
En este trabajo se utiliza la tcnica de visualizacin de imgenes para evaluar el
comportamiento de parmetros hidrodinmicos del rgimen slug. Para diferentes caudales
de una mezcla aire-agua que ascienden por una tubera a diversos ngulos de inclinacin
(entre 0-90 medidos desde la horizontal), se mide el espesor de la pelcula a ambos
lados de la burbuja, la frecuencia de aparicin de las burbujas de Taylor y la longitud y
velocidad de ascenso de las mismas.
2. Metodologa Experimental
El sistema experimental empleado est formado por un banco de tres tuberas en forma
de U de 8.3 m de largo y de 0.75, 1.25 y 1.75 de dimetro interno, cada una, soportado
por una estructura metlica que puede rotar de 0 a 90. El lquido es suministrado por una
bomba de desplazamiento positivo y para medir el caudal se utilizan dos rotmetros en
paralelo de 0.04-0.30 L/s y 5.06x10-4-1.53x10-2 L/s. El aire es proporcionado por un
compresor y su caudal es medido con otro rotmetro 0.021-0.667 L/s. Ambos flujos
fueron controlados usando vlvulas de aguja. El mezclado de los fluidos se realiza en un
conector tipo T antes de ingresar a la tubera de prueba. La visualizacin se efectu con
una cmara de video de alta velocidad (llega a filmar hasta 4500 fps) que permite
almacenar las imgenes recolectadas para su procesamiento posterior. Con la finalidad de
evitar el efecto de la curvatura de la tubera en las imgenes, se film a travs de una
celda de visualizacin llena de glicerina (Figura 2).
10
8 =0
f (slug/s)
6 =90
4
2
0
0 50 100 150 200
Ns
Fig. 5. Frecuencia de aparicin de las burbujas de Taylor en funcin
del nmero de tapones contabilizados (Ns).
3. Resultados y Discusin
250
200 y = 1.1447x + 67.955
VTB (cm/s)
VTB (cm/s)
R2 = 0.968
150 200
y = 1.1223x + 68.949
2
R = 0.987
100 150
40 60 80 100 120 60 110 160
(a) VM (cm/s) (b) VM (cm/s)
300 300
VTB (cm/s)
y = 1.347x + 59.379
VTB (cm/s)
200 250 2
R = 0.8798
100 y = 1,1734x + 73,126
200
R2 = 0,9931
0 150
50 100 150 200 60 110 160
(c) VM (cm/s) (d) VM (cm/s)
300
200
150
80
60
LTB (cm)
40
20
0
200 220 240 260 280 300
VTB (cm/s)
El carcter irregular y poco estable del flujo slug se verific en las pruebas
experimentales realizadas, ya que las longitudes de las burbujas de Taylor que se
obtuvieron, para los mismos caudales, presentan una desviacin promedio de 55% para
cada ngulo estudiado (Figura 9). Sin embargo, los resultados obtenidos en todos los
casos para LTB/D coinciden con los reportados en la literatura, pues estn contenidos en el
intervalo de 4 a 30 (Shosho y Ryan, 2000).
35
30 VSG=15,82 cm/s VSL=56,08 cm/s
VSG=47,87 cm/s VSL=54,82 cm/s
25 VSG=15,92 cm/s VSL=48,5 cm/s
20
LTB/D
15
10
5
0
0 20 40 60 80
()
50
VSG=15,92 cm/s VSL=56,08 cm/s
40 VSG=91,32 cm/s VSL=48,5 cm/s
VSG=47,87 cm/s VSL=54,82 cm/s
30
Lf /D
20
10
0
0 20 40 () 60 80 100
90
80 VSG=15,92 cm/s VSL=56,08 cm/s
VSG=91,32 cm/s VSL=54,82 cm/s
70 VSG=47,87 cm/s VSL=54,82 cm/s
60
50
LU/D
40
30
20
10
0
0 20 40() 60 80
Fig. 11. Variacin de Lu/D en funcin del ngulo de inclinacin de la tubera para
diferentes condiciones de operacin.
7
VSG=15,92 cm/s VSL=56,08 cm/s
6 VSG=47,87 cm/s VSL=54,82 cm/s
5 VSG=47,87 cm/s VSL=28,29 cm/s
VSG=15,92 cm/s VSL=30,27 cm/s
f (slug/s)
Modelo k p j
Este trabajo 0.00045 2.10 1.91
Greskovich y Shrier, 1972 0.0226 1.2 2.02
Heywood y Richardson, 1979 0.0434 1.02 2.02
3,5
Greskovich y Shrier 20%
3,0
Heywood y Richardson
2,5
2,0
f cal
1,5 -20%
1,0
0,5
0,0
0,0 1,0 2,0 3,0
f exp
2,0
Espesor de pelcula superior
1,5 Espesor de pelcula inferior
pelicula /R
1,0
0,5
0,0
0 2 4 6 8
lB/D
Fig. 14. Altura de la pelcula de lquido alrededor de una burbuja de Taylor (=50).
90
70
45
30
Fig. 15. Evolucin de la forma de la burbuja de Taylor a medida que se inclina la tubera
para VSG=91.32cm/s y VSL=48.50 cm/s.
2,00
VSG=15,92 cm/s VSL=56,08 cm/s
1,50 VSG=47,87 cm/s VSL=54,82 cm/s
inf /R VSG=91,32 cm/s VSL=48,50 cm/s
1,00
0,50
0,00
0,00 4,00 8,00 12,00
LB/D
(a)
1,50
VSG=15,92 cm/s VSL=56,08 cm/s
VSG=47,87 cm/s VSL=54,82 cm/s
1,00 VSG=91,32 cm/s VSL=48,50 cm/s
inf /R
0,50
0,00
0,00 10,00 20,00 30,00 40,00
LB/D
(b)
Fig. 17. Espesor de la pelcula inferior a lo largo de la burbuja de Taylor
(a): =0. (b): =90.
(a) (b)
(c) (d)
Fig. 18. Comparacin de los valores experimentales del espesor de la pelcula inferior y
los estimados por el modelo de Cook y Behnia (1997). (a) 0; (b) 10; (c) 20; (d) 90.
Tal y como se observa en la Figura anterior, el modelo de Cook y Behnia predice valores
de espesor de pelcula con un 9% de desviacin promedio para inclinaciones de la tubera
inferiores a 20. Por encima de este ngulo las diferencias son mayores. Esto se debe a que
la velocidad de ascenso de la burbuja de Taylor se incrementa a medida que aumenta el
ngulo de inclinacin de la tubera (Taitel y Barnea, 1990) y el modelo es muy sensible a
este parmetro.
4. Conclusiones
5. Nomenclatura
6. Referencias
Akagawa K., Sakaguchi T., (1966). Fluctuation of Void Ratio in Two-phase. JSME, 9,
104-110.
Barnea D., Sller L. (1989) Void-Fraction Measurements in Vertical Slug Flow:
Applications to Slug Characteristics and Transition, International Journal of Multiphase
Flow, 495-504.
Brill J.P., H.D. Beggs. (1991). Two Phase Flow in Pipes. 6ta. Edicin.
Cook, M., M. Benhia. (2000). Pressure Drop Calculation And Modelling Of Inclined
Intermittent Gas-Liquid Flow. Chemical Engineering Science, vol. 55.
Dukler A., Hubbard M. (1975). A model for Gas-liquid slug Flow in horizontal and near
horizontal tubes, Ind. Engng Chem. Fundam, 337-347.
Fernandes R., Semial R., Dukler A. (1983). Hydrodynamic Model For Gas-Liquid Slug
Flow in Vertical Tubes, AICHE Journal, 981-989.
Gmez, L.E., Shoham, O., Taitel, Y. (1999). Prediction of Slug Liquid Holdup: Horizontal
to Upward Vertical Flow. International Journal of Multiphase Flow, vol 26, No. 3, 517.
Gopal, M., Jepson, W. P. (1987). Development Of Digital Image Analysis Techniques For
The Study Of Velocity And Void Profiles In Slug Flow. International Journal of
Multiphase Flow, vol 23, No. 5.
Hagedorn, A. R., Brown, K. E. (1965). Experimental Study of Pressure Gradients
Occurring During Continuous Two-Phase Flow in Small Diameter Vertical Conduits, J.
Pet Tech, 17, 475-48.
Hasan A., Kabir C. (1992). Two Phase Flow in Vertical and Inclined Annuli, International
Journal of Multiphase Flow, 279-293.
Moissis R., Grifith P., (1962) Entrance Effects in a Two-Phase Slug Flow, Journal of Heat
Transfer, 29-39.
Nicklin D., Wilkes D., Davidson J., (1962). Two Phase Flow in Vertical Tubes, Trans.
Instn Chem, 61-68.
Orell, A., Rembrand, R. (1986). A Model for Gas-Liquid Slug Flow in a Vertical Tube. Ind.
Eng. Chem. Fundam., 196-206.
Polonsky S., Barnea D., Shemer L. (1999). Averaged and Time-Dependent Characteristics
of the Motion of an Elongated Bubble in a Vertical Pipe, International Journal of
Multiphase Flow, 25, 795-812.
Shosho C., Ryan M. (2001). An Experimental Study of the Motion of Long Bubbles in
Inclined Tubes, International Journal of Multiphase Flow, 56, 2191-2204.
Taitel Y., Barnea, D. (1990). Two-Phase Slug Flow. Academic Press Inc., 83-132
Taitel, Y., Sarica, C., Brill, J.P. (1999). Slug Flow Modeling For Downward Inclined Pipe
Flow: Theoretical Considerations. International Journal of Multiphase Flow, vol. 26.
van Hout R., Barnea D., Shemer L., (1998). Evolution of Statistical Parameters of Gas-
Liquid Slug Flow along Vertical Pipes, International Journal of Multiphase Flow, 1579-
1602
van Hout R., Shemer L., Barnea D., (2003). Evolution of Hydrodynamics and Statistical
Parameters of Gas-Liquid Slug Flow along Inclined Pipes, Chemical Engineering
Science, 115-133
Xiao, J.J., Shoham, O., Brill, J. P. (1990). A Comprehensive Mechanistic Model For Two-
Phase Flow. SPE 20631.
Apndice D:
Grupo G10. Departamento de Termodinmica y Fenmenos de Transferencia, Universidad Simn Bolvar. Apartado
Postal 89000, Caracas, 1080 A, Venezuela.
e-mail: [email protected], [email protected]
RESUMEN
En este trabajo se determin experimentalmente la retencin de agua y se observaron los patrones de flujo en un
sistema bifsico lquido-lquido en tuberas verticales e inclinadas, los fluidos de trabajo fueron aceite (densidad
8672 kg/m3 a 21C y viscosidad de 0,13100,0001 Pas a 23C) y agua. Se trabaj con una tubera acrlica de
0,03175m de dimetro interno. Mediante una cmara de alta velocidad se observaron los patrones de flujo para una
tubera a 90 y 45. A partir de estas observaciones se elaboraron mapas de patrones de flujo y se compararon con los
reportados en la literatura. El holdup se midi mediante el mtodo de cierre simultneo de vlvulas. Para el caso de
90 como 45 se obtuvo el comportamiento esperado, el holdup de agua aumento al aumentar Usw manteniendo Uso
constante. Se evalu la correlacin obtenida por Hasan y Kabir (1988) para estimar holdup con buenos resultados
pues el error relativo promedio de los valores de la correlacin con respecto a los experimentales fueron de 6% para
flujo vertical y de un 2% para flujo inclinado, por lo que se recomienda para predecir holdup en flujo bifsico agua-
aceite a condiciones de operacin similares a las de este trabajo.
El flujo bifsico lquido-lquido est presente en una gran variedad de industrias: qumica, nuclear, de procesamiento,
petroqumica y en especial en la industria petrolera. Cuando el crudo es transportado desde el yacimiento hasta los
cabezales de produccin, y posteriormente a las estaciones de flujo y de descarga, ste va acompaado de otros
fluidos como agua, gas natural y sedimentos slidos, formando distintas configuraciones geomtricas en la tubera,
definidas como patrones de flujo. Existen diferentes tipos de flujo multifsico, para la industria petrolera los ms
relevantes son: gas-lquido, lquido-lquido, lquido-slido, gas-slido, gas-liquido-slido y gas-liquido-lquido.
El flujo gas-lquido es el que ha sido estudiado ms extensamente, a pesar de la gran importancia que tiene el
transporte lquido-lquido ya que, se ha estimado que el 80% de las fallas ocurridas en los sistemas de produccin y
transporte de petrleo, son causados por corrosin interna, elevando as los costos de reemplazo de tuberas. Los
hidrocarburos lquidos que se transportan mediante tuberas, frecuentemente contienen agua libre y la presencia de
este fluido puede provocar corrosin severa en las paredes de la tubera.
El fenmeno de corrosin interna depende de la configuracin geomtrica del flujo de los dos lquidos. Las
inspecciones ultrasnicas de oleoductos [1] revelan que la corrosin ms severa se presenta en tuberas en las que
una baja velocidad de flujo ocasiona segregacin de fases. Por lo tanto la determinacin del patrn de flujo es una
consideracin importante en la evaluacin del riesgo de la corrosin en oleoductos. Adicionalmente, muchas de las
variables de diseo como la cada de presin, el holdup y el coeficiente de transferencia de masa y calor dependen
del patrn de flujo existente.
En el caso del flujo lquido-lquido, an no se ha llegado a una unificacin de criterios en cuanto a los patrones de
flujo desarrollados y sus transiciones. La principal causa de esto es la similitud entre las propiedades de los dos
fluidos, existiendo un diferencial finito de densidad y viscosidad entre s, por lo que no puede ser despreciada una
variable frente a otra, como en el caso gas-lquido. Por tal motivo, no es evidente cual de las fases se mueve ms
rpido en un flujo estratificado, o bien cual es el fluido ms propenso a formar una gota o a coalescer, hacindose
necesario considerar factores adicionales.
Cada autor de acuerdo a su apreciacin ha propuesto patrones de flujo distintos. Oglesby 1979 report catorce
patrones de flujo diferentes, mientras que otros como Russell et al. y Malinowsky describen slo tres o cuatro [2].
Govier et al. [3] publicaron uno de los primeros trabajos en esta rea. Empleando un aceite de viscosidad 0,0201Pas,
en un intervalo de Usw y Uso de 0,003 m/s hasta 3,048 m/s, encontraron los patrones de flujo burbuja, tapn, espuma,
rizado (o semianular), pelcula (o anular), niebla.
Zavareh et al. [4] desarrollaron un trabajo experimental empleando una tubera de 12,8 m de largo, y dimetro de
1,842 m, para un intervalo de Uso entre 0,05 m/s y 0,5 m/s y Usw entre 0,02 m/s y 0,6 m/s, obteniendo 4 patrones de
flujo diferentes: burbuja, burbuja disperso, burbuja invertido y burbuja disperso invertido (Fig. 1).
Con la tubera vertical todos los regmenes de flujo pueden ser clasificados como flujo burbuja, en el cual las
burbujas de aceite relativamente grandes son arrastradas por el agua, que se presenta aqu como la fase continua. El
flujo burbuja disperso, consiste en burbujas muy finas o pequeas de aceite que estn dispersas en el agua, en tanto
que el flujo burbuja invertido, las gotas grandes de agua son arrastradas por un flujo de aceite, que en este patrn es
la fase continua; y finalmente el flujo burbuja disperso invertido, es aquel donde gotas muy pequeas de agua se
encuentran dispersas en el flujo de aceite.
Ms recientemente, Flores et al. [5] publicaron un trabajo muy completo en esta rea. Emplearon fluidos con una
relacin o/w = 0,85, o/w = 20 y tensin interfacial de 0,029 N/m, en una tubera de 15,5 m de longitud y 0,0508
m de dimetro. Los patrones de flujo encontrados estn clasificados en dos categoras principales: flujo segregado y
flujo disperso (Fig. 2).
El flujo segregado est caracterizado por la continuidad en direccin axial de las dos fases, mientras que en flujo
disperso la continuidad ha sido interrumpida en cualquiera de las dos fases. Los patrones de flujo identificados
fueron: i) dominados por agua: dispersin de aceite en agua (D o/w), dispersin muy fina de aceite en agua (VFD
o/w) y aceite en flujo agitado de agua (o/w CF), y ii) dominados por aceite: agua en flujo agitado de aceite (w/o CF),
dispersin de agua en aceite (D w/o) y dispersin muy fina de agua en aceite (VFD w/o).
Bannwart et al. [6] realizaron un trabajo experimental con la finalidad de estudiar los patrones de flujo en tuberas
empleando agua y crudo de alta viscosidad, 0,488 Pas en una tubera de 5,43 m de longitud y 0,0284 m de dimetro.
Los patrones de flujo los clasificaron por analoga con patrones gas-lquido, empleando una letra que denota el
patrn principal, una segunda letra que indica un subnivel y tambin una tercera si es necesaria una mejor
denominacin. En la Fig. 3 se incluyen las fotografas de algunos de ellos.
Figura 3. Patrones de flujo en tuberas verticales observados por Bannwart et al. [6].
En la literatura los patrones se reportan a travs de los llamados mapas de patrones de flujo, los cuales son la
representacin bidimensional en un sistema cartesiano de los regmenes de flujo. La gran mayora de los mapas son
experimentales y la transicin de un rgimen de flujo a otro lo determina el autor de manera subjetiva, de acuerdo a
la visualizacin que haya hecho de los mismos.
El comportamiento del flujo bifsico depende de las propiedades fsicas de cada fluido (viscosidad, densidad,
tensin interfacial) y de las condiciones de operacin (holdup, velocidades superficiales, dimetro e inclinacin de la
tubera), por tal motivo cada mapa de flujo ser vlido para las condiciones de operacin y los fluidos utilizados para
su generacin. De esta manera, si se desea hacer un estudio, ya sea de diseo o de evaluacin en un sistema de
produccin, habra que realizar experimentos con los fluidos y las condiciones de operacin deseadas. Por esta razn,
surgi la necesidad de predecir el comportamiento del flujo bifsico sin tener que realizar pruebas experimentales
para cada par de fluidos o condiciones de operacin. As, nacieron las aproximaciones empricas [7], las que si bien
son sencillas de usar presentan el inconveniente de que no explican el fenmeno fsico del flujo en dos fases. La otra
gran desventaja es que su uso est limitado a condiciones similares para las que fueron desarrolladas, fuera de ellas,
las desviaciones en los valores predichos aumentan significativamente.
Debido a las limitaciones mencionadas, surgi una nueva aproximacin por modelos mecansticos. En stos, se
pretende aclarar el fenmeno fsico modelando matemticamente el mecanismo que causa el flujo bifsico. El primer
objetivo de esta aproximacin es predecir el rgimen de flujo existente en un sistema dado, para luego desarrollar un
modelo separado para cada patrn de flujo con el fin de predecir su hidrodinmica y/o transferencia de calor.
El modelo mecanstico desarrollado por Trallero [7] y posteriormente modificado por Trallero et al. [8] es el ms
utilizado actualmente para predecir la transicin entre los patrones de flujo en tuberas horizontales o ligeramente
inclinadas. Este modelo se basa en un anlisis de estabilidad para las transiciones entre los patrones segregados, y
una comparacin entre las fuerzas gravitacionales o de flotacin con las fuerzas de arrastre, para las transiciones
entre los flujos dispersos. Para realizar estos clculos se necesita conocer variables como las velocidades
superficiales y el holdup de cada fase, el dimetro de las gotas, entre otras. Sin embargo, muchas veces no se dispone
de correlaciones precisas para su estimacin, por lo que se requieren sistemas que permitan determinarlas
experimentalmente.
Para el caso de flujo vertical, la informacin es ms escasa razn por la cual en este trabajo, se determinarn los
patrones de flujo y la retencin de fases para el sistema agua-aceite.
PARTE EXPERIMENTAL
El equipo utilizado consiste en una tubera de acrlico en forma de U de 8.3 m de largo y de 1.25 in de dimetro
interno, soportada por una estructura metlica que puede rotar de 0 a 90 (Fig. 4).
El sistema de alimentacin de los fluidos consta de 3 lneas principales, una para la fase gaseosa (aire) y otras 2
para las fases lquidas (una para fluidos acuosos y otra para aceites). Cada lnea de lquido cuenta con una bomba de
desplazamiento positivo con una potencia mxima de 5,5 kW, mientras que el aire es proporcionado por un
compresor capaz de mantener la presin en 100 psig. Para la medicin del flujo volumtrico de gas se utiliza un
rotmetro de 2 a 20 pie3/min y un medidor tipo VORTEX con un intervalo de medicin de 0 a 22,78 L/s. Para las
lneas de lquido se utiliza: un rotmetro de 0 a 10 gal/min para la lnea de lquidos acuosos, y un rotmetro de 1 a 10
gal/min para la lnea de los aceites. Ambas lneas de lquido cuentan con un medidor tipo VORTEX de 0 a 13,33 L/s.
Una vez que los fluidos pasan por la tubera de prueba, llegan a la seccin de separacin. Esta est conformada por
dos tanques de 500 L de capacidad donde se deja que la mezcla agua-aceite se separe por gravedad antes de retornar
el aceite al tanque de alimentacin respectivo.
Una vez seleccionados los caudales de ambas fases y el ngulo de inclinacin de la tubera, se inicia la prueba,
esperando unos dos minutos hasta alcanzar condiciones estacionarias. El patrn de flujo se observa directamente a
travs de las paredes de la tubera y se graba con la cmara de alta velocidad a travs de la celda visualizacin
ubicada a una distancia L/D de 200, medida desde el punto de mezcla de las fluidos, con el fin de almacenar la
informacin para anlisis posteriores. Para determinar el holdup de lquido se acciona el interruptor On-Off que
permite aislar la parte del flujo en donde se van a realizar las mediciones. Posteriormente se espera a que la mezcla
bifsica se separe en sus fases y se procede a medir la altura de lquido. Se repite el procedimiento para la
determinacin del holdup hasta que la desviacin entre una medicin y el promedio de las anteriores sea menor al
3%. Posteriormente se procede a fijar otra combinacin de caudales y se repite el mismo procedimiento.
En la Tabla 1 se muestran las propiedades fsicas de los fluidos usados en este trabajo.
(kg/m3) (Pas) (N/m)
Aceite 8672 0,1310
0,02970009
Agua 997,538 0,000941
RESULTADOS Y DISCUSION
Variando los caudales de ambas fases, se realiz un total 106 experimentos. Los patrones de flujo observados
fueron cuatro, de los cuales tres son dominados por la fase acuosa (Intermitente o tapn, dispersin de aceite en agua,
Do/w, y dispersin muy fina de aceite en agua, VFDo/w) y el restante posee dos fases continuas (anular).
La base de los patrones de flujo es bastante similar al caso gas-lquido en tuberas verticales. A bajos flujos de
aceite y moderados flujos de agua el patrn de flujo consiste en gotas de aceite en agua. Pero a altos flujos de agua,
se observa dispersin de aceite en agua. El flujo intermitente consiste de una serie de gotas de aceite en forma de
campana, muy similar al flujo gas-lquido. El flujo anular consiste en un cilindro de aceite rodeado por un anillo de
agua. Cada uno de estos patrones base, puede aparecer con una dispersin de pequeas gotas de aceite en agua,
causando patrones de flujo combinados [6]. En todos los casos, se observ que las transiciones de un patrn de flujo
a otro ocurren de forma continua y no abrupta.
En las Figuras 6 (a) y (b) se pueden apreciar los patrones de flujo dispersin de aceite en agua, tapn y anular. El
patrn de flujo dispersin de agua en aceite presenta burbujas de aceite grandes separadas entre s por una distancia
considerable, y no presenta dispersin muy fina de aceite entre las mismas. El flujo intermitente o tapn (Fig. 7 (a)),
consiste en una burbuja larga de aceite en forma de bala en su extremo superior, sin embargo hacia su parte
intermedia presenta una interfase ondulada hasta que llega la cola, donde culmina con una forma irregular. El espacio
comprendido entre un tapn y otro suele estar ocupado por burbujas de aceite pequeas. Es posible ver dos tapones
seguidos sin la presencia de burbujas de aceite entre uno y otro. En el anular se observa una fase continua de aceite
en el centro de la tubera cubierta por un nulo de agua, donde la interfase del aceite es ondulada por lo que el
dimetro de la columna de aceite no es constante.
En la Figura 6 (a), para una Usw=0,24m/s y Uso=0,02m/s se aprecia un patrn de flujo de burbuja dispersa (A). Si
se aumenta el caudal de aceite, estas gotas comienzan a coalescer adquiriendo la forma de campana caracterstica del
flujo intermitente (B, C), presentndose la transicin entre estos dos patrones. Al seguir aumentando la velocidad de
aceite se observan tapones de forma similar a los mostrados en la Fig. 8, pero que difieren en longitud y cantidad de
gotas asociadas a la cola del mismo.
En la Figura 6 (b) se observa la transicin entre el patrn de flujo tapn y el anular. Este ltimo muy similar al
encontrado por Bannwart et al. [6] y que denomin Anular Ondulado Distorsionado + Burbuja Dispersa (AODBD),
que se muestra en la Fig. 4.
En lneas generales, las transiciones que se observan en las Figs. 6 (a) y (b) son bastante claras, y la tendencia es
que al ir aumentando Uso, manteniendo Usw constante, las burbujas de aceite coalescen y forman tapones. Al
aumentar ms Uso los tapones van aumentando su tamao hasta que finalmente forman una fase continua de aceite
en el centro de la tubera, rodeado por la fase acuosa.
En la Fig. 7 (b) se muestran los patrones de flujo obtenidos manteniendo la velocidad superficial de aceite
constante Uso=0,02m/s y variando la de agua. Para Usw de 0,29m/s (A) y 0,32m/s (B) se aprecian burbujas de aceite
grandes, con forma de tapn, antecedidas y precedidas por burbujas de aceite de menor tamao. A medida que se
incrementa la velocidad superficial del agua pasando por 0,35m/s (C), 0,38m/s (D) y 0,41m/s (E), se observa una
mayor turbulencia, provocando que aumente la frecuencia de aparicin en los tapones, llegando en algunos casos a
chocar entre s, lo que conlleva a su colapso y separacin en gotas de aceite ms pequeas dando paso al patrn de
flujo de gotas dispersas, provocndose as la transicin intermitente-gotas dispersas.
(a) (b)
Fig. 6. (a) Transicin dispersin de aceite en agua a intermitente o tapn. Usw=0,24m/s. Uso: (A) 0,02m/s; (B)
0,05m/s; (C) 0,08m/s; (D) 0,13m/s y (E) 0,18m/s. (b) Transicin intermitente anular Usw=0,24m/s. Uso: (A) 0,23m/s;
(B) 0,29m/s; (C) 0,34m/s y (D) 0,39m/s.
(a) (b)
Fig. 7. (a) Patrn de flujo intermitente o tapn. Uso=0,02m/s. Usw : (A) 0,10m/s, (B) 0,13m/s, (C) 0,16m/s, (D)
0,21m/s y (E) 0,24m/s. (b) Transicin intermitente dispersin de aceite en agua. Uso=0,02m/s. Usw: (A) 0,29m/s, (B)
0,32m/s, (C) 0,35m/s, (D) 0,38m/s y (E) 0,41m/s.
Todos los valores experimentales de patrones de flujo, fueron graficados en funcin de las velocidades
superficiales de ambas fases obtenindose el mapa de patrones presentado en la Fig. 8(a). El intervalo de velocidades
superficiales con que se trabaj corresponde a 0,02 m/s < Uso < 0,39 m/s y 0,16 m/s < Usw < 0,72 m/s. Como se
aprecia en la Fig. 8 (b), hay zonas bien diferenciadas donde se ubican las regiones de patrones de flujo, en la parte
inferior el flujo tapn, en la regin media la dispersin de aceite en agua, en la regin superior la dispersin muy fina
de aceite en agua, y en la regin inferior derecha el flujo anular. El patrn de flujo ms observado fue el de
dispersin muy fina de aceite en agua, seguido por el flujo tapn y la dispersin de aceite en agua. El que menos se
present fue el patrn de flujo anular, por las condiciones de operacin del equipo.
En la Figura 9 (a) y (b) se muestra la retencin de agua promedio, obtenida para cada prueba experimental para
inclinaciones d tubera de 45 y 90, en funcin de los caudales de agua y aceite. Tal y como se observa, Hw aumenta
con el caudal de agua y disminuye con el caudal de aceite. En la misma Fig. se muestra en lneas continuas los
valores arrojados por la relacin de Hasan y Kabir [11]:
USO
Ho = 1 H w = (1)
g ( w o )
1/ 4
1, 2U M + 1,53 H 2
2w
w
1 1
PS
CC
VDFo/w
CC-PS
CT
Tapn
VFDo/w
Anular
Anular + VFDo/w
Do/w
Tapn-Do/w
0,1
0,1
0,01 0,1 1
0,0 0,1 1,0
Uso (m/s) Uso (m/s)
.
(a) (b)
Fig. 8. Mapa de patrones de flujo para (a) 45 y (b) 90 de inclinacin obtenido experimentalmente.
1 1
0,8 0,8
Hw Hw
0,6 0,6
CONCLUSIONES
Para el flujo simultneo de agua y aceite de densidad 867 kg/m3 y viscosidad 0,1310 Pas, en una tubera
vertical ascendente de 1,25 in de dimetro interno, se obtuvieron los patrones de flujo: dispersin de aceite en agua
pseudo patrn (PS), dispersin de aceite en agua cocorriente (CC), dispersin muy fina de aceite en agua (VFD o/w),
transicin entre flujo CC y PS y dispersin de aceite en agua contracorriente (CT) para 45, para el caso de 90 los
patrones de flujo observados fueron intermitente, dispersin de aceite en agua (D o/w), dispersin muy fina de aceite
en agua (VFD o/w) y anular. En todos los casos, se observ que las transiciones entre patrones son suaves y
continuas.
El holdup de agua obtenido experimentalmente aumenta a medida que aumenta Usw y disminuye Uso. La
correlacin de Hasan y Kabir [11] predijo los valores de retencin medidos con un error promedio de 6%.
NOMENCLATURA
H Retencin (-)
U Velocidad superficial (m/s)
Densidad (kg/m3)
Viscosidad (Pa.s)
Tensin interfacial (N/m)
Subndices
m Referido a la mezcla
so Superficial de aceite
sw Superficial de agua
o Referido a la fase oleosa
w Referido a la fase acuosa
AGRADECIMIENTOS
Los autores agradecen al FONACIT (Proyecto 97003788) por el financiamiento de los equipos y el sistema
experimental empleado.
REFERENCIAS
[1] Oronel C., Estudio Experimental del Flujo Bifsico Lquido-Lquido en Tuberas Horizontales y Ligeramente
Inclinadas, Anteproyecto de Maestra, Universidad Simn Bolvar, Caracas, 2005.
[2] Russell, T. W., Hodgson, G. W. y Govier, G. W. Horizontal Pipeline Flow Mixtures of Oil and Water, Can. J.
Chem. Eng. 37, 9-17, 1959.
[3] Govier, G. W., G. A Sullivan y R.K. Wood, The Upward Vertical Flow of Oil-Water Mixtures, The Canadian
Journal of Chemical Engineering, 39, 67-72, 1961.
[4] Zavareh F., Shoham O. y Brill J. P. A Comprehensive Mechanistic Model for Two-Phase Flor in Pipelines.
SPE 63 Annual Technical Conference and Exhibition, paper SPE 20631, New Orleans, 1990.
[5] Flores J. G., X. T. Chem y J. P. Brill, Characterization of Oil-Water Flow Pattern in Vertical and Deviated
Wells, en SPE Paper 38810, SPE Proceedings Annual Technical Conference and Exhibition, 1; 601-610, 1997.
[6] Bannwart A., Rodrguez O. y Carvalho C. Flow Patterns in Heavy Crude Oil-Water Flow, Journal of Energy
Resources Technology, Vol. 126, 184-189, 2004.
[7] Charles, M.E., G. W Govier. y G. W. Hodson, The horizontal pipeline flow of equal density oil-water mixture,
The Canadian Journal of Chemical Engineering, 39, 27-36,1961.
[8] Trallero J. L. Oil-Water Flow Patterns in Horizontal Pipes, PhD Thesis, The University of Tulsa, 1995.
[9] Trallero J. L., Sarica, C. y Brill, J. P. A Study of Oil-Water Flow Patterns in Horizontal Ppipes, SPE Annual
Technical Conference and Exhibition, paper SPE 36609, Denver, Colorado, Octubre 6-9, 1996.
[10] Romero C., Estudio de Patrones de Flujo Bifsico lquido-gas en tuberas inclinadas de 0 a 90, Trabajo de
Grado en la Maestra de Ingeniera Qumica (a ser presentada en marzo 2007), Universidad Simn Bolvar,
Caracas, 2007.
[11] Hasan, A. R y Kabir. C. S., A new Model for Two-Phase Oil-Watern Flow: Production Log Interpretation and
Tubular calculations, Proceeding Annual Technical Conference and Exhibition, 3, 369-376, 1988.
Apndice E:
1. INTRODUCCIN
Durante la exploracin y produccin de los yacimientos petrolferos es de vital
importancia predecir la facilidad de transporte de los fluidos presentes en el yacimiento. Para
ello, se utilizan las curvas de permeabilidad relativa, kri, en funcin de la saturacin, Si, de
alguna de las fases presentes en el sistema. Dada la importancia de estas curvas, se han
desarrollado diversas metodologas (Gutirrez, 2004) para su estimacin, para las que se
requiere la determinacin exacta de la fraccin volumtrica de cada una de las fases en el
medio poroso. Esta variable se puede inferir indirectamente a partir de la medicin de los
fluidos producidos del medio poroso o se puede determinar directamente en el ncleo de
pruebas (Honarpour y Mahmood, 1986).
La visualizacin de flujo es una tcnica interna para el clculo de la saturacin que
consiste en visualizar la distribucin de cada uno de los fluidos dentro del medio poroso. En
estudios recientes (Avraam y Payatakes, 1995; Lago et al., 2000; Heny et al., 2001), se ha
utilizado esta tcnica en medios porosos artificiales, transparentes y poco profundos
(micromodelos) para calcular las saturaciones de las fases y obtener informacin estadstica
(tamao y distribucin de burbujas) de la fase gaseosa en el medio poroso. Esta metodologa
requiere el uso de cmaras de alta resolucin, iluminacin adecuada y programas de
digitalizacin de imgenes. Una de las ventajas ms importantes de la tcnica de visualizacin
es que permite establecer de manera exacta, la presencia o no de gradientes de saturacin en el
medio poroso, los cuales son consecuencia de fenmenos como el efecto de borde (retencin
de la fase mojante en la salida del sistema) o la compresibilidad de la fase gaseosa en el caso
de flujos gas-lquido.
El presente estudio tiene como objetivos medir y registrar las variables experimentales
necesarias que permitan caracterizar un micromodelo de vidrio y utilizar la tcnica de
visualizacin con cmaras de alta velocidad, para cuantificar los gradientes de saturacin de
un flujo bifsico gas-lquido en un micromodelo de vidrio.
2. MARCO TERICO
Las curvas kri-Si pueden ser afectadas por diversos factores (Geffen et al., 1951; Torsaeter
y Abtahi, 2000). En el caso particular de las curvas obtenidas experimentalmente en
laboratorios, se debe considerar con especial importancia el efecto de borde y la
compresibilidad del gas.
El efecto de borde no es ms que la retencin de la fase mojante, producida por fuerzas
capilares, justo antes de la salida del medio poroso. Esta discontinuidad en las propiedades
capilares hacia el final del sistema siempre existir cuando por l circulen dos o ms fluidos
inmiscibles y generar un gradiente de saturacin de la fase mojante en la direccin del flujo
(Richardson et al., 1952). Entre los mecanismos que se utilizan para minimizar el efecto de
borde en pruebas de laboratorio, se encuentran (Osoba et al., 1951): colocar una muestra del
mismo material justo antes y despus del ncleo donde se realiza la toma de datos
experimentales, utilizar altos flujos volumtricos y retener la fase mojante en el ncleo a
travs de fuerzas capilares dejando fluir solamente la fase no mojante, entre otros.
La compresibilidad del gas reviste especial importancia en experimentos de laboratorios
donde la presin de operacin es cercana a la presin atmosfrica (Geffen et al., 1951). El
volumen poroso ocupado por el gas se incrementa progresivamente desde la entrada y hacia la
salida del sistema debido a la disminucin de la presin, esto genera un gradiente de
saturacin del gas en la direccin del flujo. Este efecto se puede minimizar haciendo que la
cada de presin en el ncleo sea despreciable frente a la presin de operacin a la entrada del
sistema.
Los efectos de gradientes de saturacin en el medio poroso se pueden puntualizar de la
siguiente manera (Osoba et al., 1951):
El clculo de una saturacin promedio implica el clculo de valores promedios de
permeabilidades relativas a ambas fases, esto se traduce en velocidades estimadas de
produccin muy diferentes a las reales.
En un sistema donde est presente un gradiente de saturacin siempre existir un
gradiente de presin capilar. Esto se traduce en que cada una de las fases fluyentes tenga
una cada de presin diferente por unidad de longitud. Considerar despreciable la
diferencia en las cadas de presin, puede resultar en valores errneos de
permeabilidades relativas a la fase mojante.
2.1 Clculo de la saturacin en el laboratorio
Dado que es indispensable conocer en forma exacta la Si que ocupa cada una de las fases
dentro del medio poroso para la estimacin de las curvas kri-Si, se han desarrollado tcnicas
externas y tcnicas internas para el clculo de la saturacin (Honarpour y Mahmood, 1986).
En las tcnicas externas la saturacin en el ncleo es inferida indirectamente a partir de
las mediciones de los fluidos producidos. Estas tcnicas permiten reportar un valor de
saturacin promedio, ms no un perfil de saturaciones en el ncleo. Cuando el volumen poroso
es pequeo, esta tcnica no resulta confiable, dada la presencia de volumen muerto y prdida
por evaporacin, entre otros. Las tcnicas externas requieren remover el ncleo de pruebas del
montaje experimental.
En las tcnicas internas o in situ, la cantidad de cada uno de los fluidos se determina sin
perturbar la distribucin de los fluidos en el ncleo de pruebas. Por tanto, estas tcnicas son
ms exactas y confiables que las tcnicas externas (Hornapour y Mahmood, 1986).
Una de las tcnicas in situ consiste en visualizar la distribucin de cada uno de los fluidos
dentro del medio poroso. En aos recientes (Avraam y Payatakes, 1995; Lago et al., 2000;
Heny et al., 2001), cada vez es ms comn el uso de esta tcnica en estudios experimentales
para calcular las saturaciones de las fases en medios porosos transparentes. Esta tcnica
tambin se ha utilizado con el fin de obtener informacin estadstica (saturacin, tamao y
distribucin de las burbujas) de la fase gaseosa en el medio poroso. Esta metodologa requiere
el uso de cmaras de alta velocidad, iluminacin adecuada y programas de digitalizacin de
imgenes. Adicionalmente, el medio poroso debe ser transparente y poco profundo, de forma
que permita la visualizacin dentro de los canales porosos.
A pesar del error que involucra la magnificacin de las imgenes de visualizacin de los
flujos en el sistema, esta tcnica es una de las ms eficientes y exactas al momento de calcular
fracciones volumtricas y/o estudiar patrones de flujo para sistemas bifsicos en medios
porosos transparentes (Ferrigni et al., 2003).
2.2 Experimentos en micromodelos de vidrio
Los micromodelos son medios porosos creados artificialmente con el objetivo de observar
el flujo de fluidos a escala de poro. Bsicamente, consisten en dos placas de vidrio, plexigls o
cualquier material resistente y transparente, donde una de las placas tiene grabado un patrn
geomtrico, ordenado o al azar, que simula los obstculos y canales de un medio poroso
natural. Una de las limitaciones del uso de micromodelos, es que slo simulan un flujo
bidimensional, lo que permite obtener slo una representacin cualitativa ms que cuantitativa
de lo que ocurre realmente en el yacimiento. Sin embargo, son una herramienta til en la
visualizacin de los patrones de flujo y la validacin de los modelos de nucleacin,
crecimiento, coalescencia y ruptura de burbujas.
Recientemente, en diversos estudios experimentales (Bora et al., 1997; Lago et al., 2000;
Heny et al., 2001; Kamp et al., 2001) se han utilizado micromodelos en representacin de
medios porosos, con el fin de estudiar y analizar el desarrollo de los fluidos en procesos de
empuje por gas en solucin.
3. DISEO Y METODOLOGA EXPERIMENTAL
El equipo experimental (Fig. 1) consta de una tubera de cobre de de dimetro externo.
El aire proviene de una bombona de aire comprimido (BB) que cuenta con un regulador de
presin (RP) de escala 0-400 psig y apreciacin 10 psig. La fase lquida forma un circuito
cerrado y se almacena en un tanque de plstico (TT) de 38,8 l (0,0388 m3) de capacidad. La
salida del tanque est conectada a una electrobomba con rodete perifrico (B) que permite
impulsar el lquido a la tubera. El flujo de lquido se regula a travs de un rotmetro que
permite medir un caudal mximo de agua de 564 ml/min (9,4.10-6 m3/s) en escala milimtrica
(0-150 mm) con apreciacin de 1mm.
El aire y el lquido se unen en un conector tipo T. El flujo bifsico pasa por un tubo de
vidrio en forma de L que permite el ingreso del flujo al micromodelo (MM) y al transductor
de presin (T), simultneamente. El transductor de presin (T) est conectado a un indicador
de presin (I) donde se registra la cada de presin que se produce en el micromodelo. Este
indicador tiene una apreciacin de hasta 0,001% de la cada de presin mxima.
El medio poroso utilizado en este estudio es un micromodelo de 250x100 mm2 que consta
de dos placas planas de vidrio. La placa inferior, de 3 mm de espesor, tiene grabada un patrn
de obstculos cuadrados de arista 0,6 mm aproximadamente y una separacin entre ellos de
0,6 mm. El micromodelo y el patrn de obstculos se pueden observar en la Fig. 2.
Figura 2. Micromodelo de vidrio utilizado en las pruebas experimentales (3, Fig. 1).
3.1 Sistema de visualizacin
El sistema de visualizacin (Fig. 3) consta de una cmara de alta velocidad marca Kodak
Ektapro modelo 4540mx Imager, la cual permite grabar desde 30 hasta 4500 cuadros por
segundo en la modalidad de pantalla completa. Se cuenta, adems, con un lente 60mm f/2.8D
AF Micro-Nikkor, marca Nikon, que permite aumentar o disminuir la abertura focal. A este
lente se le puede adaptar un anillo de extensin PK13 para la magnificacin de la imagen.
La cmara de alta velocidad est conectada a un monitor donde se observan las imgenes
del flujo en el micromodelo. Adicionalmente, la cmara est conectada a una estacin de
trabajo donde se graban las imgenes digitales para su posterior anlisis.
4,E+05
3,E+05
3
P (Pa) = 2E+12.Q (m /s)
2
R = 0,2631
2,E+05
1,E-07 2,E-07 2,E-07 3,E-07 3,E-07
3
Flujo volumtrico, QL (m /s)
FM . u FM
N Ca = (1)
FM .
.u k
Re = (2)
Estos valores, NCa y Re, cuantifican las fuerzas viscosas y las fuerzas capilares del fluido
presente en el micromodelo. En particular, el aumento en el valor del NCa se interpreta como
un aumento de las fuerzas viscosas en relacin con las fuerzas interfaciales. Igualmente, y
como consecuencia del aumento de las fuerzas viscosas, se aprecia una disminucin en el Re.
Sin embargo, y a pesar del aumento del NCa, no se obtiene un valor lo suficientemente alto
como para que prevalezcan, de forma significativa, las fuerzas viscosas sobre las fuerzas
interfaciales. Diversos estudios (Fulcher et al., 1983) han demostrado que para valores del NCa
en el orden de 1.10-3 se favorece el desplazamiento de la fase mojante en el medio poroso.
Tabla 2. Viscosidad, tensin superficial, nmero capilar y nmero de Reynolds para el agua
destilada y la solucin diluida de glicerina en el micromodelo de vidrio.
Sustancia (cP) (mN/m) NCa.104 Re
Agua 1 72,50 62 0,160,08
Sol. diluida de glicerina 7,5 62,47 4020 0,020,01
100
pG=30 psi
Sat. de lquido, SL (%) pG=50 psi
80
pG= 70 psi
pG=90 psi
60
40
0 50 100
65
pG=30 psi
45
Sat. de gas, SG (%)
pG=50 psi
pG=70 psi
pG=90 psi
25
5
0 50 100
24
Veloc. del gas, uG (cm/s)
pG=90 psi
16
pG=70 psi
pG=50 psi
pG=30 psi
0
0 50 100
Distancia desde el ingreso de los fluidos al
micromodelo (%)
1 INTRODUCCIN
2 MARCO TEORICO
La tensin en la interfase de una gota esfrica suspendida en otra fase inmiscible es dada
por la ecuacin de Laplace-Young:
2.
Pc = Pb Pa = (1)
R
2. 2.k .
Pc = Pb Pa = + (2)
R R
Figura 1. Modelos Reolgicos: a) Estado Incial; b) Surfactante Insoluble (A/A=a/a) y c)
Surfactante Soluble (A/A>a/a).
1 A
t
= ; = dt (3)
A t 0
1 1 a
t
= = ; = dt (4)
t a t 0
= EG + d (5)
= EG = (6)
1
= EG . = = (t ) (t = 0) = .( Pc(t ).R(t ) Pc(0).R(0)) (7)
2
3 METODOLOGA EXPERIMENTAL
Se selecciono como fase dispersa, heptano, con un grado de pureza del 99+% HPLC de
SIGMA-ALDRICH, y dodecano de Fisher Scientific. Para la fase continua se seleccion
Tritn X-100 extra puro de laboratorios Aldrich preparado en mezcla con agua destilada a tres
concentraciones diferentes de Tritn X-100; dichas concentraciones fueron de 10-3 mol/l, 10-4
mol/l y de 10-5 mol/l.
(a) (b)
Figura 2. a) Vista frontal del equipo de trabajo, b) Visualizacin y grabacin del
experimento.
Las imgenes fueron obtenidas por una cmara de video (CCD-72) que est colocada
sobre un sistema de rieles que permiten el movimiento de la misma en tres direcciones, a la
cual se le adapt un lente Nikkon 60 mm con un anillo de extensin PK-13, que se conecta a
un monitor Sony (PVM-2030) (Fig. 2b).
Las imgenes temporales de la gota (Fig. 3), finalmente son digitalizadas en una
computadora con una tarjeta Targa, la cual convierte en formato digital la imagen grabada y
la almacena como un archivo del tipo *.tga, al cual se le pueden hacer modificaciones
(limpieza de la imagen), que permiten transformar las imgenes a formato *.bmp. En la Figura
4 se puede observar el diagrama general del equipo experimental.
Actuador
Jeringa con
alcano
Cmara Celda de
visualizacin
Alta Transductor de
Lmpara presin
Baja
Monitor
Solucin
surfactante
/ t (mN/m.s)
TID (m N /m ) 0.30
25 C = 10-4 mol/l
0.25
20 0.20
C = 10-3 mol/l C = 10-3 mol/l
15
0.15
10 C = 10-4 mol/l 0.10
5 C = 10-5 mol/l 0.05
0 0.00
0 2 4 6 0.0 1.0 2.0 3.0 4.0
t (s) t (s)
(a) (b)
Figura 5. a) Tensin interfacial dinmica, b) Tasa de variacin de la tensin interfacial, del
dodecano en interfase con solucin acuosa de Tritn X-100 a diferentes concentraciones.
En funcin de lo anterior se puede especular que el zigzagueo de las curvas de la Fig. 5b,
se deben a la lucha entre el instinto adsortivo del surfactante, reflejado en disminuciones de
la tasa de crecimiento de la tensin y, el efecto de dispersin molecular en la interfase (visto
como picos en la Fig. 5b), que genera la gota durante su crecimiento dentro de la solucin.
45 0.70
40 0.60
35 C = 10-5 mol/l
/t (mN/m.s)
0.50
T ID (m N /m )
30 C = 10-4 mol/l
25 0.40
C = 10-3 mol/l
20 C = 10-3 mol/l 0.30
15 0.20
10 C = 10-4 mol/l
5 C = 10-5 mol/l 0.10
0 0.00
0 1 2 3 4 0.0 1.0 2.0 3.0 4.0
t (s) t (s)
(a) (b)
Figura 6. a) Tensin interfacial dinmica, b) Tasa de variacin de la tensin interfacial, de
heptano en interfase con solucin acuosa de Tritn X-100 a diferentes concentraciones.
Se observa de la Fig. 6b que para bajas concentraciones el incremento en la TID tiene una
conducta diferente al de las concentraciones mayores. Esta conducta, se caracteriza por una
alta tasa de aumento de la TID en los instantes siguientes a la formacin de la interfase fresca
y la drstica disminucin de la misma en poco tiempo, que luego se hace prcticamente
constante. Para las concentraciones mayores se observa que la variacin es mnima en lo que
respecta al aumento de la TID, lo cual se traduce en la conducta lineal de incremento de la
tensin interfacial dinmica que se observa en la Fig. 5a, para estas concentraciones. El
anlisis utilizado en el caso del dodecano tambin es aplicable en este caso.
Los resultados obtenidos con el montaje experimental concuerdan con el resultado
obtenido en el caso del heptano por Gbel y Joppien (1997) quienes utilizaron el tensimetro
de volumen de la gota para obtener la TID del heptano en concentraciones de 10-5 y 10-4 mol/l
de solucin acuosa de Tritn X-100. Ellos obtuvieron en su experiencia resultados similares a
los reportados en este trabajo. La conducta en lo que a concentracin del medio se refiere, fue
la misma. En el caso de la concentracin menor, el rango de valores en el cual se mueve la
TID es distinto (de 22 a 45 mN/m). Esta diferencia puede tener su origen en caractersticas
inherentes a la tcnica de medicin como tambin a condiciones experimentales distintas, tales
como, la temperatura de operacin, tamao del capilar, caudal de inyeccin del alcano, entre
otros.
Es completamente factible que el caudal de inyeccin sea el factor que establece la
diferencia entre estos rangos. Durante el crecimiento de la gota imperan los procesos de
dispersin molecular en la interfase y la formacin de la monocapa en la superficie, pero a
medida que transcurre el tiempo el efecto tensoactivo del surfactante comenzar su accin (de
disminucin de la TI) opacando el efecto de dispersin molecular.
4.3 Comparacin entre el comportamiento de los dos alcanos en las distintas soluciones
surfactantes
30
70 25
20
60
15
50 Dodecano 10
Heptano 5
40
0
30 1.00E-06 1.00E-05 1.00E-04 1.00E-03 1.00E-02
1.0E-06 1.0E-05 1.0E-04 1.0E-03 1.0E-02 Concentracin de la solucin (mol/l)
Concentracin de solucin surfactante Dodecano TID reol Dodecano TID cap Dodecano T
(mol/l) Heptano TID reol Heptano TID cap Heptano TID
(a) (b)
Figura 7. (a)Porcentaje de variacin neta del valor de tensin interfacial de dodecano y
heptano en interfase con solucin acuosa de Tritn X-100 a diferentes concentraciones, (b)
Contribuciones a los valores de TID promedio.
5. CONCLUSIONES
6. REFERENCIA
Binks, B., Cho, W., Fletcher, P y Petsev, D, Stability of oil-in-water emulsions in a low interfacial tension
system, Langmuir, Vol. 16, N 3, pp. 1025 - 1034, (2000).
Bofillon, A., Sicoli, F y Langevin, D, Dynamic surface tension of ionic surfactant solutions, Journal of colloid
and interface science, Vol. 168, pp. 497 - 504, (1994).
Edwards, D., Brenner, H y Wasan, D, Interfacial transport processes and rheology, Butterworth-Heinemann,
Stoneham, M.A, pp. 1051 - 1064, (1991).
Goebel, A y Lunkenheimer, K, Interfacial tensin of the water/n-alkane interfacial, Langmuir.,Vol. 13, pp. 369
372, (1997).
Horozov, T y Joos, P, Dynamic Surface Tension of Surfactant Solutions Studied by Peaktensiometry, Journal
of Colloid and Interface Science, Vol. 173, pp. 334 - 342, (1995).
Horonov, T., Dushkin, C., Danov, K., Arnaudov, L., Velev, O., Mehretead, A y Broze, G, Effect of the surface
expansion and wettability of the capillary on the dynamic surface tension measured by the maximum bubble
pressure method, Colloids Surfaces, Vol. 113, pp. 117 - 126, (1996).
Horozov, T y Arnaudov, L, A novel fast technique for measuring dynamic surface and interfacial tension of
surfactant solutions at constant interfacial area, Journal of Colloid and Interface Science, Vol. 219, pp. 99 -
109, (1999).
Kim, Y., Koczo, K y Wasan, D, Dynamic film and interfacial tensions in emulsion and foam systems, Journal
of Colloid and Interface Science, Vol. 187, pp. 29 - 44, (1997).
Liggieri, L., Ravera, F y Passerone, A, Dynamic interfacial tension measurements by a capillary pressure
method, Jounal of Colloid and Interface Science, Vol. 169, pp. 226 - 237, (1995).
MacLeod, C y Radke, C, A growing drop technique for measuring dynamic interfacial tension, Journal of
Colloid and Interface Science, Vol. 160, pp. 435 - 448, (1993).
MacLeod, C y Radke, C, Charge Effects in the Transient Adsorption of Ionic Surfactants at Fluid Interfaces,
Langmuir, Vol. 10, pp. 3555 - 3566, (1994).
Miller, R., Joos, P y Fairnerman, V, Dynamic studies of soluble adsorption layers, Progressing Colloid
Polymeric Science, Vol. 97, pp. 188 - 193, (1994).
Nagaranjan, R y Wasan. D, Measurement of dynamic interfacial tension by an expanding drop tensiometer,
Journal of Colloid and Interface Science, Vol. 159, pp. 164 - 173, (1993).
Touhami, Y., Neale, G., Hornof, V y Khalfalah, H, A modified pendant drop method for transient and dynamic
interfacial tension measurement, Colloids and Surface A:Physicochemical and Engineering Aspects, Vol.
112, pp. 31 - 41, (1996).
Thermal conductivity and surface tension experimental values for
three alcohol binary mixtures: heptanol-octanol, heptanol-nonanol
and octanol-nonanol
ABSTRACT
Surface tension was measured for different alcohol binary mixtures (heptanol-octanol, heptanol-
nonanol and octanol-nonanol) using the pendant drop technique. The experimental equipment
assembled consisted of a stainless steel chamber, an RGB monitor, VHS, computer and software
to calculate surface tension from the pendant drop profile. The composition of the mixture was
varied between 0 and 100% v/v. The density values for pure alcohols compare satisfactorily with
those reported in the literature. A non-linear behavior was observed in the surface tension
values as a function of volumetric composition. Maximum and minimal values were observed
depending on the alcohol binary mixture. For the same alcohol mixture, thermal conductivity
were experimentally measured using coaxial cylinder equipment with an accuracy of 0.09
W/m.K. Thermal conductivity for heptanol, octanol and nonanol experimentally obtained
matched with theoretical values reported in the literature with less than 2% of deviation. Similar
behavior to surface tension was obtained with the thermal conductivity, since maximum and
minimal values were detected depending on the alcohol binary mixture. In the case of octanol-
nonanol, two maximum and minimal values for surface tension and thermal conductivity were
obtained for the same composition: 20-30% v/v and 70-80% v/v. These observations suggested
that there could be a relationship between these two properties, the capacity of a substance to
conduct heat through the liquid is related with the surface tension forces between the liquid
phase and its vapor. This nonlinear behavior on thermal conductivity could have an application
on heat transfer equipment, since it is then possible to create a mixture more conductive from
pure liquid than are less conductive.
The knowledge of the transport properties of different materials is of great importance for the
work of the chemical engineer as the above mentioned properties determine the viability of a
process, its costs and yields.
The heat transfer phenomenon by conduction is one of the most frequent in the processes the
chemical engineer designs and controls. To quantify heat transfer, it is necessary to know the
Fourier Laws proportionality constant, to measure heat flow. This constant is a thermal
property. For the majority of pure substances, the constants have been reported, but for mixtures
the information is rather scarce.
The surface tension is another phenomenon of great importance in the chemical industry.
This property defines the surface forces present in the gas-liquid interface. The study of the
surface tension started more than two centuries ago and there is considerable information
available for pure substances, which is not the case for mixtures. The fact that the two properties
have their origin in the same molecular interactions, mainly the Van der Walls forces, favored
the conjunct measurement of both properties, using different alcohol mixtures: heptanol-octanol,
heptanol-nonanol and octanol-nonanol for concentrations between 0 and 100% v/v. A heat
exchanger with concentric pipes of small annular region was used for the determination of the
thermal conductivity. The fundamental objective of this work is to compare the behavior of both
phenomena and to establish some sort of common relationship.
FUNDAMENTAL PRINCIPLES
Thermal Conductivity
The process of heat transfer is explained by the Fourier Law. After numerous experimental
tests, it was determined that the speed of heat transfer is directly proportional to the cross transfer
area and the difference in temperatures, while it is inversely proportional to the length of
transfer. The proportionality becomes equality including a coefficient which is a measure of the
behavior of the material, named material conductivity coefficient, K. The Fourier Law is then as
follows:
T
qx = K A (1)
x
Where qx is the speed of heat transfer, and K is the material thermal conductivity in units W/(m
K). Evaluating this expression when x tends to zero the following differential expression is
obtained
dT
qx = K A (2)
dx
The negative sign is due to the fact the heat is transfer in the direction of the decreasing
temperature.
Using this expression, the heat flow by unit of area can be deducted
q dT
qx = x = K (3)
A dx
It is important to highlight that this equation, in the presented form, only considers
unidirectional heat transfer, normal to the cross area (Incropera and Dewitt, 1996). Nevertheless,
this equation is applicable to this experiment, if the conduction effects are discarded in any
direction other than radial.
The surface tension is considered an extensive thermodynamic property which depends of the
interfacial area of the system. The weak Van der Walls molecular attraction forces are used to
explain the superficial and interfacial tension phenomena. The molecules in a surface are not
completely surrounded by similar molecules, and because of that, the interaction forces between
them are anisotropic and inferior in magnitude to the net attraction forces that take place between
the molecules that are in the center of the substance in question. Considering a molecule of m
mass in the center of a fluid, this embraces the same number of molecules in all directions, so all
the molecular actions experimented by the m molecule balance; but, for a molecule in the liquid
surface there will be an unbalance between the forces that act on it. Perpendicularly to the
surface, a force will exist, which is translated into a molecular pressure impossible to be
measured directly, tangentially there will be a chain of surface molecules united by their mutual
attraction, which effect will be the manifestation of a surface tension, expressed in milligrams
force divided by linear millimeter of surface.
Surface tension
A pendant drop is a technique that uses the shape of a drop and the Laplace-Young equation
to compute the surface or interfacial tension of gas-liquid or liquid-liquid systems. If there is no
gravity the shape of the drop is perfectly spherical (Rabockai, 1979). This method is one of the
most accurate techniques to predict surface tension for liquids (Garandet et al., 1994). The
accuracy of the methods is 0.1 N/m (Lin et al. 1995) and in certain cases of 0.1% (Touhami et
al. 1996). One of its advantages is that a small quantity of liquid is required (Adamson, 1967).
The first work was done by Worthington (1881) and since then this technique has been
applied successfully to calculate surface and interfacial tension from the shape of a drop
(Fergunson, 1912; Andreas et al., 1938).
Nowadays with the advances in optics, electronics and computation process, it is possible to
capture the drop shape and from its curvature estimate the surface tension with a precision
around 0.1%. The results are independent to the contact angle since the mathematical model does
not need any boundaries condition in the region near the place where the drop suspends. Another
advantage of these methods is the small quantity of liquid that it needs (Andreas et al. 1938)
The theory of the pendant drop is based in the balance between the gravitational and the
flotation forces. The mechanic equilibrium of the drop is given by the Laplace-Young equation.
Al the calculus variations that make use of this method must solve this equation, which related
the difference or change in pressure in any point over an interface of two non miscible liquids or
between a liquid and a gas phase.
The drop profile corresponds to a convolution solid, in other words, the drop is a cylindrical
surface that has a symmetry center and can be generated by the rotation of a curve in three-
dimensional space, so the Laplace-Young equation can be expressed as:
d 2z dz
2 g z dx 2 dx
= +
R0 3 1
(11)
dz 2
2 dz 2 2
1 + x 1 +
dx
dx
This is a second order differential equation for the form of a meniscus in equilibrium found
between two liquids of uniform density. Its analytical integrations are up to date unknown. To
ease the resolution of this equation, this equation is written in function of an arc of S length, in
spherical coordinates, knowing that:
dZ
= sen (12)
dS
dX
= cos (13)
dS
d sen
= 2 z (14)
ds x
Where,
gR 0 2
= (15)
Where is the difference in density between the drop and the surrounding media and is
the shape factor. These equations are subject to the boundary conditions:
The numeric integration of the Laplace-Young equations was done for the first time by
Bashforth and Adams (1883), years later was recalculated by Fordham (1948), Niederhauser and
Bartell (1948). Later on, complex methods were developed, taking advantage of computer, as
Padday and Pitt (1971), Girault et al. (1982), Rotemberg et al. (1983), Huh amd Reed (1983),
Jenning and Pallas (1988), Pallas and Harrison (1990) and Rodsrud (1990) and Lpez de Ramos
et al. (1993), who developed a simple and precise method for the determination of surface
tension using the pendant drop technique. This work adjusts the drop profile by means of fourth
order spline functions, used to compute the drop curvature and finally the surface tension.
EXPERIMENTAL SET-UP
The experimental equipment used to calculate the thermal conductivity was a co-cylindrical
cell heated at the center with an electrical resistance and an outer jacket with a water circulation
system to keep the external wall temperature constant (Fig. 1). The fluid is introduced by the
bottom of the cell, after extracting almost all the air with a vacuum pump. The total heat transfer
is equal to d x and the thickness is r, where d is the heat transfer distance, x is the length of
the heat exchange and r is the internal radial distance.
The conductivity cell is connected to a console where voltage, electrical resistance and
temperature are registered. The water for the refrigeration cycle is provided from a tank with
liquid level constant (Fig. 2).
The equations that model the heat process in the cell are very simple and come from an
energy balance in the unit.
Qe = Qc + Qi (17)
Where,
V2
Qe = (18)
R
The calibration curve made with a liquid of known thermal conductivity allows determining
the value of the incidental heat flow (Qi). From Eq. (17):
Qc = Qe Qi (19)
The conductive heat transfer can be calculated from the Fouriers Law:
T
Qc = k A (20)
r
Qc r
k= (21)
A T
For the measurement of surface tension values of the alcohol mixtures at different
concentrations an experimental set-up, similar to the one showed in figure 3, was used.
(4)
(1) (2) (5)
(6)
(3)
(7)
(9)
(8)
(10)
This experimental arrangement consists of a stainless steal cell with two optical glass
windows of 5 mm of thickness, sailed to the cell body with gaskets of 0.65 mm. The mean body
has four channels that come throughout the wall to keep the temperature of the chamber constant.
Also, there are several connections to introduce liquid and gas inside the cell. In the case of the
liquid, a syringe is used to form the pendant drop. The cell also has drainage holes to take out the
liquid of the chamber. More detailed information could be found at Lpez de Ramos (1993).
The cell is fixed to an aluminum support that it is placed on the top of an optical table (with
isolation from outside vibration sources) to avoid perturbations of the liquid drop surface and
obtain a sharp image with accurate contour.
The visualization cell also has a thermocouple and a pressure transducer connected to the
chamber in order to register and control these operation variables during the experiment. The
drop formation is made by means of a syringe with a needle of 0.1 cm of diameter. This syringes
firmly fixed at the top of the cell to avoid any movement during the drop formation, assuring the
vertical position of the needle inside the cell.
The cell is illuminated by its back using a fiber optic lamp. A filter is placed between the
source of light and the cell in order to diffuse the light and get a high quality image of the alcohol
drop (Hansen and Rdsrud, 1990).
At the front side of the cell, a video camera is placed to capture the drop images. This camera
is connected to and xyz linear translator that moves the camera to center and focus the image.
The captured image is sent to a video system to be record in digital and video formats. The
digital drop image is processed by a custom made program (Fuentes, 1998) in order to obtain its
xy profile and compute the surface tension.
EXPERIMENTAL RESULTS
Using the experimental techniques explained in the previous sections, it was possible to
calculate the surface tension and thermal conductivity of the mixtures of heptanol-octanol,
heptanol-nonanol and octanol- nonanol at different concentrations. In this section the final results
are presented.
These data can be expressed in a graphic way in order to have a better idea of the tendency of
thermal conductivity and surface tension with the concentration (figure 4).
31 0,24
Surface Tension (mN)/m
Thermal Condutivity
30
0,235
29
W/(m.K)
28 0,23
27
0,225
26
25 0,22
0 0,25 0,5 0,75 1
Volumetric % (v/v)
Figure 4: Thermal conductivity and surface tension of the binary mixture of heptanol-octanol.
On tables 4 and 5, it is shown the experimental values obtained in this work for thermal
conductivity and surface tension.
Table 4: Thermal conductivity and surface tension of heptanol-nonanol mixture.
31 0,23
Surface Tension (mN)/m
Thermal Condutivity
30 0,22
W/(m.K)
29 0,21
28 0,2
27 0,19
26 0,18
0 0,2 0,4 0,6 0,8 1
Volumetric % (v/v)
Figure 5: Thermal conductivity and surface tension for the binary mixture of heptanol-nonanol.
On tables 5 and 6, it is presented the experimental results for the binary mixture of octanol-
nonanol as a function of the nonanol concentration (% v/v).
Table 5: Thermal conductivity and surface tension of the binary mixture of octanol-nonanol
Thermal Condutivity
0,24
34
0,22
W/(m.K)
32
0,2
30
0,18
28
26 0,16
24 0,14
0 0,2 0,4 0,6 0,8 1
Volumetric % (v/v)
Surface Tension Thermal Condutivity
Figure 6: Thermal conductivity and surface tension of the binary mixture of octanol-nonanol.
DISCUSSIONS
A quick inspection of the experimental results, allows observing a similar behavior in the
surface tension and thermal conductivity properties, when the composition of the alcohol
mixtures A-B is varied from 0 (pure A) to 100% (pure B). None of these tendencies are
represented by any empirical correlation or thermodynamic equation (Yaws, 1977, Reid, 1977,
Perry and Grenn, 1997). Those correlations predict a lineal behavior for alcohol mixtures;
nevertheless the experiment results obtained in this work reveal something very different. An
extensive experimental work has been done since 2000 at the Transport Phenomena Laboratory
at Universidad Simon Bolivar (Materan and Sanabria, 2000, Saume et al., 2001, Santelli et al.,
2001). In all these experimental work, a non lineal behavior was permanently observed, not only
in thermal conductivity data, but also in surface tension experimental values.
Observing the thermal conductivity and surface tension for the pure alcohols, there is an
increment of these properties as the number of carbons increases. This behavior can be translated
to a stronger dispersion forces due to the size of the aliphatic chain. In figure 4, it can be
observed the presence of two maximum clearly differentiated around the same composition
values of 25 % and 75 % v/v of octanol (xoctanol=0.25, xoctanol=0.75). The thermal conductivity
values of all the mixture of heptanol-nonanol studied are bigger than the pure alcohols.
In figure 5, it can be observed the presence of maximum and minimal values, presenting a
non linear behavior for thermal conductivity and surface tension properties. The observed
maximum is in a point near by 25% v/v (xnonanol=0.25). A maximum value of surface tension
could indicate a short separation between the molecules (short intermolecular distance), favoring
the transport of heat by conduction along the alcohol mixture. This enhancing of the conductive
heat flow is reflected in a higher value of thermal conductivity at this point. In the alcohol
mixtures with small percentage of nonanol, probably the interactions of molecules favor the
reduction of the intermolecular distances, allowing a better and more compact molecular
arrangement and a higher conduction of the heat between them. This compact arraignment of the
molecules could be demonstrated by an increment of the surface tension of the alcohol mixture.
The experimental results obtained in this work for the binary mixture of octanol-nonanol are
reported in figure 6 as a function of the volumetric fraction of nonanol (xnonanol). As it can be
observed, thermal conductivity and surface tension of pure alcohols are lower than the mixtures
at any intermediate composition, being octanol (xnonanol = 0, k = 0.176 0.003 W/(m.K), = 27.1
0.3 mN/m) lower than the nonanol (xnonanol =1, k = 0.184 0.004 W/(m.K), = 27.9 0.1
mN/m).
A possible explanation of this non linear behavior is that the alcohols have always been
considered highly associated in the liquid phase. This association takes into the formation of
polymeric chains. If the attraction forces are big enough, the minimal potential energy of
interaction will be more negative than the mean thermal energy. The equilibrium stage of the
system will have a considerable quantity of dimmers. In alcohols, the hydroxyl groups are able to
form dipoles and hydrogen bridges.
The geometry of the molecular structure of the alcohol at the interface is controlled by the
van der Waals forces that act over the aliphatic chain (Adamson, 1960). This arrangement fixes
the orientation of the dipole of the hydroxyl group and origins an increment of the surface
potential at the interface, because the non-parallel-alignment to the interfacial plane. This
behavior could explain the presence of maximum and minimal values in the curves of surface
tension as a function of concentration due to variation of surface potential energy at the interface.
This hypothesis pretends to explain the appearance of maximum and minimal values in the
curves of thermal conductivity and surface tension experimentally obtained for the binary
mixtures of alcohol selected in this work. Nevertheless there is no more than a possible
explanation that needs to be corroborated in the near future.
CONCLUSIONS
Thermal conductivity and surface tension of binary alcohol mixtures present a non linear
relationship with concentration.
Thermal conductivity of binary alcohol mixtures presents values that are bigger than the
correspondent pure elements.
Thermal conductivity curves as a function of the mixture composition present maximum and
minimum values. The same behavior was found for surface tension of the binary solutions.
Those maximum and minimum property values tend to appear at similar composition of the
alcohol mixtures.
In the case of the measurement of thermal conductivity, the difference of temperature applied
and humidity of the air was strictly controlled during al the experiments.
It can be concluded that as the carbon chain of the alcohol increases, thermal conductivity values
of the mixtures also increases.
NOMENCLATURE
k: Thermal conductivity [W/(m.C)]
A: Cross area [m2]
T: Temperature [C]
qx: Heat flow rate in x direction [W]
qx: Heat flow per unit of area [W/m2]
V: Voltage [V]
R: Electrical resistance of the equipment H470 []
T: Temperature difference [C]
T1: Temperature of the shell [C]
T2: Temperature of the jacket [C]
r: Radial distance [m]
Qi: Incidental heat [W]
Qc: Conduction heat [W]
Qe: Electrical flow [W]
ACKNOWLEDGEMENT
The authors want to thanks the financial supports of Decanato de Investigaciones y Desarrollo of
Universidad Simon Bolivar and FONACIT. Also wants to recognize the enthusiastically help
given by Marysol Materan, Viviana Sanabria, Miguel Flores, Leopoldo Duarte, Cristina
Laveglia, Luis Aguirre, Rodolfo Santelli, Cristina Leon, Marjorie Nogueira, Carolina Pena y
Erika Saume.
REFERENCES
Adamson, A. W., Physical chemistry of surfaces. Inter Science Publishers, New York, pp. 4-6,
14, 28, 30-33, 37 (1967).
Andreas, J. M., Hauser, E. A., Tucker, W. B., Boundary tension by pendant drops, J. Physical
Chemistry., 42, pp 1001-1018, (1938).
Bashforth, F and Adam, J, An attempt to test the theories of capillary action by comparing the
theoretical and measured forms of drops of fluid, Cambridge Uni. Press, c. 107, (1883).
Fergunson, A, Photographic measurement of pendent drops, Philosophical Magazine, Vol. 23,
pp. 417 - 431, (1912).
Fordham, S, On the calculation of surface tension from measurements of pendant drops, Journal
Fluid Mechanical, Vol. 194, pp. 1 - 15, (1948).
Garandet, J. P., Vinet, B and Gros, P., Considerations on the pendant drop method: A new look
at Tates Law and Harkins correction factor, Journal of Colloid and Interface Science, 165, pp.
351-354, (1994).
Girault, H, Schiffrin, D and Smith, B, Drop image processing for surface and interfacial tension
measurement, Journal of Electroanalitical Chemical., Vol. 137, pp. 207 - 217, (1982).
Hansen, F and Rodsrud, G, Surface tension by pendant drop, Journal of Colloid and Interface
Science, Vol. 141, N 1, pp. 1 - 9, (1990).
Huh, C and Reed, R. L, A method for estimating interfacial tension and contact angles from
sessile and pendant drops shapes, Journal of Colloid and Interface Science, 91 (2), pp. 470-484,
(1983).
Incropera, F and Dewitt, D, Fundamentos de Transferencia de Calor. Prentice Hall; Mxico,
(1996).
Jennings, J. W and Pallas, N. R, An efficient method for the estimation of interfacial tension of
interfacial drop profiles, Langmuir, 4 (4), pp. 960-967, (1988).
Lin, S, An examination on the accuracy of interfacial tension measurement from pendant drop
profiles, Langmuir, Vol. 11, pp. 159 4166, (1995).
Lopez de Ramos, A., Redner, R. and Cerro, R., Surface tension from pendant drop curvature,
Langmuir, 9, 3691-3694, (1993).
Niederhauser, D. O and Bartell, F. E, Report of progress fundamental research on occurrence and
recovery of petroleum, American Petroleum Institute, Baltimore, pp. 114-146, (1948).
P.A. Hilton Ltd. Manual Experimental de Operacin y Mantenimiento de la Unidad de
Conductividades Trmicas de Lquidos y Gases. England.1989.
Padday, J. F and Pitt, A., Axisymetric meniscus profiles, Journal of Colloid and Interface
Science, 38, pp. 323-334, (1971).
Pallas, N. R and Harrison, Y, An automated drop shape apparatus and the surface tension of pure
water colloids and surfaces, 43, pp.169-194, (1990).
Parfitt, G. D., Trottman-Dickerson, A. F., Chemical kinetic surface & colloid chemistry,
Pergamon Press, Oxford, pp. 79-80, 82, 86, (1966).
Perry, R and Grenn, D, Perrys Chemical Engineers Handbook. Mc Graw Hill, Estados Unidos,
(1997).
Rabockai, T., Fisico-quimica de superficies, Monografa N 20, Secretaria-Geral da Organizao
dos Estados Americanos. Programa Regional de Desenvolvimiento Cientfico e Tecnolgico.
Serie de qumica, Washington D. C., pp. 7, 8, 18-19, (1979).
Rotenberg, Y., Boruvka, L and Newmann, A. W, Determination of surface tension and contact
angle from the shape of axisymmetric fluid interfaces, Journal of Colloid and Interface Science,
93 (1), pp. 169-183, (1983).
Touhami, Y., Neale, G., Hornof, V and Khalfalah, H, A modified pendant drop method for
transient and dynamic interfacial tension measurement, Colloids and Surface A:
Physicochemical and Engineering Aspects, Vol. 112, pp. 31 - 41, (1996).
Worthington, A, On pendant drops, Proceeding of the Royal Society of London, Vol. 32, pp.
362-377, (1881).