OPERACIONES DE FLUJO DE FLUIDOS

TEMA 2: FLUJO INTERNO DE FLUIDOS

Índice
TEMA 2: FLUJO INTERNO DE FLUIDOS .......................................................................................... 1
1. Balances de materia, energía y cantidad de movimiento ................................................. 2
1.1 Ley de conservación de la masa ................................................................................ 2
1.2 Concepto de tensor de fuerzas por unidad de superficie en un fluido. .................... 4
1.3 Ley de conservación de la cantidad de movimiento ................................................. 5
2. Distribución de velocidades en el interior de conducciones cilíndricas............................ 8
2.1 Flujo laminar.............................................................................................................. 9
2.2 Flujo turbulento....................................................................................................... 11
3. Conservación de energía mecánica. Ecuación de Bernoulli ............................................ 12
3.1 Balance de energía total. Análisis de los diferentes tipos de energía. .................... 12
3.2 Balance de energía mecánica .................................................................................. 14
3.3 Ecuación de Bernoulli para líquidos. Perfil hidráulico ............................................. 16
4. Pérdidas de energía mecánica ........................................................................................ 18
4.1 Pérdidas de energía mecánica en conducciones cilíndricas.................................... 18
4.2 Pérdida de energía mecánica de accidentes en conducciones ............................... 21

1
OPERACIONES DE FLUJO DE FLUIDOS

1. Balances de materia, energía y cantidad de movimiento

El desarrollo de una descripción analítica de un fluido en movimiento, se basa en la expresión

de las leyes físicas relacionadas con el flujo de fluidos, en una forma matemática adecuada.
Hay tres leyes físicas que se aplican a todos los flujos, independientemente de la naturaleza del
fluido que se esté considerando. Estas leyes son: la conservación de la masa, del momento
(segunda ley de Newton del movimiento) y de la energía (primera ley de la termodinámica).

De manera general, las leyes de conservación de cualquiera de las tres propiedades se

establecen aplicando un balance de la propiedad, de la forma

 velociad   velociad   velociad 

     
 entrada al  salida del    acumulación en el 
 volumen de control   volumen de control   volumen de control 
     

1.1 Ley de conservación de la masa

Consideremos un elemento de volumen de un flujo de fluido circulando en el interior de una

conducción (¡Error! No se encuentra el origen de la referencia.).

Figura 1. Elemento de volumen x, y, z, a través del que pasa el flujo de fluido.

Aplicando un balance de materia al volumen de control se tiene

 velociad 
 masa entrada   masa salida   
     acumulación 
 por unidad de tiempo   por unidad de tiempo   de masa 
 

2
OPERACIONES DE FLUJO DE FLUIDOS

Considerando el flujo de masa en dirección del eje x, perpendicular a las caras de sección y
z, el caudal másico que entra en el volumen de control es  vx  x yz y el que lo

abandona  vx  x x yz . Para el resto de las caras se pueden escribir expresiones
similares.

  
La velocidad de acumulación de masa en el elemento de volumen es x y z  .
 t 

Utilizando estas expresiones, el balance de materia quedaría de la siguiente forma

   
x y z     v x  x y z   v x  x x  y z 
 t  
  v  x z   v   x z  (1)
 y
y
y
y y 
  v z  x y   v z   x y
 z z z 

Dividiendo la ecuación por x y z, cuando tienden a cero, se obtiene

      
   v x  vy  v z  ó    v  (2)
t  x y z  t

    
donde  es el operador nabla  i  j  k  y v es el vector velocidad
 x y z 
v i  v
x y j  vz k  .

Al producto    v  se le conoce como divergencia de  v , que a veces se escribe como div. v .

La ecuación anterior se conoce como ecuación de continuidad (es equivalente a un balance

microscópico de materia total (acumulación= -(entrada-salida) expresado en kg/sm3)).

Si se considera flujo unidimensional en dirección del eje x, la ecuación de continuidad queda

   v x 
 (3)
t x

Si el sistema se encuentra en régimen estacionario, la masa que pasa a través del volumen de
control permanece constante con el tiempo, no hay acumulación, por tanto


0 (4)
t

Con lo que la ecuación de continuidad unidimensional sería

0  d  v  (5)

3
OPERACIONES DE FLUJO DE FLUIDOS

Integrando la ecuación (5) entre los puntos 1 y 2 de la conducción de la Figura 1, se obtiene

2
0   d  v    v 2   v 1 (6)
1

por lo tanto

 v 2   v 1  v  cte (7)

El producto v es la densidad de flujo total de materia o gasto másico, que denominaremos G

(kg/s m2). Multiplicando el gasto másico G, por el sección transversal de la conducción, S (m2),
obtenemos el caudal másico, m (kg/s), de fluido que circula por la conducción.

 kg 
m  G1 S1  G2 S2  1v1 S1  2v2 S2  cte   (8)
 s 

El producto v, además de representar una densidad de flujo de materia (kg/sm2) representa
también una concentración de cantidad de movimiento ( (kg. m/s)/m3).

La ecuación de continuidad también se puede escribir en coordenadas cilíndricas y esféricas.

1.2 Concepto de tensor de fuerzas por unidad de superficie en un fluido.

