# Flow Patterns Transition Law of Oil-Water Two-Phase Flow under a Wide Range of Oil Phase Viscosity Condition

Authors: Wei Wang; Wei Cheng; Kai Li; Chen Lou; Jing Gong
Journal: Journal of Applied Mathematics
 (2013)

---

## Abstract

A systematic work on the prediction of flow patterns transition of the oil-water two-phase flows is carried out under a wide range of oil phase viscosities, where four main flow regimes are considered including stratified, dispersed, core-annular, and intermittent flow. For oil with a relatively low viscosity, VKH criterion is considered for the stability of stratified flow, and critical drop size model is distinguished for the transition of o/w and w/o dispersed flow. For oil with a high viscousity, boundaries of core-annular flow are based on criteria proposed by Bannwart and Strazza et al. and neutral stability law ignoring that the velocity of the viscous phase is introduced for stratified flow. Comparisons between predictions and quantities of available data in both low and high viscosity oil-water flow from literatures show a good agreement. The framework provides extensive information about flow patterns transition of oil-water two-phase flow for industrial application.


---

## Body

## 1. Introduction

Oil-water two-phase systems abound in many applications in the petroleum industry including emulsions preparation, oil-water mixture separation, and transportation. Accurate prediction of oil-water flow characteristics, such as flow pattern transition, is important in many engineering applications. Despite the importance, however, oil-water flow has not been explored to the same extent as gas-liquid flow. The density difference between the two phases is relatively small, while the viscosity ratio encountered can extend over several orders of magnitude, both of which bring great complexities and difficulties to the study of oil-water flow transitions.Various experimental works have been proposed in the literature to enhance the understanding of oil-water flow transitions [1–4] in horizontal pipes with the flow pattern being classified as stratified, dispersed, and mixed flow. Trallero [4] conducted extensive experimental and theoretical studies on oil-water flow patterns, where segregated and dispersed flow were mainly focused and a total of six patterns were identified.For the transition law of the flow patterns, Brauner and Maron [5] first developed a two-fluid model to characterize the transition from stratified flow to other flow patterns. Trallero [4] proposed that VKH criterion could predict the transition between stratified flow and unstable wavy stratified flow. Brauner [6] suggested an approach for the transition to dispersed flow pattern. Recently, Sharma et al. [7] introduced the energy minimization concept, a principle that a system can be stabilized to its minimum total energy, into oil-water flow pattern transition.Most of the published works focus on low viscosity oil-water flows. However, when the viscosity of oil increases, new characteristics could be observed. For heavy oil-water dispersions, core-annular flow, especially oil core-water annular flow, becomes dominant. Moreover, drop entrainment as well as oil clot and slug can easily be found in such flow.For the transition of core-annular flow pattern, Bannwart [8] argued that core-annular flow is prone to occur in a pipe where the two fluids have very dissimilar viscosities but relatively close densities. In this case, however, it is proposed that velocity difference of the two phases should be small. Grassi et al. [9] found that the extension of theoretical models for low viscosity ratio to high viscosity ratio flows should not be straightforward. Rodriguez et al. [10] developed a refined pressure-loss prediction model which includes a slip ratio term that implicitly accounts for the buoyancy of the oil core.From the previous literature, it could be gained that systematic work for the prediction of flow pattern transition is still demanded. This paper aims to establish a framework for flow pattern transition prediction especially over a wide range of oil viscosity, which is expected to provide extensive support for the industrial application of oil-water flow pattern transition theories. Results are validated with a wide range of existing published data, even that of heavy crude oils.


## 2. Theoretical Models of Flow Pattern Transition

### 2.1. Flow Pattern Classification

Although large amounts of experiments on oil-water two-phase flow in horizontal pipe have been conducted by scholars, there still lacks overall consensus about how the flow patterns could be classified. This may be due to the fact that experimental conditions differ and ways to observe flow regimes are not identical. In spite of that, the common regimes without drop entrainment are mainly considered here. Combine Taitel and Dukler’s gas-liquid flow pattern classification method with oil-water pipe flow characteristics, the flow patterns of low viscosity oil-water two-phase flow are classified as stratified flow, w/o dispersed flow, o/w dispersed flow, and intermittent flow.


### 2.2. Transition of Stratified Flow

#### 2.2.1. Transition of Stratified Flow for Low Viscosity Ratio

