2.5 Masa

2.5. Difusión y reacción en un tubo.
Considera una mezcla gaseosa de A diluido en B que

está dentro de un tubo de 10 cm de largo. Dentro del tubo ocurre la reacción A → B a una
tasa igual a Ra=kCa En el extremo izquierdo del tubo la concentración de “A” se mantiene
en 1x10-5 mol/cm3 y en el extremo derecho en 0.3x10-5 mol/cm3. El coeficiente de difusión
de “A” en “B” es 0.7 cm/s.
a) Deriva las expresiones que describan el perfil de concentraciones y el flux de “A”
dentro del tubo.
b) Explica el significado del grupo

√
k L2 .
D
c) Calcula el perfil de concentraciones y el flux de A cuando k= 0.1 s-1; k= 0.01 s-1; k=
0.001 s-1; y k= 0 s-1.
d) Grafica en papel milimétrico los perfiles de concentraciones y los fluxes calculados
en el inciso anterior usando como parámetro
conclusiones.
√ k L2 . Analiza las gráficas y extrae
D
Ca=Cao= 0.3x10-5 mol/cm3 Ca=Cal=0.3x10-5 mol/cm3
z=0 z=L=10cm
Suposiciones:
 Estado estacionario
 Proceso reactivo
 P,T cte.
 Sistema diluido, predomina la difusión (D)
 Coeficiente de difusividad efectivo cte.
 Flujo unidireccional axial (en z)
 As constante
Balance de materia y energía:
El sistema presenta entrada (E) y salida (S) de flux (Ja) además de presentar una reacción
dentro del sistema, por lo que hay termino de consumo (C) del reactivo A dada por la
velocidad de reacción (Ra).
E−S−C=0
El sistema está en coordenadas cilíndricas, considerando esto, se analiza el sistema a una

escala microscópica y posteriormente se lleva a escala macroscópica.
Primero se toma un elemento diferencial de volumen de la misma geometría del sistema:
Δr
Δz
2
∆ V =2 π ∆ r ∆ z= As ∆ z
Estableciendo así en la ecuación de balance de materia para el flux difusivo del sistema.
As Ja∨¿ z− AsJa∨¿z +∆ z−RaAs ∆ z=0¿ ¿
Posteriormente se divide la expresión entre el volumen diferencial y se le aplica el límite

cuando Δz tiende a cero puesto que es un valor muy diminuto.
Ja∨¿ z +∆ z
lim Ja∨¿ z− −Ra=0 ¿ ¿
∆ z→o ∆z
Esta expresión puede modificarse aplicando la definición de la derivada.
−dJa
+ Ra=0 … Ec. (1)
dz
De acuerdo con la ley de Fick:

dCa
Ja=−D … Ec. (2)
dz
Sustituyendo (2) en (1):
D
d ( dCa
dz )
−Ra=0
dz
2
d Ca
D 2
−Ra=0 … Ec. (3)
dz
Como se indica en el enunciado:
Ra=krCa … Ec. (4)
Finalmente sustituyendo (4) en (3) se obtuvo la siguiente ecuación diferencial homogénea

de segundo orden:
2
d Ca
D 2
−krCa=0 … Ec. (5)
dz
Determinando las condiciones límites:
C.L.:
1. Ca=Cao z=0
2. Ca=Ca z=L
Por fines prácticos se puede modificar (5) en términos adimensionales multiplicando por:
2
L
D
2 2
2 d Ca L krCa
L 2
− =0… Ec. (6)
dz D
Considerando y derivando llegamos a la siguiente ecuación diferencial:

2
z dz 2 kL
z∗¿ dz∗¿ β=
L L D
2
d Ca
2 2 … Ec. (7)
d z∗¿ −β Ca=0 ¿
Una vez hechos estos cambios, ya se puede resolver (6), se utilizó el método de coeficientes
constantes:
m z∗¿=0¿
¿
m2 em z∗¿−β e
Ec. Auxiliar:
2 2
m −β =0
2 2
m =β
m1=β m2 =−β
Resolviendo (6) tenemos:

−βz∗¿ ¿
¿
Ca=K 1 e βz∗¿+ K e 2
… Ec. (8)
Considerando las siguientes identidades hiperbólicas podemos reescribir (6):

1. cosh ⁡¿
2. sen h ⁡¿
Ca=C 1 cosh ¿ ¿… Ec. (9)
Estableciendo las nuevas condiciones límite, se despejan las constantes de integración:

C.L.:
1. z∗¿ 0 Ca=Cao
2. z∗¿ 1 Ca=Cal
Cao=C 1 … Ec. (10)
Cal=C 1 cosh ( β ) +C 2 sen h(β ) … Ec. (11)
Sustituyendo (10) en (11):

