Tarea 7

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 17

UNIVERSIDAD VERACRUZANA

Nombre de la Licenciatura
INGENIERÍA QUÍMICA

Nombre del alumno


CARLOS ALBERTO MARTINEZ HERNANDEZ

Matrícula
S18015627

Número de la Tarea
TAREA 7

Unidad 9
DIFUSIÓN EN ESTADO ESTACIONARIO

Nombre del Maestro


LUCIO ÁLVAREZ PÉREZ

Fecha
03/12/2021
Fundamentos de Transferencia de Calor y
Materia:
Masa

Tarea: Tarea 7

CONVECCIÓN NATURAL Y FORZADA DE MASA EN FLUJO


LAMINAR Y TURBULENTO
Hasta ahora, se ha considerado la difusión de masa, la cual es la transferencia de masa debida
a un gradiente de concentración. A continuación, se considerará la convección de masa (o
transferencia de masa por convección), la cual es la transferencia de masa entre una superficie
y un fluido en movimiento, debido tanto a la difusión de masa como al movimiento de toda la
masa de fluido. Se mencionó con anterioridad que el movimiento del fluido mejora de manera
considerable la transferencia de calor, al quitar el fluido calentado que está cerca de la superficie
y reemplazarlo por el fluido más frío que se encuentra más alejado. De modo semejante, el
movimiento del fluido mejora de manera considerable la transferencia de masa, al quitar el fluido
con alta concentración que está cerca de la superficie y reemplazarlo por el fluido con
concentración más baja que está más alejado. En el caso límite de que no se tenga movimiento
de la masa de fluido, la convección de masa se reduce a la difusión de la misma, precisamente
como la convección se reduce a conducción. La analogía entre la convección de calor y la de
masa se cumple tanto para convección forzada como para la natural, para el flujo laminar como
para el turbulento, así como para los flujos interno y externo. Como la convección de calor, la de
masa también es compleja debido a las complicaciones asociadas con el flujo de fluidos, como
la configuración geométrica de la superficie, el régimen, la velocidad de flujo y la variación de las
propiedades del fluido y de la composición. Por lo tanto, hay que apoyarse en relaciones
experimentales para determinar la transferencia de masa. También, la convección de masa suele
analizarse en términos de masa, en lugar de sobre una base molar. Por esa razón, se presentarán
las formulaciones en términos de concentración de masa (densidad r o fracción de masa, w), en
lugar de concentración molar (densidad molar C o fracción molar y). Pero pueden obtenerse las
formulaciones sobre una base molar usando la relación 𝐶 𝜌/𝑀, donde M es la masa molar.
También, por facilidad, se restringirá la atención a la convección en fluidos que son (o pueden
tratarse como) mezclas binarias. Considérese el flujo de aire sobre la superficie libre de una masa
de agua, como un lago, en condiciones isotérmicas. Si el aire no está saturado, la concentración
del vapor de agua variará desde un máximo en la superficie del agua, en donde el aire está
siempre saturado, hasta el valor en el flujo libre, lejos de la superficie. En la convección de calor,
la región en la cual existen los gradientes de temperatura se definió como la capa límite térmica.
De manera análoga, en la convección de masa, se define la región del fluido en la cual existen
gradientes de concentración como la capa límite de concentración, como se muestra en la figura
14-38. En el flujo externo, el espesor de la capa límite de concentración, 𝛿𝐶 , para una especie A,
en un lugar especificado sobre la superficie, se define como la distancia normal y desde la
superficie, a la cual

𝜌𝐴 , 𝑠 − 𝜌𝐴
= 0.99
𝜌𝐴 , 𝑠 − 𝜌𝐴 , ∞
Fundamentos de Transferencia de Calor y
Materia:
Masa

Tarea: Tarea 7

donde 𝜌𝐴 , 𝑠𝑦𝜌𝐴 , 8 son las densidades de la especie A en la superficie (en el lado del fluido) y en
el flujo libre, respectivamente. En el flujo interno, se tiene una región de entrada de concentración,
donde se desarrolla el perfil de concentración, además de las regiones hidrodinámica y térmica
(figura 14-39). La capa límite continúa desarrollándose en la dirección del flujo hasta que su
espesor llega al centro del tubo y se juntan las capas límite. La distancia medida desde la entrada
del tubo hasta el lugar donde ocurre esta unión se llama longitud de entrada de concentración,
Lc, y la región que se encuentra más allá de ese punto se conoce como la región completamente
desarrollada, la cual se caracteriza por

𝜕 𝜌𝐴 , 𝑠 − 𝜌𝐴
( )=0
𝜕𝑥 𝜌𝐴 , 𝑠 − 𝜌𝐴,𝑏

donde 𝜌𝐴,𝑏 es la densidad media de la masa de la especie A, definida como

1
𝜌𝐴,𝑏 = ∫ 𝜌 𝑢𝑑𝐴𝐶
𝐴𝐶 𝑉𝑝𝑟𝑜𝑚 𝐴𝐶 𝐴

Por lo tanto, el perfil de diferencia de concentración, en forma adimensional, así como el


coeficiente de transferencia de masa permanecen constantes en la región completamente
desarrollada. Esto es análogo a los coeficientes de fricción y de transferencia de calor que
permanecen constantes en la región completamente desarrollada. En la convección de calor, las
magnitudes relativas de la cantidad de movimiento y la difusión de calor en las capas límite de
velocidad y térmica se expresan por el número de Prandtl, adimensional, definido como (figura
14-40)

