DUNS Y ROSS

La correlación de Duns y Ros es el resultado de un extenso estudio de laboratorio

en el que se midieron el colgamiento de líquidos y los gradientes de presión. Se
desarrollaron correlaciones para la velocidad de resbalamiento (a partir de la cual
se puede calcular el colgamiento) y el factor de fricción para cada uno de los tres
regímenes de flujo (Fig. 3.1). Los regímenes de flujo se definen como funciones de
las cantidades adimensionales Ngv, NLv, L1, L2, Ls, Lm y Nd donde:
Ls = 50 + 36*NLv
Lm =75 +84*NLv
L1 y L2 son funciones de Nd como se muestra en la Fig. 3.11. Duns y Ros
presentaron correlaciones para una velocidad de resbalamiento adimensional a
partir de la cual la velocidad de resbalamiento real y el colgamiento del líquido
pueden calcularse utilizando las siguientes ecuaciones:
1/4
ρL
S=v s ( ) σLg

Donde:
v sg v sL
v s=v g−v L = −
1−H L H L

Solución de la ecuación anterior para el colgamiento

2 1/ 2
v s−v m+ [ ( v m−v s ) +4 v s v sL ]
H L=
2 vs

El procedimiento para calcular el gradiente de presión debido al cambio de

elevación es:
1. Calcular la velocidad de resbalamiento adimensional S, utilizando la
correlación apropiada. La correlación para S es diferente para cada régimen
de flujo.
2. Resolver la ecuación Num. (Ecuación de S) para la velocidad de
resbalamiento, vs.
3. Calcular el colgamiento del líquido, HL usando la ecuación Num. (Ecuación
de HL).
4. Calcular la densidad de la mezcla, ρ s=ρ L H L + ρg ( 1−H L ).
5. Calcular el gradiente de presión debido al cambio de elevación,

( dZdp ) = gg ρ
el c
s
FLUJO BURBUJA
Limites: 0 ≤ N gv ≤ L1+ L2 N Lv

Densidad para dos fases

2
N gv
S=F 1 + F2 N Lv+ F 3 '
( 1+ N Lv )
Donde F1, F2, F3, y F4 se obtienen de la figura ___, Los cuales son
funciones del número de viscosidad del líquido, N L. F '3 puede ser obtenido
de la siguiente forma:

' F4
F 3 =F 3−
Nd

Factor de Fricción

dP f ρ V V
( )
dZ f
= m L sL m
2 gc d

De datos experimentales Duns y Ros llegaron a la siguiente ecuación para

f m:

f2
f m=f 1
f3
f 1 es obtenido del diagrama de Moody como una función del número de
Reynolds del líquido:
ρL V sL d
Nℜ = L
μL

El factor f 2es una corrección para la relación de gas-líquido en situ, y es

V sg 2/ 3
dada en la figura __ como una función del grupo f 1 ( ) N . El factor f 3es
V sL d
considerado por Duns y Ros como un factor de corrección adicional para
ambas relaciones viscosidad y gas-líquido en situ y se convierte en
importante para viscosidades mayores que 50 centistokes
aproximadamente. Este factor es calculado de la siguiente forma:
0.5
V sg
f 3=1+ f 1 ( 50 V sL )
Término de aceleración
 El término de aceleración fue considerado insignificante en el régimen de
flujo burbuja.

FLUJO SLUG
Límites: L1 + L2 N Lv ≤ N gv ≤ L s

Densidad para dos fases

N gv 0.982+ F 6'
S=(1+ F 5)
(1+ F7∗N Lv )2

F5, F6 y F7 se obtienen de la figura ___ como una función del número de

viscosidad del líquido, N L, y

F 6' =0.029 N d + F 6

Factor de fricción
El gradiente de presión debido a la fricción en este régimen de flujo es
calculado usando el procedimiento dado para el flujo burbuja.
Término de aceleración
El término de aceleración fue considerado insignificante en el régimen de
flujo slug.

