SciELO - Scientific Electronic Library Online

 
vol.22 número2Efecto del tratamiento con explosivos sobre las tensiones residuales en uniones soldadas de acero al carbonoParámetros elásticos y acústicos Para la predicción de fluidos y litología en el lago de Maracaibo. índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay articulos similaresSimilares en SciELO

Compartir


Revista de la Facultad de Ingeniería Universidad Central de Venezuela

versión impresa ISSN 0798-4065

Rev. Fac. Ing. UCV v.22 n.2 Caracas jun. 2007

 

Study of an aqueous foam flow in horizontal pipes.

NÓLIDES GUZMÁN, PH.D.*, OVADIA SHOHAM, PH.D.**, RAM MOHAN, PH.D.***

*Escuela de Ingeniería Química, Facultad de Ingeniería, Universidad Central de Venezuela. Venezuela.

**Petroleum Engineering Department, The University of Tulsa, U.S.A.

***Mechanical Engineering Department, The University of Tulsa, U.S.A.

ABSTRACT

Several sets of experimental data were acquired for an aqueous two-phase foam flow in horizontal pipes, including flow pattern maps and rheological flow behavior. The data were acquired in two pipe diameters under various inlet flow conditions, utilizing a commercial surfactant, F-450, at two concentrations. The experimental results were plotted on a Taitel and Dukler (1976) flow pattern map predicted for the experimental flow conditions. As this model does not predict accurately the transition boundaries for foam flow, modifications of the original model were carried out for matching the experimental data. Different rheological models were tested for characterizing the behavior of the studied aqueous foam. The results obtained show that this foam flow behaves as non-Newtonian flow, following the Power Law model. It was also observed that foam rheology is flow pattern dependent. Comparison between the developed models and the acquired experimental data shows a good agreement.

Keywords: foam, foam rheology, aqueous foam, two-phase foam flow, foam flow patterns.

Estudio del flujo de una espuma acuosa en tuberías horizontales.

RESUMEN

Varios grupos de datos experimentales fueron obtenidos para el flujo bifásico de espuma en tuberías horizontales; estos incluyen: mapas de patrones de flujo y comportamiento reológico del flujo. Los datos fueron obtenidos para dos diámetros de tubería a varias condiciones de flujo de entrada, utilizando el surfactante comercial F-450 a dos concentraciones. Los resultados experimentales fueron representados en el mapa de patrones de flujo de Taitel y Dukler (1976) predicho a las condiciones experimentales de flujo. Como este modelo no predice apropiadamente los límites de transición para flujo de espuma, se requirieron modificaciones del modelo original para ajustarlo a los datos experimentales. Se probaron diferentes modelos reológicos para caracterizar el comportamiento de la espuma acuosa estudiada. Los resultados muestran que este flujo de espuma se comporta como no Newtoniano, siguiendo el modelo de la Ley de Potencia. También se observó que la reología de la espuma depende del patrón de flujo. Las comparaciones entre los modelos desarrollados y los datos experimentales muestran buen ajuste.

Palabras clave: espuma, reología de espumas, espuma acuosa, flujo bifásico de espuma, patrón de flujo de espuma.

Recibido: noviembre de 2006 Recibido en forma final revisado: abril de 2007

INTRODUCTION

Foams are useful in many industrial applications, such as in paper processing, fire-fighting, and enhanced oil recovery (EOR). However, in other processes the production of foam is often an unwanted consequence. A proper design of a foam flow system requires an understanding of the physical phenomena involved. The nature of the foam, i.e., its structure, foamability, stability, and its overall rheological properties can determine both  the economical and technical successes of an industrial process.

Several models have been developed for foam characterization to relate the different effects observed on foams. Scientists who have studied foamy fluids rheology have used different models for foam characterization. Utilizing the model that better represents the foam flow behavior enables more accurate prediction of the pressure drop. This requires that all the characteristic parameters of the foam are known.

Rheology can be evaluated using different equipments. Among the most commonly used are: Couette rheometer, circulating pipe rheometer and single-pass pipe rheometer. In the present study, a single-pass pipe rheometer was used for foam characterization due to its simple design and suitability for the flow conditions that were covered in this research.

Foam flow in straight pipes has been the subject of several investigations, such as: Calvert and Nezhati (1986); Calvert (1990); Enzendorfer et al. (1995); Hanselmann and Windhab (1996), and Deshpande and Barigou (2000). However, no specific correlation between the flow conditions and the rheology of the foam has yet been established. This is the gap the present study attempts to fill.

EXPERIMENTAL WORK

This study was carried out at the Tulsa University Separation Technology Projects (TUSTP). The rheology section, installed in a two-phase flow facility, consisted of 56 ft long 1 in. and 2 in.-diameter horizontal pipes (Figure 1). Two absolute pressure transducers are placed one at each end of the rheology section, for measurement of the pressure at the inlet and at the outlet of the pipes, enabling determination of the foam flow pressure drop. Also, a visualization window is installed in each pipe for flow pattern observations.

Water and surfactant are mixed in a tee to form the surfactant solution. This solution and compressed air are combined in a mixing-tee to promote foam formation. The generated foam flows through a static mixer to ensure homogeneity of the flow. The foam then flows through the rheology section where all measurements are acquired. Downstream, the flow enters into a separation tank where the air is vented into the atmosphere and the surfactant solution is drained.

Data were acquired for various inlet flow conditions: superficial gas velocity between 20 and 70 ft/s and superficial liquid velocity between 0.25 and 0.85 ft/s, using commercial surfactant (F-450) at two concentrations of 0.02 and 0.04% vol./vol..

Foam Flow Pattern Maps

Experimental flow pattern maps for aqueous foam flow in the 2 in.-diameter horizontal pipe of the rheology section were obtained. Three different flow patterns were observed: stratified, slug and annular flow. Figures 2 and 3 show the flow pattern maps for the two surfactant concentrations, namely, Cs = 0.02% and Cs = 0.04 % vol./vol., respectively.

As can be seen, the results for the two surfactant concentrations are similar. Also, the Taitel and Dukler (1976) model (represented by the solid lines) does not predict accurately the stratified to non-stratified and the slug to annular transition boundaries for foam flow. Modifications of the original model are required for matching the experimental data.

Note that flow patterns were also observed in the 1 in.-diameter rheology pipe section. For this case, only annular flow was observed for all flow conditions.

Foam Rheology in Pipes

The experimental data for the 2 in.-diameter pipe for the Cs = 0.02% vol./vol. case are presented in figure 4 in the form of the measured pressure drop (Δp) on the y-axis vs. the superficial gas velocity (vSG) on the x-axis, whereby the parameter is the superficial liquid velocity (vSL). An error analysis for the pressure drop measurement was carried out and the corresponding uncertainty is included with the reported results in Figure 4. As expected, the pressure drop increases when either the gas or the liquid flow rates are increased. For the analysis of the rheological behavior of the foam flow, the experimental data are re-plotted in Figures 5 and 6. Figure 5 presents the experimental data in the form of the wall shear stress vs. shear rate. The shear stress is determined as tw = d/4xdp/dL = d/4xΔp/L, where Δp is the pressure drop across the rheology section pipe and L is the distance between the two pressure transducers. As can be seen from the figure, foam flow does not behave as a Newtonian fluid.

The experimental data are plotted in figure 6 in the form of natural log of wall shear stress vs. natural log of shear rate. The figure demonstrates that there are two different regions; one corresponding to stratified flow conditions while the other is for non-stratified (slug and annular) flow conditions. Within each region, the flow behaves according to the Power Law model, as the data exhibit straight lines, with different slopes for the stratified and the non-stratified cases.

Similar results are obtained for Cs = 0.04%  vol./vol., whereby the data for this case, too, demonstrates that foam flow behaves as a Power Law fluid, with different slopes for the stratified and non-stratified cases, and different slopes as compared with the lower concentration results given in Figure 6.

MODELING

Foam Rheology in Pipes

Various rheological models were tested against the experimental data of foam flow acquired in this study. The best correlation with the data was obtained with the Power Law Model, as shown in figure 6.

