La difusión transitoria sólido-líquido se rige por las leyes de Fick, la conservación de la masa y la diferencia de concentración cambiante entre las fases sólida y líquida. Dentro del sólido, el soluto se desplaza por difusión hacia la superficie; la agitación transporta ese soluto a través del líquido y mantiene el líquido a granel aproximadamente uniforme. Debido a que el volumen de líquido es finito, su concentración aumenta a medida que avanza la extracción, reduciendo la fuerza impulsora hasta alcanzar el equilibrio o una recuperación especificada.

Conclusión clave: Un modelo de extracción a escala piloto debe acoplar la difusión transitoria dentro del sólido con un balance de masa dinámico en la fase líquida. La difusión interna, el tamaño de partícula, la difusividad efectiva, el equilibrio interfacial, el volumen de líquido y la agitación determinan conjuntamente la tasa de extracción y el tiempo de lote requerido.

El modelo físico gobernante

La primera ley de Fick determina el flujo local

La primera ley de Fick relaciona el flujo de soluto con el gradiente de concentración:

[ \mathbf{N}A=-D{\mathrm{eff}}\nabla C_A ]

Aquí, (D_{\mathrm{eff}}) es la difusividad efectiva del soluto (A) a través de la matriz sólida, y (C_A) es su concentración dentro del sólido.

El signo negativo indica que el soluto se mueve desde regiones de alta concentración hacia regiones de baja concentración. En una partícula porosa, (D_{\mathrm{eff}}) refleja los efectos de la estructura de los poros, la porosidad, la tortuosidad, la saturación, la temperatura y la viscosidad del líquido.

La segunda ley de Fick describe el agotamiento transitorio

Para una difusividad constante y sin reacción dentro del sólido, el campo de concentración transitorio satisface:

[ \frac{\partial C_A}{\partial t}

D_{\mathrm{eff}}\nabla^2 C_A ]

Para una partícula esférica, una geometría común en los modelos de lixiviación, esto se convierte en:

[ \frac{\partial C_A}{\partial t}

D_{\mathrm{eff}} \left[ \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2\frac{\partial C_A}{\partial r} \right) \right] ]

La geometría es importante porque las placas planas, los cilindros y las esferas tienen diferentes longitudes de trayectoria de difusión y perfiles de concentración.

La condición inicial representa el sólido cargado

El modelo comienza con una distribución inicial de soluto, a menudo aproximada como uniforme:

[ C_A(r,0)=C_{A,0} ]

Los sólidos más complejos pueden tener un perfil de concentración radial causado por secado, impregnación, reacción o precipitación previos. Esa condición inicial puede afectar sustancialmente la tasa de extracción inicial.

La simetría restringe el centro de la partícula

Para una partícula esférica, no hay flujo neto que cruce el centro:

[ \left. \frac{\partial C_A}{\partial r} \right|_{r=0}=0 ]

Esta condición evita una singularidad artificial en el centro y refleja la simetría física de la partícula.

Cómo entra el volumen líquido finito en el modelo

El líquido a granel a menudo se trata como bien agitado

Con una agitación vigorosa, la concentración del líquido a granel puede aproximarse como espacialmente uniforme:

[ C_L=C_L(t) ]

Esta suposición elimina la necesidad de resolver un campo de concentración completo en la fase líquida. No elimina la transferencia de masa; significa que los gradientes de concentración del lado del líquido se representan a través de una resistencia interfacial o se asumen suficientemente pequeños.

El balance de líquido acopla ambas fases

El soluto liberado del sólido se acumula en el volumen de líquido finito:

[ V_L\frac{dC_L}{dt}

\dot{m}_{A,\mathrm{sólido\ a\ líquido}} ]

Para (N_p) partículas idénticas con área superficial (A_p), esto se puede escribir como:

[ V_L\frac{dC_L}{dt}

N_p A_p N_A ]

donde (N_A) es el flujo de soluto hacia afuera en la superficie de la partícula.

El balance total de soluto proporciona una importante verificación de consistencia:

[ M_{A,0}

M_A^{\mathrm{sólido}}(t) + V_L C_L(t) ]

asumiendo que no hay pérdida de soluto, degradación o eliminación por muestreo.

El volumen finito crea una fuerza impulsora decreciente

Una condición interfacial comúnmente utilizada es:

[ -D_{\mathrm{eff}} \left. \frac{\partial C_A}{\partial r} \right|_{r=R}

k_L\left(C_{A,i}-C_L\right) ]

donde (C_{A,i}) es la concentración interfacial equivalente al líquido y (k_L) es el coeficiente de transferencia de masa del lado del líquido.

Si el equilibrio local se representa mediante una relación de partición,

[ C_{A,i}=K_{\mathrm{eq}}C_{A,s} ]

entonces la concentración superficial en el sólido y la concentración del líquido están vinculadas a través de (K_{\mathrm{eq}}). A medida que aumenta (C_L), la diferencia entre la concentración de equilibrio del líquido y (C_L) disminuye, por lo que la extracción se ralentiza.

Los principios adimensionales clave

El número de Fourier mide el progreso de la difusión

El tiempo de difusión adimensional es el número de Fourier:

[ \mathrm{Fo}=\frac{D_{\mathrm{eff}}t}{R^2} ]

Aquí, (R) es el radio característico de la partícula. Los valores pequeños indican que la difusión ha penetrado solo una parte del sólido, mientras que los valores más grandes indican un agotamiento más extenso.

El tiempo de difusión característico escala como:

[ t_D\sim \frac{R^2}{D_{\mathrm{eff}}} ]

Esta dependencia al cuadrado hace que el tamaño de partícula sea especialmente importante en el diseño de plantas piloto. Reducir a la mitad el radio de la partícula puede reducir el tiempo de difusión interna característico en aproximadamente un factor de cuatro, siempre que otras limitaciones permanezcan sin cambios.

El número de Biot compara la resistencia externa e interna

Un número de Biot de transferencia de masa compara la resistencia del lado del líquido con la difusión interna del sólido:

[ \mathrm{Bi}_m=\frac{k_L R}{D_{\mathrm{eff}}} ]

Un (\mathrm{Bi}_m) grande generalmente indica que la resistencia externa es relativamente pequeña en comparación con la resistencia a la difusión interna. Un (\mathrm{Bi}_m) pequeño indica que la capa límite líquida puede controlar o influir significativamente en la tasa de extracción.

La definición exacta puede variar según la geometría y las convenciones de partición, pero la pregunta subyacente es la misma: ¿dónde ocurre la mayor parte de la caída de concentración?

Las soluciones de similitud explican el comportamiento en tiempos tempranos

Para la difusión semi-infinita dentro o fuera de un sólido, la variable de similitud es:

[ \eta=\frac{x}{2\sqrt{D_{\mathrm{eff}}t}} ]

La ecuación de difusión puede entonces producir perfiles de concentración de función de error, tales como:

[ \Theta(\eta)=\operatorname{erfc}(\eta) ]

Estas soluciones muestran que la profundidad de penetración crece con la raíz cuadrada del tiempo:

[ \delta(t)\propto \sqrt{D_{\mathrm{eff}}t} ]

Son útiles para comprender la extracción temprana de sólidos grandes o capas superficiales delgadas, aunque se requieren soluciones de partícula finita una vez que el frente de difusión llega al centro de la partícula.

El flujo temprano a menudo decae como (t^{-1/2})

Para un sólido semi-infinito inicialmente uniforme expuesto a una concentración superficial fija, el flujo superficial tiene la forma:

[ N_A\propto \left(\frac{D_{\mathrm{eff}}}{\pi t}\right)^{1/2} \Delta C ]

La dependencia (t^{-1/2}) describe la disminución característica en el flujo de difusión transitoria. En un extractor de volumen finito, la disminución suele ser más fuerte o modificada de otro modo porque la concentración del líquido cambia y la diferencia de concentración disponible disminuye.

Aplicación del modelo a plantas piloto

El tiempo de extracción es una variable de diseño acoplada

Una planta piloto debe proporcionar suficiente tiempo de contacto para que la fracción objetivo de soluto salga del sólido. Ese tiempo depende de (D_{\mathrm{eff}}), el radio de la partícula, la temperatura, las propiedades del solvente, la relación líquido-sólido, el comportamiento de equilibrio y la agitación.

Una estimación simple basada solo en (R^2/D_{\mathrm{eff}}) es útil para el dimensionamiento inicial, pero no puede predecir la recuperación final cuando la acumulación de líquido y las limitaciones de equilibrio son importantes.

La agitación controla principalmente el transporte externo

La agitación reduce la capa límite líquida que rodea a las partículas y aumenta el coeficiente de transferencia de masa del lado del líquido (k_L). También puede suspender partículas, evitar la sedimentación, mejorar la uniformidad de la temperatura y reducir los gradientes de concentración en el recipiente.

Una vez que la resistencia externa se vuelve pequeña, la agitación adicional produce mejoras decrecientes porque la difusión interna del sólido se ha convertido en la resistencia controladora.

El tamaño de partícula influye fuertemente en el rendimiento

Reducir el tamaño de partícula acorta la trayectoria de difusión característica y aumenta el área superficial externa disponible por unidad de masa. Estos cambios generalmente aceleran la extracción, pero también pueden aumentar el consumo de energía, crear finos, complicar la filtración y causar problemas de manejo o seguridad.

Por lo tanto, se necesitan pruebas piloto para identificar el tamaño de partícula práctico más pequeño en lugar de asumir que la molienda máxima es óptima.

El volumen de líquido controla la fuerza impulsora

Un mayor volumen de solvente puede mantener la concentración del líquido a granel más baja para una cantidad dada de soluto extraído. Esto mantiene una fuerza impulsora mayor y puede aumentar la recuperación final cuando el equilibrio de partición limitaría de otro modo la extracción.

La contrapartida es un mayor uso de solvente, un mayor deber de concentración aguas abajo y un volumen de equipo potencialmente mayor. Por lo tanto, la relación líquido-sólido adecuada es una variable de optimización en lugar de una constante universal.

Los datos piloto conectan la teoría con el escalado

Una campaña piloto útil mide la concentración del líquido a granel en función del tiempo mientras registra el tamaño de partícula, la temperatura, la velocidad de agitación, el volumen de solvente y la carga de sólido. Estos datos pueden ajustarse a un modelo de difusión transitoria para estimar (D_{\mathrm{eff}}), (k_L), los parámetros de equilibrio o combinaciones de ellos.

La velocidad de agitación por sí sola no es un criterio de escalado transferible. El escalado debe preservar el comportamiento de transporte controlador, utilizando cantidades relevantes como el tiempo de mezcla, la entrada de potencia, la calidad de la suspensión, el número de Reynolds o coeficientes de transferencia de masa correlacionados experimentalmente.

Comprensión de las contrapartidas

La suposición de "bien agitado" tiene límites

Un líquido puede parecer vigorosamente mezclado y, sin embargo, contener gradientes de concentración locales cerca de las partículas, las paredes del recipiente o las regiones mal suspendidas. Si el volumen de líquido es grande o la suspensión es densa, una (C_L(t)) espacialmente uniforme puede ser inexacta.

Las mediciones del tiempo de mezcla, las pruebas con trazadores o los modelos de compartimentos y de dinámica de fluidos computacional pueden usarse cuando la uniformidad a granel es incierta.

La difusión interna puede no ser el único mecanismo

La tasa de extracción aparente puede incluir lavado superficial rápido, difusión en poros, hinchamiento, rotura de partículas, disolución, reacción química, adsorción o desorción. Tratar cada desviación del comportamiento fickiano ideal como un cambio en la difusividad puede producir parámetros ajustados engañosos.

Puede ser necesario un modelo de reacción-difusión o un modelo de partícula encogible cuando la estructura sólida cambia durante la extracción.

La difusividad efectiva suele ser una propiedad aparente

(D_{\mathrm{eff}}) no es necesariamente la difusividad molecular medida en un líquido puro. A menudo incorpora la geometría de los poros, la tortuosidad, la constricción, la saturación, la partición y las interacciones con la matriz sólida.

Su valor puede cambiar con la temperatura, el contenido de humedad, la concentración de soluto o el grado de agotamiento. Por lo tanto, un coeficiente ajustado a escala de laboratorio debe validarse en condiciones piloto.

El equilibrio puede limitar la recuperación final

Aumentar la agitación no puede superar una relación de equilibrio desfavorable. Si el líquido se acerca a su concentración de equilibrio, la tasa de extracción se vuelve pequeña incluso cuando la mezcla es excelente.

La operación en contracorriente, la renovación del solvente, la extracción por etapas, el ajuste del pH o una reacción química selectiva pueden mejorar la fuerza impulsora general cuando estos enfoques son compatibles con el producto y el proceso.

