Page 371 - Defrosting for Air Source Heat Pump

Appendices 367

x(2)-1.4748.*x(3).^(1/3)./1005./1.258./0.845^(2/3).*(denspipe-
dens_air)*2.6852*0.8*2.5*((mr0-smvaw)./mr0).^1.5*7.8;
% the mass of water vaporized into the ambient air is equal with the
x(4)-0.0679*5.5*(kTri-2.0*x(3))./(kRr+3.6316e-06); % 0.755 is the
modification value of refrigerant side
x(4)-1.25*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

functionF=mystage43(x,kTw1,mr0,smvaw,i,denspipe,dens_air,kTri,
kRr,kMr,khri)
% solve mrw=x(1), mvaw=x(2), Tw=x(3), qr=x(4), Rr=x(5)
F=[x(1)-x(2); % the water left on the coil is the water vaporized
from this coil; mass conservation law
x(4)-4195.2*((mr0-smvaw-x(2))*x(3)-(mr0-smvaw)*kTw1)-1480.54.*
((i+1)^0.5-i^0.5).*2.6852*2.5*0.3*x(3)-1.4748*x(3)^(4/3)
*2.6852*2.5*0.7-2443*1000*x(2);
% the energy comes from refrigerant was used in improve the
temperature
% of retained water, and heat the ambient air, and retained water
% vaporized; energy conservation, 2.0
x(2)-1.4748.*x(3).^(1/3)./1005./1.258./0.845^(2/3).*(denspipe-
dens_air)*2.6852*0.7*2.5*((mr0-smvaw)./mr0).^1.5*7.8;
% the mass of water vaporized into the ambient air is equal with the
x(4)-0.0679*4.0*(kTri-2.0*x(3))./(kRr+3.6316e-06); % 0.755 is the
modification value of refrigerant side
x(4)-1.15*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