Question

using the exergy and anergy analysis of pulse tube cryocooler inline configuration using R407 c refrigerant producing -20K temperature under atmospheric pressure conditions....simulate using fortran 77 and discuss its effect on overall efficiency.....

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

REAL:: E=2.718281828, pie=3.141592654, dia=1.63E-3, w=0.010
c   print*, "enter the value of temperature"
   t1=40
   z=2418.4/(t1+273.15)
   y=15.06-z
   p1=1000*E**y
   print*, "Pressure=", p1/1000
   Vf1=(0.777+0.002062*t1+0.00001608*t1**2)/1000
   V1=Vf1
   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-(t1*1.715E-6)+(t1*t1*8.869E-9)
   print*, "viscosity of saturated liquid=", uf1
   d= (w*4)/(pie*dia**2)
   Vel1=d*Vf1
   print*, "velocity=", Vel1
   Re1=(Vel1*dia)/(Vf1*uf1)
   print*, "Reynolds Number=", Re1
   f1=0.33/(Re1)**0.25
   print*, "Friction Factor=", f1

   t=39
   z=2418.4/(t+273.15)
   y=15.06-z
   p=1000*E**y
   print*, "Pressure=", p/1000
   Vf=(0.777+0.002062*t+0.00001608*t*t)/1000
   print*, "Vf=", Vf
   Vg= (-4.26+94050*(t+273.15)/p)/1000
   print*, "Vg=", Vg
   hf= 200+1.172*t+0.001854*t**2
   print*, "hf=", hf
   hg=405.5+0.3636*t-0.002273*t**2
   print*, "hg=", hg
   uf=0.0002367-(t*1.715E-6)+(t*t*8.869E-9)
   print*, "uf=", uf
   ug= (11.945E-6)+(t*50.06E-9)+(t*t*0.2560E-9)
   print*, "ug=", ug
   d= (w*4)/(pie*dia**2)
   a=0.5*(Vg-Vf)*(Vg-Vf)*d**2
   print*, "a=", a
   b=1000*(hg-hf)+Vf*(Vg-Vf)*d**2
   print*, "b=", b
   c=1000*(hf-hf1)+(0.5*d*d*Vf**2)-(0.5*Vel1**2)
   print*, "c=", c
   x=(-b+ SQRT((b*b)-(4*a*c)))/(2*a)
   print*, "dryness fraction=", x
   h=hf*(1-x)+x*hg
   print*, "enthalpy=", h
   v=vf*(1-x)+vg*x
   print*, "specific volume=", v
   u=uf*(1-x)+x*ug
   print*, "viscosity=", u
   Vel=d*v
   print*, "velocity=", Vel
   Re=(Vel*d)/(V*u)
   print*, "Reynolds Number=", Re
   f=0.33/(Re)**0.25
   print*, "friction factor=", f
  
  
   stop
   end