Suscríbete a DeepL Pro para poder traducir archivos de mayor tamaño.

Más información disponible en www.DeepL.com/pro.

Materiales hoy: Proceedings xxx (xxxx) xxx

Listas de contenidos disponibles en ScienceDirect

Materiales hoy: Actas

Pág i n a w e b d e l a r e v i s ta :
w w w .e l s e vi e r.co m /l o c a te /m a tp r

Modelización CFD no permanente de un destilador solar de pendiente única

Gaurav Mittal
Departamento de Ingeniería Mecánica, Graphic Era Deemed to be University, Dehradun, Uttarakhand, India

A R T Í C U L OEN F OR A B S T R A C T O
MA
Se ha realizado un modelo CFD no permanente de una destilación solar de pendiente simple utilizando
Historia del artículo: ANSYS Fluent. En el modelo, sólo se considera la zona húmeda (mezcla de aire + vapor de agua) con
Recibido el 6 de enero de 2021 temperaturas especificadas en el nivel inferior (agua) y superior (vidrio). Las condiciones de contorno de
Recibido en versión revisada el 30 de la temperatura se especifican para una duración de 24 horas, y las simulaciones no estacionarias
enero de 2021 Aceptado el 4 de febrero también se realizan para una duración de 24 horas. Las tasas de evaporación y condensación se
de 2021
obtienen a partir del gradiente de la fracción de masa de vapor de agua en las superficies de
Disponible en línea xxxx
evaporación (inferior) y condensación (superior). Para ello se utilizan varias funciones definidas por el
usuario (UDF). El método proporciona una producción de destilado instantánea y diaria. El método no
Palabras
es costoso desde el punto de vista computacional y puede ayudar a evaluar la producción de destilado
clave:
de una variedad de configuraciones cuando se conocen las temperaturas de las superficies inferior y
Destilació
n solar superior.
CFD 2021 Elsevier Ltd. Todos los derechos reservados.
Rendimiento del destilado Selección y revisión por pares bajo la responsabilidad del comité científico de la International Confer-
ence on Technological Advancements in Materials Science and Manufacturing.

computacionalmente más exigentes. En este trabajo, como primer

1. Introducción
paso, la modelización es

Los alambiques solares de pendiente simple son Dirección de correo electrónico: gauravmittal@geu.ac.in

configuraciones muy sencillas que se utilizan para obtener agua

destilada a partir de salmuera [1-4]. Estos alambiques suelen
emplear una cubierta de vidrio en la parte superior y un depósito
de agua en la parte inferior. El vidrio es casi transparente a la
radiación solar entrante. La radiación transmitida es captada por el
depósito de agua, que calienta el agua. La opacidad del cristal a la
radiación térmica de longitud de onda más larga de la cubeta
impide la salida de la radiación térmica de la cubierta superior.
Este es el efecto invernadero. La superficie superior de cristal
pierde calor a su alrededor por convección y radiación térmica y se
mantiene a una temperatura inferior a la de la cubeta, por lo que
cumple la función de superficie de condensación. Los vapores de
agua se elevan desde el recipiente debido a la convección natural y
se condensan en la superficie superior. El destilado se desliza por
la cubierta superior y se recoge en un canal lateral.
El efecto de varias configuraciones en la productividad de los
alambiques se ha estudiado experimental y computacionalmente
[por ejemplo, 1-6]. Las simulaciones CFD ofrecen medios rentables
de exploración. Los estudios CFD se realizan principalmente para
configuraciones 2 D y a veces para 3D [1- 8]. Además, algunos
enfoques se basan en la prescripción de las temperaturas en las
superficies inferior y superior como condiciones de contorno,
mientras que los enfoques más elaborados incluyen la radiación
también en el modelado y prescinden de la especificación de las
temperaturas como condiciones de contorno. Estos enfoques son
 temperatura se basan en [6]. Las demás condiciones de contorno
