Cap5 Flujo Monofasico Tuberias
Cap5 Flujo Monofasico Tuberias
Cap5 Flujo Monofasico Tuberias
CAPITULO 5
ECUACIONES
FUNDAMENTALES PARA
EL FLUJO MONOFASICO
EN TUBERIAS
1
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 2
Contenido Pag.
5. 2 Ecuación de continuidad 7
5. 3 Ecuación de momentum 12
5. 4 Ecuación de la energía 15
5. 5 Flujo monofásico 20
2
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 3
Contenido Pag.
3
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 4
Estas leyes pueden ser aplicadas a un fluido en diferentes formas. Entre las más
comunes tenemos el enfoque de tratar al fluido como un sistema o un volumen de
control. Por definición, “un sistema es una cantidad de materia, de masa y de
identidad fijas, la cual está separada del medio que la rodea por medio de una
superficie. El sistema puede ser rígido, móvil o bien puede interactuar con su medio
ambiente. La superficie que separa al sistema de su medio ambiente puede ser real,
imaginaria, fija o móvil, con la condición de no permitir la transferencia de masa”. Por
otro lado, el volumen de control representa “un volumen en el espacio, el cual puede
estar fijo, moverse o relacionarse con su medio ambiente a través de su superficie, la
cual puede ser real o imaginaria, fija o deformable (móvil), con la condición de
permitir la transferencia de masa”.
4
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 5
Las leyes que gobiernan el flujo de fluidos a través de un volumen de control son
derivadas a partir de su aplicación a un sistema. Por ejemplo, el desarrollo de la
ecuación de continuidad se obtiene del hecho de considerar que la masa de un
sistema es constante, mientras que la ecuación de cantidad de movimiento se
obtiene de la aplicación de la ecuación del cambio de momentum aplicable a un
sistema o la segunda ley de movimiento de Newton, que establece el producto de la
masa de un sistema por su aceleración es igual a la suma de todas las fuerzas
actuantes sobre dicho sistema. De esta manera se desarrollan las expresiones de
las ecuaciones de continuidad y la ecuación de momentum lineal para un volumen
de control.
5
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 6
B = bm (5.1)
⎛ DB ⎞ d
= ∫vc ρbdV + ∫sc ρbV.ndA
⎜ Dt ⎟
⎝ ⎠ sist dt (5.2)
donde:
⎛ DB ⎞
⎜ Dt ⎟ , representa el cambio de cualquier propiedad extensiva B en el sistema.
⎝ ⎠sist
d
∫ ρbdV
dt vc , representa la tasa de cambio de B en el volumen de control.
6
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 7
Considere que la propiedad extensiva “B” del sistema está representada por su
masa “m”, esto es B = m. Luego, por definición de una propiedad específica b = 1.
Por el teorema del transporte de Reynolds se tiene:
⎛ Dm ⎞ = d
⎜ Dt ⎟ ∫vc ρdV + ∫sc ρV.ndA
⎝ ⎠sist dt (5.3)
7
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 8
⎛ Dm ⎞ = 0
⎜ Dt ⎟
⎝ ⎠sist (5.3a)
d
∫vc ρdV + ∫sc ρV.ndA = 0 (5.4)
dt
Las dimensiones en la ecuación (5.4) son dimensiones de masa (m) sobre tiempo
(t), comúnmente llamada relación de flujo de masa.
d
∫ ρdV = 0
dt vc (5.4a)
8
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 9
Luego
Lo que permite afirmar que “el flujo de masa neto a través de la superficie de control
es igual a cero”, lo que se puede interpretar como una igualdad entre el flujo de
masa que se le suministra o extrae al volumen de control. Bajo esta situación se le
considera al flujo como flujo estacionario o estable.
9
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 10
⎡ ⎛ r ⎞2 ⎤
V =V ⎢1 − ⎜ ⎟ ⎥
max ⎢ ⎜⎝ R ⎟⎠ ⎥ (5.6)
⎣ ⎦
Para tener el valor de una distribución de velocidad uniforme se hace uso del
concepto de la velocidad media o promedio en una sección transversal, la cual se
define como:
∫ V.ndA
V = A (5.7)
A
V =V /2
med max (5.7a)
( )
∫sc ρV.ndA = ρ salVsal Asal −( ρ entVent Aent ) (5.8)
10
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 11
d
ρ ∫ dV + ρ ∫sc V.ndA = 0
dt vc (5.9)
⎛ dV ⎞ + V.ndA = 0
⎜ dt ⎟ ∫
⎝ ⎠vc sc (5.9a)
Q = ∫ V.ndA
A (5.10b)
11
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 12
⎛ DmV ⎞ = ∑F
⎜ ⎟
⎝ Dt ⎠ sist (5.11)
d
∑F = ∫ ρVdV + ∫sc ρVV.ndA (5.13)
dt vc
La ecuación anterior establece que “la fuerza neta actuante sobre un volumen de
control es igual a la suma de la tasa de cambio de momento dentro del volumen de
control más el flujo neto de momento que atraviesa la superficie de control”. Las
12
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 13
fuerzas involucradas en esta ecuación están formadas por las fuerzas de cuerpo y
las fuerzas de superficies.
d
∑F = ∫ ρVdV + ∫sc ρVV.ndA (5.14)
dt vc
d
∫ ρVdV = 0
dt vc (5.15)
Luego,
13
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 14
∑F = 0 (5.17)
La ecuación (5.17) establece que “la fuerza neta actuante sobre un volumen de
control es igual a cero”.
dP
( P ) A − τdAw − ( P − dL) A = 0 (5.18)
dL
dP τ dAw
= (5.19)
dL A dL
Luego:
⎛ dP ⎞ 4τ w
⎜ ⎟ =
⎝ dL ⎠w D (5.20)
14
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 15
La ecuación (5.20) permite afirmar que para determinar las pérdidas de presión por
fricción en la pared de una tubería solo es necesario conocer la expresión del
esfuerzo de corte en la pared de la misma.
2 rτ τ τ
τ = w = w
D r R (5.21)
15
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 16
Qsist > 0
Qsist < 0
⎛ DE ⎞ = d
⎜ Dt ⎟ ∫ ρedV + ∫sc ρeV.ndA
⎝ ⎠sist dt vc (5.23)
Para la condición en la cual el contenido del sistema coincide con el contenido del
volumen de control se cumple que:
⎛•⎞ ⎛ • ⎞ d
⎜ Q ⎟ − ⎜W ⎟ = dt ∫vc ρedV + ∫sc ρeV.ndA
⎝ ⎠vc ⎝ ⎠vc (5.25)
En esta ecuación, la tasa de transferencia de calor Qvc representa todas las formas
de transferir energía en forma de calor, ya sea por radiación, convección o
conducción como consecuencia de la diferencia de la temperatura del fluido dentro
del volumen de control y la correspondiente a su medio ambiente. La transferencia
de calor se considera positiva (Qvc >0) cuando es suministrada al volumen de
control. Cuando la transferencia de calor se extrae del volumen de control (Qvc <0)
se considera negativa y cuando no existe se considera igual a cero (Qvc =0), en este
caso el proceso se considera adiabático.
16
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 17
t t +dt
La magnitud de este trabajo se determina por medio del producto escalar entre la
fuerza normal que produce el desplazamiento de la superficie de control y la
velocidad del fluido, para un elemento de fluido sobre la superficie de control, así:
δW = V • dF
σ σ (5.26)
Expresando la fuerza normal en función del esfuerzo normal y este a su vez por la
presión se tiene:
δW = V • σdA = PV • ηdA
σ (5.27)
δW = PυρV • ηdA
σ (5.28)
17
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 18
Luego, el trabajo neto como consecuencia del flujo de masa que atraviesa la
superficie de control es:
W = ∫ PυρV • ηdA
σ sc (5.29)
• • • •
W = Weje + Wτ + Wσ (5.30)
Considerando:
• • •
W s = Weje + Wτ (5.31)
⎛•⎞ ⎛ • ⎞ d
⎜ ⎟
Q − ⎜W s ⎟ = dt ∫vc ρedV + ∫sc (e + Pυ ) ρV.ndA (5.32)
⎝ ⎠vc ⎝ ⎠ vc
18
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 19
⎛ ⎛ V2 gZ ⎞⎟ ⎞⎟ ⎛ ⎛ V2 gZ ⎞⎟ ⎞⎟ ⎛ V2 gZ ⎞ ⎛ V 2 gZ ⎞
Q − W = ⎜ m⎜ u + + − ⎜ m⎜ u + + + ∑ m⎜⎜ h + 2 g + ⎟⎟ − ∑ m⎜⎜ h + + ⎟
vc vc ⎜ ⎜
⎝ ⎝
2g
c
g ⎟⎟
c ⎠⎠ f ⎝ ⎝
⎜ ⎜ 2g
c
g ⎟⎟
c ⎠ ⎠i s ⎝ c gc ⎠ e ⎝ 2 g c g c ⎟⎠
(5.33a)
Q −W =
vc vc ∑ m ( e + Pυ ) − ∑ m ( e + Pυ )
s e
(5.34)
⎛ V 2 gZ ⎞ ⎛ V 2 gZ ⎞
Q − W = ∑ m⎜⎜ h + + ⎟⎟ − ∑ m⎜⎜ h + + ⎟⎟ (5.35)
vc vc s ⎝ 2 g c g c ⎠ e ⎝ 2 g c g c ⎠
⎛ P V 2 gZ ⎞ ⎛ P V 2 gZ ⎞
Q − W = ∑ m⎜ u + +
⎜ + ⎟ − ∑ m⎜ u + ρ + 2 g + g ⎟⎟
⎟ ⎜ (5.36)
vc vc s ⎝ ρ 2 g c g c ⎠ e ⎝ c c ⎠
V2
Efecto de la energía cinética = ∫sc ρV.ndA (5.37)
2
19
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 20
V2
∫sc ρV.ndA
α= 2
V2 (5.38)
m
2
Para flujo laminar en una tubería descrito por un perfil de velocidad parabólico:
⎡ ⎛ r ⎞2 ⎤
V =V ⎢1 − ⎜ ⎟ ⎥ ⇒ α = 2
max ⎢ ⎜⎝ R ⎟⎠ ⎥ (5.39)
⎣ ⎦
Mientras que para flujo turbulento, con un perfil de velocidad definido por:
1/ n
⎡ ⎛⎜ r ⎞⎟⎤ V 2n 2
V =V 1− ⇒α =
max ⎢⎣ ⎜⎝ R ⎟⎠ ⎥⎦ V (3 + n)(3 + 2n) (5.40)
max
Aquí se puede observar que el valor del coeficiente de energía cinética es una
función del exponente utilizado para definir el perfil de velocidad para flujo turbulento,
para un valor de n=6, ∝=1.08, mientras que para n=10, ∝=1.03. En virtud de que ∝
es razonablemente próximo a uno para números de Reynold altos, con frecuencia se
asume la unidad para efectuar los cálculos de flujo de fluidos en tuberías. Sin
embargo, para flujos en desarrollo a números de Reynolds moderados, el cambio de
energía cinética puede ser importante.
20
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 21
Considerando que el participante está familiarizado con las herramientas básicas del
análisis de flujo de fluidos, como por ejemplo las ecuaciones fundamentales
(Continuidad, Energía, Balance de energía mecánica, Segunda ley de la
Termodinámica). Nos concentraremos en el análisis del flujo de fluidos en conducto
cerrado. Si tal conducto presenta un área transversal circular, se denomina tubería,
en caso contrario haremos referencia a conductos. En todas las situaciones que
analizaremos a continuación se considera que la tubería está llena de fluido como se
muestra en la siguiente figura 5.5.
Los sistemas de tuberías que hacen posible el transporte de los fluidos desde los
centros de producción a los centros de consumo, se construyen uniendo diferentes
componentes como se muestra en la figura 5.6. Además de conducir el fluido, los
diferentes componentes sirven para cambiar su dirección, lo dividen, lo aceleran o lo
desaceleran, controlan el flujo y como consecuencia de estas acciones se producen
21
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 22
La diferencia entre los flujos laminar y turbulento en una tubería fueron establecidos
por Reynolds, el cual condujo una serie de experimentos en los cuales inyectó un
colorante en el agua que fluía en una tubería de vidrio.
22
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 23
ρVD
Re =
µ (5.41)
23
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 24
El fluido en contacto con la pared de la tubería tiene una velocidad cero (debido a la
condición de no deslizamiento). El gradiente de velocidad cerca de la pared de la
tubería se asocia a un esfuerzo cortante retardador sobre el fluido, como
consecuencia de la viscosidad del fluido. Una capa viscosa (capa limite) se
desarrolla a lo largo de la tubería hasta que los efectos viscosos se manifiesten en
todo el campo de la tubería. La región comprendida entre la sección de entrada de la
tubería y la sección donde el fluido se ha desarrollado, se conoce como sección de
entrada, cuya longitud característica se le denomina longitud de entrada.
Para el caso de un flujo laminar en una tubería circular con un perfil de velocidad en
la entrada, la longitud de entrada está dada por:
Lent
= 0.06 Re
D (5.42)
24
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 25
Lent
= 4.4 Re
1/ 6
(5.43)
D
Lent
= 0.06 × 2300 = 138 ≅ 140
D
Lent
= 4.4 × 40001 / 6 = 17.5 ≅ 18
D
Lent
= 4.4 × 100001 / 6 = 20,4 ≅ 21
D
25
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 26
Lent
= 4.4 × 200001 / 6 = 22.9 ≅ 23
D
Se puede observar que se requiere una longitud de desarrollo de flujo mayor en flujo
laminar que para flujo turbulento.
La ecuación fundamental para el flujo de fluido a lo largo de una tubería puede ser
desarrollada a partir de la ecuación general de la energía aplicable a un volumen de
control como el mostrado en la figura 5.9.
∂ V 2 gZ V 2 gZ
Q vc - ∫ eρ dV + ∫ ( h + + ) ρVdA = Wvc + ∫ ( h + + ) ρVdA (5.44)
∂t 2g g 2g g
vc ent c c sal c c
Esta ecuación establece que “el flujo de energía en forma de calor menos el cambio
de energía dentro del volumen de control, más el flujo de energía suministrado a
través de la superficie de control debe ser igual el flujo de energía que se genera en
forma de trabajo más el flujo de energía extraído a través de la superficie de control”.
26
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 27
∂
∫ eρdV = 0
∂t (5.45)
vc
Obteniéndose:
V 2 gZ V 2 gZ
Q vc + ∫ (h + + ) ρVdA = Wvc + ∫ (h + + ) ρVdA (5.46)
2g g 2g g
ent c c sal c c
Bajo condiciones estacionarias esta ecuación establece que “el cambio de energía
dentro del volumen de control es igual a cero, de allí que el flujo de energía
suministrado a través de la superficie de control, más el flujo de energía en forma de
calor debe ser igual el flujo de energía extraído a través de la superficie de control
más el flujo de energía que se genera en forma de trabajo”.
⎛ V 2 gZ ⎞⎟ ⎛ V 2 gZ ⎞⎟
Qcv + ∑ ent m⎜ h + + = Wvc + ∑ sal m⎜ h + +
⎜ 2g g ⎟ ⎜ 2g g ⎟ (5.47)
⎝ c c ⎠ ⎝ c c ⎠
27
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 28
V 2 gZ V 2 gZ
∑ ent m( h + + ) ∑ sal m(h + + )
2 gc gc 2 gc gc
Q vc
⎡ ⎛ V2 gZ ⎞⎟⎤ ⎡ ⎛ V2 gZ ⎞⎟⎤
⎜
Qcv + ⎢m h + + ⎜
⎥ = Wvc + ⎢m h + + ⎥
⎢⎣ ⎜⎝ 2 g c g c ⎟⎠⎥⎦ ent ⎢⎣ ⎜⎝ 2 g c g c ⎟⎠⎥⎦ sal
(5.48)
Dividiendo la ecuación (5.48) por el flujo de masa m, para obtener una ecuación de
la energía por unidad de masa,
⎡⎛ V 2 gZ ⎞⎟⎤ ⎡⎛ V 2 gZ ⎞⎤
qcv + ⎢⎜ h + + ⎥ = wvc + ⎢⎜ h + + ⎟⎥
⎢⎣⎜⎝ g ⎟⎥ ⎜ 2g ⎟⎥ (5.49)
2g
c c ⎠⎦ ent ⎢
⎣⎝ c
g
c ⎠⎦ sal
⎡ dV gdZ ⎤
δqcv - ⎢dh + V + ⎥ = δwvc (5.50)
⎢⎣ g
c
g ⎥
c ⎦
dV gdZ
δqcv = dh + V +
gc gc (5.51)
28
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 29
Tds ≥ δq (5.53)
Tds = δq + dI (5.54)
dh − νdP = δq + dI (5.55)
De allí,
dh − νdP − dI = δq (5.56)
dV gdZ
dh − νdP − dI = dh + V + (5.57)
gc gc
Obteniéndose:
29
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 30
dV gdZ
νdP + V + + dI = 0 (5.58)
gc gc
dP V dV g dZ dI
ν + + + =0 (5.59)
dL g c dL g c dL dL
dP Vρ dV gρ dZ dI
+ + +ρ =0 (5.60)
dL g c dL g c dL dL
⎛ dP ⎞ ⎛ dP ⎞ ⎛ dP ⎞ ⎛ dP ⎞
⎜ ⎟ +⎜ ⎟ +⎜ ⎟ +⎜ ⎟ =0 (5.61)
⎝ ⎠total ⎝ ⎠ene, cin ⎝ ⎠ene, pot ⎝ ⎠ fricción
dL dL dL dL
donde:
30
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 31
⎛ dP ⎞ Vρ dV
⎜ ⎟ = (5.62)
⎝ dL ⎠ene, cin g c dL
⎛ dP ⎞ gρ dZ
⎜ ⎟ =
⎝ dL ⎠ ene, pot g c dL (5.63)
⎛ dP ⎞ dI
⎜ ⎟ = ρ
⎝ dL ⎠ fricción dL (5.64)
En vista que las pérdidas de presión por fricción en una tubería no dependen de la
orientación de las misma, para expresar el gradiente de presión como una función
de los parámetros de la tubería y del fluido se selecciona una tubería horizontal a
través de la cual circula un fluido bajo condiciones estacionarias y completamente
desarrolladas, como se esquematiza en la figura 5.11.
τdAw
dP
(PA) (P − dL) A
dL
31
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 32
∂
Fvc = ∫ VρdV + ∫ VρVdA
∂t (5.65)
vc sc
Esta ecuación establece que “la suma de todas las fuerzas de superficie y másicas
(cuerpo), que actúan sobre un volumen de control no acelerado, es igual a la suma
neta de la relación de cambio de momento dentro del volumen de control, más la
relación neta de flujo del momento que sale a través de la superficie de control”.
∂
∫ VρdV = 0
∂t (5.66)
vc
∫ VρVdA = 0 (5.67)
sc
Por lo tanto,
Fvc = 0 (5.68)
De allí, se obtiene:
32
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 33
dP
( P ) A − τdAw − ( P − dL ) A = 0 (5.69)
dL
Resultando:
dP
AdL = τdAw (5.70)
dL
dP πD 2
( )dL = τπDdL (5.71)
dL 4
dP τ
=4 (5.72)
dL D
τ ρfV 2
f = 8gc τ=
ρV 2
8gc (5.73)
dP fρV 2 ⎛ dP ⎞ fρV 2
= ⎜ ⎟ = (5.74)
dL 2 g c D ⎝ dL ⎠f 2 g c D
dP Vρ dV gρ dZ fρV 2
+ + + =0 (5.75)
dL g c dL g c dL 2 g c D
33
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 34
∆P = f ( D , L , ε , ρ µ , V ) (5.76)
h( D, L, ε , ρ µ , V , ∆P ) = 0 (5.77)
⎛ L ε ρVd ∆P ⎞
g ⎜⎜ , , , ⎟=0
2 ⎟
⎝ D D µ ρV ⎠
(5.78)
34
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 35
De allí.
∆P ⎛ L ε ρVd ⎞
= j ⎜⎜ , , ⎟⎟
ρV 2 ⎝ D D µ ⎠ (5.79)
∆P L ⎛ ε ρVd ⎞
= k⎜ , ⎟
ρV 2 D ⎜⎝ D µ ⎟⎠ (5.80)
35
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 36
El diagrama de Moody se utiliza como una ayuda para determinar el valor del
factor de fricción f para flujo turbulento. Debe conocerse los valores del número
de Reynolds y de la rugosidad relativa de la tubería. Por consiguiente, los datos
básicos requeridos son el diámetro interior de la tubería, el material de
construcción de la tubería (que permite determinar la rugosidad equivalente de la
tubería), la velocidad del flujo, el tipo de fluido y la temperatura (con la cual se
36
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 37
37
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 38
38
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 39
Tubería lisa:
Blausius:
0.316
f = 0.25 (5.81)
Re
0.5
f = 0.0056 + 0..32 (5.82)
Re
Tubería rugosa:
Colebrook - White
1 ⎛ ε 2.51 ⎞
= −2 log⎜⎜ + ⎟
⎟ (5.83)
f ⎝ 3.7 D Re f ⎠
Nikuradse:
39
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 40
1 ⎛ ε ⎞
= −2 log⎜ ⎟
f ⎝ 3.7 D ⎠ (5.84)
D
Re f = 200
ε (5.85)
1 ⎛ 2.51 ⎞
= −2 log⎜⎜ ⎟
⎟ (5.86)
f ⎝ Re f ⎠
1
= −2 log( Re f ) − 0.8 (5.87)
f
1 ⎛ ε ⎞
= −2 log⎜ ⎟
f ⎝ 3.7 D ⎠ (5.88)
40
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 41
2
⎛ ⎞
⎜ ⎟
⎜ ⎟
1
fc = ⎜ ⎟
⎜ ⎛ ε ⎞⎟ (5.89)
⎜ − 2 log⎜ ⎟⎟
2.51
+
⎜ ⎜ 3.7 D Re f ⎟⎟
⎝ ⎝ g ⎠⎠
Los valores del factor de fricción ƒg son estimados y luego el factor de fricción ƒc es
calculado hasta que los valores de ƒg y ƒc se encuentren dentro de una tolerancia
aceptable. Las expresiones de Blausius, Drew-Koo-McAdams o Nikuradse pueden
ser utilizadas para inicializar el cálculo. En cada iteración el valor calculado es
utilizado como el valor supuesto para la próxima iteración. Con este método la
convergencia de la solución es relativamente rápida.
fricción con una tolerancia en el orden de 0.00001, realice una sensibilidad para
41
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 42
1 ⎛ ε 2.51 ⎞
= −2 log⎜⎜ + ⎟
⎟ (5.90)
f ⎝ 3.7 D Re f ⎠
⎛ ε 2.51X ⎞
g ( X ) = −2 log⎜⎜ + ⎟⎟
⎝ 3.7 D Re ⎠
X = g( X )
(5.91)
g( X ) − X = 0
Este método es una aceleración del método anterior, hecho que lo hace más
conveniente; por lo general se requieren sólo tres iteraciones. Sin embargo, la
función g(x) debe cumplir las siguientes condiciones para que exista convergencia:
1.- La función g(x) debe existir en un intervalo I = [a,b], de modo que para todo x
perteneciente a I, la función g(x) esté definida y pertenezca a I.
2.- La función g(x) debe ser continua en I
3.- La función g(x) debe ser diferenciable en I. La pendiente de g(x) debe ser menor
que 1 y mayor que –1.
i +1 g( X i ) − X i
X =X −i
(5.92)
g'( X i ) − 1
42
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 43
⎡ 2.51X i ⎤
2 ⎢⎢ Re ⎥
⎥
g'( X ) = −
ln(10) ⎢ ε 2.51X i ⎥ (5.93)
⎢ 3. 7 D + R ⎥
⎣ e ⎦
En los últimos años, con la finalidad de evitar el proceso iterativo se han desarrollado
expresiones explícitas para el factor de fricción. Una de las formas de obtener
expresiones explícitas sencillas consiste en evaluar el intervalo de valores posibles
para el factor de fricción. En el diagrama de Moody, el factor de fricción puede estar
entre 0.008 y 0.08 con un promedio de 0.044, sustituyendo este valor para el factor
de fricción en el término de la derecha de la ecuación de Colebrook-White.
Colebrook - White
1 ⎛ ε 2.51 ⎞
= −2 log⎜⎜ + ⎟
⎟ (5.94)
f ⎝ 3.7 D Re f ⎠
43
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 44
1 ⎛ ε 2.51 ⎞
= −2 log⎜⎜ + ⎟
f ⎝ 3 . 7 D Re 0.044 ⎟⎠ (5.95)
1 ⎛ ε 11.966 ⎞
⎜
= −2 log⎜ + ⎟⎟
⎝ 3.7 D Re ⎠ (5.96)
f
Colebrook-White.
Un resumen de las expresiones con precisión simple, intermedia y alta que han sido
reportadas en la literatura se presenta a continuación:
Jain
1 ⎛ ε 21.25 ⎞
= 1.14 − 2 log⎜⎜ + 0.9 ⎟⎟ (5.97)
f ⎝ D Re ⎠
Zigrang- Sylvester
1 ⎛ ε 13 ⎞
= −2 log⎜⎜ + ⎟
⎟
⎝ 3.7 D Re ⎠
f (5.98)
Haaland
44
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 45
1 ⎡⎛ ε ⎞1.11 6.9 ⎤
= −1.8 log ⎢⎜ ⎟ + ⎥
f ⎣⎢⎝ 3. 7 D ⎠ R ⎥
e ⎦
(5.99)
Serghides
1 ⎛ ε 12 ⎞
= −2 log⎜⎜ + ⎟⎟ (5.100)
f ⎝ 3.7 D Re ⎠
Chen
1 ⎛ ε 2.51 ⎞⎟
= −2 log⎜ +
f ⎜ 3.7 D R f ⎟ (5.102)
⎝ e s ⎠
45
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 46
las expresiones del factor de fricción obtenido por las ecuaciones de Jain, Zingrang –
Silvestre, Haaland o Chen.
Jain
1 ⎛ ε 21.25 ⎞
= 1.14 − 2 log⎜⎜ + 0.9 ⎟⎟ (5.103)
fs ⎝ D Re ⎠
Zigrang- Sylvester
1 ⎛ ε 13 ⎞
= −2 log⎜⎜ + ⎟⎟
fs ⎝ 3.7 D Re ⎠ (5.104)
Haaland
1 ⎡⎛ ε ⎞1.11 6.9 ⎤
= −1.8 log ⎢⎜ ⎟ + ⎥
⎢⎣⎝ ⎠
3. 7 D Re ⎥⎦ (5.105)
fs
1 ⎛ ε 2.51 ⎞⎟
= −2 log⎜ +
f ⎜ 3.7 D R f ⎟ (5.106)
⎝ e s ⎠
1 ⎡⎛ ε ⎞1.11 6.9 ⎤
= −1.8 log ⎢⎜ ⎟ + ⎥ (5.107)
fs ⎣⎢⎝ 3 . 7 D ⎠ R ⎥
e ⎦
46
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 47
⎛ ε 4.518 ⎛⎜ ⎡⎛ ε ⎞
1.11
6.9 ⎤ ⎞⎟ ⎞⎟
1
= −2 log ⎜ − log ⎢⎜ ⎟ + ⎥
f ⎜ 3. 7 D Re ⎜⎝ ⎢⎣⎝ 3.7 D ⎠ Re ⎥⎦ ⎟⎠ ⎟ (5.109)
⎝ ⎠
1 ⎛ ε 2.51 ⎞
= −2 log⎜ + ⎟
f ⎜ 3 .7 D R f ⎟ (5.110)
⎝ e 1 ⎠
1 ⎛ ε 2.51 ⎞
= −2 log⎜ + ⎟
f1 ⎜ 3.7 D R f ⎟ (5.111)
⎝ e s ⎠
47
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 48
1 ⎛ ε 2.51 ⎞
= −2 log⎜ + ⎟
f ⎜ 3 .7 D R f ⎟ (5.112)
⎝ e 1 ⎠
⎛ ε 4.518 ⎛⎜ ⎡⎛ ε ⎞
1.11
6.9 ⎤ ⎞⎟ ⎞⎟
1
= −2 log ⎜ − log ⎢⎜ ⎟ + ⎥
f1 ⎜ 3 .7 D Re ⎜⎝ ⎣⎢⎝ 3.7 D ⎠ Re ⎦⎥ ⎟⎠ ⎟ (5.113)
⎝ ⎠
Resultando:
⎧ ⎫
⎪ ε 5.02 ⎡⎢ ⎛⎜ ε 4.518 ⎛⎜ ⎡⎛ ε ⎞ 6.9 ⎤ ⎞⎟ ⎞⎟⎤⎥ ⎪
1.11
1
=−2 log ⎨ − log − log ⎢⎜ ⎟ + ⎥ ⎬ (5.114)
f ⎪⎩ 3.7 D Re ⎢ ⎜ 3.7 D Re ⎜⎝ ⎣⎢⎝ 3.7 D ⎠ Re ⎦⎥ ⎟⎠ ⎟⎥ ⎪
⎣ ⎝ ⎠⎦ ⎭
Existen otros tipos de expresiones para determinar el factor de fricción, las cuales
pueden ser encontradas como referencias técnicas del tópico.
Serghides
−2
⎡ ( B − A) ⎤ 2
f = ⎢A − ⎥
⎣ C − 2 B + A ⎦ (5.115)
48
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 49
⎛ ε 12 ⎞
A = −2 log⎜⎜ + ⎟⎟ (5.116)
⎝ 3.7 D Re ⎠
⎛ ε 2.51A ⎞
B = −2 log⎜⎜ + ⎟⎟
⎝ 3.7 D Re ⎠ (5.117)
⎛ ε 2.51B ⎞
C = −2 log⎜⎜ + ⎟⎟
⎝ 3.7 D Re ⎠ (5.118)
−2
⎡ ( A − 4.781) 2 ⎤
f = ⎢4.781 − ⎥
⎣ B − 2 A + 4.781 ⎦ (5.119)
Churchill
⎡ 1 ⎤
f = 8⎢ A −
⎣ (B + C )1 / 12 ⎥⎦ (5.120)
12
⎛ 8 ⎞
A = ⎜⎜ ⎟⎟
⎝ Re ⎠ (5.121)
49
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 50
16
⎡ ⎧⎪⎛ 7 ⎞ 0.9 0.27ε ⎫⎪⎤
B = ⎢2.457 ln ⎨⎜⎜ ⎟⎟ + ⎬⎥ (5.122)
⎢ ⎪⎩⎝ e ⎠
R D ⎪⎭⎥⎦
⎣
16
⎛ 37530 ⎞
C = ⎜⎜ ⎟⎟
⎝ Re ⎠
(5.123)
Cuarta iteración
1 ⎛ ε 2.51 ⎞⎟
= −2 log⎜ +
f ⎜ 3.7 D R f ⎟ (5.124)
⎝ e 3 ⎠
Tercera iteración
1 ⎛ ε 2.51 ⎞
= −2 log⎜ + ⎟
f3 ⎜ 3 .7 D R f ⎟ (5.125)
⎝ e 2 ⎠
Segunda iteración
50
JAGP 1805 – FACILIDADES DE SUPERFICIE EN INSTALACIONES PETROLERAS 51
1 ⎛ ε 2.51 ⎞
= −2 log⎜ + ⎟
f2 ⎜ 3 .7 D R f ⎟ (5.126)
⎝ e 1 ⎠
Primera iteración
1 ⎛ ε 2.51 ⎞⎟
= −2 log⎜ +
f1 ⎜ 3.7 D R f ⎟ (5.127)
⎝ e 0 ⎠
1 ⎡⎛ ε ⎞1.11 6.9 ⎤
= −1.8 log ⎢⎜ ⎟ + ⎥
⎢⎣⎝ ⎠
3 .7 D Re ⎥⎦ (5.128)
f0
para determinar el factor de fricción. Desarrolle una hoja de cálculo que permita
Compare sus resultados para las diferentes expresiones y establezca sus propias
conclusiones.
D
Re f = 200
ε
51
Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.
Alternative Proxies: