Article

 INTRODUCIÓN

La transferencia de masa de aire a un grupo de partículas de naftaleno en un lecho fluidizado circulante (CFB)
se investiga mediante un enfoque de dinámica de fluidos computacional (CFD). Los resultados calculados
indican que la transferencia de masa de aire a partículas en el cúmulo se reduce debido al agrupamiento de
partículas y a los incrementos del tamaño y la temperatura de las partículas. Se analizan las influencias de la
porosidad del cúmulo, la velocidad del gas de entrada y la temperatura en la transferencia de masa de aire al
cúmulo.

 PLANTEAMIENTO DEL SISTEMA

Para simplificar, se supone que (1) las partículas de naftaleno en el grupo están uniformemente distribuidas
en el espacio; (2) el diámetro y la masa de las partículas de naftaleno permanecen constantes durante la
sublimación; (3) las temperaturas de la fase gaseosa y las partículas son constantes.

Figura 1. Racimo de naftaleno en un sistema CDF con flujo de aire

Consideraciones:

- Dado que el diámetro de la partícula es muy pequeño, el gradiente de temperatura dentro de la

partícula de naftaleno es insignificante.
- La transferencia de calor entre el gas y la partícula de naftaleno se desprecia debido a la pequeña
diferencia de temperatura entre ellos.
- El flujo de aire se toma como un flujo turbulento isotérmico tridimensional.
La porosidad del racimo se calcula con la ecuación:
D3C − nd3p
εg =
D3C
Tabla 1. Parámetros usados en la simulación.

Datos Valor
Número de esferas 65
Número de esferas en el plano x-y 19
Diámetro del racimo [μm] 688
Diámetro de partícula [μm] 100
Porosidad 0.8
Temperatura [K] 300
Densidad del aire [kg/m3] 1.124
Densidad del vapor de naftaleno [kg/m3] 4.967
Presión de saturación del vapor naftaleno 50.767
[Pa]
Peso molecular del naftaleno [g/mol] 128.174
Peso molecular del aire [g/mol] 28.966
Difusividad [m2/s] 6.5525e-6
Viscosidad dinámica [kg/m۰s] 1.9139e-5
Cte. Determinada empíricamente 1.44
Cte. Determinada empíricamente 0.09
Cte. Determinada empíricamente 1.3
Cte. Determinada empíricamente 1.92

Teniendo en cuenta estas suposiciones, las ecuaciones de conservación de masa, momento, especie y
energía en el modelo CFD 3D se dan a continuación:
- Ecuación de continuidad de la fase gaseosa
La ecuación de conservación de la masa del gas se puede escribir de la siguiente manera:
∂ρg ∂
+ (ρ u ) = 0
∂t ∂xj g j (1)

Donde ρg y uj son la densidad y la velocidad de la fase gaseosa.

- Ecuación del momentum de la fase gaseosa

∂ ∂ ∂p ∂ ∂uj ∂ ∂ui (2)
(ρg uj ) + (ρg ui uj ) = + (μg )+ (μg ) + ρg g
∂t ∂xj xi ∂xj ∂xi ∂xj ∂xj
 Donde μg = μ1 + μt es la viscosidad efectiva de la fase gaseosa, μt = Cμ ρg k2 /ε es viscosidad
turbulenta de la fase gaseosa.
- Ecuación de la energía cinética turbulenta en fase gaseosa
∂ ∂ ∂ μg ∂k
(ρg k) + (ρg ui k) = ( ) + G k − ρg ε (3)
∂t ∂xj ∂xj σk ∂xj

Donde el término de producción (𝐺𝑘 ) es:

∂ui ∂uj ∂ui 2 ∂uk
G k = μg ( + ) − δ
∂xj ∂xi ∂xj 3 ij ∂xk (4)

- Ecuación tasa de disipación de energía cinética turbulenta en fase gaseosa

∂ ∂ ∂ μg ∂ε ε
(ρg ε) + (ρg uj ε) = ( ) + [C1 Gk − C1 ρg ε] (5)
∂t ∂xj ∂xj σε ∂xj k

- Especies en equilibrio (aire y racimo)

La ecuación de equilibrio para las especies k (aire y el cúmulo) se puede escribir como (Fan y Zhu,
1998):
∂ ∂ ∂ μt ∂Yk
(ρg Yk) + (ρg uj Yk) = [(ρg Dk + ) ] (6)
∂t ∂xj ∂xj Sct ∂xj

Donde 𝑆𝑐𝑡 es el número turbulento de Schmidt con un valor constante de 0,7

- Condiciones iniciales y de frontera
Inicialmente, la velocidad del aire es cero y no hay naftaleno en el racimo. Se asume una velocidad
uniforme del aire de entrada. En la salida, las condiciones de frontera son:
∂ug,y ∂ug,x ∂ug,z ∂k ∂ε (7)
= = = 0; = 0; =0
∂y ∂x ∂z ∂y ∂y

