Difusion de Masa

Ecuación diferencial de la transferencia de masa
Erick Salinas Ojeda1 Viviana Carrillo2, PhD. Verónica Pinos3

Universidad de Cuenca, Facultad de Ciencias Químicas, Carrera de Ingeniería Química
Asignatura: Transferencia de Masa, Cuenca – Ecuador, Fecha de entrega: 18 – Junio – 2018
1. INTRODUCCIÓN
En la presente investigación se abordara el tema de la ecuación diferencial de
transferencia de masa, este trabajo se hace con el objeto de afianzar y ampliar los
conocimientos adquiridos en la catedra impartida sobre el tema, para el desarrollo del
presente trabajo se realizó una consulta netamente bibliográfica tomando como medio
de investigación libros y artículos de internet. Primordialmente se aborda en donde se
establece la ecuación diferencial de la transferencia de materia y a partir de la misma se
realiza el análisis para diferentes compuestos, en mezclas binarias, a diferentes
concentraciones y a diferentes volúmenes de control (rectangulares, cilíndricos y
esfericos).
2. OBJETIVO GENERAL
Analizar y estudiar la ecuación de la transferencia de masa de una especie determinada
en un volumen de control diferencial de masa, enfocándose en la evolución de la
ecuación con la concentración de los componentes y el volumen donde se contienen las
especies que se difundirán, así establecer la ecuación de continuidad para compuestos
simples y para mezclas binarias en diferentes geometrías de volumen de control.
3. DESARROLLO DEL TRABAJO

3.1. DEFINICÓN DE TRANSFERENCIA DE MASA
Se entiende por transferencia de masa al proceso en el cual una especie de un
sistema de múltiples componentes fluye desde una región de concentración alta a una de
concentración baja en el caso de la difusión y puede darse por movimiento global del fluido
(convección). http://bdigital.unal.edu.co/6542/1/9589532349_Parte1.pdf
Es importante destacar que, las leyes establecidas para este tipo de transferencia,
muestran la relación entre el flujo de las sustancias difundidas y el gradiente de
concentración que lo produce. https://es.slideshare.net/w3r0o0zoffith/48385740-
principiosdetransferenciademasadifusionmolecular
3.2. ECUACIÓN DIFERENCIAL DE LA TRANSFERENCIA DE MASA
Se establece un balance de masa con un volumen de control diferencial el cual nos

proporcionara la ecuación de continuidad para una especie determinada. En la Ilustración 1 se
muestra un volumen de control diferencial.
GRAFICO
Tenemos que considerar al volumen en las 3 dimensiones del espacio por donde fluye
una mezcla que tiene un componente A. A partir de esto se obtiene la expresión del volumen de
control:
𝜕
∬ 𝜌(𝒗 ∗ 𝒏) 𝑑𝐴 + ∭ 𝜌𝑑𝑉 = 0 (1)
𝜕𝑡 𝑐.𝑣.
𝑐.𝑠.
𝑅𝑎𝑝𝑖𝑑𝑒𝑧 𝑛𝑒𝑡𝑎 𝑑𝑒
rapidez neta de 𝑅𝑎𝑝𝑖𝑑𝑒𝑧 𝑑𝑒 𝑝𝑟𝑜𝑑𝑢𝑐𝑐𝑖ó𝑛
𝑎𝑐𝑢𝑚𝑢𝑙𝑎𝑐𝑖ó𝑛 𝑑𝑒 𝐴
( emisión de masa de A ) + ( ) − (𝑞𝑢í𝑚𝑖𝑐𝑎 𝑑𝑒 𝐴 𝑑𝑒𝑛𝑡𝑟𝑜 𝑑𝑒𝑙 ) = 0
𝑑𝑒𝑛𝑡𝑟𝑜 𝑑𝑒𝑙 𝑣𝑜𝑙𝑢𝑚𝑒𝑛
𝑑𝑒𝑙 𝑣𝑜𝑙𝑢𝑚𝑒𝑛 𝑑𝑒 𝑐𝑜𝑛𝑡𝑟𝑜𝑙 𝑣𝑜𝑙𝑢𝑚𝑒𝑛 𝑑𝑒 𝑐𝑜𝑛𝑡𝑟𝑜𝑙
𝑑𝑒 𝑐𝑜𝑛𝑡𝑟𝑜𝑙
(2)
Seguidamente se explicará cada uno de los términos:
3.2.1. RAPIDEZ NETA DE EMISIÓN DE MASA DEL VOLUMEN DE CONTROL.- Está