Recent researches of stratified flow stability are all based on momentum and continuity equation of double fluid model by introducing the influence of small disturbance on the interface. It is common that ways of treating the stability equations could be classified into two categories: ZNS (zero neutral stability)/ZRC (zero real characteristics) criterion and Kelvin-Helmholz (IKH and VKH) instability theory. Brauner and Moalem Maron [11] used ZNS/ZRC criterion to predict the transition criterion of smooth stratified flow pattern to other flow patterns. ZNS line is set as the boundary of smooth stratified and stratified wavy and ZRC line as the upper limit of smooth stratified. Trallero [4] argued that the VKH equation determines the transition between stable stratified flow and either unstable stratified with wavy interface or some other flow patterns and may be used to predict onset of entrainment, which is suitable for low viscosity oil-water systems. The equations are as follows: (1)(CV-CIV)2+ρwρoρHwHo(Uw-Uo)2-(ρw-ρo)gcosβAρSi-σAρSik2︸IKH+ρf(Uw-Uo)2CSAρ(1Aw-1Ao)︷VKH<0(2)Jμ+JU-Jg-Jσ︸IKH+Js︷VKH<0(3)CV=∂F/∂Hw∂F/∂Uso-∂F/∂Usw(4)F=-τwSwAw+τoSoAo+τiSi(1Aw+1Ao)=0(5)CIV=ρwUwHw+ρoUoHoρwHw+ρoHo(6)ρ=A(ρwAw+ρoAo),where ρw and ρo are the density of water and oil, respectively; Hw and Ho are the holdup of water and oil, respectively; Usw and Uso are superficial velocity of water and oil, respectively; Uw and Uo are real mean velocity of water and oil, respectively; τw, τo, and τi are water, oil, and interfacial shear stresses; Sw, So, and Si are water, oil, and interface wet perimeter; Aw, Ao, and A are cross-sectional area of water, oil, and the pipe; CV and CIV are interface wave velocity and critical interface wave velocity, and k is the wave number.The first item on the left of (1) is the viscous instability, depending on τw, τo, and τi; the second item is the velocity instability, related to velocity difference of the two phases; the third item is the gravity stability; the fourth item is interface tension stability, related to interface tension and interface wave number; the fifth is a correction term which stands for the influence of other interface instable factors. CS is the correction coefficient.In order to apply the stratified flow criterion, the in-situ water and oil velocities as well as cross-section parameters must be resolved firstly as in-situ water and oil velocities as well as cross-section parameters can be regarded as a function of liquid level,hw.Fairuzov [12] found that, while solving the steady state parameters of (4) for gas-liquid horizontal stratified flow, there may appear more than one solution, but, for oil-water stratified flow, there exists only one solution in the global scope. Therefore, in the current study, liquid level hw is obtained by varying the parameter from 0 to D.


#### 2.2.2. Transition of Stratified Flow for High Viscosity Ratio

In derivation of the neutral stability equation of stratified liquid-liquid flow, Brauner [13] found that the lower fluid of stratified flow is sometimes much faster, so the velocity of the upper phase may be neglected in neutral stability condition. Such an extreme situation may be of practical relevance, for example, in the transportation of high viscosity oil-water system. Ignoring the time and space variations in the upper fluid velocity yields
(7)Usw2(CrnUw-1)2<16π2DA~w3ρbA~w′(Δρg+σk2)(8)Crn=Uw(∂F/∂Uw)-(Aw/Aw′)(∂F/∂H)(∂F/∂Uw),
where Crn is the real neutral stable wave velocity.


### 2.3. Transition of Dispersed Flow

Dispersed flow exists on the condition that the continuous phase is turbulent and the turbulent kinetic energy in the whole region is sufficiently intense to break the dispersed droplets into drops smaller than the critical diameter:(9)dmax<dcr,ReC≥2100,1.82ReC-0.7<dcritD<0.1,
where D is the inner diameter of the pipe.


#### 2.3.1. Phase Inversion Point

It is important to introduce the phase inversion point. Factors that influence phase inversion vary, including two-phase density, viscosity, interfacial tension, and other physical properties and operation conditions. At the same time, temperature, oil-water system formation, and mixture container wettability also play an impact. The inversion point aims to provide a probable range of distinguishment between continuous and disperse phase where Decarre and Fabre [14] model is introduced:
(10)φinv=[1+(μoμw)1/6(ρoρw)1/6]-1,
where φinv stands for water content of inverting point.


#### 2.3.2. Maximum Diameter

In the turbulent field, the dispersed droplet is dominately affected by the turbulent inertia shear. Therefore, the maximum stable drop is the balance between turbulent kinetic energy and droplet interface energy [6]:
(11)Wecrit=τdmaxσ,(12)ρcu′22≈4σdmax.In Chen’s work [15], maximum drop size of dilute oil-water dispersed flow depends on the energy dissipation rate of the continuous phase and is directly proportional to (ε-0.4). Meanwhile, by introducing Barnea’s [16] idea of gas-liquid bubbly flow, the dispersed volume fraction is taken into consideration, which gives
(13)dmax=(1+aφs)(ρcρd)0.33(σρc)0.6ε-0.4,(14)ε=2fUm3D,
where ε is the energy dissipation rate, Um is the mixture velocity, f is the friction coefficient, and the subscripts c and d stand for continuous phase and dispersed phase, respectively.The higher the concentration of the dispersed phase, the more the energy needed by the continuous phase. Brauner [6] supposed that the incoming flow of the continuous phase should carry sufficient turbulent energy to disrupt the tendency to coalescence and to disperse the other phase. So the maximum droplet diameter can be expressed as
(15)dmax=2.22CH3/5(ρcUc2Dσ)-0.6×(εw1-εw)0.6[ρmρc(1-εw)f]-0.4D,
where CH is a constant, CH=O(1).  εw is the water holdup. ρm is the density of mixture.


