disp('*******************DESIGN OF DPHX*********************')
disp('*******************DESIGN DPHX*********************') ;
fprintf ('\n');
jch=input('1.Counter flow 2. Parallel Flow ');
disp('**Tube side parameter entry**');
fprintf ('\n');
hit= input ('Choose the tube side fluid 1. Water 2. Ethelene glycol 3. Engine oi
l');
fprintf ('\n');
mt=input('Mass flow rate of fluid ');
thi=input('Input fluid temperature ');
tho=input('Exit fluid temperature ');
fprintf ('\n');
ch=input('Type of pipe :- 1. Plain pipe
pipe 2. Outer fin');

fprintf
if ch==1('\n');
  di=input('Inner diameter of the pipe ');
  do=input('Outer diameter of the pipe ');
  Di= input ('Inner diameter of the shell');
  L=input('Length of hair pin');
  Nt=input('No: of tubes');
  kp=input('Thermal conductivity of the pipe material');

else
  di=input('Inner diameter of the pipe ');
  do=input('Outer diameter of the pipe ');
  Di= input ('Inner diameter of the shell');
  L=input('Length of heat pipe');
  hf=input('Fin height ');
  t=input('Thickness of fin ');
  Nt=input('No: of tubes ');
  Nf=input('No: of fins ');
  kf= input ('Thermal conductivity of the fin');
  kp=input('Thermal conductivity of the pipe material');
end
tt=(thi+tho)/2;
fprintf ('\n');
disp ('Tube side fluid properties at ');
disp (tt);
if hit==1
den=-0.0027*tt^2 - 0.1468*tt+ 1003;
pr=7*10^-13*tt^6 - 7*10^-10*tt^5 + 3*10^-07*tt^4 - 6*10^-05*tt^3 + 0.0071*tt^2 -
 0.4335*tt + 13.528;
kv=(7*10^-14*tt^6 - 8*10^-11*tt^5 + 3*10^-08*tt^4 - 7*10^-06*tt^3 + 0.0008*tt^2
- 0.051*tt + 1.7807)*10^-6;
k=-1*10^-10*tt^4 + 6*10^-08*tt^3 - 2*10^-05*tt^2 + 0.0026*tt + 0.552;
cph=8*10^-7*tt^4 - 0.0004*tt^3 + 0.068*tt^2 - 3.8529*tt + 4227.6;
end
if hit==2
den= (-2*10^(-6)*tt^4) + (0.0003*tt^3) - (0.0161*tt^2) - (0.5202*tt) + 1130.9 ;
pr= (3*10^-5*tt^4) - (0.0075*tt^3) + (0.7286*tt^2) - (32.173*tt) + 614.39 ;
kv= (3*10^(-12)*tt^4) - (7*10^(-10)*tt^3) + (7*10^(-8)*tt^2) - (3*10^(-6)*tt) +
(6*10^(-5));
k= (5*10^(-10)*tt^4) - (9*10^-8*tt^3) + (3*10^-6*tt^2) + (0.0003*tt) + 0.2419 ;
cph = (4.4743*tt) + 2293.6 ;
end
if hit==3
den= (-0.5933*tt) + 899.47 ;
pr= (10^(-7)*tt^6) - (6*10^(-5)*tt^5) + (0.013*tt^4) - (1.5329*tt^3) + (97.718*t
t^2) - (3266.5*tt) + 47085 ;
kv= (10^(-14)*tt^6) - (5*10^-12*tt^5) + (10^(-9)*tt^4) - (10^-7*tt^3) + (9*10^(-
 

6)*tt^2) - (0.0003*tt) + 0.0043 ;

k= (3*10^(-14)*tt^6) - (10^(-11)*tt^5) + (3*10^-9*tt^4) - (2*10^(-7)*tt^3) + (8*
10^(-6)*tt^2) - (0.0002*tt) + 0.1477 ;
cph = (4.2925*tt) + 1792.4 ;
end
dv=den*kv;
disp ('Specific heat ');
disp (cph);
disp ('Density ');
disp (den);
disp ('Prandtl Number ');
disp (pr);
disp ('Kinematic viscosity ');

disp
disp (kv);
('Thermal conductivity ');
disp (k);
Rfi= input('Fouling factor ');
jj=den;
fprintf ('\n');
if thi>tho
  Q=mt*cph*(thi-tho);
  disp('TUBE Side fluid is hot');
  j=0.3;
else
  Q=mt*cph*(tho-thi);
  disp('TUBE Side fluid is cold');
  j=0.4;
end
Act=3.14*di^2/4;
Umt=mt/(den*Act*Nt);
Re=(den*Umt*di)/dv;
disp ('********************************************');
('********************************************');
disp('Mean velocity ');
disp(Umt);
disp('Reynolds Number ');
disp(Re);
disp('Heat Transfer rate');
disp(Q);
fprintf ('\n');
if Re>2300
  Nu= 0.023*Re^0.8*pr^j;
  ft=0.184*Re^-0.2;
else
  x1=0.04*Re*di;
  x2=0.04*Re*di*pr;
  ft=64/Re;
  disp ('Hydrodynamic entry length');
  disp (x1)
  disp ('Thermal entry length');
  disp (x2);
  if x1<L && x2<L
  disp ('Hydrodynamicaly and thermaly fully developed');
  else
  if x1<L && x2>L
disp('Hydrodynamicaly developed and thermal developing');
  else
  if x1>L && x2<L
  disp ('Hydrodynamicaly under developed and thermaly developed');
  else
  disp('Hydrodynamicaly and thermaly under developed');
 

  end
  end
  end
Nu=3.66 +((0.0668*(di/L)*Re*pr)/(1+(0.04*(di/L)*Re*pr)
+((0.0668*(di/L)*Re*pr)/(1+(0.04*(di/L)*Re*pr)^0.67));
^0.67));
end
disp ('***********************************************
('***********************************************');
');
fprintf ('\n');
disp ('Nusslet Number');
disp (Nu);
hi= (Nu*k)/di;
disp ('Tube side Heat transfer coefficient');
disp (hi);
fprintf ('\n');

fprintf ('\n');
disp('**Shell side parameter entry**');
fprintf ('\n');
hot= input ('Choose the shell side fluid 1. Water 2. Ethelene glycol 3. Engine o
il');
fprintf ('\n');
ms=input('Mass flow rate of fluid ');
tci=input('Input fluid temperature ');
tco=input('Exit fluid temperature ');
if ch==1
  Acs= 3.14*(Di^2-di^2)/4;
  At= 3.14*do*L*Nt*2;
  dh= do-di;
  de= (4*Acs)/(3.14*do);
else
  Acs= 0.0012632 ;
  At=7.612 ;
  Au=0.51107;
  Af=7.101;
  Pw=(3.14*(Di+(do*Nt)))+(2*hf*Nf*Nt);
  Ph=(3.14*do*Nt)+(2*hf*Nf*Nt);
  dh=0.005 ;
  de=0.005974;
end
Ai= 0.5915 ;
ts=(tci+tco)/2;
fprintf ('\n');
disp ('Shell side fluid properties at ');
disp (ts);
if hot==1
den=-0.0027*ts^2 - 0.1468*ts+ 1003;
pr=7*10^-13*ts^6 - 7*10^-10*ts^5 + 3*10^-07*ts^4 - 6*10^-05*ts^3 + 0.0071*ts^2 -
 0.4335*ts + 13.528;
kv=(7*10^-14*ts^6 - 8*10^-11*ts^5 + 3*10^-08*ts^4 - 7*10^-06*ts^3 + 0.0008*ts^2
- 0.051*ts + 1.7807)*10^-6;
k=-1*10^-10*ts^4 + 6*10^-08*ts^3 - 2*10^-05*ts^2 + 0.0026*ts + 0.552;
end
if hot==2
den= (-2*10^(-6)*ts^4) + (0.0003*ts^3) - (0.0161*ts^2) - (0.5202*ts) + 1130.9 ;
pr= (3*10^-5*ts^4) - (0.0075*ts^3) + (0.7286*ts^2) - (32.173*ts) + 614.39 ;
kv= (3*10^(-12)*ts^4) - (7*10^(-10)*ts^3) + (7*10^(-8)*ts^2) - (3*10^(-6)*ts) +
(6*10^(-5));
k= (5*10^(-10)*ts^4) - (9*10^-8*ts^3) + (3*10^-6*ts^2) + (0.0003*ts) + 0.2419 ;
cph = (4.4743*ts) + 2293.6 ;
end
if hot==3
den= (-0.5933*ts) + 899.47 ;
 

pr= 1050;
kv= 0.001*exp(-0.041*ts) ;
k= (3*10^(-14)*ts^6) - (10^(-11)*ts^5) + (3*10^-9*ts^4) - (2*10^(-7)*ts^3) + (8*
10^(-6)*ts^2) - (0.0002*ts) + 0.1381 ;
cph = (4.2925*ts) + 1792.4 ;
end
dv=den*kv;
fprintf ('\n');
disp ('Specific heat ');
disp (cph);
disp ('Density ');
disp (den);
disp ('Prandtl Number ');

disp
disp (pr);
('Kinematic viscosity ');
disp (kv);
disp ('Thermal conductivity ');
disp (k);
fprintf ('\n');
Rfo= input('Fouling factor ');
fprintf ('\n');
if tci>tco
  Q=ms*cph*(tci-tco);
  disp('SHELL Side fluid is hot');
  j=0.3;
else
  Q=ms*cph*(tco-tci);
  disp('SHELL Side fluid is cold');
  j=0.4;
end
Ums=ms/(den*Acs);
Re=(den*Ums*dh)/dv;
disp ('**************************************');
('**************************************');
disp('Mean velocity ');
disp(Ums);
disp('Reynolds Number ');
disp(Re);
disp('Heat Transfer rate');
disp(Q);
fprintf ('\n');
if Re>2300
  Nu= 0.023*Re^0.8*pr^j;
  fs=0.184*Re^-0.2;
else
  x1=0.04*Re*di;
  x2=0.04*Re*di*pr;
  fs=64/Re;
  disp ('Hydrodynamic entry length');
  disp (x1)
  disp ('Thermal entry length');
  disp (x2);
  if x1<L && x2<L
  disp ('Hydrodynamicaly and thermaly fully developed');
  else
  if x1<L && x2>L
disp('Hydrodynamicaly developed and thermal developing');
  else
  if x1>L && x2<L
  disp ('Hydrodynamicaly under developed and thermaly developed');
  else
 

  disp('Hydrodynamicaly and thermaly under developed');

  end
  end
  end
Nu=3.66 +(0.0668*(de/L)*Re*pr)/(1+(0.04*((de/L)*Re*pr)
+(0.0668*(de/L)*Re*pr)/(1+(0.04*((de/L)*Re*pr)^0.67));
^0.67));
end
disp ('******************************************');
('******************************************');
fprintf ('\n');
disp ('Nusslet Number');
disp (Nu);
ho= (Nu*k)/de;
disp ('Shell side Heat transfer coefficient');
disp (ho);

Rw= log(do/di)/(2*3.14*kp*
log(do/di)/(2*3.14*kp*L);
if ch==1 L);
  Uo= 1/((1/ho)+Rfo+((do/2)*
1/((1/ho)+Rfo+((do/2)*(log(do/di))/kp)+(do/(di
(log(do/di))/kp)+(do/(di*hi)));
*hi)));
else
  m=((2*ho)/(t*kf))^0.5;
  Yf= 0.6835 ;
  Yo= 0.7047;
  Uo=1/((At/(Ai*hi))+((At*Rfi)/Ai)+(At*Rw)+(Rfo/Yo)+(1/(Yo*ho)));
end
if ch==1
  Ui= 1/((1/hi)+Rfi+((di/2)*
1/((1/hi)+Rfi+((di/2)*(log(do/di))/kp)+(di/(do
(log(do/di))/kp)+(di/(do*ho)));
*ho)));
else
  m=((2*ho)/(t*kf))^0.5;
  Yf= 0.6835 ;
  Yo= 0.7047;
  Ui=1/((1/hi)+Rfi+(Ai*Rw)+((Ai*Rfo)/(At*Yo))+(Ai/(At*Yo*ho)));
end
disp ('*****************************');
('*****************************');
fprintf ('\n');
disp('Overall heat transfer coefficient w.r.t outer surface');
disp (Uo);
disp('Overall heat transfer coefficient w.r.t inner surface');
disp (Ui);
if jch==1
  if thi>tho
  T1=thi-tco;
  T2=tho-tci;
  else
  T1=tci-tho;
  T2=tco-thi;
  end
  if T1==T2
  LMTD=T1;
  else
  LMTD=(T1-T2)/log(T1/T2);
  end
else
  if thi>tci
  T1=thi-tci;
  T2=tho-tco;
  else
  T1=tci-thi;
  T2=tco-tho;
  end
  LMTD= (T1-T2)/log(T1/T2);
end
A= Q/(Uo*LMTD);
 

fprintf ('\n');
fprintf ('\n');
disp ('***RESULTS***');
disp ('Total heat transfer area ');
disp (A);
Njp= A/At;
disp ('No: of hair pins ');
disp (Njp);
Nhp= round (Njp+1);
Yp= input('Enter Pump efficiency ');
Ps=(fs*2*L*den*Ums^2*Nhp)*2^(-1)*dh^(-1);
Ws = (ms*Ps)/(den*Yp);
den=jj;

Pt=(ft*2*L*den*Umt^2*Nhp)*(2^(-1)*di^(-1));
Wt = (mt*Pt)/(den*Yp) ;
disp ('Shell side pressure drop ');
disp(Ps);
disp ('Tube side pressure drop ');
disp(Pt);
disp ('Shell side pumping power ');
disp(Ws);
disp ('Tube side pumping power ');
disp(Wt);