se realiza sólo para la zona húmeda (mezcla de agua + aire)
son las siguientes. Las paredes laterales están aisladas, y se
dentro del alambique con temperaturas especificadas como
especifica un gradiente cero para las especies y ningún
condiciones de contorno. El modelo se implementa en ANSYS
deslizamiento para el momento. En las caras superior e inferior, se
Fluent y la simulación se realiza durante todo el día para obtener
supone que el aire está saturado debido a la condensación y
el rendimiento diario del destilado. Los detalles computacionales
evaporación, respectivamente. Así pues, la fracción másica de
y los resultados se describen en las siguientes secciones.
vapor de agua correspondiente al aire saturado a la temperatura
de la superficie se especifica en las superficies inferior y superior.
La especificación de las condiciones de contorno no estacionarias
2. Modelización computacional para la temperatura y la fracción de masa de vapor de agua se
realiza mediante UDF.
La geometría 2D del destilador solar de pendiente única se ha En el dominio se resuelven las ecuaciones de continuidad,
tomado de [8]. El alambique tiene una anchura de 75 cm. La momento, transporte de especies y energía. Se utiliza el modelo de
altura del lado izquierdo es de 8 cm y el ángulo de inclinación de turbulencia RNG k - s. El aire húmedo se trata como un gas ideal y
la cubierta de vidrio es de 30o . En la Fig. 1 se muestra la rejilla de se aplican las ecuaciones termofísicas.
cálculo. En la Fig. 2 se muestran los perfiles de temperatura del aire y del vapor de agua se han tomado de [9].
especificados en las superficies inferior y superior. Los perfiles de

https://doi.org/10.1016/j.matpr.2021.02.090
2214-7853/© 2021 Elsevier Ltd. Todos los derechos reservados.
Selección y revisión por pares bajo la responsabilidad del comité científico de la Conferencia Internacional sobre Avances Tecnológicos en Ciencia de Materiales y Fabricación.

Citar este artículo como: G. Mittal, Unsteady CFD modeling of a single slope solar still, Materials Today: Proceedings, https://doi.org/10.1016/j.
matpr.2021.02.090
 G. Mittal Materiales hoy: Proceedings xxx (xxxx) xxx

Dm = 2,6x10—5 101325 T 2

P 298
Se utiliza una UDF para especificar la difusividad basándose en
la ecuación anterior. Las ecuaciones de gobierno se discretizan
utilizando un esquema de viento ascendente de segundo orden
para el momento, la continuidad, las especies, la energía y la
turbulencia. El método implícito de segundo orden se utiliza para
la diferenciación temporal. El esquema de escalonamiento de
presión se utiliza para la interpolación de presión junto con el
algoritmo SIMPLE.
Los cálculos no estacionarios se realizan para una duración de
24 h. Se especifican como condiciones iniciales una temperatura
uniforme de 302 K y aire saturado en reposo a las 7 de la mañana.
Para garantizar la independencia del paso de tiempo y del tamaño
de la malla, las simulaciones se realizaron para un intervalo de
pasos de tiempo de 3 a 10 s en mallas de cálculo de 12.150 a
42.000 celdas. Para todas las distribuciones de malla, la capa límite
Fig. 1. Dominio computacional y distribución de la malla.
en la zona de aire húmedo se refinó utilizando 15 capas de
"inflación" con las celdas cercanas a la pared de un tamaño de 0,2
mm. La resolución de la capa límite es fundamental, ya que las
tasas de evaporación y condensación dependen de los gradientes
de la capa límite. La diferencia en el rendimiento diario de
destilación simulado entre el paso de tiempo y el tamaño de malla
más pequeños y el más grandes fue inferior al 1%.

3. Resultados

Las simulaciones se realizaron durante 24 horas, desde las 7 de

la mañana hasta las 7 de la mañana del día siguiente, utilizando las
temperaturas especificadas como condiciones límite. La Fig. 3
muestra el rendimiento instantáneo simulado del destilado, así
como el rendimiento diario acumulado. El rendimiento se calcula a
partir de la tasa de condensación de la siguiente manera.

