CAPITULO I

ELEMENTOS DE CINETICA

1.0 GENERALIDADES

En los procesos de obtención de un metal desde el mineral o concentrado hasta su

refinación, existen una o más etapas en las que ocurren transformaciones físico-
químicas. Solo contados metales, tales como oro o platino se encuentran nativos. En
todos los otros casos el metal debe ser químicamente extraído y purificado, por
ejemplo, en la obtención de cobre desde sulfuros u óxidos; en la obtención de arrabio
(hierro) y acero desde minerales oxidados de hierro; en la obtención de zinc desde
sulfuros o silicatos, etc.

En el caso mencionado del cobre, su extracción a partir de sus sulfuros ocurre por una
transformación inicial a alta temperatura para formar el eje fundido de cobre (Cu2S-
FeS) y luego al oxidar y reducir el eje formado a metal crudo o cobre blister (Fig. 1).

Procesos físicos (temp. amb.)

Minera 1 OPERACIONES CONCEN-

MECANICAS TRACION Relaves
Energía

Energía
Escoria Gases mecánica

1300 ºC 1250 ºC
REFINACION FORMACION
Cu Gases
PRIMARIA DEL EJE
blister escoria
(CONVERSION) (FUSION)

Chatarra Oxígeno Oxígeno Calor

Cu

Figura 1. Obtención de cobre crudo (cobre blister) desde sulfuros

Otro ejemplo típico lo representa la obtención de zinc desde sulfuros (Fig. 2)

1
 Procesos físicos

OPERACIONES CONCEN- Concentrado

Minera 1 de plomo (a
MECANICAS TRACION
Pb - Zn obtención de
Pb metálico)

Oxidante H2 SO4 Relaves

Energía Energía
Conc. Zn
Precipitación de: mecánica mecánica
Fe Aire
LIXIVIACION OXIDACION A
Co PURIFICACION 40 ºC OXIDO 950 ºC
As Gases

ELECTRO-
Electrolito Cátodos
OBTENCION
gastado de zinc

Energía eléctrica

Figura 2. Obtención de zinc desde minerales de plomo – zinc sulfurados.

En metalurgia, la mayor parte de las reacciones ocurren tanto en medio acuoso

(lixiviación, cementación, precipitación, electroobtención y electrorefinación), como a
alta temperatura (tostación, reducción, oxidación, etc.) o en medio fundido a alta
temperatura (obtención de aluminio, formación de ejes o matas fundidas de sulfuros
reducción de wustita (FeO) líquida con carbono, etc.).

El análisis inicial de la extracción, purificación o refinación de un metal debe ser

termodinámico, para establecer si la o las reacciones pueden ocurrir
termodinámicamente.

El segundo paso es establecer que ocurre durante la transformación química del o de

las substancias involucradas. Es necesario aquí conocer el o los mecanismo (s) que
participa (n), así como la velocidad de aproximación al equilibrio que se establece. Esto
corresponde a la cinética de reacción. Por otra parte es necesario conocer y
caracterizar el flujo y la mezcla entre reactivos y productos en el espacio físico en que
ocurre la transformación. Este espacio o volumen se denomina en general “reactor”,
aún cuando en metalurgia, y dependiendo del proceso, se conoce por horno, (cuba,
cubilote, rotativo, flash, CT, etc.) estanque (lixiviación en estanques o columnas,
estáticas o agitadas) convertidor (para ejes de cobre níquel), tostador (lecho fluido,
horno de piso múltiples), etc.

En el diseño de un reactor, o un proceso se comienza por analizar los aspectos

termodinámico y las condiciones bajo las cuales sería posible producir las
transformaciones deseadas, o bien, dadas ciertas condiciones de operación, verificar si
2
la transformación es posible de efectuar. Conocida las condiciones termodinámicas, se
debe entonces conocer el aspecto cinético del proceso: es decir, que equilibrio se
establece en las condiciones fijadas y con que velocidad se aproxima al equilibrio.
Finalmente, se debe conocer el patrón de flujo y de mezcla que ocurren en el interior
del reactor, e integrar la expresión analítica en la cinética para así caracterizar la
transformación global y diseñar el reactor o proceso.

1.1 CINETICA DE REACCIONES

Dependiendo de las fases que toman parte en una transformación química, las
reacciones se pueden clasificar en reacciones homogéneas y reacciones heterogéneas.

Reacciones homogéneas

Ocurren entre fases gaseosas, líquidas o sólidas. En metalurgia este tipo de reacción no
es tan frecuente como las heterogéneas, aún cuando existen varios casos de éstos:

- Gas-Gas, como ocurre en la etapa de reducción del Cu3Cl3 gaseoso a 700ºC, generado
en el proceso de segregación de minerales de cobre con el hidrógeno naciente
(atómico) en las superficies de las partículas de carbón, según:

Cu3Cl3(g) + 3H(g) 3Cu(s) + 3HCl(g)

- Líquido-líquido, como la oxireducción del óxido de cobre líquido a 1200ºC formado

durante el soplado del metal blanco, para la conversión de cobre:

2Cu2O()l + Cu2S()l ƒ 6Cuo(l) + 2SO2(g)

Otro caso de este tipo de reacciones Líquido-Líquido lo representa la refinación de un

metal fundido mediante una escoria líquida.

- Sólido-sólido, como en el caso de la formación de ferrita de zinc durante la tostación

de esfalerita (ZnS) a óxido de zinc:

ZnO(s) + 2Fe2O3(s) ƒ ZnO·2FeO(s)

Reacciones heterogénea

Son más frecuentes en metalurgia, tanto sólido-gas como sólido-líquido y líquido-gas,

como son por ejemplo:

- Sólido –Gas, como ocurre en la reducción de magnetita (Fe3O4) a wustita (FeO) en el

Alto Horno mediante el CO:

Fe3O4(s) + CO(g) ⇔ 3FeO(l) + CO2(g)

3
- Sólido-Líquido como en la lixiviación de un mineral de cobre con ácido sulfúrico:

CuSiO3 ·2H2O(s) + H2 SO4(aq) ƒ CuSO4(aq) + SiO2(s) + 3H2O(l)

O en la reducción de wustita líquida (1200ºC) con carbono sólido en el crisol del Alto
Horno:

FeO(l) + C(s) ƒ Fe(l) + CO(g)

- Líquido-Gas, como la conversión del eje de níquel o cobre a níquel o cobre metálicos:

NiS(l) + O2(g) ƒ Ni(l) + SO2(g)

Además de éstas reacciones heterogéneas, existen también algunas reacciones

catalíticas que tienen interés en metalurgia, por ejemplo, la transformación de SO2 a
SO3 mediante un catalizador (pentóxido de vanadio) según:

SO2(g) + 0.5 O2(g) ƒ SO3(g)

En ciertos casos, reacciones homogéneas en metalurgia requieren de otra fase sólida

líquida para nuclear, como es el caso de la oxidación del carbono a CO en la refinación
de acero.

Además de los factores termodinámicos, la cinética en algunos casos está asociada a la

superficie reaccional disponible. Por ejemplo, en la refinación de acero líquido con
oxígeno, la reacción de descarburación del acero según:

C(l) + 0.5O2(g) ƒ CO(g)

tiene una energía de Gibbs fuertemente negativa de ∆G1600ºC

0
= −65.9 kcal/mol y una
constante de equilibrio muy alta, sin embargo, se requiere una gran superficie de
metal líquido expuesto al oxígeno para que la reacción ocurra en forma rápida. Si el
oxígeno se aporta sólo desde la superficie del baño fundido la refinación puede
demorar varias horas, en tanto que al emplear una lanza refrigerada que inyecte el
oxígeno a alta velocidad al baño y disperse el metal líquido, la refinación no toma más
de 20 a 30 minutos. La razón se debe a que el fenómeno en su totalidad está
controlado por la transferencia de masa (oxígeno) desde la fase gaseosa hacia el baño,
como se verá mas adelante.

Otro tanto ocurre, por ejemplo, con la reacción de un sólido y un gas. Para obtener una
alta velocidad de reacción será necesario tener un sólido poroso o finamente dividido.

4
1.2 EQUILIBRIO QUIMICO

De físicoquímica se vio que para una transformación cualquiera

r1 &
A + B ¬¾r¾
¾
¾
&
® C+ D
2

las velocidades de reacción se pueden representar en la forma

n m dC A dCB
r&
1 = kr1 C A CB = - =-
dt dt

n m dCC dCB
r&
2 = kr2 C C CD = =
dt dt

donde C representa la concentración de los reactivos y productos y kr las denominadas

constantes de reacción, (o constantes de velocidad o constantes específicas de
reacción), y n y m las ordenes de reacción respecto a A, B, C o D.

Se debe tener en cuenta que los valores de kr son en condiciones muy alejadas del
equilibrio, puesto que con excepción de las reacciones independientes de las
concentraciones de los reactivos, la velocidad depende de éstas.

Para sistemas ideales, la actividad de los reactivos y productos son iguales a sus
concentraciones. Para sistemas inmiscibles o sólidos puros, por lo general también se
considera que C=a=1, y para sistemas en que a ≠ C , aA = p A poA = g ANA en que NA
es la fracción molar de la especia A y γ A su coeficiente de actividad, inferior o mayor a
1.
Suponiendo las órdenes de reacción iguales a 1 para simplificar, en el equilibrio, r1 = r2
y se puede colocar entonces

kr1 C'A C'´B = kr2 C'CCD'

en que C’ denotan las concentraciones en la condición de equilibrio, distintas a las

anteriores. De esta forma se tiene que

C'CCC' kr
' '
= 1 = Kc (1.2)
C A CB kr2

Kc es la Constante de Equilibrio de la reacción (1.1) a la temperatura especificada, en

términos de las concentraciones. Si los reactivos y productos son gases, la expresión
toma la forma:

5
 p'CpD'
KP = (1.3)
p'A pB'

Donde p’ son las presiones de equilibrio respectivas.

Para una reacción en la que participan más de 1 mol de cada compuesto, como por
ejemplo:

aA + bB ƒ cC + dD
1 dC A 1 dCC
r&= - = .... = = etc (1.4)
a dt c dt

y la constante de equilibrio en términos de las concentraciones será:

c d
(C'C ) (C'D )
KC = a b
(C'A ) (C'B )

En términos de las presiones parciales, considerando que para un mol de

A se tiene A A = 1 / V y pA V = RT , se puede colocar entonces:

Kp =
(p ) (p ) = (C RT ) (C RT )
' c
c
' d
D
'
C
C '
D

(p ) (p ) (C RT ) (C RT )
' a
A B
' b '
A
a '
B
b

Kp = (RT ) =
(C ) (C )
∆n
' C
C
' d
D

(C ) (C ) ' a
A
' b
B

en que D n = (c + d) - (a + b), es decir,

Kp = Kc (RT)D n (1.5)

Cuando no hay variaciones en el número de moles que toman parte en la

transformación, D n = 0 y Kp = Kc.

En sistemas heterogéneos en que hay, por ejemplo, líquidos y gases como la oxidación
del sulfuro de fierro del eje de cobre a oxido ferroso a 1200ºC se tiene entonces que:

FeS(l) + 3
2 O2(g) € FeO(l) + SO2(g)

6
 aFeOp SO2
K=
aFeSp O2

o bien:

γ FeO NFeO p S O2
K=
γ FeO NFeO p 3O/ 2
2

Para una reacción en que participan Q componentes y pueden ocurrir R reacciones

químicas, la relación (1.1) se puede generalizar en la forma

αi1A1 + αi2 A 2 + ... + αiQ A Q = 0 (1.1)a

o bien,

N,R

å a ijA j = 0
j= 1

(1.1)b
donde j = 1,2...,R

En este caso, los coeficientes estequeométricos a ij son positivos para los reactantes y
negativos para los productos.

Para la expresión general (1.1)a, la velocidad de reacción puede escribirse en función

de una propiedad extensiva, por ej. el volumen:

1 1 dnj
r&
i = (1.6)
V a ij dt

Para el caso en que se tiene

aA + bB ® cC + dD
1 dnA 1 d
r&= - = (C A V)
aV dt aV dt

o bien,

1æ ççV dC A + C A dV ÷
ö
r&= - ÷
÷ (1.7)
aV çè dt dt ø

7
CA representa la concentración de A en (moles/unidades vol.). Si la densidad no varía,
el volumen permanece constante y la relación (1.7) se reduce a:

1 dC A
r& = −
a dt

Se puede ver que en este caso basta con medir la variación de concentración de A con
el tiempo para obtener la velocidad de reacción.
La forma particular de la expresión de velocidad varía según se exprese en función de
las concentraciones, volumen del reactor o moles participantes. De esta forma, para la
reacción general (1.4) y un reactor de volumen V donde participan nA, nB, etc. moles de
reactantes y productos, la expresión cinética toma las forma:

1 dnC dCC 1 dC A 1 dCB

r&= = =- =-
V dt dt a dt b dt

Otras formas pueden también ser definidas, según la conveniencia de cada caso. Así,
se puede expresar también la velocidad de reacción en la formas siguientes:

1 dnC moles de C formados

r'
&= =
W dt (masa de sólido total)(tiempo)

1 dnC moles de C formados

r''
&= =
S A dt (superficie del sólido A )(tiempo)

1 dnC moles de C formados

r'''
&= =
VSA dt (volumen del sólido A )(tiempo)

En estos casos, las velocidades así calculadas quedan relacionadas entre sí en la forma:

rV
& = r'W
& = r''S
& A = r'W
& = r'''
&S
A
(1.8)

donde W es la masa total de sólido en el sistema, SA la superficie o la superficie

específica del reactante y v sA su volumen o volumen específico.

Se debe notar que cuando en un sistema abierto se forma hay un producto gaseoso,
no se alcanza el equilibrio y la reacción ocurre hasta la desaparición de los reactivos.
Este caso es relativamente frecuente en metalurgia. Por ejemplo, en la reducción de
ZnO (sólido) para producir Zn gaseoso a 900ºC según:

8
ZnO(s) + CO(g) ƒ Zn(g) + CO2(g)

La variación de energía libre de Gibbs en condiciones distintas a las de equilibrio en

una reacción, se expresa mediante la llamada Isoterma de Van’t-Hoff:

ép ·p ù
D GT = - RTlnKP + RTln ê C D ú (1.9)a
êp A ·pB ú
ë û
o bien

éC ·C ù
D GT = - RTlnKC + RTln ê C D ú (1.9) b
êC A ·CB ú
ë û

Si las presiones parciales de cada gas son las normales (1 atm), o la fase condensada se
comporta en forma ideal de tal manera que a = C = 1, se tiene que

D G0T = - RTlnKP (1.10)a

D G0T = - RTlnKC (1.10) b

es decir, la isoterma de Van’t-Hoff se puede escribir

p p 
∆G T = ∆GoT + RT ln c D  (1.10)c
pA pA 

Problema

Calcular a qué presión de CO se debe efectuar la reducción de la alúmina con carbono

a 1500ºC.

Para la reacción de reducción, según

Al2O 3( s ) + 3C ( s ) ⇔ 2Al(l ) + 3CO( g)

0
D G1500ºC = 36.09 kcal

9
Se ve que en condiciones estándar (pCO=1 atm y Al, Al2O3 y C puros) no es posible
obtener aluminio líquido. Como en condiciones distintas de las estándar se tiene que

ép3 ·a2 ù
0
D G1500ºC = D G1500ºC + RTln êê CO Al3 úú
êëaAl2O3 ·aC úû
Considerando que las fases condensadas no forman soluciones líquidas o sólidas (son
puras), a = C = 1, se tiene que:

0
D G1500ºC = D G1500ºC + RTlnp3CO

Se puede calcular entonces D Gº 1500ºC para distintas presiones de CO, como se ve en la

Fig. 1.3. Se observa que es posible producir aluminio solo cuando pCO<5x10-2

30

20

10

0
G 1500 ºC (Kcal)

-10

-20

-30
V

-40

-50

-60
10 2 10 1 10 -1 10 -2 10 -3 10 -4 10 -5 10 -6
PCO (Atm.)

Figura 1.3 Reducción de Al2O3 con carbono a 1500ºC.

1.3 CONSTANTE DE EQUILIBRIO Y CALOR DE REACCION

La constante de equilibrio está relacionada con el calor de reacción mediante la

relación llamada Isocora de Van’t Hoff. Diferenciando la expresión (1.10)b con
respecto a la temperatura a presión constante, se tiene:

é¶ (D G0 )ù
ê ú dlnKP
ê ¶ T ú = - RTlnKP - RT dT (1.11)
êë úûP

o bien, como D G0 = - RTlnKP

10
 é¶ (D G0 )ù
T êê ú = D G0 - RT2 dlnKP
ú
êë ¶ T úûP dT

Como se tiene que según la relación de Gibbs-Duhem:

é¶ (D G)ù
D G = D H - TD S = D H + T ê ú (1.12)
ê ¶T ú
ë û
y para los reactivos y productos en su estado estándar la energía libre estándar de
Gibbs vale:

é¶ (D G0 )ù
D G = D H + T êê
0 0 ú
ú
êë ¶ T úû

se puede colocar que

¶ (D G0 ) dlnKP
T = D G0 - D H0 = D G0 - RT 2
¶T dT

es decir,

dlnKP
D H0 = RT2
dT

y por lo tanto,

dlnKP D H0
= (1.13)
dT RT 2

Esta expresión correlaciona el efecto de la temperatura en la constante de equilibrio,

y permite calcular KP a cualquier temperatura si se conoce Kp y D Ho a una cierta
temperatura, para rangos moderados no mayores de 200ºC, puesto que D Ho es
también dependiente de la temperatura.

Para compuestos puros o soluciones en que uno de sus componentes tenga una alta
presión de vapor respecto a los otros, para el componente de alta presión de vapor A
se tiene

11
 A()l ƒ A(g)

y entonces Kp ≅ p oA , si la presión del otro (s) componente (s) es despreciable. En este

caso, la relación (1.13) toma la forma
dlnpA D H0 L
= 2
= V2 (1.14)
dT RT RT

donde Lv es el calor latente de vaporización del componente de alta presión de vapor.

Esta relación es la denominada ecuación de Clausius-Clapeyron, y permite calcular la
presión de vapor a distintas temperaturas, conociendo el valor de Lv a otra dada.

1.4 INFLUENCIA DE LA CONCENTRACION DE LOS REACTIVOS

Dos conceptos deben ser definidos en relación a la concentración de los reactantes:

- molecularidad (estequeometría)
- orden de reacción

Molecularidad representa el número de átomos o moléculas que toman parte en una

determinada reacción, en tanto que orden de reacción denota el valor numérico de la
potencia a que se elevan las concentraciones de los reactivos que determinan la cinética
del proceso, y es igual a la sumatoria de los ordenes individuales de los exponentes de
las concentraciones en la expresión que representa la velocidad de reacción
considerada. El orden de reacción debe ser determinado en forma experimentalmente,
y puede ser de orden cero, 1, 2, etc. o fraccionario.

El orden de reacción se puede considerar como el número mínimo de moléculas o

átomos que toman parte en una reacción total. En ciertos casos, el conocimiento de los
ordenes parciales y orden global de la reacción permite dilucidar el o los mecanismos
involucrados en la transformación.

Ambos términos, molecularidad y orden de reacción no son equivalentes y tampoco

son numéricamente iguales.

Por ejemplo, para la reacción hipotética:

A (g) + B(g) ƒ C(g)

12
es unimolecular con respecto a los reactivos. Si se supone que es de primer orden
respecto a A y orden ½ respecto a B, la expresión cinética para la velocidad de reacción
r& toma la forma:

r& = k r C A CB0.5

El orden global de reacción sería 1 + 0.5 = 1.5

1.4.1 Reacciones de orden cero

En este caso, la velocidad de reacción es independiente de las concentraciones de los

reactantes y para un sistema transciente es constante en el tiempo.

Para una reacción como:

A+ Bƒ C

Se tendrá

dC A dCB dCC
r&= - =- = = kr (1.15)
dt dt dt

La pendiente de la recta que se obtiene al graficar la variación de la concentración de A

o B con el tiempo es cero, es decir, es una línea paralela a las absisas del tiempo.

En metalurgia existen casos de reacciones de orden cero, como la oxidación del

carburo de titanio y carburo de tungsteno; oxidación del sodio y potasio y ciertas
reacciones de disolución entre metales a 485-530ºC, como el Fe-Zn y Ti-Al. Para la
oxidación del carburo de titanio la oxidación es de orden cero respecto del oxigeno y
orden 1 respecto del TiC:

TiC(s) + 2O2(g) ƒ TiO2(s) + CO2(g)

dC TiC 1 dCO2 dC TiO2 dCCO2
- =- = = = krC TiC (1.16)
dt 2 dt dt dt

En ciertos casos, existen reacciones de seudo-orden cero que aparentemente son

independientes de la concentración (y el tiempo), como en la oxidación de ciertos
metales y aleaciones.

1.4.2 Reacciones de Primer Orden

13
Este caso en metalurgia es frecuente, lo que hace que el análisis de los sistemas
reaccionantes sea relativamente simple, ya que la velocidad de reacción es
directamente proporcional a la concentración del reactivo o los reactivos. Por ejemplo,
en la descomposición del carbonilo de níquel gaseoso según:

Ni(CO)4(g) ƒ Ni(s) + 4CO(g)

El orden de reacción es 1 con respecto a la concentración del carbonilo, por lo tanto

dCNi(CO) dCNi
- 4
= = kr CNi(CO)
dt dt 4

Reacciones de primer orden son frecuentes entre un sólido y un fluido, donde el orden
es cero o cercano a cero respecto del sólido y orden uno respecto al fluido.

Para una reacción tipo

A (s) + B(f) ƒ C(s)

en que sea de orden cero respecto al sólido A y de orden 1 respecto al gas B, se tiene
que

dCC dC A dCB
=- = = kr CB
dt dt dt

en que CB es la concentración de B en el tiempo t. Para la concentración inicial CBo y

otra cualquiera CB. Integrando esta última expresión:

t=0 t=t
CB = CB0 CB = CB

CB0
- ln = - kr t (1.17)a
CB

1 C0
kr = ln B (1.18)a
t CB

o bién

14
CB = CB0e- kr t (1.18)b

Si la concentración se expresa en moles, la expresión (1.17) a toma la forma:

NB
ln = - kr t (1.17)b
NB0
De esta forma, del gráfico de ln(NB / NB0 ) versus t (ó ln(CB / CB0 ) versus t) se obtiene una
línea recta cuya pendiente vale -kr

Reacciones de primer orden en metalurgia son por ejemplo la oxidación de sulfuros de

metales a óxidos ó sulfatos, las que son de orden cercano o igual a cero respecto al
sólido y orden 1 respecto al oxigeno, como ocurre por ejemplo con las reacciones:

2CuS(s) + 3O2(g) € 2CuO(s) + 2SO2(g)

1 dCO2 1 dCCuS
- =- = kr CO2
3 dt 2 dt

CuS(s) + 2O2(g) € CuSO4(s)

1 dCO2 dCCuS
- =- = kr CO2
2 dt dt

Otro ejemplo lo representa la transferencia de azufre disuelto en el metal (arrabio) a

1600ºC a la escoria en la producción de hierro:

S(l)m ⇔ S (l)e

dC(S) dC(S)
e
=- m
= krC(S)
dt dt m

El tiempo necesario para reducir la concentración de los reactivos a la mitad se puede

deducir de la expresión (1.18)a colocando para CB = 0.5CB0 es decir, la expresión queda

1 éê CB0 ùú
kr = ln
t0.5 êë0.5CB0 úû

es decir,

15
 ln2
t0.5 = (1.19)
kr

1.4.3 Reacciones de segundo orden

En este caso, la velocidad depende de dos concentraciones que pueden ser el mismo
reactivo o dos diferentes:

(a) 2A ® C
(b) A + B → C

La reacción puede ser, por ejemplo, de segundo orden con respecto a A en el caso (a),
y de segundo orden global y de orden 1 con respecto a A y orden 1 respecto a B en el
caso (b), es decir, se puede escribir:

1 dC A dCC
- = = krC2A
2 dt dt

dC A dCB dCC
- =- = = k r C A CB
dt dt dt

Llamando CºA = CBº las concentraciones iniciales en el caso (b) y X e Y la cantidad de A y

B reaccionada después del tiempo t (si esta es mol a mol), se tiene que la velocidad de
reacción queda dada por

dY dY
dX d(C0A - X) d(CB0 - dX
Y) dX 0
== - = - =- =- = = kr (C=A - X)(C- B0 - Y) - (1.20)
dt dt dt dt dt dt dt dt

es decir,

dY dX
= 0 = k dt
(C - X)(CB - Y) (CA - X)(CB0 - Y) r
0
A
0

Esta expresión se puede resolver con las CB siguientes:

I) t=0 C A = C 0A CB = CB0
II) t=t C A = C0A - X CB = CB0 - Y

Luego, se puede escribir

16
 1 éC0 (C0 - X)ù
ln êê 0 0 ú
B A
kr = 0 0 ú (1.21)
t(C A - CB ) êC A (CB - Y)ú
ë û

Si las concentraciones iniciales son iguales C0A = CB0 y la reacción ocurre mol a mol
entre A y B, CC0A0A -- xXx== CCB0B0-- YYy la expresión (1.20) se puede colocar en la forma:
dY dX 2 2
= = (C 0A - X) = kr (CB0 - Y) (1.22)
dt dt

Luego,

1 X 1 Y
kr = = (1.23)
t C A (C A - X) t CB (CB0 - Y)
0 0 0

Cuando X = 0.5, C0A = 0.5CB0 y el tiempo de reacción para alcanzar el 50% de conversión
del reactivo A o B será:

1 1
kr = 0
= 0
tC A tCB

o bien

1 1
t0.5 = = (1.24)
krC A kr CB0
0

Si la reacción procede según

2A ® B

dX 2
= kr (C 0A - X)
dt

y valen igualmente las expresiones (1.20) a (1.24).

Existen también reacciones de orden fraccionario inferior a 1, y de orden tres y

superiores, pero por lo general ocurren entre gases, aún cuando algunas reacciones
sólido-líquido pueden ser tercer orden.

17
1.5 DETERMINACION DEL ORDEN DE REACCION

Los métodos para determinar el orden de una reacción son los métodos integral, de
vida media, diferencial y de aislamiento.

1.5.1 Método Integral

Este método es muy útil en el caso de reacciones simples y consiste en ajustar la

información experimental a ecuaciones de distinto orden hasta encontrar la que mejor
representa la transformación, lo que permite entonces calcular el orden de reacción.
Para una reacción cualquier de la forma:
A® C

la expresión general tiene la forma

d(C0A - X) dX n
- = = kr (C0A - X) (1.25)
dt dt

donde n es el orden de reacción (0, 1, etc….) y dX/dt la disminución de concentración

con el tiempo.

o bien

dNA dNC
- = = kr CnA
dt dt

Una forma de presentar los datos experimentales es graficando log (C 0A / C A ) en

función de t o bien log (N0A / N A ). Si el orden es cero, la velocidad de reacción será

independiente de t. Si es orden 1, será una recta. Una curva si es 2, etc. (Fig. 1.4).

18
 1

n=2
=
n

log (CºA / CA )
n=0
n = 0.5

0
t
Figura 1.4 Formas de las curvas experimentales para diferentes ordenes de
reacción.
Para reacciones de primer orden, la pendiente vale entonces kr/2.303.

Este método en ocasiones es difícil de aplicar si los puntos experimentales están muy
dispersos y la forma de la curva no es simple de visualizar. El método tampoco se
presta cuando los órdenes parciales son fraccionarios.

1.5.2 Método de Vida Media

Se vio que el tiempo necesario para alcanzar la mitad de la concentración inicial de los
reactivos quedaba dado por la expresión (1.17) para reacciones de orden 1:

ln2 0.69
t0.5 = = (1.26)
kr kr

Para reacciones de 2º orden, según (1.21) o (1.23) cuando las concentraciones iniciales
son iguales y la reacción procede mol a mol o es un solo reactivo inicial:
1 1
t0.5 = =
krC A kr CB0
0

Es decir,

log t0.5 = cte. (reacciones de 1er orden)

log t0.5 = cte. - logC A0 (reacciones de 2º orden)

y para una reacción de orden n será:

19
log t0.5 = cte - (n - 1)logC A0 (1.27)

Si se conocen las concentraciones para diferentes condiciones iniciales se puede

conocer el orden de reacción al graficar log t 0.5 en función de log CºA , puesto que la
pendiente da el valor de –n.

1.5.3 Método Diferencial

Este método, conocido también como método diferencial de Van’t Hoff, es apropiado
para reacciones simples, donde la velocidad de reacción se puede expresar por

r&= kr CnA

o bien:
log r&= log kr + n logC A = cte. + n logC A (1.28)

Si se conoce la velocidad de reacción a distintas concentraciones, el gráfico de log r&

versus log C A dará una línea recta y la pendiente será el orden de reacción respecto al
reactivo A.

Se pueden emplear dos métodos para calcular la velocidad de reacción: uno de ellos
consiste en graficar la concentración del reactivo (s) en función del tiempo, de manera
que las pendientes de las rectas o curvas obtenidas representan el “orden de reacción
respecto al tiempo”. El otro método consiste en graficar el logaritmo de la
concentración de los reactivos en función del logaritmo de la velocidad de reacción
para diferentes concentraciones iniciales al comienzo de la reacción (es decir, la
velocidad de reacción calculada de la tangente de las curvas obtenidas para cada
concentración). Este representa el “orden con respecto a la concentración” o el “orden
real” de la reacción.

Cuando existe más de un reactivo, se efectúan determinaciones variando la

concentración de sólo uno de ellos a la vez; el orden global de reacción es la suma de
las ordenes parciales determinadas.

Ejemplo.
Si se tiene la reacción simple homogénea

A+ B€ C

y se obtuvo la siguiente información:

20
 Test # Conc. Inicial (moles/kg) Vida media (hrs)
A B
1 2 2 2.5
2 5 5 1.0

Calcular el orden de reacción.

En cada test, la vida media será (t0.5 )1 y (t0.5 )2 a las concentraciones iniciales de (C 0A )
1

y (C 0A ) . De la relación (1.28) se puede escribir entonces

2

(log t0.5 )1 = cte. - (n - 1)log (C0A )1

y
(log t0.5 )2 = cte. - (n - 1)log (C0A )2
es decir,

(log t0.5 )1 - (log t0.5 )2

n=
log (C0A ) - log (C0A )
2 1

(log t0.5 )1 - (log t0.5 )2

n= =2
log 5 - log2

Es decir, la reacción es de segundo orden.

1.5.4 Método de Aislación

Conocido también por método de Oswald, consiste en mantener uno de los reactivos
en exceso de manera que la concentración sea aproximadamente constante al
comienzo de la reacción o durante ella.

Para una reacción del tipo

aA + bB ® cC + dD

La variación de la concentración de las especies A y B será:

1 dC A 1 dCB
- - = kr CmA CnB
a dt b dt

21
Puesto que se mantiene constante la concentración de A en toda la experiencia, la
reacción es de pseudo orden cero respecto a A y puesto que CmA ; cte, es decir,
independiente de C A que se puede escribir:

dCB
- = k'r CnB
dt

De esta forma, variando CB se puede obtener el orden n trazando el logaritmo de la

velocidad inicial de reacción versus log CB en la forma

dCB
- log = logk'r ++nC
nlog
B CB (1.29)
dt

La pendiente de la recta será entonces el orden n con respecto a B, y la intercepción

con el eje de las ordenadas log k'r . Igualmente se procede con el componente A,
empleando un gran exceso de B.

Ejemplo:

Se estudió la disolución de una aleación aluminio-titanio, sumergiendo una lámina de

la aleación en aluminio líquido a 800ºC. Las curvas que se obtuvieron se observan en
la Fig. 1.5.

12
b
a

c
d

11
e

10
Disminución espesor, (mm)

9
8
7
6
5
a : Ti 1,0000 (Frac. atómica)
4 b : Ti 0,9535
c : Ti 0,9386
3 d : Ti 0,9014
2 e : Ti 0,8850

