IMPLEMENTACION DE TERMA SOLAR DE 120 LITROS EN
CASAS RURALES DE HUAMACHUCO.
ESTRUCTURA DEL PROYECTO
1. TITULO DEL PROYECYO
2. RESUMEN
3. ANTECEDENTES
4. JUSTIFICACION
5. PROBLEMÁTICA
6. OBJETIVOS
7. MARCO TEORICO
8. VARIABLES DE DISEÑO
9. CALCULOS
10. RESULTADOS
11. CONCLUSIONES
12. REFERENCIAS BIBLIOGRAFICAS
I. ANTECEDENTES
La energía solar es una de las opciones que se están desarrollando como alternativas
a las energías provenientes de la quema de combustibles fósiles. A diferencia de los
países nórdicos, el territorio peruano, por estar mucho más próximo al Ecuador,
cuenta con sol durante la mayor parte del año. Según el Atlas Solar del Perú
elaborado por el Ministerio de Energía y Minas
En la gran mayoría de localidades del Perú, la disponibilidad de la energía solar es
bastante uniforme durante todo el año, estando casi siempre dentro de un margen de
+/- 20 % del promedio anual. es lo suficientemente alta y uniforme (comparada con
otros países) para ser considerada como una fuente energética utilizable para
fomentar el desarrollo de las comunidades
El Perú tiene una elevada radiación solar anual siendo en la sierra de
aproximadamente 5.5 a 6.5 kWh/m2; 5.0 a 6.0 kWh/m2 en la Costa y en la Selva de
aproximadamente 4.5 a 5.0 kWh/m2.
En el Perú hay tres ámbitos donde se ha desarrollado el uso de energía solar en el
Perú. El primer ámbito (y más tradicional) es el uso como fuente térmica a través de
termas de agua en zonas del sur peruano, principalmente Arequipa y Puno,
departamentos en los que existe cerca de 30 empresas dedicadas a la fabricación y
mantenimiento de estos aparatos.
Existen alrededor de 10000 termas solares instaladas en Arequipa, Ayacucho 44,
Lima 25, Puno 52, Tacna 29, Ancash 11 son los principales; de acuerdo a estudio
realizado por el MEM durante el 2004(Fuente: MINIAGRI)
ENERGIA SOLAR TERMICA
1
De ese mismo estudio se obtuvo que se habían instalado 640 cocinas solares en total
(Ancash 242, Junior 79, Cusco 39, Lima 53, Puno 128 son los principales. (Fuente:
MINIAGRI)
Si cada terma solar de uso doméstico transforma 3 kWh de energía solar al día,
5,000 termas transformarán 15 MWh, lo que corresponde (según el uso promedio
diario de instalaciones termoeléctricas) a una planta de 3 MW. (Fuente:
SOLARTEC)
Las termas solares se pueden utilizar con excelentes resultados en toda la costa y
sierra, principalmente en Arequipa, Moquegua y Tacna; en las Regiones Grau, Nor
Oriental del Marañón, Libertad y Chavín. Su rango de operación está entre los 20 –
100 °C. Las capacidades de las termas solares tanto las de fabricación local como las
importadas varían entre 50 – 2000 litros. Para el Perú puede considerarse
dependiendo de la región y de los materiales de construcción el valor de 60– 70
litros / m2 de superficie. (Fuente: SOLARTEC)
Fuente: SENAMH
ENERGIA SOLAR TERMICA
2
Figura 1: Atlas solar del territorio peruano
II. PROBLEMÁTICA
¿Cómo implementar una terma solar de 120 litros para la obtención de agua caliente en
las casas rurales de Huamachuco?
III. OBJETIVOS
1. Objetivo general
Implementación de una terma solar de 120 litros para la obtención de agua
caliente en la sierra libertina,
2. Objetivos específicos
Determinar el área efectiva de absorción de energía radiante en cada colector
solar, la energía radiante recibida en el colector solar, el calor entregado al agua
en cada colector solar.
Determinar la eficiencia instantánea de la terma solar, la eficiencia de la placa
absorbedora, el factor de eficiencia del colector.
Satisfacer una porción considerable del mercado de hogares unifamiliares en
busca de termas eficientes, económicas y ecológicas.
Fomentar una cultura eco amigable que permita sustituir parte del saturado
consumo eléctrico en favor de fuentes de energía renovable, como la solar.
IV. MARCO TEÓRICO
La energía solar, siendo una fuente inagotable y disponible, se presenta como alternativa
para reemplazar a las convencionales. En el Perú, a cada metro cuadrado de superficie
horizontal llega cada día cerca de 5 KW de energía solar; por tanto, no existe problema
en la disponibilidad de esta energía. Comparando este valor con otras formas de energía.
Conversión foto térmica, transforma energía solar en calor útil, es la técnica más simple
y difundida; consiste en colectores solares que absorben la radiación solar en una
superficie negra.
RADIACION SOLAR
La radiación solar es un tipo de energía que se propaga a través del espacio, a la
velocidad de la luz y en forma de ondas electromagnéticas; dicha radiación se distribuye
sobre un rango de longitudes de onda a la que se le reconoce como espectro de radiación
electromagnética. Dentro del espectro, la radiación solar que llega a la superficie de la
tierra está ubicada en el rango de longitud de onda comprendido entre 0,39 a 3,0 um.
ENERGIA SOLAR TERMICA
3
La parte del espectro que va de 0,40 um a 0,78 um forma el espectro visible que
denominamos comúnmente luz. El resto del espectro, que no es visible, lo forman
las radiaciones con longitudes de onda inferiores a 0,4 um denominados radiaciones
ultravioletas (UV) y con longitudes superiores a 0,75 um denominados radiaciones
infrarrojas (IR).
COMPONENTES DE LA RADIACION SOLAR
No toda la radiación solar que incide sobre la atmósfera llega hasta la superficie del
planeta tierra; los rayos solares al entrar en contacto y atravesar el obstáculo que supone
la atmósfera, sufren una serie de fenómenos físicos entre los que se destacan los
siguientes:
Reflexión. Una parte de la radiación se refleja en la parte superior de las nubes (como
en un espejo).
Difusión. Los rayos solares chocan con las moléculas gaseosas y partículas de polo
difundiéndose en todas direcciones, de forma que una parte vuelve al espacio.
Absorción. Una parte de la radiación es absorbida por las moléculas presentes en la
atmósfera.
Debido a los anteriores fenómenos físicos, un objeto situado en la superficie
terrestre recibirá los siguientes tipos de radiación.
a) Radiación directa. Aquellos rayos que alcanzan la superficie terrestre sin ser
desviados. Esta radiación en días nublados es muy pequeña y provoca sombras.
b) Radiación difusa. La que llega a la superficie terrestre después de haber cambiado
varias veces de dirección al atravesar la atmósfera. A mayor nubosidad más importante
es la radiación difusa.
c) Radiación reflejada. La parte de radiación reflejada por el terreno y otros elementos
de la superficie terrestre y que puede ser reabsorbida por otros objetos. Al cociente entre
la radiación reflejada y la incidente en la superficie terrestre se le llama albedo. A la
suma de los tres tipos de radiación recibida se le llama radiación total.
ENERGIA SOLAR TERMICA
4
COLECTOR SOLAR PLANO
En este tipo de colector, la radiación solar incide en una placa plana, por lo general de
un material con alta conductividad térmica como aluminio o cobre debido a que el
rendimiento térmico del colector rectangular depende en gran medida de la velocidad de
transferencia de calor entre la superficie de absorción y el agua (Hussain, 2005).
La placa almacena temporalmente energía térmica y la transfiere a un fluido que circula
en una tubería la cual está unida la placa. La energía almacenada en la placa es
transferida al fluido. Por lo general los colectores solares de placa plana son de aspecto
oscuro para maximizar la absorción del espectro solar (Kalogirou, 2000).
Además de la placa y las tuberías, este tipo de colectores se construyen con una cubierta
superior fabricada con un material transparente. Por lo general se utiliza vidrio de
borosilicato transparente. Este material presenta transmisividad alta a la radiación solar
pero baja a la radiación infrarroja. Esto permite que la mayor parte de los rayos solares
pasen a través de la cubierta transparente e incidan sobre la placa absorvedora de
radiación. Cuando esta última se calienta, emite radiación en la longitud de onda del
infrarrojo la cual no puede escapar en su totalidad debido a la cubierta transparente y
con esto se favorece el efecto invernadero. También se utilizan polímeros transparentes,
pero no son de uso común debido a su menor transmisividad (Mills, 1988).
Otra función de la placa es: cuando la radiación solar pasa a través de una cubierta
transparente y afecta la superficie oscura de la placa colectora, una gran parte de la
energía es absorbida por la placa y luego se transfiere por medio de los tubos al fluido
de transporte para su uso o almacenamiento (Kalogirou, 2000). También disminuye las
corrientes de aire ya que este al estar confinado no es afectado por el viento, y a la
circular más lento disminuye la magnitud de calor que se pierde por convección. Esto
aumenta la temperatura interior y favorece el “efecto invernadero” dentro del colector
como se muestra en la siguiente figura.
Figura 2: Descripción grafica de la función de una placa abservadora de radiación.
ENERGIA SOLAR TERMICA
5
Se debe tener en cuenta que los colectores de placa plana, generalmente son fijados de
forma estable en su posición y no requieren de seguimiento del sol. Los colectores
deben estar orientados directamente hacia el ecuador mirando hacia el sur en caso de
estar en el hemisferio norte, y en caso de estar en el hemisferio sur deben colocarse
mirando hacia el hemisferio norte (Kalogirou, 2000) esto hace que sean de fácil
instalación y funcionamiento. Según Nandwani (2013) un ángulo ideal de instalación
para el colector es de 20 grados con respecto a la horizontal, ya que permite un flujo
eficiente a lo largo de las tuberías, y reduce la cantidad de depósitos acumulados a lo
largo del tiempo. A diferencia de otros, el sistema de placa plana solo Cubierta
Transparente 45 es eficiente durante pocas horas al día en las que la radiación solar
incide de manera más perpendicular a la placa.
Figura 3: Componentes principales de un colector solar de placa plana y su
comportamiento ante la radiación solar.
ENERGIA SOLAR TERMICA
6
Componentes de la terma solar
Los componentes básicos de la terma solar son:
a) Placa de absorción: “Es el elemento encargado de absorber la energía disponible del
sol y transformarla en energía térmica para luego ser transferida al agua, generalmente
está hecha de un metal (cobre, acero, aluminio, etc.).”
b) Cubierta: Generalmente es un vidrio transparente montado en frente de la placa
absorbente, en la parte superior del colector, creando un espacio entre la placa y ella.
“La función de la cubierta es permitir el paso de la radiación solar absorbida por la
placa, igualmente disminuye la cantidad de radiación infrarroja que se escapa al
exterior, disminuyendo de esta forma las pérdidas del colector.” (Burbano et al., 2006,
p.86)
c) Conductos para la circulación del fluido: El colector solar de placa plana tiene que
poseer una serie de conductos por los cuales circula el fluido de trabajo (agua), el cual
recibe y transporta la energía absorbida por la placa hacia el tanque de almacenamiento.
(Orozco, 1993)
d) Aislante térmico: Para Burbano et al. (2006) es el punto básico para disminuir las
pérdidas de calor por conducción en la parte inferior y lateral del colector. Las
características que debe poseer el material utilizado para ser un buen aislante son: baja
conductividad térmica (menor de 0,040 W/m°C en el rango de 20 a 120 °C) y no debe
absorber agua.
e) Caja: “La caja es el elemento que soporta todos los componentes del colector, la cual
impide que la humedad, polvo y aire penetren por el colector y disminuyan su
eficiencia.” (Burbano et al., 2006, p.86)
f) Tanque de almacenamiento: Manrique (1984) menciona que el calor puede ser
almacenado en el tanque por circulación directa entre el colector y el tanque o el agua
calentada en el colector puede circular por un serpentín dentro del tanque, transfiriendo
así calor al agua que se encuentra dentro del mismo.
El mayor requerimiento técnico del tanque es su total aislamiento, con el fin de
conservar la temperatura requerida en el agua. La ubicación del tanque es importante, ya
que de ésta, depende la eficiencia del termosifón y el evitar el fenómeno de flujo
inverso.
ENERGIA SOLAR TERMICA
7
Figura 4: Flujo del agua a través de una terma solar
Radiación Solar
El Sol es una fuente de energía aprovechable, de ahí la importancia de
determinar la distribución solar horaria y la magnitud de esta energía disponible
ya que es utilizada en muchas aplicaciones en el campo de la energía solar
(Siqueira et al, 2011). A lo largo del presente capítulo se desarrollará la teoría
para la determinación de la radiación solar teórica, y posteriormente se expondrá
el método del día claro de Hottel (1976). Este método permite determinar la
radiación solar dependiendo de la latitud y la altura de un lugar específico, en
condiciones de atmosfera clara, con la ventaja de ser un método sencillo y su
utilidad se ve potenciada para simulaciones donde no se dispone de instrumentos
adecuados para la medición de la radiación solar (Passamai, 2000).
Constante solar
Se define como el flujo de energía solar interceptada por una superficie normal
a la dirección del Sol en Watts por metro cuadrado (W/m2 ). Dicha constante
representa la cantidad de energía extraterrestre procedente del Sol. En la
ENERGIA SOLAR TERMICA
8
superficie terrestre el valor es menor debido a que sufre diferentes procesos de
atenuación, como las nubes, la capa de ozono y la reflexión de las ondas en la
atmosfera entre otros factores (Çengel, 2007). El valor de la constante solar es
de 1 366 W/m2 con una incertidumbre de ± 1,6 W/m2 y una desviación máxima
de + 7 W/m2 (García Cole, 2012).
PARÁMETROS GEOMÉTRICOS
Para el presente estudio y también para poder abordar el cálculo de la radiación sobre la
superficie en la Tierra, varios parámetros relevantes a la radiación incidente deben ser
definidos. Esta sección está dedicada a su estudio
Ángulo acimutal βs
Este ángulo es el que forma la componente horizontal de la radiación solar directa con
el eje de rotación de la Tierra como se muestra en la Figura 5, este se mide en el plano
horizontal que pasa por el lugar de medición (Fernández, 2005). Es medido hacia el este
positivamente y oeste negativamente, (sur = cero). Este varía entre 0° y 180°.
Figura 5: Parámetros geométricos para el estudio de la radiación solar en la
tierra.
Declinación solar δs Indica la posición angular del Sol al mediodía con respecto al
plano del ecuador. Este equivale a la latitud proyectada sobre la esfera celeste
(Fernández, 2005). Este ángulo varía con las estaciones y su período es de un año, es
decir, el tiempo que tarda la Tierra en dar una vuelta completa alrededor del Sol. En la
Figura 6 se muestra el ángulo de declinación solar y cómo afecta la incidencia de la
radiación en los dos hemisferios terrestre a lo largo de un año solar.
ENERGIA SOLAR TERMICA
9
Figura 6: Declinación solar.
El valor de la declinación varíe entre ± 23°45´ a lo largo del año. La declinación
al Norte del Ecuador se considera positiva. Para calcular este valor se pueden
utilizar las siguientes ecuaciones (Vera, 2005). Fórmula de Perrin de
Brichambaut.
360
δ=sin−1 0,4 sin d −82 (1)
365 j
Fórmula de Cooper
360
δ =23,45 sin d +2824 (2)
365 j
d j : Día juliano en que se desea realizar el cálculo (1 < d j < 365)
(ejemplo :2 de febrero correspondería a d j=33)
ENERGIA SOLAR TERMICA
10
Ángulo horario ω
Este ángulo hace referencia al desplazamiento horario del sol con respecto al mediodía.
Cada hora es igual a 15° de longitud angular, tomando un valor positivo para horas
antes del mediodía y por las tardes toma un valor negativo.
Por último, es importante mencionar dos conceptos como lo son la latitud λ que indica
la posición con respecto al ecuador, donde los lugares en el hemisferio norte reciben
valores positivos y valores negativos en el hemisferio sur. La longitud del lugar L es el
ángulo que forma el plano que pasa por los polos en un lugar de la superficie de la
Tierra y el plano que pasa por el meridiano de Greenwich.
Cálculo de la radiación
A continuación, se describen las diferentes ecuaciones que modelan el cálculo teórico de
la radiación sobre un lugar específico.
Distancia de la Tierra al Sol
Es bien conocido que la órbita de la Tierra alrededor del Sol no es redonda si no que es
más bien una órbita elíptica, por lo que el radio de la órbita varía dependiendo en la
época del año. Al ser su órbita una elipse, existe un semieje mayor el cual es la mayor
distancia entre el Sol y la Tierra y un semieje menor, el cual corresponde a la menor
distancia entre el Sol y la Tierra (Vera, 2005). Se define también la unidad astronómica
(UA) como la distancia media entre el Sol y la Tierra. Esta tiene un valor de 1,496x108
km, por lo que la distancia entre el Sol y la Tierra tendría un valor mínimo de 0,983 UA
y un valor máximo de 1,017 UA (Vera, 2005).
La importancia de esta información es dejar claro el hecho de que el radio de la órbita
terrestre varía conforme la Tierra se mueve alrededor del Sol, por lo que para el cálculo
de la radiación, se debe introducir un factor de corrección para el radio de rotación
terrestre. Para este cálculo Duffie & Beckman (1974) propusieron la siguiente ecuación:
r 2 2 πd j
2
ρ=
ro( )
=1+0.033 cos
365
(3)
2
ρ : Factor de correció n debidoa la traslaciónterrestre .
ENERGIA SOLAR TERMICA
11
Relaciones trigonométricas
El colector solar se debe colocar a una cierta inclinación, 20 grados con respecto a la
horizontal según Nandwani (2013), por lo que se deben ajustar los valores de radiación
para la proyección horizontal del área del colector la cual recibe la radiación solar. En
una superficie las relaciones trigonométricas para un plano en cualquier posición, se
plantean de la siguiente forma:
cos θ s=sinδ sinλ cosψ + sinδ cosλ sinψ cos β s +cosδ cosλ cosψ cosω+cosδ sinλ sinψ cos β s cosω+ cosδ sinψ s
Donde:
ψ : Águlo que forma el terrenocon la horizontal .
λ : Latitud de la zona a evaluar .
β s : Águlo acimutal .
Para simplificar, se puede suponer que el plano del terreno coincide con el plano
horizontal, dicho de otra forma ψ = 0, los ángulos cenital s y el acimutal βs, de la
radiación directa se definen con las siguientes ecuaciones:
cos θ s=sinδ sinλ+cosδ cosλ cosω (5)
tanδ cosλ−sinλ cosω
cot β s= (6)
sinω
RADIACIÓN SOLAR DIRECTA
Para encontrar la ecuación que permita calcular la cantidad de radiación sobre una
superficie inclinada, se procede a utilizar el “Método de Hottel” (Hottel, 1976). Este
método permite calcular la radiación para condiciones de un día claro utilizando
parámetros geométricos, geográficos y características climáticas del sitio para el cual se
desea calcular la radiación. Además de ser un método sencillo en comparación con otros
métodos presentados en textos clásicos (Passamai, 2000).
ENERGIA SOLAR TERMICA
12
Este modelo define la transmitancia directa de la radiación (porcentaje de la radiación
solar extraterrestre que se ve absorbida por la atmósfera terrestre) usando un modelo de
gas negro-gris-claro que arroja la siguiente relación.
−k /cosθ
τ b =a0 + a1 e (7)
τ b :trasmitancia directa de la radeacion solar .
Los parámetros a0 , a1 y k adimensionales deben ser ajustados dado que la altitud
a la que se encuentra el punto donde se desea conocer la radiación introduce
irremediablemente un error, mismo que se corrige usando las siguientes ecuaciones
cuadráticas para una visibilidad de 23 km en un día claro (Hottel, 1976).
¿ 2
a0 =0,4237−0,008216− A (8)
¿ 2
a1=0,5055+0,00595 6,5− A (9)
k ¿ =0,2711+ 0,018582,5− A2 (10)
A : Altura sore el nivel del mar .
Finalmente, estos valores deben ser corregidos nuevamente para mostrar el efecto que
produce el tipo de clima sobre la cantidad de radiación incidente. Definiendo las
¿ ¿ ¿
siguientes proporciones: r 0=a0 a0 , r 1=a1 a 1 y r k =k k , donde los valores de
r0 , r1 y rk se resumen en la Tabla 1 para distintos tipos de clima.
Tabla 1. Factores de corrección de radiación solar directa para algunos climas.
Tipo de clima r0 r1 rk
Tropical 0,95 0,98 1,02
Verano, latitud media 0,97 0,99 1,02
Verano, sub-ártico 0,99 0,99 1,01
Invierno, latitud media 1,03 1,01 1,00
Fuente: (Hottel, 1976).
Nota: Donde los valores son para un clima estándar tropical, templado (invierno y
verano) y sub ártico.
ENERGIA SOLAR TERMICA
13
Combinando las ecuaciones 8, 9 y 10 con las correcciones presentadas en la Tabla 1
para un clima tropical y despejando los valores necesarios para determinar la
transmitancia en la ecuación de Hottel, se tienen las siguientes relaciones:
2
a0 =r 0 0,4237−0,008216− A ( 11)
2
a1=r 1 0,5055+ 0,005956,5− A (12)
2
k =r k 0,2711+0,01858 2,5− A (13)
Con esto ya se puede determinar la transmitancia solar directa utilizando la ecuación 7.
Ahora, para obtener la radiación incidente sobre la superficie terrestre, es necesario
multiplicar la radiación extraterrestre definida como constante solar (1 366 W/m2 ), por
el valor de transmitancia obtenido anteriormente.
Como se mencionó anteriormente, la órbita de la Tierra es elíptica y la distancia entre el
Sol y la Tierra varía a lo largo del año, por lo que el valor de radiación se verá afectado
por esta distancia. Se debe utilizar la ecuación 3 para corregir el valor de la radiación
(Passamai, 2000). Con esto, la radiación solar directa sobre el plano perpendicular a los
rayos solares en la superficie terrestre viene dada por:
Gb=τ b∗G∗ρ2 cos θs (14)
Donde:
Gb : Radiación solar directa sobre la superficie terrestre(W m−2).
τ b :Transmitancia directa de la radiación solar .
ρ2 : Factor de correccióndebido a latraslación terrestre .
cos θ s :Componente perpendicular a lasuperficie .
G: Radiación solar extraterrestre ( 1 366 W m−2 ) .
3.3.4. Radiación solar difusa.
Duffie & Beckaman (1974) definen la componente difusa de la radiación como “la
radiación solar recibida desde el Sol luego que su dirección ha sido cambiada por
reflexión y dispersión en la atmosfera”. Para determinar el valor de la componente
difusa de la radiación se utilizará la ecuación propuesta por Liu & Jordán (1960), la cual
define un valor de transmitancia difusa para una atmosfera clara basándose en curvas
ENERGIA SOLAR TERMICA
14
teóricas y experimentales y utilizando el método de los menores cuadrados obtienen la
siguiente ecuación:
τ d =0,2710+ 0,2939 τ b (15)
Donde
τ b :Transmitancia difusa de la radiación solar .
De forma similar a la ecuación 3.14, la componente difusa de la radiación finalmente
viene dada por:
2
Gd =τ d∗G∗ρ cos θ s (16)
Albedo
El albedo es la fracción de la radiación, que proviene luego de reflejarse contra las
paredes y otras superficies. Una superficie inclinada ve una porción del suelo, esta luz
que es reflejada por las diferentes superficies, la capacidad de las mismas de reflejar
esta energía se denomina coeficiente de albedo . En la Tabla 2 se presentan los valores
del porcentaje de albedo para algunas superficies.
Tabla 2. Abeldo de algunas superficies comunes
Fuente:
(Fernandez,2005).
El albedo es una fracción de la radicación
global o total que llega a la tierra
q rad =Gd −Gb por lo que la radiación
que recibe una superficie inclinada debido
al albedo estará dada por la siguiente
expresión :
Donde:
Ga : Radiacion solar debido al
albedo ( W . m2 )
ENERGIA SOLAR TERMICA
15
q rad : Radiación solar total. Componente difusa más componente directa (
2
W .m )
s : es el ángulo que forma el plano con la superficie (grados)
a : Albedo
Con estas ecuaciones se puede calcular la radiación solar incidente sobre una
placa horizontal, pero como el colector propuesto y la mayoría de los mismos no se
instalan en dicha posición, hay que corregir estos valores debido a la posición del
colector con respecto a la horizontal, para la componente directa Duffie & Beckaman
(1974) proponen la siguiente relación geométrica:
cos θT
Rb= (17)
cos θ S
Rb : Relación angular de la posición del colector con respecto a la horizontal para la
componente directa de la radiación.
Para la componente difusa la situación es un poco más complicada, dado que esta
depende de la distribución de la radiación sobre el cielo y esta misma normalmente no
se conoce ya que se distribuye en todas direcciones y los métodos para determinarla son
complejos (Duffie & Beckman, 1974). Debido a la dificultad para determinar una
relación sencilla como en el caso de la componente directa, Duffie & Beckaman (1974)
proponen las siguientes suposiciones para determinar la corrección angular para la
componente difusa:
Se asume que la mayoría de la radiación difusa que llega al colector viene de un
punto común muy cercano a la posición del Sol, esto es que la dispersión de la
radiación es una dispersión fundamentalmente orientada Esta aproximación es
mejor en días claros.
Se asume que la radiación difusa está uniformemente distribuida en el
firmamento y además es bastante mayor que la componente reflejada por otras
superficies (albedo) por lo que se desprecia la corrección para la radiación
difusa.
Bajo estas suposiciones, el valor total de la radiación incidente sobre un colector solar
de placa plana inclinado viene dado por la suma de la componente difusa más la
componente directa corregida por la inclinación.
q rad =Rb Gb +Gd (18)
ENERGIA SOLAR TERMICA
16
Factor de remocion de calor FR
Conociendo la distribución de la temperatura del fluido a lo largo de las tuberías, se
F
define el factor de remoción de calor (¿¿ R) como el parámetro que relaciona el
¿
calor de calentamiento que se obtiene efectivamente en el colector, con el que se
obtendría si toda la superficie del colector se encontrara a la temperatura del fluido a la
T
entrada (¿¿ fe) .
¿
τα . q rad −U T . T fe−T a
¿
A¿
mC p (T fs −T fe )
F R= ¿
A : Área efectiva del colector m2
T fe : Temperatura del fluido de salida (°C)
m Cp τα q rad ( U T −T fs −T a )
F R= 1− (20)
A UT τα q rad ( U T −T fe −T a )
Sustituyendo, tenemos los siguiente
m Cp −U T MnL
F R= 1−exp (21)
A UT m Cp
En consecuencia, el calor útil en términos del factor de remoción de calor es como
sigue:
qutil = Ae F R τα qrad −U T T fe−T a (22)
ENERGIA SOLAR TERMICA
17
CIRCULACIÓN DEL AGUA POR CONVECCIÓN NATURAL O EFECTO
BOYANTE EN EL COLECTOR SOLAR Y POSICIONAMIENTO DEL TANQUE
DE ALMACENAMIENTO RESPECTO AL COLECTOR
Convección natural
La circulación natural o circulación por termosifón es un fenómeno que se fundamenta
en la variación de la densidad del agua, que es inversamente proporcional al incremento
de la temperatura, esto quiere decir por ejemplo, que la densidad del agua a 60 ºC es
menor que la densidad de agua a 20 ºC.
Por lo tanto para una tubería con un diámetro constante, el incremento de flujo se debe a
la absorción de calor de este ya que se afecta la densidad del fluido con la temperatura
(Quiros, 1979). Existen dos fuerzas principales que afectan el flujo dentro de las
tuberías y son: La fuerza boyante en función del peso específico del fluido y la fuerza de
fricción contraria a la fuerza boyante y que se ve afectada por el diámetro de la tubería,
la velocidad de flujo, la viscosidad cinemática y la rugosidad de la tubería (Taherian et
al, 2011).
Debido a que en este tipo de colector solar aprovecha la circulación natural de agua, no
se necesita de una bomba para la trasegar el fluido desde el colector hasta el tanque de
almacenamiento. Un aspecto importante de la circulación natural es la ubicación del
depósito de almacenamiento respecto del colector solar como se observa en la Figura
24. En la Figura 7 se muestra la distribución típica del agua dentro del colector.
Figura 7. Distribución de temperaturas del agua para la convección natural.
4.5.2. Posición del tanque respecto al colector.
La geometría de la distribución del tanque y el colector para la convección natural con
las respectivas entradas y salidas de agua, es la que se presenta a continuación.
ENERGIA SOLAR TERMICA
18
Figura 8. Posición del tanque de almacenamiento respecto al colector.
M: Distancia horizontal desde el colector hasta el tanque de almacenamiento (m).
L: longitud del colector (m).
h1 : Altura desde la entrada al colector hasta la parte superior del tanque (m).
h2 : Altura del tanque de almacenamiento (m).
h3 : Altura entre el colector y el tanque de almacenamiento (m).
h4 : Altura desde la entrada y salida del colector (m).
La relación de ángulos necesarios para el análisis se basa en la siguiente Figura 24.
Figura 9. Ángulos y distancias asociadas a la posición del tanque de
almacenamiento con respecto al colector.
ENERGIA SOLAR TERMICA
19
h2+ h3
α 2=tan −1
M
h1
α 1=tan −1
L cos β + M
Altura de presión de la circulación natural
La altura de presión de circulación natural ∆ PCN es la altura resultante del
diferencial de presión generada en virtud de la diferencia de densidades del agua ∆ P
y la pérdida de carga por fricción en los tubos h p .
∆ PCN =∆ P−h p
h p : Pérdida de carga dentro de las tuberías debido a la fricción del fluido con la
tubería (m.c.a.).
∆ P : Deferencia de presiones debido a la diferencia de densidades dentro del
colector (m.c.a.).
∆ PCN : Presión de circulación natural medida en metros de columna de agua (m.c.a.).
La ecuación 4.50 se cumple para el estado transitorio, es decir, cuando el fluido está en
aceleración. Es de interés práctico analizar el caso para estado estacionario o cuando el
fluido está en reposo o a velocidad constante. Para este caso, la diferencia de presión
generada por la diferencia de densidades es igual a las pérdidas en todo el circuito. Es
decir cuando ∆ PCN =0 . Por lo que la ecuación 4.50 toma la siguiente forma.
∆ P=h p
Para calcular la presión generada por el cambio de densidad del agua, el calentador solar
se representa como se muestra en la Figura 25. Representación de columnas de fluido
para el circuito de la Figura 23 y las densidades del fluido asociadas.
Figura 10. Representación de columnas de fluido para el circuito de la Figura 23 y
las densidades del fluido asociadas.
ENERGIA SOLAR TERMICA
20
En la entrada del colector el agua tiene una densidad, mientras que a la salida el agua
tiene una densidad, además se asume una densidad promedio ρ. Aplicando la ecuación
de Bernoulli para estado estacionario se obtiene:
∆ P∗ρs= ρh2 + ρ e h3 +h 4−ρh4 −ρs h3 +h 2( 23)
ρe : Densidad del agua a la entrada del colector ( Kg m3 ).
ρs : Densidad del agua a la salida del colector ( Kg m3 ).
ρ : Densidad promedio del agua dentro del colector ( Kg m3 ).
Además:
ρ s+ ρe
ρ= (24)
2
Reemplazando la ecuación 4.53 en la 4.52:
h2 +h4
ρe − ρs +h3
2
∆ P= (25)
ρs
INCLINACION DEL COLECTOR SOLAR
Según Nandwani (2013) un ángulo ideal de instalación para el colector es de 20 grados
con respecto a la horizontal, ya que permite un flujo eficiente a lo largo de las tuberías,
y reduce la cantidad de depósitos acumulados a lo largo del tiempo. A diferencia de
otros, el sistema de placa plana solo Cubierta Transparente 45º es eficiente durante
pocas horas al día en las que la radiación solar incide de manera más perpendicular a la
placa. Guinebault & Rozis (1997), mencionan que el ángulo de inclinación y la
orientación del colector cumplen un rol fundamental para la correcta captación de la
energía solar. En general, se sigue las recomendaciones señaladas por la siguiente tabla
y se muestra a continuación.
Tabla 3: Orientación del colector respecto a la latitud del lugar para una mejor
captación de la radiación solar.
Latitud del lugar Inclinación
De 0° a 10° i = 10°
De 10° a 20° i = lat.
De 20° a 35° i = lat.+10°
> de 35° i = lat.+15°
ENERGIA SOLAR TERMICA
21
Balance de energía en la terma solar
El balance de energía tiene como objetivo determinar el calor útil transmitido al agua
para establecer la eficiencia de la terma solar.
El calor útil � u𝑡𝑖�, empleado para calentar el agua del depósito de almacenamiento, es
una parte de la radiación global incidente sobre el colector solar 𝐻𝑇𝐴𝑐�; la otra parte es
calor disipado al ambiente
Por el colector solar �𝑐�.
Por el depósito de almacenamiento �𝑑�.
Por las conexiones �𝑐�ε
Figura 11: Balance de energía en el calentador solar.
Balance de energía en el colector solar
La radiación global que incide sobre el colector solar 𝐻𝑇𝐴𝑐�, disminuye su valor como
resultado de las características propias de reflexión y absorción de la cubierta del
colector; este valor, es corregido por un factor (𝜏�).
Esta radiación corregida 𝐻𝑇𝐴𝑐�(𝜏�) es empleada como calor para calentar el agua �𝑐��;
la parte restante constituye el almacenamiento de energía en el colector mismo 𝑑�/𝑑𝑡
que usualmente es despreciable; así como, las pérdidas de calor hacia el ambiente por
conducción, convección y re-radiación �𝑝�.
ENERGIA SOLAR TERMICA
22
dU
H T A cs ( τα )=q cal +q perdidas + (26)
dt
Donde:
A cs : Área efectiva del colector
q cal : Calor útil transportado al fluido
q perdidas : Perdidas que experimenta el colector
dU
: Cambio de energía interna almacenada en el colector
dt
dU
=0
dt
Además
q cs ( τα )=q perd
Para determinar el valor del calor de calentamiento del agua en el colector solar es
primeramente debemos hallar las pérdidas de calor en el colector mediante la siguiente
ecuación:
q perd =U L A cs ( T P−T a ) (27)
Entonces el calor útil quedaría expresado de la siguiente manera
q cal= Acs ( H T ( τα ) −U L ( T P−T a ) ) (28)
Donde
H T : Radiación solar sobre la superficie horizontal
τα : Producto transmitancia-absortancia
U L : Coeficiente global de perdidas
T P Temperatura promedio de la placa del colector
T a :Temperatura ambiente
Bien, podemos ingresar el factor F’ en la ecuación anterior, obteniendo lo siguiente
' A cs
q cal=F ( H T ( τα ) −U L ( T F−T a ) ) (29)
ENERGIA SOLAR TERMICA
23
Donde:
F’: Factor de eficiencia del colector
T F : Temperatura promedio del fluido
También se puede expresar el calor de calentamiento de agua de la siguiente
manera:
q cal=F ' A ( H T ( τα ) −U L ( T Fi−T a ) ) (30)
cs
Donde :
FR: Factor de remoción de calor
Tfi : Temperatura del fluido al ingreso del colector
Coeficiente global de perdidas
Para determinar el coeficiente global de transferencia de calor UL se requiere determinar
las conductancias del inferior, superior, lateral en el colector solar.
En la siguiente figura se muestra el circuito térmico de un colector solar plano que nos
ayuda a observar los fenómenos físicos que se producen en el colector solar y de apartir
de este circuito plantear las ecuaciones de las conductancias.
Figura 12: Circuito térmico de un colector solar plano.
Conductancia en la parte inferior
ENERGIA SOLAR TERMICA
24
La conductancia se evalúa conociendo la resistencia térmica conductiva �1 y la
resistencia térmica convectiva �2.
1
U fo = (31)
R1 + R 2
Donde �1 está dado por la conductividad del aislante 𝐾�, y la longitud transversal del
mismo �; �2 está dado por el coeficiente convectivo de transferencia de calor entre la
placa de absorción y el ambiente, el cual en colectores bien aislados es despreciable.
l
R1= (32)
Ka
l
R2= (34)
hfo
Si �2 << �1 , entonces �2 se desprecia; por lo que la conductancia en la parte inferior
sería:
1
U fo = (35)
R1
Conductancia en la parte superior
Esta conductancia se evalúa determinando las resistencias térmicas �3 , �4 y �5. El
calor se transfiere entre la placa de absorción y la cubierta interior por convección y
radiación en forma paralela; esto mismo, también ocurre entre las dos cubiertas, interior
y exterior. Los mismos mecanismos de transferencia son igualmente importantes en la
disipación de calor por parte de la cubierta exterior hacia el ambiente.
1
U ¿= (36)
R 3+ R 4 + R 5
R3 es la resistencia térmica entre la placa de absorción y la cubierta interior
El valor de R3 está determinado por la siguiente ecuación
T P −T C 2
R3= (37)
(hc 2+ hR 2 )(T P−T C 1)
Donde:
0.310
( T P−T C 2)
h́C 2={1−0.0018 ( T P −10 } 1.14 0.070
(38)
l
ENERGIA SOLAR TERMICA
25
2 2
σ (T P +T C 2 )(T P + T C 2 )
h R 2= (39)
1 1
+ −1
εP εP
Donde:
�: Constante de Boltzman 5.67�10−8 �/(�2𝐾4)
𝑇𝑝: Temperatura promedio de la placa de absorción
: Temperatura promedio de la cubierta exterior
T C 2 : Temperatura promedio de la cubierta interior
h́C 2 : Coeficiente de transferencia de calor por convección entre la placa de absorción
y la cubierta interior
ℎ�2: Coeficiente de transferencia de calor por radiación entre la placa de absorción y la
cubierta interior
�: Distancia entre placa y cubierta
𝜀�2: Emitancia de la cubierta interior
𝜀𝑝: Emitancia de la placa de absorción en el infrarrojo
R4, la resistencia térmica entre la cubierta exterior y la cubierta interior podemos
hallarla mediante la siguiente expresión
T C 2−T C 1
R4 = (40)
( h´c 1 +h R 1)(T C 2−T C 1)
2 2
σ (T C 1+ T C 2 )(T C 1 +T C2 )
h R 2= ( 41)
1 1
+ −1
ε C 2 ε C1
0.310
( T C 2−T C 1 )
h́C 1={1−0.0018 ( T C 2−10 } 1.14 0.070
( 42)
l
𝜀�1: Emitancia de la cubierta exterior
h́C 1 : Coeficiente de transferencia de calor por convección entre las cubiertas exterior
e interior
h R 1 : Coeficiente de transferencia de calor por radiación entre las cubiertas exterior e
interior
ENERGIA SOLAR TERMICA
26
�: Distancia entre cubiertas
R5: Resistencia térmica entre la cubierta exterior y el ambiente
T C 1−T a
R5= (43)
( h´c ∞ +hfir )(T C 1−T a )
Donde :
h´c ∞=5.7 +3.8 V
hfir =σ ε C 1 ( T C 1 +T a ) ( T C 12 +T a2 ) (44)
Ta. Temperatura ambiental
h´c ∞ : Coeficiente de transferencia de calor por convección entre la cubierta exterior y
el ambiente
hfir : Coeficiente de transferencia de calor por radiación entre la cubierta exterior y el
ambiente.
El proceso de cálculo de �3, �4 y �5 es iterativo, para lo cual se utiliza las siguientes
ecuaciones
U ¿ ( T p−T a )
T C 2 =T P− (45)
h́C 1 + h́C 2
U ¿ ( T p−T a )
T C 1=T C 2− (46)
h́C 1 + h́C 2
Conductancia en la parte lateral
Las pérdidas de calor a través de los lados del colector se evalúan de manera similar,
aplicando la siguiente ecuación
k ' a hP
U l= ( 47)
l ' AC
Dónde:
�: Perímetro del colector solar
ℎ: Altura del colector solar
�’: Espesor del aislante por los lados
𝑘�’: Conductividad del aislante lateral
ENERGIA SOLAR TERMICA
27
Finalmente, la evaluación del coeficiente total de transferencia de calor para todo el
colector está representada por la sumatoria de las tres conductancias indicadas
anteriormente
U L=U fo +U ¿ +U l ( 48)
Convección entre la superficie y un fluido
El proceso de transferencia de calor por convección entre una superficie y un fluido, se
basa en la aplicación de la Ley de Newton y las formulas experimentales de “Hauzhten
y Mc Demortt” en donde se corrige la velocidad de viento tomando en cuenta la
rugosidad de la superficie de la pared.
Coeficiente de película:
El coeficiente de película o de convección, representado como �, cuantifica la influencia
de las propiedades del fluido, de la superficie y del flujo cuando se produce
transferencia de calor por convección
Figura 13: Coeficiente de película en una superficie plana
La transferencia de calor entre la superficie de un sólido en contacto con el fluido se da
por combinación de conducción y transporte de masa, por lo tanto, establecemos si el
movimiento de fluidos es régimen laminar o turbulento.
En el flujo laminar el fluido se mueve en capas, cada partícula de este sigue una
trayectoria uniforme y continua. Contrario a este flujo, el movimiento de partículas del
flujo turbulento tiene presencia de remolinos y no siguen una trayectoria uniforme
debido a su desorden.
Cuando el fluido fluye con movimiento laminar sobre una superficie de temperatura
diferente, el calor se transfiere únicamente por conducción molecular. En el flujo
turbulento la conducción sufre una modificación debido a los remolinos que acarrean
masas de fluido, por lo que se incrementará la rapidez de flujo de calor por convección.
El valor del coeficiente f se incrementa al aumentar la rugosidad de las paredes y crece
también con la velocidad del fluido. Sin embargo, existen fórmulas experimentales,
ENERGIA SOLAR TERMICA
28
debidas a “Hauzhten y Mc Demortt” en donde se corrige por velocidad del viento
tomando en cuenta la rugosidad de la pared.
Para superficies muy lisas f =1.4+ 0.28V
Para madera y yeso f =1.6+0.3 V
Para superficies rugosas f =2.1+0.5 V
Para aire quieto f =1.65
Donde V esta dado en unidades de milla por hora
V. VARIABLES DE DISEÑO DEL COLECTOR SOLAR PLANO
Tabla 4: Medidas correspondientes a las respectivas características del colector.
CARACTERISTICAS DEL MEDIDAS
COLECTOR
ANCHO Bc=1.00 m
LONGITUD O LARGO Lc=1.50 m
2
ÁREA DEL COLECTOR Ac=1.5 m
ALTURA DEL COLECTOR Zc=0.145m
ESPESOR DEL AISLANTE Za=0.04 m
ESPACIAMIENTO ENTRE PLACA Y l=0.0254 m
CUBIERTA
DIÁMETRO EXTERIOR DEL TUBO 1
De= pulg=12.70 mm
2
DIÁMETRO INTERIOR DEL TUBO Di=11.1 mm
DISTANCIA ENTRE TUBOS W =0.100 m
ESPESOR DE PLACA g=0.508 mm
ESPESOR DE TUBO e t=7.62 ×10−4 mm
ANGULO DE INCLINACIÓN ^
S=10 °
Tabla 5: Otros datos a tener en cuenta para el diseño de la terma.
OTROS DATOS PARA EL DISEÑO MEDIDAS
Emisividad de la placa en infrarrojo Gp=0.90=α p
(pintura negra)
Emisividad del vidrio en infrarrojo ε v =0.88
Conductividad térmica del aislante (… Ka=0.04 W /mK
de vidrio)
Conductividad térmica del cobre K=386 W / mK
Temperatura ambiente Ta=11.42°=284.42 K
Temperatura promedio de placa (valor Tp=60 ° =333 K
asumido)
Velocidad del viento v =7.7 Km/h
VI. CALCULOS, ANALISIS Y DESARROLLO
ENERGIA SOLAR TERMICA
29
A. DETERMINACIÓN DEL ÁNGULO DE INCLINACIÓN:
Latitud: 7° 47’ 59”s
Longitud: 78° 0.4’ 00” 0
Altitud sobre el nivel del mar: 3341 m.s.n.m
Del criterio de Guinebault y Rozis, la inclinación del colector solar será de i=10°
Tabla 6: Tabla climática del tiempo de Humachuco.
B. TEMPERATURAS PROMEDIOS
12.4 +11.9+12+11.7 +11+10.3+10.4+ 11.1+11.9+11.9 +12
Ta=Tm=
12
Ta=11.42 ° C=284.42 K
19+18.3+18.5+18.5+18.9+19.4 +19.7+19.4+19.5+ 19.4+19.5+19.5
Tmax=
12
Tmax=19.13 ° C=292.13 K
5.6+5.6+5.6 +5+3.1+1.2+1.1+ 1.4+2.8+ 4.4+ 4.3+4.6
Tmax=
12
Tmin=3.75 ° C=276.75 K
C. CÁLCULO DEL COEFICIENTE GLOBAL DE PÉRDIDAS (U L )
ENERGIA SOLAR TERMICA
30
Asumiendo inicialmente la temperatura del vidrio igual a 40 °C tenemos que
Tv=313 K
Coeficiente de transferencia de calor por radiación entre la placa y el vidrio :
2 2
σ (T p+ T v )(T p +T v )
hr p−v =
1 1
+ −1
ε p εv
−8 2 2
5.67 ×10 (333+313)(333 +313 )
hr p−v =
1 1
+ −1
0.90 0.88
hr p−v =6.132W /m2 K
2. Coeficiente de transferencia de calor por radiación entre el vidrio y el amiente.
hr =σε (T v +T a)(T 2v +T 2a)
v−a
hr =(5.67∗10−8 )(0.88)(333+292.13)(3332 +292.132 )
v−a
W
hr =6.121 2
m k
v−a
3. Coeficiente de transferencia de calor por convección entre la placa y el vidrio.
h́ ¿ Δ T 0.310
=1−0.0018 ( T́ −10 ) , pero h́ =1.14
h́¿ l 0.070
0.310
Δτ
h́ p−v =[1−0.0018 ( T́ −10 ) ]1.14 0.070
l
W
h́ p−v =2.460 2
m k
4. Coeficiente de transferencia de calor por convección entre el vidrio y el
ambiente.
h́v−a =5.7+3.8 V m
9.8 ( 3.6 ) +6.9 ( 8.4 )
V m=
12
ENERGIA SOLAR TERMICA
31
h́v−a =5.7+3.8(3.158)
km m
V m=7.77 =2.158
h s
W
h́v−a =13.90
m2 k
5. Coeficiente de transferencia de calor superior.
−1
1 1
U ¿=
( hr p−v
+
+ h́ p −v hr + h́v−a
v− p
)
−1
1 1
U ¿=( +
6.132+2.460 6.121+13.90 )
W
U ¿ =6.012
m2 k
6. Verificado la temperatura resumida pero al vidrio.
U ¿ (T p−T a )
T v =T p−
h́ p−v +hr p−v
6.012(60−19.13)
T v =60−
2.460+6.132
T v =31.40° C
7. Calcular el U ¿ para el nuevo valor obtenido T v =31.40° C
T v =31.40° C=304.4 K
Coeficiente de transferencia de calor por radiación entre la placa y el
vidrio.
σ ( T p +T v )(T 2v +T 2a )
hr =
p−v
1 1
+ −1
ε p εv
−8 2 2
(5.67∗10 )(333+304.4)(333 +304.4 )
hr = p−v
1 1
+ −1
0.90 0.88
ENERGIA SOLAR TERMICA
32
W
hr =5.90
m2 k
p−v
Coeficiente de transferencia de calor por convección entre la placa y el vidrio.
ΔT 0.310
h́ p−v =[1−0.0018 ( T p−10 ) ]1.14
l 0.070
(60−31.40)0.310
h́ p−v =[1−0.0018 ( 60−10 ) ]1.14
2.540.070
W
h́ p−v =2.748
m2 k
Coeficiente de transferencia de calor por radiación entre el vidrio y el ambiente.
hr =σε (T v +T a)(T 2v +T 2a)
v−a
hr =(5.67∗10−8 )(0.88)(304.4+292.13)(304.4 2+ 292.132)
v−a
W
hr =5.298 2
m k
v−a
Coeficiente de transferencia de calor por convección entre el vidrio y el
ambiente.
W
h́v−a =13.90 2
m k
Luego:
−1
1 1
U ¿=
( hr p−v
+
+ h́ p −v hr + h́v−a
v− p
)
−1
1 1
U ¿= ( +
5.90+2.748 5.298+ 13.90 )
W
U ¿ =5.962 2
m k
8. Verificando la temperatura del vidrio ( T v ) para el nuevo U ¿
ENERGIA SOLAR TERMICA
33
U ¿ (T p−T a )
T v =T p−
h́ p−v +hr p−v
5.962(60−19.13)
T v =60−
2.748+5.90
°
T v =31.82 C
Como el valor de T U es próximo al valor anterior, concluimos que U ¿ es
correcto.
9. Calculo de U fondo
K a 0.4
U fondo = = =1
Z a 0.4
10. Cálculo de U lados
K A ZC PC 0.04( 0.145)( 4.2) 2
U lados= = =0.405W /m k
Z A AC 0.04 (1.5)
11. Calculo de U L
U L=U ¿ +U fondo +U lados
U L=5.962+1+0.405=7.367 W /m2 k
D. CALCULO DE LA EFICIENCIA DEL COLECTOR
Asumiendo la radiación instantánea es igual a 5000W/ m2
La transmitancia del vidrio en el espectro solar es igual a 0.88
1. El valor de Tα
Tα
Tα=
1−( ( 1−α ) pd)
La reflectancia difusa para una cubierta es pd=0.16
0.88(0.9)
Tα= =0.8
1−( ( 1−0.9 ) 0.16)
ENERGIA SOLAR TERMICA
34
2. Calculo de Qutil (calor de calentamiento de agua)
Qutil=H T A C ( Tα )−U L A C (T P −T a )
Qutil=5000( 1.5) ( 0.8 ) −7.367(1.5)(333−292.13)
Q util=5548.37 W
3. Calculo de la eficiencia instantánea de la terma solar
Q util 5548.37 3096.10
η= = = =0.74
H T AC 5000 ( 1.5 ) 4200
E. Calculo de la eficiencia de placa absorbedora
1. Calculo de m
m=
√ √
UL
Kδ
=
7.367
−3
386(0.508 x 10 )
=6.278 m−1
m(W −D)
2. Calculo de
2
m(W −D) 6.261 [ 0.1−(12.7 x 10 ) ]
−3
= =0.273
2 2
Finalmente, calculamos F:
F=
tanh [
m(W −D)
2
=
]
tanh(0.273)
=0.974
m(W −D) 0.273
2
F. Calculo del factor de eficiencia del colector ( F' )
Tomamos las siguientes propiedades termofísicas del agua a una temperatura de 50ºC
Tabla 7: Propiedad Termofísica vs Magnitud
Propiedad termofísica Magnitud
Viscosidad dinámica ( ν ) 0.658x 10−6 m2 / s
Dilatación térmica ( α ) m2
Conductividad térmica ( κ ) 1.533x 10−7
s
Dilatación volumétrica ( β ) 0.6395x 10−6 W /mk
ENERGIA SOLAR TERMICA
35
0.18x 10−6 −1
k /s
1. Para convección libre , calculamos el número de Ra
gβ ΔT D i3 9.8 ( 0.18 x 10−3 ) (60−50)(11.1 x 10−3)3
Ra=G r Pr = =
να (0.658 x 10−6)(1.533 x 10−7)
donde: Δ T =60 ºC−50 ºC=10ºC
Ra=239166.02
2. Calculo del número de Nusselt
1 1
N U =9 R a =0.47 ( 239166.02 ) 4
4
N U =10.39
3. A partir del número de Nusselt, hallamos el coeficiente conectivo entre el tubo y
el fluido
hf D i
NU=
k
Despejamos:
N k 10.39 (0.6395)
hf = U =
Di 11.1 x 10
−3
W
hf =598.6
m2 k
4. Hallamos el factor de eficiencia ( F' )
1
UL
F' =
w
[ 1
+
1
U L ( D e −(w−De ) ) F 3.14 DC hf ]
ENERGIA SOLAR TERMICA
36
1
' 6.74
F=
0.14
[ 1
6.74 [ 12.7 x 10 + ( 0.1−12.7 x 10 ) (0.97) ]
−3 −3
+
1
−3
3.14 (11 x 10 )(598.6) ]
' 0.148
F=
0.14 [
1
+
1
0.674 20.676 ]
F' =0.69
VII. SELECCIÓN DE EQUIPOS, COSTOS Y PRESUPUESTOS
Tabla 8: Comparación de los beneficios económicos de los calentadores solares de
agua respecto otros tipos de calentadores
CALENTADOR TERMA TERMA DUCHA TERMAS
SOLAR ELECTRICA ELECTRICA GAS (GLP)
DESCRIPCION 100 Lt 100 Lt - 6 Lt
Inversión inicial S/. 1,500 S/. 1,200 S/. 50 S/. 650
del equipo
Tiempo de vida 15 10 2 3
(años)
Salida de agua 100 litros 100 litros - 100 litros
caliente diaria
60 °C
Numero de días S/. 0.00 S/. S/. S/.
usado 0.402/KWh 0.402/KWh 0.0502/MJ
anualmente
ENERGIA SOLAR TERMICA
37
Costo de energía 0 S/. 902 S/. 665 S/. 383
Costo de energia S/. 75 S/. 150 - S/. 150
por año
Inversión en S/. 1,500 S/. 2,400 S/. 375 S/. 3,250
equipos en 15
años
Total de gastos S/. 2,625 S/. 18,180 S/. 10,350 S/. 11,245
en 15 años S/.
Riesgo No Tiene Tiene Tiene
Polución No Tiene Tiene Tiene
PLANO DE LA INSTALACION DE LA TERMA SOLAR
ENERGIA SOLAR TERMICA
38
Figura 14: plano de la instalación de una terma solar
PRESUPUESTO GENERAL
Tabla 9. Costos por fabricación propia
Descripción Medidas Cantidad P. unit. Total
ENERGIA SOLAR TERMICA
39
Tanque Solar 104.05
Plancha galvanizada 0.75mmx2.40X1.20m 2.88m2 1 1.00Uni 60.00 60.00
Niples galvanizados de 8 cm x 1/2” con ½” 2.00Uni 2.50 5.00
rosca
Barillas de estaño 6.00Uni 5.50 33.00
Pasta de soldar 0.25Uni 5.00 1.25
Pintura (225ml) 1.00Uni 2.50 2.50
Union Te con rosca Hydro (PP) ½ 1Uni 1.00 1.00
Tubo Hydro ½” 0.40m 3.40 1.30
Armazón de madera 15.00
Clavo 3" (Kg) 0.13Kg 4.00 0.50
Clavo 2” (kg) 0.1Kg 4.00 0.50
Listón 2”x3"x3m 2.00uni 7.00 14.00
Recubrimiento de tanque 111.30
Tubo Rectangular 1"x2"x0.75mm x6.0m 0.45m 25.80 1.94
Plancha galvanizada 0.25mmx0.90x2.40 1Unid 28.00 28.00
Vidrio Semidoble 3mm (m2) 1.00 0.65 25.97 17.00
Vidrio Doble (6mm) 1.00 0.75 49.33 37.00
Silicona 1.00uni 8.00 8.00
Tecnopor (2") 2.88 1.44m2 19.40 9.70
Lana de vidrio (1x1.2m) o Tecnopor 1.20 0.19m2 10.00 1.62
Tornillos 18.00unid 0.20 3.60
Varillas livianas (m) 3.70m 1.20 4.44
Total Terma Solar sin Conexiones 230.35
Tabla 10: Costos para un tanque con terma en el teccho
Conexión Balde / Tanque con Terma en el Techo 55.50
ENERGIA SOLAR TERMICA
40
Tuberia Hydro (PP) ½” (6m) 3.0m 20.00 10.0
Unión Universal 1/2" (PP) 2.00 3.50 7.00
Codos de 1/2" con rosca (PP) 6.00 1.00 6.00
Balde 19Lit o mas grande 1.00 6.00 6.00
Bomba de Water / Válvula de ingreso 1.00 15.00 15.00
Unión T de 1/2” con rosca (PP) 1.00 1.00 1.00
Cinta Teflón 2.00 1.50 3.00
Listón 2” x 3” x 3.0m 1.0 7.00 7.00
Clavos de 3” y4” 0.1Kg 5 0.50
Conexión a Red 93.50
Llaves de paso 1/2" 1.00 8.50 8.50
Llaves de Ducha ½” 2 15.00 30.00
Tuberia (PP) ½” C10 * 1Unid 20.00 20.00
Tuberia (PVC) ½” C10 para roscas * 1Unid 13.00 13.00
Difusor de ducha (cabeza) * 1Unid 10.00 10.00
Codos de 1/2" con rosca * 3.00 1.00 3.00
Union T de 1/2” con rosca * 2.00 1.00 2.00
Union Universal 1/2" 2.00 3.50 7.00
Total Conexión a la red y instalación en el baño 149.00
SELECCIÓN DE EQUIPOS:
ENERGIA SOLAR TERMICA
41
Figura 15: Terma solar plana de 120 litros- Termoinox
Tabla 11: Descripción de una placa plana - Termoinox.
Modelo Placa plana
Marca Termoinox
Procedencia Perú
Capacidad 120 L
Peso en vacío 62 Kg
Peso en bruto 182 Kg
Medidas de 0.48 m. x 1.14 m
tanque
Medidas de panel 1.0 m. x 1.5 m. x 0.08 m
Entrada de agua Ingreso de agua fría. Salida de agua caliente. Termosifón.
Número de puntos 2
simultáneos
Número de 03 a 04 Duchas personas / día
personas
Material Termotanque de acero inoxidable
Potencia 2000 W
Regulador de Si
presión de
temperatura
Aislamiento Poliuretano Rígido Ecológico
Tipo de Red pública/ Reservorio de agua o mixta
instalación
Almacenaje Hasta 24 horas
caliente
Válvula Si
ENERGIA SOLAR TERMICA
42
antiretorno
Color Aluminio
Garantía 5 años (Termotanque y colector solar), 1 año (Sistema auxiliar
eléctrico) y 06 meses (Accesorios, tuberías y válvulas).
Características Mayor tiempo de vida útil, el colector posee una carcasa de
aluminio anodizado pesado y extra fuerte, Ideal para zonas de
humedad elevada y zonas cerca al mar. Fondo de colector solar
metálico de alta resistencia mecánica y embate de las lluvias.
Irrompible, colector con vidrio solar prismático templado de
4.00mm de espesor con textura anti-reflectiva. Resiste impactos
de 300 kg/m2; es decir de hasta 4 personas/m2 sin romperse.
Incluye Instalación gratuita.
Instalación Se requiere: Arranques: 01 arranque de agua fría de 3/4" y 01
arranque de agua caliente de 3/4" (puntos no más de 2.5 m de
radio). Espacio requerido: 4.5 m x 4.5 m. Este espacio debe
estar limpio, libre de sombra y resistente al precio en bruto.
Arranque sistema eléctrico: 01 llave termo magnética de 2x 32
Amperios en tablero general y tendido de cables # 12 (ida y
vuelta) hasta el lugar de instalación.
Observaciones La instalación es gratis en cualquier punto de venta de Sodimac
- ciudad. Fuera de esta zona, solo se aplicará una tarifa por
traslado del instalador. Si el cliente decide hacer la instalación
por su cuenta, ésta pierde toda la garantía de fábrica.
Detalles Superficie selectiva, posee aletas de aluminio con tratamiento
adicionales superficial de Oxidación Anódica de Cromo Negro. Fuerte
aislamiento, con lana de vidrio y/o poliuretano rígido ecológico
de 30 mm de espesor. Menos pérdidas de calor. 100% de
contacto entre placa selectiva y tubos de elevación de cobre
mediante proceso de soldadura por ultrasonido. Eficiencia
energética medía 60.2%. Producción energética de 84.7 kw-
h/mes-m2
Grupo Termas solares
Categoría Baños
Familia Baños y cocinas
COSTOS DE LA INSTALACION DE LA TERMA SOLAR POR FABRICACION
PROPIA
Tabla 12: costos de la instalación de la terma solar por fabricacion propia
Costo Terma Solar sin Conexiones S/.230.35
ENERGIA SOLAR TERMICA
43
Costo de la conexión a la red y S/.149.00
instalación en el baño
Costo Final S/.379.35
COSTOS DE LA INSTALACION COMPRANDO LA TERMA SOLAR
Tabla 13: costos de la instalación comprando la terma solar.
Costo Terma Solar sin Conexiones S/.2730.00
Costo de la conexión a la red y S/.149.00
instalación en el baño
Costo Final S/.2879.00
VIII. CONCLUSIONES
1) Se asemejo las dimensiones del colector solar plano con la que se
dispone en el mercado obteniendo una área efectiva A efec =1,5 .
2) El calor entregado al agua en el colector es Qutil=5548,37 W
3) La eficiencia instantánea de la terma solar es del 74%
4) La eficiencia de la placa absorbedora es del 97%
5) El factor de eficiencia del colector F’=0,69
6) Vemos que es más económico fabricar nuestro colector solar plano
implementando a la conexión de la red que adquiriendo una terma solar
de tienda se ahorra aproximadamente el 86,84%
IX. REFERENCIAS BIBLIOGRAFICAS
ENERGIA SOLAR TERMICA
44