Vind med en intelligent varmepumpe - Aalborg Universitet

%a l den e n e r g i der s k a l bruges t i l at holde huset varmt om dagen
Qtot = Qnat+Qdag+c *varmtvand *25*(1− nablavand ) ;
%energimængden der s k a l bruges t i l at varme den ønskede mængde vand
L a g e r s t r = Qtot/c /(TBbund+0.5*dTB−TRind )
%Beregner l a g e r e t s t r . i n t e t semikolon f o r så b l i v e r værdien
%a f l a g e r s t r smidt ud i matlab vinduet
l a g e r e n e r g i = mB* c *(TBbund+0.5*dTB−TRind )
%Beregner mængden a f e n e r g i t i l rådighed i l a g e r e t i n t e t semikolon
%f o r så b l i v e r værdien a f l a g e r s t r smidt ud i matlab vinduet
time = 0 : dt : round ( t s l u t −dt ) ;
%hvor mange t i d s i n d e l i n g e r der er i a l t , og t i d e n s l u t t e r .
Selve programmet:
NivVaerdier %henter værdierne f r a f i l e n Nivvardier
ArrayB = z e r o s ( 1 , round (dTB/ dl ) ) ;
%Talrække over temperaturen i beholderens l a g [C]
ArrayR = z e r o s ( 1 , round ( t s l u t / dt ) ) ;
%Talrække over temperaturen a f hvert l a g i varmerøret [C]
ArrayH = z e r o s ( 1 , round ( t s l u t / dt ) ) ;
%Talrække over e n e r g i b a l a n c e n f o r hust ( f o r meget e l l e r
%l i d t e n e r g i i f h t . behovet
f o r i i = 1 : round ( t s l u t / dt ) %Temperaturen i r ø r e r n e
ArrayR ( i i )=TRind ; %s t a r t e r ved d e r e s indgangstemperatur
end
f o r j j = 0 : round (dTB/ dl −1)% temperaturen a f de f o r s k e l l i g e l a g
ArrayB ( j j +1)=TBbund+j j * dl ; %s t a r t e r a l l e ved temperaturen i
end %bunden + den f o r s k e l , der er mellem bund og l a g e t s p l a c e r i n g
ArrayBgemt = ArrayB ; %Vi gemmer talrækken t i l s e n e r e sammenligning
%da den v i s e r l a g e r e t l a g d e l i n g e r ved s t a r t
f o r i i = 1 : t s l u t / dt %i i =s t a r t : i n t e r v a l : s l u t
i f i i == round (12*3600/ dt ) %h v i s t i d e n er l i g med et halvt døgn ,
Qhustemp=Qhusdag+Qvarmtvand ; %så s k i f t t i l dag
end
ArrayH ( i i )=Qhustemp ; %s t a r t e r med at opvarmningsbehovet er f o r
%natten samt varmt vand
f o r j j = 1 :dTB/ dl
Qoverf = kV*AV*( ArrayB ( j j )−ArrayR ( i i ) ) / lV * dt ;
%beregn det der b l i v e r o v e r f ø r t t i l r ø r e t
ArrayB ( j j )=ArrayB ( j j )−Qoverf /(mBl* c ) ;
%træk det der blev o v e r f ø r t t i l r ø r e t f r a beholderen
ArrayR ( i i )=ArrayR ( i i )+Qoverf /(mRl* c ) ;
%Læg det der b l i v e r o v e r f ø r t t i l r ø r e t t i l r ø r e t s
%temperatur
Qtab = kT*ATB*( ArrayB ( j j )−Ta)/ lT * dt ;
%Beregn varmetabet f o r vandtanken
Qtab = Qtab * 1 . 3 3 ; %Pga . manglende top og bund − f o r d e l e r