FLUJO NIEBLA
Limites: N gv ≥ Lm
Densidad para dos fases
Duns y Ros asumieron que con la alta velocidad del gas en el régimen de
flujo niebla la velocidad de resbalamiento fue cero. Por lo tanto, la densidad
de la mezcla es la densidad sin resbalamiento y es calculado de la siguiente
forma:
V sL V sg
ρn =ρL λL + ρg λ g= ρL + ρg
Vm Vm

Factor de fricción
En esta región, el término de fricción es basado solo en la fase gas.
Entones:
dP f ρ V2
( ) = g sg
dZ f 2 gc d

Ya que no hay resbalamiento, el factor de fricción es dado en el diagrama

de Moody, pero como una función del número de Reynolds del gas:
ρ g V sg d
N ℜ=
μg

Duns y Ros notaron que la rugosidad de la pared para el flujo niebla es

afectado por la película del líquido en la pared del tubo. Las ondas de la
película de la pared causan un arrastre en el gas. Este proceso es
gobernado por una forma de número de Weber:

ρ g V 2sg e
N we =
σL

Y es un poco afectado por la viscosidad del líquido. Esta influencia fue

explicada haciendo N we una función de un número adimensional que
contiene la viscosidad del líquido,

μ 2L
N μ=
ρg σ L e

La relación funcional es mostrada en la fig _____ donde las coordenadas

son N we vs N we N μ
El valor de la rugosidad puede ser muy pequeña pero ϵ /d nunca se
convertirá más pequeña que el valor para la tubería en sí. En la zona de
transición para el flujo slug, ϵ /d puede aproximarse a 0.5. entre estos
límites, ϵ /d puede obtenerse de la ecuaciones siguientes las cuales son
desarrolladas de la Fig ___.

ε 0.0749 σ L
N we N μ ≤ 0.005 ; =
d ρg V 2sg d

ε 0.3717 σ L 0.302
N we N μ > 0.005; = 2 ( N we N μ )
d ρ g V sg d

Donde:
σ L= Tensión interfacial gas-líquido, dinas/cm

ρ g= Densidad del gas, lbm/ft3

V sg = velocidad superficial del gas, ft/sec

d = Diámetro del tubo, ft

 Valores de f para el régimen de flujo niebla puede ser calculada para ϵ /d>
0.05 de la siguiente forma:

1 ϵ 1.73
f=
[[ 4 log 10 (
0.27 ε
d
)
2
+0.067∗(
] d
) ∗4
]
Como la altura de la onda en las paredes del tuve incrementa, el área actual
a través del cual el gas puede fluir es disminuido, desde el diámetro abierto
hasta que el flujo de gas es d−ϵ. Duns y Ros sugestionan que la predicción
V sg d 2
de fricción baja puede ser referida por sustitución de d−ϵ por d y
por
(d −ε )2
V sg en todo el cálculo del gradiente de presión. En este caso la
determinación de la rugosidad, ϵ , es iterativo.
Término de la aceleración.

dP ρ n V sg V m dp
( )
dZ acc
=
gc P dZ

Si definimos
ρn V sg V m
Ek =
gc P

Entonces el gradiente de presión total puede ser calculado como:

dp dp
dp
=
( dZ ) +( )
dZ el f
dZ 1−Ek

FLUJO DE TRANSICIÓN
Límites: Ls < N gv < Lm

En la zona de transición entre el flujo slug y el flujo niebla, Duns y

Ros sugestionan interpolación lineal entre las fronteras del régimen
de flujos Ls y Lm, para obtener el gradiente de presión. Este valor
entre N gv falla entre Ls y Lm, el gradiente de presión puede ser
calculado usando ambos, la correlación de flujo slug y la correlación
de flujo niebla. El gradiente de presión en la zona de transición es
entonces calculado como:
dp dp
dZ
=A
dZ ( ) slug
+B ( dpdZ ) niebla