Fig. 2. Condiciones límite de temperatura especificadas.

Rendimiento instantáneo de destilado (kg/m2 - hr) = m_ c

× 3600/Superficie de la cuenca

Las tasas de evaporación y condensación durante cada iteración El rendimiento instantáneo aumenta a medida que aumenta la
se determinan utilizando UDF que integra el gradiente de la temperatura de la superficie inferior, así como la diferencia de
fracción de masa de vapor sobre la superficie inferior y superior, temperatura entre la superficie superior y la inferior. Esto se debe
respectivamente. a un aumento de la fracción másica de vapor de agua, así como a la
dA
corriente de convección. El rendimiento instantáneo de destilado
alcanza su punto máximo en torno a las 15 horas y disminuye a
partir de entonces a medida que disminuye la temperatura.
El campo de velocidad simulado y las líneas de corriente, a 2 pm
se muestran en la Fig. 4. Tres estructuras vorticales, generadas
debido a con-
Velocidad de evaporación, ∫ .-→ vección, puede verse en las líneas de corriente. Una característica
interesante es
m_ e = - qDm �Y donde la integración es
→ que la superficie de la cubierta de vidrio está atravesada por dos
sobre la superficie inferior. ∫ dA donde la integración es
= qDm �Y .- sobre vórtices que giran en sentidos opuestos. Esto puede ser perjudicial
Velocidad de condensación, para el rendimiento del destilado porque el giro en el sentido de
m_ c las agujas del reloj del vórtice del lado izquierdo puede
la superficie superior. potencialmente interrumpir el flujo del destilado por la pendiente
de
La difusividad molecular en la zona de aire húmedo se la cubierta de cristal y puede inducir la caída de gotas en el agua
calcula como

2
G. Mittal Materiales hoy: Proceedings xxx (xxxx) xxx

Fig. 3. Rendimiento de destilado por hora y acumulado a lo largo de todo un día.

3
G. Mittal Materiales hoy: Proceedings xxx (xxxx) xxx

Fig. 4. Resultados simulados a 2 pm. Resultados simulados a 2 pm. (a) campo de velocidades (b) líneas de
corriente.

Fig. 5. Resultados simulados a 2 pm. Resultados simulados a 2 pm. (a) campo de temperatura (b) campo de
densidad (c) fracción másica de vapor de agua.

4
G. Mittal Materiales hoy: Proceedings xxx (xxxx) xxx

Fig. 6. Fracción másica de vapor de agua a las 2 pm a lo largo de (a) la línea vertical que pasa por el centro del alambique (b) la línea horizontal a una altura de 4 cm desde el
fondo (c) vista ampliada a lo largo de la línea horizontal cerca de la pared lateral izquierda (d) vista ampliada a lo largo de la línea horizontal cerca de la pared lateral
derecha.

escribieron varias funciones definidas por el usuario. El modelo

piscina. Tales efectos no pueden examinarse en el marco del
utiliza condiciones de contorno inestables de temperatura en la
presente modelo CFD.
superficie superior e inferior como entrada y proporciona el
En la Fig. 5 se muestran el campo de temperatura simulado, el
rendimiento del destilado para un día entero. Aparte de la
campo de densidad y la fracción másica de vapor de agua a las 2
implementación del modelo, las cuestiones de malla y la
pm. Se observa que el aire húmedo en la superficie inferior tiene
independencia del paso de tiempo fueron
menor densidad y mayor fracción másica de agua que en la
superficie superior. El aire de baja densidad asciende desde la
superficie inferior y transporta vapor con él, que a su vez se
condensa en la superficie superior. En la Fig. 6(a) se muestra un
gráfico de la fracción másica de vapor de agua a lo largo de una
línea vertical que pasa por el centro del alambique a las 2 pm. En la
figura se puede observar un fuerte gradiente en la fracción de
masa de vapor cerca de las superficies inferior y superior. Estos
gradientes son responsables de la evaporación y la condensación,
respectivamente. La fracción de masa de vapor a lo largo de una
línea horizontal a una altura de 4 cm desde la superficie inferior se
muestra en la Fig. 6(b). Una mayor fracción de masa cerca de la
pared izquierda se debe al flujo ascendente de aire a lo largo de
esta superficie, que transporta aire húmedo desde la superficie
inferior. Las Fig. 6(c) y (d) muestran vistas ampliadas cerca de las
paredes izquierda y derecha a lo largo de la línea horizontal. Se
incluyen para reforzar que el gradiente cero de la fracción de masa
de las especies se cumple y que la capa límite está bien resuelta.

