OPS/CEPIS/03___

UNATSABAR

UNIDAD DE APOYO TCNICO

PARA EL SANEAMIENTO BSICO DEL REA RURAL

TEORA PARA EL DISEO DE

CALENTADORES SOLARES DE AGUA

Centro Panamericano de Ingeniera Sanitaria y Ciencias del Ambiente

rea de Desarrollo Sostenible y Salud Ambiental
Organizacin Panamericana de la Salud
Oficina Sanitaria Panamericana Oficina Regional de la
Organizacin Mundial de la Salud

Auspiciado por:

Agencia Suiza para el

Desarrollo y la Cooperacin

UNIDAD DE APOYO TCNICO PARA

EL SANEAMIENTO BSICO DEL REA RURAL

TEORA PARA EL DISEO DE

CALENTADORES SOLARES DE AGUA
Sixto Guevara Vsquez
UNATSABAR OPS/CEPIS

Centro Panamericano de Ingeniera Sanitaria y Ciencias del Ambiente

rea de Desarrollo Sostenible y Salud Ambiental
Organizacin Panamericana de la Salud
Oficina Sanitaria Panamericana Oficina Regional de la
Organizacin Mundial de la Salud
Auspiciado por la Agencia Suiza para el Desarrollo y la Cooperacin
Lima, 2003

OPS/CEPIS/03___
UNATSABAR

Tabla de contenido
Pgina
1. Introduccin ..................................................................................................... 3
2. Conceptos previos ........................................................................................... 3
2.1 Radiacin solar ....................................................................................... 3
2.2 Balance de energa en el calentador solar .............................................. 4
3. Balance de energa en colector solar ............................................................... 6
3.1 Coeficiente total de transferencia de calor ............................................. 7
3.1.1 Conductancia en la parte inferior, Ufo ........................................ 8
3.1.2 Conductancia en la parte superior Usup ..................................... 8
3.1.3 Conductancia en la parte lateral Ul. .......................................... 10
3.2 Temperatura a lo largo del tubo y el factor de
eficiencia del colector .......................................................................... 11
3.3 Distribucin de la temperatura del agua a lo largo del tubo ................ 15
3.4 Factor de remocin de calor ................................................................. 16
4. Prdidas de calor en el depsito de almacenamiento y las conexiones .......... 17
4.1 Depsito de almacenamiento ............................................................... 17
4.2 Conexiones ........................................................................................... 17
5. Circulacin natural del agua en el calentador solar ....................................... 17
5.1 Ubicacin del depsito de almacenamiento respecto al
colector solar ........................................................................................ 18
5.2 Altura de presin de la circulacin natural .......................................... 18
5.2.1 Altura de presin debido al cambio de densidad
del agua en el colector solar .................................................... 18
5.3 Prdidas de carga ................................................................................ 19
6. Desarrollo de la OPS/CEPIS. ......................................................................... 20
7. Referencia bibliogrfica ................................................................................ 20
Anexo 1. Simbologas empleadas ................................................................. 21

-2-

OPS/CEPIS/03___
UNATSABAR

TEORA PARA EL DISEO DE CALENTADORES SOLARES DE AGUA

1.

Introduccin

En este documento se presentan brevemente las formulaciones acerca de la teora de

calentadores solares de agua que describen los fenmenos fsicos que se producen
en el aprovechamiento de la energa solar.
La recopilacin de la informacin presentada est enfocada en un modelo
de calentador solar con colector plano tipo rejilla y de circulacin natural. Antes, se
analizan algunos conceptos previos sobre radiacin solar y balance de energa en el
calentador solar con el objetivo de ayudar a la comprensin de las formulaciones
tratadas en este documento.
El contenido principal del documento gira en torno al anlisis de tres puntos:
i) balance de energa en el colector solar plano, ii) prdidas de calor en las conexiones y en
el depsito de almacenamiento, y iii) circulacin natural del lquido en el calentador solar.
En
el anlisis se tiene en cuenta los parmetros ms importantes que influyen en el diseo del
calentador solar, tales como: radiacin solar, temperatura ambiental, material de los tubos,
distancia entre tubos, nmero de cubiertas de vidrio, material del aislante, espesor
del aislante, ubicacin del depsito de almacenamiento respecto al colector y
volumen del depsito de almacenamiento.
Este es un documento tcnico, realizado gracias al apoyo de la Agencia Suiza para
el Desarrollo y la Cooperacin (COSUDE).
2.

Conceptos previos

2.1

Radiacin solar

La energa solar, como recurso energtico terrestre, est constituida simplemente

por la porcin de luz que emite el Sol y que es interceptada por la Tierra.
La intensidad de la radiacin solar en el borde exterior de la
atmsfera, considerando la distancia promedio entre la Tierra y el Sol, se llama constante
2
solar, y su valor medio es 1353 W/m , la cual vara en un 0,2% en un perodo de
30 aos. La intensidad de energa real disponible en la superficie terrestre es menor
2
que la constante solar, siendo alrededor de 1000 W/m , debido a la absorcin y a la
dispersin de la radiacin que origina la interaccin de los fotones con la atmsfera. Esta
porcin de energa
se conoce como radiacin directa.
Otra parte de la energa solar que llega a la superficie de la tierra se
denomina radiacin difusa que es aquella energa solar reflejada por la atmsfera terrestre,
en especial por las nubes.

-3-

OPS/CEPIS/03___
UNATSABAR

Adems, a nivel del suelo se tiene la radiacin reflejada que es parte de la energa
reflejada por los objetos terrestres. Por ejemplo, la proveniente de una pared blanca,
un charco de agua o un lago, etc.
Radiacin total es la suma de las radiaciones directa, difusa y reflejada que se
reciben sobre una superficie.
De otra parte, en el diseo de calentadores solares, la radiacin que interesa es
la que llega a la parte superior de la superficie horizontal del colector solar. En este caso,
se puede considerar la falta de radiacin reflejada y por lo tanto, la energa que
recibe el colector solar se le denomina como radiacin global.
De este modo, la radiacin global es la suma de las radiaciones directa y difusa.
Un caso particular, pero de mucho inters en el estudio de calentadores solares, es
la radiacin total sobre una superficie horizontal expuesto hacia el sol. En este caso
puede considerarse que no existe radiacin reflejada y se conoce tambin como radiacin
global. Por tanto, la radiacin global es la suma de la directa ms la difusa, esta radiacin
global es
la que se aprovecha en colectores planos.

Depsito de
almacenamiento
Conexin de
agua caliente

Colector solar

Conexin de
agua fra

Figura 1. Partes del calentador solar

2.2

Balance de energa en el calentador solar

El balance de energa tiene como objetivo determinar el calor til transmitido

al agua para establecer la eficiencia del calentador solar. El calor til qtil, empleado
para calentar el agua del depsito de almacenamiento, es una parte de la radiacin
global incidente sobre el colector solar HTAe; la otra parte es calor disipado al ambiente:
i) por el colector solar qcs, ii) por el depsito de almacenamiento qda y iii) por las
conexiones qcx. En
la figura 1 se presentan las partes del calentador solar y en la figura 2 se muestra un
diagrama del balance de energa.

-4-

OPS/CEPIS/03___
UNATSABAR

Este balance de energa se expresa analticamente como:

H T Ae qtil qcs qda

(1)

qcx

El calor til introducido en la ecuacin 1 est representado por el flujo de agua m

en el calentador y las temperaturas de entrada y salida del depsito de almacenamiento Tfe,
Tfs.

qtil m C
p

T fs

fe

Donde:

(2)

Cp = capacidad calorfica del agua

La eficiencia del calentador se define, entonces, como el porcentaje de la

radiacin solar global incidente sobre la superficie efectiva del colector y que es
aprovechado como calor til.

calor aprovechado
radiacin solar incidente

til

C
H A
T

(3)
e

En la prctica la eficiencia se determina experimentalmente en un intervalo finito de

tiempo; para fines de diseo el perodo de evaluacin t puede ser de un da o un tiempo
mayor. As, los valores puntuales obtenidos durante el periodo de evaluacin son integrados
mediante la siguiente ecuacin:
t

t
0

dtq

til

A e dt

(4)

-5-

OPS/CEPIS/03___
UNATSABAR

HT Ae
H T A e (1-) :
H TA e()

q cal

q til :

Prdidas de energa por

reflexin y absorcin
de las cubiertas

q pe :

Prdidas de calor en el colector

solar

q da :

Prdidas de calor en el
depsito de almacenamiento

q cx :

Prdidas de calor en las

conexiones colector-depsito

C alor til para el calentam iento del agua

en el depsito de alm acenam iento

Figura 2.-Diagrama del balance de energa en el calentador solar

3.

Balance de energa en el colector solar

La radiacin global que incide sobre el colector solar HTAe, disminuye su

valor como resultado de las caractersticas propias de reflexin y absorcin de la
cubierta del colector; de ah que, este valor, es corregido por un factor (). Parte
de esta radiacin, HTAe(), es empleada como calor para calentar el agua qcal; la parte
restante constituye el almacenamiento de energa en el colector mismo dU/dt, que
usualmente es despreciable; as como, las prdidas de calor hacia el ambiente por
conduccin, conveccin y re-radiacin qpe.

Donde:

H T Ae ()
q

cal

dU
q pe
dt

(5)

dU/dt0
Por lo que, la ecuacin 5 se transforma en:
(6)

H T Ae () qcal q
pe

Relacionando las ecuaciones 6 y 1, se obtiene:

q cs ( ) q

(7)

pe

qcal
q
() til q da qcx

(8)

-6-

OPS/CEPIS/03___
UNATSABAR

Para la determinacin del calor de calentamiento del agua en el colector solar

es
necesario determinar las prdidas de calor en el colector mediante la siguiente ecuacin:
q pe U L Acs(TP Ta
)

(9)

Donde:
UL =
Acs =
Tp =
Ta =

coeficiente total de transferencia de calor

rea del colector solar
temperatura promedio de la placa de absorcin
temperatura ambiental

El coeficiente total de transferencia de calor UL se determina a partir de la

temperatura promedio de la placa de absorcin, la cual a su vez, requiere la determinacin
de la temperatura del tubo y del agua a lo largo de todo el tubo y del factor de remocin de
calor.
3.1

Coeficiente total de transferencia de calor UL

La evaluacin numrica del coeficiente total de transferencia de calor UL requiere

determinar las conductancias inferior, superior y lateral en el colector solar.
En la figura 3 se muestra el circuito trmico de un colector solar plano que ayuda
a visualizar los fenmenos fsicos que se producen en el colector solar y as, a partir de
este circuito se plantean las ecuaciones de las conductancias.

Temperatura ambiente, Ta

Reflexin
HT

R5
Absorcin

Tc1

Cubierta
exterior

R4
Tc2
S
Calor til
(qtil)

Cubierta
interior

R3
Tp

Ta
Placa de
absorcin

R1
Tf

Aislante
fondo

1
UL

S
qtil
Tp

R2
Temperatura Ambiente, Ta
(i) circuito en detalle

(ii) circuito equivalente

Figura 3. Circuito trmico de un colector solar plano con dos cubiertas

-7-

OPS/CEPIS/03___
UNATSABAR

3.1.1 Conductancia en la parte inferior, Ufo

La conductancia se evala conociendo la resistencia trmica conductiva R1 y la
resistencia trmica convectiva R2.
U fo

1
R1 R2

R1 est dado por la conductividad del aislante ka, y la longitud transversal del
mismo l; R2 est dado por el coeficiente convectivo de transferencia de calor entre la placa
de absorcin y el ambiente, el cual en colectores bien aislados es despreciable.
R1 = l/ka
R2 = 1/hfo
Si R2 R1 entonces R2 se desprecia; por lo que la conductancia en la parte
inferior sera:
ka
1
U fo
R1
l

(10)

3.1.2 Conductancia en la parte superior Usup

Esta conductancia se evala determinando las resistencias trmicas R3, R4, y R5, las
cuales se muestran en la figura 3. El calor se transfiere entre la placa de absorcin
y la cubierta interior por conveccin y radiacin en forma paralela; esto mismo, tambin
ocurre entre las dos cubiertas, interior y exterior. Los mismos mecanismos de
transferencia son igualmente importantes en la disipacin de calor por parte de la
cubierta exterior hacia el ambiente.
1
U sup
R3 R4 R5

(11)

R3: Resistencia trmica entre la placa de absorcin y la cubierta interior:

R3

Donde:

(h C 2
)

