%% input data for box type heat
exchanger T1=45; %C inlet temp of air
T2=12; %C outlet temp of air
ml=0.9448; %kg/sec massflowrate of
air Dl=12939; %kg/m3 density of air
Cpl=0.25; %kJ/kg C specific heat of
air
kl=0.02575; %w/m C thermal conductivity of
air P=500; %mmH2Og operating pressure of
air Vl=0.02; %cP viscosity of air
hod=2000; %kcal/hrm2 C dirt factor for shell
side Co=120; %condensation load as sensible
load M=29; %kg/kmol molar mass of air
kwater=0.628; %w/m C thermal conductivity of
water t1=6; %C inlet temp of water
t2=9; %C outlet temp of water
Cpw=1; %kcal/hrm C specific heat of water
%%specification of air cooled heat
exchanger fi=19.05; %mm Fin inside
diameter
fd=38.1; %mm Fin outside
diameter do=19.05; %mm tube
o.d. di=15.748; %mm tube i.d.
L=1.2; %m tube length
Tav=7.5; %C t average
teprature th=1.65; %mm tube
wall thickness F=7; %mm fins
per inch
tb=43; %mm tube pitch
hid=4000; %kcal/hrm2 C dirt factor for tube side
Km=205.87; %kcal/hrm C thermal conductivity of fin
metal Pr=0.699; %Prandalt number
FI=0.97; %Efficiency of
fin Ft=0.98; %Friction
factor Np=28; %Number of
passes
%% heat load & energy balance
phit=2.2*ml*Cpl*(T1-T2); %kg/s heat
load
mw= phit/(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
%% eavaluation of air side HTC
Af=3.14/4*(fd^2-fi^2) *2*F*(1000/25.4)*(10^-6); %m2 fin area per unit
length Abl=((3.14*do)- (3.14*do*0.889*10^(-3) *7*(1000/25.4*10^-
3)))*(10^-3); %m2 bare tube area
PP=2*(fd-do*10^-3)*(7*(1000/25.4*10^-3))+2*(1-0.889*10^(-3)
*7*(1000/25.4)); %m
Projected Perimeter
%% eavaluation of air side flow area
as=((1.21*1.2)-(np*do*10^-3*L))-np*((fd-do)*10^-3)*0.899*(10^-
3)*1.2*7*(1000/25.4); %m2 Air side flow area
Gs=ml/as; %kg/m2sec air side mass velocity
Res=(0.03657*Gs/Vl*1000);
jh=(0.04536*2632.21)^Pr; %Tube side factor for heat transfer
cofficeint hf=jh*(kl/0.03657)*(Pr)^0.33;% efficency
omega=0.97
;
AF=0.565;
AO=0.0717;
AI=3.14*di*10^-3*1.2*1;
Ddhf=((omega*AF)+(AO))*((hf)/AI);
disp("Ddhf");
%% eavaluation of dirty fin side
HTC Dhf = 1/((1/hf)+(1/hod));
%% Calculation of water side HTC
at=Np*3.14/4*di^2*(10^-6); %m2 tube side flow
area AT=mw/at; %kg/m2s tube side mass
velocity
Re=di*AT*(10^-3)/(0.73*10^-3); %tube side Reynolds
number Pr= 0.73*10^(-3)*1*3600/0.5355; %Tube side
Prandalt number
hi=(0.023*(Re^(0.8))*(Pr^(0.33))*(0.5355))/di*(10^3); %Tube side heat
transfer cofficient
%% Evaluation of Overall HTC
Ui=1/(1/Ddhf+1/hi+1/hid); %Kcal/hrm2 C Overall Heat Transfer Cofficient
Areq=(phit)*3600/(Ui*MTD); %m2 area required
Apro=1.1*Areq; %m^2 area
provided Ao=313.26; %m^2 total
outside area
Ex=(1-(Areq/Apro))*100;%excess
area TT=18*28;% total number of
fin tube aAF=Af*TT;% total fin
area aAbl=Abl*TT;% total bare
tube area aAI=AI*TT;
TOtal=aAF+aAbl;%total outsife area
%% Evaluation of air side pressure
drop WL=1.21*1.2;
x=37.239;
Vfv=WL*x*(10^-3)*(Np*(3.14/4)*(do^2)*L*(10^-6))-(Np*(3.14/4)*(fd^2-
do^2)*(10^- 6)*(10^-6)*7*(1000/25.4)*(0.889*10^-3));%Net free volume
Dx=37.239; %mm
Dev=0.00874; %m Volumetric equivalent
diameter Rep= Dev*AT/Vl; %Reynolds number for
Pressure Drop Jf= 0.9355;
dPs=(Jf*(Gs^2)*0.67)/(Dev*1.29)*((Dev/0.043)^(0.4))*((0.043/0.043)^0.6);%
air side pressure drop
%% Result
disp('Design Calculation of Box type
Exchanger') fprintf('Heat Load = %4.4f kg\
n',phit)
fprintf('Mass Flowrate of Water = %1.4f kg/sec\n',mw)
fprintf('Mean Temperature Difference = %2.4f Degree Celcius\
n',MTD) fprintf('Fin area per unit length= %1.5f m^2\n',Af)
fprintf('bare tube area= %1.4f m^2\n',Abl)
fprintf("Projected Perimeter=%1.5f m\n" ,PP)
fprintf('Air side flow area = %1.4f m^2\n',as)
fprintf('Air side mass velocity = %1.3f kg/m2sec \
n',Gs) fprintf('Air side Reynolds number = %d\
n',round(Res))
fprintf('Tube side Factor for heat transfer cofficient = %d\
n',round(jh)) fprintf("efficency= %2.2f kcal/hrm^2 C\n" ,hf)
fprintf('dirty fin side Heat Transfer coefficient = %2.3f kcal/hrm^2 C\
n',hf') fprintf('Air side Heat Transfer coefficient = %4.2f W/m^2 C \
n',Af)
disp('Air side HeatTransfer Cofficient')
fprintf('Tube side flow area = %1.6f m^2\
n',at)
fprintf('Tube side mass velocity = %4.3f kg/m2sec \
n',AT) fprintf('Tube side Reynolds number = %d\
n',round(Re)) fprintf('Tube side Prandel number =
%1.2f \n',Pr)
fprintf('Tube side Heat Transfer coefficient = %4.2f
Kcal/hr*m^2*C \n',hi) disp('Tube side pressure drop')
fprintf('Overall Heat Transfer coefficient = %4.2f W/m^2 C \
n',Ui) fprintf('Area Required = %2.2f m^2\n',Areq)
fprintf('Area Provided = %2.2f m^2\n',Apro)
fprintf('excess area = %d \n',round(Ex))
fprintf('Air side net free volume = %1.5f m^3 \
n',Vfv)
fprintf('Reynolds number for pressure drop = %d\
n',round(Rep)) fprintf('Jf= %1.4f \n' ,Jf)
fprintf('Air side pressure drop = %2.3f kPa \
n',dPs) disp('Air side pressure drop')
OUTPUT
Design Calculation of Box type Exchanger
Heat Load = 17.1481 kgMass Flowrate of Water = 5.7160
kg/sec Mean Temperature Difference = 16.4085 Degree
Celcius
Fin area per unit length= 0.47106
m^2 bare tube area= 0.0598 m^2
Projected Perimeter=22.49950
m Air side flow area = 0.6533
m^2
Air side mass velocity = 1.446
kg/m2sec Air side Reynolds number =
2644
Tube side Factor for heat transfer cofficient
= 28 efficency= 17.71 kcal/hrm^2 C
dirty fin side Heat Transfer coefficient = 17.707
kcal/hrm^2 C Air side Heat Transfer coefficient = 0.47
W/m^2 C
Air side HeatTransfer Cofficient
Tube side flow area = 0.005451
m^2
Tube side mass velocity = 1048.617
kg/m2sec Tube side Reynolds number =
22621
Tube side Prandel number = 4.91
Tube side Heat Transfer coefficient = 4025.90
Kcal/hr*m^2*C Tube side pressure drop
Overall Heat Transfer coefficient = 169.33
W/m^2 C Area Required = 22.22 m^2
Area Provided = 24.44
m^2 excess area = 9
Air side net free volume = 0.00052
m^3 Reynolds number for pressure
drop = 458 Jf= 0.9355
Air side pressure drop = 61.467 kPa
