%%PROBLEM_4: input data for fixed tubesheet type shell and tybe heat
exchanger
T1=65; %C inlet temp of lube oil
T2=45; %C outlet temp of lube oil
ml=9.6459; %kg/sec massflowrate of lube oil
Dl=869; %kg/m3 density of lube oil
Cpl=2.1413; %kJ/kg C specific heat of lube oil
kl=0.13; %w/m C thermal conductivity of lube oil
Vl=15; %cP viscosity of lube oil
kwater=0.628; %w/m C thermal conductivity of water
t1=35; %C inlet temp of water
t2=39; %C outlet temp of water
Cpw=4.1868; %kJ/kg C specific heat of water
Vw=0.73; %cP viscosity of water
Dw=993.32; %kg/m3 density of water
%%specification of fixed tubesheet heat exchanger
Di=418; %mm Shell inside diameter
do=15.875; %mm tube o.d.
di=13.3858; %mm tube i.d.
L=3.048; %m tube length
Np=4; %no. of tube side passes
ft=0.95; %LMTD correction
Bs=83.6; %mm baffle spacing
AUo=400; %w/m2 C assumed overall heat transfer coefficient
%% heat load & energy balance
Q=ml*Cpl*(T1-T2); %kw heat load
mw=Q/(Cpw*(t2-t1)); %kg/sec mass flowrate of water
%% mean temp difference
dT1=T1-t2;
dT2=T2-t1;
LMTD=(dT1-dT2)/log((dT1)/(dT2)); %C
MTD=LMTD*ft; %C
Apro=(Q*1000)/(AUo*MTD); %m2 Area provided
Nt=Apro/(pi*L*do/1000); %no. of tubes
%% eavaluation of tube side HTC
at=(Nt/Np)*(pi/4*(di/1000)^2); %m2 tube side flow area
Gt=mw/at; %kg/m2sec tubeside mass velocity
ut=Gt/Dw; %m/sec tube side linear velocity
Ret=(di/1000*Gt)/(Vw/1000);
Prt=(Cpw*Vw/1000*1000)/kwater;
hi=(0.023*Ret^0.8*Prt^0.33*kwater)/(di/1000); %w/m2 C tube side HTC
%% eavaluation of shell side HTC
Pt=1.25*do; %mm tube pitch
As=((Pt-do)*Di*Bs)/(Pt*1000*1000); %m2 Shell side flow area
Gs=ml/As; %kg/m2sec shellside mass velocity
us=Gs/Dl; %m/sec shell side linear velocity
de=(1.1*(Pt^2-(0.907*do^2)))/do; %mm equivalent diameter
Res=(de/1000*Gs)/(Vl/1000);
Prs=(Cpl*Vl/1000*1000)/kl;
jh=10^((0.48856*log10(Res))-2.82393);
ho=(jh*Res*Prs^(1/3)*kl)/(de/1000); %w/m2 C shell side HTC
%% overall HTC
hid=4000; %w/m2 C tube side dirt coefficient
hod=3000; %w/m2 C shell side dirt coefficient
kw=50; %w/m C thermal conductivity of wall
Uo=((ho^-1)+(hod^-1)+((do*log(do/di))/(2*kw*1000))+(do/(di*hid))+(do/
(di*hi)))^-1; %w/m2 C overall HTC
Areq=(Q*1000)/(Uo*MTD); %m2 area required
excessarea=((Apro/Areq)-1)*100; % percentage excess area
%% shell side pressure drop
AdPs=80; %kPa allowable pressure drop
Jf=0.069;
dPs=8*Jf*(Di/de)*(L*1000/Bs)*Dl*us^2*1/(2*1000); %kPa shell side pressure
drop must be less than allowable pressure drop
NBs=240; %mm baffle spacing
NAs=((Pt-do)*Di*NBs)/(Pt*1000*1000); %m2 Shell side flow area
NGs=ml/NAs; %kg/m2sec shellside mass velocity
Nus=NGs/Dl; %m/sec shell side linear velocity
Nde=(1.1*(Pt^2-(0.907*do^2)))/do; %mm equivalent diameter
NRes=(Nde/1000*NGs)/(Vl/1000);
Prs=(Cpl*Vl/1000*1000)/kl;
Njh=10^((0.48856*log10(NRes))-2.82393);
Nho=(Njh*NRes*Prs^(1/3)*kl)/(Nde/1000); %w/m2 C shell side HTC
NUo=((Nho^-1)+(hod^-1)+((do*log(do/di))/(2*kw*1000))+(do/(di*hid))+(do/
(di*hi)))^-1; %w/m2 C overall HTC
NAreq=(Q*1000)/(NUo*MTD); %m2 area required
Nexcessarea=((Apro/NAreq)-1)*100; % percentage excess area
NJf=0.065;
NdPs=8*NJf*(Di/Nde)*(L*1000/NBs)*Dl*Nus^2*1/(2*1000); %kPa shell side
pressure drop must be less than allowable pressure drop
%% tube side pressure drop
AdPt=80; %kPa allowable pressure drop
Jf=0.0037;
dPt=4*((8*Jf*(L*1000/di)*1)+2.5)*Dw*ut^2/(1000*2); %kPa shell side pressure
drop must be less than allowable pressure drop
%% Result
disp('Design Calculation of Shell & Tube Exchanger with No Phase Change')
fprintf('Heat Duty = %4.2f kW\n',Q)
fprintf('Mass Flowrate of Water = %8.3f kg/sec\n',mw)
fprintf('Mean Temperature Difference = %4.2f Degree Celsius\n',MTD)
fprintf('Area Provided = %8.3f m^2\n',Apro)
fprintf('number of tubes = %d\n',round(Nt))
fprintf('Tube side flow area = %8.3f m^2\n',at)
fprintf('Tube side mass velocity = %4.2f kg/m2sec \n',Gt)
fprintf('Tube side linear velocity = %4.2f m/sec \n',ut)
fprintf('Tube side Reynolds number = %d\n',round(Ret))
fprintf('Tube side Prandtl number = %4.2f \n',Prt)
fprintf('Tube side Heat Transfer coefficient = %4.2f W/m^2 C \n',hi)
fprintf('Shell side flow area = %8.3f m^2\n',As)
fprintf('Shell side mass velocity = %4.2f kg/m2sec \n',Gs)
fprintf('Shell side linear velocity = %4.2f m/sec \n',us)
fprintf('Shell side Reynolds number = %d\n',round(Res))
fprintf('Shell side Prandtl number = %4.2f \n',Prs)
fprintf('Shell side Heat Transfer coefficient = %4.2f W/m^2 C \n',ho)
fprintf('Overall Heat Transfer coefficient = %4.2f W/m^2 C \n',Uo)
fprintf('Area Required = %8.3f m^2\n',Areq)
fprintf('percentage excess area = %d % \n',round(excessarea))
disp('percentage heat transfer area is too high but we will first check
shell side pressure drop')
fprintf('Allowable Shell side pressure drop = %4.2f kPa \n',AdPs)
fprintf('Shell side pressure drop = %4.2f kPa \n',dPs)
disp('Shell side pressure drop is too high than allowable shell side
pressure drop')
disp('to decrease it we have to increase Baffle spacing to 200 mm')
fprintf('Shell side flow area = %8.3f m^2\n',NAs)
fprintf('Shell side mass velocity = %4.2f kg/m2sec \n',NGs)
fprintf('Shell side linear velocity = %4.2f m/sec \n',Nus)
fprintf('Shell side Reynolds number = %d\n',round(NRes))
fprintf('Shell side Prandtl number = %4.2f \n',Prs)
fprintf('Shell side Heat Transfer coefficient = %4.2f W/m^2 C \n',Nho)
fprintf('Overall Heat Transfer coefficient = %4.2f W/m^2 C \n',NUo)
fprintf('Area Required = %8.3f m^2\n',NAreq)
fprintf('percentage excess area = %d % \n',round(Nexcessarea))
fprintf('Shell side pressure drop = %4.2f kPa \n',NdPs)
disp('Shell side pressure drop is less than allowable shell side pressure
drop which is acceptable')
fprintf('Tube Shell side pressure drop = %4.2f kPa \n',AdPt)
fprintf('Tube side pressure drop = %4.2f kPa \n',dPt)
disp('Tube side pressure drop is less than allowable tube side pressure
drop which is acceptable')
Design Calculation of Shell & Tube Exchanger with No Phase Change
Heat Duty = 413.10 kW
Mass Flowrate of Water = 24.667 kg/sec
Mean Temperature Difference = 15.91 Degree Celsius
Area Provided = 64.921 m^2
number of tubes = 427
Tube side flow area = 0.015 m^2
Tube side mass velocity = 1641.67 kg/m2sec
Tube side linear velocity = 1.65 m/sec
Tube side Reynolds number = 30103
Tube side Prandtl number = 4.87
Tube side Heat Transfer coefficient = 6961.77 W/m^2 C
Shell side flow area = 0.007 m^2
Shell side mass velocity = 1380.16 kg/m2sec
Shell side linear velocity = 1.59 m/sec
Shell side Reynolds number = 1053
Shell side Prandtl number = 247.07
Shell side Heat Transfer coefficient = 3373.98 W/m^2 C
Overall Heat Transfer coefficient = 889.97 W/m^2 C
Area Required = 29.179 m^2
percentage excess area = 122 percentage heat transfer area is too high but we will first check shell
side pressure drop
Allowable Shell side pressure drop = 80.00 kPa
Shell side pressure drop = 805.48 kPa
Shell side pressure drop is too high than allowable shell side pressure drop
to decrease it we have to increase Baffle spacing to 200 mm
Shell side flow area = 0.020 m^2
Shell side mass velocity = 480.76 kg/m2sec
Shell side linear velocity = 0.55 m/sec
Shell side Reynolds number = 367
Shell side Prandtl number = 247.07
Shell side Heat Transfer coefficient = 702.06 W/m^2 C
Overall Heat Transfer coefficient = 444.12 W/m^2 C
Area Required = 58.471 m^2
percentage excess area = 11 Shell side pressure drop = 32.07 kPa
Shell side pressure drop is less than allowable shell side pressure drop which is acceptable
Tube Shell side pressure drop = 80.00 kPa
Tube side pressure drop = 50.14 kPa
Tube side pressure drop is less than allowable tube side pressure drop which is acceptable
