Page 365 - Defrosting for Air Source Heat Pump

Appendices 361

A
x00=real(A);
fval
exit
2
hair(j,i)=1.4748.*Tri(j,1).^(1/3); % W/(K˙m )
qair(j,i)=1.4748.*Tri(j,1).^(4/3).*2.6852*2.5*2; % W
s_qair(j,i)=sum(qair(:,i))*5; % W
2
hd(j,i)=hair(j,i)/1005./1.258./0.845^(2/3); % W/(K˙m )
qvap(j,i)=mvaw(j,i)*2443*1000; % W
smvaw(j,i)=5.*sum(mvaw(:,i)); % kg
qm(j,i)=334000.*mf(j,i); % W
qvap(j,i)=mvaw(j,i)*2443*1000; % W
s_qvap(j,i)=sum(qvap(:,i))*5; % W
watertray(j,i)=0; % kg/s
swatertray(j,i)=sum(watertray(:,i)); % kg
hro(j,i)=44518+1170.36*Tro(j,i)+1.68674*Tro(j,i)^2+5.2703/
1000*Tro(j,i)^3;
2
% W/(K˙m )
qr2(j,i)=kMr*(khri-hro(j,i)); % W
s_qr2(j,i)=sum(qr2(:,i))*5; % W
end
end
end
% here is the end of stage 4 for Circuit 2: water layer evaporating
stage
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

for i=3
for j=18:45
% for the 18*5 seconds for the 1st circuit;
kmw2=mf(j-1,i-2)+mf(j-1,i-1)+0.0192*2.0; % kg/s
kmw3=kmw2+mf(j-1,i)+0.0192; % kg/s
ksmrw=sfrost(17,i); % kg
kTw1=Tw(j-1,i); % °C
kTri=Tri(j,i); % °C
2
kRr=Rr(j,i); % (K˙m )/W
kMr=Mr(j,i); % kg/s
khri=hri(j,i); % kJ/kg
% all the input parameters in the function listed here
x0=[0.0042 0.0042 0.335 1200 0.001];%mf=x(1), mr=x(2),
Tw=x(3); qr=x(4); Tro=x(5) the values of debugging
options=optimset(’display’,’off’,’MaxIter’,100000,
’MaxFunEvals’,20000); % number
[A,fval,exit]=fsolve(@(x) mystage33(x,kmw2,kmw3,ksmrw,
kTw1,i,kTri,kRr,kMr,khri),x0,options);