Cieplne właściwości dynamiczne grzejnika podłogowego

Cieplne właściwości dynamiczne grzejnika podłogowego
Zatem dla prostego układu zaprezentowanego na rysunku 6.5. bilans ciepła w węźle
o współrzędnych i,j będzie miał postać:
gdzie:
λ
⎛
⎞ ⎛
⎞
⎜t'
−t'
⎟ + λ ⎜t'
−t'
⎟ +
i −1;
j ⎝ i −1;
j i;
j ⎠ i;
j −1⎝
i;
j −1
i;
j ⎠
+ λ
⎛
⎞ ⎛
⎜t'
−t'
⎟ + λ ⎜t'
−t'
i + 1;
j ⎝ i + 1;
j i;
j ⎠ i;
j + 1⎝
i;
j + 1 i;
j
= V ⋅ ρ ⋅c
⋅
⎛
⎜t'
−t''
p ⎝ i;
j i;
j
λ – współczynnik przewodzenia ciepła pomiędzy sąsiadującymi węzłami, [W/mK],
t’ – temperatura w n-tym kroku czasowym [K],
t’’ – temperatura w n-1 kroku czasowym, [K].
65
⎞
⎟
⎠
⎞
⎟
⎠
=
(6.24)
Opracowano autorski program obliczeniowy do rozwiązania nieustalonej wymiany
ciepła w wodnych grzejnikach podłogowych. Program ten dotyczy obliczania nieustalonej
wymiany ciepła za pomocą metody bilansów elementarnych.
Program w sposób automatyczny generuje siatkę podziałów elementarnych, badając
kryteria stabilności rozwiązań. Gęstość podziału siatki jest kompromisem zawartym
pomiędzy stabilnością obliczeń a ich szybkością. Zaawansowane funkcje edycji parametrów
grzejnika podłogowego pozwalają na dokładny opis symulowanego układu.
Przykładowe obliczenia, wraz z wizualizacją wyników obliczeń wykonanych za
pomocą autorskiego programu przedstawiono na rysunku 6.6.