GRUPO DE ENERGÉTICA

ETS de Ingenieros Industriales

Universidad de Málaga

Aplicación de la energía solar al

calentamiento de piscinas

Fernando Domínguez Muñoz

Málaga, Noviembre de 2006

1
§1. INTRODUCCIÓN

Una de las aplicaciones más extendidas y rentables de la energía solar térmica de

baja temperatura es el calentamiento de piscinas. Para mantener la temperatura del agua
de una piscina en condiciones de confort, en torno a 26ºC, se precisa de calefacción
durante todo el año debido a que los flujos de calor dominantes (evaporación y
renovación de agua) siempre representan pérdidas de calor para el agua (§4).

La energía solar es particularmente interesante en el caso de las piscinas debido al

reducido nivel térmico requerido y a la coincidencia temporal entre la disponibilidad de
radiación solar y el horario de uso de la piscina. La zona de operación habitual de los
colectores solares en esta aplicación se muestra en la figura 1. A modo de ejemplo,
consideremos una temperatura del agua de 26ºC y una temperatura de impulsión
requerida de 40ºC: la temperatura media (aritmética) del fluido en el colector1 será de
unos 33ºC, muy próxima a la temperatura ambiente durante las horas de operación. En
estas condiciones, las pérdidas de calor son pequeñas y los captadores funcionarán en la
región de alta eficiencia de la curva característica.

0.8

Región de
funcionamiento
para piscinas

0.21 (Tm-Tamb)/I

Figura 1: Recta de eficiencia para un colector acristalado simple típico. Cuando la

carga es una piscina, la temperatura de trabajo del colector es próxima a la del ambiente,
lo que minimiza las pérdidas de calor

Otras características distintivas de este tipo de instalaciones son: (1) no es necesario

utilizar un sistema de acumulación específico, ya que la gran inercia térmica del agua

1
Recuérdese que la temperatura del fluido evoluciona exponencialmente a lo largo de los tubos de un
captador solar. Por este motivo, la temperatura media del fluido no es la media aritmética entre entrada y
salida. En nuestro caso la media aritmética se utiliza a título orientativo.

2
del vaso puede realizar esta función y, en conexión con esto, (2) la piscina demanda
energía durante todas las horas del día, y en especial durante las horas de ocupación,
parte de las cuales coinciden con las horas de sol. Por este motivo carga y producción
solar están más acompasadas en el tiempo que, por ejemplo, en las aplicaciones de
A.C.S. doméstica. Compárense las figuras 2 y 3.

Captadores

A consumo

Calentador
Pérdidas instantáneo
tanque

Intercambiador
de calor
De red

Figura 2: Instalación doméstica para la preparación de A.C.S. El acumulador permite

disponer de agua caliente solar durante horas sin radiación.

Captadores

Pérdidas
Piscina

HXR solar
HXR apoyo

Caldera

Figura 3: Instalaciones solar y de apoyo para calentamiento de una piscina. Las

pérdidas (carga) son continuas durante todo el día y la inercia del vaso puede realizar el
papel de sistema de acumulación.

En el esquema de la figura 3, la instalación solar cede toda la energía captada al

3
agua del vaso, venciendo parcial o totalmente las pérdidas de calor de éste. Cuando la
energía solar captada sea nula o menor que la demandada, el sistema auxiliar aportará la
diferencia para mantener la temperatura del agua en su valor de consigna. Cuando, por
el contrario, el aporte solar supere a la demanda, el exceso de energía hará aumentar la
temperatura (energía interna) del agua por encima de la consigna, lo que provocará la
desconexión del sistema auxiliar hasta que el agua vuelva de nuevo a enfriarse por
debajo de la consigna. La piscina estará actuando entonces como un acumulador.

Estos comportamientos se ilustran en la figura 4 para dos días consecutivos. La

temperatura de consigna se ha fijado en 25ºC. Durante la noche del primer día, el
sistema auxiliar está combatiendo todas las pérdidas de calor del agua, a fin de mantener
la temperatura de consigna. Tras la salida del Sol, el aporte solar va aumentando y el
auxiliar disminuyendo en la misma cuantía. La temperatura del agua se mantiene
constante hasta que llega un momento en que el aporte solar supera la demanda y la
temperatura del agua del vaso empieza a aumentar por encima de la de consigna.

Temperatura del agua

Energía solar cedida
Energía auxiliar cedida
26.0
500
25.8

400
25.6
Temperatura [ºC]
Energía [MJ]

300
25.4

200
25.2

100 25.0

0 24.8
0 10 20 30 40 50 60
Horas

Figura 4: Evolución de la temperatura del agua y de los aportes solar y auxiliar

Esto provoca la desconexión del sistema auxiliar (línea verde a cero). Tras la puesta de
Sol (línea azul a cero), la temperatura del agua comienza a disminuir, pero el sistema
auxiliar no volverá a funcionar hasta que no baje de 25ºC. Al día siguiente se repite un

4
ciclo similar. Este ejemplo ilustra claramente cómo la instalación solar entra en el
balance de energía del vaso, reduciendo el aporte auxiliar, y cómo el agua del vaso
puede actuar como medio de almacenamiento en esa tarea.

Las piscinas se clasifican en cubiertas o interiores cuando se encuentran en el

interior de un edificio, y descubiertas o exteriores cuando están a la intemperie. En cada
caso los flujos de calor que afectan al agua y los requerimientos que impone la
normativa son diferentes. En §2 y §3 se describe la tipología de las instalaciones usuales
en piscinas exteriores e interiores respectivamente. En §4 se definen y modelan los
flujos de calor y de masa que afectan al agua del vaso en cada caso. Conocido el valor
de estos flujos, puede estimarse la carga térmica punta y la demanda térmica del vaso.
En §5 se revisan algunos programas de ordenador disponibles para calcular la fracción
solar, es decir, la fracción de la demanda total del vaso que puede cubrirse con un
sistema solar concreto (un sistema de cierta área de captación, tipo de colector…).

§2. INSTALACIONES SOLARES PARA PISCINAS EXTERIORES

En las piscinas exteriores la masa de agua está en contacto directo con la atmósfera
a través de la superficie libre del agua. En España está expresamente prohibido el uso de
fuentes de energía convencionales (electricidad, gas, gasoil,..) para calentar piscinas
exteriores. Únicamente pueden utilizarse fuentes de calor residual (de procesos
industriales, calor de condensación de equipos de climatización u otros,…) o fuentes
gratuitas (energía solar).

Nos centraremos en el uso de la energía solar. En la figura 5 se esquematiza una

instalación típica. El sistema solar se integra como un elemento más en la cadena, figura
6. En este caso se ha utilizado la bomba del sistema de filtrado y tratamiento químico
del agua para hacer circular el fluido a través de los captadores solares. Como el caudal
que mueve el sistema de filtrado es mucho mayor que el requerido por la instalación
solar, únicamente parte del caudal total se deriva hacia la instalación solar (línea roja).
Si la bomba del sistema de filtrado no pudiera vencer la pérdida de carga adicional que
supone la instalación solar, deberá instalarse una bomba independiente para circular el
fluido por los captadores.

Nótese que el agua caliente (impulsión) debe introducirse por la parte inferior de la

5
piscina, y el agua de retorno (la que va al sistema de calentamiento) debe tomarse de la
parte superior. Haciendo esto se evita la formación de un estrato de agua caliente en la
parte superior, que aumentaría las pérdidas de calor al aire.

Bomba
Desin-
Filtro fección

Figura 5: Piscina exterior

Captadores

IN -1
Controlador IN-2
OUT

Filtro Desf.
Bomba

Calentador
auxiliar

PROHIBIDO

Figura 6: Instalación solar simple para piscina exterior

En este esquema, el sistema de control (todo/nada) acciona una válvula de tres vías
que envía o no el agua hacia los captadores; esta decisión se toma en función del nivel
térmico de los captadores con respecto al agua del vaso. Como las válvulas de tres vías
motorizadas suelen ser costosas, puede utilizarse una válvula simple que bloquee el
camino directo hacia la piscina y obligue al agua a circular a través del sistema solar,
figura 7.

Lo habitual es que una piscina exterior únicamente se utilice durante el verano,

6
cuando la temperatura ambiente y la radiación son elevadas. En estas condiciones, el
captador descubierto (sin cubierta ni aislamiento) puede ser la mejor opción por coste y
prestaciones. Recuérdese (figura 1, §1) que la temperatura de operación del colector
durante el verano no tiene que ser muy superior a la temperatura ambiente.

Pueden existir otros condicionantes que desaconsejen utilizar colectores descubiertos.

Por ejemplo, una limitación en la superficie disponible para instalar el campo de
colectores: en este caso, para combatir la demanda se necesitan niveles de temperaturas
mayores a la salida (menos flujo a mayor temperatura), por lo que habrá que utilizar
colectores cubiertos. También puede ocurrir que se requiera un suministro solar a mayor
temperatura, por ejemplo, para atender a una red de agua caliente sanitaria con la misma
instalación solar.

CAPTADORES

IN -1
IN-2
OUT

Bomba
Filtro

Calentador
auxiliar
PROHIBIDO

Figura 7: Instalación solar simple para piscina exterior

El mejor dispositivo de ahorro energético en una piscina exterior es la denominada

“manta térmica”. Se trata de una lámina de material plástico que se coloca cubriendo la
superficie del agua en los momentos en que la piscina no se utiliza, figura 8. Su efecto
se comprenderá mejor tras la lectura del apartado §4, pero básicamente suprime la
evaporación (perdida de calor y agua a reponer) y reduce las pérdidas convectivas y
radiantes. También reduce apreciablemente la pérdida de agentes desinfectantes y el
ensuciamiento.

7
 Otro aspecto a tener en cuenta es que, si la manta es traslucida a la radiación solar y
opaca a la radiación de larga, y la piscina recibe suficiente radiación solar, el conjunto
manta-piscina se comportará como un colector solar con almacenamiento, ya que la
manta crea el efecto invernadero y reduce las pérdidas de calor en la superficie. Es
interesante evaluar si esta técnica pasiva de calentamiento es suficiente como para
alcanzar condiciones de baño aceptables durante el verano. Dependerá de las horas de
uso previsto y de la cantidad de radiación solar incidente.

Figura 8: Ejemplos de mantas térmicas

§3. INSTALACIONES SOLARES PARA PISCINAS INTERIORES

Las piscinas interiores están situadas dentro de un edificio, normalmente una gran
instalación deportiva, figura 9. En este caso el Reglamento de Instalaciones Térmicas en
la Edificación (R.I.T.E.) permite utilizar energías convencionales para acondicionar
tanto el vaso como el aire del local. La intensa evaporación y la elevada inercia térmica
del vaso dan lugar a unas condiciones termo-higrométricas singulares en el aire del
local, que precisa de unidades de tratamiento especiales. Las condiciones del aire
interior (temperatura y humedad relativa) están, al menos en parte, controladas: lo
normal es disponer de un sistema de calentamiento y deshumidificación, junto con un
free-cooling entálpico. Este último elemento permite, entre otros, retirar carga latente

8
(humedad) del local utilizando aire exterior cuando las condiciones exteriores sean
favorables para ello.

Figura 9: Piscina interior

Al permitir el uso de energías convencionales, el RITE limita la temperatura a la

que puede calentarse el vaso según el tipo de actividad que se realice, tabla 1.

Recreo 23 ºC
Competición 24 ºC
Chapoteo 24 ºC ± 1 ºC
Enseñanza 25 ºC
Entrenamiento 26 ºC

Tabla 1: Temperatura del agua

En piscinas de uso polivalente se suelen utilizar 26ºC; en piscinas para personas

mayores se pueden alcanzar hasta 28ºC.

La temperatura del aire debe ser del orden de 2ºC mayor que la del agua (con un
límite de 29 ºC) a fin de evitar la sensación de frío. La humedad relativa estará en torno
al 60%, con mínimo 55% y máximo 70%.

La gran mayoría de piscinas interiores las encontraremos en instalaciones

deportivas públicas o comunitarias. En estas instalaciones aparecen tres consumos:

9
 1. Calentamiento del vaso (~ 50%)
2. Agua caliente sanitaria para duchas (~ 20%)
3. Tratamiento de aire (~ 30%)

La figura 10 muestra un posible esquema de principio del sistema auxiliar de

calentamiento con bombas de calor. También es muy común utilizar calderas. En la
figura 10 se denomina colector a un recipiente al que se vierte y del que se toma el
fluido de trabajo. Por ejemplo, las bombas de calor vierten el fluido caliente al colector
de impulsión, del que después cada bomba asociada a cada consumo toma el caudal que
necesite.

Climatizadora

Vaso

A.C.S.

Colector de retorno Colector de impulsión (CALOR)

Bombas de calor

Figura 10: Esquema de principio del sistema convencional

Obviando la rama que va a la climatizadora, el sistema auxiliar con sus cargas se

muestra en la figura 11. Nótese que se utilizan depósitos de acumulación en el circuito
de agua caliente sanitaria. Estos depósitos permiten hacer frente a las puntas de
consumo (gran número de duchas simultáneas) con una máquina de menor potencia
instantánea que la que se necesitaría para hacer un calentamiento instantáneo. Los
depósitos también permiten que la máquina trabaje en un régimen más estable, lo que
mejora su eficiencia.

La instalación convencional debe dimensionarse como si la solar no existiera, a fin

10
de asegurar la cobertura de la demanda en cualquier momento. No obstante, la
concepción del sistema debe ser global ya que ambos sistemas, solar y auxiliar,
interaccionarán y deberán adecuarse el uno al otro. Es evidente que el aporte solar hará
que el sistema auxiliar trabaje a carga parcial, a veces muy reducida, durante muchas
horas del año. Para que esto no resulte en pérdida de eficiencia o deterioro del sistema
auxiliar, éste deberá estar proyectado para trabajar a carga parcial sin mayores
problemas. Una solución inmediata es no utilizar una bomba de calor de gran potencia,
sino varias bombas de calor más pequeñas, como se ilustró en la figura 10 dibujando
dos bombas de calor. Esto proporciona varios escalones de parcialización y, por tanto,
máquinas trabajando individualmente en mejores condiciones.

PISCINA
Filtro y
tratamiento

Agua para DUCHAS

Agua de
red

Acumuladores ACS-apoyo
Colector de retorno Colector de impulsión (CALOR)

Bombas de calor

Figura 11: Esquema de principio del sistema convencional completo

La siguiente cuestión es cómo incorporar el sistema solar al sistema convencional

de la figura 11. En este tipo de instalaciones es frecuente diseñar el sistema solar para

11
que atienda tanto al vaso como a la preparación de agua caliente sanitaria. Una
posibilidad se esquematiza en la figura 12.

12
 Rebosadero

Al sistema auxiliar
(Fig. 11- colector de impulsión,

DUCHAS
salida a piscina) Tratami
Captadores ento
Vd

Vm
HXR
vaso Bomba de
recirculación

HXR Al sistema auxiliar

A.C.S. (Fig. 11- colector de impulsión,
salida a duchas)
Agua de RED

Acumulación de apoyo
Acumulación solar en paralelo
en serie

Figura 12: Esquema de principio del sistema solar acoplado al sistema de la figura 11

13
 La instalación solar se conecta a los consumos a través de dos intercambiadores de
calor, “HXR vaso” y “HXR ACS” en la figura 12. La válvula motorizada “Vd” envía el
fluido calentado por los colectores a uno de esos intercambiadores de calor, al que
decida el sistema de control en cada momento.

• En el circuito del vaso, el intercambiador solar precalienta el agua de la piscina

antes de que entre en el intercambiador de calor del sistema auxiliar

• En el circuito de agua caliente sanitaria es necesario instalar depósitos de

acumulación solar, diferentes de los depósitos de acumulación del sistema
auxiliar. Como se comentó al estudiar las instalaciones de agua caliente
sanitaria, en ningún caso un sistema auxiliar deberá estar actuando sobre el
volumen de acumulación asignado a la parte solar.

La manta térmica también es útil en piscinas interiores (figura 13), aunque no

siempre se utiliza debido a que es más complicada de manejar en piscinas de gran
tamaño. Para colocarla y retirarla suelen utilizarse sistemas automáticos. En este caso y
por motivos de seguridad, no debe permitirse que la manta cubra parcialmente la piscina
mientras está en uso, es decir, la piscina debe estar o completamente cubierta o
completamente descubierta.

Figura 13: Manta térmica en una piscina interior

En una piscina interior, la manta reduce principalmente las pérdidas por

evaporación, es decir, la carga latente en la climatizadora. No parece ser una buena idea
el desconectar la climatizadora del local durante las horas en que la piscina no se utiliza,
pues el aire se saturará de humedad y aparecerán condensaciones en muchos lugares,
que con el tiempo dañarán el mobiliario y la propia estructura del edificio. Si la
climatizadora se mantiene en funcionamiento durante todo el día, la manta térmica
reducirá apreciablemente la carga en la sección de deshumidificación.

§4. FLUJOS DE CALOR Y MASA EN PISCINAS INTERIORES Y

EXTERIORES

La demanda térmica del vaso vendrá determinada por el balance de flujos de calor
que afectan al agua. Para el propósito de nuestros modelos, la masa de agua puede
considerarse isoterma (modelo de un nodo) y con propiedades constantes. También se
supondrá que el volumen de agua contenido en el vaso no varía.

En la figura 14 se esquematizan los flujos de calor y masa intercambiados entre el

agua y su entorno.
Int. larga
Evaporación. Radiación
Convección solar
Ganancias
internas

Renovación
(suele ser carga en el
equipo)
Conducción

Figura 14: Flujos de calor y masa en una piscina

Las pérdidas de calor más importantes están asociadas con la evaporación hacia el
ambiente y con la renovación, por motivos de higiene, de parte del agua del vaso con
agua fresca. La siguiente pérdida en importancia viene del intercambio de radiación de
onda larga entre la superficie del agua y las superficies con las que tiene contacto visual,
que serán las del local en el caso de una piscina interior, y el cielo en el caso de una
exterior. Las pérdidas por conducción son normalmente muy reducidas, menores del 5%
del total del resto de pérdidas.

15
 §4.1 Pérdidas por evaporación (Qev)

Este flujo, siempre pérdida y cuantitativamente el más importante, es el resultado

de tres fenómenos:

• DIFUSIÓN: justo por encima de la superficie libre del agua existe

permanentemente una fina capa de aire saturado (interfase) que, por difusión
molecular, intercambia vapor de agua con el aire. Éste es el único mecanismo de
transporte que se manifiesta en una piscina desocupada con el aire en calma

• ADVECCIÓN: en piscinas exteriores el viento renueva continuamente la masa de

aire por encima de la interfase, lo que aumenta considerablemente el nivel de
evaporación. En este caso, la evaporación es combinación de difusión y advección

• AGITACIÓN: en una piscina ocupada, el movimiento de los bañistas produce

ondas en la superficie del agua (aumento de la interfase), salpicaduras,
evaporación desde sus cuerpos,… Todos estos factores incrementan la
evaporación.

La masa total evaporada (me) puede escribirse como

m& e = ER ⋅ m& o , (1)

siendo mo la masa de agua evaporada desde la piscina desocupada y ER un coeficiente

de corrección función del número de bañistas. En la práctica, la única forma de predecir
estos términos es utilizando correlaciones experimentales.

La evaporación desde una superficie de agua en calma (sin bañistas) es

directamente proporcional a la diferencia de presiones parciales de vapor entre la
superficie del agua y el aire del local

( )
m& o = A ⋅ ( k1 + k2 v ) psat (Tag ) − ϕ ai psat (Tai ) h fg , (2)

donde v es la velocidad del aire, psat(Tag) la presión parcial de vapor del aire saturado a
la temperatura del agua, φai la humedad relativa del aire, psat(Tai) la presión parcial de
vapor del aire saturado a la temperatura del local (o exterior), y hfg el calor latente del
agua. En ocasiones la ecuación (2) se expresa, de manera equivalente, en términos de
diferencia de contenidos de humedad.

16
 La idea del proceso representado por la ecuación (2) se ilustra en la figura 15: el
agua cambia a fase vapor para incorporarse a la interfase, y después cambia su
temperatura hasta la del aire. En el modelo, este cambio de temperatura no está incluido
en (2), sino en el término convectivo.

Me
Aire (26ºC, 60%)
Cambio
de temperatura
Interfase (24ºC, 100%)

Cambio
de fase
Agua, T=24ºC

Figura 15: Transferencia de masa del agua al aire

La velocidad del aire sólo es importante en piscinas exteriores. En la práctica, su

valor está sujeto a gran incertidumbre debido a que depende de la topografía del entorno
de la piscina (exterior). Normalmente se conoce de la velocidad de la corriente de aire
no perturbada, es decir, medida en campo abierto (típicamente en un aeropuerto) a cierta
altura del suelo (~10m). Ponderando este valor puede estimarse la velocidad del viento
en las cercanías de la superficie del agua:

v = ξ ⋅ vaeropuerto . (3)

El coeficiente de resguardo ξ tiene en cuenta la reducción de la velocidad de viento con

la altura y con la presencia de obstáculos. Valores típicos son:

ξ = 0.3 lugar no resguardado

ξ = 0.2 moderadamente resguardado
ξ = 0.1 bien resguardado
ξ = 0.15 para hacer estimaciones

Debe mencionarse que algunos investigadores cuestionan la forma funcional (2) a

favor del ajuste no lineal

(
m& o = A ⋅ ( k1 + k2 v m ) psat (Tag ) − ϕ ai psat (Tai ) )
n
h fg , (4)

17
con 0<m<1 y 0<n<1. El exponente para la velocidad da cuenta de la transición
turbulenta a cierta distancia del borde de ataque de la corriente de aire. El exponente en
la fuerza impulsora manifiesta una desviación en la analogía transferencia de calor-masa
vinculada al movimiento que la masa evaporada induce en dirección normal a la
superficie del agua, fenómeno que ensancha la capa límite de densidad de vapor y que,
consecuentemente, reduce el gradiente en su interior. No existen ajustes de (3)
específicos para piscinas.

Veremos ahora cuatro de las muchas correlaciones que han sido propuestas a lo
largo del tiempo. Para piscinas interiores y exteriores, Smith [1] propone utilizar los
siguientes coeficientes en la expresión (2)

( k1 , k2 ) = ( 0.2314 , 0.231056 ) , (5)

valores que han sido obtenidos a partir de ensayos y que corrigen los que históricamente
se han recomendado en el conocido libro ASHRAE Applications [2]. También según
estos autores, el término de corrección por ocupación se calcula como

ER = 1.04 + 4.2735 ⋅ N A . (6)

Así, la correlación de Smith para piscinas es

( )
m& e = A ⋅ (1.04 + 4.2735 ⋅ N A ) ⋅ ( 0.2314 + 0.23105v ) psat (Tag ) − ϕ ai psat (Tai ) h fg , (7)

La velocidad del viento únicamente es importante en piscinas exteriores, de manera

que para piscinas interiores (7) se reduce a

( )
m& e = 0.2314 ⋅ A ⋅ (1.04 + 4.2735 ⋅ N A ) ⋅ psat (Tag ) − ϕai psat (Tai ) h fg . (8)

El significado y las unidades de las variables son los siguientes:

me Masa de agua evaporada [kg / h]

A Área de la superficie libre del agua (m2)

N Número de bañistas

v Velocidad del viento por encima de la superficie del agua [m/s]

Psat(Tag) Presión parcial de vapor en saturación a la temperatura del agua [Pa]

Psat(Tai) Presión parcial de vapor a la temperatura del aire [Pa]

φai Humedad relativa del aire %/100

18
 hfg Calor latente de cambio de fase del agua ≈ 2448 kJ/kg

La siguiente correlación sólo es aplicable a piscinas interiores. Shah [3] revisó toda
la literatura disponible sobre este tema y publicó la siguiente correlación a partir de un
ajuste de los datos publicados por otros autores:

⋅ (Wag − Wai ) ,
1/ 3
m& o = k ⋅ A ⋅ ρ ag ⋅ ρ ai − ρ ag (9)

donde
⎧⎪35 ρ ai − ρ ag > 0.02
k=⎨ . (10)
⎪⎩40 ρ ai − ρ ag ≤ 0.02

Para considerar el efecto de los ocupantes propone la siguiente corrección ER,

⎧ 1 + 14.85 ⋅ N A > 45
⎪ A N
⎪
ER = ⎨1.2 + 5.85 ⋅ N 4.5 ≤ A ≤ 45 , (11)
A N
⎪
⎪ 2.5 A < 45
⎩ N

donde

mo Masa de agua evaporada en la piscina desocupada [kg / h]

ρag Densidad del agua [kg / m3]

ρai Densidad del aire [kg / m3]

Humedad absoluta del aire en saturación a la temperatura del agua
Wag
[kg agua / kg aire seco]
Humedad absoluta del aire en las condiciones del local [kg agua / kg
Wai
aire seco]
pag Presión parcial de vapor en saturación a la temperatura del agua [Pa]

pai Presión parcial de vapor saturado a la temperatura del aire [Pa]

φai Humedad relativa del aire %/100

N Número de bañistas

A Área de la superficie libre del agua [m2]

La correlación de Pizzetti se encuentra con cierta frecuencia en proyectos en

19
España, y en principio sólo debería utilizarse para piscinas interiores debido a que no
incluye la dependencia con la velocidad del viento

m& e = A ⋅ (16 + 133 ⋅ N ) ⋅ (Wag − Wai ) , (12)

donde las variables tienen el mismo significado y unidades que las comentadas
anteriormente.

Otra correlación es la de Biasin, utilizada en el componente no estándar de

