ECUACIONES FUNDAMENTALES DE UN FLUJO

José Agüera Soriano 2012 1

ECUACIONES FUNDAMENTALES DE UN FLUJO

• ECUACIÓN DE CONTINUIDAD
• ECUACIÓN DE LA ENERGÍA
• ECUACIÓN CANTIDAD DE MOVIMIENTO
• APLICACIONES ECUACIÓN DE LA ENERGÍA
• APLICACIONES ECUACIÓN DE LA CANTIDAD
DE MOVIMINTO

José Agüera Soriano 2012 2

INTRODUCCIÓN
Son tres las ecuaciones fundamentales de un flujo:

• ecuación de continuidad (conservación de la masa)

• ecuación de la energía (conservación de la energía)
• ecuación de la cantidad de movimiento
(conservación de la cantidad de movimiento).

José Agüera Soriano 2012 3

ECUACIÓN DE CONTINUIDAD
En un flujo permanente,
m 1  m 2  m
V2
1  Q1   2  Q2    Q S2
2

V1

S1
volumen
de control
1

José Agüera Soriano 2012 4

ECUACIÓN DE CONTINUIDAD
En un flujo permanente,
m 1  m 2  m
V2
1  Q1   2  Q2    Q S2
2
a) para gases,
m  1  V1  S1   2  V2  S 2
V1

S1
volumen
de control
1

José Agüera Soriano 2012 5

ECUACIÓN DE CONTINUIDAD
En un flujo permanente,
m 1  m 2  m
V2
1  Q1   2  Q2    Q S2
2
a) para gases,
m  1  V1  S1   2  V2  S 2
V1
b) para líquidos ( 1   2   )
S1
volumen
Q  V1  S1  V2  S 2 de control
1

José Agüera Soriano 2012 6

Energía de un líquido en reposo (recordatorio)
p1 p2
E  g  z1   g  z2
 
La energía de un líquido en reposo viene dada por la suma
de dos términos: energía de presión y energía de posición.

José Agüera Soriano 2012 7

Energía de un líquido en reposo (recordatorio)
p1 p2
E  g  z1   g  z2
 
La energía de un líquido en reposo viene dada por la suma
de dos términos: energía de presión y energía de posición.

En metros de columna de líquido,

p1 p2
H  z1   z2
 

José Agüera Soriano 2012 8

Energía de un líquido en movimiento
A las dos energías anteriores hay que sumar la energía
cinética del flujo:

V2 p
E   g  z J/kg en S.I.
2 

José Agüera Soriano 2012 9

Energía de un líquido en movimiento
A las dos energías anteriores hay que sumar la energía
cinética del flujo:

V2 p
E   g  z J/kg en S.I.
2 

En metros de columna de líquido,

V2 p
H  z m
2g 

José Agüera Soriano 2012 10

Ecuación energía en conducciones de líquidos

H1 = H2 + Hr12

SLL plano de carga inicial (PC)

V 21 / 2g
A'
p 1 / línea pie H r 12
zométrica
(LP) línea d
1 e energí
a (LE)
A
LP V 22 / 2g
B'
H1

z1 p 2 /
B H2
2

V z2
plano de referencia

José Agüera Soriano 2012 11

Ecuación energía en conducciones de líquidos

V12 p1 V22 p2
H1 = H2 + Hr12   z1    z 2  H r12
2g  2g 

SLL plano de carga inicial (PC)

V 21 / 2g
A'
p 1 / línea pie H r 12
zométrica
(LP) línea d
1 e energí
a (LE)
A
LP V 22 / 2g
B'
H1

z1 p 2 /
B H2
2

V z2
plano de referencia

José Agüera Soriano 2012 12

V 2 2 g = altura cinética (desnivel entre LE y LP)
p  = altura de presión (desnivel entre LP y conducto)
z = altura de posición respecto plano de referencia

SLL plano de carga inicial (PC)

V 21 / 2g
A'
p 1 / línea pie H r 12
zométrica
(LP) línea d
1 e energí
a (LE)
A
LP V 22 / 2g
B'
H1

z1 p 2 /
B H2
2

V z2
plano de referencia

José Agüera Soriano 2012 13

V 2 2 g = altura cinética (desnivel entre LE y LP)
p  = altura de presión (desnivel entre LP y conducto)
z = altura de posición respecto plano de referencia
p   z = altura piezométrica

SLL plano de carga inicial (PC)

V 21 / 2g
A'
p 1 / línea pie H r 12
zométrica
(LP) línea d
1 e energí
a (LE)
A
LP V 22 / 2g
B'
H1

