PHT Assignment#3 Final Excel Sheet
PHT Assignment#3 Final Excel Sheet
Correction factor for end effects (Ψ) For Inner Pipe (Ψp)
1.02346201139361
And by defination
Nu=(hDfc)/k Euqation9
Similarly for Isothermal Conditions (means wall temp.=mean temp of fluid) Equation9 becomes
Nu'=h'Dfc/k
Rearranging gives
h'=Nu'k/Dfc Equation11
For Inner Pipe
h'
h'
h'i
279.010051696219
For Isothermal condition H.T coefficient based on outer surface (of inner pipe)
U'=h'/2
For Inner Pipe
U'
U'i
139.505025848109
Now As
Wall temperature (θw)=(h'iotc+h'oTc)/(h'io+h'o) Euqation12
where
Now we are able to calculate Viscosity Number (ɸ) and then Correct film coefficient,Correct H.T coefficient.
And then using Kern's model i.e. Equation 2 we can find Area (Ao) from which Length of H.E can be determined and then fin
For Inner Pipe OR for (cf)
Viscosity Number (ɸ) ɸp
1.15110188421556
Now using ho and hio we can find I/Uco using Equation 5 i.e. 1/Uco=1/ho + 1/hio
1/Uco 0.00871221774324467
Now using Equation 7 i.e. Rsio= Rsi * Ao/Ai = Rsio * do/di Calculate Rsio
Rsio 0.00120289855072464
Now calculate overall dirt factor based on outer surface of inner pipe
RDo=Rsio + Rso 0.00220289855072464
Now calculate Wall resistance based on outer surface of inner pipe (Rwo)
Rwo=(do/2kw)*(ln(do/di))
As
Ao=∏doL
Therefore
L=Ao/∏do 207.354446789159
L=NHPLHP
NHP=L/LHP 6.47982646216123
and
Head loss due to skin friction (hsf)=Kf*hk
Head loss due to form friction (hff)=KL*hk
Total Head Loss ∆h= hsf+hff
Friction Parameter (Kf)=(fd*Lfc)/Dfc
where
Darcy friction factor (fd) for turbulent flow=0.01
K'L
For Annulus
1
Units
30 K or °F
40 K or °F
For (hf) in (W/m.K) For (cf) in (Btu/hr.ft. °F) For (hf) in (Btu/hr.ft. °F)
0.168916578204 0.082859970575756 0.097639640580347
Value of pi 3.141592654
Note: ɸ is taken as 1
For Annulus
179.30807906805
Euqation13
Euqation14
Equation15
Units (°F)
Units (K)
For Annulus OR (for (hf)) For Inner Pipe OR (for (cf)) For Annulus OR (for (hf))
μwa OR μwhf in (cp) μwp OR μwcf in (lbm/ft.hr) μwa OR μwhf in (lbm/ft.hr)
2.3754819062615 1.14515588637368 5.74866621315282
hr.oF.ft2/Btu
hr.oF.ft2/Btu
hr.oF.ft2/Btu
(Btu/hr.ft. °F)
hr.oF.ft2/Btu
hr.oF.ft2/Btu
Btu/hr.ft^2. °F
ft^2
ft
192 ft
Units
ft
8189.79851341796 2100<Re<10000
Transition zone
n friction (hsf)=Kf*hk
m friction (hff)=KL*hk
For Annulus
fda
0.0379932548459893
For Annulus
Kfa
215.077295246092
For Annulus
hsfa (ft)
83.161820980174
KL K'L=0 if U<3ft/s
6 K'L=1 if U>3ft/s
For Annulus
hffa (ft)
2.31996094850514
For Annulus
∆ha (ft)
85.4817819286791
For Annulus
∆P (Psia)
36.8571665875581
K
K
Ua(ft/s)
4.98806654450414
OK
Turbulent flow
Property Name
Specific Heat Capacity (Cp)
Specific Heat Capacity (Cp)
Density (ρ)
Density (ρ)
Dynamic Viscosity ( μ)
Dynamic Viscosity ( μ)
Thermal Conductivity (k)
Thermal Conductivity (k)
Our calculations for Heat Exchange will be based on outer surface of inner pipe
For the calculations of No.of hairpins we have to calculate Length of exchanger as
For the calculations of Length of exchanger as we have to find Area (Ao)
And for the calculation of outer area we shall use Kern's model
Basic Design Equation or Kern's model
where
LMTD=(∆θ1-∆θ2)/ln(∆θ1/∆θ2)
As we know
Law of Heat exchange
Qdot
Fulid Allocation
Note:Now as information regarding fluid allocation is not given in problem so we have to decide which fluid should be place in
From properties value calculations it can be seen that viscosity of cold fluids is comparatively low, but mass flow rate of cold f
And More viscous fluid should be placed outside tube/inner pipe.
Now from Equation 2
Inorder to calculate (Ao) using this equation we need Qdot,UDO and LMTD, out of which Qdot and LMTD have been
As
where
1st term in Equation 4 is
and
Empirical Equation
Cp= A+BT+CT^2+DT^3
Cp= A+BT+CT^2+DT^3
ρ=AB-(1-T/C)^D
ρ=AB-(1-T/C)^D
log μ=A+(B/T)+CT+DT2
log μ=A+(B/T)+CT+DT2
logk=A+B(1-T/C)2/7
logk=A+B(1-T/C)2/7
UDO
Ao
LMTD
34.7605949678221
Q̇ absorbed = Q̇ evolved
ṁcfCpcf∆t = ṁhfCphf∆T
ṁcfCpcf∆t
9224.17338677156
n is not given in problem so we have to decide which fluid should be place inside and which fluid should be placed outside the inner pipe.T
hat viscosity of cold fluids is comparatively low, but mass flow rate of cold fluid i.e. Benzene is higher than that of hot fluid i.e. Aniline. Th
ube/inner pipe.
Qdot=UDoAo(LMTD)
n we need Qdot,UDO and LMTD, out of which Qdot and LMTD have been calculated hence only UDO is to be calculated
1/UDO=1/Uco + RDo + Rwo
Uco
1/Uco=1/ho + 1/hio
ho
hio
hio= hi * Ai/Ao = hi * di/do
RDo
RDo=Rsio + Rso
Rso
Rsio
Rsio= RSi * Ao/Ai = RSi * do/di
Rwo
Rwo=(do/2kw)*(ln(do/di))
Specified Units A B
T in (K) and Cp in (J/gmol.K) -31.662 1.3043
T in (K) and Cp in (J/gmol.K) 63.288 0.9896
T in (K) and ρ in (g/mL) 0.3009 0.2677
T in (K) and ρ in (g/mL) 0.3119 0.25
T in (K) and μ in (cp) -7.4005 1181.5
T in (K) and μ in (cp) -13.8625 2510.9
T in (K) and k in (W/mK) -1.6846 1.052
T in (K) and k in (W/mK) -1.3485 0.6888
Euqation1
Euqation2
overall H.T/design coefficient based on outer surface of inner pipe under dirty conditions
H.T area based on outer surface of inner pipe
Mean temperature difference or overall driving force
K or °F
Euqation3
255197.454155659
lbm/hr
be place inside and which fluid should be placed outside the inner pipe.This decision can be made on the basis of property values,flow rates
e of cold fluid i.e. Benzene is higher than that of hot fluid i.e. Aniline. Therefore Cold fluid should be placed inside the inner pipe and Hot f
have been calculated hence only UDO is to be calculated
Euqation4
overall H.T coefficient based on outer surface of inner pipe under clean conditions
Euqation5
film coefficinet of outside fluid based on outer surface of inner pipe
film coefficinet of inside fluid based on outer surface of inner pipe
Euqation6
μ (cP) at θmax
μ (cP) at θmin
Properties Calculations at Caloric Temperature
Specific Heat Capacity
Cp
Thermal Conductivity
k
Density
ρ
Dynamic Viscosity
μ
From given NPS of inner and outerpipe note Inner and Outer Diameters of both pipies
Reference PHT by Kern 1st Edn. Pg844
Flow Channel
Inner Pipe
Outer Pipe
Name of Parameter
And by defination
Nu=(hDfc)/k
Similarly for Isothermal Conditions (means wall temp.=mean temp of fluid) Equation9 becomes
Nu'=h'Dfc/k
Rearranging gives
h'=Nu'k/Dfc
h'
For Isothermal condition H.T coefficient based on outer surface (of inner pipe)
U'=h'/2
U'
Now As
Wall temperature (θw)=(h'iotc+h'oTc)/(h'io+h'o)
Now calculate dynamic viscosity (μ) at wall temperature (θw) using Empirical Correlation
μw
Now we are able to calculate Viscosity Number (ɸ) and then Correct film coefficient,Correct H.T coefficient.
And then using Kern's model i.e. Equation 2 we can find Area (Ao) from which Length of H.E can be determined and then fin
Viscosity Number (ɸ)
h=h'*(ɸ^(n))=(ɸ)^(0.14)
Now using ho and hio we can find I/Uco using Equation 5 i.e. 1/Uco=1/ho + 1/hio
1/Uco
Now using Equation 7 i.e. Rsio= Rsi * Ao/Ai = Rsio * do/di Calculate Rsio
Rsio
Now calculate overall dirt factor based on outer surface of inner pipe
RDo=Rsio + Rso
Now calculate Wall resistance based on outer surface of inner pipe (Rwo)
Rwo=(do/2kw)*(ln(do/di))
As
Ao=∏doL
Therefore
L=Ao/∏do
L=NHPLHP
NHP=L/LHP
and
where
Dynamic Pressure Pd
For Annulus
Hence Devide annulus side fluid flow into two parallel streams.
H.T Calculations
Benzene (cf)
At 91.4°F (306.1457793K) calooric temperature (inside Inner Pipe)
Note:As wall temperature (θw) is unkwon yet,therefore viscosity number can't be calculated right now
Hence take ɸ=(μ/μw)=1 i.e. isothermal condition and after finding wall temperature use
Nu'=0.023(Re^0.8)(Pr)^1/3*(ɸ)^n*(1+x^-2/3)
(Seder-Tate Equation for Turbulent flow)
389.553139256751
As Nu'=h'Dfc/k
Hance h'=Nu'k/Dfc
Viscosity Number Correction factor for turbulent and transition flow is (ɸ^(n))=(ɸ)^(0.14)
1.01409165186426
Now Correct Nusselt Number (Nu) for Non-Isothermal Condition during Turbulent or Transition flow
Nu=Nu'*(ɸ^(n))
395.042586477786
Now Correct Fil Coefficient (h) for Non-Isothermal Condition during Turbulent or Transition flow
hi=hi'*(ɸ^(n))
(Btu/hr.ft^2. °F)
284.162151659408
overall H.T coefficient based on outer surface of inner pipe under clean conditions (Uco)
Uco= 1/(1/ho + 1/hio + Rwo)
where Rwo=(do/2kw)*(ln(do/di))
Hence, Uco
Over all Design coefficient based on outer surface of inner pipe Under dirty/fouling condition (Udo)
1/Udo=1/Uco+Rdo
Hence Udo=1/((1/Uco)+Rdo))
Designe Equation
Qdot=UDoAo(LMTD)
Rearranging gives Ao= Qdot/Udo*(LMTD)
Ao
Dynamic Pressure Pd
For Annulus
Ado=Qdot/Udo∆ΘLM
ACo=Qdot/UCo∆ΘLM
Ao=∏doLHPNHP
∆θ1=(T1-t2) 30
∆θ2=(T2-t1) 40
De=(Di^2-do^2)/do
0.076149046184739
Euqation8
end effects during turbulent flow
Euqation9
of inner pipe)
Euqation12
where tc is caloric temp. of (cf)
Tc is caloric temp. of (hf)
h'io is film coefficient of inner fluid based on outer surfac
h'o is film coefficient of outer fluid based on outer surfac
Equation15
1.33333333333333
0.00619078575495324
165.530383958111
1.02145320183428
161.530383958111
0.523426887906908
Fc>0.5
Therefore θc>ΘAM
91.4056132744145 tc>tAM
306.153118485786 K
126.171344395345 Tc>TAM
325.46741355297 K
co=1/ho + 1/hio
0.0086782150766756 hr.oF.ft2/Btu
0.00120289855072464 hr.oF.ft2/Btu
0.00220289855072464 hr.oF.ft2/Btu
0.00135930235510152 hr.oF.ft2/Btu
0.0122404159825018 hr.oF.ft2/Btu
81.696569906574 Btu/hr.ft^2. °F
89.9491911089788 ft^2
206.976265883021 ft
6.4680083088444
6
Lfc 192
De'=Di-do Units
0.0339166666666667 ft
Re'=ρhf*De'Uhf/μhf 8314.79940743734
Equation 16
where
∆Psf=ρghsf
∆Pff=ρghff
32.174
K'L KL
1 6
Aniline (hf)
At 126.16°F (325.4612976K) calooric temperature (in annulus)
Nu'=0.116((Re^(2/3))-125)((Pr)^(1/3))*(ɸ)^n*(1+x^-2/3))
Hausen Equation for Transition Flow
111.073557464373
°F K
104.590888504193 313.478271391218
y were already calculated in sheet1 on the basis of which we are performing our calculations.Results don't change appreciably as it can be s
μw (cP)
2.51937189640719
y number= ɸ= (μ/μw)
0.746206721288496
on flow is (ɸ^(n))=(ɸ)^(0.14)
0.959843176647664
0.00135930235510152 hr.oF.ft2/Btu
77.4446716677012 (Btu/hr.ft^2. °F)
0.00220289855072464 hr.oF.ft2/Btu
255.589123814076 ft
7.98716011918989
8 Inner pipe
4 Annulus
Lfc 256
Lfc 128
K'L KL
0 0
K'L KL
1 4
nd each consists of 4 hairpins. Cold fluid i.e. Benzene will flow (inside inner pipe) in series through both H.E, while hot fluid i.e. An
111.075706411296 ft^2
94.8876174594104 ft^2
111.254267853653 ft^2
Fraction Percentage
0.0016075652194909 0.16075652194909
0.172484575252868 17.2484575252868
306.1457793
325.4612976
Units
lbm/hr
O
F 288.705555555556
O
F 322.038888888889
O
F 288.705555555556
O
F 338.705555555556
O
F 310.927777777778
O
F 283.15
hr.oF.ft2/Btu
hr.oF.ft2/Btu
psia
ft
ft
in
in
O
F 305.372222222222
F
O
324.816666666667
g/gmol
g/gmol
Units
K or °F
K or °F
Value of pi 3.141592654
Units
ft
ft^2
Note:Ap>Aa
lbm/ft^2.hr
Minor Checks
ft/hr Up(ft/s)
4.95273398434034
3≤U≤8 ft/s
Re>10000
t of inner fluid based on outer surface (of inner pipe) under isothermal conditions
of outer fluid based on outer surface (of inner pipe) under isothermal conditions
Units (°F)
Units (°F)
ft
2100<Re<10000
Transition zone
Units
ft
K'L=0 if U<3ft/s
K'L=1 if U>3ft/s
OK
n't change appreciably as it can be seen from the table of Rough work.
Lfc
ft
ft
h both H.E, while hot fluid i.e. Aniline will flow (on annulus side) in two parallel streams through H.E means half aniline flow will p
Allowable limits
Slightly below the allowable limit (OK) 5≤∆d≤20%
Slightly below the allowable limit (OK) 20≤∆c≤40%
K
K
K
K
K
K
K
K
mp for calculations.
For (cf) in Btu/lbm.°F
0.425732726225996
Inner Diameter (ID)
(ft)
0.115
0.17225
Ua(ft/s)
4.99274548104421
OK
Turbulent flow
K
K
Property Name
Specific Heat Capacity (Cp)
Specific Heat Capacity (Cp)
Density (ρ)
Density (ρ)
Dynamic Viscosity ( μ)
Dynamic Viscosity ( μ)
Thermal Conductivity (k)
Thermal Conductivity (k)
Our calculations for Heat Exchange will be based on outer surface of inner pipe
For the calculations of No.of hairpins we have to calculate Length of exchanger as
For the calculations of Length of exchanger as we have to find Area (Ao)
And for the calculation of outer area we shall use Kern's model
Basic Design Equation or Kern's model
where
LMTD=(∆θ1-∆θ2)/ln(∆θ1/∆θ2)
As we know
Law of Heat exchange
Qdot
Fulid Allocation
Note:Now as information regarding fluid allocation is not given in problem so we have to decide which fluid should be place in
From properties value calculations it can be seen that viscosity of cold fluids is comparatively low, but mass flow rate of cold f
And More viscous fluid should be placed outside tube/inner pipe.
Now from Equation 2
Inorder to calculate (Ao) using this equation we need Qdot,UDO and LMTD, out of which Qdot and LMTD have been
As
where
1st term in Equation 4 is
and
Empirical Equation
Cp= A+BT+CT^2+DT^3
Cp= A+BT+CT^2+DT^3
ρ=AB-(1-T/C)^D
ρ=AB-(1-T/C)^D
log μ=A+(B/T)+CT+DT2
log μ=A+(B/T)+CT+DT2
logk=A+B(1-T/C)2/7
logk=A+B(1-T/C)2/7
UDO
Ao
LMTD
34.7605949678221
Q̇ absorbed = Q̇ evolved
ṁcfCpcf∆t = ṁhfCphf∆T
ṁcfCpcf∆t
9227.55425884813
n is not given in problem so we have to decide which fluid should be place inside and which fluid should be placed outside the inner pipe.T
hat viscosity of cold fluids is comparatively low, but mass flow rate of cold fluid i.e. Benzene is higher than that of hot fluid i.e. Aniline. Th
ube/inner pipe.
Qdot=UDoAo(LMTD)
n we need Qdot,UDO and LMTD, out of which Qdot and LMTD have been calculated hence only UDO is to be calculated
1/UDO=1/Uco + RDo + Rwo
Uco
1/Uco=1/ho + 1/hio
ho
hio
hio= hi * Ai/Ao = hi * di/do
RDo
RDo=Rsio + Rso
Rso
Rsio
Rsio= RSi * Ao/Ai = RSi * do/di
Rwo
Rwo=(do/2kw)*(ln(do/di))
304.8058330547
324.344675693732
Specified Units A B
T in (K) and Cp in (J/gmol.K) -31.662 1.3043
T in (K) and Cp in (J/gmol.K) 63.288 0.9896
T in (K) and ρ in (g/mL) 0.3009 0.2677
T in (K) and ρ in (g/mL) 0.3119 0.25
T in (K) and μ in (cp) -7.4005 1181.5
T in (K) and μ in (cp) -13.8625 2510.9
T in (K) and k in (W/mK) -1.6846 1.052
T in (K) and k in (W/mK) -1.3485 0.6888
Euqation1
Euqation2
overall H.T/design coefficient based on outer surface of inner pipe under dirty conditions
H.T area based on outer surface of inner pipe
Mean temperature difference or overall driving force
K or °F
Euqation3
255439.635735598
lbm/hr
be place inside and which fluid should be placed outside the inner pipe.This decision can be made on the basis of property values,flow rates
e of cold fluid i.e. Benzene is higher than that of hot fluid i.e. Aniline. Therefore Cold fluid should be placed inside the inner pipe and Hot f
have been calculated hence only UDO is to be calculated
Euqation4
overall H.T coefficient based on outer surface of inner pipe under clean conditions
Euqation5
film coefficinet of outside fluid based on outer surface of inner pipe
film coefficinet of inside fluid based on outer surface of inner pipe
Euqation6