TRNSYS “type344” [4] para piscinas exteriores

( )
m& o = A ⋅ ( 0.1526 + 0.20347 ⋅ v 0.5 ) ⋅ psat (Tag ) − ϕai psat (Tai ) h fg , (13)

donde v [m/s] es la velocidad del viento medida a 3 metros por encima de la superficie
de la piscina. Este autor no facilita una expresión analítica para considerar el efecto de
los ocupantes. De acuerdo con el manual del componente type344, deben utilizarse los
siguientes valores para tal fin

ER = 1 para piscinas desocupadas

ER = 1.5 para piscinas privadas con actividad moderada
ER = 2 para piscinas públicas y actividad normal
ER = 2.8 para piscinas de recreo (mucha agitación del agua)
ER = 3.5 para piscinas con olas artificiales

Ejemplo: Calcule la masa de agua evaporada a lo largo de un día en una piscina interior
de 425 m2. Las condiciones de diseño son: temperatura del agua 25ºC, temperatura del
aire 27ºC y humedad relativa 60%. La piscina se utiliza durante 11 horas al día, con una
ocupación media de 20 bañistas por hora.

En primer lugar se harán todos los cálculos psicrométricos necesarios. Pueden utilizarse
tablas o expresiones analíticas [4], el diagrama psicrométrico o algún programa de
ordenador:

Psat(Tag) = presión_parcial (25ºC, 100%HR) = 3166.9 [Pa]

Psat(Tai) = presión_parcial (27ºC, 100%HR) = 3564.3 [Pa]
ρag = densidad (25ºC, 100% HR) = 1.147 [kg / m3]
ρai = densidad (27ºC, 60% HR) = 1.151 [kg / m3]
Wag = humedad_absoluta (25ºC, 100% HR) = 0.02007 [kg agua / kg aire seco]
Wai = humedad_absoluta (27ºC, 60% HR) = 0.01341 [kg agua / kg aire seco]
hfg = 2448 kJ / kg

Para piscinas interiores podemos utilizar las correlaciones (8), (9) con (11) y (12).

SMITH: La masa evaporada en la piscina desocupada es

0.2314 ⋅ ( 3166.9 − 0.6 ⋅ 3564.3)

m& o = 425 ⋅ = 41.31 [kg / h]
2448

20
 Durante las horas de uso, la masa evaporada aumentará en la cuantía dada por (6)
debido a la agitación provocada por los bañistas

ER = 1.04 + 4.2735 ⋅ 20 425 = 1.241

m& e = ER ⋅ m& o = 1.241⋅ 41.31 = 50.15 [kg / h]

La piscina se utiliza durante 11 horas al día, de manera que la masa total evaporada en
un día es
M = 41.31 ⋅ ( 24 − 11) + 50.15 ⋅ 11 = 1088.68 [kg / dia]

SHAH: La diferencia de densidades es 0.004 < 0.02, de manera que según la ecuación
(10) k=40. La masa evaporada en la piscina en calma se calcula mediante la ecuación
(9)

⋅ ( 0.02007 − 0.01341) = 21.28 [kg / h]

1/ 3
mo = 40 ⋅ 425 ⋅ 1.147 ⋅ 1.1514 − 1.147

Para calcular el efecto de la ocupación se utiliza la ecuación (11). Ya que A/N = 425/20
= 21.25 < 45, se tiene
20
ER = 1.2 + 5.85 ⋅ = 1.475 ,
425
de manera que

m& e = 1.475 ⋅ 21.28 = 31.388 [kg / h]

La masa de agua evaporada en un día será

M = 21.28 ⋅ ( 24 − 11) + 31.388 ⋅ 11 = 621.908 [kg / dia]

PIZZETTI: La masa evaporada en la piscina desocupada se calcula haciendo N=0 en

(12),
m& o = 425 ⋅ (16 + 133 ⋅ 0 ) ⋅ ( 0.02007 − 0.01341) = 45.288 [kg / h]

Durante las horas de ocupación

⎛ 20 ⎞
m& e = 425 ⋅ ⎜16 + 133 ⋅ ⎟ ⋅ ( 0.02007 − 0.01341) = 63 [kg / h]
⎝ 425 ⎠
Para un día completo

M = 45.288 ⋅ ( 24 − 11) + 63 ⋅ 11 = 1281.744 [kg / dia]

▄

El ejemplo anterior pone de manifiesto la gran dispersión en los resultados de las

diferentes correlaciones; por ejemplo la de Pizzetti predice el doble de evaporación que
la de Shah bajo las condiciones estudiadas. Si probamos otras correlaciones disponibles
en la bibliografía comprobaremos que este problema persiste.

Profundizando un poco más en la comparativa, en la figura 16 se muestran las

predicciones de cada correlación para mo al variar la temperatura del agua manteniendo

21
las condiciones en el aire constantes en 27ºC con 60% HR.

100

Smith
80 Shah
Flujo másico evaporado [kg/h]

Pizzetti

60

40

20

0
20 22 24 26 28 30
Temperatura del agua [ºC]

Figura 16: Evolución de mo con Tagua

La figura 17 muestra la evolución de me con el número de usuarios cuando se

mantienen constantes las condiciones en el agua (25ºC) y en el aire (27ºC, 60%HR)