En un sólido rígido en reposo, al aplicar una fuerza en una determinada dirección se puede
producir un movimiento de traslación/rotación de éste dependiendo de la dirección, pero el
sólido continua siendo rígido. En un fluido la situación es diferente, y puede deformarse y
adaptarse al recipiente que lo contiene. Las fuerzas que se pueden transmitir dentro del fluido
pueden ser de dos tipos: las fuerzas de presión que se transmiten perpendicularmente al plano
considerado y las fuerzas de rozamiento debido a la diferente velocidad en ambas caras de la
sección considerada, pero puede ser normal o tangencial según la cara considerada. Este
conjunto de fuerzas se expresa mediante un tensor de esfuerzos, que se puede expresar con la
siguiente matriz:

𝑝 + 𝜏𝑥𝑥 𝜏𝑥𝑦 𝜏𝑥𝑧

( 𝜏𝑦𝑥 𝑝 + 𝜏𝑦𝑦 𝜏𝑦𝑧 )
𝜏𝑧𝑥 𝜏𝑧𝑦 𝑝 + 𝜏𝑧𝑧

Para un elemento de volumen cúbico, tal como se expresa en la siguiente figura:

4
OPERACIONES DE FLUJO DE FLUIDOS

Figura 2. Elementos del tensor de esfuerzos sobre un cubo cuando el fluido se mueve en la dirección x.

xx, yy y zz son los componentes normales o perpendiculares a la cara que actúan, igual que la
presión pero debido a diferentes causas: los tensores porque hay un gradiente de cantidad de
movimiento (sería cero para un fluido en reposo) y la presión por el empuje de las moléculas
debido a su situación energética (tendrá un valor aunque el fluido esté en reposo).

Los otros componentes del tensor son tangenciales.

Los componentes del tensor tienen un doble sentido, aunque están relacionados:

ab: Es un tensor, fuerza transmitida por unidad de superficie, donde el primer subíndice “a”
representa el plano sobre el que actúa, y el segundo subíndice “b” la dirección de la velocidad
o de la cantidad de movimiento. Las unidades serían las mismas que las de la presión (Pa).

ab: Representa una densidad de flujo de cantidad de movimiento, donde el primer subíndice
“a” representa la dirección en la que se transmite la cantidad de movimiento, y el segundo
subíndice “b” la dirección de la velocidad o de la cantidad de movimiento. Las unidades serían
(kg. m/s)/(m2.s) que coincide con las unidades de la presión (Pa).

Realmente ambos conceptos son equivalentes, y el orden de los subíndices no es importante

porque se cumple que ab = ba.

1.3 Ley de conservación de la cantidad de movimiento

Para obtener la ecuación de movimiento habrá que aplicar un balance de cantidad de

movimiento sobre un elemento de fluido

 velocidad   momento   momento 

       Fuerzas externas 
 acumulación 
  entrada 
  salida    sobre el fluido 
 momento   por unidad de tiempo   por unidad de tiempo   
     

Consideremos una conducción cilíndrica (

5
OPERACIONES DE FLUJO DE FLUIDOS

Figura 3) y un elemento volumétrico de fluido como el de la Figura 4, donde se considera flujo

