UNIVERSIDAD NACIONAL DEL CENTRO DEL PERÚ

FACULTAD DE INGENIERÍA QUÍMICA

BALANCE
MACROSCÓPICO
FOLLETO

Ing. Dr. ORLANDO A. VILCA MORENO

2018
 2

CONTENIDO

INTRODUCCIÓN
I. BALANCE MACROSCÓPICO EN SISTEMAS ISOTÉRMICOS

1.1. BALANCE MACROSCÓPICO DE MATERIA

MÉTODOS DE CÁLCULO

EJEMPLO DE APLICACIÓN

1.2. BALANCE MACROSCÓPICO DE CANTIDAD DE MOVIMIENTO

1.3. BALANCE MACROSCÓPICO DE ENERGÍA MECÁNICA (ECUACIÓN DE

BERNOULLI)
II. BALANCE MACROSCÓPICO DE ENERGÍA EN SISTEMAS NO ISOTÉRMICOS

2.1. BALANCE MACROSCÓPICO DE ENERGÍA

2.2. E JEMPLO DE APLICACIÓN

III. BALANCE MACROSCÓPICO EN SISTEMAS DE VARIOS COMPONENTES

3.1. BALANCE MACROSCÓPICOS DE MATERIA

3.2. BALANCE MACROSCÓPICOS DE CANTIDAD DE MOVIMIENTO

3.3. BALANCE MACROSCÓPICOS DE ENERGÍA

3.4. BALANCE MACROSCÓPICOS DE ENERGÍA MECÁNICA

BIBLIOGRAFÍA

2
 UNCP – FIQ Ing. Dr. O. Vilca

INTRODUCCIÓN

Todas las operaciones básicas, físicas y químicas de la ingeniería química implica el

transporte de una o varias de las tres propiedades fundamentales de cantidad de
movimiento, energía y materia. Los tres fenómenos de transporte se producen casi siempre
simultáneamente en todas las operaciones básicas, tanto físicas y químicas, pero la
importancia relativa de los mismos varía en cada caso, siendo frecuente que sólo uno de
ellos, por su mayor lentitud, determine el tamaño de equipo necesario para el desarrollo de
dichas operaciones.

Los balances macroscópicos son muy utilizados en el análisis de sistemas ingenieriles de

flujo. Los balances se aplican descartando los términos que resultan despreciables en un
determinado problema. Para saber que términos pueden despreciarse se requiere cierta
intuición y en algunos casos se necesitan observaciones experimentales acerca del
comportamiento del flujo.

3
FENÓMENOS DE TRANSPORTE – BALANCE MACROSCÓPICO
 4

BALANCE MACROSCÓPICO EN SISTEMAS ISOTÉRMICOS

BALANCE MACROSCÓPICO DE MATERIA

Las ecuaciones correspondientes al balance macroscópico de materia constituyen una de

las herramientas matemáticas más útiles de la Ingeniería Química, donde su utilización
resulta imprescindible para el estudio de cualquier proceso u operación unitaria y existen un
gran número de problemas que pueden resolverse mediante su conveniente aplicación.
Considérese, por ejemplo, un sistema como el expuesto en la figura BMM-1, en la que se
muestra el llenado de un tanque agitado que contiene, inicialmente, una disolución de
cualquier compuesto A. Dada esta situación, puede resultar de interés el tiempo necesario
para que la concentración de A alcance un determinado valor. Si al tanque entra A puro, la
concentración aumentará, mientras que entra solo disolvente, deberá disminuir. Una
complicación adicional se puede presentar si existe una corriente de salida del tanque, o
mas aún, si A experimenta una reacción química. Las respuestas a estas cuestiones y a
otras de esta índole, que relacionan corrientes de entrada y de salida en sistemas
cualesquiera, pasa por la formulación de adecuados balances macroscópicos de materia. El
objetivo final de la realización de un balance macroscópico de materia es llegar a conocer
los caudales y composiciones de las distintas corrientes de entrada y salida de un sistema y
ls cantidades totales y composiciones medias que están en el interior del mismo en un
momento dado.

Disol. A

Disol.A
Disol. A

Fig. BMM-1: Llenado de un tanque con cierta disolución de un compuesto A

Los balances macroscópicos de materia, son una generalización de la ley de conservación