120

Smith
Shah
100 Pizzetti
Flujo másico evaporado [kg/h]

80

60

40

20
10 20 30 40 50 60 70
Número de bañistas

Figura 17: Efecto del número de usuarios

22
 Las discrepancias son considerables. Probablemente la correlación de Pizzetti
sobrestima la masa evaporada, especialmente cuando la ocupación es elevada. En
cuanto a la de Shah, aunque su autor afirma que es la que mejor ajusta los resultados de
distintas fuentes, vemos que predice valores bastante por debajo de la correlación de
Smith, que se obtuvo a partir de experimentos cuidadosos y bien documentados
realizados en cinco piscinas. Desde el punto de vista del diseñador y teniendo en cuenta
que la incertidumbre persiste, los resultados de Shah serían comprometedores si
resultara que subestimaran la carga. Los resultados de Smith pueden considerarse una
solución intermedia, hasta que se publiquen resultados más concluyentes.

§4.2 Pérdidas por radiación de onda larga

El intercambio de radiación de onda larga entre la superficie de la piscina y las

paredes interiores del local en el caso de una piscina interior, o el cielo en el caso de una
exterior, se expresa mediante la ley de Stefan-Boltzmann,

QOL = A ⋅ ε ⋅ σ ⋅ (Tag4 − Tx4 ) , (14)

donde ε es la emisividad del agua (~0.9), σ la constante de Stefan-Boltzmann

W
σ = 5.67 ⋅10−8 ,
m2 K 4
Tag [K] la temperatura del agua y Tx [K] la temperatura de las superficies del local o del
cielo, según corresponda. El modelo para calcular la temperatura de cielo dependerá de
la implementación particular de cada programa El lector interesado puede consultar
algunos modelos en las referencias [6], [7] y [8].

§4.3 Pérdidas por convección

El flujo de calor por convección se calcula mediante la Ley de Newton,

Qcv = Av ⋅ hcv (Tag − Tai ) . (15)

En piscinas interiores la convección vendrá determinada casi exclusivamente por la

diferencia de temperaturas entre el agua y el aire, mientras que en piscinas exteriores
también incluirá la velocidad del viento. Una correlación simple y muy utilizada es la
siguiente [8]:
hcv = ( 3.1 + 4.1 ⋅ v ) (Tag − Tai ) . (16)

23
 §4.4 Pérdidas por renovación

Por motivos de higiene, diariamente debe reemplazarse cierta cantidad de agua del
vaso por agua fresca. Típicamente se renueva un 5% del volumen del vaso cada día.
Este término suele representar una pérdida de calor porque el agua de la red está más
fría que la del vaso,
ξ ⋅V
Qrenv = ⋅ c p ⋅ ρ ⋅ (Tagua − Tred ) , (17)
24

donde ζ es el tanto por uno de volumen del vaso a renovar diariamente (por ejemplo,
0.05), V [m3] es el volumen del vaso, cp [kJ/kg K] el calor específico del agua, ρ[kg/m3]
la densidad del agua y Tred[ºC] la temperatura del agua de red.

§4.5 Ganancias por radiación

La masa de agua puede ganar energía procedente de fuentes de radiación de alta

temperatura, como es el caso del Sol o, eventualmente, ciertos tipos de lámpara.

En piscinas interiores existirá aporte solar si hay superficies vidriadas al exterior en

el local. En nuestras latitudes estas superficies suelen orientarse al sur, a fin de permitir
la entrada de radiación durante el invierno, cuando la altura solar es baja, y bloquearla
durante el verano (figura 9). Los modelos de piscina interior no suelen considerar la
ganancia de radiación solar: se trata de una hipótesis conservadora que, además,
simplifica considerablemente el problema.

En el caso de piscinas exteriores expuestas al Sol, no puede despreciarse la

ganancia solar sin cometer errores considerables en los cálculos. En la figura 18 se
ilustra el proceso que experimenta un rayo de Sol en el interior del vaso. La secuencia
es: reflexión – refracción – absorción y transmisión en el agua – absorción y reflexión
en la pared del vaso – …. El proceso continúa hasta que el rayo se desvanece.

La expresión para calcular la radiación solar absorbida por el conjunto agua + vaso
(recipiente) tendrá la forma
Gsolar = A ⋅ ⎡⎣(1 − rb ) Rb + (1 − rd ) Rd ⎤⎦ , (18)

24
 Figura 18: Seguimiento de un rayo solar en el interior del vaso

donde A[m2] es el área de la interfase, rb la reflectividad total (considerando todas las

reflexiones) de la interfase para la radiación directa, Rb [kJ/h m2] la radiación directa
incidente sobre superficie horizontal durante el periodo de cálculo considerado, rd la
reflectividad total de la interfase para la radiación difusa y Rb [kJ/h m2] la radiación
difusa incidente sobre la horizontal durante el periodo de cálculo considerado. El
cálculo de estos coeficientes depende de la implementación particular de cada
programa.

§4.6 Pérdidas por conducción

Las pérdidas por conducción tienen lugar a través de las paredes del vaso hacia el
medio circundante, que puede ser el terreno, un sótano, un vacío sanitario, etc… Según
la bibliografía, las pérdidas de calor por conducción no superan el 5% de la suma del
resto de pérdidas de calor del vaso, de manera que pueden calcularse como