unidimensional en la dirección del eje z (normalmente en fluidos el símbolo z se emplea para
la cota de altura, pero también es el símbolo utilizado en coordenadas cilíndricas para
determinar la dirección axial.

Figura 3. Elemento de fluido.

r
L
r

V=2rLr

Figura 4. Elemento de volumen

El elemento de volumen tiene como tres dimensiones 2r, r y L, y por tanto su volumen es
2rrL. En principio, los tensores que aparecen son los correspondientes a las coordenadas
cilíndricas, y actuará también la presión. Para este caso considerado, el movimiento tiene lugar
en la dirección z y la transmisión de cantidad de movimiento en la dirección r, siendo el tensor
de transmisión de cantidad de movimiento rz, y el área de transmisión de cantidad de
movimiento es 2rL. Por otra parte, hay que considerar en la cara z=0 y z=L, actuarán las
fuerzas de presión y zz. Se va a considerar que zz vale cero

El balance de flujo de cantidad de movimiento aplicado al elemento de volumen vendrá dado

por


  vz 2 r rL   2 rL rz r   2 rL rz r r  2 rr  vz2 1  2 rr  vz2  2
t (9)
2 r r ( p1  p2 )  g z  2 r rL

Se asume también que vz es la misma en la sección 1 que en la sección 2, por lo que la

ecuación anterior se simplifica a

6
OPERACIONES DE FLUJO DE FLUIDOS


  vz 2 rrL   2 rL rz r   2 rL rz r r  2 rr ( p1  p2 )  g z  2 rrL (10)
t

El valor de gz es el componente de la gravedad que hace aumentar la cantidad de movimiento.

Dividiendo la expresión anterior por 2rrL, y tomando límites cuando t y r tienden a cero
se obtiene

vz p  p1 1 
  2   r rz    g z (11)
t L r r

La ecuación anterior es un balance microscópico, obtenido al dividir por el elemento de

volumen, y por tanto con unidades de flujo de cantidad de movimiento por unidad de volumen
(kg m/s)/(s m3).

Si no se hubiese hecho ninguna simplificación y se hubiese elegido un diferencial de volumen

rrz (Figura 5), el balance microscópico de cantidad de movimiento tendría tres
ecuaciones para cada una de las tres coordenadas r,  (ó ) y z. La componente z tendría la
siguiente forma:

 vz v v v v  p  1  1  
  vr z   z  vz z       r rz    r  zz    g z (12)
 t r r  z  z  r r r  z 

Figura 5. Diferencial de volumen de fluido

Se puede observar que si se considera que vz sólo depende de r, y que vz y v son nulos, y
teniendo en cuenta que bajo estas consideraciones r y zz son nulos al no existir gradientes de
velocidad, la ecuación (12) se convierte en la ecuación (11). Si además se considera régimen
estacionario, la ecuación (11) se convierte en

p2  p1 1 
0   r rz    g z (13)
L r r

ecuación que se utilizará posteriormente.

7
OPERACIONES DE FLUJO DE FLUIDOS

Las ecuaciones anteriores son válidas tanto para régimen laminar como turbulento, ya que no
se ha introducido ninguna ley de velocidad.

2. Distribución de velocidades en el interior de conducciones cilíndricas

Consideremos un fluido incompresible ( cte) que circule en estado estacionario por el

interior de una conducción cilíndrica hacia arriba. En este caso, el componente gz será
negativo, y por tanto;

gz  g cos (14)

La ecuación de movimiento (componente z) quedará

1 d(r rz ) p p   p  p  gL cos 

   2 1   g cos    2 1  (15)
r dr  L   L 

como

L cos  h2  h1 (h equivale a la cota de altura, que en otras ocasiones se expresa como z) (16)

se obtiene

1 d(r rz )  p  p  g  h2  h1  
  2 1  (17)
r dr  L 

Y si se define p+ = p + gh, la ecuación (19) se transforma en:

1 d (r rz )  p  p1   p1  p2 

  2  
r dr  L   L 
(18)
 p1  p2 
  rdr  d (r rz )
 L 

Integrando la ecuación (18)

 p1  p2  C2
 rz   r  (19)
 2L  r

Si r = 0, como rz, C2=0, con lo que

 p1  p2 
 rz   r (20)
 2L 

es decir, el esfuerzo cortante aumenta linealmente desde el centro de la conducción hasta la

pared de la misma (Figura 6-(a)).

8
OPERACIONES DE FLUJO DE FLUIDOS

Obsérvese, una vez más, que el esfuerzo cortante tiene las dimensiones de una presión, y que,
en este caso particular, la proporcionalidad entre el esfuerzo cortante y el radio es,
precisamente, la pérdida de presión con contribución gravitatoria que experimenta el fluido
por unidad de longitud. La ecuación (22) pone de manifiesto, en definitiva, la estrecha relación
que existe entre la variación de presión que experimenta un fluido que circula por una
conducción y el esfuerzo cortante resultante de la interacción del fluido con la pared de la
conducción.

2.1 Flujo laminar

En el caso de que el fluido circule a número de Reynolds menor de 2100-2400, el régimen es

laminar y el transporte de cantidad de movimiento rz tiene lugar por un mecanismo
molecular. Aplicando la ley de Newton en coordenadas cilíndricas se tiene:

 vz   dvz 
 rz         (21)
 r   dr 

Sustituyendo esta última expresión en la ecuación (20), se obtiene

 p  p2 
dv z    1  rdr (22)
 2 L 

e integrando, con la condición: r = R, vz= 0, se obtiene:

(p1  p2 )R2   r  

2

vz  1     (23)
4 L   R  

Como puede apreciarse, la relación vz = vz(r) es una parábola (Figura 6-(b)).

rz

Figura 6. Variación del esfuerzo cortante (a) y la velocidad (b) con el radio.

Es evidente que la velocidad máxima tiene lugar en el centro de la conducción (r=0), y su valor
viene dado por:

(p1  p2 )R2

vz max  (24)
4 L

La velocidad media será:

9
OPERACIONES DE FLUJO DE FLUIDOS

(vz )m 
 v dS   v 2 rdr  (p
S z z

1  p2 )R2
(25)
 dS R 2
8 L
S

es decir, la mitad de la velocidad máxima

En la Figura 6 se aprecia perfectamente, comparando el perfil de velocidad con el del esfuerzo

cortante, que la velocidad máxima tiene lugar donde el esfuerzo cortante es nulo, y que se
anula allí donde el esfuerzo cortante es máximo, es decir, en la pared.

Por último, el caudal volumétrico vendrá expresado por:

 (p1 - p2 )R 4
Q  (v z )m S  (26)
8L

relación que recibe el nombre de ecuación de Hagen-Poiseuille, y que indica, que el caudal
volumétrico es directamente proporcional a la pérdida de presión por unidad de longitud y al
radio a la cuarta potencia, e inversamente proporcional a la viscosidad.

Esta ley, que resume el comportamiento de un fluido newtoniano que fluye laminarmente por
una conducción cilíndrica, fue puesta de manifiesto simultánea e independientemente por
Gotthilf Hagen (1797-1884), ingeniero alemán que dedicó gran parte de su vida a trabajos de
hidráulica en ríos y puertos, y por Jean-Louis-Marie Poiseuille (1799-1869), médico francés que
contribuyó valiosamente al conocimiento de la circulación sanguínea. De campos e intereses
aparentemente alejados se confluyó en una ley física que, lógicamente había de tener una
aplicación general.

Conviene hacer hincapié en que esta última ecuación, así como el perfil de velocidad en el que
se basa (ec. (23)), son válidas si se cumplen estrictamente todas las hipótesis que se han
venido haciendo o han quedado implícitas, y que se resumen a continuación:

 Medio continuo. Es decir, que el tamaño de las moléculas, o el recorrido libre medio si
se trata de un gas, sea despreciable en comparación con el diámetro de la conducción.
Esta hipótesis podría no cumplirse, y de hecho así ocurre a veces, en lo que respecta,
por ejemplo, al flujo de materia en los poros de un catalizador sólido.
 Estado estacionario.
 Régimen laminar.
 Fluido Newtoniano.
 Densidad constante. Esta hipótesis implica que las conclusiones obtenidas sean más
aplicables a líquidos que a gases.
 Efectos finales despreciables. En este sentido, la longitud de entrada para que se
formen los perfiles parabólicos es Le= (0.035) D (Re), y si esta longitud no fuera
despreciable frente a la longitud total de la conducción habría que efectuar una
corrección.
 No hay deslizamiento en la pared. El caso opuesto sería el de un flujo ideal cuya
viscosidad fuera nula, que se deslizaría sobre la pared sin rozamiento, es decir, sin que
hubieran fuerzas de cizalla.

10
OPERACIONES DE FLUJO DE FLUIDOS

El estudio realizado en este apartado pone de manifiesto una pauta para resolver problemas
de flujo laminar que se resume a continuación:

 Establecer el modelo con las hipótesis pertinentes.

 Aplicar la ecuación de continuidad, la ecuación de movimiento y la ley de Newton.
 Delimitar las condiciones de contorno e integrar las ecuaciones diferenciales
resultantes.

En el caso aquí discutido la densidad y la viscosidad del fluido se han considerado constantes.
En los problemas en los que tales hipótesis no fueran admisibles habría que tener en cuenta,
además, la ecuación de estado del fluido compresible y una ecuación de viscosidad. Si el fluido
no fuera newtoniano, se sustituiría la ley de Newton por otras expresiones.

2.2 Flujo turbulento

El flujo turbulento es mucho más frecuente que el laminar, sin embargo no se pueden expresar
con la misma precisión los flujos de materia, cantidad de movimiento y energía debidos a la
turbulencia. La integración de las ecuaciones diferenciales es más compleja. Se han deducido
diversas ecuaciones empíricas, basadas en las distintas teorías propuestas para la turbulencia,
para representar el perfil de velocidades en el caso de un flujo turbulento estacionario de un
fluido newtoniano por una conducción cilíndrica. Ninguna de las teorías propuestas es
completamente satisfactoria y las ecuaciones indicadas solo serán aproximadas, aunque
algunas representan el perfil de velocidades con bastante precisión.

Desde un punto de vista cualitativo, puede decirse que al aumentar el número de Reynolds,
desde valores inferiores al indicado como valor crítico (2100), el perfil parabólico de
velocidades típico del flujo laminar (ecuación (23)), se achata progresivamente llegando a
convertirse en una superficie prácticamente plana, perpendicular a las paredes del tubo.
Solamente en la porción de fluido próxima a la pared existe un perfil de velocidades ya que la
velocidad del fluido en contacto con la pared es 0 (Figura 7).

Figura 7. Forma general del perfil de velocidad para el flujo turbulento.

11
OPERACIONES DE FLUJO DE FLUIDOS

Dado el perfil prácticamente plano de velocidades en flujo turbulento, se puede considerar

que la velocidad máxima prácticamente coincide con la velocidad media a través de la sección
del tubo

vz max  vz m (27)

3. Conservación de energía mecánica. Ecuación de Bernoulli

3.1 Balance de energía total. Análisis de los diferentes tipos de energía.

Considérese el sistema de flujo representado en la ¡Error! No se encuentra el origen de la

eferencia., consistente en una conducción de sección circular (que para mayor generalidad se
ha supuesto variable con la posición, pero no con el tiempo), por la que fluye, en régimen
estacionario, un fluido con caudal másico m (kg/s). Se supone que existe una máquina (bomba
o turbina) que intercambia energía con el fluido y que al mismo tiempo este puede
intercambiar también energía calorífica con el exterior a través de las paredes de la
conducción

2
1

P2
2
P1 V2
1
W, Q
V1

Figura 8. Sistema de flujo estacionario.

Aplicando un balance macroscópico de energía total entre las secciones 1 y 2:

 V2 V2   J 
 p1 1  p2 2  g(z1  z2 )  W   Q  (u2  u1 )   2  1   kg  (28)
 22 21   

Los términos que aparecen en la ecuación (28), son:

 p1 1  p 2 2 : trabajo que las fuerzas de presión externas ejercen sobre la unidad de
masa de fluido al pasar desde la sección S1 a la S2.

 g( z1  z2 ) : trabajo que la fuerza de la gravedad ejerce sobre la unidad de masa de

fluido al pasar desde la sección S1 a la S2.

12
OPERACIONES DE FLUJO DE FLUIDOS


 W : trabajo específico que la máquina ejerce sobre la unidad de masa de fluido entre
las secciones S1 y S2; es positivo o negativo si el trabajo hace aumentar o disminuir la
energía del sistema.

 Q : calor específico intercambiado con los alrededores por la unidad de masa de fluido
al pasar desde la sección S1 a la S2; será positivo o negativo según sea ganado o perdido
por el mismo.

 (u2-u1): variación de energía interna de la unidad de masa de fluido al pasar desde la

sección S1 a la S2.

 V22   V12 
    : variación de energía cinética de la unidad de masa de fluido al pasar
 2 2   21 
de la sección S1 a la S2. El valor de  depende del régimen de circulación y considera que
no todo el fluido tiene la misma velocidad que la velocidad media a través de toda la
sección. Para régimen laminar (Re=VD/ < 2100) el valor de  es de 0.5, y para
régimen muy turbulento (Re < 104) el valor de  es de 0.96, pero se suele tomar el valor
de 1.

Centramos ahora la atención en la masa mT existente entre las secciones 1 y 2 de la Figura 8.

Supongamos que transcurrido un instante de tiempo d, esta porción de materia se desplaza
desde 1-2 hasta 1’-2’, durante el que se desarrolla la transformación representada en la Figura
9.

2 d 2
1 1

2’ 2’
1’ 1’

Sistema a Sistema b

Figura 9. Situación de la porción mT desde 1-2 a 1’-2’ tras un tiempo d.

Según el balance de materia (ecuación (8)), la masa de fluido que atraviesa cualquier sección
del sistema tras el instante d, será md. Estableciendo diferencias entre los estados físicos a y
b, es fácil distinguir que el cambio existente entre ambos es el transporte de la porción de
materia md de la porción 1-1’ en el estado a, a la porción 2-2’ en el estado b. Ambas
situaciones energéticas están relacionadas por el primer principio de la termodinámica, y para
obtenerlo, bastará multiplicar por tal masa todos los términos del balance de energía
(ecuación (28)):

ˆ   md Q  md (u2  u1 )  md  V2  V1 

2 2
md (p1 1  p2 2 )  g(z1  z2 )  W  J (29)
 22 21 

Obteniéndose

13
OPERACIONES DE FLUJO DE FLUIDOS

Wea b  Q' a b  (U b  U a )  (E cb  E ca ) (30)

es decir:

Weab  md (p1 1  p2 2 )  g(z1  z2 )  W

ˆ   trabajo de las fuerzas exteriores
 (31)

Q' ab  mdQ (32)

(U b  U a )  md (u 2  u 1 ) (33)

 V22 V12 
( E cb  E ca )  md   (34)
 2 2 21 

3.2 Balance de energía mecánica

La ecuación de conservación de la energía mecánica se puede establecer aplicando una de las

leyes fundamentales de la mecánica: el teorema de las fuerzas vivas. Según este teorema, el
incremento de energía cinética en un sistema material en movimiento es igual a la suma de
 
los trabajos realizados por las fuerzas exteriores sobre el sistema Weab , más la suma de

 
los trabajos realizados por las fuerzas interiores que actúan sobre el mismo Wiab , es decir:

(E cb  E ca )  Wea b  Wia b (35)

Todos los términos de las ecuación de conservación (35) están evaluados (ecuaciones (30) a
(34)) excepto el trabajo de las fuerzas internas Wi a b . De él nos ocupamos a continuación.

El término Wi a b , está integrado por dos tipos de fuerzas:

 Trabajo debido a las fuerzas ejercidas mutuamente entre los diversos elementos del
fluido, normalmente a sus superficies de separación o fuerzas de presión interiores
(Wpi).
 Trabajo debido a las fuerzas ejercidas mutuamente entre los diversos elementos del
fluido, tangencialmente a sus superficies de separación o fuerzas rasantes de
viscosidad, siempre negativas (-W).

Por tanto:

Wia b  Wpi  W (36)

El trabajo realizado por las fuerzas de presión interiores (Wpi) resultará de restar del total
(Wpt), debido a todas las fuerzas de presión interiores y exteriores, el que corresponde
solamente a las exteriores (Wpe), que según la ecuación (31) es:

Wpe  md (p1 1  p 2 2 ) (37)

14
OPERACIONES DE FLUJO DE FLUIDOS

es decir:

Wpi  Wpt  Wpe  Wpt  md( p11  p 2  2 ) (38)

r
A fin de evaluar Wpt considérese un pequeño elemento de volumen del fluido con base dS ,
sobre la sección transversal h-h’. El elemento diferencial en h alcanza h’ tras el instante d. Las
generatrices de dicho elemento son paralelas a la velocidad local del fluido V y que terminan
en la sección transversal h' muy próxima a h (Figura 10). Dado que las generatrices son
paralelas a la velocidad del fluido, las fuerzas normales a la superficie lateral del pequeño
cilindro no desarrollarán trabajo, debiendo considerarse solamente el trabajo realizado por las
r
fuerzas de presión normales a sus dos bases dS localizadas sobre los planos h y h'. Durante el
intervalo de tiempo d, tal trabajo, producto de la fuerza por el espacio, será:
 

dWpt = Fuerza × distancia = ( p  p')dS • Vd    (39)

e integrando para abarcar toda la sección transversal S entre h y h':

   
S (p  p' )dS • Vd  (p  p' )dSV • dS  (p  p' )dS
h h h
dQ v (p  p' )md (40)

Figura 10. Esquema para el cálculo del trabajo realizado por todas las fuerzas de presión (internas y
externas).

La integral del segundo miembro de esta igualdad representa el caudal volumétrico a través de
la sección h (caudal volumétrico=m), como se muestra en la última igualdad de la ecuación
(40).

Finalmente, integrando sobre la transformación de a hasta b del volumen de control:

p2
Wp t   (p  p' )md    mddp  md  dp (41)
a b a b
p1

expresión del denominado trabajo técnico.

De las expresiones (38) y (41) se deduce:

Wpi  md    dp  ( P1 1  P2 2 )   md  pd

P2 v2
(42)
 P1  v1

15
OPERACIONES DE FLUJO DE FLUIDOS

expresión del trabajo de expansión o compresión realizado por las fuerzas interiores, que será
nulo en el caso de los fluidos incompresibles, como los líquidos, al ser constante su volumen
específico.

Si se representa por F el trabajo desarrollado por las fuerzas rasantes de viscosidad sobre la
unidad de masa de fluido entre las secciones S1 y S2, es decir, si se supone: W=md(F), las
ecuaciones (36) y (42) conducen a la expresión:

Wi a b  md  pd  F

v2
(43)
 v1 

De las ecuaciones (31), (34), (35) y (43) se deduce la ecuación de conservación de la energía
mecánica:

 V22 V12  P2  J 
    g(z2  z1 )  P1  dp  F  W
ˆ
 kg  (44)
2
 2 2 1   

Esta ecuación representa el balance de energía cinética o mecánica en régimen estacionario,

denominada también ecuación de Bernoulli para los fluidos reales, puesto que si el fluido
fuera incompresible (=constante), y sin viscosidad (F=0), y no hubiera máquina (W=0), se
simplificaría a la deducida por dicho autor para fluidos ideales.

3.3 Ecuación de Bernoulli para líquidos. Perfil hidráulico

En el Sistema Internacional de unidades, todos los términos de la ecuación (44) están

expresados en J/kg. En el caso de líquidos, como su comportamiento es prácticamente
incompresible, la forma de la ecuación de Bernoulli (44) al integrar el término de presión es:

 V22 V2  P P  J 
  1   g(z2  z1 )  2 1  F  W
ˆ
 kg  (45)
2
 2 2 1    

Si se dividen todos los términos de la ecuación de Bernoulli para líquidos por la aceleración de
la gravedad, la ecuación está expresada en unidades de energía por unidad de peso de fluido
(J/N, dimensionalmente m), con un sentido físico idéntico al anterior. Esta forma de expresar
el balance es el clásico de la Hidráulica, recibiendo los distintos términos el nombre de cargas.

1  V22 V12  P2  P1 J 
    (z2  z1 )   hf  H   m (46)
g  22 21  g N 

cuyos términos reciben el nombre, respectivamente, de cargas cinéticas o de velocidad, cargas

estáticas o potenciales o simplemente cotas, cargas de presión, pérdida de carga (hf = F/g) y
carga o altura de la bomba (H = W/g).

En muchos manuales de Hidráulica, se presenta el perfil hidráulico para conducciones de agua,

que no es más que una representación de las cargas que tiene un fluido a lo largo del sistema
hidráulico. Tal como se observa en la Figura 11, se representa la variación de la:

16
OPERACIONES DE FLUJO DE FLUIDOS

- línea de carga estática, que representa el nivel de energía más alto de la tubería, o
un nivel de referencia como por ejemplo el del depósito del cual se parte. Esta línea
es siempre horizontal y también se le conoce como plano inicial o línea de carga
inicial.
- línea de energía, suma de las cargas estáticas debidas a la cota, a la energía de
presión y la cinética. Su diferencia con respecto a la estática representa la pérdida de
carga de rozamiento respecto a ese punto.
- línea piezométrica, que representaría la carga manométrica, suma de la carga de
presión y la cota (sería la altura a la que estaría un piezómetro o manómetro de rama
abierta). Puede darse algún caso donde la línea piezométrica cruce la trayectoria
(cuando la carga de la presión atmosférica no se ha considerado), y en este caso se
podría deducir que la presión es inferior a la atmosférica.

Estas cargas se representan en un gráfico (Figura 11) frente a la trayectoria de la conducción,

cuya distancia al plano de referencia (z=0) representa la posición de la conducción.

Figura 11. Perfil hidráulico de un sistema de flujo de fluidos.

Otra forma de escribir el balance de energía mecánica para líquidos entre una sección 1 y una
sección cualesquiera, donde no hay ninguna bomba ni turbina es:

1  V12  p1 1  V 
2
p 1
 z    z   F (47)
g  21  g g  2  g g
1
 

El primer término representa la carga total en el extremo 1, y el segundo término las cargas en
cualquier sección, donde F representa la pérdida de energía mecánica entre el extremo 1 y la
sección considerada. Otra forma de expresarlo sería:

17
OPERACIONES DE FLUJO DE FLUIDOS

1  V 2  1  V12  p p 1
z1  z       1  F (48)
g  2  g  21   g  g g

4. Pérdidas de energía mecánica

4.1 Pérdidas de energía mecánica en conducciones cilíndricas

Considérese un fluido incompresible que circula por una conducción cilíndrica de sección
constante, pudiendo ser tanto en régimen laminar como turbulento (Figura 12).

P
2R 2

z2
x

P1
r=0
z1

Figura 12. Circulación de un fluido incompresible por una conducción cilíndrica.

Aplicando el balance de cantidad de movimiento a los extremos 1 y 2:

p1 R2  p2 R2   R 2L g cos   rx r R 2 RL  0

(49)
 (p   gz2 )  (p1   gz1 ) 
 rx r 0
  2 R
 2L 

El factor de fricción de Fanning (f) se define como:

 rx
f r R
(50)
1 2
V
2

Aplicando el balance de energía mecánica (ecuación (45)) sobre los mismos puntos (1 y 2) se
obtiene:

p  p1
g ( z 2  z1 )  2  F  0 (51)


de donde:

F  g (z2  z1 )   p2  p1  (52)

18
OPERACIONES DE FLUJO DE FLUIDOS

Por tanto se deduce que:

F 1
 rz r 0  R  f V 2 (53)
2L 2

de donde:

L
F  2 fV2 (54)
D
que es la conocida ecuación de Fanning para una conducción cilíndrica, que relaciona las
pérdidas de energía mecánica con las condiciones de circulación, independientemente de la
inclinación de la tubería. Esta ecuación es completamente equivalente a la denominada
ecuación de Darcy-Weissbach, excepto que el factor de fricción de Darcy es igual a fD = 4f, y
por tanto:

V2 L
F  f D (55)
2 D

Este factor de fricción de Darcy está definido para que en la ecuación de pérdidas aparezca que
V 2 
las pérdidas son proporcionales a la energía cinética   , y suele ser más utilizado en
 2 
ingeniería hidráulica, mientras que en ingeniería química se utiliza el de Fanning.

La ecuación (49) es resoluble analíticamente en régimen laminar al conocerse las componentes

del tensor  en función de la ley de Newton de la viscosidad (ecuación (21)). En este caso, dada
la definición del factor de fricción, sustituyendo en la ecuación (25), se obtiene

16
f
Re

Para el régimen de circulación turbulento, el factor de fricción se puede determinar mediante

la gráfica de Moody (Figura 13) o por ecuaciones semiempíricas, en función del número de
Reynolds y la rugosidad de las conducciones, tales como la ecuación de Colebrook-White, que
no es más que el ajuste matemático de la gráfica de Moody a una ecuación:

1   1.256 
 4 log    (56)
f  3.7D Re f 

A partir de la ecuación (56), P. K. Swamee y A. K. Jain, desarrollaron una ecuación que permite el
cálculo directo del factor de fricción para flujo turbulento.

0.0625 (57)
f  2
  
  1 5.74  
log   0.9  
  3.7  D  Re  
    
   

19
OPERACIONES DE FLUJO DE FLUIDOS

La rugosidad  representa la altura media de los desniveles en la superficie y los valores que
aparecen en muchas tablas son para conducciones nuevas o limpias (es posible que con el
tiempo se incrementen los desniveles superficiales). En la Tabla 1 se muestran las rugosidades
de algunos materiales típicos de conducciones.

Tabla 1. Rugosidades de varios materiales

En la utilización de las gráficas de Moody y ecuaciones semiempíricas, hay que cerciorarse de si

expresan el factor de Fanning o el de Darcy.

Figura 13. Gráfica de Moody para la determinación del factor de fricción.

Es muy frecuente expresar la ecuación de Fanning en función del caudal volumétrico Q, en vez
Q
de en función de la velocidad. Sustituyendo V por :
S

20
OPERACIONES DE FLUJO DE FLUIDOS

2
 4Q  L 32fQ L
2
F  2f    (58)
 D  D  2 D5
2

Para la circulación de agua o disoluciones acuosas con densidades y viscosidades próximas a la

del agua a través de tuberías, se pueden utilizar las ecuaciones empíricas, que son de la forma

F Qn
J  m (59)
gL D

donde el parámetro J representa la pendiente hidráulica (pérdida de carga por rozamiento y

por unidad de longitud L de la conducción) en función del caudal Q (m3/s) y el diámetro D. Las
más usadas son las siguientes:

Q1.75
i) para régimen turbulento liso (de la ecuación de Blausius) J  0.008 (60)
D 4.75

ii) para régimen de transición

Q1.78
-para tubos de amianto-cemento J  0.00098 (61)
D 4.78

iii) para régimen de circulación turbulento

10.62 Q1.85
- fórmula de Hazen-Williams J (62)
C1.85 D 4.87

donde C depende del material de la tubería: 140 para tuberías rectas y muy lisas (plásticos),
135 para tuberías de amianto-cemento, 130 para hormigón armado, 120 para acero nuevo,
110 para acero usado, 100 para tuberías de fundición nuevas, 90 para tuberías de fundición
con algunos años de servicio y 80 en tuberías de fundición en malas condiciones. Puede
demostrarse que esta expresión es similar a la de Fanning.

-fórmula de Manning. Esta ecuación tiene una expresión general con independencia de la
geometría de la conducción, y el resultado de aplicárselo a una sección cilíndrica es:

Q2
J  10.3n 2 (63)
D 5.33

El coeficiente n toma valores constantes en función del material de construcción de la

conducción, por lo que se considera que es independiente del Reynolds: 0.08 para plástico liso
(PVC), 0.09 para amianto-cemento, 0.01 para hormigón armado, 0.012 para acero nuevo,
0.013-0.014 para hormigón en masa y 0.016-0.018 para plástico corrugado.

4.2 Pérdida de energía mecánica de accidentes en conducciones

En las conducciones cilíndricas rectas se ha visto que el fluido pierde energía mecánica por
rozamiento y fricción en las paredes de la conducción, cuyo valor viene dado por la ecuación
de Fanning. Estas pérdidas tienen lugar teniendo en cuenta que las líneas de flujo del fluido

21
OPERACIONES DE FLUJO DE FLUIDOS

siguen una trayectoria totalmente recta. En todos los sistemas de conducciones es preciso que
haya una serie de accesorios (accidentes), que hacen variar las líneas de flujo del fluido,
distorsionando la uniformidad de la corriente, aumentando la turbulencia donde están estos
accesorios. Esto implica que, necesariamente, se produzcan nuevas pérdidas de energía
mecánica, que serán mayores cuanto mayor sea la perturbación introducida en ese punto. Las
pérdidas de energía mecánica provocadas por estos accesorios se suelen denominar pérdidas
menores a causa de que su valor, para cada accidente, suele ser pequeño. Sin embargo, la
suma de estas pequeñas pérdidas puede adquirir importancia y suponer un elevado porcentaje
de la pérdida total.

Las pérdidas menores se suelen calcular por alguno de los dos procedimientos empíricos
siguientes:

a) Determinando la longitud de tubo recto de igual diámetro, que produciría una pérdida
de carga análoga a la del accidente. Esta longitud se denomina "longitud equivalente"
del accidente. Basta luego añadir esta longitud equivalente a la longitud de tramos
rectos en la ecuación de Fanning para calcular la pérdida de energía correspondiente al
sistema.
b) Si el régimen es turbulento, como es lo más frecuente, se puede calcular según la
ecuación
V2
F  K (64)
2