de la masa a sistemas abiertos, esto es, sistemas con posibles entradas y/o posibles salidas
de materia al exterior. Su utilidad en el campo de la ciencia y tecnología es muy amplia, y en
4
 UNCP – FIQ Ing. Dr. O. Vilca

general, su complejidad matemática es escasa. Así, en el estado estacionario las

ecuaciones correspondientes son algebraicas, normalmente lineales, y si se trata de régimen
no estacionario, son diferenciales ordinarias de primer orden, siendo la magnitud tiempo la
variable independiente.

EJEMPLO 1:
Como etapa previa a la depuración de cierta corriente de gas residual con un 10% molar de
H2S, se debe rebajar su contenido de gas tóxico hasta el 0,1%. Para ello se ha decidido
emplear una columna de absorción, utilizándose agua como disolvente:
a) Determinar la cantidad de H2S absorbida por cada segundo por el agua, si el caudal
del gas es de 100 m3/h, medido en condiciones normales.
b) Si la columna de absorción trabaja a la presión atmosférica con agua impurificada con
un contenido en H2S de 0,001 g/L de agua. La solubilidad del agua en el gas es de
0,10 mol/L.atm; determinar el caudal mínimo de agua a emplear.
SOLUCIÓN:
El problema consiste, para la primera parte, en calcular las cantidades de H2S en la
corriente gaseosa a la entrada y a la salida del sistema (Fig. BMM-2). Para ello se
considerará el caudal del gas inerte (no se absorbe) que debe ser el mismo a la entrada y
a la salida del absorbedor:

H2O H2S +
Inerte

H2O +
H2S H2S +
Inerte
Fig. BMM-2: Columna de absorción de H2S

a) ENTRADA DE GAS INERTE = SALIDA DE GAS INERTE

5
FENÓMENOS DE TRANSPORTE – BALANCE MACROSCÓPICO
 6

m 3 tot. 90 m 3 inerte 1 kmol inerte 1h

BASE: 100 x 3
x 3
x = 1,1161x10 −3 kmol inerte / s
h 100 m tot. 22,4 m inerte 3600 s

kmol inerte 10 kmol H 2 S

ENTRADA: 1,1161x10 −3 x = 1,24 x10 − 4 kmol H 2 S / s
s 90 kmol inerte

kmol inerte 0,10 kmol H 2 S

SALIDA: 1,1161x10 −3 x = 1,11722 x10 −6 kmol H 2 S / s
s 99,9 kmol inerte

Por lo tanto, la cantidad de H2S absorbido por segundo será:

ENTRADA – SALIDA = 1,22883x10 −4 kmol H 2 S Absorbido / s

b) El caudal mínimo de agua a emplear, corresponderá a la saturación de H2S en la

corriente de salida. Si se realiza un balance de H2S en la fase líquida de la columna
de absorción (Fig. BMM-3), se tendrá:

w1

w3

w2

Fig. BMM-3: Balance de H2S en la fase líquida de la columna de absorción

w1, H 2 S + w3, H 2 S − w2, H 2 S = 0 (A)

Donde: wm , j : Es el caudal del componente j que entra o sale del sistema a través de

la corriente m-ésima.
El caudal de H2S, proveniente de la corriente gaseosa es de 1,22883x10 −4 kmol H 2 S / s
y es igual en este caso al caudal global w3.

6
 UNCP – FIQ Ing. Dr. O. Vilca

El caudal molar de H2S que entra con la corriente de agua, es conveniente expresarlo

en función del caudal másico de agua w1, H 2 S , puesto que los datos de
concentración a la entrada del absorbedor vienen dados en unidades másicas.
kmolH 2 S kg H 2 O 1 m 3 H 2 O 0,001kg H 2 S 1kmol H 2 S
( w1, H 2 S ) = w1, H 2O x x x
s s 1000 kg H 2 O 1 m3 H 2O 34 kg H 2 S

10 −6
= ( w1, H 2O )
34
kmol H 2 S kg H 2 O 1m 3 H 2 O 0,1kmol H 2 S
( w2, H 2 S ) = w2, H 2O x x x 0,1 atm
s s 1000 kg H 2 O atm. m 3 H 2 O

= (w2,H 2O ) 10 −5

Como: w1, H 2O = w2, H 2O

Reemplazamos estos valores en la ecuación (A):

