Unidad 4 Analisis de Datos Cineticos11

Descargar como xlsx, pdf o txt
Descargar como xlsx, pdf o txt
Está en la página 1de 87

Para la Reaccion en Fase Liquida A + ⅟ 2 B → C, se obtivieron los siguientes datos de conversion en un reactor discontinuo

concentraciones iniciales de A y B fueron 0,1 y 0,05 gmol L -1, respectivamente. Usando el metodo integral, encuentre el orde
reaccion, la constante de velocidad de reaccion aparente y la energia de activacion.

Tiempo XA Pasos de Resolucion del Problema


XA T (K)
(min) @60min 0. Consideraciones
0 0 300 0.86 1. Ecuacion de Diseño
20 0.67 310 0.957 2. Modelo Cinetico
40 0.79 320 0.987 3. Estequiometria
60 0.86 4.
5. Combinar
Aplicar Metodo de Analisis de Datos
50000 0.9999 → Reaccion irreversible a este
"tiempo infinito" aun no alcanza
un equilibrio estable
CA0 (gmol/L)
0. Consideraciones CB0 (gmol/L)
Reaccion Homogenea (Fase Liquida) k' 0.1008
Reactor Discontinuo (Bacth) EA 100.33335
Mezcla Perfecta n 2
Metodo de Analisis Cinetico (Integral)

1. Ecuacion de Diseño 2. Modelo Cinetico


𝐹_𝐴0−𝐹_𝐴+∫1▒ 〖𝑟 _𝐴 𝑑𝑉 〗 = 〖𝑑𝑛〗 _𝐴/𝑑𝑡 𝑟_𝐴=−𝑘𝐶_𝐴^(𝛽_𝐴 ) 𝐶_𝐵^(𝛽_𝐵 ) (b)

rAV = dnA/dt (a) 3. Estequiometria


𝜉=( 〖−𝑛〗 _𝐴0 𝑋_𝐴)/𝑎=𝑛_𝐴0 𝑋_𝐴
𝑛_𝑗=𝑛_𝑗0+𝜐_𝑗 𝜉

A: nA = nA0 - nA0XA
B: nB = ⅟2 nA0 - ⅟2 nA0XA
C: nC = nA0XA
T: nT =

4. Combinar (a), (b), (c) y (d)

𝐴= 〖−𝑘 𝐶〗 _𝐴^(𝛽_𝐴 ) 𝐶_𝐵^(𝛽_𝐵 )=−𝑘 〖 𝐶〗 _𝐴0^(𝛽_𝐴 ) (⅟2 𝐶_𝐴0 )^(𝛽_𝐵 ) (1−𝑋_𝐴 )^(𝛽_𝐴 ) (1−𝑋_𝐴 )^(𝛽_𝐵
con)

𝐴/𝑑𝑡=−𝐶_𝐴0^(−1) (−𝑘)(𝐶_𝐴0^(𝛽_𝐴 ) ) (⅟2 𝐶_𝐴0 )^(𝛽_𝐵 ) (1−𝑋_𝐴 )^(𝛽_𝐴+𝛽_𝐵 )


〖𝑑𝑋〗 _𝐴/𝑑𝑡=𝑘′(1−𝑋_𝐴 )^𝑛 (f) 𝑘^′=𝑘 (𝐶_𝐴0^((𝑛−1)))/2^

5. Aplicar Metodo de Analisis de Datos


∫24_0^(𝑋_𝐴)▒ 〖 (𝑑𝑋_𝐴)/(1−𝑋_𝐴 )^𝑛 =𝑘′∫24_0^𝑡▒𝑑𝑡 〗
@n=0
Separando variables e integrando

XA = k't

XA vs t
1 Como la grafica de XA vs t no
f(x) = 0.0135 x + 0.175 Recta : n ≠ 0
0.9 R² = 0.78051391862955
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 10 20 30 40 50 60

@n=1 ∫24_0^(𝑋_𝐴)▒ 〖 (𝑑𝑋_𝐴)/(1−𝑋_𝐴 )^1 =𝑘′∫24_0^𝑡▒𝑑𝑡 〗

Separando variables e integrando

Ln( 1 - XA ) = k't
Ln( 1 - XA) vs t
0

Tiempo f(x) = − 0.0317516184643078 x − 0.206307253360545


XA Ln (1 - XA)
(min) -0.5 R² = 0.934106143180977
0 0 0
20 0.67 -1.108663 -1
40 0.79 -1.560648
60 0.86 -1.966113 -1.5

-2

-2.5
0 10 20 30 40 50

∫24_0^(𝑋_𝐴)▒ 〖 (𝑑𝑋_𝐴)/(1−𝑋_𝐴 )^2 =𝑘′∫24_0^𝑡▒𝑑𝑡 〗


@n=2 ∫24_0^(𝑋_𝐴)▒ 〖 (𝑑𝑋_𝐴)/(1−𝑋_𝐴 )^2 =𝑘′∫24_0^𝑡▒𝑑𝑡 〗

Separando variables e integrando


1/(1−𝑋_𝐴 )=1+𝑘^′ 𝑡

1/(1-XA) vs t
Tiempo 8
XA 1/(1-XA)
(min)
7
0 0 1 f(x) = 0.100800865800866 x + 0.95974025974026
6 R² = 0.996289762024474
20 0.67 3.030303
40 0.79 4.7619048 5
60 0.86 7.1428571
4

k' 0.1008 2

0
0 10 20 30 40 50 60

Ln (k') = Ln (A') + (EA/R)(1/T)

T (K) XA @60min k' 1/T vs ln (k`)


300 0.86 0.102381 0.5
310 0.957 0.3709302
0 f(x) = − 12068.1883361396 x + 37.9448676786456
320 0.987 1.2653846 R² = 0.999976949952974

-0.5
1/T (K )
-1
ln (k')
-1
0.00333333 -2.2790545956
0.00322581 -0.9917412865 -1.5
0.003125 0.2353761197
-2

-2.5
0.0031 0.00315 0.0032 0.00325 0.0033
Ln (A') 37.945
A' (min-1) 3.02E+16
-EA/R -12068
EA (kJ/mol) 100.333352
on en un reactor discontinuo. Las
do integral, encuentre el orden de

de Resolucion del Problema


raciones
n de Diseño
Cinetico
metria
ar
Metodo de Analisis de Datos

0.10
0.05

𝐴0 𝑋_𝐴)/𝑎=𝑛_𝐴0 𝑋_𝐴 𝐶_𝑗=𝑛_𝑗/𝑉 θB 0.5


V = V0

nA0 - nA0XA = nA0 ( 1 - XA ) → CA = nA/V = nA0( 1 - XA )/V0


⅟2 nA0 - ⅟2 nA0XA = ⅟2nA0 ( 1 - XA) → CB = nB/V = ⅟2 nA0 ( 1 - XA)/V0
nA0XA =
=

βA + βB = n

CA = CA0( 1 - XA )
dCA/dt = -CA0 dXA/dt
-CA0-1 dCA/dt = dXA/dt
𝑘^′=𝑘 (𝐶_𝐴0^((𝑛−1)))/2^(𝛽_𝐵 )

𝐴 )^𝑛 =𝑘′∫24_0^𝑡▒𝑑𝑡 〗

Como la grafica de XA vs t no es una linea


Recta : n ≠ 0

A) vs t

360545
Como la grafica de Ln(1 - XA) vs t no es una
linea Recta : n ≠ 1

40 50 60
vs t

Como la grafica de 1/(1 - XA) vs t es una


linea Recta : n = 2

40 50 60 70