en la que K es una constante característica de cada accidente y V es la velocidad del fluido.

Para régimen laminar, Crane propone la siguiente relación empírica:

L Re  L 
     (65)
 D l 1000  D t

donde el subíndice “l” se refiere a la longitud equivalente (expresada en diámetros de tubo)

para régimen laminar, si Re < 1000 y el “t” a la longitud equivalente (expresada en diámetros
de tubo) para régimen turbulento.

En la Tabla 2 aparece una relación de valores de K para distintos accidentes, así como algunas
longitudes equivalentes, expresadas en número de diámetros.

22
OPERACIONES DE FLUJO DE FLUIDOS

Tabla 2. Valores de K de la ecuación (64) para distintos accidentes y algunas longitudes equivalentes.

Accidente K L.E.* Velocidad

 D 12 
2
2
1  2  1
 D2  medida en el
tubo estrecho 
Ensanchamiento
V1

2
1


Estrechamiento " "

Entradas

-encañonada 0.78 la del tubo

-cantos vivos 0.50 "


-ligeramente
redondeada
0.23 "

-bien redondeada 

0.23 "

* (número de diámetros)

23
OPERACIONES DE FLUJO DE FLUIDOS

Tabla 2: Continuación.

Accidente K L.E* Velocidad

Salidas

