PRIMERA LEY DE LA TERMODINÁMICA: ECUACIÓN DE ENERGÍA

PROBLEMA TEÓRICO:
Considerando la situación ordinaria de un volumen de control para el cual los
flujos de salida y entrada son permanentes y unidimensionales sin pérdidas
debidas a la fricción. (Ver figura). Realizar un balance de energía.
dQ
dt

P1, V1

1, A1 P2, V2

 2 , A2
dw
dt

SOLUCIÓN: Estableciendo el volumen de control, tenemos:

Q
dt

P1, V1, 1, A1

1
P2, V2, A2, 2,

2
WS
dt

Aplicando la ecuación de balance de energía, al volumen de control

considerado, tenemos:
Q WS  P  W
    e   . (v.n)dA   e...dv  (1.1)
dt dt c. s
  t c.v dt
De acuerdo con las condiciones especificadas en el enunciado, la ecuación que
corresponde al balance de energía total, se transforma:
Q WS  P  W
   e   . (v.n)dA , donde,
t c.v
 e.  ..dv  0 , 0.
dt dt c . s
  dt
Si analizamos con atención la integral de superficie, veremos que el producto
.(v.n) .dA es la rapidez de flujo de masa y el signo de éste producto indica si el
flujo de masa ocurre hacia adentro o hacia fuera del volumen de control,
dependiendo del signo del producto escalar de (v.n) .

El factor por el cual está multiplicada la rapidez de masa, o sea e + P/,

representa los tipos de energía que pueden entrar o salir del volumen de
control por unidad de masa de fluido. La energía especifica total, e, se puede
considerarla como la expresión que incluya las contribuciones de energía:
cinética, potencial e interna, de tal manera que:
P v2 P
(e  )  gy  u
 2 
Ya que la masa únicamente entra al volumen de control en la sección (1) y sale
en la (2), la integral de superficie se transforma en,
 P  v22 P2   v12 P1 
C.S   
 e    .( v.n ).dA  
 2
 gy2  u 2 


2 
(  2 v2 A2 )  
 2
 gy1  u1 
1
 ( 1v1 A1 )

En un problema anterior se dedujo que en esta misma situación, para el

V1
1 V2
A1 2
A2

1
2

flujo permanente unidimensional una vez establecido el volumen de control, por

las líneas punteadas en la figura, se tiene:
Considerando el volumen de control indicado en la figura, aplicando la
expresión integral de la conservación de la masa a un flujo permanente,
tenemos:

 C .S
 (v.n )dA 
t C .V
 dv  0 (1-2)

Como la masa atraviesa la superficie de control solamente en las posiciones (1)

y (2) y es permanente, la expresión esta dada por:

C .S
 (v.n )dA    (v.n )dA    (v.n )dA  0 ,
A1 A2 t C .V
 dv  0
El valor absoluto del producto escalar (v.n ) es igual a la magnitud de la
velocidad en cada una de las integrales ya que los vectores velocidad así como
los vectores normales dirigidos hacia afuera, son colineales, tanto en (1) como
en (2). En (2) ambos vectores tienen el mismo sentido, por lo que este producto
es positivo, como debe ser para un flujo hacia afuera de masa. En (1), donde la
masa fluye hacia el volumen de control, ambos vectores tienen sentidos
opuestos, por lo que el signo es negativo. Ahora se puede expresar la ecuación
de continuidad en forma escalar.

 C .S
 (v.n)dA    vdA    vdA  0
A1 A2

La integración produce el resultado:

1v1 A1 = 2 v2 A2 (1-3)
que nos es familiar que corresponde al equilibrio de masa, esta ecuación puede
ser escrita de la siguiente forma:
1v1 A1  2 v2 A2  m
Sustituyendo, todos estos valores en la ecuación (1-1), la expresión que
corresponde al balance de energía, en éste problema, será.

 Q  ws  v2 P   v2 P
   2  gy2  u2  2  ( 2 v2 A2 )   2  gy1  u1  1  ( 1v1 A1 ) (1-4)
dt dt 2 2  2 1 

Si se divide cada uno de los términos de la expresión anterior entre la rapidez

de flujo de masa, se obtendrá la ecuación:
q  ws  v22 P   v2 P
   gy2  u2  2    1  gy1  u1  1  (1-5)
m 2 2   2 1 
Ó la ecuación más familiar:
q ws v22  v12
   g ( y2  y1 )  h2  h1 (1-6)
