Skip to content

Теоретические основы тепловых расчетов ​

1. Законы сохранения и фундаментальные принципы ​

Теплотехнические расчёты сушильных агрегатов и печей в программном комплексе «UrFU HeatBalance» базируются на совместном решении уравнений материального (баланс массы сухого вещества и влаги) и энергетического (первый закон термодинамики) балансов.

Для установившегося (стационарного) или интегрально осреднённого квазистационарного технологического процесса суммарный приход энергии в рабочий объём агрегата строго равен суммарному расходу:

∑Qприход=∑Qрасход\sum Q_{\text{приход}} = \sum Q_{\text{расход}}

Невязка баланса ΔQ\Delta Q, возникающая вследствие погрешностей исходных измерений или неучтённых факторов, определяется выражением:

ΔQ=∑Qприход−∑Qрасход\Delta Q = \sum Q_{\text{приход}} - \sum Q_{\text{расход}}

В верифицированных инженерных расчётах относительная величина невязки не должна превышать ±3%\pm 3\%:

δ=∣ΔQ∣∑Qприход⋅100%≤3%\delta = \frac{|\Delta Q|}{\sum Q_{\text{приход}}} \cdot 100\% \le 3\%

2. Теория и уравнения горения газообразного топлива ​

В качестве основного энергетического топлива для камерных, конвейерных печей и сушильных барабанов в системе принят природный газ стандартного сетевого состава.

2.1. Низшая теплота сгорания топлива ​

Низшая теплота сгорания сухого газового топлива при нормальных физических условиях (T0=273,15 KT_0 = 273{,}15 \text{ K}, P0=101,325 кПаP_0 = 101{,}325 \text{ кПа}) рассчитывается по формуле Менделеева:

Qнр=127,5⋅CO+107,9⋅H2+234,1⋅H2S+358,2⋅CH4+637,5⋅C2H6+912,7⋅C3H8+1186,5⋅C4H10+1461,1⋅C5H12Q_н^р = 127{,}5 \cdot \text{CO} + 107{,}9 \cdot \text{H}_2 + 234{,}1 \cdot \text{H}_2\text{S} + 358{,}2 \cdot \text{CH}_4 + 637{,}5 \cdot \text{C}_2\text{H}_6 + 912{,}7 \cdot \text{C}_3\text{H}_8 + 1186{,}5 \cdot \text{C}_4\text{H}_{10} + 1461{,}1 \cdot \text{C}_5\text{H}_{12}

где процентное содержание каждого компонента газовой смеси подставляется в объёмных долях (% об.\% \text{ об.}). Для типового природного газа уральского региона в системе заложена эталонная расчётная величина:

Qнр=37222,40 кДж/м3(≈8890 ккал/м3)Q_н^р = 37222{,}40 \text{ кДж/м}^3 \quad \left(\approx 8890 \text{ ккал/м}^3\right)

2.2. Потребность в кислороде и воздухе для горения ​

Теоретически необходимый объём кислорода для полного сгорания 1 м31 \text{ м}^3 газа:

VO20=0,01⋅[0,5⋅CO+0,5⋅H2+1,5⋅H2S+∑(m+n4)⋅CmHn−O2]V_{\text{O}_2}^0 = 0{,}01 \cdot \left[0{,}5 \cdot \text{CO} + 0{,}5 \cdot \text{H}_2 + 1{,}5 \cdot \text{H}_2\text{S} + \sum\left(m + \frac{n}{4}\right) \cdot \text{C}_m\text{H}_n - \text{O}_2\right]

Для базового состава природного газа расчётное значение составляет:

VO20=2,073879 м3/м3 газаV_{\text{O}_2}^0 = 2{,}073879 \text{ м}^3/\text{м}^3 \text{ газа}

Теоретический объём сухого воздуха при содержании кислорода 21% об.21\% \text{ об.}:

L0=VO200,21=4,762⋅VO20=9,871665 м3/м3 газаL_0 = \frac{V_{\text{O}_2}^0}{0{,}21} = 4{,}762 \cdot V_{\text{O}_2}^0 = 9{,}871665 \text{ м}^3/\text{м}^3 \text{ газа}

Действительный объём воздуха, подаваемого на горение, с учётом коэффициента избытка воздуха α\alpha:

Lα=α⋅L0L_\alpha = \alpha \cdot L_0

2.3. Объёмы продуктов полного сгорания ​

При теоретическом горении без избытка воздуха (α=1,0\alpha = 1{,}0) объёмы составляющих дымовых газов определяются стехиометрическими уравнениями:

VCO20=0,01⋅(CO2+CO+∑m⋅CmHn)=1,067409 м3/м3V_{\text{CO}_2}^0 = 0{,}01 \cdot \left(\text{CO}_2 + \text{CO} + \sum m \cdot \text{C}_m\text{H}_n\right) = 1{,}067409 \text{ м}^3/\text{м}^3
VH2O0=0,01⋅(H2O+H2+∑n2⋅CmHn)+0,00124⋅dв⋅L0=2,033991 м3/м3V_{\text{H}_2\text{O}}^0 = 0{,}01 \cdot \left(\text{H}_2\text{O} + \text{H}_2 + \sum \frac{n}{2} \cdot \text{C}_m\text{H}_n\right) + 0{,}00124 \cdot d_{\text{в}} \cdot L_0 = 2{,}033991 \text{ м}^3/\text{м}^3

где dвd_{\text{в}} — влагосодержание атмосферного воздуха (в расчётах принято dв=14 г/кг сухого воздухаd_{\text{в}} = 14 \text{ г/кг сухого воздуха}).

Объём азота:

VN20=0,01⋅N2+0,79⋅L0=0,03+3,762⋅VO20=7,83193 м3/м3V_{\text{N}_2}^0 = 0{,}01 \cdot \text{N}_2 + 0{,}79 \cdot L_0 = 0{,}03 + 3{,}762 \cdot V_{\text{O}_2}^0 = 7{,}83193 \text{ м}^3/\text{м}^3

Полный теоретический объём продуктов сгорания:

V0=VCO20+VH2O0+VN20=10,93333 м3/м3V_0 = V_{\text{CO}_2}^0 + V_{\text{H}_2\text{O}}^0 + V_{\text{N}_2}^0 = 10{,}93333 \text{ м}^3/\text{м}^3

2.4. Объёмы продуктов сгорания при α>1,0\alpha > 1{,}0 ​

При подаче избыточного воздуха появляются избыточный кислород и дополнительный азот:

VO2изб=(α−1)⋅VO20V_{\text{O}_2}^{\text{изб}} = (\alpha - 1) \cdot V_{\text{O}_2}^0
VN2α=VN20+(α−1)⋅3,762⋅VO20V_{\text{N}_2}^\alpha = V_{\text{N}_2}^0 + (\alpha - 1) \cdot 3{,}762 \cdot V_{\text{O}_2}^0
VH2Oα=VH2O0+(α−1)⋅0,00124⋅dв⋅L0V_{\text{H}_2\text{O}}^\alpha = V_{\text{H}_2\text{O}}^0 + (\alpha - 1) \cdot 0{,}00124 \cdot d_{\text{в}} \cdot L_0

Суммарный объём влажных дымовых газов:

Vα=VCO20+VH2Oα+VN2α+VO2избV_\alpha = V_{\text{CO}_2}^0 + V_{\text{H}_2\text{O}}^\alpha + V_{\text{N}_2}^\alpha + V_{\text{O}_2}^{\text{изб}}

3. Баланс массы влаги при сушке материалов ​

Исходная влажность материала задаётся в долях или процентах на влажную основу (WW):

W=mвлагиmвлажного материала⋅100%W = \frac{m_{\text{влаги}}}{m_{\text{влажного материала}}} \cdot 100\%

Влагосодержание материала на сухую основу (WcW^c):

Wc=W100−W⋅100%[кг влагикг сухой массы]W^c = \frac{W}{100 - W} \cdot 100\% \quad \left[\frac{\text{кг влаги}}{\text{кг сухой массы}}\right]

Масса влаги, удаляемой из материала в процессе сушки:

Wисп=Gсм⋅(W1c−W2c100)=Gсм⋅(W1100−W1−W2100−W2)[кг]W_{\text{исп}} = G_{\text{см}} \cdot \left(\frac{W_1^c - W_2^c}{100}\right) = G_{\text{см}} \cdot \left(\frac{W_1}{100 - W_1} - \frac{W_2}{100 - W_2}\right) \quad \left[\text{кг}\right]

где:

  • GсмG_{\text{см}} — масса абсолютно сухого материала (кг);
  • W1,W2W_1, W_2 — начальная и конечная влажность материала на влажную основу (%).

4. Термодинамические энтальпии дымовых газов и воздуха ​

Энтальпия продуктов сгорания Ig(t)I_g(t) при температуре газов tt (°C) рассчитывается через объёмные изобарные теплоёмкости компонентов:

Ig(t)=VCO2⋅(c⋅t)CO2+VH2O⋅(c⋅t)H2O+VN2⋅(c⋅t)N2+VO2⋅(c⋅t)O2[кДжм3 газа]I_g(t) = V_{\text{CO}_2} \cdot (c \cdot t)_{\text{CO}_2} + V_{\text{H}_2\text{O}} \cdot (c \cdot t)_{\text{H}_2\text{O}} + V_{\text{N}_2} \cdot (c \cdot t)_{\text{N}_2} + V_{\text{O}_2} \cdot (c \cdot t)_{\text{O}_2} \quad \left[\frac{\text{кДж}}{\text{м}^3 \text{ газа}}\right]

В алгоритмах системы для экспресс-расчёта в широком диапазоне температур (100…1200 ∘C100 \dots 1200 \text{ }^\circ\text{C}) используются непрерывные степенные аппроксимации термодинамических таблиц ВТИ:

Iгазов=(t6,7851)10,6906[кДжм3]I_{\text{газов}} = \left(\frac{t}{6{,}7851}\right)^{\frac{1}{0{,}6906}} \quad \left[\frac{\text{кДж}}{\text{м}^3}\right]
Iвозд=1,0286⋅t1,0501[кДжм3]I_{\text{возд}} = 1{,}0286 \cdot t^{1{,}0501} \quad \left[\frac{\text{кДж}}{\text{м}^3}\right]

5. Классификация статей теплового баланса ​

Все составляющие теплового баланса агрегатов приводятся к единой энергетической шкале (кДж за цикл или кДж/ч).

5.1. Статьи прихода тепла (∑Qприход\sum Q_{\text{приход}}) ​

  1. QтоплQ_{\text{топл}} — химическое тепло сгорания топлива:
    Qтопл=B⋅QнрQ_{\text{топл}} = B \cdot Q_н^р
  2. QвоздQ_{\text{возд}} — физическое тепло подогретого воздуха, подаваемого на горение:
    Qвозд=B⋅Lα⋅cвозд⋅tвоздQ_{\text{возд}} = B \cdot L_\alpha \cdot c_{\text{возд}} \cdot t_{\text{возд}}
  3. QэлQ_{\text{эл}} — электрическая энергия, подводимая к нагревателям (для электропечей):
    Qэл=Pэл⋅3600⋅ηтэнQ_{\text{эл}} = P_{\text{эл}} \cdot 3600 \cdot \eta_{\text{тэн}}

5.2. Статьи расхода тепла (∑Qрасход\sum Q_{\text{расход}}) ​

  1. Q1Q_1 — полезное тепло на сушку и нагрев материала:
    Q1=Qисп+Qнагрев.сух+Qнагрев.водыQ_1 = Q_{\text{исп}} + Q_{\text{нагрев.сух}} + Q_{\text{нагрев.воды}}
  2. Q2Q_2 — потери тепла с уходящими дымовыми газами:
    Q2=B⋅Vα⋅Ig(tух)Q_2 = B \cdot V_\alpha \cdot I_g(t_{\text{ух}})
  3. Q3Q_3 — потери от химической неполноты сгорания топлива:
    Q3=q3⋅B⋅Qнр(q3≈0,01…0,02)Q_3 = q_3 \cdot B \cdot Q_н^р \quad (q_3 \approx 0{,}01 \dots 0{,}02)
  4. Q5Q_5 — потери тепла в окружающую среду через футеровку:
    Q5=∑jKj⋅Fj⋅Δtj⋅τQ_5 = \sum_{j} K_j \cdot F_j \cdot \Delta t_j \cdot \tau
  5. Q6Q_6 — теплоаккумуляция футеровкой и конструктивными элементами печи при разогреве.

6. Критерии энергетической эффективности ​

6.1. Термический коэффициент полезного действия (КПД) ​

η=Q1∑Qприход⋅100%\eta = \frac{Q_1}{\sum Q_{\text{приход}}} \cdot 100\%

6.2. Удельный расход теплоты на испарение 1 кг влаги ​

qисп=∑QприходWисп[кДжкг влаги]q_{\text{исп}} = \frac{\sum Q_{\text{приход}}}{W_{\text{исп}}} \quad \left[\frac{\text{кДж}}{\text{кг влаги}}\right]

6.3. Удельный расход условного топлива ​

Условное топливо определяется эквивалентной теплотой сгорания Qусл=29308 кДж/кгQ_{\text{усл}} = 29308 \text{ кДж/кг} (7000 ккал/кг7000 \text{ ккал/кг}):

bуд=B⋅QнрQусл⋅Wисп⋅1000[кг у.т.т влаги]b_{\text{уд}} = \frac{B \cdot Q_н^р}{Q_{\text{усл}} \cdot W_{\text{исп}}} \cdot 1000 \quad \left[\frac{\text{кг у.т.}}{\text{т влаги}}\right]

Разработано в рамках проекта моделирования тепловых балансов и промышленного мониторинга