-encañonada 1.00 la del tubo

-cantos vivos 1.00 "

-redondeada 1.00 "

codo de 45º

-standard 0.35 la del tubo

-radio medio (33,

34, 39, 40, 41)
0.30 15.0 "

-gran curvatura (33)

0.20 10.0 "

codo de 90º

-standard 0.75 la del tubo

-radio medio (33,

34, 40, 41, 42, 43)
0.75 37.5 "

-gran curvatura (33,

34, 39, 40)
0.45 22.5 "

-pequeña curvatura 1.30 "

* (número de diámetros)

24
OPERACIONES DE FLUJO DE FLUIDOS

Tabla 2: Continuación.

Accidente K L.E* Velocidad

codo de 90º

-A escuadra 1.20 60.0 "

"U" de retorno (180º) (33, 1.50 75.0 la del tubo

34, 40)

Uniones y acoplamientos 0.04 2.0 la del tubo

(33, 41)

T standard

-con la bifurcación la del tubo

cerrada 
0.40

-usada como codo 

 ó
1.00 "



-con división del la de la corriente

caudal principal
1.00

-Entrada rama la del tubo

perpendicular (39)
1.80 90.0

-Salida rama 1.20 60.0

perpendicular (39)

unión roscada 0.04 la del tubo

manguito de unión 0.04 la del tubo