determinada por la masa que se difunde en cada una de las superficies
perpendiculares al vector de flujo (𝒏𝒂 = 𝝆𝑨 𝒗𝑨 ) de la dirección que se tome; por
ejemplo en la dirección X sería: 𝑛𝐴,𝑥 ∆𝑦∆𝑧|𝑥 .
Las expresiones resultantes en los 3 ejes son:

 Dirección x: 𝑛𝐴,𝑥 ∆𝑦∆𝑧|𝑥+∆𝑥 − 𝑛𝐴,𝑥 ∆𝑦∆𝑧|𝑥 (3)
 Dirección y: 𝑛𝐴,𝑦 ∆𝑥∆𝑧|𝑦+∆𝑦 − 𝑛𝐴,𝑦 ∆𝑥∆𝑧|𝑦 (4)
 Dirección z: 𝑛𝐴,𝑧 ∆𝑥∆𝑦|𝑧+∆𝑧 − 𝑛𝐴,𝑧 ∆𝑥∆𝑦|𝑧 (5)
3.2.2. RAPIDEZ DE ACUMULACIÓN DE A.- Este término está dado por la siguiente
expresión:
𝜕𝜌𝐴
∆𝑥∆𝑦∆𝑧 (6)
𝜕𝑡
3.2.3. RAPIDEZ DE PRODUCCIÓN QUÍMICA DE A: Dentro del volumen de control se

forma el componente A con una rapidez de reacción 𝑟𝐴 ,la cual es directamente
proporcional a la masa que reacciona e inversamente proporcional al volumen que
la contiene en un tiempo establecido.
𝑟𝐴 ∆𝑥∆𝑦∆𝑧 (7)
Reemplazando los términos antes descritos en la ecuación 2 obtenemos la siguiente

expresión:
𝑛𝐴,𝑥 ∆𝑦∆𝑧|𝑥+∆𝑥 − 𝑛𝐴,𝑥 ∆𝑦∆𝑧|𝑥 + 𝑛𝐴,𝑦 ∆𝑥∆𝑧|𝑦+∆𝑦 − 𝑛𝐴,𝑦 ∆𝑥∆𝑧|𝑦 + 𝑛𝐴,𝑧 ∆𝑥∆𝑦|𝑧+∆𝑧
𝜕𝜌𝐴
− 𝑛𝐴,𝑧 ∆𝑥∆𝑦|𝑧 + ∆𝑥∆𝑦∆𝑧 − 𝑟𝐴 ∆𝑥∆𝑦∆𝑧 = 0 (8)
𝜕𝑡
La siguiente expresión se obtiene, dividiendo para el volumen de control y realizando
las correspondientes simplificaciones:
𝑛𝐴,𝑥 |𝑥+∆𝑥 − 𝑛𝐴,𝑥 |𝑥 𝑛𝐴,𝑦 |𝑦+∆𝑦 − 𝑛𝐴,𝑦 |𝑦 𝑛𝐴,𝑧 |𝑧+∆𝑧 − 𝑛𝐴,𝑧 |𝑧 𝜕𝜌𝐴
+ + + − 𝑟𝐴 = 0 (9)
∆𝑥 ∆𝑦 ∆𝑧 𝜕𝑡
Con aproximaciones a cero de ∆𝑋, ∆𝑌 𝑦 ∆𝑍 se obtiene la ecuación de continuidad para

