Inge 2
Inge 2
Inge 2
TRANSFERENCIA DE MASA EN
MEDIO SOLIDO
Curso: Ingeniera de Alimentos II
Docente: Ing. Isabel Berrocal Martinez
Alumnos:
INTRODUCCION
Los fenmenos de transporte tienen lugar en aquellos procesos, conocidos como procesos de
transferencia, en los que se establece el movimiento de una propiedad ( masa, momentum o
energa) en una o varias direcciones bajo la accin de una fuerza impulsora. Al movimiento de una
propiedad se le llama flujo.
Los procesos de transferencia de masa son importantes ya que la mayora de los procesos
qumicos requieren de la purificacin inicial de las materias primas o de la separacin final
de productos y subproductos. Para esto en general, se utilizan las operaciones de transferencia de
masa.
Con frecuencia, el costo principal de un proceso deriva de las separaciones ( Transferencia de
masa). Los costos por separacin o purificacin dependen directamente de la relacin entre la
concentracin inicial y final de las sustancias separadas; s esta relacin es elevada, tambin
sern los costos de produccin.
Algunos de los ejemplos del papel que juega la transferencia de masa en los procesos industriales
son: la remocin de materiales contaminantes de las corrientes de descarga de los gases y aguas
contaminadas, la difusin de neutrones dentro de los reactores nucleares, la difusin de
sustancias al interior de poros de carbn activado, la rapidez de las reacciones qumicas
catalizadas
y
biolgicas
as
como
el
acondicionamiento
del aire,
etc.
En la industria farmacutica tambin ocurren procesos de transferencia de masa tal como la
disolucin de un frmaco, la transferencia de nutrientes y medicamento a la sangre, etc.
La ley de Fick es el modelo matemtico que describe la transferencia molecular de masa,
en sistemas o procesos donde puede ocurrir solo difusin o bien difusin mas conveccin. En
este trabajo, una idea central ser el clculo de los coeficientes de transferencia de masa para
diferentes
sistemas(
estados
de
agregacin
de
la materia ).
TRANSFERENCIA DE MASA
Clasificacin
general
de
la
transferencia
de
masa.
El mecanismo de transferencia de masa, depende de la dinmica del sistema en que se lleva
acabo.
Hay dos modos de transferencia de masa:
a.
molecular: La masa puede transferirse por medio del movimiento molecular fortuito en los
fluidos ( movimiento individual de las molculas ), debido a una diferencia de concentraciones.
La difusin molecular puede ocurrir en sistemas de fluidos estancados o en fluidos que se
estn moviendo.
b.
convectiva: La masa puede transferirse debido al movimiento global del fluido. Puede
ocurrir que el movimiento se efecte en rgimen laminar o turbulento. El flujo turbulento resulta
del movimiento de grandes grupos de molculas y es influenciado por las caractersticas
dinmicas del flujo. Tales como densidad, viscosidad, etc.
Usualmente, ambos mecanismos actan simultneamente. Sin embargo, uno puede ser
cuantitativamente dominante y por lo tanto, para el anlisis de un problema en particular, es
necesario considerar solo a dicho mecanismo. La transferencia de masa en slidos porosos,
lquidos y gases sigue el mismo principio, descrito por la ley de Fick.
El mecanismo real de transporte difiere en gran medida entre gases, lquidos y slidos, debido a
las diferencias sustanciales en la estructura molecular de estos 3 estados fsicos.
Gases: los gases contienen relativamente pocas molculas por unidad de volumen. Cada
molcula tiene pocas vecinas o cercanas con las cuales pueda interactuar y las fuerzas
moleculares son relativamente dbiles; las molculas de un gas tienen la libertad de
moverse a distancias considerables antes de tener colisiones con otras molculas.
El comportamiento ideal de los gases es explicado por la teora cintica de los gases.
Lquidos: los lquidos contienen una concentracin de molculas mayor por unidad de
volumen, de manera que cada molcula tiene varias vecinas con las cuales puede
interactuar y las fuerzas intermoleculares son mayores. Como resultado, el movimiento
molecular se restringe ms en un lquido. La migracin de molculas desde una regin
hacia otra ocurre pero a una velocidad menor que en el caso de los gases. Las molculas
de un lquido vibran de un lado a otro, sufriendo con frecuencia colisiones con las
molculas vecinas.
Slidos: En los slidos , las molculas se encuentran ms unidas que en los lquidos; el
movimiento molecular tiene mayores restricciones. En muchos slidos, las fuerzas
intermoleculares son suficientemente grandes para mantener a las molculas en
una distribucin fija que se conoce como red cristalina.
Y neto =
Segn la ecuacin (1), para que la densidad de flujo Y sea positiva, el gradiente dG /dz tiene que
ser negativo.
Los coeficientes de transferencia de masa tienen mucha importancia, por que al regular la rapidez
con la cual se alcanza el equilibrio, controlan el tiempo que se necesita para la difusin.
Los coeficientes de rapidez para los diferentes componentes en una fase dada difieren entre si en
mayor grado bajo condiciones en donde prevalece la difusin molcular. En condiciones de
turbulencia, en que la difusin molecular carece relativamente de importancia, los coeficientes de
transferencia se vuelven mas parecidos para todos los componentes. 1
= - C DAB
Donde
C: es la concentracin de A y B en mol Kg de (A + B) / m3
Xa: es la fraccin mol de A en la mezcla de A y B
JAZ: es el flujo de masa en molKg/(seg m2)
S c es constante, tenemos que cA = cxA
cdxA = d ( cxA ) = dcA
Entonces, para una concentracin total constante
= - C DAB
(2)
gradiente de concentracin puede considerarse por consiguiente como una fuerza impulsora. La
magnitud numrica de la difusividad indica la facilidad con que el componente A se transfiere en la
mezcla. Si la difusividad tiene un valor elevado, entonces hay mucha facilidad para el transporte de
masa.. El flujo del componente A se mide con relacin a la velocidad molar promedio de todos los
componentes.
El signo negativo hace hincapi que la difusin ocurre en el sentido del decremento en
concentracin, y el gradiente es negativo, pero el flujo de masa debe ser positivo. La difusividad es
una caracterstica de un componente y su entorno (temperatura, presin, concentracin; ya sea en
solucin lquida, gaseosa o slida y la naturaleza de los otros componentes)
Hasta ahora se ha considerado la ley de Fick para la difusin en un fluido estacionario; es decir , no
ha habido un movimiento neto ( o flujo convectivo ) de la totalidad de la mezcla A y B. El flujo
especfico de difusin JAZ se debe en este caso al gradiente de concentracin. La velocidad a la
cual los moles de A pasan por un punto fijo hacia la derecha, lo cual se tomar como flujo positivo.
Este flujo puede transformarse en una velocidad de difusin de A hacia la derecha por medio de la
expresin.3
Donde
n Ad: es la velocidad de difusin de A en m/seg.
Considerando ahora lo que sucede cuando la totalidad del fluido se mueve con un flujo general o
convectivo hacia la derecha. La velocidad molar promedio de la totalidad del fluido con respecto a
un punto estacionario es n M m/seg. El componente A sigue difundindose hacia la derecha, pero
ahora su velocidad de difusin n Ad se mide con respecto al fluido en movimiento. Para un
observador estacionario, A se desplaza con mayor rapidez que la fase total, pues su velocidad de
difusin n Ad se aade a la fase total n M. Expresada matemticamente, la velocidad de A con
respecto al punto estacionario es la suma de la velocidad de difusin y de la velocidad convectiva o
promedio2.
n A = n Ad + n M (4)
Expresndolo esquemticamente:
nA
n Ad n M
Multiplicando la ecuacin (4) por cA
cAn A = cAn Ad + cAn M (5)
Cada uno de estos 3 componentes es un flujo especfico. El primer trmino cAn A puede
representarse con el flujo NA en mol kg A / seg. m2. Este es el flujo total de A con respecto al punto
estacionario. El segundo trmino es JAZ ,esto es, el flujo especfico de difusin con respecto al
fluido en movimiento. El tercer trmino es el flujo convectivo de A con respecto al punto
estacionario.
Por
consiguiente
,
la
ecuacin
(5)
se
transforma
en:
NA = JAZ + cAn M (6)
Sea N el flujo convectivo total de la corriente general con respecto al punto estacionario. Entonces:
NA = cn M = NA + NB (7)
Despejando n M
n M = NA + NB / c (8)
Sustituyendo la ecuacin (8) en la ecuacin (6)
NA = JAZ +
( NA + NB ) (9)
Puesto que JAZ es la ley de Fick, por lo tanto la ecuacin (9) se transforma en la expresin general
para difusin mas conveccin2.
(10)
-DAB C
La cual describe la difusin a travs de una superficie fija en el espacio; en esta ecuacin, los
efectos del flujo global y el de la difusin molecular estn representados por el primer y segundo
trmino respectivamente.
Desde el punto de vista matemtico, esta ecuacin posee una estructura vectorial, y la direccin
del flujo global por unidad de rea, o sea, el primer trmino coincide con la direccin del gradiente .
El signo negativo del segundo trmino solo indica una disminucin de la concentracin, dada por
xA en la direccin del gradiente.
Una vez analizada la ley de Fick, se observa la necesidad de disponer de valores numricos del
parmetro difusividad. En las siguientes secciones se discutirn sus diversos mtodos de clculo.
Difusividad
de
gases
La difusividad, o coeficiente de difusin es una propiedad del sistema que depende de la
temperatura , presin y de la naturaleza de los componentes. Las expresiones para calcular la
difusividad cuando no se cuenta con datos experimentales, estn basadas en la teora cintica de
los
gases.
Hirschfelder, Bird y Spotz, utilizando el potencial de Lennard Jones para evaluar la influencia de las
fuerzas intermoleculares, encontraron una ecuacin adecuada al coeficiente de difusin
correspondiente a parejas gaseosas de molculas no polares, no reactivas a temperaturas y
presiones moderadas. Conocida como la ecuacin de Chapman-Enskog1,2,5
s AB = 1/2 ( s A + s B ) (12)
e AB = ( e Ae B )1/2 (13)
W D se calcula en funcin de KT/e AB donde K es la constante de Boltzmann y e AB es la energa
de interaccin molecular correspondiente al sistema binario AB
Hay tablas y apndices que tabulan estos valores. En ausencia de datos experimentales, los
valores de los componentes puros se pueden calcular a partir de las siguientes relaciones
empricas.
s = 2.44
1/3 (16)
Donde:
Vb = volmen molecular en el punto normal de ebullicin, en cm3 / g mol
Vc = volmen molecular crtico, en cm3 / g mol
Tc = temperatura crtica en grados kelvin
Pc = presin crtica en atmsferas
lugar en la red para ocupar una vacante cercana (creando un nuevo sitio vaco
en su posicin original en la red). Se presenta un reflujo de tomos y vacantes.
Difusin intersticial: Un tomo se mueve de un intersticio a otro. Este
mecanismo no requiere de vacantes para llevarse acabo. En ocasiones un
tomo sustitucional deja su lugar en la red normal y se traslada a un intersticio
muy reducido.
Difusin intersticial desajustada: Es poco comn, debido a que el tomo no se
ajusta o acomoda fcilmente en el intersticio, que es ms pequeo.
Intercambio simple: Puede darse el intercambio simple entre tomos o por
medio del mecanismo cclico( desplazamiento circular).
Aponte, 2005). Tambin es comn utilizar modelos fenomenolgicos muy simples, slo vlidos
para condiciones operativas similares a las experimentales (Spiazzi y Mascheroni, 1997).
METODOLOGIA
DESHIDRATACION OSMOTICA
Ambos tipos de frutas fueron seleccionadas y lavadas en agua corriente para remocin de
suciedad superficial. Luego, fueron peladas manualmente y cortadas en rodajas de 1 cm de
espesor, utilizando cuchillo de acero inoxidable. Las muestras fueron pesadas y colocadas en
vasos
de
precipitado
con
una
solucin
deshidratante
a
30
C.
La relacin muestra: solucin utilizada fue de 1:10 (p/p) para garantizar una concentracin de
solutos prcticamente todo el proceso.
Luego del tiempo de inmersin requerido, las muestras fueron drenadas, colocadas sobre papel
absorbente para remocin de solucin en exceso y pesadas nuevamente para la determinacin de
prdida de peso en relacin a la masa inicial del producto. La evolucin de la transferencia de
masa
durante el proceso fue monitoreada a travs de la variacin en el tiempo de la prdida de peso
(WR), contenido de slidos totales (TS) y de slidos solubles (Brix). La determinacin de estos
parmetros fue realizada por duplicado, utilizando las siguientes expresiones:
Reduccin de Peso (WR): fue calculada en una balanza semi analtica a partir de la
WR ( )=
(mimf )
x 100
mi
..(1)
Dnde: WR: prdida de peso de producto deshidratado (%); mi: masa inicial de la muestra fresca y
mf: masa de muestra deshidratada osmticamente a tiempo t.
Slidos Totales (TS): fueron determinados en una estufa de vaco a 70 0C (A.O.A.C., 2002)
hasta alcanzar un peso constante.
ms
( mo
) x 100
TS ( )=
. (2)
A partir de estos valores experimentales, se evaluaron los parmetros cinticos del proceso de DO:
WL ( ) =
m wo mwf
x 100 . (3)
mwo
siendo: WL = prdida de agua del producto deshidratado (%); mwo = masa inicial de agua de la
muestra (g) y mwf = masa final de agua de la muestra (g)
MODELOS UTILIZADOS
El primer modelo utilizado para ambos productos, est basado en la Ley de Difusin de Fick para
estado no estacionario, que calcula la cantidad de agua que abandona la fruta y los slidos
solubles
que
difunden
hacia
la
fruta,
como
una
funcin
del
tiempo.
Las siguientes suposiciones fueron usadas en el desarrollo del modelo: las rodajas de kiwi y pera
de 1 cm de espesor se consideraron placas infinitas; las concentraciones iniciales de agua y
slidos solubles en ambas frutas fueron uniformes; el proceso es isotrmico (30 0C); el coeficiente
de difusin aparente es constante (D f(C)); flujos simultneos en contra corriente: difusin de
agua desde la fruta y difusin de azcar hacia la fruta son solamente considerados, otros
mecanismos de transferencia son despreciados; se considera que un film de solucin de azcar
est adherido a la superficie de la fruta como una capa lmite (se asume que este film est a la
concentracin de equilibrio y que el proceso procede directamente por difusin); se desprecia el
encogimiento. Este fenmeno, en difusin unidimensional puede describirse con la Ley de Fick
para
estado
no
estacionario
dado
por
las
siguientes
expresiones:
PARA AGUA
C
C
=Dw
t
x2
. (5)
C: concentracin;
D: coeficiente de difusin efectivo; t: tiempo; x: distancia desde el centro de la placa y los
suprandices: w: agua y s: slidos solubles, siendo sacarosa en este caso.
Condiciones iniciales:
C=CO en t=O -L/2 < x< +L/2
Condiciones de contorno:
C = C1 para t > 0; x = L/2
donde
CO y Cl son las concentraciones inicial y global (promedio en todo el volumen de solucin),
respectivamente. La solucin de las ecuaciones (5) y (6) promediada en el volumen, considerando
que hay algo de soluto inicialmente presente en el slido, puede escribirse como (Crank, 1975):
2
W W 8
1
2
w
= 2
exp
D
(
2
n+
1
)
t
W iW n=0 ( 2 n+1 )2
4 l2
2
SS 8
1
2
s
= 2
exp
D
(
2
n+1
)
t
2
S oS n=0 ( 2n+ 1 )2
4l
(7)
.. (8)
Donde:
Dwy Ds: coeficientes de difusin efectivos para el agua y el azcar respectivamente;
I: medio espesor de la rodaja de fruta = l/2;
W: contenido medio de agua en la fruta a tiempo t;
Wi: contenido medio de agua en la fruta a tiempo t = 0 (contenido inicial de humedad);
W: contenido medio de agua en la fruta a tiempo t = (concentracin de equilibrio);
S: contenido medio de sacarosa en la fruta a tiempo t; SO: contenido medio de sacarosa en la fruta
a tiempo t = 0 (concentracin inicial) y
S: contenido medio de sacarosa en la fruta a tiempo t = (concentracin de equilibrio).
El segundo modelo utilizado, desarrollado por Azuara y col, (1992), se us slo para determinar las
concentraciones del film en la capa lmite en la superficie de cada fruta, asumiendo que est a la
concentracin de equilibrio de la solucin durante el proceso de smosis. En este modelo, los
autores realizaron un balance de masa sobre el movimiento de agua dentro del alimento,
obtenindose la siguiente expresin de la velocidad de prdida de agua WL como una funcin del
tiempo:
WL=
S1 tW l
1+ S 1 t
dnde: WL: fraccin de la prdida de agua en el equilibrio; SI: constante del modelo relacionado
con la velocidad de prdida de agua y t: tiempo.
La ecuacin anterior se linealiza mediante la siguiente expresin:
t
1
t
=
+
WL S1 tW L W L
De forma similar, para la Ganancia de Slidos SG se obtiene:
t
1
t
=
+
SG S 2 S G S G
donde:
SG: fraccin de la ganancia de slidos solubles en el equilibrio y S2: constante del modelo
relacionado con la velocidad de ganancia de slidos.
Utilizando los valores calculados de WL y SG obtenidos de los datos experimentales para cada
tiempo de proceso estudiado, la prdida de agua en el equilibrio (WL) y la constante SI fueron
estimados de la pendiente e intercepcin del grfico WIL vs t, mientras que la ganancia de slidos
en el equilibrio (SG' y la constante S2 fueron determinados de la pendiente e intercepcin del
grfico t/SG vs t.
C provistas por el software Excel, eligiendo aquella que tuviese mayor coeficiente de correlacin.
ANALISIS Y RESULTADOS
En las Figuras 1-4 se observan las evoluciones de los parmetros: WL y SG de las rodajas de pera
y kiwi, respectivamente, durante las experiencias de DO, utilizando sacarosa como agente
De acuerdo a las Figuras 1 y 3, se observa una mayor velocidad de Prdida de Agua WL en ambas
frutas durante las 4 primeras horas de proceso en comparacin con la prdida de agua alcanzada
en las 72 horas de deshidratacin, principalmente en los ensayos realizados con sacarosa a 600
Brix, consecuencia de una mayor fuerza osmtica de deshidratacin existente entre el alimento y la
solucin hipertnica al inicio del proceso. Adems, en ambos ensayos, se observ una tendencia a
la estabilizacin hacia las 72 horas de proceso, fundamentalmente debido a una aproximacin a la
condicin de equilibrio. El aumento de la concentracin de solucin osmtica favorece la prdida
de agua en el producto debido al mayor gradiente de presin osmtica en la interface
producto/solucin, proporcionando una mayor velocidad de transferencia de masa.
Con relacin a la Ganancia de Slidos SG, en las Figuras 2 y 4, se observa que un incremento en
la concentracin de sacarosa aumenta significativamente el ingreso de slidos en ambos
productos. El efecto de la concentracin de la solucin osmtica sobre la transferencia de masa
est ms marcado para el caso de la pera, que, en el kiwi, donde se podra deducir que pasadas
las 72 horas la solucin ms concentrada (60 Brix) limit la entrada de slidos al producto por la
formacin
de
una
costra
superficial
de
azcar
alrededor
de
la
fruta.
Ambos resultados de WL y SG muestran las mismas tendencias tanto para la pera como para el
kiwi, pero en este ltimo caso la velocidad de ambos fenmenos es mucho mayor. Este
comportamiento puede atribuirse a las diferencias en la estructura celular de cada fruta.
(aproximacin para tiempos largos). Se observa que los coeficientes son claramente ms altos
para la solucin de 60 Brix, para sacarosa respecto de agua y para kiwi comparados con los de
pera.
Tabla 1. Valores de los coeficientes de difusin aparentes durante la DO de peras
Fruta y condiciones
Kiwi 40brix
Kiwi 60brix
Pera 40brix
Pera 60brix
D Agua (m2/s)
6.19x10-11(R2=0.867)
1.18x10-10(R2=0.973)
7.32x10-11(R2=0.987)
9.85x10-11(R2=0.955)
D Sacarosa (m2/s)
3.01x10-10(R2=0.957)
3.91x10-10(R2=0.836)
7.88x10-11(R2=0.992)
1.18x10-10(R2=0.882)
CONCLUSIONES
Las rodajas de peras sufren menor prdida de agua y ganancia de slidos que las de kiwis
en
iguales
condiciones
de
proceso
en
soluciones
de
sacarosa.
Tanto la solucin analtica de difusin como las regresiones empricas son adecuadas para
simular el proceso de DO.