Se emplea una condición de frontera sin deslizamiento para el flujo de gas en la superficie de la
partícula (Fan y Zhu, 1998):
u=k=ε=0 (8)

Se especifica una temperatura constante para la superficie de la partícula de naftaleno 𝑇𝑤 .

La fracción másica del naftaleno en la interface se estimó a partir de la presión de vapor del naftaleno,
que es la condición límite para la Ec. (6) en la superficie de las partículas. La presión de vapor del
naftaleno se determina usando la ecuación (Soguin, 1958; Weast, 1973):
3729.4 (9)
lg Pn = 13.564 −
Tw
 Concentración de naftaleno
Pn (10)
Yw =
Rn Tw

La densidad de vapor de naftaleno se estima usando la presión de vapor del naftaleno y la ecuación
de estado de gas ideal. La correlación empírica para el coeficiente de difusión de Reid, Prausnitz y
Poling (1987), presentan:
1.75 (11)
−6 (
Tg
D = 5.13x10 )
273.15

 PROCEDIMIENTO DE SIMULACIÓN

- Las ecuaciones (1) y (6) se resuelven utilizando el método de volumen finito con las condiciones límite
apropiadas.
- La solución de la ecuación de balance de cantidad de movimiento de la fase gaseosa necesita adivinar o
evaluar el gradiente de presión.
- Se ha utilizado el método semi-implícito para ecuaciones vinculadas a la presión (SIMPLE) para la
solución de la ecuación del momento de la fase gaseosa, adivinando el gradiente de presión y logrando la
convergencia a través de la iteración (Patankar, 1981).

 RESULTADOS Y ANÁLISIS

- Efecto de la porosidad y Reynolds

Se modelan dos velocidades de entrada respectivas de 0.783 y 1.0 m/s, correspondientes a los números
de Reynolds de 85.7 y 109. Con el aumento de la porosidad del racimo, se incrementa la relación de flujo
de aire. La relación del flujo de aire que pasa a través del cúmulo aumenta exponencialmente después de
que la porosidad del cúmulo sea superior a 0.8.

Figura 2. Relación entre el flujo de gas a través del cúmulo y el flujo total de gas en función de la porosidad:
u=1.0, 2.0, 3.0 y 5,0 m/s
- Efecto de la velocidad inicial
A baja velocidad del aire de entrada, las concentraciones de vapor de naftaleno no muestran cambios
aparentes en el grupo ya que el aire no puede ingresar al grupo. Sin embargo, con el aumento de la
velocidad del aire de entrada, puede pasar más aire a través del grupo y provocar una disminución de
la concentración de vapor de naftaleno. Sin embargo, se puede esperar la rotura de racimos cuando la
velocidad del aire es lo suficientemente alta.

Figura 3. Concentraciones molares promediadas en el tiempo calculadas de vapor de naftaleno en el plano x-y

La siguiente ecuación se utiliza para el cálculo de la tasa de transferencia de masa (Bird, Stewart y
Lightfoot, 1960):

(12)

ΔYn es la diferencia media logarítmica de la concentración másica de vapor de naftaleno:

(13)

El flujo másico del vapor de naftaleno es:

Donde el coeficiente de transferencia de masa es calculado con la Ec. (14).

(14)

El número de Sherwood correspondiente se define como:

 El número de Sherwood para la partícula aislada es tres veces mayor que el del grupo, lo que indica
que el grupo afecta negativamente a la transferencia de masa de un CFB.

Figura 4. Número de Sherwood del grupo y la partícula aislada en función de la velocidad del aire de entrada.

- Efecto de la temperatura

Las concentraciones de vapor de naftaleno están fuertemente relacionadas con la temperatura de las
partículas. En general, las temperaturas más altas hacen que se evapore más naftaleno en la corriente
de aire. La alta temperatura del aire, hace que el aire sea más viscoso y da como resultado un menor
flujo de aire a través del grupo, por lo tanto, los coeficientes de transferencia de masa hidrodinámica
disminuyen con el incremento de temperatura con una alimentación de aire constante al sistema.

Figura 5. Distribución del coeficiente de transferencia de masa en función de la temperatura.

- CONCLUSIONES

- Los resultados numéricos indican que la tasa de transferencia de masa del cúmulo de partículas en un
CFB es menor que la de las partículas dispersas en las mismas condiciones de operación.
- Los coeficientes de transferencia de masa del gas al cúmulo aumentan con el aumento de la porosidad
del cúmulo y la velocidad de aire de entrada.
Nota: Se debe tener en cuenta que los análisis se basan en una suposición isotérmica con el cúmulo de
partículas en estado estacionario en la corriente de aire y la ruptura del cúmulo no se simula. También existe
la suposición de que el cúmulo tiene la forma de un objeto circular. Esta es una buena base para el inicio de
las simulaciones, pero ciertamente es una aproximación.