m m 2
P
Donde la suma de la energía interna y el flujo de energía, o sea u  ha sido

reemplazada por la entalpía h, que es igual a la suma de estas cantidades, por
P
definición: h = u  .


OBSERVACIÓN I.
Aplicando la ecuación (1-1) de balance de energía total a un volumen de
control (como el de la figura, en el cual el flujo sea permanente, incompresible y
no viscoso, en el cual no ocurra ninguna transferencia de calor ni cambio en la
energía interna) con las consideraciones adecuadas, obtendremos el resultado
siguiente:
Q  ws   w
dt
 0,
dt
 0, 
t C .V
e dv  0 ,
dt
0

P P P
C .S
 (e  )(v.n )dA    (e  )(v.n )dA    (e  )(v .n )dA  0
 A 2  A 1

 v22 P2   v12 P
=   gy2   ( 2 v2 A2 )    gy1  1  ( 1v1 A1 )  0
2 2  2 1 
Así, la expresión para la primera ley se transforma en:
 v22 P2   v12 P1 
  gy   2 2 2   gy1   ( 1v1 A1 )  0
(  v A ) 
2
2
2  2 1 

Como el flujo es permanente, el flujo de masa es, 1v1 A1  2 v2 A2 así, la

ecuación se simplifica, a:
v12 P1 v2 P
gy1    gy2  2  2 (1-7)
2 1 2 2
Dividiendo la ecuación (27-5) entre “g”, tenemos:
v12 P v2 P
y1   1  y2  2  2
2 g 1 g 2 g 2 g (1-8)

Cualquiera de las ecuaciones anteriores se llama ecuación de Bernoulli. Debe

tenerse en cuenta que cada uno de los términos de la ecuación (1-8), tiene la
unidad de longitud. Con frecuencia a las cantidades se les llama “carga” debido
a la elevación, a la velocidad y a la presión. Estos términos, tanto individual
como conjuntamente indican las cantidades que se pueden convertir
directamente para producir energía mecánica.

OBSERVACIÓN II.
Aplicando la ecuación (1-1) de balance de energía total al análisis del flujo en el
tubo. En este caso el volumen de control es el tubo, junto con el líquido que
contiene. Si se aplica a este volumen de control la ecuación (1-1) y se obtiene
la ecuación (1-5) la cual puede ser escrita de la siguiente manera:
q ws v22  v12 P P
   g ( y2  y1 )  2  1  u2  u1
m m 2 2 1
Si el sistema considerado no intercambia calor, q/m=0, el cambio de energía
interna puede escribirse: ghL, que es la pérdida de carga debida a la fricción,
así la ecuación puede ser escrita:
v12 P1 v2 P
y1    wP  y2  2  2  hL (1-9)
2g  2g 
Donde, ws/m=wp
P
A la expresión y  , se le conoce con el nombre de altura piezométrica.

v2
P
A la expresión y   , se le conoce con el nombre de energía total.
 2g
La ecuación (1-9) puede ser escrita y agrupada convenientemente, así:
 P1 v12   P2 v22 
 y1    y     wp  hL (1-10)
  2g   2  2g 
La expresión del lado izquierdo se llama energía disponible y la expresión del
lado derecho energía requerida.
El término wP es la energía absorbida por el fluido del medio ambiente (a través
de motor eléctrico, motor a combustión).
El término hL, es la suma de todas las pérdidas de carga por fricción e incluye a
las pérdidas locales (válvulas, conexiones, accesorios) y a la pérdida de carga
continua (tubería).
v2
La expresión matemática para la pérdida de carga local es: h L= k donde k
2g
es un factor de resistencia y su valor depende del tipo de conexión ó
accesorio.
La expresión matemática para la pérdida de carga continúa esta dado por la
ecuación de DARCY:
L v2
hL  f
D 2g
La ecuación de DARCY resulta del estudio del flujo de un fluido a través de
una tubería por el método del análisis dimensional.

PROBLEMA Nº 2
En una tubería horizontal de 4” de diámetro nominal nº 40 de cédula, circula
vapor de agua, se halla instalado suficientemente distanciados dos
termómetros y dos manómetros en los cuales se observan las siguientes
lecturas:
Sección 1: 103 psia, 330ºF
Sección 2: 67 psia, 300ºF
Considerando que la tubería esta totalmente aislada, determine la cantidad de
vapor que circula por la tubería.

1, v1,A 2, v2,A