𝑣 Difusividad de la cantidad de movimiento


𝑁ú𝑚𝑒𝑟𝑜 𝑑𝑒 𝑃𝑟𝑎𝑛𝑑𝑡𝑙: 𝑃𝑟 = =
∝ Difusividad térmica

La cantidad correspondiente en la convección de masa es el número de Schmidt, definido como

𝑣 Difusividad de la cantidad de movimiento


𝑁ú𝑚𝑒𝑟𝑜 𝑑𝑒 𝑆𝑐ℎ𝑚𝑖𝑑𝑡: 𝑆𝑐 = =
𝐷𝐴𝐵 Difusividad de masa

el cual representa las magnitudes relativas de la cantidad de movimiento molecular y la difusión


de masa en las capas límite de velocidad y de concentración, respectivamente. En el flujo laminar,
el crecimiento relativo de las capas límite de velocidad y térmica lo rige el número de Prandtl, en
tanto que el crecimiento relativo de las capas límite de velocidad y de concentración lo rige el
número de Schmidt. Un número de Prandtl cercano a la unidad (Pr 1) indica que la cantidad de
movimiento y la transferencia de calor por difusión son comparables, y las capas límite de
velocidad y térmica casi coinciden entre sí. Un número de Schmidt cercano a la unidad (Sc1)
Fundamentos de Transferencia de Calor y
Materia:
Masa

Tarea: Tarea 7

indica que la cantidad de movimiento y la transferencia de masa por difusión son comparables, y
las capas límite de velocidad y de concentración casi coinciden entre sí.

La transferencia de masa debida a la convección consiste en la transferencia entre un fluido en


movimiento y una superficie entre dos fluidos en movimiento, relativamente no miscibles. Este
modo de transferencia depende, tanto de las propiedades de transferencia como de las
características dinámicas del fluido que está fluyendo. Se debe hacer una distinción entre los
tipos de flujo, ya como se hizo en el caso de la transferencia convectiva de calor. Cuando una
bomba u otro aparato semejante ocasiona el movimiento del fluido, el proceso se llama
convección forzada. Si el movimiento del fluido se debe a una diferencia de densidades, que
puede haber surgido como resultado de una concentración o una diferencia de temperatura, el
proceso se llama convección libre o natural. De igual manera, como se hizo en el caso de la
transferencia de masa molecular, la transferencia convectiva de masa se lleva a cabo en la
dirección de una concentración que disminuye.

El coeficiente de transferencia convectiva de masa, 𝑘𝐶, estará en función del flujo de masa y la
diferencia de concentración desde que inicie hasta el final de la trayectoria que presente la
difusión. Es decir, dicho coeficiente incluirá características de ambas regiones, flujo laminar y
turbulento, del fluido en cualesquiera proporciones en las que se encuentren. Siempre existirá
una capa, a veces extremadamente delgada, cercana a la superficie, donde el flujo laminar y las
partículas próximas a la frontera sólida se encuentran en reposo. Como esto siempre se cumple,
el mecanismo de transferencia de masa entre una superficie y un fluido debe incluir una
transferencia de masa molecular a través de las capas estancada y laminar de fluido.

La resistencia que controla la transferencia de convectiva de masa es, a menudo, el resultado de


esta “película” de fluido y el coeficiente 𝑘𝐶 se llama, de acuerdo con esto, coeficiente de la
película. Si el flujo de fluido es laminar, entonces toda la transferencia entre la superficie y el fluido
en movimiento se llevará a cabo por medios moleculares. Por otra parte, el flujo de fluido es
turbulento, habrá un movimiento físico de paquetes de materia a lo largo de las líneas de flujo
transportadas por los remolinos presentes en el flujo turbulento. Así como en el caso de la
transferencia de calor, una gran rapidez de transferencia de calor asocia condiciones turbulentas.
La distinción entre los flujos laminar y turbulento será de vital importancia en cualquier situación
convectiva.

Parámetros importantes en la transferencia convectiva de masa

A menudo se utilizan parámetros adimensionales para relacionar los datos relativos a la


transferencia convectiva. En la transferencia de momento, nos hemos topado con los números
de Reynolds y de Euler. En la relación de datos de transferencia convectiva de calor, son
importantes los números de Prandtl y Nusselt. Algunos de los mismos parámetros junto con
Fundamentos de Transferencia de Calor y
Materia:
Masa

Tarea: Tarea 7

ciertas razones adimensionales recién definidas, resultaran de utilidad en la correlación de los


datos correspondientes a la transferencia de masa. El número de Schmidt tiene una importancia
en la transferencia de masa convectiva de masa, análoga a la del número de Prandtl, en la
transferencia convectiva de calor. Por otra parte, tenemos el número de Lewis, este aparece
cuando un proceso consta de la transferencia convectiva simultanea de masa y energía. Los
números de Schmidt y Lewis son combinaciones de las propiedades de los fluidos, y, por lo tanto,
cada uno de ellos se puede considerar como una propiedad del sistema en difusión.

Análisis dimensional de la transferencia convectiva de masa

El análisis adimensional predice los diversos parámetros adimensionales que resultan útiles en
la correlación de los datos experimentales. Hay dos procesos importantes de transferencia de
masa: la transferencia de masa hacia una corriente que fluye en condiciones de convección
forzada y de transferencia de masa hacia una fase que se mueve en condiciones de condiciones
de convección natural.