z1 p 2 /
B H2
2

V z2
plano de referencia

José Agüera Soriano 2012 14

Ecuación de la energía en máquinas de flujo líquido
El trabajo (cedido o recibido) que atraviesa los límites de un
sistema abierto se llama genéricamente, trabajo técnico, Wt
Cuando de turbomáquinas se trata, como ocurre en la figura,
es más intuitivo llamarlo trabajo interior en el eje.

turbina

Wt

José Agüera Soriano 2012 15

Ecuación energía en máquinas de flujo líquido
H1 = H2 + |Hr12| + Ht

H r 12 H r 12

Ht
Ht
Ht Ht
H1
H2
turbina bomba
1 2
H2 H1
z1
2 z1 1
plano de referencia plano de referencia

José Agüera Soriano 2012 16

Ecuación energía en máquinas de flujo líquido
H1 = H2 + |Hr12| + Ht
V12  V22 p1  p2
Ht    z1  z 2  H r
2g 

H r 12 H r 12

Ht
Ht
Ht Ht
H1
H2
turbina bomba
1 2
H2 H1
z1
2 z1 1
plano de referencia plano de referencia

José Agüera Soriano 2012 17

Ecuación energía en máquinas de flujo líquido
H1 = H2 + |Hr12| + Ht
V12  V22 p1  p2
Ht    z1  z 2  H r
2g 

H r 12 H r 12

Ht
Ht
Ht Ht
H1
H2
turbina bomba
1 2
H2 H1
z1
2 z1 1
plano de referencia plano de referencia

José Agüera Soriano 2012 18

Representación gráfica
H1  H2 = |Hr12| + Ht
He
He

Hm
Ht Hm
Ht

Ht
Ht
Hr
TURBINA
H2 H1 Hr
H1 H2
BOMBA

José Agüera Soriano 2012 19

Representación gráfica
H1  H2 = |Hr12| + Ht
He
He

Hm
Ht Hm
Ht

Ht
Ht
Hr
TURBINA
H2 H1 Hr
H1 H2
BOMBA

José Agüera Soriano 2012 20

Potencia de un flujo
Turbina de reacción
Caudal : Qm s 
3

3
  Q kg/s (caudal másico)
Densidad :  kg m 

José Agüera Soriano 2012 21

Potencia de un flujo
Turbina de reacción
Caudal : Qm s 
3

3
  Q kg/s (caudal másico)
Densidad :  kg m 
Altura : Hm 
2
g  H m 2
s 2
(J/kg)
Gravedad : g m s 

José Agüera Soriano 2012 22

Potencia de un flujo
Turbina de reacción
Caudal : Qm s 
3

3
  Q kg/s (caudal másico)
Densidad :  kg m 
Altura : Hm 
2
g  H m 2
s 2
(J/kg)
Gravedad : g m s 

P    g  Q  H J/s (W)

José Agüera Soriano 2012 23

EJERCICIO
3
La energía de un caudal de agua de 60 m /s, es
H = 50 m. Calcúlese su potencia.
Solución
P    Q  H  1000  9,81 60  50 
 29430 103 W  29430 kW (40000 CV)

José Agüera Soriano 2012 24

 ECUACIÓN DE LA CANTIDAD DE MOVIMIENTO
Cuando a lo largo de un volumen de control, la velocidad
del flujo varía, es porque actúan fuerzas sobre él que lo
aceleran:

F  m  a
B
B
B V1 C
B C 
V1 F V2 V1
A C V2
C 
V1 F V2 D
F
D
A V2
A D
D
F

A
José Agüera Soriano 2012 25
 B
B B' V1
C
C 
V1 F V2
A V2
m1 D
F
D

V1 
A
(a) C (b)
F B B'
C'
m2
A A' m1
V1 
D V2 F
C
C'
m2
D' A A'

El impulso (F  dt ) sobre la masa del volumen de control

D V2
(c) D'

provocará una variación de su cantidad de movimiento [d ( m  V )]:


F  dt  d (m  V )  dp

José Agüera Soriano 2012 26

 B
B B' V1
C
C 
V1 F V2
A V2
m1 D
F
D

V1 
A
(a) C (b)
F B B'
C'
m2
A A' m1
V1 
D V2 F
C
C'
m2
D' A A'

El impulso (F  dt ) sobre la masa del volumen de control

D V2
(c) D'

provocará una variación de su cantidad de movimiento [d ( m  V )]:


F  dt  d (m  V )  dp

Esta variación dp del sistema es la corresponde al instante
(t + dt), menos la que tenía en el instante t:
  
