Artículos de Investigación Científica y Tecnológica

Límites termodinámicos a la productividad de los destiladores solares pasivos

Thermodinamical productivity limitations for passive solar stills

Limites termodinâmicos á produtividade de equipamentos de destilação solar pasivi

Henry Alberto Salinas-Freire ^{1 *}

Osney Pérez-Ones ²

Susana Rodríguez-Muñoz ²

¹ Universidad Técnica de Ambato, Campus Huachi, Av. Los Chasquis y Río Payamino, Ambato, Ecuador

² Universidad Tecnológica de La Habana José Antonio Echeverría, Ave. 114 N° 11901 e/119 y 127. Marianao, Ciudad Habana, Cuba.

* ha.salinas@uta.edu.ec

Cita: Salinas Freire HA, Pérez Ones O, Rodríguez Muñoz S. Límites termodinámicos a la productividad de los destiladores solares pasivos. rev.ion. 2019;32(1):7-20.

Creative Commons: https://creativecommons.org/licenses/by/4.0/

https://doi.org/10.18273/revion.v32n1-2019001

Fecha recepción: 26 de febrero de 2018

Fecha aceptación: 1 de febrero de 2019

Resumen

La desalinización de agua de mar se ha estudiado con interés debido a la escasez de agua dulce para consumo humano. La destilación solar es un método antiguo; la productividad, consumo energético del proceso y el costo del agua desalinizada así obtenida dependen de la eficiencia que se alcance en cada una de las etapas de estos sistemas. La limitada capacidad de absorber la radiación solar y transformarla en calor útil para la evaporación, la interacción con el medio circundante, y las pérdidas de calor limitan la eficiencia global del proceso térmico y la productividad. Como la energía proviene de la radiación solar, la productividad máxima de este proceso estará limitada por la magnitud de la radiación solar total disponible en una zona del planeta debido a su ubicación geográfica, época del año y condiciones climáticas locales. Los procesos a partir de esa energía estarán limitados termodinámicamente por los coeficientes de transferencia térmica que se alcancen en el equipo, por el máximo que pueda alcanzar el calor de evaporación, mientras que las pérdidas al ambiente por convección y radiación sean mínimas. Se presentan análisis comparativos de varios modelos propuestos, datos reportados de destiladores, datos reportados de radiación solar que alcanzan valores promedio de hasta 7,2 – 7,4 kWh/m² en algunas regiones del planeta y se hacen estimaciones para productividad de estos equipos que alcanzan entre 6,7 y 6,9 kg/m²día y la eficiencia máxima teórica se encuentra alrededor de 0,16 de la radiación solar total.

Palabras clave: Desalinización; Destilador Solar Pasivo; Productividad Destilador Solar; Eficiencia Destilador Solar.

Abstract

Seawater desalination is a process that has been widely studied with great interest due the growing water scarcity for human use. Solar distillation is a very ancient method, the productivity, energy requirements and the consequent cost of the desalinated water depends at least of the efficiency that can be reached in every single step of those systems. Due the limited solar irradiance absorption and to water evaporation transformation capacity, due to glass and water radiation absorbance and the ambient interaction, heat losses limit the overall thermal process efficiency and consequently productivity. Since in passive solar stills all energy comes from solar irradiation, this process maximum productivity will be limited for total solar irradiance that may be available in a defined zone of the world defined by their geographical localization, season of the year and local climatic conditions. All thermal process with this energy will be thermodynamically limited by heat transfer coefficient values reached in the dispositive, by the maximum value that can reach the evaporative heat transfer coefficient, and heat radiation and convective losses coefficients being minimum. In this paper, comparative analysis is presented from several process proposed models, device reported data; solar irradiation reported data reaching about 7.2 – 7.4 kWh/m² and estimations for productivity intervals for this devices area about 6.7 – 6.9 kg/m²d; efficiency calculated with available reported data are about 0.16 of the global solar irradiance.

Keywords: Desalination; Passive Solar Still; Solar Still Productivity; Solar Still Efficiency.

Resumo

A dessalinização da água do mar tem sido estudada com interesse por causa da escassez de água doce para consumo humano. Destilação solar é um método antigo, a produtividade, consumo de energia do processo e o custo da água dessalinizada assim obtido dependem da eficiência em cada fase de estes sistemas. A capacidade limitada de absorver radiação solar e transformando-a em calor útil para a evaporação, a interação com o meio circundante e as perdas de calor, limita a eficiência global do processo térmico e produtividade. Como a energia vem de radiação solar, a máxima produtividade deste processo é limitado pelo tamanho da radiação solar total disponível em uma área do planeta devido à sua localização geográfica, a época do ano e das condições climáticas locais. Os processos a partir dessa energia estão limitadados termodinamicamente pelo coeficientes de transferência térmico que são alcançados no equipamento, pôr o valor máximo que pode atingir o calor de evaporação, enquanto que as perdas de calor para o ambiente por convecção e radiação são mínimas. São apresentados, análise comparativa dos diferentes modelos propostos, dados reportados de destiladores, dados de irradiação solar relatados que atingem valores médios até 7,2 - 7,4 kWh/m² em algumas regiões do mundo. São feitas estimações de coeficientes de transferência em estas etapas e os intervalos de produtividade para estes dispositivos que atingem até 6,7 - 6,9 kg/m²d e eficiência é entre 0,16 da radiação solar total disponível.

