Balances de energías en sistemas reactivos

Andrés F. Tovar Parada. 1​​ (​[email protected]​), Camilo A. Lizarazo C.​1

(​[email protected]​), David Camilo Guerrero G. 1​​ (​[email protected]​),
1. Programa curricular: Ingeniería Química, Departamento de Ingeniería Química y Ambiental, Facultad de
Ingeniería, Universidad Nacional de Colombia-Sede Bogotá.

Balances de sistemas reactivos

Problema 9.5-1
Cien moles de N H 3 por segundo y 200 moles de O2 por segundo son alimentadas a un reactor
donde el amoniaco se consume por completo. El gas producido emerge a 300℃. Calcular la
velocidad del flujo de calor de la operación, suponiendo una presión aproximada de 1 atm.

Inicialmente se plantea el diagrama del problema ordenando la información dada en el

enunciado:

Posteriormente se realiza el balance de energía para el sistema del reactor por medio de la
primera ley:

Q = H3 − H1 − H2

Realizando el despeje trivial se obtiene la entalpía de la mezcla del flujo (3) en términos de las
entalpías residuales y de gas ideal:

H 3 = H 3R + Σni,3 · h*i,3
De esta manera se realiza el balance de energía completo, realizando el despeje para cada
corriente:

Q = H 3R − H 2R − H 1R + n · h*H + nN O ;3 · h*N O ;3 + nO ;3 · h*O ;3 − n · h*N H −n · h*O ;2

H 2 O;3 2 O;3 2 2 N H 3 ;1 3 ;1 O2 ;2 2
Posteriormente se realiza el balance molar utilizando la reacción balanceada del reactor:

Sustancia Moles de entrada Moles de salida

N H 3 (gas) A A − 4ξ = 0

O2 (gas) B B − 5ξ

N O (gas) 0 4ξ

H 2 O (vap) 0 6ξ

Puesto que en el enunciado se menciona que el amoniaco se consume totalmente resulta sencillo
determinar el avance de la reacción como:

4ξ = A

Sustituyendo las moles de entrada y salida en el balance de energía se tiene:

Q = H 3R − H 2R − H 1R + 6ξ · h*H + 4ξ · h*N O ;3 + (B − 5ξ) · h*O ;3 − 4ξ · h*N H − B · h*O ;2

2 O;3 2 3 ;1 2

Agrupando los avances se tiene:

Q = H 3R − H 2R − H 1R + (6 · h*H + 4 · h*N O ;3 + 5 · h*O ;3 − 4 · h*N H ;1 ) · ξ + B · (h*O ;3 − h*O ;2 )

2 O;3 2 3 2 2

En la literatura es posible encontrar la entalpía de reacción para la oxidación de amoniaco

(H oR = ) [1].FELDER, sin embargo dicha información sólo se suministra para una temperatura
de 298.15K, por lo que es necesario dejar expresadas cada una de las entalpías que describen la
reacción en términos de dicha temperatura:

573,15
h*H = h*H + ∫ cp,H O dT
2 O;3 2 O;o
298,15 2

De esta manera se tiene para las otras sustancias:

573,15
h*N O ;3 = h*N O ;o + ∫ cp,N O dT
298,15
 573,15
h*
O2 ;3
= h*
O2 ;o
+ ∫ cp,O dT
298,15 2

Puesto que los flujos (2) y (1) ya se encuentran a temperatura de 298.15 K no es necesario
corregir sus entalpías. Nuevamente despejando el balance de energía se tiene:

Q = H 3R − H 2R − H 1R + (6 · h* + 4 · h* + 5 · h* − 4 · h* ) · ξ + B · (h* − h* )+
H 2 O;o N O ;o O2 ;o N H 3 ;1 O2 ;3 O2 ;2
573,15 573,15
(6 · ∫ cp,H O dT + 4 · ∫ cp,N O dT ) · ξ
298,15 2 298,15

De esta manera se tiene:

H oR = (6 · h* + 4 · h* + 5 · h* − 4 · h* KJ
) =− 904, 74 mol
H 2 O;o N O ;o O2 ;o N H 3 ;1

Las integrales se calcularon por medio de Excel (hoja 1) utilizando las siguientes expresiones de
los Cp:

c = 33, 46 × 10−3 + 0, 688 × 10−5 T + 0, 7604 × 10−8 T 2 − 3, 593 × 10−12 T 3

p,H 2 O

c = 29, 5 × 10−3 + 0, 8188 × 10−5 T − 0, 2925 × 10−8 T 2 + 0, 3652 × 10−12 T 3

p,N O

Adicionalmente se tiene el cambio de entalpía del oxígeno que pasa de la corriente (2) a la
corriente (3):

573,15
h*O ;3 − h*O ;2 = ∫ cp,O dT = 8, 474 mol
KJ
2 2 298,15 2

Finalmente se tiene el cálculo de las entalpías residuales para cada corriente. Tanto la corriente
(1) como la (2) están conformadas por sustancias puras por lo que fácilmente se pueden calcular
por medio de tablas [2]:

Iniciando con la temperatura reducida para cada corriente:

T
TR = Tc

T R1 = 0, 73 T R2 = 1, 93
Posteriormente se calcula la presión reducida:

P
PR = PR

P R1 = 0, 0089 P R2 = 0, 002

De esta manera utilizando las tablas de Byung y Kesler [5], para las entalpías residuales por
medio de las correlaciones de pitzer y colaboradores, se determinó la entalpía residual:

HR R
H 0 R
H 1
RT c
= [ RT c
] + ω[ RT c
]

H 1R =− 4, 12 Kg
KJ
H 2R =− 0, 2219 Kg
KJ

Puesto que el flujo (3) corresponde a una mezcla de gases, la tercera entalpía residual fue
calculada por medio de un programa en Excel (hoja 2). Dicho programa utiliza una ecuación
virial truncada en el segundo coeficiente que se muestra a continuación:

Donde:

De esta manera se tiene la entalpía residual del flujo (3):

H 3R =− 0, 0002961 Kg
KJ

Finalmente se tiene el balance de energía como:

 KJ KJ KJ KJ KJ
Q =− 0, 0002961 Kg · F 3 + 0, 2219 Kg · F 2 + 4, 12 Kg · F 1 + ξ · (− 904, 74 mol ) + ξ · (91, 23 mol )
mol KJ
+ 200 s
· (8, 474 mol )
Con el balance molar:

A = 100mol/s ; B = 200mol/s y ξ = 25 mol/s

Luego los flujos en kilogramos:

Kg Kg Kg
F 1 = 1, 7 S
F 2 = 6, 4 S
F 3 = 8, 1 S

Donde el calor que debe fluir es:

Q =− 18635, 68 KW

Al comparar el resultado obtenido por el de la literatura [1], se tiene:

19700−18635,68
%Error = 19700
= 0.054

Con un error porcentual del 54%, es posible establecer que aunque el modelo ideal propuesto por
Felder [1], posea errores considerables, puede ser de gran ayuda para realizar cálculos a
condiciones cercanas a la idealidad.

problema 9.5-2

se oxida metano con aire para producir formaldehído en un reactor continuo. La combustión del
metano para formar CO2 es una reacción que compite con la anterior.
1. CH 4 (g) + O 2 → H CHO (g) + H 2 O (V )
2. C H 4 (g) + 2 O 2 → CO 2 + 2H 2 O (V )

El siguiente es el diagrama de flujo del proceso para la base supuesta de 100 mol de alimentación
de metano al reactor.

Solución:
Siendo un proceso isobárico sin realización de trabajo, la expresión de la primera ley de la
termodinámica se reduce a :

Q = ΔH
descomponiendo ΔH en corrientes:

Q = Hd − Ha − Hb

donde la corriente ​a ​es la alimentación de metano, la corriente ​b e​ s la alimentación de oxígeno, la

corriente ​b’ e​ s la alimentación de nitrógeno molecular y la corriente ​d e​ s la corriente de salida.
siendo el nitrógeno molecular inerte en el proceso, la corriente de ingreso de aire no será
considerada como una solución de gases sino dos corrientes individuales. por el contrario, la
corriente de salida si será considerada una solución de vapores y gases.
Con el fin de modelar el comportamiento real del proceso, se expresa el balance de energía en
función de entalpías residuales:

Q = H dR − H aR − H R
b −H R
b′ − a · h *CH 4 ;a − b · h *O 2 ;b − b′ · h *N 2 ;b′ + H *
d

Para escribir el balance de la corriente d es necesario primero escribir el balance de masa en

función de la coordenada de reacción:

Sustancia Moles de entrada Moles de salida

 C H 4 (g) A A−ξ 1−ξ 2

O2 (gas) B B − ξ 1 − 2ξ 2

N 2 (gas) B’ B’

HCHO(gas) - ξ1

CO2 (gas) - ξ2

H 2 O (vap) - 2ξ 2

de esta forma, expresamos H * como:

d

H * = (a − ξ 1 − ξ 2 ) · h *CH + (b − ξ 1 − 2ξ 2 ) · h *O + b′ · h *N + ζ 1 h *HCHO;d + ζ 2 h *CO + (ξ 1 − 2ξ 2 ) · h *H

d 4 ;d 2 ;d 2 ;d 2 ;d 2 O;d

Reorganizando y factorizan las coordenadas de reacción:

Q = H dR − H aR − H R
b −H R
b′ − a · (h *CH − h *CH ) − b · (h *O − h *O ;d ) − b′ · (h *N − h *N ;d )
4 ;a 4 ;d 2 ;b 2 2 ;b′ 2

+ ξ 1 (− h *CH ;d − h *O ;d + h *HCHO;d + h *H O;d )

4 2 2
* * * *
+ ξ 2 (− h CH ;d − 2h O ;d + h CO ;d + 2h H O;d )
4 2 2 2

Sin embargo, estos términos de entalpía se encuentran en las condiciones de la corriente de

salida, por tanto deberán ser expresados en su forma integral más el calor de formación en el
estado de referencia, al realizar esta operación y factorizar los calores de formación y las
integrales obtenemos:
d d d d
+ ξ 1 (− h *CH 4 ;25 − h *O2 ;25 + h *HCHO;25 + h *H 2 O;25 ) + ξ 1 ( - ∫ cp,CH 4 dT − ∫ cp,O dT + ∫ cp,HCHO dT + ∫ cp,H O dT )
2
ref ref 2 ref ref
d d d d
+ ξ 2 (− h *
CH 4 ;25 − 2h *
O2 ;25 +h *
CO 2 ;25 + 2h *
H 2 O;25 ) + ξ 2( − ∫ cp,CH dT − 2 ∫ cp,O dT + ∫ cp,CO dT + 2 ∫ cp,H O dT )
4 2 2
ref ref 2 ref ref
siendo los factores de ξ 1 y ξ 2 el calor de reacción estándar de cada reacción. En resumen :

Q = H dR − H aR − H R
b −H R
b′ − a · (h *CH 4 ;a − h *CH 4 ;d ) − b · (h *O2 ;b − h *O2 ;d ) − b′ · (h *N 2 ;b′ − h *N 2 ;d )
d d d d
o
+ ξ 1 · (ΔH r1 ) + ξ 1 · (− ∫ cp,CH dT − ∫ cp,O dT + ∫ cp,HCHO dT + ∫ cp,H O dT )
4 2
ref ref 2 ref ref
d d d d
o
+ ξ 2 · (ΔH r2 ) + ξ 2 · (− ∫ cp,CH 4 dT − 2 ∫ cp,O dT + ∫ cp,HCHO dT + 2 ∫ cp,H O dT ) 2
ref ref 2 ref ref

Para el cálculo de las respectivas integrales se reemplaza el c p por su respectivo polinomio en

función de la temperatura:
 c = 1, 2 + 3, 25 × 10−3 T + 0, 75 × 10−6 T 2 − 0, 71 × 10−9 T 3
p,CH 4

cp,O = 0, 88 − 0, 0001 × 10−3 T − 0, 54 × 10−6 T 2 − 0, 33 × 10−9 T 3

2

c = 1, 79 + 0, 107 × 10−3 T + 0, 586 × 10−6 T 2 − 0, 20 × 10−9 T 3

p,H 2 O

cp,HCHO = 33.26 J/mol · K (Varía menos de 0,02 en el intervalo de temperatura).

Al realizar los cálculos de todas las entalpías, la expresión resultante:

Q = H dR − H aR − H bR − H bR′ − 18065,6585 KJ
Utilizando el esquema de cálculo de entalpía residual del problema 9,5-1 (cálculos en el archivo
excel), se obtiene:

H aR =− 0, 982983514
H bR =− 0, 222904388
H bR′ =− 0, 01014888
H dR =− 0, 000119284

la computar los datos obtenemos el calor total:

Q =− 18064, 44258 KJ

problema 9,5-3

la deshidrogenación de etanol para formar acetaldehído:

1. C 2 H 5 OH (v) → C H 3 CHO (v) + H 2(v)

se lleva a cabo en un reactor adiabático continuo. Se alimenta vapor de etanol al reactor 400ºC y
se logra una conversión de 30%. Calcule la temperatura del producto.
los balance de materia conducen a la información que se incluye en el siguiente diagrama de
flujo:

Solución:
siendo un proceso adiabático e isobárico, el balance de energía se reduce a el balance de
entalpías:
ΔH = H b − H a
Donde ​a ​ representa la corriente de entrada de etanol, y ​b​ representará la corriente de salida.
Realizando el balance de masa:

sustancia Moles de​ e​ ntrada Moles de​ s​ alida

C 2 H 5 OH A A- ζ

C H 3 CHO - ζ

H 2 - ζ

Expresando el balance en términos de entalpías residuales:

ΔH = H bR − H aR − a · h *C H OH;a + H *
b
2 5

donde H * en términos de la coordenada de reacción:

b

H *b : (a − ζ ) · h *C + ζ · h *CH + ζ · h *H
2 H 5 OH;b 3 CHO;b 2 ;b

por tanto, la expresión resultante:

ΔH = H bR − H aR − a · h *C H + (a − ζ ) · h *C + ζ · h *CH + ζ · h *H
2 5 OH;a 2 H 5 OH;b 3 CHO;b 2 ;b

debido a que la corriente de entrada no se encuentra a la temperatura de referencia, es necesario

reescribirla en términos de la entalpia de formacion de referencia. este proceso se realiza
también para la corriente de salida, con el fin de hacer explícito el término del calor estándar de
la reacción.

a T
ΔH = H bR − H R
a − a ( h *C 2 H 5 OH;25 + ∫ cp,C 2 H 5 OH
dT ) + a(h *C 2 H 5 OH;25 + ∫ cp,C 2 H 5 OH
dT )
ref ref
T T T
− ζ (h *
C 2 H 5 OH;25 + ∫ cp,C 2 H 5 OH
dT ) + ζ (h *
CH 3 CHO;25 + ∫ cp,CH CHO dT ) + ζ (h
3
*
H 2 ;25 + ∫ cp,H dT ) 2
ref ref ref
Reorganizando:
a T
0 = ΔH = H bR − H R
a -a · ∫ cp,C 2 H 5 OH
dT + a · ∫ cp,C 2 H 5 OH
+ ζ ( h *CH 3 CHO;25 + h *H 2 ;25 − h *C 2 H 5 OH;25 ) +
ref ref
T T T
ζ( ∫ cp,CH CHO dT + ∫ cp,H dT - ∫ cp,C
3 2 2 H 5 OH
dT )
ref ref ref
lo cual es equivalente a :
 a T
0 = ΔH = H bR −H R
a -a · ∫ cp,C 2 H 5 OH
dT + a · ∫ cp,C 2 H 5 OH
+ ζ ( h *CH + h *H − h *C )+
3 CHO;25 2 ;25 2 H 5 OH;25
ref ref
T T T
ζ( ∫ cp,CH CHO dT + ∫ cp,H dT - ∫ cp,C
3 2 2 H 5 OH
dT )
ref ref ref

Problema 9.5-4

La reacción de deshidrogenación del ejemplo 9.5-3 se realiza introduciendo la alimentación a

300º C. La alimentación contiene 90.0mol% y el balance del acetaldehído entra al reactor a
razón de 150 mol/s.. Para evitar que la temperatura descienda demasiado y entonces la velocidad
de reacción disminuya a un nivel bajo inaceptable, se transfiere calor al reactor. Cuando la
velocidad de adición de calor es de 2440 kW, la temperatura de salida es 253ºC. Calcule la
fracción de conversión del etanol que se logra en el reactor.

Primero realizamos un diagrama del reactor y las corrientes de entrada y salida de este:

Después procedemos a plantear el balance de materia:

-Balance de C:
(0, 900) (150) (2) + (0, 100) (150) (2) = 2n1 + 2n2

Dividimos todos los términos de la ecuación por 2

150 (0, 900 + 0, 100) = n1 + n2

150 = n1 + n2

-Balance de H:
(0.900) (150) (6) + (0, 100) (150) (4) = 6n1 + 4n2 + 2n3

Dividimos todos los términos de la ecuación por 2

150 [(0, 900) (3) + (0, 100)(2)] = 3n1 + 2n2 + 1n3

150 (2, 9) = 3n1 + 2n2 + 1n3

 435 = 3n1 + 2n2 + 1n3

-Balance general:

Sustancia Entrada (mol/s) Salida (mol/s)

C2H3OH A A-ε
CH3CHO B B+ ε
H2 - ε

Del balance general se deduce:

B + ε = n2

n2 = B + n3

ε = n3

reemplazando esta última ecuación en las 2 anteriores tenemos

n1 = A − ε

n1 = A − n3

y también calculamos los valores de A y B

A = (0, 900) (150) = 135

B = 150 − 135 = 15

Ahora procedemos con el balance de energía del reactor

Q = H S − H E = H SR + H S* − H RE − H E*

donde:

H S* = ∑ ni h*i = (A − n3 ) hC2H5OH
* * *
253ºC + (B + n3 ) hCH3CHO 253ºC + n3 hH2 253ºC

H E* = AhC2H5OH
* *
300ºC + B hCH3CHO 300ºC

Reemplazando estos términos tenemos:

Q = H SR − H RE + (A − n3 ) hC2H5OH
* * * * *
253ºC + (B + n3 ) hCH3CHO 253ºC + n3 hH2 253ºC − AhC2H5OH 300ºC − B hCH3CHO 300ºC
R R * * * * * * *
Q = H S − H E + (A) (hC2H5OH 253ºC − hC2H5OH 300ºC ) + (B) (hCH3CHO 253ºC − hCH3CHO 300ºC ) + n3 (hH2 253ºC + hC2H5OH 253ºC − hCH3CHO 253ºC )

( )
300ºC 300ºC
Q = H SR − H RE − (A) ∫ CpC2H5OH dT − (B) ( ∫ *
CpCH3CHO dT ) + n3 (hH2 * *
25ºC + hC2H5OH 25ºC − hCH3CHO 25ºC
253ºC 253ºC

253ºC
+ ∫ (CpCH3CHO + CpH2 − CpC2H5OH )dT )
25ºC

Los Cp se pueden calcular usando los siguientes polinomios:

 CpC2H5OH = (61, 34x103 ) + (15, 72x105 )T − (8, 749x108 )T 2 + (19, 83x1012 )T 3

CpH2 = (28, 84x103 ) + (0, 00765 x105 )T − (0, 3288 x108 )T 2 − (0, 898 x1012 )T 3

Mientras que el Cp del acetaldehído se puede tomar como una constante para este intervalo de
temperatura, y su valor sería de 57,3 J/mol k

lo siguiente es calcular los cambios de entalpías molares y por último despejar el valor de n3 de la
ecuación; se debe tener en cuenta que el cálculo de la entalpía residual de la corriente de salida
depende de este mismo valor, por lo tanto se debe proceder con un método iterativo con el fin de hallar
n3

( )
300ºC 300ºC
n3 = [Q − H RS + H RE + (A) ∫ CpC2H5OH dT + (B) ( ∫ CpCH3CHO dT )}/(h*H2 25ºC + h*C2H5OH 25ºC − h*CH3CHO 25ºC
253ºC 253ºC

253ºC
+ ∫ (CpCH3CHO + CpH2 − CpC2H5OH )dT )
25ºC

al reemplazar estos valores tenemos que el valor de n3 es

n3 = 43 mol
s =ε

Problema 9,5-6

Una solución acuosa al 10,0% por peso de H2SO4 a 40ºC se neutralizara con una solución acuosa de
NaOH al 20,0% por peso a 25ºC en un reactor continuo. ¿A qué velocidad en kJ/kg de solución de
H2SO4 se debe retirar calor si la solución que se produjo emerge a 35ºC?

Para una base de cálculo de 1 kg de solución de H2SO4 el balance de materia queda de la siguiente
forma

Sustancia Entrada Salida

H2SO4 1,019 mol 0

NaOH 2,038 mol 0

H20 68,115 mol 2ε

Na2SO4 ε

El balance de energía para este reactor quedaría de la siguiente manera

Q = △H E + nC hCid − nF hFid + nG hGid

Q = △H E + εhNida2SO4 25ºC + (2ε + 68, 115)hH2O
id id id id id
35ºC − 50hH2O 40ºC − εhH2SO4 40ºC − 18, 115hH2O 40ºC − 2εhN aOH 25ºC

donde:
40ºc
id id
hH2SO4 40ºC = hH2SO4 25ºC + ∫ C pH2SO4 dT
25ºC

35ºc
hNida2SO4 35ºC = hNida2SO4 25ºC + ∫ C pN a2SO4 dT
25ºC

35ºc
id id
hH2O 35ºC = hH2O 25ºC + ∫ C pH2O dT
25ºC

hNidaOH 25ºC = hNidaOH 25ºC

35ºc 40ºc
Q = △H E + εhNida2SO4 25ºC + (2ε + 68, 115)(hH2O
id
25ºC +
id
∫ C pH2O dT ) − 50(hH2O 25ºC + ∫ C pH2O dT )
25ºC 25ºC
40ºc
id id id
− ε(hH2SO4 25ºC + ∫ C pH2SO4 dT ) − 18, 115hH2O 25ºC − 2εhN aOH 25ºC
25ºC

Problema 9.6-2

En un horno se quema metano líquido con un 100% de aire en exceso. Cual es la temperatura
más alta que tendría que soportar las paredes del horno?. El metanol entra a 25℃ y el aire a 100
℃.

Inicialmente se plantea el diagrama del problema ordenando la información dada en el

enunciado:

De esta manera se tiene el balance de energía para el sistema de temperatura de flama adiabática:

0 = H3 − H1 − H2
Realizando el despeje trivial se descompone la entalpía real en términos de las entalpías
residuales y de gas ideal:

H = H R + Σni · h*i

De esta manera se realiza el balance de energía completo teniendo en cuenta que el flujo (1) solo
posee metanol puro líquido:

0 = H 3R − H 2R + n · h* + n · h* + n · h* + n · h* −n · ho
H 2 O;3 H 2 O;3 CO2 ;3 CO2 ;3 O2 ;3 O2 ;3 N 2 ;3 N 2 ;3 CH 3 OH;1 CH 3 OH;1

−n · h* − n · h*
O2 ;2 O2 ;2 N 2 ;2 N 2 ;2

Posteriormente se realizan los balances de materia utilizando la reacción balanceada del reactor y
los datos del enunciado:

Sustancia Moles de entrada Moles de salida

C H 3 OH (L) A A−ξ =0

O2 (gas) B B − 32 ξ

N 2 (gas) 3,76B 3,76B

CO2 (gas) 0 ξ

H 2 O (vap) 0 2ξ

Puesto que del enunciado se puede inferir que el metanol se consume totalmente, resulta sencillo
determinar el avance de la reacción como:

ξ=A

Sustituyendo las moles de entrada y salida en el balance de energía se tiene:

0 = H 3R − H 2R + 2ξ · h*H + ξ · h*CO ;3 + (B − 32 ξ) · h*O ;3 + 3, 76B · h*N − ξ · hoCH

2 O;3 2 2 2 ;3 3 OH;1

−B· h* − 3, 76B · h*
O2 ;2 N 2 ;2
Agrupando los avances se tiene:

