MODELAMIENTO DEL REACTOR

 
En este sistema se lleva a cabo simultáneamente las reacciones elementales
mostradas: 

 
Estas cuatro reacciones pueden resumir el proceso de gasificación de carbón que
se lleva a cabo al interior del lecho fluidizado burbujeante. En la primer reacción se
considera la gasificación del carbón con vapor de agua, en la segunda considera la
gasificación de carbón con dióxido de carbono (Boudouard reaction) y finalmente
en la tercer reacción se presenta el equilibrio entre los reactivos y productos
principales involucrados en el proceso (water-gas shift reaction). [!] 
 

 
 
 
 
 
Los balances de materia y energía en el lecho fluidizado se llevaron a cabo
haciendo uso de la Figura ** ya que esta aclara fácilmente las corrientes que entran
y salen de este sistema. 
  
se realizaron balances de material por componente, con el fin de conocer la
concentración de cada componente a la salida del sistema. Inicialmente se definen
las velocidades de las reacciones representadas por (37), (38), y (39) suponiendo
un mecanismo de reacción elemental. Las expresiones (54) a (57) representan
estas velocidades expresadas en moles/segundo
Velocidades de reacción
r rxn 1 ( t )=K 1∗C char∗P H O ( t )
2

r rxn 2 ( t )=K 2∗C char∗PCO ( t ) 2

r rxn 3 ( t ) =K 3∗C char∗P H O (t )∗P CO ( t )

2