Palavras-chave: Dessalinização; Destilador Solar Passiva; a Produtividade Destilador Solar; a Eficiência Destilador Solar.

Introducción

La desalinización del agua de mar se conoce desde la antigüedad ^{(1),( 2)}, se ha reseñado ^{(1),( 3),( 4)} que Tales de Mileto y Aristóteles lo mencionan en sus escritos y describen cómo se habrían utilizado dispositivos primitivos basados en hervir el agua en vasijas de barro o calderos y recogerla en esponjas en los barcos para conseguir agua dulce durante viajes largos. Desde entonces se han estudiado y desarrollado varios métodos y sistemas de desalinización que han pasado de los prototipos a instalaciones de gran capacidad y esto ha permitido conocer costos y eficiencias de las etapas en cada proceso ^{(5),(6),(7), (8),( 9)}.

Entre los sistemas de desalinización empleados actualmente, se estudia con mucho interés aquellos basados en energías renovables por sobre aquellos que emplean energía de combustibles fósiles ^{(10), (11), (12), (13), (14), (15), (16)}. La desalinización solar o SD consiste en colocar agua de mar en un recipiente donde a través de una película transparente, generalmente vidrio, se emplea la radiación solar para calentar la solución salina y evaporar una parte del agua que pasa al aire de la cámara y posteriormente condensa al entrar en contacto con una superficie más fría desde donde se colecta ^(17),(18). Estos equipos no emplean otra fuente de energía y por lo tanto su rendimiento está en función de la radiación solar que se pueda captar, y de la eficiencia de los procesos al interior del equipo.

Por la forma de operación de los desalinizadores solares pasivos, todas las etapas de transferencia de masa y energía tienen límites, los cuales se analizan en este trabajo y se comparan con datos reportados, con el fin de determinar la productividad máxima teórica esperada.

Marco teórico

Los desalinizadores solares pasivos reciben toda la energía para su operación a partir de la radiación solar incidente, ésta se transfiere desde la cubierta de vidrio hacia el fondo del destilador y desde ahí, transformada en energía térmica, hacia el agua, las paredes, el fondo y la cubierta del equipo. Una parte de esa energía evapora agua y la restante se pierde al ambiente. Para su análisis se emplean varios modelos de transferencia de calor y masa ^(19),(20), así como correlaciones para las propiedades termodinámicas y de transporte del agua y del vapor ⁽²¹⁾. Los modelos matemáticos se usan para analizar el comportamiento del desalinizador en las mejores condiciones de operación para la temperatura del agua, del vidrio y la radiación  solar y se comparan los resultados. Para el desarrollo del presente trabajo se considera un destilador solar pasivo de una sola cubierta representado mediante el esquema de la Figura 1.

Figura 1. Esquema del destilador solar pasivo.

En el destilador, se considera que no hay pérdidas de masa por fugas del equipo, que no hay gradientes de temperatura en el vidrio, el agua o el material aislante y que las variaciones debidas a la temperatura en la capacidad calorífica de la base del destilador del vidrio y del aislante son despreciables ^{(21), (22)}.

Radiación solar

La radiación solar que alcanza la superficie de la tierra no es monocromática ⁽²³⁾, abarca una gama de longitudes de onda que van del infrarrojo hasta el ultravioleta, pasando por el espectro visible. En los desalinizadores solares esta radiación I_h llega hasta la superficie del vidrio donde una fracción α_b es absorbida, una fracción τ_v es transmitida y una fracción R_v es reflejada. La fracción transmitida atraviesa el aire húmedo con pérdidas despreciables y alcanza el agua en el fondo del destilador, donde nuevamente la radiación total se divide entre lo absorbido α_w, reflejado R_w y transmitido τ_w. Un esquema de este proceso para el vidrio se muestra en la Figura 2, la fracción transmitida tiene el mismo compartimiento en el aire, aunque con pérdidas despreciables, al igual que en el agua del fondo del destilador.

Figura 2. Esquema del proceso de incidencia de la radiación solar en el vidrio.