el componente A.
𝜕 𝜕 𝜕 𝜕𝜌𝐴
𝑛 + 𝑛 + 𝑛 + − 𝑟𝐴 = 0 (10)
𝜕𝑥 𝐴,𝑥 𝜕𝑦 𝐴,𝑦 𝜕𝑧 𝐴,𝑧 𝜕𝑡
Dónde:
𝑛𝐴 = 𝜌𝐴 𝑣𝐴 𝑑𝑒𝑛𝑠𝑖𝑑𝑎𝑑 𝑑𝑒 𝑓𝑙𝑢𝑗𝑜 𝑑𝑒 𝑚𝑎𝑠𝑎
𝑁𝐴 = 𝐶𝐴 𝑣𝐴 𝑑𝑒𝑛𝑠𝑖𝑑𝑎𝑑 𝑑𝑒 𝑓𝑙𝑢𝑗𝑜 𝑑𝑒 𝑚𝑜𝑙𝑎𝑟
Para expresar en función de base molar se realiza el mismo procedimiento, resultando lo

siguiente:
𝜕 𝜕 𝜕 𝜕𝐶𝐴
𝑁𝐴,𝑥 + 𝑁𝐴,𝑦 + 𝑁𝐴,𝑧 + − 𝑟𝐴 = 0 (11)
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡
Puesto que los vectores de flujo descritos en la ecuación anterior son los componentes
rectangulares del vector global de flujo de masa (𝑛𝐴 ) puede expresarse de esta forma:
𝜕𝜌𝐴
∇. 𝑛𝐴 + − 𝑟𝐴 = 0 (12)
𝜕𝑡
El análisis para un componente B se realizará de la misma forma en la que se

desarrolló para el componente A.
Por lo tanto, para una mezcla binaria se tendrán las siguientes ecuaciones:
𝜕(𝜌𝐴 + 𝜌𝐵 )
∇ ∙ (𝑛𝐴 + 𝑛𝐵 ) + − (𝑟𝐴 + 𝑟𝐵 ) = 0 (13)
𝜕𝑡
𝑛𝐴 + 𝑛𝐵 = 𝜌𝐴 ∗ 𝑣𝐴 + 𝜌𝐵 ∗ 𝑣𝐵 = 𝜌𝑣 (14)
𝜌𝐴 + 𝜌𝐵 = 𝜌 (15)
𝑟𝐴 = −𝑟𝐵 (16)
Remplazando las expresiones (14), (15) y (16) en la ecuación (13), obtenemos la

ecuación de continuidad para la mezcla:
𝜕𝜌
∇ ∙ 𝜌𝑣 + =0 (17)
𝜕𝑡
Esta ecuación puede definirse en función de su derivada:

𝐷𝜌
+ 𝜌∇ ∙ 𝒗 = 0 (18)
𝐷𝑡
Además tanto para los componentes A y B la ecuación de continuidad puede ser

expresada en función de su derivada:
𝜌𝐷𝑤𝐴
+ ∇ ∙ 𝒋𝑨 − 𝑟𝐴 = 0 (19)
𝐷𝑡
Se puede expresar las ecuaciones en función de unidades molares tanto para los
componentes A y B así, como también para la mezcla, teniendo las siguientes expresiones:
COMPONENTE A COMPONENTE B PARA LA MEZCLA

𝜕𝑐𝐴 𝜕𝑐𝐵 𝜕(𝑐𝐴 + 𝑐𝐵 )
∇ ∙ 𝑁𝐴 + − 𝑅𝐴 ∇ ∙ 𝑁𝐵 + − 𝑅𝐵 ∇ ∙ (𝑁𝐴 + 𝑁𝐵 ) +
𝜕𝑡 𝜕𝑡 𝜕𝑡
=0 (20) =0 (21) − (𝑅𝐴 + 𝑅𝐵 ) = 0 (22)
Para la mezcla binaria de A y B, se tiene
𝑁𝐴 + 𝑁𝐵 = 𝑐𝐴 𝑣𝐴 + 𝑐𝐵 𝑣𝐵 = 𝑐𝑉 (23)
𝑐𝐴 + 𝑐𝐵 = 𝑐 (24)
Cuando se tiene ecuaciones reversibles se considera que 𝑅𝐴 = −𝑅𝐵 por lo que la