#### 2.3.3. Critical Diameter

Due to the gravity or buoyancy effect, droplets will be separated from the mainstream, migrating to the wall. Considering the equilibrium of gravity and turbulent forces only, the droplet of critical dimension can be obtained through(16)dcbD=3ρcfUc28|ρc-ρd|Dg.By assuming the buoyancy force is strong enough to overcome interfacial tension and break a dispersed drop, the critical size can be expressed through the balance between interfacial tension and buoyancy forces:(17)dcσD=[0.4σ|ρc-ρd|gD2]0.5.The practical critical diameter is the minimum between the two conditions including oil-water two-fluid system and operation condition:(18)dcritD=Min(dcσD,dcbD).


### 2.4. Transition of Core-Annular Flow

Joseph et al. [17] used the hydrodynamic stability theory and concluded that a fully developed laminar-laminar core flow having the same density could stay stable when more viscous fluid is placed at the core and occupies most of the cross-section. Considering the stability of a viscous liquid flowing in the core, which was surrounded by a turbulent annulus, Bannwart [8] extended the previous model by using effective viscosity:
(19)μo,eff>μw,eff,μo>μw+0.0005ρwUswD,
where μeff is the effective viscosity which is the sum of the absolute molecular viscosity with the turbulent viscosity. To obtain the predicted transition boundary, it is necessary to find the critical Usw. Through a two-fluid model, Brauner [13] proposed an expression for calculating the core diameter, Dc, as a function of Usw and Uso, given the setup conditions and fluid properties. Considering a laminar core-turbulent annular flow, the expression that includes the influence of drop entrainment is given as
(20)DcD={+ϕ+G(ϕ)F(ϕ)[μoμwρwρb(νwνb)2G(ϕ)1.8F(ϕ)]1/2)-1ϕ×(χϕ[μoμwρwρb(νwνb)2G(ϕ)1.8F(ϕ)]1/2ccccccccccccz2[μoμwρwρb(νwνb)2G(ϕ)1.8F(ϕ)]1/2+ϕ+G(ϕ)F(ϕ)[μoμwρwρb(νwνb)2G(ϕ)1.8F(ϕ)]1/2)-1}1/2,
where ϕ=Uso/Usw and χ is the Lockhart-Martinelli parameter which is fully based on superficial variables:
(21)F(ϕ)=αwϕ-(1-αw)ϕ(αw+αc-1),(22)G(ϕ)=αc-(1-αc)ϕαw+αc-1,
where αw and αc stand for water volume fraction in the annular and oil volume fraction in the core. Once the oil superficial velocity Uso and Ho=Ho are fixed, the critical water velocity corresponding to the transition from core-annular flow is obtained.


## 2.1. Flow Pattern Classification

Although large amounts of experiments on oil-water two-phase flow in horizontal pipe have been conducted by scholars, there still lacks overall consensus about how the flow patterns could be classified. This may be due to the fact that experimental conditions differ and ways to observe flow regimes are not identical. In spite of that, the common regimes without drop entrainment are mainly considered here. Combine Taitel and Dukler’s gas-liquid flow pattern classification method with oil-water pipe flow characteristics, the flow patterns of low viscosity oil-water two-phase flow are classified as stratified flow, w/o dispersed flow, o/w dispersed flow, and intermittent flow.


## 2.2. Transition of Stratified Flow

### 2.2.1. Transition of Stratified Flow for Low Viscosity Ratio