Cal=Cao cosh ( β ) +C 2 senh ( β )
Cal−Cao cosh ( β )
C 2= … Ec. (12)
senh ( β )
Una vez conocidas las constantes de integración, se sustituye (10) y (12) en (9):
Ca=Cao cosh ( β )+
( Cal−Caocosh ( β )
senh ( β )
senh (β)
)
Ca=Cao cosh ¿ ¿… Ec. (13)
Se procede a sustituir los valores de β y z* en (13) respectivamente:
( (√ )
)
k L2
Cal−Cao cosh
(√ ) D
√
2 2
kL kL
Ca=Cao cosh + senh( )
(√ )
D 2 D
kL
senh
D
(√ )
( (√ ) (√
) )
2 2 2
kL z Cal kL kL z
Ca=Cao cosh ( )+ −Cao cotanh senh ( )
(√ )
D L 2 D D L … Ec.
kL
senh
D
(14)
Para fines prácticos, se le denominará “γ” a todo el término en rojo:
(√ ) (√ )
2 2
kL z k L z … Ec. (15)
C a=Cao cosh ( ) + γ senh ( )
D L D L
Reescrita en números adimensionales:

Ca=Cao cosh ¿ ¿… Ec. (16)
Una vez determinado el perfil de concentración, podemos determinar el flux difusivo

derivando (14) respecto a z y sustituyendo en la Ley de Fick:
dC a= (√ ( )) (√ ) (√ ( )) (√ )
k L2 1
D L
γ cosh
k L2 z
D L
( )−
k L2 1
D L
Cao senh
k L2 z
D L
( ) dz
Ja=−D
((√ ( )) (√ ( )) ( √ ( )) (√ ( )))
k L2 1
D L
γ cosh
k L2 z
D L
−
k L2 1
D L
Cao senh
k L2 z
D L
Ja=−D ( √ ( ))( (√ ( ))
k L2 1
D L
γ cosh (√ ( )))
k L2 z
D L
−Cao senh
k L2 z
D L
… Ec. (17)
Reescrita en números adimensionales:
Ja=−D ( βL ) ¿… Ec. (18)

A partir de los modelos matemáticos determinados laboriosamente (16 y 18) y los datos
proporcionados por el problema, se puede calcular, en Excel, valores para el perfil de
concentración y el flux difusivo a distintas constantes velocidad de reacción en distintas
longitudes del sistema; posteriormente estos valores serán graficados.
Tabla 1. Datos proporcionados
Cao Cal
L (cm) D (cm2/s) k (s-1)
(mol/cm3) (mol/cm3)
1.00E-05 3.00E-06 10 0.7 0 0.001 0.01 0.1
k (s-1) 0 k (s-1) 0.001 k (s-1) 0.01 k (s-1) 0.1

z (cm) β 0 β 0.377964473 β 1.195228609 β 3.77964473
Ca Ja Ca Ja Ca Ja Ca Ja
(mol/cm3) (mol/cm2s) (mol/cm ) (mol/cm2s)
3
(mol/cm ) (mol/cm2s)
3
(mol/cm ) (mol/cm2s)
3
0 #¡DIV/0! #¡DIV/0! 1.00E-05 5.28E-07 1.00E-05 8.38E-07 1.00E-05 2.61E-06

1 #¡DIV/0! #¡DIV/0! 9.25E-06 5.18E-07 8.87E-06 7.44E-07 6.90E-06 1.78E-06
2 #¡DIV/0! #¡DIV/0! 8.52E-06 5.09E-07 7.87E-06 6.60E-07 4.80E-06 1.20E-06
3 #¡DIV/0! #¡DIV/0! 7.80E-06 5.01E-07 6.98E-06 5.86E-07 3.39E-06 7.94E-07
4 #¡DIV/0! #¡DIV/0! 7.09E-06 4.94E-07 6.19E-06 5.20E-07 2.48E-06 5.04E-07
5 #¡DIV/0! #¡DIV/0! 6.39E-06 4.87E-07 5.49E-06 4.62E-07 1.92E-06 2.86E-07
6 #¡DIV/0! #¡DIV/0! 5.69E-06 4.81E-07 4.87E-06 4.10E-07 1.64E-06 1.10E-07
7 #¡DIV/0! #¡DIV/0! 5.01E-06 4.76E-07 4.32E-06 3.64E-07 1.60E-06 -4.96E-08
8 #¡DIV/0! #¡DIV/0! 4.33E-06 4.71E-07 3.82E-06 3.24E-07 1.79E-06 -2.17E-07
9 #¡DIV/0! #¡DIV/0! 3.66E-06 4.67E-07 3.39E-06 2.88E-07 2.23E-06 -4.15E-07
10 #¡DIV/0! #¡DIV/0! 3.00E-06 4.64E-07 3.00E-06 2.56E-07 3.00E-06 -6.74E-07
γ #¡DIV/0! -2.00E-05 -1.00E-05 -9.87E-06
Tabla 2. Memoria de cálculo
Perfil de concentración (Ca)-Longitud (z)
1.20E-05
1.00E-05
8.00E-06
Ca (mol/cm3)
k=0
6.00E-06 k=0.001
k=0.01
4.00E-06 k=0.1
2.00E-06
0.00E+00
0 2 4 6 8 10 12
z (cm)
Fig. 1: Gráfica de perfil de concentración contra longitud
Flux difusivo (Ja)-Longitud (z)

