Page 369 - Defrosting for Air Source Heat Pump

Appendices 365

x(4)-0.0679*(kTri-2.*x(3))./(kRr+3.6316e-06); % 0.755 is the mod-
ification value of refrigerant side
x(4)-0.32*kMr*(khri-44518-1170.36*x(5)-1.68674*x(5)^2-5.2703/
1000*x(5)^3)] % 0.40 is to modify the Mr. -5.2703/1000*kTro^3, W
end

function F=mystage31(x,kmw1,ksmrw,kTw1,i,kTri,kRr,kMr,khri)
% solve mf=x(1), mrw=x(2), Tw=x(3), qr=x(4), Rr=x(5)
F=[x(1)-x(2)-kmw1; % mass conservation, frost changes to water, kg/s
x(4)-334000.*x(1)-4195.2.*ksmrw*(x(3)-kTw1)-1.4748.*x(3).^(4/3).
*2.6852*2.5*0.55; % the energy in refrigerant used in frost melt-
ing, retained water temperature improve, heating the air, W
% 2.6852*2.5*4.0*0.75 is the actual area of contact fins surface
and the ambient air
1.4748.*x(3).^(4/3).*2.6852*2.5*0.55+334000.*x(1)-1480.54.*((i
+1)^0.5-i^0.5).*2.6852*2.5*x(3)*0.45; % the heat transfer from
water layer; 0.85 stand the area is not 2.6852*2.5*4.0;
x(4)-0.0679*1.6*(kTri-2.*x(3))./(kRr+3.6316e-06); % 0.755 is the
modification value of refrigerant side
x(4)-0.42*kMr*(khri-44518-1170.36*x(5)-1.68674*x(5)^2-5.2703/
1000*x(5)^3)] % 0.40 is to modify the Mr.-5.2703/1000*kTro^3, W
end

function F=mystage32(x,kmw1,kmw2,ksmrw,kTw1,i,kTri,kRr,kMr,khri)
% solve mf=x(1), mrw=x(2), Tw=x(3), qr=x(4), Rr=x(5)
F=[x(1)+kmw1-x(2)-kmw2; % mass conservation, frost changes to
water, kg/s
x(4)-334000.*x(1)-4195.2.*(ksmrw+kmw1-kmw2)*(x(3)-kTw1)-1.4748.
*x(3).^(4/3).*2.6852*2.5*0.50; % the energy in refrigerant used in
frost melting, retained water temperature improve, heating the air
% 2.6852*2.5*4.0*0.75 is the actual area of contact fins surface
and the ambient air
1.4748.*x(3).^(4/3)*2.6852*2.5*0.50+334000.*x(1)-1480.54.*((i+1)
^0.5-i^0.5).*2.6852*2.5*x(3)*0.50; % the heat transfer from water
layer; 0.85 stand for the area, which is not 2.6852*2.5*4.0;
x(4)-0.0679*1.3*(kTri-2*x(3))./(kRr+3.6316e-06); % 0.755 is the
modification value of refrigerant side
x(4)-0.375*kMr*(khri-44518-1170.36*x(5)-1.68674*x(5)^2-5.2703/
1000*x(5)^3)] % 0.30 is to modify the Mr.-5.2703/1000*kTro^3, W
end

function F=mystage33(x,kmw2,kmw3,ksmrw,kTw1,i,kTri,kRr,kMr,khri)
% solve mf=x(1), mrw=x(2), Tw=x(3), qr=x(4), Rr=x(5)
F=[x(1)+kmw2-x(2)-kmw3; % mass conservation, frost changes to
water, kg/s
x(4)-334000.*x(1)-4195.2.*(ksmrw+kmw2-kmw3).*(x(3)-kTw1)-1.4748.
*x(3).^(4/3).*2.6852*2.5*0.45; % the energy in refrigerant used in
frost melting, retained water temperature improve, heating the air