Libro Fluidos Practico
Libro Fluidos Practico
dmELECT
DMELECT, S.L.
C/General Alvear, 4, 3º B
04800 Albox (Almería)
Tlfno: 950 120757; Fax: 950 120891
[Link]
e-mail:info@[Link]
INDICE
1.2.2. Bombas.
1.2.3. Válvulas.
Figura 1.4. Intercambiador contracorriente. Figura 1.5. Intercambiador equicorriente. Figura 1.6. Radiador.
En la figura anterior pueden establecerse las siguientes relaciones, en un punto x con dA:
dq = U ⋅ dA ⋅ ∆Tx (1.1)
dq = −Ch ⋅ dTh
dq = Cc ⋅ dTc
Siendo:
−1 1
dTh − dTc = dq ⋅ −
Ch Cc
−1 1 1
− =
Ch Cc C
De donde:
dq = C ⋅ d (∆T ) (1.2)
U ⋅ dA ⋅ ∆Tx = C ⋅ d (∆T )
U ⋅ dA d (∆T )
=
C ∆Tx
Integrando:
U⋅A ∆T 2
= ln
C ∆T 1
Y empleando (1.3):
∆T 2 − ∆T 1
q = U ⋅ A⋅
∆T 2
ln
∆T 1
O también:
∆T 1 − ∆T 2
q = U ⋅ A⋅ (1.4)
∆T 1
ln
∆T 2
∆T 1 − ∆T 2
∆Tm = (1.5)
∆T 1
ln
∆T 2
P v²
+z+ (1.6)
ρ⋅g 2⋅ g
Siendo:
v max = 2
v max ²
hc =
2 ⋅ 9,81
hc = 0,204mca
P
H= +z (1.7)
ρ⋅g
La relación entre la pérdida de carga (h) para una rama que va del nudo i al j y el caudal de
paso(Q), puede escribirse:
Hi − Hj = hij
∑ (Qij − Di ) = 0 (1.9)
j
Para todo i, siendo i los nudos de caudal de nuestra red, y siendo Di el caudal de demanda en el nudo
i, el cual por convenio se toma positivo cuando entra en el nudo. En conducciones cerradas a
presión no abiertas a la atmósfera lógicamente Di = 0
Dados los valores de las alturas en los nudos de altura conocida (en las redes de calefacción
son los nudos donde se fija la presión de llenado y en redes abiertas a la atmósfera son los
embalses y depósitos), se trata de encontrar una solución para las alturas Hj en los restantes
nudos y para los caudales Qij de todas las ramas o líneas, que satisfagan las ecuaciones 1.8 y
1.9.
- Tuberías, Darcy-Weisbach.
L v²
hf = f ⋅ ⋅ (1.10)
D 2⋅ g
Siendo:
v ⋅ π ⋅ D² 4000 ⋅ Q
Q= o v= (1.11)
4000 π ⋅ D²
Siendo:
v = Velocidad en m/s.
D = Diámetro interior en mm.
Q = Caudal en l/s.
10 9 ⋅ 8 ⋅ f ⋅ L
hf = ⋅ Q² (1.12)
π ² ⋅ g ⋅ D5
Siendo:
Por último, el factor de fricción "f" se calcula de diferentes formas según el tipo de flujo, que a
su vez depende del Re:
64
f = (1.13)
Re
0,25
f = 2
(1.14)
ε 5,74
log + 0,9
3,7 ⋅ D Re
- 2000 < Re < 4000: Se emplea una interpolación cúbica del diagrama de Moody:
- Tuberías, Hazen-Williams.
En el cálculo de las pérdidas por fricción o lineales también se utilizan fórmulas empíricas
como la de Hazen-Williams (sólo para agua), de amplia aceptación en los paises anglosajones.
12.171 ⋅ 10 9 ⋅ L 1,852
hf = ⋅Q (1.16)
C 1,852 ⋅ D 4,871
Siendo:
v²
hl = k ⋅ (1.17)
2⋅ g
Siendo:
10 6 ⋅ 8 ⋅ k
hl = ⋅ Q² (1.18)
π ² ⋅ g ⋅ D4
- Resumen Tuberías.
Fórmula Coef. Resistencia (r) Exp. Caudal (n) Coef. Pérdidas menores (m)
10 9 ⋅ 8 ⋅ f ⋅ L 10 6 ⋅ 8 ⋅ k
Darcy-Weisbach 2
π ² ⋅ g ⋅ D5 π ² ⋅ g ⋅ D4
12.171 ⋅ 10 9 ⋅ L 10 6 ⋅ 8 ⋅ k
Hazen-Williams 1,852
C 1,852 ⋅ D 4,871 π ² ⋅ g ⋅ D4
Tipo Accesorio k
No obstante, las pérdidas localizadas pueden introducirse como un porcentaje de las lineales (10 % -
20 %).
1.2.2. Bombas.
Son elementos que comunican energía al fluido elevando su altura piezométrica:
Qij
n
hij = −ω ² ⋅ h0 − rb ⋅ (1.19)
ω
Siendo:
hij ⋅ Qij
P = 9,81 ⋅ (1.20)
η
Siendo:
P = Potencia en W.
η = Rendimiento en tanto por uno.
1.2.3. Válvulas.
Son elementos que limitan o regulan la presión o el caudal en algún punto de la red.
- Válvulas de operación.
Permiten el mantenimiento de la red y actúan todo o nada (on_off). Son las válvulas de
corte o seccionamiento.
10 9 ⋅ 8 ⋅ f ⋅ L 10 6 ⋅ 8 ⋅ k
π ² ⋅ g ⋅ D5 π ² ⋅ g ⋅ D4
f = 0,02 Tipo de válvula k0 (Coef. pérdidas abierta)
L = 2·D/1000
Compuerta 0,1-0,3
Mariposa 0,2-0,6
Bola 0,1
Asiento 3,0-10,0
- Válvulas de regulación.
Permiten el control de una variable hidráulica, normalmente el caudal. Son las Válvulas de 2
vías o 3 vías (mezcladoras), que reciben órdenes a distancia mediante señal eléctrica o
termostática.
- Válvulas de protección.
Protegen los elementos de la red. Son las válvulas de alivio de presión o de seguridad,
válvulas de retención, etc.
VPG Las válvulas de propósito general se comportan como una tubería, permitiendo modelar
elementos cuya pérdida de carga varía con el caudal (separadores de aire, interacumuladores,
intercambiadores de calor independientes, turbinas, etc).
DATOS
CUESTIONES
1. Diseñar los emisores teniendo en cuenta las condiciones del local: metros para el PC 600 y
número de elementos en los DUBAL 60, teniendo en cuenta que el PV 1500/300 es fijo.
2. Diseñar las tuberías, circulador y demás elementos de la instalación considerando una vmax= 2
m/s y hu = 40 mmca/m (pérdida unitaria).
3. Analizar el equilibrado de la red y calcular el vaso de expansión.
elemento con una diferencia de Tª media entre el emisor y el ambiente de 50 ºC, según
consideraciones de UNE-EN 442 y que debe proporcionar el fabricante). Según esto:
i = 0.....5 i = 0.....5
781 908,303
103,9 120,836
103,9 120,836
Pce50 = en kcal/h ⇒ Pce50 = Pce50 · 1,163 ⇒ Pce50 = en W
500 581,5
103,9 120,836
103,9 120,836
En el vector anterior el primer elemento es el modelo PC 600 -panel horizontal de acero-, el
cuarto es el modelo PV 1500/300 -panel vertical de acero- (se trata de un radiador fijo, como los de
toallero, sin posibilidad de diseñar elementos), el resto son del modelo DUBAL 60 (FA).
Por otra parte, n representa el exponente de la curva característica del emisor, necesario ya
que Pce = Pce50(∆t/50)n.
1,3
1,35
1,35
n=
1,3
1,35
1,35
te = 75 St = 10 ts = te - St ts = 65 ta = 21
Dado que el sistema bitubo conecta los radiadores en paralelo se tendrá la misma Tª de entrada y salida
en todos nuestros emisores térmicos, luego:
Las ecuaciones anteriores, puesto que ta la consideramos cte, también se podrán expresar:
te + ts te − ts
∆t1 = − ta ∆t 2 =
2 ∆te
ln
∆ts
Resultando:
Al ser a ≥ 0,7 tomamos por tanto la diferencia de Tª media (∆t1 = 49), con lo cual:
884,758
117,585
ni 117,585
∆t
Pcei = Pce50 i ⋅ Pce =
50 566,427
117,585
117,585
Dado que el cuarto elemento es un radiador fijo, la potencia del local menos la potencia de ese radiador
(6520 W - 566,427 W) se la repartirán el resto de radiadores, y así diseñaremos los elementos de cada
uno.
1190,715 1,346
1190,715 10,126
1190,715 Pcri 10,126
Pcr = ⇒ nei = ⇒ ne =
566,427 Pcei 1
1190,715 10,126
1190,715 10,126
Si analizamos los resultados llegamos a la siguiente conclusión: el primer modelo -PC 600-
tiene como unidad el m (1000 mm), así que el superior en base de datos de nuestro fabricante será
1,35 m (1350 mm), el cuarto es nuestro panel vertical fijo -PV 1500/300-, y los demás -DUBAL
60(FA)- son de aluminio de elementos verticales y se necesitarán 11 (el entero superior). No
obstante, de esta forma se mayora mucho, y más en este caso donde ne esta más cerca del entero
inferior que del superior. DmELECT SL tiene una forma ingeniosa de resolver esta situación: a la Pc
(6520 W, potencia calefacción local) le resta los elementos fijos y del resto (configurables) va uno a
uno aproximando al entero superior que exista en la base de datos del fabricante, de este manera lo
que va quedando de potencia de calefacción a satisfacer lo va repartiendo entre los siguientes y así
sucesivamente hasta obtener una optimización en el material necesario para aportar la potencia
de calefacción en cada local; lógicamente esto lo hace local a local, pues cada uno tiene sus
necesidades térmicas, de esta forma quedará:
1,35 1194,42
11 1293,43
10 1157,85
ne = ⇒ Pcri = Pcei · nei ⇒ Pcr = ⇒ Pclocal = ∑ Pcr j = 6581,82 W
1 566,43 i
10 1175,85
10 1175,85
No obstante lo expuesto, el proyectista puede imponer cualquier número de elementos y por
tanto cualquier potencia en los radiadores.
En primer lugar debemos determinar el caudal que debe circular por cada radiador para
proporcionar su potencia calorifica, asignada en función del salto térmico:
0,0285
0,0309
Pcri 0,0281
Qri = ⇒ Qr = en l/s (estando Pcr en W)
4186 ⋅ St 0,0135
0,0281
0,0281
Por otra parte, el caudal de la bomba será:
1/ Para el dimensionado de tuberías y válvulas debemos conocer el caudal (Qij) por cada una de
ellas, esto supone que por cada radiador pase el caudal necesario para aportar la potencia térmica;
esto sólo es posible si dichos emisores disponen de detentores o válvulas limitadoras de caudal
(VLQ) taradas al caudal necesario (sin estos elementos de regulación, al ser un sistema en paralelo
pasaría más caudal por los recorridos de menor resistencia hidráulica). Esto se conoce como el
equilibrado hidráulico en un circuito cerrado, recirculación.
2/ El diámetro interior óptimo será el mayor entre los siguientes: Dmín, Dvmáx, y Dhumax.
Criterio Dmín
i = 0.....6
15 16,1
20 16
20 16
Dn min = 20 Di min =16
20 16
20 16
20 16
Donde Dnmin y Dimin representan el diámetro nominal e interior mínimo respectivamente. El primer
elemento del vector corresponde a una válvula, su diámetro nominal (valor medio) es inferior al interior y
se designa en pulgadas (1/2").
Criterio Dvmax
El siguiente vector -Ql- corresponde a una válvula de corte y 6 tuberías que transportan el
caudal de los radiadores de nuestro ejemplo. En realidad, este proceso de cálculo hay que plantearlo
para las n tuberías de nuestra red, pero por comodidad de explicación lo hacemos para esta muestra.
0,1572
0,1572
0,1287
Ql = 0,0978
0,0697
0,0562
0,0281
Criterio Dhumax
16,1
18,189
16,875
Di =16
16
16
16
Los diámetros interiores y nominales normalizados, así como la designación UNE serán:
Finalmente, con estos diámetros podemos verificar la pérdida de carga y la velocidad real en tuberías.
El primer término del vector de pérdidas (hf) representa la pérdida lineal en la válvula, para la cual se
considera fijo f = 0,02 y L = 2D/1000; lógicamente éste término en una válvula es menos importante
que el debido al coeficiente de pérdidas localizadas, m. Según esto:
k = 0,5
10 6 ⋅ 8 ⋅ k
hv = 0,001 + 2 ⋅ 0,1572 2 ⇒ hv = 0,016
π ⋅ 9,81 ⋅ 16,1 4
Conclusión
1/ Puede observarse que en la justificación manual se obtienen los mismos resultados que con el
Software dmELECT, no obstante hay que señalar que para no perder claridad en la exposición hemos
tomado una pequeña muestra de tuberías y válvulas; así mismo hemos asumido directamente el
vector de factores de fricción " f " para cada tubería, cuando sabemos que depende del Re y éste del
diámetro (el cual buscamos), por tanto habrá que iterar hasta que no se modifique ni " f " ni " Di ".
Qb = 0,1572 l/s
hb = 1,5 mca
Por tanto:
Fig. 1.9. Curva de la bomba del ejemplo (generada por el software de DMELECT)
Otro punto fundamental en los vasos de expansión es que, para su correcto funcionamiento,
la presión inicial en frío en el lado del gas deberá ser igual a la diferencia de cotas entre el punto
mas alto y la posición del vaso más 5 mca (0,5 bar) - 10 mca (1 bar) por encima de la presión
atmosférica, para evitar entrada de aire en la red (UNE 100155).
Vu = Vinstal ⋅ Ce
Siendo:
1 d
β= ⋅ V ⇒ dV = β ⋅ dT ⋅ V ⇒ Ce = β ⋅ dT
V dT
( )
Ce1 = − 1,75 + 0,064t + 0,0036t 2 ⋅ 10 −3 para Tª ≤ 70 ºC
Ce2 = (− 33,48 + 0,738t ) ⋅ 10 −3 para Tª > 70 ºC
Para calcular Vinstal debemos acudir a la medición del software dmELECT SL para ver los elementos
con los que contamos, si no corremos el riesgo de olvidar alguno:
l = 0...3
5,94
i = 0...1 j = 0...1 k = 0...1
14,6 16 30,8 16,1 3,96
G = Dit = Lreal = Div = Ut =
10,8
0 20 32,39 0
7,8
Por lo tanto:
Por último, hay que señalar que estos 1,302 litros se introducirán en el vaso disminuyendo el
volumen de la cámara que ocupa el gas, aumentando por tanto su presión y con ello la de toda
nuestra red hidráulica, la cual viene limitada por nuestra válvula de seguridad Pvs. Si consideramos
que el gas evoluciona isotérmicamente, se tendrá (ver Figura 1.10):
P min
P max⋅ V min = P min⋅ V max ⇒ V min = V max⋅ ⇒ Vu = V max − V min
P max
P max
= Cp
P max − P min
Donde Cp es el Coeficiente o factor de presión. Así mismo las presiones son absolutas o
termodinámicas.
En nuestro ejemplo:
Por lo tanto:
El volumen total nominal será el superior de base de datos al Vtc (Volumen total de cálculo), en este caso
12 litros. Como puede apreciarse los cálculos coinciden con los desarrollados con el software de
dmELECT.
De donde:
Finalmente:
i = 0...5 i = 0...5
781 908,303
103,9 120,836
103,9 120,836
Pce50 = en kcal/h ⇒ Pce50(W)=Pce50·1,163 ⇒ Pce50 = en W
500 581,5
103,9 120,836
103,9 120,836
En el vector anterior el primer elemento es el modelo PC 600 -panel horizontal de acero-, el
cuarto es el modelo PV 1500/300 -panel vertical de acero- (se trata de un radiador fijo, como los de
toallero, sin posibilidad de diseñar elementos), el resto son del modelo DUBAL60 (FA).
Por otra parte, n representa el exponente de la curva característica del emisor, necesario ya
que Pce = Pce50(∆t/50)n:
1,3
1,35
1,35
n=
1,3
1,35
1,35
Pcri Pcri
tei +1 = tei − ts i = tei −
4186 ⋅ Q 4186 ⋅ Qri
75
75 69,795
73,178
71,357 73,178 68,624
71,357 66,803
te = 69,535 te = ts =
68,643 69,535 66,987
68,643 64,089
66,822
66,822 62,268
65
54 48,795
52,178 47,624
50,357 45,803
∆te = ∆ts =
48,535 45,987
47,643 43,089
45,822 41,268
Si a ≥ 0,7 el valor buscado será ∆t1, caso contrario ∆t2. Por lo tanto:
Análisis
1. La potencia emitida por elemento Pce cambia de radiador a radiador aunque sean del mismo
modelo, ya que al estar en serie la te y ts es distinta en cada uno de ellos.
2. Existen muchas formas de optimizar la elección de elementos por radiador para no
sobredimensionar la instalación, pero en todas se necesita un proceso de iteración por varias
razones. A la potencia de calefacción del local o espacio térmico hay que restar la del radiador fijo, a
continuación fijamos el ne del primer radiador del anillo al inmediato superior y el resto de potencia
del mismo local se reparte entre los demás radiadores, a continuación volvemos a calcular te y ts
(pues dependen de Pcr) y así sucesivamente hasta que el sistema converja, y siempre comprobando
que la suma de la potencia emitida por los radiadores del local es mayor o igual que la necesaria.
3. El resto de elementos (tuberías, válvulas, bomba y vaso de expansión) se calculan exactamente
igual al bitubo.
4. El caudal total es el mismo en el sistema bitubo y monotubo, pues combaten la misma carga
térmica del local.
5. Los caudales en los radiadores del sistema monotubo son mayores que los del bitubo (salto
térmico menor).
6. Las pérdidas hidráulicas en el sistema monotubo (serie) son mayores que en el bitubo (paralelo).
7. Precisa menos tubería el sistema monotubo que el bitubo.
Figura 1.15. Vista en 3D generada por el Software de DMELECT. Colector y distribución de tuberías.
DATOS
1. Consideramos un local con unas necesidades de calefacción de Pc = 5700 W y ta =21 ºC, estando
el material del pavimento o solera emisora constituido por baldosa de 2 cm.
CUESTIONES
1. Diseñar los emisores - diámetro de tuberías - teniendo en cuenta una vmax = 2 m/s y hu = 40
mmca/m, así como la temperatura de impulsión.
i = 0 ..... 2 h = 11 W/m²·ºC St = 5 ºC
1900 22,107
P = 1900 S = 22,107
1900 22,107
El primer vector señala la carga térmica a combatir por cada unidad terminal y el segundo el
área correspondiente.
28,81
tsoi = ∆Tsa i + tai ⇒ tso = 28,81
28,81
Siendo:
Por otra parte, considerando Rse como la resistencia de la solera emisora compuesta por tubo PEX
(0,006) más la solera de hormigón y el pavimento:
0,062 5,329
P ⋅ Rsei
Rse = 0,062 ∆Tmas i = i ⇒ ∆Tmas = 5,329
0,062 Si 5,329
Siendo:
St
tmai = ∆Tmasi + tsoi tia i = tmai +
2
34,142 36,642
tma = 34.142 tia = 36.642
34,142 36,642
Siendo:
Para obtener los caudales por cada unidad terminal (Qi, l/s) se utilizará:
0,0908
Pi
Qi = ⇒ Q = 0,0908
4186 ⋅ St 0,0908
0, 2
1012 ⋅ 8 ⋅ f ⋅ Q 2
Dc = 2
π ⋅ 9,81 ⋅ hu max
Siendo:
Figura 1.18. Vista en 3D generada por el Software de DMELECT. Fan Coils Cassette 2 tubos.
DATOS
1. Consideramos un local con unas necesidades de refrigeración de Pft = 10909 W, Pfs = 5773 W y
ta =24 ºC.
CUESTIONES
1. Diseñar los Fan_coils y diámetro de tuberías teniendo en cuenta una vmax = 2 m/s y hu = 40
mmca/m.
Una vez establecido el diseño térmico (elección de la potencia total de las UT) estamos en
disposición de realizar el dimensionado hidráulico de nuestras unidades terminales de fan_coils el
cuál sigue el mismo criterio utilizado en suelo radiante.
11460
i = 0 ..... 1 St = 5 ºC Pf := Siendo Pfi la potencia total frigorífica en W.
5730
Pi 0,548
Qi = ⇒ Q = Siendo Qi los caudales l/s por cada unidad terminal.
4186 ⋅ St 0, 274
0, 2
0,032 1012 ⋅ 8 ⋅ f i ⋅ (Qli 2 )
f = hu max = 40mmca / m Dcf i = 2
0,036 π ⋅ 9,81 ⋅ hu max
Siendo:
Finalmente en este ejemplo tan sencillo (pero representativo del modelo) se han seleccionado
estas dos tuberías, la primera para el caudal total y la segunda para alimentar una unidad terminal, el
material es Multicapa Pol/Al/PE_X y la denominación será 40x3,5 y 32x3 mm.
En primer lugar para identificar el problema señalaremos lo que indica el CTE sección HS_4,
suministro de agua del DB_HS apartado [Link], en el punto 3. "Tanto en instalaciones individuales
como en instalaciones de producción centralizada, la red de distribución de ACS debe estar
dotada de una red de retorno cuando la longitud de la tubería de ida al punto de consumo más
alejado sea igual o mayor que 15 m".
También es muy interesante lo que señala en el punto 6 del mismo apartado." En los
montantes debe realizarse el retorno desde su parte superior y por debajo de su última
derivación particular. En la base de dichos montantes se dispondrán válvulas de asiento para
regular y equilibrar hidráulicamente el retorno".
DATOS
CUESTIONES
Para determinar el caudal que circulará por el circuito de retorno, se estimará que en el grifo
más alejado, la pérdida de temperatura sea como máximo de 3ºC desde la salida del acumulador o
intercambiador en su caso.
Este es el método exacto y para resolverlo consideramos la red de ACS como un emisor
térmico (como una unidad terminal semejante al suelo radiante).
Siendo:
Planteando la ecuación de Fourier para conducción del calor y la ley de Newton para una
sección de tubería:
q d
= − λ ⋅ T Ecuación de Fourier q = h ⋅ S ⋅ dT Ecuación de Newton
2 ⋅π ⋅ L ⋅ r dr
Integrando y desarrollando se obtiene la transmitancia por la superficie para una tubería
circular:
2πL
US = (1.7.2)
req ln reso
ln
ret + req + 1
λq λso he ⋅ reso
Siendo:
Aplicando la ecuación (1.7.2) por unidad de longitud (L) a nuestro ejemplo se tiene:
det
λq = 0,05 λso = 1,975 he = 11 det = 22 ret = ret = 11
2
eq = 9 mm req = ret + eq req = 20 mm eso = 50 + 20
2π
USL = USL = 0,4577 W / m K
req ln reso
ln
ret + req + 1
λq λso he ⋅
reso
1000
Podemos simplificar el problema considerando sólo la coquilla de la tubería ya que esta es la
que proporciona mayor aislamiento:
2π
USL = USL = 0,5255 W / m K
req
ln
ret
λq
L = 20 m q = qL ⋅ L q = 367,845 W
Dado que esta potencia calorífica debe ser aportada por un caudal de agua con un salto
térmico de 3ºC:
q
Q(l / s ) = Q = 0,0293 l/s
4186 ⋅ 3
Así mismo el DB_HS4 contempla la posibilidad de estimar el caudal de retorno según reglas
empíricas de la siguiente forma:
½ 140
¾ 300
1 600
1¼ 1100
1½ 1800
2 3300
Aplicando este segundo criterio a nuestro caso, caudal instalado en ACS 0,295 l/s, el 10 %
será 0,0295 l/s, valor muy semejante al obtenido por el método exacto.
Una vez efectuado el diseño de la red de retorno, el cálculo hidráulico es exactamente igual al
efectuado en los casos anteriores típico de una red cerrada no abierta a la atmósfera.
Análisis
1. En este tipo de redes de retorno es interesante comprobar con el programa dmELECT, el hecho de
que si conectamos la primera derivación particular y la última a la red de retorno se produciría un
cortocircuito hidráulico, no retornando agua desde la última derivación particular.
2. También se puede comprobar que si no ponemos una válvula limitadora de caudal que modeliza el
equilibrado hidráulico se establecerá el caudal que imponga la resistencia hidráulica.
La última red cerrada no abierta a la atmósfera que vamos a estudiar es la constituida por
captadores solares térmicos para dar cumplimiento a la sección HE_4 (Contribución solar mínima
de agua caliente sanitaria) del DB_HE (ahorro de energía) del CTE.
Figura 1.23. Vista en 3D generada por el Software de DMELECT. Instalación Solar Térmica.
DATOS
2 La instalación la situamos en Albox (Almería), 37º de Latitud Norte y 420 metros de altitud. La
ubicación de los captadores en el edificio corresponde al caso general, sin superposición ni
integración. La fuente de energía del apoyo será electricidad (efecto Joule).
4 Las características hidráulicas del captador nos dice que el caudal de diseño recomendado para
2
el uso previsto es de 50 l / hm .
CUESTIONES
2 Calcular y comprobar el diámetro de las tuberías principales, teniendo en cuenta, vmax = 2 m/s y
hu = 40mmca/m.
En primer lugar vamos a expresar los datos necesarios para el cálculo en la localidad de
Albox (Almería).
i = 0 ..11 ∆z = 404 Ptub = 1,05 Rend cubierta = 0,96 Rend suciedad = 0,98
Acap = 2,23
12 0,0066
12 0,0066
13 0,0066 0
0 9,334
14 0,0033 1 9,334
16 0,0033 2 10,334
3 12,667
18 0,0033 4 14,667
tagr = Br = tagri = tagri - ∆z ⋅ Bri tagr = 5 16,667
20 0,0033 6 18,667
21 0,0033 7 19,667
8 17,667
19 0,0033 9 14,334
17 0,0066 10 11,334
11 9,334
14 0,0066
12 0,0066
12,4 0,01
13 0,01
14,4 0,01
0
16,1 0,005 0 8,36
18,7 0,005 1 8,96
2 10,36
22,3 0,005 3 14,08
tamb = Ba = tambi = tambi - ∆z ⋅ Bai tamb = 4 16,68
25,5 0,005 5 20,28
26 0,005 6 23,48
7 23,98
24,1 0,005 8 22,08
20,1 0,01 9 16,06
10 12,16
16,2 0,01 11 9,26
13,3 0,01
Siendo:
Con esta información, ya estamos preparados para abordar los cálculos; en primer lugar
determinaremos el volumen de almacenamiento de cada mes (para obtener el más desfavorable), ya
que almacenamos agua a 50ºC.
0
60 − tagri 0 622,95
D50º i = D60º i ⋅ D50º =
1 622,95
50 − tagri 2 626,05
3 633,93
El primer dato que obtenemos es que el volumen de 4 641,51
acumulación para el mes más desfavorable es de 664,84 l/dia; elegimos 5 650
un acumulador solar (serpentín sumergido), con función de depósito de 6 659,58
consumo (circuito abierto), marca Ferroli, modelo Inoxunit de 750 litros. 7 664,84
8 654,64
9 640,19
10 629,31
11 622,95
El siguiente paso es calcular la demanda de energía necesaria para nuestro uso mes a mes:
1
DEacsi = D60º i ⋅4186 ⋅ (60 − tagri ) ⋅ Ptub ⋅ ni ⋅ Mj
10 6
DEacsaño = ∑i DEacsi 0
0 3451,76
1 3117,72
DEacsaño = 37145,25 Mj/año DEacs = 2 3383,63 Mj/mes
3 3120,65
4 3088,42
5 2856,94
6 2815,91
7 2747,79
8 2791,01
9 3111,12
10 3208,55
11 3451,76
Ahora calculamos la energía que llega a nuestros captadores y que penetra a través del
cristal, según la localidad de ubicación:
El próximo paso será calcular el rendimiento medio anual de nuestro captador, para
determinar qué parte de la energía que llega a Rad 37 º se convierte en útil para calentar el agua.
50 − tambi
η i = 0,788 − 3.955 ⋅ 0
Gi 0 0,222
1 0,329
2 0,459
∑η i 3 0,524
ηmedioanual = i
η= 4 0,581
12 5 0,614
6 0,637
7 0,621
8 0,582
ηmedioanual = 0,464 9 0,456
10 0,332
11 0,21
Según el DB_HE sección HE_4, nuestra zona climática es 5 (recibe una radiación media
anual sobre superficie horizontal mayor de 5 kWh (18 Mj); dado que el apoyo será por efecto Joule la
contribución solar mínima anual será del 70 %.
En el apartado del rendimiento ya hemos hecho una primera aproximación, pues hemos
dejado sólo el término independiente (ganancia del captador) y el primer término de pérdidas y así
mismo hemos dejado las pérdidas en función de la temperatura de entrada (la de acumulación). Por
último, para diseñar el número de captadores necesario y la fracción solar en tanto por uno
mes a mes, sin superar en ningún mes el 100 %, se debe emplear un método iterativo, probando
captador a captador hasta cumplir con lo especificado; el método f_chart, es reconocido
internacionalmente y es el utilizado por dmELECT en su programa SOLTE; su justificación puede
verse en el Anexo de Cálculos del fichero Curso_Fluidos.IED que se acompaña; dado que su
desarrollo y justificación es algo tediosa y nos aparta del carácter didáctico de esta obra, aquí
emplearemos un método muy simple basado en los valores medios anuales pero suficiente para
comprender el principio de funcionamiento.
DEacsaño ⋅ 0,7
nc = nc = 3.936 tomando el entero superior:
Radcaptador 37 º⋅ηmedioanual
Radcaptador 37º⋅ηmedioanual ⋅ nc
f = f = 0,711 Fracción solar = 71,1%
DEacsaño
Vac
Vac = 750 = 84,08
Sc
A continuación pueden verse los diagramas de demanda y contribución solar mes a mes.
De los datos del captador solar plano obtenemos su longitud a = 2 metros; así mismo los
datos para el cálculo serán:
φ = 36,9 Latitud en º
β ⋅π
asen
β ⋅π 180
dcap = a ⋅ cos + dcap = 4,29 metros.
180 tan α ⋅ π
180
h
dcap _ ob = dcap _ ob = 2,24 metros.
α ⋅π
tan
180
Ce = 0,06
Para calcular el Vinsta debemos acudir a la medición del software dmELECT SL para ver los
elementos con que contamos, si no, corremos el riesgo de olvidar alguno.
π ⋅ ( Dit ) ⋅ Lreal 2
11 ⋅ π ⋅ ( Div k ) 3 ⋅ 2
Vinsta = (∑ G ) + ∑
+ ∑ + ∑ Ut l
j j
i
i
4000
j k 4 ⋅ 10 6
l
Vinsta = 32,09 litros
Vvap = 1,1 ⋅ ∑ Gi Vvap = 7,48 litros
i
Vres = 0,03 Vinsta Vres = 0,963 Vres = if (Vres > 3,Vres,3) Vres = 3
Pllenado
Pllenado = Pllenado = 1,5 bar P min = Pllenado + 1 P min = 2,5 bar
10
P max 1 = 0,9 Pvs + 1 P max 1 = 5,5 bar P max 2 = Pvs + 0,65 P max 2 = 5,65 bar
P max
Cp = Cp = 1,8333 Vtc = Vu ⋅ Cp Vtc = 22,92 litros
P max − P min
El volumen total nominal será el superior de base de datos al Vtc (volumen total de cálculo),
en este caso 25 litros. Como puede apreciarse, los cálculos coinciden con los desarrollados con el
software dmELECT SL.
50 ⋅ Ac
i = 0.. 1 Ac = 2,23 m 2 Qc = Qc = 0,031 l/s
3600
Dado que se tiene una batería con dos captadores en paralelo y un total de 2 baterías (4
captadores):
0, 2
0,06194 0,031 1012 ⋅ 8 ⋅ f i ⋅ (Qtsoltei ) 2
Qtsolte = f = hu max = 40mca Dcsoltei =
π ⋅ 9,81 ⋅ hu max
2
0,1239 0,043
Siendo:
Finalmente, en este ejemplo tan sencillo (pero representativo del modelo) se han
seleccionado estas dos tuberías, la primera para el caudal total y la segunda para alimentar una
batería de dos captadores, el material es Cu y la denominación será 20x22 mm.
Este tipo de redes presentan características distintas a las cerradas del capítulo 1, ya que
pueden intercambiar con el medio, materia (dotaciones_caudales) y también energía (trabajo de
expansión al estar en carga, calor, etc). Para comprender su dimensionamiento vamos a definir
algunos conceptos:
Acometida: Tubería que enlaza la red de distribución con la instalación general del edificio.
Montante o Ascendente: Tubería que une verticalmente el distribuidor principal (tubo de
alimentación) con las instalaciones interiores particulares o las derivaciones particulares en edificios
de propiedad múltiple.
Caudal instantáneo minimo, Qmin (l/s): Caudal instantáneo que se debe suministrar a cada uno de
los aparatos sanitarios.
Caudal total instalado, Qt (l/s): Suma de los caudales instantáneos mínimos de todos los aparatos
instalados.
Caudal simultáneo o caudal de cálculo, Qc (l/s): Caudal que se produce por el funcionamiento
lógico simultáneo de aparatos de consumo. Veremos que lo obtenemos según UNE 149201.
Conviene señalar en este punto una curiosidad de las redes de fontanería y es que en ellas no se
cumplen las leyes de Kirchhoff (en los modelos de cálculo), concretamente el principio de
continuidad; la razón es que las normas establecen simultaneidades por ramas y no por nudos, de ahí
que el caudal que entra en un nudo a través de ramas es distinto del que sale; esto hace que los
modelos de cálculo por ordenador obedezcan a sistemas en árbol (desde un tronco hacia las ramas).
Presión estática, Pe: Presión en cualquier punto de la red sin consumo alguno.
Presión dinámica, Pd: Presión en cualquier punto de la red cuando hay consumo.
Presión mínima de funcionamiento, Pmin: Presión mínima necesaria para el funcionamiento de un
aparato, a saber:
Presión de suministro, Ps: Presión manométrica mínima necesaria para un correcto funcionamiento
de la instalación.
Presión de trabajo, Pt: Valor de la presión manométrica interna máxima para la que se ha diseñado
el conducto, considerando un uso continuado de 50 años.
Presión nominal, Pn: Número convencional que coincide con la presión máxima de trabajo a 20ºC.
Presión de arranque, Pb: Presión manométrica a la que se dispone el arranque del grupo de
sobreelevación.
Presión de paro, Pp: Presión manométrica a la que se dispone el paro del grupo de sobreelevación.
Presión de consigna n (en velocidad variable), Pc: Presión manométrica que se mantiene
constante a la salida de la bomba para el caudal de cálculo.
Número de arranques/hora, máximo n: Máximo número de arranques que se le debe permitir a una
bomba para un adecuado mantenimiento.
2
Pa = N/m , Pascal es la unidad del Sistema Internacional.
1 bar = 100000 Pa = 0,1 Mpa = 100 kPa = 10,2 mca
2
1 atm = 1,01325 bar = 101325 Pa = 1013,25 HPa = 1,033 Kg/cm = 10,33 mca = 760 Torr
(1 Torr = mmcHg)
2
1 Kg/cm = 1 atm técnica = 10 mca 1 mbar = 10,2 mmca = 100 Pa
DATOS
2 La instalación será realizada con Cu (tubería lisa), vmáx = 2 m/s; cáculo de caudales simultáneos
de cálculo según UNE 149201, tipo de edificio “Centros Comerciales”. Número de arranques/hora
máximos permitidos al grupo de sobreelevación = 20. Diferencia entre la presión de paro y la de
arranque (d) 20 mca. Acumulador hidráulico de membrana. El depósito de alimentación (para no
aspirar directamente de las tuberías) se dimensionará para un tiempo previsto de utilización de 20
minutos.
CUESTIONES
1 Calcular el caudal simultáneo o de cálculo por cada rama, así como su diámetro, pérdidas y
velocidad cumpliendo con la sección HS_4.
2 Calcular todos los elementos complementarios, contador, grupo de presión con calderín
presurizado, bomba, depósito auxiliar de alimentación, etc.
i = 0..11 j = 0..3
0.15
0.1
0.1
0.1
0.1
0.065
0.1 0.065
Qf min = Qft = ∑ Qf min 1 Qacs min = Qacst = ∑ Qacs min j
0.1 0.065
i
j
0.1 0.1
0.1
0.1
0.25
0.2
Qft = 1.5 l/s Qacst = 0.295 l/s
Hemos determinado el caudal total instalado en agua fría y agua caliente correspondiene a
los tramos principales que transportan toda el agua. Para el cálculo de tuberías vamos a elegir cuatro
( el resto se calculan igual ), tres representan el tramo principal y las derivaciones a los locales ( baño
y cafetería ) de agua fría y la otra corresponde al total de ACS.
1.5 2.22
1.05 2.89
Qtinsta = ltub =
0.45 7.24
0.295 2.29
Una vez determinado el caudal total instalado por cada rama en estudio y su longitud real,
vamos a determinar el caudal simultáneo o de cálculo según UNE 149201, para centros comerciales:
0.7349
0.5952
Qc = 0.698(Qtinsta l )
0.5
− 0.12 Qc = l/s
0.3482
0.2591
El diseño en fontanería es más sencillo que en el resto de redes ( al estar muy normalizado),
de hecho, ya vimos en la tabla anterior los diámetros nominales de los aparatos ( para tubería lisa y
rugosa) así como sus caudales instantáneos mínimos apartado 4.3.1 del DB HS_4, así mismo en este
apartado punto 2 en función del tipo de tramo se tiene:
21.63
0.5
4000 ⋅ Qcl 19.47
Dcv max l = Dcv max =
π ⋅ v max 14.89
12.84
Estos diámetros lógicamente son interiores, cogiendo normalizados en Cu.
26
20
D int v max = Los mínimos interiores según tablas anteriores serán:
16
13
26 26
20 20
D int min = de donde: Dc int =
20 20
20 20
Con estos diámetros interiores normalizados (mayor del criterio de mínimo y de velocidad), ya
podemos calcular el resto de parámetros:
vl := if (l < 3, vf , vc)
1.1 × 10 −6 32715.7
1.1 × 10 −6 4 ⋅ Qc1 34449.1
v= Re l = Re =
−6
1.1 × 10 π ⋅ Dc int l ⋅ vl 20153.8
6.6 × 10 −7 24993.2
Puede observarse lo turbulento (Re > 4000 ) del régimen hidráulico incluso para estos
caudales tan pequeños.
0.025
0.25 0.026
fl = f = Donde f es el factor de fricción.
ε 5.74
2
0.028
+
log 0.9
0.027
3.7 ⋅ Dc int l (Re l )
0.252 1.38
0.81 1.89
h= mca v= m/s
0.764 1.11
0.128 0.82
DEPÓSITO DE ALIMENTACIÓN.
V = Qc ⋅ t ⋅ 60 V = 881.9 l
CALDERÍN HIDRÁULICO.
A = 23 Hp = 0 hp = 3
3600 Vu Vu 3600 Qs 2
= + Vu = ⋅ Qs −
Nc Qe − Qs Qs Nc Qe
Qe
derivando respecto a Qs para encontrar el máximo Vu para un Nc ,se tiene para Qs =
2
900
Vu = ⋅ Qe
Nc
Considerando una evolución isoterma del aire en el calderín al igual que en los vasos de
expansión y definiendo de igual forma el coeficiente de presión Cp, siendo aquí la presión máxima la
de paro Pp y la mínima la de arranque Pa en mca:
Pp − Pa
Vu = V ⋅ Pp − Pa = d
Pp
Pa = A + Hp + hp Pa = 26 mca
Siendo:
POTENCIA DE LA BOMBA.
Pa ⋅ Qe ⋅ 9.81
η = 65 P= P = 288.4 Watios
65
100
Nota: El circuito de retorno del ACS ya se calculó en su apartado correspondiente en redes no
abiertas a la atmósfera. También es muy interesante el estudio de fuerzas internas en las tuberías
( presión de trabajo y presión nominal ), sobre todo en las tuberías plásticas donde se definen
conceptos como serie S = σ / Pt y SDR = Dn / e ( e es el espesor) SDR = 2 S + 1 . Por último, el
usuario puede ver cómo se dimensiona un equipo de apoyo en ACS según la tipología elegida; estos
últimos apartados pueden consultarse en el Libro Curso de Instalaciones Térmicas y de fluidos de
dmELECT.