Kinetics of Salt Creeping on a Free Surface: From Nucleation to Saturation

Título: Cinética del "Arrastre" de Sal en una Superficie Libre: Desde la Nucleación hasta la Saturación

1. El Problema

El fenómeno del "arrastre" de sal (salt creeping) consiste en el crecimiento espontáneo de cristales de cloruro de sodio (NaCl) que ascienden por superficies verticales parcialmente sumergidas en una solución salina, desafiando la gravedad. Este proceso es crucial tanto por sus implicaciones negativas (degradación y corrosión acelerada de estructuras en entornos costeros) como por sus aplicaciones positivas (desalinización impulsada por energía solar).

A pesar de ser estudiado desde hace más de un siglo (Washburn, 1900s), la dinámica completa de este fenómeno, especialmente la transición desde la nucleación inicial hasta la saturación final, no está totalmente comprendida. Los desafíos principales incluyen:

• La complejidad de los mecanismos microscópicos que impulsan el ascenso (capilaridad vs. adhesión).
• La falta de un modelo unificado que describa las tres etapas cinéticas observadas macroscópicamente.
• La influencia de parámetros ambientales como la humedad relativa (HR) y la temperatura en la morfología y velocidad de crecimiento.

2. Metodología

Los autores combinaron un enfoque experimental riguroso con modelado teórico analítico y simulaciones numéricas.

• Configuración Experimental:

  • Se utilizó una varilla de vidrio parcialmente sumergida en una solución saturada de NaCl dentro de un tanque controlado.
  • Se varió la humedad relativa (25%, 26%, 30%, 32%) y se mantuvo la temperatura ambiente.
  • Se midió el crecimiento de la altura y el espesor del depósito cristalino mediante imágenes tomadas cada 90 segundos y análisis de imagen (ImageJ).
  • Se realizaron observaciones microscópicas en portaobjetos planos para analizar la distribución de tamaños de cristal.
  • Se midió la pérdida de masa de una solución de referencia para estimar el flujo de evaporación ($J \approx 0.56 \times 10^{-4} \, \text{kg} \cdot \text{m}^{-2} \cdot \text{s}^{-1}$).

• Modelo Teórico:

  • Se desarrolló un modelo analítico basado en la conservación de flujo en la interfaz cristal-aire.
  • Se asumió que la evaporación de una molécula de agua induce la cristalización de una fracción de moléculas de NaCl.
  • El flujo de evaporación se modeló con tres parámetros clave: el flujo base ($\Phi_0$), un factor de obstrucción por acreción ($\xi$) y una altura característica ($h_\infty$) donde la gravedad limita el transporte capilar.
  • Se aplicó el teorema de transporte de Reynolds para derivar ecuaciones diferenciales para la altura $h(t)$ y el espesor $e(z,t)$.

• Simulaciones Numéricas:

  • Se empleó un autómata celular bidimensional para simular el crecimiento mesoscópico.
  • El modelo incorporó la evolución de la concentración de sal, la evaporación y una probabilidad estocástica de cristalización basada en un factor de Boltzmann, considerando efectos de tensión superficial (relación Ostwald-Freundlich).

3. Contribuciones Clave

El trabajo presenta una descripción unificada de la cinética del arrastre de sal, identificando y modelando matemáticamente tres regímenes cinéticos distintos:

1. Regimen Inicial (Crecimiento Exponencial): Dominado por la nucleación y la expansión libre del cristal. La altura crece exponencialmente mientras el espesor permanece constante.
2. Regimen Intermedio (Crecimiento Lineal): Una fase de larga duración donde la altura aumenta linealmente con el tiempo. El espesor sigue siendo constante, pero la altura se vuelve comparable a la altura característica $h_\infty$.
3. Regimen de Saturación (Estancamiento de Altura y Crecimiento Logarítmico): La altura deja de crecer (satura) debido a la gravedad que impide el ascenso de la salmuera, mientras que el espesor del depósito aumenta logarítmicamente con el tiempo.

Además, el estudio valida que la geometría del sustrato (cilíndrico vs. plano) no altera el comportamiento macroscópico global, y demuestra que las simulaciones pueden reproducir la distribución de tamaños de cristal observada experimentalmente.

4. Resultados

• Validación Experimental del Modelo:

  • Los datos experimentales de altura y ancho se ajustaron con éxito a las ecuaciones derivadas teóricamente (Ecuaciones 5, 7 y 9).
  • Regimen 1: El tiempo característico $\tau_1$ predicho teóricamente ($\approx 4$ horas) coincide en orden de magnitud con la observación experimental ($\approx 2.4$ horas).
  • Regimen 2: Se observó un crecimiento lineal de la altura. La pendiente experimental ($0.8 \pm 0.013$ cm/h) es del mismo orden que la predicción teórica ($3$ cm/h), aunque la discrepancia se atribuye a fluctuaciones en la tasa de evaporación y ruido experimental.
  • Regimen 3: El ensanchamiento lateral del depósito sigue una tendencia logarítmica, confirmando la predicción de que, una vez saturada la altura, el crecimiento se vuelve lateral.

• Observaciones Microscópicas y Simulaciones:

  • La distribución de tamaños de cristal obtenida en simulaciones (promedio $\approx 3.9 \, \mu m$) coincide satisfactoriamente con las mediciones experimentales ($\approx 3.1 \, \mu m$).
  • Efecto de la Temperatura: Las simulaciones revelaron que un aumento de temperatura reduce el tamaño promedio de los cristales y la prevalencia de cristales grandes. Esto sugiere que la temperatura afecta tanto la cinética como la morfología final, mientras que la humedad relativa afecta principalmente la velocidad de cristalización.
  • Mecanismo Propuesto: Cristales más pequeños (inducidos por mayor temperatura) crean espacios capilares más estrechos. Según la ley de Jurin, esto permitiría un ascenso mayor, pero al mismo tiempo aumenta la resistencia al flujo radial (aumentando $\xi$), limitando el ensanchamiento lateral.

• Limitaciones Identificadas:

  • La medición de masa no logró modelar con precisión el primer régimen debido a la dificultad de distinguir entre condiciones iniciales macroscópicas y mesoscópicas (mojado de la varilla).
  • La tomografía de coherencia óptica (OCT) no pudo confirmar la existencia de un capilar principal continuo entre la varilla y el cristal debido a la difusión de la luz en el medio cristalino grueso.

5. Significado e Impacto

Este estudio proporciona una comprensión fundamental de la dinámica de crecimiento de sales en superficies libres, cerrando la brecha entre la observación macroscópica y los mecanismos microscópicos.

• Científico: Establece un marco teórico robusto que explica la transición entre regímenes de crecimiento, desafiando modelos previos basados únicamente en la ley de Darcy (que fallan al no capturar la naturaleza fuera de equilibrio y la reducción de permeabilidad por cristalización).
• Práctico:
  • Protección de Infraestructuras: Al entender los tiempos característicos y los factores que limitan la altura máxima ($h_\infty$), se pueden desarrollar mejores estrategias para mitigar la corrosión en edificios costeros.
  • Tecnología de Desalinización: La capacidad de predecir la morfología y la tasa de crecimiento de los cristales es vital para optimizar sistemas de desalinización solar que dependen de la evaporación y cristalización controlada.
  • Dirección Futura: El trabajo sugiere que la temperatura es un parámetro más crítico que la humedad relativa para controlar la morfología final de los depósitos, abriendo nuevas vías de investigación para el control de procesos de cristalización.

En conclusión, el artículo logra una descripción unificada del fenómeno de "salt creeping" desde la nucleación hasta la saturación, validando un modelo teórico de tres etapas mediante experimentos controlados y simulaciones numéricas.