3.modelado Matemático y Dimensionado de Captadores
3.modelado Matemático y Dimensionado de Captadores
3.modelado Matemático y Dimensionado de Captadores
Clima
Medio
Ambiente
Costos
de la
instalacin
Carga
Trmica
Suelo /
Geologa
Dimensionado
Instalacin
Emplazamiento
Bomba de
Calor
20
Dimensionado de instalaciones
3.2. Propiedades trmicas del suelo
Las propiedades ms importantes del terreno que influyen en un sistema geotrmico de
manera directa son la conductividad trmica, la resistividad trmica de las superficies de
intercambio y la difusividad trmica.
La ley de Fourier explica la transferencia de calor por conduccin. De esta forma, la tasa de
transferencia de calor por conduccin en una direccin dada, es proporcional al rea normal a
la direccin del flujo de calor y al gradiente de temperatura en esa direccin.
(3.1)
donde:
es el salto trmico, K.
K
Otra propiedad muy importante es el calor especfico (J/KgK), que se define como la
cantidad de calor que hay que suministrar a la unidad de masa del sistema considerado para
elevar su temperatura en una unidad (kelvin o grado Celsius) a partir de una temperatura
dada, en general, el valor del calor especfico depende de dicha temperatura inicial.
Atendiendo a esta definicin el suelo se calentar tanto ms cuanto menor sea su calor
especfico. El suelo tiene elementos constitutivos con diferente calor especfico, lo que hace que
la textura del suelo juegue un papel importante en este sentido.
La humedad supone un contenido en agua del suelo. El calor especfico del agua es del
orden de cinco veces el calor especfico de los constituyentes del suelo, por lo que el agua se
calentara ms lentamente que el suelo. Y como resultado de esto, un suelo con alto contenido en
agua se calentara ms lentamente que uno con bajo nivel de humedad.
La difusividad trmica (m2/s) se define como el ratio entre la conductividad y el
almacenamiento trmico del terreno, relacionndose con el calor especifico en la siguiente
expresin:
(3.2)
p
donde:
21
Dimensionado de instalaciones
Al producto de la densidad y el calor especfico se le llama capacidad trmica y representa
la capacidad del material para almacenar calor en condiciones de transicin en su entorno.
Valores de conductividad y difusividad para suelos tpicos se muestran a continuacin:
Conductividad Trmica
(W/mK)
0,77
1,67
1,11
0,91
2,50
1,67
Tipo de suelo
Arena
Limo
Arcilla
Marga
Arena saturada
Arcilla saturada
Difusividad
Trmica (cm2/s)
0,0045
0,0054
0,0049
0,0093
0,006
Un valor que se debe tener en cuenta, no siendo tanto una propiedad del terreno pero si un
parmetro trmico que influye directamente en un sistema geotrmico, es la resistividad trmica
entre el fluido caloportador y las paredes del pozo de intercambio. Para una tasa de transferencia
de calor concreta q (W/m), la diferencia de temperaturas entre el fluido caloportador y las
paredes del pozo de intercambio se define como:
f
donde:
(3.3)
3.2.1. Ensayos de medida. TP02 Sonda Non steady state para la medida de la
conductividad trmica
Para la obtencin de la conductividad trmica existen ensayos que permiten su clculo, la
difusividad trmica en cambio se puede calcular a partir de la conductividad.
TP02 es una sonda de alta calidad que ofrece la posibilidad de realizar una prctica y rpida
medida de la conductividad trmica del medio en el que se inserta. sta trabaja en el marco de
los estndares ASTM (American Society for Testing and Materials) D 5334-92, D 5930-97 y
IEEE (Institute of Electrical and Electronics Engineers) 442-1981.
El mtodo de medida NSSP (Non-steady-state probe) [3] posee la fundamental ventaja de
ser rpida y absoluta mientras las dimensiones de la muestra no sean crticas.
22
Dimensionado de instalaciones
La sonda TP02 Non-Steady-State consiste en una aguja con 2 enlaces de termopar (uno de
los cuales acta como referencia) y un cable que acta como foco caliente. Se inserta en el
medio que se investiga. En la base, se monta un sensor de temperatura. Las ventajas de este
diseo: una precisin ptima independiente de la temperatura media, una sensibilidad mnima a
gradientes trmicos, un sensor de alta estabilidad y la posibilidad de utilizar cables y conectores
normales.
El principio del NSSP se basa en la propiedad de la lnea de salida: tras un corto periodo
transitorio el au ento de te peratura, , solo depende de la fuente de calor,
, y de la
donde:
ln t
(3.4)
tiempo, s.
constante.
En base a esto la conductividad trmica puede ser calculada midiendo la fuente de calor y
siguiendo la evolucin de la temperatura en funcin del tiempo (2 minutos en un ensayo tpico
TP02).
23
Dimensionado de instalaciones
trmica elevada y otro material con conductividad trmica baja. Despus del periodo transitorio
observa un comportamiento lineal de la tendencia. En este periodo estacionario la pendiente de
la curva es inversamente proporcional a la conductividad trmica K.
ASTM D 5334-92 y D 5930-97
IEEE Estndar 442-1981
ango
0.1 a 6 W/mK
Sensi ilidad
Termopar K, ANSI MC96.1-1982
Rango de temperatura (sensor completo)
-55 a +180 C
Precisin (20 C)
(3% + 0.02) W/mK
Dependencia de la precisin con la temperatura
0.02 %/K (adicional)
Duracin del ciclo de medicin
200 s (tpico)
Necesidades de potencia
3 V, 1 W (mx.)
Sensor de temperatura (T)
Pt1000 Clase B, IEC 751:1983
Proteccin de la aguja y la base
IP 68
Proteccin total sensor:
IP 67
Mtodo
Ao
Objetivo
1948
Kavanaugh
1985
24
Dimensionado de instalaciones
1987
Eskilson
1965
Esta teora fue la base para los modelos de Kalman, Penrod e Ingersoll-Plass. La ecuacin
integral bsica usada en este modelo es:
e
(3.5)
25
Dimensionado de instalaciones
donde:
(3.6)
t
donde:
es el tiempo transcurrido, s.
Ingersoll y Plass en 1948 formularon el siguiente modelo como solucin al modelo integral
bsico:
,
s
donde:
log
(3.7)
es el tiempo transcurrido, s.
evolucin de la temperatura del fluido a lo largo del tiempo, K.
s
q
siendo:
3.3.2.
donde:
y,t
cos
tt
y
s
(3.8)
es la profundidad, m.
26
Dimensionado de instalaciones
es la amplitud de la oscilacin anual de la temperatura exterior, K.
t
es el tiempo transcurrido, s.
y,t
(3.9)
(3.10)
Estas son las temperaturas mxima y mnima utilizadas en el epgrafe 3.4.1.1 y 3.4.1.2 para
el clculo de la longitud del intercambiador enterrado.
27
Dimensionado de instalaciones
Desarroll el modelo considerando una sola tubera aislada rodeada por un suelo infinito
con las propiedades constantes. Asume conduccin pura en un contacto perfecto entre el terreno
y la tubera. Como ya se ha apuntado, desprecia las interferencias de pozos adyacentes.
La solucin del modelo cilndrico para un flujo de calor constante se basa en el trabajo de
.p
(3.11)
t
r
r
r
donde:
p
r0
es la distancia a la tubera, m.
es el tiempo, s.
28
Dimensionado de instalaciones
Finalmente, estas respuestas muestran curvas de T frente a t, si se adimensionalizan, se
obtendr lo que se denomina g-function. Se puede ver un ejemplo a continuacin, donde se ve
como varia la g-function en funcin del tiempo:
Una vez conocida la respuesta a una funcin escaln, la respuesta a cualquier excitacin
puede ser determinada por la descomposicin en una serie de funciones escaln unitario (figura
3.5).
Luego, utilizando los factores de respuesta (g-functions) a las funciones de cada escaln
unidad puede ser determinada la respuesta global.
29
Dimensionado de instalaciones
t r
,
ts
g
donde:
eff
(3.12)
ts
, horas.
eff
tn tn r
,
ts
eff
(3.13)
(3.14)
eff
eff
donde:
(3.15)
a
30
Dimensionado de instalaciones
rp
a
ln
ln
rp
(3.16)
rp r r
rp
donde:
r
ln
rp
ln
ln
rp
(3.17)
donde:
(3.18)
(3.19)
(3.20)
3.4.Dimensionado de captadores
Se va a explicar el mtodo de IGSHPA para el dimensionado de los captadores en las
estaciones de invierno (calefaccin) y verano (refrigeracin).
donde:
Bomba de calor
Pa Pc
Local a climatizar.
(3.21)
31
Dimensionado de instalaciones
Pa
Pc
Si este flujo de calor se expresa por unidad de longitud de intercambiador enterrado se tiene:
donde:
qc
c
(3.22)
qc
qc
donde:
(3.23)
32
Dimensionado de instalaciones
La ecuacin (3.23) se obtiene aplicando la ley de Fourier en el sistema representado en la
figura 3.7. RT se calcula como:
t
donde:
u,c
(3.24)
RS
Fu,c
La resistencia trmica del suelo (Rs) es una funcin del tiempo desde cuando el terreno cede
o absorbe calor, del radio afectado (distancia del punto del terreno hasta la tubera) y la
difusividad trmica del mismo (definida en el epgrafe 3.2)
El factor de utilizacin (Fu,c) es la fraccin de tiempo en el cual la bomba de calor
permanece en funcionamiento. Su clculo se deriva de las simulaciones que determinan las
cargas trmicas. En cada edificio, y dependiendo de cada clima y sistema de climatizacin, el
factor de utilizacin varia.
La resistencia trmica de la tubera se puede expresar como:
donde:
ln
Rt
kt
Di
De
(3.25)
Pc
Pa
Pc
(3.26)
(3.27)
s
33
Dimensionado de instalaciones
Pa
Pc
t
f
u,c
(3.28)
Si adems se toman las temperaturas mximas para obtener el caso extremo, TMAX del fluido
y TH (high temperature) del suelo, y se expresan las potencias en funcin del COP de la bomba
de calor, se obtienen las expresiones que describe el mtodo IGSHPA:
Pc
(3.29)
Pa
es el coefficient of performance en rgimen de refrigeracin de la bomba
Pc
Pc
donde:
de calor, adimensional.
Pa Pc Pc
Pc
c
donde:
Pc
Pc
Pc
Pc
(3.30)
u,c
(3.31)
La temperatura mxima del suelo se obtendr en epgrafes posteriores a partir del modelo
transitorio del sistema, con el modelo IGSHPA de Kusada & Achenbach.
donde:
Ph - Pa
(3.32)
34
Dimensionado de instalaciones
Pa
Ph
Si este flujo de calor se expresa por unidad de longitud de intercambiador enterrado se tiene:
donde:
qh
h
(3.33)
qh
En el intercambiador se tiene,
qh
donde:
qh
(3.34)
(3.35)
Ph
Pa
s
(3.36)
f
35
Dimensionado de instalaciones
y sustituyendo (3.24) en (3.36), se obtiene la longitud necesaria de intercambiador enterrado,
(3.37), donde Rs y Fu se definieron en el epgrafe de refrigeracin,
Ph
Pa
u,h
(3.37)
Si adems se toman las temperaturas mnimas para obtener el caso extremo, T MIN del fluido
y TL (low temperature) del suelo, y se expresan las potencias en funcin del COP de la bomba
de calor, se obtienen las expresiones que describe el mtodo IGSHPA:
Ph
Pa
Ph
donde:
Ph
(3.38)
Ph Pa Ph
Ph
h
donde:
Ph
Ph
(3.39)
u,h
(3.40)
La temperatura mnima del suelo se obtendr en epgrafes posteriores a partir del modelo
transitorio del sistema, con el modelo IGSHPA de Kusada & Achenbach.
La longitud que se emplear en el sistema ser la mayor de las obtenidas en regmenes de
verano e invierno, para cumplir as las especificaciones en las dos estaciones.
36