Page 370 - Defrosting for Air Source Heat Pump
P. 370

366                                                          Appendices

            1.4748.*x(3).^(4/3).*2.6852*2.5*0.45+334000.*x(1)-1480.54.*((i
            +1)^0.5-i^0.5).*2.6852*2.5*0.55*x(3);  %  the  heat  transfer  from
            water layer
               % 2.6852*2.5*4.0*0.75 is the actual area of contact fins surface
            and the ambient air;
            x(4)-0.0679*1.0*(kTri-2*x(3))./(kRr+3.6316e-06); % 0.755 is the
            modification value of refrigerant side
            x(4)-0.33*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
            functionF=mystage41(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.1*x(3)-1.4748*x(3)^(4/3)
            *2.6852*2.5*0.9-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*2.5*0.9*((mr0-smvaw)./mr0).^1.5*7.8;
            % the mass of water vaporized into the ambient air is equal with
            the
            x(4)-0.0679*7.0*(kTri-2.0*x(3))./(kRr+3.6316e-06); % 0.755 is the
            modification value of refrigerant side
            x(4)-1.35*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=mystage42(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.2*x(3)-1.4748*x(3)^(4/3)
            *2.6852*2.5*0.8-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
   365   366   367   368   369   370   371   372   373   374   375