Vind med en intelligent varmepumpe - Aalborg Universitet

Kildekode til Matlab
program og Excel program
C
Bilag C findes også på den vedlagte cd, hvor værdierne ligger i filen NivVærdier.m og selve
behandlingen af dem ligger i filen Niveauer.m Modellen for vandtanken med omrøring
findes under model1lag.m og her optræder alle værdierne igen i starten, da den har et
anderledes defineret tidsinterval. Tiden i den første model bliver defineret efter, hvor lang
tid vandet er i hvert lag, hvorimod den i model for et lag er 0.1 sek. Hvis tidsskridtene
bliver for store eller flowet for langsomt fungerer programmet ikke. Hvis man vil ind og
ændre i koden for at se hvad en bedre varmeisolering af huset, kh, for eksempel betyder,
så opfordres det til at man kun ændrer i værdierne i filen NivVærdier.m, gemmer og så
køre filen Niveauer.m for at se, hvad det betyder. Selve koden står ude til venstre, og alt
der kommer efter % er forklaringer til, hvad der sker på den enkelte linje kode. Til sidst i
bilaget findes dokumentationen for Excelkoden
Alle de værdier der optræder i programmet
%værdier f o r s e l v e tanken :
kT = 0 . 1 4 ; %Varmekonduktivitet f o r beholder , a l t s å tab ud a f
%tanken [W/(m*C) ]
dTB =25; %F o r s k e l i temp . mellem ø v e r s t e og n e d e r s t e l a g
%i beholder [C]
TBbund = 4 0 ; %bundtemperaturen [C]
lT = 0 . 0 5 ; %Tykkelse a f beholderens væg [m]
dl = 1 ; %Antal grader mellem hvert l a g i varmebeholder [C]
dTBl = dl /dTB; %Forholdet mellem grad l a g t y k k e l s e og temperatur−
%f o r s k e l l e n mellem top og bund hver d e l a f l a g e t udgør
h = 1 . 3 ; %Højde a f varmebeholder [m]
mB = 5 00; %Masse a f vand i beholder [ kg ]
AD = (mB/( pi *1000*h ))ˆ(1/2)/0.5+ lT ; %diameter a f bunden [m]
AT = AD* pi *h ; %O v e r f l a d e a r e a l a f beholderen [mˆ 2 ]
%værdier f o r r ø r
kV = 401; %Varmekonduktivitet f o r varmeudvekslingen mellem
%vand og r ø r [W/(m*C) ]