(k`)

8676786456

0.00325 0.0033 0.00335


nA0( 1 - XA )/V0 = CA0( 1 - XA ) (c)
⅟2 nA0 ( 1 - XA)/V0 = ⅟2 CA0 ( 1 - XA) (d)
Determine el orden de reaccion de la descomposicion en fase gaseosa de Peroxido de di-terbutilo, utilizando el metodo difere
a partir de los datos presentados. Esta reaccion se efectuo en el laboratorio en un sistema Isotermico en el que se registro la P
Total en diversos momentos durante la reaccion. Inicialmente solo habia peroxido puro en el recipiente de reaccion.

Tiempo Presion Pasos de Resolucion de


(min) Total (Torr) A (g) Þ B (g) + 2C (g) 0. Consideraciones
0 7.5 1. Ecuacion de Diseño
5 12.5 2. Modelo Cinetico
10 15.8 3. Estequiometria
15 17.9 4.
5. Combinar
Aplicar Metodo de Analisis de
20 19.4

0. Consideraciones
Reaccion Homogenea (Fase Gaseosa)
Reactor Discontinuo (Bacth)
Mezcla Perfecta
Metodo de Analisis Cinetico (Diferencial)
DT = 0
PT vs t
PA0 = PT0

1. Ecuacion de Diseño 2. Modelo Cinetico


𝐹_𝐴0−𝐹_𝐴+∫1▒ 〖𝑟 _𝐴 𝑑𝑉 〗 = 〖𝑑𝑛〗 _𝐴/𝑑𝑡 𝑟_𝐴=−𝑘𝐶_𝐴^(𝛽_𝐴 )

rAV = dnA/dt (a) 3. Estequiometria CA = ƒ(PT)

𝑛_𝑗=𝑛_𝑗0+𝜐_𝑗 𝜉

𝐶_(𝐴 )= ( 〖 3𝑃 〗 _𝐴0−𝑃_𝑇)/2𝑅𝑇
A: nA =
(𝑑𝐶_𝐴)/𝑑𝑡=𝑑(3𝑃_𝐴0−𝑃_𝑇 )/(2𝑅𝑇 𝑑𝑡) B: nB =
C: nC =
(𝑑𝐶_𝐴)/𝑑𝑡=(−1/2𝑅𝑇) (𝑑𝑃_𝑇)/𝑑𝑡 T: nT =
(d)

Inicialmente solo habia (A) 𝜉=(𝑛_𝑇 〖−𝑛〗 _𝐴0)/


nA0 = nT0

4. Combinar (a), (b), (c) y (d) (𝑑𝑃_𝑇)/𝑑𝑡=𝑘(2𝑅𝑇)^(1−𝛽_𝐴 ) ( 〖 3𝑃 〗 _𝐴0−𝑃_𝑇 )^(𝛽_𝐴 )


〖𝑑𝑛〗 _𝐴/𝑑𝑡=𝑟_𝐴 𝑉 → 〖𝑑𝐶〗 _𝐴/𝑑𝑡=𝑟_𝐴

(𝑑𝑃_𝑇)/𝑑𝑡=𝑘′( 〖 3𝑃 〗 _𝐴0−𝑃_𝑇 )^(𝛽_𝐴 )


(−1/2𝑅𝑇) (𝑑𝑃_𝑇)/𝑑𝑡=−𝑘(( 〖 3𝑃 〗 _𝐴0−𝑃_𝑇)/2𝑅𝑇)^(𝛽_𝐴 )

con 𝑘^′=𝑘(2𝑅𝑇)^(1−𝛽_𝐴 )

𝑃_𝑇)/𝑑𝑡=𝑘(2𝑅𝑇/(2𝑅𝑇)^(𝛽_𝐴 ) ) ( 〖 3𝑃 〗 _𝐴0−𝑃_𝑇 )^(𝛽_𝐴 )

5. Aplicar Metodo de Analisis de Datos

𝐿𝑛((𝑑𝑃_𝑇)/𝑑𝑡)=𝐿𝑛(𝑘^′ )+𝛽_𝐴 𝑙𝑛( 〖 3𝑃 〗 _𝐴0−𝑃_𝑇 )

Tiempo Presion dPT/dt


n (min) Total (Torr) ΔPT/dt (3PA0 - PT)
Grafico Numerico Polinomio
0 0 7.5 1.2 1.17 1.2102 15
1 5 12.5 1.000 0.82 0.83 0.8087 10
2 10 15.8 0.660 0.53 0.54 0.5122 6.7
3 15 17.9 0.420 0.35 0.36 0.3207 4.6
4 20 19.4 0.300 0.26 0.24 0.2342 3.1

x y x y
0 1 0 1.2 Metodo Diferencial Grafico
5 1.4
1 2.5 1
5 0 7.5 0.66
5 1 12.5 0.42
5 1.2 0.66 17.5 0.3
f(x) = 0.00198635712238464 x² − 0.0868027515047292 x + 1.20112353109774
10 0.66
R² = 0.999784209615179 22.5 0.25
10 0
1
10 0.66
10 0.42
15 0.8 0.42
15 0
ΔPT/Δt

15 0.42
0.6
15 0.3
20 0.3
20 0.4 0
20 0.3

0.2
0.4

0.2

0
0 5 10 15 20

tiempo (min)

Modelo Diferencial Polinomico


25

20
f(x) = 0.00073333333333334 x³ − 0.0454285714285716 x² + 1.2102380952381 x + 7.49857142857143
R² = 0.999998416007198

15

10

0
0 5 10 15 20
tilo, utilizando el metodo diferencial,
ermico en el que se registro la Presion
ecipiente de reaccion.

Pasos de Resolucion del Problema


. Consideraciones
. Ecuacion de Diseño
. Modelo Cinetico
. Estequiometria
.. Combinar
Aplicar Metodo de Analisis de Datos

(b)

𝐶_𝑗=𝑛_𝑗/𝑉 V = V0

nA0 - x = nA0 - ½(nT - nA0) = nA0 + ½nA0 -½nT


nB0 + x ((3𝑃_𝐴0)/𝑅𝑇−𝑃_𝑇/𝑅𝑇)/2
CA (3nA0 - nT)/2V =
nC0 + 2x
nT0 + 2x = nA0 + 2x

𝜉=(𝑛_𝑇 〖−𝑛〗 _𝐴0)/2

Metodo De Analisis de Datos Cinetico (Diferencial)


0.4000

0.2000
1.07098943274328 x x−−2.56796787839409
1.01832339573626
f(x) = 0.994952734773053 2.70579902428351
2.51396026367995
0.996116451217882
R² = 0.998233600564242
Metodo De Analisis de Datos Cinetico (Diferencial)
0.4000

0.2000
1.07098943274328 x x−−2.56796787839409
1.01832339573626
f(x) = 0.994952734773053 2.70579902428351
2.51396026367995
0.0000 0.996116451217882
0.998233600564242
R² = 0.996466095953561

-0.2000

-0.4000
Ln (dPT/dt)
-0.6000

-0.8000

-1.0000

-1.2000

-1.4000

-1.6000
1.000 1.200 1.400 1.600 1.800 2.000 2.200

Ln ( 3PA0 - PT)

Ln (Graf) Ln (Num) Ln (Pol) Ln (3PA0 - PT)


βA ≡ 1.00
0.1823 0.1570 0.1908 2.7081 Ln (k') -2.7058
-0.1985 -0.1863 -0.2123 2.3026 k' 0.0668
-0.6349 -0.6162 -0.6690 1.9021
-1.0498 -1.0217 -1.1372 1.5261
-1.3471 -1.4271 -1.4516 1.1314

Metodo Nu
Punto inicial

Puntos Intermedios

( 〖𝑑𝑃〗 _𝑇/𝑑𝑡)_𝑡𝑖=(𝑃

Punto Final
( 〖𝑑𝑃〗 _𝑇/𝑑𝑡)_𝑡𝑓=(𝑃
20 25

Ecuacion Polinomica

PT a4 t3 a3 t2 a2 t a1
0.0007 -0.0454 1.2102 7.4986
Grado del polinomio n 3
Derivada del Polinomio

dPT/dt 3a4 t2 2a3 t a2


0.0021 -0.0908 1.2102

20 25
= /2nA0 - ½nT
3

−𝑃_𝑇/𝑅𝑇)/2 (3𝑃_𝐴0−𝑃_𝑇)/2𝑅𝑇
= (c)

s Cinetico (Diferencial)
s Cinetico (Diferencial)

2.000 2.200 2.400 2.600 2.800

A0 - PT)

Metodo Numerico
en terminos de PT Δt 5

𝑑𝑃〗 _𝑇/𝑑𝑡)_𝑡𝑖=(𝑃_(𝑇𝑖+1)−𝑃_(𝑇𝑖−1))/2Δ𝑡

𝑑𝑃〗 _𝑇/𝑑𝑡)_𝑡𝑓=(𝑃_(𝑇(𝑓−2))−4𝑃_𝑇(𝑓−1) +3𝑃_𝑇𝑓)/2Δ𝑡


Se han encontrado los siguientes datos para la
conversion de Acido Hidroxivalerico (A) a
Valerolactona (B)
A (aq) Þ B (aq) + C (aq)
Tiempo Pasos de Resolucion del Problema
CA (mol/L)
(min) 0. Consideraciones
0 19.04 1. Ecuacion de Diseño
48 17.6 2. Modelo Cinetico
76 16.9 3. Estequiometria
124 15.8 4.
5. Combinar
Aplicar Metodo de Analisis de Datos
204 14.41
238 13.94
289 13.37

0. Consideraciones
Reaccion Homogenea (Fase Liquida)
Reactor Discontinuo (Bacth)
Mezcla Perfecta
Metodo de Analisis Cinetico (Integral)

1. Ecuacion de Diseño 2. Modelo Cinetico


𝐹_𝐴0−𝐹_𝐴+∫1▒ 〖𝑟 _𝐴 𝑑𝑉 〗 = 〖𝑑𝑛〗 _𝐴/𝑑𝑡 𝑟_𝐴=−𝑘𝐶_𝐴^(𝛽_𝐴 )

rAV = dnA/dt (a) 3. Estequiometria

4. Combinar (a), (b), y (c) No aplica, los datos ya estan en terminos de la


variable a considerar
〖𝑑𝐶〗 _𝐴/𝑑𝑡=𝑟_𝐴= 〖−𝑘 𝐶〗 _𝐴^𝑛 con βA = n

∫_(𝐶_𝐴0)^(𝐶_𝐴)▒ 〖𝑑𝐶〗 _𝐴/(𝐶_𝐴^𝑛 )=−𝑘∫_0^𝑡▒𝑑𝑡


〖𝑑𝐶〗 _𝐴/(𝐶_𝐴^𝑛 )=−𝑘𝑑𝑡

5. Aplicar Metodo de Analisis de Datos (Integral)


@n=0
Separando variables e integrando

CA - CA0 = -kt
CA0 - CA= kt
(CA0 - CA) vs t
Tiempo (CA0 - CA) 6
f(x) = 0.0193633656616857 x + 0.466180716744246
(min) (mol/L) R² = 0.97713990753186
5
0 0
CA0 - CA (mol/L)

48 1.44 4

1
6
f(x) = 0.0193633656616857 x + 0.466180716744246
R² = 0.97713990753186
5

CA0 - CA (mol/L)
4
76 2.14 3
124 3.24
204 4.63 2
238 5.1 1
289 5.67
0
0 50 100 150 200 250 300
t (min)

@n=1
Separando variables e integrando
-(ln(CA)-ln(CA0) = kt
ln(CA0/CA) = kt
ln(CA0/CA) vs t
Tiempo 0.4
(min) ln (CA0/CA)
0.35 f(x) = 0.00122080101172592 x + 0.0190239021765612
0 0 R² = 0.989027713330823
0.3
48 0.0786431273 0.25
Ln (CA0/CA)

76 0.1192284074
0.2
124 0.1865320893
0.15
204 0.2786196194
238 0.311779624 0.1

289 0.3535286382 0.05


0
0 50 100 150 200 250 300
t (min)

@n=2

Separando variables e integrando

(CA)-1 - (CA0)-1 = kt

Tiempo ((CA)-1 - (CA0)-1) vs t


(CA)-1 - (CA0)-1
(min) 0.025
0 0
f(x) = 7.75789369906874E-05 x + 0.000590238381415203
48 0.0042971734 0.02 R² = 0.996438906388666
(CA0)-1 -(CA)-1 (L/mol)

76 0.0066505892
0.015
124 0.0107701308
204 0.0168752442
0.01
238 0.0192150031
289 0.0222733072 0.005

0
0 50 100 150 200 250 300
t (min)
(CA0)-1 -(C
0.01

0.005

0
0 50 100 150 200 250 300
t (min)

@n=3

Separando variables e integrando ((CA)-2 -(CA0)-2) vs t


0.0016
⅟2(C A
-2
-C A0 )= kt
-2
0.0014
f(x) = 4.96838808471404E-06 x + 8.99651049012396E-07

((CA)-2 -(CA0)-2) L2/mol2


0.0012 R² = 0.999535304910649
Tiempo 0.001
⅟2(CA-2 - CA0-2)
(min)
0.0008
0 0
0.0006
48 0.0002349247
0.0004
76 0.0003714108
0.0002
124 0.000623656
0
204 0.0010286918 0 50 100 150 200 250 300 350
238 0.0011937995
t (min)
289 0.0014178667
(b)

en terminos de la

〖𝑑𝐶〗 _𝐴/𝑑𝑡= 〖−𝑘 𝐶〗 _𝐴^𝑛

5. Aplicar Metodo de Analisis de

Tiempo
(min)

0
48
Como la grafica de (CA0 - CA)vs t no es una linea 76
Recta : n ≠ 0 124
204
238
289

0
48
48
250 300 48
48
76
76
76
76
124
124
124
124
Como la grafica de Ln(CA0/CA) vs t no es una linea 204
612 Recta : n ≠ 1 204
204
204
238
238
238
238
289
289
00 250 300
289

Como la grafica de (CA-1 - CA0-1) vs t no es una linea


Recta (tiene fluctuaciones por encima y por debajo
1415203 de la linea de tendencia)(argumentamos no
uniformidad de los datos) : n ≠ 2

00 250 300
00 250 300

2396E-07
Como la grafica de (CA-2 - CA0-2) vs t es una linea
Recta : n = 3

k (L2/mol2 min) 5.00E-06

250 300 350


0.000
2.550 2.600
-0.500
/𝑑𝑡= 〖−𝑘 𝐶〗 _𝐴^𝑛 𝐿𝑛( 〖−𝑑𝐶〗_𝐴/𝑑𝑡)=𝐿𝑛(𝑘)+𝑛 𝐿𝑛(𝐶_𝐴 ) -1.000
-1.500
-2.000

. Aplicar Metodo de Analisis de Datos (Diferencial) -2.500


-3.000

-dCA/dt Ln(dCA/dt) -3.500 f(x) =


CA (mol/L) ΔCA/Δt Ln(CA) R² =
-4.000
M. Graf M. Poly Ln(M. Graf) Ln(M. Poly) -4.500
19.04 0.034 3.22E-02 -3.381 -3.436 2.947 -5.000
17.6 0.03 0.0262 2.73E-02 -3.642 -3.601 2.868
16.9 0.025 0.024 2.46E-02 -3.730 -3.703 2.827
15.8 0.02291667 0.021 2.05E-02 -3.863 -3.887 2.760
14.41 0.017375 0.0143 1.47E-02 -4.247 -4.220 2.668
13.94 0.01382353 0.0129 1.26E-02 -4.351 -4.371 2.635
13.37 0.01117647 0.01 1.00E-02 -4.605 -4.605 2.593

0.030000 0 0.034
0.030000 24 0.030000 Metodo Grafico
0.000000 0.036000
62 0.025000
0.030000 0.034000 100 0.022917
0.025000 0.032000 164 0.017375
0.025000 0.030000 221 0.013824
0.000000 263.5
0.028000 0.011176
0.025000 0.026000 289 0.01
0.022917 0.024000
0.022000
0.022917
0.020000
0.000000 0.018000
0.022917 0.016000
0.017375 0.014000
0.017375 0.012000
0.010000
0.000000
0.008000
0.017375
0.006000
0.013824 0.004000
0.013824 0.002000
0.000000 0.000000
0 50 100 150 200 250 300 350
0.013824
0.011176
0.011176
0.000000
0.011176
Chart Title
0.000
2.550 2.600 2.650 2.700 2.750 2.800 2.850 2.900 2.950 3.000
-0.500
-1.000
-1.500
-2.000
-2.500
-3.000
-3.500 f(x) = 3.30799629002703
3.26894104589214 x − 13.0936362608222
12.9864091959965
R² = 0.981027040039635
0.977455662043945
-4.000
-4.500
-5.000

M. Graf M. Poly
ln(k) -13.094 ln(k) -12.986
k 2.05754E-06 k 2.292197E-06
n 3.3 n 3.26
Metodo Polinomico
20

19
f(x) = − 3.52912944569097E-08 x³ + 5.37052008477574E-05 x² − 0.0321872391973788 x + 19.0375605556
R² = 0.999984727772582
18

17
CA (mol/L)

16

15

14

13

12
0 50 100 150 200 250

t (min)

Ecuacion Polinomica

CA a 4 t3 a3 t 2 a2 t a1
-3.53E-08 5.37E-05 -3.22E-02 1.90E+01
Grado del polinomio n 2
Derivada del Polinomio

dCA/dt 3a4 t2 2a3 t a2


-1.06E-07 0.0001074 -0.0322
1872391973788 x + 19.0375605556972

200 250 300


V. N. Prusakov y V. B. Sokolov (Kinetics and Catalysis, 12, 24 (1971).) han presentado los siguientes datos para la descomposic
fase gaseosa de KrF2 a 70°C en un reactor rígido. Inicialmente solo está presente KrF2 en el reactor.

KrF2 → Kr + F2

Tiempo Presion Pasos de Resolucion de


(min) Total (Torr) A (g) Þ B (g) + C (g) 0. Consideraciones
0 40 1. Ecuacion de Diseño
3.5 50 2. Modelo Cinetico
6.25 56 3. Estequiometria
11.17 64 4.
5. Combinar
Aplicar Metodo de Analisis de
16.6 70

0. Consideraciones
Reaccion Homogenea (Fase Gaseosa)
Reactor Discontinuo (Bacth)
Mezcla Perfecta
Metodo de Analisis Cinetico (Diferencial)
DT = 0
PT vs t
PA0 = PT0

1. Ecuacion de Diseño 2. Modelo Cinetico


𝐹_𝐴0−𝐹_𝐴+∫1▒ 〖𝑟 _𝐴 𝑑𝑉 〗 = 〖𝑑𝑛〗 _𝐴/𝑑𝑡 𝑟_𝐴=−𝑘𝐶_𝐴^(𝛽_𝐴 )

rAV = dnA/dt (a) 3. Estequiometria CA = ƒ(PT)

𝑛_𝑗=𝑛_𝑗0+𝜐_𝑗 𝜉

𝐶_(𝐴 )= ( 〖 2𝑃 〗 _𝐴0−𝑃_𝑇)/𝑅𝑇
A: nA =
(𝑑𝐶_𝐴)/𝑑𝑡=𝑑(2𝑃_𝐴0−𝑃_𝑇 )/(𝑅𝑇 𝑑𝑡) B: nB =
C: nC =
(𝑑𝐶_𝐴)/𝑑𝑡=(−1/𝑅𝑇) (𝑑𝑃_𝑇)/𝑑𝑡 T: nT =
(d)

Inicialmente solo habia (A)


nA0 = nT0

4. Combinar (a), (b), (c) y (d) (𝑑𝑃_𝑇)/𝑑𝑡=𝑘(𝑅𝑇)^(1−𝛽_𝐴 ) ( 〖 2𝑃 〗 _𝐴0−𝑃_𝑇 )^(𝛽_𝐴 )


〖𝑑𝑛〗 _𝐴/𝑑𝑡=𝑟_𝐴 𝑉 → 〖𝑑𝐶〗 _𝐴/𝑑𝑡=𝑟_𝐴

(𝑑𝑃_𝑇)/𝑑𝑡=𝑘′( 〖 2𝑃 〗 _𝐴0−𝑃_𝑇 )^(𝛽_𝐴 )


−1/𝑅𝑇) (𝑑𝑃_𝑇)/𝑑𝑡=−𝑘(( 〖 2𝑃 〗 _𝐴0−𝑃_𝑇)/𝑅𝑇)^(𝛽_𝐴 )

con 𝑘^′=𝑘(𝑅𝑇)^(1−𝛽_𝐴 )

𝑃_𝑇)/𝑑𝑡=𝑘(𝑅𝑇/(𝑅𝑇)^(𝛽_𝐴 ) ) ( 〖 2𝑃 〗 _𝐴0−𝑃_𝑇 )^(𝛽_𝐴 )

5. Aplicar Metodo de Analisis de Datos

𝐿𝑛((𝑑𝑃_𝑇)/𝑑𝑡)=𝐿𝑛(𝑘^′ )+𝛽_𝐴 𝑙𝑛( 〖 2𝑃 〗 _𝐴0−𝑃_𝑇 )

dPT/dt Ln(dPT/dt)
Tiempo Presion
n ΔPT/Δt (2PA0 - PT)
(min) Total (Torr)
M. Graf M. Poly Ln(M. Graf)
0 0 40 40 3.2621 3.2481 1.182
1 3.5 50 2.857 30 2.4725 2.4687 0.905
2 6.25 56 2.182 24 1.9690 1.9697 0.678
3 11.17 64 1.626 16 1.3247 1.3260 0.281
4 16.6 70 1.105 10 0.9959 0.9865 -0.004

0 2.85714286 0
3.5 2.85714286 Histograma
1.75 2.85714286Presion Vs Tiempo
3.5
4 0 4.875 2.18181818
3.5 2.85714286 8.71 1.62601626
3.5 2.18181818
3.5 13.885 1.10497238
6.25 2.18181818 16.6
f(x) = 0.00679806905881176 x² − 0.249389111584565 x + 3.26209225508529
6.25 0
3 R² = 0.999272898238328
6.25 2.18181818
Presion Total (Torr)

2.5
6.25 1.62601626
11.17 1.62601626
2
11.17 0
11.17 1.62601626
1.5
11.17 1.10497238
16.6
1 1.10497238
16.6 0
16.6 1.10497238
0.5

0
0 2 4 6 8 10 12 14 16 18

Tiempo (min)
1

0.5

0
0 2 4 6 8 10 12 14 16 18

Tiempo (min)
ntes datos para la descomposición en
nte KrF2 en el reactor.

Pasos de Resolucion del Problema


. Consideraciones
. Ecuacion de Diseño
. Modelo Cinetico
. Estequiometria
.. Combinar
Aplicar Metodo de Analisis de Datos

(b)

𝐶_𝑗=𝑛_𝑗/𝑉 V = V0

nA0 - x = nA0 - (nT - nA0) = nA0 + nA0 -nT


nB0 + x (2𝑃_𝐴0)/𝑅𝑇−𝑃_𝑇/𝑅𝑇
CA (2nA0 - nT)/V =
nC0 + x
nT0 + x = nA0 + x

𝜉=𝑛_𝑇 〖−𝑛〗 _𝐴0


Analisis De Datos Cinetico Me
1.400

1.200
f(x) = 0.870217953030766
0.87343305146018 x x−−2.07463254530205
2.06196729482314
R² = 0.989306675186603
0.99093798289599
1.000

0.800

Ln (dPT/dt)
0.600

0.400

0.200

Ln(dPT/dt)
Ln(2PA0 - PT) 0.000
2.000 2.200 2.400 2.600 2.800
Ln(M. Poly)
1.178 3.689 -0.200
0.904 3.401
ln(2PA0 - P
0.678 3.178
0.282 2.773
-0.014 2.303

PT vs t
75

70
f(x) = 0.00217324794604577 x³ − 0.122897568987384 x² + 3.24812843906857 x + 40.0084456250831
R² = 0.999994402417439
65
Presion Total (Torr)

60

55

50

45

40

35
0 2 4 6 8 10 12 14 16

Tiempo (min)
45

40

35
0 2 4 6 8 10 12 14 16

Tiempo (min)

Ecuacion Polinomica

PT a 4 t3 a 3 t2 a2 t a1
0.0022 -0.1229 3.2481 40.008
Grado del polinomio n 3
Derivada del Polinomio

dPT/dt 3a4 t2 2a3 t a2


0.0066 -0.2458 3.2481
= 2nA0 - nT
(2𝑃_𝐴0−𝑃_𝑇)/𝑅𝑇
= (c)
s De Datos Cinetico Metodo Diferencial

0766
018 x x−−2.07463254530205
2.06196729482314
603
99

2.600 2.800 3.000 3.200 3.400 3.600 3.800

ln(2PA0 - PT)

M. Graf M. Poly
ln(k') -2.035 ln(k') -2.0746
k' 0.13068048 k' 0.12560666
n 0.86 n 0.87

k 0.04290121 k 0.04465031

𝑘^′=𝑘(𝑅𝑇)^(1−𝛽_𝐴 )

14 16 18
14 16 18
Se reportaron los siguientes datos para la descomposición en fase gaseosa a volumen constante de éter dimetílico a 504°
A→B+C+D
Inicialmente solo había el reactivo. Determine los parámetros cinéticos

Presion
Tiempo (s)
Total (Torr) A (g) Þ B (g) + C (g) + D (g)
0 310.3
390 408.0 T (ºC) 504
777 488.0 T(K) 777.15
1195 562.0
3155 700.0
∞ 931.0
P0 V0 = n0RT Pf Vf = nfRT
0. Consideraciones V0 = Vf
Reaccion Homogenea (Fase Gaseosa) P0/n0 = RT Pf/nf = RT
Reactor Discontinuo (Bacth)
Mezcla Perfecta P0/n0 = Pf/nf
Metodo de Analisis Cinetico (Diferencial) P0 = Pf (n0/nf)
DT = 0
PT vs t
PA0 = PT0

1. Ecuacion de Diseño 2. Modelo Cinetico


𝐹_𝐴0−𝐹_𝐴+∫1▒ 〖𝑟 _𝐴 𝑑𝑉 〗 = 〖𝑑𝑛〗 _𝐴/𝑑𝑡 𝑟_𝐴=−𝑘𝐶_𝐴^(𝛽_𝐴 )

rAV = dnA/dt (a) 3. Estequiometria CA = ƒ(PT)

𝑛_𝑗=𝑛_𝑗0+𝜐_𝑗 𝜉

𝐶_(𝐴 )= ( 〖 3𝑃 〗 _𝐴0−𝑃_𝑇)/2𝑅𝑇
A: nA =
(𝑑𝐶_𝐴)/𝑑𝑡=𝑑(3𝑃_𝐴0−𝑃_𝑇 )/(2𝑅𝑇 𝑑𝑡) B: nB =
C: nC =
(𝑑𝐶_𝐴)/𝑑𝑡=(−1/2𝑅𝑇) (𝑑𝑃_𝑇)/𝑑𝑡 D: nD =
(d)
T: nT =
Inicialmente solo habia (A)
𝜉=(𝑛_𝑇 〖−𝑛〗 _𝐴0)/
nA0 = nT0
4. Combinar (a), (b), (c) y (d) (𝑑𝑃_𝑇)/𝑑𝑡=𝑘(2𝑅𝑇)^(1−𝛽_𝐴 ) ( 〖 3𝑃 〗 _𝐴0−𝑃_𝑇 )^(𝛽_𝐴 )

〖𝑑𝑛〗 _𝐴/𝑑𝑡=𝑟_𝐴 𝑉 → 〖𝑑𝐶〗 _𝐴/𝑑𝑡=𝑟_𝐴

(𝑑𝑃_𝑇)/𝑑𝑡=𝑘′( 〖 3𝑃 〗 _𝐴0−𝑃_𝑇 )^(𝛽_𝐴 )


−1/2𝑅𝑇) (𝑑𝑃_𝑇)/𝑑𝑡=−𝑘(( 〖 3𝑃 〗 _𝐴0−𝑃_𝑇)/2𝑅𝑇)^(𝛽_𝐴 )

con 𝑘^′=𝑘(2𝑅𝑇)^(1−𝛽_𝐴 )

𝑃_𝑇)/𝑑𝑡=𝑘(2𝑅𝑇/(2𝑅𝑇)^(𝛽_𝐴 ) ) ( 〖 3𝑃 〗 _𝐴0−𝑃_𝑇 )^(𝛽_𝐴 )

5. Aplicar Metodo de Analisis de Datos

𝐿𝑛((𝑑𝑃_𝑇)/𝑑𝑡)=𝐿𝑛(𝑘^′ )+𝛽_𝐴 𝑙𝑛( 〖 3𝑃 〗 _𝐴0−𝑃_𝑇 )

Tiempo (s) Presion ΔPT/dt dPT/dt (3PA0 - PT) Ln (Graf) Ln (Pol)


Total (Torr) Grafico Polinomio
0 310.3 0.2710 0.2664 620.6667 -1.3056 -1.3228
390 408.0 0.2504 0.2285 0.2292 523.0000 -1.4760 -1.4731
777 488.0 0.2067 0.1896 0.1931 443.0000 -1.6628 -1.6444
1195 562.0 0.1770 0.1511 0.1551 369.0000 -1.8898 -1.8637
3155 700.0 0.0704 0.0200 0.0104 231.0000 -3.9121 -4.5636

x y
0 0.25042735 0 Metodo
0.2691 Grafico
390 0.35 0.25042735 195 0.25042735
390 0 583.5 0.20671835
390 0.3 0.25042735 986 0.17703349
390 0.20671835 2175 0.07040816
777 0.25 0.20671835
f(x) = 3155x² − 0.000112870135425844
1.06406989900394E-08 0.023 x + 0.271136211303716
R² = 0.998131531236532
777 0
777 0.2 0.20671835
777
ΔPT/Δt

0.17703349
11950.15 0.17703349
1195 0
1195 0.1 0.17703349
1195 0.07040816
31550.05 0.07040816

0
0 500 1000 1500 2000 2500 3000 3500

tiempo (s)
0.1

0.05
3155 0
3155 0 0.07040816
0 500 1000 1500 2000 2500 3000 3500

tiempo (s)
stante de éter dimetílico a 504°C

néticos

Pasos de Resolucion del Problema


0. Consideraciones
1. Ecuacion de Diseño
2. Modelo Cinetico
3. Estequiometria
4.
5. Combinar
Aplicar Metodo de Analisis de Datos

n0 1
nf 3

(b)

𝐶_𝑗=𝑛_𝑗/𝑉 V = V0

nA0 - x = nA0 - ½(nT - nA0) = nA0 + ½nA0 -½nT


nB0 + x ((3𝑃_𝐴0)/𝑅𝑇−𝑃_𝑇/𝑅𝑇)/2
CA (3nA0 - nT)/2V =
nC0 + x
nD0 + x
nT0 + 2x = nA0 + 2x
𝜉=(𝑛_𝑇 〖−𝑛〗 _𝐴0)/2
Analisis de Datos Cinetico Metodo
0.0000
5.2000 5.4000 5.6000 5.8000 6.0
-0.5000

-1.0000 f(x) = 3.29886113102107 x − 22.0571741927336


f(x)
R² ==0.859586149683926
2.65214160653302 x − 18.0348943713445
-1.5000 R² = 0.896302330431025

Ln (dPT/dt) -2.0000

-2.5000

-3.0000

-3.5000

-4.0000

-4.5000

-5.0000
Ln( 3PA0 - PT)

M. Graf M. Poly

Ln (3PA0 - PT) ln(k') -18.035 ln(k') -17.298


k' 1.4706151E-08 k' 3.0730813E-08
6.4308 n 2.652 n 2.5327
6.2596
6.0936 k 0.0236018637 k 0.01757486008
5.9108
5.4424

PT vs t
750.0 Ecua

700.0
PT
f(x) = 9.19072323162799E-10 x³ − 4.82227985914439E-05 x² + 0.266429365486779 x
+ 310.569511686638 Grado del polinomio
650.0 R² = 0.999981488692278
Deriva
600.0
dPT/dt
550.0

500.0

450.0

400.0

350.0

3500 300.0

250.0
-300 200 700 1200 1700 2200 2700 3200
450.0

400.0

350.0

3500 300.0

250.0
-300 200 700 1200 1700 2200 2700 3200
= /2nA0 - ½nT
3

(3𝑃_𝐴0−𝑃_𝑇)/2𝑅𝑇
= (c)
os Cinetico Metodo Diferencial
5.8000 6.0000 6.2000 6.4000 6.6000

741927336
943713445

Ln( 3PA0 - PT)

Ecuacion Polinomica
a 4 t3 a 3 t2 a2 t a1
9.19E-10 -4.82E-05 2.66E-01 3.10E+02
rado del polinomio n 3
Derivada del Polinomio
3a4 t2 2a3 t a2
2.757E-09 -0.00009644 0.2664
Gas Phase Reaction
A → 2B

Volumen
Tiempo (s)
(dm3) A (g) Þ 2B (g)
0 0.200
30 0.251 T (ºC) 40
60 0.276 T(K) 313.15

120 0.302 P (kPa) 303.39

240 0.322 yA 0.85


yI 0.15

0. Consideraciones
Reaccion Homogenea (Fase Gaseosa) (Gas id.)
Reactor Discontinuo (Bacth) (Tubo Capilar)
Mezcla Perfecta
Metodo de Analisis Cinetico (Diferencial)
DT = 0
DP = 0
VT vs t

1. Ecuacion de Diseño 2. Modelo Cinetico


𝐹_𝐴0−𝐹_𝐴+∫1▒ 〖𝑟 _𝐴 𝑑𝑉 〗 = 〖𝑑𝑛〗 _𝐴/𝑑𝑡 𝑟_𝐴=−𝑘𝐶_𝐴^(𝛽_𝐴 )

rA= (1/V)dnA/dt (a)

3. Estequiometria
𝐶_𝑗=𝑛_𝑗/𝑉
𝑛_𝑗=𝑛_𝑗0+𝜐_𝑗 𝜉

A: nA = nA0 - ξ = yA0nT0 - nT + nT0

I: nI = nI0
Gas ideal
𝑉/𝑉_0 =(𝑍/𝑍_0 )(𝑛_𝑇/𝑛_𝑇0 )(𝑇/𝑇_0 )(𝑃_0/𝑃)
B: nB = nB0 + 2ξ

T: nT = nT0 + ξ 𝑉=)𝑉_0 (𝑛_𝑇/𝑛_𝑇0 )


𝑉/𝑉_0 =(𝑛_𝑇/𝑛_𝑇0

ξ = nT -nT0 𝑛_𝑇0=(𝑃𝑉_0)/𝑅𝑇 𝑛_𝑇=𝑃𝑉/𝑅𝑇

𝑃/𝑅𝑇 (𝑉_0/𝑉 (𝑦_𝐴0+1)−1)


A: CA = nA/V =

𝑛_𝐴=𝑃/𝑅𝑇 (𝑉_0 (𝑦_𝐴0+1)−𝑉)

(𝑑𝑛_𝐴)/𝑑𝑡=−𝑃/𝑅𝑇 𝑑𝑉/𝑑𝑡 d)

4. Combinar a), b), c) y d)


1/𝑉 (−𝑃/𝑅𝑇 𝑑𝑉/𝑑𝑡)=−𝑘(𝑃/𝑅𝑇 (𝑉_0/𝑉 (𝑦_𝐴0+1)−1))^𝛼
1/𝑉 〖𝑑𝑛〗 _𝐴/𝑑𝑡=−𝑘𝐶_𝐴^(𝛽_𝐴 )

1/𝑉 𝑑𝑉/𝑑𝑡=𝑘 𝑅𝑇/𝑃 (𝑃/𝑅𝑇)^𝛼 (𝑉_0/𝑉 (𝑦_𝐴0+1)−1)^𝛼 1/𝑉 𝑑𝑉/𝑑𝑡=𝑘 〖𝑅𝑇〗 ^(1−𝛼) 𝑃^(𝛼−1) (𝑉_0/𝑉 (𝑦_𝐴0+1)−1)^𝛼

1/𝑉 𝑑𝑉/𝑑𝑡=𝑘′(𝑉_0/𝑉 (𝑦_𝐴0+1)−1)^𝛼 con 𝑘^′=𝑘 〖𝑅𝑇〗 ^(1−𝛼) 𝑃^(𝛼−1)

Aplicando Logaritmo en ambos lados

𝐿𝑛(1/𝑉 𝑑𝑉/𝑑𝑡)=𝐿𝑛(𝑘^′ )+𝛼𝐿𝑛(𝑉_0/𝑉 (𝑦_𝐴0+1)−1)

Volumen dVT/dt Ln (V-1dV/dt)


Tiempo (s) 1/V (dm-3) ΔV/Δt
(dm3) Grafico Polinomio Grafico Polinomio
0 0.200 5 0.0023 0.001931 -4.46540824 -4.64027936
30 0.251 3.98406375 0.0017 0.0012 0.00126603 -5.34313138 -5.28957008
60 0.276 3.62318841 0.00083333 0.00065 0.0007435 -6.05118378 -5.916782
120 0.302 3.31125828 0.00043333 0.00034 0.00012582 -6.78923668 -7.78336177
240 0.322 3.10559006 0.00016667 0.00005 0.00059986 -8.77028382 -6.28560386
0 0.0017 0 0.0023
30 0.0017 15 0.0017 Histograma V vs t
30 0
0.0025 45 0.00083333
30 0.0017 90 0.00043333
0.0023
30 0.00083333 180 0.00016667
60 0.0021
0.00083333 240 0.00005
60 0
0.0019
60 0.00083333
60 0.0017
0.00043333
120 0.00043333
0.0015
120 0
0.0013
120 0.00043333
120 0.00016667
0.0011
240 0.00016667
0.0009
240 0
240 0.00016667
0.0007

0.0005

0.0003

0.0000999999999999998

-0.0001 0 50 100 150 200

-0.0003
Pasos de Resolucion del Problema
0. Consideraciones
1. Ecuacion de Diseño
2. Modelo Cinetico
3. Estequiometria
4.
5. Combinar
Aplicar Metodo de Analisis de Datos

Reactores Intermitentes de Volumen Variable


Mas complejos que el reactor intermitente de flujo constante
Microprocesos → Reactor → Tubo Capilar con Particula Movil
Transcurso de la reaccion se sigue por el movimiento de la particula con resp
tiempo
Mas simple que intentar medir la composicion de la mezcla

V0 → Volumen Inicial Reactor


V → Volumen en el tiempo t

Reactores Intermitentes de Volumen Variable


ΔT = 0
ΔP = 0
Estequiometria Unica

(b)

= nT0(yA0 +1) - nT

(𝑃𝑉_0)/𝑅𝑇 (𝑦_𝐴0+1)−𝑃𝑉/𝑅𝑇
=
𝑃/𝑅𝑇 (𝑉_0 (𝑦_𝐴0+1)−𝑉)
=

𝑛_𝑇=𝑃𝑉/𝑅𝑇

/𝑉 (𝑦_𝐴0+1)−1)
c)

0/𝑉 (𝑦_𝐴0+1)−1)^𝛼

Analisis de Datos Cineticos Modelo Diferencial

f(x) = 0.563847222262097 x + 2.32341986518456


R² = 0.995140449792354
f(x) = 0.388556315020855 x + 1.42617732906916
R² = 0.880449407338046
ln(V-1dV/dt)

Ln(V0/V (yA0+1)-1)

-0.162518929497775 -8 -7.5 -7 -6.5 -6 -5.5 -5 -4.5


-0.746329446020254
Ln(V0/V (yA0+1)-1)
-1.07710608344715
-1.49091931219876
-1.90335053463652
Chart Title
0.330

f(x) = 2.63796727689145E-08 x³ − 1.22746364534801E-05 x² + 0.00193145206833262 x + 0.200837865311309


R² = 0.998611816181913
0.310

0.290

0.270

0.250

0.230

0.210

200 250

0.190
0 50 100 150 200
n Variable
nstante
a Movil
de la particula con respecto al

mezcla

n Variable
delo Diferencial
0 Metodo: Grafico Metodo: Polinomico
ln(k´) 2.3234 ln(k´) 1.4262
-0.2
x + 2.32341986518456 k' 10.2103305 k' 4.16285027
-0.4 n 0.56 n 0.38
-0.6
k( 3.96527108 k 1.09795076
-0.8

-1

-1.2

-1.4

-1.6
-5.5 -5 -4.5 -4

)-1)
Ecuacion Polinomica

VT a4 t
3
a 3 t2 a2 t a1
2.64E-08 -1.23E-05 1.93E-03 2.01E-01
Grado del polinomio n 3
06833262 x + 0.200837865311309 Derivada del Polinomio

dVT/dt 3a4 t2 2a3 t a2


7.914E-08 -0.00002454 0.001931

0 200 250
La Penicilina G (Reactivo A) reacciona irreversiblemente con Hidroxilamina NH2OH (Reactivo B) para formar
Acido Hioxamico. Para determinar el orden de reaccion global, se mezclaron concentraciones iguales de A y
B en un tanque de 250ml. Las muestras se retiraban cada 10 min y se adicionaban a una solucion de cloruro
Ferrico. El Acido Hidroxamico forma un complejo colorido con Fe(III). Con la ayuda de un colorimetro, se
obtuvo como una funcion del tiempo la concentracion del complejo colorido y por lo tanto la concentracion
de acido hidroxamico. La absorbancia es directamente proporcional a la concentracion del acido
hidroxamico.

Tiempo (min) Ab
0 0
A+BÞC
10 0.337
20 0.433
30 0.495
40 0.539
50 0.561
∞ 0.685

0. Consideraciones
Reaccion Homogenea (Fase Liquida)
Reactor Discontinuo (Bacth)
Mezcla Perfecta
Metodo de Analisis Cinetico (Diferencial)
Ci = ƒ(Abs)

1. Ecuacion de Diseño 2. Modelo Cinetico


𝐹_𝐴0−𝐹_𝐴+∫1▒ 〖𝑟 _𝐴 𝑑𝑉 〗 = 〖𝑑𝑛〗 _𝐴/𝑑𝑡 𝑟_𝐴=−𝑘𝐶_𝐴^𝛼 𝐶_𝐵^𝛽

rA = dCA/dt (a)

3. Estequiometria Ci = ƒ(Abs)

𝐶_𝑃=𝛾𝐴_𝑏

A: nA = nA0 - ξ = nA0 - nP

B: nB = nB0 - ξ = nB0 - nP

P: nP = nP0 + ξ = nP = ξ

nP0 = 0 Cj = nj/V
A: CA = nA/V = (nA0 - nP)/V

B: CB = nB/V = (nB0 - nP)/V

CA = CA0 - γAb → dCA/dt = -γdAb/dt

Evaluando la ecuacion en los limites

@t = 0 → Ab = 0 CP = 0

CA = CA0 - γ(0) → CA = CA0

@t = ∞ → Ab = Ab∞ CP = CP∞ CA = 0

0 = CA0 - γAb∞ → CA0 = γAb∞ → γ = CA0/Ab∞

4. Combinar

〖𝑑𝐶〗 _𝐴/𝑑𝑡=𝑟_𝐴=−𝑘𝐶_𝐴^𝛼 𝐶_𝐵^𝛽

−𝛾 〖𝑑𝐴〗 _𝑏/𝑑𝑡=−𝑘(𝐶_𝐴0−𝛾𝐴_𝑏 )^𝛼 (𝐶_𝐴0−𝛾𝐴_𝑏 )^𝛽 con n=α+β

𝛾 〖𝑑𝐴〗 _𝑏/𝑑𝑡=𝑘(𝐶_𝐴0−𝛾𝐴_𝑏 )^𝑛

𝐴_𝑏∞ 〖𝑑𝐴〗 _𝑏/𝑑𝑡=𝑘(( 〖𝐴 _𝑏∞ 𝐶 〗 _𝐴0−𝐶_𝐴0 𝐴_𝑏)/𝐴_𝑏∞ )^𝑛

𝐴〗 _𝑏/𝑑𝑡=𝑘 𝐴_𝑏∞/𝐶_𝐴0 (𝐶_𝐴0/𝐴_𝑏∞ )^𝑛 (𝐴_𝑏∞−𝐴_𝑏 )^𝑛

〖𝑑𝐴〗 _𝑏/𝑑𝑡=𝑘(𝐴_𝑏∞ )^(1−𝑛) (𝐶_𝐴0 )^(𝑛−1) (𝐴_𝑏∞−𝐴_𝑏 )^𝑛

〖𝑑𝐴〗 _𝑏/𝑑𝑡=𝑘′(𝐴_𝑏∞−𝐴_𝑏 )^𝑛 con 𝑘^′=𝑘(𝐴_𝑏∞ )^(1−𝑛) (𝐶_𝐴0 )^(𝑛−1)

ln⁡( 〖𝑑𝐴〗 _𝑏/𝑑𝑡)=ln⁡(𝑘′)+𝑛 ln⁡(𝐴_𝑏∞−𝐴_𝑏 )


5. Aplicar Metodo de Analisis de Datos

dAb/dt Ln (dAb/dt)

Tiempo (min) Ab ΔAb/Δt Grafico -


Numerico Polinomico Ln (Ab∞ Ab) Grafico
0 0 0.055 0.04575 0.0529 -0.37833644 -2.90042209
10 0.337 0.0337 0.02 0.02165 0.0182636 -1.0555528 -3.91202301
20 0.433 0.0096 0.0064 0.0079 0.0058008 -1.37832619 -5.05145729
30 0.495 0.0062 0.0051 0.0053 0.0049612 -1.66073121 -5.27851474
40 0.539 0.0044 0.002 0.0033 0.0051944 -1.92414866 -6.2146081
50 0.561 0.0022 0.0015 0.0011 -0.00405 -2.08747371 -6.50229017

x y x y
0 Histograma
0.0337 Absorbancia
0 0.055vs tiempo Aj
100.055 0.0337 5 0.0337 0.6
10 0.05 0 15 0.0096
f(x) = − 0.000000439583333 x⁴ + 0.0
100.045 0.0337 25 0.0062 + 0.00067857143
0.5 R² = 0.999470448735885
10 0.0096 35 0.0044
0.04
20 0.0096 45 0.0022
0.035 0.4
20 0 50 0.0015
Absorbancia

20 0.03 0.0096
200.025 0.0062 Δt 10 0.3
30 0.02 0.0062 1/(2Δt) 0.05
300.015 0
30 0.0062 0.2
0.01
30 0.0044
400.005 0.0044 0.1
40 0 0 0 10 20 30 40 50
40 0.0044
tiempo 0
40 0.0022 0 10 20
50 0.0022
50 0
50 0.0022
Pasos de Resolucion del Problema
0. Consideraciones
1. Ecuacion de Diseño
2. Modelo Cinetico
3. Estequiometria
4.
5. Combinar
Aplicar Metodo de Analisis de Datos

(b)
= CA0 - CP = CA0 - γAb

= CB0 - CP = CA0 - γAb

Alimentacion Equimolar CA0 = CB0

-2.2 -2 -1.8 -1.6 -1.4 -1.2 -1 -0.8 -0.6

f(x) = 2.15609953339323
2.04946121580272 x − 1.92762335254554
2.02262284715541
R² = 0.977646523845277
0.936089908014811
Ln (dAb/dt)

Numerico Polinomico Metodo: Grafico Metodo: Polinomico


-3.08456349 -2.93935194 ln(k´) -1.9276 ln(k´) -2.0226
-3.83274982 -4.00284527 k' 0.14549697 k' 0.13231101
-4.84089252 -5.14975944 n 2.2 n 2
-5.24004846 -5.30610763
-5.71383281 -5.26017416
-6.8124451 Err:502

Ajuste Polinomico Ecuacion Polinomica

Ab a 5 t4
-4.40E-07
(x) = − 0.000000439583333 x⁴ + 0.0000545416667 x³ − 0.00246229167 x² + 0.05290119048 x
+ 0.00067857143 Grado del polinomio
R² = 0.999470448735885
Derivada del Polinomio

dAb/dt 4a5 t3
-1.7584E-06

10 20 30 40 50 60
0
-1 -0.8 -0.6 -0.4 -0.2

-1

-2

-3

-4

-5

-6

-7

-8
-6

-7

-8

Ecuacion Polinomica
a 4 t3 a 3 t2 a2 t a1
5.45E-05 -2.46E-03 5.29E-02 6.79E-04
n 3
Derivada del Polinomio
3a4 t2 2a3 t a2
0.00016362 -0.004924 0.0529
Tiempo CA x103 CB0 (mol/dm3) 0.5
(min) (mol/dm3)
0 5.00E-02 Pasos de Resolucion del Problema
50 3.80E-02 1. Postular Una ley de Velocidad
100 3.06E-02 2. Procesar datos en terminos de la variable medida
150 2.56E-02 3. Buscar simplificaciones
200 2.22E-02 4.
5. Aplicar
Aplicar el Algortimo
Metodo CRE de Datos
de Analisis
250 1.95E-02
300 1.74E-02

1. Postular Una ley de Velocidad

−𝑟_𝐴=𝑘𝐶_𝐴^𝛼 𝐶_𝐵^𝛽

2. Procesar datos en terminos de la variable medida


Variable medida CA

3. Buscar simplificaciones
Como CB0 es 10 veces CA0, CB ≡ CB0

−𝑟_𝐴=𝑘𝐶_𝐴^𝛼 𝐶_𝐵0^𝛽

−𝑟_𝐴=𝑘 〖′𝐶〗 _𝐴^𝛼 con k' = kCβB0

4. Aplicar el Algortimo CRE

Balance de moles Combinando el Balance de moles con la Ley de velocidad


Reactor Intermitente: 〖𝑑𝑁〗 _𝐴/𝑑𝑡=𝑟_𝐴 𝑉

1/𝑉 〖𝑑𝑁〗 _𝐴/𝑑𝑡=𝑟_𝐴 (𝑑𝐶_𝐴)/𝑑𝑡=𝑟_𝐴

−𝑟_𝐴=𝑘 〖′𝐶〗 _𝐴^𝛼


Ley de Velocidad: −𝑟_𝐴=𝑘 〖′𝐶〗 _𝐴^𝛼
(𝑑 𝑁_𝐴⁄𝑉)/𝑑𝑡=𝑟_𝐴 −(𝑑𝐶_𝐴)/𝑑𝑡=𝑘′𝐶_𝐴^𝛼
Estequiometria: Liquido
V = V0
𝐶_𝐴=𝑁_𝐴/𝑉_0

Tomando Logaritmo natural en ambos lados


𝐿𝑛(−(𝑑𝐶_𝐴)/𝑑𝑡)=𝐿𝑛(𝑘^′ )+𝛼𝐿𝑛(𝐶_𝐴 )

5. Aplicar Metodo de Analisis de Datos

Metodo Grafico

Tiempo CA x103 -dCA/dt


-ΔCA/Δt
(min) (mol/dm3) Grafico Numerico Polinomico Ln(CA) Ln(Graf)
0 5.00E-02 3.00E-04 2.86E-04 2.98E-04 -2.99573227 -8.11172808
50 3.80E-02 2.40E-04 1.82E-04 1.94E-04 1.88E-04 -3.27016912 -8.61150387
100 3.06E-02 1.48E-04 1.20E-04 1.24E-04 1.19E-04 -3.48675527 -9.02801882
150 2.56E-02 1.00E-04 7.00E-05 8.40E-05 8.01E-05 -3.66516293 -9.56701532
200 2.22E-02 6.80E-05 5.80E-05 6.10E-05 6.03E-05 -3.80766299 -9.75506755
250 1.95E-02 5.40E-05 4.90E-05 4.80E-05 4.84E-05 -3.93734081 -9.92369026
300 1.74E-02 4.20E-05 4.50E-05 3.60E-05 3.34E-05 -4.05128507 -10.0088481
x y
0 0.00024 Histograma
0 CA0.0003
Vs t Concentracio
50 0.00035 0.00024 25 0.00024 5.50E-02
50 0 75 0.000148
5.00E-02
50 0.0003 0.00024 125 0.0001 f(x) = 3.6969697E-12 x⁴ − 3.48484848E-09 x³ + 0.00
R² = 0.999986967245954
4.50E-02
50 0.000148 175 0.000068
Concentracion de A (mol/dm3)

0.00025
100 0.000148 225 0.000054 4.00E-02
100 0 275 0.000042
0.0002 3.50E-02
100 0.000148 300 0.000038
-ΔCA/Δt

1000.00015 0.0001 3.00E-02


150 0.0001 Δt 50
2.50E-02
150 0.0001 0
150 0.0001 2.00E-02
1500.00005 6.80E-05
1.50E-02
200 6.80E-05
200 0 1.00E-02
0 50 0 100 150 200 250 300 350 0 50 100
200 0.000068
tiempo (min) Tie
200 0.000054
250 0.000054
250 0
250 0.000054
250 0.000042
300 0.000042
300 0
300 0.000042
moles con la Ley de velocidad

(𝑑𝐶_𝐴)/𝑑𝑡=𝑟_𝐴
𝑑𝐶_𝐴)/𝑑𝑡=𝑘′𝐶_𝐴^𝛼

Analisis de Datos Cinetico Modelo Diferencial

en ambos lados
^′ )+𝛼𝐿𝑛(𝐶_𝐴 )

f(x)
f(x) == 2.04906498733202
1.99585927365838 xx −− 1.92312963509051
1.90041551744449 2.09203136396756
2.44122699604109
R² = 0.995466122696827
0.994239216517507
R² = 0.983020939714896
Ln(-dCA/dt)

Ln (Numr) Ln (Polyn)
-8.15951875 -8.11892055
-8.5476524 -8.58019955
-8.99522899 -9.03717729
-9.38469376 -9.4320537
-4.3 -4.1 -3.9 -3.7 -3.5 -3.3 -3.1 -2
-9.70463669 -9.71581368
-9.94430955 -9.93523625 Ln (CA)
-10.2319916 -10.3078533
Ecuacion Polinomica
Concentracion vs Tiempo Orden Polinomio 4
2 CA a4t4 a3t3
3.70E-12 -3.48E-09
2
697E-12 x⁴ − 3.48484848E-09 x³ + 0.00000134348485 x² − 0.00029784632 x + 0.04999025974
dCA 4a4t3 3a3t2
6967245954
2 1.48E-11 -1.05E-08
2

2
0 50 100 150 200 250 300

Tiempo (min)
Diferencial
-7 Pendiente α ≈ 2
ln (k') -1.9231
-7.5 k' 0.14615318

-8

-8.5 Se predijo que la reaccion era de orden 1 para B, asi que podemos
despejar como CB ≈ CB0 :
-9
k' = kCβB0
-9.5
k = k'/CB0
-10 k (dm3/mol)2/min 0.29230637

-10.5 por lo que la Ley de Velocidad para la Reaccion estudiada seria

-11
-3.3 -3.1 -2.9
−𝑟_𝐴=[0,2923( 〖𝑑𝑚〗 ^3∕𝑚𝑜𝑙)^2/𝑚𝑖𝑛] 𝐶_𝐴^2 𝐶_𝐵

Ecuacion Polinomica

a2t2 a1t a0
1.34E-06 -2.98E-04 5.00E-02
2a2t a1
2.69E-06 -2.98E-04
tiempo CA 1/CA CB0 (mol/dm3) 0.5
(min) (mol/dm3) (dm3/mol)
0 0.05 20.00 Pasos de Resolucion del Problema
50 0.038 26.32 1. Postular Una ley de Velocidad
100 0.0306 32.68 2. Procesar datos en terminos de la variable medida
150 0.0256 39.06 3. Buscar simplificaciones
200 0.0222 45.05 4.
5. Aplicar
Aplicar el Algortimo
Metodo CRE de Datos
de Analisis
250 0.0195 51.28
300 0.0174 57.47

1. Postular Una ley de Velocidad

−𝑟_𝐴=𝑘𝐶_𝐴^𝛼 𝐶_𝐵^𝛽

2. Procesar datos en terminos de la variable medida


Variable medida CA

3. Buscar simplificaciones
Como CB0 es 10 veces CA0, CB ≡ CB0
−𝑟_𝐴=𝑘′𝐶_𝐴^𝛼 con k' = kCβB

4. Aplicar el Algortimo CRE

Balance de moles Combinando el Balance de moles con la Ley de velocidad


Reactor Intermitente: 〖𝑑𝑁〗 _𝐴/𝑑𝑡=𝑟_𝐴 𝑉

1/𝑉 〖𝑑𝑁〗 _𝐴/𝑑𝑡=𝑟_𝐴 (𝑑𝐶_𝐴)/𝑑𝑡=𝑟_𝐴

Ley de Velocidad: −𝑟_𝐴=𝑘 〖′𝐶〗 _𝐴^𝛼


(𝑑 𝑁_𝐴⁄𝑉)/𝑑𝑡=𝑟_𝐴 −(𝑑𝐶_𝐴)/𝑑𝑡=𝑘′𝐶_𝐴^𝛼
Estequiometria: Liquido
V = V0
𝐶_𝐴=𝑁_𝐴/𝑉_0 −(𝑑𝐶_𝐴)/𝑑𝑡=𝑘′𝐶_𝐴^2
−(𝑑𝐶_𝐴)/(𝐶_𝐴^2 )=𝑘′𝑑𝑡
1/CA vs t
70.00
4_(𝐶_𝐴0)^(𝐶_𝐴)▒(𝑑𝐶_𝐴)/(𝐶_𝐴^2 )=𝑘′∫24_0^𝑡▒𝑑𝑡
60.00
f(x) = 0.124794016573669 x + 20.1175244754766
50.00 R² = 0.999924169389169
(𝐶_𝐴)^(𝐶_𝐴0)▒(𝑑𝐶_𝐴)/(𝐶_𝐴^2 )=𝑘′∫24_0^𝑡▒𝑑𝑡

1/CA (dm3/mol)
40.00

30.00

20.00
1/𝐶_𝐴 −1/𝐶_𝐴0 =𝑘^′ 𝑡1/𝐶_𝐴 = 〖 1/𝐶_𝐴0 +𝑘 〗 ^′ 𝑡
10.00

0.00
0 50 100 150 200 250
tiempo (min)

k' (dm3/mol-min) 0.1248

k' = kCβB0
k = k'/CB0
k (dm3/mol)2/min 0.2496

−(𝑑𝐶_𝐴)/𝑑𝑡=𝑘′𝐶_𝐴^𝛼 α 2.0234285
k' 0.13590687

tiempo CA CAC
Error % (Cam-Cac)
(min) (mol/dm3) (mol/dm3)
0 0.05 0.05 1.92593E-32 0
50 0.038 0.03800331 7.585864E-09 -3.3097E-06
100 0.0306 0.03068299 7.355321E-06 -8.2989E-05
150 0.0256 0.02574611 3.25762E-05 -0.00014611
200 0.0222 0.02218933 2.309796E-07 1.06694E-05
250 0.0195 0.01950365 3.512011E-08 -3.6544E-06
300 0.0174 0.01740323 3.451605E-08 -3.2327E-06
4.023972E-05
moles con la Ley de velocidad

(𝑑𝐶_𝐴)/𝑑𝑡=𝑟_𝐴

𝑑𝐶_𝐴)/𝑑𝑡=𝑘′𝐶_𝐴^𝛼
α 2

𝐶_𝐴)/𝑑𝑡=𝑘′𝐶_𝐴^2
/CA vs t

175244754766

150 200 250 300 350


tiempo (min)

a0 20.1177787
a1 0.12479276
ca = a0 + a1t
tiempo CA 1/CA CAC (Cam-Cac)2
(min) (mol/dm3) (mol/dm3)
0 0.05 20 20.1177787 0.01387181
50 0.038 26.3157895 26.3574169 0.00173284
100 0.0306 32.6797386 32.597055 0.00683656
150 0.0256 39.0625 38.8366932 0.05098869
200 0.0222 45.045045 45.0763314 0.00097884
250 0.0195 51.2820513 51.3159696 0.00115046
300 0.0174 57.4712644 57.5556078 0.00711382
Suma 0.08267302
Tiempo α 2.04198525
Cam Cac Error % (Cam-Cac)
(min)
k 0.14532468
0 5.00E-02 5.00E-02 1.93E-32 -6.94E-18
50 3.80E-02 3.79E-02 4.02E-06 7.62E-05
100 3.06E-02 3.06E-02 4.67E-08 -6.61E-06
150 2.56E-02 2.57E-02 1.23E-05 -8.99E-05
200 2.22E-02 2.22E-02 4.12E-06 4.51E-05
250 1.95E-02 1.95E-02 3.38E-07 1.13E-05
300 1.74E-02 1.74E-02 7.14E-08 -4.65E-06
Suma %E2 2.09E-05

Tiempo α 2.0400515
CA tc Error % (tm-tc)
(min) k 0.14431005
0 5.00E-02 0.00E+00 0.00E+00 0.00E+00
50 3.80E-02 4.96E+01 9.52E+01 3.71E-01
100 3.06E-02 1.00E+02 1.70E+01 -1.26E-01
150 2.56E-02 1.51E+02 2.10E+03 -1.17E+00
200 2.22E-02 1.99E+02 9.28E+02 6.76E-01
250 1.95E-02 2.50E+02 1.08E+02 2.03E-01
300 1.74E-02 3.00E+02 5.42E+01 -1.28E-01
Suma %E2 3.30E+03
El proceso de obtención de sulfato de sodio, ocurre en varias etapas de reacción. En la etapa final se purifica el s
solución del producto final, teniendo en cuenta que la máxima concentración permisible en el producto (expresa
aporta el Na para el proceso, y se ha encontrado que en este Reactor 3 no está ocurriendo la eliminación requer
concentración de impurezas, por lo que se pide al Ingeniero de Procesos hacer los correctivos pertinentes.
Para eliminar las impurezas hasta el valor permitido, el Ing. de Procesos propone colocar un reactor CSTR de 5m
soda (NaOH, reactivo B) para que se produzca la reacción con el (NH4)2SO4 (reactivo A). De esta reacción se prod
apruebe su instalación. Por lo tanto, se decide estudiar en el Laboratorio de la empresa, en un beaker, la reacció
a la concentración requerida.
Para este análisis se toma una muestra líquida de la planta, a la salida del Reactor 3 y se lleva al Laboratorio. Los
DATOS INICIALES PAR CADA ENSAYO:
Cantidad de Solución: 300 mL Densidad: 1
Cantidad de NaOH al 42% agregada: 2,52 g Na2SO4%: 20.2
Contenido de NH3 en la muestra inicial: 966 ppm

ENSAYOS Y RESULTADOS:_x000B_Transcurrido el tiempo para tomar cada muestra, se retira una pequeña porción, se deti
agregando un ácido, se enfría y se determina el NH3 residual en la muestra líquida:

Analice los datos reportados por el Laboratorio para:


a ) Presentar una tabla de la concentración de cada reactivo y la conversión (C A, CB, X) en función del tiempo.
b ) Establecer si es una reacción de primer orden con respecto a cada especie como lo informa un estudio publicado en la
c ) literatura y cuál sería el valor de k.
d ) Definir si es posible instalar en la planta el CSTR de 5m 3 (Reactor 4) y cuál sería la conversión final para una carga de 65
e ) Calcular la dependencia con la temperatura. Esto es: determinar Ea y factor pre-exponencial (o de frecuencia)
d ) Si se decide operar a 90°C el Reactor 4, para incrementar la conversión, ¿Cuál sería el tiempo de residencia requerido p
conversión permisible de NH3 cuando se utiliza un PFR después del Reactor 3 vs. cuando se emplea un CSTR, para una carg

Resumen
Reaccion Homogenea (Fase Liquida)
Reactor Discontinuo (Bacth 4 Reactor) A + 2B Þ 2C + 2D + E
Reactivo A (NH4)2SO4
Reactivo B (NaOH)
Alimentacion de A expresado en NH 3
Determinar:
CA, CB, X
n, k, k0, EA
VCSTR, XFinal
VCSTR, VPFR Para 65 gpm
CA, CB, X en funcion del tiempo

〖𝑐𝑚〗 ^3 𝑠𝑜𝑙)×(1000 〖𝑐𝑚〗 ^3 𝑠𝑜𝑙)/(1 𝐿 𝑠𝑜𝑙)×(1 𝑚𝑜𝑙 〖𝑁𝐻〗 _3)/(17,031𝑔 〖𝑁𝐻〗 _3 )×(1 𝑚𝑜𝑙 𝑁)/(1 𝑚𝑜𝑙 〖𝑁𝐻〗 _3 )×(1 𝑚𝑜𝑙 ( 〖𝑁𝐻

=2,52 𝑔 𝑁𝑎𝑂𝐻×0,42×(1 𝑚𝑜𝑙 𝑁𝑎𝑂𝐻)/(39,997 𝑔 𝑁𝑎𝑂𝐻)×1/(300 𝑚𝐿 𝑠𝑜𝑙)×(1000 𝑚𝑙)/(1 𝐿 𝑠𝑜𝑙)=0,0882 (𝑚𝑜𝑙 𝑁𝑎𝑂𝐻)/(𝐿 𝑠𝑜𝑙)

60º C
No. Tiempo (min)
NH3 ppm g NH3 CA (mol/L) CB (mol/L) XA
0 0 966 0.966 0.03403206 0.0882 0.000
1 5 829 0.829 0.02920557 0.0786 0.142
2 10 766 0.766 0.02698608 0.0741 0.207
3 15 707 0.707 0.02490752 0.0700 0.268
4 20 647 0.647 0.02279373 0.0657 0.330
5 25 615 0.615 0.02166637 0.0635 0.363
6 30 567 0.567 0.01997534 0.0601 0.413

0. Consideraciones
Reaccion Homogenea (Fase Liquida)
Reactor Discontinuo (Bacth 4 Reactor)
Mezcla Perfecta
Densidad Constante (V = V0)

1. Ecuacion de Diseño
A es el reactivo limitante
𝐹_𝐴0−𝐹_𝐴+∫1▒ 〖𝑟 _𝐴 𝑑𝑉 〗 = 〖𝑑𝑛〗 _𝐴/𝑑𝑡

rA V = dnA/dt como V = V0

rA = dCA/dt a)

2. Modelo Cinetico

𝑟_𝐴=−𝑘𝐶_𝐴^𝛼 𝐶_𝐵^𝛽 b)

3. Estequiometria CB = f(CA)
𝐶_𝑗=𝑛_𝑗/𝑉
𝑛_𝑗=𝑛_𝑗0+𝜐_𝑗 𝜉

A: nA = nA0 - ξ
A: nA = nA0 - ξ

B: nB = nB0 - 2ξ

A: CA = (nA0 - ξ)/V = 𝐶_𝐴0−𝜉 ̇

𝐶_𝐵0−2𝜉 ̇
B: CB = (nB0 - 2ξ)/V =

4. Combinar a), b) y c)

〖𝑑𝐶〗 _𝐴/𝑑𝑡=𝑟_𝐴=−𝑘𝐶_𝐴^𝛼 𝐶_𝐵^𝛽

〖𝑑𝐶〗 _𝐴/𝑑𝑡=−𝑘(𝐶_𝐴 )^𝛼 (𝐶_𝐵0−2(𝐶_𝐴0−𝐶_𝐴 ))^𝛽

〖𝑑𝐶〗 _𝐴/𝑑𝑡=−𝑘(𝐶_𝐴 )^𝛼 (𝐶_𝐵0−2𝐶_𝐴0+2𝐶_𝐴con


)^𝛽 m = CB0 - 2CA0

〖𝑑𝐶〗 _𝐴/𝑑𝑡=−𝑘(𝐶_𝐴 )^𝛼 (𝑚+2𝐶_𝐴 )^𝛽

〖𝑑𝐶〗 _𝐴/((𝐶_𝐴 )^𝛼 (𝑚+2𝐶_𝐴 )^𝛽 )=−𝑘𝑑𝑡 si n = 2 entonces α = β = 1

resolviendo por fracciones parciales

1/(𝐶_𝐴 )(𝑚+2𝐶_𝐴 ) =𝐴/𝐶_𝐴 +𝐵/((𝑚+2𝐶_𝐴 ) )

{█(2𝐴+𝐵=0@𝐴𝑚=1)┤ 𝐴=1/𝑚 𝐵=−2(1/𝑚)=−2/𝑚

𝐶_𝐴)▒ 〖 (𝑑𝐶_𝐴)/𝐶_𝐴 −2/𝑚 ∫25_(𝐶_𝐴0)^(𝐶_𝐴)▒ 〖𝑑𝐶〗 _𝐴/(𝑚+2𝐶_𝐴 )=−𝑘𝑡 〗

1/𝑚 [ln⁡〖𝐶 _𝐴−ln⁡〖𝐶 _𝐴0 〗 〗 ]−1/𝑚 [ln⁡〖 (𝑚+2𝐶_𝐴 )− 〗 ln⁡(𝑚+2𝐶_𝐴0 ) ]=−𝑘𝑡

𝐶_𝐴/𝐶_𝐴0 ) ]−[ln⁡(((𝐶_𝐵0−2𝐶_𝐴0+2𝐶_𝐴 ))/((𝐶_𝐵0−2𝐶_𝐴0+2𝐶_𝐴0 ) )) ]=−𝑚𝑘𝑡

[ln⁡(𝐶_𝐴/𝐶_𝐴0 ) ]−[ln⁡(((𝐶_𝐵0−2(𝐶_𝐴0−𝐶_𝐴 )))/((𝐶_𝐵0 ) )) ]=−𝑚𝑘𝑡

[ln⁡(𝐶_𝐴/𝐶_𝐴0 ) ]−[ln⁡(𝐶_𝐵/𝐶_𝐵0 ) ]=−𝑚𝑘𝑡

[ln⁡(𝐶_𝐵/𝐶_𝐵0 ) ]−[ln⁡(𝐶_𝐴/𝐶_𝐴0 ) ]=𝑚𝑘𝑡


ln⁡((𝐶_𝐵 𝐶_𝐴0)/(𝐶_𝐵0 𝐶_𝐴 ))=𝑚𝑘𝑡

60º C
No. Tiempo (min)
NH3 ppm g NH3 CA (mol/L) CB (mol/L) XA
0 0 966 0.966 0.03403206 0.0882 0.000
1 5 829 0.829 0.02920557 0.0786 0.142
2 10 766 0.766 0.02698608 0.0741 0.207
3 15 707 0.707 0.02490752 0.0700 0.268
4 20 647 0.647 0.02279373 0.0657 0.330
5 25 615 0.615 0.02166637 0.0635 0.363
6 30 567 0.567 0.01997534 0.0601 0.413

80º C
NH3 ppm g NH3 CA (mol/L) CB (mol/L) XA
966 0.966 0.034032 0.0882 0.000
415 0.415 0.014620 0.0494 0.570
282 0.282 0.009935 0.0400 0.708
151 0.151 0.005320 0.0308 0.844
107 0.107 0.003770 0.0277 0.889
65 0.065 0.002290 0.0247 0.933
33 0.033 0.001163 0.0225 0.966

c) Definir si el Volumen de 5m3

𝐹_𝐴0−𝐹_𝐴+∫1▒ 〖𝑟 _𝐴 𝑑𝑉 〗 =0 v0 (gpm) 65
v0 (L/min) 246.1
V(m )
3
5
〖𝜐 ̇ _0 𝐶 〗 _𝐴0− 〖𝜐 ̇ _0 𝐶 〗 _𝐴+𝑟_𝐴 𝑉=0 V(L) 5000
〖𝜐 ̇ _0 𝐶 〗 _𝐴0− 〖𝜐 ̇ _0 𝐶 〗 _𝐴+(−0,2383𝐶_𝐴 𝐶_𝐵 )𝑉=0
𝐶_𝐴0=(0,966𝑔 𝑁𝐻_3)/(1000 𝑔 𝑠𝑜𝑙)×(1,2 𝑔 𝑠𝑜𝑙)/(1 〖𝑐𝑚〗 ^3 𝑠𝑜𝑙)×(1000 〖𝑐𝑚〗 ^3 𝑠𝑜𝑙)/(1 𝐿 𝑠𝑜𝑙)×(1 𝑚𝑜𝑙 〖𝑁𝐻〗 _3)/
CA (mol/L) ƒ(CA)
0.01949566 0
0,034 (𝑚𝑜𝑙 ( 〖𝑁𝐻〗 _3 )_2 〖𝑆𝑂〗 _4)/(𝐿 𝑠𝑜𝑙)×(2 𝑚𝑜𝑙 𝑁)/(1 𝑚𝑜𝑙 ( 〖𝑁𝐻〗 _3 )_2 〖𝑆𝑂〗 _4 )×(1 𝑚𝑜𝑙 〖𝑁𝐻〗 _3)/(1 𝑚𝑜𝑙 𝑁
NH3 ppm 572.86158413

d) Calcular EA y k(90ºC)

ln⁡〖 (𝑘_80/𝑘_60 )=−𝐸_𝐴/𝑅 [1/𝑇2−1/𝑇1] 〗

EA (J/mol) f(EA) EA (kJ/mol) RT


126608.066 -1.617328E-11 126.608 3019.2291
k90 1.4637416E-19

〖𝜐 ̇ _0 𝐶 〗 _𝐴0−𝜐 ̇_0 (𝐶_𝐴0 (1−𝑋_𝐴 ))+𝑟_𝐴 𝑉=0


v0 (gpm) 65
v0 (L/min) 246.1
Xa f(Xa)
0.42713842 0
ión. En la etapa final se purifica el sulfato de sodio obtenido, utilizando un CSTR a 80°C (Reactor 3). El objetivo de esta etap
n permisible en el producto (expresado como NH3) es de 20 ppm. Este proceso opera normalmente desde hace varios años e
tá ocurriendo la eliminación requerida del Sulfato de Amonio, ya que el NH3 en el producto final está entre 900 y 1200 ppm
er los correctivos pertinentes.
one colocar un reactor CSTR de 5m3 (Reactor 4) que se encuentra disponible en el Taller de Mantenimiento, y conectarlo en
eactivo A). De esta reacción se produce NH3 gaseoso, sulfato de sodio y agua. Para poder realizar este cambio en la planta (i
a empresa, en un beaker, la reacción entre los reactivos A y B a 60°C y 80°C, para determinar la expresión cinética, y poste
actor 3 y se lleva al Laboratorio. Los datos son los siguientes (Reacción tipo A + 2B → 2C + 2D +E) :

Densidad: 1,2 g/cm3


Na2SO4%: 20.2

Tiempo Contenido de NH3 (ppm)


No.
(min) 60º C
se retira una pequeña porción, se detiene la reacción 1 5 829
2 10 766
3 15 707
4 20 647
X) en función del tiempo. 5 25 615
lo informa un estudio publicado en la 6 30 567
conversión final para una carga de 65 gpm.
exponencial (o de frecuencia)
ría el tiempo de residencia requerido para alcanzar la
ando se emplea un CSTR, para una carga de 65gpm?

Datos Iniciales
Sistema Abierto, Fase Liquida
A + 2B Þ 2C + 2D + E Volumen de la solucion : 300mL
Densidad: 1,2 g/cm3
NH3 Inicial : 966 ppm
NaOH inicial: 2,52g al 42%
𝑙 𝑁)/(1 𝑚𝑜𝑙 〖𝑁𝐻〗 _3 )×(1 𝑚𝑜𝑙 ( 〖𝑁𝐻〗 _3 )_2 〖𝑆𝑂〗 _4)/(2 𝑚𝑜𝑙 𝑁)=0,034 (𝑚𝑜𝑙 ( 〖𝑁𝐻〗 _3 )_2 〖𝑆𝑂〗 _4)/(𝐿 𝑠𝑜𝑙)

𝑠𝑜𝑙)=0,0882 (𝑚𝑜𝑙 𝑁𝑎𝑂𝐻)/(𝐿 𝑠𝑜𝑙)


CB = CB0 -2CA0XA

80º C
NH3 ppm g NH3 CA (mol/L) CB (mol/L) XA
966 0.966 0.034032 0.0882 0.000
415 0.415 0.014620 0.0494 0.570
282 0.282 0.009935 0.0400 0.708
151 0.151 0.005320 0.0308 0.844
107 0.107 0.003770 0.0277 0.889
65 0.065 0.002290 0.0247 0.933
33 0.033 0.001163 0.0225 0.966
con 𝜉 ̇=𝐶_𝐴0−𝐶_𝐴
𝐶_𝐴0−𝜉 ̇

𝐶_𝐵0−2𝜉 ̇ 𝐶_𝐵0−2(𝐶_𝐴0−𝐶_𝐴 ) c)
=

60ºC
0.16
0.14 f(x) = 0.00476261635989551 x + 0.00763430420052433
R² = 0.991724437763569
0.12
CBCA0)/(CB0CA))

0.1
0.08
0.06
60ºC
0.16
0.14 f(x) = 0.00476261635989551 x + 0.00763430420052433
R² = 0.991724437763569
0.12

Ln((CBCA0)/(CB0CA))
0.1
0.08
0.06
0.04
0.02
0
Ln((CBCA0)/(CB0CA)) 0 5 10 15 20 25 30 35
0 tiempo (min)
0.0370432864881477
0.0579131266504783
0.0803396930803208
0.106690335928805 80ºC
0.122509561051502 2.5
0.149018843993445
2
Ln((CBCA0)/(CB0CA))

f(x) = 0.0639414094605404 x − 0.103980810100537


1.5 R² = 0.971082353616897
Ln((CBCA0)/(CB0CA))
0 1
0.264815523226556
0.440757870129343 0.5
0.803129515213391
1.04142423071643 0
0 5 10 15 20 25 30 35
1.42680353976011 tiempo (min)
2.00905164360715

^3 𝑠𝑜𝑙)/(1 𝐿 𝑠𝑜𝑙)×(1 𝑚𝑜𝑙 〖𝑁𝐻〗 _3)/(17,031𝑔 〖𝑁𝐻〗 _3 )×(1 𝑚𝑜𝑙 𝑁)/(1 𝑚𝑜𝑙 〖𝑁𝐻〗 _3 )×(1 𝑚𝑜𝑙 ( 〖𝑁𝐻〗 _3 )_2 〖𝑆𝑂〗 _4)/(2 𝑚𝑜𝑙 𝑁

〖𝑆𝑂〗 _4 )×(1 𝑚𝑜𝑙 〖𝑁𝐻〗 _3)/(1 𝑚𝑜𝑙 𝑁)×(17,031𝑔 〖𝑁𝐻〗 _3)/(1 𝑚𝑜𝑙 〖𝑁𝐻〗 _3 )×(1 𝐿 𝑠𝑜𝑙)/(1000 〖𝑐𝑚〗 ^3 𝑠𝑜𝑙)×(1 〖𝑐𝑚〗 ^3 𝑠𝑜𝑙)/(1
or 3). El objetivo de esta etapa es eliminar impurezas de sulfato de amonio [(NH4)2SO4] que están presente en la
ente desde hace varios años en la empresa, sin embargo, recientemente se cambió el proveedor de materia prima que
al está entre 900 y 1200 ppm. Esto implica que el producto (sulfato de sodio) no puede ser comercializado por la alta

antenimiento, y conectarlo en serie con el reactor existente (Reactor 3). Para esto se requiere adicionar al Reactor 4,
ar este cambio en la planta (i.e., conectar en serie el reactor de 5m 3) se necesita justificarlo ante la Gerencia para que
la expresión cinética, y posteriormente, determinar si el volumen de 5m 3 servirá para llevar el NH3 en el producto final
+E) :

𝟐𝑵𝑯_𝟑+𝑯_𝟐 〖𝑺𝑶〗 _𝟒→( 〖𝑵𝑯〗 _𝟒 )_𝟐 〖𝑺𝑶〗

Contenido de NH3 (ppm)


80º C Pasos de Resolucion del Problema
415 0. Consideraciones
282 1. Ecuacion de Diseño
151 2. Modelo Cinetico
107 3. Estequiometria
65 4.
5. Combinar
Aplicar Metodo de Analisis de Datos
33
mk (60ºC) 0.0048
m (60ºC) 0.0201425
k(60ºC) 0.23830213 L/(min - mol)

−𝑟_𝐴=0,2383𝐶_𝐴 𝐶_𝐵
30 35

mk (80ºC) 0.0639
m (80ºC) 0.0201425
k(80ºC) 3.17239713 L/(min - mol)

−𝑟_𝐴=3,1723𝐶_𝐴 𝐶_𝐵

30 35

〖𝑁𝐻〗 _3 )_2 〖𝑆𝑂〗 _4)/(2 𝑚𝑜𝑙 𝑁)=0,034 (𝑚𝑜𝑙 ( 〖𝑁𝐻〗 _3 )_2 〖𝑆𝑂〗 _4)/(𝐿 𝑠𝑜𝑙)

𝑚〗 ^3 𝑠𝑜𝑙)×(1 〖𝑐𝑚〗 ^3 𝑠𝑜𝑙)/(1,2 𝑔 𝑠𝑜𝑙)×(1000 𝑔 𝑠𝑜𝑙)/(0,966𝑔 𝑁𝐻_3 )=𝑝𝑝𝑚


presente en la
de materia prima que
cializado por la alta

cionar al Reactor 4,
a Gerencia para que
3 en el producto final

〗 _𝟒→( 〖𝑵𝑯〗 _𝟒 )_𝟐 〖𝑺𝑶〗 _𝟒

También podría gustarte