Recent researches of stratified flow stability are all based on momentum and continuity equation of double fluid model by introducing the influence of small disturbance on the interface. It is common that ways of treating the stability equations could be classified into two categories: ZNS (zero neutral stability)/ZRC (zero real characteristics) criterion and Kelvin-Helmholz (IKH and VKH) instability theory. Brauner and Moalem Maron [11] used ZNS/ZRC criterion to predict the transition criterion of smooth stratified flow pattern to other flow patterns. ZNS line is set as the boundary of smooth stratified and stratified wavy and ZRC line as the upper limit of smooth stratified. Trallero [4] argued that the VKH equation determines the transition between stable stratified flow and either unstable stratified with wavy interface or some other flow patterns and may be used to predict onset of entrainment, which is suitable for low viscosity oil-water systems. The equations are as follows: (1)(CV-CIV)2+ρwρoρHwHo(Uw-Uo)2-(ρw-ρo)gcosβAρSi-σAρSik2︸IKH+ρf(Uw-Uo)2CSAρ(1Aw-1Ao)︷VKH<0(2)Jμ+JU-Jg-Jσ︸IKH+Js︷VKH<0(3)CV=∂F/∂Hw∂F/∂Uso-∂F/∂Usw(4)F=-τwSwAw+τoSoAo+τiSi(1Aw+1Ao)=0(5)CIV=ρwUwHw+ρoUoHoρwHw+ρoHo(6)ρ=A(ρwAw+ρoAo),where ρw and ρo are the density of water and oil, respectively; Hw and Ho are the holdup of water and oil, respectively; Usw and Uso are superficial velocity of water and oil, respectively; Uw and Uo are real mean velocity of water and oil, respectively; τw, τo, and τi are water, oil, and interfacial shear stresses; Sw, So, and Si are water, oil, and interface wet perimeter; Aw, Ao, and A are cross-sectional area of water, oil, and the pipe; CV and CIV are interface wave velocity and critical interface wave velocity, and k is the wave number.The first item on the left of (1) is the viscous instability, depending on τw, τo, and τi; the second item is the velocity instability, related to velocity difference of the two phases; the third item is the gravity stability; the fourth item is interface tension stability, related to interface tension and interface wave number; the fifth is a correction term which stands for the influence of other interface instable factors. CS is the correction coefficient.In order to apply the stratified flow criterion, the in-situ water and oil velocities as well as cross-section parameters must be resolved firstly as in-situ water and oil velocities as well as cross-section parameters can be regarded as a function of liquid level,hw.Fairuzov [12] found that, while solving the steady state parameters of (4) for gas-liquid horizontal stratified flow, there may appear more than one solution, but, for oil-water stratified flow, there exists only one solution in the global scope. Therefore, in the current study, liquid level hw is obtained by varying the parameter from 0 to D.


### 2.2.2. Transition of Stratified Flow for High Viscosity Ratio

In derivation of the neutral stability equation of stratified liquid-liquid flow, Brauner [13] found that the lower fluid of stratified flow is sometimes much faster, so the velocity of the upper phase may be neglected in neutral stability condition. Such an extreme situation may be of practical relevance, for example, in the transportation of high viscosity oil-water system. Ignoring the time and space variations in the upper fluid velocity yields
(7)Usw2(CrnUw-1)2<16π2DA~w3ρbA~w′(Δρg+σk2)(8)Crn=Uw(∂F/∂Uw)-(Aw/Aw′)(∂F/∂H)(∂F/∂Uw),
where Crn is the real neutral stable wave velocity.


## 2.2.1. Transition of Stratified Flow for Low Viscosity Ratio

Recent researches of stratified flow stability are all based on momentum and continuity equation of double fluid model by introducing the influence of small disturbance on the interface. It is common that ways of treating the stability equations could be classified into two categories: ZNS (zero neutral stability)/ZRC (zero real characteristics) criterion and Kelvin-Helmholz (IKH and VKH) instability theory. Brauner and Moalem Maron [11] used ZNS/ZRC criterion to predict the transition criterion of smooth stratified flow pattern to other flow patterns. ZNS line is set as the boundary of smooth stratified and stratified wavy and ZRC line as the upper limit of smooth stratified. Trallero [4] argued that the VKH equation determines the transition between stable stratified flow and either unstable stratified with wavy interface or some other flow patterns and may be used to predict onset of entrainment, which is suitable for low viscosity oil-water systems. The equations are as follows: (1)(CV-CIV)2+ρwρoρHwHo(Uw-Uo)2-(ρw-ρo)gcosβAρSi-σAρSik2︸IKH+ρf(Uw-Uo)2CSAρ(1Aw-1Ao)︷VKH<0(2)Jμ+JU-Jg-Jσ︸IKH+Js︷VKH<0(3)CV=∂F/∂Hw∂F/∂Uso-∂F/∂Usw(4)F=-τwSwAw+τoSoAo+τiSi(1Aw+1Ao)=0(5)CIV=ρwUwHw+ρoUoHoρwHw+ρoHo(6)ρ=A(ρwAw+ρoAo),where ρw and ρo are the density of water and oil, respectively; Hw and Ho are the holdup of water and oil, respectively; Usw and Uso are superficial velocity of water and oil, respectively; Uw and Uo are real mean velocity of water and oil, respectively; τw, τo, and τi are water, oil, and interfacial shear stresses; Sw, So, and Si are water, oil, and interface wet perimeter; Aw, Ao, and A are cross-sectional area of water, oil, and the pipe; CV and CIV are interface wave velocity and critical interface wave velocity, and k is the wave number.The first item on the left of (1) is the viscous instability, depending on τw, τo, and τi; the second item is the velocity instability, related to velocity difference of the two phases; the third item is the gravity stability; the fourth item is interface tension stability, related to interface tension and interface wave number; the fifth is a correction term which stands for the influence of other interface instable factors. CS is the correction coefficient.In order to apply the stratified flow criterion, the in-situ water and oil velocities as well as cross-section parameters must be resolved firstly as in-situ water and oil velocities as well as cross-section parameters can be regarded as a function of liquid level,hw.Fairuzov [12] found that, while solving the steady state parameters of (4) for gas-liquid horizontal stratified flow, there may appear more than one solution, but, for oil-water stratified flow, there exists only one solution in the global scope. Therefore, in the current study, liquid level hw is obtained by varying the parameter from 0 to D.