10 −6
( w1,H 2O ) + 1,22883x10 −4 = ( w1,H 2O )10 −5
34
De donde se tiene: w1,H 2O = 12,3263 kg / s
Este valor corresponde al valor mínimo del agua, puesto que, de emplear menor cantidad,
antes de absorber la totalidad del H2S ( 1,22883x10 −4 kmol H 2 S / s ), el agua estaría ya
saturada.

EJEMPLO Nº 2

DESCARGA DE UN LÍQUIDO POR GRAVEDAD DE UN DEPÓSITO A TRAVÉS DE UN

TUBO

DETERMINAR:
• La altura de un líquido en cualquier instante.
• El tiempo requerido para vaciar completamente el depósito.
CASOS:
• El depósito puede ser cilíndrico, esférico, cono invertido, etc.
• El tubo puede estar en forma vertical, horizontal o inclinado.
• El fluido puede ser newtoniano o no newtoniano.
• El régimen de salida del fluido puede ser, laminar, transiente, o
turbulento.

Aplicar el caso de un DEPÓSITO CILÍNDRICO con tubo de salida vertical en régimen

laminar.

VARIABLES:
7
FENÓMENOS DE TRANSPORTE – BALANCE MACROSCÓPICO
 8

• Diámetro del recipiente [D]

• Diámetro del tubo de salida [d]
• Altura del líquido en el recipiente [h]
• Densidad del fluido []
• Tiempo [t]
Para t =0, la altura del líquido en el recipiente será [H].

D
D
Area = A
1 Volumen
Descendido

H H

Area = B

Volumen
Evacuado

Figura BMM-4: Depósito cilíndrico con tubo de salida vertical

FORMULACIÓN DEL MODELO

SUPOSICIONES (O HIPÓTESIS)
• Depósito en forma cilíndrica.
• Líquido se asume que es agua (fluido incompresible)
• El sistema es isotérmico
• Los diámetros del depósito y del tubo de salida son constantes.
• El régimen de salida es laminar.
• No hay generación o producción
• No hay flujo de entrada al sistema

Por la Ley General de CONSERVACIÓN:

ENTRADA + PRODUCCIÓN = SALIDA + CONSUMO + ACUMULACIÓN (I)

De acuerdo a las suposiciones realizadas, la ecuación (I), será:

(B Z )2 + ACUMULACION = 0 (1)

 z = z , m + z ,c
(Macroscópico)
Donde :  z , m = 0
Por tan to :  z =  z ,c (2)

8
 UNCP – FIQ Ing. Dr. O. Vilca

El término acumulación, desde el punto de vista integral es:


t v
Acumulación = dv

D 2
Pero : dv = dh = Adh (3)
4
h
d
dt H
Acumulación = Adh

Reemplazando en la ecuación (1):

h
d
Bv z +  Adh = 0 ; A,  : Constantes
dt H
dh
A = − B  v z  ; donde < v z > velocidad media (4)
dt

Realizando un Balance de Energía Mecánica (Ecuación de Bernoulli) entre el nivel superior

del depósito y la salida del tubo:

La ecuación de Bernoulli es:

2 P
^ ^
1 1
  vz2  + gh +  dP + W + Ev = 0 (5)
2 P1

Donde:
^
W : Velocidad a la que el sistema realiza trabajo mecánico.
^
Ev : Pérdida por fricción

v12 P v2 P ^ ^
+ gh1 + 1 = 2 + gh2 + 2 + W + Ev (6)
2  2 

Como:
P1 = P2 : No hay variación de la presión (presión atmosférica).
V1≈ 0 : Suponemos que el diámetro del depósito es grande.
^
W = 0 : No hay trabajo mecánico.
^
Ev : No hay pérdida por fricción

 vz2 
h= →  vt  = 2 gh (7)
2g

Donde:
<vt> = vt: Velocidad teórica, ya que se ha despreciado las pérdidas entre dos puntos.

9
FENÓMENOS DE TRANSPORTE – BALANCE MACROSCÓPICO
 10

COEFICIENTE DE VELOCIAD:

vr
cv = → vr = cv vt (8)
vt

Donde: cv: Coeficiente de velocidad

vr: Velocidad real

Reemplazando la ecuación (8) en la ecuación (7):

vr = cv 2 gh (9)

COEFICIENTE DE CONTRACCIÓN (cc)

B0
cc =
B
Donde: Bo: Área de contracción del orificio

DESCARGA REAL DEL ORIFICIO (Qr)

Qr = Bo v r
Qr = cc cv B 2 gh