T1, P1, h1 T2, P2, h2

4” Nº 40
c.v.

1 2

SOLUCIÓN:
- Aplicando la ecuación (1.2) de balance de masa al volumen de control
considerado se obtiene la ecuación (3):
1v1 A1  2 v2 A2   vA  m (3)

- Aplicando la ecuación (1.1) al volumen de control considerado

obtenemos la ecuación (1.4), así el balance de energía es:
q ws v22  v12
   g ( y2  y1 )  h2  h1 (1.6)
m m 2
Para el problema considerado se tiene: q/m=0; ws/m=0; y2=y1, analizando las
ecuaciones (1.2) y (1.6), las incógnitas P1, P2, h1, h2 y A (tubería Nº 40) pueden
ser obtenidas de tablas de vapor y de un manual, así tenemos:
Sección 1: 330ºF, 103 psia, 1188 BTU/lbm, 4.31 ft3/lbm
Sección 2: 300ºF, 67 psia, 1180 BTU/lbm, 6.47 ft3/lbm
Sustituyendo en la ecuación (28.1) y sabiendo que V=1/, se tiene:
q ws v22  v12
   g ( y2  y1 )  h2  h1
m m 2
Sustituyendo, este valor en la ecuación (28.2), obtenemos:
 4.31 2  2
1  ( 6.47 )  v
 778(1180  1188)  0
2(32.2)
Resolviendo, v2=853 ft/s, para A40=12,7 in2
La masa circulante por la tubería será:
v2 A 853 12.7
m  2 v2 A   x lb / s  11, 63 lbm / s
V2 6.47 144
PROBLEMA Nº 3
Una masa de aire a 500ºC y 5 atm. Se expansiona en una turbina hasta la
presión de una atm. El flujo de aire es de 50 kg/h a través de una tubería de
3/8”, y en la turbina se realiza un trabajo de 1 kw-h a la vez que pierde la
cantidad de calor de 700 kcal/h. Para la masa molecular media del aire puede
tomarse 29, y su calor especifico, prácticamente independiente de la
temperatura y la presión, es igual a 0.24 kcal/kg ºC. Calcúlese la temperatura y
la velocidad del aire a la salida de la tubería.
SOLUCIÓN
c.v
ws
dt
3/8”
w(aire)

Q 1 2

dt v1 V2
A A
P1 P2
1 2

- Aplicando la ecuación (1.1) de balance de masa al volumen de control

considerado, obtenemos la ecuación (1.3):
 1v1 A1  2 v2 A2  m ó 1v1  2 v2 (1.2)

- Aplicando la ecuación (1.1) de balance de energía al volumen de control

considerado, obtenemos la ecuación (1.6):
q ws v22  v12
    g ( y2  y1 )  h2  h1 (1.6)
m m 2
Analizando las ecuaciones (1.2) y (1.6) se concluye que se pueden calcular P1,
v1, q/m, ws/m, con los datos que se conocen, m,A,P 1,P2, M,R,T y con ayuda de
las relaciones adicionales para un gas ideal:
PM=RT, dQ=dh=cpdT
Podemos calcular, 2=f2( T2 ); v2=f1(T1) y h= f(T2), así:
h2 T2
1  PM
1 / RT1 , 2  P2 M / RT2 , h1
dh   cP dT , h2  h1  CP (T2  T1 )
T1

Y v1= m/1A= mRT1/AMP1, sustituyendo valores de 1 y 2 en (1.2) se tiene:

v2=1v1/2 = m/2A = mRT2/ AMP2 (1)
Sustituyendo valores en la ecuación (1.6) se obtiene:
(mRT2 / AMP2 )2  (mRT1 / AMP1 ) 2
(q / m  ws / m)   Cp(T2  T1 ) (2)
2
en la ecuación (2) la única incógnita es T2.
Para hallar la solución de la ecuación (2), calculamos los valores conocidos,
así:
q/m = (700kcal/h)/(50 kg/h)=14 kcal/kgm
ws/m=(1kw-h/h)(860kcal/1kw-h)/(50kg/h) = 17.2 kcal/kgm.
mRT1/AMP1 = (50x0.082x773)/(3600x1.23x10-4x29x5)=49.36(m/s)
mRT2/AMP2 = (50x0.082xT2)/(3600x1.23x10-4x29x1)=0.32 T2(m/s)
Cp(T2-T1)=0.24(T2-773) kcal/kg.
La ecuación (2) no es homogénea respecto a las unidades, por tanto el primer
termino del lado derecho debemos introducir el valor de la constante g y las
unidades mecánicas convertirlas en térmicas, así:
(0.32T2 ) 2  (49.36) 2 1
 (14  17.2)  x  0.24(T2  773)(kcal / kg)