Transferencia a una fase cuyo movimiento se debe a la convección natural

Las corrientes de convección natural se crearán si existen cualquier variación en la densidad en


una fase liquida o en una fase gaseosa. La variación en la densidad se puede deber a diferencias
de temperatura o a diferencias de concentración, relativamente grandes. En el caso en que la
convección natural incluye transferencia de masa desde una pared vertical plana hasta un fluido
adyacente, las variables diferirán de las utilizadas en el análisis de la convección forzada.

Consideraciones acerca del flujo turbulento

En la mayor parte de las aplicaciones prácticas, el flujo, en la corriente principal, es turbulento en


lugar de ser laminar. Aunque muchos investigadores han hecho grandes aportaciones en la
comprensión del flujo turbulento, hasta el momento nadie ha tenido éxito en la predicción de los
coeficientes de transferencia convectiva. Es un esfuerzo para caracterizar este tipo de
movimiento.

Modelos para explicar la transferencia convectiva de masa

Teoría de película ficticia. Se basa en la presencia de una película ficticia de fluido en flujo laminar,
próxima a la frontera, que ofrece una resistencia a la transferencia de masa, igual a la que existe
verdaderamente en todo el fluido en movimiento. En otras palabras, se supone que toda la
resistencia a la transferencia existe en una película ficticia en la que la transferencia se realiza
totalmente por medio de difusión molecular. El espesor de la película, 𝛿, debe extenderse más
allá de la subcapa laminar para incluir la resistencia equivalente que se encuentra cuando la
Fundamentos de Transferencia de Calor y
Materia:
Masa

Tarea: Tarea 7

concentración cambia dentro de la capa de amortiguamiento y el núcleo turbulento. La teoría de


la penetración explica la transferencia de masa en la fase liquida durante la absorción del gas.
Se ha aplicado al flujo turbulento y también cuando la componente en difusión solo penetra a
corta distancia en la fase de interés a causa de su desaparición rápida por medio de una reacción
química o del tiempo, relativamente corto, de contacto.

Determinación de coeficientes fenomenológicos y de interfase

En términos generales, estos coeficientes dependen de la presión, temperatura y concentración


de un determinado sistema, por lo que no influye la fuerza impulsora en el cálculo de estos. El
coeficiente fenomenológico predominante en el flujo de transferencia de materias es el de
difusión. En algunas ocasiones el coeficiente de difusión no depende de la densidad, sino del
recorrido medio y de la velocidad media que estén presentes. En este caso se puede predecir
una dependencia del coeficiente con la presión y la temperatura, con esto nos referimos a que si
se hace un aumento en la temperatura existirá un aumento en la velocidad media, como
consecuencia la difusión será más rápida. Si tenemos más de una fase, en cada fase existirá un
gradiente de concentración ocasionando la transferencia de masa.

La interfase está en equilibrio siendo la resistencia muy baja para que exista transferencia de
materia en dado caso de no haber reacción. Si nuestro sistema presenta reacción química
debemos tener en consideración la resistencia presente de la interfase. La determinación de estos
coeficientes no es fácil siendo el método para determinarlo experimental. Para hacer uso de estos
coeficientes, se deberá hacer un análisis correcto del sistema a tratar, tener fundamentos los
conocimientos básicos presentes en el fenómeno. Si se tiene cuidado en el análisis y uso de los
coeficientes fenomenológicos nuestra modelación será correcta.

Conclusión

Es muy importante el estudio del principio de transferencia de masa por convección forzada, estos
parámetros son fundamentales para describir la transferencia de masa por convección y el
modelo propuesto para explicar el mecanismo de transferencia de masa por convección. Se ha
visto que la transferencia de masa por convección está estrechamente relacionada con la
dinámica del fluido en movimiento, especialmente la dinámica del fluido cerca del límite. Porque
el mecanismo de transmisión y la transmisión de par, energía y masa tienen grandes similitudes.
No hace falta decir que el mismo método que se utilizó originalmente para analizar el coeficiente
de transferencia de calor por convección se puede utilizar para evaluar el coeficiente de
transferencia de masa por convección. La transferencia de masa por convección puede existir en
flujos laminares y turbulentos. Este estudio introdujo muy brevemente la analogía de la
transferencia de masa por convección.
Fundamentos de Transferencia de Calor y
Materia:
Masa

Tarea: Tarea 7

Problemas

Problema 1. Una partícula de carbón pulverizado se quema en aire a 1145 k. el proceso se


encuentra limitado por la difusión del contraflujo de oxígeno a monóxido de carbono, CO, formado
en la superficie de la partícula. Suponga que el carbón está formado por carbono puro con una
densidad de 1280 𝑘𝑔/𝑚 3 y que la partícula es de forma esférica y tiene un diámetro inicial de
0.015 cm. Hay aire (21 % de oxígeno) a varios diámetros de la esfera. En las condiciones de
proceso de combustión, la difusividad del oxígeno en la mezcla gaseosa se puede suponer de
10^-4 𝑚 2 /𝑠𝑒𝑔. Determine el tiempo necesario para reducir el diámetro de la partícula a 0.005 cm.

Estado pseudo estacionario:

𝐶 + 𝑂 = 𝐶𝑂

𝜕𝐶𝐴
+ ∇ ∗ 𝑁𝐴 − 𝑅𝐴 = 0
𝜕𝑡