0 = H 3R − H 2R + (2 · h* + h* − 3
2
· h* − ho )·ξ
H 2 O;3 CO2 ;3 O2 ;3 CH 3 OH;1
+ B · (h* − h* + 3, 76B · h* − 3, 76B · h* )
O2 ;3 O2 ;2 N 2 ;3 N 2 ;2

En la literatura es posible encontrar la entalpía de combustión para el metanol (H oc = )

[1].FELDER, sin embargo dicha información sólo se suministra para una temperatura de
298.15K, por lo que es necesario dejar expresadas cada una de las entalpías que describen la
reacción en términos de dicha temperatura. Por otro lado, como se menciona en la descripción
del problema, la temperatura adiabática de salida del flujo (3) debe dejarse expresada puesto que
es la incógnita que se desea determinar:

T ad
h*
H 2 O;3
= h*
H 2 O;o
+ ∫ cp,H O dT
298,15 2

T ad
h*
CO2 ;3
= h*
CO2 ;o
+ ∫ cp,CO dT
298,15 2

T ad
h*
O2 ;3
= h*
O2 ;o
+ ∫ cp,O dT
298,15 2

Y para los flujos de oxígeno y nitrógeno que no reaccionan:

T ad
h*O ;3 − h*O ;2 = ∫ cp,O dT
2 2 100 2

T ad
h*N
;3
− h*N
;2
= ∫ cp,N dT
2 2 100 2

Continuando con el balance de materia:

T ad
0 = H 3R − H 2R + (2 · h*H + h*CO ;3 − 3
2
· h*O ;3 − hoCH )·ξ+B·( ∫ cp,O dT )
O;3 OH;1
2 2 2 3 100 2

T ad T ad T ad T ad
+ (2 · ∫ cp,H O dT + ∫ cp,CO dT − 32 · ∫ cp,O dT ) · ξ + 3, 76B · ( ∫ cp,N dT )
298,15 2 298,15 2 298,15 2 100 2

De esta manera se tiene:

 H oc = (2 · h* + h* − 3
2
· h* − ho KJ
) = − 726, 6 mol
H 2 O;3 CO2 ;3 O2 ;3 CH 3 OH;1

Finalmente los balances másicos se muestran a continuación:

Con una base de cálculo de 1 mol de metanol líquido en la corriente (1):

A = 1mol

2· 32 mol O2
B = 1mol CH 3 OH · = 3 mol O2
1mol CH 3 OH
Luego las moles de nitrógeno son:

N 2 (gas) = 3.76 · B = 11.28 mol

Y las moles de salida:

CO2 (gas) = ξ = 1 mol

H 2 O (vap) = 2ξ = 2 mol

Puesto que los flujos (2) y (3) corresponden a una mezcla de gases, sus entalpías fueron
calculadas por medio del programa en Excel (hoja --), así mismo el límite superior de las
integrales de los Cp, que corresponde a la temperatura de flama adiabática, fue determinado
utilizando la función solver del Excel (hoja---):

T ad = 1256, 669K

Al comparar el resultado obtenido por el de la literatura [1], se tiene:

1256,669−1256,63
%Error = 1256,669
= 0.003

Con un error porcentual del 0,003%, se confirma que para este caso el modelo ideal propuesto
por Felder [1], no posee errores considerables, por lo que es válida la suposición que desarrolla
el ejemplo del texto.

Bibliografía
[1]. ​Richard M, Felder. Ronald W, Rousseau. Principios elementales de los procesos químicos.
3a Edición. Editorial Limusa WILEY. México.:1978

[2]​. Byung, Lee. Kesler. A generalized thermodynamic correlation based on three-Parameter

Corresponding states. Mobil research and development corporation. Princeton, New jersey.:
Estados Unidos.1875

[3]​. Borgnakke. Sonntag.Fundamentals of Thermodynamics.Primera edición.WILEY.

University of Michigan.:USA.2009

[4]​ ​G.Rayner-Canham. Química Inorgánica Descriptiva. 2° Edición. Editorial Prentice Hall. México.
2000.

[5] ​P. Atkins. Physical Chemistry. 6th Edition. W. H. Freeman and Company. New York. 1998.