TP TC 2
hR 2 )(TP TC1

(T
T P
hR 2

)(T T
C2

(12)
2

C2

1
1
1

-8-

OPS/CEPIS/03___
UNATSABAR

=
Tp =
TC1 =
TC2 =
h C2 =
hR2 =

C2 =
P =

constante de Boltzman
temperatura promedio de la placa de absorcin
temperatura promedio de la cubierta exterior
temperatura promedio de la cubierta interior
coeficiente de transferencia de calor entre la placa de absorcin y la
cubierta interior
coeficiente de transferencia de calor por radiacin entre la placa de
absorcin y la cubierta interior
emitancia de la cubierta interior
emitancia de la placa de absorcin en el infrarrojo

R4: Resistencia trmica entre la cubierta exterior y la cubierta interior:

R4

TC TC
2

(hC1 hR1 )(TC 2 TC1

)

(13)

Donde:
(
T
hR1

T )
T )
C 1 (T C 2
C1
C2
1
1
1

C1

C1 =

emitancia de la cubierta exterior

h C1

coeficiente de transferencia de calor entre las cubiertas exterior e

interior
coeficiente de transferencia de calor por radiacin entre las
cubiertas exterior e interior

hR1 =

R5: Resistencia trmica entre la cubierta exterior y el ambiente:

R5

TC 1 Ta
(hc h fir )(TC1 Ta )

(14)

Donde:

fir (TC1

T )
(T

h
C1

C1

T )
a

Ta = temperatura ambiental
hc= coeficiente de transferencia de calor por conveccin entre la
cubierta y el ambiente
hfir = coeficiente de transferencia de calor por radiacin entre las
cubiertas exterior e interior.
El proceso de clculo de R3, R4 y R5 es iterativo, adems requiere de datos como
las temperaturas de las cubiertas y los coeficientes convectivos de transferencia de calor,
que por lo general son difciles de evaluar; esto hace que el clculo de estas
resistencias trmicas sea complicado.

-9-

OPS/CEPIS/03___
UNATSABAR

La ecuacin propuesta por Hottel y Wertz para el clculo de la conductancia por la

parte superior, simplifica el procedimiento indicado anteriormente. Esta ecuacin se utiliza
para temperaturas de la placa de absorcin entre 25 y 225 C, con un margen de
2
error aproximado de 0,3 W/m .

T p
a
p
a

U sup
e
2
N

f 1
hv
1
C T p Ta

p 0,00591Nhv
N
T N f

(15)
Donde:
hv 5,7 3,8vv

f 1
0,04h

0,0005hv

1 0,058N
2

C 520 1 0,000051

de =70

para

0 < < 70; para 70< < 90 usar un valor

e 0,43 1 100T
p

= ngulo de inclinacin del colector

N = nmero de cubiertas
h = coeficiente de transferencia de calor por conveccin entre el aire y el
v
colector
c = emitancia de las cubiertas en el infrarrojo
Tp = temperatura promedio de la placa de absorcin
3.1.3 Conductancia en la parte lateral Ul.
Las prdidas de calor a travs de los lados del colector se evalan de manera
similar, aplicando la siguiente ecuacin:
qlados U l Ac (T p Ta ) k a '
hP

l
'

(T p Ta
)

Donde:
U l

k a ' hP
l'A c

P = permetro del colector solar

h = altura del colector solar
l = espesor del aislante por los lados
ka = conductividad del aislante lateral

(16)

- 10 -

OPS/CEPIS/03___
UNATSABAR

Finalmente, la evaluacin del coeficiente total de transferencia de calor para todo el

colector est representada por la sumatoria de las tres conductancias indicadas
anteriormente; es decir, inferior, superior y lateral.
UL = Ufo + Usup +Ul
3.2

(17)

Temperatura a lo largo del tubo y factor de eficiencia del colector

Un balance de energa en la placa de absorcin (ver figura 4) muestra que, para un

plano perpendicular al flujo, se cumple la siguiente ecuacin:
(q" (1)) x (q" (1)) x x Sx(1) U L (T Ta )x(1)
0

Donde:
S
HT

x
T

= HT ()
= radiacin solar global
= espesor de la placa de absorcin
= elemento del ancho de la placa de absorcin
= temperatura del elemento (x)

W-D

W-D
2

x
(i) Configuracin y dimensiones

(ii) Volumen de control

U L x(1)( T-T a)

S x(1)

(q"(1))x
x

(q"(1))x

x
(iii) Elemento diferencial

Figura 4.- Balance de energa en la placa absorcin

- 11 -

(18)

OPS/CEPIS/03___
UNATSABAR

Dividiendo la ecuacin 18 entre x y haciendo x 0, se obtiene que:

dq"
S U L (T Ta ) 0
dx

(19)

Donde:
q" -k

dT
dx

k = coeficiente de conductividad trmica de la placa de absorcin

q = calor por unidad de rea que fluye por conduccin de la placa
de absorcin hacia el tubo
Reemplazando el valor de q en la ecuacin 18, se obtiene la siguiente ecuacin:
2

d T UL
S
(T Ta )

2
k
UL
dx

(20)

Esta ecuacin diferencial requiere de dos condiciones de frontera, ver figura 6.

a)
dT
0
dx

b)

(21)
x 0

Si la temperatura de la placa en la vecindad del tubo es Tb:

T
T

(22)

W-D
b x
2

Resolviendo la ecuacin 20 con las condiciones de frontera establecidas se obtiene

que:
T Ta S/U
L

Tb
Ta

S/U

cosh(cx)
coshc(W D)/ 2

(23)

Donde se define el parmetro c de la siguiente manera:

2

U
k

La ecuacin 23 permite conocer la temperatura T en la placa de absorcin a lo largo

de x para cualquier posicin y a lo largo del flujo. De esta ecuacin se concluye
que, el mximo de temperatura en la placa ocurre entre los dos tubos, para x = 0,
tal como se observa en la figura 5.

- 12 -

OPS/CEPIS/03___
UNATSABAR

De la ecuacin 19 se deduce el calor por unidad de longitud que fluye por

conduccin a lo largo de la placa de absorcin hacia uno de los tubos qcond, esto es:
dT
q' cond q" WD q"WD 2k
dx W D
2
2

(24)

Derivando la ecuacin (23) se obtiene el dT/dx, con lo cual se obtiene:

q' cond . (W D)S UL (Tb T )
2
a

tanh c(W D) /

(25)

c(W D) / 2

y=C

T(x,y)
Tubo

y
x
(ii) Distribucin de temperatura en toda la placa

Tubo
x

(i) Perfil de temperatura en un plano y=C

Figura 5.- Distribucin de la temperatura en la placa de absorcin

Se define F como la eficiencia geomtrica de una superficie plana. As:

tanh c(W D) / 2
F
c(W D) / 2

(26)

La eficiencia geomtrica F vara de acuerdo al parmetro c(W-D)/2, segn se

muestra en la figura 6.
Haciendo uso de la definicin para la eficiencia de una superficie plana, la ecuacin
25 se convierte en:
q' cond (W D) F S
U

(T Ta )

(27)

Debido a que el tubo recibe radiacin solar por la parte superior.

q'rad DS U (T T )
L
a

(28)

Sumando las ecuaciones (27) y (28) se obtiene:

q' cal q' q' (W D)F D S
cond
rad
U

(Tb T )

(29)

- 13 -

OPS/CEPIS/03___
UNATSABAR

Este calor de calentamiento por unidad de longitud qcal, es transferido al agua. As,
este calor de calentamiento se presenta en trminos de la resistencia trmica de conveccin
que presenta la interfase tubo-agua, mediante la siguiente ecuacin.
(30)

q' cal Di h f (Tb T f

)

Donde:
Tf = temperatura del fluido en cualquier posicin
Di = dimetro interior del tubo
hf = coeficiente de transferencia de calor por conveccin
en la interfase tubo-fluido
W-D

1,0

0,9

0,8

W-D

0,7

0,6
0,0

0,5

1,0

1,5
c(W-D)/2

Figura 6.- Variacin de la eficiencia geomtrica en funcin del parmetro c(W-D)/2

Las ecuaciones 29 y 30 se combinan para determinar el calor de calentamiento del

agua en trminos de la temperatura ambiente y de la temperatura del agua. Despejando Tb
en la ecuacin 30 y sustituyendo en la ecuacin 29 se obtiene que:
q' cal WF'[S U L (T f Ta
)]

(31)

F se conoce como el factor de eficiencia del colector y esta dado por la relacin:
F '

1
UL

1
1

U L D (W D)F Di h f

(32)

Fsicamente, el numerador corresponde a la resistencia trmica entre la

superficie
del colector y el aire ambiente, mientras que el denominador constituye la resistencia
trmica entre el fluido y el ambiente. El factor de eficiencia del colector F es una funcin
de F, hf y UL dependen ligeramente de la temperatura, por ello, F se considera un
parmetro de diseo.

- 14 -

OPS/CEPIS/03___
UNATSABAR

3.3

Distribucin de la temperatura del agua a lo largo del tubo

Un balance de energa en un elemento diferencial de agua dentro del tubo, tal como
se muestra en la figura 7, indica que:
m

C pT
f
n

C pT y y q'cal y 0
y
f
qcaly
n

C pT f
y
n

Flujo de agua

(33)
m
C pT f
n

y y

Figura 7.- Balance de energa en un elemento diferencial de agua dentro del tubo
m /n

obtiene:

es el flujo de masa de agua por cada tubo, haciendo y 0, de la ecuacin 33 se

mC
p

dT f
nWF ' S U L (T Ta )
dy
f
0

(34)

La ecuacin se tiene la siguiente condicin de frontera:

Tf

y 0

T fe

Donde:
Tfe = temperatura del fluido a la entrada del colector
Cp = capacidad calorfica del agua
m = flujo de masa de agua en el colector
Resolviendo la ecuacin 34 empleando la condicin de frontera mencionada:
T f Ta S / U
L

T fe T a S / U L

exp( U L nWF ' y / mC P

)

(35)

La ecuacin 35 permite calcular la temperatura del fluido en cualquier posicin y.

Si los tubos del colector tienen una longitud L, la ecuacin 35 se convierte en:
T fs Ta S / U L
T fe Ta S / U L

)
exp(U
L nWF ' L / mC
P

(36)

- 15 -

OPS/CEPIS/03___
UNATSABAR

Una condicin lmite de inters es cuando L o bien, m 0. En estas

condiciones:
T fs,max Ta S /
L

(37)

Mediante esta ecuacin se determina la temperatura mxima que puede alcanzar un

colector plano. El valor de la temperatura mxima es importante en la seleccin de
los materiales que constituyen un colector. Si bien es cierto que la eficiencia trmica
del colector es igual a cero a esta temperatura, su importancia fundamental radica en el
hecho que, en condiciones de estancamiento (o flujo igual a cero), la temperatura
que puede alcanzarse en el interior del colector sobrepasara el lmite de operacin de
algunas de sus partes: aislante, empaques, etc.
3.4

Factor de remocin de calor

Conociendo la distribucin de la temperatura en el fluido, se define el factor

de remocin de calor FR como el parmetro que relaciona el calor de calentamiento
que se obtiene efectivamente en el colector, con el que se obtendra si toda la
superficie del colector se encontrara a la temperatura del fluido a la entrada Tfe.
FP
C
R

(T fs T fe )

Ae S U (T Ta )

(38)

fe

El rea efectiva del colector Ae,, se define como el rea de la placa de absorcin
expuesta a la radiacin solar, por lo tanto el rea efectiva tiene la siguiente ecuacin:
Ae WnL

(39)

S
(
mC
U/ L
T

F R P 1
Ae L
S/ L
U (T
U

fs
fe

T )

Ta )

(40)

Sustituyendo la ecuacin (36) en (40) se obtiene:

FR

mC P

exp(U
AeU L

F A e / mC P )

'

(41)

En consecuencia, el calor til en trminos del factor de remocin de calor es como

sigue:
qtil Ae FR [S U L (T fe Ta
)]

(42)

- 16 -

OPS/CEPIS/03___
UNATSABAR

4.

Prdidas de calor en el depsito de almacenamiento y las conexiones

4.1

Depsito de almacenamiento
Las prdidas de calor en depsito de almacenamiento se expresan como:
q
U
da

Uda =
Ada =
Tfda =
4.2

A
da

da

T T
fda

(43)

coeficiente global de prdidas de calor en el depsito de

almacenamiento
rea externa del depsito de almacenamiento
temperatura promedio del agua dentro del depsito de
almacenamiento

Conexiones

Por otro lado, el balance de energa en las conexiones se divide en prdidas en la

conexin de agua fra y en la conexin de agua caliente, as:
q
U
cx

Ucx =
Acx =
Tfcx =
5.

A
cx

cx

T T
fcx

(44)

coeficiente global de prdidas de calor en las conexiones

rea externa de las conexiones
temperatura promedio del agua dentro de las conexiones

Circulacin natural del agua en el calentador solar

La circulacin natural o circulacin por termosifn es un fenmeno que

se fundamenta en la variacin 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. Si se colocan en un tubo en U dos volmenes
iguales de agua a las temperaturas mencionadas, el nivel de la columna de agua
ms caliente es mayor respecto a la de agua ms fra (ver figura 8), esto se explica
porque la columna de agua a 20 C tiene mayor masa que la columna de agua a 60 C.

Columna de
agua caliente

Columna de
agua Fra
P: Altura de presin
generada por la
diferencia de densidades

Figura 8. Representacin de la altura de presin generada

por agua a temperaturas diferentes

- 17 -

OPS/CEPIS/03___
UNATSABAR

Debido a que en estos tipos de calentadores solares se aprovecha la circulacin

natural de agua, no se necesita de una bomba para la circulacin del agua. Un
aspecto importante de la circulacin natural es la ubicacin del depsito de
almacenamiento respecto del colector solar.
5.1

Ubicacin del depsito de almacenamiento respecto al colector solar

En la figura 9 se muestra la disposicin de las entradas y salidas de agua tanto en el

colector como en el depsito de almacenamiento, as como las alturas o niveles en que
se encuentran ubicadas dichas entradas y salidas. A partir de esta figura, se definen
los ngulos 1 y 2 de la siguiente manera:
h h
1 arctang 3 4
M

(45)

h1

2 arctang
cos

M
L

(46)

B
h3
L

h4

A
C

h1

h2

2
O

1
M

Lcos+ M

h3+h4
C
h2
D

Figura 9. Ubicacin del depsito de almacenamiento respecto al colector solar

5.2

Altura de presin de la circulacin natural

La altura de presin de circulacin natural PCN es la altura resultante de la

diferencia de la presin generada por la diferencia de densidades del agua P y la prdida
de carga por friccin en los tubos hp.
(47)

PC P h p
N

5.2.1 Altura de presin debido al cambio de densidad del agua en el colector solar
Para calcular la presin generada por el cambio de densidad del agua, el calentador
solar se representa como se muestra en la figura 10.

- 18 -

OPS/CEPIS/03___
UNATSABAR

Conexin de
h h
agua caliente 3+ 4
Colector
solar

Depsito de agua
h3 caliente

h2

Conexin
de agua fra

Figura 10. Representacin en columnas de agua de diferentes

densidades en el calentador solar

La entrada del colector el agua tiene una densidad, e , mientras que a la salida
el agua tiene una densidad s, adems se asume una densidad promedio .
Aplicando la ecuacin de Bernoulli se obtiene:

Adems,

P h
3h
s

h h
2h 2

h
s

(48)

s e

2

(49)

Reemplazando la ecuacin 49 en la ecuacin 48, se obtiene:

5.3

h2 h3
h4
2

e s

(50)

Prdidas de carga

Para determinar las prdidas de carga se determinan el rgimen del flujo de agua, la
longitud de tubos equivalente Lt y el coeficiente de friccin del flujo.
El rgimen del flujo se determina con el nmero de Reynolds:
Df
Re

v
h

(51)

vf = velocidad del agua dentro de los tubos

Dh = dimetro hidrulico
= peso especfico del agua
Por lo general, el rgimen del flujo de agua es laminar, por lo que el coeficiente de
friccin f se determina de la siguiente manera:
f 64 /
Re

(52)

- 19 -

OPS/CEPIS/03___
UNATSABAR

La longitud de tubos equivalentes Lt se determina mediante la siguiente ecuacin,

ver figura 9:
h1
Lt L
sen

h3 h4

sen 1

(53)

Finalmente, las prdidas de carga se determinan por:

t

hp

L f v

(54)

2Dh
g

Para determinar la velocidad del agua en un calentador solar de circulacin natural

se requiere un clculo iterativo, por lo cual se asume una velocidad del agua en el colector
solar, luego, se determina la presin de circulacin natural y las prdidas de carga,
finalmente se calcula la velocidad del agua, hasta que el valor asumido coincida
con el valor calculado.
6.

Desarrollo de la CEPIS/OPS

A partir de estas formulaciones de desarroll un programa de clculo el cual

nos permiti disear dos modelos de colectores solares denominados CS1 y CS2.
Ambos modelos son muy similares, variando solamente en el nmero y material de los
tubos que conforman la placa de absorcin. La eficiencia global de estos modelos se
presenta en el cuadro 1.
Cuadro N 1. Eficiencias globales en los sistemas puntual y continuo

Eficiencia global diaria

sistema puntual (gsb)
Eficiencia global diaria
sistema continuo (gsc)

Modelo CS1

Modelo CS2

43%

42%

53%

52%

7.

Referencia bibliogrfica

ATECYR. Aplicaciones de la energa solar a baja temperatura. Editorial Index,

1984
BAITSELL George. Uso directo de la energa solar. H. Blume ediciones, 1982
BECKMAN William A. Proyecto de sistemas trmico-solares. Editorial Index,
1984
DUFFIE John A. Solar engineering of thermal processes. John Wiley & Sons, 1976
HUNT Daniel V. Diccionario de energa. Publicaciones Marcombo, 1976
MANRIQUE J.A. Energa solar, fundamentos y aplicaciones fototrmicas.
Industria editorial fototrmicas, reg. 723, 1984
McCARTNEY Kevin. Agua saliente solar H. Blume ediciones, 1980
PITTS Donald. Teora y problemas de transferencia de calor. Coleccin Schaum,
1977
RAU Hans. Energa solar. Marcombo ediciones, 1977
VALERA, Anbal. Energa solar. Lima, 1993

- 20 -

OPS/CEPIS/03___
UNATSABAR

Anexo 1
Simbologas empleadas
Acs =
Acx =
Ae =
2
c
=
Cp =
Di =
Dh =
dU/dt =
F
=
h
=
h1 =
h2
h3

=
=

h4

hc=
hC1
hf
hfir
hfo

=
=
=
=

hR1 =
hR2 =
hv

mt
M

=
=
=
=
=
=
=
=
=
=
=

n
N
P

=
=
=

hC 2

HT
k
ka
l
l
Lt
m

rea del colector

rea externa de las conexiones
rea efectiva del colector
Constante para determinar el factor de eficiencia geomtrico
Capacidad calorfica del agua
Dimetro interior de los tubos de la rejilla colectora
Dimetro hidrulico de los tubos de la rejilla colectora
Incremento de energa interna
Eficiencia geomtrica
Altura del colector
Distancia vertical de la base del colector a la entrada de agua caliente del
depsito
Distancia vertical de la base del colector a la salida de agua caliente del colector
Distancia vertical de la salida de agua fra del depsito a la entrada de agua
caliente del depsito de almacenamiento
Distancia vertical de la salida de agua caliente del colector a la salida de agua fra
del depsito de almacenamiento
Coeficiente de transferencia de calor por conveccin entre la cubierta y le aire
ambiente
Coeficiente de transferencia de calor entre las 2 cubiertas
Coeficiente de transferencia de calor por conveccin en la interfase tubo fluido
Coeficiente de transferencia por radiacin entre la primera y la segunda cubierta
Coeficiente de transferencia de calor por conveccin entre el fondo del colector y
el ambiente
Coeficiente de transferencia por radiacin entre la cubierta superior y la cubierta
interior
Coeficiente de transferencia por radiacin entre la placa de absorcin y la
cubierta interior
Coeficiente de transferencia de calor por conveccin entre aire y el colector
Coeficiente de transferencia de calor entre la placa y la segunda cubierta
Radiacin solar total incidente por unidad de rea
Conductividad de los tubos
Conductividad trmica del aislante
Espesor del aislante en el fondo
Espesor del aislante por los lados
Longitud de los tubos del calentador solar
Flujo de masa de agua en el colector
Flujo de masa por cada tubo
Distancia horizontal entre la salida de agua fra del depsito de almacenamiento
y la salida de agua caliente del colector
Nmero de tubos de la placa de absorcin
Nmero de cubiertas
Permetro del marco del colector

- 21 -

OPS/CEPIS/03___
UNATSABAR

qcal
qcs
qcx
qda
qprd
qtil
R1
R2
Re
S
Ta
Tb
TC1
TC2
Tf
Tfcx
Tfe
Tfs
TP
Ucx
Uda
Ufo
UL
Ul
Usup
vf
vv
1

=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=

=
P =
PCN =
C1 =
C2 =
P =
=
f
c =

=
e
=
=
s
=

Calor entregado al agua en el colector solar

Calor disipado al ambiente por el colector
Calor disipado al ambiente por las conexiones
Calor disipado al ambiente por el depsito de almacenamiento
Prdidas de calor en el colector
Calor til entregado al agua en el calentador solar
Resistencia a la conduccin a travs del aislante
Resistencia a la conveccin hacia el ambiente
Nmero de Reynolds
Radiacin solar incidente sobre la placa de absorcin
Temperatura ambiente
Temperatura de la placa en la vecindad del tubo
Temperatura de la cubierta exterior
Temperatura de la cubierta interior
Temperatura de fluido en cualquier posicin
Temperatura promedio del agua dentro de las conexiones
Temperatura del fluido a la entrada del colector
Temperatura del fluido a la salida del colector
Temperatura promedio de la placa de absorcin
Coeficiente global de prdidas de calor en las conexiones
Coeficiente global de prdidas de calor en el depsito de almacenamiento
Coeficiente de prdidas de calor inferior
Coeficiente de prdidas de calor global
Coeficiente de prdidas de calor lateral
Coeficiente de prdidas de calor superior
Velocidad del fluido en el calentador solar
Velocidad del viento
ngulo de elevacin relativa del depsito de almacenamiento, respecto dela
salida de agua caliente del colector y la entrada de agua caliente del depsito
ngulo de elevacin relativa del depsito de almacenamiento, respecto de l
entrada de agua fra del colector y la entrada de agua caliente del depsito
ngulo de inclinacin del colector
Espesor de la placa de absorcin
Diferencia de presin por termosifn
Diferencia de presin de circulacin natural
Emitancia de la cubierta exterior en el infrarrojo
Emitancia de la cubierta interior en el infrarrojo
Emitancia de la placa de absorcin en el infrarrojo
Peso especfico del fluido
Eficiencia del calentador solar
Constante de Boltzman
Promedio de densidades del fluido
Densidad del fluido a la entrada del colector
Densidad del fluido a la salida del colector
Producto transmitancia absortancia

- 22 -