1 𝑑
(𝑟 2 𝑁𝐴𝑟 ) = 0
𝑟 2 𝑑𝑟

𝑟 2 𝑁𝐴𝑟 = 𝐶𝑡𝑒.

𝑑𝑦𝑂
𝑁𝑂𝑟 = −𝐶𝐷𝑂𝑚 + 𝑦𝑂 (𝑁𝑂𝑟 + 𝑁𝐶𝑂𝑟 )
𝑑𝑟

𝑁𝐶𝑂_ 1
= − & 𝑁𝐶𝑂 = −1𝑁𝑂𝑟
𝑁𝑂 1

𝑑𝑦𝑂
[𝑁𝑂𝑟 = −𝐶𝐷𝑂𝑚 ] 4𝜋𝑟 2
𝑑𝑟

𝑑𝑦𝑂
𝑊𝐴𝑟 = −4𝜋𝑟 2 𝐶𝐷𝑂𝑚
𝑑𝑟
𝑅2 𝑌 𝑂2
𝑑𝑟
𝑊𝐴𝑟 ∫ = −4𝜋𝐶𝐷𝑂𝑚 ∫ 𝑑𝑦𝑂
𝑅1 𝑟2 𝑌𝑂 1

1 1
𝑊𝐴𝑟 ( − ) = −4𝜋𝐶𝐷𝑂𝑚 (𝑌𝑂2 − 𝑌𝑂2 )
𝑅1 𝑅2
−4𝜋𝐶𝐷𝑂𝑚 (𝑌𝑂2 − 𝑌𝑂1 )
𝑊𝑂𝑟 =
1 1
( − )
𝑅1 𝑅2
Fundamentos de Transferencia de Calor y
Materia:
Masa

Tarea: Tarea 7

𝑊𝑂𝑟 = −4𝜋𝐶𝐷𝑂𝑚 𝑅1 (−𝑌𝑂1 )

𝑊𝐶𝑟 = −[−4𝜋𝐶𝐷𝑂𝑚 𝑅1 (−𝑌𝑂1 )]


𝑚
𝑑𝑚𝑜𝑙 𝑑 (𝑚𝑜𝑙 ) 1 𝑑𝑚 1 𝑑 (𝑝𝑐 𝑉) 𝑃𝑐 𝑑𝑉
= = = =
𝑑𝑡 𝑑𝑡 𝑀𝑐 𝑑𝑡 𝑀𝑐 𝑑𝑡 𝑀𝑐 𝑑𝑡

𝑃𝑐 𝑑 4 3 4𝜋𝑃𝑐 𝑑𝑟 3 4𝜋𝑃𝑐 3𝑟 2 𝑑𝑟
( 𝜋𝑟 ) = =
𝑀𝑐 𝑑𝑡 3 3𝑀𝑐 𝑑𝑡 3𝑀𝑐

4𝜋𝑃𝑐 𝑟 2 𝑑𝑟
𝑀𝑐 𝑑𝑡

4𝜋𝑃𝑐 𝑟 2 𝑑𝑟
= [4𝜋𝐶𝐷𝑂𝑚 𝑟(−𝑌𝑂1 )]
𝑀𝑐 𝑑𝑡

4𝜋𝑃𝑐 𝑟=𝑅2 𝑡=𝑡


∫ 𝑟 𝑑𝑟 = 4𝜋𝐶𝐷𝑂𝑚 𝑟(−𝑦𝑂1 ) ∫ 𝑑𝑡
𝑀𝑐 𝑟=𝑅1 𝑡=0

𝑃𝐶 (𝑅22 − 𝑅12 ) = 𝑀𝐶 𝐶𝐷𝑂𝑚 𝑟(−𝑦𝑂1 )𝑡

𝑃𝐶 (𝑅22 − 𝑅12 )
𝑡=
𝑀𝐶 𝐶𝐷𝑂𝑚 𝑟(−𝑦𝑂1 )

𝑃𝐶 (𝑅22 − 𝑅12 )
𝑡=
𝑀𝐶 𝑃𝐷𝑂𝑚 𝑟(−𝑦𝑂1 )

2
𝑘𝑔 pa m3 1.5𝑥10−4 𝑚
1280 (8.314 𝑚𝑜𝑙 𝑘 ) ( ) (1145 𝑘)
𝑚3 2
𝑡=
𝑔 𝑚2 1 𝑘𝑔
16 (101.325𝑥103 𝑝𝑎) (10−4 ) (1.5𝑥10−4 𝑚 ) ( )
𝑚𝑜𝑙 𝑠 1000𝑔

𝑡 = 2812.51 𝑠𝑒𝑔
Fundamentos de Transferencia de Calor y
Materia:
Masa

Tarea: Tarea 7

Problema 2 En la ruptura catalítica de los aceites de hidrocarburo los gases de los hidrocarburos
pesados se difunden en la superficie catalizadora donde se rompen, esto es, se descomponen.
El producto se vuelve a difundir en la corriente de gas. Una investigación cinética indicó que el
compuesto se rompe de acuerdo con la reacción siguiente:

𝐻 → 3𝐿

La reacción se efectúa en una partícula catalizadora esférica con una rapidez tan grande que la
difusión en la película estancada que circunda a la partícula controla la rapidez de la reacción.

Encuentre una expresión correspondiente a la rapidez de descomposición en función de las