Como: cd = cc cv : Coeficiente de descarga.

Bibliografía:
cd = 0.61: Entrada de borde afilado
cd = 0.80: Para tubo corto ligeramente redondeado.
cd = 0.98: Para orificio bien redondeado

Bv z = Qr = cd B 2 gh
(10)

Separando variables:

h t
dh
A = −cd B 2 g  dt
H
h 0

Integrando: 2 A[(h)1 / 2 ]hH = −cd B 2 g t

cd B 2 g t
H− h=
2A

El área para un cilindro:

10
 UNCP – FIQ Ing. Dr. O. Vilca

D 2
A=
4
 d2
B=
4

ALTURA DEL LÍQUIDO:

2
c 2gt  d 
h= H− d   (11)
2 D

TIEMPO DE VACIADO:

( ) ( )
2
2A 2 D
t= H− h =   H− h (12)
cd B 2 g cd 2 g  d 

Ejemplo Nº 3

Determinar la potencia necesaria para bombear agua en el sistema de la Figura Nº 2. El

agua cuya densidad es 1.0 g/cm3 y cuya viscosidad es de un centipoise, se eleva al depósito
superior a razón de 0.34 m3/min. Toda la conducción esta formada por una tubería de acero
circular lisa de tamaño nominal 4 pulgadas y código 40 ST.

24 m

12 m

2
36 m

90 m
mm
1,5 m m
1

Figura BMM-5: Flujo en una conducción con pérdidas por fricción provocadas por accesorios

De la ecuación de Bernoulli:

11
FENÓMENOS DE TRANSPORTE – BALANCE MACROSCÓPICO
 12

2 P
^ ^
1 1
  vz2  + gh +  dP + W + Ev = 0
2 P1

Haciendo un balance del punto “1” al punto “2”:

v22 P v2 P ^ ^
+ gh2 + 2 = 1 + gh1 + 1 + W + Ev (1)
2  2 

1. Cálculo de Ev:

^ 1 L  1 
Ev = i   vz2  f  + i   vz2  ev 
2 Rh i 2 i
1 L 
  2  v
i
2
z  f  : Pérdida de fricción por tuberías
Rh i
(2)

1 
  2  v
i
2
z  ev  : Pérdida de fricción por accesorios
i

1.1. Cálculo de Pérdida por fricción por las tuberías:

1 L  1 4L  2  v 2 
  2  v
i
2
z  f  = i  vz2 
2
f =
D i
f iLi =
 Rh i D
2(0.69) 2
= (0.0049)(1.5 + 90 + 36 + 24 + 12) = 7.46 m2 / s 2
0.10226

1.2. Cálculo de pérdida por fricción por accesorios:

1  v 
2

i 2 z v 
  v 2
 e  =
2
iev (3)
i

1.2.1. Calculo del Nº de Reynolds (Re)

 vd
Re =

De tablas: d = 4.026 pulg. = 0.10226 m

Dato: Q = 0.34 m3/min; como Q = <v>A,

m3  1min 
4(0.34 )  
4Q min  60 s  = 0.69m / s
 v = =
d 2  (0.10226)m 2

Por lo tanto: Re = 7.056x104 (Flujo Turbulento)

1.3. Cálculo del factor de fricción “f”

12
 UNCP – FIQ Ing. Dr. O. Vilca

i) Se puede determinar de la Figura 6.2-2 (R. Bird) el Factor de fricción vs. Número
de Reynolds
Donde f = 0.0049
ii) También se puede determinar de la siguiente relación:
0.0791
f = , para tuberías “hidráulicamente lisas”. Reemplazando estos valores se
Re1/ 4
tiene:
f= 0.004853

Como se puede observar los valores son muy cercanos.

1.4. Cálculo de ei (Pérdida de fricción en los accesorios):

Se determina a partir de los factores de pérdida por fricción utilizados en el sistema
para flujo turbulento. Se puede obtener estos valores a partir de la tabla 7.4-1 (R.
Bird).
i) Contracción brusca : 0.45
ii) 3 Codos de 90º : 1.50
iii) Ensanchamiento brusco : 1.00
Total : 2.95
1.5. Reemplazando en la ecuación (2):

Ev = 7.46 + 0.702 = 8.162 m2/s2

1.6. De la ecuación (1), se tiene:

<v1> = <v2> ; P1 = P2 ;
g (37.5 − 12) + W + 8.162 = 0
^