## 2.2.2. Transition of Stratified Flow for High Viscosity Ratio

In derivation of the neutral stability equation of stratified liquid-liquid flow, Brauner [13] found that the lower fluid of stratified flow is sometimes much faster, so the velocity of the upper phase may be neglected in neutral stability condition. Such an extreme situation may be of practical relevance, for example, in the transportation of high viscosity oil-water system. Ignoring the time and space variations in the upper fluid velocity yields
(7)Usw2(CrnUw-1)2<16π2DA~w3ρbA~w′(Δρg+σk2)(8)Crn=Uw(∂F/∂Uw)-(Aw/Aw′)(∂F/∂H)(∂F/∂Uw),
where Crn is the real neutral stable wave velocity.


## 2.3. Transition of Dispersed Flow

Dispersed flow exists on the condition that the continuous phase is turbulent and the turbulent kinetic energy in the whole region is sufficiently intense to break the dispersed droplets into drops smaller than the critical diameter:(9)dmax<dcr,ReC≥2100,1.82ReC-0.7<dcritD<0.1,
where D is the inner diameter of the pipe.


### 2.3.1. Phase Inversion Point

It is important to introduce the phase inversion point. Factors that influence phase inversion vary, including two-phase density, viscosity, interfacial tension, and other physical properties and operation conditions. At the same time, temperature, oil-water system formation, and mixture container wettability also play an impact. The inversion point aims to provide a probable range of distinguishment between continuous and disperse phase where Decarre and Fabre [14] model is introduced:
(10)φinv=[1+(μoμw)1/6(ρoρw)1/6]-1,
where φinv stands for water content of inverting point.


### 2.3.2. Maximum Diameter

In the turbulent field, the dispersed droplet is dominately affected by the turbulent inertia shear. Therefore, the maximum stable drop is the balance between turbulent kinetic energy and droplet interface energy [6]:
(11)Wecrit=τdmaxσ,(12)ρcu′22≈4σdmax.In Chen’s work [15], maximum drop size of dilute oil-water dispersed flow depends on the energy dissipation rate of the continuous phase and is directly proportional to (ε-0.4). Meanwhile, by introducing Barnea’s [16] idea of gas-liquid bubbly flow, the dispersed volume fraction is taken into consideration, which gives
(13)dmax=(1+aφs)(ρcρd)0.33(σρc)0.6ε-0.4,(14)ε=2fUm3D,
where ε is the energy dissipation rate, Um is the mixture velocity, f is the friction coefficient, and the subscripts c and d stand for continuous phase and dispersed phase, respectively.The higher the concentration of the dispersed phase, the more the energy needed by the continuous phase. Brauner [6] supposed that the incoming flow of the continuous phase should carry sufficient turbulent energy to disrupt the tendency to coalescence and to disperse the other phase. So the maximum droplet diameter can be expressed as
(15)dmax=2.22CH3/5(ρcUc2Dσ)-0.6×(εw1-εw)0.6[ρmρc(1-εw)f]-0.4D,
where CH is a constant, CH=O(1).  εw is the water holdup. ρm is the density of mixture.


### 2.3.3. Critical Diameter

Due to the gravity or buoyancy effect, droplets will be separated from the mainstream, migrating to the wall. Considering the equilibrium of gravity and turbulent forces only, the droplet of critical dimension can be obtained through(16)dcbD=3ρcfUc28|ρc-ρd|Dg.By assuming the buoyancy force is strong enough to overcome interfacial tension and break a dispersed drop, the critical size can be expressed through the balance between interfacial tension and buoyancy forces:(17)dcσD=[0.4σ|ρc-ρd|gD2]0.5.The practical critical diameter is the minimum between the two conditions including oil-water two-fluid system and operation condition:(18)dcritD=Min(dcσD,dcbD).


## 2.3.1. Phase Inversion Point

It is important to introduce the phase inversion point. Factors that influence phase inversion vary, including two-phase density, viscosity, interfacial tension, and other physical properties and operation conditions. At the same time, temperature, oil-water system formation, and mixture container wettability also play an impact. The inversion point aims to provide a probable range of distinguishment between continuous and disperse phase where Decarre and Fabre [14] model is introduced:
(10)φinv=[1+(μoμw)1/6(ρoρw)1/6]-1,
where φinv stands for water content of inverting point.