dp  pA'B'C'D'  pABCD 
   
 ( pA'B'CD  pCDD'C' )  ( pABB'A'  pA'B'CD )
José Agüera Soriano 2012 27
 B
B B' V1
C
C 
V1 F V2
A V2
m1 D
F
D

V1 
A
(a) C (b)
F B B'
C'
m2
A A' m1
V1 
D V2 F
C
C'
m2
Por ser el régimen permanente D' A A'
V2
D
   (c) D'

F  dt  dp  pCDD'C'  pABB'A' 
 m2 V 2  m1 V 1  m 2  dt V 2  m 1  dt V 1

José Agüera Soriano 2012 28

 B
B B' V1
C
C 
V1 F V2
A V2
m1 D
F
D

V1 
A
(a) C (b)
F B B'
C'
m2
A A' m1
V1 
D V2 F
C
C'
m2
Por ser el régimen permanente D' A A'
V2
D
   (c) D'

F  dt  dp  pCDD'C'  pABB'A' 
 m2 V 2  m1 V 1  m 2  dt V 2  m 1  dt V 1
F  m 2  V 2  m 1  V 1

F  m  (V 2  V 1 )
válida para líquidos y para gases

José Agüera Soriano 2012 29

APLICACIONES ECUACIÓN DE LA ENERGÍA
Salida por un orificio
pa
SLL

p V2
 2g =
H

José Agüera Soriano 2012 30

APLICACIONES ECUACIÓN DE LA ENERGÍA
Salida por un orificio
pa
SLL

V2
H
p V2 2g
 2g =
H

V  2 g  H

José Agüera Soriano 2012 31

Flujo en tuberías con salida libre
Con el mismo diámetro y el mismo desnivel entre el
extremo 2 y la SLL, se cumple para cualquier longitud,
V2
H  Hr 
2g

José Agüera Soriano 2012 32

 V2
H  Hr 
2g

José Agüera Soriano 2012 33

 V2
H  Hr 
2g
Mayor longitud L de la tubería origina:
- más pérdida de carga Hr
- menos velocidad V del flujo en la tubería
- menos caudal Q
- menos pérdida de carga unitaria J (Hr/L)
José Agüera Soriano 2012 34
Salida por tobera
Se trata de colocar al final de la conducción una reducción de
sección (tobera).
Cuanto más se reduzca la sección S de salida, menor será el
V2
2g
LE
LP
LE

caudal Q que circula, menor la velocidad V y menor Hr V 2i

LP
2g
pB V S2
H
 2g =
pi

V
S VS
B

SLL A' plano de carga inicial

línea piezom Hr
A étrica con
tobera B'
líne L
a pi
ezo V
mé H1
tric pB V S2
as in t
obe  = H = 2g
ra (
Q
má
x)

S
B
VS

línea piezométrica entre

Q = 0 y Q = Qmáx
tobera
José Agüera Soriano 2012 35
En B, la energía H está prácticamente en forma de presión
(pB/). A lo largo de la tobera la velocidad aumenta y la presión
disminuye. A la salida S, pS = pa

LE
V2 LE
2g LP
V 2i

LP
2g
pB V S2
H
 2g =
pi

V
S VS
B

una tobera es un transformador de energía potencial

(energía de presión en líquidos) en energía cinética.
José Agüera Soriano 2012 36
EJERCICIO
Calcúlese la sección de salida de la tobera al final de una
conducción, de 10 km y 1 m de diámetro, origine la pérdida
de 40 m en un salto de 400 m. Tómese,
L V2
H r  0,02  
D 2g
Solución
Velocidad en la conducción
2 2
L V V Hr
H r  0,02   ;   0,005  H r
D 2 g 2 g 0,02 10000

V 2 2 g  0,2 m ; V  1,981 m s

José Agüera Soriano 2012 37

Caudal
Q    r 2  V    0,5 2  1,981  1,556 m 3 s

José Agüera Soriano 2012 38

Caudal
Q    r 2  V    0,5 2  1,981  1,556 m 3 s

Velocidad de salida
H  VS2 2 g  360 m; VS  84,04 m s

José Agüera Soriano 2012 39

Caudal
Q    r 2  V    0,5 2  1,981  1,556 m 3 s

Velocidad de salida
H  VS2 2 g  360 m; VS  84,04 m s

Diámetro de salida
SS  Q VS  1,556 84,04; DS  0,154 m

El diámetro de la tobera pasa de 1 m a 0,15 m

José Agüera Soriano 2012 40

 Conducción hidroeléctrica de Villarino