propiedades de la fase gaseosa, de la concentración de los compuestos 𝐻 y 𝐿 en la fase gaseosa,
global, del diámetro de la película catalizadora y del espesor de la película estancada que
circunda al catalizador.

𝜕𝐶𝐴
+ 𝛻 ∙ 𝑁𝐴 − 𝑅𝐴 = 0
𝜕𝑡

Laplaciano para coordenadas esféricas

1 𝜕 1 𝜕 1 𝜕𝑁𝐴 𝜙
2
(𝑟 2 𝑁𝐴 𝑟) + (𝑁𝐴 𝜃𝑆𝑒𝑛𝜃) + − 𝑅𝐴 = 0
𝑟 𝜕𝑡 𝑟𝑆𝑒𝑛𝜃 𝜕𝜃 𝑟𝑆𝑒𝑛𝜃 𝜕𝜙

Eliminamos términos

𝜕𝐶𝐴 1 𝜕 1 𝜕 1 𝜕𝑁𝐴 𝜙
+ 2 (𝑟 2 𝑁𝐴 𝑟) + (𝑁𝐴 𝜃𝑆𝑒𝑛𝜃) + − 𝑅𝐴 = 0
𝜕𝑡 𝑟 𝜕𝑡 𝑟𝑆𝑒𝑛𝜃 𝜕𝜃 𝑟𝑆𝑒𝑛𝜃 𝜕𝜙

El resultado fue la siguiente ecuación diferencial

1 𝑑
(𝑟 2 𝑁𝐻 𝑟) = 0
𝑟 2 𝑑𝑡

Integramos la ecuación

∫ 𝑑(𝑟 2 𝑁𝐻 𝑟) = 0
Fundamentos de Transferencia de Calor y
Materia:
Masa

Tarea: Tarea 7

𝑟 2 𝑁𝐻 𝑟 = 𝐶

Aplicamos la ley de Fick

𝑑𝑦𝐻
𝑁𝐻 𝑟 = −𝐶𝐷𝐻𝐿 + 𝑦𝐻 (𝑁𝐻 𝑟 + 𝑁𝐿 𝑟)
𝑑𝑟

Hacemos el balance estequiométrico

𝑁𝐻 𝑟 −1
=
𝑁𝐿 𝑟 3

𝑁𝐿 𝑟 = −3𝑁𝐻 𝑟

Sustituimos en la Ley de Fick

𝑑𝑦𝐻
𝑁𝐻 𝑟 = −𝐶𝐷𝐻𝐿 + 𝑦𝐻 (𝑁𝐻 𝑟 − 3𝑁3 𝑟)
𝑑𝑟

𝑑𝑦𝐻
𝑁𝐻 𝑟 = −𝐶𝐷𝐻𝐿 − 2𝑦𝐻 𝑁𝐻 𝑟
𝑑𝑟

𝑑𝑦𝐻
𝑁𝐻 𝑟 + 2𝑦𝐻 𝑁𝐻 𝑟 = −𝐶𝐷𝐻𝐿
𝑑𝑟

𝑑𝑦𝐻
𝑁𝐻 𝑟(1 + 2𝑦𝐻 ) = −𝐶𝐷𝐻𝐿
𝑑𝑟

𝑑𝑦𝐻
−𝐶𝐷𝐻𝐿
𝑁𝐻 𝑟 = 𝑑𝑟
(1 + 2𝑦𝐻 )

Se multiplica por el área de una esfera, 4𝜋𝑟 2 , ya que la ecuación diferencial tiene un 𝑟 2 que no
tenemos, resultando la siguiente expresión

𝑑𝑦
−4𝜋𝑟 2 𝐶𝐷𝐻𝐿 𝑑𝑟𝐻
2
4𝜋𝑟 𝑁𝐻 𝑟 =
(1 + 2𝑦𝐻 )

El término 4𝜋𝑟 2 𝑁𝐻 𝑟 lo consideraremos como una constante 𝑊𝐻 𝑟

−4𝜋𝑟 2 𝐶𝐷𝐻𝐿 𝑑𝑦𝐻


𝑊𝐻 𝑟 =
(1 + 2𝑦𝐻 ) 𝑑𝑟
Fundamentos de Transferencia de Calor y
Materia:
Masa

Tarea: Tarea 7

Procedemos a integrar la ecuación resultante

𝑅2
𝑑𝑟 −4𝜋𝑟 2 𝐶𝐷𝐻𝐿 𝐻2 2𝑑𝑦𝐻
𝑊𝐻 𝑟 ∫ = ∫
𝑅1 𝑟2 2 𝐻1 1 + 2𝑦𝐻

El resultado de la integral es el siguiente

1 1 1 + 2𝑑𝑦𝐻2
𝑊𝐻 𝑟 ( − ) = −2𝜋𝐶𝐷𝐻𝐿 ln [ ]
𝑅1 𝑅2 1 + 2𝑑𝑦𝐻 1

Perfil de flujo

2𝜋𝐶𝐷𝐻𝐿 1 + 2𝑑𝑦𝐻2
𝑊𝐻 𝑟 = − ln [ ]
1 1 1 + 2𝑑𝑦𝐻 1
( − )
𝑅1 𝑅2

Para determinar el perfil de concentraciones, sustituimos nuestro valor 𝑁𝐻𝑟 en la ecuación


diferencial sin integrar