^
− W = 9.81(37.5 − 12) + 8.162 = 258.317m 2 / s 2
Por tanto la bomba comunica un trabajo de 258.32 m2/s2

1.7. Cálculo de la Potencia de la Bomba (W)

^
− W = W w ; Donde : w = flujo másico
m3 1 min Kg
w = Q = 0.34( )( )(1000 3 ) = 5.667 Kg / s
min 60 s m

Por lo tanto:

Kg f Kg Kg f .m
W = 26.33 (5.667) = 149.212  2 Hp
Kg.m s s

13
FENÓMENOS DE TRANSPORTE – BALANCE MACROSCÓPICO
 14

I.2. MÉTODOS DE CÁLCULO

• Integración de la ecuación de la continuidad
Planteamiento en un volumen de control macroscópico

Fig. BMM-6: Sistema macroscópico de flujo

Consideramos el sistema de la Fig. 1, en el que existe una sola entrada de fluido (de sección
S1, en el plano “1”) y una sola salida (de sección S2 en el plano “2”).
Estos planos de referencia se toman perpendiculares a las paredes del tubo.
Hacemos dos suposiciones que para la mayor parte de problemas no son muy restrictivas

SUPOSICIONES
a) En los planos “1” y “2” la velocidad de tiempo ajustado es paralela a las paredes de la
conducción.
b) b) La densidad y las demás propiedades físicas no varían a lo largo de la sección
transversal de los planos “1” y “2”
Aplicando al sistema de la expresión de la conservación de la materia se obtiene:

VELOCIDAD DE VELOCIDAD DE VELOCIDAD DE

ENTRADA DE SALIDA DE ACUMULACION (BMM-1)
MATERIA MATERIA DE MATERIA

14
 UNCP – FIQ Ing. Dr. O. Vilca

− − d
1  v1  S1 −  2  v2  S 2 = mtot (BMM-2)
dt
Donde: mtot : masa total del fluido.

−
Si: w = 1  v  S VC
 w = w2 − w1
El balance macroscópico de
materia en estado no estacionario <v1> <v2>

es:
d
mtot = −w (BMM-2.1)
dt 1 2

En estado estacionario, la masa total del fluido en el sistema no varía con el tiempo, y por lo
tanto el balance macroscópico de materia en estado estacionario es:

 w = w2 − w1 = 0 ; w1 = w2 (BMM-3)

II. BALANCE MACROSCÓPICO DE CANTIDAD DE MOVIMIENTO

Consideramos de nuevo el sistema de la Fig. BMM-6, haciendo las dos mismas

suposiciones anteriores. Aplicando la expresión de la conservación de la cantidad de

movimiento a este sistema, se obtiene la siguiente ecuación vectorial:

VELOCIDAD VELOCIDAD SUMATORIA FUERZAS VELOCIDAD DE

ENTRADA DE SALIDA DE QUE ACTÚAN SOBRE ACUMULACIÓN
MOMENTO MOMENTO EL SISTEMA DE MOMENTO
− −
1  v12  S1 −  2  v 22  S 2 + P1 S1 − P2 S 2  − F + mtot .g  =
d
Ptot (BMM- 4)
dt
− −
Donde: 1  v12  S1 −  2  v22  S 2 : Cantidad de movimiento global del fluido
P1S1 − P2 S2  : Fuerza de presión
F : Fuerza neta de las superficies sólidas que actúan sobre el sólido.
mtot .g: Fuerza de la gravedad que actúa sobre la masa total del fluido
15
FENÓMENOS DE TRANSPORTE – BALANCE MACROSCÓPICO
 16

Ptot : Cantidad de movimiento total del fluido

La ecuación anterior también se puede escribir de la siguiente forma:

 −2 
d v  
Ptot = − − w + PS  − F + mtot .g (BMM- 5)
dt  v 
 
El balance macroscópico de cantidad de movimiento en estado no estacionario es:
 −2 
v  
F = − − w + PS  + mtot .g (BMM- 6)
 v 
 
Este resultado es útil para calcular las fuerzas que actúan sobre piezas o aparatos, tales
como paletas de turbinas y curvas de tuberías, también es útil, juntamente con otros
balances macroscópicos, para describir el funcionamiento de eyectores de vapor de agua y
otros aparatos en las que hay mezcla de corrientes de fluidos.
− −
La relación  v 2  /  v  pueden evaluarse si se determinan experimentalmente los perfiles
de velocidad o es posible calcular por los métodos propuestos anteriormente. Sin embargo,
en la mayor parte de los problemas ingenieriles de flujo turbulento, donde el perfil de
− − −
velocidad es aproximadamente plana, se puede substituir  v 2  /  v  por  v  .