## 2.3.2. Maximum Diameter

In the turbulent field, the dispersed droplet is dominately affected by the turbulent inertia shear. Therefore, the maximum stable drop is the balance between turbulent kinetic energy and droplet interface energy [6]:
(11)Wecrit=τdmaxσ,(12)ρcu′22≈4σdmax.In Chen’s work [15], maximum drop size of dilute oil-water dispersed flow depends on the energy dissipation rate of the continuous phase and is directly proportional to (ε-0.4). Meanwhile, by introducing Barnea’s [16] idea of gas-liquid bubbly flow, the dispersed volume fraction is taken into consideration, which gives
(13)dmax=(1+aφs)(ρcρd)0.33(σρc)0.6ε-0.4,(14)ε=2fUm3D,
where ε is the energy dissipation rate, Um is the mixture velocity, f is the friction coefficient, and the subscripts c and d stand for continuous phase and dispersed phase, respectively.The higher the concentration of the dispersed phase, the more the energy needed by the continuous phase. Brauner [6] supposed that the incoming flow of the continuous phase should carry sufficient turbulent energy to disrupt the tendency to coalescence and to disperse the other phase. So the maximum droplet diameter can be expressed as
(15)dmax=2.22CH3/5(ρcUc2Dσ)-0.6×(εw1-εw)0.6[ρmρc(1-εw)f]-0.4D,
where CH is a constant, CH=O(1).  εw is the water holdup. ρm is the density of mixture.


## 2.3.3. Critical Diameter

Due to the gravity or buoyancy effect, droplets will be separated from the mainstream, migrating to the wall. Considering the equilibrium of gravity and turbulent forces only, the droplet of critical dimension can be obtained through(16)dcbD=3ρcfUc28|ρc-ρd|Dg.By assuming the buoyancy force is strong enough to overcome interfacial tension and break a dispersed drop, the critical size can be expressed through the balance between interfacial tension and buoyancy forces:(17)dcσD=[0.4σ|ρc-ρd|gD2]0.5.The practical critical diameter is the minimum between the two conditions including oil-water two-fluid system and operation condition:(18)dcritD=Min(dcσD,dcbD).


## 2.4. Transition of Core-Annular Flow

Joseph et al. [17] used the hydrodynamic stability theory and concluded that a fully developed laminar-laminar core flow having the same density could stay stable when more viscous fluid is placed at the core and occupies most of the cross-section. Considering the stability of a viscous liquid flowing in the core, which was surrounded by a turbulent annulus, Bannwart [8] extended the previous model by using effective viscosity:
(19)μo,eff>μw,eff,μo>μw+0.0005ρwUswD,
where μeff is the effective viscosity which is the sum of the absolute molecular viscosity with the turbulent viscosity. To obtain the predicted transition boundary, it is necessary to find the critical Usw. Through a two-fluid model, Brauner [13] proposed an expression for calculating the core diameter, Dc, as a function of Usw and Uso, given the setup conditions and fluid properties. Considering a laminar core-turbulent annular flow, the expression that includes the influence of drop entrainment is given as
(20)DcD={+ϕ+G(ϕ)F(ϕ)[μoμwρwρb(νwνb)2G(ϕ)1.8F(ϕ)]1/2)-1ϕ×(χϕ[μoμwρwρb(νwνb)2G(ϕ)1.8F(ϕ)]1/2ccccccccccccz2[μoμwρwρb(νwνb)2G(ϕ)1.8F(ϕ)]1/2+ϕ+G(ϕ)F(ϕ)[μoμwρwρb(νwνb)2G(ϕ)1.8F(ϕ)]1/2)-1}1/2,
where ϕ=Uso/Usw and χ is the Lockhart-Martinelli parameter which is fully based on superficial variables:
(21)F(ϕ)=αwϕ-(1-αw)ϕ(αw+αc-1),(22)G(ϕ)=αc-(1-αc)ϕαw+αc-1,
where αw and αc stand for water volume fraction in the annular and oil volume fraction in the core. Once the oil superficial velocity Uso and Ho=Ho are fixed, the critical water velocity corresponding to the transition from core-annular flow is obtained.


## 3. Results and Discussion

### 3.1. Model Diagram