válvula de compuerta (39,

41, 44)

-abierta 0.17 8.5 la del tubo

-3/4 abierta 0.90 45.0 "

-1/2 abierta 4.50 225.0 "

-1/4 abierta 24.0 1200.0 "

* (número de diámetros)

25
OPERACIONES DE FLUJO DE FLUIDOS

Tabla 2: Continuación.

Accidente K L.E* Velocidad

válvula de diafragma(45)

-abierta 2.30 115.0 la del tubo

-3/4 abierta 2.60 130.0 "

-1/2 abierta 4.30 215.0 "

-1/4 abierta 21.0 1050.0 "

válvula de asiento(41,42)

cóncavo -abierta 6.00 300.0 la del tubo

-1/2 abierta 9.50 475.0 "

complejo -abierta 6.00 300.0 "

-1/2 abierta 8.50 425.0 "

plano -abierta 9.00 450.0 "

-3/4 abierta 13.0 650.0 "

-1/2 abierta 36.0 1800.0 "

-1/4 abierta 112.0 5600.0 "

válvula angular abierta la del tubo

(39, 41)
2.00 100.0

válvula en Y abierta (39,

44)
3.00 150.0 "

válvula de retención (39,

41, 44) (abierta)

-de charnela 2.00 100.0 "

-de disco 10.0 500.0 "

-de bola 70.0 3500.0 "

* (número de diámetros)

26
OPERACIONES DE FLUJO DE FLUIDOS

Tabla 2: Continuación.

Accidente K L.E* Velocidad

válvula de bola o macho

(29)

=5º
0.05 2.5
10º
0.29 14.5 la del tubo 
20º
1.56 78.0
40º
17.3 865.0
60º
206.0 10300

válvula de mariposa (29)

=5º 0.24 12.0

10º 0.52 26.0


20º 1.54 77.0 la del tubo

40º 10.8 540.0

60º 118.0 5900.0

válvula de fondo (41) 15.0 750.0 la del tubo

caudalímetros mecánicos

-de disco 7.00 7.00 la del tubo

-de pistón 15.0 15.0 "

-rotatorio 10.0 10.0 "

-de turbina 6.00 6.00 "

* (número de diámetros)

27