II.1. BALANCE MACROSCÓPICO DE ENERGÍA MECÁNICA (ECUACIÓN DE

BERNOULLI)
Para poder obtener el siguiente balance macroscópico de energía mecánica en estado no
estacionario para flujo isotérmico, es necesario integrar la ecuación de energía mecánica
(obtenida anteriormente):

 1 2  1 2 
 v  = − . v v  − (.Pv) − P(−.v) − (.[ .v]) − ( .v)  (v.g )
t  2   2 
sobre el volumen del sistema de la Figura BMM-4, la cual resulta:

 −
 
   3
  
(Ect + Ept + At ) = −  + Ep + G  w − W − Ev
d 1 v
dt  2  v   
  
En esta expresión, Ett, Ept y At, son respectivamente, la energía cinética, la energía
potencial y la energía libre (Helmholtz) totales del sistema de flujo.
W es la velocidad a la que el sistema realiza trabajo mecánico sobre los alrededores;
16
 UNCP – FIQ Ing. Dr. O. Vilca

Ev es la pérdida por fricción, es decir, la velocidad con la que la energía mecánica se

convierte irreversiblemente en energía calorífica.
  
La magnitud de G = H − T S es la energía libre (o energía libre de Gibbs) por unidad
de masa.
Para el flujo de un fluido de densidad constante, la ecuación anterior se transforma en:

 −
 
   3
 P 
d
(Ect + Ept ) = −  1 v
+ Ep +  w − W − Ev
dt 2 v  
  
En aquellos sistemas en los que el comportamiento de flujo es independiente del tiempo,
desaparece la derivada con respecto al tiempo. Puesto que w1=w2=w, se puede dividir toda
la ecuación por la velocidad de flujo de masa e introducir la notación:
 W  Ev
W= y Ev =
w w
P
 2
1
Además para sistemas isotérmicos se puede substituir: G =  )dP . De esta forma se
(
P1

obtiene el BALANCE MACROSCÓPICO DE ENERGÍA MECÁNICA en estado estacionario:
−
1  v3 
P
2
1  
 + Ep +  ( )dP + W − Ev = 0
2 v P1

La ecuación anterior puede evaluarse si se conoce la ecuación de estado (para sistemas
isotérmicos). Dos casos especiales son los siguientes:

a) GAS IDEAL ISOTÉRMICO

P2 P2
1 RT RT P
 ( )dP =  (
P1
 P1
M
)dP =
M
Ln( 2 )
P1

b) FLUIDO INCOMPRESIBLE
P2
1 1
 (  )dP =  ( P
P1
2 − P1 )

En la práctica ingenieril se admiten con frecuencia estos dos casos límite.

III. BALANCE MACROSCÓPICO DE ENERGÍA EN SISTEMAS NO ISOTÉRMICOS

III.1. BALANCE MACROSCÓPICO DE ENERGÍA

En el tema anterior se ha tratado de balances macroscópicos de materia, energía
mecánica y cantidad de movimiento, para el caso de sistemas que están a temperatura
17
FENÓMENOS DE TRANSPORTE – BALANCE MACROSCÓPICO
 18

constante.

En realidad, esta restricción es un poco artificial, puesto que en los sistemas reales de
flujo siempre se convierte algo de energía mecánica en energía calorífica.
Lo que realmente se supone en el tema de “isotérmicos” es que, o bien el calor que se
produce es demasiado pequeño para dar lugar a una variación apreciable de las
propiedades del fluido, o que se elimina inmediatamente a través de las paredes del sistema
que limitan el fluido.
Para el caso de un sistema no isotérmico, existen balances macroscópicos que describen
las relaciones entre las condiciones de la corriente a la entrada y a la salida.

Fig. BMM-7: Sistema macroscópico de energía

VELOCIDAD DE ACU- VELOCIDAD DE EN- VELOCIDAD DE SA- VELOCIDAD NETA DE