The oil-water flow pattern transition prediction is based on sound results of fluid properties, certain pipelines, and operation parameters. Thus the criterion can be integrated as given in Figure1.Figure 1
Diagram of flow pattern transition program.Figure2(a) indicates how the diameter affects the stability of the stratified flow and dispersed flow. The diameters are 1, 2, and 3 inches, respectively. Viscosity ratio of oil to water (μ~) is 29.6 and interface tension is 0.036 N/m. Prediction indicates that the scope of stratified flow becomes broader when the diameter increases, but the scope of o/w and w/o dispersed flow becomes narrower. It can be concluded that, for oil-water two-phase flows, stratified flow pattern is more likely to form in a large-diameter pipe, as the effect of gravity stable term (1) is more obvious, which, to a certain extent, hinders the growth of interface wave and stabilizes the two-phase flow while, in smaller pipes, dispersed oil-water flow is relatively common, due to the fact that flow energy dissipation rate and turbulent kinetic energy are large though mixing velocity is unchanged. In addition, under the same diameter and internal phase concentration, the velocity needed to maintain a complete o/w dispersed flow pattern is bigger than that of a complete w/o dispersed flow pattern.Prediction of flow pattern transition.
(a)
Pipe diameter
(b)
Oil phase viscosityFigure2(b) shows the influence of viscosity on stability of stratified and dispersed flow. The viscosities are 10, 50, and 80 mPa·s, respectively, while pipe diameter is 2 inch, and the interfacial intension is 0.036 N/m. Prediction indicates that the scope of stratified flow becomes narrower with larger oil viscosity while the scope of w/o dispersed flow becomes broader, and the sensitivity of o/w dispersed flow boundary is weaker.The more viscous the oil, the stronger the influence of viscous instability and the smaller the superficial velocity needed to maintain the two phases stratified. From (9) and (12), dispersed viscosity plays no impact on maximum droplet diameter and critical droplet diameter. Therefore, for o/w dispersed flow, the model does not provide the trend of boundary changes with oil viscosity. This is due to the present model neglecting the effect of dispersed drop viscosity.


### 3.2. Discussions

Models are validated against available experimental data. Two cases of low viscosity oil (viscosity ratio of oil/aqueous phase <100) and two higher ones (viscosity ratio >100) are selected and compared. Comparisons with Oglesby [1] and Trallero [4] results are chosen as a representation of low viscosity.Figure3(a) is a comparison between predicted results and Trallero [4] experimental data (viscosity ratio 29.6), which indicates a good agreement between the predicted results of o/w and w/o transition boundaries and the experimental data. The scope of stratified flow with mixed interface (stratified wavy flow) in experiment is very large; the proper explanation is that the low viscosities of both phases make it easier for drop entrainment to take place. Actually this flow pattern is a transitional flow pattern from stratified flow to dispersed flow. The predicted transition boundary of stratified flow includes some stratified flows with mixed interface.Comparison of experimental data (symbols) and prediction (lines) of oil-water two phase flow pattern transition.
(a)
(b)Figure3(b) compares the predicted boundaries with the experimental data from Oglesby [1] (viscosity ratio 62.8). The prediction of w/o boundary fits well with experimental data. The prediction of o/w dispersed flow boundary is narrower than experimental results. This is partly because the sensitivity of the o/w dispersed flow criterion is weak (discussed above). Another explanation is that oil and water were mixed before entering the experimental system. The premix process makes o/w dispersed flow scope larger and stratified flow scope narrower. In addition, the predicted scope of stratified flow is narrow, and this is because the viscous instability becomes stronger.Figure4(a) is a comparison between predicted results and Yao and Gong [18] experimental data (viscosity ratio 325.4), while Figure 4(b) compares the predicted boundaries and experimental data from Wang et al. [19] (viscosity ratio 620). Core-annular, stratified, w/o dispersed, and intermittent flow are considered.Comparison between prediction (lines) and experimental data (symbols) of oil-water flow.
(a)
(b)Both figures show that the stratified flow transition criterion fail to predict the accurate region. This may be caused by drop entrainment in viscous oil. And pure single phase in certain region of the pipe cross-section rarely exists. It should be noted that (7) does not take into account the drop entrainment in both phases, whose regime actually is Dw/o and Do/w. In addition, the core-annular/Do/w transition boundary agrees well with experimental data. This boundary is the lowest oil holdup requirement, beyond which the oil core will be broken and finally dispersed in the water. As water drop is easily entrained into the viscous oil phase, the flow regime of Dw/o becomes dominant when phase inversion occurs. Phase inversion point can be perfectly regarded as a criterion for Dw/o transition boundary.


## 3.1. Model Diagram