La fracción transmitida por el agua, llega hasta el fondo del destilador, generalmente pintado de color negro, donde la mayor parte se absorbe, una pequeña fracción se refleja y prácticamente nada se transmite; allí calienta la superficie del fondo del destilador, y este calor empieza a transmitirse hacia el agua y hacia el material aislante en el fondo y en las paredes del destilador. El agua es calentada y del total de la energía térmica generada que se transfiere al vidrio, una parte lo hace por convección natural, otra parte por radiación y otra parte evapora el agua de la solución salina y se transfiere en forma de calor de evaporación; esta es la única fracción de todo calor involucrado en el proceso que sirve para obtener destilado, se transfiere al vidrio donde se condensa el vapor por la diferencia de temperatura, y de allí se disipa al ambiente ⁽²⁴⁾. Este proceso se muestra en la Figura 3.

Figura 3. Transferencia de calor desde el agua caliente en el fondo del destilador.

La radiación solar varía de un lugar a otro dependiendo de factores como la ubicación geográfica, la época del año, la nubosidad y la hora del día. Los valores se pueden obtener por medición directa o a través de publicaciones de instituciones como el Banco Mundial. Es posible encontrar datos útiles para la mayoría de regiones ⁽²⁵⁾.

En el proceso de absorción de la radiación solar, Cooper ⁽¹⁸⁾ y otros autores ⁽²⁶⁾ han determinado que solamente entre el 85 y el 95% de la radiación I_h incidente alcanza el fondo del destilador.

Balance de energía

De acuerdo con el diagrama de la Figura 1 se plantean los balances de energía así ⁽²²⁾:

Para el recipiente que contiene el agua en el fondo del destilador:

Donde:

I_h  = radiación incidente en el vidrio  [W/m²],

τ_v = transmisividad del vidrio, adimensional

τ_w = transmisividad del agua, adimensional

α_b = absortividad del recipiente o fondo del destilador, adimensional

A_b = área de la base del destilador [m²],

A_f = área del fondo del destilador [m²],

h_c,b-w = coeficiente de transferencia de calor por convección del fondo del destilador al agua [W/ m²°C],

T_p = temperatura de la base del destilador [°C],

T_f = temperatura del fondo del destilador [°C],

T_w = temperatura del agua [°C],

U_f = coeficiente total de transferencia de calor del fondo del destilador al ambiente [W/m²°C], y

T_a = temperatura ambiente [°C].

Balance de Energía para el agua en el recipiente del destilador:

Donde:

A  = área del agua del destilador [m]

m_w = masa del agua en el destilador [kg],

c_w = capacidad calorífica del agua [J/kg°C], dT_w/dt = variación de la temperatura del agua respecto del tiempo [°C/s],

h_t,w-v = coeficiente total de transferencia de calor desde el agua al vidrio [W/m²°C], y

T_vi = temperatura de la cara interna del vidrio [°C].

El coeficiente total de transferencia de calor del agua al vidrio, definido por Sampathkumar ⁽²⁷⁾ es

Donde:

h_c,w-v = coeficiente de transferencia de calor por convección del agua al vidrio [W/m²°C],

h_r,w-v = coeficiente de transferencia de calor por radiación del agua al vidrio [W/m²°C],

h_e,w-v = coeficiente de transferencia de calor por evaporación del agua al vidrio [W/m²°C].

Balance de Energía para el vidrio de la cubierta

Donde:

α_v = absorbancia del vidrio, adimensional A_v = área del vidrio [m²], h_r,v-c = coeficiente de transferencia de calor por radiación del vidrio al cielo [W/m²°C],

T_vo = temperatura de la cara externa del vidrio [°C], T_c = temperatura del cielo [°C], h_c,v-a = coeficiente de transferencia de calor por convección del vidrio al ambiente [W/m²°C].

Correlaciones para las propiedades del agua y del aire

Para evaluar los coeficientes es necesario conocer las propiedades del agua y el aire en las condiciones de operación. Estas propiedades como densidad, viscosidad y presión de saturación de vapor, varían en función de la temperatura, por lo que se emplean correlaciones para su cálculo.

Presión de saturación o vapor del agua. Existen varias ecuaciones que relacionan la presión de vapor del agua con la temperatura, de las cuales se evalúan tres. La correlación de Dunkle ⁽²⁸⁾, válida para temperaturas inferiores a 70°C, permite calcular la presión de vapor P del agua en Pa a la temperatura T (en °C) y se expresa en la Ecuación 5.

La ecuación de Antoine ⁽²⁹⁾ también correlaciona presiones de vapor con temperatura para substancias en un intervalo mayor. En el caso del agua, la ecuación que es válida entre 0 y 200°C es:

Sharma y Mullic (21) presentan una ecuación de Keenan y Keyes (1936) para calcular la presión de vapor del agua, la cual es válida entre 10 y 150°C, mediante:

Donde:

,

 X = 647,27 – T,

 T = temperatura [°C],

a = 3,2437814,

b = 5,86826x10^-3, c = 1,1702379x10^-8, y

d = 2,1878462 x10^-3.

Transferencia de calor interna

Comprende la energía transferida desde la superficie del agua hacia la superficie interna del vidrio que ocurre principalmente por radiación, convección y evaporación ⁽²⁷⁾.

Calor transferido por convección. Se define el calor de convección q_c,w-v en W/m² transferido desde la superficie del agua hacia el vidrio (30) por la Ecuación 8.

Cooper ⁽¹⁸⁾, Sharma y Mullick ⁽²¹⁾ y Sampathkumar ⁽²⁷⁾ establecen el flujo de calor desde la superficie del agua hacia el vidrio usando las ecuaciones de Dunkle:

Según Sampathkumar ⁽²⁷⁾, el coeficiente de transferencia de calor de convección se puede expresar por:

Donde

P_w = presión de vapor del agua a la temperatura

T_w, [Pa], y

P_vi = presión de vapor del agua a la temperatura de la pared interna del vidrio, [Pa].

Estas correlaciones son válidas para temperaturas de operación alrededor de 50°C y un valor de ΔT’ alrededor de 17°C, además son independientes del volumen de la cámara y válidas solamente para el flujo de calor hacia arriba en el espacio cerrado entre las superficies de evaporación y condensación.

Calor transferido por radiación. La Ecuación 12 define el calor transferido por radiación entre la superficie del agua y la superficie interna del vidrio ⁽³⁰⁾:

Donde:

q_r,w-v = calor transferido por radiación, [W/m²],

El coeficiente de transferencia de calor por radiación se define en la Ecuación 13 empleando la constante de Stefan-Woltzman σ:

Donde:

ԑ_ef = emisividad efectiva, adimensional, y σ = 5,67x10^-8 [W/m²K⁴].

Cooper ⁽¹⁸⁾ y Sharma-Mullick ⁽²¹⁾ establecen como un valor estimado de 0,9; reemplazándolo en la ecuación 13 y ésta en la 12 se obtiene:

Para obtener datos más exactos o casos particulares, la emisividad efectiva se calcula a partir de las emisividades del vidrio ԑ_v y del agua ԑ_w, mediante la Ecuación 15 (27).

Calor transferido por evaporación. El calor q_e,w-v en W/m² transferido por evaporación desde el agua hacia el vidrio, como función del calor de convección, según Cooper ⁽¹⁸⁾ es:

Los valores de los coeficientes están relacionados de acuerdo con Setoodeh (31) mediante la Ecuación 17, donde se asume que el calor de vaporización del agua es constante.

El coeficiente de transferencia de calor por evaporación según Sharma y Mullick (21) es definido por la Ecuación 18. Esta ecuación incluye las variaciones del calor de vaporización del agua con la temperatura.

Coeficiente total de transferencia entre el agua y el vidrio. De acuerdo con el balance expresado en la Ecuación 2, el coeficiente total de transferencia se puede calcular por la Ecuación 3 ⁽²⁷⁾.

Transferencia de calor externa. El calor se pierde en el destilador solar desde la superficie del vidrio hacia el ambiente por convección natural y radiación, y por las paredes y el fondo del destilador solar, que tienen un material aislante, por radiación y convección (27). Se analizará únicamente el calor perdido por el vidrio al ambiente por convección natural, y por radiación; asumiendo que en comparación con aquellos, el calor perdido por las paredes y el fondo es despreciable.

Calor transferido por convección. Las pérdidas de calor por convección q_c,v-a en W/m² del destilador se pueden expresar en función de la velocidad del viento. Aboul-Einein ⁽²²⁾, Shukla ⁽³²⁾] y Sampathkumar ⁽²⁷⁾ emplean para el calor perdido y el coeficiente de transferencia entre el vidrio y el ambiente una correlación propuesta por Duffie and Beckman (1980):

Donde:

V = velocidad del viento [m/s].

El coeficiente de transferencia de calor por convección perdido por el destilador solar desde el vidrio hacia la atmósfera, o hacia el cielo, según Sharma y Mullic ⁽²¹⁾ se define en la Ecuación 21:

Donde

h_c,v-c = coeficiente de transferencia de calor por convección [W/m²], y

h_w = coeficiente de transferencia de calor debido al viento [W/m²].

Coeficiente de transferencia de calor por radiación entre el vidrio y el ambiente. El coeficiente de transferencia de calor de radiación transferido entre el vidrio y el ambiente definido por Mullick ⁽²¹⁾ es:

Donde:

ε_g= emisividad del vidrio, adimensional T_c = temperatura del cielo, [K].

Coeficiente de pérdidas totales de calor entre el vidrio y el ambiente. Este coeficiente comprende la suma de las pérdidas por convección y por radiación, y se expresa mediante la Ecuación 23:

Donde:

h_t,gs = coeficiente total de transferencia de calor entre el vidrio y el ambiente, [W/m²K], Reemplazando las definiciones de la ecuación 21 y 22 en la ecuación 23 se obtiene:

Algunos autores ^{(1, 27, 33)} emplean para el cálculo del coeficiente de pérdidas totales desde el vidrio al ambiente en función de la velocidad del viento una correlación que se muestra en la Ecuación 25.

Indicadores de rendimiento

Productividad empleando el calor de evaporación. La tasa de evaporación de agua instantánea  en kg/s, de acuerdo con Mowla ⁽³⁴⁾, está dada por:

La evaporación M_ew para un período de tiempo t es:

Si el período de tiempo es una hora, en segundos, se tiene ^{(20, 22, 32)}:

Para Setoodeh ⁽³¹⁾ y Sampathkumar ⁽²⁷⁾ la masa de agua evaporada m_ew en un período de tiempo t se puede calcular por la Ecuación 29.

Sampathkumar ⁽²⁷⁾ define la producción diaria, incluyendo las horas de la noche, como:

Productividad empleando los números adimensionales Pr y Gr. Este método ^(19),(27), emplea los números adimensionales de Prandtl y Grashoff, como se muestra en la Ecuación 31.

Donde: k = conductividad térmica del aire húmedo

[W/ m°C], y

d = longitud media característica entre las superficies de evaporación y de condensación [m], C = constante, adimensional Gr = número de Grashoff, adimensional y Pr = número de Prandtl, adimensional.

Los números adimensionales de Grashoff y Prandtl son:

Donde:

β = factor del coeficiente volumétrico de expansión

térmica del aire [K^-1], ρ = densidad del vapor [kg/m³], y μ = viscosidad del aire húmedo [Pa.s].

Setoodeh ⁽³¹⁾ también presenta una correlación para ciertas propiedades que se evalúan a T_v que es la temperatura promedio entre T_w y T_g:

Donde para T < 70°C el calor latente de vaporización del agua h_{f g} es:

Para la viscosidad:

La capacidad calorífica del aire:

Reemplazando las definiciones y despejando , se establece ⁽¹⁹⁾ que la productividad m_ew se puede  calcular empleando los números de Grashoff y Pandtl mediante la Ecuación 41:

Los valores de C y n se calculan, como propone de cada destilador solar y del intervalo de temperatura Tiwari ⁽³³⁾, mediante regresiones a partir de datos del agua en que operan, por lo que se deben experimentales; en su trabajo, Sampathkumar ⁽²⁷⁾ determinar en cada caso particular. En la Tabla 1 se establece que C y n dependen del diseño específico presentan algunos valores reportados para C y n.

Tabla 1. Valores de C y n reportados y su rango de validez

Productividad empleando la radiación solar incidente. P. I. Cooper ⁽¹⁸⁾ en su trabajo presenta una correlación lineal, válida para niveles de radiación I_h entre 0,4 y 1,39 kW/m² y para una temperatura ambiente de 30°C, para calcular la productividad instantánea la cual se muestra en la Ecuación 42.

Eficiencia

A partir del agua evaporada. Para calcular la eficiencia η_t de un destilador solar pasivo a partir del calor de vaporización la masa de agua condensada y la radiación incidente medida, se emplea la Ecuación 43 ⁽²⁰⁾:

Donde:

A_s = superficie de la base del destilador, [m], y ∆t = período de tiempo, [s].

Eficiencia máxima teórica a partir de la radiación solar. Cooper ⁽¹⁸⁾ presenta una ecuación para hallar la eficiencia máxima teórica diaria que se muestra en la Ecuación 45. Empleando esta ecuación se obtienen datos más altos que los hallados empleando la ecuación 44 para el mismo caso.

Donde:

η_{o =} eficiencia máxima teórica diaria, adimensional, θ_{s =} horas de luz de un día medidas desde la salida del sol hasta el ocaso, [h].

Para los cálculos se consideran los datos ⁽³⁵⁾ de un destilador solar con los valores presentados en la Tabla 2 ⁽²²⁾:

Donde x_a y x_v son el espesor del agua y del vidrio respectivamente y k_v la conductividad térmica del vidrio.

Tabla 2. Parámetros del destilador solar pasivo.

Resultados

Desviaciones de las correlaciones para estimar la presión de saturación de vapor del agua. La presión de vapor del agua entre 0 y 100°C calculada mediante las correlaciones de Antoine, Sharma y Dunkle difieren entre sí y con respecto a los datos de Lemmon, McLinden y Friend reportados en Perry ⁽³⁶⁾; su desviación se muestra en la Figura 4.

Figura 4. Desviación porcentual de la presión de vapor del agua calculada según las correlaciones de Dunkle, Antoine y Sharma-Mullick (Keenan-Keyes).

Se aprecia que la correlación de Dunkle es la que presenta mayores desviaciones, entre +2,7 y -3,1%, mientras que los menores valores son para la de Antoine con -0,4 y 0,7%.

Transferencia de calor interna

Calor de convección

En la Figura 5 se muestra el comportamiento del coeficiente de transferencia de calor por convección en función de la temperatura del agua T_w y la temperatura de la cara interna del vidrio T_vi. Se aprecia que se incrementa a medida que aumenta Tw y disminuye Tvi.

Figura 5. Coeficiente de transporte de calor por convección h  en función de las temperaturas delagua T c,w-v w  y del vidrio T vi

Calor transferido por radiación

El comportamiento del coeficiente de transferencia de calor por radiación con respecto a las temperaturas del agua T w  y del vidrio T  se aprecia en la Figura 6. El coeficiente puede alcanzar 9,0 W/m 2 v °C cuando la diferencia de temperatura es de 50°C, valores superiores, cuando ambas temperaturas son iguales o la temperatura del vidrio es mayor que la del agua, carecen de sentido físico. La magnitud del coeficiente aumenta a medida que se incrementa la temperatura del agua y disminuye la temperatura del vidrio.

Figura 6. Coeficiente de transporte de calor por radiación h_r,w-v en función de las temperaturas del agua y del vidrio.

Calor de evaporación

El comportamiento del coeficiente de transferencia de calor de evaporación calculado a partir de la ecuación 17 en función de las temperaturas del agua y del vidrio se muestra en la Figura 7. El valor que alcanza este coeficiente es mucho mayor en magnitud que el calor de convección y el calor de radiación, alrededor de 120 W/m²°C para diferencias de temperatura de 60°C, el segmento donde la temperatura del vidrio es mayor que la del agua, carece de sentido físico.

Figura 7. Coeficiente de transporte de calor por evaporación propuesto por Mowla y Tiwari ^h_e,w-v en función de las temperaturas del agua T_w y la cara interna del vidrio T_vi.

Los valores calculados con las Ecuaciones 17 y 18 varían muy poco. El comportamiento del coeficiente de transferencia de calor por evaporación entre el agua y la superficie interna del vidrio calculado mediante la Ecuación 18 se muestra en la Figura 8.

Figura 8. Coeficiente de transporte de calor de evaporación h_,e,w-g considerando el calor de vaporización variable en función de las temperaturas del agua y del vidrio.

Coeficiente total de transferencia entre el agua y el vidrio

Corresponde a la suma del calor transferido por radiación, convección y evaporación, se calcula con la ecuación 3 y se muestra en la Figura 9. Este coeficiente puede alcanzar hasta 135 W/m²°C para diferencias de temperatura de 50 °C.

Figura 9. Coeficiente total de transferencia de calor entre el agua y el vidrio h_t,w-v en función de las temperaturas del agua T_w y la cara interna del vidrio T_vi.

Indicadores de rendimiento

Productividad empleando el calor de evaporación

El comportamiento de la productividad diaria calculada a partir de los coeficientes de transferencia por convección y por evaporación, por la Ecuación 30, se observa en la Figura 10 en función de ^T_w y ^T_vi.

Figura 10. Productividad diaria del destilador calculada empleando los coeficientes de transferencia, en función de las temperaturas del agua T_w y la cara interna del vidrio T_vi.

Productividad empleando números adimensionales

En la Figura 11 se observa el comportamiento de la productividad en función de las temperaturas del agua y del vidrio, donde se han empleado 0,0112 y 0,4088 para C y n respectivamente. Se debe notar que en este caso, los valores de productividad son menores que aquellos mostrados en la Figura 10 en las mismas condiciones de temperatura calculadas empleando el calor de evaporación. Los valores máximos reportados en la práctica alcanzados para la productividad son alrededor de 0,550 kg/m²h ⁽³³⁾. Se aprecia que la productividad aumenta con el incremento de la temperatura del agua y la disminución de la temperatura del vidrio.

Figura 11. Productividad de destilado calculada empleando los números adimensionales, en función de las temperaturas del agua T_w y del vidrio T_vi.

Magnitudes de los coeficientes de transferencia

Empleando los modelos de transferencia descritos, se puede calcular un rango que pueden alcanzar los valores de los coeficientes de transferencia de calor, que se muestra en la Tabla 3.