Donde:
 Lm −N gv N −L
A= , B= gv s =1−A
Lm−Ls Lm−Ls

Mayor precisión fue reclama si el cálculo de la densidad del gas

usada en el gradiente de presión del flujo niebla fue modificada
como:

' ρg N gv
ρ g=
Lm

Donde ρ g es la densidad del gas calculada a las condiciones de

presión y temperaturas dadas. Esta modificación cuenta para
algunos líquidos arrastrados en el gas.
EJERCICIO

Dado:
ft
V sg =4.09
s
ε
=0.000 6
d
μo =18Cp

ft
V sL =2.65
s
dynas
σ L =30
cm
μ g=0.018 Cp

d=0.249 f t
N d =41.3 4

N Lv=6.0 2

P=720 psia
N gv=9.2 9

T =128 ° F
lbm
ρ L=56.6
ft 3
 lbm
ρ g=2.84
ft 3
N L=0.0 8

Despreciando la aceleración, calcule el gradiente de flujo en estas condiciones.

1. Calcular los límites del régimen de flujo y determinar el régimen de flujo.

De la Figura 1

L1=1.53 L2=0.88
L1 + L2∗Nde
Fig. 1 Número 6.02 )=6.83
régimen de (flujo
Lv =1.53+0.88∗

Ls =50+36∗N Lv =50+36∗( 60.2 )=266.72

Lm=75+ 84∗N Lv 0.75=397.8

Ya que L1 + L2∗N Lv < N gv ≤ Ls , el régimen de flujo es SLUG

2. Determinar la densidad de dos fases.

De la Figura 2 para N L=0.08
F 5=0.065 F6 =1.95 F 7=0.041

F ´ 6=0.029∗N d + F 6=0.029∗( 41.34 ) +1.95=3.15

N gv 0.982+ F ´ 6 ( 9.29 )0.982+3.15

S= ( 1+ F 5 )
( 1+ F 7∗N Lv )
2
=( 1+0.065 )
[ ( 1+0.041∗6.02 )2 ] =8.29

S σ 1 8.29 30 14 ft
V s=
1.938 ρ L ( )
∗ L 4= ∗
1.938 56.6
=3.64
sg( )
2 0.5
V s −V m + [ ( V m −V s ) + 4∗V s∗V sL ]
H L=
2∗V s
 2 0.5
3.64−6.741+ [ ( 6.74−3.64 ) + 4∗(3.64)∗¿ 2.65 ¿ ]
H L= =0.527
2∗( 3.64 )

3. Determinar la densidad de la mezcla

lbm
ρ s=ρ L∗H L + ρg ( 1−H L ) =56.6 ( 0.527 ) +2.84 ( 1−0.527 )=31.7
ft 3

4. Determinar el gradiente de fricción

1488∗ρ L∗V sL∗d 1488∗56.6∗2.65∗0.249
N ℜ= = =3087
μL 18
De la Figura 3, f 1=0.044

Fig. 3 Diagrama de Moody

2 2
V sg 4.09
( )
f 2∗
V sL
∗ N d 3 =0.044∗( )
2.65
∗( 41.34 ) 3 =0.81
 De la figura 4, f 2=0.87

Fig. 4 Corrección de fricción por burbuja

0.5 0.5
V sg 4.09
f 3=1+ f 1 (
50∗V sL ) =1+0.044 (
50∗2.65 ) =1.008

f 2 0.044∗0.87
f m=f 1 = =0.038
f3 1.008

dP f m∗ρ L∗V sL∗V m lb f

( )
dZ f
=
2∗g c∗d
=2.40 3
ft
5. Determinar el gradiente de presión total
lbf
( dPdZ ) =( dPdZ ) +( dPdZ ) =31.17+2.40=33.572 .40 ft
T el f
3

( dPdZ ) =0.233 psift

T