ecuación de continuidad en base molar es:
𝜕𝑐
∇ ∙ cV + − (𝑅𝐴 + 𝑅𝐵 ) = 0 (25)
𝜕𝑡
3.3. SEGUNDA LEY DE FICK
Esta ley describe el comportamiento que tienen las especies que se difunden en un
sistema respecto a su acumulación o pérdida, debido a que su flujo y gradiente de concentración
varían con el tiempo. Generalmente los procesos de difusión se dan en estado no estacionario,
en este caso resulta conveniente la aplicación de esta ley, expresada de la siguiente manera:
𝜕𝑐𝐴
= 𝐷𝐴𝐵 ∗ ∇2 ∗ 𝑐𝐴 (26)
𝜕𝑡
3.3.1. PARA DIFUSIÓN DE SÓLIDOS O PARA SISTEMAS DONDE NO INTERVIENE
LA CONVECCIÓN: se considera la difusión únicamente en la dirección Z y además
que 𝐷𝐴𝐵 no depende de la concentración, resultando la siguiente expresión:
𝜕𝑐𝐴 𝜕 2 𝑐𝐴
( ) = 𝐷𝐴𝐵 ∗ ( 2 ) (27)
𝜕𝑡 𝑧 𝜕𝑧 𝑡
3.3.2. DIFUSIÓN EN ESTADO ESTACIONARIO: Gradiente de concentración y flujo

constante, por lo tanto, se tienen las siguientes expresiones:
𝜕𝑐𝐴 𝜕 2 𝑐𝐴
=0 =0
𝜕𝑡 𝜕𝑧 2
𝜕𝑐𝐴
= 𝑐𝑡𝑒 → 𝐽 = 𝑐𝑡𝑒
𝜕𝑥
3.3.3. EN FUNCIÓN DE COORDENADAS RECTANGULARES
𝜕𝑐𝐴 𝜕 2 𝑐𝐴 𝜕 2 𝑐𝐴 𝜕 2 𝑐𝐴
= 𝐷𝐴𝐵 ∗ [ 2 + + ] (28)
𝜕𝑡 𝜕𝑥 𝜕𝑦 2 𝜕𝑧 2
3.3.4. EN FUNCIÓN DE COORDENADAS CILÍNDRICAS
𝜕𝑐𝐴 𝜕 2 𝑐𝐴 1 𝜕𝑐𝐴 1 𝜕 2 𝑐𝐴 𝜕 2 𝑐𝐴
= 𝐷𝐴𝐵 ∗ [ 2 + + + ] (29)
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝑟 𝑟 2 𝜕𝜃 2 𝜕𝑧 2
3.3.5. EN FUNCIÓN DE COORDENADAS ESFÉRICAS
𝜕𝑐𝐴 1 𝜕 𝜕𝑐𝐴 1 𝜕 𝜕𝑐𝐴 1 𝜕 2 𝑐𝐴

= 𝐷𝐴𝐵 ∗ [ 2 (𝑟 2 )+ 2 (𝑠𝑒𝑛𝜃 )+ 2 ] (30)
𝜕𝑡 𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝑠𝑒𝑛𝜃 𝜕𝜃 𝜕𝜃 𝑟 𝑠𝑒𝑛𝜃 𝜕∅2
Una vez descritas las ecuaciones de la Segunda Ley de Fick en los tres tipos de
coordenadas, se puede realizar el análisis y obtener la ecuación diferencial de transferencia de
masa para cada una de ellas. Siendo las ecuaciones correspondientes, las indicadas a
continuación:
3.3.6. COORDENADAS RECTANGULARES
𝜕𝑐𝐴 𝜕𝑁𝐴,𝑥 𝜕𝑁𝐴,𝑦 𝜕𝑁𝐴,𝑧

+[ + + ] = 𝑅𝐴 (31)
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
3.3.7. COORDENADAS CILÍNDRICAS

