Universidad

Michoacana de
San Nicolás de
Hidalgo

Facultad de Ingeniería Química

Ingeniería de Reactores Homogéneos

Tarea 3

Reactores Isotérmicos de Flujo Pistón

Profesor: Ing. Luis Nieto Lemus

Oros Ferreira Luis Fernando

1804857F

Morelia Michoacán a 13 de mayo de 2022

Datos:

Para un reactor discontinúo tendremos las siguientes características:

V=cte

T=100°C=cte

Ahora hablando del reactor en flujo pistón tendremos:

T=100°C=373.15K

P=1atm

Fao=100mol/hr

Za=0.95

Yao=0.8

𝒏𝒊 = 𝟎. 𝟐𝒏𝑻𝒐

Solución:

𝟏 𝒅𝒏𝑨 𝒅𝑪𝒂
𝒓𝑨 = − ( ) ∗ = −( ) = 𝒌𝟏𝑪𝑨
𝒗 𝒅𝒕 𝒅𝒕
𝒕𝑹 𝑪𝑨
𝒅𝑪𝑨
𝒌𝟏 ∫ 𝒅𝒕 = − ∫ ( )
𝟎 𝑪𝑨𝒐 𝑪𝑨

𝑪𝑨𝒐
𝒌𝟏 ∗ 𝒕𝑹 = [𝐥𝐧(𝑪𝑨)]𝑪𝑨𝒐
𝑪𝑨 = 𝐥𝐧 ( )
𝑪𝑨
𝑪𝒐𝒏𝒔𝒊𝒅𝒆𝒓𝒂𝒏𝒅𝒐 𝒍𝒂 𝒐𝒑𝒆𝒓𝒂𝒄𝒊ó𝒏 𝒑𝒂𝒓𝒂 𝒈𝒂𝒔𝒆𝒔 𝒊𝒅𝒆𝒂𝒍𝒆𝒔:

𝑪𝑨𝒐 𝑷𝑨𝒐/𝑹𝑻 𝑷𝑨𝒐

𝒌𝟏 ∗ 𝒕𝑹== 𝐥𝐧 ( 𝑪𝑨 ) = 𝐥𝐧 ( 𝑷𝑨/𝑹𝑻 ) = 𝐥𝐧 ( 𝑷𝑨 )
𝟏 𝑷𝑨𝒐
𝒕𝑹 = ∗ 𝐥𝐧 ( )
𝒌𝟏 𝑷𝑨

Hacemos ajuste por mínimos cuadrados con Excel para ajustar la ecuación de la forma
lineal:
 𝒚 = 𝒂 + 𝒃𝒙
𝒚 = 𝒕𝑹
𝑷𝑨𝒐
𝒙 = 𝐥𝐧 ( )
𝑷𝑨
𝒂=𝟎
𝟏
𝒃= = 𝟏𝟎𝟐. 𝟐𝟒𝟏𝟏𝟕𝟏𝒔𝒆𝒈
𝒌𝟏
𝒓 = 𝟎. 𝟗𝟗𝟗𝟗𝟓𝟕𝟎𝟔
Corroborando que la reacción es de segundo orden determinamos k1:

𝟏
𝒌𝟏 = = 𝟗. 𝟕𝟖𝟎𝟕𝟗𝟔 ∗ 𝟏𝟎−𝟑 𝒔−𝟏
𝟏𝟎𝟐. 𝟐𝟒𝟏𝟏𝟕𝟏𝒔

Pasamos a la ecuación de diseño de un reactor en flujo pistón:

𝒁𝑨
𝑽𝑹 𝒅𝒁𝑨
=∫
𝑭𝑨𝒐 𝟎 𝒓𝑨
𝒁𝑨
𝒅𝒁𝑨
𝑽𝑹 = 𝑭𝑨𝒐 ∫
𝟎 𝒓𝑨
𝒁𝑨
𝒅𝒁𝑨
𝑽𝑹 = 𝑭𝑨𝒐 ∗ 𝒀𝑨𝒐 ∫
𝟎 𝒓𝑨
𝒅𝑽𝑹 𝑭𝑨𝒐 𝑭𝑻𝒐 ∗ 𝒀𝑨𝒐
= =
𝒅𝒁𝑨 𝒓𝑨 𝒓𝑨

Teniendo una cinética de primer orden:

𝟏 𝒅𝒏𝑨 𝒏𝑨 𝒏𝑨 ∗ 𝑷
𝒓𝑨 = − ( ) ∗ = 𝒌𝟏 ∗ 𝑪𝑨 = 𝒌𝟏 ( ) = 𝒌𝟏 ( )
𝒗 𝒅𝒕 𝑽 𝒏𝑻 ∗ 𝑹𝑻
𝒌𝟏 ∗ 𝑷 𝒏𝑨
𝒓𝑨 = ∗
𝑹𝑻 𝒏𝑻

Utilizamos la estequiometría de la reacción:

𝒏𝑨 = 𝒏𝑨𝒐(𝟏 − 𝒁𝒂) = 𝒚𝑨𝒐 ∗ 𝒏𝑻𝒐(𝟏 − 𝒁𝒂)

