Example 9.6-1
Example 9.6-1
Example 9.6-1
The inlet humidity ratio (a,in) and the dewpoint associated with the incoming air (Tdp, the
temperature at which the relative humidity reaches 100% and therefore water will begin to
condense) are computed:
omega_a_in=HumRat(AirH2O,P=p_a,T=T_a_in,R=RH_a_in)
T_dp=Temperature(AirH2O,P=p_a,w=omega_a_in,R=1)
T_dp_C=converttemp(K,C,T_dp)
which leads to Tdp = 295.3 K (22.19C). Because the entering water temperature is less than Tdp,
water will condense onto the coil surface provided that the cooling coil is sufficiently large.
The water-side thermal resistance is evaluated. The inner diameter of the tube is:
Din = Dout 2 th
and the total length of tube is:
Rin =
1
Din Ltube hw
The area-specific fouling resistance for closed-loop treated water ( Rf ,in ) is obtained using the
FoulingFactor function in EES and used to compute the fouling resistance on the water side:
R f ,in =
Rf ,in
Din Ltube
"Fouling resistance"
R``_f_in=FoulingFactor('Closed-loop treated water') "fouling factor on inner surface of tube"
R_f_in=R``_f_in/(pi*D_in*L_tube)
"fouling resistance on inner surface of tube"
The conductivity of copper (km) is obtained at the average temperature using the internal property
routine in EES and used to compute the conduction resistance through the tube:
Rcond
D
ln out
D
= in
2 km Ltube
"Conduction resistance"
k_m=k_('Copper',T_avg)
R_cond=ln(D_out/D_in)/(2*pi*k_m*L_tube)
"tube conductivity"
"tube resistance"
where a,in is the density of the air at the inlet. The heat exchanger geometry corresponds to
compact heat exchanger core fc_tubes_s80-38T. The compact heat exchanger function
CHX_h_finned_tube is used to evaluate the dry heat transfer coefficient on the air side ( ha ):
"Compact heat exchanger correlation"
rho_a_in=density(Air,T=T_a_in,P=P_a)
"density of air at inlet"
m_dot_a=V_dot_a*rho_a_in
"mass flow rate of air"
TypeHX$='fc_tubes_s80-38T'
"heat exchanger identifier name"
Call CHX_h_finned_tube(TypeHX$, m_dot_a, W*H, 'Air',T_avg, P_a:h_bar_a)
Dout
As , fin ,tot = 2
H
L
N
N
t
row
t
col
,
,
p fin
4
p fin
The dry coil/wet coil analysis divides the heat exchanger into a separate dry coil and wet coil. In
the dry coil, the temperature of the moist air is reduced to the dew point temperature; the exit
temperature from the dry coil (Ta,dc,out) is therefore the dew point temperature.
T_a_dc_out=T_dp
The fraction of the total heat exchange surface required by the dry coil section is determined.
The specific heat capacity of the moist-air on a per mass of dry air basis ( ca ) is calculated using
EES' internal psychrometric function.
c``_a=cP(AirH2O,P=p_a,T=T_a_in,R=RH_a_in) "specific heat capacity of moist air in dry coil"
The specific heat capacity of the cooling water (cw) is calculated using EES' internal property
routine:
c_w=cP(Water,P=p_w,T=T_avg)
The capacitance rate of the air and the water in the dry-coil are:
C a , dc = m a ca
C w = m w cw
C_dot_a_dc=m_dot_a*c``_a
C_dot_w=m_dot_w*c_w
Provided that the cooling coil is sufficiently large so that the dew point temperature is reached,
the rate of heat transfer in the dry coil section of the heat exchanger is given by:
dc =
q_dot_dc=C_dot_a_dc*(T_a_in-T_a_dc_out)
C_dot_min_dc=Min(C_dot_a_dc,C_dot_w)
q_dot_dc_max=C_dot_min_dc*(T_a_in-T_w_in)
eff_dc=q_dot_dc/q_dot_dc_max
qdc
qdc ,max
"heat transfer rate in dry coil"
"minimum heat capacity rate in dry coil"
"maximum possible heat transfer rate in dry coil"
"effectiveness of dry coil"
The number of transfer units in the dry coil section (NTUdc) is obtained using the -NTU solution
for a cross-flow heat exchanger (with the air unmixed and the water mixed) accessed using the
HX function in EES. The conductance required in the dry-coil section is:
UAdc = NTU dc C min ,dc
NTU_dc=HX('crossflow_one_unmixed', eff_dc, C_dot_a_dc, C_dot_w, 'NTU')
"number of transfer units in dry coil"
UA_dc=NTU_dc*C_dot_min_dc
"conductance in dry coil"
The total resistance to heat transfer in the dry coil is the inverse of the conductance:
Rdc =
1
UAdc
and is made up of the resistance to convection on the water side, a fouling resistance, resistance
to conduction through the tube, and a convection resistance on the air side (assuming that the fin
efficiency is 100%):
Rdc =
1
ha Atot Fdc
where Fdc is the fraction of the total heat exchanger required by the dry coil.
R_dc=1/UA_dc
R_dc=(R_in+R_f_in+R_cond)/F_dc+1/(h_bar_a*A_tot*F_dc)
which leads to Fdc = 0.329 (i.e., 33% of the coil is dry, according to this analysis). The
remainder of the heat exchanger coil is wet:
Fwc = 1 Fdc
"Wet coil analysis"
F_wc=1-F_dc
The wet coil is analyzed using the concepts of a saturation specific heat capacity and wet coil
conductance. The saturation specific heat capacity is used to account for the latent heat
associated with the condensing water vapor and is calculated using Eq. (9-139):
ia (Tdc ,out , p,RH = 1) ia (Ta ,out , p,RH = 1)
ca,sat =
(Tdc ,out Ta ,out )
Initially, a value for Ta,out is assumed; this value is subsequently commented out in order to close
the solution:
T_a_out=T_w_in
"an initial assumption for the saturation specific heat capacity"
c``_a_sat=(enthalpy(AirH2O,T=T_a_dc_out,R=1,P=p_a)-&
enthalpy(AirH2O,T=T_a_out,R=1,P=p_a))/(T_a_dc_out-T_a_out)
"saturation specific heat capacity"
and the minimum capacitance rate in the wet coil ( C min , wc ) is computed using the Min function:
C_dot_a_wc=m_dot_a*c``_a_sat
C_dot_min_wc=Min(C_dot_a_wc,C_dot_w)
The conductance of the wet coil is computed by augmenting the air side heat transfer coefficient
by the ratio ca, sat / ca in order to account for the enhancement of the energy transfer rate due to
mass transfer. The thermal resistance of the wet coil section is:
Rwc =
ca
ha ca, sat Atot Fwc
1
Rwc
The effectiveness of the wet coil (wc) is obtained using the -NTU solution for a cross-flow heat
exchanger, obtained using the HX function in EES:
NTU_wc=UA_wc/C_dot_min_wc
"number of transfer units in the wet coil"
eff_wc=HX('crossflow_one_unmixed', NTU_wc, C_dot_a_wc, C_dot_w, 'epsilon')
"effectiveness of wet coil"
The heat transfer rate in the wet coil section is obtained from the effectiveness:
qwc = wc C min , wc (Tdc ,out Tw,in )
q_dot_wc=eff_wc*C_dot_min_wc*(T_a_dc_out-T_w_in)
The temperature of the air leaving the wet coil is obtained using an energy balance on the airside:
Ta ,out = Tdc ,out
qwc
Ca , wc
Before calculating the value of Ta,out, the guess values for the problem are updated (select Update
Guesses from the Calculate menu) and the initial guess for Ta,out is commented out:
{T_a_out=T_w_in}
"an initial assumption for the saturation specific heat capacity"
T_a_out=T_a_dc_out-q_dot_wc/C_dot_a_wc
"temperature of air leaving wet coil"
T_a_out_C=converttemp(K,C,T_a_out)
"in C"
The total heat transfer in the cooling coil predicted by the wet coil/dry coil method is the sum of
the heat transferred in the wet and dry sections:
qwc / dc = qdc + q wc