UNIVERSIDAD AUTONOMA DE SAN LUIS POTOSI
Facultad de Ciencias Químicas
Fenómenos de Transporte II
Cuarto Examen Parcial
Profesor:
Erik Cesar Herrera Hernández
Equipo No 5:
Alviso Camacho Alex
Carrizalez Castro Adán
Delgado Alejandre Aránzazu
Falcón González Luis Felipe
Gaitán Domínguez Pablo Cesar
Hernández Hernández Nohemí
Fecha: 3/12/2021
Ejemplo 18.3 Bird 2da edición: DIFUIÓN CON UNA REACCION HETEROGENEA
A continuación consideraremos un modelo sencillo para un reactor catalítico, como el que se muestra, donde
se lleva a cabo una reacción 2A + B. Un ejemplo de una reacción de este tipo sea la dimerización catalizada
por un sólido de CH3CH = CH2. Se supone que cada partícula catalizadora está rodeada de una película gaseosa
estancada a través de la cual A debe difundirse para alcanzar la superficie del catalizador, como se muestra
en la figura. En la superficie del catalizador se supone que la reacción 2A + B ocurre instantáneamente, y que
luego el producto B se difunde retirándose a través de la película gaseosa hacia la corriente turbulenta
principal compuesta por A y B. Se desea obtener una expresión para la velocidad local de conversión de A a B
cuando se conocen el espesor efectivo de la película de gas y las concentraciones XAO y XBo en la corriente
principal. Se supone que la película de gas es isotérmica, aunque en muchas reacciones cataliticas no puede
despreciarse el calor que se genera en h reacción.
Diagrama:
Suposiciones:
Se lleva a cabo una reacción
heterogénea instantánea
Estado Estacionario
La película de gas es isotérmica
Transferencia solo en una
dirección
Balance general
𝐸−𝑆±𝐺 =𝐴 → 𝐸−𝑆 = 0
𝑆𝑁𝐴𝑍 |𝑍 − 𝑆𝑁𝐴𝑍 |𝑍+∆𝑍 = 0
Dividiendo entre el volumen y tomando su límite cuando ΔZ→0
𝑁𝐴𝑍 |𝑍 − 𝑁𝐴𝑍 |𝑍+∆𝑍
lim
∆𝑧→0 ∆𝑧
Aplicando la definición de derivada
𝒅𝑵𝑨𝒁
− =𝟎
𝒅𝒛
Como la reacción es del tipo 2A→B por la estequiometria se tiene:
1
𝑁𝐵𝑧 = − 𝑁𝐴𝑧
2
Sustituyendo este valor de la estequiometria en la ecuación constitutiva:
𝜕𝜒𝐴
𝑁𝐴𝑧 = 𝜒𝐴 (𝑁𝐴𝑧 + 𝑁𝐵𝑧 ) − 𝑐𝐷𝐴𝐵
ð𝑧
𝑐𝐷𝐴𝐵 𝑑𝜒𝐴
𝑁𝐴𝑧 = −
1
1 − 2 𝜒𝐴 𝑑𝑧
Insertando la ecuación de NAZ en el balance de masa:
𝑑 1 𝑑𝜒𝐴
( )=0
𝑑𝑧 1 − 1 𝜒 𝑑𝑧
2 𝐴
Integrando dos veces para z:
1 𝑑𝜒𝐴
∫𝑑( ) = ∫ 𝑑𝑧
1
1 − 2 𝜒𝐴 𝑑𝑧
1 𝑑𝜒𝐴
= 𝐶1
1
1 − 2 𝜒𝐴 𝑑𝑧
𝑑𝜒𝐴
∫ = 𝐶1 ∫ 𝑑𝑧
1
1 − 2 𝜒𝐴
1
−2 ln (1 − 𝜒𝐴 ) = 𝐶1 𝑧 + 𝐶2
2
Otra forma de encontrar las constantes de integración es la siguiente:
−(2𝑙𝑛𝐾1 )𝑧 − (2𝑙𝑛𝐾2 ) = 𝐶1 𝑧 + 𝐶2
Condiciones Frontera
𝑪𝑭𝟏 𝒛 = 𝟎 𝝌𝑨 = 𝝌𝑨𝟎
1
−2 ln (1 − 𝜒𝐴0 ) = 𝐶1 (0) + 𝐶2
2
1
−2 ln (1 − 𝜒𝐴0 ) = 𝐶2
2
𝑪𝑭𝟐 𝒛 = 𝜹 𝝌𝑨 = 𝟎
1 1
−2 ln (1 − (0)) = 𝐶1 𝛿 + (−2 ln (1 − 𝜒𝐴0 ))
2 2
1
2 𝑙𝑛 (1 − 2 𝜒𝐴0 )
= 𝐶1
𝛿
Sustituimos C1 y C2 en la ecuación general:
1
1 2 𝑙𝑛 (1 − 2 𝑥𝐴0 ) 1
−2 ln (1 − 𝑥𝐴0 ) = 𝑧 − 2 𝑙𝑛 (1 − 𝑥𝐴0 )
2 𝛿 2
1 1
1 2 𝑙𝑛 (1 − 𝑥𝐴0 ) 2𝑙𝑛(1 − 2 𝑥𝐴0 )
ln (1 − 𝑥𝐴0 ) = 2 𝑧−
2 −2𝛿 −2
1
1 𝑙𝑛 (1 − 𝑥𝐴0 ) 1
ln (1 − 𝑥𝐴0 ) = − 2 𝑧 + 𝑙𝑛(1 − 𝑥𝐴0 )
2 𝛿 2
1 1 𝑧
ln (1 − 𝑥𝐴0 ) = 𝑙𝑛 (1 − 𝑥𝐴0 ) [1 − ]
2 2 𝛿
Utilizando las ecuaciones de frontera se llega al resultado final:
𝑧
[1− ]
1 1 𝛿
(1 − 𝑥𝐴0 ) = (1 − 𝑥𝐴0 )
2 2
De esta forma, la densidad de flujo a través de la membrana utilizando y cambiando valores en la ecuación constitutiva
queda como:
2𝑐𝐷𝐴𝐵 1
𝑁𝐴𝑍 = 𝑙𝑛 (1 − 𝑥𝐴0 )
𝛿 2
Discusión problema 1
Para el análisis de este sistema debe recurrirse a un modelo muy simplificado. A tal efecto se imagina que cada partícula
catalítica está rodeada de una película gaseosa estancada a través de la cual difunde el reactante A hasta alcanzar la
superficie catalítica. El producto B difunde en sentido contrario hasta alcanzar la corriente turbulenta de gas compuesto
por A y B.
La ecuación anterior [𝑁𝐴𝑧 ] representa la velocidad de dimerización por unidad de área de superficie catalítica. Debe
observase que, aunque la reacción es instantánea en la superficie del catalizador, la conversión de A en B se realiza a
velocidad “finita” ya que la difusión en la película es la que aporta el reactante en la superficie, luego, la reacción está
“controlada por la difusión”, o sea la “difusión controla el fenómeno”.
Otras limitaciones al modelo son las desigualdades en la distribución del flujo global, los efectos térmicos, la perdida de
actividad del catalizador, etcétera.
Problema 2
2.- Se tiene un sólido hecho de un material ligeramente soluble que se pone en contacto con una gran cantidad
de agua. La rapidez de disolución del sólido está controlada por la difusión del soluto en una película
estancada de fluido de espesor R1-R2. La concentración del soluto en la superficie del cilindro es CA1 y fuera
de la película es CA2. La solución es diluida y la disolución ocurre en estado estacionario.
Análisis del problema:
La geometría del sistema es cilíndrica, debido a que el sólido está disolviéndose, la concentración del material en R1 es
mayor que R2 por lo tanto CA1 > CA2.
La difusión se dará en dirección radial (desde la superficie del cilindro r=R1 hasta la superficie delimitada por r=R2).
Consideraciones:
Estado estacionario
Sin reacción química (no generación)
Propiedades constantes
a. Formule un balance de masa sobre el soluto A en un elemento diferencial del volumen
represntativo. ¿Es constante el área de difusión?
Comenzamos realizando un balance de masa sobre el soluto en el elemento diferencial de volumen, para nuestro
problema será el de un cilíndro 2𝜋𝑟∆𝑟𝐿.
El balance de masa general puede enunciarse como:
𝐹𝑙𝑢𝑗𝑜 𝑚𝑜𝑙𝑎𝑟 𝑑𝑒 𝐴 𝑞𝑢𝑒 𝐹𝑙𝑢𝑗𝑜 𝑚𝑜𝑙𝑎𝑟 𝑑𝑒 𝐴 𝑞𝑢𝑒 𝑅𝑎𝑝𝑖𝑑𝑒𝑧 𝑑𝑒 𝑔𝑒𝑛𝑒𝑟𝑎𝑐𝑖ó𝑛 𝑅𝑎𝑝𝑖𝑑𝑒𝑧 𝑑𝑒 𝑎𝑐𝑢𝑚𝑢𝑙𝑎𝑐𝑖ó𝑛
(𝑒𝑛𝑡𝑟𝑎 𝑝𝑜𝑟 𝑑𝑖𝑓𝑢𝑠𝑖ó𝑛 𝑒𝑛) − ( 𝑠𝑎𝑙𝑒 𝑝𝑜𝑟 𝑑𝑖𝑓𝑢𝑠𝑖ó𝑛 𝑒𝑛 ) + (𝑑𝑒 𝐴 𝑝𝑜𝑟 𝑟𝑒𝑎𝑐𝑐𝑖ó𝑛 𝑞𝑢í𝑚𝑖𝑐𝑎) = ( 𝑑𝑒 𝐴 𝑒𝑛 𝑒𝑙 𝑣𝑜𝑙𝑢𝑚𝑒𝑛 )
𝑙𝑎 𝑠𝑢𝑝𝑒𝑟𝑓𝑖𝑐𝑖𝑒 𝑟 𝑙𝑎 𝑠𝑢𝑝𝑒𝑟𝑓𝑖𝑐𝑖𝑒 𝑟 + ∆𝑟 𝑒𝑛 𝑒𝑙 𝑣𝑜𝑙𝑢𝑚𝑒𝑛 2𝜋𝑟∆𝑟𝐿 2𝜋𝑟∆𝑟𝐿
𝐸−𝑆± 𝐺 =𝐴
Como en el sistema no ocurre reacción química alguna y se encuentra en estado estacionario, los términos de generación
y de acumulación de soluto A dentro del elemento diferencial de volumen son cero, pordemos escribir el balance de
materia como:
𝑆𝑁𝐴𝑟|𝑟 − 𝑆𝑁𝐴𝑟|𝑟+∆𝑟 = 0
Donde S es el área de la sección transversal
𝑁𝐴𝑟|𝑟 (2𝜋𝑟𝐿) − 𝑁𝐴𝑟|𝑟+∆𝑟 (2𝜋(𝑟 + ∆𝑟)𝐿 = 0
(𝑁𝐴𝑟 ∙ 2𝜋𝑟𝐿)|𝑟 − (𝑁𝐴𝑟 ∙ 2𝜋𝑟𝐿)|𝑟+∆𝑟 = 0
Dividir por el elemento diferencial de volumen:
(𝑁𝐴𝑟 ∙ 2𝜋𝑟𝐿)|𝑟 − (𝑁𝐴𝑟 ∙ 2𝜋𝑟𝐿)|𝑟+∆𝑟
=0
2𝜋𝑟∆𝑟𝐿
1 (𝑁𝐴𝑟 ∙ 𝑟)|𝑟+∆𝑟 − (𝑁𝐴𝑟 ∙ 𝑟)|𝑟
− [ ]=0
𝑟 ∆𝑟
El término 𝑁𝐴𝑟 ∙ 𝑟 evaluado en 𝑟 + ∆𝑟 es diferente del que se está evaluando en 𝑟, ya que ahora el área a través de la cual
ocurre la difusión es variable, es una función de la radio.
Si tomamos el límite cuando ∆𝑟 → 0
1 (𝑁𝐴𝑟 ∙ 𝑟)|𝑟+∆𝑟 − (𝑁𝐴𝑟 ∙ 𝑟)|𝑟 1 𝑑(𝑟 ∙ 𝑁𝐴𝑟 )
lim − [ ]=− =0
∆𝑟→0 𝑟 ∆𝑟 𝑟 𝑑𝑟
𝑑(𝑟 ∙ 𝑁𝐴𝑟 ) 1
= 0 (− )
𝑑𝑟 𝑟
𝒅
𝒓𝑵 = 𝟎
𝒅𝒓 𝑨𝒓
b. Derive la ecuación diferencial del flux difusivo
Integrar la ecuación del balancde de masa:
𝑑
𝑟 ∙ 𝑁𝐴𝑟 = 0
𝑑𝑟
∫ 𝑑(𝑟 ∙ 𝑁𝐴𝑟 ) = ∫ 0𝑑𝑟
𝑪𝟏
(𝑟 ∙ 𝑁𝐴𝑟 ) = 𝐶1 → 𝑵𝑨𝒓 =
𝒓
Nos indica que en Nos indica que el flux difusivo 𝑁𝐴𝑟 es
coordenadas cilindricas el una función inversa de la variable
producto (𝑟 ∙ 𝑁𝐴𝑟 ) es cte independiente, y no una constante
c. Derive la ecuación diferencial cuya solución proporcionará el perfil de concentraciones
𝜕𝐶𝐴
Sustituyendo la ley de Fick en coordenadas cilíndricas para difusión radial 𝐽𝐴,𝑟 = −𝐷𝐴𝐵 𝜕𝑟
en la ecuación
diferencial del flux difusivo tenemos:
𝑑𝐶𝐴 𝐶1
𝑁𝐴𝑟 = −𝐷𝐴𝐵 =
𝑑𝑟 𝑟
Integrando:
𝐶𝐴2 𝑟2
𝑑𝑟
−𝐷𝐴𝐵 ∫ 𝑑𝐶𝐴 = 𝐶1 ∫
𝐶𝐴1 𝑟1 𝑟
Solución general:
𝑪𝑨 = 𝑪𝟏 ′ 𝑳𝒏 𝒓 + 𝑪𝟐 ′
En donde:
𝐶1 𝐶2
𝐶1 ′ = ; 𝐶2 ′ =
−𝐷𝐴𝐵 −𝐷𝐴𝐵
d. Escribe y explique las condiciones de frontera:
𝑟 = 𝑅1 ; 𝐶𝐴 = 𝐶𝐴1
𝑟 = 𝑅2 ; 𝐶𝐴 = 𝐶𝐴2
Tenemos en nuestro sistema un sólido poco soluble lo cual la concentración del soluto en la
superficie será menor que la concentración fuera de la pelicula, es decir, en R1 la
concentración será mayor que en el espesor R2.
e. Derive la expresion del perfil de concentraciones.
Sustituyendo las condiciones frontera
CFI r=R1 CA=CA1 CFII r=R2 CA=CA2
𝐶𝐴1 = 𝐶′1 𝑙𝑛𝑅1 + 𝐶′2 𝐶𝐴2 = 𝐶′1 𝑙𝑛𝑅2 + 𝐶′2
*Despejar 𝐶′2 *Despejar 𝐶′2
𝐶′2 = 𝐶𝐴1 − 𝐶′1 𝑙𝑛𝑅1 𝐶′2 = 𝐶𝐴2 − 𝐶 ′1 𝑙𝑛𝑅2
Igualar las expresiones obtenidas
𝐶𝐴1 − 𝐶 ′1 𝑙𝑛𝑅1 = 𝐶𝐴2 − 𝐶 ′1 𝑙𝑛𝑅2
𝐶𝐴1 − 𝐶𝐴2 = 𝐶 ′1 𝑙𝑛𝑅1 − 𝐶 ′1 𝑙𝑛𝑅2
𝐶𝐴1 − 𝐶𝐴2 = 𝐶 ′1 (𝑙𝑛𝑅1 − 𝑙𝑛𝑅2)
𝐶𝐴1 − 𝐶𝐴2
𝐶 ′1 =
𝑅1
(𝑙𝑛 𝑅2)
Sustituir 𝐶 ′1 en 𝐶′2
𝐶′2 = 𝐶𝐴1 − 𝐶′1 𝑙𝑛𝑅1
𝐶𝐴1 − 𝐶𝐴2
𝐶′2 = 𝐶𝐴1 − 𝑙𝑛𝑅1
𝑅1
(𝑙𝑛 𝑅2)
Sustituir 𝐶 ′1 y 𝐶′2 en la solución general
𝐶𝐴1 − 𝐶𝐴2 𝐶𝐴1 − 𝐶𝐴2
𝐶𝐴 = 𝑙𝑛𝑟 + 𝐶𝐴1 − 𝑙𝑛𝑅1
𝑅1 𝑅1
(𝑙𝑛 𝑅2) (𝑙𝑛 𝑅2)
Reacomodando
𝒓
𝒍𝒏
𝑪𝑨 = 𝑪𝑨𝟏 − (𝑪𝑨𝟏 − 𝑪𝑨𝟐 ) 𝑹𝟏
𝑹𝟐
𝒍𝒏 𝑹𝟏
f. Derive la expresión del flux difusivo.
Para poder obtener la expresión final del flux difusivo debemos comenzar derivando la expresión del perfil de
concentraciones:
𝑑𝐶𝐴 𝑑 𝐿𝑛(𝑟/𝑅1 )
= [𝐶𝐴1 − (𝐶𝐴1 − 𝐶𝐴2 ) ]
𝑑𝑟 𝑑𝑟 𝐿𝑛(𝑅2 /𝑅1 )
𝐿𝑛(𝑟/𝑅1 )
𝑑𝐶𝐴1 𝑑(𝐶𝐴1 − 𝐶𝐴2 ) (𝐿𝑛(𝑅2 /𝑅1 ))
−
𝑑𝑟 𝑑𝑟
(𝐶𝐴1 − 𝐶𝐴2 ) 𝑑(𝐿𝑛(𝑟/𝑅1 ))
− ∙
𝐿𝑛(𝑅2 /𝑅1 ) 𝑑𝑟
−𝐶𝐴1 − 𝐶𝐴2 1
∙
𝐿𝑛(𝑅2 /𝑅1 ) 𝑟
𝑑𝐶𝐴 −𝐶𝐴1 − 𝐶𝐴2 1
= ∙
𝑑𝑟 𝐿𝑛(𝑅2 /𝑅1 ) 𝑟
Sustiruir la expresión anterior en la ecuación del fluix difusivo
𝑑𝐶𝐴 −𝐶𝐴1 − 𝐶𝐴2 1
𝑁𝐴𝑟 = −𝐷𝐴𝐵 → −𝐷𝐴𝐵 ( ∙ )
𝑑𝑟 𝑅
𝐿𝑛 (𝑅2 ) 𝑟
1
𝑪𝑨𝟏 − 𝑪𝑨𝟐 𝟏
𝑵𝑨𝒓 = 𝑫𝑨𝑩 ∙
𝑹
𝑳𝒏 (𝑹𝟐 ) 𝒓
𝟏
Discusión problema 2
La transferencia de materia de un soluto de una fase a otra incluye el transporte desde el seno de una fase a la
superficie limite o interface, pero esto también se puede desde la interface al seno del segundo fluido, esta teoría se
basa en la idea de que en la zona de contacto entre las dos fases se forma una película de fluido o capa límite de
transferencia de materia (Doran, 1998).
Se formuló un balance de materia para un sistema de forma cilíndrica en estado estacionario en unidades de
concentración, no tenemos reacción química, por tanto, sabemos que la generación y acumulación en este caso serán
nulas, se llevó a cabo el estudio sobre un incremento o una sección del sistema ∆𝑟 estableciendo así que la cantidad que
entra de A en el plano r será igual a la cantidad de A que sale en el plano 𝑟 + ∆𝑟.
Para nuestro sistema sabemos que es una solución diluida por lo que nuestra ecuación del flux másico se simplificara ya
que ‘’Xa’’ tiende a 0.
La rapidez de la disolución del solido se controla por la difusión del soluto en una película estancada y este se disolverá
al agregarlo en una solución diluida asumiendo que se presenta en la fase acuosa poca cantidad del sólido.
Para R1-R2 se contiene una especie A en un disolvente B, podemos ver conforme nos desplazamos en el radio se
presentará una difusión que se dará poco a poco o lentamente debido a que es ligeramente soluble entonces en una
distancia R1 tenemos una concentración “CA1” que en este punto se presenta una mayor concentración del sólido que en
una distancia R2 donde tenemos una concentración “CA2”.
Por último, hay que recordar que la difusión es el movimiento de la materia de una región de alta concentración a una
región de baja concentración, uno de los factores que afecta la difusión es el tamaño y forma del soluto, pero también
hay otros factores como la viscosidad del disolvente y naturaleza de la partícula.
Bibliografía
1.- Referencia: R. Byron Bird Warren E. Stewart Edwin N. Lightfoot. (2006). Fenómenos de transporte. México: Limusa
Wiley.
2.- Referencia: Ing. Jorge E. Robín Ing. Marcela N. Kaminsky Juan M. Dominguez. (2008). ANÁLISIS ENVOLVENTE PARA
DIFUSIÓN. 30/11/2021, de UNIVERSIDAD TECNOLOGICA NACIONAL FACULTAD REGIONAL ROSARIO DEPARTAMENTO DE
INGENIERIA QUIMICA Sitio web:
https://www.frro.utn.edu.ar/repositorio/catedras/quimica/3_anio/fenomenos_transporte/Unidad_tematica_8.pdf
3.- Referencia: Diego Dolores. (2015). Difusión a través de una película de gas. 29/11/2021, de Academia.edu Sitio web:
https://www.academia.edu/39732137/Difusi%C3%B3n_a_trav%C3%A9s_de_una_pel%C3%ADcula_de_gas#:~:text=Difu
si%C3%B3n%20a%20trav%C3%A9s%20de%20una%20pel%C3%ADcula%20de%20gas%20estancada%20La,transferencia
%20de%20materia%20entre%20fases.&text=Esta%20teor%C3%ADa%20se%20basa%20en,materia%20(Doran%2C%201
998).