1
0 1 2 3 4 5 6 7 8 9 10
Tiempo, (hrs)
Figura 1.5 Resultados de la disolución de una placa de aleación de aluminio-titanio
en aluminio líquido a 800ºC.

22
Como en todos los casos se obtienen rectas experimentales, la velocidad de disolución
se puede calcular entonces directamente de la pendiente de las rectas, como se indica
en la Tabla 1.1.

Tabla 1.1 Valores deducidos de la Fig. 1.4

% atómico de Log % r&, velocidad disol. (*)
log r&
Ti en aleación atómico Ti (cm/seg)
100.0 2.0000 4.201x10-6 -5.377
95.35 1.9793 3.527x10-6 -5.453
93.86 1.9725 2.944x10-6 -5.531
90.14 1.9549 2.500x10-6 -5.602
88.50 1.9420 2.222x10-6 -5.653
( )
* Velocidad lineal, por cada cara de la plancha de aleación.

Se puede calcular ahora la contribución del aluminio al orden total de reacción. La

velocidad de reacción (disolución) es proporcional a las concentraciones de Ti y Al en la
aleación, es decir,
dC Ti dC
= = kr CnTiCmAl (1.30)
dt dt

en que C Ti y C Al son las concentraciones atómicas del Ti y A y dCTi/dt la velocidad de

disolución de la aleación Ti – Al. Como la disolución de la aleación es pequeña y la
variación de Al en la solución líquida (Al + Ti) – Al se mantiene virtualmente constante
con el tiempo, CmAl ; cte y entonces la expresión (1.30) se reduce a

dC Ti
= krCnTi
dt

o bien,

ædC ö
log çç Ti ÷ = logkr + nlogC Ti (1.31)
çè dt ÷
÷
ø

De esta forma, graficando log (dCTi / dt ) en función de log CTi (el componente aislado
en este caso), la contribución del titanio al orden global de reacción de disolución será
la pendiente de la recta, como se observa en la Fig. 1.6.

23
 -5,7

-5,6
log (dC/dt)

-5,5

-5,4

-5,3
1,90 1,92 1,94 1,96 1,98 2,00
log (% atómico Ti)

Figura 1.6. Determinación del orden de reacción respecto al titanio

De la Fig. 1.6 se encuentra que:

- 5.725 - (5.376)
n= ; 4.98 ; 5
1.93 - 2.00
Es decir, en un sistema diluido en titanio en que la concentración del aluminio es
aproximadamente constante, a 800ºC la velocidad de reacción se puede expresar en la
forma:

dC Ti
= krC 5Ti
dt

En sistemas heterogéneos que incluyen la desaparición de una fase o formación de

otra, a menudo hay varias etapas consecutivas que pueden ocurrir:

- Transporte de las moléculas desde el seno del fluido hacia la superficie reaccional
(que puede ser a través del film que rodea la partícula o cuerpo reaccionante).
- Absorción del reaccionante en la superficie.
- Reacción química en la interfase o volumen reaccional
- Nucleación de la nueva fase
- Transporte del producto desde la interfase hacia el seno del fluido, a través de un
film que rodea las partículas (si éste es un fluido).

La transformación completa será controlada por la etapa más lenta entre ella, o la
combinación de etapas más lentas, si hay dos o más de magnitudes similares.

24
1.6 INFLUENCIA DE LA TEMPERATURA EN LA VELOCIDAD DE REACCION

La temperatura ejerce una fuerte influencia en la velocidad de reacción, casi sin

excepción. Para reacción homogéneas y heterogéneas, se ha encontrado que la
velocidad de reacción se duplica aproximadamente por cada 10 grados de aumento de
temperatura. La relación puede expresarse mediante la relación fenomenológica:

E
-
kr = k 0 e RT
(1.32)

Es expresión se conoce por ecuación de Arrhenius, y en ella kr es la constante

específica de reacción, ko el llamado factor de frecuencia (o factor de colisiones entre
átomos o moléculas) y E la energía de activación de la reacción. La ecuación de
Arrhenius se puede aplicar tanto a reacciones homogéneas como heterogéneas. Se ve
que tanto kr como ko tienen iguales unidades de medida.

Esta ecuación puede ser deducida también de la isocora de Van’t Hoff y la ley de
acción de masas. La razón de ello está en que puesto que las fuerzas intermoleculares
actúan sólo en un rango muy pequeño de espacio (corta distancia), sólo las moléculas
más cercanas capaces de entrar en colisión pueden reaccionar, de forma que la
colisión es un requisito para que ocurra la reacción. Sin embargo, no todas las
moléculas pueden reaccionar ya que deben además se activas.

Este concepto de moléculas activas y energía de activación se puede explicar mejor

mediante la ley estadística de distribución de Maxwell. De esta forma, el exponencial
exp(E-/RT) se puede interpretar como la fracción de moléculas o átomos que tienen
una energía igual a o mayor a E. Físicamente, esto se puede representar como una
cierta barrera de dimensiones E que los reactivos A y B, con cierto valor inicial de
energía, deben sobrepasar para llegar a un nivel inferior de energía, en la forma de
productos C y D, como se observa en la Fig. 1.7.

25
 Energía Total del sistema
E

A B
 H,  G ó  F
C D

Figura 1.7 Representación gráfica del concepto de energía de activación.

Se puede ver en la Fig. 1.7 que, a menos que las moléculas o átomos de A y B posean
en la colisión una energía igual o superior a E, no podrán superar esta barrera
energética y llegar al estado C y D. Se observa además que la diferencia entre los
estados inicial y final, permite conocer D H, D G ó D F , pero no E.

1.6.1 Determinación de la energía de activación

Es posible deducir en forma teórica el valor de la energía de activación, sin embargo en

reacciones simples se emplea por lo general el método experimental. Para ello, se
determinan las velocidades iniciales de reacción (a concentración aproximadamente
constante) para diferentes temperaturas a partir de determinaciones de velocidad de
reacción en función del tiempo para las temperaturas consideradas y se calcula la
constante específica de velocidad (constante de reacción) en cada caso. De la ecuación
de Arrhenius se tiene que

E
lnkr = lnk 0 - (1.33)
RT

De esta forma, conociendo k r a diferentes temperaturas, al graficar ln kr versus l/T, la

pendiente de la recta obtenida representa la energía de activación E/R.

Ejemplo:

El caso visto de disolución de la aleación Ti-Al en aluminio fundido, se hicieron

determinaciones de la velocidad de disolución a distintas temperaturas (Fig. 1.8). En
este caso la velocidad de reacción es constante con el tiempo y por lo tanto la
constante de reacción se puede calcular de las rectas obtenidas. Así se tiene que:

26
 5

Disminución total espesor, (mm)

4

C
0º
90
ºC
0
3 85
C
0º
80
2
ºC
750

1
700 ºC

0
0 1 2 3 4 5 6 7 8 9 10
Tiempo, (hrs)
Figura 1.8 Resultados experimentales de disolución de la aleación Ti.Al en aluminio
líquido.

Tabla 1.3 Valores de la constante de velocidad

Ti ºK ln kr 1/T (ºK-1)
973 -14.228 10.280x10-4
1023 -14.773 9.775x10-4
1073 -15.243 9.320x10-4
1123 -15.630 8.905x10-4
1173 -16.052 8.525x10-4

Esto se ha graficado en la Fig. 1.9

-14
In k r

-15

-16

8,2 8,4 8,6 8,8 9,0 9,2 9,4 9,6 9,8 10,0 10,2 10,4
I / T (ºK -1)

Figura 1.9 Expresión de Arrhenius para la disolución de la aleación Al – Ti en aluminio

27
De esta forma, como la pendiente de la recta vale E/R se tiene que

é - 14.228 - (- 16.052) ù
E = 1.9875´ ê ú
ê10.24 ´ 10- 4 - 8.52´ 10- 4 ú
ë û

cal kcal
E ; 21077 ; 21.1
mol mol

Por otra parte, la constante de frecuencia (o el factor de frecuencia) ko vale:

k
k0 = E
e RT

6.62´ 10- 7
k0 = - 2107
e 1.987·973

æcm ö
k 0 = 3.583´ 10- 2 çç ÷ ÷
çèseg ÷
ø

1.7 TEORIAS DE REACCION

Existen por lo menos dos teorías de importancia en cinética de reacciones: la teoría de

colisión y la teoría de reacción absoluta o de estados de transición.

1.7.1 Teoría de Colisión

Esta teoría se adapta bien a los gases y con restricciones a los líquidos debido a la
relativa movilidad de sus átomos y moléculas, y considera que para que se efectúe la
reacción entre los componentes A y B, debe ocurrir una colisión entre éstas. La
velocidad de reacción en tal caso queda dada por la expresión derivada de la ley de
distribución de energía de Maxwell-Boltzmann.

De la teoría cinética de los gases se tiene que el número de colisiones entre A y B por
unidad de volumen y tiempo está dado por la relación:

0.5
é æM + MB ö÷ù
Z AB = s n n êê4p kT çç A
2
÷ú (1.34)a
AB A B çè M ·M ø÷
÷úú
êë A B û

28
 0.5
é8p RT ù
2
o bien Z AB = s n n ê ú
AB A B (1.34)b
êë m úû

en que:
s AB Diámetro medio molecular de A y B, (cm)
nA y nB Concentración de A y B, (moléculas/cm3)
MA y MB Masa de las moléculas de A y B, (gr/moléculas)
k Constante de Boltzmann = 1.3 x 10-16 erg/(ºK)
T Temperatura absoluta, ºK
m “Masa reducida” de las moléculas dado por m= éëMA ·MB (MA + MB )ùû

De la ley de distribución de energía de Maxwell-Boltzmann se tiene a su vez que la

fracción de moléculas dnu en el rango de velocidad entre u y (u + du) vale:

2
32 çæ Mu ö÷
÷
æ M ÷ ö ççèç- ÷
÷
RT ø÷ 2
dnu = 4p çç e u du
çè2p RT ÷
÷
ø

Si se consideran sólo colisiones directas coplanares, la fracción de moléculas con

energía entre E y E + dE es

æ Eö
1 ççèç- RT ÷÷÷ø
dnE = e dE
RT

y entre un cierto valor E = E y ¥ será

¥ ¥ æ Eö
1 ççèç- RT ø÷÷÷
ò E ò RT e dE = exp (−E / RT)
dn = (1.35)
nE E

De esta forma, combinando (1.34)b y (1.35) se tiene la velocidad de reacción en

términos del número de colisiones que producen una reacción química:

0.5
æ ö æçççè- RTE ÷÷÷öø
2 ç 8p RT ÷
r = nAnB s AB ç
& ÷ e
÷ (1.36)
çè m ÷ ø

Como la velocidad de reacción es proporcional al número de moléculas presentes, se

puede escribir

r&= krnAnB
29
en que kr está expresado en (m3/seg), por ejemplo. De esta forma, se puede colocar

0.5
æ8p RT ö÷ æçççè- RTE ÷÷÷öø
kr = s çç
2
÷ e (1.37)
çè m ø÷
AB
÷

Si se compara esta expresión con la de Arrhenius, se puede observar que, para

Arrhenius:

æ E÷ ö
ççç- ÷
è RT ø÷
kr = k 0 e

se puede ver que

0.5
æ8p RT ö÷
k 0 = s çç
2
AB ÷÷
çè m ø÷

En el caso de reacciones entre las mismas moléculas, tal como 2A ® C , la expresión

(1.34)b toma la forma

0.5
2 2æ 8p RT ö÷
Z AA = s AA nA çç (1.34)c
çè MA ø÷ ÷

Además de ello, otro factor se debe considerar para que la velocidad predicha sea igual
a la experimental, y se refiere a la orientación que deben tener las moléculas para que
se produzca la reacción, es decir, las moléculas (o átomos) para que efectivamente
reaccionen deben tener una cantidad mínima de energía igual o mayor que un cierto
valor E (energía de activación) y poseer una orientación preferencial caracterizada por
un factor esteárico P, menor de 1. De esta forma, la velocidad real de reacción queda
dada por la reacción:

-12 æ Eö
æ æM + MB ö÷ö÷ ççç- ÷
÷
r&= nAnB s ççç4p kT çç A
2
÷÷ Pe è RT ø÷
(1.38)
èç
AB çè M ·M ÷
A
÷÷
÷
B øø

æ 0.5
E÷ ö
æ ö
2 ç8p RT ÷
ççç- ÷
è RT ø÷
o bien, r = nAnB s AB ç
& ÷ Pe
÷
çè m ø÷

1.7.2 Teoría de la velocidad absoluta

30
Conocida también como de los complejos activados o de los estados de transición, fue
desarrollada para evitar los problemas derivados de la teoría de colisión que puede
predecir velocidades de reacción por exceso o defecto.

En la teoría de velocidad absoluta, se considera que las reacciones ocurren igualmente

por colisión entre moléculas o átomos, los que forman compuestos intermedios
(complejos activados) que a su vez se descomponen continuamente para dar origen a
los productos finales de la reacción, es decir, para una reacción A + B ® C + D ,
ocurrirá según:

A + B ¾k*¾
® A * + B* ® C + D (1.39)

En este caso, K* representa la constante de equilibrio entre los reactivos y el (o los)

complejos activados. De esta forma, para estos últimos se puede escribir que:

D G* = - RTlnK * (1.40)

o bien

éC* ùéC* ù
ê A úê B ú
D G* = - RTln ë ûë û (1.41)
éC A ùéCB ù
ëê ûëúê ûú
es decir,
D G*
éC* ùéC* ù= éC ùéC ùe- RT
(1.42)
úê B ûú ëê A úê
ëê A ûë ûë B ûú

Esta última expresión permite entonces calcular la concentración de los complejos

activados a una cierta temperatura conociendo D G * y las concentraciones de A y B.
Gráficamente, se puede observar en la Fig.1.10.

En términos de superficies de energía potencial, para una reacción tipo

AB + C ® A + BC se puede observar los niveles mínimos y máximos de energía
potencial en la Fig. 1.11.

31
 A*+ B*

Energía total del sistema

E

A+B
V V V
H, G ó F
C+D

Coordenada de reacción
Figura 1.10 Representación de una reacción mediante la formación de complejos
activados.

A*+B*+C*

AB+C
dxz

dyz A+BC

dxy

Fig. 1.11. Representación esquemática de la energía potencial de un sistema de 3

átomos en configuración lineal en función de las distintas inter atómicas. La línea
segmentada representa la reacción AB + C ƒ A + BC .

Se observa en la expresión (1.42) que el miembro derecho de esta última expresión se

asemeja a la relación de Arrhenius.
La velocidad global de reacción en un instante cualquiera se considera igual al
producto de la concentración de las moléculas de los complejos activados A* y B* en
un nivel máximo de energía (E = D G*) por el factor de frecuencia para cruzar esta
barrera energética. Se puede demostrar que el factor de frecuencia vale kT h , en que
k es la constante de Boltzmann y h la constante de Planck, igual a 6.63 x 10-27 (erg x
seg).

La velocidad de reacción de formación de los complejos queda expresada entonces en

la forma:

kT kT
r&
A*,B* = [C A* ][CB* ]= k *[C A ][CB ] (1.43)
h h
32
Reemplazando (1.43) en (1.42), y puesto que para los complejos activados se tiene
que ∆G* = ∆H* −T∆S * , se llega a la expresión

æD S*ö æ D H* ö÷
kT ççç ÷ çç-
÷
è RT ø÷ èç RT ø÷
÷
r&
A*,B* = [C A ]
[CB ]e e (1.44)
h

Como ∆S * varía muy poco con T, se puede considerar constante y entonces:

æ D H* ÷
ö
ækT ö÷ ççç-
è RT ÷
÷
ç
rA*,B* = ç ÷ ø
& ÷C [C A ][CB ]e
çè h ø÷
(1.45)

Estas últimas dos expresiones muestran que la velocidad de reacción depende de una
constante universal, común a todas las reacciones; de la variación de entropía para la
formación de complejos activados y de la variación de entalpía para la formación de
estos compuestos, así como de la concentración de los reactantes. Se debe notar que
aún cuando ambas teorías permiten predecir las velocidades de reacción y las energías
involucradas, a menudo existe gran discrepancia entre los valores calculados y los
experimentales, lo que refleja la mayor complejidad de las reacciones sobre lo
expuesto, sin embargo, las teorías permiten predecir ordenes de magnitud, de utilidad
al predecir el comportamiento de un determinado sistema.

BIBLIOGRAFÍA CAPITULO I

- J.J. Moore, “Chemical Metallurgy”, Butterworth, (1981)

- J. Mackowiak, “Physical Chemistry for Metallurgists”, G. Allen and Unwind Ltda.,

(1965).

33
 CAPITULO II
TRANSFERENCIA DE MASA

2.0 INTRODUCCION

Transferencia de masa se refiere al transporte o transferencia de una o más especies

(átomos o moléculas) en un medio. Esta puede ser en un fluido (gas, líquido, metal o
escoria fundida) o bien en un sólido (difusión en estado sólido).

El transporte o transferencia de masa entre un fluido (gas o líquido) y otro fluido

puede ocurrir por tres mecanismos:

- difusión molecular
- microturbulencia ( eddy difusión)
- Convección (natural y forzada)

El transporte de masa entre un fluido (gas o líquido) y un sólido puede ocurrir por dos
mecanismos diferentes o combinación de ambos:

- difusión interpartícula
- transporte en el film de fluido que rodea la partícula (difusión molecular).

El transporte de masa de un sólido en otro, ocurre solo a través de mecanismos de

difusión en estado sólido y están regidos por las leyes de difusión de Fick.

Para evitar ambigüedades, en adelante se denominará transporte o transferencia de

masa a los fenómenos que ocurren en un fluído y una interfase, tal como transporte
de oxígeno desde el seno de un gas hasta la superficie de una partícula de carbón en
combustión y se llamará a difusión al fenómeno asociado al transporte de una (o más)
especies de un fluído en otro por difusión molecular; de un fluído hacia (o desde) el
interior de un sólido, o de un sólido en otro (difusión en estado sólido).

Para los fluídos, el mecanismo intrínseco de transferencia de materia está asociado

tanto al movimiento de los átomos o moléculas (difusión molecular) como a la
producida por un efecto dinámico debido a microturbulencia del fluido (difusión
convectiva) causada por efectos térmicos, de diferencia de densidad, etc. En la mayor
parte de los casos, solo una de éstas tiene un efecto importante y la otra puede ser
despreciada, sin embargo, como se verá más adelante aún en fluídos con un alto grado
de agitación, la difusión a través del film estático (capa límite de concentración) toma
parte en el fenómeno en su totalidad.

2.1 GENERALIDADES SOBRE CONCENTRACIONES Y FLUJOS

Para un componente A de un sistema multicomponente, se pueden definir las

concentraciones, velocidades y flujos de la siguiente forma:

34
DEFINICIONES EN TRANSFERENCIA DE MASA

a) Concentraciones

- Concentración molar:
C

/
M
A

A
= r
n
T

r = å r
i
1

- Fracción en peso:
W

/
A

= r å r= r r

para gases:
n AV
C

m
o
l
e
s
/
c
m
3
A

r
R
T

( )
A

= =

C T = concentración total (moles/cm3)

C A = concentración de A (moles/cm3)
r A = densidad de A (gr/cm3)
r T = densidad total (gr/cm3)
WA = fracción en peso de A (-)
nA = nº moles A (-)
V = volumen total (cm3)
pA = presión parcial de A (atm)
P = presión total (atm)
X A = fracción molar de A (-)
MA =peso molecular A (gr/mol)

- Concentración molar total:

C T = nT / V = P / RT

- Fracción molar:

X A = C A / C T (líquidos y sólidos)

X A = C A / C T = (pA / RT) / (P / RT )= p A / P (gases)

Para una mezcla de gases binarios A-B se cumple que

35
P p A pB
= +
RT RT RT

r A = C AMA
WA / MA X AMA
XA = y WA =
WA / MA + WB / MB X AMA + XBMB

b) Velocidades:

- Velocidad másica media:

En sistemas multicomponentes, cada especie se mueve con una cierta velocidad ui y la

u
mezcla lo hace con cierta velocidad másica media definida en la forma:

n n

å r iui å r iui
i= 1 i= 1 r AuA + r BuB
u= n
= =
rT rT
å r
i= 1

En que uA y uB son las velocidades másicas absolutas de las especies A y B respecto de

coordenadas fijas.
U

- Velocidad molar media es a su vez:

n n

å CiUi å CiUi
i= 1 i= 1 C AUA + CBUB
U= n
= =
CT CT
å Ci
i= 1

Donde Ui, UA, y UB son las velocidades molares absolutas y (UA - U) y (uA - u) son las
velocidades molares y másicas relativas de A y B, respectivamente.

c) Flujos

El flujo molar o másico de una especie es una cantidad vectorial que define la
cantidad de masa transferida por unidad de tiempo y área normal al vector y se
puede expresar como el producto del coeficiente de difusión molecular por la
gradiente de concentración.

La relación fundamental es la relación fenomenológica de Fick, que define la difusión

de una especie en condiciones isotérmicas e isobáricas. Si la difusión es
unidireccional en z para un sistema de dos componentes AB la expresión tiene la
forma:

36
&
JA z = - DAB (dC A / dz) (ML- 2 t- 1 ) (2.1)

o bien,

jA z = - C T DAB (dX A / dz) (ML- 2 t- 1 )

J

U
&A es el flujo molar relativo a la velocidad molar media para una gradiente (dCA/dz)
en z. DAB es un coeficiente de proporcionalidad o coeficiente de difusión molecular
para la especie A difundiendo en B con dimensiones L2 t-1.

El flujo másico de la especie A que difunde en dirección Z respecto a la velocidad

másica media será a su vez:

dWA dr dX
jA z = - r T DAB = - DAB A = - C TDAB A (2.2)
dz dz dz

Para un sistema binario en que la velocidad media molar es constante en la dirección z

el flujo molar se puede expresar en la forma:

(
jA z = C A UA z - U
z
) (M2 T- 1L- 2 ) (2.3)

Es decir,

dX A
C AUA z = - C T DAB + CA U z (2.4)
dz

Como se tiene que

C AUA z + CBUB z
U z=
CT

o bien,

C A U = X A (C AUA z + CBUB z )
z

Se puede escribir que:

dX A
C AUA z = - C T DAB + X A (C AUA z - CBUB z ) (2.5)
dz

37
Como UA Z
y UB Z son velocidades relativas a la coordenada fija z, las cantidades C AUA Z

y CBUB Z son flujos también molares relativos a estas coordenadas fijas, de manera que
se pueden definir un nuevo flujo molar en la forma:

N&A Z = C AUA Z

N&B Z = CBUB Z

o en general,

N&A = C AUA
N&B = CBUB

De esta forma, reemplazando en la expresión (2.5) se tiene que

dX
(
N&A z = - C T DAB A + X A N&A z + N&B z
dZ
) (2.6)

y en general para el flujo molar de una especie A en un sistema binario A-B:

N&A = - C T DABÑ X A + X A (N&A + N&B ) (2.7)

Se puede notar que el flujo molar N&A es el resultado de dos cantidades vectoriales:
- C T DABÑ X A que representa el flujo molar causado por la gradiente de concentración
y X A (N&A + N&B )= C A U que representa al flujo molar debido a la convección del fluido.

Si la especie A depende de un sistema multicomponente, la expresión (2.7) toma la

forma:

n
N&A = - C T DANÑ X A + X A å Ni (2.7)a
i= 1

Donde DAN es el coeficiente de difusión de la especie A en la mezcla multicomponente

con otros fluidos N.

Por otra parte se puede también definir un flujo másico n&A de las especie A con
respecto a coordenadas espaciales fijas para una sistema binario en forma similar a
N&A :

n&A = - r T DABÑ WA + WA (n&A + n&B ) (2.8)

38
En que n&A = r AUA

En condiciones isotérmicas e isobáricas, se tiene entonces que

n&A = - DA B Ñ r A + WA (n&A + n&B ) (2.9)

Las cuatro expresiones que denotan los flujos (jA &

JA , N&A y n&A ) son, expresiones
equivalentes de la ley de Fick y el coeficiente de difusión DAB es igual en todos los
casos. Cualquiera de ellas puede ser empleada, sin embargo, en fenómenos de
transferencia de masa en que ocurren reacciones químicas generalmente se emplean
j
y
N

los flujos molares &.

En la Tabla 2.1 se encuentran las formas equivalentes de las ecuaciones de flujo para
sistemas binarios.

Tabla 2.1. Formas equivalentes de las ecuaciones de flujo másico en sistemas

binarios.
Flujo Gradiente Ecuación de Fick Condiciones
n&A Ñ WA (*)
n&A = - r T DABÑ WA + WA (n&A + n&B ) r = cte.
Ñr A n&A = - DABÑ r A + WA (n&A + n&B )

N&A Ñ XA N&A = - C T DABÑ X A + X A (N&A + N&B ) CT=cte.

Ñ CA N&A = - DABÑ C A + X A (N&A + N&B )

jA Ñ WA jA = - r T DABÑ r A r = cte.
jA = - DABÑ r A CT=cte.
A

Ñr

JA Ñ XA JA = - C T DABÑ X A
Ñ CA JA = - DABÑ C A
(*) dWA dWA dWA
Ñ WA = + + ,etc
dx dy dz

El transporte de materia puede también ser producido por condiciones físicas

impuestas al sistema distintas a la vorticidad (microturbulencia o “eddy diffusion”) y a
la difusión molecular. Por ejemplo, en un sistema multicomponente, al aplicar un
gradiente de temperatura se produce la transferencia de masa desde un punto a otro
(efecto Soret). El efecto de la fuerza centrífuga también puede producir la
transferencia de masa, tal como se emplea para separar algunos isótopos gaseosos del
uranio. Otro tanto ocurre con el transporte de iones bajo un campo eléctrico
(electroobtención de metales) y bajo una gradiente de tensión superficial (efecto
Marangoni).

39
2.2 EL COEFICIENTE DE DIFUSION

El coeficiente de proporcionalidad DAB o coeficiente de difusión o difusividad se

puede obtener de la expresión (2.1) vista:

&
JA Z M l L2
DAB = = 2 =
dC A / dz L t M / L3 1L T

y tiene por consiguiente iguales dimensiones que otros coeficientes de

proporcionalidad como la viscosidad cinemática v y la difusividad térmica a . Este
coeficiente de difusión depende por lo tanto de la temperatura, presión y composición
del sistema. Su similitud con otros coeficientes de proporcionalidad es evidente:

du
Newton : τ = v (esfuerzo de corte)
dx

dT
Fourier : q = a (flujo de calor)
dx

dC
Fick : J = D (flujo de masa)
dx

Los valores de los coeficientes de difusión varían considerablemente de una especie a

otra y en ausencia de datos experimentales se emplean relaciones semi empíricas para
estimar sus valores.

El coeficiente de difusividad varía aún dentro de sistemas homogéneos. Por ej., para
gases fluctúa entre 0.3 a 13 (cm2/seg); para líquidos entre 10-5 a 10-3 (cm2/seg) y entre
10-2 a 10-5 (cm2/seg) para metales y escorias fundidos, en tanto que para sólidos varía
ampliamente entre 10-8 a 10-30 (cm2/seg), dependiendo de la temperatura y del
sistema en particular.

2.2.1 Difusión en gases

Para difusión en gases en un sistema binario no reaccional, se puede emplear la

relación de Andrusov, entre varias propuestas:

é T1.78 ùé1 + (M + M )0.5 ù

ê
- 2
DAB = 6.06 x 10 ê 0.33 úê A B ú (2.10)
0.33 úê 0.5 ú
P ˆ
êë ( A
v + ˆ
v B )úê
(M M
A B ) úû
ûë

En que T está en ºK y P en mmHg. v̂ son los volúmenes molares, que para gases de
interés en metalurgia tienen los valores indicados en la Tabla 2.1.

40
Tabla 2.1 Valores de los volúmenes molares de algunos gases
v̂ Aire Cl2 CO CO2 H2 H2O N2 O2 SO2
3
(cm /gr-mol) 29.9 48.4 30.7 34.0 14.3 18.9 31.2 25.6 44.8

Otra correlación que puede ser empleada es la de Chapmann – Enskog, válida para
gases monoatómicos:

0.5
1.8583x10- 3 T 3/2 æ 1 1 ÷ö
DAB = çç + ÷
Ps AB 2WAB çèMA MB ÷
÷
ø

Donde s AB = 0.5(s A + s B ) y es el diámetro de colisión entre las especies A y B. Ω es la

integral de colisión para A-B a la temperatura adimensional T* para el potencial de
Lennard – Jones T* = (KB/ε)ABT, donde KB es la constante de Boltzmann (1.38 x 10-16
erg-,molécula-1 ºK-1), ε un parámetro característico de la función potencial de Lennard-
22 22 ù 0.5
Jones e igual a su vez a éêéê(ee//KK)AA2Ag(ee//KK)BBù en que e / K = 1.15Te (ºK)
ëë ûúûú
y e / K = 1.92Tf (ºK)) válida para metales gaseosos (Te=ebullición, Tf=fusión). W se
puede calcular de la Figura 2.1

10.0
T* = (K/ ) T

7.5

5.0

2.5

0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5
AB
Figura 2.1 Variación de la integral de colisión con la temperatura reducida para el
potencial de Lennard-Jones.

2.2.2 Difusión en líquidos

Para líquidos no es posible aplicar la teoría cinética de los gases y las expresiones son
semi empíricas o empíricas. En el caso de electrolitos, esta puede ser de origen iónico
o molecular y representa un campo separado de análisis en los Cursos de
Electroquímica. Para el caso de no electrolitos, existen dos teorías para explicar el
transporte de masa en líquidos. Una es la Teoría de Eyring y la otra es la Teoría
Hidrodinámica.

41
En la Teoría de Eyring, el fluido se considera ideal con estructura cristalina y una cierta
distribución de vacancias El transporte (difusión) de una especie se interpreta
entonces como un proceso unimolecular en el cual hay “saltos” de la especie que es
transportado dentro de estas vacancias.
En la teoría hidrodinámica se postula que el coeficiente de difusión en el líquido está
relacionado con su movilidad, esto es, la velocidad neta de ésta bajo la influencia de
una fuerza impulsora unitaria. Una ecuación derivada de esta teoría es la Stokes-
Einstein:

KT (2.11)
DAB =
6 p rA mB

Donde DAB es el coeficiente de difusión de A en una solución de diluida de B, K la

constante de Boltzman, T la temperatura absoluta, rA el radio atómico o molecular de
la especie A y µB la viscosidad de la especie B (solvente).

Una relación similar a la anterior es la de Sutherland:

1
ö3 (2.12)
K TæçNa ÷
DAB = ç ÷
2p çè vˆ A ø÷
÷

En que Na es el Nº de Avogadro (60.23x1023 moléculas/gr-mol) y v̂ A el volumen molar

de A.

Existen otras correlaciones basadas en diferentes teorías, todas las cuales predicen el
coeficiente de difusión con un cierto grado de precisión, dependiendo de las
características de la especie que difunde y el medio. (ver por ej., R.C. Reid y T.K.
Sherwood, “The Properties of Gases and Liquids”, Mc Graw Hill, NY, 1958).

Es importante notar que, independiente del sistema o especies que difunden, los
valores del coeficientes de difusión en líquidos fluctúan entre 10-2 y 10-5 (cm2/seg) con
energías de activación del orden de 1-4 Kcal/mol y un comportamiento similar al
predicho por la relación de Arrhenius en función de la temperatura.

Algunos valores se encuentran en la Tabla 2.1 y fig.2.2

42
 10-3

10 -4
Sb
Cd en B C
en i Ce en
P n 1 2 .5
Bi en b B . 5% %
D, (cm2/seg)

Sb ie nP C C
b

D, (cm2 /seg)
Si
Sn
en 10-4 P
P b Mn
S Ti
S en
0.6 %
10 -5 S
Ni
1.0 1.1 1.2 1.3 1.4 1.5
10 3/ T, (ºK-1 ) S
en
0.9
%
S
10-5
Fe líq. saturado en C .50 .55 .60 .65 .70
Fe liq. puro 3 -1
10 / T, (ºK )

Figura 2.2 Valores del coeficientes de difusión de metales líquidos.

Tabla 2.2 Valores del coeficiente de difusión para algunos metales líquidos en
escorias fundidas y en otros metales líquidos.
Sistema Temperatura D
ºK (cm /seg) x 10-5
2

Ca en 40 CaO – 20 Al2O3 – 40 SiO2 1350 –1450 0.067

Si en 40 CaO-20 Al2O3 – 40 SiO2 1365 – 1460 0.010
0 en 40 CaO-20 Al2O3 – 40 SiO2 1370 – 1520 0.40
Fe en 43 CaO-22 Al2O3 – 35SiO2 1500 0.35
Fe en 61 FeO-39 SiO2 1250 – 1305 9.6
P en 40 CaO-21 Al2O3-39SiO2 1300 – 1500 0.2
Fe en Fe -29S 1958 – 1254 2020
Cu en Cu-19.8 S 11.60 – 1256 69.3
Cu en 32 Fe-Cu 1160 – 1245 56.4
O2 en H2O 25 2.5

2.2.3 Difusión en Sólidos

Se deben distinguir dos clases de difusión en el interior de un sólido:

- difusión de un fluido (gas o líquido) en los poros hacia (o desde) el interior de

un sólido poroso.
- Difusión de un sólido en otro.

a) Difusión de un sólido en otro (difusión superficial)

La difusión de un sólido en otro sigue las leyes de Fick y su estudio con aplicación a la
metalurgia física se ve en detalle en los cursos de Metalurgia Física.

43
La difusión de un sólido en otro o difusión superficial se produce de acuerdo a cuatro
mecanismos independientes: difusión por vacancia; difusión intersticial; reemplazo y
por e intercambio directo (Fig. 2.2).

Vacancia Intersticial

Reemplazo Intercambio directo

Figura 2.3 Mecanismos de difusión en un sólido cristalino.

Para los sólidos, la temperatura tiene considerable influencia en los coeficientes de

difusión, encontrándose experimentalmente que ésta sigue una correlación
fenomenológica tipo Arrhenius:

-
E
(2.13)
D = D0 ×e RT

Donde D0 es un equivalente del factor de frecuencia visto anteriormente para

reacciones químicas. Valores para el coeficiente de difusión de algunos metales en
otros metales se encuentran en la Fig. 2.4

44
 10-9

C e % M .5
Ce
co n e co n
nF

n Fe n
1 2
10-9

%
Bi

Cr e
en C e
N
en

F n Fe

Be
Fe

nF
e
10-9
D, (cm 2 /seg)

10-9

Zn
20

en
5. 6
%

Al
15

Cu Cu
Zn
%

35
%

Sn
-8

e
%
Al

n
0%
-9
-8

10-9
4 .4
5%
10 % Cu

-6
Cu
90
%

5%
Cu

en n Cu
Ni e n C

Cu

Al
en

Cu
-

e
Cu

10-10
u

0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3

103 / T, (ºK -1 )

Figura 2.4 Coeficientes de difusión de algunos metales y aleaciones en otros

metales.

b) Difusión de un fluído en un sólido poroso

La difusión de un fluido (gas o líquido) hacia (o desde) el interior de un sólido puede

ocurrir debido a la porosidad del sólido (difusión en poros) o por difusión en estado
sólido (superficial) o difusión de Fick. Normalmente esta última es despreciable frente
a la difusión en poros, puesto que para la mayor parte de los sólidos (con excepción de
metales, aleaciones y algunos compuestos) estos tienen una cierta porosidad macro o
microscópica.

En la difusión de un fluido en los poros de un sólido, ésta puede ocurrir mediante la

difusión de Fick, difusión de Knudsen o superficial. Si los poros del sólido son
relativamente grandes y el gas es denso, ésta ocurre según las leyes de Fick. Si los
poros son pequeños e inferiores al camino medio libre recorrido por las moléculas que
difunden, esta ocurre según la difusión de Knudsen, la que se puede expresar por la
relación

DAB- K = 1.33 (8 R T / p MA )0.5 K0 (2.14)

45
K0 un parámetro característico del sólido y vale Ko=ε dap /4), donde ε es la porosidad
macroscópica del sólido y dap el diámetro medio de los poros. Cuando ambas
difusividades (Fick y Knudsen) tienen relevancia, la difusividad resultante se denomina
“coeficiente de difusión efectivo” y vale

DAB e (2.15)
DABef =
t

ε es la porosidad (fracción hueca) del sólido y τ el factor de tortuosidad de los poros

del sólido que varía entre 1 y 10. Los valores de ε y τ deben ser medidos en forma
experimental. En particular, τ es difícil de medir.

2.3 DERIVACION DE LAS ECUACIONES GENERALES DE CONSERVACION DE MASA.

2.3.1 Ecuaciones Generales

Si se considera un elemento diferencial de volumen ΔXΔYΔZ y una especie A que

difunde (es transportada) dentro de éste y dentro de él hay generación o consumo de
A, se puede hacer un balance general en la forma:

(entrada neta de A por difusión) + (entrada de A por convección) + (velocidad de

generación o desapariciñon de A) = (velocidad de acumulaciñon de A).

O bien:

é flujo neto de A ù (2.16)

éflujo neto de A ù é flujo neto de gene- ù ê ú
ê ú+ ê ú= êacumulación (o desap.) deú
êque entra al vol.ú êración (o desap.) de A en el vol.ú ê ú
ë û ë û ê A en el vol. ú
ë û

El flujo másico transferido en el volumen, en forma vectorial es igual a:

v v
nA = r AUA

y para cada dirección en coordenadas cartesianas será:

n&A X D yD z = r AU&A X D yD z

e igual para los ejes y y z . El flujo neto de salida del elemento será entonces para cada
eje igual a:

n&A D yD z X + D X - n&A X D yD z

46
y en igual forma para los ejes y y z. La velocidad de acumulación de A en el volumen
será a su vez

¶rA
D xD yD z
¶t

y la generación de A (o consumo) en el volumen es

r̂A D xD yD z

sustituyendo en la ecuación (2.18) se tiene

n&A X + D X - n&A X n&A Y+ D Y - n&A Y n&A Z + D Z - n&A Z ¶rA

+ + + ± r̂A = 0
Dx Dy Dz ¶t

Evaluando al límite para Δx, Δy y Δz→0

¶ ¶ ¶ ¶r (2.17a)
n&A ,X + n&A ,Y + n&A,Z + A ± rˆA = 0
¶x ¶y ¶z ¶t

Esta relación representa la ecuación de continuidad del componente A. Como n&A ,X ,

v
n&A ,Y y n&A ,Z son los componentes rectangulares del vector de flujo másico nA la
expresión anterior se puede colocar como la divergencia del flujo másico de A:

v ¶r (2.17b)
Ñ ×nA + A ± rˆA = 0
¶t

Para la especie B en un sistema binario igualmente se puede demostrar que

v ¶r (2.17) c
Ñ ×nB + B ± rˆB = 0
¶t

Para un sistema binario se tiene de (2.17b) y (2.17c) que

v v ¶ (r A + r B )
Ñ ×(nA + nB )+ ± (rˆA + rˆB )= 0 (2.17) d
¶t

Y como en un sistema binario se tiene también que:

v v v v r
nA + nB = r AUA + r BUB = r T uT

y
47
rA + rB = rT

Si se considera que se crea A y consume B según A ƒ B, se puede escribir que

v ¶r (2.18)a
Ñ ×r T u + T = 0
¶t

Esta representa la ecuación de continuidad de una mezcla binaria, la que también

puede escribirse en términos de la derivada sustantiva o material de ρT cómo

D rT v (2.18)b
+ r T Ñ ×u = 0
Dt

Donde

Dr T ¶ r T ¶r ¶r ¶r
= + uX T + uY T + uZ T
Dt ¶t ¶X ¶Y ¶Z

Empleando las definiciones de flujo vistas, se puede colocar igualmente la ecuación de

continuidad para una especie A en la forma

D WA v (2.19)
rT + Ñ ×jA ± rˆA = 0
Dt

El mismo desarrollo puede hacerse en términos de moles. Llamando r& A la velocidad

molar de generación o consumo de la especie A por unidad de volumen y r&

B para B, se

puede escribir:

v ¶ CA (2.20)a
Ñ ×NA + ± ˆrA = 0
¶t

v ¶ CB (2.20)b
Ñ ×NB + ± ˆrB = 0
¶t

y para la mezcla binaria AB,

v v ¶ (C A + CB ) (2.21)a
Ñ ×(NA + NB )+ + (rˆA + rˆB )= 0
¶t

Como se tiene para una mezcla binaria que

v v v v v
NA + NB = C A uA + CBuB = C T uT

48
y que

C A + CB = C T

y si se considera que la reacción ocurre según A € B , es decir, - r&

A = - rB
&

v v ¶ (C A + CB ) (2.21)b
Ñ ×(NA + NB )+ =0
¶t

que es similar a la expresión (2.18)a. En general, se puede también expresar la

ecuación de continuidad en términos de la concentración de las especies A y B en la
forma:

v ¶ CT (2.22)
Ñ ×C T U + ± (rˆA + rˆB )= 0
¶t

2.3.2. Formas particulares de la ecuación general de transferencia de masa.

En ciertos casos, es necesario conocer el perfil de concentración de una especie que es

transportada. Recordando que el flujo másico en un sistema binario puede ser
expresado en la forma:
v v v
nA = - r T DABÑ WA + WA (nA + nB ) (2.23)a

o bien
v v
nA = - r T DABÑ WA + r A u

Substituyendo esta última expresión en la expresión general (2.17)b de conservación

para la especie A

v ¶r
Ñ ×nA + A ± r&
A = 0 (2.23)b
¶t

entonces

Ñ g nA = Ñ g - r T DAB Ñ WA + Ñ g r A u (2.23)c

se tiene

v ¶r
- Ñ ×r T DABÑ WA + Ñ ×r A u + A ± r&
A = 0 (2.24)
¶t

49
Esta expresión permite entonces calcular el perfil de concentración de la especie A en
un sistema binario. En igual forma, puede deducirse la expresión (2.24) empleando el
flujo molar:
v
NA = - C T DABÑ X A + X A (N&A + N&B )

o bien
v v
NA = - C TDABÑ X A + C AU

De igual manera, reemplazando la expresión (2.7) en la (2.17)b se tiene para el perfil

de concentración de la especie A:

v ¶ CA
- Ñ ×C TDABÑ X A + Ñ ×C AU + ± r&
A = 0 (2.25)
¶t

Estas expresiones pueden simplificarse, ya que en su forma general son difíciles de

emplear. Para ello se pueden hacer las siguientes consideraciones:

a) La densidad total del sistema y el coeficiente de difusión son constantes.

De esta forma, por ej., la expresión (2.24) queda:

v v ¶r
- DABÑ 2 r A + r A Ñ ×u + u ×Ñ r A + A ± r&
A = 0
¶t

v
Dividiendo por MA y arreglando términos, para u constante:

v ¶ CA
u ×Ñ C A + = DABÑ 2C A ± r&
A = 0 (2.26)
¶t

b) Si además no hay generación ni desaparición de la especie A (no hay reacción

química).

La expresión anterior toma la forma:

v ¶ CA
u ×Ñ C A + = DABÑ 2C A (2.27)a
¶t

¶ CA v
el término + u ×Ñ C A es la derivada material de CA y entonces esta última relación
¶t
se puede colocar:

50
D CA
= DABÑ 2C A (2.27)b
Dt

Esta ecuación tiene la misma forma que la ecuación de transferencia de calor

¶T
= a Ñ 2T
¶t

lo que indica la similitud entre los fenómenos de transferencia de masa y transferencia

de calor.

c) No hay movimiento convectivo del fluido (u = 0) (además de todas las

condiciones anteriores impuestas).

La relación (2.27)b se reduce a la forma:

¶ CA
= DABÑ 2C A (2.28)
¶t

Esta es la ecuación conocida como la segunda ley de difusión de Fick, y se aplica en

forma rigurosa en los casos de difusión en sólidos; para líquidos estáticos de
relativamente alta viscosidad y sistemas de líquidos y gases binarios en que los flujos
N&A y N&B son iguales pero de sentido opuesto (contra-difusión equimolecular), con las
restricciones impuestas anteriormente.

Esta última expresión también es idéntica a la ecuación de conducción de calor de

Fourier:

¶T
= a Ñ 2T (2.29)
¶t

d) El sistema está en régimen permanente (¶ C A / ¶ t = 0) con densidad de total

del sistema constante y difusividad constante.

La expresión (2.27) se reduce a:

v
u ×Ñ C A = DABÑ 2C A (2.30)
y
D
T

A
B

e) El sistema, además de permanente, r son constantes, no ocurre

reacción química y el fluído está estático

La expresión (2.27) queda solo en la forma:

51
Ñ 2C A = 0 (2.31)

Esta es la ecuación de Laplace en términos de la concentración y es idéntica a la

ecuación de Laplace para transferencia de calor:

Ñ 2T = 0

2.3.3 Transporte de masa por difusión molecular en régimen permanente

Recordando la ecuación general de conservación:

r ¶ CA
Ñ g NA + ± r&
A = 0 (2.31)a
¶t

C
/
t
A
Dependiendo de las condiciones impuestas, los términos de acumulación ¶ ¶ y
generación o consumo de la especie A pueden ser cero. Varios casos se pueden
presentar:

a) Transporte en régimen permanente unidireccional sin reacción química

En este caso, tanto la concentración como el flujo másico son funciones de solo una
coordenada. Para un sistema binario se tiene entonces que:

dX
(
N&A z = - C T DAB A + X A N&A z + N&B z
dz
) (2.32)

A su vez se pueden tener varios casos de esta condición:

- Transporte a través de un gas estático

Si se considera un líquido (metal líquido, líquido inorgánico, etc.) y un tubo de pequeña

sección (Fig. 2.5) (celda de difusión de Arnold) el líquido A se mantiene a T y P
constantes, en tanto que el gas B es insoluble en el líquido A y químicamente inerte
con respecto a éste. El líquido A se vaporiza y su gas difunden a lo largo del tubo y
mezcla con el gas B que fluye sobre él.

52
 C ºB º
CA
Gas B Gas A+B
(Portador) z = z2
NA v
v
z z+ z
z

Gas A
NA
z

z = z1
C iB CiA
C
Líquido A

Figura 2.5 Celda de difusión de Arnold

Para una sección S del tubo, se tiene en régimen permanente que

S N&A z+ D Z - S N&A z = 0

Dividendo por S∆z y evaluando al límite ∆z→0

dN&A z
=0
dz

Es decir, el flujo molar de A en dirección z es conservativo, lo que también puede

deducirse de la expresión general (2.22) a:

r ¶ CA
Ñ gNA + ± r&
A = 0 (2.32)a
¶t

Como ¶ C A / ¶ t = 0 y r&
A =0

¶ N&A ¶ N&A y
¶ N&A z
x
+ + =0
¶x ¶y ¶z

Para flujo unidireccional el primer y segundo términos de esta expresión son cero,
entonces

¶ N&A z
=0
¶z

53
Puesto que el flujo de B es cero (estático) en el interior del tubo puesto que N&B Z = 0 y
el flujo molar de A a su vez en régimen permanente se puede expresar en la forma
vista anteriormente.

dX
(
N&A z = - C T DAB A + X A N&A z + N&B z
dZ
) (2.32)b

luego,

C D ædX ö
N&A z = - T AB çç A ÷ ÷ (2.33)
1 - X A çè dz ÷
ø

NA
Esta ecuación se puede resolver considerando que & es constante a lo largo de z y

z
que DAB también lo es, es decir, es independiente de CA.

Las CB son en este caso:

z = z1 z=z2
xA=xA1 xA=xA2

Resolviendo:

C D é1 - X A ù
N&A Z = T AB ln êê 2 ú
(2.34)
z2 - z1 ëê1 - X A1 úûú

Se puede definir la concentración logarítmica media de XB en la forma:

XB2 - XB1
XB,ln = (2.35)
ln(XB2 XB1 )

O bien, para una mezcla binaria A-B:

(1 - X A1 )- (1 - X A2 ) X A1 - X A2
XB,ln = =
ln éë(1 - X A2 )/ (1 - X A1 )ùû ln éë(1 - X A2 )/ (1 - X A1 )ùû

Reemplazando esta última expresión en la (2.34), se tiene:

C D (X A1 - X A2 )
N&A Z = T AB (2.36)
(Z2 - Z1 ) XB,1n

Esta relación también se puede escribir en términos de las presiones parciales de A y B.

Para una gas ideal, se tiene
54
C A = nA / V = PA / RT
X A = pA / P

de donde se tiene que

DABP (pA1 - p A2 )
N&A Z = (2.37)
R T (z2 - z1 ) pB,ln

Las expresiones (2.36 y (2.37) representan las ecuaciones de difusión de un fluido o gas
en otro inerte, estanco.

- Contradifusión equimolecular en régimen permanente

En este caso, la transferencia de las especies (sistema binario) es igual y opuesta

N&A Z = N&B Z y la ecuación general de conservación se reduce a:

Ñ g NA = 0 (2.38)

Para transporte unidireccional se tiene:

&
dNA z
=0 (2.39)
dz

es decir N&A = cons tante a lo largo de z

Como se tiene que

dC
dZ
(
N&A Z = - DAB A + X A N&A Z + N&B Z ) y N&A Z
= - N&B z
(2.40)
dC
N&A Z = - DAB A
dZ

Ecuación que puede resolverse para

z = z1 z = z2
CA = CA1 CA = CA2

Luego,

æ  D ö
N& A z == ç AB ÷ ÷(C A1 - C A2 )
÷ (2.41)
è Zz22 − Zz1 ÷
ø

55
ecuación que también es válida para un gas ideal, considerando que

C A = nA / V = pA / RT

DAB
N&A z = (pA1 - pA2 ) (2.42)
R T (z2 - z1 )

El perfil de concentración de la especie A o B puede obtenerse reemplazando la

expresión del flujo de A (ó B) (2.41) en la relación que describe al flujo de A en
dirección Z (flujo permanente):

&
dNA z
=0 (2.43)
dz

Es decir,

d2CA
= 0
d2 z

Ecuación diferencial que al resolverla se tiene:

CA = C1z+C2
Las ctes. pueden evaluarse a su vez de

z = z1 z = z2
CA = CA1 CA =CA2

De esta manera, el perfil de la especie A tiene la forma lineal:

C A - C A1 z - z1
= (2.44)
C A1 - C A2 z1 - z2

b) Transferencia unidireccional con reacción química

Este caso se presenta en forma relativamente frecuente en metalurgia y puede ocurrir

tanto en fase homogénea como heterogénea. En estos casos, se debe notar que en la
ecuación general de conservación para una especie:

¶ CA
Ñ g NA + ± r&
A = 0 (2.45)
¶t

El término r&
a aparece solo en reacciones en fase homogénea. En fases heterogéneos

no aparece, puesto que no ocurre dentro del volumen de control desde el cual se

56
desarrolló la expresión (2.22)a, pero puede aparecer en las condiciones de borde
impuestas (interfase reaccional).

Se pueden presentar varios casos, sin embargo el de mayor interés lo representan las
reacciones heterogéneas.

- Transporte y reacción química en fase heterogénea con variación de la interfase

reaccional.

Este caso ocurre cuando la reacción química es muy rápida y la transformación global
depende de la velocidad con que el reactivo (fluído) llegue a la interfase (o el producto
sea eliminado de ésta). Este tipo de reacción es común en sólidos que reaccionan con
gases a alta temperatura, en este caso, se trata de control por transporte de la especie
en el gas desde el seno de éste hacia la interfase reaccional, o viceversa.

Un caso típico lo representa la combustión de carbón con aire en ausencia de flujo

convectivo. Este es un caso hipotético, puesto que industrialmente este proceso no
sería práctico ya que la combustión sería relativamente lenta.

Si se considera la reacción global

2C(s) + O2(g) = 2CO(g)

El oxígeno es transportado desde el seno del fluido (aire) hacia la interfase reaccional,
del carbón la cual al combustionarse decrece en tamaño (Fig. 2.6).

Figura 2.6. Combustión de carbón con oxígeno en ausencia de flujo convectivo.

Si se consideran partículas esféricas, flujo unidireccional y que no hay reacción entre el

oxígeno y CO, el balance para el oxígeno en un sistema polar será:

57
(N&
O2 r 4 p r2
r
)- (N&
O2 r + D r 4 p r2
r+ D r
)= 0 (2.46)

Dividiendo por un elemento de volumen 4p r2D r y evaluando al límite:

d 2
dr
(r NO2 r )= 0 (2.47)

Esta expresión indica que r2N&O2 es constante para cualquier valor del radio. Si se
r

comparan las expresiones (2.44) con la última (2.47) se observa que son diferentes,
debido al distinto sistema de coordenadas empleado en este último caso (esféricas):

dN&
=0 (2.48)
dz Z

d 2 &
dr
(
r NO2 = 0
r
) (2.49)

Esta última expresión puede también deducirse directamente de la ecuación de

continuidad tal como la (2.26):

¶ CA
- Ñ g C T DABÑ X A + Ñ g C AU + ± r&
A = 0 (2.50)
¶t

la que en coordenadas esféricas y en ausencia de flujo convectivo con r T y DAB

constantes, será:

∂C A  1 ∂ 2 & 1 ∂ & 1 ∂N& A φ 

+  2 (r NA r ) + (NA senθ) +  − rA = 0 (2.51)
∂t  r ∂ r senθ ∂θ r senθ ∂φ 

Para régimen permanente y sin generación ni consumo de oxigeno en el volumen de

control, flujo unidireccional (los flujos en θ y φ son nulos) se tiene que:

d 2
dr
( r NCO r ) = 0 (2.52)

Por estequeometría, en la interfase se tiene que

- N&CO r = 2N&O2 r

válido en cualquier punto radial. Por otra parte, para un sistema pseudo binario O2-
aire se tiene:

58
 dX O2
- N&O2 r = - C TDO 2 - Aire
dr
(
+ X 02 N&O2 r + N&N2
r
)
Como el nitrógeno no difunde, se puede simplificar a un binario real:

dX
r dr r
(
N&O2 = - C T DO2 - Aire O2 + X O2 N&O2 + N&CO r ) (2.53)

Ecuación diferencial que puede resolverse con las CB

r = ro r=∞
x O2 = x O2 x O2 = 0.21
ro

Es más útil sin embargo conocer la cantidad total de masa transferida, la que puede
& en términos de flujo,
calcularse definiendo la velocidad de transferencia de masa MO2

en un radio dado:

&O = 4p r 2N&O
M (Mt- 1 )
2 2 r

& al igual que N& es constante a lo largo de r. Reemplazando esta última expresión
MO2 O2

en (2.52), se tiene:

2
&O = - 4p r C TDO2 - aire ççdX O2 ÷
æ ö
M ÷ (2.54)
2
1 + X O2 è dr ø÷
ç

o bien,

dr æ dX O2 ö÷
Mo2 = - = - 4 p C D ç
ç ÷ (2.55)
r2
T O2 - aire
çè1 + X O2 ø÷
÷

Resolviendo para las CB vistas anteriormente,

é1 + X O2 ù
&O = 4p C TDO - aire ro ln ê
M
ro ú
(2.56)
2 2 ê1 + 0.21 ú
êë úû

De esta forma, mediante la expresión (2.56) se puede calcular la velocidad a la cual es

transferido O2 hacia el sólido, la velocidad de combustión del carbón y la cantidad de
energía generada.

Otros modelos más complejos y más reales se han propuesto para la combustión de
carbón y oxidación de sulfuros que se verán más adelante y que consideran el flujo
convectivo y la reacción química interfacial.

59
2.4 TRANSFERENCIA DE MASA CONVECTIVA

En la mayoria de los casos reales en que hay involucrado un fluído la transferencia de

masa no ocurre solo por difusión molecular, sino que también debido al efecto
convectivo del fluido causado por una gradiente de densidad, temperatura
(transferencia convectiva natural) o un efecto mecánico impuesto al fluído
(transferencia convectiva forzada) tal como una agitación mecánica.

La ecuación de velocidad de transferencia (flujo) se puede expresar en forma similar a

la de enfriamiento de un cuerpo (ecuación de Newton):

N&A Z = hmD C A (2.57)

Donde N&A es el flujo convectivo unitario de la especie A ,∇CA la gradiente de

concentración y hm el coeficiente de transferencia de masa. Este coeficiente incluye
tanto los efectos convectivos naturales y forzados como el coeficiente de difusión
molecular y depende de la geometría del sistema, propiedades del fluído, régimen de
flujo y gradiente de concentración.

Cuando la transferencia de masa la representa por ej. un sólido disolviéndose

(lixiviación) o un gas reaccionando con un fluído, de mecánica de fluídos se ha visto
que existe un film que rodea las partículas (interfase reaccional) donde el flujo es
laminar y no turbulento. El mecanismo de transferencia de masa a través de esta capa
estática o cuasi estática debe entonces involucrar transferencia por difusión molecular.
En este caso, el flujo de la especie A desde el seno del fluido (C oA ) hasta la superficie
reaccional (C sA ) se expresa en la forma:

N&A Z = hm (CSA - C0A ) (2.58)

La resistencia a la transferencia de masa es entonces el resultado de la suma de las

resistencias de la transferencia desde el seno del fluido hasta el film de fluido y del film
de fluído. Si el fluido es turbulento, no hay resistencia en el fluido y toda esta reside
en el film de fluido por lo que el coeficiente de transferencia de masa también se le
denomina coeficiente de film.

La determinación del coeficiente de transferencia de masa no es simple y existen

varias formas de cálculo o estimación de éste:

- Análisis dimesional
- Análisis exacto de la ecuación de capa límite
- Solución aproximada de la ecuación de capa límite
- Analogía entre momento, energía y transferencia de masa.
- Teorías del film y de renovación de superficie.

60
Dependiendo si la transferencia de masa ocurre en flujo convectivo natural o en flujo
convectivo forazado, las expresiones para el coeficiente de transferencia de masa se
pueden adquirir en función de los números adimensionales de Sherwood, Freude,
Lewis, Grashoff, Schmidt, Reynolds y Péclet:

- Convección natural Sh=f (Fr, Sc)

- Convección forzada Sh=f (Re, Sc)
- Convección natural y forzada Sh=f (Fr, Re, Sc)

Lur Fuerzas inerciales

R e= = (Reynolds)
m fuerzas viscosas

m difusividad convectiva
S c= = (Schmidt)
rD difusividad molecular

hmL transferencia global de masa

S h= = (Sherwood) (Nusselt para Transf. de masa)
D difusividad molecular

g r 2L3 (r 0 - r s ) æfuerzas ascensionales ö÷

G r= = Reçç ÷ (Grashoff)
m2 çè fuerzas viscosas ø÷

Lu flujo total masa

P e= = (Péclet)
D flujo difusivo masa

m2 fuerzas inerciales
F r= = (Freude)
gL fuerzas gravitacionales

k difusividad térmica
L e= = (Lewis)
r CP D difusividad másica

Donde:

L = Dimensión característica del sistema (diámetro, longitud, etc.)

u = Velocidad relativa del fluído
µ = Viscosidad del fluído
ρ o , ρs = Densidad del fluído en su seno y superficie del sólido
D = Coeficiente de difusión
k = Conductividad térmica del fluído
ν = Viscosidad cinemática

61
La razón entre el flujo total de masa transferido y el flujo difusivo de masa puede
expresarse mediante el número de Péclet definido como Pe=Lu/D en que D es un
coeficiente de difusión equivalente, de iguales dimensiones que D.

2.4.1 Solución exacta de la capa límite de concentración

Debido a la complejidad que representa la geometría de un cuerpo irregular, la

solución analítica de la capa límite se ha efectuado sólo para la capa formada sobre un
cuerpo plano en que el flujo es laminar y paralelo a la superficie. La solución para la
capa de concentración es idéntica a la solución de Blasius para la capa límite
hidrodinámica y la capa térmica, de manera que su deducción se hace en forma
resumida.

Para la transferencia de momento en 2-D en régimen permanente de un fluído

incompresible, la ecuación de continuidad tomaba la forma de:

¶ vx ¶ vy
+ =0 (2.59)
¶x ¶y

La ecuación de movimiento en dirección x a su vez, para viscosidad y presión constante

tenía la forma:

¶ vx ¶v ¶ 2v
vx + v y y = n 2x (2.60)
¶x ¶y ¶y

Como se vió en transferencia de calor, la ecuación equivalente a esta última y que

describe la transferencia de energía térmica es un sistema incompresible en 2-D con
flujo isobárico y difusividad térmica constante tenía la forma:

¶T ¶T ¶ 2T
vx + vy =a 2 (2.61)
¶x ¶y ¶y

Donde α es la difusividad térmica, de igual significado que ν y D. Una expresión

análoga puede ser derivada para la transferencia de masa en el film de concentración
(Fig. 2.3).

62
 A
oo
CA

Capa límite
y
s
CA CA = CA (y)

Figura 2.3 Capa límite de concentración para un fluido que se mueve paralelo a una
superficie plana.

Si no existe reacción química dentro de la capa límite y la segunda derivada de CA

respecto al eje x se puede despreciar frente a la primera derivada respecto al eje, en
régimen permanente, fluído incompresible, difusividad constante y flujo 2-D, la
expresión para la ecuación de concentración se puede escribir, al igual que (2.60):

¶ CA ¶C ¶ 2C A
vx + v y A = DAB (2.62)
¶x ¶y ¶ y2

Las condiciones de borde para las tres ecuaciones vistas (2.60), (2.61) y (2.62) son:

- Momento

y= 0 y= ¥

v x - v xS v x - v xS
=0 =1
v ¥x - v xS v ¥x - v xS

- Calor

y= 0 y= ¥

T - TS T - TS
=0 =1
T¥ - TS T¥ - TS

- Concentración

y= 0 y= ¥

63
 C A - CSA C A - CSA
=0 =1
C¥A - C SA C¥A - CSA

La similitud de las tres ecuaciones diferenciales y sus condiciones de borde sugiere que
la solución de Blasius puede ser aplicada también a la capa límite de concentración.
Para la transferencia de calor, la solución fue aplicada cuando el número de Prandtl = 1
(Pr = Cp µ/k). Igual caso será válido en la transferencia de masa cuando número de
Schmidt = 1.

Blasius obtuvo una solución para la capa límite de velocidad definiendo la función
potencial ψ. De esta forma, las ecuaciones diferenciales parciales se reducen a una
ecuación diferencial ordinaria al transformar las variables independientes x e y en η y
la variable dependiente ψ (x,y) a f(η), donde f es un parámetro de similitud respecto a
η. De esta forma, se tiene para la transferencia de masa

dC A æ0.332 0.5 ö÷
= (C¥A - CSA )çç Rex ÷ (2.63)
dy y= 0
çè x ø÷

La solución de Blasius no considera velocidad en y sobre la superficie (y=0), lo que

aplicado a la solución de concentración implica que la velocidad a la cual la especie que
difunde es tal que no altera el perfil de la capa límite de concentración.

Cuando la velocidad en y es cero, Vys 0 el término convectivo en la ecuación de Fick

es también cero (en dirección y). De esta forma, el flujo desde la superficie plana hacia
el fluído (o viceversa) a través de la capa límite se puede expresar como:

dC
N&A y = - DAB A (2.64)
dy y= 0

De (2.62) y (2.63) entonces:

æ0.332 0.5 ÷
ö S
N&A y = DAB çç Rex ÷ (C - C¥A ) (2.65)
çè x ø A
÷

Definiendo el flujo en términos del coeficiente de transferencia de masa hm:

N&A y = hm (CSA - C¥A )

se tiene:

64
 DAB
hm =
x
(0.332Re0.5
x )

o bien,

hmx
Sh = = 0.332Re0.5
x (Sh = Nu para transferencia de masa) (2.66)
DAB

Esta última expresión es válida sólo para Sc 1 y valores bajos del coeficiente de
transferencia de masa desde o hacia la superficie plana. En la Fig. 2.7 se han trazado
los perfiles de concentración para este caso, graficando el parámetro η en la forma:

y C SA - C A
h= (Re0.5
x ) versus
x CSA - C¥A

1.0
. 5
-2 0
0.8
25

6
CAs - CA 0.6 0.
0.

s
0.

Vy
(Re x) 0.5
CAs - oo
CA 0.4 V oo

0.2
0.0
0 2 4 8 6 10 12
y 0.5
= x (Re x )
Figura 2.7 Perfiles de concentración para la transferencia de masa a través de un film
de concentración sobre una superficie plana en régimen laminar.

En la Figura 2.7 la solución gráfica se ha trazado para valores (+) y (-) del parámetro de
capa límite ( Vys / V ∞ ) (Re x )0.5 . Los valores (+) son válidos para la transferencia de masa
desde la superficie plana hacia la capa límite y los (-) para el caso inverso.

Cuando el valor se hace cero, la transferencia de masa decrece hasta no tener

influencia en el perfil de concentración. La pendiente de la curva vale 0.332 para y = 0.

2.4.2 Concepto de coeficiente de transferencia de masa

En la mayor parte de los casos en ingeniería más importante que conocer el perfil de
concentración de una especie que es transferida es conocer el flujo de masa desde el

65
fluído al sólido (o viceversa), el que en general puede expresarse por ejemplo,
mediante la relación (2.7) para sistemas binarios:

dX
N&A Z = C T DAB A + X A N&A Z + N&B Z
dz
( )
Sin embargo, no siempre es posible encontrar el flujo de masa debido a lo complejo
del problema o cuando no existe un perfil de concentración definido (flujo turbulento).
En estos casos, es más conveniente emplear el coeficiente de transferencia de masa
definido en la forma:

N&A = hmD C A (2.67)

Como se vió anteriormente hm es un coeficiente de proporcionalidad con dimensiones

(flujo molar/gradiente molar) ó (flujo másico/gradiente másico). Cuando hm varía en
un sistema de punto a punto y en el tiempo, es conveniente emplear un valor medio
del coeficiente de transferencia integrando para una superficie y tiempo definido:

S t

ò ò h dtdS m
0 0
hm = S t
(2.68)
ò dSò dt
0 0

Las derivaciones de las ecuaciones analíticas para el cálculo del coeficiente de

transferencia de masa son, desde el punto de vista del rigor analítico las únicas
exactas, sin embargo en la mayor parte de los casos la solución de las ecuaciones
diferenciales es muy compleja como para ser de utilidad. En tales casos, se puede
emplear teorías simplificadas que representan con suficiente aproximación la realidad.
Dos de ellas son la Teoría de renovación de superficie y Teoría del film.

a) Teoría del film

En esta teoría se postula que, independientemente del régimen de flujo (laminar o

turbulento), toda la resistencia a la transferencia de masa reside en un film estático
sobre el sólido de espesor δm . Esto también es válido para la transferencia de masa
entre dos fluidos no miscibles entre sí, como por ej., en extracción por solventes entre
la fase orgánica y la fase acuosa o entre un metal fundido y la escoria.

Cuando difunde una especie A en B en un sistema binario estacionario de fluido B en

régimen permanente, unidireccional:

N&A Z = N&A Z= 0
= N&A Z= L

66
y como

dX
(
N&A Z = - C TDAB A + X A N&A Z + N&B Z
dz
)
Si la especie B está estacionaria o el flujo es contradifusión equimolecular:

N&B Z = 0

C D ædX ö
N&A Z = - T AB çç A ÷ ÷
1 - X A çè dz ÷
ø

æN& ö dX A
çç A Z= 0 ÷ ÷
çç ÷
÷dz = -
çè C TDAB ÷
ø 1 - XA

Resolviendo para las CB:

ïìï z = 0 z = L = dm
í
ïïî X A = X SA X A = X 0A

En que XsA es la concentración en la superficie y X ºA en el seno del fluido. De esta

forma, se tiene:

æN& ö÷ æ1 - X 0A ö÷
ççç A Z= 0 ÷ d
÷m = ln çç ÷ (2.69)
÷
ççè C TDAB ø÷ çè1 - X SA ø÷
÷

o bien,

C TDAB æ S ö
ççC A - C A ÷
N&A = ÷
÷ (2.70)
Z= 0
dn çèç CB,ln ÷ ø

Si se compara esta última expresión con la expresión anteriormente vista para el

coeficiente de transferencia de masa según:

N&A = hm (CSA - C0A )

se observa que hm es un coeficiente de proporcionalidad igual a

C TDAB
hm = (2.71)
dnCB,ln

67
En que δm es el espesor de la capa límite de concentración.

Si la especie que difunde (A en este caso) está muy diluida, C T ; CB,ln puesto que
CsA y CºA << C T y la expresión se simplifica a:

DAB
hm ;
dm

Aunque esta expresión es simple de emplear, δm depende también de otras variables

(u, S, T, µ, etc.) así como de la geometría del sistema. A pesar de ello y debido a la
simplicidad del modelo, se usa con frecuencia ya que es posible estimar con cierta
precisión el valor de δm . Para ello, existen relaciones empíricas desarrolladas para
calcular el espesor de la capa de difusión o de concentración δm . Por ejemplo:

- Para flujo convectivo laminar forzado sobre una superficie la expresión tiene la
forma:

0.5
æL ö 0.167 0.33
÷
dm @3ççç ÷
÷ n DAB (2.72)
èu ÷
0ø

uo es la velocidad en el seno del fluído paralelo o por el exterior de la superficie, ν es la

viscosidad cinemática del fluído y L la longitud de la superficie sobre la cual se forma el
film de espesor δm .

- Para flujo convectivo forzado turbulento sobre una superficie, la expresión toma
la forma:

0.5
æL0.1 ö 0.57 0.33
dm @3ççç 0.9 ÷
÷ n DAB
÷ (2.73)
èu ÷ ø
o

L es la dimensión característica del sistema.

b) Teoría de renovación de superficie

Se conoce también con el nombre de Teoría de Higbie y aún cuando el planteamiento

físico puede parecer ideal, ha resultado bastante aproximada en muchos casos.

Si se considera un flujo turbulento de un fluído A en contacto con una superficie B (que

puede ser un sólido o dos fluídos inmiscibles), se puede considerar que continuamente
elementos del fluído (cantidades discretas de éste) entran en contacto con la
superficie por un corto tiempo durante el cual ocurre difusión en régimen transciente

68
desde o hacia el elemento de fluido y el sólido (o líquido y gas o dos líquidos
inmiscibles, etc.) (Fig. 2.8).

A(f) A(f) A(f)

Seno del fluído
CAº = cte. º
CA º
CA
C

Elemento B
de fluído

(Sólido o líquido) B(s) (Sólido o líquido) B(s) (Sólido o líquido) B(s)

t<0 t=0 t>0

Figura 2.8 Representación del modelo de Higbie o de renovación de superficie.

Si se considera que el tiempo de contacto entre el elemento de fluido y el sólido (u

otro fluido) es suficientemente corto, el problema se resuelve considerando el
elemento del fluido A como un cuerpo semi infinito, es decir, la cara del volumen de
fluido A en contacto con el seno del fluído no varía de composición durante el tiempo
de contacto. En estas condiciones y si se considera que se transfiere la especie sólida B
al fluido A y que la concentración de B en A es diluida y el flujo es unidireccional
normal a la superficie de B, para un elemento de fluido se tenía en general de (2.29)
que

¶ CB
= DABÑ 2CB (2.74)
¶t

es decir para flujo unidireccional

¶ CB ¶ 2CB
= DAB (2.75)
¶t ¶ z2

para 0 ≤ z ≤ ∞ y 0 ≤ t ≤ tc donde tc es el tiempo total de contacto del elemento de

fluído con el sólido.

Las CB son en este caso:

I) z= 0 CB = CB0 t£ 0
II) z® ¥ CB = CB0 t¹ 0
s
III) z= 0 CB = C B t> 0

La solución de esta ecuación diferencial es idéntica a la encontrada para el flujo de

calor transciente en un cuerpo semi infinito y tiene la forma general:

CB = M + N cfer éêz 2 DAB t ùú

ë û
69
 Para las BC impuestas en este caso,

CB - CB0 ìï z üïï
= cfer íï (2.76)
ý
CBs - CB0 ïï 2 DAB t ïï
î þ

donde:

ïì z ïüï 2 z/2 DAB t

fer ïí
2
ý = 0.5 ò e- h dh (cfer = 1 - fer)
ïï 2 DAB t ïï p 0
î þ

Donde η es función de DAB y t. La función error llamada también integral error de

Gauss, se encuentra tabulada. Por otra parte, el flujo difusivo desde la interfase
reaccional (despreciando el convectivo) vale:

¶C DAB
N&A = - DAB B = CB - CB (2.77)
¶z z= 0 pt

Diferenciando la función error de la expresión (2.76) y reemplazando en la (2.77) se

tiene:

0.5
CB & æD ö
NA = - DAB = (CBs - CB0 )çç AB ÷
NA =
¶z çè p t ÷
÷
ø

La cantidad total de la especie A transferida en el tiempo tc por unidad de área

interfacial será entonces:

0.5
tc æD ö tc dt
ò N&Bdt = (CBs - CB0 )çç AB ÷ ò (2.78)
0 çè p ø÷
÷ 0 t0.5

luego:

0.5
æD ö
N&A = 2(C - C )çç AB ÷
s 0
÷ (2.79)
çè p t ÷
B B
÷
ø c

Esta última relación da el valor promedio de la especie A transferida en el tiempo de

contacto tc. Si se compara esta última correlación con la vista para hm:

N&B = hm (CBs - CB0 )

se puede observar que para concentraciones bajas de B en A el coeficiente de

transferencia de masa vale:

70
 0.5
æD ö
hm = 2çç AB ÷ ÷ (2.80)
çè p t ÷
c
÷
ø

Se observa además que un parámetro de importancia en la teoría de Higbie es el

tiempo de contacto tc el que puede estimarse a partir del flujo del fluido, o en ciertos
casos medirse directamente, como por ejemplo, el tiempo de permanencia de las
burbujas de gas (CO) en un reactor de refinación de acero; de las burbujas de gas en un
lecho fluidizado, de las gotas de orgánico en el acuoso en extracción por solvente, etc.
En estos casos, tc ; db / ub en que db y ub son el diámetro de las burbujas o gotas y ub
su velocidad relativa al fluido.

Se puede observar que la teoría de renovación de superficie no considera la existencia

de un film, el que para geometrías irregulares del sólido o velocidades altas es
virtualmente imposible de estimar.

2.4.3 Intercambio de masa en interfases fluido-fluido

Cuando existen dos fases líquidas o fundidas inmiscibles tal como escoria-metal
fundido, eje-escoria fundidos, o dos líquidos inmiscibles entre sí como orgánico –
acuoso en extracción por solvente, una o más especies pueden difundir desde una de
ellas a la otra a través de la interfase. Considerando que una especie A que se
encuentra en la Fase 1 se transfiere hacia la Fase 2 (Fig. 2.9):

º
CA
Fase 1 1
Fase 1 Perfil de conc.
de la especie A
A1 en la fase 1
*
CA m A1
Interfase 1
*
CA m A2
2

A2 Perfil de conc.
de la especie A
en la fase 1 º
CA Fase 2
Fase 2 2

Figura 2.9 Flujo de una especie entre dos fases inmiscibles.

En un caso como éste C*A1 no tiene por que ser igual a C*A2 . En régimen permanente
esto se puede expresar en la forma:

N&A = hm1 (C0A1 - C*A1 )= hm2 (C0A2 - C*A2 ) (2.81)

En que C*A y CoA denotan las concentraciones de A en la interfase y el seno de las fases,
respectivamente y hm1 y hm2 son las coeficientes de transferencia de masa de A en las
71
fases 1 y 2. Como en la interfase debe existir equilibrio entre la especie, se tiene que
K=( C*A1 / C*A2 ) en que la constante de equilibrio K se denomina también coeficiente de
partición de A entre las fases 1 y 2.

Arreglando la expresión (2.80) se puede eliminar las concentraciones en la interfase,

las que no se conocen. Se puede colocar entonces que:

KC0 - C0A2
N&A = A1 (2.82)
K 1
+
hm1 hm2

En esta expresión, el término (KC 0A1 − C0A2 ) representa el gradiente de concentración

global de transferencia de masa y (K / hm1 + 1 / hm2 ) la resistencia global a la
transferencia de masa entre las fases y 1 y 2.

Cuando el sistema es multicomponente, por ejemplo existen las especies A, B y C en

las fases 1 y 2 de manera que la especie A existente en la fase 1 reacciona en la
interfase con la especie B existente en la fase 2 para generar la especie C según:

A1 + B2 ® C1

y la especie C difunde en la fase 1, la constante de equilibrio para la interfase se puede

entonces expresar en la forma

C*C1
K= (2.83)
C*A1 ×C*B2

C* denota las concentraciones de equilibrio en la interfase entre las fases 1 y 2. Existen

entonces tres coeficientes de transferencia de masa: (hmA)1, (hmB)2, (hmc)1. Si se supone
régimen permanente y reacciones equimoleculares, los flujos son iguales, entonces:

N&A = N&B = N&C

es decir,

N&A1 = (hmA ) (C 0A1 - C*A1 )

1

N&B2 = (hmB ) (CB02 - CB*2 )

2

N&C1 = (hmC ) (CC01 - C*C1 )

1

72
De esta forma se pueden eliminar las concentraciones en la interfase que son
desconocidas y colocar las expresiones en función la constante de equilibrio K
resultando así:
ïïì æ
1 çç 0 (hmB )2 (hmB )2 0 ÷÷ö
& & & ï 0
NA1 = NB2 = NC1 = (hmA ) í C A1 - çC A1 - - CB ÷
2 ççè ÷-
1ï
ïîï K(hmC ) K(hmA ) 2 ÷ ÷
ø
1 1

é (h )
ê mB 2 (hmB )2 0 0 (hmB )2 0 (hmB )2 ùúüïïï
ê + CB - C A1 + CA + úý (2.84)
êK(hmC ) K (hmA ) 2 K (hmC ) 1 K (hmA ) úïï
ë 1 1 1 1 ûïþ

Un ejemplo de este tipo de reacción lo representa la transferencia de manganeso

desde el baño de acero a la escoria durante el período de refinación del acero a 1550-
1600ºC, donde el manganeso disuelto en el acero reacciona en la interfase con el FeO
de la escoria para generar MnO (que pasa a la escoria) y Fe que se integra al acero:

Mn(m) + FeO(e) → MnO(e) + Fe(l)

2.4.4 Correlaciones de cálculo en transferencia de masa convectiva

Debido a lo limitado de las soluciones analíticas en el caso de geometrías y/o

condiciones de flujo no laminar, para cálculos de ingeniería se han desarrollado
expresiones semi-analíticas y empíricas para estimar los valores del coeficiente de
transferencia o las gradientes de concentración en casos específicos de geometría y
patrones de flujo. Para algunos de estos casos se tiene:

a) Transferencia de masa en tubos y ductos

- Flujo laminar

Para la concentración de salida de un ducto en el que hay transferencia desde la pared

hacia el fluido o viceversa, se puede emplear la correlación gráfica de Litton y
Sherwood, basada en valores experimentales (Fig. 2.10)
1

10 -1 Pe r
fil p
l an o
-2
- º 10 Pe
A A rf il
par
*- º
A A
ab
ó lic
10 -3 o

10 -4

10 -5
1 10 102 103 104 105 106 107 108 109
D
ReSc ( AB )
Z 4
Figura 2.10 Transferencia de masa en flujo convectivo laminar en tubos y ductos
(desde o hacia el fluido) (Litton & Sherwood).
73
Para flujo parabólico se puede emplear la expresión de Lévêque:

r A - r 0A - 23
* 0
= 5.5(n&A r TDABz) (2.85)
rA - rA

Para flujo laminar (zona lineal) se puede usar a su vez la expresión:

r A - r 0A - 12
* 0
= 4.0(n&A r TDABz) (2.86)
rA - rA

Donde z es la distancia axial desde la entrada del ducto o desde el punto de formación
del perfil plano.

- Flujo turbulento

Se puede emplear la correlación experimental de Gilliland y Sherwood:

0.83 0.44
hL æLu r ö æ m ÷ ö
Sh = m = 0.023çç o T ÷
÷ çç ÷
÷ (2.87)
DAB çè m ÷÷
ø çè r TDAB ÷
ø

Donde L es la dimensión característica del sistema (por ej. el diámetro del ducto), uo la
velocidad media del fluido y µ su viscosidad.

Para gases en iguales condiciones el número de Sherwood se expresa en la forma:

hmL æççp A,ln ÷

ö
Sh = ÷
DAB çè P ÷ ø

Donde P es la presión total del sistema y pA ,ln la presión logarítmica media de A.

Otra expresión similar es la relación

Sh=0.023Re0.8 Sc0.33 (2.88)

- Partículas planas en régimen convectivo natural

Se puede aplicar la relación de Wilke:

Sh = 0.673 Gr 0.25 Sc0.25 (2.89)

Válida para los rangos de 500 < Sc 800,000 y 104 < Gr < 109

74
- Partículas o placas planas en régimen convectivo forzado.

Se puede emplear la expresión:

Sh = 0.037Re0.8 Sc0.33 (2.90)

- Partículas o cuerpos esféricos en régimen convectivo forzado.

Se puede aplicar la expresión de Ranz y Marshall:

Sh = 2.0 + 0.6Re0.5 Sc0.33 (2.91)

O bien las expresiones de Rowe y Claxton:

aire : Sh = 2.0 + 0.69Re0.5 Sc0.33 (2.92)

agua : Sh = 2.0 + 0.79Re0.5 Sc0.33 (2.93)

BIBLIOGRAFÍA CAPITULO II

1. J. Welty, C. Wicks and R. Wilson, “Fundamentals of momentum, Heat and mass

transfer”, Willy, (1967).

2. D.R. Poirier and G.H. Geiger, “Transport phenomena in materials processing”, 2nd
Ed., TMS, (1999).

3. H.Y. Sohn and M. Wadsworth, Ed., “Rate processes of Extractive Metallurgy”,

Plenum Press, (1979).

75
 CAPITULO III
REACCIONES HETEROGENEAS

3.0 GENERALIDADES SOBRE REACCIONES HETEROGENEAS (NO CATALITICAS)

A diferencia de las reacciones homogéneas (sólido-sólido, líquido-líquido o gas-gas), en

las reacciones heterogéneas y como en las sólido-gas, líquido-gas o sólido-líquido,
pueden ocurrir otros tipos de mecanismos de control además de la reacción química,
ya sea en forma individual o en conjunto: estos son los fenómenos de transporte de
masa en fase homogénea en el fluído (ya vista en capítulos anteriores), la difusión de la
(s) especies (s) reaccionantes a través del reactivo o de la capa de producto que se
pueda formar (o producto gaseoso o líquido que fluye hacia o desde el exterior a
través de esa capa) y la reacción química intrínseca que ocurre en la superficie
reaccional o volumen reaccional.

Las reacciones heterogéneas sólido- fluido o líquido-gas por lo tanto son generalmente
más complejas que las homogéneas, puesto que en estas últimas se hace abstracción
del mecanismo íntimo molecular y se trata como un medio continuo, no así en las
sólido-fluido en las que los mecanismos de trasferencia puden ser los mecanismos de
control de la transformación del sistema en su totalidad. En este capítulo se dará
énfasis a las reacciones sólido-fluido, aún cuando su análisis también es aplicable a
reacciones líquido-gas o a dos líquidos inmiscibles entre sí.

3.1 MECANISMOS EN REACCIONES SOLIDO- FLUIDO

El sólido reaccionante pueden ser partículas individuales o aglomerado de

micropartículas (pellets). Además el sólido puede ser denso o poroso inicialmente;
puede formar un producto sólido denso o poroso que crece sobre el sólido
reaccionante, o formar un producto sólido que se desprende, se vaporiza o cambia de
fase, etc. Dependiendo entonces de las características físicas de los productos y
reactivos, se debe aplicar el modelo más apropiado.
Independiente del tipo de sólido, sin embargo, tres pueden ser los mecanismos que
controlan la trasformación en forma independiente o en combinación:

- Transporte de los reactivos y/o productos líquidos o gaseosos através del film
de fluido (film de concentración) que rodea la partícula.

- Difusión de los reactivos y/o productos líquidos o gaseosos a través del

producto sólido formado.

- Reacción química en la interfase reaccional o volumen reaccional

Existe un cuarto mecanismo que también puede controlar la transformación y es

la transferencia de calor desde o hacia el volumen o interfase reaccional. Este
76
 control es significativo solo en reacciones a alta temperatura, generalmente por
sobre los 800-1000°C y puede bajo ciertas condiciones ser el mecanismo de
control principal.
Para el caso más simple de un sólido denso A que reacciona con un fluido B líquido o
gaseoso para formar un sólido poroso C y otro fluido D (líquido o gaseoso) por ej.,
según la reacción del sólido A con el gas B:

A(s) + B(g) = C(s) + D(g) (3.1)

Si se supone que el producto C se forma progresivamente sobre el sólido A como una

capa de espesor creciente que se mueve hacia el interior del sólido inicial, los tres
primeros mecanismos antes mencionados pueden tomar parte en la transformación
(Figura 1)

º = cte.
B(f)

A (s) C(s)
m
B (f)
D(f)
Reacción
química en Difusión de Transporte de
la interfase B y D en la capa B y D en el film
de producto de fluído

Figura 3.1 Mecanismos de control posibles en reacciones sólido-fluido.

En ciertos casos, el control por transferencia de masa a través del film es poco
significativo, en particular, cuando la agitación del fluido es intensa como por ejemplo
en la lixiviación de un sólido por agitación, en las reacciones en lecho fluido sólido-gas,
en el flujo de gas en hornos de cuba a alto Reynolds, etc. En casos como éstos, los
controles pueden reducirse a dos: la difusión en la capa de sólido y/o reacción química
en la superficie reaccional. Cuando la difusión controla la transformación, la velocidad
de reacción química es muy alta comparada con la reacción en régimen puramente
difusional, y viceversa.

En la Fig. 3.2 se esquematica el efecto de dos tipos de control (químico y difusional) en

la velocidad global de reacción (rg).

77
 rq
rd
a) b)

rd
rq

rg rg
Ci
Ci C
C

Figura 3.2. Controles extremos: a) Control difusional, la velocidad de reacción en

régimen difusión << velocidad de reacción en régimen químico; b) control químico, la
velocidad reacción << velocidad de reacción en régimen de difusión.
( rg = vel. reacción global).

En todos los casos sin embargo las tres resistencias cinéticas actúan en serie y por
consiguiente la resistencia cinética total puede expresarse como:

1 1 1 1
= + +
RT rq rd rt

En que RT es la resistencia cinética total y rq, rd y rt las resistencias individuales de la

reacción química, difusión en la capa de producto y transporte de masa a través del
film de concentración, respectivamente.

Bajo ciertas condiciones se puede pasar de un régimen de control a otro durante la

transformación, por ejemplo, al oxidar MoS2 a MoO3 en aire a 560 – 600°C, en la
primera etapa y hasta cerca de un 90% del MoS2 transformado a MoO3, el régimen es
químico, seguido de uno mixto y finalmente por uno puramente difusional por sobre el
95-97% de transformación.

La definición del o los mecanismos depende también si el sólido es inicialmente poroso

o si se trata de un lecho de partículas de una cierta porosidad. Para una partícula
individual depende del tamaño de los poros del sólido que reacciona y/o del producto
que se forma sobre éste. Como se vió en el capitulo anterior, cuando el tamaño de los
poros es inferior al camino libre medio recorrido de las moléculas que difunden, la
difusión molecular debe ser reemplazada por la difusión de Knudsen. En la mayor
parte de los casos, sin embargo, una aproximación es suficiente, ya que el término
convectivo y la gradiante de presión interna en el sólido son pequeñas frente al
término puramente difusivo, es decir, se puede escribir para el reactivo gaseoso o
líquido B:
78
nB = −DBDef ∇CB (3.2)

Si los poros son muy pequeños, o la difusión de Knudsen es mucho menor que la
molecular DBDK << DBDef es decir DBDef DBDK y la expresión vista se transforma en:

nB −DBDK ∇CB (3.3)

El coeficiente efectivo de difusión de una especie líquida o gaseosa en un sólido se

puede estimar mediante la expresión de Bosanquet:

1 τ 1 1 
=  +  (3.4)
DBDef 
ε  DBD DBDK 

En esta expresión τ es el factor de tortuosidad de los poros y ε la porosidad media del

sólido. El valor de ε es relativamente simple de determinar, no así el de τ que tiene
un valor entre 2 y 10. Para sólidos tales como hematita, óxidos de metales no ferrosos,
sulfatos, etc. τ 4 – 6. Esta expresión es similar a la (2.15) anteriormente vista.

Para materiales a granel, tal como concentrados de minerales, el coeficiente de

difusión efectivo interpartícula puede estimarse de la Figura 3.3, conociendo el valor
de la porosidad global del lecho y el coeficiente de difusión molecular.

1.0

0.5
DAB ef / DAB

0.1
0.1 0.5 1.0
(cm 2 /cm 3 )
Figura 3.3 Relación entre la superficie específica ε de un lecho de partículas y el
coeficiente efectivo de difusión.

79
A diferencia de los fenómenos de transferencia sólido-fluido y de difusión, los
fenómenos de adsorción previos a la reacción química han sido poco estudiados en
sistemas metalúrgicos. Por lo general, se hace abstracción de ellos y las expresiones
cinéticas se consideran como reacciones simples, sin etapas intermedias. En la mayor
parte de los casos en metalurgia esta aproximación ha probado ser suficiente.

Para una reacción simple de un sólido A con un gas B que forman un sólido C según:

aA(s) + B(g) = C(s) (3.5)

La velocidad de reacción para el gas B y el sólido A se puede escribir como:

r&B = r&A = (1 a ) r&A = (1 a ) kr CmA CnB (3.6)

En la mayor parte de los sistemas heterogéneos sólido – fluido el orden de reacción

respecto del sólido es cero o cercano a cero, (esto no es siempre válido en sistemas en
que se forman aleaciones o soluciones sólidas o líquidas) en tanto que el orden con
respecto al gas o líquido en muchos casos es uno o cercano a uno. En tal caso, es válida
la simplificación:

r&B = (1 a ) r&A = ( 1 a ) kr CB (3.7)

Si se forma otro gas D en la transformación:

aA(S) + B(g) = C(s) + D(g)

(3.8)

La expresión ( 5 ) toma la forma

r&B = 1 a r&A = 1 a kr ( CmB − CnD / K ) (3.9)

Si tanto m como n son unitarios y el sistema es muy diluido o la reacción es

irreversible, CnD / K → 0 y la expresión (3.7) es suficiente para representar la cinética de
reacción.

Existen entonces distintas condiciones bajo los cuales pueden reaccionar un sólido y un
fluido, puesto que depende de las condiciones físicas de porosidad y densidad del
sólido reaccionante y del o los productos formados, la geometría de éstos, variación de
volumen, cambio físico, etc., de manera que es posible aplicar un modelo más
apropiado para cada caso específico. Sin embargo, existen dos condiciones extremas
características del sólido reaccionante bajo las cuales puede ocurrir una
transformación:

80
- La reacción ocurre en un frente reaccional (superficie) bien delimitada (Modelo
de núcleo no reaccionado o núcleo recesivo o topoquímico)

- La reacción ocurre en todo el volumen de la partícula (modelo de reacción

generalizada o reacción volumétrica)

El primer de estos modelos heterogéneos requiere de una velocidad de reacción alta y

un sólido reaccional denso, en tanto que el segundo requiere de un sólido reaccional
inicial poroso y una velocidad de reacción relativamente lenta.

Otros modelos particulares toman en cuenta otros factores físicos como la

fragmentación del sólido, el agrietamiento del producto, etc., en tanto que otros
modelos se pueden aplicar a la transformación de un conjunto de partículas (pellets)
tanto a nivel macro como microscopico (partículas individuales).

Para decidir cuál es el modelo más apropiado a emplear en cada caso, se debe iniciar
con ensayos de laboratorio bajo diferentes condiciones de temperatura, tiempo
reacción y condiciones de flujo del fluido y concentración de los reactantes, y
completar el estudio con cortes micro y macrográficos de partículas parcialmente
transformadas para así aproximar de la mejor forma el modelo cinético a los
resultados de laboratorio obtenidos.

En este capítulo se verán los casos extremos (reacción superficial y reacción

volumétrica), ya que en la mayor parte de las reacciones heterogéneas en metalurgia,
tanto en los casos sólido-liquido a baja o moderada temperatura como sólido-gas a
alta temperatura estos siguen aproximadamente el modelo de núcleo recesivo
llamando también modelo de núcleo sin reaccionar o modelo topoquímico.

3.1.1 Modelo heterogéneo de núcleo no reaccionado o núcleo recesivo

Este modelo de reacción requiere que el sólido que reacciona sea denso (caso común
en metalurgia) de tal forma que la reacción ocurre en un frente reaccional
relativamente bien definido.

Dentro de este modelo se pueden presentar dos sub-modelos:

- El producto sólido formado se desprende del núcleo sin reaccionar, o bien se

forma un producto líquido o gaseoso.

- El producto sólido formado crece en forma homotetica sobre el núcleo no

reaccionado.

La forma general de una expresión que permite predecir el grado de avance de la

transformación cuando hay un control mixto puede obtenerse haciendo ciertas
simplificaciones:

81
 - Las partículas iniciales son de forma geométrica regular (esfera, cilindro o
placa plana).

- La transformación es homotética es decir, las partículas no cambian de

forma ni de tamaño.

- El sistema está en estado seudo-estacionario, es decir, la velocidad de

reducción del núcleo es pequeña comparada con la velocidad de difusión
del fluido reaccionante a través de la capa de sólido formada sobre el
núcleo sin reaccionar, es decir, -dr/dt<<DBDef.

- La velocidad de difusión a través de la capa de producto es constante en el

tiempo, es decir, independiente de rn (coeficiente de difusión constante).

- El sólido inicial es denso y el reactivo no difunde a través de éste.

- La reacción es irreversible, de primer orden respecto al fluido y orden cero

respecto al sólido.

En estas condiciones, los flujos de reactivos y productos (liquido o gases) son

conservativos y las velocidades de transferencia en el film, difusión en la capa de
producto sólido y de reacción química en la interfase son iguales.

Las ecuaciones que expresan la transformación global se pueden colocar en función de

la desaparición de B con el tiempo o la variación del radio de sólido A o la velocidad
másica de reacción del sólido A con el tiempo. Para una reacción general como la
considerada anteriormente:

A) NUCLEO RECESIVO SIN FORMACION DE CAPA DE PRODUCTO

Este modelo se puede aplicar en los casos en que un sólido A reacciona con un fluido B
para formar otro fluido, o bien se forma otro sólido que se desprende del núcleo. Un
ejemplo de ello es la combustión de carbón con desprendimiento de la ceniza formada
y generación de productos gaseosos (CO y CO2):

aA (s) + B(f) = C(f) (3.10)

A (s) + B(f) = C(s) (C se desprende del núcleo de A)

(3.11)

Para simplificar el modelo se puede considerar que el sistema es isotérmico y en

régimen seudo-permanente. De esta forma, en la superficie externa de la partícula de
sólido A se tiene que el flujo de B desde el seno del fluido vale:

.
N B = (1 a ) hm ( CBo - CnB ) (3.12)

82
El subíndice n se refiere a la concentración de B en interfase reaccional sobre el núcleo
no reaccionado y CBo es la concentración del fluido B en el seno de éste. Al inicio de la
reacción rBn = rBs en que rBs (ó rBo ) es el tamaño inicial de la partícula.

En régimen permanente el flujo NBo es igual a la velocidad de reacción química en la

interfase reaccional, suponiendo una reacción irreversible de orden cero respecto del
sólido y orden uno respecto del fluido. Por lo tanto:

N& B = 1
a krCnB (por unidad de superficie) (3.13)

de donde:

 h  0o
CnB =  m  CCBB
 hm + kr 

Entonces:

 hk  1 1
NB =  m r  CCBo ==CC Bo// +  (3.14)
 hm + kr   kr hmm 

En esta expresión, si kr << hm se tiene que Cnb −~ CBo (el control por transferencia en el
film es despreciable) y la transformación puede representarse con buena aproximación
sólo por la expresión para el control por reacción química:

NB (1 a ) kr CBo (3.15)

Lo contrario ocurre si el control es por transferencia en el film de fluido, ya que kr >>

hm y CnB → CB∗ de manera que:

N& B (1 a ) hm (CCBoB − CC*BB

donde CB∗ es la concentración de equilibrio para la reacción química. Si CB∗ → 0 esta

última expresión se puede colocar como:

N& B ( 1 a ) hm CBo

83
Cuando tanto el control por transferencia en el film como la reacción química
participan en la transformación, vale la expresión (14).

Por otra parte, si la partícula del sólido A se considera una esfera de radio r, la
variación de tamaño de la partícula con el tiempo está relacionada con el flujo de B
según:

ρ& dr
N& = A (3.16)
a dt

En que ρ& A es la densidad molar del sólido A. De esta forma, se tiene de (3.14) y (3.16)
que

dr   11 11 
− = = (a /ρρ& A ) C Bo / ++  (3.17)
dt  hhm kk r 

En la expresión (3.17) todos los términos son independientes del radio de la partícula
excepto hm el que se puede estimar mediante alguna de las expresiones vistas
anteriormente en el curso.

Para los casos de control por reacción química y control por transferencia, la variación
del radio de una partícula con el tiempo se puede deducir de estas expresiones y
esquematizar simplificadamente como en la Figura 3.4.

Film de concentración
del fluído

A A

CºB
CBº
CsB = CBn
a) b)

CsB = CBn C*B

C*B
0 rn m 0 rn m

Figura 3.4 Modelo de núcleo recesivo (núcleo sin reaccionar) sin formación de capa
de producto sólido adherido al núcleo o formación de productos líquidos o gaseosos:
a) Control por transferencia en el film de fluido
b) Control por reacción química en la interfase reaccional

84
a) Control por reacción química

En este caso, si se considera que la reacción es irreversible de primer orden con

respecto al fluido y cero respecto al sólido, se puede escribir para una reacción como la
(3.10) ó (3.11):

1 dnA ρ& A dr adnB

− = =− = ak r C Bo (3.18)
s dt a dt sdt

En esta expresión la constante de velocidad está expresada en cm/seg., s es la

superficie reaccional del sólido A (superficie externa inicial) y ρ& A la densidad molar de
A. De esta forma para una partícula esférica de radio inicial ro la velocidad de reacción
de la especie A se puede escribir:

dr
= a kr CBo / ρ& A (3.19)
dt

Las condiciones de borde son t = 0, r = ro y t = t, r = r. De esta forma, se tiene que para

alcanzar un radio cualquiera r el tiempo t requerido será:

ρ& A
t= (ro − r ) (3.20)
ak r CBo

y el tiempo total τ de transformación para r = 0 será:

ρ Aro
τ= (3.21)
ak r CBo

La fracción de sólido XA transformado cuando el núcleo tiene un cierto tamaño r se

puede expresar como una relación entre los volúmenes inicial del sólido y cuando el
sólido alcanza un radio r, es decir:

VA − VA 4 πro3 − 4 3 πr 3
XA = o r
= 3

VA 4 πr 3
o 3 o

X A = 1 − (r / ro )3 (3.22)

85
b) Control por transferencia en el film

La forma de la expresión para el coeficiente de transferencia de masa hm depende de la

relación empleada. Por ejemplo, para partículas pequeñas en régimen de Stokes, la
expresión de Ranz y Marshall antes vista se reduce a

2 DBD DBD
hm = (3.23)
d r

Como el control es sólo por transferencia de masa en el film, reemplazando la

expresión (3.23) en la expresión (3.17) y resolviendo para las CB t = o, r = ro y t = t, r = r
se encuentra que el tiempo para alcanzar un radio r del núcleo sin transformar será:

t=
ρ& Aro2
2aC BoDBD
[
1 − (r / ro )2 ] (3.24)

y el tiempo total de transformación τ para r = 0 valdrá:

ρ& A ro2
τ= (3.25)
2a CB° DBD

La fracción transformada de A para un valor cualquiera r del radio del núcleo es igual
que para control químico, es decir,

X A = 1 − (r / ro )3 (3.26)

B) NÚCLEO RECESIVO CON FORMACIÓN DE CAPA DE SÓLIDO DENSO

Este caso ocurre con relativamente mayor frecuencia en metalurgia y puede tomar
tres formas límites dependiendo del tipo de control particular: por transferencia en el
film de fluido, químico en la interfase reaccional o difusivo a través del producto sólido
formado, o bien una combinación de estos. (Fig. 3.5).

Este modelo requiere que el (o los) productos formados se adhieran al núcleo de sólido
reaccionante formando una capa de producto de espesor creciente.

86
 A Film A A
fluido
C C C

CBº CºB
CBº
CsB n CsB
CB
a) b) c)

CsB
CnB CnB
C*B C*B C*B
0 rn rs m 0 rn rs m 0 rn rs m

Figura 3.5. Formas de control límite en el modelo de núcleo recesivo con

formación de capa de producto sólido:
a) Control por transferencia en el film b) Control por reacción química
en la interfase c) Control por difusión en la capa de producto sólido formado.

Si se tiene un control mixto (como el de la Figura 3.6) todas las resistencias son
significativas en la transformación:

rn rs
A (s) Film
fluido
C(s)
º
CB

s
CB

n
CB

C*B
0 rn rs m

Figura 3.6. Mecanismos mixtos de control en una partícula densa con formación de
capa sólida sobre el núcleo recesivo. (transporte en el film de fluido + difusión en la
capa de producto sólido + reacción química en la interfase reaccional).

Un ejemplo de esta condición la representa la sulfatación de óxido cúprico con SO2 a

650°C en presencia de oxigeno, la que se puede representar mediante un modelo de
estas características. La reacción que ocurre es

C u O(s) + S O2(g ) + 0.5 O2(g) = C u S O4(s) (3.27)

87
Experiencias efectuadas con pellets esféricos de 1 cm diámetro (formados por un gran
número de partículas individuales) demostró que el modelo de núcleo recesivo de CuO
con formación de capa de producto sólido (CuSO4) representaba bien el fenómeno, a
nivel macroscópico como puede verse en la micrografía de un corte de pellet
parcialmente reaccionado (Fig. 7).

Figura 3.7. Corte transversal a través de un pellet de óxido cúprico parcialmente

reaccionado con SO2.

Para una reacción general como:

aA(s) + B(f) = C(s) + D(f)

si se expresa esta transformación en función de la variación de la concentración del

fuido B con el tiempo (variación del numero de moles de B), se puede escribir para
cada mecanismo individual que:

- Transferencia en el film de fluido:

= 4 π rs2 hm (CCBoB0−−CCBsBs)
dnB
− (3.28)
dt

Donde CBs es la concentración del gas en la superficie externa del sólido

- Difusión en la capa de sólido:

dnB dC
− = 4 π rn2 DBDef B r =rn (3.29)
dt dr

- Reacción química en la interfase reaccional:

dnB
− = 4 π rn2 kr C Bn (3.30)
dt

88
Donde C Bn es la concentración del fluido B en la interfase reaccional del núcleo no
reaccionado de sólido A.
Para evaluar la gradiante de CB en la expresión (3.29), se puede hacer un balance de
materia para un elemento ∆r del sólido formado C considerando el proceso semi-
estacionario, es decir, la difusión de B en C, es independiente de drn/dt. Como el flujo
en conservativo (no hay reacción en el interior del producto sólido C), se puede
escribir:

æ dC ö çç4 p r2 dC B ÷
æ ö
DBDef çç4 p r2 B ÷ ÷ - DBDef ÷ =0 (3.31)
çè dr ø÷ r çè dr ø÷ r+ D r

De donde se tiene que para D r ® 0 :

d 2 dCB 
 r DBDef =0 (3.32)
dr  dt 

Esta ecuación puede igualmente deducirse del balance general de masa para una
esfera:

∂ 1 ∂ ∂C
( εC B ) = 2  DBDef r2 B  − r&B ρ& A = 0 (3.33)
∂t r ∂r  ∂t 

Si se considera régimen seudo permanente, la reacción sólo ocurre para r=rn y DBDef no
varía con rn el término del miembro izquierdo y el 2º término del miembro derecho
son cero, por lo que se llega igualmente a la expresión (3.32).

Esta ecuación diferencial (3.32) puede resolverse directamente por doble integración
con las CB:

I) r = rs C B = C Bs

II) r = rn C B = C Bn

De esta forma se obtiene:

 1 − (rn / r ) 
CB - C nB = ( C sB − C nB )   (3.34)
 1 − (rn / rs ) 

Para obtener ahora la expresión para el gradiante de CB con r se puede diferenciar esta
última expresión con respecto al radio y luego evaluar para r = rn:

dC B C BS − C nB
= (3.35)
rn (1 − rn / rs )
r =rn
dt

89
Esta expresión puede ahora ser reemplazada en la expresión (3.29) para obtener así:
dn  C s − C nB 
− B = 4 π rn DBDef  B  (3.36)
dt  1 − rn / rs 

Se puede obtener entonces Cnb en términos de CBo y rn eliminando CBS y (dnB/dt) en las
ecuaciones (3.28) y (3.30), obteniendo de esta forma:

CBo
C nB = (3.37)
  rn  2 k r kr  rn 
1 +   + rn  1 − 
  r  hm DBDef  rs 

La velocidad global de transformación de la partícula (en términos del consumo de B o

de A) puede entonces expresarse reemplazando el valor de CnB en una expresión
adecuada, por ejemplo en la reacción para control por reacción química.
Reemplazando entonces (3.37) en (3.30) se tiene que la velocidad de reacción del
fluido B vale:

dnB 4 πrn2 CBo k r

− = (3.38)
dt   r 2 k kr  rn 
1 +  n  r + r n  1 − 
  rs  hm DBDef  rs 

En general, es más importante conocer como varía el sólido A en función del tiempo.
La variación de A con t se puede expresar en función de la variación del radio del
núcleo del sólido A con el tiempo. Por estequeometría de la reacción se tiene:

dnB 1 dnA 1 ρA dvn 1 pA d 4 3

= = = ( πrn )
dt a dt a PMA dt a PMA dt 3
(3.39)

Entonces para una partícula esférica:

dnB 1 dnA 4 π rn2ρA  drn 

= =   (3.40)
dt a dt a PMA  dt 

Reemplazando esta última relación en la expresión para control químico (3.30) se tiene
entonces que:

drn a PMA kr n
=− CB (3.41)
dt ρA

90
Reemplazando la expresión (3.37) encontrada para CnB en esta última expresión, la
variación del radio del núcleo de sólido A no reaccionado con el tiempo toma la forma:

drn aPMAkr CBo / ρ A

− = (3.42)
dt   r 2 k k k r  r 
1 +  n  r + r + r n  1 − n 
  rs  hm hm DBDef  rs 

o bien

drn aCBo / ρ& A

− = (3.43)
dt  r 2 1  r − r  k r r 1 
 n  +  s n  r n n + 
 s  m  DBDef  DBDef rs k r 
r h

en que ρ& A es la densidad molar de A.

Las expresiones diferenciales (3.36), (3.41) ó (3.42) dan la solución para rn=f(t) en
o
términos de CB y t. La solución de esas ecuaciones diferenciales dependerá de cómo
varía CBo con t. Por ejemplo, en un reactor perfectamente mezclado CBo es constante,
o
(como un reactor de lecho fluidizado o un estanque de lixiviación), en cambio CB varía
con t para reactores de flujo pistón, como en el caso de un reactor de reducción
vertical (cuba de alto horno, sistema Midrex etc.) o en un sistema de lixiviación por
percolación en columnas, etc.

Cuando la concentración del fluido en el seno de éste es constante, las ecuaciones

diferenciales (3.42) ó (3.43) se pueden resolver directamente con las CB siguientes:

I) t = 0, r = ro = rs

II) t = t, r = rn

De esta forma, se tiene:

a k r o   rn  k r  rn 
r 2
kr
− C B ∫ 1 +   + rn  1 + drn (3.44)
ρ& A ro   rs  hm DBD ef  rs 


91
Luego:

é æ é æ ö3 ù ì 2
ö ïïü 1 æ rn ÷ öùú
ρ& 1 1
ê1 ç 1 ro ö÷
÷ê rn ÷ ú ro ïïï ærn ÷
t =t =A ror o 0
ê çç - ç
÷ê1 - çç ÷ ú+
ç
í 1 - çç ÷ ç
÷ ýï + k çç1 - r ÷ ú (3.45)
a a CBCB ÷ê èro ÷
ê3 èçhm DBDef ÷
ø ÷
ø ú 2DBDef ïï èro ÷
øï rè
÷
÷
o øú
ëê ë û ïî þï ûú
El tiempo de transformación total τ cuando rn = 0 será entonces:

ρ& Aro  1 r 1
τ= 
o 
+ o +  (3.46)
aC B  3hm 6DBD ef k r 

Esta última relación es general, y de ella se pueden deducir las otras vistas
anteriormente. Cuando el régimen es químico kr << hm y kr << DBDef se tiene:

ro ρ& A
τq = (3.47)
CBo k r

Para régimen difusional (DBdef<<kr y DBdef<<hm):

ro2 ρ& A
τd = (3.48)
CBo k r

y para régimen de transporte en el film hm<<DBdef y hm << kr:

ro ρ& A
τt = (3.49)
3 a C Bohm

La fracción convertida del sólido A cuando r = rn será entonces igual a:

3
r 
x A = 1 −  n  = 1 − (rn / ro ) 3 (3.50)
 rs 

En la deducción de estas expresiones se consideró la reacción irreversible, con un

tamaño constante de la partícula. Si estas condiciones no se cumplen, aún es posible
deducir expresiones analíticas, pero las soluciones son más complejas. Considerando
una distribución de tamaño discreta de fracciones de tamaño desde d1 a di para cada
fracción de tamaño la fracción del sólido transformada será:

3
r 
X Ai = 1 −  n  (3.51)
 rs i

y el tiempo total τ requerido para cada fracción de tamaño en llegar a rni serán
respectivamente:
92
 ρ& A roi  1 r 1
τi = º 
+ oi +  (3.52)
a CB  3 hm 6 DefBD kr 

ρ& A roi  1  1 1    rn      r 3 1  r  
3
roi r
ti = º
  + +  1 −     + oi 1 −   +  1 −  
n n
(3.53)
a C B  3  hm DB Def kr    ro    2 DBDef r kr  ro i 
  i    i  

Luego, la fracción promedio transformada del sólido A para un tiempo de reacción t

será la sumatoria de la fracción transformada para cada fracción de tamaño de
partícula por su fracción en peso Wi correspondiente a la distribución de tamaño
considerada:

i   r 3 
X A = ∑ 1 −  n  Wi (3.54)
n=1  r
  s i 

Si el tiempo de reacción t es igual, por ej., al tiempo total de transformación τ de la

fracción mayor, obviamente todas las otras fracciones menores estarán
completamente transformadas. En la práctica, un reactor opera siempre con t < τ para
las fracciones mayores, de manera que X A < 1.

Para cualquier condición sin embargo, en una reacción heterogénea se observa que el
tiempo total de transformación vale:

ttotal = tfilm + treacción + tdifusión

Y para el tiempo total de conversión, éste será la suma de los tiempos individuales
totales:

t total = t film + t reacción + t difusión

Para la combinación de resistencia (transferencia en el film, difusional en la capa y

reacción química), la ecuación general tiene la forma vista (3.22.a) o (3.22.b).

Para otras formas geométricas (esfera, cilindro y placa plana), en la Tabla 3.1 se
resumen las expresiones derivadas para casos límite de control individual, para
partículas de tamaño constante y partículas de tamaño decreciente.

La mayor dificultad que se presenta con los modelos de núcleo no reaccionado con
superficie reaccional es que sólo en contados casos el sólido inicial es tan denso y la
reacción química lo suficientemente rápida o el producto sólido formado lo
suficientemente poroso como para que se forme un frente claramente definido. En la
mayoría de los casos, evidencia directa de partículas parcialmente transformadas

93
 indican que más que una superficie reaccional existe un frente de pequeño espesor
con una cierta gradiente de concentración, semejante al frente formado en el modelo
de reacción generalizado en que la velocidad de reacción química es alta.

Tabla 3.1 Expresiones para la conversión de partículas de acuerdo al modelo de núcleo recesivo
Control por
Control por Difusión Control por Reacción
Transferencia
en la capa de sólido química
en el film
Placa plana t t = XA t t = X 2A t t = XA
e  e = espesor 
Partículas tamaño constante

XA = 1−   t = rˆ A ahmCCBBo o
t = rˆ AL2 2aDefBDCCB0B o
t = rˆ AL ak rCCBB
L  L = longitud 
Cilindro 0.5
2 t t = XA t t = X A + (1 - X A )ln(1 - X A ) t t = 1 - (1 - X A )
 rn 
X A = 1 -   t = rˆ Ar0 2ahmCCBBo τ = ρˆ Ar 02 4aDefBDCCBo0 t = rˆ Ar0 akrCCBBo
 ro 
Esfera 2 1
3 t t = XA t t = 1 - 3(1 - X A ) 3 + 2(1 - X A ) t t = 1 - (1 - X A ) 3
 rn 
  t = rˆ Ar0 3ahmCCBo t = rˆ Ar 02 6aDefBD CCBo0
XA = 1- t = rˆ Ar0 akrCCBo
 ro 
2 1
Esfera
Part. disminuyen tamaño con

t t = 1 - (1 - X A ) 3 t t = 1 - (1 - X A ) 3
sólido que se desprende

Part. pequeña (No se aplica)

o
reg. Stokes t = rˆ Ar 02 2aDCC B0 t = rˆ Ar0 ak r C BoB

0.5
t t = 1 - (1 - X A ) 1
Esfera t t = 1 - (1 - X A ) 3
32 o (No se aplica)
Part. grande t =r /C
0 BB t = rˆ Ar0 ak r C B0
K = cte.

Para las partículas que disminuyen de tamaño (esfera pequeña), el coeficiente de difusión Dji es el
molecular en el film (y no el efectivo Defkj), de tal manera que se asemeja a uno difusional puro.

3.1.2. Modelo de transformación generalizado o modelo de reacción volumétrico

En este caso, a diferencia del modelo de reacción superficial, la transformación ocurre

simultáneamente en toda la partícula, lo que implica que ésta debe ser inicialmente
porosa (Def grande) y la reacción química relativamente lenta (kr pequeño) y que el
producto formado sea igualmente poroso o que este sea otro fluido, como en caso
de la combustión de coke que genera CO2 o CO.

En el caso extremo, si kr < < DefBD , prácticamente no existe gradiente de

concentración de B en el sólido y la reacción ocurre efectivamente uniformemente en
todo el volúmen de la partícula. Por el contrario, si kr>>DefBD todo el fluido B que
penetra al interior de la partícula se consume de inmediato y la reacción se restringe a
un frente estrecho que avanza hacia el interior de la partícula, es decir, representa el
caso del modelo de núcleo recesivo. En la práctica siempre existe una cierta gradiente
de concentración del fluido hacia el interior del sólido, así como también una gradiente

94
de concentración del sólido reaccionante, (Figura 3.6) puesto que Def no es infinito y/o
k → 0.

Figura 3.8 Transformación de un sólido poroso en todo su volumen cuando hay

formación de producto sólido.

A) GRADIENTE DE CONCENTRACION DE UN FLUIDO QUE REACCIONA EN EL

INTERIOR DE UN SÓLIDO POROSO.

Para el flujo unidireccional de un fluido B en el interior de un sólido poroso (Figura

3.9), la variación del flujo de B causado por la reacción química dentro del sólido se
puede expresar en la forma:

dN&
r&B = − B
∆x
dx

A su vez, el flujo de fluido en el interior del sólido poroso queda dado por la primera
ley de Fick:

& = −D S dCB
NB efBD
dx

Es decir:

& = −D dy dx dCB
NB efBD
dx

Si se considera el coeficiente efectivo de difusión constante con la temperatura y la

composición, se puede escribir que:
95
 &B d2 CB
−
dN
= D efBD (dy dx )
dx dx 2

Si ocurre una reacción química en el interior del sólido en que se alcanza CB = CB*
donde CB* > 0,

−
dN&
B
= DefBD
d2 CB − CB* (
(dy dx ) ) (3.55)
dx dx 2

Por otra parte, la expresión cinética se puede escribir en general en la forma vista

− r&B = −r&A = 1
4 k r S A (CB − CB*

donde C* es la concentración de B en el equilibrio.

NB V
y
N B - dN B
GB G B - dG B
V
z

V
x
Sólido poroso
Figura 3.9 Variación del flujo de un gas reactante en un elemento de volumen de
sólido poroso.

Llamado Ŝ A la superficie específica del sólido A, en (cm2/cm3) se tiene entonces que:

( )
r&B = ( 1a )Ŝ A k r CB − CB* dz dy dz

kr tiene las dimensiones de (cm/seg) en este caso. Reemplazando en (3.55), se tiene

que

d2 (CB - CB* )
DefBD 2
= (1 a ) Sˆ A kr (CB - C*B ) (3.56)
dx

O bien

d2 (CB - CB* ) CB - C*B

- =0
dx 2 H

en que

96
 DefBD
H=
( a ) Sˆ A kr
1

La ecuación (3.56) se puede colocar en general como:

(
∇ 2 CB − CB* − ) 1
H
( )
CB − CB* = 0 (3.57)

Se pueden presentar varios casos, según la geometría de la partícula:

a) Partícula esférica

En este caso el flujo es radial y por simetría esférica y asimutal la expresión anterior
toma la forma:

d2 (CB - CB* ) 2 d(CB - C*B ) 1

2
+ - (CB - C*B )= 0 (3.58)
dr r dr H

Las CB son en este caso,

Cs = CBo CB = C*B

r = ro r=0

La solución de esta última expresión con estas CB tiene la forma

K1 K
CB - C*B = senh(r H)+ 2 cosh (r H)
r r

Las constantes K1 y K2 se evalúan de las CB, resultando así:

ro (CBo − C *B )
k1 =
senh(ro / H))

K2 = 0

La solución general para la partícula esférica toma entonces la forma:

CB − C *B ro senh (r / H)
= (3.59)
C Bo − C*B rsenh (ro / H)

b) Partícula plana:

97
Llamando eo el espesor en el que ocurre transformación (simultánea desde ambas
caras) y despreciando la superficie correspondiente a los bordes de ésta, la solución de
la ecuación (3.57) tiene la forma:

CB - C*B = K3 exp(- x / H) + K4 exp(+ x / H)

con las CB:

x=e x=0
dCB
CB = C*B =0
dx

En estas condiciones, las constantes de esta expresión valen:

C Bo − C *B
K3 = K 4 =
cosh(e o / h))

y la expresión general toma entonces la forma

CB − C *B cosh(x / H)
= (3.60)
C B − CB cosh(e o / H)
o *

La forma de los perfiles de concentración del fluído en el interior de un sólido poroso

en función del parámetro H para una esfera y partícula plana se puede ver en la Figura
3.10.

ro
eo /2

H =1 H =1
0.5 0. 5
0.1 0.1
5 5
0.02 0.0
- -2
3
3

10
10 -

10
10 -

0 ro 0 eo /2
Figura 3.10 Perfiles de concentración del fluido reactante en el interior de un sólido
poroso, cuando la reacción es generalizada en el total del volumen de la partícula.

3.1.3 Modelos generales de aglomerados de partículas (pellets)

98
Existe una diversidad de modelos propuestos. Algunos simplifican el problema al
hacer abstracción de la geometría de los poros, otros consideran poros cilíndricos
dispuestos al azar, ó simétricos, grietas, etc. (Fig. 3.11).

a) b) c)

Figura 3.11 Modelo de poros

Un caso especial lo representa el modelo de pellets compuesto de un agregado de

partículas esféricas individuales, las cuales a su vez pueden seguir el modelo de núcleo
recesivo o modelo de reacción generalizado (Figura 3.12).

Figura 3.12 Modelos de pellet de agregados de partículas esféricas

3.2 FENOMENOS DE NUCLEACION

Fenómenos de este tipo ocurren por ej. en reacciones sólido-gas o cuando el sólido por
efecto de la temperatura sufre una transformación cristalina o bien cuando hay
descomposiciones térmicas. Esta transformación puede o no generar a su vez un
producto sólido, líquido o gaseoso, según:

A (s) € C(s)
A (s) € C(s) + D(g)

Existen dos modelos para explicar el fenómeno, ambos postulan el crecimiento del
producto sólido formado desde determinados puntos iniciales o gérmenes de
nucleación. Estos modelos no están restringidos solo a descomposiciones térmicas, y
se pueden aplicar a otros fenómenos de reacciones químicas.
En general, se considera que las reacciones sólido-gas en particular se inician en los
sólidos en algunos puntos de este los cuales tienen mayor energía superficial. En estos

99
puntos crecen los gérmenes iniciales de la reacción que luego se propagan a la
totalidad de la masa o de la superficie.

3.2.1 Modelo de Prout – Tomkin

En este modelo se considera que al inicio de la transformación se forman un cierto

número inicial de núcleos de crecimiento los cuales van creciendo lentamente a lo
largo de los planos de fracturas de la partícula y en la superficie, es decir, el
crecimiento de los núcleos ocurre en dos dimensiones (Fig. 3.13).

Núcleos C
iniciales de C
A A
t=0 t>0

Figura 3.13 Modelo de Prout- Tomkin de crecimiento de gérmenes de nucleación.

Puesto que al inicio la interfase entre reactivo y producto es pequeña, el crecimiento

inicial es lento, para luego acelerarse y finalmente decrecer nuevamente cuando los
planos del producto se entrecurzan. De ésta forma, la transformación es autocatalítica
y se puede expresar en la forma:

ln éëX A (1 - X A )ùû= Mt - N (3.61)

En que M es la pendiente de la curva al graficar ln éëX A (1 - X A )ùû versus t y N se

obtienen del valor de ti que es la intersección de la recta obtenida para
ln éëX A (1 - X A )ùû= 0 (Fig. 3.13).
ln X / (1 - XA )

0
ti
t

Figura 3.14 Modelo de Prout-Tomkin

3.2.2. Modelo de Avrami

100
Este modelo postula que el crecimiento del producto ocurre a partir de ciertos
gérmenes (núcleos) iniciales los que se forman continuamente. El modelo se puede
aplicar tanto a crecimiento lineal (1 dimensión) como a superficial (2 dimensiones) o
volumétrico (3 dimensiones) (Fig. 3.13).

t=0 t>0 t >> 0

Núcleos
iniciales

Figura 3.15 Modelo de Avrami

El crecimiento de la fase formada a partir de los núcleos iniciales se puede expresar en

la forma:

N
X A = 1 - e- Mt (3.62)a

o bien,

ln éëX A (1 - X A )ùû= MtN (3.62)b

é æ 1 öù
Gráficamente, M y N pueden determinarse graficando ln êêlnçç ÷
÷ú versus ln t (Fig.
çè1 - X ÷÷úú
ëê A øû

3.16).

N
(1 -1XA)

0
ln ln

t=M Int

Figura 3.16 Modelo de Avrami

En la expresión analítica de Avramí, N depende de cómo y cuando se forman los

núcleos del producto, como se observa en la Tabla 3.2.
Tabla 3.2 Valores de N en el modelo de Avrami
Producto crece en Producto crece en Producto crece en
101
 forma forma 2 3 dimensiones
unidimensional dimensiones (esferas)
(agujas) (placas planas)
Todos los núcleos
se crean N=1 N=2 N=3
inicialmente.
Se crean núcleos
continuamente, a
N=2 N=3 N=4
velocidad
constante

Ambos modelos, de Prout – Tomkin y de Avrami son similares. Más aún, predicen
conversiones similares en condiciones iguales.

3.3 DETERMINACIÓN DE LAS ETAPAS QUE CONTROLAN UNA TRANSFORMACIÓN

La determinación de la o las etapas que participan en la transformación de un sólido

debe hacerse a partir de datos experimentales. El tamaño de las partículas, efecto de
la temperatura y tiempo requerido para la transformación pueden dar una indicación
de o los mecanismos involucrados.

a) Tamaño de partículas

De las expresiones vistas se pueden observar que para partículas de igual tamaño, el
tiempo de reacción es proporcional al tamaño inicial de estas (ro por ejemplo para
esferas). De esta forma, se tiene que para partículas esféricas:

- Control por transferencia en el film: t µ ron (n varía entre 1.0 y 2.0,

decreciente al aumentar el valor de Re).

- Control por reacción química: t µ ro

- Control por difusión a través del producto sólido: t µ ro2

Por consiguiente, haciendo experiencias con distintos tamaños uniformes de

partículas, se puede deducir el tipo de control prevalenciente.

b) Temperatura

Tanto la difusión del fluido a través de la capa de producto sólido formado como la
transferencia en el film varían poco con la temperatura, (en un rango moderado de
temperatura), no así la cinética química, lo que es más sensible a una variación de
temperatura. De la ecuación de Arrhenius se observa que la constante específica de
reacción kr es función exponencial de la temperatura (expresión 1.26):

102
k r = k ro exp− (e / RT)

En que E es la energía de activación y k ro el factor de frecuencia.

De esta forma, al variar la temperatura si quien controla la transformación es la

reacción química se debe cumplir la relación de Arrhenius, es decir, se debe obtener
una recta al graficar la kr versus 1/T, además de obtener un valor de la Energía de
Activación entre 30 y 100 kcal/mol.

c) Tiempo de reacción

Una forma directa de determinar la (s) etapas (s) de control se puede obtener al trazar
las curvas de avance de la transformación en función del tiempo adimensional
reducido (Fig. 3.17). Se observa que si el control es químico, la transformación es
lineal con el tiempo para el modelo de núcleo no reaccionado. Si el control es
difusional en la capa de sólido formado, ésta es rápida al comienzo (cuando no hay
capa formada o es delgada) para decrecer luego a un mínimo y acelerarse nuevamente
hacia el término de la transformación debido a la rápida disminución de la superficie
reaccional respecto de la velocidad de difusión del reactivo hacia la superficie
reaccional. Finalmente si el control es por transferencia, ésta aumentará
progresivamente con el aumento de la transformación debido a la disminución del
tamaño de las partículas o del núcleo sin reaccionar.
1.0
Control dif. capa
prod. (part. plana)
0.999

0.99
Control reacc. quim.
(part. plana y esfera) 0.97
= ro ó e o
e

0.9
XA
r

0.5
0.8
Control dif. capa
prod. (esfera)
0.6

0.4
Control trans. film
(part. plana y esfera)
0.2

0.0 0.0
0.0 0.5 1.0
= t/

Figura 3.17 Formas de las curvas de avance de la transformación de un sólido de

tamaño constante.

Si se grafican las resistencias a la transformación correspondiente a los tres

mecanismos, (ver por ejemplo ecuación 3.52), se tiene para el tiempo total de
transformación:

103
 ro ρˆ A  1  1  ro  1  1 
τ=    +   +   (3.63)
aCBo  3  hm  6  D efBD   kr 

Se puede ver que en régimen químico el control es constante; decreciente para el

control por transferencia en el film, y con valor máximo y para el control difusional. En
la Figura 3.18 se observan los valores de cada resistencia para un pellet denso de
hematita, reducida en presencia de hidrógeno a 800ºC.

3.5
Resistencia total
3.0
Resistencia (seg/cm)

2.5

2.0
Capa sólido
1.5 Transf. film

1.0
Reacción química
0.5

0.0
0.0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0
=t/
Figura 3.18 Variación de las resistencias individuales en la transformación de una
partícula densa (pellet de hematita reducido con H2 a 800ºC).

Como la resistencia total es superior al valor de las resistencias individuales, la

velocidad de reacción global debe ser igual o inferior al valor individual de cada
mecanismo, como se observa en la Figura 3.19, en la que se ha trazado la velocidad de
reacción hipotética de un sólido, y en que la velocidad de reacción global se compara
con las que se obtendrían si solo un mecanismo controlara la transformación. Tanto la
velocidad de transformación en régimen difusional como en régimen de transferencia
son menos afectadas por la temperatura, no sí en régimen químico.