1 𝑑 2
(𝑟 𝑁𝐻 𝑟) = 0
𝑟 2 𝑑𝑟

1 𝑑 2 −𝐶𝐷𝐻𝐿 𝑑𝑦𝐻
[𝑟 ( )] = 0
𝑟 2 𝑑𝑟 1 + 2𝑦𝐻 𝑑𝑟

Integramos la ecuación

−𝐶𝐷𝐻𝐿 𝑑𝑦𝐻
∫ 𝒅 [𝑟2 ( )] = 0
1 + 2𝑦𝐻 𝑑𝑟

−𝐶𝐷𝐻𝐿 𝑑𝑦𝐻
𝑟2 ( ) = 𝐶1
1 + 2𝑦𝐻 𝑑𝑟

1 𝑑𝑦𝐻
𝑟2 ( ) = 𝐶1
1 + 2𝑦𝐻 𝑑𝑟

1 2𝑑𝑦𝐻 𝑑𝑟
∫ = 𝐶1 ∫ 2
2 1 + 2𝑦𝐻 𝑟

𝐶1
ln(1 + 2𝑦𝐻 ) = 2 + 2𝐶2
𝑟
Fundamentos de Transferencia de Calor y
Materia:
Masa

Tarea: Tarea 7

Condiciones limites

C. L. I r = R1 yH = y𝐻1

C. L. II r = R2 yH = y𝐻2

𝐶1
ln(1 + 2𝑦𝐻1 ) = 2 + 2𝐶2
𝑅1

𝐶1
ln(1 + 2𝑦𝐻2 ) = 2 + 2𝐶2
𝑅2

1 + 2𝑦𝐻 1 1
1
ln [ ] = 2𝐶1 ( − )
1 + 2𝑦𝐻 𝑅1 𝑅2
2

1 + 2𝑦 𝐻
ln [1 + 2𝑦 1 ]
𝐻2
∴ 𝐶1 = 1 1
2 (𝑅 − 𝑅 )
1 2

1 + 2𝑦𝐻
ln [1 + 2𝑦 1 ]
ln(1 + 2𝑦𝐻1 ) 𝐻2
∴ 𝐶2 = − 1 1
2 2 (𝑅 − 𝑅 ) 𝑅1
1 2

Sustituyendo los valores de 𝐶1 𝑦 𝐶2 en la ecuación original, tenemos lo siguiente

1 + 2𝑦 𝐻
1 + 2𝑦 𝐻
ln [1 + 2𝑦 1 ] ln [1 + 2𝑦 1 ]
𝐻2 2 𝐻
ln(1 + 2𝑦𝐻 ) = 1 1
+ ln(1 + 2𝑦𝐻1 ) − 1 1
(𝑅 − 𝑅 ) (𝑅 − 𝑅 ) 𝑅1
1 2 1 2

Simplificamos

1 1
1 + 2𝑦𝐻 1 + 2𝑦𝐻 (𝑟 − 𝑅 )
1 1
ln [ ] = ln [ ] 1 1
1 + 2𝑦𝐻 1 + 2𝑦𝐻 ( − )
1 2
𝑅1 𝑅2
Fundamentos de Transferencia de Calor y
Materia:
Masa

Tarea: Tarea 7

Perfil de concentraciones
1 1
( − )
𝑟 𝑅1
1 + 2𝑦𝐻 1+ 2𝑦𝐻1 ( 1 − 1 )
[ ]=[ ] 𝑅1 𝑅2
1 + 2𝑦𝐻1 1 + 2𝑦𝐻2

PROBLEMA 3: El compuesto hace transforma en el compuesto B en los sitios activos dentro de


un catalizador esférico. la rapidez de desaparición de A en la relación con la concentración De A
en el catalizador, Por medio de la reacción:
𝑅𝐴 = 𝑘𝑎𝐶𝐴
donde a es el área superficial del catalizador por unidad de volumen del mismo. Encuentra una
expresión que corresponda:
a) al perfil de concentraciones dentro del catalizador
b) el flujo molar de la superficie r=R
en función de un coeficiente de difusión efectiva de la concentración de a sobre la superficie de
la constante de rapidez de reacción y el radio del catalizador, R.

𝑑𝐶𝐴 1 𝑑(𝑟 2 𝑁𝐴𝑟 ) 1 𝑑(𝑁𝐴𝛳 𝑠𝑖𝑛 𝛳) 1 𝑑𝑁𝐴𝜙


+ 2 + + − 𝑅𝐴 = 0
𝑑𝑡 𝑟 𝑑𝑟 𝑟 𝑠𝑖𝑛 𝛳 𝑑𝑟 𝑟 𝑠𝑖𝑛 𝛳 𝑑𝜙

1 𝑑(𝑟 2 𝑁𝐴𝑟 )
− 𝑅𝐴 = 0 𝑅𝐴 = 𝐾𝑎𝐶𝐴
𝑟2 𝑑𝑟

1 𝑑(𝑟 2 𝑁𝐴𝑟 )
− 𝐾𝑎𝐶𝐴 = 0
𝑟2 𝑑𝑟

𝑑𝐶𝐴 𝑑𝐶𝐴
𝑁𝐴𝑟 = −𝐷𝐴𝐵 + 𝑦𝐴 (𝑁𝐴𝑟 + 𝑁𝐵𝑟 ) ⟹ 𝑁𝐴𝑟 = −𝐷𝐴𝐵
𝑑𝑟 𝑑𝑟