Qcd = 0.05 ⋅ ( Qev + Qcv + Qlw + Qrenv ) , (19)

Para determinarlas de forma analítica es necesario conocer la temperatura del

medio circundante y la composición de las paredes del vaso. En buena medida las
paredes se encontrarán en régimen permanente debido a que las oscilaciones de
temperatura a ambos lados (agua y exterior) del muro deben ser pequeñas. En la figura
19 se tiene el esquema eléctrico equivalente del problema, con el flujo radiante de corta
absorbido por la cara interna del vaso (§4.5) impuesto en el nodo correspondiente a la
superficie interior.

El flujo por conducción vale

Qcd = Ahi (Tag − Tvi ) , (20)

25
 Iabs,v
1/hi Upared 1/he

Tag Tvi Tve Ts

Qcd

Figura 19: Efecto del número de usuarios

siendo la temperatura en la cara interior del vaso

1 ⎡ I abs ,v ⎤
Tvi = ⎢ + U p Ts + hi Tag ⎥ , (21)
U p + hi ⎣ Av ⎦

y el coeficiente global de transferencia

1
Up = , (22)
ej 1
∑
capas k j
+
he

donde hi [W/m2K] es el coeficiente de convección en la cara interior de la pared del

vaso, he [W/m2K] el coeficiente convectivo-radiante en la cara exterior de la pared del
vaso, ej [m] el espesor de la capa “j” del vaso y kj [W/m K] la conductividad de dicha
capa.

§5. PROGRAMAS DE CÁLCULO

Las ideas expuestas en los apartados anteriores pueden elaborarse algo más y
combinarse con los modelos matemáticos de radiación solar y de la instalación solar
para implementar un programa informático que calcule la fracción solar de una
instalación solar conectada a una piscina. La descripción completa de estos algoritmos
está fuera del alcance de este tema y, en todo caso, su conocimiento detallado no es
necesario en la práctica.

Los siguientes programas son muy sencillos de manejar y están disponibles

gratuitamente en Internet:

1. RETSCREEN: Se trata de un programa Excel creado y distribuido por el

Departamento de Recursos Naturales de Canadá. Puede descargarse a través de Internet
en la dirección dada en [8]. Modela piscinas interiores y exteriores. Para calcular la

26
fracción solar utiliza el método de utilizabilidad de Klein, cuya descripción
pormenorizada puede encontrarse en la referencia [6]. La base de datos meteorológicos
del programa incluye localidades de España y es ampliable por el usuario. Los datos
meteorológicos requeridos son valores mensuales, fáciles de conseguir para muchas
localidades.

2. PACSOL: este programa resuelve el sistema mostrado en la figura 12, es decir, una
instalación solar que proporciona suministro simultáneo a una piscina interior (no
modela piscinas exteriores) y a un sistema de preparación de agua caliente sanitaria. El
programa está distribuido gratuitamente por la Junta de Andalucía en la dirección de
Internet [9]. Este programa utiliza el simulador TRNSYS como motor de cálculo, y
resuelve detalladamente los flujos de calor y masa en la instalación a lo largo de un año
completo. Por este motivo, la información meteorológica que requiere es cuantiosa:
radiación y temperatura exterior para cada hora del año típico para la localidad
considerada. En este momento el mayor problema para su uso es que únicamente
dispone de datos meteorológicos para las capitales de provincia de Andalucía. En
España no existe aún una colección de años típicos accesible y reconocida a nivel
nacional. En todo caso, en el manual del programa se explica cómo incluir una nueva
localidad cuando se dispone de la información necesaria.

BIBLIOGRAFÍA
[1] Charles C. Smith, George O.G. Löf, Randy W. Jones, “Rates of Evaporation from
Swimming Pools in Active Use”, ASHRAE Transctions104 (1A) (1999) 514-523.
[2] ASHRAE, 1995. “AHSRAE Handbook – HVAC applications”, p.4.7. American
Society of Heating, Refrigerating and Air-Conditioning Engineers, Atlanta (USA)
[3] M. Mohammed Shah, “Prediction of evaporation from occupied indoor pools”,
Energy and Buildings 35(2003) 707-713
[4] TYPE344 manual, disponible en http://sel.me.wisc.edu/trnsys/components/type144-
manual.pdf a 6 de noviembre de 2006.
[5] ASHRAE, 2005. “AHSRAE Handbook – HVAC fundamentals”, American Society
of Heating, Refrigerating and Air-Conditioning Engineers, Atlanta (USA)

27
[6] Duffie J.A., Beckman W.A., “Solar engineering of thermal processes” Third
Edition. John Wiley and Sons, 2006
[7] Klein S.A., Duffie J.A., Beckman W.A., “TRNSYS, a transient system simulation
program. User’s Manual (version 16)”. Solar Energy Laboratory - University of
Wisconsin-Madison
[8] RETSCREEN, disponible en http://www.retscreen.net/es/home.php a 6 de
noviembre de 2006
[9] PACSOL, disponible en
http://www.agenciaandaluzadelaenergia.es/agenciadelaenergia/nav/com/contenido.jsp?p
ag=/contenidos/aplicaciones/Piscinas_cubiertas_mas_acs
a 6 de noviembre de 2006

28