Universidad Nacional Del Santa: Facultad de Ingeniería
Universidad Nacional Del Santa: Facultad de Ingeniería
FACULTAD DE INGENIERÍA
INGENIERO EN ENERGÍA
AUTOR :
NUEVO CHIMBOTE-PERU
2018
UNIVERSIDAD NACIONAL DEL SANTA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICA PROFESIONAL DE INGENIERÍA EN ENERGIA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICA PROFESIONAL DE INGENIERÍA EN ENERGIA
A mis hermanos:
Elmer, Jaimito, Mayreli, Isaura, Anali.
Quienes me brindaron todo su apoyo
para progresar en este tiempo.
A todas las personas que contribuyeron en mi formación, por que sin ellos hubiera sido
imposible llegar hasta donde estoy, ya que su influencia fue necesaria para llegar mis
objetivos.
El autor
ÍNDICE
DEDICATORIA
AGRADECIMIENTO
RESUMEN
ABSTRACT
CAPITULO I
I. INTRODUCCIÓN 1
1.1. Antecedentes 2
1.2. Formulación del problema 5
1.3. Objetivos 6
1.3.1. Objetivo General 6
CAPITULO III 31
III. MATERIALES Y MÉTODO 32
3.2. Método 32
3.2.1. Determinación de la radiación efectiva en Nuevo Chimbote 32
CAPITULO IV 65
IV. RESULTADOS Y DISCUSIÓN 66
4.1. Irradiación efectiva orientación del eje este-oeste y norte- sur 66
4.2. Dimensionamiento del concentrador solar cilíndrico parabólico 68
4.3. Características generales en el diseño del concentrador solar CSCP 69
4.4. Temperatura logarítmica en el intercambiador de calor (ºC) 71
V. CONCLUSIONES Y RECOMENDACIONES 72
5.1. Conclusiones 72
5.2. Recomendaciones 73
VI. REFERENCIAS BIBLIOGRÁFICAS Y VIRTUALES 74
VII. ANEXOS 77
INDICE DE FIGURAS
Figura 1. Longitud de onda 11
Figura 2. Tres clases de radiación solar. 12
Figura 3. Colector plano. 14
Figura 4. Colectores Cilíndricos Parabólicos. 15
Figura 5. Principio de funcionamiento y componentes de un CCP. 18
Figura 6. Esquema del seguimiento solar a un eje de un CCP. 20
Figura 7. Intercambiador de calor a contra flujo. 23
Figura 8. Balance energético de un CCP. 29
Figura 9. Espectro de Radiación Solar. 33
Figura 10. Relaciones angulares terrestres. 34
Figura 11. Esquema solar frente al vector normal a la superficie 35
Figura 12. Mapa del territorio de la loza sosal de la Universidad Nacional del Santa 37
Figura 13. Perfil de la superficie captadora 44
Figura 14. Esquema de relaciones y distancias dentro del perfil del CSCP 47
Figura 15. Proyección del disco solar sobre el absorbedor 49
Figura 16. Método volumétrico para hallar el flujo másico 51
Figura 17. Medición del gradiente de temperatura 52
Figura 18. Perdidas por conducción, convección y radiación en la tubería 54
Figura 19. Intercambiador de calor a contraflujo 60
Figura 20. Medición de temperatura en el intercambiador de calor 61
Figura 21. Sistema de seguimiento solar 62
Figura 22. Concentrador solar orientado hacia el este 63
Figura 23. CS orientado perpendicularmente a la posición del sol (12:00 horas) 63
Figura 24. Concentrador solar orientado hacia el oeste 64
Figura 25. CSCP acoplado al intercambiador de calor a contra flujo 71
INDICE DE GRAFICAS
INDICE DE TABLAS
Dicho trabajo se realizó teniendo en cuenta las condiciones de radiación solar, y otros
parámetros meteorológicos de la ciudad de Nuevo Chimbote, obteniendo 1128.69 W de
energía incidente, 1.712 𝑚2 de área eficiente de la superficie reflectante, 0.55 de eficiencia
óptica y 71 % de rendimiento global del concentrador solar cilíndrico parabólico.
El concentrador solar cilíndrico parabólico permitió obtener la energía térmica necesaria (50-
70°C) para poder suplir las necesidades energéticas requeridas por el intercambiador de calor
a contra flujo obteniendo un gradiente de temperatura de 44.3°C valor requerido y aceptable,
y así poder realizar las prácticas de laboratorio del curso de Transferencia de Calor y Masa de
la Escuela Profesional de Ingeniería en Energía de la Universidad Nacional del Santa.
The aim of this research work was to build and implement a cylindrical parabolic solar
concentrator for the preheating of water required in the counter-flow heat exchanger of the
thermo-fluid laboratory of the National University of Santa.
The research work was done take into account the conditions of solar radiation, and other
meteorological parameters of the city of Nuevo Chimbote, obtaining 1128.39 W of incident
energy, 1,712 𝑚2 of efficient area of the reflecting surface, 0.55 of optical efficiency and 71%
of global performance of the cylindrical parabolic solar concentrator.
The cylindrical parabolic solar concentrator made it possible to obtain the necessary thermal
energy (50-70 °C) to be able to supply the energy needs required by the heat exchanger in the
counter-flow, obtaining a temperature gradient of 44.3 °C, a value required and acceptable by
it, and so to be able to made the laboratory practices of the court of Heat and Mass Transfer
the professional career of Engineering in Energy at National University of Santa.
Es claro que las fuentes de energías limpias poseen un elevado potencial económico, siendo
en lo más mínimo el impacto con el medio ambiente, sabiendo que hasta la actualidad aún
no se han explotado ni comercializado a gran escala. Lo anterior se da por el elevado costo
de inversión en desarrollar estas tecnologías, puesto que, dichos sistemas no han tenido
mucha aceptación debido a no ser muy continuo y siendo necesario el almacenamiento de
energía, teniendo como consecuencia optar por las fuentes convencionales, sin embargo, en
las tecnologías con fuentes de energías limpias, teniendo un bajo impacto ambiental, lo que
se debe considerar es el costo de operación. De esta forma, el coste de los sistemas con
fuentes de energías limpias se da por el ahorro en el consumo de combustible. Puesto que la
energía renovable por si es gratuita, siendo el desafío diseñar métodos de captación y
aprovechamiento que conviertan dicha energía limpia comparada con las fuentes de
combustibles fósiles, puesto que se quiere tener sistemas con muy alta eficiencia de
conversión. Teniendo otro problema en la discontinuidad, siendo por el ciclo día – noche ó
también por las condiciones atmosféricas, tales como la lluvia o días nublados.
1
El siguiente trabajo de investigación tiene como objetivo de construir e implementar un
concentrador solar cilíndrico parabólico (CSCP) para el pre-calentamiento de agua requerida
`por un intercambiador de calor a contra flujo, es decir generar un gradiente de temperatura
de 40ºC.
1.1. Antecedentes
2
“Evaluación Experimental de un Concentrador Solar Cilíndrico Parabólico”,
(UNIVERSIDAD NACIONAL JORGE BASADRE GROHMANN DE TACNA –
Perú).
Realizado por CARLOS POLO BRAVO, ELISBAN J. SACARI SACARI Y JORGE
ARMANDO CHOQUE CHACOLLA, el 18 de octubre del 2011.
El colector solar cilíndrico parabólico (CCP), fue diseñado y construido por estudiantes
de la Especialidad de Energías Renovables de la Escuela Académica Profesional de Física
Aplicada.
Dicho trabajo se evaluó teniendo en cuenta las condiciones meteorológicas de radiación
solar de la ciudad de Tacna, y conlleva a lo siguiente: definir la superficie reflectante más
eficiente, perdidas de calor según el ángulo de inclinación del colector, determinar la
eficiencia óptica, el coeficiente global de transferencia de calor y la eficiencia térmica.
3
movimiento gradual de los dos ejes del concentrador solar, haciendo uso de una
programación, el cual se realizó en el programa Labview, también se hizo uso de una
tarjeta Microcontroladora Arduino y una verificación de la posición proveniente de dos
sensores adaptados mediante instrumentos diseñados para fijación de estos a los ejes del
concentrador solar.
El mando de motores se diseña en forma de un subprograma, el cual es usado por el
programa principal de mando, entretanto en la etapa de potencia de los motores se realiza
haciendo uso de una segunda tarjeta Microcontroladora el cual se usa al mismo tiempo
para la obtención de datos y monitoreo de esta manera las variables del sistema el cual
contienen un HMI diseñada y esquematizada en el entorno Labview.
4
1.2. Formulación del problema
Actualmente la Escuela Académica Profesional de Ingeniería en Energía (EAPIE), tiene
incluido en su currícula el curso de TRANFERENCIA DE CALOR Y MASA. El
Laboratorio de Termo-Fluidos de la EAPIE cuenta con un intercambiador de calor a
contra flujo, donde los alumnos realizan sus prácticas de laboratorio, para poder analizar
y conocer el fenómeno de transferencia de calor. Para realizar dicha práctica y observar
los fenómenos que se dan en el intercambiador de calor, se requiere de agua caliente
aproximadamente 60 °C y agua fría a T ambiente (20 °C).
Para lograr las condiciones de temperatura requerida por el intercambiador de calor, los
estudiantes hacen uso de los recursos disponibles en dicho laboratorio, en la actualidad
únicamente se cuenta con una cocina a gas para calentar el agua, siendo este un método
eficaz, pero no eficiente para esta práctica, puesto que demanda tiempo, coste, una
combustión que emite CO2 y contamina el medio ambiente.
5
1.3. Objetivos
7
CAPITULO II
8
II. MARCO TEORICO
2.1. Energía solar
La energía solar es una energía limpia, que se obtiene a partir de la radiación
electromagnética procedente del Sol. La radiación solar que incide sobre la Tierra ha
sido empleada por el ser humano desde la Antigüedad, haciendo uso de diferentes
tecnologías que han ido evolucionando en el tiempo (Vásquez & Zuñiga, 2015). En la
actualidad, el calor y la luz del Sol puede usar por medio de múltiples captadores como
células fotovoltaicas, helióstatos o concentradores térmicos, logrando transformarse en
energía térmica o eléctrica. Es conocida con fuentes de energías renovables o energías
limpias, que podrían apoyar en resolver muchos de los problemas más inminente que
enfrenta la humanidad.
Abril (2010), indica que la energía solar térmica es la energía que transforma la radiación
solar en calor, y cede este calor a un líquido, comunmente agua. Se puede aprovechar el
calor del líquido calentado para conseguir agua caliente sanitaria, calefacción de
viviendas, calefacción en las piscinas, etc. El propósito prinicipal de cualquier sistema
térmico solar es lograr la energía proveniente del sol para convertirla y transportarla de
la forma más económica y eficaz posible. La utilización del aprovechamiento solar
pueden ser directas, en modo de calor, o bien indirectas, empleando este calor para hallar
trabajo mecánico en un eje y finalmente energía eléctrica (Paredes, 2012). A partir de
1973, la energía solar ha sido objeto de diversas investigaciones, y sin obligación de
técnicas complicadas, cada vez se conoce más estudios que son sumamente prácticas y
atractivas en el calentamiento de agua.
9
El sol a través de reacciones nucleares que se dan en su interior genera su propia energía.
En estas reacciones nucleares los átomos de hidrogeno (elemento abundante del sol), se
combinan entre si y se obtiene átomos de hielo. Es importante mencionar que una
pequeña parte de la masa de los átomos de hidrogeno, se transforman en energía según
la famosa ecuación de Einstein, E = m.𝑐 2 , donde E = cantidad de energía liberada cuando
desaparece la masa m; C = velocidad de la luz. Esta energía fluye desde el interior hasta
la superficie (fotosfera) y es irradiada cada segundo una energía de 4x1020 J en todas las
direcciones al espacio, lo que genera una potencia de 4x1023 kW (Salgado, 2013). En
un solo segundo, el sol irradia más energía que la consumida por todo el género humano
desde sus albores hasta nuestros días.
Constante solar
Sobre la atmosfera terrestre, llega una cantidad de radiación solar constante. Debido a la
radiación que emite el sol y sus condiciones geométricas que presenta respecto de la
tierra. Lo explicado define lo denominado constante solar. La constante solar (I cs), es la
cantidad de energía que recibe una superficie de área unitaria que se encuentra
perpendicular a la radiación, esta energía es emitida por el sol en unidad de tiempo.
(Castañeda & Hernández, 2013). El valor normal o estándar que tiene la constante solar
(Ics) propuesto por Thekaekara & Drumnond es de 1.37 kW/m2 h ó 428 BTU/ft2 . Este
valor ha variado en los últimos años, por las técnicas de medición empleadas, lo cual no
implica que haya variado la magnitud de la energía del Sol.
La constante solar, está distribuida según distancias de onda. Sin embargo, la radiación
solar adquiere solamente valores significativos de longitud de onda, entre 0.2 y 0.4 µm
como se muestra en la figura 1.
10
Figura 1. Longitud de onda. Fuente: Castañeda &
Hernández (2013)
11
- Radiación directa: Es aquella radiación solar que incide en la
superficie de la tierra sin tener cambios en su trayectoria lineal desde
el disco solar. Esta radiación está caracterizada por proyectar una
sombra definida de los objetos opacos que la interceptan.
12
Conocer la radiación solar directa y difusa es necesario para el análisis y diseño de
algunos sistemas solares. Por ejemplo, el valor de la radiación directa normal es
necesario para el cálculo de la radiación solar sobre superficies inclinadas, así como para
establecer el desempeño de colectores concentradores. (Solartronic, 2003).
13
2.5.1. En función de la temperatura del fluido
En esta definición se basa en las temperaturas que puede alcanzar el fluido, por
consiguiente, las distintas aplicaciones que pueden aprovecharse. También, el rango de
operación del campo solar está en función del tipo de fluido que utilice en dicho
concentrador solar, siendo de mucha influencia para la obtención del rendimiento del
ciclo de potencia. (Paredes, 2012).
14
la producción de energía eléctrica. Asimismo, se conocen como colectores de
concentración. Estando incluido los Concentradores Cilíndricos Parabólicos,
siendo la mitad de un cilindro el área que recibe toda la energía del sol. (Fig. 4.),
se hace uso de la energía solar directa y siempre se debe realizar una orientación
hacia el sol, haciendo uso de un mecanismo apropiado.
Á𝑟𝑒𝑎 𝑑𝑒 𝑎𝑝𝑒𝑟𝑡𝑢𝑟𝑎
𝐶=
Á𝑟𝑒𝑎 𝑑𝑒 𝑎𝑏𝑠𝑜𝑟𝑐𝑖ó𝑛
15
2.5.2.1.Colectores no concentrados.
Los colectores solares sin concentración (C=1), se usan para temperaturas bajas.
Debido a que las superficies de apertura y absorción son iguales, lo que significa
que todo el elemento del subsistema solar tiene el mismo emplazamiento físico.
Este tipo de concentradores son ampliamente utilizados en aplicaciones térmicas
domésticas, pues su costo y rango de temperatura permite una buena adaptación a
las necesidades del usuario (Castañeda & Hernández, 2013).
16
Tabla 1. Clasificación en función del factor de concentración. Fuente: Paredes.
Los sistemas cilindro parabólicos Este tipo de sistemas permiten obtener razones
de concentración en un rango de 30-90. Este sistema cuenta con espejos cilíndricos
cuya sección transversal es una parábola, permitiendo que la radiación solar se
centre en el eje central focal.
Los sistemas de torre central este tipo de sistemas permiten obtener factores de
concentración entre 200 y 1000. Están compuesto por espejos cuasi-planos
(helióstatos), que se distribuyen en una superficie horizontal, inclinados de tal
forma que la radiación solar se refleje hacia la parte superior de la torre, donde
generalmente se coloca el receptor.
17
Los discos parabólicos son sistemas que consiguen razones de concentración
altas, entre 1000 y 5000. Consisten en espejos paraboloide de revolución, que giran
siempre orientados al sol.
La concentración óptica de radiación solar implica que la superficie del tubo absorbedor
sea menor al área de abertura del captador, permitiendo de esta forma que se reduzca
significativamente las pérdidas térmicas del colector, debido a que estas pérdidas están en
función de la superficie del absorbedor, que es el componente “caliente” del captador.
(Echeverría, 2011). Debido a que solo la radiación solar directa es aprovechable por el
captador, es necesario que éste gire a lo largo del día en dirección al sol.
18
2.6.1. Componentes
Los elementos principales de un CCP son (Fig. 5): el reflector cilindro parabólico, tubo
absorbedor o receptor, sistema de seguimiento solar y fluido de trabajo
19
2.6.1.3. Sistema de Seguimiento Solar
Los CCP necesitan un sistema de seguimiento solar para poder seguir al sol durante
el dia, esto exige que el sistema permita cambiar de posicion al CCP conforme el
sol se mueva, y de esta forma concentrar la radiacion solar sobre el tubo receptor.
Los sistemas de seguimiento solar consisten en un dispositivo que permita girar
los reflectores cilindricos alrededor de un eje (Fig. 6). (Zarza, 2002). La
orientación puede ser Norte-Sur o Este-Oeste.(Paredes, 2012).
Existen CCP que giran en dos ejes, sin embargo se ha demostrado que este tipo de
captadores son poco rentables por que demandan de mayor mantenimiento (por su
complejidad mecánica) a pesar de que el flujo útil de radiación solar sobre el área
de abertura sea mayor a diferencia de los CCP con sistemas de seguimiento de un
solo eje, las perdidas térmicas son mayores debido a que la longitud de tuberías
pasivas dentro del captador es mayor, por tanto la potencia térmica útil que
suministra es igual o incluso menor a los que suministra un CCP con sistemas de
seguimiento en un solo eje.
20
2.6.1.4. Fluidos de Trabajo para los CSCP
Los CCP usan fluidos de transferencia de calor que, circulan en el tubo receptor, y
absorben en forma de energía térmica la radiación solar, transportándola hasta el
bloque de potencia. La función de los fluidos de trabajo es transportar el calor
producido a tanques de almacenamiento térmicos, que posteriormente son usados
en sistemas de calentamiento o enfriamiento. Los fluidos de trabajo pueden ser:
Agua, glicolpropileno, glicol-etileno, sales fundidas, aceites térmicos y aire
(Velázquez, 2014). La elección del tipo de fluido depende de las temperaturas que
se desee obtener. Por ejemplo, si se desea obtener temperaturas mayores a 175 °C
el agua desmineralizada como fluido de trabajo es ideal, ya que la presión de
trabajo no es excesiva. Por otro lado, la aplicación de aceites sintéticos como fluido
de trabajo permite obtener temperaturas más altas (175 ºC< T <400 ºC).
21
2.8. Intercambiador de calor
Los intercambiadores de calor son dispositivos que estan diseñados de forma que
permiten transferir calor entre fluidos, que esten separados por una barrera o en
contacto (Cuadrado, 2010). Según Palomino (2010), el proceso de transferencia de
calor se realiza en fluidos que presentan diferentes temperatruras, pues el objetivo de
esta transeferencia de calor es variar dichas temperaturas hasta alcanzar un equilibrio
termico. Esta transferencia comprende convección en cada fluido y conducción
mediante la pared que los separa. Los intercambiadores de calor según Cengel (2004)
son utilizados para: Calentar fluidos frios mediante fluidos que tienen mayor
tempertura, reducir la temperatura de un fluido a traves de fluidos con temperturas
bajas, llevar al punto de ebullicón a un fluido, asi mismo condensar un fluido que se
encuentre en estado gaseoso mediante un fluido frío.
Flujo paralelo
En este tipo de intercambiadores el flujo interno y externo fluyen en la misma
dirección. Estos fluidos ingresan al intercambiador por el mismo extremo y tienen una
diferencia de temperatura significativa. La temperatura de dichos fluidos se
aproximan, un fluido disminuye de temperatura y el otro aumenta, hasta alcanzar el
equilibrio térmico.
Flujo cruzado
En este caso uno de los fluidos circula de manera perpendicular al otro, un fluido circula
a traves de los tubos internos mientras el otro fluido pasa alrededor de estos. Estos
intercambiadores son usados comunmente cuando un fluido presenta un cambio de
fase.
22
Contraflujo
En los intercambiadores de calor a contraflujo los fluidos van hacia una misma
dirección, pero uno de estos fluidos tiene sentido contrario (Fig. 7). El ingreso de los
fluidos se da por diferentes extremos. Esto permite que al salir un fluido de temperatura
menor por el extremo donde entra el fluido de mayor temperatura, la tempratura del
primer fluido se aproxime a la temperatura del fluido de entrada. Este tipo de
intercambidores son considerados como los mas eficientes.
q = U.A.ΔTml (2.1)
Siendo: q = calor total (W), U = coeficiente global de transferencia de calor (W/𝑚2 . °𝐶),
A = área de transferencia de calor (𝑚2 ) y ΔTml = temperatura media logarítmica (°C).
𝛥𝑇2− 𝛥𝑇1
ΔTml = 𝛥𝑇1 (2.2)
ln(𝛥𝑇2)
23
2.9. Diseño de un concentrador cilíndrico parabólico
2.9.1. Características generales
Un CSCP consiste de un arreglo de espejos, formados de lámina de acero inoxidable con
acabado espejo, diseñado en forma de cilindro que presenta generatriz parabólica. Este
mecanismo tiene un solo grado de libertad, con orientación Este- Oeste y un sistema de
seguimiento Norte – Sur. Estos CCP tienen un tubo de cobre en la línea focal, por donde
pasa el fluido termo portador, que absorbe la radiación solar y posteriormente la
transforma en energía calorífica permitiendo calentar agua a altas temperaturas.
𝐴𝑐
𝐶= (2.3)
𝐴𝑎
El ángulo máximo que pueden formar dos rayos en un plano transversal a la apertura de
un concentrador se denomina ángulo de aceptancia, de tal forma que dichos rayos
reflejados interceptan en el tubo receptor en su reflexión. El valor de la razón de
concentración se encuentra alrededor de 220 teóricamente, sin embargo, en la práctica es
20. El ángulo de aceptancia de los CCP actuales son inferiores a 1°, esto se toma como
referencia para el error máximo que se puede presentar en el sistema de seguimiento solar.
Por otro lado, la temperatura ideal de trabajo se encuentra dentro del rango 150-400°C
(Paredes, 2012).
24
Pérdidas ópticas: Estas pérdidas de radiación se debe a que los elementos del
concentrador no son perfectos. Esto provoca que solo una cantidad de la radiación
solar directa llegue a la apertura del concentrador y alcance al fluido que se
encuentra en el interior del tubo receptor (Balado, 2011). Los cuatro parámetros
que afectan a las pérdidas ópticas son:
25
Estos valores permanecen constantes, siempre y cuando no se degraden los
materiales del CSCP, debido a que no depende ni la radiación solar ni de la
temperatura del fluido. Sin embargo, las partículas del polvo o suciedad que se
encuentren sobre las superficies del tubo receptor y de los espejos del CSCP,
producto de tener una instalación al aire libre, si influye en el rendimiento óptico.
Por tal motivo lo habitual es proporcionar los parámetros ópticos con un grado del
100% de limpieza (Paredes, 2012).
Donde:
Para un CSCP con tubo receptor de vacío es de unos 4 (W/𝑚2 . °𝐶) para
temperaturas de entorno a los 350°C.
26
Perdidas geométricas: Estas pérdidas provocan la disminución del área efectiva
de captación solar. Un CSCP tiene pérdidas debido a su posición relativa de los
concentradores entre sí, las pérdidas inherentes a cada concentrador son motivo
del ángulo de incidencia, lo que provoca que en los extremos del concentrador
exista una pérdida de superficie reflexiva útil. El área que el concentrador pierde
se define de la siguiente manera:
Fm = F + (F . 𝐿𝑟 2/48 . 𝐹 2 )
Siendo:
27
Rendimiento óptico, Este rendimiento consiste en considerar para un ángulo de
incidencia nulo, los efectos de la reflectividad de los espejos, la absortividad del
tubo receptor, la transmisividad de la cubierta de vidrio y el factor de intercepción
del concentrador. Donde se debe tener en cuenta el grado de ensuciamiento de la
instalación a través del llamado factor de ensuciamiento (o limpieza).
Rendimiento térmico, ήth este rendimiento considera todas las pérdidas que
existe en el CSCP.
Rendimiento global, ήglobal Considera todas las pérdidas. Que puede definirse
como el cociente entre la potencia térmica útil que entrega al fluido térmico, y el
recurso solar disponible en ese instante, es decir, la potencia Qsol.
Donde:
28
Siendo:
El cociente entre la energía térmica útil suministrada por el CSCP y la energía solar
que incide sobre el concentrador se denomina rendimiento global (Fig. 8).
Qútil
ήglobal = Qsol (2.9)
29
Para un grado de limpieza del 100% habitualmente el valor del rendimiento óptico
pico es 0.75, K, depende directamente del ángulo de incidencia, siendo K = 1 para
φ = 0°, y K = 0 para φ = 90°. El valor de K es una función K=K (φ) que se determina
experimentalmente. La radiación solar directa y la temperatura de trabajo del
fluido influyen en el rendimiento térmico (Balado, 2011).
30
CAPITULO III
31
III. MATERIALES Y MÉTODO
3.1. Materiales
- 01 unid. Plancha de aluminio anodizado de (1.5 m de largo/1.2 m de ancho).
- 1.5 m de tubo de cobre de ½”.
- 01unid. Cubierta de vidrio (1.5 m de largo/1.2 m de ancho).
- 10 m. Manguera de ½”.
- 6 unid. Niples de ½”x 8mm.
- 10 und de acoples.
- 01 gal. De pegamento de PVC.
- 01 gal. Pintura negra para cobre.
- 1 und. Tanque Rotoplas de 750 Lts.
- 10 und. Válvulas de ½”.
- 15 und. De codo de 90 °.
- 10 und. Termocúplas.
3.2. Método
3.2.1. Determinación de la radiación efectiva en Nuevo Chimbote
3.2.1.1.Estimación de la radiación directa para procesos de cálculos
El sol, como cualquier estrella, es un reactor de fusión en el que cada segundo se
generan 600 millones de toneladas de helio a través del ciclo protón-protón. Estos
procesos de fusión liberan una gran cantidad de potencia calorífica, la cual se
evalúa en 3.86𝑥1023 Wth. A pesar de que la estructura solar es muy compleja,
para el aprovechamiento de su energía se adopta un modelo simplificado, donde
se toma el sol como a un cuerpo negro que radia energía a una temperatura de
5.780K, ya que su distribución espectral es muy parecida a la de dicho cuerpo
negro.
32
Figura 9. Espectro de Radiación Solar. Fuente: Paredes, 2012.
33
Declinación (δ): es la ubicación angular del sol en el medio día solar, en
referencia al plano del Ecuador. Sabiéndose que su valor está definido y varía
entre ± 23º15´.
Ángulo horario (τ): esta defino por ser el movimiento angular del sol en relación
al mediodía, siendo este el instante del día donde el sol se encuentra más alto en
el cielo. Siendo cada hora correspondiente a 15° de longitud.
34
Para desarrollar el proyecto, es necesario definir la radiación directa, sabiendo
que los Concentradores Solares Cilíndricos Parabólicos únicamente captan este
tipo de radiación, por ende, las demás radiaciones (difusa y reflejada) se
descartan. Para minimizar los cálculos, se plasmó una serie de prototipos
estadísticos el cual permitirán encontrar varios valores de la radiación directa I o
teniendo en cuenta las siguientes ecuaciones:
−1
((3.8).cos(𝜃𝑠−1.6))
Io (media) = 1230.𝑒 , siendo para condiciones normales.
−1
(6.cos(𝜃𝑠−1))
Io (máxima) = 1220. 𝑒 , siendo para cielos muy puros.
−1
((2.3).cos(𝜃𝑠−3))
Io (mínima) = 1260. 𝑒 , siendo para zonas industriales.
35
El ángulo ϴ es muy importante en esta clase de concentrador ya que es un factor
multiplicativo de la radiación directa, otorgando como producto la energía por
unidad de superficie sobre el Concentrador Solar Cilíndrico Parabólico.
Giro del colector sobre el eje horizontal Este-Oeste, teniendo un control diario:
lo vectores n y s tienen igual dirección al mediodía de cada día.
cosϴ = (senδ . senλ) + (cosδ . cosλ. cos τ) (3.3)
36
3.2.1.2. Cálculos
(284+𝑛)
δ = 23.45º . [360.( )] , siendo “n” el nº de día natural. (3.7)
365
37
Tabla 2. Declinación media mensual en grados sexagesimal y en radianes.
Declinación
Meses Día del año sexagesimal radianes
Enero 17 -20.92 -0.37
Febrero 47 -12.95 -0.23
Marzo 75 -2.42 -0.04
Abril 105 9.41 0.16
Mayo 135 18.79 0.33
Junio 162 23.09 0.40
Julio 198 21.18 0.37
Agosto 228 13.45 0.23
Septiembre 258 2.22 0.04
Octubre 288 -9.60 -0.17
Noviembre 318 -18.91 -0.33
Diciembre 344 -23.05 -0.40
Ángulo horario
Tiempo
τ en seg. τ en rad.
(h)
6 -90 -1.57
7 -75 -1.31
8 -60 -1.05
9 -45 -0.79
10 -30 -0.52
11 -15 -0.26
12 0 0.00
13 15 0.26
14 30 0.52
15 45 0.79
16 60 1.05
17 75 1.31
18 90 1.57
38
Para hallar el valor de la radiación directa (Io) se tuvo en cuenta la región de
Nuevo Chimbote en condiciones normales, y se tiene lo siguiente:
−1
((3.8) cos(𝜃𝑠−1.6))
Io (media) = 1230.𝑒 (3.8)
ENERO
Angulo horario Angulo Declinación
Tiempo latitud cos ϴ Io media
seg. horario rad. rad.
07:00 75 1.31 -0.16 -0.37 0.295 489.622
08:00 60 1.05 -0.16 -0.37 0.516 744.376
09:00 45 0.79 -0.16 -0.37 0.705 856.265
10:00 30 0.52 -0.16 -0.37 0.856 911.970
11:00 15 0.26 -0.16 -0.37 0.947 938.466
12:00 0 0 -0.16 -0.37 0.978 945.306
13:00 15 0.26 -0.16 -0.37 0.947 938.466
14:00 30 0.52 -0.16 -0.37 0.856 911.970
15:00 45 0.79 -0.16 -0.37 0.705 856.265
16:00 60 1.05 -0.16 -0.37 0.516 744.376
17:00 75 1.31 -0.16 -0.37 0.295 489.622
39
Para calcular el rendimiento óptico se usó la siguiente fórmula:
Rendimiento Óptico
ρ: reflectividad, superficie reflectora. 0.76
τ: trasmisividad del cristal absorbedor. 1
α: Absortividad del tubo del absorbedor. 0.95
d: reflectancia espectral. 0.80
f: error óptico de reflexión. x
ɣ: factor de interceptación. 0.96
Rend.opt= 0.554496
Por tanto, se pudo hallar la Irradiación Efectiva (Iefec) como el producto de Io por
el cosϴ y por el rendimiento óptico.
Enero
40
En la gráfica 1 se presenta la curva que forma la radiación a lo largo de las
horas de luz (07:00 – 17:00 h) definido para el mes de enero.
500
400
300
200
100
0
6 7 8 9 10 11 12 13 14 15 16 17 18
1100
1000
900
800
700
600 Io media
500 Io máxima
400 I efectiva
300
200
100
0
6 7 8 2.
Gráfica 9 Comparativa
10 11 12 entre
13 la14irradiación
15 16 17 18 máxima y
normal,
efectiva para el mes de enero.
41
Finalmente, en la gráfica 3 se muestra la irradiación efectiva mensual para todo
un año.
42
En el estudio de la radiación solar, se vio por conveniente realizar una orientación
horizontal N – S, con su sistema de seguimiento constante.
Para llevar a cabo los cálculos en la orientación del concentrador solar N-S (ver
Anexo 2), se tuvo en cuenta la ecuación (3.5).
Por tanto, para poder tener un único valor de Irradiación efectiva para una
orientación Norte - Sur, se realizó una media aritmética de los diferentes
parámetros que conforman el cálculo de la misma (ver Anexo 2).
43
3.3.2. Dimensionamiento del concentrador solar cilíndrico parabólico
Para determinar la abertura, siendo el lado recto “L”, se usó la ecuación del arco
de una curva, el cual permitió conocer la longitud de lado recto de la parábola.
44
𝑏
S = ∫𝑎 √1 + (𝑓´(𝑥))2 dx
Donde:
𝑋2
Se sabe que f(x) depende del valor del foco, pues f(x) = y = , por lo que:
4∗𝐹
𝑋2
F = 0.50 m → f(x) = = f´(x) = x → Xmáx. =0.43
4∗0.5
𝑋2
F = 0.30 m → f(x) = = f´(x) = 1.67x → Xmáx. =0.39
4∗0.3
𝑋2
F = 0.20 m → f(x) = = f´(x) = 2.5 x → Xmáx. =0.39
4∗0.2
45
Una vez encontrado los valores máximos de la parábola, se procedió a realizar
una gráfica para escoger con que distancia focal nos quedaremos para el proceso
de dimensionamiento del concentrador solar cilíndrico parabólico.
Tal como se aprecia en la gráfica 5, se tiene que la abertura para un foco con 50
cm de distancia nos da un valor muy elevado para lo que se necesita en el
proyecto, puesto que la altura máxima sobre el vértice es de 9cm, por
consiguiente, este valor del foco queda descartado. En un foco de 20 cm se tiene
una altura máxima (Ymáx) y es la misma, viendo la gráfica se puede ver que la
parábola es más favorable para el propósito de este trabajo, ya que el foco y Ymáx
se ajustan perfectamente al mismo plano Z. Por último, en un foco de 30 cm se
cuenta con la relación distancia focal – altura máxima, que es casi el doble la
una de la otra. Descrito lo anterior se tomó la decisión que la distancia focal será
de 30 cm siendo el que más se asemejó a los propósitos del trabajo de
investigación. Las razones porque no se escogió el foco de 20 cm se explica más
adelante.
46
La parábola definitiva tiene las dimensiones siguientes: Foco = 0.30 metros, lado
curvo = 1.2 metros y lado recto = 0.78 metros. Y se muestra a continuación:
47
A continuación, se calcula el ángulo de abertura (𝛼)
𝛼 𝑋𝑚á𝑥
Tg ( 2 ) = (3.11)
𝐹−𝑌𝑚á𝑥
𝛼 𝑋𝑚á𝑥 0.39
Tg ( 2 ) = = = 2.29
𝐹−𝑌𝑚á𝑥 0.30−0.13
α = 2.32 rad.
Y para la distancia ρ mediante la fórmula (3.12):
ρ = √𝑊 2 + (𝐹 − 𝑌)2 (3.12)
ρ = 0.425 m
3.3.2.3. Cálculo del diámetro del tubo receptor
El sol no es una fuente solar luminosa puntual. Visto desde la superficie de la
tierra, el disco solar que conforma corresponde a un ángulo de apertura de 32´,
por tanto, la radiación solar directa sobre la superficie terrestre no está formada
por rayos perfectamente paralelos entre sí, sino que se distribuyen sobre un cono
de direcciones se semiángulo igual a 16´ alrededor de la línea que une el punto
de observación con el centro del disco solar. Así, no todos los rayos de este cono
alcanzarán un receptor de tamaño arbitrario.
48
Figura 15. Proyección del disco solar sobre el absorbedor.
El diámetro mínimo que debe tener el tubo receptor está en función de la distancia
focal a través del valor de ρ, puesto que, para un foco ubicado a 30 cm sobre el
vértice del colector, el diámetro mínimo del tubo se determinó haciendo uso de
la siguiente ecuación:
2.𝜌.𝑡𝑔16´
Φmin = 𝑎 (3.13)
cos(2 )
2.𝜌.𝑡𝑔16´ 2.(0.425).𝑡𝑔16´
Φmin = 𝑎 = 133
cos( 2) cos( 2 )
Se = Lr . Le = Lr . Fm . tan (ϴ)
(3.14)
Fm = F + (F . 𝐿𝑟 2/48 . 𝐹 2 ) (3.15)
49
Haciendo uso de las fórmulas (3.15) y (3.14), se procede a cálcular las perdidas
geométricas en el concentrador solar y se desarrolla a cotinuación:
Sc = 1.712 𝒎𝟐
Donde, se tienen todos los valores para reemplazar en dicha formula, siendo
indespensable el valor de la irradiación solar incidendete real (Ver anexo 5).
Reemplazando los valores tenemos:
Qsol = (1.712).(693.98).(0.95) = 1128.69 W
Finalmente :
Qsol = 1128.69 W
50
3.3.4. Flujo másico a la entrada y salida en el concentrador solar (kg/s)
ṁ = 𝜌. 𝑄 (3.17)
Donde:
𝑘𝑔
(1000 𝑠 ).(70𝑥10−3 𝑙) 1𝑚3
ṁ = 𝜌. 𝑄 = . 1000 𝑙
(15 𝑠)
Finalmente:
ṁ = 4.67x𝟏𝟎−𝟑 𝒌𝒈/𝒔
51
3.3.5. Gradiente de temperatura del concentrador solar (°C)
Figura 17. Medición del gradiente de temperatura. A) Ingreso y salida del flujo de agua. B)
Instrumentos para la medición de temperatura.
ΔT = Ts –Ts (3.18)
Qútil
ήglobal =
Qsol
52
Dónde: Q sol = energía solar incidente sobre el concentrador solar cilíndrico
parabólico, este valor ya se ha calculado anteriormente en el enunciado 3.3.
Convección en el interior a través del agua para el caso que tenemos suponemos
una temperatura interior de la tubería e igual a la temperatura en la pared interior
del tubo Ti = T1.
53
Dónde: h=coeficiente conectivo (W/°C 𝑚2 ), A=superficie de la tubería (𝑚2 ),
T2= temperatura al exterior de la tubería (°C) y Ta = temperatura ambiente (°C).
𝐾(𝑇1−𝑇2)
Q= (3.22)
𝑑
54
3.3.6.1. Cálculos
Una vez presentado la tubería de cobre con sus temperaturas se procede a
calcular las pérdidas de calor en la tubería, y se procede de la siguiente manera:
1
R i = R conv,1 = (3.23)
h1 A
r
Ln ( 2 )
r1
R1 = R TUBO = (3.24)
2πk1 L
1
R o = R conv,2 = (3.25)
h2 A3
Dónde: V=velocidad del fluido (m/s), Q=caudal del fluido (m3 /𝑠), D=diámetro
interior de la tubería (m).
55
Reemplazando los valores en la ecuación 3.26, se tiene lo siguiente.
4. (4.67x10−6 )
V= = 0.043 𝑚/𝑠
π(0.0117)2
Ahora se procede a calcular el número de Reynolds, haciendo uso de la
siguiente ecuación.
ρ. V. D
Re = (3.27)
μ
Donde,
56
Habiendo determinado el Nuset, ya se puede hallar el valor de h 1 y se realiza de
la siguiente manera:
Nu ∗ k
h= (3.29)
D
Reemplazando los valores encontrados se tiene lo siguiente.
(16.43)(0.02693)
h1 = = 37.87w/m2ºC
0.0117
encuentra en la tabla A-15 del libro de Yunus A. Cengel (Ver anexo 4).
Reemplazando los valores en la ecuación 3.30 se tiene lo siguiente:
(3.65). (0.01127)
Re = = 2661.02
1.742x10−5
Una vez obtenido el valor del número de Reynolds y sabiendo que es un flujo
laminar, se usa la ecuación de Crurchil y Berstein para hallar el Nuset (Nu) y se
detalla a continuación dicha ecuación:
1 1
0.3 + 0.62 (Re 2 ) (Pr 3 ) 5 4
Nu = 2 1 . {(1 + (Re /282000)8 )}5 (3.31)
{(1 + (0.4/Pr )3 )}4
Reemplazando los valores y se obtiene lo siguiente:
1 1
0.3 + 0.62 (2661.022 ) (0.72433 ) 5 4
Nu = 2 1 . {(1 + (2661.02/282000)8 )}5
{(1 + (0.4/0.7243)3 )}4
𝐍𝐮 = 𝟐𝟔. 𝟔𝟐
57
Habiendo determinado el Nuset, ya se puede hallar el valor de h 1 y se hace uso
de la ecuación 3.29 y se obtiene lo siguiente:
(26.62)(0.02693)
h2 = = 56.447 W/m2ºC
0.0127
Resistencia total:
RTotal = Ri + R1 + Ro
ºC ºC ºC
R Total = 0.7146 + 3.3588𝑥10−5 + 0.4429
W W W
ºC
R Total = 1.1575
W
q cd (66.3 − 22)°C
= = 38.27W
L ºC
1.1575
W
58
Habiendo hallado los valores de la energía en el concentrador solar y las pérdidas
de calor existentes en la tubería de cobre, se reemplazó en la ecuación de energía
útil y se obtuvo lo siguiente:
Qutil = Qc – Q perd.
Qutil = 864.56 - 57.41
Qutil = 807.15 W
Qútil
ήglobal = Qsol
807.15
ήglobal = 1128.69
ήglobal = 71 %
59
3.3.7. Temperatura media logarítmica en el intercambiador de calor a contra
flujo (ºc)
60
A B
61
Se reemplazó en la ecuación de la temperatura media logarítmica
ΔTml = 21.7 °C
A B
C D
62
Este sistema de seguimiento solar nos permitió ajustar el concentrador solar en
cualquier posición, dependiendo de la hora y de la posición que se encuentra el
sol. Se resume en tres posiciones más importantes a lo largo del día:
Primera posición del concentrador: orientado hacia el este, esto fue posible en
horas de la mañana, sabiendo que el sol sale por el este como se muestra en la
figura 22.
Segunda posición del concentrador: orientado perpendicular a la posición del sol, esto
es al medio día solar (12:00 h), a esta hora se tiene la mayor cantidad de radiación,
haciendo más eficiente al concentrador solar cilíndrico parabólico (Fig. 23).
63
Tercera posición del sol: orientado hacia el oeste, siendo esto posible en horas de la
tarde, sabiendo que el sol a esas horas se oculta por el oeste (Fig. 24)
64
CAPITULO IV
65
IV. RESULTADOS Y DISCUSIÓN
4.1. Irradiación efectiva teórico con orientación del eje este-oeste y norte- sur
En la tabla 8, se muestra la irradiación efectiva con orientación en el eje E-O y N-S
de la ciudad de Nuevo Chimbote, expresada en (W/𝑚2).
Los valores de irradiación efectiva reportados por Almanza et al. (2013), fueron
900 W/𝑚2 siendo estos valores superiores a los obtenidos en este trabajo de investigación
debido a su ubicación geográfica.
66
A continuación, se presenta un cuadro comparativo entre los valores teóricos y reales de
la radiación solar medidas para los meses de abril y mayo del 2018 (ver anexo 5).
Tabla 8.1. Cuadro comparativo entre los datos reales y teóricos de la radiación solar, medida
para los meses de abril y mayo del 2018.
Al ver en la tabla 8.1 se puede apreciar una comparación entre las Irradiaciones Solares
reales y teóricas. Por el lado de los valores obtenidos de Irradiación Solar real, se
realizaron las mediaciones para los meses de abril – mayo, siendo el día 1, 15 y 30 de
abril – mayo que se llevaron a cabo las mediciones de Irradiación solar incidente en el
Concentrador Solar Cilíndrico Parabólico. Todo lo descrito se muestra en los cuadros
existentes en el anexo 5.
La elección de la mejor irradiación efectiva para el cálculo de la energía incidente en el
Concentrador Solar Cilíndrico Parabólico, se hizo en base al promedio de los valores
obtenidos de los datos tomados en la práctica, en este caso fue 693.98 W/𝑚2 con una
orientación del eje de Norte – Sur, tal como se muestra en los cuadros y en la gráfica 8,
se hizo la comparación de dichos valores solamente para los meses de abril - mayo. Cabe
resaltar que dichos valores de Irradiación fueron obtenidos, mediante el uso de un
piranómetro.
67
.
Gráfica 8. Comparación de las irradiaciones solares reales y teóricas para los meses de
abril-mayo.
68
Utilizando las ecuaciones más importantes del modelado termodinámico del
concentrador solar cilíndrico parabólico presentados en la sección 3.2, se obtuvo los
valores de foco = 30 cm, habiéndose escogido entre tres valores tal como se muestra en
la gráfica 5, donde la abertura para un foco de 50 cm de distancia fue demasiado grande
para las necesidades térmicas del proyecto, además de que la altura máxima sobre el
vértice es de tan solo 9 cm, quedando por lo pronto descartado. Por otro lado, para un
foco de 20 cm, la altura Ymax es prácticamente la misma, siendo estructuralmente este el
caso más favorable, ya que el foco y la altura máxima se pueden ajustar al mismo plano
Z. Sin embargo, al calcular el diámetro de la tubería con esta distancia de 20 cm se
necesitaría una tubería con un diámetro de 13.8 cm, siendo este diámetro muy excesivo
dificultando la captación de energía, por consecuencia esta distancia también queda
descartado. Dicho lo anterior se decidió que la distancia focal adecuada es de 30 cm,
obteniendo con esta distancia focal un diámetro de 1cm.
En la tabla 9, el valor del ángulo de abertura obtenido fue 2.32 rad, y una distancia de
0.425 metros, siendo estos valores los que definieron la posición del tubo de cobre
respecto al concentrador solar cilíndrico parabólico, teniendo en cuenta que la distancia
focal es 30 cm. Así mismo, el diámetro el tubo receptor obtenido fue 1 cm, siendo este
un diámetro aceptable para el diseño dado que el presente trabajo de investigación es a
pequeña escala, teniéndose un flujo másico pequeño. Finalmente, el valor reportado para
las perdidas geométricas fue 0.088 𝑚2 , implicando esto una reducción aproximada del 5
% de la superficie total de reflexión.
En la tabla 10, se presenta las características generales del diseño del concentrador solar
cilíndrico parabólico, los cuales fueron obtenidos mediante las ecuaciones presentadas
en el capítulo 3 y se tiene lo siguiente:
69
Tabla 10. Resultado de cálculos en el concentrador solar cilíndrico parabólico.
Según los datos obtenidos, la energía solar incidente reportada es 1128.69 W, en base a
estos resultados obtenidos es posible asegurar que el diseño del concentrador solar es
factible termodinámicamente, puesto que como se evidencia en el enunciado 3.3, dicho
valor es directamente proporcional de la irradiación efectiva, el área del concentrador y
el ángulo de incidencia. Se puede notar que el valor de la energía incidente no es muy
alto, debido a que se presenta un trabajo a pequeña escala con la finalidad de generar
agua caliente y abastecer a un intercambiador de calor a contra flujo.
Teniendo los valores de energía solar incidente y del calor útil se determinó que el
concentrador solar cilíndrico parabólico tiene un rendimiento de 71 %, que es adecuado
en este tipo de investigación y se pudo demostrar en la práctica.
70
4.4. Temperatura media logarítmica en el intercambiador de calor a contra flujo (ºc)
Según los datos obtenidos en el enunciado 3.7, se tiene 21.7°C de temperatura media
logarítmica en el intercambiador de calor, temperatura requerida para realizar las
prácticas de laboratorio de transferencia de calor por parte de los alumnos de la Escuela
Académico Profesional de Ingeniería en Energía, siendo este el principal objetico en el
trabajo de investigación.
71
V. CONCLUSIONES Y RECOMENDACIONES
5.1. Conclusiones
Se determinó el potencial de energía incidente 1128.69 W en el concentrador solar
cilíndrico parabólico, habiendo hecho uso de la irradiación efectiva 693.98 W/𝑚2
orientado el eje Norte-Sur, siendo esta orientación donde se aprovecha la mayor cantidad
de radiación solar debido a que el concentrador cuenta con un sistema manual de
seguimiento solar.
Se evaluó el área de captación del concentrador solar teniendo una plancha de aluminio
de 1.8 𝑚2 , obteniéndose como el área de captación útil 1.712 𝑚2 debido a las perdidas
existentes en la construcción del concentrador, siendo vidrio de alta calidad el material
reflectante.
72
5.2. Recomendaciones
Se recomienda implementar una bomba en el sistema instalado que consta de un
concentrador solar, el intercambiador de calor y un tanque de abastecimiento de agua,
esto para que se cumpla el proceso de recirculación del agua y así hacer más eficiente el
concentrador solar.
73
VI. REFERENCIAS BIBLIOGRÁFICAS Y VIRTUALES
Almanza Gómez, G., & Cabarcas Toscano, J. (2013). Diseño de un sistema de secado
industrial basado en la concentración de energía solar usando un colector tipo
cilíndrico parabólico acoplado a un intercambiador de calor de tubo coraza.
(Tesis de pregrado). Universidad de Cartagena. Colombia.
Almanza, R. & Muñoz, F. (2003). Ingeniería de la energía solar. En serie ingeniería (Edit.),
México: cromo color.
Bagarozza Bastidas, E., & Del Risco De León, J. (2002). Diseño de un colector solar
cilíndrico parabólico y construcción de un prototipo. (Tesis de maestría).
Corporación Universitaria Tecnológica Del Bolívar. Cartagena- Colombia.
74
Daghero, J., Garnica, J., Buitrago, A., Dubini, D., Lorenzo, C., Manero, C., . . . Sayago, U.
(2015). Concentrador cilíndrico parabólico aplicado a la generación de vapor de
uso industrial. evaluación óptica y térmica de un prototipo. Energías Renovables
y Medio Ambiente, 19(2), 45-54.
Kostyuk, A., & Frolov, V. (1988). Steam and gas turbines. Moscú, Rusia: Editorial Mir.
Montes Pita, M. (2008). Análisis y propuestas de sistemas solares de alta exergía que
emplean agua como fluido calorífico. (Tesis doctoral). Universidad Politécnica
de Madrid. Madrid.
Polo Bravo, C., Sacari Sacari, E., & Choque Chacolla, J. (2011). Preliminary results with
solar concentrating systems for water vapor generation. Ciencia & desarrollo.
13(1), 11-20.
75
Salgado Conrado, L. (2013). Modelado y control de un concentrador solar cilíndrico
parabólico con un mecanismo por desplazamiento de agua. (Tesis de doctorado).
Instituto Politécnico Nacional. México.
Sánchez, C., Piedra Castillo, D., & Mendoza, I. (2008). Validación y ajuste de modelos de
radiación solar directa para la ciudad de Bogotá a partir de datos
experimentales tomados en la Universidad Distrital Francisco José de Caldas.
(Tesis posgrado). Universidad Distrital. Colombia.
Vásquez Chigne, L., & Zuñiga Anticona, B. (2015). Proyecto de pre factibilidad para la
implementación de energía solar fotovoltaica y térmica en el campamento
minero Comihuasa. (Tesis de pregrado). Universidad Peruana de Ciencias
Aplicadas. Lima. Perú.
Venegas, E., Jaramillo, O., Castrejón, R., Aguilar, J., & Montemayor, F. (2012). Desing,
construction and testing of a parabolic trough solar concentrator for hot water
and low enthalpy steam generation. Renewable sustainable energy. 2(4). 303-
320.
76
VII. ANEXOS
Enero
Angulo Angulo
Tiempo Latitud Declinación Iefec I efec
horario horario cos ϴ Io media Io máxima Io mínima I efec media
(h) (°) (δ) máxima mínima
(τ seg) (τ rad)
07:00 75 1.31 -0.44 -0.36 0.29 489.62 666.44 308.55 79.43 108.108 50.05
08:00 60 1.05 -0.44 -0.36 0.52 744.38 882.72 567.12 211.08 250.314 160.82
09:00 45 0.79 -0.44 -0.36 0.71 856.26 967.69 701.91 332.23 375.459 272.34
10:00 30 0.52 -0.44 -0.36 0.86 911.97 1008.33 773.48 429.54 474.925 364.31
11:00 15 0.26 -0.44 -0.36 0.95 938.47 1027.45 807.85 488.85 535.202 420.81
12:00 0 0.00 -0.44 -0.36 0.98 945.31 1032.68 815.24 508.50 555.497 438.53
13:00 15 0.00 -0.44 -0.36 0.95 938.47 1027.45 807.85 488.85 535.202 420.81
14:00 30 0.00 -0.44 -0.36 0.86 911.97 1008.33 773.48 429.54 474.925 364.31
15:00 45 0.00 -0.44 -0.36 0.71 856.26 967.69 701.91 332.23 375.459 272.34
16:00 60 0.00 -0.44 -0.36 0.52 744.38 882.72 567.12 211.08 250.314 160.82
17:00 75 0.00 -0.44 -0.36 0.29 489.62 666.44 308.55 79.43 108.108 50.05
77
ANEXO 2. IRRADIACIÓN EFECTIVA ANUAL- ORIENTACION ESTE-OESTE / NORTE-SUR
Enero
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 -0.37 0.295 489.622 0.55 79.426
08:00 60 1.05 -0.16 -0.37 0.516 744.376 0.55 211.084
09:00 45 0.79 -0.16 -0.37 0.705 856.265 0.55 332.226
10:00 30 0.52 -0.16 -0.37 0.856 911.970 0.55 429.541
11:00 15 0.26 -0.16 -0.37 0.947 938.466 0.55 488.849
12:00 0 0 -0.16 -0.37 0.978 945.306 0.55 508.496
13:00 15 0.26 -0.16 -0.37 0.947 938.466 0.55 488.849
14:00 30 0.52 -0.16 -0.37 0.856 911.970 0.55 429.541
15:00 45 0.79 -0.16 -0.37 0.705 856.265 0.55 332.226
16:00 60 1.05 -0.16 -0.37 0.516 744.376 0.55 211.084
17:00 75 1.31 -0.16 -0.37 0.295 489.622 0.55 79.426
0.69244591 802.43 326.432
78
Febrero
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 -0.23 0.28417416 489.6218002 0.55 76.52583
08:00 60 1.05 -0.16 -0.23 0.5146012 744.3759507 0.55 210.6812
09:00 45 0.79 -0.16 -0.23 0.71287822 856.2648251 0.55 335.7269
10:00 30 0.52 -0.16 -0.23 0.87049483 911.9698369 0.55 436.6258
11:00 15 0.26 -0.16 -0.23 0.96524404 938.4661737 0.55 498.2169
12:00 0 0 -0.16 -0.23 0.997551 945.3062152 0.55 518.6451
13:00 15 0.26 -0.16 -0.23 0.96524404 938.4661737 0.55 498.2169
14:00 30 0.52 -0.16 -0.23 0.87049483 911.9698369 0.55 436.6258
15:00 45 0.79 -0.16 -0.23 0.71287822 856.2648251 0.55 335.7269
16:00 60 1.05 -0.16 -0.23 0.5146012 744.3759507 0.55 210.6812
17:00 75 1.31 -0.16 -0.23 0.28417416 489.6218002 0.55 76.52583
Marzo
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 -0.04 0.260724 489.6218002 0.55 70.21088
08:00 60 1.05 -0.16 -0.04 0.49719382 744.3759507 0.55 203.5545
09:00 45 0.79 -0.16 -0.04 0.70067052 856.2648251 0.55 329.9777
10:00 30 0.52 -0.16 -0.04 0.8624205 911.9698369 0.55 432.5758
11:00 15 0.26 -0.16 -0.04 0.95965445 938.4661737 0.55 495.3318
12:00 0 0 -0.16 -0.04 0.99280864 945.3062152 0.55 516.1795
13:00 15 0.26 -0.16 -0.04 0.95965445 938.4661737 0.55 495.3318
14:00 30 0.52 -0.16 -0.04 0.8624205 911.9698369 0.55 432.5758
15:00 45 0.79 -0.16 -0.04 0.70067052 856.2648251 0.55 329.9777
16:00 60 1.05 -0.16 -0.04 0.49719382 744.3759507 0.55 203.5545
79
17:00 75 1.31 -0.16 -0.04 0.260724 489.6218002 0.55 70.21088
Abril
Tiempo Angulo horario Angulo horario (τ Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 0.16 0.22592292 489.6218002 0.55 60.83923
08:00 60 1.05 -0.16 0.16 0.45955926 744.3759507 0.55 188.1467
09:00 45 0.79 -0.16 0.16 0.66059782 856.2648251 0.55 311.1057
10:00 30 0.52 -0.16 0.16 0.82040965 911.9698369 0.55 411.5039
11:00 15 0.26 -0.16 0.16 0.9164785 938.4661737 0.55 473.0462
12:00 0 0 -0.16 0.16 0.94923542 945.3062152 0.55 493.525
13:00 15 0.26 -0.16 0.16 0.9164785 938.4661737 0.55 473.0462
14:00 30 0.52 -0.16 0.16 0.82040965 911.9698369 0.55 411.5039
15:00 45 0.79 -0.16 0.16 0.66059782 856.2648251 0.55 311.1057
16:00 60 1.05 -0.16 0.16 0.45955926 744.3759507 0.55 188.1467
17:00 75 1.31 -0.16 0.16 0.22592292 489.6218002 0.55 60.83923
Mayo
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 0.33 0.18919536 489.6218002 0.55 50.94879
08:00 60 1.05 -0.16 0.33 0.41308491 744.3759507 0.55 169.1198
09:00 45 0.79 -0.16 0.33 0.60573658 856.2648251 0.55 285.269
10:00 30 0.52 -0.16 0.33 0.75888142 911.9698369 0.55 380.6423
11:00 15 0.26 -0.16 0.33 0.85094248 938.4661737 0.55 439.2194
12:00 0 0 -0.16 0.33 0.88233286 945.3062152 0.55 458.7411
13:00 15 0.26 -0.16 0.33 0.85094248 938.4661737 0.55 439.2194
14:00 30 0.52 -0.16 0.33 0.75888142 911.9698369 0.55 380.6423
80
15:00 45 0.79 -0.16 0.33 0.60573658 856.2648251 0.55 285.269
16:00 60 1.05 -0.16 0.33 0.41308491 744.3759507 0.55 169.1198
17:00 75 1.31 -0.16 0.33 0.18919536 489.6218002 0.55 50.94879
Junio
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 0.40 0.17242071 489.6218002 0.55 46.43152
08:00 60 1.05 -0.16 0.40 0.3903982 744.3759507 0.55 159.8317
09:00 45 0.79 -0.16 0.40 0.57796268 856.2648251 0.55 272.189
10:00 30 0.52 -0.16 0.40 0.72706355 911.9698369 0.55 364.683
11:00 15 0.26 -0.16 0.40 0.81669363 938.4661737 0.55 421.5416
12:00 0 0 -0.16 0.40 0.84725511 945.3062152 0.55 440.5035
13:00 15 0.26 -0.16 0.40 0.81669363 938.4661737 0.55 421.5416
14:00 30 0.52 -0.16 0.40 0.72706355 911.9698369 0.55 364.683
15:00 45 0.79 -0.16 0.40 0.57796268 856.2648251 0.55 272.189
16:00 60 1.05 -0.16 0.40 0.3903982 744.3759507 0.55 159.8317
17:00 75 1.31 -0.16 0.40 0.17242071 489.6218002 0.55 46.43152
Julio
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 0.37 0.17971815 489.6218002 0.55 48.39666
08:00 60 1.05 -0.16 0.37 0.40036192 744.3759507 0.55 163.9109
09:00 45 0.79 -0.16 0.37 0.59022067 856.2648251 0.55 277.9619
10:00 30 0.52 -0.16 0.37 0.74114533 911.9698369 0.55 371.7462
11:00 15 0.26 -0.16 0.37 0.83187177 938.4661737 0.55 429.3759
12:00 0 0 -0.16 0.37 0.86280707 945.3062152 0.55 448.5893
81
13:00 15 0.26 -0.16 0.37 0.83187177 938.4661737 0.55 429.3759
14:00 30 0.52 -0.16 0.37 0.74114533 911.9698369 0.55 371.7462
15:00 45 0.79 -0.16 0.37 0.59022067 856.2648251 0.55 277.9619
16:00 60 1.05 -0.16 0.37 0.40036192 744.3759507 0.55 163.9109
17:00 75 1.31 -0.16 0.37 0.17971815 489.6218002 0.55 48.39666
Agosto
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 0.23 0.21153222 489.6218002 0.55 56.96393
08:00 60 1.05 -0.16 0.23 0.44195926 744.3759507 0.55 180.9411
09:00 45 0.79 -0.16 0.23 0.64023628 856.2648251 0.55 301.5165
10:00 30 0.52 -0.16 0.23 0.79785289 911.9698369 0.55 400.1898
11:00 15 0.26 -0.16 0.23 0.8926021 938.4661737 0.55 460.7223
12:00 0 0 -0.16 0.23 0.92490906 945.3062152 0.55 480.8773
13:00 15 0.26 -0.16 0.23 0.8926021 938.4661737 0.55 460.7223
14:00 30 0.52 -0.16 0.23 0.79785289 911.9698369 0.55 400.1898
15:00 45 0.79 -0.16 0.23 0.64023628 856.2648251 0.55 301.5165
16:00 60 1.05 -0.16 0.23 0.44195926 744.3759507 0.55 180.9411
17:00 75 1.31 -0.16 0.23 0.21153222 489.6218002 0.55 56.96393
Septiembre
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 0.04 0.24798194 489.6218002 0.55 66.77955
08:00 60 1.05 -0.16 0.04 0.48445176 744.3759507 0.55 198.3378
09:00 45 0.79 -0.16 0.04 0.68792846 856.2648251 0.55 323.9769
10:00 30 0.52 -0.16 0.04 0.84967845 911.9698369 0.55 426.1846
82
11:00 15 0.26 -0.16 0.04 0.94691239 938.4661737 0.55 488.7549
12:00 0 0 -0.16 0.04 0.98006658 945.3062152 0.55 509.5547
13:00 15 0.26 -0.16 0.04 0.94691239 938.4661737 0.55 488.7549
14:00 30 0.52 -0.16 0.04 0.84967845 911.9698369 0.55 426.1846
15:00 45 0.79 -0.16 0.04 0.68792846 856.2648251 0.55 323.9769
16:00 60 1.05 -0.16 0.04 0.48445176 744.3759507 0.55 198.3378
17:00 75 1.31 -0.16 0.04 0.24798194 489.6218002 0.55 66.77955
Octubre
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 -0.17 0.27784092 489.6218002 0.55 74.82034
08:00 60 1.05 -0.16 -0.17 0.51108855 744.3759507 0.55 209.2431
09:00 45 0.79 -0.16 -0.17 0.71179263 856.2648251 0.55 335.2156
10:00 30 0.52 -0.16 -0.17 0.87133857 911.9698369 0.55 437.049
11:00 15 0.26 -0.16 -0.17 0.96724758 938.4661737 0.55 499.251
12:00 0 0 -0.16 -0.17 0.99995 945.3062152 0.55 519.8924
13:00 15 0.26 -0.16 -0.17 0.96724758 938.4661737 0.55 499.251
14:00 30 0.52 -0.16 -0.17 0.87133857 911.9698369 0.55 437.049
15:00 45 0.79 -0.16 -0.17 0.71179263 856.2648251 0.55 335.2156
16:00 60 1.05 -0.16 -0.17 0.51108855 744.3759507 0.55 209.2431
17:00 75 1.31 -0.16 -0.17 0.27784092 489.6218002 0.55 74.82034
Noviembre
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 -0.33 0.29244726 489.6218002 0.55 78.75371
08:00 60 1.05 -0.16 -0.33 0.51633682 744.3759507 0.55 211.3918
83
09:00 45 0.79 -0.16 -0.33 0.70898849 856.2648251 0.55 333.895
10:00 30 0.52 -0.16 -0.33 0.86213332 911.9698369 0.55 432.4318
11:00 15 0.26 -0.16 -0.33 0.95419439 938.4661737 0.55 492.5135
12:00 0 0 -0.16 -0.33 0.98558477 945.3062152 0.55 512.4237
13:00 15 0.26 -0.16 -0.33 0.95419439 938.4661737 0.55 492.5135
14:00 30 0.52 -0.16 -0.33 0.86213332 911.9698369 0.55 432.4318
15:00 45 0.79 -0.16 -0.33 0.70898849 856.2648251 0.55 333.895
16:00 60 1.05 -0.16 -0.33 0.51633682 744.3759507 0.55 211.3918
17:00 75 1.31 -0.16 -0.33 0.29244726 489.6218002 0.55 78.75371
Diciembre
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 -0.40 0.29650358 489.6218002 0.55 79.84604
08:00 60 1.05 -0.16 -0.40 0.51448107 744.3759507 0.55 210.632
09:00 45 0.79 -0.16 -0.40 0.70204554 856.2648251 0.55 330.6253
10:00 30 0.52 -0.16 -0.40 0.85114641 911.9698369 0.55 426.9209
11:00 15 0.26 -0.16 -0.40 0.9407765 938.4661737 0.55 485.5878
12:00 0 0 -0.16 -0.40 0.97133797 945.3062152 0.55 505.0165
13:00 15 0.26 -0.16 -0.40 0.9407765 938.4661737 0.55 485.5878
14:00 30 0.52 -0.16 -0.40 0.85114641 911.9698369 0.55 426.9209
15:00 45 0.79 -0.16 -0.40 0.70204554 856.2648251 0.55 330.6253
16:00 60 1.05 -0.16 -0.40 0.51448107 744.3759507 0.55 210.632
17:00 75 1.31 -0.16 -0.40 0.29650358 489.6218002 0.55 79.84604
84
CUADRO RESUMEN PROMEDIOS
Meses del año cos ϴ Io Iefec ή óptica
ENERO 0.69 802.43 326.43 0.55
FEBRERO 0.70 802.43 330.38 0.55
MARZO 0.69 802.43 325.41 0.55
ABRIL 0.65 802.43 307.53 0.55 I efectiva final
MAYO 0.59 802.43 282.65 0.55 288.48 W/𝒎𝟐
JUNIO 0.57 802.43 269.99 0.55
JULIO 0.58 802.43 275.58 0.55
AGOSTO 0.63 802.43 298.32 0.55
SEPTIEMBRE 0.67 802.43 319.78 0.55
OCTUBRE 0.70 802.43 330.10 0.55
NOVIEMBRE 0.70 802.43 328.22 0.55
DICIEMBRE 0.69 802.43 324.75 0.55
Enero
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 -0.37 0.94785684 489.6218 0.55 255.250255
08:00 60 1.05 -0.16 -0.37 0.95909363 744.375951 0.55 392.659428
09:00 45 0.79 -0.16 -0.37 0.96761001 856.264825 0.55 455.691727
10:00 30 0.52 -0.16 -0.37 0.97363921 911.969837 0.55 488.361274
11:00 15 0.26 -0.16 -0.37 0.9769535 938.466174 0.55 504.260799
12:00 0 0 -0.16 -0.37 0.97803091 945.306215 0.55 508.496286
13:00 15 0.26 -0.16 -0.37 0.9769535 938.466174 0.55 504.260799
14:00 30 0.52 -0.16 -0.37 0.97363921 911.969837 0.55 488.361274
15:00 45 0.79 -0.16 -0.37 0.96761001 856.264825 0.55 455.691727
85
16:00 60 1.05 -0.16 -0.37 0.95909363 744.375951 0.55 392.659428
17:00 75 1.31 -0.16 -0.37 0.94785684 489.6218 0.55 255.250255
0.96621248 802.427581 427.358478
Febrero
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 -0.23 0.98272586 489.6218 0.55 264.640204
08:00 60 1.05 -0.16 -0.23 0.98900516 744.375951 0.55 404.90541
09:00 45 0.79 -0.16 -0.23 0.99326285 856.264825 0.55 467.772822
10:00 30 0.52 -0.16 -0.23 0.99590126 911.969837 0.55 499.527551
11:00 15 0.26 -0.16 -0.23 0.99717175 938.466174 0.55 514.696576
12:00 0 0 -0.16 -0.23 0.997551 945.306215 0.55 518.645138
13:00 15 0.26 -0.16 -0.23 0.99717175 938.466174 0.55 514.696576
14:00 30 0.52 -0.16 -0.23 0.99590126 911.969837 0.55 499.527551
15:00 45 0.79 -0.16 -0.23 0.99326285 856.264825 0.55 467.772822
16:00 60 1.05 -0.16 -0.23 0.98900516 744.375951 0.55 404.90541
17:00 75 1.31 -0.16 -0.23 0.98272586 489.6218 0.55 264.640204
Marzo
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 -0.04 0.99999877 489.6218 0.55 269.291659
08:00 60 1.05 -0.16 -0.04 0.99921045 744.375951 0.55 409.083524
09:00 45 0.79 -0.16 -0.04 0.99736353 856.264825 0.55 469.70402
10:00 30 0.52 -0.16 -0.04 0.99512019 911.969837 0.55 499.135776
11:00 15 0.26 -0.16 -0.04 0.99343917 938.466174 0.55 512.76998
12:00 0 0 -0.16 -0.04 0.99280864 945.306215 0.55 516.179496
86
13:00 15 0.26 -0.16 -0.04 0.99343917 938.466174 0.55 512.76998
14:00 30 0.52 -0.16 -0.04 0.99512019 911.969837 0.55 499.135776
15:00 45 0.79 -0.16 -0.04 0.99736353 856.264825 0.55 469.70402
16:00 60 1.05 -0.16 -0.04 0.99921045 744.375951 0.55 409.083524
17:00 75 1.31 -0.16 -0.04 0.99999877 489.6218 0.55 269.291659
Abril
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 0.16 0.98023457 489.6218 0.55 263.969319
08:00 60 1.05 -0.16 0.16 0.97186396 744.375951 0.55 397.887686
09:00 45 0.79 -0.16 0.16 0.9634227 856.264825 0.55 453.719734
10:00 30 0.52 -0.16 0.16 0.95587407 911.969837 0.55 479.450575
11:00 15 0.26 -0.16 0.16 0.95097093 938.466174 0.55 490.849729
12:00 0 0 -0.16 0.16 0.94923542 945.306215 0.55 493.524977
13:00 15 0.26 -0.16 0.16 0.95097093 938.466174 0.55 490.849729
14:00 30 0.52 -0.16 0.16 0.95587407 911.969837 0.55 479.450575
15:00 45 0.79 -0.16 0.16 0.9634227 856.264825 0.55 453.719734
16:00 60 1.05 -0.16 0.16 0.97186396 744.375951 0.55 397.887686
17:00 75 1.31 -0.16 0.16 0.98023457 489.6218 0.55 263.969319
87
Mayo
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 0.33 0.93342687 489.6218 0.55 251.364378
08:00 60 1.05 -0.16 0.33 0.91872458 744.375951 0.55 376.132065
09:00 45 0.79 -0.16 0.33 0.90472837 856.264825 0.55 426.077896
10:00 30 0.52 -0.16 0.33 0.89267378 911.969837 0.55 447.75036
11:00 15 0.26 -0.16 0.33 0.88501625 938.466174 0.55 456.806799
12:00 0 0 -0.16 0.33 0.88233286 945.306215 0.55 458.741104
13:00 15 0.26 -0.16 0.33 0.88501625 938.466174 0.55 456.806799
14:00 30 0.52 -0.16 0.33 0.89267378 911.969837 0.55 447.75036
15:00 45 0.79 -0.16 0.33 0.90472837 856.264825 0.55 426.077896
16:00 60 1.05 -0.16 0.33 0.91872458 744.375951 0.55 376.132065
17:00 75 1.31 -0.16 0.33 0.93342687 489.6218 0.55 251.364378
Junio
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 0.40 0.90646461 489.6218 0.55 244.103659
08:00 60 1.05 -0.16 0.40 0.88923078 744.375951 0.55 364.057103
09:00 45 0.79 -0.16 0.40 0.87299568 856.264825 0.55 411.133524
10:00 30 0.52 -0.16 0.40 0.85911073 911.969837 0.55 430.915692
11:00 15 0.26 -0.16 0.40 0.8503273 938.466174 0.55 438.901873
12:00 0 0 -0.16 0.40 0.84725511 945.306215 0.55 440.503537
13:00 15 0.26 -0.16 0.40 0.8503273 938.466174 0.55 438.901873
14:00 30 0.52 -0.16 0.40 0.85911073 911.969837 0.55 430.915692
15:00 45 0.79 -0.16 0.40 0.87299568 856.264825 0.55 411.133524
16:00 60 1.05 -0.16 0.40 0.88923078 744.375951 0.55 364.057103
88
17:00 75 1.31 -0.16 0.40 0.90646461 489.6218 0.55 244.103659
Julio
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 0.37 0.91855344 489.6218 0.55 247.359084
08:00 60 1.05 -0.16 0.37 0.90239768 744.375951 0.55 369.447723
09:00 45 0.79 -0.16 0.37 0.88711758 856.264825 0.55 417.784169
10:00 30 0.52 -0.16 0.37 0.87401461 911.969837 0.55 438.39123
11:00 15 0.26 -0.16 0.37 0.86571282 938.466174 0.55 446.843209
12:00 0 0 -0.16 0.37 0.86280707 945.306215 0.55 448.589287
13:00 15 0.26 -0.16 0.37 0.86571282 938.466174 0.55 446.843209
14:00 30 0.52 -0.16 0.37 0.87401461 911.969837 0.55 438.39123
15:00 45 0.79 -0.16 0.37 0.88711758 856.264825 0.55 417.784169
16:00 60 1.05 -0.16 0.37 0.90239768 744.375951 0.55 369.447723
17:00 75 1.31 -0.16 0.37 0.91855344 489.6218 0.55 247.359084
Agosto
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 0.23 0.96423081 489.6218 0.55 259.659633
08:00 60 1.05 -0.16 0.23 0.95322862 744.375951 0.55 390.258251
09:00 45 0.79 -0.16 0.23 0.94248513 856.264825 0.55 443.859277
10:00 30 0.52 -0.16 0.23 0.93307412 911.969837 0.55 468.014497
11:00 15 0.26 -0.16 0.23 0.92703503 938.466174 0.55 478.49506
12:00 0 0 -0.16 0.23 0.92490906 945.306215 0.55 480.877256
13:00 15 0.26 -0.16 0.23 0.92703503 938.466174 0.55 478.49506
14:00 30 0.52 -0.16 0.23 0.93307412 911.969837 0.55 468.014497
89
15:00 45 0.79 -0.16 0.23 0.94248513 856.264825 0.55 443.859277
16:00 60 1.05 -0.16 0.23 0.95322862 744.375951 0.55 390.258251
17:00 75 1.31 -0.16 0.23 0.96423081 489.6218 0.55 259.659633
Septiembre
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 0.04 0.99675252 489.6218 0.55 268.417469
08:00 60 1.05 -0.16 0.04 0.99293168 744.375951 0.55 406.512957
09:00 45 0.79 -0.16 0.04 0.98845354 856.264825 0.55 465.5079
10:00 30 0.52 -0.16 0.04 0.98409782 911.969837 0.55 493.607141
11:00 15 0.26 -0.16 0.04 0.98113597 938.466174 0.55 506.419606
12:00 0 0 -0.16 0.04 0.98006658 945.306215 0.55 509.554665
13:00 15 0.26 -0.16 0.04 0.98113597 938.466174 0.55 506.419606
14:00 30 0.52 -0.16 0.04 0.98409782 911.969837 0.55 493.607141
15:00 45 0.79 -0.16 0.04 0.98845354 856.264825 0.55 465.5079
16:00 60 1.05 -0.16 0.04 0.99293168 744.375951 0.55 406.512957
17:00 75 1.31 -0.16 0.04 0.99675252 489.6218 0.55 268.417469
Octubre
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 -0.17 0.99196235 489.6218 0.55 267.127514
08:00 60 1.05 -0.16 -0.17 0.99604127 744.375951 0.55 407.786041
09:00 45 0.79 -0.16 -0.17 0.99840246 856.264825 0.55 470.193298
10:00 30 0.52 -0.16 -0.17 0.99952695 911.969837 0.55 501.346137
11:00 15 0.26 -0.16 -0.17 0.99988329 938.466174 0.55 516.096157
90
12:00 0 0 -0.16 -0.17 0.99995 945.306215 0.55 519.892423
13:00 15 0.26 -0.16 -0.17 0.99988329 938.466174 0.55 516.096157
14:00 30 0.52 -0.16 -0.17 0.99952695 911.969837 0.55 501.346137
15:00 45 0.79 -0.16 -0.17 0.99840246 856.264825 0.55 470.193298
16:00 60 1.05 -0.16 -0.17 0.99604127 744.375951 0.55 407.786041
17:00 75 1.31 -0.16 -0.17 0.99196235 489.6218 0.55 267.127514
Noviembre
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 -0.33 0.95969591 489.6218 0.55 258.438421
08:00 60 1.05 -0.16 -0.33 0.96954599 744.375951 0.55 396.938694
09:00 45 0.79 -0.16 -0.33 0.97687323 856.264825 0.55 460.054201
10:00 30 0.52 -0.16 -0.33 0.9819569 911.969837 0.55 492.533289
11:00 15 0.26 -0.16 -0.33 0.98470178 938.466174 0.55 508.260122
12:00 0 0 -0.16 -0.33 0.98558477 945.306215 0.55 512.423673
13:00 15 0.26 -0.16 -0.33 0.98470178 938.466174 0.55 508.260122
14:00 30 0.52 -0.16 -0.33 0.9819569 911.969837 0.55 492.533289
15:00 45 0.79 -0.16 -0.33 0.97687323 856.264825 0.55 460.054201
16:00 60 1.05 -0.16 -0.33 0.96954599 744.375951 0.55 396.938694
17:00 75 1.31 -0.16 -0.33 0.95969591 489.6218 0.55 258.438421
91
Diciembre
Tiempo Angulo horario Angulo horario Latitud Io media Eficiencia I efec
Declinación (δ) cos ϴ
(h) (τ seg) (τ rad) (°) (W/𝒎𝟐 ) óptica (ή) (W/𝒎𝟐 )
07:00 75 1.31 -0.16 -0.40 0.93801043 489.6218 0.55 252.598695
08:00 60 1.05 -0.16 -0.40 0.95026911 744.375951 0.55 389.046611
09:00 45 0.79 -0.16 -0.40 0.95966064 856.264825 0.55 451.948009
10:00 30 0.52 -0.16 -0.40 0.96638505 911.969837 0.55 484.722709
11:00 15 0.26 -0.16 -0.40 0.97011775 938.466174 0.55 500.732479
12:00 0 0 -0.16 -0.40 0.97133797 945.306215 0.55 505.016504
13:00 15 0.26 -0.16 -0.40 0.97011775 938.466174 0.55 500.732479
14:00 30 0.52 -0.16 -0.40 0.96638505 911.969837 0.55 484.722709
15:00 45 0.79 -0.16 -0.40 0.95966064 856.264825 0.55 451.948009
16:00 60 1.05 -0.16 -0.40 0.95026911 744.375951 0.55 389.046611
17:00 75 1.31 -0.16 -0.40 0.93801043 489.6218 0.55 252.598695
92
ANEXO 3. CÁLCULOS PARA HALLAR VOS VALORES MAXIMOS DE LA
PARÁBOLA.
𝑏
S = ∫𝑎 √1 + (𝑓´(𝑥))2 dx
𝑋2
Este valor depende del foco “F” que se escoja, pues f(x) = y = , por lo que:
4∗𝐹
𝑋2
F = 0.50 m → f(x) = = f´(x) = x → Xmáx. =0.43
4(0.5)
𝑏
S = ∫𝑎 √1 + (𝑥)2 dx, integrando y luego reemplazando en f(x) se obtiene Xmáx.
𝑋2
F = 0.30 m → f(x) = = f´(x) = 1.67x → Xmáx. =0.39
4.(0.3)
𝑏
S = ∫𝑎 √1 + (1.67𝑥)2 dx, integrando y luego reemplazando en f(x) se obtiene
Xmáx.
𝑋2
F = 0.20 m → f(x) = = f´(x) = 2.5 x → Xmáx. =0.39
4.(0.2)
𝑏
S = ∫𝑎 √1 + (2.5𝑥)2 dx, integrando y luego reemplazando en f(x) se obtiene Xmáx.
93
ANEXO 4: Tabla para hallar las propiedades del agua saturada.
94
ANEXO 5. Mediciones de la radiación en los meses de abril y mayo el 2018.
95
ANEXO 6.
a) Cuadros de datos para obtener el flujo másico (realizado el 17 de abril del 2018).
96
19 12:30 22 69.21 87
20 12:35 22 69.48 83
21 12:40 22 70.04 87
22 12:45 22 72.4 83
23 12:50 22 72.84 86
24 12:55 22 74.04 89
25 13:00 22 74.21 94
PROMEDIO 22 66.69 91
97
Datos tomados el 19 de abril del 2018.
Temperatura de Salida
promedio (°C)
Día 1 Día 2 Día 3
66.69 64.86 67.34
Promedio 66.2967
98
c) Cuadros de datos para obtener la temperatura media logarítmica en el
intercambiador de calor a contra flujo.
Datos tomados el 17 de abril del 2018.
Item Tiempo Tc,entr (°C) Tc,sal (°C) Tf,entr (°C) Tf,sal (°C)
(h)
1 11.00 56.21 37.11 22 38.45
2 11.05 56.32 37.17 22 38.50
3 11.10 56.45 36.35 22 38.71
4 11.15 57.42 36.24 22 39.81
5 11.20 57.52 37.18 22 39.94
6 11.25 58.38 37.28 22 39.54
7 11.30 58.45 37.25 22 39.84
8 11.35 59.21 35.4 22 39.61
9 11.40 59.36 36.96 22 40.01
10 11.45 60.32 37.52 22 40.21
11 11.50 68.7 45.5 22 42.98
12 11.55 69.08 45.58 22 43.15
13 12:00 75.24 46.1 22 43.65
14 12:05 75.26 46.18 22 43.75
15 12:10 74.98 45.98 22 42.14
16 12:15 73.89 45.9 22 43.34
17 12:20 74.21 45.75 22 43.45
18 12:25 74.01 45.61 22 43.00
19 12:30 69.21 44.21 22 43.05
20 12:35 69.48 44.36 22 41.98
21 12:40 70.04 44.41 22 42.45
22 12:45 72.4 44.89 22 42.14
23 12:50 72.84 44.91 22 42.32
24 12:55 74.04 45.71 22 42.36
25 13:00 74.21 45.86 22 42.40
PROMEDIO 66.69 41.3 22 42.12
99
Datos tomados el 18 de abril del 2018.
100
Datos tomados el 19 de abril del 2018.
101