스테파노님 안녕하세요. ^^
항상 명쾌한 답변에 너무 감사드립니다. 구정명절이 얼마 얼마남지 않았는데 새해 복많이 받으시길
바랍니다.
Flue Gas의 Enthalpy를 계산하는 방법에 대해 궁굼하여 질문을 드려봅니다.
Flue Gas의 조성은 다음과 같습니다.
N2 70.448 Mole%
O2 1.119 Mole%
Ar 0.832 Mole%
CO2 9.035 Mole%
H2O 18.565 Mole%
기준온도는 14 degC에서 주어진 온도 152 degC에서의 Enthapy를 구하여 90 T/H을 곱해 나중에
Heat Flow를 알려고 합니다.
SWEC이라는 Vendor에서의 계산은 아래의 계산식에 의해 14 degC와 152 degC에서의 2지점 enthalpy를
줍니다.)
공식을 보시면 아시겠지만 H = ∫S Cp dt로 적분식이고 온도에 따라 Enthalpy는 선형의 모습을 갖습니다.
하지만 Flue Gas 조성을 고려하면 약 59 degC에서 H2O가 응축이 시작됩니다.
저는 Dew Point 온도이하인 14 degC와 152 degC에서의 2지점 enthalpy 값의 차이를 Enthalpy로 보는
것이 이상하다고 생각이 들었습니다.
엑셀 = GasTXH(152, 14,70.448,1.119,0.832,9.035,18.565) -------> 36.860 kcal/kg
ES_FlueGas에서 계산하면 36.762 kcal/kg 입니다.
그래서 HYSYS를 이용하여 위 조성을 꾸미고 Enthalpy를 14 degC와 152 degC 두 지점의 값을 계산한 후
빼 보았습니다.
그랬더니 103.836 kcal/kg이 나왔습니다.
14 degC Enthalpy = -761.677 kcal/kg
152 degC Enthalpy = -657.842 kcal/kg
△ Enthalpy = 103.836 kcal/kg
36.860 kcal/kg과 103.836 kcal/kg 두개의 값중에 어느 것이 더 정확할까요?
Furnace Efficiency를 계산함에 있어 이부분이 혼란스러웠습니다.
정확한 계산 가이드를 주시면 대단히 감사하겠습니다. ^^
================================================================================
Function GasTXH(temp_given, temp_base, mol_n2, mol_o2, mol_arg, mol_co2, mol_h2o)
Dim coeff_a(1 To 5), coeff_b(1 To 5), coeff_c(1 To 5), coeff_d(1 To 5)
Dim mw(1 To 5), wt(1 To 5)
Dim each_ha, each_hb, each_hc, each_hd, each_h
Dim comp, sum_mol, sum_mw, sum_h
coeff_a(1) = 6.921525928892
coeff_a(2) = 6.905598931009
coeff_a(3) = 4.964694915254
coeff_a(4) = 8.44682086499
coeff_a(5) = 7.967187267909
coeff_b(1) = 0.000264324135
coeff_b(2) = 0.00156659667
coeff_b(3) = 0#
coeff_b(4) = 0.00575849121
coeff_b(5) = 0.000876904142
coeff_c(1) = 0.000000433542006
coeff_c(2) = -0.000000453431059
coeff_c(3) = 0#
coeff_c(4) = -0.00000215927489
coeff_c(5) = 0.000000575161121
coeff_d(1) = -1.17644032E-10
coeff_d(2) = 5.37208504E-11
coeff_d(3) = 0#
coeff_d(4) = 3.05898491E-10
coeff_d(5) = -1.47239368E-10
mw(1) = 28.0134
mw(2) = 31.9988
mw(3) = 39.948
mw(4) = 44.01
mw(5) = 18.0152
temp_given = temp_given * 1.8 + 32
temp_base = temp_base * 1.8 + 32
sum_mol = mol_n2 + mol_o2 + mol_arg + mol_co2 + mol_h2o
mol_n2 = mol_n2 / sum_mol
mol_o2 = mol_o2 / sum_mol
mol_arg = mol_arg / sum_mol
mol_co2 = mol_co2 / sum_mol
mol_h2o = mol_h2o / sum_mol
sum_mw = mw(1) * mol_n2 + mw(2) * mol_o2 + mw(3) * mol_arg + mw(4) * mol_co2 + mw(5) * mol_h2o
wt(1) = mw(1) * mol_n2 / sum_mw
wt(2) = mw(2) * mol_o2 / sum_mw
wt(3) = mw(3) * mol_arg / sum_mw
wt(4) = mw(4) * mol_co2 / sum_mw
wt(5) = mw(5) * mol_h2o / sum_mw
If temp_given < 0 Then GasTXH = -7000
If temp_base < 0 Then GasTXH = -7000
If temp_base > temp_given Then GasTXH = -7000
sum_h = 0
For comp = 1 To 5
each_ha = coeff_a(comp) * (temp_given - temp_base)
each_hb = coeff_b(comp) / 2 * (temp_given ^ 2 - temp_base ^ 2)
each_hc = coeff_c(comp) / 3 * (temp_given ^ 3 - temp_base ^ 3)
each_hd = coeff_d(comp) / 4 * (temp_given ^ 4 - temp_base ^ 4)
each_h = (each_ha + each_hb + each_hc + each_hd) / mw(comp) * wt(comp) * 0.556
sum_h = sum_h + each_h
Next comp
GasTXH = sum_h
End Function
================================================================================