𝜕𝑐𝐴 1𝜕 1 𝜕𝑁𝐴,𝜃 𝜕𝑁𝐴,𝑧
+[ (𝑟𝑁𝐴,𝑟 ) + + ] = 𝑅𝐴 (32)
𝜕𝑡 𝑟 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧
3.3.8. COORDENADAS ESFÉRICAS
𝜕𝑐𝐴 1 𝜕 1 𝜕 1 𝜕𝑁𝐴,𝜃
+ [ 2 (𝑟 2 𝑁𝐴,𝑟 ) + (𝑁𝐴,𝜃 𝑠𝑒𝑛𝜃) + ] = 𝑅𝐴 (33)
𝜕𝑡 𝑟 𝜕𝑟 𝑟𝑠𝑒𝑛𝜃 𝜕𝜃 𝑟𝑠𝑒𝑛𝜃 𝜕∅
Estas ecuaciones pueden ser reestructuradas en sistemas en los cuales 𝜌 y 𝐷𝐴𝐵 sean
constantes, para lo cual se igualará la ecuación diferencial de transferencia de masa de cada
una de las coordenadas con las expresiones correspondientes de la Segunda Ley de Fick,
dichas ecuaciones son las siguientes:
3.3.9. COORDENADAS RECTANGULARES
𝜕𝑐𝐴 𝜕𝑐𝐴 𝜕𝑐𝐴 𝜕𝑐𝐴 𝜕 2 𝑐𝐴 𝜕 2 𝑐𝐴 𝜕 2 𝑐𝐴

+ [𝑣𝑥 + 𝑣𝑦 + 𝑣𝑧 ] = 𝐷𝐴𝐵 ∗ [ 2 + + ] + 𝑅𝐴 (34)
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 2 𝜕𝑧 2
3.3.10. COORDENADAS CILÍNDRICAS
𝜕𝑐𝐴 𝜕𝑐𝐴 1 𝜕𝑐𝐴 𝜕𝑐𝐴 1 𝜕 𝜕𝑐𝐴 1 𝜕 2 𝑐𝐴 𝜕 2 𝑐𝐴

+ [𝑣𝑟 + 𝑣𝜃 + 𝑣𝑧 ] = 𝐷𝐴𝐵 ∗ [ ( )+ 2 + ] + 𝑅𝐴 (35)
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧 𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝜃 2 𝜕𝑧 2
FORMAS ESPECIALES DE LA ECUACIÓN DIFERENCIAL DE LA

TRANSFERENCIA DE MASA
Para analizar perfiles de concentración en un sistema que se difunde se remplazan las siguientes
expresiones:
𝑛𝐴 = −𝜌𝐷𝐴𝐵 ∇𝑤𝐴 + 𝜌𝐴 𝑉 (36)

𝑁𝐴 = −𝑐𝐷𝐴𝐵 ∇𝑦𝐴 + 𝑐𝐴 𝑉 (37)
En la ecuación (12) se reemplazará con la ecuación (36 y 37), obteniendo lo siguiente:

𝜕𝜌𝐴
∇ ∙ 𝜌𝐷𝐴𝐵 ∇𝑤𝐴 + ∇ ∙ 𝜌𝐴 𝑣 + − 𝑟𝐴 = 0 (38)
𝜕𝑡
En la ecuación (20) se reemplazará con la ecuación (23), obteniendo lo siguiente:
𝜕𝑐𝐴
−∇ ∙ 𝑐𝐷𝐴𝐵 ∇𝑦𝐴 + ∇ ∙ 𝑐𝐴 𝑉 + − 𝑅𝐴 = 0 (39)
𝜕𝑡
4. CONCLUSIONES
Al analizar la ecuación diferencial de la transferencia de masa en un volumen de control
establecido en caso general un cubo, se establece que la velocidad de difusión es
directamente proporcional a la masa que se difunde en un volumen determinado
respecto al cambio de sección establecido por el tiempo en que se difunde la masa
gracias a esto y llevando a cada uno de los ejes de coordenadas (∆X,∆Y y ∆Z) a cero se
obtiene la ecuación de continuidad (Ecuación (10)); esta ecuación se puede expresar en
base molar (Ecuación (11)), en esta forma se puede modificar la ecuación de
continuidad de la transferencia de masa para diferentes concentraciones, diferentes
sustancias e incluso para mezclas binarias.
Para analizar la ecuación de continuidad de la transferencia de masa para distintos
volúmenes de control como pueden ser formas rectangulares, cilindricas y esféricas se
parte de la segunda ley de Fick en donde se analiza el comportamiento que tienen las
especies que se difunden en un sistema respecto a su acumulación o pérdida, debido a
que su flujo y gradiente de concentración varían con el tiempo; a partir de este análisis
podemos establecer las ecuaciones ((34) y (35)) de continuidad para la transferencia de
masa en las formas antes expuestas.
Existen diferentes formas especiales de esta ecuación pero la que nos pareció de mayor
importancia es la que relaciona los perfiles de concentración de transferencia de masa la
cual nos permite determinar la transferencia de masa en un rango determinado de flujo
másico y molar.
5. BIBLIOGRAFÍA.