MULACIÓN DE ENER- TRADA DE ENERGÍA LIDA DE ENERGÍA ADICIÓN DE CALOR
GÍA INTERNA, CINÉTI- INTERNA, CINÉTICA INTERNA, CINÉTICA AL SISTEMA PROCE-
CA Y POTENCIAL Y POTENCIAL Y POTENCIAL DENTE DE LOS ALRED

VELOCIDAD NETA DE
TRABAJO PRODUCIDO (BMM-7)
POR EL SISTEMA SOBRE
LOS ALREDEDORES
− − − −
d
(U t + Ec + Ep) = 1  v1  U1 S1 −  2  v2  U 2 S 2 + 1  v1  Ec1 S1 −  2  v2  Ec2 S 2
dt
− − − −
+ 1  v1  Ep1 S1 −  2  v2  Ep2 S 2 + Q − W + P1  v1  S1 − P2  v2  S 2
(BMM-8)

18
 UNCP – FIQ Ing. Dr. O. Vilca

−
Si Et = U t + Ec + Ep y como w =   v  S

(Et ) = − U + Ec + Ep + PV w + Q − W

d −
(BMM-9)
dt   
Si el sistema de flujo se encuentra en estado estacionario:
d
(Et ) = 0 ;  w = w2 − w1 = 0 ; w1 = w2 = w
dt
Dividiendo entre w la ecuación (BM-10):

 −
  
U + Ec + Ep + P V  = Q −W (BMM-10)
 
Donde:

Q : Calor suministrado por unidad de masa

W : Trabajo producido por unidad de masa
La magnitud: H = U + PV
H = U + PV (BMM-11)
 
 H + Ec + Ep = Q+ W (BMM-12)
En los cálculos ingenieriles resulta útil el uso de las dos formas de H :
a) GASES IDEALES
PV = nRT
T2
R
T2

c p − cv = R H =  c p dT =   − 1 dT (BMM-13)
T1
M T1

cp
=
cv

b) LíQUIDOS INCOMPRESIBLES
 : Constante
T2
H =  c p dT +
1
(P2 − P1 ) (BMM-14)
T1

c p = cv

Para muchos sistemas de ingeniería química, los efectos de la energía cinética,

energía potencial y del trabajo son pequeños, de forma que:

H = Q

19
FENÓMENOS DE TRANSPORTE – BALANCE MACROSCÓPICO
 20

EJEMPLO DE APLICACIÓN

Por el tubo interior de un cambiador de calor se introduce 100 kg/h de aire seco a 150
°C y 2 atmósferas absolutas de presión y una velocidad de 30 m/s. El aire sale del
intercambiador a -18 °C y 1 atmósfera; tres metros encima del punto de entrada del
cambiador. Calcular la velocidad de de eliminación de energía a través de la pared del
tubo. Supóngase que el flujo es turbulento, que el gas se comporta como ideal y
utilizar para el calor específico del aire:

cP = 6,39 + 1,76x10−3T − 2,65x10−7 T 2 ;

c P = T = K
kcal
kmol K

SOLUCION

De acuerdo al problema se tiene el siguiente gráfico:

P2=1 atm.
T2=-18°C
Salida de
aire
2
Entrada del
liquido frio

3m

Salida del líquido

frio

1
Entrada de
aire seco
m=100 kg/h
T1=150°C
P1=2 atm.
<v>=30 m/s

De la ecuación BMM-12:
 
H + Ec + Ep = Q+ W

20
 UNCP – FIQ Ing. Dr. O. Vilca


Donde el término W no existe debido a que no hay trabajo del sistema, por tanto la ecuación
queda de la siguiente forma:

H + Ec + Ep = Q

i) CALCULO DE H :
Integrando desde una T1=150 °C (423,15 K) hasta T2=-18 °C (255,15 K), se
obtiene el siguiente valor para la variación entálpica:
H = −40,2957 kcal / kg
ii) CALCULO DE Ec :
1
Ec = (v 22  − v 12  )
2
Del balance macroscópico de materia en estado estacionario:
w2 − w1 = 0 → w2 = w1 → 2 v2  S2 = 1 v1 S1
Como S2 = S1 , entonces despejando v2  se tiene
1 v1 
v2  =
2
Considerando un gas ideal como se plantea en el problema, se tendrá:
−
n P PM
PV = nRT → = → =
V RT RT
De donde reemplazando en la ecuación se tiene:
P1T2
 v2  = v1 
P2T1
Por tanto, reemplazando en la ecuación de la energía cinética, se tiene:

1  P1T2  2   v 2   PT  2 
2

