Ingeniería

ISSN: 1665-529X
ingenieria@[Link]
Universidad Autónoma de Yucatán
México

Cetina-Quiñones, A. J.; Bassam, A.; Hernández-Chan, G.; Hernández Benítez, J.;

Hernández Reyes, I.; Lugo Chávez, D.
Modelación térmica de un colector solar de canal parabólico mediante el método de
elementos finitos
Ingeniería, vol. 21, núm. 1, 2017, pp. 1-12
Universidad Autónoma de Yucatán
Mérida, México

Disponible en: [Link]

Cómo citar el artículo

Número completo
Sistema de Información Científica
Más información del artículo Red de Revistas Científicas de América Latina, el Caribe, España y Portugal
Página de la revista en [Link] Proyecto académico sin fines de lucro, desarrollado bajo la iniciativa de acceso abierto
 Modelación térmica de un colector solar de canal parabólico
mediante el método de elementos finitos
Cetina-Quiñones, A. J.1, Bassam, A.2, *, Hernández-Chan, G.3, Hernández Benítez, J.3, Hernández Reyes, I.4,
Lugo Chávez, D.4

Fecha de recepción: 23 de abril de 2017 – Fecha de aprobación: 16 de junio de 2017

Resumen
El presente trabajo tiene como propósito presentar el modelo de un Colector Solar de Canal Parabólico (CCP)
empleando el Método de Elementos Finitos con el fin de predecir el comportamiento térmico del fluido de trabajo a
lo largo del tubo receptor del colector. Así mismo, se estima la eficiencia térmica con base a las ecuaciones
gobernantes que intervienen en los procesos de transferencia de calor. Para validar los resultados del modelo, se
realizó una simulación térmica del fluido empleando el software SolidWorks. El error máximo obtenido de la
comparación de la modelación con la simulación fue del 7.6% a un caudal de 1 l/min. De acuerdo a los resultados
obtenidos de los errores estadísticos, el método puede predecir eficazmente la temperatura del fluido en caudales
altos. El modelo desarrollado puede ser útil como herramienta de diseño, en la optimización del tiempo empleado en
las simulaciones generadas por el software y en la minimización de los costos de manufactura referentes a Colectores
Solares de Canal Parabólico.

Palabras clave: energía solar, simulación numérica, sistemas fototérmicos, eficiencia térmica

Abstract
The purpose of the present work is to present a model of a parabolic trough solar collector (PTC) using the finite
element method in order to predict the thermal behavior of the working fluid along the receiver tube of the collector.
Likewise, the thermal efficiency is estimated bases on the governing equations involved in heat transfer processes.
To validate the results of the modeling, a thermal simulation of the fluid is performed using the SolidWorks
software. The maximum error obtained from the modeling comparison with the simulation was 7.6% at a flow rate of
1 l/min. According to the results obtained from statistical errors, the method can effectively predict the temperature
of the fluid at high flow rates. The developed model can be useful as a design tool, as well as for comparison with the
experimental data, in the optimization of the time employed in the simulations generated by the software and in the
manufacturing costs of references to the parabolic trough solar collectors

Keywords: Solar energy, numerical simulation, photothermal systems, thermal efficiency

______________________
1
Estudiante de Licenciatura de Ingeniería en Energías Renovables, Facultad de Ingeniería-UADY, Av. Industrias
No Contaminantes por Anillo Periférico Norte s/n Apdo. Postal 150 Cordemex, Mérida, Yucatán, México
2
Facultad de Ingeniería-UADY, Av. Industrias No Contaminantes por Anillo Periférico Norte s/n Apdo. Postal 150
Cordemex, Mérida, Yucatán, México
3
Instituto Tecnológico de Mérida, Av. Tecnológico Km 4.5 S/N C.P 97118, Mérida, Yucatán México.
4
División de Ingeniería Electromecánica del Instituto Tecnológico Superior de Centla, Calle Ejido S/n., Col. Siglo
XXI, Frontera, Centla, Tabasco. C.P. 86750.
*
Autor de correspondencia: baali@[Link]
Nota: Este artículo de investigación es parte de Ingeniería Revista Académica de la Facultad de Ingeniería,
Universidad Autónoma de Yucatán, Vol. 21, No. 1, 2017, ISSN 2448-8364.
 Cetina-Quiones et al. / Ingeniería 21-1 (2017) 1-12

INTRODUCCIÓN receptores con cubierta evacuado. Estudios más