Los sistemas multicomponentes requieren atención adicional

Para la extracción de un solo soluto diluido, un modelo fickiano o pseudobinario suele ser adecuado. En sistemas concentrados o multicomponentes, las interacciones de las especies, la adsorción competitiva, la termodinámica no ideal y la difusión acoplada pueden requerir descripciones de transporte de Stefan-Maxwell u otras descripciones multicomponentes.

El modelo debe coincidir con el rango de composición y la precisión requerida para la decisión de la planta piloto.

Cómo aplicar esto a su proyecto

Se debe desarrollar un modelo de extracción transitoria en torno a propiedades medibles y validarse con datos de concentración frente al tiempo. Las siguientes prioridades ayudan a alinear el modelo con el objetivo de ingeniería:

Si su enfoque principal es el tiempo de extracción: Estime (R^2/D_{\mathrm{eff}}), mida las distribuciones de tamaño de partícula y ajuste el modelo completo de difusión sólida transitoria para determinar el tiempo de lote para la recuperación requerida.
Si su enfoque principal es el rendimiento de la agitación: Varíe la velocidad del impulsor o las condiciones de fluidización mientras mide la tasa de extracción, la calidad de la suspensión y el tiempo de mezcla para separar la resistencia del lado del líquido de la difusión interna.
Si su enfoque principal es el consumo de solvente: Utilice el balance de líquido de volumen finito y la relación de equilibrio para seleccionar la relación líquido-sólido o el número de etapas de extracción que logre la recuperación con un uso aceptable de solvente.
Si su enfoque principal es el escalado de piloto a producción: Determine si la difusión interna, la transferencia de masa externa, la mezcla o el equilibrio controlan el proceso, luego preserve el comportamiento adimensional controlador y valídelo a mayor escala.
Si su enfoque principal es la fiabilidad del modelo: Verifique la conservación total de la masa, mida la temperatura y la viscosidad, caracterice la porosidad y el tamaño de partícula, y pruebe si (D_{\mathrm{eff}}) y (k_L) siguen siendo válidos en el rango operativo.

La tarea central de ingeniería es identificar la resistencia controladora y diseñar la planta piloto de modo que sus mediciones distingan los efectos de difusión, mezcla y equilibrio.

Tabla resumen:

Principio	Rol en el diseño de extracción
Primera ley de Fick	Calcula el flujo de soluto a partir de gradientes de concentración.
Segunda ley de Fick	Predice el agotamiento transitorio de soluto dentro de las partículas.
Balance de líquido finito	Rastrea el aumento de la concentración a granel y la disminución de la fuerza impulsora.
Número de Fourier	Relaciona el progreso de la difusión con el tiempo y el tamaño de partícula.
Número de Biot de transferencia de masa	Compara la resistencia externa del líquido con la difusión interna.
Radio de partícula	Afecta fuertemente el tiempo de extracción, aproximadamente proporcional a R^2/D_eff.
Agitación	Aumenta la transferencia de masa del lado del líquido, la suspensión y la uniformidad de la temperatura.
Relación de equilibrio	Determina la recuperación máxima alcanzable.
Datos piloto	Valida D_eff, k_L, parámetros de equilibrio y suposiciones de escalado.

Mejore su planta piloto de extracción

¿Necesita determinar el tiempo de extracción, el volumen de solvente, los requisitos de agitación o las condiciones de escalado? Contacte a LABPARK para discutir una unidad piloto de operaciones unitarias educativa o vocacional adaptada a su proyecto.

LABPARK ayuda a universidades, institutos de investigación y empresas a evaluar la extracción sólido-líquido utilizando equipos prácticos a escala piloto para ingeniería química, bioprocesos y biotecnología, y aplicaciones de tratamiento de agua y medio ambiente. Nuestros sistemas admiten estudios controlados de tamaño de partícula, temperatura, relación líquido-sólido, agitación, mezcla, transferencia de masa y equilibrio para que su equipo pueda identificar la resistencia controladora y generar datos fiables de concentración frente al tiempo.

Trabaje con LABPARK para conectar la teoría de la difusión transitoria con el rendimiento medible del proceso, mejorar la eficiencia del solvente y la energía, validar los parámetros de escalado y seleccionar las condiciones operativas que logren su recuperación objetivo.