Lab Problem: 2
Design of Shell & Tube Heat Exchanger
%Design of Shell & Tube Heat Exchanger (For No Phase Change)
%Flow rate in kg/sec
m=input('Enter the value of flow rate')
% Pressure in Kpa
PdT=80;
PdS=80;
%Unit of specific heat of oil is kj/kg c
Cp=2.14;
%Inlet temperature of lube oil in celcius
T1=65;
%Outlet temperature of lube oil in celcius
T2=45;
%Total heat duty
Q=m*Cp*(T1-T2)
%Unit of heat duty in KW
%Take the heat duty balance and to get the flow rate of water
%Specifc heat of water in kj/kg*celcius
Cpw=4.18;
%Inlet temperature of water in celcius
t1=32;
%Outlet temperature of water in celcius
t2=36;
%Flow rate in kg/sec
Mw=Q/(Cpw*(t2-t1))
%Calculate the LMTD
%Correction factor is given
Ft=0.97;
Tlm= 19.097;
%Tlm=((T1-t1)-(T2-t2))/(2.303*log((T1-t1)/(T2-t2)))
Tm=Ft*Tlm;
%Overall heat transfer coefficient is given in W/m^2*c
U=400;
%Area provided in m^2
Ap=(Q*1000)/(U*Tm);
%Data is given for calculation of the no of tubes
%Outside diameter of tube in mm
do=19.05;
%Tube length in m
L=3;
disp('Number of tubes')
Nt=Ap*1000/(3.14*do*L)
fprintf(' Number of tubes Nt = %d\n',round(Nt))
%Calculate the bundle diameter in m
k1=0.175;
n1=2.285;
db= (do*((Nt/k1)^(1/n1)))/1000
%Shell ID is given in m
Ds=db + 0.0125
% Calculate the tube side flow area in m^2
%Number of passes are given
Np=4;
di=0.0157;
At=(Nt*3.14*(di)^2)/(Np*4)
%Tube side mass velocity in kg/m^2*s
Gt=Mw/At;
%Tube side linear velocity in m/s
%Density of water is given in kg/m^3
Density=1000;
disp('Tube side velocity')
Ut=Gt/Density
%Viscocity in cP is given
Viscosity1=0.67;
Re=(di*Gt*1000)/Viscosity1
fprintf('Tube side Reynolds number = %d\n',round(Re))
%Thermal conductivity of water in W/m*c
K=0.628;
Pr=(Cpw*Viscosity1)/K
%Calculation of tube side HTC
hi=0.023*((Re)^0.8)*((Pr)^0.33) *(K/di)
disp('Tube side heat transfer coefficient')
%hi=Nu*K*1000/di
%Shell side calculations
%Shell side flow area in m^2
Pt=1.25*do;
%Baffle spacing is given in m
Bs=0.16;
As=((Pt-do)/Pt)*Ds*Bs
%Calculate the equivalent diameter of shell side
de=((1.1/do)*(((Pt)^2)-0.907*(do)^2))/1000
%Shell side mass velocity
Gs=m/As
%Viscosity of oil in cP
Viscosity2=15;
Res=(Gs*de*1000)/Viscosity2
fprintf('Shell side Reynolds number = %d\n',round(Res))
%Thermal conductivity of oil in W/m*c is given
k=0.13;
pr=Cp*Viscosity2/k
%Shell side HTC in W/m^2*c
disp('Shell side heat transfer coefficient')
ho=(0.025*(Res)*(pr)^0.33)*(k/de)
%given in kcal/hrmc
Hid=5000;
Hod =3000;
Kw=45;
Do=19.05;
Di=15.57;
Uoo=((1/ho)+ (1/Hod) + (0.001*Do*2.303*log(Do/Di)/(2*Kw)) +
((Do/Di)*(1/hi)) + ((Do/Di)*(1/Hid)))
Uo = 1/Uoo
%Uo in W/m2c
%Required area
Ar=Q*1000/(Uo*Tlm)
if(Ap<Ar)
disp('Required condition not satisfied')
else
if(Ap>Ar)
disp('Condition is satisfied')
end
end
Excessarea =((Ap/Ar)-1)*100
MATLAB Command Window
Enter the value of flow rate
9
m = 9
Q = 385.2000
Mw = 23.0383
Number of tubes
Nt = 289.6967
Number of tubes Nt = 290
db = 0.4882
Ds = 0.5007
At = 0.0140
Tube side velocity
Ut = 1.6440
Re = 3.8523e+04
Tube side Reynolds number = 38523
Pr = 4.4596
hi = 7.0248e+03
Tube side heat transfer coefficient
As = 0.0160
de = 0.0137
Gs = 561.7028
Res = 514.3701
Shell side Reynolds number = 514
pr = 246.9231
Shell side heat transfer coefficient
ho = 749.6250
Uo = 457.76
Ar = 44.0637
Condition is satisfied
Excessarea = 17.9801