Scribd
Upload
24 views5 pages

Matlab Problem 5

Heat exchang problem

Uploaded by

Harvey Specter
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
24 views5 pages

Matlab Problem 5

Heat exchang problem

Uploaded by

Harvey Specter
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 5

%% 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

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy