Diffusioosmosis of electrolyte solutions in axisymmetric channels

Resumen Técnico: Difusioósmosis de Soluciones Electrolíticas en Canales Axisimétricos

Planteamiento del Problema
El artículo aborda el desafío teórico de describir el flujo de difusioósmosis en estado estacionario de soluciones electrolíticas en canales axisimétricos largos (por ejemplo, poros cilíndricos, cónicos o periódicamente modulados) sujetos a gradientes de concentración y presión. Si bien la difusioósmosis en rendijas planas con solutos neutros está bien comprendida, la teoría para canales cargados que contienen electrolitos sigue estando poco desarrollada. En tales sistemas, el flujo es impulsado no solo por gradientes de concentración de solutos, sino también por campos eléctricos tangenciales no uniformes (contribuciones electroosmóticas) derivados de la difusión diferencial de cationes y aniones. Los tratamientos analíticos previos a menudo asumieron gradientes de concentración constantes o se centraron únicamente en geometrías planas, sin tener en cuenta el acoplamiento no lineal entre el flujo de fluido, el flujo iónico y las secciones transversales variables del canal, presentes en medios porosos y dispositivos microfluídicos realistas.

Metodología
Los autores desarrollan una teoría analítica no lineal basada en las ecuaciones acopladas de Nernst-Planck, Poisson y Stokes. El modelo considera un canal axisimétrico cargado de longitud $L$ y radio local $R(X)$ que conecta reservorios de diferentes salinidades ($C_0$ y $C_1$) y presiones hidrostáticas ($P_0$ y $P_1$).

Características metodológicas clave:

• Aproximación de Doble Capa Delgada: La teoría asume que las capas difusas electrostáticas son delgadas en comparación con el radio local del canal ($\lambda_D \ll R$), lo que permite tratar el canal como electro-neutro en la región "bulk" con condiciones de contorno de deslizamiento (slip) en las paredes.
• Formulación Adimensional: Las ecuaciones gobernantes se adimensionalizan para derivar expresiones de forma cerrada para el flujo de fluido global ($Q$) y el flujo iónico ($J$).
• Solución Auto-consistente: A diferencia de las teorías lineales, los autores resuelven el perfil de concentración $c_m(x)$ y el potencial $\psi_m(x)$ de manera auto-consistente, reconociendo que las fuerzas impulsoras locales dependen del flujo de fluido global y viceversa.
• Modelado Geométrico: Las secciones transversales variables se modelan mediante una familia de curvas $r(x) = (1 + A\omega \sin(\omega x))^{-1/2}$, que abarca geometrías en expansión, contracción y moduladas sinusoidalmente.
• Aproximación de Cilindro: Para interpretar los resultados de las secciones transversales variables, los autores introducen una aproximación de "cilindro efectivo", relacionando la resistividad hidrodinámica de formas complejas con un radio cilíndrico equivalente.

Contribuciones Clave y Resultados

1. Tasa de Flujo en Cilindros vs. Rendijas: Los autores demuestran que la magnitud de la tasa de flujo de difusioósmosis pura ($Q_{DO}$) en un canal cilíndrico es idéntica a la de una rendija plana de espesor equivalente. Sin embargo, la respuesta a una caída de presión aplicada ($\Delta p$) difiere cuantitativamente; la contribución impulsada por presión al flujo total es menor en cilindros que en rendijas debido a factores geométricos.
2. Acoplamiento No Lineal: El estudio confirma que, para caídas de concentración grandes, la relación entre la tasa de flujo y las fuerzas impulsoras es no lineal. El flujo de fluido global $Q$ y el flujo iónico $J$ están acoplados a través del perfil de concentración local en el bulk, lo que impide una representación de matriz lineal simple de la movilidad.
3. Efecto de la Geometría en el Flujo:
   • Canales en Expansión: Los canales que se ensanchan en la dirección del flujo (hacia el reservorio salino) aceleran el flujo de difusioósmosis.
   • Canales en Contracción: Los canales que se estrechan en la dirección del flujo retardan el flujo de difusioósmosis.
   • Linealidad de $\Delta p(Q)$: Se encuentra que la relación entre la caída de presión y la tasa de flujo es casi lineal para todas las geometrías probadas, aunque la pendiente (resistividad hidrodinámica) es altamente sensible a la forma del canal.
4. Sintonización del Flujo Iónico: El artículo deriva una ecuación que vincula el flujo iónico con la tasa de flujo total. Muestra que, si bien el signo del flujo de difusioósmosis es generalmente hacia el reservorio de agua dulce (igualando concentraciones), el flujo iónico neto puede sintonizarse para ser positivo o negativo dependiendo de la geometría del canal y la caída de presión aplicada. Específicamente, la geometría puede desplazar el punto donde el flujo iónico neto es cero.
5. Validez de la Aproximación de Cilindro: Los autores proponen una aproximación analítica simple donde un canal complejo es reemplazado por un cilindro imaginario con un radio efectivo $R_{eff}$. Este radio se determina por la resistividad hidrodinámica ($W$) del canal. Esta aproximación predice con precisión las curvas $\Delta p(Q)$ y la caída de presión de flujo cero, proporcionando una herramienta práctica para interpretar datos.

Significancia y Reivindicaciones
El artículo afirma proporcionar un marco para interpretar datos experimentales y numéricos sobre la difusioósmosis en geometrías axisimétricas, que son más representativas de materiales reales (por ejemplo, membranas porosas, poros biológicos) que las rendijas planas.

Los autores enfatizan que su análisis revela cómo la geometría del canal actúa como un parámetro de control para las tasas de flujo de fluido, ya sea potenciando o retardando el flujo en comparación con un cilindro estándar. Concluyen que, si bien el flujo de difusioósmosis en canales gruesos tiende inherentemente a igualar las concentraciones de sal (limitando su uso directo para la desalinización en este régimen), la capacidad de sintonizar el flujo iónico y la tasa de flujo mediante la forma del canal y la caída de presión ofrece un potencial significativo para el diseño de dispositivos micro y nanofluídicos. El trabajo extiende los hallazgos previos de las rendijas planas a los sistemas axisimétricos, ofreciendo una base teórica más realista para aplicaciones en tecnologías de separación, recolección de energía y biosensado.