Question

I need a mathematical model governing the mechanism of air source heat pump with split evaporator connected with bryair dehumidifier for water heating and space cooling application....draw a schematic arrangement using AutoCAD and develop a mathematical model in fortran 77...

please don't write irrelevant answers...if u don't know it's ok skip but don't write wrong answers

Answer 1

C PROGRAM FOR CALCULATION OF LENGTH OF CAPILLARY TUBE(R-22)
   IMPLICIT NONE
   integer::t1,t2,Re,Re1,Re2
   Real::p1,p2,Vf1,Vf2,Vg1,Vg2,V1,V2,hf1,hf2,hg1,hg2,h1,
*h2,uf1,uf2,u1,u2,ug1,ug2,Area,Vel1,Vel2,x1,x2,x,d,a,b,c,fm,f1,f2,
*inc_length1,inc_length2,cum_length1,cum_length2,Vm,
*e=2.718281828,pie=3.141592654,dia=1.63E-3,w=0.010,tc=40,te=5,
* t,p,h,V,u,inc_length, cum_length
   Area=(pie*dia**2)/4
   print*, "Area=", Area
   d=(4*w)/(pie*dia**2)
   print*, "d=", d

   t1=40
   p1=1000*E**(15.06-(2418.4/(t1+273.15)))
   p=p1/1000
   print*, "Pressure=", p1/1000
   Vf1=(0.777+0.002062*t1+0.00001608*t1**2)/1000
   PRINT*, "specific volume of saturated liquid=", Vf1
   hf1= 200+1.172*t1+0.001854*t1**2
   print*, "enthalpy of saturated liquid=", hf1
   uf1=0.0002367-(1.715E-6*t1)+(8.869E-9*t1**2)
   print*, "viscosity of saturated liquid=", uf1
   V1=Vf1
   h1=hf1
   u1=uf1
   Vel1=d*V1
   print*, "velocity=", Vel1
   Re1=(Vel1*dia)/(V1*u1)
   print*, "Reynolds Number=", Re1
   f1=0.33/(Re1)**0.25
   print*, "Friction Factor=", f1

   inc_length1=0
   cum_length1=0
  
      
   do 10 t2=39,5,-1
   print*, "t2=",t2  
   p2=1000*E**(15.06-(2418.4/(t2+273.15)))
   p=p2/1000
C   print*, "Pressure=", p2/1000
   Vf2=(0.777+0.002062*t2+0.00001608*t2**2)/1000
c   PRINT*, "specific volume of saturated liquid=", Vf2
   Vg2=(-4.26+(94050*(t2+273.15))/p2)/1000
c   print*, "Specific Volume of saturated vapour=", Vg2
   hf2= 200+(1.172*t2)+(0.001854*t2**2)
c   print*, "enthalpy of saturated liquid=", hf2
   hg2=405.5+(0.3636*t2)-(0.002273*t2**2)
c   print*, "enthalpy of saturated vapour=", hg2
   uf2=(0.0002367)-(1.715E-6*t2)+(8.869E-9*t2**2)
c   print*, "viscosity of saturated liquid=", uf2
   ug2=(11.945E-6)+(50.06E-9*t2)+(0.2560E-9*t2**2)
c   print*, "viscosity of saturated vapour=", ug2
c
   a=0.5*(Vg2-Vf2)*(Vg2-Vf2)*d**2
c   print*, "a=", a
   b=1000.0*(hg2-hf2)+(Vf2*(Vg2-Vf2)*d**2)
c   print*, "b=", b
   c=1000.0*(hf2-h1)+((0.5*d**2)*(Vf2**2))-(0.5*Vel1**2)
c   print*, "c=", c
   x1=(-b+ SQRT(b**2-4*a*c))/(2*a)
   x2=(-b- SQRT(b**2-4*a*c))/(2*a)
   x=max(x1,x2)
   print*, "dryness fraction=", x

   h2=hf2*(1-x)+x*hg2
   print*, "enthalpy=", h2
   V2=Vf2*(1-x)+x*Vg2
   print*, "specific volume=", V2
   u2=uf2*(1-x)+x*ug2
   print*, "viscosity=", u2
   Vel2=d*V2
   print*, "velocity=", Vel2
   Re2=(Vel2*dia)/(V2*u2)
   print*, "Reynolds Number=", Re2
   f2=0.33/(Re2)**0.25
   print*, "Friction Factor=", f2

   fm=(f1+f2)/2
   print*, "fm=", fm
   Vm=0.5*(Vel1+Vel2)
   print*, "Vm=", Vm
   inc_length2=(2*Area*dia*((p1-p2)-(w*(V2-V1))))/(fm*Vm*w)
   print*, "inc_length=", inc_length2
   cum_length2=cum_length1+inc_length2
   print*, "cum_length=", cum_length2

   p1=p2
   Vf1=Vf2
   Vg1=Vg2
   V1=V2
   hf1=hf2
   hg1=hg2
   h1=h2
   uf1=uf2
   ug1=ug2
   u1=u2
   Vel1=Vel2
   Re1=Re2
   f1=f2
   inc_length1=inc_length2
   cum_length1=cum_length2
   OPEN (UNIT=1, FILE='RESULTS1', STATUS='OLD')
   OPEN (UNIT=2, FILE='RESULTS2', STATUS='OLD')
   WRITE (2,*) t1,p,x,V1,h1,Vel1,inc_length1, cum_length1
   CLOSE(2)
   WRITE (1,*) t2,p,x,V2,h2,Vel2,inc_length2, cum_length2
10   continue
   close(1)
  
  
   stop
   end