De los estudios mencionados (variación del flujo de fluido, variación de la temperatura

y variación del tamaño del sólido), la interpretación de ellos, con excepción de casos
en que ocurren cambios de estado como fusión, volatilización, condensación, etc., se
puede resumir en la Tabla 3.3.

104
 o
mic
f. film
Trans

Quí
Difusión sólido

r Curva experimental

0
T
Figura 3.19 Velocidad de reacción global e individual en cada mecanismo de control.

Tabla 3.3 Test para determinar la etapa de control de reacciones sólido fluido.
Resultados esperado para
Test Control Control Control
Químico difusional transferencia
Aumento del flujo Inalterado Inalterado Aumenta la
de fluido velocidad de
reacción

Variación de la E = 30-100 E = 5-20 kcal/mol E = 20-50 kcal/mol

temperatura Kcal/mol Puede ser
válida la Válida una
expresión de expresión tipo
Arrhenius Arrehenius

Variación del a) t α tamaño (sólido denso) t α (tamaño 2) t α (tamaño)n

tamaño de las b) t ∼ constante (sólido denso) 1<n<2
partículas

BIBLIOGRAFÍA CAPITULO III

1. J. Szekely, J. W. Evans and H. Y. Sohn, “Gas – Solid reactions”, Academic Press,

(1976).

2. L- Coudurier, D. W. Hopkins and I. Wilkomirsky, “Fundamentals of the Metallurgical

Processes” 2nd Ed, Pergamon Press, (1986).

105
 CAPITULO IV

DISEÑO DE REACTORES
4.0 CONCEPTOS EN EL DISEÑO DE REACTORES

4.1 REACTORES IDEALES

El volumen físico en que ocurre una transformación o reacción se denomina en general

reactor y corresponde al espacio útil o neto y que puede o no ser igual al volumen
total del reactor. Por otra parte reactor es una palabra genérica que define el equipo
en que ocurre la transformación o reacción química y puede tener cualquier geometría
o forma.

En este volumen neto o volumen de reactor VR las reacciones que ocurren entre los
reactivos requieren de una mezcla de ellos, así como de un tiempo de permanencia y
condiciones de concentración, presión y temperatura determinadas en su interior.

Las propiedades o condiciones intensivas de concentración y temperatura son

independientes de la geometría o características del reactor, sin embargo, el tiempo
de retención y patrón de flujo en su interior no lo son.

Existen tres tipos extremos de reactores en cuanto a mezcla interna y tiempo de

permanencia de los reactivos en su interior:

- Reactor discontinuo (batch)

- Reactor de flujo pistón (sin mezcla axial)
- Reactor de mezcla perfecta o mezcla total

Tipos de reactores empleados en metalurgia se esquematizan en la Figura 4.1.

En reactores que operan a alta temperatura tal como un convertidor al oxigeno para
refinar acero (1600oC) o un horno de fusión flash (1300ºC), la cinética es muy rápida y
por tanto la mezcla entre sólido (arrabio a refinar o concentrado a fundir) y el gas no es
primordial en el diseño del horno, no así en un estanque agitado de lixiviación donde la
aproximación al equilibrio es lenta y la mezcla de éste es más importante.

En cuanto al tiempo de reacción (o tiempo de residencia en el reactor), se pueden

clasificar los reactores en dos grandes grupos:

- reactores discontínuos o por carga (batch)

- reactores continuos

El primer tipo de reactor al igual que en un reactor batch, la reacción ocurre en un

tiempo bien definido t y la concentración de los reactivos decrece desde su valor inicial
hasta la de salida, que puede ser la de equilibrio.
106
 S G G
S
G
G+S

G G G G
S
Flash
G
S
G
S

G G

S Rotativo
S G
G
Lecho móvil G
en contracorriente G
S
S
G G
S S

G+S
G G+S
S S

S
Reactor ciclónico (fundido) Pisos múltiples
S S

G+S S L G

S S S

G+S
G
Lecho fluído
(S en MP, G en FP G
S+G Sinter
en ciertos casos)
Estanque agitado (flujo cruzado
(S y en MP) S en lecho estático)
G

S G

L L
(fundido) (fundido, L
L
escoria) (fundido) (fundido)
Reverbero
Conox refinación Horno
acero con oxígeno eléctrico

Figura 4.1 Tipos de reactores metalúrgicos sólido – líquido, sólido – gas, líquido –
líquido y sólido – líquido –gas.

107
El segundo tipo de reactor se puede subdividir en dos dependiendo
fundamentalmente de la mezcla que ocurra en el interior del reactor: si no hay
mezcla en absoluto (flujo pistón) todos los reactivos en la salida habrán permanecido
en el reactor un tiempo preciso igual a t de manera que este tipo de reactor es
equivalente al reactor batch. En el otro extremo de comportamiento está el reactor
de mezcla total o mezcla perfecta en que cada elemento que entra al reactor se
mezcla instantáneamente y totalmente con el resto de la masa. Esto implica que
aparecerán reactivos sin reaccionar en la descarga del reactor puesto que todos los
elementos en el interior de éste tienen igual posibilidad de salir en la descarga.

Para los reactores contínuos, los balances de masa y energía (macroscópicos) se

pueden escribir en general en la forma vista anteriormente para las ecuaciones de
continuidad:

- Masa: (4.1)

 Tasa de acumulació n   Flujo másico total   Flujomásic o total   Tasa de reación o 

       
 de los reactivos en  =  de entrada al  −  de salida  +  desaparici ón de los 
 el reactor   reactor   del reactor   reactivos en el reactor 
       

- Energía:

 
 Acumulación de   Total energía   Generación o consumo   Total energía que 
       
 energía en  =  que entra  +  de energía − sale del reactor, 
 el reactor       incluyendo pérdidas 
   al reactor   en el reactor   
 de energía alambiente 
(4.2)

4.2 REACTORES DISCONTINUOS (BATCH)

Para reactores discontinuos también llamados reactores batch o reactores por carga,
los reactivos se cargan en este el que luego se presuriza, y/o aumenta la temperatura
hasta alcanzar las condiciones requeridas y luego se hacen reaccionar durante un
tiempo definido para luego descargarlos y comenzar otro ciclo de operación.

La ecuación general de conservación de masa para una especie A que reacciona para
generar otra B será en un reactor batch es por lo tanto:

 Tasa de acumulación   Tasa de desaparición 

   
 de B en el reactor  =  de A en el reactor 
   
   
(4.3)

Esta última expresión puede escribirse también para el reactivo A en la forma:

108
 dNA dX A
− = −NoA = −R
& =R
A
&
B (4.4)
dt dt

Donde NoA indica el número total inicial de moles del reactivo A en el reactor, XA la
fracción molar de A reaccionada y RA la tasa global de reacción de A en el reactor.

Por otra parte, la tasa global de reacción de A en el reactor se puede expresar como:

−R
& = ( −r& )V
A A R (4.4)

Donde VR es el volumen útil o volumen neto del reactor ocupado por reactivos y
productos en la reacción y −r&A es la velocidad de reacción por unidad de volumen de
reactor. De esta forma, se tiene:

dX A
− NoA = (− r&A )VR (4.5)
dt

ecuación que se puede resolver con las CB:

t = t1 t = t2

X A = X1A X A = X2A

Por lo tanto:

X 2A dX A
t 2 − t 1 = NºA ∫ 1 (4.6)
XA (− r&A )VR
La expresión (4.6) representa la forma más general de diseño de un reactor batch, y
puede simplificarse si VR = cte. Por otra parte −r&A por lo general varía con t y sólo en
contados casos es independiente del tiempo. De esta forma; la expresión (4.6) se
puede escribir como:

X 2A
NoA dX A
t = t 2 − t1 =
VR ∫ − r& (4.7)
X1A A

Por otra parte, si al comienzo de la reacción X1A = 0 es decir, no hay transformación

previa del reactante A, como frecuentemente ocurre, se tiene que:

XA
NoA dX A
t=
VR ∫ (− r& )
o A
(4.8)

109
Si el volumen neto de reacción VR es contante (densidad global del sistema constante)
y como X A = 1 – C A / CoA y NoA / VR = CoA la expresión (4.8) puede escribirse
entonces en la forma:

XA C
dX A A
dC A
t = C oA ∫ = −∫ (4.9)
o
(− r&A ) CoA (− r&A )

Si el volumen neto de reacción VR varía con el tiempo, por ejemplo en forma lineal, se
puede escribir:

VR = VRo (1 + ξX A ) (4.10)

Donde VRo es el volumen inicial (útil) del reactor y ξ el coeficiente de variación de

volumen de la mezcla reactivos más productos (cero inicialmente). El valor de ξ puede
ser + o – dependiendo si el volumen total se contrae o expande durante la reacción. De
esta forma, la ecuación (4.9) se puede escribir:

XA
dX A
t =N ∫ (− r& )(1 + ξX )
o
A (4.11)a
o A A

En todas estas expresiones, el término cinético −r&A debe ser representado en función
de kr y CA. El valor de −r&A es un valor experimental y como se vió anteriormente debe
ser determinado en forma experimental, tanto para el valor de la constante específica
de velocidad kr como para el orden de reacción de CA.

4.3 REACTORES IDEALES CONTINUOS

4.3.1 REACTORES DE FLUJO PISTÓN

El reactor de flujo pistón, también llamado reactor tubular, considera que no existe
mezcla axial (longitudinal) en éste y que no hay gradiente radial en el reactor. De esta
forma, el reactor se puede imaginar como un volumen reaccional en el cual entra en
un extremo un frente continuo de reactivos ya mezclados los que avanzan
ordenadamente a lo largo del reactor reaccionando entre sí hasta emerger en la
descarga, sin mezcla axial interna y después de permanecer todos y cada uno de los
elementos de reactivos un tiempo t único e igual para cada elemento en el reactor.
(Fig. 4.2). Se observa que un reactor de flujo pistón perfecto es equivalente a un
reactor batch puesto que el tiempo de permanencia de todas las entidades en su
interior es uno solo y para el reactor FP dependerá solo de la tasa de alimentación de
A y B y del volumen del reactor.

110
 dXA

FA
º uo
Productos
FB
º
dVR
L

Figura 4.2 Esquema de un reactor de flujo pistón ideal.

Reactores de este tipo existen sólo en forma aproximada puesto que siempre existe un
cierto grado de mezcla de los reactivos axial dentro del reactor debido a la difusión
molecular de reactivos y productos y perfil de velocidad en el interior del reactor, sin
embargo, en algunos casos es posible aproximar un reactor a uno de tipo pistón, como
es por ejemplo el flujo del sólido descendente en la cuba del alto horno.

Si se llama FA la tasa molar de alimentación de A que entra a un elemento de volumen

dVR de reactor (Fig. 4.2), se puede establecer un balance de masa para un elemento de
volumen de reactor dVR en la forma siguiente:

FA − (FA + dFA ) = ( −r&A ) dVR (4.12)

Como FA = FAo (1 − X A ) en que XA es la fracción convertida de A, se tiene que:

dFA = −FAo ( dX A ) (4.13)

Donde FAo es la tasa molar de alimentación del reactante A a la entrada del reactor,
considerada constante.

De (4.12) y (4.13) se tiene entonces:

dVR dX A
= (4.14)
FAº − r&A

Esta última expresión puede ser integrada para el total del volumen útil del reactor.
Las CB pueden ser dos valores X1A y X2A cualquiera. De esta manera, la solución de la
expresión (4.14) toma la forma general:

VR X2A
dX A
∫V =F ∫ (− r& )
o
R A (4.15)
o X1A A

111
VR X2A dX

FAo ∫X1A (− r&A )

= A
(4.15)a

Si la reacción es completa e irreversible (t = τ) es decir, los reactivos se convierten

totalmente en productos y el reactante A entra al reactor sin haber reaccionado, se
puede resolver la expresión (4.15)a para t = 0 X1A = 0 y t = τ X 2A = 1 . De esta
forma, se tiene

1
VR dX A
=∫ (4.15)b
FA 0 (− r&A )
o

Y en general,

VR X A dX
=∫ A
(4.15)c
º
FA 0 (− r&A )
Estas últimas dos expresiones son útiles para diseñar un reactor de flujo pistón (o
cercano a este comportamiento) y requiere conocer la relación entre la fracción
convertida de A y el recíproco de la velocidad de reacción −r&A . Si el sistema es
isotérmico, esta relación generalmente es simple, pero no así si el sistema no es
isotérmico.

Si se conoce el volumen VR de un reactor y el flujo de alimentación FAo , para un grado

de conversión XA es posible calcular −r&A lo que puede hacerse para distintos valores de
XA. De esta forma es posible establecerse una ecuación para la velocidad de reacción
−r&A en términos de XA y tener la ecuación (4.15) en forma apropiada para el diseño del
reactor.

La expresión (4.15) no contiene el tiempo de reacción o tiempo de residencia t en el

reactor. Si no existe contracción o expansión del volumen entre reactivos y
productos (o se puede despreciar), el tiempo de retención en un reactor de flujo
pistón se puede escribir para el total de la alimentación en la forma:

VR WR L
t= o
= o = (4.16)
v F uo

Donde:

t = Tiempo de retención o reacción en el reactor (t)

uo = Velocidad espacial (lineal) en el reactor (Lt-1)
L = Longitud del reactor (L)
vo = Tasa de alimentación volumétrica al reactor (L3t-1)

112
WR = Masa total en el reactor en régimen permanente (M)
Fo = Tasa de alimentación molar o másica al reactor (Mt-1)

La expresión (4.15) puede entonces también colocarse en términos del tiempo de

reacción anteriormente definido, por ejemplo, para el componente A:

VR
Como t = , se puede escribir:
vo

t V V
o
= o R o = oR
C A v C A FA

Reemplazando en la ec. (4.15) se tiene:

X2A dX A
t = CoA ∫ 1 (4.17)a
XA (− r&A )

Y si X oA = 0 en la entrada del reactor:

XA dX A
t = C oA ∫ (4.17)b
o (− r&A )

Se puede observar que ésta última expresión es igual a la expresión (4.9) para un
reactor batch. Esto es lógico puesto que en ambos reactores todos los elementos o
entidades en su interior permanecen un tiempo único igual a t.

Si la temperatura es constante a lo largo del reactor, el término kr es constante en la

expresión de velocidad r&A y puede ser sacado de la integral y resolver en forma
analítica las expresiones (4.15) o (4.17). En forma más general, la integral en estas
expresiones puede ser resuelta en forma gráfica o numérica, como se muestra en la
Figura 4.3. En ella, la forma de la curva debe ser determinada previamente en forma
experimental para distintos valores de XA de manera que la integral bajo la curva
valdrá VR / FAo .

113
 Area = t = VR / FºA Area = t = CºA VR / FºA
(-1/ rA )

(-1/ rA )
0 XAº XA 0 CA CAº

Figura 4.3 Determinación gráfica de la integral de las ecuaciones (4.15) y (4.17) para
un reactor de flujo pistón.

La expresión (4.17) puede también ser escrita en la forma:

CA dC A
t=∫0 (4.18)
CA (− r&A )
Si la reacción es de primer orden con respecto al reactante A y orden cero respecto de
los otros reactantes, se puede colocar entonces:

CA dC A 1
t=∫ = ln (C 0A / C A ) (4.19)
C0A k r C A kr

La definición del tiempo de residencia en el reactor en la forma descrita en la ecuación

(4.16) es válida sólo si no hay variación del volumen útil del reactor (contracción o
aumento de volumen), de tal manera que el tiempo se puede definir como el espacio-
tiempo y en rigor el tiempo de reacción o residencia en el reactor se define como la
integral:

VR
dVR
t= ∫
o vo
(4.20)

En esta expresión, v o es la tasa volumétrica de alimentación de los reactivos en el

interior del reactor, la cual es variable punto a punto y que vale a su vez:

v = v º (1 + ξ X A ) (4.21)

Donde ξ es el coeficiente de variación de volumen de la suma de los reactivos y

productos en el reactor y v o la tasa volumétrica de alimentación de los reactivos en la
entrada del reactor.
114
Si los reactantes no varían de densidad, o lo que es lo mismo, si el total de volumen de
reactantes más productos es constante, ξ =0, v o es constante y dVR = 0, siendo valida
entonces la definición del tiempo de residencia dado por la relación (4.16).

4.3.2 REACTORES PERFECTAMENTE MEZCLADOS O MEZCLA TOTAL

En los reactores de mezcla perfecta (MP) también denominados de retromezcla, a

diferencia del reactor de flujo pistón, la mezcla es total (perfecta), lo que implica que
cada elemento o entidad de fluido o sólido que entra al reactor se mezcla
completamente con el resto en un tiempo infinitamente corto, de tal manera que la
composición de la descarga es por definición igual a la composición en el interior del
reactor. Por esta razón, también se les denomina reactores de retromezcla total o de
mezcla completa.

Puesto que por definición todas las entidades (fluido, átomos, partículas, etc.) tienen
igual posibilidad estadística de salir por la descarga, habrá una distribución normal de
tiempos de residencia o de reacción en el reactor: algunos elementos permanecerán
un tiempo muy corto en tanto que otros estarán un largo tiempo en él y una fracción
mayoritaria permanecerá un tiempo cercano al tiempo medio de residencia o
espacio-tiempo, definido ahora como:

WR VR
t= = 0 (4.22)
F0 v

Para reactores perfectamente mezclados t = L / uo no tiene significado físico ya que L

no es un parámetro caractrístico y uo es un campo de velocidad y no un valor único.

Este tiempo medio o espacio tiempo sin embargo, no tiene un significado físico
exacto como en el caso de un reactor batch o de flujo pistón, puesto que cada entidad
en el interior del reactor tiene un valor de tiempo de reacción (residencia) real
diferente a las demás entidades.

Frecuencia
F Aº F Bº
C Aº C Bº
Reactor
-r A MP Desc. CB, CB
-r B VR FA = FºA (1-XA )
CA WR
CB FB = FºB (1-XB)
t>t t t>t

Figura 4.4. El reactor ideal perfectamente mezclado o de mezcla total

115
Reactores perfectamente mezclados ideales no existen, puesto que se requiere de un
tiempo finito para alcanzar la mezcla total de sus elementos en su interior, pero
existen mucho casos en que debido a la intensa agitación en el interior del reactor, la
mezcla puede aproximarse a una mezcla ideal, como es por ejemplo un estanque
agitado de lixiviación de partículas finas, un reactor del lecho fluidizado (condición
para el sólido), etc.

En un reactor perfectamente mezclado en régimen permanente se puede hacer un

balance de masa para un reactante cualquiera A para el total del volumen útil del
reactor en la forma siguiente:

FAº = FA + (− r&A )VR (4.23)

Debe notarse que r&A está evaluada para las condiciones internas del reactor (o
descarga) puesto que son idénticas, pero no para FAo . La ecuación (4.23) también se
puede escribir en la forma

FA = FAo (1 − X A ) + (− rA )VR (4.24)

en que X oA es la fracción transformada del reactivo A en la entrada, normalmente = 0

De esta forma, se puede escribir:

VR X A − X oA
= (4.25)
FAo (− r&A )

Si X oA = 0 (el reactante entra al reactor sin ningún grado de conversión):

VR X
= A (4.26)
FA (− r&A )
o

Estas dos expresiones (4.25) y (4.26) pueden ser empleadas directamente para el
diseño de reactores perfectamente mezclados, si se conocen tres de los términos.
Normalmente, el término ( −r&A ) se debe determinar en forma experimental de un
estudio de laboratorio o planta piloto; XA es el que económicamente dicte el grado de
conversión deseado y FAo es un parámetro de diseño, normalmente el que dimensiona
la capacidad de diseño de la planta. Esto permite conocer VR que es el tamaño del
reactor requerido.

Adicional a ello, esta ecuación permite conocer −r&A directamente de un reactor dado si
se conoce VR y FAo para distintos valores de la fracción transformada de XA. Se observa

116
que ambas expresiones son independientes de la geometría del reactor y dependen
solo del volumen de éste.

Estas expresiones pueden también ser expresadas en función del tiempo medio de
residencia o espacio-tiempo t anteriormente definido en la forma
t = VR / v º = WR / Fº .

Puesto que se tiene:

VR VR X A − X oA
= = (4.27)
FAo C oA v o (− r&A )

Se puede escribir entonces:

X 2A

∫ dX
o
C
C o (X − X oA )
A A
V X1A
t = Ro = A A = (4.28)
v (− r&A ) (− r&A )

o bien, si X oA = 0

X 2A

∫ dX
o
C A A
V X1A C oA X A
t = Ro = = (4.29)
v (− r&A ) (− r&A )

Si se analizan las expresiones de diseño para un reactor batch y uno de flujo pistón y se
comparan con las de un reactor de mezcla perfecta (ec. 4.9 y 4.17a y ec. 4.29) se
observa que el término r&A en el caso del reactor de mezcla perfecta está fuera de la
integral y está evaluado en la descarga del reactor, (puesto que ésta es igual en
composición que la mezcla en el interior del reactor), a diferencia del reactor batch
en que −r&A varia con el tiempo, o bien − r&A varía punto a punto en el interior de éste
como para un reactor de flujo pistón.

En igual forma para un reactor de mezcla total cuando hay variación de volumen en el
reactor se debe expresar el tiempo medio en la forma vista anteriormente.

VR

∫ dV R

t= o
(4.30)
v o (1 + ξ A X A )

Es importante notar que en el interior de un reactor de mezcla perfecta hay una

distribución de tiempos de residencia de todas las entidades (como se verá más en
detalle más adelante del curso) de tal forma que puesto que todas las entidades en la

117
descarga del reactor tienen entre sí distintos tiempos de permanencia en el reactor,
habrá algunas totalmente convertidas en tanto que otras sólo estarán parcialmente
convertidas y algunas muy poco.

De esta forma, la composición de un determinado producto en la descarga del

reactor de MP es el promedio estadístico de conversión de todas las entidades de
acuerdo con su grado de conversión individual y su ponderación en peso o volúmen.
Para efectos de cálculo, sin embargo, el valor del tiempo medio de reacción es valor de
referencia y se emplea para diseño.

En igual forma que para el reactor de FP, es posible calcular en forma gráfica la razón
VR / FAo y por lo tanto el volumen del reactor requerido (Fig. 4.5). Para el reactor de MP
esta razón la representa el rectángulo de área (X 1
A − X2A ) x ( 1 −r&A ) puesto que r&A es
constante y no varía punto a punto, como ocurre con un reactor FP (Fig. 4.5)

Cinética
del reactor
FP descarga
y reactor MP (vR / FºA ) MP
(-1/ rA )

i ca
ét
ci n (vR / FºA ) FP
rva
Cu
Cinética
reactor
FP (entrada)
0 1
XA 2
XA
XA

Figura 4.5. Determinación del área VR / FAo para un reactor de mezcla perfecta (MP) y
uno de flujo pistón (FP).

Se observa que para el área VR / FAo de ambos reactores (FP y MP) se tiene siempre que

 VR  V 
 o  >  oR  (4.30)
 FA MP  FA FP

Para un flujo molar FAo dado (capacidad de una planta, por ejemplo), la relación de
volúmenes de reactor en cada caso será entonces:

(VR )MP > (VR )FP y (t )

MP ()
> t FP (4.31)a

Esto significa que para obtener un cierto grado de transformación XA en un reactor, si

este opera en régimen de flujo pistón se minimiza su volumen y por lo tanto la
inversión, en tanto que la inversión es mayor para el reactor que opera en régimen de

118
mezcla perfecta debido a que requiere un tamaño de reactor significativamente mayor
para tener el mismo grado de conversión que el de flujo pistón, como se observa en la
Figura 4.5.

Se debe tener presente, sin embargo, que este solo aspecto, (tamaño del reactor) no
es suficiente como para definir el tipo de reactor más apropiado. Factores tales como
características de los reactivos, tamaño de las partículas sólidas, densidad, generación
o consumo de calor, forma y requerimientos de agitación de la masa reactante
simplicidad de diseño, simplicidad de construcción y de operación así como
confiabilidad del diseño empleado juegan también un papel importante en la selección
y diseño definitivo del reactor a emplear.

4.3.3 COMPARACIÓN ENTRE REACTORES DE MEZCLA PERFECTA Y FLUJO PISTÓN.

Se pueden visualizar las condiciones de operación de un reactor de FP y otro de MP al

graficar la velocidad global de reacción en el reactor en función de la concentración del
reactivo en el reactor (Figura 4.6).

Area = t = CºA V R / FºA

Cinétic a de (reactor MP)
salida de
los reactores
FP y MP
(-1/ rA )

Cu
rv
a
ci
né
ti c
a

Cinética de
entrada de
reactor FP
CA º
CA
CA
Area = t = CºA VR / FºA
(reactor FP)

Figura 4.6 Comparación entre la cinética de un reactor de mezcla total y otro de flujo
pistón.

En la Figura 4.6 se puede ver que la composición de CA en la descarga del reactor FP es

igual a la interna (o descarga) del reactor de MP. Esto indica que la cinética global es
más lenta en reactor MP y explica entonces porque es necesario emplear un reactor
de mayor tamaño (mayor tiempo de retención) si se desea obtener el mismo grado de
conversión en un reactor MP respecto de un reactor de FP.

Se puede también observar que cuando el grado de transformación no es alto, la

razón entre los volúmenes del reactor de MP respecto al de FP crece lentamente.
Para una reacción de un sólido A que reacciona con un fluido B, si la reacción es de
primer orden respecto a B y orden cero respecto al sólido, se tiene la curva de la Figura
4.7.

119
 7

FP
/ (VR)
MP
4

= (VR )

1
0 0.5 1.0
XA

Figura 4.7 Razón entre los volúmenes de un reactor MP y FP para una reacción de
orden 1 respecto al fluido y orden cero respecto del sólido A.

Se observa en la Figura 4.7 que para valores de XA sobre 0.7, aproximadamente, la

razón de volúmenes crece rápidamente. Por ejemplo, para XA =0.9, (VR)MP/(VR)FP=4.6,
mientras que para XA=0.95 la razón vale 6.9 y para X=0.99 esta vale cerca de 14, es
decir, se requiere un reactor casi 14 veces mayor si éste opera un régimen de mezcla
total a que si opera en régimen sin mezcla para una conversión de 99%.

Si se emplean reactores de flujo pistón en serie la conversión total en la descarga del

último de ellos es aditiva, al igual que el tiempo medio total de retención:

VR1 VR2 VRn VRT

t T = t1 + t2 + ... + tn = 0
+ 0
+ ... + 0
= 0
(v o = cte)
v 1 v 2 v n v

Esto también es cierto para una serie de reactores de mezcla perfecta en serie, sin
embargo, al emplear reactores de mezcla perfecta en serie es posible reducir el
volumen total y si se extrapola hasta un número infinito de reactores en serie, estos
se comportarán como un solo reactor de flujo pistón ideal, como se puede observar
en la Figura 4.8 para 1 reactor de MP, 2, 3 y 5 reactores de MP en serie. Se puede
observar en la Figura 4.8 que si emplean 5 reactores de MP en serie el volumen total
de estos es sólo cerca del 55% del volumen que si se emplea un solo reactor MP y
aproximadamente 12% mayor que si se emplea un reactor FP. Sin embargo, se debe
observar que a medida que aumenta el número de reactores de MP en serie también
aumenta el costo de capital, el que puede crecer rápidamente. Por ej., si es un tren de
reactores agitados en serie, se requieren 5 agitadores, 5 bombas, etc. lo cual puede ser
más alto en costo de capital que un solo reactor.

120
 1 Reactor MP 2 Reactores MP
1
A 1 = VR / FºA

A 2 = VR2 / FºA 1
A 1 = VR / FºA

-1/ rA

-1/ rA
XA XA

3 Reactores MP 5 Reactores MP
1
A 1 = VR / FºA

2
A 2 = VR / FºA

-1/ rA
-1/ rA

3
A 3 = VR / FºA

XA XA

Figura 4.8. Reactores de mezcla perfecta en serie comparados con un reactor de flujo
pistón.

Si todos los reactores tienen igual volumen, para los 5 reactores de la Figura 4.8 se
tiene que los tiempos medios de retención en cada uno son iguales y entonces:

t1 = t2 = t3 = t4 = t5
t T = 5t1
VRT = 5VR1

Se puede observar también de la Figura 4.8 que a medida que aumenta el número de
reactores en serie el volumen total de estos se aproxima al de un reactor de flujo
 ∞ 
pistón y cuando n→ → ∞,  ∑ VR  = ( VR )FP .
 1 MP
En general, si se considera una reacción irreversible (K muy grande) de un fluido A +
con un sólido B con orden n para el fluido reactivo A y orden cero para el reactivo
sólido B y además se considera que X oA = 0 y ξ tiene un valor cualquiera, positivo o
negativo, se puede escribir entonces en general que:

k (C o ) (1 − X A )n
n

(−rA ) = r A
(1 + ξ A X A )

Para esta reacción, si el reactor opera como reactor flujo pistón, de la ecuación (4.17)
se tiene que:

121
 VR 
 o  = ∫
XA (1 + ξA XA )
n

dX A (4.32)
 FA FP
o
( )
k r C oA (1 − X A )
n n

En tanto que para el reactor operando como reactor de mezcla perfecta, de la

ecuación (4.28) se puede escribir:

n
 VR  XA  1 + ξA XA 
 o  =   (4.33)
( )
 FA MP k r C oA
n
 1− XA 

De (4.32) y (4.33) se tiene entonces que la razón ψ entre volúmenes de los reactores
DP y MP vale:

 VR (CoA ) 
( )
n

 FA MP
o
XA 1+ ξ A X A n
1− X A (VR )MP
ψ= = = (4.34)
 VR (CoA ) 
n X A (1+ ξX )n
(VR )FP
 FA FP
o ∫ o
A
(1− X A ) dX A

El valor del parámetro Ψ se ha calculado para distintos valores de n y ξ y valores

constantes de FAo a la entrada del reactor (FP ó MP). Esto se encuentra graficado en la
Figura 4.9. Esta figura es válida para un reactante cualquiera (para efectos de orden
de reacción) puesto que se ha supuesto el orden de reacción igual a cero para el otro
(u otros), por ejemplo:

A →B →C r&A = k r CnA
A +B → C+D r&A = k r CnA CmB (m = 0)
A +B +C → C+D r&A = k r CnA CmB CpC (m, p = 0 )

Si para ambos reactores de flujo pistón y mezcla perfecta se tiene que las tasas de
alimentación másica y concentración de entrada de los reactantes a los reactores son
iguales, (FAo ) = (FAo ) y ( CoA ) = ( CoA ) . Se puede entonces establecer una
MP FP MP FP

relación entre los volúmenes de los dos reactores en la forma de:

( VR )MP = ψ ( VR )FP

Como se puede observar en la Figura 4.9 (al igual que en la Figura 4.8), la razón Ψcrece
rápidamente cuando la conversión es alta.

122
 10 2

n = orden de reacción
n=2
=2

n=1 =1
=2
=0
=1

10 =0
= -0.5
n = 0.5 = -0.67
=1

n = 0.25
=1
1
10 -2 10 -1 1
1 - XA = CA / CºA

Figura 4.9. Valores del parámetro Ψ en función del orden de reacción y coeficientes
de variación de volumen del reactor.

4.3.4 REACTORES PERFECTAMENTE MEZCLADO EN SERIE

Como se vió anteriormente, una forma de obtener una reducción del volumen total
de un reactor de mezcla perfecta es empleando varios reactores más pequeños en
serie, cada uno perfectamente mezclado y normalmente de igual volumen entre sí.

Reactores perfectamente mezclados en serie se emplean en forma relativamente

extensa, en particular, en sistemas de lixiviación sólido – líquido para maximizar la
extracción con una menor inversión de capital y operación.

Para reactores en serie, para un reactor cualquiera i se puede escribir, de acuerdo a la

relación (4.27):

VRi X iA − X iA−1
= (4.36)
FAº (− rA )i

Donde ( −r&A )i es la velocidad de reacción en el reactor i constante para ese reactor. En