La contaminación ambiental por gases de efecto recientes como el de Uzgoren (2016), en el cual
invernadero es un problema que se ha presentado a lo evaluó un CCP en estado transitorio mediante un
largo de estas últimas décadas. Estos gases se modelo matemático desarrollado por el método del
producen por diferentes factores, dentro de los que volumen finito, validándolo con datos
destacan la combustión de carburantes fósiles, el experimentales, obteniendo errores no significativos.
vapor de agua, el uso de fertilizantes en los campos Finalmente, Jaramillo (2016) realizó un análisis de la
agrícolas, el venteo del gas natural, entre otros. El eficiencia térmica de un CCP mediante la inserción
aumento considerable de estos gases, como el dióxido de cintas retorcidas dentro del tubo receptor,
de carbono (CO2), dióxido de azufre (SO2), y óxido obteniendo un aumento del 3.5 % en la eficiencia
nitroso (NOx) llevan a devastadoras consecuencias térmica bajo condiciones controladas.
que repercuten directamente en el cambio climático
(Kalogirou, 2004). Por esta razón, existe una La caracterización de un CCP es de suma
tendencia mundial actual centrada en el uso de importancia, ya que permite conocer los diferentes
fuentes de energía alternativas como medida para fenómenos que se presentan en éste. Uno de estos, es
mitigar las emisiones de estos gases nocivos y otros el comportamiento de la temperatura a lo largo del
contaminantes a la atmósfera (May et al. 2017). tubo receptor, desde la entrada hasta la salida de éste.
Esto permite observar el fenómeno de la trasferencia
La energía solar es una de las fuentes energéticas de calor del tubo receptor hacia el fluido, así como
renovables que atraen más atención debido a su también predecir el cambio de fase de dicho fluido,
abundancia, la limpieza y el hecho de que no genera garantizando la obtención de vapor para la generación
ninguna contaminación (May et al. 2016). Dentro de de electricidad. Del mismo modo, la modelación
esta fuente de energía se encuentra la tecnología de permite caracterizar un CCP bajo diferentes
los concentradores solares la cual es una de las más condiciones sin necesidad de recurrir a la
desarrolladas actualmente y empleadas para la experimentación y tener un concepto general del
generación de energía eléctrica (Fernández-García et comportamiento térmico del fluido.
al., 2010). Los colectores solares de tipo concentrador
son capaces de producir altas temperaturas (más de El objetivo del presente trabajo, es modelar un CCP
400 °C). Esto los convierte en una tecnología factible mediante la primera ley de la termodinámica
y prometedora para la desalinización solar, empleando el Método de Elementos Finitos (MEF).
aplicaciones en química solar, refrigeración solar Este método permite obtener un perfil de
(absorción y adsorción), producción de hidrógeno temperaturas del fluido, mediante la discretización del
solar y para Plantas Solares Concentradas (CSP) tubo receptor en un número finito de elementos
(Tzivanidis et al. 2015). llamados nodos. De igual manera se presenta la
validación del modelo mediante una simulación de la
Un CCP capta la energía radiante del Sol y la temperatura del flujo implementado en el software
concentra en un tubo receptor localizado en el eje SolidWorks®.
focal de la parábola por el cual se hace circular un
fluido de trabajo, transformando la radiación solar en METODOLOGÍA
energía térmica útil. El análisis de transferencia de Un CCP es un sistema de concentración solar
calor de estos colectores es importante para el cálculo conformado por una lámina reflectora la cual capta
de las pérdidas térmicas y el dimensionamiento de la los rayos solares y los concentra en un tubo receptor
planta de energía solar durante el diseño preliminar situado en el eje focal de la parábola. En la Fig. 1 se
(Padilla et al. 2011). muestra el diseño y los componentes básicos de un
CCP.
Dada la importancia del análisis de transferencia de
calor en los CCP, se han desarrollado diferentes Dos de los parámetros importantes que caracterizan
modelos a lo largo de los últimos años. Nwosu (2009) un CCP son la temperatura en la salida del tubo
desarrolló un modelo térmico empleando el método receptor y la eficiencia térmica. El primero permite
de elementos finitos a un tubo receptor con cubierta describir el calor ganado por el fluido de trabajo, lo
de vidrio evacuado, en el cual obtuvo la mejor que da pie al desarrollo de prototipos que puedan
ganancia a flujos de calor altos. Por otro lado, Liang optimizarse en procesos que requieran temperaturas
(2015) comparó diferentes modelos de transferencia específicas. Por otra parte, la eficiencia térmica,
de calor en CCP por el cual concluyó la innecesaridad permite obtener la razón de la energía útil respecto a
de contemplar la conducción de calor en tubos la energía de entrada al sistema.
 Cetina-Quiones et al. / Ingeniería 21-1 (2017) 1-12

Tubo receptor

Soportes del CCP

Lamina reflectora

Figura 1. Componentes básicos de un CCP

Análisis térmico del CCP intercepción, 𝐴! es el factor geométrico del colector y

En un sistema de concentración solar, uno de los 𝜃 es el ángulo de incidencia de los rayos solares
parámetros importantes es la eficiencia óptica la cual (Duffie y Beckman, 2013). Para el análisis térmico
se calcula mediante la Ec. 1. del CCP se considera al flujo unifásico, siendo agua
dicho fluido de trabajo. De igual manera, se
𝜂! = 𝜌𝜏𝛼𝛾 1 − 𝐴! 𝑡𝑎𝑛 𝜃 𝑐𝑜𝑠 𝜃 (1) consideran despreciables las pérdidas por conducción
ya que se suponen los soportes del tubo receptor
aislados. En la Fig. 2 se presenta el modelo del tubo
donde 𝜌 es la reflectividad de la lámina reflectora, 𝜏
receptor y la representación de las pérdidas térmicas
es la transmisividad de la envolvente del tubo, 𝛼 es la
como una red de resistencias.
absortividad del tubo receptor, 𝛾 es el factor de

Figura 2. Modelo del tubo receptor:

(a) nomenclatura
(b) red de resistencias térmicas para la sección transversal

Para estimar las pérdidas de transferencia de calor por donde el coeficiente ℎ! representa la pérdida de
convección y radiación en el receptor, se debe hallar transferencia de calor por convección y ℎ! la pérdida
el coeficiente global de pérdidas térmicas térmica por radiación. ℎ! se calcula mediante la Ec. 3:
𝑈! calculado como:
𝑘!
ℎ! = 𝑁𝑢! (3)
𝑈! = ℎ! + ℎ! (2) 𝐷!

siendo 𝑘! la conductividad térmica de aire, 𝐷! el diámetro externo del tubo receptor y 𝑁𝑢! el número
 Cetina-Quiones et al. / Ingeniería 21-1 (2017) 1-12

de Nusselt el cual se puede estimar mediante la Ghajar, 2011):

ecuación de Churchill y Bernstein como (Cengel y

! ! ! !
0.62ℜ!! ! 𝑃𝑟!! ! ℜ! (4)
𝑁𝑢! = 0.3 + 1+
1 + 0.4 𝑃𝑟! ! ! ! ! 282000

donde 𝑃𝑟! es el número de Prandtl para el aire y ℜ! 𝜂! 𝐺! 𝐶!

es el número de Reynolds del aire, el cual se calcula 𝑇! = 𝑇!,! + (7)
ℎ!
mediante la Ec. 5:
siendo 𝑇!,! la temperatura máxima del fluido en la
𝑉𝐷! salida, 𝐺! la radiación solar directa, 𝐶! el factor de
ℜ! = (5)
𝜐 concentración y ℎ! el coeficiente de transferencia de
calor por convección en el interior del tubo receptor,
siendo 𝑉 la velocidad del viento y υ la viscosidad el cual se halla como:
cinemática del aire
𝑘!
Por otra parte, ℎ! se estima mediante la Ec. 6: ℎ! = 𝑁𝑢! (8)
𝐷!
ℎ! = 4𝜎𝜀! 𝑇!! (6)
de aquí, 𝑘! es la conductividad térmica de agua, 𝐷! es
el diámetro interno del tubo y 𝑁𝑢! el número de
donde 𝜎 es la constante de Stefan-Boltzmann, 𝜀! es la Nusselt del agua el cual depende del tipo de flujo, y
emisividad y 𝑇! la temperatura de tubo receptor, se calcula como (Gnielinski, 2013):
calculado como:
𝟒. 𝟑𝟔
𝟏 − 𝝃 𝑵𝒖𝒍𝒂𝒎,𝟐𝟑𝟎𝟎 + 𝝃𝑵𝒖𝒕𝒖𝒓𝒃,𝟒𝟎𝟎𝟎
𝒇𝒇 𝟐 𝕽𝒘 < 𝟐𝟑𝟎𝟎
𝕽𝒘 − 𝟏𝟎𝟎𝟎 𝟏 + 𝒅𝟑 𝑷𝒓𝒘 𝑷𝒓 𝟎.𝟏𝟏 𝟐𝟑𝟎𝟎 < 𝕽𝒘 < 𝟒𝟎𝟎𝟎
𝑵𝒖𝒘 = 𝟖 𝒘
𝕽𝒘 > 𝟒𝟎𝟎𝟎 (9)
𝟐 𝑷𝒓𝒗
𝒇𝒇
𝟏 + 𝟏𝟐. 𝟕 𝑷𝒓𝟑𝒘 − 𝟏
𝟖

De la Ec. 9, 𝑃𝑟! es el número de Prandtl del agua, 𝑑 ℜ! − 2300

𝜉= (10)
es un factor calculado como la razón del diámetro 4000 − 2300
interno de tubo entre su longitud total, 𝐷! es el
diámetro interno del tubo receptor y 𝜉 es un factor el De igual modo, 𝑓! es el factor de fricción y se estima
cual se calcula mediante la Ec. 10. mediante la ecuación de Chen (Fernández-García et
al., 2015):

𝟏
𝒇𝒇 =
𝟏 𝒓𝒓 𝟓. 𝟎𝟒𝟓𝟐 𝟏 𝒓𝒓 𝟏.𝟏𝟎𝟗𝟖 𝟓. 𝟖𝟓𝟎𝟔
𝟐 (11)
−𝟐𝒍𝒐𝒈 − 𝒍𝒐𝒈 +
𝟑. 𝟕𝟎𝟔𝟓 𝑫𝒊 𝕽𝒘 𝟐. 𝟖𝟐𝟓𝟕 𝑫𝒊 𝕽𝟎.𝟖𝟗𝟖𝟏
𝒘

siendo 𝑟! la rugosidad de material del cual está hecho 4𝑚

el tubo. ℜ! = (12)
𝜋𝐷! 𝜇!

Por otra parte, 𝑁𝑢!"#,!"## = 4.36, 𝑁𝑢!"#$,!""" se donde 𝑚 es el flujo másico y 𝜇! la viscosidad
calcula mediante la ecuación de flujo turbulento dinámica del agua.
(ℜ! > 4000), tomando ℜ! = 4000, 𝑃𝑟! = 0.7,
!"
𝑑 = 0.0001 y ! = 1. Posteriormente se calcula el factor de remoción 𝐹!
!"!
como (Kalogirou, 2013):
El número de Reynolds del agua (ℜ! ) se calcula
como: !!! !!! ! ! !!
𝐹! = 1 − 𝑒𝑥𝑝 (13)
!! !! !!!
 Cetina-Quiones et al. / Ingeniería 21-1 (2017) 1-12

siendo 𝐴! es el área del tubo receptor,𝑐! el calor El Método de Elementos Finitos

específico de agua y 𝐹 ! es el factor de eficiencia del El método de los elementos finitos (MEF) es un
colector calculado como: método numérico para la solución de problemas de
ingeniería comúnmente empleado en aquellos que
! involucran un alto grado de complejidad matemática,
! !!
𝐹 = ! ! ! ! (14) así como de físico-matemáticas, con aplicaciones en
! ! ! ! !" !
!! !! !! !!! !! problemas de transferencia de calor, flujo de fluidos,
transporte de masa, etc. (Baskharone, 2013).
de la cual, 𝑘! es la conductividad térmica del tubo
receptor. El MEF involucra diversos métodos numéricos que se
La ecuación que relaciona la ganancia de energía de emplean para obtener resultados aproximados a valor
un CCP es: real. En este trabajo se empleó el método de Galerkin
el cual es un método con un error de alrededor del 7%
𝑚𝑐! 𝑇! − 𝑇! = 𝐹! 𝜂! 𝐴! 𝐺! − 𝐴! 𝑈! 𝑇! − 𝑇! (15) del valor real, y con un alto grado de precisión
(Nithiarasu, 2016). El método de Galerkin se enuncia
siendo 𝑇! la temperatura inicial del fluido y 𝐴! el área como:
de apertura de la lámina reflectora. Por lo tanto, la
𝒙𝒋
temperatura de salida del fluido queda como:
𝑵𝒊,𝒋 𝒙 𝑹 𝒙: 𝑻 𝒅𝑽 = 𝟎 (18)
𝒙𝒊
!! !! !! !! !!! !! !! !!!
𝑇! = 𝑇! + (16)
!!!
Donde:
donde 𝑇! es la temperatura ambiente. Finalmente, la 𝑥! : es el primer nodo del sistema.
eficiencia térmica del CCP se halla mediante la 𝑥! : es la longitud nodal del sistema.
!
primera ley de la termodinámica como: 𝑁!,! 𝑥 : es la función de prueba nodal, 𝑁! = 1 −
!!
!
!! !! !!!
𝑁! = , de la cual 𝐿! es la longitud nodal.
!!
𝜂! = 𝐹! 𝜂! − (17)
!! !! 𝑅 𝑥: 𝑇 : es la función de residuo, la cual mide la
diferencia entre la solución de prueba y la solución
real, siendo ésta la Ec. 15.
𝑑𝑉: es la diferencial de volumen.

De la Ec. 15, 𝑇! − 𝑇! se reescribe como diferencial.

De modo que aplicando el método de Galerkin a
dicha ecuación queda:

𝑳𝒏 (19)
𝒅𝑻 𝒙
𝑵𝒊,𝒋 𝒙 𝒅𝑽 = 𝑭𝑹 𝜼𝒐 𝑮𝒃 𝑾𝒂 − 𝑼𝑳 𝝅𝑫𝒆 𝑻𝒊 − 𝑻𝒂 𝑵𝒊,𝒋 𝒙 𝒅𝑽
𝒅𝒙 𝟎
𝑳𝒏
𝒎𝒄𝒑
𝟎

donde 𝑊! es la longitud de apertura de la lámina. nodo j multiplicada por su longitud nodal hasta dicho
punto se expresa como 𝑇 𝑥 = 𝑁! 𝑥 𝑇! + 𝑁! 𝑥 𝑇! , la
Considerando que 𝑁!,! se descompone en 𝑁! y 𝑁! y Ec. 19 se reescribe en forma matricial como:
que la temperatura en el nodo i multiplicada por su
respectiva longitud nodal, más la temperatura en el

𝒅𝑵𝒊 𝒙 𝒅𝑵𝒋 𝒙
𝑳𝒏 𝑵𝒊 𝒙 𝑵𝒊 𝑻𝒊 𝑳𝒏
𝑵𝒊 𝒙
𝒅𝒙 𝒅𝒙 (20)
𝒎𝒄𝒑 𝑻𝒋 𝒅𝒙 = 𝑭𝑹 𝜼𝒐 𝑮𝒃 𝑾𝒂 − 𝑼𝑳 𝝅𝑫𝒆 𝑻𝒊 − 𝑻𝒂 𝒅𝒙
𝟎 𝒅𝑵𝒊 𝒙 𝒅𝑵𝒋 𝒙 𝟎
𝑵𝒋 𝒙
𝑵𝒋 𝒙 𝑵𝒋
𝒅𝒙 𝒅𝒙
 Cetina-Quiones et al. / Ingeniería 21-1 (2017) 1-12

Integrando y simplificando analíticamente los términos de la Ec. 20, queda finalmente como:

𝒎𝒄𝒑 −𝟏 𝟏 𝑻𝒊 𝑭𝑹 𝜼𝒐 𝑮𝒃 𝑾𝒂 − 𝑼𝑳 𝝅𝑫𝒆 𝑻𝒊 − 𝑻𝒂 𝑳𝒏 𝟏 (21)

𝑻 =
𝟐 −𝟏 𝟏 𝒋 𝟐 𝟏

La Ec. 21 representa la temperatura del fluido a lo receptor 𝐿! y después el número de nodos 𝑁! ,

largo del tubo receptor simplificada para 2 nodos. esto con el fin de calcular la distancia nodal 𝐿! .
Para obtener un perfil de temperaturas es necesario Posteriormente se ensambla la matriz generada y el
trabajar con una mayor cantidad de nodos. vector generado. Se incluyen las condiciones de
frontera para ambos y se calcula el perfil de
Desarrollo computacional temperaturas resolviendo la matriz ensamblada.
El MEF fue desarrollado en el entorno de
programación MATLAB. El diagrama de flujo
implementado para la programación se presenta en la
Fig. 3. Primero se inicializó la longitud del tubo
INICIO

Variables iniciales
𝐿 &,𝐺+ , 𝑉, 𝑄, 𝑇0 , 𝜑2 , 𝐷"

Inicializar número de nodos

𝑁"

Calcular distancia nodal

𝐿&
𝐿" =
𝑁" − 1

Discretizar el vector
Discretizar la matriz columna

𝐿 " = 𝐿 " + ∆𝐿"

∑𝐿 " = 𝐿 & no
∑𝐿 " = 𝐿 &
no
si si
Condición Ensamblar el vector
Ensamblar la matriz
de frontera columna

Calcular perfil de
temperaturas

FIN
Figura 3. Diagrama de flujo del MEF

Simulación térmica con el software SolidWorks
 Cetina-Quiones et al. / Ingeniería 21-1 (2017) 1-12

Para la validación del modelo del CCP se empleó el receptor, el cual permitió observar los fenómenos
software SolidWorks, mediante su herramienta presentados en éste. Para la validación se siguieron
FlowSimulation, pudo estimar el comportamiento de cuatro fases presentados en la Fig. 4.
la temperatura del fluido de trabajo a lo largo del tubo

Fase 2:
Inclusión de Fase 3: Fase 4:
Fase 1: Diseño parámetros Simulación Validación
nominales
Figura 4. Fases implementadas para la validación del modelo térmico del CCP

Fase 1: Se diseña el CCP en el software como se
observa en la Fig. 1. En la Fig. 5 se muestran los diferentes resultados
Fase 2: Se incluyen los parámetros de entrada con los obtenidos de los efectos de la variación de cada uno
que se realizará la simulación. Dichos parámetros se de los parámetros mencionados anteriormente sobre
explican en la sección de resultados. la temperatura del fluido. La Fig. 5(a), presenta el
Fase 3: Se realiza la simulación considerando los efecto de la irradiancia y se observa una mayor
parámetros incluidos en la fase 2, así como las ganancia térmica a irradiancias altas debido a que el
condiciones de frontera que intervienen en el flujo de calor hacia el tubo receptor es mayor en estas
comportamiento térmico de la temperatura del fluido. condiciones. De la misma manera, la Fig. 5(b)
Fase 4: En esta última fase se comparan los resultados presenta la variación de la velocidad del viento y se
obtenidos de la simulación del software con los observa un decremento de la ganancia térmica del
obtenidos por la modelación mediante el MEF, con lo fluido cuando este parámetro aumenta, debido a que
cual se valida el modelo desarrollado. ésta es menor a velocidades de viento altas por el
aumento de las pérdidas convectivas exteriores. La
RESULTADOS Y DISCUSIÓN Fig. 5(c) presenta la variación del caudal en la entrada
En esta sección se presentan los resultados obtenidos del tubo receptor y muestra una menor ganancia
de la modelación desarrollada empleando el MEF. Se térmica por parte del fluido cuando los caudales son
realizó un análisis térmico de los efectos de siete altos debido a que se tienen mayores pérdidas
parámetros incluidos en la modelación con el fin de convectivas en el interior del tubo receptor. Por otra
observar el comportamiento del fluido a lo largo de parte, en la Fig. 5(d) se refleja la variación de la
tubo receptor, así como el comportamiento de la temperatura inicial del fluido del tubo receptor del
eficiencia térmica del CCP. Finalmente se comparan CCP, de la cual se obtuvo un mayor gradiente de
los resultados y los errores obtenidos estadísticamente ganancia térmica cuando este parámetro aumenta, por
del modelo desarrollado a través del MEF, con la el hecho de tener un alto calor útil en la entrada del
simulación realizada empleando el software tubo receptor. La Fig. 5(e) presenta la variación del
SolidWorks a través de su herramienta ángulo de borde, en la cual se obtuvo la mayor
FlowSimulation. ganancia térmica para el ángulo de borde de 90°,
debido a que el aumento del área receptora de la
Efectos de los parámetros de entrada sobre la lámina produce un mayor factor de concentración y
temperatura del fluido por ende un mayor flujo de calor concentrado hacia el
Para la modelación se trabajó con siete parámetros tubo receptor. Finalmente, la Fig. 5(f) representa la
nominales los cuales se variaron para observar el variación del diámetro nominal y su efecto sobre la
comportamiento de la temperatura del agua a lo largo temperatura del fluido. El aumento de este parámetro
del tubo receptor, así como de la eficiencia térmica conlleva a obtener mayores pérdidas convectivas
del CCP. Estos parámetros fueron la irradiancia solar debido a que el área del tubo receptor es mayor. Sin
directa (𝐺! ) que se varió desde 500 hasta 1000 embargo, el factor de concentración es mayor cuando
𝑊 𝑚 ! , la velocidad del viento (𝑉) desde 1 hasta 6 el diámetro nominal disminuye. Por tal motivo, se
𝑚 𝑠, la temperatura ambiente (𝑇! ) de 25 a 35 °C, el debe encontrar un diámetro nominal el cual sea el
caudal (𝑄) de 1 hasta 6 𝑙 𝑚𝑖𝑛, la temperatura inicial óptimo para la minimización de las pérdidas
del fluido (𝑇! ) de 35 a 60 °C, el ángulo de borde del convectivas y maximice la ganancia térmica del
CCP (𝜑! ) 30 a 90° y el diámetro nominal del tubo fluido. Para este trabajo este diámetro fue de 1
receptor (𝐷! ), al cual se le asignaron los valores de ½, pulgada.
¾, 1, 1 ½ y 2 pulgadas respectivamente.
 Cetina-Quiones et al. / Ingeniería 21-1 (2017) 1-12

43 52
V= 1 m/s
42
Gb= 500 W/m2
(a) 50
V= 2 m/s
(b)
Gb= 600 W/m2

Temperatura del fluido [°C]

Temperatura del fluido [°C]

41 48 V= 3 m/s
Gb= 700 W/m2
V= 4 m/s
40 Gb= 800 W/m2 46 V= 5 m/s
Gb= 900 W/m2 44
V= 6 m/s
39
Gb= 1000 W/m2
42
38
40
37
38
36
36
35
0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0
Posición del fluido [m] Posición del fluido [m]
50 70
Ti=35 °C
Q=1 l/min
48 Q=2 l/min (c) 65 (d) Ti=40 °C
Ti=45 °C

Temperatura del fluido [°C]

Temperatura del fluido [°C]

Q=3 l/min Ti=50 °C

46 60
Q=4 l/min Ti=55 °C
44 Q=5 l/min 55 Ti=60 °C
Q=6 l/min
42 50

40 45

38 40

36 35

0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0
Posición del fluido [m] Posición del fluido [m]

42 42
Dn=1/2"
41
ϕr=30°
(e) 41
(f)
ϕr=40° Dn=3/4"
Temperatura del fluido [°C]

Temperatura del fluido [°C]

ϕr=50° Dn=1"
40 40
ϕr=60°
Dn=1 1/2"
ϕr=70°
39 39 Dn=2"
ϕr=80°
38 ϕr=90° 38

37 37

36 36

35 35
0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0
Posición del fluido [m] Posición del fluido [m]
Figura 1. Efectos de los parámetros nominales sobre la temperatura del fluido:
5(a) efecto de la irradiancia,
5(b) efecto de la velocidad del viento,
5(c) efecto del caudal,
5(d) efecto de la temperatura inicial de fluido,
5(e) efecto del ángulo de borde,
5(f) efecto del diámetro nominal.

En cuanto al comportamiento de la eficiencia, se valor designado a cada parámetro nominal, como se

incluyeron siete casos, a cada uno se le asignó un muestra en la Tabla 1.
 Cetina-Quiones et al. / Ingeniería 21-1 (2017) 1-12

Tabla 1. Casos para la modelación de la eficiencia térmica del CCP

Caso Irradiancia Velocidad Temperatura Caudal Temperatura Ángulo de Diámetro
(𝐺! ) del viento Ambiente (𝑇! ) (Q) inicial (𝑇! ) borde (𝜑! ) nominal
(𝑉) (𝐷! )

I 500 1 25 1 35 30 ½
II 600 2 27 2 40 40 ¾
III 700 3 29 3 45 50 1
IV 800 4 31 4 50 60 1½
V 900 5 33 5 55 70 2
VI 1000 6 35 6 60 80 -
VII - - - - 90 -

En la Fig. 6 se presentan las gráficas de la eficiencia temperatura inicial son inversamente proporcionales,
térmica 𝜂! para cada caso de la Tabla 1. Se observa lo que significa que el aumento de dichos parámetros
que los parámetros tales como la irradiancia, el produce menor ganancia térmica del flujo. De igual
caudal, la temperatura ambiente y el ángulo de borde, manera, se observa que en el comportamiento de la
son directamente proporcionales a la eficiencia eficiencia para la variación del diámetro nominal se
térmica. Dichos parámetros influyen del mismo modo obtuvo un diámetro óptimo de 1 pulgada. Para
que la temperatura del fluido a lo largo de tubo diámetros mayores a éste, la eficiencia térmica
receptor presentados en la Fig. 5, es decir, el aumento decrece. Finalmente, la eficiencia térmica más alta se
de cualquiera de estos parámetros involucrados obtuvo en la variación de la temperatura ambiente, a
genera mayor ganancia térmica por parte del fluido de los 35 °C alcanzando un valor del 72 %.
trabajo. Por otro lado, la velocidad del viento y la

0.75

0.70

0.65
Eficiencia térmica (-)

0.60

0.55
Irradiancia
0.50
Velocidad del viento
Caudal
0.45 Temperatura ambiente
Temperatura inicial
0.40 Ángulo de borde
Diámetro
0.35
I II III IV V VI VII
Casos
Figura 2. Comportamiento de la eficiencia térmica para cada caso

 Cetina-Quiones et al. / Ingeniería 21-1 (2017) 1-12

Validación del modelo térmico empleados en la modelación térmica con el MEF.

Para la validación del modelo se empleó el software Esta temperatura concuerda con la temperatura
SolidWorks como herramienta de simulación. En esta obtenida en SolidWorks debido a que esa irradiancia
sección, se realizaron varias simulaciones fue uno de los parámetros nominales en la
considerando las mismas variables de entrada que se simulación, manteniendo los demás parámetros
consideraron en la simulación desarrollada con el invariantes, tanto en la modelación como en la
MEF del modelo térmico del CCP, como ejemplo de simulación. Los tiempos de cálculo correspondientes
esta simulación con el software SolidWorks. la Fig. a la simulación con el software SolidWorks fueron
7(a) muestra la temperatura en la salida del tubo aproximadamente de 45 a 58 min. Por otra parte, en
receptor como una gráfica de corte transversal. En la Fig. 7(b) se presenta la gráfica de corte del tubo
dicha gráfica se puede observar que la sección roja en receptor que representa el comportamiento de la
la salida corresponde a una temperatura de temperatura a lo largo de éste. Se observa el cambio
aproximadamente 42 °C, por otro parte, de la Fig. de dicha temperatura del fluido conforme transcurre a
5(a), la temperatura aproximada que corresponde a la lo largo de tubo receptor, en donde el color rojo
irradiancia de 900 W/m2 para el tubo no evacuado es representa la mayor temperatura.
de 41.6 °C considerando los mismos parámetros

(b)

Lamina reflectora
Tubo receptor

Dirección del flujo

Figura 3. Gráfica de la distribución térmica de la temperatura del fluido:

(a) salida del tubo receptor,
(b) a lo largo del tubo receptor

Finalmente, en la Tabla 2 se muestran los valores obtenido fue de 7.6 % para el mínimo caudal
obtenidos de la modelación y de la simulación, con empleado, y el error mínimo fue del 0.5 % para el
sus respectivos errores, para los cuales el caudal fue caudal más alto. Por lo tanto, los resultados del
el parámetro a variar, debido a que éste es uno de los modelo desarrollado se aproximan a la simulación de
que generan mayores gradientes en la temperatura de SolidWorks y puede predecir eficazmente el
salida del fluido, como se puede observar en la Fig. comportamiento de la temperatura del fluido en
5(c). De igual manera se observa el error máximo caudales altos.

Tabla 2. Comparación de los resultados obtenidos por la modelación y la simulación para cada variación del caudal
Caudal MEF SolidWorks Error
1 l/min 53.70 °C 49.61 °C 7.6 %
2 l/min 44.55 °C 42.41 °C 4.8 %
3 l/min 41.65 °C 41.33 °C 0.7 %
4 l/min 40.06 °C 39.74 °C 0.8 %
5 l/min 39.07 °C 38.82 °C 0.6 %
6 l/min 38.40 °C 38.20 °C 0.5 %

CONCLUSIONES las ecuaciones del modelo reportadas se empleó el

En el presente trabajo se desarrolló un modelo físico Método de Elementos Finitos (MEF) como
para la transferencia de calor en un Concentrador herramienta de modelación. El modelo físico
Solar de Canal Parabólico (CCP). En la solución de desarrollado es capaz de calcular la distribución de
 Cetina-Quiones et al. / Ingeniería 21-1 (2017) 1-12

temperaturas del fluido de trabajo a lo largo del tubo desarrollado reporta un error máximo del 7.6% entre
receptor del CCP. De igual manera, el modelo el modelo teórico y el software SolidWorks. El
permite la estimación de la eficiencia del CCP modelo desarrollado en este trabajo es una
mediante un análisis térmico basado en la primera ley herramienta útil para predecir el comportamiento de
de la termodinámica. La distribución de temperaturas la temperatura del fluido de trabajo bajo diferentes
a lo largo del tubo receptor del CCP se obtuvo en condiciones ambientales y de operación. La ventaja
función de siete parámetros de entrada (irradiancia de este modelo, en comparación con el software
solar, velocidad del viento, temperatura ambiente, SolidWorks, son los cortos tiempos de cálculo,
caudal, temperatura inicial del fluido, ángulo de borde debido a que, con el software, este tiempo es de
del CCP y diámetro nominal del tubo receptor), los aproximadamente 45 a 58 minutos, mientras que
cuales se variaron para estudiar el comportamiento de empleando MATLAB, dicho cálculo no supera los 30
la temperatura del fluido en función de dichos segundos. Como trabajos a futuro se desea comparar
parámetros. Los resultados obtenidos en la los resultados obtenidos de la modelación térmica con
modelación con MEF se validaron mediante la resultados experimentales del CCP, de tal manera que
simulación del fluido empleando el software pueda validarse el modelo desarrollado y pueda
SolidWorks a través de su herramienta emplearse en el dimensionamiento y la
FlowSimulation. La validación del modelo caracterización de los CCP.

REFERENCIAS
Baskharone, E. A. (2013). The finite element method with heat transfer and fluid mechanics applications. Cambridge
University Press.

Cengel,Y.A. y Ghajar,A.(2011). Transferencia de calor y masa. Editorial: Editorial McGraw–Hill.

Duffie, J. A. y Beckman, W. A. (2013). Solar engineering of thermal processes (Vol. 3). New York:Wiley.

Edenburn, M. W. (1976). Performance analysis of a cylindrical parabolic focusing collector and comparison with
experimental results. Solar Energy, 18(5), 437-444.

Eiamsa-Ard, S., Thianpong, C., y Eiamsa-Ard, P. (2010). Turbulent heat transfer enhancement by counter/co-
swirling flow in a tube fitted with twin twisted tapes. Experimental Thermal and Fluid Science, 34(1), 53-62.

Fernández-García, A., Rojas, E., Pérez, M., Silva, R., Hernández-Escobedo, Q., y Manzano-Agugliaro, F. (2015). A
parabolic-trough collector for cleaner industrial process heat. Journal of Cleaner Production, 89, 272-285.

Fernandez-Garcia, A., Zarza, E., Valenzuela, L., y Pérez, M. (2010). Parabolic-trough solar collectors and their
applications. Renewable and Sustainable Energy Reviews, 14(7), 1695-1721.

Gnielinski, V. (2013). On heat transfer in tubes. International Journal of Heat and Mass Transfer, 63, 134-140.

Jaramillo, O. A., Borunda, M., Velazquez-Lucho, K. M., y Robles, M. (2016). Parabolic trough solar collector for
low enthalpy processes: An analysis of the efficiency enhancement by using twisted tape inserts. Renewable
Energy, 93, 125-141.

Kalogirou, S. A. (2013). Solar energy engineering: processes and systems. Academic Press.

Kalogirou, S. A. (2004). Solar thermal collectors and applications. Progress in energy and combustion science, 30(3),
231-295.
Liang, H., You, S., y Zhang, H. (2015). Comparison of different heat transfer models for parabolic trough solar
collectors. Applied Energy, 148, 105-114.

May, O., Ricalde, L. J., Ali, B., López, E. O., Venegas-Reyes, E., y Jaramillo, O. A. (2016). Neural Network Inverse
Modeling for Optimization. In Artificial Neural Networks-Models and Applications. InTech.

May Tzuc, O., Bassam, A., Escalante Soberanis, M. A., Venegas-Reyes, E., Jaramillo, O. A., Ricalde, L. J., ... y El
Hamzaoui, Y. (2017). Modeling and optimization of a solar parabolic trough concentrator system using inverse
 Cetina-Quiones et al. / Ingeniería 21-1 (2017) 1-12

artificial neural network. Journal of Renewable and Sustainable Energy, 9(1), 013701.

Nithiarasu, P., Lewis, R. W., y Seetharamu, K. N. (2016). Fundamentals of the finite element method for heat and
mass transfer. John Wiley & Sons.

Nwosu, N. P. (2009). Finite-element analysis of an absorber in an evacuated solar tube heat exchanger employing the
Galerkin method. International Journal of Sustainable Energy, 28(4), 247-255.

Padilla, R. V., Demirkaya, G., Goswami, D. Y., Stefanakos, E., y Rahman, M. M. (2011). Heat transfer analysis of
parabolic trough solar receiver. Applied Energy, 88(12), 5097-5110.

Tzivanidis, C., Bellos, E., Korres, D., Antonopoulos, K. A., y Mitsopoulos, G. (2015). Thermal and optical
efficiency investigation of a parabolic trough collector. Case Studies in Thermal Engineering, 6, 226-237.

Uzgoren, E. (2016). One-dimensional transient thermal model for parabolic trough collectors using closed-form
solution of fluid flow. International Exergy, Energy ans Environment Symposium (IEEES-8).

_____________________________

Este documento debe citarse como: Cetina-Quiñones, A. J., Bassam, A., Hernández-Chan, G., Hernández Benítez, J.,
Hernández Reyes, I., Lugo Chávez, D. (2017). Modelación térmica de un colector solar de canal parabólico
mediante el método de elementos finitos. Ingeniería, Revista Académica de la FI-UADY, 21-1, pp. 1-12, ISSN
2448-8364.