Artículo 979 Análisis
Artículo 979 Análisis
Artículo 979 Análisis
Resumen
Este artículo muestra la simulación de un sistema de concentración solar partiendo de la
restauración de un disco parabólico de 1.80 m de diámetro y 0.675 m de distancia focal, el cual se
encuentra en el Instituto Tecnológico de Apizaco (ITA), Tlaxcala, México. En este estudio se
construyó una base para el disco parabólico y se recubrió el colector con aluminio de alta
reflectividad (Prominox ASTM 240), posteriormente el sistema se simuló en el software de trazado
de rayos “Tonatiuh” con la finalidad de analizar la distribución del flujo de la radiación en el punto
focal del sistema y con la imagen focal obtener de primera instancia las dimensiones del diámetro
del receptor. Para un error de 4.39 mrad se obtuvo una potencia máxima de 1829.03 W con un
radio del receptor de 0.06 m. Con el error de 6.7 mrad se obtuvo una potencia máxima de 1818.52
W que corresponde al radio del receptor de 0.075 m. Se observó que a medida que el radio del
receptor es mayor, la potencia aumenta y el desbordamiento disminuye, pero al tener un receptor
de mayores dimensiones el factor de sombra se eleva, el cual afecta el desempeño óptico del
sistema.
Palabras clave: Concentrador solar, disco parabólico, energía termosolar, simulación, trazado de
rayos.
76
Moisés Sánchez Espejel y Vicente Flores Lara / Ingeniería Vol. 25-2 (2021) 75-88.
Para el diseño de los sistemas de concentración diseño de un sistema de bajo costo para
solar y su análisis óptico, se utilizan softwares procesos térmicos de media temperatura.
como Tonatiuh o SolTrace, los cuales se basan Rodríguez et al. (2014) utilizaron el trazado de
en el método de trazado de rayos de Monte- rayos para diseñar un colector solar térmico que
Carlo (MCRT). En estos programas el usuario permita generar energía uniformemente
puede crear el concentrador solar de acuerdo durante el año.
con la geometría deseada, así como modificar
sus propiedades ópticas, la inclinación, forma En este trabajo se busca evaluar la capacidad de
solar, número de rayos del sol, entre otros concentración y la distribución geométrica de
factores (National Renewable Energy la energía concentrada en el receptor de un
Laboratory, 2020). sistema de disco parabólico de 1.80 m de
diámetro mediante las simulaciones de trazado
Las simulaciones de trazado de rayos son muy de rayos que nos demuestren como los errores
efectivas para el diseño y la optimización de ópticos y de seguimiento influyen en el
parámetros de los sistemas de concentración rendimiento del concentrador, con el propósito
solar, ya que en ellas se pueden analizar las de conocer la potencia y las dimensiones del
alteraciones en los sistemas y evaluarlas receptor que mejor aprovechen la energía solar.
fácilmente. Los softwares de trazado de rayos Posteriormente con ayuda de los resultados de
se han utilizado para simular el fenómeno de este trabajo, realizar una propuesta de un
recepción y reflexión en las superficies, receptor de cavidad y de un sistema de
incluyendo la forma solar, inclinación, seguimiento solar, que en conjunto formen un
propiedades ópticas, número de rayos, entre concentrador solar de disco parabólico
otros factores. funcional para su estudio e investigación en el
ITA.
Diversos autores han realizado simulaciones de
trazado de rayos para sistemas de
concentración solar con el objetivo de predecir
su comportamiento y mejorar su rendimiento, 2. Marco Teórico
tal es el caso de Supa et al. (2014), quienes 2.1 Método de Trazado de rayos de Monte-
realizaron un estudio de trazado de rayos para Carlo
observar cómo afecta el ángulo de borde de un Los softwares que se basan en el método
concentrador de disco parabólico en su imagen MCRT utilizan los principios de la geometría
focal. Pavlovic y Stefanovic (2015) simularon óptica, así como el uso de cálculos estadísticos
un disco parabólico para analizar los aspectos para la simulación de los sistemas de
geométricos de su imagen focal, concentración solar, generando rayos a partir
posteriormente obtuvieron el flujo de la de una fuente simulada, permitiendo observar
irradiancia lo que les permitió establecer los la interacción de los rayos en las superficies.
parámetros de diseño de su receptor. Terrón et Este tipo de programas son muy útiles para el
al. (2018) efectuaron un análisis de un CPC, diseño y análisis de los concentradores solares
mediante una simulación del sistema, (Giovinazzo et al., 2014).
observaron la interacción de los rayos del Sol
con el concentrador y obtuvieron la energía El código de MCRT trabaja con un número
disponible en el receptor en un tiempo específico de rayos, cada uno se traza a través
determinado. Pavlovic et al. (2014) de un sistema definido, mientras se va
examinaron el flujo en el receptor de un encontrando varias interacciones ópticas.
concentrador de disco parabólico mediante un Algunas de estas interacciones son de
análisis óptico de trazado de rayos para el naturaleza probabilística y otras son
deterministas. Tal código tiene la ventaja de
77
Moisés Sánchez Espejel y Vicente Flores Lara / Ingeniería Vol. 25-2 (2021) 75-88.
replicar la interacción de fotones reales y, por Para el cálculo del factor de distribución 𝑄𝑖𝑗 se
lo tanto, puede proporcionar resultados emplea la siguiente ecuación:
precisos para sistemas complejos que no se
pueden modelar de otra manera. La precisión
aumenta con el número de rayos trazados, pero 𝑄𝑖𝑗 = 𝑄𝑖 𝐷𝑖𝑗 (3)
un mayor número de rayos significa más
tiempo de procesamiento (Wendelin et al.,
2013). 2.2 Parámetros ópticos-geométricos
78
Moisés Sánchez Espejel y Vicente Flores Lara / Ingeniería Vol. 25-2 (2021) 75-88.
79
Moisés Sánchez Espejel y Vicente Flores Lara / Ingeniería Vol. 25-2 (2021) 75-88.
por la empresa “antenas de México”, con los utilizó PTR, placas y tubos de acero para la
parámetros principales que se muestran en la construcción de la base del disco y se usaron
tabla 1. chumaceras de piso, así como de pared para que
el sistema pudiera tener movimiento sobre sus
2 ejes. El recubrimiento del disco fue hecho con
Para este estudio solo se disponía de un disco aluminio de alta reflectividad (Prominox
parabólico, por lo que fue necesario la ASTM 240), el cual tiene un valor de 0.85 de
construcción de una base que fuera capaz de reflectividad. En la figura 3 se muestra el disco
mover el concentrador sobre 2 de sus ejes. Se montado en la base y recubierto con el material
reflejante.
Figura 3. Concentrador de disco parabólico, colocado en su base y recubierto con aluminio de alta
reflectividad.
80
Moisés Sánchez Espejel y Vicente Flores Lara / Ingeniería Vol. 25-2 (2021) 75-88.
3.2 Simulación del sistema Tonatiuh permite elegir dos formas solares:
Pillbox y Buie, para este análisis se eligió la
La simulación se llevó a cabo empleando el distribución Buie, ya que es más cercana a la
software Tonatiuh, con el fin de analizar la forma real del sol. Considerando que en este
distribución del flujo de la radiación en el punto caso se está analizando un disco restaurado, el
focal del sistema y con la imagen focal obtener cual no es perfecto, se consideran los valores
de primera instancia el diámetro del receptor. máximos para los errores ópticos previamente
Se estableció un sistema óptico, el cual está mencionados, resultando en un error total de
compuesto de tres elementos: Sol,
𝜎 = 4.39 𝑚𝑟𝑎𝑑. Los errores del sistema de
concentrador y receptor.
seguimiento no serán considerados, ya que en
De acuerdo con los datos proporcionados por el este documento únicamente se está analizando
Sistema Meteorológico Nacional (2020), para el disco parabólico y su receptor. Los
el estado de Tlaxcala, en el mes de mayo de parámetros principales de entrada para realizar
2020 se tuvo un promedio aproximado de la simulación se muestran en la tabla 2.
850 𝑊/𝑚2 de radiación solar entre las 16:00 y El disco parabólico se caracterizó, tomando en
las 17:00 horas. Se eligió la elevación de 90°, cuenta sus propiedades ópticas y dimensiones.
ya que se considera que el concentrador contará El radio del receptor comenzó con un valor
con un sistema de seguimiento solar, el cual aleatorio y posteriormente se fue variando para
buscará que los rayos incidan analizar el comportamiento del flujo. Un
perpendicularmente en la superficie del disco. trazado de 1000 rayos fue necesario para
Se utilizó el valor de 0.85 para la reflectividad, comprobar que los rayos llegaran a su objetivo
la cual corresponde al material Prominox (figura 4).
ASTM 240.
81
Moisés Sánchez Espejel y Vicente Flores Lara / Ingeniería Vol. 25-2 (2021) 75-88.
Se ejecutó el análisis del flujo en el receptor. se puede notar que la distribución es simétrica
Para obtener un resultado más preciso en este a lo largo del eje X y del eje Y. También se
proceso, se utilizó la máxima cantidad de rayos puede observar que el flujo máximo se
(7500000). La gráfica de la distribución de encuentra en el centro del receptor y disminuye
flujo incidente para un receptor con radio de progresivamente cuando se aleja del centro.
0.025 m se muestra en la figura 5. Visiblemente
Las figuras 6 y 7 muestran las gráficas en donde nuevamente se verifica su simetría en los dos
se percibe de manera más detallada el planos y se puede observar que la mayor parte
comportamiento del flujo en el receptor en los del flujo se encuentra en un diámetro
sectores vertical y horizontal. En ellas aproximado de 0.016 m.
82
Moisés Sánchez Espejel y Vicente Flores Lara / Ingeniería Vol. 25-2 (2021) 75-88.
83
Moisés Sánchez Espejel y Vicente Flores Lara / Ingeniería Vol. 25-2 (2021) 75-88.
Figura 8. Resultados de la simulación del disco parabólico con un error de 4.39 mrad.
84
Moisés Sánchez Espejel y Vicente Flores Lara / Ingeniería Vol. 25-2 (2021) 75-88.
Figura 10. Resultados de la simulación del disco parabólico con un error de 4.39 mrad.
1818.52 W con un radio del receptor de 0.075 sombra aumenta y la relación de concentración
m. De acuerdo con la potencia obtenida en la disminuye.
simulación del disco parabólico, el posible uso
del sistema se enfocaría para fines domésticos, Con esta investigación fue posible conocer el
ya que diversos autores han obtenido resultados comportamiento del concentrador de disco
similares, tal es el caso de André-Charles et al. parabólico de acuerdo con sus propiedades
(2019), los cuales obtuvieron una potencia de óptico-geométricas y se pudo observar la
1660 W utilizando un concentrador solar para variación de la potencia y la imagen focal
proveer de agua caliente a una casa habitación. dependiendo de los errores ópticos y de
Y Arteconi et al. (2019), quienes consiguieron seguimiento. Con los resultados es posible
una potencia de 2000 W utilizando la energía continuar la investigación con la propuesta de
termosolar para el suministro de agua caliente un receptor de cavidad partiendo de la imagen
de un complejo residencial de apartamentos. focal y potencia obtenidos en la simulación,
posteriormente el diseño de un sistema de
Por otro lado, se puede concluir que los errores seguimiento solar preciso que presente el
ópticos no afectan demasiado a la potencia ya menor error posible.
que solo hubo una diferencia del 0.57% entre
las potencias de ambas pruebas, pero las
dimensiones de la imagen focal si sufren Agradecimientos
cambios más drásticos, ya que hubo una
diferencia del 20% entre los radios del receptor Al Consejo Nacional de Ciencia y Tecnología
de ambas pruebas. Esto se debe a que, al (CONACYT) por el apoyo económico para
aumentar el error, los rayos del Sol tienden a cursar un posgrado y realizar una estancia de
tener una mayor dispersión sobre la superficie investigación en el extranjero
del receptor, por lo tanto, es necesario una
Asimismo, al departamento de posgrado del
mayor área para obtener una potencia más alta.
Instituto Tecnológico de Apizaco, por brindar
Se debe tener en cuenta que entre más grandes
el apoyo necesario para la realización de este
sean las dimensiones del receptor, el factor de
proyecto
Referencias
Abas, N., Kalair, A. and Khan, N. (2015). Review of fossil and future energy technologies. Futures,
69, 31-49. doi: 10.1016/j.futures.2015.03.003
Ahmed, A., Tyurina, Y. G., Smailova, Z. P., Kurilova, A. A. and Shulus, A. A. (2019). Government
Policy and Financing Options for Solar Energy: World Prospects. International Journal of Energy
Economics and Policy, 9, 131-145. doi: 10.32479/ijeep.8212
Lentz, A., Rincón, E. y Estrada, V. (2012). Presente y Futuro de la Energía Solar en México. XV
Congreso Ibérico y X Congreso Iberoamericano de Energía Solar, Galicia, España.
Secretaria de Energía (SENER). Prospectiva de Energías Renovables 2012-2026.
Muhammad-Sukki, F., Ramirez-Iniguez, R., McMeekin, S.G., Stewart, B.G. and Clive, B. (2010).
Solar Concentrators. IJAS, 1.
86
Moisés Sánchez Espejel y Vicente Flores Lara / Ingeniería Vol. 25-2 (2021) 75-88.
Reddy, G.C.M., Rao, G. A., Prasad, V.V.S.H and Mukkoti, V. V. (2017). Design of concentrated
parabolic dish collector for vapor absorption refrigeration system. International Journal of Civil
Engineering and Technology, 8, 714-724.
Yosra, K., Nizar, D., Aribia, B. and Abdallah H. (2013). Design of a solar dish concentrator
according to the needed energy for a given application. 14th international conference on Sciences
and Techniques of Automatic control & computer engineering, Sousse, Tunisia.
National Renewable Energy Laboratory (2020), SolTrace. [Online], Available:
https://www.nrel.gov/csp/soltrace.html
Supa, B. A., Zainudina, M. F., Zanariah, T., Alia, S., Bakara, R. A. and Minga, G. L. (2014). Effect
of rim angle to the flux distribution diameter in solar parabolic dish collector. 2nd International
Conference on Sustainable Energy Engineering and Application, ICSEEA.
Pavlovic, S. R. and Stefanovic, V. P. (2015). Ray Tracing Study of Optical Characteristics of the
Solar Image in the Receiver for a Thermal Solar Parabolic Dish Collector. Journal of Solar Energy,
2015, Article ID 326536
Terrón, M., Peña, M. I., Carrillom, J. G., Ayala, U. D. and Flores, V. (2018). Solar Ray Tracing
Analysis to Determine Energy Availability in a CPC Designed for Use as a Residential Water
Heater. Energies, 11.
Pavlovic, S., Vasiljevic, D. and Stefanovic, V. (2014). Optical Design of a Solar Parabolic Thermal
Concentrator Based on Trapezoidal Reflective Petals. Conference paper. doi:
10.2298/TSCI160213104Y
Rodriguez, D., Rosengarten, G., Belmonte, J. F., Izquierd,o M., Molina A. and Almendros, J. A.
(2014). Ray Tracing of a Solar Collector Designed for Uniform Yearly Production. Energy
Procedia, 57, 2221-2230. doi: 10.1016/j.egypro.2014.10.189
Giovinazzo, C., Bonfiglio, L., Gomes, J. and Karlsson, B. (2014). Ray Tracing Modelling of an
Asymmetric Concentrating PVT. ISES, EuroSun 2014.
Wendelin, T., Dobos, A. and Lewandowski, A. (2013). SolTrace: A Ray- Tracing Code for
Complex Solar Optical Systems”, Technical Report, NREL/TP-5500-59163
Mahan, J. R. (2019). The Monte Carlo Ray-Trace Method in Radiation Heat Transfer and
Applied Optics. Virginia, USA, ISBN 9781119518525
Duffie, J. A. and Beckman, W. A. (2013). Solar Engineering of Thermal Processes. Solar Energy
laboratoy, University of Wisconsin, USA.
Dincer, I., Zamfirescu C. (2014). Advanced Power Generations Systems. 369-453. doi:
https://doi.org/10.1016/C2009-0-63866-6
Stine, B. W. and Greyer, M. (2001). Power from the sun. [Online], Available:
http://www.powerfromthesun.net/book.html
Sistema Meteorológico Nacional. (2020)., Gobierno de México, [Online], Available:
https://smn.conagua.gob.mx/es/estaciones-meteorologicas-automaticas
87
Moisés Sánchez Espejel y Vicente Flores Lara / Ingeniería Vol. 25-2 (2021) 75-88.
Blázquez, R., Carballo, J., Cadiz, P., Frasquet, M., Silva, M., Fontela, P. and Ballesteros, J. C.
(2015). Optical test of the DS1 prototype concentrating Surface. Energy Procedia, 69, 41-49. doi:
10.1016/j.egypro.2015.03.006
André-Charles, M. A., Gael, L., Norman, H., Gabriel, H. and Arthur, L. (2019). Low temperature
and power solar energy conversion for domestic use. Energy Procedia, 161, 454-463. doi:
10.1016/j.egypro.2019.02.072
Arteconi, A., Zotto, L., Tascioni, R., Mahkamov, K., Underwood, C., Cabeza, L. F., Gracia, A.,
Pili, P., Mintsa, A. C., Bartolini, C. M., Gimbertani, T., Botargues, T., Halimic, E. and Cioccolanti,
L. (2019). Simulation analysis of an innovative micro-solar 2kW Organic Rankine Cycle plant
coupled with a multi-apartment building for domestic hot water supply. Energy Procedia, 158,
2225-2230. doi: 10.1016/j.egypro.2019.01.168
88