The oil-water flow pattern transition prediction is based on sound results of fluid properties, certain pipelines, and operation parameters. Thus the criterion can be integrated as given in Figure1.Figure 1
Diagram of flow pattern transition program.Figure2(a) indicates how the diameter affects the stability of the stratified flow and dispersed flow. The diameters are 1, 2, and 3 inches, respectively. Viscosity ratio of oil to water (μ~) is 29.6 and interface tension is 0.036 N/m. Prediction indicates that the scope of stratified flow becomes broader when the diameter increases, but the scope of o/w and w/o dispersed flow becomes narrower. It can be concluded that, for oil-water two-phase flows, stratified flow pattern is more likely to form in a large-diameter pipe, as the effect of gravity stable term (1) is more obvious, which, to a certain extent, hinders the growth of interface wave and stabilizes the two-phase flow while, in smaller pipes, dispersed oil-water flow is relatively common, due to the fact that flow energy dissipation rate and turbulent kinetic energy are large though mixing velocity is unchanged. In addition, under the same diameter and internal phase concentration, the velocity needed to maintain a complete o/w dispersed flow pattern is bigger than that of a complete w/o dispersed flow pattern.Prediction of flow pattern transition.
(a)
Pipe diameter
(b)
Oil phase viscosityFigure2(b) shows the influence of viscosity on stability of stratified and dispersed flow. The viscosities are 10, 50, and 80 mPa·s, respectively, while pipe diameter is 2 inch, and the interfacial intension is 0.036 N/m. Prediction indicates that the scope of stratified flow becomes narrower with larger oil viscosity while the scope of w/o dispersed flow becomes broader, and the sensitivity of o/w dispersed flow boundary is weaker.The more viscous the oil, the stronger the influence of viscous instability and the smaller the superficial velocity needed to maintain the two phases stratified. From (9) and (12), dispersed viscosity plays no impact on maximum droplet diameter and critical droplet diameter. Therefore, for o/w dispersed flow, the model does not provide the trend of boundary changes with oil viscosity. This is due to the present model neglecting the effect of dispersed drop viscosity.


## 3.2. Discussions

Models are validated against available experimental data. Two cases of low viscosity oil (viscosity ratio of oil/aqueous phase <100) and two higher ones (viscosity ratio >100) are selected and compared. Comparisons with Oglesby [1] and Trallero [4] results are chosen as a representation of low viscosity.Figure3(a) is a comparison between predicted results and Trallero [4] experimental data (viscosity ratio 29.6), which indicates a good agreement between the predicted results of o/w and w/o transition boundaries and the experimental data. The scope of stratified flow with mixed interface (stratified wavy flow) in experiment is very large; the proper explanation is that the low viscosities of both phases make it easier for drop entrainment to take place. Actually this flow pattern is a transitional flow pattern from stratified flow to dispersed flow. The predicted transition boundary of stratified flow includes some stratified flows with mixed interface.Comparison of experimental data (symbols) and prediction (lines) of oil-water two phase flow pattern transition.
(a)
(b)Figure3(b) compares the predicted boundaries with the experimental data from Oglesby [1] (viscosity ratio 62.8). The prediction of w/o boundary fits well with experimental data. The prediction of o/w dispersed flow boundary is narrower than experimental results. This is partly because the sensitivity of the o/w dispersed flow criterion is weak (discussed above). Another explanation is that oil and water were mixed before entering the experimental system. The premix process makes o/w dispersed flow scope larger and stratified flow scope narrower. In addition, the predicted scope of stratified flow is narrow, and this is because the viscous instability becomes stronger.Figure4(a) is a comparison between predicted results and Yao and Gong [18] experimental data (viscosity ratio 325.4), while Figure 4(b) compares the predicted boundaries and experimental data from Wang et al. [19] (viscosity ratio 620). Core-annular, stratified, w/o dispersed, and intermittent flow are considered.Comparison between prediction (lines) and experimental data (symbols) of oil-water flow.
(a)
(b)Both figures show that the stratified flow transition criterion fail to predict the accurate region. This may be caused by drop entrainment in viscous oil. And pure single phase in certain region of the pipe cross-section rarely exists. It should be noted that (7) does not take into account the drop entrainment in both phases, whose regime actually is Dw/o and Do/w. In addition, the core-annular/Do/w transition boundary agrees well with experimental data. This boundary is the lowest oil holdup requirement, beyond which the oil core will be broken and finally dispersed in the water. As water drop is easily entrained into the viscous oil phase, the flow regime of Dw/o becomes dominant when phase inversion occurs. Phase inversion point can be perfectly regarded as a criterion for Dw/o transition boundary.


## 4. Conclusions

In current work, theoretical models concerning systems with both low and high viscosity ratio are introduced for criterion of flow pattern transition of oil-water horizontal flow. The main flow regimes considered include water stratified flow, core-annular flow, w/o dispersed flow, o/w dispersed flow, and intermittent flow. A systematic work is carried out for flow pattern transition with oil viscosity varying in a wide range.It is concluded that stability of the oil-water stratified flow in horizontal pipe is strongly related to oil viscosity, gravity, and interfacial intension. For viscous oils, the influence of shear stress becomes much more obvious and can be characterised by ignoring the velocity of the viscous phase while, in dispersed flow, as the viscosity increases, oil droplets become harder to breakup, meaning that the ability of droplets to recover deformation becomes stronger. It is more difficult for o/w dispersed flow to be formed in viscosity oil. Additionally, core-annular flow tends to occur in a pipe where the two fluids have much different viscosities but relatively close densities. Drop entrainment occurs easily in core-annular flow and must be taken into consideration in the transition criterion.



---

*Source: 291217-2013-12-10.xml*