Resumen: Contexto: Dentro de las tecnologías de procesamiento con membranas, la evaporación osmótica es una alternativa promisoria para la transformación de frutas exóticas, generando productos concentrados que pueden ser usados en la alimentación diaria, siendo más fáciles de consumir, disminuyendo gastos de transporte e incrementando la vida útil. Método: En este trabajo de investigación se estudió y desarrollo una estrategia integral para la modelación y simulación multifísica de los fenómenos de transferencia de masa y movimiento en el proceso de evaporación osmótica, por medio del software Comsol® y Matlab®, usando un enfoque de geometría axial en dos dimensiones como simplificación del módulo real y el método de elementos finitos para la solución numérica, las simulaciones fueron validadas experimentalmente en un sistema de evaporación osmótica a escala de laboratorio. Resultados: Los modelos utilizados y las simulaciones generadas fueron estadísticamente significativas (p<0,05) en predecir el comportamiento del flux, teniendo en cuenta el efecto de la velocidad y temperatura de alimentación, junto con la velocidad de la salmuera, obteniéndose correlaciones mayores al 96% entre los datos experimentales y los calculados. Conclusiones: Se encontró que para las condiciones estudiadas el modelo difusional Knudsen es el más adecuado en describir la transferencia del vapor de agua a través de la membrana hidrófoba; las simulaciones desarrolladas describen de forma adecuada el proceso de evaporación osmótica, convirtiéndose en una herramienta para el desarrollo más rápido y económico de esta tecnología
Palabras clave:difusión, membrana hidrófoba, transferencia de calor, transferencia de masa. difusión, membrana hidrófoba, transferencia de calor, transferencia de masa. .
Abstract: Context: Within the processing technologies with membranes, osmotic evaporation is a promising alternative for the transformation of exotic fruits, generating concentrated products that can be used in the daily diet, being easier to consume, reducing transportation costs and increasing shelf life. Method: In this research, it was studied and developed a comprehensive strategy for multiphysics modeling and simulation of mass and momentum transfer phenomena in the process of osmotic evaporation through Comsol® and Matlab® software. It was used an axial geometry approach in two dimensions as simplifications of real module and the finite element method for the numerical solution. The simulations were validated experimentally in an osmotic evaporation system of laboratory scale. Results: The models used and the generated simulations were statistically significant (p <0,05) in predicting the flux behavior, taking into account the effect of flow and temperature feed together with the brine flow, being obtained correlations above 96% between experimental and calculated data. Conclusions: It was found that for the conditions studied the Knudsen diffusion model is most suitable to describe the transfer of water vapor through the hydrophobic membrane. Simulations developed adequately describe the process of osmotic evaporation, becoming a tool for faster economic development of this technology
Keywords: diffusion, hydrophobic membrane, heat transfer, mass transfer..
Investigación
Modelación y simulación computacional del proceso de evaporación osmótica
Modeling and computational simulation of osmotic evaporation process
Universidad Distrital Francisco José de Caldas
Recepción: 10 Octubre 2014
Aprobación: 15 Febrero 2016
Las tecnologías de membranas tienen ventajas sobre otros procesos de conservación, siendo la más importante, que las características del producto se mantienen mejor, pues se utilizan temperaturas bajas, además existen requerimientos de energía bajos, reducidas exigencias de espacio y flexibilidad para operar. Dentro de esta gama de operaciones con membranas, la concentración por evaporación osmótica (EO) tiene potencial para ser usada en las industrias cosmética, alimentaria, farmacéutica y química, esta técnica ha generado mucho interés, en el campo de los concentrados de jugos de frutas, vegetales, café instantáneo, leche, té y otros materiales sensibles al calor, (Shaw et al., 2001; Vaillant et al., 2012; Valdés, Romero, Saavedra, Plaza, & Bubnovich, 2013), debido a que se realiza a temperatura ambiental, presión atmosférica, y condiciones isotérmicas (Bailey, Barbe, Hogan, Johnson, & Sheng, 2012; Nii, Jebson, & Cussler, 2002) disminuyendo el pardeamiento no enzimático, pérdida de aromas, degradación de color y sabor.
En la concentración osmótica el transporte de materia (figura 1a) se divide en tres fases: las etapas inicial y final son el movimiento de agua, desde el producto diluido a través de la interface y desde la zona de condensación hacia la solución concentrada, la etapa intermedia es cuando el vapor pasa a través de la membrana (Romero, 2003; Ratchaneekorn Thanedgunbaworn, Jiraratananon, & Nguyen, 2009). El diferencial de presión en la membrana, se produce con soluciones salinas como el MgSO4, MgCl2, CaCl2, NaCl, y líquidos como el glicerol, que tienen alta solubilidad, baja actividad de agua y alta tensión superficial. El sistema básico para describir la transferencia de masa, está dado en la ecuación (1) donde se relaciona el flux (N) y el diferencial de presiones, por medio de una constante (K) que es la permeabilidad de la membrana (Alves & Coelhoso, 2004; Cassano et al., 2003).
(1)
(2)Como las condiciones existentes en la interface no siempre están disponibles, una representación más completa del proceso está dada por la ecuación , donde el coeficiente global K integra varias resistencias a la transferencia de masa (figura 1b).
Figura 1. a)
Principio del proceso
elaboración propia.
Figura 1. b)
Perfil de concentraciones y resistencia a la transferencia de materia en la EO.
elaboración propia.
Puesto que en los poros de la membrana existe aire, sólo dos formas de difusión participan en la transferencia de agua, difusión Knudsen y molecular, según la teoría cinética de los gases (R. Thanedgunbaworn, Jiraratananon, & Nguyen, 2007a). El número de Knudsen (Kn) definido por la ecuación se usa para estimar cuál de los mecanismos de difusión es el predominante, comparando la trayectoria media libre (l) de la molécula (ecuación 4) que se difunde con el radio de los poros en la membrana (Qtaishat, Matsuura, Kruczek, & Khayet, 2008).
(3)
(4)Para un tamaño de poro relativamente pequeño, Kn ³ 10, las moléculas tienden a colisionar frecuentemente con las paredes del poro y el modelo Knudsen representado por las ecuaciones y (6) es el más conveniente.
(5)
(6)Cuando los poros son grandes, Kn £ 0.01, las moléculas de gas colisionan más frecuentemente entre sí, la difusión molecular es considerada predominante (Celere & Gostoli, 2002) y el flujo de vapor puede ser descrito por las ecuaciones , (8) y (9) donde PAlm es la presión media logarítmica dentro de los poros y D (coeficiente difusional) es función de la temperatura y la presión .
(7)
(8)
(9)Finalmente, en la región de transición, 0.01 < Kn < 10, los dos fenómenos se combinan y en este caso la mejor opción es utilizar un modelo mixto como el de la ecuación (10) con una permeabilidad del tipo planteada por la ecuación (11), en la cual se incluye el término PAlm que tiene en cuenta el efecto del aire presente en los poros.
(10)
(11)El flux de vapor en las zonas de polarización del lado de la alimentación y del permeado, se expresa por las ecuaciones y , algunas investigaciones desprecian esta polarización cuando la alimentación es solo agua, y en condiciones donde se cuantifican otros parámetros de la membrana (Mathilde Courel, Dornier, Herry, Rios, & Reynes, 2000; M. Courel, Tronel-Peyroz, Rios, Dornier, & Reynes, 2001).
(12)
(13)Al tener en cuenta la resistencia por polarización, las concentraciones en las capas límites son calculadas por el modelo presentado en las ecuaciones y , donde no se tiene en cuenta el cambio de la distribución de la capa en el módulo.
(14)
(15)En módulos de fibras huecas, se aplican modelos semi empíricos para estimar el coeficiente ka, con lo cual el flujo de alimentación al interior de la fibra se asimila al de una tubería, con la ecuación de Sieder – Tate y las respectivas para transferencia de materia, ecuaciones y (Martínez-Díez, Florido-Díaz, & Vázquez-González, 2000).
(16)
(17)Varios estudios han generado modelos para estimar el coeficiente (kp) al exterior de las fibras (carcasa) (Gawronski & Wrzesinska, 2000; Lipnizki & Field, 2001; Wu & Chen, 2000). Recientemente (R. Thanedgunbaworn, Jiraratananon, & Nguyen, 2007b) se ha desarrollado una nueva ecuación , con la cual se logra describir mejor el fenómeno, hallando que el número de Reynolds es función del empaquetamiento (θ), tal como se presenta en la ecuación (19).
(18)
(19)El módulo de fibras huecas (MFH) constituye la parte principal del sistema a modelar, cada fibra se compone de un cilindro con radio (r) y la longitud (z) rodeado por una membrana porosa que actúa de forma independiente, la alimentación (agua) fluye con un comportamiento laminar por dentro de la fibra (lumen), siendo introducido a z = 0, mientras que la solución de extracción (salmuera) es pasada a través de la carcasa del módulo, provocando la remoción de agua en forma de vapor desde la alimentación y que por subsecuente difusión a través de la membrana es absorbida por la salmuera. La geometría del módulo fue simplificada tal como se muestra en la figura 2, con el fin de facilitar el análisis y reducir el consumo de recursos computacionales, por lo que se utilizó un enfoque de simetría axial, donde se aprecian tres subdominios lumen, membrana y carcaza, que pueden ser identificados también como agua, membrana y salmuera.
Figura 2.
Geometría general del módulo de fibras huecas (Izq), dominios y límites para la simetría de la fibra (Der).
elaboración propia.
El modelo de la superficie libre de Happel (Chen, Gao, Jin, & Zheng, 2011) fue usado para estimar el radio r3 en el lado de la salmuera, de acuerdo a este enfoque, se considera que cada fibra en el espacio de la carcasa está rodeada por un fluido envolvente, en el cual no existe transferencia de momento, masa o calor en la superficie externa. Las fibras se consideran homogéneamente distribuidas y el flujo es puramente axial; el modulo es considerado como una serie de celdas de superficies libres, cada celda tiene solo una fibra en el centro y está rodeada de fluido uniforme, matemáticamente la superficie libre es definida por las ecuaciones (20) y (21).
(20)
(21)Donde Φ es la fracción volumétrica de vacío y puede ser calculada usando la ecuación , para cual n es el número de fibras y rim es el radio interior del módulo.
La ecuación (22) representa el estado estacionario para el transporte de vapor de agua en el lado del lumen (C1) y puede ser escrita como:
(22)Donde Dw es el coeficiente de auto difusión del agua (Holz, Heil, & Sacco, 2000) y es descrito por la ecuación (23):
(23)La ecuación es la expresión para el transporte de agua (C2) al interior de la membrana, el cual es considerado debido solamente a la difusión.
(24)En los poros de la membrana existe aire que viene mezclado con la alimentación y solo la difusión Knudsen (Dk) y molecular (Dm), están involucradas en la transferencia de vapor; para tener encuentra estas dos contribuciones y obtener el valor de Dwa que representa el coeficiente de difusión efectivo del vapor de agua en el aire retenido dentro de los poros, se hizo uso de la aproximación de Bonsaquet representada por las ecuaciones (26) y (27) (Rezakazemi, Shirazian, & Ashrafizadeh, 2012):
(25)
(26)
(27)En las anteriores ecuaciones, ε, τ, dp, son la tortuosidad, porosidad y diámetro de poro en la membrana, Mw es la masa molar del agua, R constante de los gases, T, temperatura absoluta. La ecuación (28) calcula la transferencia del vapor de agua en estado estacionario para el lado de la carcasa (C3) del contactor en coordenadas cilíndricas, se obtiene utilizando la ley de difusión de Fick:
(28)La difusión del agua en soluciones salinas puede ser vista como la auto difusión del agua retardada por la presencia de los iones en solución, el coeficiente de auto difusión del agua presenta un máximo valor en el límite de la dilución extrema y este tiende a cero en el límite de la sal anhidra (Suárez-Iglesias, Medina, Pizarro, & Bueno, 2008). El coeficiente de difusión del agua en la salmuera (Dws) puede ser calculado con la ecuación , donde Dw es determinado con la expresión y xs es la fracción de sal en la solución.
(29)A la entrada de una fibra, se asume que el flujo está completamente desarrollado bajo un régimen laminar y se especifica un flujo parabólico según la ecuación (30) y a la salida de la fibra, la presión es la atmosférica, según se especifica en la ecuación (31).
(30)
(31)La línea central interna de la fibra es una simetría axial, definida por la ecuación (32)
(32).En la superficie de la membrana, la velocidad axial (vz) es establecida como cero o condición de no deslizamiento, especificada por la ecuación (33). A la entrada de una fibra, la concentración inicial de agua en la alimentación (34) es CJ0, la cual corresponde a la inicial.
(33)
(34)La condición de frontera para el flujo convectivo asume que la masa que pasa por la frontera z = L es dominada por la convección. En otras palabras, se supone que el flujo de masa debido a la difusión a través de este límite es cero (35), la línea central de la fibra es una simetría axial (36):
(35)
(36)En la superficie de la membrana para la concentración (C1), se establece la siguiente condición especial de flux (37), que es descrita más en detalle en las condiciones límite de la membrana.
(37)Debido a que existen discontinuidades del perfil de concentración en los límites entre los líquidos y la membrana, se deben utilizar tres variables separadas para describir la concentración en las respectivas fases; para obtener un flux continuo a través de los límites, se aplica una condición de frontera especial llamada el método “stiff-spring”, en lugar de definir las condiciones de concentración tipo Dirichlet de acuerdo con el coeficiente de partición (K), lo que destruiría la continuidad del flujo; por medio de las ecuaciones (38) y (39) podemos definir condiciones de flux continuo que, al mismo tiempo, fuerzan las concentraciones a los valores deseados:
(38)
(39)Donde los coeficientes ka y kp se calculan a partir de las ecuaciones y . Finalmente, con la ecuación (40) se establece que no hay transferencia de masa en ambos bordes de la membrana:
(40)Para el flujo de salmuera sobre la fibra, el modelo de superficie libre de Happel se usa para caracterizar el perfil de velocidad fuera de las fibras (41). Aunque el flujo en el módulo real de fibras huecas no está absolutamente de acuerdo con el modelo de Happel, este se ha utilizado ampliamente para contactores de membrana (Zhang, Huang, Chi, & Pei, 2012).
(41)A la salida de la carcasa, la presión es la atmosférica y se especifica por la ecuación (42), la línea central entre las fibras se consideró una simetría axial (43).
(42)
(43)En la superficie de la membrana del lado de la salmuera, la velocidad axial (vz) es establecida como cero y siguiendo la ecuación (44), en la entrada de la carcasa, se toma la concentración de agua (C3) presente en el tanque de salmuera, convenientemente convertida en unidades de presión de vapor (45).
(44)
(45).La condición de frontera para el flujo convectivo asume que la masa que pasa por la frontera z = L es dominada por la convección, se supone que el flujo de masa debido a la difusión a través de este límite es cero (46), la línea central de la fibra es una simetría axial (47).
(46)
(47)En la superficie de la membrana (lado salmuera) para la concentración (C3), se establece la también la condición especial de flux, que se describió en las condiciones de la membrana (48).
(48)Los modelos matemáticos que describen las propiedades físicas de los fluidos involucrados en el proceso de evaporación osmótica se resumen en la tabla 1.
Las ecuaciones del modelo se resolvieron numéricamente utilizando un enfoque multifisico (Mendieta Reyes, Montaña Parra, & Noguera Vega, 2012; Puin Ávila, Noguera Vega, & Caicedo Morera, 2015) por medio del software COMSOL Multiphysics® 4.2b, conjuntamente con Matlab® por medio de la herramienta LiveLink®, este software emplea el método de elementos finitos (FEM) por sus siglas en inglés, se usó el solucionador PARDISO, la ejecución y resolución se llevó a cabo en un computador personal Intel I7 (4 GHz - 16 GB de RAM), el tiempo de cálculo fue de aproximadamente 20 minutos para cada experimento, el número máximo de iteraciones se fijó en 250 y la tolerancia relativa era 1E- 6, un factor de escala de 50 fue empleado en la dirección z debido a la gran diferencia entre el radio y la longitud de las membranas, con el propósito de corroborar el funcionamiento del modelo se compararon estadísticamente (Statgraphics Centurion XVI) los resultados generados por el software contra datos experimentales de evaporación osmótica en los que se usó agua destilada como alimentación, evaluando por triplicado las condiciones de proceso que se presentan en la tabla 2.
Propiedades físicas de los fluidos
elaboración propia.
Condiciones de proceso usadas en las simulaciones y validación experimental.
elaboración propia.
El sistema donde se evaporo osmóticamente el agua, se presenta en la figura 3, fue utilizado un módulo de fibras huecas, que se describe en la tabla 3, se usaron bombas peristálticas con una velocidad de flujo variable para la solución salmuera y el agua, la salmuera fue mantenida constante a concentración del 44% p/v de cloruro de calcio (CaCl2) por medio de saturación con cloruro nuevo y actividad de agua 0,33 (Aqualab 3TE–Decagon), la relación alimentación/salmuera fue de 1:7 (l), la temperatura del agua se ajustó por medio de un baño termostático, el flux es la pérdida de agua medida cada 10 minutos.
Características del módulo empleado en la experimentación.
elaboración propia.
Figura 3.
Diagrama del sistema de evaporación osmótica.
Elaboración propia
Se estableció cuál debería ser el modelo de difusión a ser utilizado en las posteriores simulaciones, para lo cual se desarrolló un ensayo (triplicado) bajo las siguientes condiciones: flujo “agua” 1 l/min, temperatura 25 °C, []sal: 44 %, Flujo salmuera 0,2 – 1,6 l/min, el tamaño de malla fue el sugerido por Comsol® y se estableció en “normal”. Al efectuar el ANOVA para comparar los 3 modelos utilizados y los datos experimentales se encontraron diferencias significativas (p<0,05), a un nivel del 95 % de confianza, al realizar la prueba de rangos múltiples (LSD) se encontró que las predicciones generadas por Comsol usando el mecanismo difusional tipo Knudsen, no presentaron diferencias significativas. En la figura 4 se presentan los flux predichos con los modelos contra los valores experimentales, el mejor ajuste presentado por el modelo Knudsen, esta explicado por el pequeño tamaño de poro que tiene la membrana (0,51 um), lo cual produce Kn ³ 10, en donde las moléculas tienden a colisionar frecuentemente con las paredes del poro siendo el modelo Knudsen más conveniente (Liu & Wei, 2014).
Figura 4
Comparación de los flux predichos por Comsol® usando diferentes modelos de difusividad.
elaboración propia
Las membranas utilizadas en el proceso de EO no hacen ninguna contribución a la selectividad, su única función es la estabilización de la interfaz entre los dos medios de contacto. Además, la propia membrana está presente como una barrera al flujo de vapor, el efecto de las características de la membrana sobre la transferencia de masa y calor ha sido estudiado por varios investigadores (Domingues, Ramos, Cardoso, & Reis, 2014; Viguié, Savart, Duru, Rouch, & Remigy, 2013; Yang, Yu, Yan, Tu, & Dahlquist, 2013), para tamaños de poro entre 0,05 y 0,2 um, el flux no muestra ser dependiente de este parámetro, en membranas con poros entre 0,2 y 0,45, se informa de ligeros aumentos en el flujo de vapor, debido principalmente a la dominancia del mecanismos difusional Knudsen. Para membranas con poros mayores a 1um, se ha observado un incremento significativo en el flux de evaporación a medida que aumenta el tamaño de poro, explicado esto por la marcada influencia de la difusión molecular (Naim, Khulbe, Ismail, & Matsuura, 2013; Rahbari-Sisakht, Ismail, Rana, & Matsuura, 2013).
El número de elementos finitos fue el segundo criterio analizado para comprobar los cálculos, la figura 5, presenta la comparación de los flux predichos con el modelo contra los valores experimentales del flux obtenidos para unas condiciones operativas de flujo “agua” 1 l/min, temperatura 25 °C, []sal: 44 %, Flujo salmuera 0,2 – 1,6 l/min y difusión Knudsen, se evidencio que usar un número reducido de elementos, produce resultados de flux bajos, a medida que se refina la malla los resultados se acercan cada vez más a los valores reales.
Figura 5
Efecto del tamaño de malla y número de elementos finitos usados para el cálculo del flux
elaboración propia.
Tal como lo evidenciaron los resultados del análisis estadístico ANOVA y la comparación de medias, la mejor correlación se obtuvo al usar una densidad de malla “normal”, esta no presento diferencias significativas frente a los datos experimentales, esta malla consistía de 4244 elementos finitos, con lo cual la diferencia entre el flux calculado y el experimental vario entre el 2 y 4 %, con este nivel de mallado el tiempo de cálculo fue de 5 minutos, seguir aumentando el número de elementos más allá de este límite, no mostro evidencia de mejorar significativamente la calidad de los resultados, por el contrario ocurrió una ligera sobre estimación de los mismos y se aumentaron los tiempos de convergencia
Se evaluaron estadísticamente cada curva de datos del flux experimental (flujo “agua” 0,25 – 1,5 l/min, temperatura 25 °C, []sal: 44 %, Flujo salmuera 0,2 – 1,6 l/min y difusión Knudsen) vs calculados, por medio de una prueba de t para comparar las medias, a un nivel del 95 % de confianza (p≤0,05), teniendo como hipótesis nula: mediaexp = mediacal y la hipótesis alternativa: mediaexp <> mediacal.
Figura 6.
Comparación del flux experimental y los calculados por Comsol® usando diferentes velocidades de alimentación.
elaboración propia.
Los resultados mostrarón que no existierón diferencias significativas entre los datos experimentales y los calculados por Comsol®, para ninguna de las 4 velocidades de alimentación, esto confirma la bondad de ajuste del modelo utilizado. Como se puede apreciar en la figura 6, el flux de evaporación aumenta proporcionalmente con la velocidad de la alimentación, este incremento no es indefinido estando limitado por las características estructurales y operativas de la membrana (Domingues et al., 2014).
Para el efecto de la temperatura sobre el flux de evaporación, dado en la figura 7, bajo condiciones operativas de: flujo “agua” 1 l/min, temperatura 25 - 40 °C, []sal: 44 %, flujo salmuera 0,2 – 1,6 l/min y difusión Knudsen; es claro que el flux es exponencialmente dependiente de la temperatura, este tipo de relación era esperada pues la presión de vapor del agua pura presenta este mismo comportamiento (Meng, Ye, Mansouri, & Chen, 2014; Zhao, Feron, Xie, Zhang, & Hoang, 2014).
Figura 7.
Comparación del flux experimental y los calculados por Comsol® usando diferentes temperaturas de alimentación.
elaboración propia.
Matemáticamente, lo anterior puede ser descrito para cualquier par de concentraciones de agua (j) y salmuera (s), con actividades de agua awj y aws, se tiene que el diferencial de presiones de vapor (Pv) estaría dado por la ecuación (49).
Donde A, B y C son constantes de la ecuación de Antoine, y T es la temperatura de la solución, esta expresión claramente muestra la relación exponencial que se mencionó anteriormente, entre la temperatura y la fuerza impulsora del proceso (Kujawski et al., 2013). Al comparar estadísticamente a un nivel del 95% de confianza , los resultados experimentales con los predichos para las 4 temperaturas evaluadas, no se encontraron diferencias significativas, los coeficientes de correlación fueron superiores al 97% en todos los casos, este buen nivel de ajuste por parte del modelo, puede ser explicado por el hecho de haber utilizado ecuaciones que describen con alto detalle, los fenómenos de transferencia y especialmente las propiedades físicas del agua y la salmuera.
Este parámetro se mostró consistentemente a través de todas las verificaciones realizadas, incrementar la velocidad de la salmuera, provoca un mayor flux de evaporación, la figura 8 presenta las tres zonas identificadas en el comportamiento del flux, en la zona I este crecimiento es el más rápido, para velocidades entre los 0,1 y 0,4 l/min con un incremento de 15% en el flux, la segunda zona comprendida desde 0,4 hasta 1,0 l/min presenta una fase de estabilización llegándose a un máximo de eficiencia cuando el flujo alcanza los 1,0 l/min, pasado este límite se establece un estancamiento o meseta (zona III), donde los incrementos en el flux son muy bajos (<0,9%).
Figura 8.
Efecto del flujo de salmuera sobre el comportamiento del flux de evaporación.
elaboración propia.
La principal explicación de este crecimiento progresivo del flux, es el aumento en el coeficiente de transferencia de materia en el lado de la carcasa (Manawi, Khraisheh, Fard, Benyahia, & Adham, 2014), por el aumento en la turbulencia y también el posible efecto de limpieza sobre la membrana que se ejerce al aumentar la velocidad, retirando depósitos del cloruro de calcio, que pudieran estar acumulados debido a la intrincada geometría formada por el empaquetamiento de las fibras en el módulo, este tipo de colmatación de la membrana no se evidencio durante los experimentos realizados.
El modelo difusional Knudsen fue el más adecuado en describir la transferencia del vapor de agua a través de la membrana, explicado esto por las características estructurales de la membrana que hacen predominar este modelo sobre la difusión molecular, siempre ha de verificarse de forma experimental la selección del modelo difusivo, porque este puede cambiar rápidamente de un tipo a otro con pequeños cambios en el proceso o la membrana. El uso de una geometría axial en dos dimensiones como simplificación del módulo real de evaporación osmótica es una decisión adecuada que permite reducir de forma significativa los recursos de cómputo necesarios para resolver las ecuaciones diferenciales sobre los tres dominios establecidos (agua, membrana, salmuera). Los modelos utilizados y las simulaciones generadas por Comsol® y Matlab® fueron adecuadas fenomenológicamente y estadísticamente significativas (p<0,05) en predecir el comportamiento del flux obteniéndose correlaciones mayores al 96% entre los datos experimentales y los calculados.
Colciencias. Departamento Administrativo de Ciencia, Tecnología e Innovación.
A Área (m-2)
a Actividad
C Concentración molar soluto (mol l-1)
D Coeficiente de difusión (m-2 s-1)
dh Diámetro hidráulico (m)
d Diámetro (m)
dp Diámetro de poro (m)
K Coeficiente transf masa (kg m-2 h-1 Pa-1)
kb Constante Boltzmann (1.3807x10-23 J K-1)
L Longitud (m)
M Peso molecular (kg mol-1)
m Masa (kg)
N Flux vapor, masa (kg m-2 h-1), (mol m-2 s-1).
P Presión (Pa)
P* Presión vapor saturado (Pa)
PAlm Presión aire media logarítmica
Q Flux de calor (W m-2)
R Constante de los gases (8.314J K-1 mol-1)
r Radio de poro (m)
T Temperatura (°C, K)
Fracción masa (p/p %)
Símbolos
Porosidad
Espesor (m)
Diferencia
Coeficiente de actividad
Trayectoria media libre (m)
Viscosidad dinámica (Pa s)
Tortuosidad
Valor asintótico
Densidad empaquetamiento
Densidad (kg m-3)
Diámetro medio de colisión
Números adimensionales
Gz Graetz
Kn Knudsen
Nu Nusselt
Pr Prandtl
Re Reynolds
Sc Schmidt
Sh Sherwood
Figura 1. a)
Principio del proceso
elaboración propia.
Figura 1. b)
Perfil de concentraciones y resistencia a la transferencia de materia en la EO.
elaboración propia.
Figura 2.
Geometría general del módulo de fibras huecas (Izq), dominios y límites para la simetría de la fibra (Der).
elaboración propia.
Propiedades físicas de los fluidos
elaboración propia.
Condiciones de proceso usadas en las simulaciones y validación experimental.
elaboración propia.
Características del módulo empleado en la experimentación.
elaboración propia.
Figura 3.
Diagrama del sistema de evaporación osmótica.
Elaboración propia
Figura 4
Comparación de los flux predichos por Comsol® usando diferentes modelos de difusividad.
elaboración propia
Figura 5
Efecto del tamaño de malla y número de elementos finitos usados para el cálculo del flux
elaboración propia.
Figura 6.
Comparación del flux experimental y los calculados por Comsol® usando diferentes velocidades de alimentación.
elaboración propia.
Figura 7.
Comparación del flux experimental y los calculados por Comsol® usando diferentes temperaturas de alimentación.
elaboración propia.
Figura 8.
Efecto del flujo de salmuera sobre el comportamiento del flux de evaporación.
elaboración propia.