2(9.81) 427
Resolviendo la ecuación se obtiene: T2 = 624ºK
Con esta temperatura calculamos la velocidad de salida del aire, así:
v2  0.32T2  0.32(624)  200 m / s
UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU
FACULTAD DE INGENIERIA QUIMICA
TRABAJO INDIVIDUAL Y/O GRUPAL DE LA ASIGNATURA DE FLUJO DE FLUIDOS
APELLIDOS Y NOMBRES:………………………………..CICLO:…SECCION:…
INDICACIONES GENERALES:
- El orden, la limpieza, la claridad en las explicaciones así como los resultados numéricos se tomaran en
cuenta en la calificación.
1.- Se necesita transportar 50 m3/h de etanol desde un depósito situado en la planta baja de una fabrica, hasta un reactor situado 20 m
sobre el depósito (en sentido vertical). La conducción se ha de efectuar a través de tubería de 4”, y la instalación tiene una longitud
de 45 m con 4 codos cerrados y tres válvulas de asiento. Calcúlese:
a) La potencia de la bomba a instalar si el rendimiento del grupo motor-bomba es del 65 %.
b) El coste de bombeo si el kilovatio-hora cuesta 0,40 nuevo sol.
DATOS: Densidad= 789 kg/m3 , Viscosidad = 1,194x10-3 kg/m.s (Rpta. 5,33 CV ; 1,57 ptas/h).
2.- Encuéntrese la descarga en el sistema de la figura para las bombas A y B en
paralelo o en serie con una longitud de 1800 m de tubería de hierro forjado nuevo,
D= 300 mm y H = 25 m y determínese la potencia necesaria para operar las
bombas.
BOMBA A BOMBA B
Hp (m) Qp ( l/s) e% Hp (m) Qp ( l/s) e %
21,3 0 0 24,4 0 0 H
18,3 56,6 59 21,3 74 54 ,
H 2O 200 C
16,8 72,5 70 18,3 112 70
15,2 85,8 76 15,2 140 80 P
13,7 97,7 78 12,2 161 73
12,2 108,0 76,3 9,1 174 60
10,7 116,0 72 6,1 177 40
9,1 127,0 65
7,6 130,0 56,5
6,1 134,0 42

3.- Determinar los gastos y las perdidas de carga en cada uno de los ramales del sistema de tuberías mostrado en la figura.

400 l/s (1) 14”

A 1000 m
C=120

16” C=100

(2)
800 m
(5) D 10” C=120
(4)
12”
1800 m 500 m
C= 120 1600 m
8” (3) B
Figura : 1 C=130

4.-Si el sistema de tuberías mostrado en la figura 1, se dispone en serie ¿Cuál seria la mejor distribución?

5.- Si el sistema de la figura 1, se dispusiera en paralelo, ¿Cuál seria el gasto en cada una de las tuberías?.

6.- Si en el sistema de tuberías mostrado en la figura 1, todas las tuberías tienen C= 120. ¿Qué diferencias existen con respecto al
sistema inicial?

Factores de conversión: 1HP=33 000 ft-lb/min = 76,04 kg.m/s ; 1CV = 75 kg.m/s = 0,736 Kw ; 1 CVh = 270 000 kg.m
1 Kwh = 3,6x106 julios = 860 Kcal.
Ing.P.V.G.Y
1.- SOLUCION:

Q wS  P  w
    e  . (v .n )dA   e  dv  27.1)
dt dt C .S
  t C .V dt
De acuerdo con las condiciones especificadas en el
enunciado, la ecuación que corresponde al balance de
energía total, se transforma:
Q wS  P  w
    e  . (v.n )dA, donde,  e  dv  0,  0.
dt dt C .S
  t C .V dt
Si analizamos con atención la integral de superficie,
veremos que el producto  (v.n)dA es la rapidez de flujo de
masa y el sigo de éste producto indica si el flujo de masa
y el signo de éste producto indica si el flujo de masa
ocurre hacia adentro o hacia fuera del volumen de control,
dependiendo del sentido de (v.n ) .
El factor por el cual está multiplicada la rapidez de masa,
o sea e + P/, representa los tipos de energía que pueden
entrar o salir del total, e , se puede considerarla como la
expresión que incluya las contribuciones de energía:
cinética, potencial e interna, de tal manera que:
P v2 P
(e  )  gy  u
 2 
