CRISTALIZACIÓN INDUSTRIAL

TALLER N° 3.11

10 DE MAYO DE 2016
SONIA CORTÉS HIDALGO
Alumna Magíster en Ciencias de la Ingeniería, Mención Proceso de Minerales
Profesora: Dra. María Elisa Taboada
 En un cristalizador que se supone MSMPR se han obtenido los siguientes
datos de distribución de tamaño de cristal obtenidos al tamizar (tamices Tyler):
Datos:
Tiempos de residencia, τ=33 min
Densidad del Cristal=ρ=1690 kg/m3
Factor de forma, kv=0.56
Volumen de muestra=350 ml

Abertura tamiz (mm) Masa en tamiz en g

0.71 0.00
0.50 18.10
0.36 55.40
0.25 309.10
0.18 121.30
0.13 41.64
0.09 10.70
0.06 3.10
0.00 0.50
Tabla 1: datos del problema

a) Dibuje la distribución de tamaño en % Retenido y en % acumulado

versus tamaño. Calcule el tamaño medio, el tamaño moda y el
coeficiente de variación CV.

diámetro abertura tamiz masa en tamiz % %

malla (mm) (gr) retención acumulado
24-32 0.605 18.1 3.233 3.233
32-42 0.43 55.4 9.896 13.129
42-60 0.305 309.1 55.212 68.341
60-80 0.215 121.3 21.667 90.008
80-115 0.155 41.64 7.438 97.446
115-170 0.11 10.7 1.911 99.357
170-250 0.075 3.1 0.554 99.911
250 0.03 0.5 0.089 100.000
Tabla 2: datos de porcentaje de Retención y Acumulado

1
La tabla 2 muestra los valores de % retenido y acumulado, para lo cual, se calculan
los valores medio de abertura de tamiz y de masa de tamiz, con los datos del
problema (Tabla 1).

 Para encontrar el % de Retención se calcula por la masa de cada tamiz sobre

la suma de masas en los diferentes tamices, en este caso tiene un valor de
559,84 gramos.

 Por otra parte, para él % Acumulado se sumaron los % de retención para

cada apertura de tamiz, es decir, se sumaron los % de Retención desde la
apertura de tamiz inicial hasta la deseada.

Distribución de tamaño de cristales

60.000 100.000
90.000
50.000 80.000
70.000

% Acumulado
40.000
% Retenido

60.000
30.000 50.000
40.000
20.000 30.000
10.000 20.000
10.000
0.000 0.000
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
tamaño en mm

Gráfica 1: Distribución de Tamaño de Cristales (DTC), % Retención versus Tamaño, en el eje

primario. Distribución de Tamaño de Cristales (DTC), % Acumulación versus Tamaño, en el
eje secundario.

De gráfica 1, se obtienen los datos de tamaño dominante, tamaño medio, L 16%, L84%.
Lo cual lleva a tener los siguientes valores y cálculos, que también se puede
apreciar en la gráfica 2, determinadas con los puntos en color morado:

2
 Distribución de Tamaño de Cristales
60 120

50 100

% Acumulado
40 80
% Retenido

30 60

20 40

10 20

0 0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
L (mm)

Gráfica 2: Distribución de Tamaño de Cristales (DTC), % Retención y % Acumulado. Datos

necesarios para el problema

 Tamaño moda o dominante:

𝐿𝐷 = 0.305 𝑚𝑚
 Tamaño medio:
𝐿𝑀 = 0.344𝑚𝑚

 Coeficiente de variación

100 ∗ (𝐿16% − 𝐿84% )

𝐶𝑉 =
2 ∗ 𝐿50%

𝐿16% = 0.420 𝑚𝑚
𝐿50% = 0.344 𝑚𝑚
𝐿84% = 0.246 𝑚𝑚

100 ∗ (0.424 − 0.239)

𝐶𝑉 = = 25.3%
2 ∗ 0.347

3
 b) Determine B, G y el coeficiente de correlación (R2) de ln(n) v/s L

Se debe encontrar el número de partículas por litro de volumen para cada diámetro
de malla, con la ecuación que sigue:

𝑚 # 𝑑𝑒 𝑐𝑟𝑖𝑠𝑡𝑎𝑙𝑒𝑠
𝑛= 3
=
𝐾𝑣 ∗ 𝑉 ∗ 𝐿 ∗ ∆𝐿 ∗ 𝜌𝑐 𝐿𝑖𝑡𝑟𝑜𝑠 ∗ 𝑚𝑚

Todos estos cálculos y datos necesarios para formar la gráfica 3 de densidad de

población versus longitud se muestran en la Tabla 3, ésta gráfica se emplea para
encontrar la velocidad de crecimiento y de nucleación de cristales.

 La velocidad de crecimiento se calcula mediante la fórmula:

1
𝑝𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 = −
𝐺 ∗ 𝜏0

 La velocidad de nucleación se calcula mediante la fórmula:

𝐵0 = 𝐺 ∗ 𝑛0

masa en tamiz (gr) Kv ΔL L (mm) n ln(n)

18.1 0.56 0.21 0.605 1175033.86 13.9768075
55.4 0.56 0.14 0.43 15025661 16.52527
309.1 0.56 0.11 0.305 298995069 19.5159376
121.3 0.56 0.07 0.215 526386333 20.081546
41.64 0.56 0.05 0.155 675154002 20.3304514
10.7 0.56 0.04 0.11 606740458 20.2236117
3.1 0.56 0.03 0.075 739458615 20.4214289
0.5 0.56 0.06 0.03 931777489 20.6526046
Tabla 3: datos para LN(n)

4
 Densidad de población
25

20

15
y = -11.686x + 21.778
Ln(n)

R² = 0.9024
10

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
L, mm

Gráfica 3: Densidad de población frente a Longitud

Desde el gráfico 3, utilizando un ajuste lineal, que queda expresado como y = -

11.686x + 21.778. Donde la intersección con la longitud igual a cero tiene un valor
de 21.778 y una pendiente de -11.686. También se encontró el coeficiente de
correlación que es R² = 0.9024.

Por lo tanto, el cálculo de velocidad de crecimiento y nucleación quedará de la

siguiente forma:
1 1
𝐺= − =−
𝑝𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 ∗ 𝜏0 −11.686 ∗ 0.528

𝒎𝒎
𝑮 = 𝟎, 𝟏𝟓𝟓𝟔
𝒉

𝐵0 = 𝐺 ∗ 𝑛0 = 0.04518 ∗ 𝑒 21.778

𝒏°
𝑩𝟎 = 𝟒𝟒𝟔, 𝟕𝟐𝟏𝟐 ∗ 𝟏𝟎𝟔
𝑳∗𝒉

5
 c) Determine la densidad de suspensión MT y el número total de cristales
NT
Para calcular la densidad de suspensión y el número total de cristales, se hace de
acuerdo a las ecuaciones de momentos, que siguen a continuación:

𝑀𝑇 = 6 ∗ 𝑘𝑣 ∗ 𝜌 ∗ 𝑛0 ∗ (𝐺 ∗ 𝜏)4

𝑁𝑇 = 𝑛0 ∗ 𝐺 ∗ 𝜏

Luego queda:

𝑀𝑇 = 6 ∗ 0.56 ∗ 0.00169 ∗ 𝑒 21.778 ∗ (0.1556 ∗ 0.55)4

𝒈 𝒄𝒓𝒊𝒔𝒕𝒂𝒍𝒆𝒔
𝑴𝑻 = 𝟖𝟕𝟒. 𝟐𝟑𝟑𝟗
𝑳

𝑁𝑇 = 𝑒 21.778 ∗ 0.1556 ∗ 0.55

#𝒄𝒓𝒊𝒔𝒕
𝑵𝑻 = 𝟐𝟒. 𝟓𝟔𝟗𝟔 ∗ 𝟏𝟎𝟕
𝑳

6
 d) Luego con los datos cinéticos, encuentre la DTC y el tamaño medio y
el tamaño moda y compare con los experimentales

𝐿
𝑋=
𝐺∗𝜏

Tamaño acumulado:

𝑀 𝑋2 𝑋3
𝑀(𝑋) = = 1 − [1 + 𝑋 + + ] ∗ 𝑒 −𝑋
𝑀𝑇 2 6

Tamaño retenido:
𝑑𝑀(𝑋) 1 3
= 𝑋 ∗ 𝑒 −𝑋
𝑑𝐿 6

M(X) dM(X)/dL
L (mm) X
acumulado Retenido
0.605 7.07003 0.92181294 0.05007691
0.43 5.02498 0.73846312 0.13897308
0.305 3.56423 0.47716154 0.21370786
0.215 2.51249 0.24509707 0.21428967
0.155 1.81133 0.11053565 0.16187913
0.11 1.28546 0.04165827 0.09789404
0.075 0.87645 0.01232716 0.04670816
0.03 0.35058 0.00047628 0.00505772
Tabla 4: datos de % acumulado y Retenido por ecuaciones

7
 Distribución de tamaño de cristales
0.25 1
0.9
0.2 0.8
dM(X)/dL Retenido

0.7

M(x) Acumulado
0.15 0.6
0.5
0.1 0.4
0.3
0.05 0.2
0.1
0 0
0 1 2 3 4 5 6 7 8
X

Gráfica 4: Distribución de tamaño de cristales de acuerdo a Ecuaciones

El tamaño moda y medio se calcula de acuerdo a las siguientes ecuaciones:

𝐿𝐷 = 3 ∗ 𝐺 ∗ 𝜏