1 𝑑 𝑑𝐶𝐴
2
(−𝑟 2 𝐷𝐴𝐵 ) − 𝐾𝑎𝐶𝐴 = 0
𝑟 𝑑𝑟 𝑑𝑟

𝐷𝐴𝐵 𝑑 2 𝑑𝐶𝐴
− (𝑟 ) − 𝐾𝑎𝐶𝐴 = 0
𝑟 2 𝑑𝑟 𝑑𝑟
Fundamentos de Transferencia de Calor y
Materia:
Masa

Tarea: Tarea 7

𝐹 = 𝐶𝐴 𝑟

𝑑𝐹 𝑑𝐶𝐴
= 𝐶𝐴 + 𝑟
𝑑𝑟 𝑑𝑟

𝑑2 𝐹 𝑑𝐶𝐴 𝑑2 𝐶𝐴 𝑑𝐶𝐴 𝑑𝐶𝐴 𝑑2 𝐶𝐴


= + 𝑟 + = 2 + 𝑟
𝑑𝑟 2 𝑑𝑟 𝑑𝑟 2 𝑑𝑟 𝑑𝑟 𝑑𝑟 2

𝐷𝐴𝐵 2 𝑑2 𝐶𝐴 𝑑𝐶𝐴
− 2
[𝑟 2
+ 2𝑟 ] − 𝐾𝑎𝐶𝐴 = 0
𝑟 𝑑𝑟 𝑑𝑟

𝑑2 𝐹 𝑑𝐶𝐴 𝑑2 𝐶𝐴 𝑑2 𝐹 𝑑𝐶𝐴 2
𝑑2 𝐶𝐴
𝑟( = 2 + 𝑟 ) ⟹ 𝑟 = 2𝑟 + 𝑟
𝑑𝑟 2 𝑑𝑟 𝑑𝑟 2 𝑑𝑟 2 𝑑𝑟 𝑑𝑟 2

𝐷𝐴𝐵 𝑑2 𝐹
− [ 𝑟 2 ] − 𝐾𝑎𝐶𝐴 = 0
𝑟2 𝑑𝑟

𝐷𝐴𝐵 𝑑2 𝐹 𝑟
[− 2
− 𝐾𝑎𝐶𝐴 = 0] −
𝑟 𝑑𝑟 𝐷𝐴𝐵

𝑑2 𝐹 𝐾𝑎 𝐾𝑎
+ 𝐶 𝑟=0 𝑏=
𝑑𝑟 2 𝐷𝐴𝐵 𝐴 𝐷𝐴𝐵

𝑑2 𝐹
+ 𝑏𝐹 = 0
𝑑𝑟 2

𝐷2 + 𝑏 = 0 𝐷2 = ±√−𝑏 = ±√𝑏𝑖

𝐹 = 𝑐1 cos √𝑏 𝑟 + 𝑐2 sin √𝑏 𝑟

𝐶𝐿1 𝑟 = 0 𝐶𝐴 =? 𝐹=0

𝐶𝐿2 𝑟 = 𝑅 𝐶𝐴 = 𝐶𝐴𝑆 𝐹 = 𝐶𝐴𝑠 𝑅

0 = 𝑐1 cos √𝑏 (0) + 𝑐2 sin √𝑏 (0) ⟹ 𝑐1 = 0

𝐹 = 𝑐2 sin √𝑏 𝑟 𝐶𝐿2

𝐶𝐴𝑠 𝑅 = 𝑐2 sin √𝑏 𝑅
Fundamentos de Transferencia de Calor y
Materia:
Masa

Tarea: Tarea 7

𝐶𝐴𝑠 𝑅
𝑐2 =
sin √𝑏 𝑅

𝐶𝐴𝑠 𝑅
𝐹= sin √𝑏 𝑟
sin √𝑏 𝑅

Sustituyendo

𝐶𝐴𝑠 𝑅
𝐶𝐴 𝑟 = sin √𝑏 𝑟
sin √𝑏 𝑅

𝐶𝐴 𝑅 sin √𝑏 𝑟
=( )
𝐶𝐴𝑠 𝑟 sin √𝑏 𝑅

𝑑𝐶𝐴
[𝑁𝐴𝑟 = −𝐷𝐴𝐵 ] 4𝜋𝑟 2
𝑑𝑟

𝑑𝐶𝐴
𝜔𝐴𝑟 = −4𝜋𝑟 2 𝐷𝐴𝐵
𝑑𝑟

𝑑 𝐶𝐴 𝑅 𝑟√𝑏 cos √𝑏 𝑟 − sin √𝑏𝑟


( )= ( )
𝑑𝑟 𝐶𝐴𝑠 sin √𝑏 𝑅 𝑟2

𝑑𝐶𝐴 𝐶𝐴𝑠 𝑅[𝑟√𝑏 cos √𝑏 𝑟 − sin √𝑏𝑟]


=
𝑑𝑟 𝑟 2 sin √𝑏 𝑅

−4𝜋𝐷𝐴𝐵 𝐶𝐴𝑠 𝑅[𝑟√𝑏 cos √𝑏 𝑟 − sin √𝑏𝑟]


𝜔𝐴𝑟 =
sin √𝑏 𝑅