igual forma, se puede escribir para el reactor i de acuerdo a la relación (4.28) que:

VRi
ti = º =
(
C ºA X iA − X iA−1 ) (4.37)
v (− r&A )i
Estas dos expresiones (4.36) y (4.37) son aplicables a cualquier tipo y número de
reactores en serie de igual tamaño.

Puesto que para el reactor de mezcla perfecta i se tiene también que:

123
C iA = C oA (1 − X iA ), X iA = 1 − C iA / C oA y X iA-1 = 1 − C iA−1 / C oA (4.38)

VRi CiA−1 − CiA

ti = =
vo (− r&A )i

y para el tiempo medio total en la serie de reactores:

n
t total = ∑ t i
i=1

4.3.5 ANÁLISIS GRÁFICO DE REACTORES PERFECTAMENTE AGITADOS EN SERIE

Frecuentemente se encuentra el problema de conocer la composición (y por lo tanto

la eficiencia del sistema) con la que salen los productos en reactores de mezcla
perfecta en serie, o bien y más comunmente determinar cual es el número óptimo de
reactores de mezcla perfecta en serie a emplear para alcanzar un determinado valor
de conversión.

Como se indicó anteriormente, la forma de la curva cinética de la reacción puede ser

obtenida de información de laboratorio o planta piloto. En general, se puede obtener
una expresión de la forma:

(− r&A ) = f (C A )
Donde f(CA) es una función de una forma cualquiera conocida. Por ejemplo, si se tienen
tres reactores perfectamente mezclados en serie pero de volúmenes diferentes
VR1 ,VR2 y VR3 ,(Fig. 4.10) y se conoce la concentración CoA en la entrada del reactor 1, para
conocer la concentración de A en la descarga del reactor V3 se puede colocar la
expresión (4.38) en forma general para un reactor cualquier en serie:

i−1
− CiA
(− rA )i = C A (4.39)
ti

Al graficar − (r&A )i versus CiA en esta expresión se obtiene una recta de pendiente
( −1 / ti ) de origen en la abscisa en CiA y que corta la curva cinética en un punto el cual
su vertical a la vez corta la absisa en CiA−1

Para el sistema de tres reactores perfectamente mezclados en serie de distintos

volúmenes de la Fig. 4.10, para el primer reactor se tiene que:

C ºA − C1A
(− rA )i = (4.40)
t 1i

124
 FºA
CºA , º

3
CA
CºA = C1A 2
CA
1 2 V3R t3
VR VR
t1 t2 3
CA
º

1
CA
2
CA

Figura 4.10 Tres reactores perfectamente mezclados en serie de diferentes

volúmenes.

Para un valor dado de CoA esta última expresión representa una recta de pendiente
( −1 / t1 ) de origen en la abscisa en CoA y que corta la curva cinética en el punto A que
representa la cinética en el reactor Nº1 y cuyo valor C1A es la composición en el interior
del reactor 1 y por lo tanto es la composición también de alimentación al reactor Nº2.
De esta forma, ser repite el procedimiento pero ahora con el valor C1A como origen en
las absisas. La pendiente de la recta será ( −1 / t2 ) , puesto que:

(− r&A )2 = C A − C A
1 2

t2

La cual corta la curva cinética en B y cuya vertical corresponde a la composición del

reactor 2 e igual a C 2A y que es la composición de alimentación al reactor 3.

Lo mismo será válido para el reactor Nº3 y cualquier otro número de reactores en
serie. Gráficamente, esto puede observarse en la Figura 4.11.

125
 Curva cinética de
reactor FP ó batch
Cinética A
reactor 1

Cinética B
reactor 2
Cinética C
(- rA )
reactor 3
-(1/ t 1 )

-(1/ t 2 )

-(1/ t 3 )

0 3
CA 2
CA 1
CA CAº
CA

Figura 4.11 Determinación de la composición de salida de un sistema de tres

reactores perfectamente mezclados en serie de distintos volúmenes.

Para trazar las rectas CoA − A, C1A − B y C2A − C, el valor del tiempo medio de reacción t
se puede obtener en cada caso de las relaciones:

t 1 = VR1 / v o , t 2 = VR2 / v o , t 3 = VR3 / v o

Donde vo se considera constante (ξ −

~ 0 ) y los volúmenes de los reactores VR1 , VR2 y
VR3 conocidos.

Cuando los volúmenes de los reactores son todos iguales, que es el caso común,
VR1 = VR2 = VR3 se tiene que t1 = t2 = t3 y la pendiente de las rectas es la misma, lo que
simplifica la determinación de las rectas.

Si no se conoce el número de reactores a emplear y solo se conoce el valor de la

composición de entrada CoA y la composición final CA que se desea obtener en la
descarga de estos (ó XA), se puede ver que cualquier número de reactores de mezcla
perfecta en serie desde uno hasta n puede satisfacer la condición de descarga CA. Si
se fija a priori el número de reactores a emplear, el problema se limita entonces a
determinar el volumen de ellos (se consideran todos de igual volumen).

Por ejemplo, en la Figura 4.12, se tiene una curva cinética conocida y los valores de
entrada CoA y el valor deseado de salida CA. Se puede ver que si se fija n = 3 (línea

126
segmentada), se debe ajustar al valor de la pendiente ( −1 / t ) para que se cumpla:
AB + BC + CD = CoA - C A

1
CA H´
I´
J´ E
CA2 K´
F
L´
(-r Aº ) M´
3
CA

G·

D A
CA G´ CF´ E´ D´ B C´ B´ CºA
CA
Figura 4.12 Determinación gráfica del volumen de cada reactor perfectamente
mezclado en serie (3 reactor ------ y 7 reactores ____).

Se observa que tgα = BE / BA = FC / BC = etc.....y como se tiene que

V 1
C −Cº 1
t1 = R
= A A
, se puede calcular el volumen de cada reactor puesto que se
vº (− r&A )
conoce la pendiente de cada recta ( −1 / t ) y la tasa volumétrica de alimentación vo
que se supone constante (ξ=0). Ajustando en forma gráfica el valor de α o de
( −1 / t ) se obtiene la condición requerida de AB + BC + CD = CoA − CA .
Si se requieren, por ejemplo, 7 reactores (línea llena, Fig. 4.12), se repite el
procedimiento, ajustando ahora el valor de β de tal forma que se cumpla:

AB`+ B`C`+ C`D`+ D`E`+ E`F`+ F`G`+ G`D = C oA - C A

En igual forma, para reactores de igual volumen, el volumen de cada reactor se puede
calcular a partir de la pendiente de cada recta ( −1 / t ) y de la relación ti = VRi / v o .

4.3.6 ASPECTOS ECONÓMICOS DEL DISEÑO DE REACTORES

Si bien el rendimiento máximo de un reactor es una de las condiciones básicas para el

diseño, éste valor debe ser compatibilizado con los costos de capital y operación
asociados a ello.

127
Por ejemplo, si se considera un proceso en el cual se extrae oro mediante cianuración
en un reactor agitado, y se supone que la extracción depende solo del tiempo medio
de residencia de las partículas en el interior del reactor, un mayor tiempo de retención
significa una mayor extracción de oro y por ende mayores entradas por ventas del
metal, sin embargo, también aumenta el costo de capital y de operación al aumentar
el tamaño o número de reactores de tal forma que habrá un punto de equilibrio
pasado el cual un aumento del tiempo de residencia (y por consiguiente del tamaño o
número de reactores) no producirá mayor incremento del aumento de entrada por
ventas debido a la mayor extracción del metal. Esto se puede observar en la Figura
4.15 en la cual se han trazado curvas típicas de rendimiento versus costo total, (costo
de capital + costo de operación).

Esto puede también aplicarse a un sistema de reactores en serie, donde al aumentar el

número de reactores aumenta el rendimiento del proceso, pero también aumenta el
costo de capital y de operación de la planta.

$ Aumento de rendimiento
$ Aumento del costo total

(entradas incrementales)
(Costos incrementales)

(capital operación)

iento
d e rendim
ent o
A um l
o tota
t
cos
e n to
Aum
V

Eficiencia (escala arbitraria)

Figura 4.13 Curvas hipotéticas típicas de la relación entre rendimiento y costo total.

Otra condición que se encuentra frecuentemente es la de considerar cuantos reactores

son adecuados para una determinada extracción o rendimiento dado. Por ejemplo, si
se desea extraer 95% de cobre de un mineral oxidado por lixiviación por agitación, se
requiere para un tiempo medio de retención del mineral de 6 hrs. Para una capacidad
de 35 TPH de mineral de cobre se requiere un reactor de MP de 7 m. de diámetro y 15
m. de alto que vale 120 millones de pesos completo. Si se emplean reactores 2 en
serie, estos tienen 2.2 m de diámetro c/u y 5 m de alto y valen 36 millones c/u, es decir
72 millones ambos, con lo cual el costo de capital decrece en 48 millones respecto de
emplear solo un reactor. Si este análisis se continua, para 5 reactores en serie, cada
uno vale 21 millones y el total 126 millones, es decir, ahora es más caro emplear 5
128
reactores de MP en serie que uno solo, de tal forma que existe un número óptimo que
es el más económico, como se observa en la Figura 4.14.

Como no se conoce a priori el volumen de cada reactor (puesto que no se conoce el

número óptimo) se debe hacer el cálculo desde n = 1 hasta n = n empleando el método
gráfico por ejemplo, o bien, si se conoce vo y ti = t / n en que n = 1, 2, etc.
Generalmente, 3 reactores en serie son suficientes y económicos, pero hay casos
como en el proceso de cianuración de minerales de oro en que se emplean 5 y hasta 7
reactores MP en serie debido a la lenta cinética de cianuración.

200
(Costo capital + opreación)

150
costo total $ x 10 6

100

50

0
0 1 2 3 4 5 6 7
Número de reactores en serie

Figura 4.14. Determinación del número más económico de reactores en serie

(valores hipotéticos).

4.4 DESVIACION DE LAS CONDICIONES DE FLUJO IDEALES

En la práctica, los reactores se comportan sólo en forma aproximada como flujo

pistón ideal o mezcla perfecta ideal y a veces se alejan por completo de ambos
patrones de flujo debido a la velocidad no uniforme del fluido y/o sólido; a
fluctuaciones de velocidades debido a difusión convectiva o molecular; a corto circuito
del fluido y/o sólido; canalizaciones y recirculación interna o externa del reactor o
volúmenes estáticos (muertos) dentro de éste. Estas condiciones no ideales requieren
por lo tanto caracterizar en forma cuantitativa el comportamiento de las entidades en
el interior del reactor para predecir el comportamiento de este y su rendimiento.

4.4.1 FUNCIONES DE DISTRIBUCIÓN DE EDADES: MODELOS POBLACIONES

Si se considera que la alimentación a un reactor y su masa interna está formada de

una población de entidades discretas tales como moléculas, elementos de fluido,
partículas, etc., una forma de caracterizarla es mediante las denominadas funciones
129
de distribución de edades o funciones de distribución de tiempos de residencia (DTR)
de las entidades en el interior del reactor o en su descarga.

Si un sistema se considera entonces constituido de entidades discretas, se puede

establecer relaciones que permitan conocer la distribución del tiempo que las
entidades permanecen en el reactor y las interrelaciones entre las distintas entidades.

Dos funciones caracterizan el historial de las entidades en un reactor: estas son las
funciones de distribución de los tiempos de residencia o de edades (DTR) en el
interior del reactor y en su salida (descarga).

Para conocer la DTR en un reactor se emplea la técnica de trazadores en la cual se

introduce en la entrada del reactor (generalmente) una perturbación o señal
(trazador) de una determinada forma matemática, haciendo luego un análisis del
trazador en trasciente en la salida del reactor (o en algún punto del interior de éste)

Señal del Detección

trazador del trazador

ALIMENT. REACTOR PRODUCTOS

Señal de salida
Señal de entrada

Figura 4.15 Forma de emplear un trazador en un reactor cerrado (sin mezcla de los
reactivos antes o después de entrar o salir del reactor).

Este trazador puede ser un producto químico, un colorante, un elemento radioactivo

etc. el que no debe alterar ni el flujo ni la composición de los reactantes en el reactor.

Se designa por I(t) o simplemente por I la función frecuencia de distribución del

tiempo de residencia en el interior del reactor, con unidades (fracción de
edades/unidad de tiempo). La fracción de entidades de edades entre t y (t + ∆ t) en el
reactor será I(t) ∆ t .

Como I(t) es una función continua, se puede normalizar haciendo la sumatoria de

todas las fracciones de edades en el interior del reactor desde las que salen del reactor
de inmediato con t=0 hasta las que permanecen un tiempo = ∞, de tal forma que se
tiene:

∞
∫o
l (t) dt = 1 (4.42)

Donde t = 0 se refiere a un tiempo inicial, arbitrariamente elegido que se puede ser el

tiempo de entrada al reactor. La ecuación (4.48) implica que todos los elementos de

130
fluido del reactor tienen una edad entre 0 e ∞. De esta forma, las fracciones de
fluido en el reactor de edades inferiores a una cierta edad t1 y mayor que t1 serán
respectivamente:

t1 ∞ t1

∫ l (t) dt
o
y ∫
t1
l (t) dt = 1 - ∫ l (t) dt
o
(4.43)

Esto se puede observar en la Figura 4.16 para una función arbitraria cualquiera de I(t).

La función de distribución de tiempo de residencia o función de distribución de

frecuencia de edades en la descarga (salida) de un reactor se designa por E(t) ó E, con
unidades fracción de edades/unidad de tiempo. En igual forma que para la
distribución interna, para el intervalo t y (t + ∆ t) la función vale E (t) ∆ t y puede ser
normalizada:

∫
o
E (t) dt = 1 (4.44)

Al igual que para I, las fracciones de fluido de edades menores a t1 y mayor que t1 serán
respectivamente:

t1 ∞ t1

∫ l (t) dt
o
y ∫
t1
E (t) dt = 1 - ∫ E (t) dt
o
(4.45)

Esto se observa en la Figura 4.16 para una distribución arbitraria cualquiera de la

función E(t) en la descarga de un reactor.

Las funciones E(t) e (I(t) no son iguales, de hecho, sólo lo son cuando el reactor opera
como un reactor perfectamente mezclado, puesto que en tal caso la composición
interna es igual a la composición de la descarga del reactor.

La detección y análisis de las funciones E(t) I(t) puede hacerse mediante análisis
químico, radioactivo, DRX u otra forma desde muestras tomadas en el interior o la
descarga del reactor. Como esto último es más sencillo, se prefiere por lo general
analizar la función E(t) y no I(t).

131
 I (t) E (t)
t1 t1
t1 t1
I(t)dt E(t)dt
1- I(t)dt = 1- E(t)dt =
0 0
0 0

0
0

0
0
I(t)dt = E(t)dt
t1 t1

0 t1 0 t1
t t

Figura 4.16 Función de distribución de edades de residencia en el interior de un

reactor I(t) y en la salida de éste E(t). (Forma típicas).

Por otra parte, para el reactor en su totalidad se tiene que el promedio de la

distribución de tiempo de residencia en el interior del reactor es igual a:

∞
VR WR
t = ∫ t l ( t ) dt = = (4.50)
0
vº Fº

Donde t se definió como el tiempo medio para reactores MP o bien como el espacio
– tiempo en el caso de reactores sin mezcla o flujo pistón (FP) donde t = t .

El concepto espacio – tiempo representa el tiempo que toma en pasar un volumen de

alimentación en el interior del reactor. Si la densidad del sistema permanece
constante (o aproximadamente constante) para un reactor FP será t = L / uo .

4.4.2 DETERMINACIÓN DE LAS FUNCIONES DE DISTRIBUCIÓN DE TIEMPOS DE

RESIDENCIA

Para determinar las funciones de distribución se emplea la técnica de marcación de

entidades mediante el uso de trazadores, los que como se mencionó pueden ser las
mismas entidades que forman la alimentación al reactor o bien otras de propiedades
físicas iguales o similares a las de fluido, partículas o elementos que componen la
alimentación al reactor.

Por ejemplo, si se desea conocer la DTR del eje en un reactor (horno) de fundición de
cobre, se puede emplear un metal soluble en el eje que tenga peso específico similar al
eje y sea soluble en él y se conozca la concentración inicial en él (cero o finita), un
compuesto podría ser en este caso NiS. En algunos casos, el trazador puede ser el
mismo elemento o compuesto que se alimenta al reactor pero en la forma de un
isótopo de éste de corta vida media, el cual se puede detectar en cantidades muy
bajas, o bien un colorante u otro.

132
El trazador se introduce al reactor como una señal o perturbación en la alimentación.
La forma de la señal o perturbación empleada puede ser de cualquier tipo, por
ejemplo una función impulso (instantáneo) o una función escalón o sinusoidal, etc.,
de ellas, las más simples son las impulso (o función Delta-Dirac) y la escalón continua
(step function).
Las respuestas a estas perturbaciones iniciales inyectadas en la entrada (u otro punto)
del reactor simultáneamente con la alimentación, se analiza la salida de éste (u otro
punto). Si se ha introducido una señal de concentración inicial CoT del trazador la
función respuesta adimensional en la salida, denominada F(θ θ) para la función de
entrada escalón y C(θ θ) para la señal impulso o Delta-Dirac, representa la respuesta
del sistema en términos de la mezcla que tiene lugar en el reactor.

A) Señal de entrada contínua o señal escalón (step function)

Al emplear este tipo de señal (perturbación) del trazador, para un tiempo t = 0 se

comienza a alimentar el trazador junto con la alimentación a la entrada del reactor, la
que se mantiene constante durante la operación normal de este, analizando
simultáneamente la concentración del trazador en la salida del reactor. Se define la
función de distribución de edades en la descarga del reactor E(t) en la forma
θ) para una señal de entrada escalón como:
adimensional F(θ

CT
F(θ) = (4.51)
CoT

Donde el tiempo adimensional θ queda definido como

t
θ= (4.52)
t

CoT es la concentración inicial del trazador (constante) en la entrada del reactor y CT la

concentración de éste en la salida, para un tiempo cualquiera t.

La forma de la función de descarga en la salida del reactor, como respuesta a una

función escalón de entrada, puede tener diferentes formas límites dependiendo si el
reactor se comporta como flujo pistón, mezcla perfecta, flujo mixto intermedio, tiene
cortocircuitos; regiones estancadas, etc.

Si el reactor se comporta como flujo tipo pistón, el trazador no se detecta hasta que
t = L/uo y este se detectará en la salida con una concentración constante e igual a CoT ,
lo que significa que no hay ningún tipo de mezcla interna. Casos que se aproximan a
éste lo representan cañerías de gran longitud (en flujo laminar); la zona del etalaje de
un alto horno u hornos de cuba (Midrex, Wiberg) o columnas de relleno de lixiviación.
En la figura 4.17 se esquematiza este tipo de función y su respuesta.

133
 E(t) 1
F( )

0 0 1
t

Figura 4.17 Respuesta en la descarga de un reactor flujo pistón para una señal
continua escalón.

Si el reactor se comporta como un sistema de mezcla perfecta, cada entidad

(elemento de fluido, partícula, etc.) que entra al reactor se mezcla en forma
instantánea y completamente con el resto del volumen VR del reactor y por lo tanto
con la masa WR en el interior de este.

Lo mismo ocurre por lo tanto con el trazador, de manera que idealmente, si se

introduce el trazador desde un tiempo inicial t = 0, para t = 0 también se detectará
éste en la salida del reactor (aun cuando muy diluido) puesto que se ha mezclado con
el total de la masa en el interior del reactor en forma instantánea. En estas
condiciones, el balance de masa para el trazador se puede expresar en la forma:

(Tasa de variación de la concentración del trazador en el reactor) = (Concentración del

trazador en la entrada) – (Concentración de trazador en la salida).

Es decir, para un reactor de volumen VR y tasas volumétricas de entrada y salida

iguales a vo y v, respectivamente, el balance de masa se expresa en la forma siguiente:

dC  vº   v 
= CoT   − C T  
dt  VR   VR 

Si vo = v = constante (ξ = 0), se tiene

o o
dC v (C T - C T )
=
dt VR

dC T vo 1
o
= dt = dt (4.53)
C T - C T VR t

134
Puesto que E(t) = CT/ CoT para θ = t/ t para la señal escalón la forma adimensional será
igualmente F(θ) = CT/ CoT . La solución analítica de la ecuación (4.53) (que tiene la forma
de una expresión cinética de primer orden) para las condiciones de borde t = 0, CT = 0
y t = t, CT = CT (ó θ = 0, F(θ) = 0 y θ = 1, F(θ) = 1) es la siguiente:

E(t ) = 1 − exp − (t / t )
(4.54)
F(θ) = 1 − exp( −θ)

Esto se representa en la Figura 18 donde se puede observar que aun cuando el

trazador se agrega desde t = 0 a una concentración constante CoT por un tiempo
indefinido, la concentración aumenta solo lentamente puesto que cada entidad de
trazador que entrar al reactor se mezcla por completo e instantáneamente con el resto
de la carga. Esto explica el porque del aumento tan lento en el tiempo.

F( )
E(t)
1

0 0 1
t
Figura 4.18 Respuesta para una señal continua de entrada (señal escalón) en un
reactor perfectamente mezclado.

En ciertas circunstancias, es conveniente emplear la función recíproca de F(θ), es decir,

la extinción de la señal y no el aumento. En este caso, la función extinción R(θ) queda
definida por R(θ) = 1- F(θ).

Si el reactor tiene zonas muertas o cortocircuitos, o presenta una combinación de flujo

pistón y mezcla perfecta, la señal de salida será distinta a la esperada a los
comportamientos extremos flujo pistón y mezcla perfecta.

B) Señal de entrada impulso

En muchos casos no es conveniente emplear una señal de entrada continua como la

vista anteriormente porque ésta puede contaminar el producto o ser de costo alto. En
tal caso, es más adecuado agregar el trazador como una señal instantánea o impulso
llamada también señal Delta - Dirac, mediante la adición instantánea de una cierta
cantidad de trazador, generalmente pequeña, en la entrada del reactor o en otro
punto de éste.

135
Así como se definió la función adimensional de distribución de tiempos de residencia
F(θ) para una señal escalón, se define igualmente la función de distribución
θ) en la salida del reactor para una señal delta en la forma:
adimensional C(θ

CT CT
C(θ) = º
= (4.56)
C T WT / VR

Donde CT es la concentración del trazador en la salida del reactor, CoT la concentración

inicial de éste WT la cantidad de trazador introducido y VR el volumen del reactor. En
igual forma, se tiene entonces que:

CT
C(θ)dθ = dθ = 1 (4.57)
C Tº

y puesto que θ = t / t , dθ = dt / t

C(θ)dθ = C( t )dt
1
t

CoT se puede calcular entonces conociendo WT y VR, o bien a partir del área bajo la
curva respuesta obtenida.

Para los casos extremos, se tiene que para flujo pistón la totalidad del trazador se
detecta en la salida del reactor para t / t = 0 = 1 y durante un tiempo infinitamente
corto puesto que la mezcla es nula. (Figura 4.18).

Si el reactor opera en mezcla perfecta, (la mezcla es instantánea) la concentración

mayor se alcanzará para t = 0 (inicial) (Fig. 4.20 a). Para un reactor de mezcla total, se
puede hacer un balance de materia para el trazador y puesto que no hay acumulación
en el reactor de mezcla perfecta, se puede escribir que:

dC T v o
- = CT (4.58)
dt VR

O bien,

dC T v o
- = dt (4.59)
CT VR

136
 C( )
E(t)

0 0
t 1

Figura 4.19. Respuesta de un reactor de flujo pistón a una señal impulso.

La solución de esta ecuación para las CB de t = 0, CT = 0 (θ = 0, CoT = 0) y t = t, CT = CT

(θ = 0, CT = CT) es entonces:

CT
E( t ) = = exp− ( t vº / VR ) = exp - (t Fº /WR ) = exp - (t / t) (4.60)a
CTº

y en forma adimensional:

C(θ) = exp− θ (4.60)b

Como

E(θ ) = C(θ ) = t E( t )
1 (4.60)c
E( t ) = exp− ( t / t )
t

Esto se puede observar en la Figura 20 para un reactor perfectamente mezclado en

que la concentración máxima se alcanza para θ = 0 (t = 0) es decir, al inicio de la
inyección del trazador.

C( )
E(t)
1

0 0 1
t
Figura 4.20. Respuesta de un reactor perfectamente mezclado a una señal impulso
137
Si se tienen n reactores perfectamente agitados en serie de igual tamaño, la función de
distribución de tiempos de residencia en la descarga para una perturbación de entrada
delta vale:

n(nθ)n−1 −nθ
C(θ) = E(θ) = e (4.61)
(n − 1)!

O bien:

nn ( t / t i )n−1 −nt / ti
C( t ) = E( t ) = e (4.62)
t i (n − 1)!

Para cada reactor se tiene que t i = VRi / v o y para la totalidad tT = n t i considerando

ξ ∼ 0.

4.4.3 PROPIEDADES ESTADISTICAS DE LA DTR

Se pueden definir propiedades estadísticas de las distribuciones (poblaciones) de

entidades en un reactor o en su descarga las cuales son útiles para su análisis de
comportamiento. De esta forma, se tiene que la distribución de una población de
entidades se puede caracterizar mediante el tiempo medio de residencia o el
denominado centroide de distribución, definido como:

t=
∫
o
tCdt
(4.63)
∞
∫o
Cdt

En que C es la concentración de cualquier entidad, por ejemplo, el trazador en el

tiempo t. Para valores discretos de t se puede también colocar para t :

t∼
∑ t C ∆t i i i
(4.64)
∑ C ∆t i i

La dispersión de la distribución por otra parte, se puede caracterizar mediante la

función varianza estadística o cuadrado de la desviación estándar, definida en la
forma:

∫ (t − t ) C
2
dt ∞
∫ t C dt − t
2

σ =
2 0
∞
= 0
∞
2
(4.65)
∫C dt ∫ C dt
0
0

La varianza representa el cuadrado de la amplitud de la distribución en (tiempo)2.

Para valores discretos de t la varianza se puede aproximar a:
138
 ∑ (t − t ) C ∑ t C ∆t
2
∆t i 2

σ ≈ = − t2
2 i i i i
(4.66)
∑ C ∆t i i ∑ C ∆t i i

Cuando se aplica a funciones normalizadas de DTR, las relaciones anteriores toman las
formas siguientes:

∞ ∑ t E(t) = t E(t ) ∆t
t = ∫ tE( t ) dt ≈ ∑i i
(4.67)
∑ E( t)
0 i i
i

∞
σ2 = ∫0
t 2E( t )dt − t 2 (4.68)

σ 2
≈
∑ t E( t) 2
i i
− t 2 = ∑ t i2 E( t )i ∆t i − t 2 (4.69)
∑ E(t ) i

y por lo tanto, E(θ)= tCi / ∑ Ci ∆t i

En igual forma que para las funciones adimensionales, se puede definir la varianza
adimensional σ 2 (θ ) en la forma σ 2 (θ ) = σ 2 / t 2 .

Finalmente, puesto que tC = tE = t también se tiene que θC = θE = θ .

La varianza es una propiedad aditiva, si hay n reactores en serie cada uno

caracterizado por su varianza σi se tiene que:

σ 2 total = σ12 + σ 22 + ...... + σn2 (4.70)

4.4.4 CÁLCULO DE LA CONVERSIÓN EN UN REACTOR MEDIANTE EL EMPLEO DE

TRAZADORES.

Si la respuesta a una señal de un trazador en un reactor es lineal con la concentración

del reactivo que interesa, se puede calcular la conversión directamente de la función
de tiempos de residencia determinada con el trazador.

En este caso, se tiene que la concentración media del reactivo que se desea conocer en
la salida del reactor es la sumatoria del producto de la concentración del reactivo que
existe en un elemento de edad entre t y (t + dt) por la fracción de entidades en la
salida que tienen una edad también entre t y (t+dt), es decir, E (t)dt, de tal forma que
se tiene para un reactivo A que

139
 ∞
CA = ∫ C eA E(t)dt (4.71)
0

Donde CeA es la concentración de A en un elemento de fluido con edad entre t y (t+∆t).

Si la reacción que ocurre en el reactor es irreversible y de primer orden respecto de A

y orden cero respecto de los otros reactivos, la expresión cinética para cada elemento
que contiene el reactivo A queda dado por:

dC A
− r&A = = krCA
dt

Cuya solución para las CB de t = 0, C A = CoA y t = t, C A = CeA es:

CeA = CoA exp – (kr t)

CoA es la concentración de A en la entrada del reactor y t su tiempo de retención. Por lo

tanto, la concentración promedio en la salida del reactor será entonces:

∞
CA = C ∫ exp− (k t ) E( t )dt
º
A r (4.72)
0

La forma de la función E(t) dependerá de la mezcla que ocurra en el reactor. Por

ejemplo, si el reactor opera en régimen de mezcla perfecta, la expresión anterior
tomará la forma

∞
CA = C ºA ∫ [exp− (k r t )]  exp− (t / t ) dt
1 
0 t 
(4.73)

Si el reactor opera como flujo pistón, la expresión (4.70) se reduce a

∞
CA = C ºA ∫ [exp− (k r t )] dt
0

Si el trazador tiene una respuesta lineal con la concentración de la especie A,

conociendo CT se puede entonces conocer CA .

140
4.4.5 INTERPRETACIÓN DE SEÑALES TÍPICAS DE TRAZADORES

En los reactores industriales las respuestas a las señales de entrada de trazadores no

necesariamente tienen las expresiones idales extremas de flujo pistón o mezcla
perfecta, aún cuando en muchos casos se aproximan bastante a una de ellas. En otros
casos, sin embargo, las respuestas pueden diferir considerablemente de la condición
ideal.

En la Figura 4.21 se indican respuestas típicas a señales de entrada escalón y delta en

reactores reales.

F( )
(a) C( ) (b)
FP
FP
1 1

P
M

MP

0 1 0 1
Figura 4.21 Respuestas típicas en reactores industriales
(a) señal de entrada escalón en reactores MP y FP ……ideal ___real
(b) señal de entrada impulso en reactores MP y FP …..ideal ___real

En cualquier caso, sin embargo, el área bajo o sobre la curva normalizada (señal
impulso o escalón, respectivamente) sigue siendo unitaria. También se puede deducir
de las curvas de las funciones de distribución E(t) o C(θ) que el grado de mezcla queda
relacionado con la altura y extensión de la dispersión de la curva obtenida. A mayor
grado de dispersión, mayor extensión de la curva y menor altura de la señal, como se
verá más adelante en los modelos que representan el grado de mezcla en un reactor,
de acuerdo con el grado de dispersión caracterizado por un parámetro llamado
coeficiente de dispersión.

141
 F( ) C( )
1
1

0 1 0 1

Figura 4.22 Reactor diseñado para mezcla perfecta cortocircuito: parte de los
reactivos salen del reactor sin reaccionar). El área oscura es proporcional al volumen
de cortocircuito (señal escalón y señal delta).

La forma de la señal de salida del reactor puede también dar información en cuanto a
canalización, cortocircuito y zonas estancadas que puede tener un determinado
reactor. En las figuras 4.21 a 4.24 se muestran respuestas típicas a señales de entrada
escalón y delta cuando existen estos problemas en reactores diseñados para
comportarse según un flujo de mezcla perfecta o flujo pistón.

F( ) C( )
1
1

0 1 0 1

Figura 4.23 Reactores que muestran un grado de cortocircuito: parte de los reactivos
salen del reactor sin reaccionar. (Reactor diseñado para flujo pistón). La zona oscura
es proporcional del volumen de cortocircuito.

142
 C C C

Cortocircuito Recirculación
Flujo esperado
temprano interna

Tiempo
medio

t t t
(a) (b) (c)

C C

Zona estancada Zona estancada

(respuesta tardía) (respuesta tardía)

t t
(d) (e)

Figura 4.24. Reactores que muestran un grado de zonas muertas (estancadas). Gran
parte de los reactivos salen antes o después de t del reactor.

4.5 MODELOS DE MEZCLA EN REACTORES