Tabla 3. Valores que alcanzan los coeficientes de transferencia de calor en los destiladores solares pasivos

En la tabla 3 se aprecia que en condiciones de bajas temperaturas y con una baja velocidad del viento, el total de energía incidente solamente se puede aprovechar en un 37 % sin considerar las pérdidas por reflexión del material transparente que varían entre 10 y 15 % de la radiación incidente ⁽¹⁸⁾, y en condiciones de alta temperatura y alta velocidad de viento este calor evaporativo puede alcanzar hasta un 62 % del total. Esto es en condiciones ideales, considerando por ejemplo que todo el tiempo el cristal del evaporador está en posición perpendicular respecto a la radiación incidente, lo cual no se cumple en los evaporadores solares pasivos. La velocidad del viento afecta de dos maneras el proceso; enfriando la cara externa del vidrio del evaporador y ayudando en la condensación, y al mismo tiempo enfriando el resto del equipo, por lo que una alta velocidad incrementa las pérdidas por convección.

Productividades teóricas y experimentales

Se ha calculado la productividad máxima teórica a partir de la radiación solar mediante la ecuación 44 y se muestra junto con los datos obtenidos experimentalmente en varios ensayos en la Figura 12.

Las máximas productividades teóricas se ubican entre 0,42 y 0,75 kg/h mientras que las experimentales alcanzan máximos entre 0,15 y 0,30 kg/h. De acuerdo con el atlas publicado por el Banco Mundial ⁽²⁵⁾, hay zonas de alta radiación solar que alcanzan un promedio anual entre 7,2 y 7,4 kW.h/m², empleando este modelo ⁽¹⁸⁾, el valor de productividad de 6,7 a 6,9 kg/m²día representa un máximo teórico en condiciones ideales.

Se ha calculado la eficiencia máxima teórica empleando la Ecuación 43 del modelo de radiación solar y se presenta en la Figura 13 junto con la eficiencia obtenida experimentalmente en varios ensayos.

Figura 12. Comparación de productividades máxima teórica y experimental para cuatro casos. a) Shukla ^(32),(32), c) ⁽³²⁾ y d) Dev ⁽³⁷⁾.

Figura 13. Comparación de eficiencia máxima teórica y experimental para cuatro ensayos. a) Shukla ⁽³²⁾, b) ⁽³²⁾,⁽³²⁾ y d) Dev ⁽³⁷⁾ Se aprecia que la eficiencia máxima teórica alcanza valores entre 0,25 y 0,63 mientras que los datos experimentales se ubican entre 0,01 y 0,45.

Conclusiones

1. En un destilador solar pasivo siempre estarán presentes los mecanismos de transferencia por convección y radiación al interior, así como las pérdidas por reflexión y convección al exterior, por lo tanto el calor aprovechable máximo teórico para evaporar el agua oscilará entre 33 y 56 % de la radiación solar total en condiciones ideales.

2. El rendimiento del destilador solar calculado mediante el modelo de los coeficientes de transferencia de calor, y el calculado empleando las correlaciones de los números adimensionales proporcionan datos que representan con precisión el fenómeno, mientras que las correlaciones basadas en la radiación solar presentan los máximos teóricos para productividad y eficiencia, que son coherentemente superiores a los medidos en todos los casos.

3. La productividad de los destiladores solares pasivos en función de la irradiación solar disponible, tiene límites termodinámicos definidos para su productividad entre 6,7 y 6,9 kg/m²día, los cuales están determinados por la intensidad de la radiación solar incidente, los materiales de construcción, las condiciones locales de temperatura y la velocidad del viento.

Referencias bibliográficas

1. Dwivedi V, Tiwari G. Experimental validation of thermal model of a double slope active solar still under natural circulation mode. Desalination. 2010;250(1):49-55

2. Hernández H, Rubalcaba E, Hermosillo J-J. Improvement of a MEH desalination unitby means of heat recovery. Energy Procedia. 2014;57:2781-6

3. Moya EZ, editor Desalinización del agua del mar mediante energías renovables. Actas del I y II seminario del agua; 1997: Instituto de Estudios Almerienses

4. Tiwari G.N. SHN, Tripathi R. Present status of solar distillation. Solar energy. 2003;75:7

5. Moore BA, Martinson E, Raviv D. Waste to water: a low energy water distillation method. Desalination. 2008;220(1-3):502-5

6. Youssef P, Al-Dadah R, Mahmoud S. Comparative analysis of desalination technologies. Energy Procedia. 2014;61:2604- 7

7. Al-Weshahi MA, Tian G, Anderson A. Performance enhancement of MSF desalination by recovering stage heat from distillate water using internal heat exchanger. Energy Procedia. 2014;61:381-4