For the Power Law Model, the shear stress vs. shear rate is given by:

(1)

where Kf is the consistency index and n is the flow behavior index.

For this study, Kf is determined from the experimental data as follows:

(2)

where the wall shear stress can be determined from the experimental pressure gradient data, namely, tw = d/4xdp/dL

The flow behavior index, n, can be obtained from the slope of the curve, namely, ln(tW) vs. ln(8vM/d), as presented in figure 6. As only a slight difference occurs between n for stratified and non-stratified flow, an average n value is obtained for each surfactant concentration: n = 1.18 for Cs = 0.02%, and n = 1.18 for Cs = 0.04%

The following equations were obtained by correlating the experimental data, using the form of the Law of the Wall. The consistency index for the stratified flow regime is given by:

Kfs = as In(ReM) + bs,  (3)

which is valid for the range:

4x103 < ReSL < 2x104 and 103 < ReM < 106

Similarly, for non-stratified flow:

KfNS = aNS In(ReM) + bNS + ReMf, (4)

valid for the range 8x103 < ReSL < 3X104 and 103 < ReM < 107.

The coefficients ai and bi for both the stratified and nonstratified flow conditions can be determined from Eqs. 5 and 6, and the exponent f =-3.508 is a constant, for all cases:

ai = Ai . In(ReSL) + A2,  (5)

bi = Bi . In(ReSL) + B2, = (6)

The coefficients Ai and Bi are flow pattern dependent and are given in Table 1 for both Cs = 0.02% and Cs = 0.04%.

The superficial liquid Reynolds number and the mixture Reynolds number are defined, respectively, by

(7)

  (8)

The definition of the mixture Reynolds number is similar to the one presented by Briceño and Joseph (2003) for selflubrication transport of aqueous foams.

Flow Pattern Prediction

Stratified to non-Stratified Transition

Figures 2 and 3 show that Taitel and Dukler (1976) model underpredicts this transition boundary, namely, the predicted transition occurs at lower superficial velocities, as compared with the experimental data.

For the transition between stratified to non-stratified flow, Taitel and Dukler (1976) used a smooth interface approach, namely, fI = fG . For stratified wavy flow, a constant interfacial friction factor was proposed, i.e., fI = 0.142 (Shoham and Taitel, 1984). In foam flow, a higher interfacial friction factor occurs due to the rough foamy interface. In this study, the following interfacial friction factor is proposed for foam flow conditions.

fI = 0.02  (9)

Slug to Annular Transition

As can be seen in figures 2 and 3, this transition boundary is also underpredicted by the original model. The original Taitel and Dukler (1976) model proposed that this transition occurs when the stratified flow configuration is unstable and the equilibrium liquid level is hL/d = 0.5 This criterion was later modified by Barnea et al. (1980), accounting for an average liquid holdup in the slug, namely, 0.7, yielding hL/d = 0.5x0.7 = 0.35 as the criterion for this transition. However, for foam flow lower values of liquid holdup may occur in the slug. In this study, a value of 0.26 is assumed for the liquid holdup in the slug (corresponding to rhombohedra packing of gas bubbles), resulting in the following criterion for the transition from slug to annular under foam flow conditions.

   (10)

RESULTS

Foam Rheology in Pipes

Figure 7 shows a comparison between the experimental data and the predictions of the extension of the Power Law Model for foam flow for the Cs = 0.02% surfactant concentration case. The y-axis is the pressure drop across the test section, Δp in psi. The x-axis and the parameter are chosen based on the developed correlation as ReM and ReSL, respectively. The solid lines represent the predictions of the proposed correlation for both stratified and non-stratified flow conditions. As can be seen, good agreement is observed between the experimental data and correlation predictions.

To demonstrate the accuracy of the proposed foam rheology correlation, a comparison between the measured pressure drop, ΔpEXP, and the predicted pressure drop, ΔpCOR, was carried out. The measured and predicted pressure drop values agree within % 10 ± for both surfactant concentration cases.

Flow Pattern Prediction