r rxn 3 } left (t right ) = {K} rsub {3∗C char∗PH O ( t )∗PCO ( t )

2 2

Balance por componentes

Balance de C char :

d C char ( t ) 1
= [F M (t)−FVTotal−out ( t )∗Cchar ( t )−r rxn 1 (t)−r rxn 2 (t)]
dt V lecho char
 d C char ( t ) 1
dt
= [ F M ( t ) −FVTotal−out ( t )∗Cchar ( t )−K 1 ( t )∗Cchar∗P H O ( t )−K 2∗C char∗PCO ( t ) ]
V lecho char 2 2

Linealizando :

d C D char ( t ) MD VD D D D D
=a 1 F char ( t )+ a2 F Total−out ( t ) +a 3 C char ( t ) +a 4 T ( t ) +a5 P H O ( t )+ a6 P CO ( t )
dt 2 2

Valor de las constantes :

1 −C ´char
a 1= a 2=
V lecho V lecho
´
− F VTotal−out
a 3=
V lecho
EA EA
−K 10 ¿C char∗PH´ O EA 1 T́
−K 20 ¿ C char∗P´CO EA 2
1 2

a 4= *( 2 ∗e )
2
*( ´ 2 ∗e T́ ) 2

V lecho T́ V lecho T
EA 1 EA 2
T́ T́
−K 10 ¿ C char∗e −K 20 ¿ C char∗e
a 5= a 6=
V lecho V lecho

´ K 10 ¿ C char∗P ´H O EA 1 EAT́
d C D char ( t ) C ´char F VTotal−out
[ K 20 ¿ C c 1
1
dt
=
V lecho
MD
F char ( t )−
V lecho
VD
∗F Total−out ( t ) −
V lecho
D
C char ( t )−
V lecho
2
∗ 2 ∗e +
T´ (
Vl )
EA
V lecho d C Dchar ( t ) C´char
[
K 10 ¿ C char∗ P ´H O EA 1 K 1
1
F V ´
Total−out
dt
+C Dchar ( t )= V
F ´ Total−out
F MD
char ( t )− V
´
F Total−out
VD
∗F Total−out ( t )− V ´
F Total−out
∗ 2 ∗e T́ +
T´
2

( )
Ecuación canonica :

d C D char ( t ) D MD VD D D D
𝜏1* +C char ( t )=K 1 F char ( t ) + K 2∗F Total−out ( t ) + K 3 T (t ) + K 4 P H O ( t )+ K 5 P CO ( t )
dt 2 2

Aplicandola transformadade Laplace:

K1 K2 K3 K4 K5
C D char ( S ) = F M Dchar ( S ) + ∗FV DTotal−out ( S ) + T D (S )+ P D H O ( S )+ P DCO ( S )
τ 1 S+1 τ 1 S +1 τ 1 S+ 1 τ 1 S+1 τ 1 S+1
2 2
 𝜏1 V lecho
F ´
V
Total−out

K1 1
F ´
V
Total−out

K2 −C ´char
FV ´
Total−out

K3 −K 10 ¿ C char∗P ´H O EA 1 EAT́ K 20 ¿ Cchar∗P´CO EA 2 EAT́

1

´
F VTotal−out
2

T´ (
∗ 2 ∗e −
´
FVTotal−out ) 2
∗ 2 ∗e
T´ (

K4 EA 1
T́
−K 10 ¿ C char∗e
FV ´ Total−out

K5 EA 2
T́
−K 20 ¿ C char∗e
FV ´ Total−out

Balance de C H O :
2

d CH O (t ) 1
F M ( t )−F VTotal ( t )∗C H O ( t )−r rxn1 ( t )+ r rxn3 ( t )−r rxn 3 ” ( t ) ]
V lecho [ H O
2
=
dt 2 out 2

d CH O (t )
2 1
= ¿
dt V lecho

Linealizando :
d CH O (t )
2
=b1 F M DH O (t ) +b 2 FV DTotal−out ( t )+ b3 C D H O ( t )+ b4 T D ( t )+ b5 PD H O ( t ) +b 6 P DCO ( t ) +b7 PD CO ( t )
dt 2 2 2 2

Valor de las constantes

1 −C ´H O
b 1= b 2= 2

V lecho V lecho
´
− F VTotal−out
b 3=
V lecho
EA EA
−K 10 ¿C char∗ PH´ O EA 1 T́
+ K 30 ¿ C char∗ PCO∗P
1 ´ H O EA
3
3

b 4= *( 2 ∗e )
2
*( ´ 2 ∗e T́ ¿ 2

V lecho T́ V lecho T
−K 30 } {*C} rsub {char} * acute {{P} rsub {{CO} rsub {2}} * {P} rsub {{H} rsub {2} O}}} over {{V} rsub {lecho} ¿ ¿ *(
EA 3 }} over {acute {{T} ^ {2}}} * {e} ^ {{{EA} rsub {3
¿)
T́
EA 1
T́
−K 10 ¿ C char∗e +
b 5=
V lecho
EA 3
T́
K 30 ¿ C char∗P´CO∗e K 30 } {*C} rsub {char} * {acute {acute {{P} rsub {CO2}}} *e} ^ {{{EA} rsub {3
− ¿ ¿
V lecho T́ V lecho
EA 3

K 30 ¿ C char∗P ´´H O∗e T́

b 6= 2

V lecho
−K 30} {*C} rsub {char} * acute {acute {{P} rsub {{H} rsub {2} O}}} * {e} ^ {{{EA} rsub {3 ¿
b 7= ¿
T́ V lecho

d CH O (t ) MD
−C H´ O V D ´
F VTotal−out K 10 ¿ C char∗ P ´H O
1
2
= F (t ) F 2
( t )− C D H O ( t )− * 2

dt V lecho H O V lecho
2 Total−out V lecho V lecho 2

EA ´ HO
EA 1 T́
K 30 ¿ C char∗PCO∗P
1

∗e + *( 2

T́ 2 V lecho
EA 3
EA 3 T́
∗e ¿−K 30 } {*C} rsub {char} * acute {{P} rsub {{CO} rsub {2}} * {P} rsub {{H} rsub {2} O}}} over {{V} rsub {lecho} ¿ ¿
T´ 2
 EA 1
EA 3 }} over {acute {{T} ^ {2}}} {e} ^ {{{EA} rsub {3 T́
*( ¿)T D ( t )− K 10 ¿ C char∗e +
T́ V lecho
EA 3

K 30 ¿ C char∗P´CO∗e T́
K 30 } {*C} rsub {char} * {acute {acute {{P} rsub {CO2}}} *e} ^ {{{EA} rsub {3 ¿ D
K 30 ¿ Cchar∗ PH´´ O∗e
¿ P ( t )+ 2

V lecho T́ V lecho H O 2
V lecho

´
d C H O ( t ) FVTotal−out MD
−C H´ O V D K 10 ¿C char∗ PH´ O
D 1
2
+ C H O ( t )= F (t ) F ( t )−
2
*( 2

dt V lecho V lecho
2
H O V lecho Total−out
2
V lecho
EA 1 EA
EA 1 ´ H O EA
+ K 30 ¿ C char∗ PCO∗P 3

∗e T́
) 2
*( 23 ∗e T́ ) -
T́ 2 V lecho T´
K 30 } {*C} rsub {char} * acute {{P} rsub {{CO} rsub {2}} * {P} rsub {{H} rsub {2} O}}} over {{V} rsub {lecho} ¿ ¿*(
EA 1
EA 3 }} over {acute {{T} ^ {2}}} {e} ^ {{{EA} rsub {3 T́
¿)T D ( t )− K 10 ¿ C char∗e +
T́ V lecho
EA3
´
K 30 ¿ C char∗ P´CO∗e T́
K 30 } {*C} rsub {char} * {acute {acute {{P} rsub {CO2}}} *e} ^ {{{EA} rsub {3 ¿ D
K 30 ¿C char∗P H´ O∗
− ¿ P ( t )+ 2

V lecho T́ V lecho H O 2
V lecho

Ecuación canonica :
𝜏2*
d CDH O (t )
+C D H O ( t )=K 6 F M char
2 D
( t )+ K 7∗F V DTotal−out ( t ) + K 8 T D ( t ) + K 9 P D H O ( t ) + K 10 P DCO (t ) + K 11 P DCO ( t )
dt 2 2 2

Aplicandola transformadade Laplace:

D D
D K6 MD K7 VD K8 D K9 D K 10 K 11
C H2 O (t )= F H O (t )+ ∗F Total−out ( t ) + T ( t )+ P ( t )+ P (t )+ P
τ 2 S +1 τ 2 S+1
2
τ 2 S+ 1 τ 2 S+1 H O τ 2 S+1 CO 2
τ 2 S+1 CO

𝜏2 V lecho
FV ´
Total−out

K6 1
F ´
V
Total−out
 K7 C H´ O
2

FV ´
Total−out

K8 −K 10 ¿ C char∗P ´H O EA 1 EAT́ 1

*( 2 ∗e )
2

´
F VTotal−out T́

+ K 30 ¿ C char∗PCO∗P ´ H O EA EA 3
3
*( 2 ∗e T́ )-
2

V ´
F Total−out T´
K 30 } {*C} rsub {char} * acute {{P} rsub {{CO} rsub {2}} * {P} rsub {{H} rsub {2} O}}} over {acu

EA 3 }} over {acute {{T} ^ {2}}} {e} ^ {{{EA} rsub {3

*( ¿)
T́

K9 EA 1
T́
−K 10 ¿ C char∗e +
FV ´ Total−out

EA 3
T́
K 30 ¿ C char∗P´CO∗e K 30 } {*C} rsub {char} * {acute {acute {{P} rsub {CO2}}} *e
−
FV ´ Total−out
T́

K 10 EA3

K 30 ¿ C char∗P ´´H O∗e 2

T́

F ´V
Total−out

K 11 −K 30 } {*C} rsub {char} * acute {acute {{P} rsub {{H} rsub {2} O}}} * {e} ^ {{{EA} rsub {3
¿
T́ V

Balance de C CO :
2

d C CO ( t ) 1 M
2
= [ F CO ( t )−F VTotal−out ( t )∗C CO ( t ) −r rxn 3 ( t )−r rxn3 ” ( t ) ]
dt V lecho 2 2
 d C CO ( t )
2 1
= ¿
dt V lecho

Linealizando :
d C CO ( t )
=b1 F MCO
D
2
( t ) +b 2 F V DTotal−out ( t ) +b3 C DCO ( t ) +b4 T D ( t ) +b5 PD H O ( t ) +b 6 P DCO ( t ) +b7 P DCO ( t )
dt 2 2 2 2

Valor de las constantes :

1 −C´CO
c 1= c 2= 2

V lecho V lecho
´
− F VTotal−out
c=
V lecho

´ H O EA EA 3
−K 30 ¿ C char∗PCO∗P 3 T́
c4 = *( ´ 2 ∗e 2
)
V lecho T
−K 30 } {*C} rsub {char} * acute {{P} rsub {{CO} rsub {2}} * {P} rsub {{H} rsub {2} O}}} over {{V} rsub {lecho} ¿ ¿ *(
EA 3 }} over {acute {{T} ^ {2}}} * {e} ^ {{{EA} rsub {3
¿)
T́
EA 3
T́
−K 30 ¿C char∗P´CO∗e K 30 } {*C} rsub {char} * {acute {acute {{P} rsub {CO2}}} *e} ^ {{{EA} rsub {3
c 5= − ¿ ¿
V lecho T́ V lecho
EA 3

−K 30 ¿ Cchar ∗PH´´ O∗e T́

c 6= 2

V lecho
−K 30 } {*C} rsub {char} * acute {acute {{P} rsub {{H} rsub {2} O}}} * {e} ^ {{{EA} rsub {3 ¿
c 7= ¿
T́ V lecho

−C´CO V D ´ EA
d C CO ( t ) 1 MD
F VTotal−out K 30 ¿ C char∗ PCO∗P
´ H O EA K 3

dt
2
=
V lecho
F
H O
(t )
2
V lecho
F
Total−out
2
(t)−
V lecho
D
C CO ( t )−
V lecho 2
∗ 2 ∗e T́ − 30 } {
T´
3 2

( )
´
d C CO ( t ) F VTotal−out 1 MD −C´CO V D K 30 ¿ C char∗PCO∗P´ H O EA 3 EA K 3

dt
2
−
V lecho
D
C CO ( t )=
V lecho
2
F
H O
(t )
V lecho
F
Total−out
2
( t )− 2

V lecho
∗ 2 ∗e T́ − 30 } {
T´
2

( )
 Ecuación canonica :
𝜏3*
d C D CO ( t )
+C DCO ( t )=K 12 F M Dchar ( t ) + K 13∗F V Total−out
2 D
( t ) + K 14 T D (t ) + K 15 P D H O ( t ) + K 16 P DCO ( t ) + K 17 P DCO ( t )
dt 2 2 2

Aplicandola transformadade Laplace:

D D
D K 12 MD K 13 VD K 14 D K 15 D K 16 K 17
C CO 2 (t )= F H O(t )+ ∗F Total−out ( t ) + T (t )+ P H O(t )+ P (t)+ P
τ 3 S+1 2
τ 3 S+ 1 τ 3 S +1 τ 3 S+1 τ 3 S+ 1 CO τ 3 S+1 CO
2

𝜏3 V lecho
F ´
V
Total−out
K 12 1
F ´V
Total−out
K 13 C H´ O 2

FV ´
Total−out
K 14 −K 30 ¿ C char∗ PCO∗P
´ H O EA EA
K 3

V lecho T́
3 2

(
∗ 2 ∗e T́ − 30 } {*C} rsub {char} * acute {{P
)
K 15 EA 3
T́
−K 30 ¿ C char∗P´CO∗e K 30 } {*C} rsub {char} * {acute {acute {{P} rsub {CO2}}
−
V lecho T́
K 16 EA 3

−K 30 ¿ C char∗P ´´H O∗e 2

T́

V lecho

K 17 −K 30 } {*C} rsub {char} * acute {acute {{P} rsub {{H} rsub {2} O}}} * {e} ^ {{{EA} rsub {3
¿
T́ V

d C CO ( t ) 1
4).
dt
= [ F M (t )−F VTota l ( t )∗C CO ( t )+ r rxn 1 ( t )+ 2∗r rxn 2 ( t ) −rrxn 3 ( t ) +rrxn 3 ” ( t ) ]
V lecho CO out
d C CO ( t ) 1
= ¿
dt V lecho

d CH (t ) 1
5.) 2
= [ F MH ( t )−FVTotal−out ( t )∗C H ( t ) +r rxn 1 ( t ) +r rxn 3 ( t ) +r rxn 3 ” ( t ) ]
dt V lecho 2 2

d CH (t )
2 1
= ¿
dt V lecho

Balance de energía:
Rxn 3
dT ( t )
dt
1
= lecho
C v mlecho [
−FVOut ( t )
( ∑
Componentes
Cip C i ∗( T ( t ) −T carbon −¿ ) − ∑ r i ∆ H i +Q out
) Rxn 1 ]
C lecho
v mlecho
dT ( t )
dt
[
= −F VOut ( t ) ( C char CO 2 H 2O
p ∗Cchar ( t )+ C p ∗C CO2 ( t ) +C p ∗C H 2O ( t )+ CCO H2
p ∗C CO ( t )+ C p ∗C H 2 (t) )∗( T ( t )−

d T D (t )
C lecho =Z 1 F V Total−out
D
v mlecho (t ) + z 2 C D char (t ) + Z3 C DCO ( t ) + Z 4 C DH 2 O ( t ) + Z 5 C D CO ( t )+ Z 6 C DH ( t ) + Z 7 T D ( t ) +Z 8
dt 2 2

Z1 =−¿)*(T́ −T carbon−¿) Z2 =−¿)*(T́ −T carbon−¿)

Z3 =−¿)*(T́ −T carbon−¿) Z 4=−¿ )*(T́ −T carbon−¿)

Z5 =−¿)*(T́ −T carbon−¿) Z6 =−¿)*(T́ −T carbon−¿)

Z7 =−¿)
EA 1 EA 2
EA 1 EA 2 E
−∆ H 1∗K 10 ¿ C char∗P ´H O∗ 2
( T́ 2
∗e T́
)
−∆ H 2∗K 20 ¿C char∗P CO2
´ ∗
( T´ 2
∗e T́
) −∆ 3∗K 30 ¿ C char∗P´CO∗P ´H O∗
2
( T
EA 1 EA 2
T́ T́
Z 8=−∆ H 1∗K 10 ¿ C char∗e Z 9=−∆ H 2∗K 20 ¿ C char∗e
EA 2
T́
Z10=−∆ H 3∗P ´H O∗K 30 ¿ Cchar∗e 2
 d T D (t ) ´ D
C lecho =−( ∑ C ip C i )∗( T́ −T carbon−¿ )F V Total−out ( t )−( F´VOut C char ´V CO
D
mlecho p )∗( T́ −T carbon−¿ ) C char ( t )− ( F Out C p )∗( T
2
v
dt
)
EA 1 EA 2
EA1 EA 2
+ ∆ H 1∗K 10 ¿ C char∗ P ´H O∗ 2
( T́ 2
∗e T́
) + ∆ H 2∗K 20 ¿C char ∗ PCO2
´ ∗
( T´ 2
∗e T́
) + ∆ H 3∗K 30 ¿ Cchar ∗ P´CO∗ P ´H O∗
2
(

Ecuaciones adicionales:

dP reactor
∗V lecho
dt
=F Mtota l −F Mtota l + K 1∗C char∗PH O ( t )∗∆ H 1 + K 2∗C char ∗P CO ( t )∗∆ H 2 + K 3∗C char ∗PH O ( t )∗PCO
RT ingas out gases 2 2 2

dP reactor
∗V lecho
dt
=a F f DF ∗a ¿ f DF ∗a P f DP ∗aP f DP ∗a P f DP
V V V
RT Out Out ¿ H 20 H 20 CO 2 CO2 CO CO

a F =1
V
Out

a ¿=−1
EA 1 EA 2
T́ T́
a P =−∆ H 1∗K 10 ¿ C char∗e
H 20
a P =−∆ H 2∗K 20 ¿ C char∗e
CO 2
 EA 2
T́
aP CO
=−∆ H 3∗P ´H O∗K 30 ¿ Cchar∗e
2

Símbolo Nombre Unidad

F Mi Flujo molar de entrada al lecho mol. s−1
de cada uno de los componentes
i
Ci Concentración del componente i
al interior del lecho fluidizado

F VOut Flujo volumétrico total de salida

V lecho Volumen del lecho
T Temperatura del lecho
fluidizado
C ip Capacidad calor’ıfica de cada
componente
T carbon−¿ Temperatura de entrada del
carbón

∆ Hi Calores de reacci’on

C lecho
v
Capacidad calor’ıfica del lecho
mlecho Masa del lecho
Preactor Presi’on del reactor
r rxni Representa la cin’etica de
reacci’on i