8. Darwish M. Qatar water problem and solar desalination. Desalination and Water Treatment. 2014;52(7-9):1250-62

9. Khayet M. Solar desalination by membrane distillation: Dispersion in energy consumption analysis and water production costs (a review). Desalination. 2013;308:89-101

10. Ma Q, Yi C, Lu H, Xie L, Fan J, He W. A conceptual demonstration and theoretical design of a novel “super-gravity” vacuum flash process for seawater desalination. Desalination. 2015;371:67-77

11. El-Sebaii A, El-Bialy E. Advanced designs of solar desalination systems: A review. Renewable and Sustainable Energy Reviews. 2015;49:1198-212

12. Asiedu-Boateng P, Nyarko B, Yamoah S, Ameyaw F, Tuffour-Acheampong K. Comparison of the Cost of Co-Production of Power and Desalinated Water from Different Power Cycles. Energy and Power Engineering. 2013;5(01):26.

13. Hamed OA, Kosaka H, Bamardouf KH, AlShail K, Al-Ghamdi AS. Concentrating solar power for seawater thermal desalination. Desalination. 2016;396:70-8

14. Gabriel KJ, Linke P, El-Halwagi MM. Optimization of multi-effect distillation process using a linear enthalpy model. Desalination. 2015;365:261-76

15. Compain P. Solar energy for water desalination. Procedia Engineering. 2012;46:220-7

16. Abdelmoez W, Mahmoud MS, Farrag TE. Water desalination using humidification/ dehumidification (HDH) technique powered by solar energy: a detailed review. Desalination and Water Treatment. 2014;52(25-27):4622- 40

17. Abdallah SB, Frikha N, Gabsi S. Study of the performances of different configurations of seawater desalination with a solar membrane distillation. Desalination and Water Treatment. 2014;52(13-15):2362-71.

18. Cooper P. The maximum efficiency of single-effect solar stills. Solar Energy. 1973;15(3):205IN1215-214217

19. Kumar S, Tiwari G. Estimation of convective mass transfer in solar distillation systems. Solar energy. 1996;57(6):459-64

20. Goosen MF, Sablani SS, Shayya WH, Paton C, Al-Hinai H. Thermodynamic and economic considerations in solar desalination. Desalination. 2000;129(1):63-89

21. Mullick S. Estimation of heat-transfer coefficients, the upward heat flow, and evaporation in a solar still. Journal of solar energy engineering. 1991;113:37

22. Aboul-Enein S, El-Sebaii A, El-Bialy E. Investigation of a single-basin solar still with deep basins. Renewable Energy. 1998;14(14):299-305

23. Ho-ming Yeh N-tM. Energy balances for upward-type,double-effect solar stills. Energy. 1990;15(12):9

24. Imad Al-Hayek OB. The effect of using different designs of solar stills on water distillation. Desalination. 2004;169:7

25. ESMAP, cartographer World GHI Global horizontal Irradiance Poster map Map]. www. worldbankgroup.com: Solargis; 2017

26. Yeh H-M, Ma N-T. Energy balances for upward-type, double-effect solar stills. Energy. 1990;15(12):1161-9

27. Sampathkumar K, Arjunan T, Pitchandi P, Senthilkumar P. Active solar distillation—a detailed review. Renewable and Sustainable Energy Reviews. 2010;14(6):1503-26

28. Dunkle R. Solar water distillation: the roof type still and the multiple effect diffusor. Int Dev in Heat Transfer;(). 1961

29. Smith JM, Van Ness HC, Abbott MM, García CR. Introducción a la termodinámica en Ingeniería Química. 7ma ed. Mexico: Interamericana M-H; 2007

30. Cengel Y. Transferencia de calor y masa. 3a ed. Mexico: McGraw Hill; 2007.

31. Setoodeh N, Rahimi R, Ameri A. Modeling and determination of heat transfer coefficient in a basin solar still using CFD. Desalination. 2011;268(1):103-10

32. Shukla S, Sorayan V. Thermal modeling of solar stills: an experimental validation. Renewable Energy. 2005;30(5):683-99

33. Tiwari G, Shukla S, Singh I. Computer modeling of passive/active solar stills by using inner glass temperature. Desalination. 2003;154(2):171-85

34. Mowla D, Karimi G. Mathematical modelling of solar stills in Iran. Solar Energy. 1995;55(5):389- 93

35. Dev R, Tiwari GN. Characteristic equation of a passive solar still. Desalination. 2009;245(13):246-65

36. Perry R. GD. Perry’s Chemical Engineers Handbook. 8 ed. Hill M, editor. New York: McGraw Hill; 2008 2008. 2735 p

37. Dev R, Tiwari G. Characteristic equation of a passive solar still. Desalination. 2009;245(13):246-65.