Figure 8 presents a comparison between the predictions of the original Taitel and Dukler (1976) model, the modified model based on Eqs. 9 and 10, and the experimental data. As can be observed, the predictions of the modified model agree closely with the experimental data.

CONCLUSIONS

The power law model was modified for aqueous foam flow.

A correlation was developed for the prediction of the consistency index, Kf, which is flow pattern dependent, namely for stratified flow and non-stratified flow conditions.

The Taitel and Dukler (1976) model was modified for flow pattern prediction in foam flow. The modifications of the model include an interfacial friction factor of fI = 0.02 for prediction of the equilibrium stratified flow and hL/d = 0.13 as the criterion of transition boundary between slug and annular flow.

Comparison between the developed correlation and model modifications with the acquired experimental data shows good agreement.

ACKNOWLEDGEMENTS

The authors acknowledge the financial support of Universidad Central de Venezuela and Tulsa University Separation Technology Projects (TUSTP).

NOMENCLATURE

A, a = coefficient (-)

B, b = coefficient (-)

CS = surfactant concentration (%)

d = diameter (in)

f = rheology exponent (-)

fI = interfacial friction factor (-)

hL = height (in)

Kf = consistency index (psia.sn)

n = flow behavior index

pAV = average pressure (psia)

Re = Reynolds number

T = temperature (°F)

v = velocity (ft/s)

Greek symbols:

Δp = pressure gradient (psi)

mL = liquid viscosity (lbm/ft x s)

q = inclination angle (°)

rL = liquid density (lb/ft3)

tW = wall shear stress (psi)

s = surface tension (dynes/cm)

Subscripts:

EXP = experimental

M = mixture

NS = non-stratified

S = stratified

SG = superficial gas

SL = superficial liquid

Abbreviations:

A = annular flow

mod = model

dat = experimental data

SF = stratified flow

SL = slug flow

REFERENCES

1. BARNEA, D., SHOHAM, O., TAITEL, Y. AND DUKLER, A.E. (1980). «Flow Pattern Transition for Gas-Liquid Flow in Horizontal and Inclined Pipes, Comparison of ( S , T = 85ºF). Experimental Data with Theory», Int. J. Multiphase Flow, 6, pp. 217-225.        [ Links ]

2. BRICEÑO, M.I. AND JOSEPH, D.D. (2003). «Self-Lubricated Transport of Aqueous Foams in Horizontal Conduits», International Journal of Multiphase Flow, 29, pp. 1817–1831.        [ Links ]

3. CALVERT, J.R. (1990). «Pressure Drop for Foam Flow through Pipes», International Journal of Heat and Fluid Flow, 11, No. 3, pp. 236-241.        [ Links ]

4. CALVERT, J.R. AND NEZHATI, K. (1986). «A Rheological Model for a Liquid-Gas Foam», International Journal of Heat and Fluid Flow, 7, No. 3, pp. 164-168.        [ Links ]

5. DESHPANDE, N.S. AND BARIGOU, M. (2000). «The Flow of Gas-Liquid Foams in Vertical Pipes», Chemical Engineering Science, 55, pp. 4297-4309.        [ Links ]

6. ENZENDORFER, C., HARRIS, R.A., VALKO, P., ECONOMIDES, M.J., FOKKER, P.A. AND DAVIES, D.D. (1995). «Pipe Viscometry of Foams» Journal of Rheology, 39, No. 2, pp. 345-358.        [ Links ]

7. HANSELMANN, W. AND WINDHAB, E. (1996). «Foam Flow in Pipes», Journal of Applied Rheology, 6, No. 6, pp. 253-260.        [ Links ]

8. SHOHAM, O. AND TAITEL, Y. (1984). «Stratified Turbulent-Turbulent Gas Liquid Flow in Horizontal and Inclined Pipes», AIChE J., 30, pp. 377-385.        [ Links ]

9. TAITEL, Y. AND DUKLER, A.E. (1976). «A Model for Predicting Flow Regime Transition in Horizontal and Near Horizontal Gas-Liquid Flow», AIChE J., 22, no. 1, pp. 47-55.        [ Links ]