Difusion de Masa

Cargado por

Copyright:

Formatos disponibles

Difusion de Masa

Cargado por

Información del documento

Descripción original:

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

Difusion de Masa

Cargado por

Copyright:

Formatos disponibles

Ecuación diferencial de la transferencia de masa

Erick Salinas Ojeda1 Viviana Carrillo2, PhD. Verónica Pinos3

3. DESARROLLO DEL TRABAJO

3.2. ECUACIÓN DIFERENCIAL DE LA TRANSFERENCIA DE MASA

Se establece un balance de masa con un volumen de control diferencial el cual nos

Seguidamente se explicará cada uno de los términos:

3.2.1. RAPIDEZ NETA DE EMISIÓN DE MASA DEL VOLUMEN DE CONTROL.- Está

Las expresiones resultantes en los 3 ejes son:

3.2.3. RAPIDEZ DE PRODUCCIÓN QUÍMICA DE A: Dentro del volumen de control se

Reemplazando los términos antes descritos en la ecuación 2 obtenemos la siguiente

Con aproximaciones a cero de ∆𝑋, ∆𝑌 𝑦 ∆𝑍 se obtiene la ecuación de continuidad para

Para expresar en función de base molar se realiza el mismo procedimiento, resultando lo

El análisis para un componente B se realizará de la misma forma en la que se

Remplazando las expresiones (14), (15) y (16) en la ecuación (13), obtenemos la

Esta ecuación puede definirse en función de su derivada:

Además tanto para los componentes A y B la ecuación de continuidad puede ser

COMPONENTE A COMPONENTE B PARA LA MEZCLA

Cuando se tiene ecuaciones reversibles se considera que 𝑅𝐴 = −𝑅𝐵 por lo que la

3.3.2. DIFUSIÓN EN ESTADO ESTACIONARIO: Gradiente de concentración y flujo

3.3.4. EN FUNCIÓN DE COORDENADAS CILÍNDRICAS

3.3.5. EN FUNCIÓN DE COORDENADAS ESFÉRICAS

𝜕𝑐𝐴 1 𝜕 𝜕𝑐𝐴 1 𝜕 𝜕𝑐𝐴 1 𝜕 2 𝑐𝐴

3.3.6. COORDENADAS RECTANGULARES

𝜕𝑐𝐴 𝜕𝑁𝐴,𝑥 𝜕𝑁𝐴,𝑦 𝜕𝑁𝐴,𝑧

3.3.7. COORDENADAS CILÍNDRICAS

3.3.9. COORDENADAS RECTANGULARES

𝜕𝑐𝐴 𝜕𝑐𝐴 𝜕𝑐𝐴 𝜕𝑐𝐴 𝜕 2 𝑐𝐴 𝜕 2 𝑐𝐴 𝜕 2 𝑐𝐴

3.3.10. COORDENADAS CILÍNDRICAS

𝜕𝑐𝐴 𝜕𝑐𝐴 1 𝜕𝑐𝐴 𝜕𝑐𝐴 1 𝜕 𝜕𝑐𝐴 1 𝜕 2 𝑐𝐴 𝜕 2 𝑐𝐴

FORMAS ESPECIALES DE LA ECUACIÓN DIFERENCIAL DE LA

𝑛𝐴 = −𝜌𝐷𝐴𝐵 ∇𝑤𝐴 + 𝜌𝐴 𝑉 (36)

En la ecuación (12) se reemplazará con la ecuación (36 y 37), obteniendo lo siguiente:

También podría gustarte