𝟏
𝒏𝑹 = ∗ (𝒀𝑨𝒐 ∗ 𝒏𝑻𝒐 ∗ 𝒁𝒂)
𝟐
𝟏
𝒏𝑺 = ∗ (𝒚𝑨𝒐 ∗ 𝒏𝑻𝒐 ∗ 𝒁𝒂)
𝟐
𝟏 𝟏
𝒏𝑻 = 𝒚𝑨𝒐 ∗ 𝒏𝑻𝒐(𝟏 − 𝒁𝒂) + ∗ (𝒀𝑨𝒐 ∗ 𝒏𝑻𝒐 ∗ 𝒁𝒂) + ∗ (𝒀𝑨𝒐 ∗ 𝒏𝑻𝒐 ∗ 𝒁𝒂)
𝟐 𝟐
𝒏𝑻 = 𝒏𝑻𝒐(𝒚𝑨𝒐 + 𝟎. 𝟐)
𝒌𝟏 ∗ 𝑷 𝒏𝑨 𝒌𝟏 ∗ 𝑷
𝒓𝑨 = ∗ = ∗ (𝒚𝑨𝒐 ∗ 𝒏𝑻𝒐(𝟏 − 𝒁𝒂))/(𝒏𝑻𝒐(𝒚𝑨𝒐 + 𝟎. 𝟐))
𝑹𝑻 𝒏𝑻 𝑹𝑻
Sustituimos lo obtenido en la ecuación de diseño:
𝒁𝑨 𝒁𝑨
𝒅𝒁𝑨 𝒅𝒁𝑨
𝑽𝑹 = 𝑭𝑨𝒐 ∫ = 𝑭𝑨𝒐 ∫
𝟎 𝒓𝑨 𝟎 𝒌𝟏 ∗ 𝑷 ∗ 𝒚𝑨𝒐 ∗ 𝒏𝑻𝒐(𝟏 − 𝒁𝒂)
𝑹𝑻 𝒏𝑻𝒐(𝒚𝑨𝒐 + 𝟎. 𝟐)
𝒁𝑨
𝑭𝑨𝒐 ∗ 𝑹𝑻 𝒅𝒁𝑨
𝑽𝑹 = ∫
𝒌𝟏 ∗ 𝒚𝑨𝒐 ∗ 𝑷 𝟎 𝒚𝑨𝒐 ∗ 𝒏𝑻𝒐(𝟏 − 𝒁𝒂)
𝒏𝑻𝒐(𝒚𝑨𝒐 + 𝟎. 𝟐)
 𝑭𝑨𝒐 ∗ 𝑹𝑻
𝑽𝑹 = ∗ (𝒚𝑨𝒐[𝐥𝐧(𝟏 − 𝒁𝑨)]𝟎𝒁𝑨 + 𝟎. 𝟐[𝐥𝐧(𝟏 − 𝒁𝑨)]𝟎𝒁𝑨 )
𝒌𝟏 ∗ 𝒀𝑨𝒐 ∗ 𝑷
𝑭𝑨𝒐 ∗ 𝑹𝑻
𝑽𝑹 = ∗ (𝒚𝑨𝒐[𝐥𝐧(𝟏)/(𝟏 − 𝒁𝑨)] + 𝟎. 𝟐[𝐥𝐧(𝟏)/(𝟏 − 𝒁𝑨)])
𝒌𝟏 ∗ 𝒀𝑨𝒐 ∗ 𝑷
𝑭𝑨𝒐 ∗ 𝑹𝑻
𝑽𝑹 = ∗ ([𝐥𝐧(𝟏)/(𝟏 − 𝒁𝑨)] ∗ (𝒚𝑨𝒐 + 𝟎. 𝟐)
𝒌𝟏 ∗ 𝒀𝑨𝒐 ∗ 𝑷
𝟏𝟎𝟎𝒎𝒐𝒍 𝟏𝒉𝒓 𝟎. 𝟎𝟖𝟐𝟎𝟔𝒂𝒕𝒎𝑳𝒕
∗ (𝟑𝟔𝟎𝟎𝒔𝒆𝒈)( )(𝟑𝟕𝟑. 𝟏𝟓𝑲)
𝒉𝒓 𝒎𝒐𝒍𝑲
𝑽𝑹 =
(𝟗. 𝟕𝟖𝟎𝟕𝟗𝟔 ∗ 𝟏𝟎−𝟑 ∗ 𝟎. 𝟖 ∗ 𝟏𝒂𝒕𝒎
∗ ([𝐥𝐧(𝟏)/(𝟏 − 𝟎. 𝟗𝟓)] ∗ (𝟎. 𝟖 + 𝟎. 𝟐)

𝑽𝑹 = 𝟏𝟒𝟐𝟖𝟔𝟒. 𝟔𝟗𝟐𝑳𝒕

Datos:

𝑃 = 2.5𝑏𝑎𝑟

4𝑘𝑚𝑜𝑙
𝐹𝐴𝑜 =
ℎ𝑟

𝑇 = 609.5𝐾

𝑘1 = 107ℎ𝑟 −1 𝑒𝑛 𝑠𝑒𝑛𝑡𝑖𝑑𝑜 𝑑𝑖𝑟𝑒𝑐𝑡𝑜

𝑘𝑝 = 2.5𝑏𝑎𝑟

𝑍𝐴 = 0.3

Solución:

De acuerdo a la ecuación de diseño del reactor de flujo pistón:

𝑍𝐴
𝑉𝑅 𝑑𝑍𝐴
=∫
𝐹𝐴𝑜 0 𝑟𝐴

Despejamos VR:
 𝑍𝐴
𝑑𝑍𝐴
𝑉𝑅 = 𝐹𝐴𝑜 ∫
0 𝑟𝐴

Siendo una reacción cinética elemental se tiene coincidencia de los coeficientes

estequiométricos con el orden de reacción:

𝑟𝐴 = 𝑘1𝐶𝐴 − 𝑘2𝐶𝐵 2

Asumiendo comportamiento de gas ideal:

𝑃𝐴 𝑃𝐵 2
𝑟𝐴 = 𝑘1𝐶𝐴 − 𝑘2𝐶𝐵 = 𝑟𝐴 = 𝑘1 ( ) − 𝑘2 ( ) = 𝑘𝑝𝑃𝐴 − 𝑘𝑝′ 𝑃𝐵 2
2
𝑅𝑇 𝑅𝑇

𝐾𝑃 = 𝑘𝑝/𝑘𝑝′

Siendo:

𝑘2
𝑘𝑝′ = ( )
(𝑅𝑇)2

0.844617𝑚𝑜𝑙 0.08314𝐿𝑡𝑏𝑎𝑟 2168.838823𝐿𝑡

𝑘2 = 𝑘𝑝′ ∗ (𝑅𝑇)2 = ( 2
) [( ) (609.5𝐾)] =
𝑏𝑎𝑟 𝐿𝑡ℎ𝑟 𝑚𝑜𝑙𝐾 𝑚𝑜𝑙ℎ𝑟

Pasamos a la cinética de la reacción:

𝑛𝐴
2
𝑛𝐵 2 𝑛𝐴 ′ 2
𝑛𝐵 2
𝑟𝐴 = 𝑘1𝐶𝐴 − 𝑘2𝐶𝐵 = 𝑘1 ( ) − 𝑘2 ( ) = 𝑘𝑝𝑃 ( ) − 𝑘𝑝 𝑃 ( )
𝑉 𝑉 𝑛𝑇 𝑛𝑇

𝑈𝑡𝑖𝑙𝑖𝑧𝑎𝑛𝑑𝑜 𝑙𝑎 𝑒𝑠𝑡𝑒𝑞𝑢𝑖𝑜𝑚𝑒𝑡𝑟í𝑎 𝑝𝑎𝑟𝑎 𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑟 𝑙𝑜𝑠 𝑚𝑜𝑙𝑒𝑠:

𝑛𝐴 = 𝑛𝐴𝑜(1 − 𝑍𝐴)

𝑛𝐵 = 𝑛𝐵𝑜 + 2𝑛𝐴𝑜 ∗ 𝑍𝐴 = 2𝑛𝐴𝑜 ∗ 𝑍𝐴

𝑛𝑇 = 𝑛𝐴𝑜(1 − 𝑍𝐴) + 2𝑛𝐴𝑜 ∗ 𝑍𝐴 = 𝑛𝐴𝑜(1 + 𝑍𝐴)

Sustituimos las ecuaciones para calcular los moles en la ecuación de velocidad de

reacción:

𝑛𝐴 ′ 2
𝑛𝐵 2 𝑛𝐴𝑜(1 − 𝑍𝐴) ′ 2
2𝑛𝐴𝑜 ∗ 𝑍𝐴 2
𝑟𝐴 = 𝑘𝑝𝑃 ( ) − 𝑘𝑝 𝑃 ( ) = 𝑘𝑝𝑃 ( ) − 𝑘𝑝 𝑃 ( )
𝑛𝑇 𝑛𝑇 𝑛𝐴𝑜(1 + 𝑍𝐴) 𝑛𝐴𝑜(1 + 𝑍𝐴)

𝑘1 1 − 𝑍𝐴 𝑍𝐴 2
𝑟𝐴 = ∗𝑃( ) − 4𝑘𝑝′ ∗ 𝑃2 ( )
𝑅𝑇 1 + 𝑍𝐴 1 + 𝑍𝐴

Sustituimos rA en la ecuación de diseño:

 𝑍𝐴
𝑑𝑍𝐴
𝑉𝑅 = 𝐹𝐴𝑜 ∫ 2 = 𝐹𝐴𝑜 ∗ 𝐼
0 𝑘1 1 − 𝑍𝐴 ′ ∗ 𝑃2 ( 𝑍𝐴 )
∗ 𝑃 ( ) − 4𝑘𝑝
𝑅𝑇 1 + 𝑍𝐴 1 + 𝑍𝐴

Evaluamos para:

𝑍𝐴 (1
+ 𝑍𝐴)2
𝐼=∫
0 𝑎 − 𝑏𝑍𝐴2

Hacemos cambio de variable para:

𝑢 = 1 + 𝑍𝐴 → 𝑍𝐴 = 𝑢 − 1

𝑑𝑢 = 𝑑𝑍𝐴

𝑍𝐴 (1 𝑍𝐴
+ 𝑍𝐴)2 (𝑢)2
𝐼=∫ =∫ 𝑑𝑢
0 𝑎 − 𝑏𝑍𝐴2 0 𝑎 − 𝑏(𝑢 − 1)
2

Y ahora:

𝑤 =𝑢−1→𝑢 =𝑤+1

𝑑𝑧 = 𝑑𝑢

𝑍𝐴 𝑍𝐴 (𝑤
(𝑢)2 + 1)2 1 𝑍𝐴 𝑍𝐴
2𝑤
𝐼=∫ 2
𝑑𝑢 = ∫ 2
𝑑𝑤 = − ( ) ∫ 𝑑𝑤 + ∫ 2
𝑑𝑤
0 𝑎 − 𝑏(𝑢 − 1) 0 𝑎 − 𝑏(𝑤) 𝑏 0 0 𝑎 − 𝑏𝑤

𝑍𝐴
𝑎 1
+ ( + 1) ∫ 2
𝑑𝑤
𝑏 0 𝑎 − 𝑏𝑤

𝐼𝑛𝑡𝑒𝑔𝑟𝑎𝑛𝑑𝑜 𝑦 𝑒𝑣𝑎𝑙𝑢𝑎𝑛𝑑𝑜 𝑙í𝑚𝑖𝑡𝑒𝑠 𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟𝑒𝑠 𝑝𝑜𝑟 𝑐𝑎𝑚𝑏𝑖𝑜 𝑑𝑒 𝑣𝑎𝑟𝑖𝑎𝑏𝑙𝑒:

𝑎 𝑎
1 1 𝑎 + 1 𝑍𝐴 + √
𝑏
𝐼 = − 𝑍𝐴 + ln ( 2
)+𝑏 (ln )
𝑏 𝑏 𝑎𝑏𝑍𝐴 𝑎 𝑎
2𝑏√ 𝑍𝐴 − √
𝑏 [ 𝑏]

𝑘1𝑃 (107ℎ𝑟 −1 )(2.5𝑏𝑎𝑟) 5.278859𝑚𝑜𝑙

𝑎= = =
𝑅𝑇 0.08314𝐿𝑡𝑏𝑎𝑟 𝐿𝑡ℎ𝑟
( ) (609.5𝐾)
𝑚𝑜𝑙𝐾

𝑘1𝑃 ′ 2
(107ℎ𝑟 −1 )(2.5𝑏𝑎𝑟) 0.844617𝑚𝑜𝑙
𝑏= + 4𝑘𝑝 𝑃 = + 4( ) (2.5𝑏𝑎𝑟 2 )
𝑅𝑇 0.08314𝐿𝑡𝑏𝑎𝑟 𝑏𝑎𝑟 2 𝐿𝑡ℎ𝑟
( ) (609.5𝐾)
𝑚𝑜𝑙𝐾

26.39428393𝑚𝑜𝑙
𝑏=
𝐿𝑡ℎ𝑟
Evaluando para I:

1
𝐼=− 0.3
26.39428393𝑚𝑜𝑙
𝐿𝑡ℎ𝑟
5.278859𝑚𝑜𝑙
1 𝐿𝑡ℎ𝑟
+ ln ( )
26.39428393𝑚𝑜𝑙 5.278859𝑚𝑜𝑙 26.39428393𝑚𝑜𝑙 2
∗ (0.3)
𝐿𝑡ℎ𝑟 𝐿𝑡ℎ𝑟 𝐿𝑡ℎ𝑟
5.278859𝑚𝑜𝑙 5.278859𝑚𝑜𝑙
𝐿𝑡ℎ𝑟 0.3 + √ 𝐿𝑡ℎ𝑟
26.39428393𝑚𝑜𝑙 + 1 26.39428393𝑚𝑜𝑙
𝐿𝑡ℎ𝑟 𝐿𝑡ℎ𝑟
+ (ln )
5.278859𝑚𝑜𝑙 5.278859𝑚𝑜𝑙
26.39428393𝑚𝑜𝑙 𝐿𝑡ℎ𝑟 𝐿𝑡ℎ𝑟
2 √ 0.3 − √ 26.39428393𝑚𝑜𝑙
𝐿𝑡ℎ𝑟 26.39428393𝑚𝑜𝑙
𝐿𝑡ℎ𝑟 [ 𝐿𝑡ℎ𝑟 ]

0.0939119𝐿𝑡ℎ𝑟
𝐼=
𝑚𝑜𝑙

Despejando este valor en la ecuación de diseño siguiente:

4000𝑚𝑜𝑙 0.0939119𝐿𝑡ℎ𝑟
𝑉𝑅 = 𝐹𝐴𝑜 ∗ 𝐼 = ( )( ) = 377.00319𝐿𝑡
ℎ𝑟 𝑚𝑜𝑙

Datos:

𝑇 𝑦 𝑃 = 𝑐𝑡𝑒𝑠

0.1𝑚𝑜𝑙
𝐹𝐴𝑜 =
𝑠𝑒𝑔

𝑇 = 777.15𝐾

𝑃 = 1𝑏𝑎𝑟

𝑍𝐴 = 0.6

Solución:

Iniciamos por la ecuación de diseño para un reactor de flujo pistón:

 𝑍𝐴
𝑉𝑅 𝑑𝑍𝐴
=∫
𝐹𝐴𝑜 0 𝑟𝐴

Y despejamos para el Volumen del reactor:

𝑍𝐴
𝑑𝑍𝐴
𝑉𝑅 = 𝐹𝐴𝑜 ∫
0 𝑟𝐴

1 𝑑𝑛𝐴 𝑛𝐴 𝑛𝐴𝑃
𝑟𝐴 = − ( ) ∗ ( ) = 4.3 ∗ 10−4 ∗ 𝐶𝐴 = 4.3 ∗ 10−4 ( ) = 4.3 ∗ 10−4 ( )
𝑉 𝑑𝑡 𝑉 𝑛𝑇𝑅𝑇

Y separamos los términos para cada operar más adelante con ellos:

(4.3 ∗ 10−4 )𝑃 𝑛𝐴
𝑟𝐴 = ∗
𝑅𝑇 𝑛𝑇

Y nos vamos a la estequiometría de la reacción:

𝑛𝐴 = 𝑛𝐴𝑜(1 − 𝑍𝐴)

𝑛𝐵 = 𝑛𝐴𝑜 ∗ 𝑍𝐴

𝑛𝐶 = 𝑛𝐴𝑜 ∗ 𝑍𝐴

𝑛𝐷 = 𝑛𝐴𝑜 ∗ 𝑍𝐴

𝑛𝑇 = 𝑛𝐴𝑜(1 − 𝑍𝐴) + 𝑛𝐴𝑜 ∗ 𝑍𝐴 + 𝑛𝐴𝑜 ∗ 𝑍𝐴

𝑛𝑇 = 𝑛𝐴𝑜 − 𝑛𝐴𝑜 ∗ 𝑍𝐴 + 3𝑛𝐴𝑜 ∗ 𝑍𝐴 = 𝑛𝐴𝑜 + 2𝑛𝐴𝑜 ∗ 𝑍𝐴 = 𝑛𝐴𝑜(1 + 2𝑍𝐴)

Ahora sustituimos los términos de moles en la ecuación de velocidad de reacción:

(4.3 ∗ 10−4 )𝑃 𝑛𝐴𝑜(1 − 𝑍𝐴)

𝑟𝐴 = ∗( )
𝑅𝑇 𝑛𝐴𝑜(1 + 2𝑍𝐴)

Y este término de velocidad de reacción se sustituye en la ecuación de Volumen de

reacción:

𝑍𝐴
𝑑𝑍𝐴
𝑉𝑅 = 𝐹𝐴𝑜 ∫
0
(4.3 ∗10−4 )𝑃 𝑛𝐴𝑜(1 − 𝑍𝐴)
∗( )
𝑅𝑇 𝑛𝐴𝑜(1 + 2𝑍𝐴)

𝑆𝑒𝑝𝑎𝑟𝑎𝑚𝑜𝑠 𝑙𝑎𝑠 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒𝑠:

𝑍𝐴
𝐹𝐴𝑜(𝑅𝑇) 1 + 2𝑍𝐴 𝐹𝐴𝑜(𝑅𝑇)
𝑉𝑅 = −4
∫ 𝑑𝑍𝐴 = 𝐼
(4.3 ∗ 10 )𝑃 0 1 − 𝑍𝐴 (4.3 ∗ 10−4 )𝑃
Evaluando y separando la integral para I:

𝑍𝐴 𝑍𝐴
1 𝑍𝐴
𝐼=∫ 𝑑𝑍𝐴 + 2 ∫ 𝑑𝑍𝐴
0 1 − 𝑍𝐴 0 1 − 𝑍𝐴

Hacemos cambio de variable:

𝑢 = 1 − 𝑍𝐴 → 𝑍𝐴 = 1 − 𝑢

𝑑𝑢 = −𝑑𝑍𝐴

𝑍𝐴 𝑧𝑎 𝑍𝐴
1 1
𝐼 = −∫ ( ) 𝑑𝑢 − 2 {∫ ( ) 𝑑𝑢 − ∫ 𝑑𝑢}
0 𝑢 0 𝑢 0

𝑍𝐴 𝑍𝐴
1
𝐼 = −3 ∫ ( ) 𝑑𝑢 + 2 ∫ 𝑑𝑢 = 3[ln(1 − 𝑍𝐴)]0𝑍𝐴 + 2[1 − 𝑍𝐴]0𝑍𝐴
0 𝑢 0

1
𝐼 = 3 ln ( ) − 2𝑍𝐴
1 − 𝑍𝐴

Ahora sustituimos en la ecuación para VR respecto a la variable I:

𝐹𝐴𝑜(𝑅𝑇) 1
𝑉𝑅 = (3 ln ( ) − 2𝑍𝐴)
(4.3 ∗ 10−4 )𝑃 1 − 𝑍𝐴

Y sustituimos todos los valores ya conocidos:

0.1𝑚𝑜𝑙 (777.15𝐾)0.08314𝐿𝑡𝑏𝑎𝑟
( ( )
𝐿𝑡 ) 𝑚𝑜𝑙𝐾 1
𝑉𝑅 = (3 ln ( ) − 2(0.6))
(4.3 ∗ 10−4 )(1𝑏𝑎𝑟) 1 − 0.6

𝑉𝑅 = 23276.102𝐿𝑡

Datos:

𝑇 𝑦 𝑃 = 𝑐𝑡𝑒𝑠
 100𝑚𝑜𝑙
𝐶𝐴𝑜 =
𝐿𝑡

100𝑚𝑜𝑙
𝐶𝑖𝑜 =
𝐿𝑡

−1
𝑘10.25 min∗

𝜃 = 349.43021𝑠𝑒𝑔

Solución:

𝑍𝐴
𝑉𝑅 𝑑𝑍𝐴
=∫
𝐹𝐴𝑜 0 𝑟𝐴

𝑍𝐴
𝑉𝑅 𝑑𝑍𝐴
=∫
𝐹𝑇𝑜𝑌𝐴𝑜 0 𝑟𝐴

𝑍𝐴
𝑉𝑅 𝑑𝑍𝐴
= 𝑌𝐴𝑜 ∫
𝐹𝑇𝑜 0 𝑟𝐴

𝑛𝐴 𝑛𝐴𝑜𝑃
𝑟𝐴 = 𝑘1𝐶𝐴 = 𝑘1 ( ) = 𝑘1 ( )
𝑉 𝑛𝑇𝑅𝑇

𝑘1𝑃 𝑛𝐴
𝑟𝐴 = ( )( )
𝑅𝑇 𝑛𝑇

Hacemos la relación de moles con la estequiometría:

𝑛𝐴 = 𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜(1 − 𝑍𝐴)

𝑛𝑅 = 2𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜 ∗ 𝑍𝐴

𝑛𝑖𝑜 = 0.5𝑛𝑇𝑜

𝑛𝑇 = 𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜(1 − 𝑍𝐴) + 2𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜 ∗ 𝑍𝐴 + 0.5𝑛𝑇𝑜

𝑛𝑇 = 𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜 − 𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜 ∗ 𝑍𝐴 + 2𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜 ∗ 𝑍𝐴 + 0.5𝑛𝑇𝑜

𝑛𝑇 = 𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜 + 𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜 ∗ 𝑍𝐴 + 0.5𝑛𝑇𝑜 = 𝑛𝑇𝑜(𝑌𝐴𝑜 + 𝑌𝐴𝑜 ∗ 𝑍𝐴 + 0.5)

Sustituimos en la ecuación de rA:

𝑘1𝑃 𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜(1 − 𝑍𝐴)

𝑟𝐴 = ( )( )
𝑅𝑇 𝑛𝑇𝑜(𝑌𝐴𝑜 + 𝑌𝐴𝑜 ∗ 𝑍𝐴 + 0.5)

Y esta relación la aplicamos a la ecuación de diseño del reactor:

 𝑍𝐴
𝑉𝑅 𝑑𝑍𝐴
= 𝑌𝐴𝑜 ∫
𝐹𝑇𝑜 0 (𝑘1𝑃 ) ( 𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜(1 − 𝑍𝐴)
)
𝑅𝑇 𝑛𝑇𝑜(𝑌𝐴𝑜 + 𝑌𝐴𝑜 ∗ 𝑍𝐴 + 0.5)

𝑉𝑅 𝑌𝐴𝑜 ∗ 𝑅𝑇 𝑍𝐴 𝑌𝐴𝑜 + 𝑌𝐴𝑜𝑍𝐴 + 0.5

= ∫ 𝑑𝑍𝐴
𝐹𝑇𝑜 𝑘1 ∗ 𝑃𝑌𝐴𝑜 0 1 − 𝑍𝐴

𝑍𝐴 𝑍𝐴
𝑉𝑅 1 𝑍𝐴
= 1/𝑘1((𝑌𝐴𝑜 + 0.5) ∫ 𝑑𝑍𝐴 + 𝑌𝐴𝑜 ∫ ( 𝑑𝑍𝐴)
𝑅𝑇 0 1 − 𝑍𝐴 1 − 𝑍𝐴
𝐹𝑇𝑜 ( 𝑃 ) 0

𝐼𝑔𝑢𝑎𝑙𝑎𝑛𝑑𝑜 𝑎 𝜃 𝑝𝑎𝑟𝑎 ℎ𝑎𝑐𝑒𝑟 𝑐𝑎𝑚𝑏𝑖𝑜 𝑑𝑒 𝑣𝑎𝑟𝑖𝑎𝑏𝑙𝑒:

𝑍𝐴 𝑍𝐴
1 𝑍𝐴
𝜃 = 1/𝑘1((𝑌𝐴𝑜 + 0.5) ∫ 𝑑𝑍𝐴 + 𝑌𝐴𝑜 ∫ ( 𝑑𝑍𝐴)
0 1 − 𝑍𝐴 0 1 − 𝑍𝐴

Donde:

𝑢 = 1 − 𝑍𝐴 → 𝑍𝐴 = 1 − 𝑢

𝑑𝑢 = −𝑑𝑍𝐴

𝑍𝐴 𝑍𝐴
1 1−𝑢
𝜃 = 1/𝑘1((𝑌𝐴𝑜 + 0.5) ∫ 𝑑𝑢 + 𝑌𝐴𝑜 ∫ ( 𝑑𝑢)
0 𝑢 0 𝑢

𝑍𝐴 𝑍𝐴
1
𝜃 = 1/𝑘1((−2𝑌𝐴𝑜 − 0.5) ∫ 𝑑𝑢 + 𝑌𝐴𝑜 ∫ 𝑑𝑢
0 𝑢 0

𝜃 = (1/𝑘1)((2𝑌𝐴𝑜 + 0.5)[𝑙𝑛𝑢]0𝑍𝐴 + 𝑌𝐴𝑜[𝑢]0𝑍𝐴 )

1 1
𝜃=( ) ((2𝑌𝐴𝑜 + 0.5) ∗ ln ( ) − 𝑌𝐴𝑜 ∗ 𝑍𝐴)
𝑘1 1 − 𝑍𝐴

Llegando a la ecuación:

1 1
𝑓(𝑍𝐴) = ( ) ((2𝑌𝐴𝑜 + 0.5) ∗ ln ( ) − 𝑌𝐴𝑜 ∗ 𝑍𝐴) − 𝜃
𝑘1 1 − 𝑍𝐴

Donde:

1 1 1𝑚𝑖𝑛
𝑓(𝑍𝐴) = ( ) ((2(0.5) + 0.5) ∗ ln ( ) − (0.5) ∗ 𝑍𝐴) − (349.43021) ( )
0.25−1 1 − 𝑍𝐴 60𝑠𝑒𝑔

1 1 1𝑚𝑖𝑛
𝑓(𝑍𝐴) = ( ) (ln ( ) − (0.5) ∗ 𝑍𝐴) − (349.43021) ( )
0.25−1 1 − 𝑍𝐴 60𝑠𝑒𝑔
Resolviendo por Polymath:

Donde ZA=0.7

Considerando gas ideal:

𝑛𝐴𝑜 𝑛𝐴𝑜𝑃 𝑃
𝐶𝐴𝑜 = = = 𝑌𝐴𝑜 ∗ ( )
𝑉𝑜 𝑛𝑇𝑜𝑅𝑇 𝑅𝑇

𝑃 𝐶𝐴𝑜
=
𝑅𝑇 𝑌𝐴𝑜

Y para CA:

𝑛𝐴 𝑛𝐴 ∗ 𝑃
𝐶𝐴 = =
𝑉 𝑛𝑇 ∗ 𝑅𝑇

Y para los moles:

𝑛𝐴 𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜(1 − 𝑍𝐴)

=
𝑛𝑇 𝑛𝑇𝑜(𝑌𝐴𝑜 + 𝑌𝐴𝑜 ∗ 𝑍𝐴 + 0.5)

Eliminando nTo y sustituyendo en CA:

𝑛𝐴 𝑌𝐴𝑜(1 − 𝑍𝐴)
=
𝑛𝑇 (𝑌𝐴𝑜 + 𝑌𝐴𝑜 ∗ 𝑍𝐴 + 0.5)

𝑛𝐴 ∗ 𝑃 𝑌𝐴𝑜(1 − 𝑍𝐴) 𝐶𝐴𝑜 𝐶𝐴𝑜(1 − 𝑍𝐴)

𝐶𝐴 = = ∗( )=
𝑛𝑇 ∗ 𝑅𝑇 (𝑌𝐴𝑜 + 𝑌𝐴𝑜 ∗ 𝑍𝐴 + 0.5) 𝑌𝐴𝑜 (𝑌𝐴𝑜 + 𝑌𝐴𝑜 ∗ 𝑍𝐴 + 0.5)

100𝑚𝑜𝑙 22.22222𝑚𝑜𝑙
𝐶𝐴 = (( )(1 − 0.7))/(0.5 + (0.5(0.7) + 0.5) =
𝐿𝑡 𝐿𝑡
Datos:

𝑇 𝑦 𝑃 = 𝑐𝑡𝑒𝑠

100𝑚𝑜𝑙
𝐶𝐴𝑜 =
𝐿𝑡

200𝑚𝑜𝑙
𝐶𝑖𝑜 =
𝐿𝑡

0.00525𝐿𝑡
𝑘2 =
𝑚𝑜𝑙𝑚𝑖𝑛

40.2299𝑚𝑜𝑙
𝐶𝐴 =
𝐿𝑡

𝜃 =¿ ? 𝑒𝑛 𝑠𝑒𝑔

Solución:

𝑛𝐴𝑜 𝑛𝐴𝑜(𝑃)
𝐶𝐴𝑜 = ( )=
𝑉𝑜 𝑛𝑇𝑜(𝑅𝑇)

Para la Yao:

𝑛𝐴𝑜 𝑛𝐴𝑜 (𝑛𝐴𝑜 ∗ 𝑃)/(𝑛𝑇𝑜 ∗ 𝑅𝑇)

𝑌𝐴𝑜 = ( )=( )=( )
𝑛𝑇𝑜 𝑛𝐴𝑜 + 𝑛𝑖𝑜 (𝑛𝐴𝑜 ∗ 𝑃)/(𝑛𝑇𝑜 ∗ 𝑅𝑇) + (𝑛𝑖𝑜 ∗ 𝑃)/(𝑛𝑇𝑜 ∗ 𝑅𝑇)

100𝑚𝑜𝑙
𝐶𝐴𝑜 𝐿𝑡 1
𝑌𝐴𝑜 = ( )=( )=
𝐶𝐴𝑜 + 𝐶𝑖𝑜 100𝑚𝑜𝑙 200𝑚𝑜𝑙 3
𝐿𝑡 + 𝐿𝑡

Teniendo en cuenta que para las fracciones iniciales:

1
𝑌𝑖𝑜 = 1 − 𝑌𝐴𝑜 = 1 − = 2/3
3

Basándonos en la ecuación de diseño para un reactor de flujo en pistón:

 𝑍𝐴
𝑉𝑅 𝑑𝑍𝐴
=∫
𝐹𝐴𝑜 0 𝑟𝐴

𝑍𝐴
𝑉𝑅 𝑑𝑍𝐴
=∫
𝐹𝑇𝑜𝑌𝐴𝑜 0 𝑟𝐴

𝑍𝐴
𝑉𝑅 𝑑𝑍𝐴
= 𝑌𝐴𝑜 ∫
𝐹𝑇𝑜 0 𝑟𝐴

Siendo una reacción irreversible de segundo orden:

2
𝑛𝐴 2 𝑛𝐴𝑃 2
𝑟𝐴 = 𝑘2𝐶𝐴 = 𝑘2 ( ) = 𝑘2 ( )
𝑉 𝑛𝑇 ∗ 𝑅𝑇

𝑃 2 𝑛𝐴 2
𝑟𝐴 = 𝑘2 ( ) ( )
𝑅𝑇 𝑛𝑇

𝑃𝑎𝑠𝑎𝑚𝑜𝑠 𝑎𝑙 𝑐á𝑙𝑐𝑢𝑙𝑜 𝑑𝑒 𝑚𝑜𝑙𝑒𝑠:

𝑛𝐴 = 𝑛𝐴𝑜(1 − 𝑍𝐴) = 𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜(1 − 𝑍𝐴)

4
𝑛𝑅 = ( ) ∗ (𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜 ∗ 𝑍𝐴)
10

2
𝑛𝑖𝑜 = ∗ 𝑛𝑇𝑜
3

4 2
𝑛𝑇 = 𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜(1 − 𝑍𝐴) + ( ) ∗ (𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜 ∗ 𝑍𝐴) + ∗ 𝑛𝑇𝑜
10 3
3 2
= 𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜 − ( ) 𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜 ∗ 𝑍𝐴 + ∗ 𝑛𝑇𝑜
5 3

3 2
𝑛𝑇 = 𝑛𝑇𝑜 (𝑌𝐴𝑜 − ∗ 𝑌𝐴𝑜 ∗ 𝑍𝐴 + )
5 3

Sustituimos estas ecuaciones en la ec. de velocidad de reacción:

2
𝑌𝐴𝑜 ∗ 𝑛𝑇𝑜(1 − 𝑍𝐴) 𝑃 2
𝑟𝐴 = 𝑘2 ( ) ∗( )
3 2 𝑅𝑇
𝑛𝑇𝑜 (𝑌𝐴𝑜 − ∗ 𝑌𝐴𝑜 ∗ 𝑍𝐴 + 3)
5

2
2
𝑌𝐴𝑜𝑃 1 − 𝑍𝐴
𝑟𝐴 = 𝑘2 ( ) ( )
𝑅𝑇 2 3
(𝑌𝐴𝑜 + 3) − ∗ 𝑌𝐴𝑜 ∗ 𝑍𝐴
5

Ahora sustituimos esta expresión en la ecuación de diseño:

 𝑍𝐴
𝑉𝑅 𝑑𝑍𝐴
= 𝑌𝐴𝑜 ∫ 2
𝐹𝑇𝑜 0 2
𝑌𝐴𝑜𝑃 1 − 𝑍𝐴
𝑘2 ( ) ( 2 3 )
𝑅𝑇
(𝑌𝐴𝑜 + 3) − ∗ 𝑌𝐴𝑜 ∗ 𝑍𝐴
5

2
2 3
𝑧𝑎 ((𝑌𝐴𝑜 + 3) − 5 ∗ 𝑌𝐴𝑜 ∗ 𝑍𝐴)
𝑉𝑅 1 𝑅𝑇 2
=( )( ) ∫ ∗ 𝑑𝑍𝐴
𝐹𝑇𝑜 𝑘2 ∗ 𝑌𝐴𝑜 𝑃 0 (1 − 𝑍𝐴)2

Teniendo el cambio de variable:

𝑢 = 1 − 𝑍𝐴 → 𝑍𝐴 = 1 − 𝑢

𝑑𝑢 = −𝑑𝑍𝐴

2
2 3
𝑧𝑎 ((𝑌𝐴𝑜 + 3) − 5 ∗ 𝑌𝐴𝑜(1 − 𝑢))
𝑉𝑅 1 𝑅𝑇 2
=( )( ) ∫ ∗ 𝑑𝑢
𝐹𝑇𝑜 𝑘2 ∗ 𝑌𝐴𝑜 𝑃 0 (𝑢)2

2
𝑧𝑎 ((𝑌𝐴𝑜
2 3 3
𝑉𝑅 1 𝑅𝑇 2 + 3 − 𝑌𝐴𝑜) ∗ ( )𝑌𝐴𝑜(𝑢))
=( )( ) ∫ 5 5 ∗ 𝑑𝑢
𝐹𝑇𝑜 𝑘2 ∗ 𝑌𝐴𝑜 𝑃 0 (𝑢)2

𝑉𝑅 1 𝑅𝑇 2 𝑧𝑎 ((𝑎 + 𝑏𝑢))2
=( )( ) ∫ ∗ 𝑑𝑢
𝐹𝑇𝑜 𝑘2 ∗ 𝑌𝐴𝑜 𝑃 0 (𝑢)2
1 𝑅𝑇 2 𝑧𝑎 𝑎2 + 2𝑎𝑏𝑢 + 𝑏 2 𝑢2
= −( )( ) ∫ ∗ 𝑑𝑢
𝑘2 ∗ 𝑌𝐴𝑜 𝑃 0 (𝑢)2

𝑉𝑅 1 𝑅𝑇 2 𝑧𝑎
1 𝑍𝐴
1 𝑍𝐴
= −( ) ( ) (𝑎2 ∫ 2
∗ 𝑑𝑢 + 2𝑎𝑏 ∫ 𝑑𝑢 + 𝑏 2
∫ 𝑑𝑢
𝐹𝑇𝑜 𝑘2 ∗ 𝑌𝐴𝑜 𝑃 0 (𝑢) 0 𝑢 0

𝑉𝑅 1 𝑅𝑇 2 1 𝑍𝐴
=− ∗ ( ) (−𝑎 [ ] + 2𝑎𝑏[𝑙𝑛𝑢]0𝑍𝐴 + 𝑏 2 [𝑢]0𝑍𝐴
2
𝐹𝑇𝑜 𝑘2𝑌𝐴𝑜 𝑃 𝑢 0

𝑆𝑢𝑠𝑡𝑖𝑡𝑢𝑖𝑚𝑜𝑠 𝑑𝑒 𝑛𝑢𝑒𝑣𝑜 𝑙𝑎 𝑣𝑎𝑟𝑖𝑎𝑏𝑙𝑒 1 − 𝑍𝐴 𝑒𝑛 𝑙𝑎 𝑣𝑎𝑟𝑖𝑎𝑏𝑙𝑒 𝑢:

𝑉𝑅 1 𝑅𝑇 2 2 1 1
= ∗ ( ) (𝑎 ( − 1) + 2𝑎𝑏 ∗ ln ( ) + 𝑏 2 𝑍𝐴
𝐹𝑇𝑜 𝑘2𝑌𝐴𝑜 𝑃 1 − 𝑍𝐴 1 − 𝑍𝐴

Con esto tendremos:

2 3 1 2 3 1 1 4
𝑎 = 𝑌𝐴𝑜 + − ∗ 𝑌𝐴𝑜 = + − ∗ ( ) = 1 − =
3 5 3 3 5 3 5 5

3 3 1 1
𝑏= ∗ 𝑌𝐴𝑜 = ∗ ( ) =
5 5 3 5
Y ahora con las fórmulas para las concentraciones:

𝑛𝐴 𝑛𝐴 ∗ 𝑃
𝐶𝐴 = =
𝑉 𝑛𝑇 ∗ 𝑅𝑇

𝑛𝐴𝑜 𝑛𝐴𝑜𝑃 𝑃
𝐶𝐴𝑜 = = = 𝑌𝐴𝑜 ∗ ( )
𝑉𝑜 𝑛𝑇𝑜𝑅𝑇 𝑅𝑇

Y con los moles:

𝑛𝐴 3 2
= (𝑌𝐴𝑜(1 − 𝑍𝐴))/(𝑌𝐴𝑜 − ( ) 𝑌𝐴𝑜 ∗ 𝑍𝐴 + )
𝑛𝑇 5 3

(𝑌𝐴𝑜(1 − 𝑍𝐴)) ∗ 𝐶𝐴𝑜

𝐶𝐴 =
3 2
(𝑌𝐴𝑜 − ( ) 𝑌𝐴𝑜 ∗ 𝑍𝐴 + 3) ∗ 𝑌𝐴𝑜
5

𝐶𝐴 ∗ 𝑌𝐴𝑜 3 2
= (𝑌𝐴𝑜(1 − 𝑍𝐴))/𝑌𝐴𝑜 − ( ) 𝑌𝐴𝑜 ∗ 𝑍𝐴 +
𝐶𝐴𝑜 5 3

𝐶𝐴 3 2
∗ (𝑌𝐴𝑜 − ( ) 𝑌𝐴𝑜 ∗ 𝑍𝐴 + ) = 1 − 𝑍𝐴
𝐶𝐴𝑜 5 3

𝐶𝐴 2 3 𝐶𝐴𝑌𝐴𝑜
(𝑌𝐴𝑜 + ) − 1 = ∗ ( ) 𝑍𝐴 − 𝑍𝐴
𝐶𝐴𝑜 3 5 𝐶𝐴𝑜

40.2299𝑚𝑜𝑙
𝐿𝑡 1 2
𝐶𝐴 2 (3 + 3) − 1
100𝑚𝑜𝑙
(𝑌𝐴𝑜 + 3) − 1
𝑍𝐴 = 𝐶𝐴𝑜 = 𝐿𝑡 = 0.64999986
3 𝐶𝐴𝑌𝐴𝑜 40.2299𝑚𝑜𝑙 1
∗ ( 𝐶𝐴𝑜 ) − 1 3 ( ) (3)
5 ∗( 𝐿𝑡 )−1
5 100𝑚𝑜𝑙
𝐿𝑡

Sustituyendo en la ec. de diseño:

𝑉𝑅 1 4 2 1 4 1
2 = 0.00525𝐿𝑡 ∗( ) ( − 1) + 2 ( ) ( )
1
𝑅𝑇
𝐹𝑇𝑜 ∗ ( 𝑃 ) ( ) (3) 5 1 − 0.64999986 5 5
𝑚𝑜𝑙𝑚𝑖𝑛
1 1 2
∗ ln ( ) + ( ) (0.64999986)
1 − 0.64999986 5

𝑉𝑅 886.007796(𝑚𝑜𝑙 ∗ min∗)
=
𝑅𝑇 2 𝐿𝑡
𝐹𝑇𝑜 ∗ ( 𝑃 )

𝑉𝑅 𝑉𝑅 𝑅𝑇 𝑉𝑅 𝑌𝐴𝑜
𝜃= ∗ = 2 ( )= 2 ( )
𝑅𝑇 𝑅𝑇 𝑃 𝑅𝑇 𝐶𝐴𝑜
𝐹𝑇𝑜 ∗ ( 𝑃 ) 𝐹𝑇𝑜 ∗ ( 𝑃 ) 𝐹𝑇𝑜 ∗ ( 𝑃 )
 1
𝑉𝑅 𝑌𝐴𝑜 886.00796(𝑚𝑜𝑙 ∗ min∗) 3
𝜃= 2 (𝐶𝐴𝑜 ) = ( 𝐿𝑡
)(
100𝑚𝑜𝑙
)
𝑅𝑇
𝐹𝑇𝑜 ∗ ( 𝑃 ) 𝐿𝑡
60𝑠𝑒𝑔
= 2.953359 min ( ) = 176.9365𝑠𝑒𝑔
1 min∗