Ec =    v 1  −  v 1   = 1  1 2  − 1

2

2  P2T1   2  P2T1  
 
Reemplazando datos, en ésta última ecuación, tendremos el valor de la energía
cinética, que será necesario llevar a las mismas unidades de energía, que la de la entalpía,
para así de esta forma comparar los valores de cada una de ellas.
Ec = 0,04882 kcal / kg
iii) CÁLCULO DE Ep :
g (h2 − h1 )
Ep = = 0,00703 kcal / kg
gc

21
FENÓMENOS DE TRANSPORTE – BALANCE MACROSCÓPICO
 22


iv) CÁLCULO DE Q

− Q = H − Ec − Ep

kcal  kg 
− Q = (40,2957 − 0,04882 − 0,00703) = 40,23985 1000 
kg  h 

− Q = 4023,985 kcal / h

IV. BALANCE MACROSCÓPICO EN SISTEMAS DE VARIOS COMPONENTES

IV.1. BALANCE MACROSCÓPICOS DE MATERIA
La ley de conservación de la materia de la especie química i en un sistema
macroscópico de flujo, establece:
d
mi ,tot = −wi + wim + ri ,tot i = 1,2,3,..., n (BMV-1)
dt
Donde:

mi ,tot : La masa total de i existente en el sistema en un instante.

− wi = wi 2 − wi1 = i 2 v2  S2 − i1 v1 S1 : Es la diferencia entre la velocidad de flujo
de materia de la especie i en los planos <1> y <2>.

wim : es la velocidad de flujo de materia de la especie i a través de las superficies

limitantes del sistema debido a la transferencia de materia.

wim es positivo cuando se adiciona materia al sistema.

IV.2. BALANCE MACROSCÓPICO DE CANTIDAD DE MOVIMIENTO

El planteamiento de la ley de conservación de cantidad de movimiento para una mezcla de
fluidos, cuya única fuerza externa que actúa sobre las especies es debido a la gravedad, es:

d   v2  
P = − w + PS  + F ( m ) − F + mtot g (BMV-2)
dt  v 

IV.3. BALANCE MACROSCÓPICO DE ENERGIA

La ley de conservación de la energía para una mezcla de fluidos se expresa de la
siguiente forma:

22
 UNCP – FIQ Ing. Dr. O. Vilca

d  
1   v3   
Etot = − U + P V +  +   w + Q ( m ) + Q − W (BMV-3)
dt  2 v   
Cuando se producen reacciones químicas se pueden absorber o desprender una
considerable cantidad de calor. Este calor de reacción se tiene automáticamente en cuenta
al calcular las entalpías de las corrientes d entrada y salida.

IV.4. BALANCE MACROSCÓPICO DE ENERGÍA MECÁNICA

d  1   v3  P 
( Ktot + tot ) = −   +  +  w + B ( m ) − W − Ev (BMV-4)
dt 2   v   

23
FENÓMENOS DE TRANSPORTE – BALANCE MACROSCÓPICO
 24

BIBLIOGRAFÍA:

1. R. B. BIRD Y OTROS, “Transport Phenomena”, Editorial Reverté, S.A, España, 2002

2. J. COSTA LOPEZ Y OTROS, “Curso de Ingeniería Química”, Editorial Reverté, S.A.
España, 1994
3. W. J. BEEK, K. M. K. MUTTZALL. “Transport phenomena”, Editorial Jhon Wiley & Sons Ltd.
1975.
4. MIRANDA OSCAR, CAMPOS DANTE, “Problemas de mecánica de fluidos e hidráulica”.
Impreso en el Perú. 1991.
5. INCROPERA FRANK P. "Fundamentos de Transferencia de Calor". Editorial Prentice Hall
1999.
6. TREYBAL ROBERT. 1980. "OPERACIONES DE TRANSFERENCIA DE MASA". Segunda
Edición. McGraw Hill, México.
7. J. R. WELTY, C.E. WICKS Y R.E WILSON, “Fundamentos de la Transferencia de
Momento, Calor y Masa”, 1ra Reimpresión, Editorial LIMUSA, México, 1982.
8. A. FOUST “Principio de Operaciones Unitarias”, Editorial CECCSA, México, 1986.
9. CHRISTIE J. GEANKOPLIS, “Procesos de Transporte y Operaciones Unitarias”, 1ra
Edición, CECCSA, México, 1986.

24