INSTITUTO POLITÉCNICO NACIONAL
ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E INDUSTRIAS
EXTRACTIVAS
Departamento de Ingeniería Química Petrolera
Ingeniería de Reactores Heterogéneos
Proyecto de Segundo Periodo
SOLUCIÓN NUMÉRICA PARA LA REACCIÓN &
DIFUSIÓN INTRAPELLET
Integrantes:
- López Vázquez Bianca
- Sosa Perea Ana Gabriela
Grupo: 3PM61
Profesor: Edgar Ramírez Jiménez
31 Marzo 2022
- Ecuaciones del modelo
El balance de materia para una partícula catalizadora en la que la difusión
de masa y la reacción tienen lugar en ella conduce a la siguiente forma
dimensionada:
𝑑 2 𝑥 2 2 𝑑𝑥
+ ( ) − 𝜙12 𝑥 = 0 … … . . (1)
𝑑𝜉 2 𝜉 𝑑𝜉
Condiciones de frontera:
𝑑𝑥
𝜉=0 =0
𝑑𝜉
𝜉=1 𝑥=1
La ecuación (1) puede resolverse mediante un cambio de variable, es decir,
𝑦 = 𝑥𝜉, por lo que la ec. (1) se reduce a:
𝑑2𝑦
− 𝜙12 𝑦 = 0 … … … … … (2)
𝑑𝜉 2
En consecuencia, la ec. (2) se evalúa ahora dentro de las condiciones de
contorno anteriores, lo que conduce al perfil de concentración
adimensional dado en la ec. (3)
1 sinh(𝜙1 𝜉)
𝑥= … … … … . (3)
𝜉 𝑠𝑖𝑛ℎ𝜙1
- Solución Numérica
En esta tarea se le pide que resuelva numéricamente la ec. (2)
discreteándola mediante un enfoque de diferencias finitas, lo que dará
lugar a un sistema de ecuaciones algebraicas lineales con estructura
tridiagonal. Desarrolle un programa utilizando el lenguaje de programación
Fortran para resolver dicho sistema de ecuaciones (el algoritmo de Thomas
es el método recomendado para estos casos) y,
a) Determinar el perfil de concentración adimensional aproximado para
diferentes valores del módulo de Thiele (𝜙1 ), así como para diferente número
de puntos.
b) Compara la solución discreteada con respecto a la solución analítica
para cada uno de los casos analizados.
c) Discute tus resultados.
- Solución
Para resolver el problema, hemos decidido adoptar la siguiente estrategia
de solución.
Para poder encontrar un perfil de concentración se comienza con la
discretización de la ecuación para el balance de materiales de la partícula
catalítica:
𝑑2𝑦
− 𝜙12 𝑦 = 0
𝑑𝜉 2
𝑦𝑖+1 − 2𝑦𝑖 + 𝑦𝑖−1
− 𝜙12 𝑖 𝑦𝑖 = 0
∆𝜉 2
Haciendo un cambio de variable
𝑦 = 𝑥𝜉
𝑥𝑖+1 𝜉𝑖+1 − 2𝑥𝑖 𝜉𝑖 + 𝑥𝑖−1 𝜉𝑖−1
− 𝜙12 𝑖 𝑥𝑖 𝜉𝑖 = 0
∆𝜉 2
Esta discretización nos daría como resultado una expresión general para
encontrar nuestro perfil de concentraciones, siendo:
𝑥𝑖+1 𝜉𝑖+1 − 𝑥𝑖 𝜉𝑖 (2 + 𝜙12 ∆𝜉 2 ) + 𝑥𝑖−1 𝜉𝑖−1 = 0
Evaluando para cada valor de 𝜉 diferentes ecuaciones con diferentes
variables(x), por por ejemplo, para 𝜉2 tenemos:
0.8𝑥3 − 𝑥2 (0.9)(2 + 𝜙12 × 0.12 ) + 1 = 0
Es importante señalar que se decidió tomar un intervalo de 10 puntos de 0.1
pasos entre ellos, además de considerar el 𝜙𝑖2 con 5 valores (0.1,0.5,1,10,100),
con el objetivo de analizar el comportamiento de la eficiencia entre rangos
grandes y pequeños de alguna manera.
También es importante saber que tenemos información en dos puntos,
𝑑𝑥
cuando ξ es 1 y cuando ξ es 0, el primer punto x=1 y el último 𝑑𝜉 =0, este último
punto es importante resolver el problema y que al realizar la discretización
de este punto con un salto hacia atrás, se puede acceder a:
𝑥𝑖 − 𝑥𝑖−1
=0
∆𝜉
𝑥𝑖 − 𝑥𝑖−1 = 0
Una vez encontradas las ecuaciones, podemos construir una matriz y un
vector de resultados de la siguiente manera:
𝐴 𝐴 0 0 0 0 0 0 0 −1
𝐴 𝐴 𝐴 0 0 0 0 0 0 0
0 𝐴 𝐴 𝐴 0 0 0 0 0 0
0 0 𝐴 𝐴 𝐴 0 0 0 0 0
0 0 0 𝐴 𝐴 𝐴 0 0 0 0
0 0 0 0 𝐴 𝐴 𝐴 0 0 0
0 0 0 0 0 𝐴 𝐴 𝐴 0 0
0 0 0 0 0 0 𝐴 𝐴 𝐴 0
0 0 0 0 0 0 0 𝐴 𝐴 0
(𝑑𝑥)∆𝜉
- Resultados
En la gráfica 1 se muestran los valores de x con respecto a ξ, a partir de la
solución numérica, se muestran los perfiles de reconversión según la
coordenada radial sin dimensiones
Cuando el módulo de Thiele es de concentración de flujo dentro del esporo,
será más alto que el área de la superficie y, por lo tanto, la resistencia al
transporte interno puede despreciarse en por el contrario, cuando el módulo
de Thiel es alto, los gradientes de concentración no se pueden descuidar, lo
que indica que la resistencia al transporte interno es muy importante.
Grafica 1. Solución Discreteada
De acuerdo con la gráfica 2, se observan los valores de x con respecto a ξ,
los perfiles de conversión los perfiles de conversión de este acto se muestran
en función de la coordenada radial sin dimensión del catalizador, para
diferentes valores del módulo de Thiele, como se observa, cuando el módulo
de Thiele es de muy baja concentración en el interior del poro, es
prácticamente igual a la superficie, por lo tanto, puede ser ignorada la
existencia de un transporte interno.
Grafica 2. Solución Analítica
Por el contrario, cuando el módulo de Thiele es alto, los gradientes de
concentración son importantes y no se pueden descuidar lo que indica que
la resistencia al transporte interno es significativa.
Al comparar ambos gráficos se puede analizar que los valores calculados
en analítica, con respecto a la numérica, existe cierta correlación, que será
mayor cuando se hay un más intervalos.
- Conclusión
Al realizar los cálculos y obtener los gráficos, se puede concluir que los
gradientes de concentración y la resistencia al transporte interno son
significantes y no se pueden descuidar, lo que ocurre cuando el módulo de
Thiele es alto.
Se puede pensar que hay un comportamiento constante cuando el módulo
de Thiele es bajo, mientras que, por el contrario, si tenemos un valor de
Tielémodulo grande, el comportamiento no será constante.
Este módulo demuestra que existe una relación entre la velocidad de
reacción y la difusión en las esporas.
Hay un mayor factor de eficacia y para los valores elevados.
- Referencia bibliografica
Smith, J. M. (1991). Ingeniería de la Cinética Química (Sexta ed.).
CECSA.
[Link]
/[Link]
Hill, C. (2014). Introduction to Chemical Engineering Kinetics and
Reactor Design. (Segunda ed.) WILEY.
[Link]
gineering_Kinetics_and_Reactor_Design_2nd_ed_2014_
- Anexos
Código utilizado para resolver el problema se muestra a continuación:
Program Thiele
implicit none
real*8 x,e,deltae,fi,ex,til,xi,xf,a,b,c,d,m,at
integer n,i,j,k
dimension x(1000),e(1000),m(1000,1000),a(1000),b(1000),c(1000)
dimension d(1000),tili(1000),at(1000,1000)
open(unit=20,file='[Link]',action='write')
write(*,*)'Evaluar x cuando e=1' read(*,*)xi
write(*,*)' Evaluar x cuando e=0' read(*,*)xf
write(*,*)'Ingresar valor de intervalos' read(*,*)n
til(1)=0.1
til(2)=0.5
til(3)=1
til(4)=10
til(5)=100
e(1)=1
e(n+1)=0
deltae=(e(1)-e(n+1))/n
do i=2,n
e(i)=e(i-1)-deltae
end do
do i=1,n+1
write(*,*)e(i)
end do
do i=3,n+2 at(i,1)=e(i-1)
end do
at(2,1)=1 at(n+2,1)=0
do k=2,6
do i=1,n-1
do j=1,n
if([Link].j)then
m(i,j)=-e(i+1)*(2+til(k-1)*(deltae**2))
else
if([Link].i-1)then
m(i,j)=e(j+1)
else
if([Link].i+1)then m(i,j)=e(j+1)
else
m(i,j)=0
endif
endif
endif
end do
end do
m(n,n)=1
m(n,n-1)=-1
d(1)=-e(1)*xi
d(n)=xf*deltae
do i=2,n-1 d(i)=0
end do
do i=1,n b(i)=m(i,i)
end do
a(1)=0
do i=2,n
a(i)=m(i,i-1)
end do
c(n)=0
do i=1,n-1
c(i)=m(i,i+1)
end do
call TH (a,b,c,d,x,n)
do i=3,n+2
at(i,k)=x(i-2)
end do
end do
do i=2,6
at(2,i)=1
end do
do i=2,6
at(1,i)=til(i-1)
end do
do i=1,n+2
write(20,*)(at(i,k),k=1,6)
end do
pause
end
subroutine TH(a,b,c,d,x,n)
implicit none
real*8 a,b,c,d,x integer n,i
dimension x(1000),d(1000),a(1000),b(1000),c(1000)
do i=1,n-1
if(b(i).eq.0)then
write(*,*)'No es la solución'
goto 6
else
b(i+1)=b(i+1)-a(i+1)*c(i)/b(i)
d(i+1)=d(i+1)-a(i+1)*d(i)/b(i)
end if
end do
if(b(n).eq.0) then
write(*,*)'No es la solución'
goto 6
else x(n)=d(n)/b(n)
end if
i=n
10 i=i-1
if([Link].1)goto 6
x(i)=(d(i)-c(i)*x(i+1))/b(i)
goto 10
6 continue
return
end