L  15000 m
H  402 m
D  7,5 m; H r  40 m

De haber sido :
D  7,0 m; H r  60 m
D  8,0 m; H r  30 m

José Agüera Soriano 2012 41

Tubo de Pitot
2
v p
h 
2g 

Mide las alturas de presión y de velocidad en el punto M

José Agüera Soriano 2012 42

Fuerza sobre un conducto corto

F  F p  Fr
(Fr insignificante) B
Fr F

V p 1· S 1
Fp
F
C G
A
1 1
p
D 2·S
2
2 2

José Agüera Soriano 2012 43

Fuerza sobre un conducto corto
Valoración indirecta de la fuerza F
Las fuerzas sobre el volumen de control entre 1 y 2 son:
p1  S1  fuerza sobre la sección 1
B
Fr F

V p 1· S 1
Fp
F

p 2  S 2  fuerza sobre la sección 2

C G
A
1 1
p
D 2·S
2
2 2

 F  fuerza que ejerce la pared

G  fuerza de gravedad

José Agüera Soriano 2012 44

Fuerza sobre un conducto corto
Valoración indirecta de la fuerza F
Las fuerzas sobre el volumen de control entre 1 y 2 son:
p1  S1  fuerza sobre la sección 1
B
Fr F

V p 1· S 1
Fp
F

p 2  S 2  fuerza sobre la sección 2

C G
A
1 1
p
D 2·S
2
2 2

 F  fuerza que ejerce la pared

G  fuerza de gravedad

F  p1  S1  p 2  S 2  ( F )  G

José Agüera Soriano 2012 45

 F  p1  S1  p 2  S 2  ( F )  G
Por otra parte,
F  m  (V 2  V 1 )

José Agüera Soriano 2012 46

 F  p1  S1  p 2  S 2  ( F )  G
Por otra parte,
F  m  (V 2  V 1 )
Igualando y despejando F

F  p1  S1  p 2  S 2  G  m  (V 2  V 1 )

José Agüera Soriano 2012 47

 F  p1  S1  p 2  S 2  ( F )  G
Por otra parte,
F  m  (V 2  V 1 )
Igualando y despejando F

F  p1  S1  p 2  S 2  G  m  (V 2  V 1 )

En conductos cortos G  0, en cuyo caso,

F  p1  S1  p 2  S 2  m  (V 2  V 1 )

José Agüera Soriano 2012 48

 F  p1  S1  p 2  S 2  ( F )  G
Por otra parte,
F  m  (V 2  V 1 )
Igualando y despejando F

F  p1  S1  p 2  S 2  G  m  (V 2  V 1 )

En conductos cortos G  0, en cuyo caso,

F  p1  S1  p 2  S 2  m  (V 2  V 1 )

Cuando no hay flujo,

F  p1  S1  p 2  S 2

José Agüera Soriano 2012 49

Generalmente la fuerza será mayor cuando no hay flujo.

Las presiones a sustituir en las fórmulas anteriores son la

diferencia entre la interior y la que actúa exteriormente
sobre el conducto. Si ésta fuera mayor, los términos de
presión cambiarían de signo y con ello la fuerza F.

Cuando la presión exterior es la atmosférica, las presiones

a sustituir son las relativas.

José Agüera Soriano 2012 50

EJERCICIO
Calcúlese la fuerza sobre el codo de una conducción
hidráulica, de 90o y de 2,8 m de diámetro,
a) cuando fluyen 38,75 m3/s y la presión manométrica
es de 8 bar,
b) cuando el flujo se anula bruscamente, momento
en el que se prevé que la presión en el codo se eleve
a 13 bar a causa de un golpe de ariete.
Solución a)
Velocidad media
Q 38,75
V   6,29 m/s
S  1,4 2

José Agüera Soriano 2012 51

Fuerza sobre el codo
F  p1  S1  p 2  S 2  m  (V 2  V 1 )

Fx  p1  S1  0    Q  V1 
 8  10 5  π  1,4 2  1000  38,75  6,29 
 4926029  243737 
 5169766 N  5169,8 kN

José Agüera Soriano 2012 52

Fuerza sobre el codo
F  p1  S1  p 2  S 2  m  (V 2  V 1 )

Fx  p1  S1  0    Q  V1 
 8  10 5  π  1,4 2  1000  38,75  6,29 
 4926029  243737 
 5169766 N  5169,8 kN

Fy  0  p2  S 2    Q  V2 
 8 105   1,42  1000  38,75  6,29 
 4926029  243737  5169766 N 
 5169,8 kN

José Agüera Soriano 2012 53