Problema 4. Calcule el tiempo requerido para la evaporación total de una esfera de naftalina
cuyo diámetro inicial es de 1 cm cuando está suspendida en una cantidad infinita de aire a 318
K. Si la temperatura de la superficie de la naftalina se supone igual a 318 k, su presión de vapor
de 1.06 𝑥 104 𝑃𝑎 y su densidad de 1.14 x 103 kg/m3. En estas condiciones, la difusividad del gas
de naftalina en aire es de 6.9 𝑥 10 − 7 𝑚2/𝑠𝑒𝑔.

𝐴𝑖𝑟𝑒

𝑇 = 318 𝐾
Fundamentos de Transferencia de Calor y
Materia:
Masa

Tarea: Tarea 7

𝑚2
𝐷𝐴𝐵 = 6.9𝑥10−7
𝑠
𝑟1 = 0.5 𝑐𝑚
𝑃𝐴 = 1.06𝑥104 𝑃𝑎
𝑘𝑔
𝑑𝑒𝑛𝑠𝑖𝑑𝑎𝑑 = 1.14𝑥103
𝑚3 Naftaleno

Con la ecuación de Fick:

𝑑𝑦𝐴
𝑁𝐴𝑍 = −𝐶𝐷𝐴𝐵
𝑑𝑟

𝑑𝑦𝐴
[𝑁𝐴𝑍 = −𝐶𝐷𝐴𝐵 ]𝐴
𝑑𝑟

𝑑𝑦𝐴
[𝑁𝐴𝑍 = −𝐶𝐷𝐴𝐵 ] (4𝜋𝑟 2 )
𝑑𝑟

𝑑𝑟
𝑁𝐴𝑍 = −𝐶𝐷𝐴𝐵 4𝜋 𝑑𝑦𝐴
𝑟2
𝑟2 𝑦𝐴2
𝑑𝑟
𝑁𝐴𝑍 ∫ 2 = −𝐶𝐷𝐴𝐵 4𝜋 ∫ 𝑑𝑦𝐴
𝑟1 𝑟 𝑦𝐴1

1 1
𝑁𝐴𝑍 ( − ) = −𝐶𝐷𝐴𝐵 4𝜋 (𝑦𝐴2 − 𝑦𝐴1 )
𝑟2 𝑟1

𝐶𝐷𝐴𝐵 4𝜋(𝑦𝐴1 − 𝑦𝐴2 )


𝑁𝐴𝑍 =
1 1
𝑟1 − 𝑟2

𝑦𝐴2 = 0; 𝑟2 = ∞

𝑁𝐴𝑍 = 𝐶𝐷𝐴𝐵 4𝜋𝑟1 𝑦𝐴1 (1)

Como hay variación del tamaño de la esfera, entonces usamos un modelo de régimen
pseudoestacionario.

𝑑𝑟
𝑁𝐴𝑍 = −𝐶𝐴 (4𝜋𝑟 2 )
𝑑𝑡

𝑁𝐴𝑍 𝑑𝑡 = −4𝜋𝐶𝐴 𝑟 2 𝑑𝑟
Fundamentos de Transferencia de Calor y
Materia:
Masa

Tarea: Tarea 7

𝑡 𝑟2
𝑁𝐴𝑍 ∫ 𝑑𝑡 = −4𝜋𝐶𝐴 ∫ 𝑟 2 𝑑𝑟
0 𝑟1

(𝑟2 3 − 𝑟1 3 )
𝑁𝐴𝑍 𝑡 = −4𝜋𝐶𝐴 𝑟2 = 0
3

4𝜋𝐶𝐴 𝑟1 3
𝑁𝐴𝑍 =
3𝑡

4𝜋𝐶𝐴 𝑟1 3 4𝜋𝐶𝐴 𝑟1 3
𝐶𝐷𝐴𝐵 4𝜋𝑟1 𝑦𝐴1 = ∴𝑡=
3𝑡 3𝐶𝐷𝐴𝐵 4𝜋𝑟1 𝑦𝐴1

𝑃𝐴 1.06𝑥104 𝑃𝑎
𝑦𝐴1 = = = 0.1046
𝑃 101 325 𝑃𝑎

𝑘𝑔
𝜌𝐴 1.14𝑥103 3
𝐶𝐴 = = 𝑚 = 8.858 𝑘𝑚𝑜𝑙
𝑀𝑀𝐴 𝑘𝑔 𝑚3
128.7
𝑚𝑜𝑙

𝑛 𝑃 1 𝑎𝑡𝑚 𝑚𝑜𝑙
𝐶= = = 3 = 38.35 3
𝑉 𝑅𝑇 𝑎𝑡𝑚 ∗ 𝐿 1𝑚 𝑚
(0.082 ) (318 𝐾) 𝑥 3
𝑚𝑜𝑙 ∗ 𝐾 10 𝑚𝑜𝑙

𝑘𝑚𝑜𝑙 2
(8.858
3 )(0.005 𝑚)
𝑡= 𝑚 = 26530.7 𝑠 = 7.4 ℎ
𝑚𝑜𝑙 𝑚2 1 𝑘𝑚𝑜𝑙
3(38.35 3 )(6.9𝑥10−7 𝑠 )(0.1046)( 3 )
𝑚 10 𝑚𝑜𝑙

Bibliografía
Bird, R., Stewart, W., & Lightfoot, E. (2006). Fenómenos de transporte. Limusa-Wiley.

Çengel, J. A. (s.f.). Transferencia de calor. McGraw-Hil.

También podría gustarte