4. Conclusiones

Se implementó un modelo CFD para modelar una destilación

solar de pendiente simple utilizando ANSYS Fluent. Para ello se

5
abordado. Además, la simulación de muestras muestra que el
modelo
G. Mittal describe muy bien las características del flujo en el Materiales hoy: Proceedings xxx (xxxx) xxx
interior del alambique. Es una herramienta valiosa para explorar
varias configuraciones que difieren en los parámetros
geométricos, la disposición de los deflectores, etc., e intentar
optimizar el destilador para obtener un alto rendimiento.

Declaración de contribución a la autoría de CRediT

Gaurav Mittal: Conceptualización, Metodología, Software,

Investigación, Redacción - revisión y edición.

Declaración de conflicto de intereses

Los autores declaran que no tienen intereses económicos ni

relaciones personales que pudieran influir en el trabajo
presentado en este artículo.

Referencias

[1] H. Taghvaei, H. Taghvaei, K. Jafarpur, M. Feilizadeh, M.R. Karimi Estahbanati,

Experimental investigation of the effect of solar collecting area on the
performance of active solar stills with different brine depths, Desalination 358
(2015) 76-83.
[2] M.M. Rahman, H.F. Öztop, A. Ahsan, M.A. Kalam, Y. Varol, Double-diffusive
natural convection in a triangular solar collector, Int. Commun. Heat Mass
Transf. 39 (2) (2012) 264-269.
[3] N. Rahbar, J.A. Esfahani, E. Fotouhi-Bafghi, Estimation of convective heat
transfer coefficient and water-productivity in a tubular solar still - CFD
simulation and theoretical analysis, Sol. Energy 113 (2015) 313-323.
[4] M. Edalatpour, A. Kianifar, S. Ghiami, Effect of blade installation on heat transfer
and fluid flow within a single slope solar still, Int. Commun. Heat Mass Transfer
66 (2015) 63-70.

6
G. Mittal Materiales hoy: Proceedings xxx (xxxx) xxx

[5] R. Alvarado-Juárez, J. Xamán, G. Álvarez, I. Hernández-López, Numerical study of [7] M. Keshtkar, M. Eslami, K. Jafarpur, Effect of design parameters on performance
heat and mass transfer in a solar still device: Efecto de la cubierta de vidrio, of passive basin solar stills considering instantaneous ambient conditions: A
Desalación 359 (2015) 200-211. transient CFD modelling, Sol. Energy 201 (2020) 884-907.
[6] M. Keshtkar, M. Eslami, K. Jafarpur, A novel procedure for transient CFD [8] M. Feilizadeh, M. Soltanieh, M.R. Karimi Estahbanati, K. Jafarpur, S.-S.
modeling of basin solar stills: Coupling of species and energy equations, Ashrafmansouri, Optimización de las dimensiones geométricas de los
Desalination 481 (2020) 114350. alambiques solares tipo cuenca de una sola pendiente, Desalinización 424
(2017) 159-168.
[9] T.L. Bergman, A.S. Lavine, F.P. Incropera, D.P. Dewitt, Fundamental of Heat and
Mass Transfer, Wiley, 2011.