0.000003
0.0000025
0.000002
Ja (mol/cm2s)
k=0
0.0000015
k=0.001
0.000001
k=0.01
0.0000005 k=0.1
0
0 2 4 6 8 10 12
-0.0000005
-0.000001
z (cm)
Fig. 2: Flux difusivo contra longitud
Como se puede apreciar tanto en la memoria de cálculo como en el gráfico, cuando la

constante velocidad de reacción es cero, quiere decir que no hay reacción química dentro
del sistema por lo que los valores del perfil de concentración y flux difusivo son
indeterminados puesto que los modelos determinados (16 y 18) solo son útiles cuando si
hay reacción en el sistema. Adicionalmente, se observa que el perfil de concentración
disminuye a medida que aumenta la constante de velocidad para posteriormente aumentar,
en contraste con el flux difusivo decrece más rápido conforme aumenta la constante de
velocidad hasta llegar a valores negativos, este comportamiento sucede porque “A” se
consume más rápido sin embargo nunca se consume, después de cierto punto se vuelve a
formar el reactivo, esta hipótesis tiene fundamento en que entre más pequeño sea el valor
de la constante de velocidad el flux difusivo se comporta de manera más uniforme porque
“A” reacciona de manera más lenta; a raíz de estas observaciones, se puede concluir que la
reacción dentro del sistema es reversible. A lo que se concluye que, a mayor velocidad de
reacción, el transporte de materia del reactivo “A” será más lento.
Finalmente el término adimensional β, también conocido como número de Hatta, es un
parámetro que relaciona la velocidad de reacción con la difusión a través de una película
líquida, también representa la resistencia que presenta la difusión y la velocidad de reacción
para el transporte de materia en un sistema reactivo homogéneo, sí la constante de
velocidad disminuye mientras que la difusividad aumenta el transporte de materia se llevara
a cabo de manera mas sencilla que por el contrario si aumenta la constante de velocidad y
disminuye la velocidad.

2.5 Masa

Cargado por

Copyright:

Formatos disponibles

2.5 Masa

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

2.5 Masa

Cargado por

Copyright:

Formatos disponibles

2.5. Difusión y reacción en un tubo.

Considera una mezcla gaseosa de A diluido en B que

b) Explica el significado del grupo

en el inciso anterior usando como parámetro

Ca=Cao= 0.3x10-5 mol/cm3 Ca=Cal=0.3x10-5 mol/cm3

El sistema está en coordenadas cilíndricas, considerando esto, se analiza el sistema a una

Posteriormente se divide la expresión entre el volumen diferencial y se le aplica el límite

De acuerdo con la ley de Fick:

Finalmente sustituyendo (4) en (3) se obtuvo la siguiente ecuación diferencial homogénea

Considerando y derivando llegamos a la siguiente ecuación diferencial:

Resolviendo (6) tenemos:

Considerando las siguientes identidades hiperbólicas podemos reescribir (6):

Estableciendo las nuevas condiciones límite, se despejan las constantes de integración:

Cal=C 1 cosh ( β ) +C 2 sen h(β ) … Ec. (11)

Sustituyendo (10) en (11):

Se procede a sustituir los valores de β y z* en (13) respectivamente:

Reescrita en números adimensionales:

Una vez determinado el perfil de concentración, podemos determinar el flux difusivo

Reescrita en números adimensionales:

Ja=−D ( βL ) ¿… Ec. (18)

k (s-1) 0 k (s-1) 0.001 k (s-1) 0.01 k (s-1) 0.1

0 #¡DIV/0! #¡DIV/0! 1.00E-05 5.28E-07 1.00E-05 8.38E-07 1.00E-05 2.61E-06

Fig. 1: Gráfica de perfil de concentración contra longitud

Flux difusivo (Ja)-Longitud (z)

Fig. 2: Flux difusivo contra longitud

Como se puede apreciar tanto en la memoria de cálculo como en el gráfico, cuando la

También podría gustarte