Desde el punto de vista del número de fases presentes, los modelos de mezcla en un
reactor se pueden clasificar en homogéneos y heterogéneos según si están
compuestos de una o más fases. Existen varios modelos de mezcla homogéneos, sin
embargo dos son los más empleados en sistemas metalúrgicos:

- Modelo de dispersión, en el que se supone que la mezcla entre las fases o

entidades reaccionantes ocurre por medio de dos mecanismos: convectivo y
difusivo. Esto representa un ensamble de un modelo de difusión con un modelo
de flujo pistón.

- Modelo de compartimentos o celdas en series, en que se supone el reactor

formado por una serie equivalente de compartimentos o celdas, cada una de
ellas perfectamente mezcladas.

Para los modelos heterogéneos (por ejemplo sólido-gas) y reactores de mezcla total no
hay un modelo que permita predecir con exactitud las características de flujo y mezcla
del sólido y del gas o los tiempos de residencia debido a la complejidad del problema.
Algunos de los modelos propuestos consideran un flujo de dos fases y otros un flujo de
tres fases.

Debido a estas limitaciones, para sistemas heterogéneos se consideran por lo general

ambas fases separadamente, cada una definida por un flujo particular y relacionadas
143
entre sí por coeficientes de transferencia de masa y calor. Para sistemas complejos,
como en lechos fluidizados sólido-gas, se consideran flujos de dos y tres fases, con
coeficientes individuales entre ellos.

Si las partículas o entidades son muy pequeñas y estas están suspendidas

uniformemente en un fluido (por ej. un concentrado fino que se lixivia con ácido en un
reactor agitado), se puede considerar el fluido o pulpa como un medio continuo y
emplear un modelo homogéneo.

Existen varios modelos de mezcla homogéneos, sin embargo dos son los más
empleados en sistemas metalúrgicos:

- Modelo de dispersión, en el que se supone que la mezcla entre las fases o

entidades reaccionantes ocurre por medio de dos mecanismos: convectivo y
difusivo. Esto representa un ensamble de un modelo de difusión con un modelo
de flujo pistón.

- Modelo de compartimentos o celdas en series, en que se supone el reactor

formado por una serie equivalente de compartimentos o celdas, cada una de
ellas perfectamente mezcladas.

4.5.1 MODELO DE DISPERSIÓN

Este modelo puede aplicarse para reactores metalúrgicos tales como reactores
continuos de fusión - conversión de cobre, flujo de la solución de lixiviación en
columnas de relleno o en bateas flujo del sólid y del gas en la cuba de alto horno;
columnas de absorción de gases, etc. y en general en reactores en que se puede
caracterizar un patrón de flujo direccional, sin embargo, no se puede aplicar a
reactores que tienen un comportamiento cercano a mezcla total, tal como lechos
fluidizados, lixiviación continua en estanques agitados, etc.

El modelo de dispersión considera que la mezcla ocurre a lo largo del reactor el cual
tiene un flujo tipo pistón y en el que la difusión produce la mezcla axial y
radialmente, es decir, es un modelo en el que se impone un modelo difusivo en uno
convectivo de flujo pistón.

Si se considera un fluido compuesto de entidades discretas tal como moléculas,

elementos de fluido o partículas macro o microscópicas que forma una pulpa
homogénea que se mueve unidireccionalmente con velocidad media u y que forman
un seudo – continuo entre el sólido suspendido y el fluido, se puede escribir para un
componente cualquiera A presente en el fluido o la pulpa (ver ecuación 2.26, por
ejemplo):

∂C A
= ∇ ⋅ (DC A ) − u ⋅ ∇ C A ± r&A
∂t tér min o
difusivo
tér min o
convectivo

144
Si se considera un reactor cilíndrico (simetría radial) en el que ocurre dispersión radial
y axial de los reactivos, se tiene que:

∂C A ∂ 2C A  ∂ 2C A 1 ∂C A  ∂C A &
= Da + Dr  +  − u ± rA (4.77)
∂t ∂x 2  ∂r
2
r ∂r  ∂x

Da y Dr son los coeficientes de dispersión axial y radial, respectivamente, de iguales

dimensiones y significado que el coeficiente de difusión, los que se pueden considerar
constantes, independientes de la concentración y posición en el reactor.

Cuando la dispersión radial es despreciable (en particular, cuando la relación

longitud/diámetro de un reactor es grande, aprox. > 5), o no hay gradiente de
concentración radial en el reactor por ser esta cero 0 o ∞ , la ecuación (4.77) se
reduce a:

∂C A ∂ 2C A ∂C A &
= Da −u ± rA (4.78)
∂t ∂x 2
∂x

Este submodelo se denomina modelo de flujo pistón con dispersión axial o

simplemente modelo de dispersión y se ha aplicado con buenos resultados en
metalurgia, en especial cuando la razón longitud/diámetro del reactor mayor a 5 y el
flujo turbulento.

Si para caracterizar el patrón de flujo en un reactor que se tiene un comportamiento

según este modelo, se puede emplear un trazador de concentración inicial de entrada
C oT . Las condiciones iniciales y de borde de la expresión (4.78) pueden ser entonces:

t = 0 t > 0
x = 0 x = 0
CT = 0 CT = CT
o

Otra CB dependerá de la señal empleada:

Señal impulso Señal escalón

t → ∞ t → ∞
x = L x = L
CT → 0 CT = CT
∂C T ∂C T
= 0 = 0
∂x ∂x
x = 0 x = L

La expresión (4.78) se puede colocar en forma adimensional definiendo una distancia

adimensional Z y un tiempo adimensional θ en la forma siguiente:

Z = ( u t + x) L

145
θ = t / t = tu / L

De esta manera, se puede escribir para la ecuación normalizada (4.78):

∂C A Da ∂ 2C A ∂C A
= − (4.79)
∂θ u L ∂z 2 ∂z

La solución de esta ecuación diferencial parcial no lineal toma distintas formas de

acuerdo con las condiciones de borde empleadas, es decir, si se emplea el una señal
F(θ) o C (θ) (señal escalón o delta).

Para una señal escalón, la solución de la ecuación (4.79) es:

  1− θ  
F(θ) = C T / CTº = 21 1 − fer  0.5 0.5  
  2θ (Da / u L )  

Para una señal impulso será a su vez:

C(θ) =
CT
=
C
=
1
C T WR / VR 2[πθ(Da / u L )]
º
2
[ ]
exp− (1 − θ) / 4θ(D a / u L )

En ambos casos se observa que, la función de distribución de tiempos de residencia

depende sólo del parámetro (Da / uL ) que representa el inverso del número de
Péclet para transferencia de masa (Pe-1), también llamado Número de Dispersión del
Reactor (RDN) o Número de Levenspiel (Le):

Da
= RDN = Pe- 1 = Le
uL

El RDN representa la razón entre la cantidad de masa transferida por dispersión

(difusión) a la transferida mediante flujo convectivo (flujo pistón). Conociendo el
RDN de un reactor es posible entonces estimar el grado de mezcla que ocurre en éste
con respecto a reactores ideales de FP o MP, así como se puede calcular el valor de
Da.

146
 1.0

0.8

F( ) 0.6 MP
RDM =

0
0
0.4
FP

-2
0 .2 RDN = 0

10
0.2

-3
0

x
2 .5 2 x 1
0 0.5 1.0 1.5 2.0

θ) calculada del modelo de dispersión.

Figura 4.25. Valores del RDN para la función F(θ

En las figuras 4.25 y 4.26 se indican valores del Pe-1 (o RDN) para señales escalón e
impulso para reactores que operan en flujo pistón, mezcla total y valores intermedios.
Se observa en la Figura 4.25 que el RDN decrece rápidamente aún para valores
relativamente grandes del grado de mezcla en el reactor. Por ej., 2 x 10-3 indica una
pequeña dispersión, en tanto que 2 x 10-2 indica una alta dispersión.

El diagrama F(θ) es más sensible que C(θ) a pequeñas variaciones en la forma de la

curva, por lo que es preferible emplear por lo general este último.

Para la señal delta se observa en la Figura 4.26 que para valores altos del RDN la curva
C(θ) tiende a condiciones de mezcla perfecta y para valores del RDN pequeños a flujo
pistón. Para desviaciones pequeñas de la condición de flujo pistón (Da / u L ) ≤ 10 −3 la
dispersión es aproximadamente simétrica con respecto al valor de t y la curva C(θ) se
puede representar mediante la distribución normal de Gauss:

C(θ) =
1
2[π(D a / u L )]
0.5
[
exp− (1 − θ) (D a / u L )
2
]

147
 2.0
0.01
- Baja dispersión:
RDN < 0.02
- Mediana dispersión:
RDN = Le = 0 RDN < 005 - 0.1
- Alta dispersión:
0.02 RDN > 0.2
1.5
- Cercano a mezcla total:
RDN > 1
0.05

E( )
RDN = Le = 0.1

0
0
10
1.0 5 2
0.2
1 0.5

0.5
Mixed
flow

0.0
0.0 0.5 1..0 1.0 2.0
-
=t/t

θ) y valores del número de RDN calculados del modelo de

Figura 4.26. Función C(θ
dispersión.

Las curvas de las Figura 4.25 y 4.26 fueron calculadas numéricamente para reactores
cerrados, es decir sin mezcla antes de la entrada o después de la descarga, que es lo
normal en la gran mayoría de reactores metalúrgicos y químicos.

Para reactores cerrados, la varianza adimensional queda relacionada con el RDN

mediante la expresión:

 Da   Da    u L 
2 2
σ
σ (θ) =   = 2Le − 2Le [1 − exp− (Le )] = 2
2 2
 − 2  1 − exp−   (4.80)
t  uL   uL   D
 a 

Para pequeñas desviaciones de flujo pistón, la expresión (4.80) se simplifica a:

D 
σ 2 (θ) ≅ 2Le = 2 a  (4.80)a
 uL 

Se puede observar que el número de Levenspiel (Le) o número de dispersión del

reactor (RDN) (inverso del Peclet) da una medida directa del grado de mezcla del
reactor por simple comparación con valores intermedios de éste entre sus extremos
de flujo pistón y mezcla perfecta.

En la práctica, con los datos experimentales obtenidos de CT versus t se compara la

curva la cual, para valores de Da dados en la expresión (4.83), se aproxima más a los
148
valores experimentales determinados. El valor de la curva experimental que más se
aproxima a la analítica representa el valor de Da. Se debe tener presente que en este
modelo (ecuación 4.83) se ha supuesto sólo difusión axial, es decir, el reactor tiene una
razón L/D >5 y el flujo es esencialmente unidireccional.

4.5.2 MODELO DE CELDAS DE MEZCLA PERFECTA EN SERIE

En este modelo, llamado también de celdas o estanques perfectamente mezclados en

serie, se considera que la mezcla se produce como si el reactor estuviera formado de
un número hipotético de celdas o reactores en serie de igual volumen y geometría
cada una de ellos perfectamente mezclados y que forman el volúmen total del
reactor (Figura 4.27)

n i
VR = VR
f=1

Co Cn
0
0

0
0

0
0

0
0
1 1 2 2 i i n n
C VR C VR C VR C VR

Figura 4.27. Modelo de reactor formado por un número hipotético de celdas en serie.

Si se considera un reactor de volumen total VR formado de una serie hipotética de

reactores en serie de volúmenes iguales VRi con una tasa volumétrica vo constante de
alimentación al reactor y considerando ξ = 0, el balance de materia en un reactor
cualquiera i se puede expresar en la forma:

dCi
VRi = v oCi- 1 - v oCi (4.81)
dt

La solución de esta ecuación diferencial se puede hacer para condiciones en que la

perturbación (trazador) es tipo delta o escalón. Para un impulso delta, las CB serán:

VR
Ci = 0 Ci = CTo
vo
t = 0 t > 0

149
en que CTo es la concentración promedio del trazador si se distribuyera
instantáneamente y en forma uniforme en todo el volumen VR del reactor, es decir,
una condición de mezcla perfecta total para todo el reactor.

La ecuación 4.84 se puede resolver mediante integral de convolución o bien mediante

transformada de Laplace, obteniendo así para el último compartimento n del reactor:

n
V  1  v º   vº 
C = Cº R
n
  n  t n−1 exp -  nt  (4.82)
vº  (n − 1)  VR  V
 R 

Puesto que

t = VR / v o y θ q= =t /t/t t

Se tiene que:

C n(nθ)n−1
E(θ)n = C(θ)n = = exp− (nθ) (4.83)
Cº (n − 1)!

En la Figura 4.28 se encuentran trazadas las curvas de respuesta para n = 1 hasta n =

∞. Se observa que a medida que crece n la respuesta se acerca a la de un reactor de
FP.

n=
0
0

2.0
E( ) = C( )

n=20
1.5

n=10
n=3
n=5
1.0 n=2
n=1

0.5

0 0.5 1.0 1.5

Figura 4.28. Respuesta a una señal delta para reactores perfectamente mezclados en
serie.
150
En igual forma, si se emplea una función escalón, la solución de la ecuación (4.81)
toma la forma ahora de:


F(θ)n = 1 − e nθ 1 + nθ + (nθ) + ....... +
1 2 1
(nθ)n−1  (4.84)
 2 (n − 1) 

En todos los casos, se tiene un tiempo medio de residencia en cada reactor igual a
ti = VRi / v o y un tiempo medio de residencia para el total igual a tT = VR / v o .

En la Figura 4.29 se han trazado las curvas que representan el % de entidades con edad
mayor que tT para distintos números de reactores perfectamente agitados en serie
para una perturbación (señal) escalón.

100
20 10 5
% de entidades en la descarga

3 2
1
80
_
de edad mayor que t

n=3
60
n=2
n=1
40
n=
0
0

20 n=5
n=10
n=20
0 0.5 1.0 1.5 2.0 2.5
t/t=

Figura 4.29. Distribución de las edades en la descarga que tienen una edad > t .

4.6 CONVERSION DE PARTICULAS SOLIDAS EN UN REACTOR.

4.6.1 CONVERSIÓN DE PARTÍCULAS SÓLIDAS DE IGUAL TAMAÑO

Cuando en las reacciones sólido-fluido estas son dos fases separadas bien
diferenciadas, sus respectivas expresiones para la transformación en el tiempo deben
ser resueltas en forma separadas por tener diferentes patrones de flujo, como en la
cuba y etalaje del alto horno en que el sólido y el gas tienen patrones de flujo
diferentes.

Sin embargo, si el sólido es muy fino y se encuentra uniformemente disperso en el

fluido, se puede hacer abstracción del carácter particulado y considerar un patrón de
flujo del sólido igual al del fluido, con una cierta concentración del reactivo sólido por

151
unidad de volumen de mezcla y de esta forma calcular la conversión del sólido en la
forma vista. Por ejemplo, para la reacción:

aA ( s ) + bB ( g) ⇔ cC ( s ) + dD( g)

1 dnA dnB a
- = = kr CnA CmB
ˆS dt dt b

Si n = 0 y m = 1

dn A a
− = ŝ k r CB
dt b

En que ŝ es la superficie del sólido, L2 y kr la constante de velocidad, L/t.

Dependiendo del mecanismo que prevalezca a nivel de partícula individual, se tendrá

una expresión para el tiempo que se necesita para la conversión completa del sólido:
transporte en el film, difusión en la capa de sólido o reacción química, o combinación
de ellas. Si se desea conocer el tiempo de reacción necesario para alcanzar un cierto
grado de conversión del sólido A, se debe conocer la función de transformación
individual de las partículas.

a) Control por reacción química

Para una esfera densa que se reduce de tamaño, se tenía de acuerdo al modelo de
núcleo sin reaccionar:

= 1 − (1 − X A )
t 0.33

es decir, para partículas individuales y de igual tamaño se tiene que:

3
 t
1 − X A = 1 − 
 τ

Puesto que para un reactor perfectamente mezclado se tiene que la distribución de

tiempos de residencia es igual en el interior del reactor que en la descarga e igual a:

1
E(t)= exp- (t / t )
t

Para la totalidad de la población de partículas se tendrá que la fracción transformada

promedio en la descarga (igual a la interna del reactor) valdrá:

152
 ∫ (1 − t / τ) exp− (t / t )dt
τ 3 1
1− XA = (4.88)
0 t

Integrando, se tiene:

2 3
t t
[
X A = 3  − 6  1 − exp(− τ / t ) ] (4.89)
τ τ

Expandiendo esta última expresión como una serie para valores de t / t > 1 se tiene:

2 3
1 τ 1  τ 1 τ
1− XA =   −   +   (4.90)
4  t  20  t  120  t 

si t / t > 5 , se puede aproximar a:

1 − X A ≅ 0.25τ / t (4.91)

b) Control difusional a través de capa de producto sólido (reactor de mezcla

perfecta)

Para partículas esféricas con control difusional a través de la capa de productos

formados, se obtiene en igual forma que:

3 4
τ τ τ  τ
1 − X A = 0.2 − 0.045 + 0.009  − 0.001  (4.92)
t t t t

c) Control por transporte en el film de fluido que rodea la partícula (reactor de

mezcla perfecta).

Considerándo igualmente partículas esféricas, para control por transporte en el film se

tiene en igual forma:

∫ (1 − t / τ) E(t ) dt
τ
1− XA =
0

y entonces:

2 3
1 τ  1  τ   1  τ 
1− XA = −    +    (4.93)
2 t  3  t   4  t 

En el caso más general, cuando la alimentación tiene una distribución de tamaño, esta
se puede expresar mediante una relación como la de Gaudin-Shuhmann u otra. De
153
esta forma, se pueden integrar en forma numérica las expresiones resultantes, como
se verá más adelante.

4.6.2 CONVERSIÓN EN UN REACTOR SEGÚN EL MODELO DE DISPERSIÓN

Como se indicó anteriormente, éste modelo considera que la mezcla ocurre de

acuerdo a un flujo pistón en el que existe mezcla axial difusiva. Haciendo un balance
para el reactivo (fluido) en una sección longitudinal de un reactor que opere de
acuerdo a este modelo, se tiene que para una reacción de primer orden respecto al
fluido B y orden cero respecto del sólido:

dC A d2C A
uo − Da = krC A (4.94)
dy dy 2

Las CB de la expresión (4.94) son:

dC A dC A
i) y=0 = uo C A + Da
dy dy
dC A
ii) y =L =0
dy
La primera CB específica que en la entrada del reactor la velocidad de entrada del
reactivo A al reactor es igual a la velocidad combinada de convección y difusión a la
cual éste es transportado al interior del reactor. La segunda CB indica que para y = L
no hay mezcla más allá de L, es decir, el reactor es cerrado en el sentido que no hay
mezcla antes de la entrada a este o más allá de su descarga.

La expresión (4.94) se puede colocar en función de una distancia adimensional Z y

tiempo medio t en la forma:

z= y /L
t = L / uo = VR / v o

De esta forma se tiene la ecuación diferencial normalizada:

 Da  d2 C A dC A
  2
− − kr t CA = 0 (4.95)
 u L  dz dz

Se observa en esta última expresión que a lo largo del reactor la concentración del
reactivo A está controlada por el grupo cinético k r , t C A y por el grupo adimensional
D / uoL ya visto, que es el inverso del número de Pèclet o RDN o Nº de Levenspiel. La
solución general de la ecuación (4.95) toma entonces la forma:

154
   2

0.5
  2
 2

0. 5

u 
 u  4k   u  
 u  4k 
C A = C1 exp 0.5  o +  o  + .... r   + C 2 exp 0.5  o  +  o  + r  
 Da  D a  Da    Da   D a  Da  
       

C1 y C2 se obtienen de las CB obteniendo así la solución particular:

  u L 
4φ exp 21  o 
CA
= 1 − XB =   D a 
( ) { ( )}
(4.96)
o
CA (1 − φ) exp 2 Da − (1 − φ)2 exp−
2 φ uoL φ
2
uoL
Da

en que φ es un parámetro igual a:

D 
φ = 1 + 4k r t  a 
 uoL 

Esta expresión se ha resuelto gráficamente en la Figura 4.30 como la razón entre el

volumen del reactor operando según el modelo difusivo y el que tendría si el flujo
fuera flujo pistón, para reacciones de primer orden.

30
Vol. reactor operando rég. difusivo
Vol. reactor operando flujo pistón

20 kr t = 50

10 16 - =
Da/uL (MP)
0
0

20 4
5
4
1
3 10
2 5
0.25
0.25 x 10 -2 2
=

1
1
10 -2 - = 0 (FP) 10 -1 10
Da/uL 1-XB = C B / CºB

Figura 4.30. Razón de volúmenes de reactores para modelo difusivo y modelo flujo
pistón. (Reacciones de primer orden respecto del fluido, orden cero respecto de
otros reactivos).

Se puede observar de la Figura 4.29 que el reactor de flujo pistón tiene ventajas
evidentes en cuanto a menor tamaño (volumen). Además se observa que un aumento

155
de la mezcla difusiva produce un aumento del tamaño del reactor necesario para
obtener igual nivel de conversión.

Para pequeñas desviaciones del flujo pistón, la expresión (4.95) se aproxima a:

 2  D 
= exp− k r t + (kr t )  a 
CA
o
CA   uoL 
(4.97)

o bien

CA  k r2 σ 2 
≅ exp − k r t +  (4.98)
C oA  2 

4.6.3 DISEÑO DE REACTORES CON DISTRIBUCIÓN DE TAMAÑO DE PARTÍCULAS

Los modelos cinéticos vistos anteriormente son aplicables a partículas de tamaño

uniforme, lo cual ocurre en el caso de pellets o bien minerales o concentrados con
tamaño de partícula uniforme, sin embargo, no lo son para concentrados o minerales
finamente divididos los que tienen una distribución amplia de tamaños.

La distribución de tamaño de partículas se puede expresar según alguna de las formas

conocidas, por ej., de Rosin-Rammbler o Gaudin- Schuhmann. En ésta última se tiene
que:

(
ψ(di ) = di / d̂ )
m
(4.100)

Donde:

0.7 < m < 1

ψ (di ) es la fracción acumulativa en peso de menor tamaño que di y d̂ el tamaño

máximo de las partículas en la distribución. m es la pendiente de la recta en el gráfico
logarítmico de % material acumulativo menor que un cierto tamaño vs. su
correspondiente tamaño. Las rectas muestran casi siempre una desviación en el rango
superior de tamaño, que representan normalmente menos del 10-15% del total del
sólido. La fracción de sólido ∆ψ (di ) de tamaño entre di y di + ∆di se obtiene
diferenciando la expresión (4.100):

m(di)m − 1∆di
∆ψ (di ) = = (4.101)
d̂m

Donde en una distribución normal de tamaño se tiene que:

156
∑ ∆Ψ(d ) = 1 i (4.102)

Para el cálculo y para simplificar el problema, las partículas se pueden aproximar a una
geometría regular tal como esferas. Hay casos sin embargo, en que esto no se cumple
y se debe emplear otra forma geométrica más apropiada, tal como plana o cilíndrica.

a) Reactores Batch y Flujo Pistón con distribución de tamaño

Tanto para el modelo de núcleo recesivo con controles individuales por transferencia,
químico o difusional en la capa de producto, así como para el modelo de reacción
generalizado en todo el volumen con control por transferencia en el film o control
químico, para reactores batch y de flujo pistón la expresión general para la fracción
promedio convertida en la descarga se puede entonces escribir en forma de una
sumatoria, como:

X A = ∑ [X A (di )][∆ψ (di )] (4.103)

X A (di ) es la función de transformación individual de una partícula de diámetro di.

Esta expresión se puede resolver en forma numérica para diferentes condiciones de
control y ambos modelo (núcleo recesivo y reacción generalizada), para reactores
ideales de flujo pistón y de mezcla perfecta. Los resultados obtenidos se observan en
las Figuras 4.31 y 4.32. En estas figuras τ es el tiempo total de transformación para
X A = 1 . En estas expresiones además se han considerado partículas esféricas y
ausencia de control mixto o control por transferencia de calor.

1.00 1.00
t Tamaño uniforme t
REACTORES DE FLUJO
PISTON Y BATCH
0.80 0.80
(b)
0.60 0.60
Tamaño uniforme
m=
0.40 1.0 0.40
m=
0.7
(a)
0.20 0.20
REACTORES DE FLUJO
m = 1.0
PISTON Y BATCH
m = 0.7
0.00 0.00
1.00 .95 .90 .85 .80 .75 .70 1.00 .95 .90 .85 .80 .75 .70
XA XA
Fi
Figura 4.31 Reactores de flujo pistón y batch

a) Fracción transformada en la descarga válido para:

- modelo de reacción volumétrico con control por:
- Transferencia en el film
- Químico
- Modelo de núcleo recesivo por control por transferencia en el film,
(partículas de tamaño constante).
b) Fracción transformada en la descarga para:
- el modelo de núcleo recesivo con control por difusión en la capa de producto
sólido formado.

157
 1.00
t
REACTORES DE FLUJO 0.4
PISTON Y BATCH
0.80 Def. t
REACTORES DE FLUJO
(a) r2
PISTON Y BATCH
0.3
Tamaño uniforme
0.60
(b)
0.2
0.40
Tamaño uniforme
0.1
m = 1.0
0.20 m m=
m=0 =0 1.0
.7 .7
0.0
0.00 1.0 0.9 0.8 0.7 0.6 0.5
1.00 .95 .90 .85 .80 .75 .70 XA
XA

Figura 4.32 Reactores de flujo pistón y batch

a) Fracción transformada en la descarga válido para:

- el modelo de núcleo recesivo con:
- Control Químico y control por transferencia (partículas que < tamaño, Re>50).
b) Fracción transformada en la descarga válido para:
- modelo de reacción generalizado con control por difusión en el producto sólido.

c) Reactores de Mezcla Perfecta con Distribución de Tamaño de Partículas

En el caso de reactores que operan como mezcla perfecta, la función de distribución

de tiempos de residencia en la descarga (o en el interior del reactor) queda dada por la
expresión:

1
C(t) = E(t) = e- t / t
t

Como E(t) es una función continua, la sumatoria se puede expresar como una integral
y la fracción media del material convertido en el descarga del reactor para un tamaño
de partículas queda dada por la expresión:

τ (di )
X A (di ) = ∫ X A (di ) E( t ) dt (4.104)
0

Si la distribución de tamaño de partículas siguen la correlación de Gaudin-Schuhmann

se tenía que para el rango de tamaño

m(di )m−1 ∆di

∆ψ (di ) =
d̂m

De esta forma se tiene que la fracción promedio del reactante B en la descarga del
reactor será entonces igual a
158
  m(di )m−1 
X A = ∑ X A (di ) m  (4.105)
 d̂ 

o bien,

τ 
X A = ∑ ∆ψ (di ) ∫ X A (di ) E ( t )dt  (4.106)
0 

Para resolver esta expresión se puede considerar que cada fracción de tamaño tiene
una distribución de tiempo de residencia en el reactor igual, sin embargo, si ésta se
conoce para cada fracción de tamaño, pueden emplearse los valores individuales,
(caso en que ésta se aleja de la condición de mezcla perfecta).

Los valores numéricos calculados para t / t en función de la fracción de sólido

transformado en la descarga o el interior de reactor se encuentran graficados en las
Fig. 4.33 y 4.34 para valores del parámetro de distribución de tamaño m entre 0.7 y 1,
para los modelos de reacción superficial (núcleo no reaccionado) y volumétrico
(reacción generalizada en todo en volumen).

20.0
20.0
t REACTORES DE t
MEZCLA PERFECTA
16.0 REACTORES DE
(a) 16.0 MEZCLA PERFECTA

(b)
12.0 12.0

8.0 8.0

Tamaño uniforme
4.0 m 4.0 Tamaño uniforme
m =1 m
=0 .0 m =1
.7 = .0
0.
7
0.0 0.0
1.00 .95 .90 .85 .80 .75 .70 1.00 .95 .90 .85 .80 .75 .70
XA XA
Figura 4.33. Reactores de Mezcla Perfecta:

a) Fracción transformada en la descarga válido para:

- modelo de reacción generalizada con:
control por transferencia y 2) control químico
- modelo de núcleo recesivo con control por transferencia, (partículas de tamaño
constante).
b) Fracción transformada de sólido en la descarga válido para:
- modelo de núcleo recesivo con:
control químico y control por transferencia, partículas de tamaño decreciente.

159
 20.0
t 0.4
REACTORES DE
16.0 MEZCLA PERFECTA Def. t
REACTORES DE
r2
(a) 0.3
MEZCLA PERFECTA

12.0 (b)

0.2
8.0 Tamaño uniforme

0.1 m
4.0 m =1
Tamaño uniforme =0 .0
.7
m

m = 1.0
=

0.0
0.

1.0 0.9 0.8 0.7 0.6 0.5

7

0.0
1.00 .95 .90 .85 .80 .75 .70 XA
XA

Figura 4.34. Reactores de Mezcla Perfecta

a) Fracción transformada en la descarga válido para modelo de núcleo recesivo con

control difusional en la capa de producto formado.

b) Fracción transformada en la descarga válido para el modelo de reacción generalizada

con control difusional.

c) Reactores de mezcla perfecta en serie

Cuando la reacción se efectúa en más de un reactor, con las propiedades de flujo

establecidas para el modelo de compartimentos o reactores en serie se tenía que la
concentración en la descarga del reactor n quedaba expresada por:

n−1
nn nn−1  t   nt 
E(t )n =   exp−  
t i (n − 1)!  t i   ti 

La expresión (4.104) toma entonces la forma:

 τ( di ) 
X A (n) = ∑ ∆ψ (di) ∫ X A (di ) E( t )n dt  (4.107)
 0 

La solución numérica como t / t versus XB se encuentra trazada en las Figs. 4.35 y 4.36
para n=1,2 y 3. Se observa que para conversiones altas (>95%), no se obtiene en gran
aumento de la conversión para n>3.

160
 20.0
20.0
t t
REACTORES DE REACTORES DE
MEZCLA PERFECTA MEZCLA PERFECTA
16.0
EN SERIE (m=1) 16.0 EN SERIE (m=1)

12.0 (a) (b)

12.0

8.0 8.0

4.0 n= 4.0
2 1
2 n=
1
3 3
0.0 0.0
1.00 .95 .90 .85 .80 .75 .70 1.00 .95 .90 .85 .80 .75 .70
X A X A

Figura 4.35. Reactores de mezcla perfecta en serie:

a) Fracción transformada en la descarga para distinto número de reactores en serie,

válido para:
- Modelo de reacción generalizada con:
control por transferencia y control químico
b) Modelo de núcleo recesivo con control por transferencia (partículas de tamaño
constante).
c) Modelo de núcleo recesivo con:
control químico y control por transferencia, (partículas de tamaño decreciente, Rep
<50).

1.00
t
REACTORES DE 0.4
MEZCLA PERFECTA
Def. t
0.80
EN SERIE (m=1) r2 REACTORES DE
0.3 MEZCLA PERFECTA
EN SERIE (m=1)
0.60 (a)
0.2 (b)
0.40
2 n
3

=
0.1 1
0

0.20
0

n=
2 1
3
0.0
0.00 1.0 0.9 0.8 0.7 0.6 0.5
1.00 .95 .90 .85 .80 .75 .70 XA
XA

Figura 4.36. Reactores de Mezcla Perfecta en serie

a) Fracción transformada en la descarga parea distinto número de reactores en serie, válido para
el modelo de núcleo recesivo con control difusional en la capa de producto sólido formado.
b) Fracción transformada en la descarga para distinto número de reactores en serie, válido
para el modelo de reacción generalizada control difusional.

Todos estos gráficos pueden ser empleados indistintamente en sistema heterogéneos

sólido-líquido como sólido-gas, dependiendo sólo del modelo mecanismo empleado.
161
BIBLIOGRAFÍA CAPITULO IV

- O. Levenspiel, “Chemical Reactin Engineering”, 3rd Ed., J. Wiley and Sons,

(1999).
- R.I.L. Guthrie, “Engineering in Process Metallurgy”, Claredon Press, (1989).
- G. F. Froment and K.B. Bischoff, “Chemical Reactor Analysis and Design”, J.
Wiley and Sons, (1979).
- N.J. Themelis, “Transport and Chemical Rate Phenomena”, Gordon and Breach,
(1995).

162