Ya que la masa únicamente entra al volumen de control en la
sección (1) y sale en la (2), la integral de superficie se
transforma en,

 P  v22 P2   v12 P1 
C.S   
e  . ( v .n ) dA  
2
 gy 2  u 2   2 2 2   gy1  u1   ( 1v1 A1 )
2 
(  v A ) 
2 1 

En el problema (18) encontramos que en esta misma

situación, para el equilibrio de masa, se obtuvo la
ecuación (18.2) cuya ecuación,
1v1 A1  2 v2 A2  m
Sustituyendo, todos restos valores en la ecuación (27.1),
la expresión que corresponde al balance de energía, en éste
problema, será.

 Q  ws  v2 P   v2 P
   2  gy2  u2  2  ( 2 v2 A2 )   1  gy1  u1  1  ( 1v1 A1 ) (27.2)
dt dt 2 2  2 1 

Si se divide cada uno de los términos de la expresión

anterior entre la rapidez de flujo de masa, se obtendrá la
ecuación:

q  ws  v22 P2   v12 P
   gy2  u2      gy1  u1  1  (27.3)
m 2 2   2 1 

q ws v22  v12 P P u u q w
  y2  y1   2  1  2 1 ,donde,  0 , s  Hm
mg mg 2g 2 g 1 g g mg mg
Así, tenemos:
v2 P v2 P L v2
y1  1  1  ( H m )  y2  2  2  hL , donde, u2  u1  hL g , hL  f
2g  2g  D 2g
Aplicando a las condiciones del problema se tiene:
(o aplicando la ecuación de Bernoulli entre los puntos (1) y ( 2) )

l v2
0 + 0 + 0 + (-Hm) = H + 0 + 0 + f (1)
D 2g
 D2 4Q 8l
Se tiene: Q = v , v , hL  f 2 5 Q 2
4 D 2
 D g

Sustituyendo en (1) :

8l 8l
Hm  H  f Q 2 ; c1  2 5
 Dg
2 5
 Dg
 H m  H  c1Q2

DATOS:  = 784 kg ;   1,194 x103 kg / m3 ; LT= 45 m ; DN = 4”

m3
8l
a) Cálculo de hL  f Q2
 D g 5 2

1) Determinación de la longitud equivalente (Le) de los accesorios.

Accesorio Nº accesorios Le de c/u (m) TOTAL (m)

Codo cerrado 04 7,7 30,8
Válvula de asiento 03 40,0 120,0
TOTAL (Le longitud equivalente de los accesorios) 150,8 (m)
En consecuencia la longitud total será: LTotal = LT + Le = 45 + 150,8 = 195,8 m

2) Determinación de “f”
 
2.1) Determinación de la rugosidad relativa   :
D
 
Del grafico:   = 0,00043 , para DN = 4”, acero comercial.
D
Calculando: Di = 4,026” = 10,226 cm , # cédula = 40 ;  = 0,0046 cm.
 0, 0046
  0, 00045
D 10, 226
2.2) Determinación del número de Reynolds, (Re)

Dv  4 1 72  784 
Re    1,14 x105
   0,10226  1,194 x10 3

del Diagrama de Moody, se tiene : fD = 0,0205
asi tenemos:
8l
 H m  H  f 2 5 Q2
 D g
8 195,8 (1 72) 2
wP   H m  20 m  0, 0205 m
 2  0,10226   9,80 
5

-Hm = 20 + 5,73 = 25,73 m

3) Potencia:
wP   H m  m  .  kg m3  .Q  m3 s  .
1C.V
75 kg.m s
 Hm   Q
wP  C.V ;   Eficiencia
75 
Para el caso:
 25, 73 784 1 72 
wP  C.V  5, 78 C.V
75  0, 65
4) Cálculo del coste/h :
0, 736 Kw 1h S / .0, 40
Coste/h = 5,783 C.V x x x = S/. 1,703 /h
1C.V 1h 1Kw  h
Considerando que la Planta trabaje 8000 horas al año, se tiene:
S / .1, 703 8000 h
Coste/año = x = S/. 13 620,4 / año
h Año

Coste / año = S/. 13 620,4 /año .