Solución b)
F  p1  S1  p 2  S 2

Fx  p1  S1  0  13 105   1,42 
 8004797 N  8004,8 kN

José Agüera Soriano 2012 54

Solución b)
F  p1  S1  p 2  S 2

Fx  p1  S1  0  13 105   1,42 
 8004797 N  8004,8 kN

Fy  0  p2  S 2  13 105   1,4 2 
 8004797 N  8004,8 kN

José Agüera Soriano 2012 55

Solución b)
F  p1  S1  p 2  S 2

Fx  p1  S1  0  13 105   1,42 
 8004797 N  8004,8 kN

Fy  0  p2  S 2  13 105   1,4 2 
 8004797 N  8004,8 kN

El anclaje del codo hay que calcularlo para esta última

situación más desfavorable.
José Agüera Soriano 2012 56
José Agüera Soriano 2012 57
 Figuras no incluidas en las diapositivas
y v2
Fr 2
p1· S
1
h
v1 v
G z1 z2 dS
1
Fr
p 2· S 2 S
v2 V z1 v1

 z2 x
V
Figura 3-8
plano de referencia
Figura 3-5
 Fx
M 1

2
M v V1 aguja V2
dS
S h
x

m
Figura 3-6 Ejercicio 3-6.1
Figura 3-14
José Agüera Soriano 2012 58
 SLL

SLL

V2
V1 p 2· S2 h2
p 1· S 1 G
h1 Fr  0
1 Fp 2
x
Figura 3-17

u
p' p'
p'
p' dif Pt _0
p 2· S2 >
uso c1 = 0 bomba
r tobera
u F w2 F
w1 A tobera c2
u c2 2
c1 = 0 1
u = w1 2
p'
p'
x
p' p'

Ejercicio 3-7.2
Figura 3-18

José Agüera Soriano 2012 59

 p' 1' 2' p'
1 p'
u
1' 2'
p'
p'
2
p1' p2'
p'
w1
volumen de control
w2
·D 2
F S=
u
D
u c2 4
w
u = w1 p'
c1 = 0 u c D
p'
p'

p'

w, p

p= p p2'
w2 p' 1' 2' p'
(x )
p1 = p' w p 2 = p'
w1
w = w (x ) p1' p1' p2'
1 1' 2' 2 x

Figura 3-19 Figura 3-20

José Agüera Soriano 2012 60

 F

A B C D
Figura 3-21

turbocompresor cámara de combustión turbina tobera

1 1' 2
bulbo de
difusor 2' regulación
3 4
Figura 3-22

José Agüera Soriano 2012 61

 1 2 3 4 5 6 7

Figura 3-23

aire
aire
aire

gases de
combustión
aire
aire
aire

Figura 3-24

José Agüera Soriano 2012 62

 oxidante reductor tobera

Figura 3-25

p
1

0,6

1 2 3 u/ w2
Figura 3-26

José Agüera Soriano 2012 63

 u/ w2

m
0,3 m

1
0,
B 0,15 m
A
2
0,5 m
0,4 m
0,75 m
A
B 1
m

0,1 m
2

h
0,

0,2 m
0,2 m 0,3 m
Hg
Problema 3.3 Hg
Hg
Problema 3.4 Problema 3.5
B

B
h =1,5 m

m
A 15 m

5
,0
=0
D
A
4,2 m

C
Problema 3.7
Problema 3.6
José Agüera Soriano 2012 64
 D 60 m
2

3 15 m A
6m
3m
B C
1 3m
bomba bomba Problema 3.10
Problema 3.8

1
1m
4 SLL (mínima)
Pt 2
turbina
D 1=1,25 m 1 3

0,9 m Problema 3..12

2 A
D 2=1,5 m

h = 6,6 m
40 m
B
tubo de
descarga
SLL 3 C

Problema 3.13
Problema 3.11
José Agüera Soriano 2012 65
 y 1
1 2
V1 F
V1 V2 45º
2 V2 2 difusor tobera
bomba
1
Problema 3.15 x Problema 3.18
Problema 3.17

1
1 2 1 1
2 2 2
V

(a) placa o válvula (b) tobera en régimen (c) tobera en régimen Problema 3.20
cerrada de máxima caudal de caudal nulo

Problema 3.19

José Agüera Soriano 2012 66

 2

3 y V2
30º

1 V1
y 60º
x

x
Problema 3.22 Problema 3.24

x
2
V2
y
1 60º V1
F = Fy F

3 x
V2
Problema 3.25
Problema 3.26
José Agüera Soriano 2012 67