𝐿𝑀 = 3.67 ∗ 𝐺 ∗ 𝜏

 Tamaño moda o dominante:

𝐿𝐷 = 0.257 𝑚𝑚
 Tamaño medio:
𝐿𝑀 = 0.314 𝑚𝑚

Nota: los puntos morados en la gráfica 4 me permiten observar los valores de X

para poder calcular L16 y L84: X16=2.1 y X84=6.07.

8
  Coeficiente de variación
𝐿 =𝑋∗𝐺∗𝜏

𝐿16% = 2.1 ∗ 0.1556 ∗ 0.55 = 0.1797 𝑚𝑚

𝐿84% = 6.07 ∗ 0.1556 ∗ 0.55 = 0.519 𝑚𝑚

𝐶𝑉 = 54%

experimental Ecuaciones
LM 0.305 0.257 mm
LD 0.344 0.314 mm
Cv 25.3% 54%
Tabla 5: tabla comparativa entre datos experimentales y las ecuaciones

% Acumulado Cristales
120 100

100 80
% acumulado Ecs.
% acumulado Exp

80
60
60
40
40

20 20

0 0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
L (mm)

Acumulado cristal exp Acumulado cristal Ecs

Gráfica 5: comparación de la curva de %Acumulado, de acuerdo a datos experimentales y

ecuaciones

9
 % Retenido
60 25

50 20

% Retenido Ecs.
% Retenido exp 40
15
30
10
20

10 5

0 0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
L (mm)

Retenido Exp. retenido Ecs.

Gráfica 6: comparación de la curva de %Retenido, de acuerdo a datos experimentales y

ecuaciones

e) En base a la forma de la curva ln(n) v/s L ¿considera posible un

alejamiento de la idealidad del MSMPR? Si es así, especifique el tipo
de fenómeno posible

En base a la forma de la curva ln(n) v/s L, se considera posible un alejamiento leve

de la idealidad del MSMPR.

Si tomamos en cuenta el comportamiento del grafico 3, donde se nota que los datos
están un poco dispersos, asumiendo un comportamiento lineal, el ajuste no resulta
preciso, dado que el coeficiente de correlación R2 es de 0.9024, que, si bien su dato
es aceptable, no es muy confiable, sin embargo, podemos atribuirlo a la retirada de
cristales de tamaño elevado, Elutriación (Proceso de separación de partículas más
ligeras de las más pesadas por medio de una corriente dirigida hacia arriba de gas
o líquido) o mezcla no adecuadas. Es por esto que el fenómeno que existe con los
datos entregados se puede atribuir a MSMPR con clasificación de producto.

10
 f) En un cristalizador batch se siembran 500 gramos de cristales con
tamaño de 100 mallas Tyler. Si después de 30 minutos se obtiene la
misma densidad de suspensión que en el MSMPR, determine el
volumen del cristalizador batch. Suponga que G es constante y que se
cumplen las condiciones para aplicar la ley L de Mc Cabe.

Datos
Masa cristal sembrados, Ms:500 gr
𝑔 𝑐𝑟𝑖𝑠𝑡𝑎𝑙𝑒𝑠
Densidad de suspensión, :𝑀𝑇 = 874.2339 𝐿

Tamaño malla (100#): Ls=0.147 mm

Tiempo 𝜏 = 30 𝑚𝑖𝑛 = 0.5 ℎ𝑟
𝑚𝑚
Velocidad de crecimiento: 𝐺 = 0,1556 ℎ

De acuerdo a la ley de ΔL Mc cabe, masa de producto total:

𝐿𝑝 3
𝑀𝑝 = ( ) 𝑀𝑠
𝐿𝑠

Con la ecuación de aumento de tamaño de los cristales:

Δ𝑙 = 𝐿𝑝 − 𝐿𝑠

Pero: Δ𝐿 = 𝐺 ∗ 𝜏
mm
Por lo tanto, Δ𝐿 = 0,1556 h
∗ 0.5 ℎ = 0.0778 𝑚𝑚

Como se tiene el valor del tamaño de la siembra de los cristales (Ls), se puede
determinar el valor del tamaño del producto:

𝐿𝑝 = Δ𝐿 + 𝐿𝑠 = 0.0778 + 0.147 = 0.2248 𝑚𝑚

11
Reemplazando los datos en la ley ΔL Mc cabe, se tiene que la masa de producto
es:

0.2248 3
𝑀𝑝 = ( ) ∗ 500 = 1788.16 𝑔𝑟.
0.147

𝑀𝑝
Utilizando la ecuación de densidad de suspensión: 𝑀𝑇 = 𝑉
, se despeja el volumen:

𝑀𝑝 1788.17 𝑔𝑟
𝑉= =
𝑀𝑇 874.2339 𝑔𝑟
𝐿

𝑽 = 𝟐, 𝟎𝟒𝟓𝟒 𝑳

12