Теплоэнергетика, 2021, № 9, стр. 64-70

Определение температуры адиабатных поверхностей в помещениях с лучистым отоплением

А. О. Шацков ***

Донбасская национальная академия строительства и архитектуры
86123 г. Макеевка, Донецкая область, ул. Державина, д. 2, Украина

* E-mail: Reggie.88@mail.ru
** E-mail: aoshatskov@donnasa.ru

Поступила в редакцию 06.08.2020
После доработки 10.12.2020
Принята к публикации 23.12.2020

Полный текст (PDF)

Аннотация

Теплоснабжение с одновременным энергоресурсосбережением является одной из важнейших задач настоящего времени. Высокую эффективность показывают системы низкотемпературных электрических отопительных приборов инфракрасного излучения. Системы лучистого отопления имеют определенные преимущества перед традиционными системами отопления: возможность обогрева как всего помещения, так и отдельных зон; отсутствие массивных трубопроводов и воздуховодов и затрат на транспортировку теплоносителя; отсутствие циркуляции пыли и загрязняющих веществ в обогреваемом помещении; возможность снижения температуры воздуха в обогреваемом помещении на 1–3°С, что обусловливается особенностями процессов теплообмена при использовании радиационного отопления и допускается строительными нормами. Тепловой комфорт в помещениях с лучистыми отопительными приборами создается благодаря повышению радиационной температуры внутренних поверхностей ограждающих конструкций при снижении температуры воздуха в помещении. При этом необходимо добиться соблюдения первого условия комфортности. Для того чтобы определить, выполняются ли эти требования и соответствуют ли они действующим строительным нормам, нужно разработать такую методику расчета теплообмена в помещении при лучистом отоплении, которая позволит учесть также и все многообразие конфигураций отопительных приборов. Важнейшими параметрами будут размеры и температура поверхности приборов, их количество и расположение в обслуживаемом помещении. Также необходимо учитывать геометрические параметры помещений, материал и излучательную способность внутренней поверхности ограждающих конструкций. В статье рассматриваются отопительные приборы для жилых и общественных зданий. В связи с небольшой высотой таких помещений для выполнения второго условия комфортности применяются приборы с температурой поверхности до 100°С. Разработанная методика имеет некоторые допущения, которые дают погрешность не более 5%.

Ключевые слова: температура поверхности, лучистый теплообмен, тепловой поток, излучательная способность, коэффициент лучеиспускания, низкотемпературное лучистое отопление жилых и общественных зданий, угловой коэффициент излучения, коэффициент конвективного теплообмена

В специализированной литературе указаны преимущества использования техники инфракрасного излучения для отопления жилых и общественных зданий [13], при этом отмечается, что такие системы отопления позволяют не только создать комфортные условия для человека [4, 5], но и снизить потребление энергии для отопления [46]. Возможности и перспективы использования лучистого отопления широко изучаются зарубежным исследователями, что находит свое отражение в материалах, публикуемых в специализированных изданиях [711].

Тем не менее, несмотря на преимущества использования лучистого отопления, отопительные приборы инфракрасного излучения не находят широкого применения на практике, в частности, из-за отсутствия единой нормативной методики расчета и проектирования таких систем. Некоторые российские исследователи тщательно изучали вопросы лучистого теплообмена в помещениях с лучистыми отопительными приборами. При этом для расчета теплообмена применялись зональные методы вычислений [6, 1215]. Следует отметить, что подобные методы расчета превосходно показывают себя при проектировании систем газового лучистого отопления производственных зданий. В жилых и общественных зданиях не допускается установка таких отопительных приборов, так как температура их поверхности превышает 100°С.

В справочной литературе приведены методики расчета теплообмена в помещениях с панельно-лучистым отоплением [1618]. Однако такие методики либо разрабатываются для систем отопления конкретных производителей [16], либо в них принимаются значительные упрощения из-за отсутствия способов расчета средних диффузных угловых коэффициентов [17, 18].

При эксплуатации зданий основным параметром служит тепловой режим помещений, от которого зависят ощущение теплового комфорта людей, нормальное протекание производственных процессов, состояние и долговечность как конструкции здания, так и его оборудования. Тепловая обстановка в помещении определяется совместным действием таких факторов, как температура, подвижность и влажность воздуха, наличие струйных течений, распределение параметров воздуха в плане и по высоте помещения, а также радиационным излучением окружающих поверхностей, которое зависит от их температуры, конфигурации и радиационных свойств [1921]. Под действием конвективного и лучистого теплообмена и процессов массопереноса температуры воздуха и поверхностей в помещении взаимосвязаны и оказывают воздействие одна на другую [22, 23].

Для изучения формирования микроклимата, его динамики и способов воздействия на него нужно знать законы теплообмена в помещении [22, 23]. Температуры поверхностей в помещении неодинаковы. Обычно зимой и летом наружные ограждения и приборы систем отопления-охлаждения бывают более нагретыми или охлажденными по сравнению с внутренними стенами, которые имеют температуру, близкую к температуре воздуха в помещении. Между поверхностями происходит теплообмен излучением, подчиняющийся общим физическим закономерностям, пользоваться которыми в инженерных расчетах сложно. Лучистый теплообмен в помещении происходит в условиях ограниченного диапазона температур, определенных радиационных свойств поверхностей и их геометрического расположения [22, 23].

Теплообмен в помещении рассчитывают при установившемся режиме передачи тепла и тепловые потоки от греющей панели в помещение и из помещения наружу считают равными [17, 18].

ПОСТАНОВКА ЗАДАЧИ РАСЧЕТА СИСТЕМЫ ОТОПЛЕНИЯ

Задачей расчета является определение температуры поверхности внутренних ограждений с учетом как лучистого теплообмена между отопительной панелью и остальными взаимно параллельными и перпендикулярными поверхностями, так и конвективного теплообмена между воздухом и ограждениями. Знание этих температур позволяет проверить соблюдение комфортной температурной обстановки, уточнить теплопотери помещения и тепловую мощность отопительной панели.

В помещении обычной конструкции имеются поверхности шести видов, участвующие в теплообмене: наружные стены, окна, пол, потолок, внутренние стены, отопительная панель. Эти поверхности можно сгруппировать следующим образом: наружные ограждения (источник тепловых потерь), внутренние ограждения (адиабатные) и отопительная панель (источник теплопоступлений).

Следует также отметить, что жилое помещение представляет собой замкнутую систему из N поверхностей конечных размеров. Для упрощения решения задачи теплообмена принимаются следующие допущения:

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

каждая поверхность системы представляет собой серое тело;

излучение, испускаемое каждой поверхностью, диффузно распределяется в полупространстве;

поверхностная плотность потока эффективного излучения одинакова во всех точках каждой из поверхностей системы; при этом угловые коэффициенты излучения не зависят от поверхностного распределения потоков излучения и их значений.

ОПРЕДЕЛЕНИЕ ТЕМПЕРАТУРЫ АДИАБАТНЫХ ПОВЕРХНОСТЕЙ

Температура адиабатных поверхностей служит важным показателем работы системы лучистого отопления, так как температура воздуха внутри помещения при работе приборов инфракрасного излучения ниже, чем при работе конвективных отопительных приборов. Поэтому выполнения первого условия комфортности, определяющего такую область сочетаний температуры воздуха внутри помещения и радиационной температуры помещения, при которых человек, находясь в центре рабочей зоны, не испытывает ни перегрева, ни переохлаждения [20, 21], можно достичь повышением температуры только внутренних поверхностей ограждающих конструкций. Без методики расчета этой температуры невозможно определить, насколько эффективно работает система лучистого отопления.

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

(1)
${{Q}_{{\text{п}}}} = \sigma {{\varepsilon }_{{e\,{\text{п}}}}}{{\varepsilon }_{{e\,{\text{а}}}}}\left[ {{{{\left( {\frac{{T_{{\text{п}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{{\text{а}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right]{{\varphi }_{{{\text{п}} - {\text{а}}}}}{{F}_{{\text{а}}}},$
где Qп – тепловой поток от поверхности отопительной панели, Вт; $\sigma $ – постоянная Стефана – Больцмана, равная примерно 5.67 × 10–8 Вт/(м2 · К4); ${{\varepsilon }_{{e\,{\text{п}}}}}$ – излучательная способность излучающей поверхности отопительной панели; ${{\varepsilon }_{{e\,{\text{а}}}}}$ – излучательная способность внутренней (адиабатной) поверхности рассчитываемого внутреннего ограждения, численно равная излучательной способности внутренней поверхности наружного ограждения; $T_{{\text{п}}}^{{\text{л}}}$ – температура поверхности отопительной панели, К; $T_{{\text{а}}}^{{\text{л}}}$ – температура адиабатных поверхностей, К; φп–а – средний диффузный угловой коэффициент излучения с поверхности панели на адиабатные поверхности; Fa – суммарная площадь адиабатных поверхностей, м2; индекс “л” обозначает лучистую составляющую теплообмена в помещении.

Поскольку любое жилое помещение представляет собой замкнутую систему лученепрозрачных поверхностей, то при расчете теплообмена можно применять свойство замыкаемости угловых коэффициентов [2224], которое, в силу закона сохранения энергии для замкнутой системы, можно записать следующим образом:

$\sum\limits_{j = 1}^n {{{\varphi }_{{i - j}}}} = 1.$

Здесь ${{\varphi }_{{i\, - j}}}$ – угловой коэффициент излучения с i‑й поверхности на все j-е поверхности, участвующие в теплообмене.

Используя это свойство, расчет угловых коэффициентов можно значительно упростить и записать выражение (1) так:

${{Q}_{{\text{п}}}} = \sigma {{\varepsilon }_{{e\,{\text{п}}}}}{{\varepsilon }_{{e\,{\text{а}}}}}\left[ {{{{\left( {\frac{{T_{{\text{п}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{{\text{а}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right]\left( {1 - {{\varphi }_{{{\text{п}} - x}}}} \right){{F}_{{\text{а}}}},$
где ${{{{\varphi }}}_{{{\text{п}}--х}}}$ – угловой коэффициент излучения с поверхности отопительной панели на поверхность наружного ограждения.

Таким образом, для расчета лучистого теплообмена в помещении достаточно определить угловой коэффициент излучения на поверхность наружного ограждения, что значительно упрощает вычисления.

Температура адиабатных поверхностей несколько выше температуры внутренней поверхности наружной ограждающей конструкции, поэтому результирующий тепловой поток между ними ${{Q}_{x}}$ отличен от нуля и вычисляется по формуле

${{Q}_{x}} = \sigma \varepsilon _{{e\,{\text{а}}}}^{2}\left[ {{{{\left( {\frac{{T_{{\text{а}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{x}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right]{{F}_{x}},$
где Fх – площадь наружной ограждающей конструкции, м2; $T_{x}^{{\text{л}}}$ – температура внутренней поверхности наружной ограждающей конструкции, К.

Конвективный теплообмен между отопительной панелью и воздухом внутри помещения определяется по выражению

(2)
${{Q}_{{{\text{вн}}}}} = {{\alpha }_{{\text{к}}}}\left( {T_{{\text{п}}}^{{\text{к}}} - {{{{\theta }}}_{{{\text{вн}}}}}} \right){{F}_{{\text{п}}}},$
где Qвн – конвективный тепловой поток от поверхности отопительной панели; αк – коэффициент конвективного теплообмена, Вт/(м2 · К); $T_{{\text{п}}}^{{\text{к}}}$ – температура поверхности панели, К; θвн – температура воздуха внутри помещения, К; Fп – площадь панели, м2.

При температуре воздуха внутри помещения 20°С коэффициент конвективного теплообмена αк для вертикальных поверхностей равен [25]

${{\alpha }_{{\text{к}}}} = 1.16{{\left( {\Delta {{T}^{{\text{к}}}}} \right)}^{{\frac{1}{3}}}},$
где $\Delta {{T}^{{\text{к}}}}$ = $T_{{\text{п}}}^{{\text{к}}} - {{{{\theta }}}_{{{\text{вн}}}}}$ – разность температур поверхности панели и воздуха внутри помещения, К.

Тогда формулу (2) можно представить следующим образом:

${{Q}_{{{\text{вн}}}}} = 1.16{{\left( {\Delta {{T}^{{\text{к}}}}} \right)}^{{\frac{4}{3}}}}{{F}_{{\text{п}}}}.$

Потери тепла на нагрев приточного воздуха определяются по выражению (написание формулы приведено в соответствии с [25])

$\begin{gathered} 0.278{{G}_{{{\text{пр}}}}}\sigma \left( {{{{{\theta }}}_{{{\text{вн}}}}} - {{T}_{{\text{н}}}}} \right) = \sum\limits_{i = 1}^N {\int\limits_{{{A}_{i}}} {{{\alpha }_{{{\text{кd}}{{A}_{i}}}}}} } \left( {{{T}_{{{\text{d}}{{A}_{i}}}}} - {{T}_{{\text{н}}}}} \right){\text{d}}{{A}_{i}} = \\ = \sum\limits_{i = 1}^N {{{\alpha }_{{{{{\text{к}}}_{i}}}}}\left( {{{T}_{{\text{а}}}} - {{{{\theta }}}_{{{\text{вн}}}}}} \right){{F}_{{{\text{a}}\,i}}}} + 1.16{{\left( {T_{{\text{п}}}^{{\text{к}}} - {{{{\theta }}}_{{{\text{вн}}}}}} \right)}^{{\frac{4}{3}}}}{{F}_{{\text{п}}}}, \\ \end{gathered} $
где Gпр – массовый расход приточного воздуха, кг/с; Tн – температура наружного воздуха, К; Ai – площадь i-й поверхности помещения, м2; ${{T}_{{{\text{d}}{{A}_{i}}}}}$ – температура i-й поверхности помещения, К; ${{\alpha }_{{{{{\text{к}}}_{i}}}}}$ – коэффициент конвективного теплообмена между i-й адиабатной поверхностью помещения и воздухом внутри помещения, Вт/(м2 · К); ${{T}_{{\text{a}}}}$ – температура адиабатных поверхностей помещения, К; ${{F}_{{{\text{a}}\;i}}}$ – площадь i-й адиабатной поверхности, м2.

Коэффициент конвективного теплообмена вычисляется по формулам из [25] для различных поверхностей:

вертикальных

${{\alpha }_{{\text{к}}}} = 1.39{{\left( {\frac{{\Delta {{T}^{{\text{к}}}}}}{l}} \right)}^{{\frac{1}{4}}}};$

нагретых горизонтальных, обращенных вниз,

(3)
${{\alpha }_{{\text{к}}}} = 1.89{{\left( {\frac{{\Delta {{T}^{{\text{к}}}}}}{l}} \right)}^{{\frac{1}{4}}}};$

нагретых горизонтальных, обращенных вверх,

(4)
${{\alpha }_{{\text{к}}}} = 0.971{{\left( {\frac{{\Delta {{T}^{{\text{к}}}}}}{l}} \right)}^{{\frac{1}{4}}}},$
где l – определяющий размер поверхности в направлении потока воздуха, м.

Геометрические размеры пола и потолка в большинстве жилых и общественных зданий одинаковы. При установке отопительных приборов на стенах формулы (3) и (4) можно представить в виде

${{\alpha }_{{\text{к}}}} = 2.861{{\left( {\frac{{\Delta {{T}^{{\text{к}}}}}}{l}} \right)}^{{\frac{1}{4}}}}.$

Тепловой баланс для адиабатных поверхностей определяется по формуле

${{Q}_{{\text{п}}}} - {{Q}_{х}} - {{Q}_{{\text{а}}}} = 0$
(где ${{Q}_{{\text{а}}}}$ – суммарный конвективный тепловой поток от всех адиабатных поверхностей помещения, Вт) или

(5)
$\begin{gathered} \sigma {{\varepsilon }_{{e\,{\text{п}}}}}{{\varepsilon }_{{e\,{\text{а}}}}}\left[ {{{{\left( {\frac{{T_{{\text{п}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{{\text{а}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right]\left( {1 - {{\varphi }_{{{\text{п}} - x}}}} \right){{F}_{{\text{п}}}} - \sigma \varepsilon _{{e\,{\text{а}}}}^{2} \times \\ \times \,\,\left[ {{{{\left( {\frac{{T_{{\text{а}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{x}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right]{{F}_{{\text{а}}}} - \sum {{{\alpha }_{{{{{\text{к}}}_{i}}}}}\left( {T_{{\text{а}}}^{{\text{к}}} - {{{{\theta }}}_{{{\text{вн}}}}}} \right)} {{F}_{i}} = 0. \\ \end{gathered} $

В выражении (5) входными являются значения следующих параметров:

отопительной панели ${{\varepsilon }_{{e\,{\text{п}}}}},$ $T_{{\text{п}}}^{{\text{л}}},$ Fп, которые определяются по паспорту прибора;

наружного ограждения Fх, ${{\varepsilon }_{{e\,{\text{а}}}}},$ а также температуры его внутренней поверхности $T_{x}^{{\text{л}}},$ численно равной температуре воздуха внутри помещения [26];

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

Перед началом расчета необходимо вычислить коэффициент конвективного теплообмена αк для всех типов поверхностей, а также средний диффузный угловой коэффициент излучения с поверхности панели на наружное ограждение. После этого выражение (5) можно представить как уравнение 4-й степени для определения температуры адиабатных поверхностей.

РАСЧЕТ ТЕМПЕРАТУРЫ С УЧЕТОМ ТЕПЛООБМЕНА С ЗАПОЛНЕНИЕМ СВЕТОВОГО ПРОЕМА

В некоторых жилых и общественных зданиях значительную часть наружного ограждения составляют заполнения световых проемов. Так как зачастую возле окон создается зона теплового дискомфорта даже при относительно высоких показателях термического сопротивления заполнения проема, как показано на рис. 1, то иногда может потребоваться расчет температурного режима помещения с учетом теплообмена через окна. Тогда тепловые потоки через наружное ограждение следует разделять на теплообмен через заполнение светового проема и через остальную часть наружной стены. Возникающие при этом тепловые потоки схематически изображены на рис. 2. При установке отопительной панели на наружном ограждении теплообмен излучением между панелью и наружным ограждением будет равен нулю (поз. 1 и 2 на указанной схеме).

Рис. 1.

Зона теплового дискомфорта. а – радиатор; б – конвектор; в – панель в перегородке; г – подоконная панель; д – границы зон комфорта в помещении: при размещении радиаторов (І) и конвекторов (ІІ) под окнами, нагревательных элементов по периметру внутренних стен (ІІІ) и в междуэтажных перекрытиях (ІV), в стеновой панели (V); А – зона комфорта; Б – зона дискомфорта; С – стена помещения в плане с заливкой

Рис. 2.

Граф тепловых потоков в помещении. I – отопительная панель; II – внутренние ограждения; III – заполнение светового проема; IV – наружное ограждение без учета заполнения светового проема; V – воздух внутри помещения; VI – наружный воздух; 1, 2, 3 – теплообмен излучением между отопительной панелью и остальными поверхностями; 4 – конвективный теплообмен между панелью и воздухом внутри помещения; 5, 6 – теплообмен излучением между внутренними и наружными ограждениями; 7 – конвективный теплообмен между внутренними ограждениями и воздухом в помещении; 8, 9 – конвективный теплообмен между наружными ограждениями и воздухом внутри помещения; 10, 11 – теплопередача через наружное ограждение

Расчет температуры наружного ограждения значительно усложняется, так как, по сравнению с выражением (5), в линейном уравнении появляются два неизвестных – температуры поверхности окна и остальной части ограждения. При этом температура окна будет ниже температуры остальной части ограждения, следовательно, необходимо учитывать лучистый теплообмен между окном и внутренними стенами. Поэтому для расчета требуется составить систему из трех уравнений, в которой первое уравнение – это тепловой баланс оконного проема, второе – тепловой баланс адиабатных поверхностей, третье – тепловой баланс внутренней поверхности наружного ограждения без учета заполнения светового проема.

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

(6)
$\left. \begin{gathered} \sigma {{\varepsilon }_{{e\,{\text{п}}}}}{{\varepsilon }_{{e\,{\text{о}}}}}\left[ {{{{\left( {\frac{{T_{{\text{п}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{{\text{о}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right]{{\varphi }_{{{\text{п}} - {\text{о}}}}}{{F}_{{\text{п}}}} + \sigma {{\varepsilon }_{{e\,{\text{о}}}}}{{\varepsilon }_{{e\,{\text{а}}}}} \times \hfill \\ \times \,\,\left[ {{{{\left( {\frac{{T_{{\text{а}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{{\text{о}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right]\sum {{{\varphi }_{{i - {\text{о}}}}}{{F}_{{{\text{a}}\,i}}}} - \hfill \\ - \,\,{{k}_{{\text{о}}}}{{F}_{{\text{о}}}}\left( {T_{{\text{о}}}^{{\text{л}}} - T_{{\text{н}}}^{{\text{л}}}} \right) + {{F}_{{\text{о}}}}{{\alpha }_{{\text{к}}}}\left( {T_{{{\text{вн}}}}^{{\text{л}}} - T_{{\text{о}}}^{{\text{л}}}} \right) = 0; \hfill \\ \sigma {{\varepsilon }_{{e\,{\text{п}}}}}{{\varepsilon }_{{e\,{\text{а}}}}}\left[ {{{{\left( {\frac{{T_{{\text{п}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{{\text{а}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right] \times \hfill \\ \times \,\,\left[ {1 - \left( {{{\varphi }_{{{\text{п}} - х}}} + {{\varphi }_{{{\text{п--о}}}}}} \right)} \right]{{F}_{{\text{п}}}} - \sigma {{\varepsilon }_{{e\,{\text{а}}}}}{{\varepsilon }_{{e\,{\text{о}}}}} \times \hfill \\ \times \,\,\left[ {{{{\left( {\frac{{T_{{\text{а}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{{\text{о}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right]\sum {{{\varphi }_{{i - {\text{о}}}}}{{F}_{{{\text{a}}\,i}}}} - \sigma \varepsilon _{{e\,{\text{а}}}}^{2} \times \hfill \\ \times \,\,\left[ {{{{\left( {\frac{{T_{{\text{а}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{x}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right]\sum {{{\varphi }_{{j - х}}}{{F}_{{{\text{a}}\,i}}}} - \hfill \\ - \,\,{{F}_{{\text{а}}}}{{\alpha }_{{\text{к}}}}\left( {T_{{\text{а}}}^{{\text{л}}} - T_{{{\text{вн}}}}^{{\text{л}}}} \right) = 0; \hfill \\ \sigma {{\varepsilon }_{{e\,{\text{п}}}}}{{\varepsilon }_{{e\,{\text{а}}}}}\left[ {{{{\left( {\frac{{T_{{\text{п}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{x}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right]{{\varphi }_{{{\text{п}} - х}}}{{F}_{{\text{п}}}} + \sigma \varepsilon _{{e\,{\text{а}}}}^{2} \times \hfill \\ \times \,\,\left[ {{{{\left( {\frac{{T_{{\text{а}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{x}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right]\sum {{{\varphi }_{{i - {\text{о}}}}}{{F}_{{{\text{a}}\,i}}}} - \hfill \\ - \,\,{{k}_{x}}F_{x}^{'}\left( {T_{x}^{{\text{л}}} - T_{{\text{н}}}^{{\text{л}}}} \right) = 0, \hfill \\ \end{gathered} \right\}$
где $T_{{\text{о}}}^{{\text{л}}}$ – температура поверхности окна, К; Fо – площадь окна, м2; $F_{x}^{'}$ – площадь поверхности наружного ограждения без учета площади заполнения светового проема, м2; φi–o – угловой коэффициент излучения с i-й поверхности внутреннего ограждения на окно; ${{F}_{{{\text{a}}\,i}}}$ – площадь i-й поверхности внутреннего ограждения, м2; φjx – угловой коэффициент излучения с j-й поверхности внутреннего ограждения на внутреннюю поверхность наружного ограждения; ${{F}_{{{\text{a}}\,j}}}$ – площадь j-й поверхности внутреннего ограждения, м2; φп–о – угловой коэффициент излучения с поверхности отопительной панели на заполнение светового проема; φп‒х – угловой коэффициент излучения с поверхности отопительной панели на внутреннюю поверхность наружного ограждения; ${{\varepsilon }_{{e\,{\text{о}}}}}$ – излучательная способность поверхности заполнения светового проема; kо – коэффициент теплопередачи через заполнение светового проема, Вт/(м2 · К); kх – коэффициент теплопередачи через наружное ограждение, Вт/(м2 · К).

При установке отопительной панели на наружном ограждении имеем

(7)
$\left. \begin{gathered} \sigma {{\varepsilon }_{{e\,{\text{о}}}}}{{\varepsilon }_{{e\,{\text{а}}}}}\left[ {{{{\left( {\frac{{T_{{\text{а}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{{\text{о}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right] \times \hfill \\ \times \,\,\sum {{{\varphi }_{i}}{{F}_{i}} - {{k}_{{\text{о}}}}{{F}_{{\text{о}}}}\left( {T_{{\text{о}}}^{{\text{л}}} - T_{{\text{н}}}^{{\text{л}}}} \right)} + \hfill \\ + \,\,{{F}_{{\text{о}}}}{{\alpha }_{{\text{о}}}}\left( {T_{{{\text{вн}}}}^{{\text{л}}} - T_{{\text{о}}}^{{\text{л}}}} \right) = 0; \hfill \\ \sigma {{\varepsilon }_{{e\,{\text{п}}}}}{{\varepsilon }_{{e\,{\text{а}}}}}\left[ {{{{\left( {\frac{{T_{{\text{п}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{{\text{а}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right]{{F}_{{\text{п}}}} - \hfill \\ - \,\,\sigma {{\varepsilon }_{{e\,{\text{о}}}}}{{\varepsilon }_{{e\,{\text{а}}}}}\left[ {{{{\left( {\frac{{T_{{\text{а}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{{\text{о}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right]\sum {{{\varphi }_{{{\text{a}} - i}}}{{F}_{{{\text{a}}\,i}}}} - \hfill \\ - \,\,\sigma \varepsilon _{{e\,{\text{а}}}}^{2}\left[ {{{{\left( {\frac{{T_{{\text{а}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{x}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right]\sum {{{\varphi }_{{{\text{a}} - j}}}{{F}_{{{\text{a}}\,j}}}} - \hfill \\ - \,\,{{F}_{{{\text{вн}}}}}{{\alpha }_{{\text{к}}}}\left( {T_{{\text{а}}}^{{\text{л}}} - T_{{{\text{вн}}}}^{{\text{л}}}} \right) = 0; \hfill \\ \sigma \varepsilon _{{e\,{\text{а}}}}^{2}\left[ {{{{\left( {\frac{{T_{{\text{а}}}^{{\text{л}}}}}{{100}}} \right)}}^{4}} - {{{\left( {\frac{{T_{x}^{{\text{л}}}}}{{100}}} \right)}}^{4}}} \right] \times \hfill \\ \times \,\,\sum {{{\varphi }_{{a - j}}}{{F}_{{{\text{a}}\,j}}} - {{k}_{x}}} F_{x}^{'}\left( {T_{{{\text{вн}}}}^{{\text{л}}} - T_{{\text{н}}}^{{\text{л}}}} \right) = 0. \hfill \\ \end{gathered} \right\}$

Данные системы уравнений применимы для подавляющего большинства возможных конфигураций взаимного размещения отопительных панелей, наружного ограждения и заполнения светового проема. Для помещений с уникальными архитектурно-планировочными решениями может потребоваться корректировка представленных выражений.

ВЫВОДЫ

1. Несмотря на кажущуюся громоздкость систем уравнений, представленная методика определения температуры адиабатных поверхностей удобна в использовании, поскольку практически все переменные в уравнениях являются исходными параметрами при проектировании систем отопления.

2. При выведении полученных зависимостей использовалась температура внутренней поверхности наружного ограждения, являющаяся вместе с температурой воздуха нормируемой величиной. Тем не менее, следует разработать методику определения фактических значений этих температур. Тогда с помощью всех зависимостей можно будет проводить тщательный анализ эффективности работы систем низкотемпературного лучистого отопления.

Список литературы

  1. Краснощёк Б.Н. Прирученное солнце: водяные потолочные панели лучистого отопления “Теплопанель” // Энергосбережение. 2017. № 1. С. 20–22.

  2. Волошин С.А. Лучистые климатические системы. Особенности проектирования // АВОК: Вентиляция, отопление, кондиционирование воздуха, теплоснабжение и строительная теплофизика. 2013. № 2. С. 88–92. https://www.abok.ru/for_spec/articles. php?nid=5473

  3. Шумилов Р.Н., Толстова Ю.И., Поммер А.А. Лучистое отопление – мифы и реальность // Сантехника, отопление, кондиционирование. 2006. № 1. С. 56–58.

  4. Мордухович И.И., Гольченко Н.Ф., Ткаченко В.Б. Лучистое отопление // Вестник Сумского нац. аграр. ун-та. 2015. № 10. С. 86–100.

  5. Чередник А.Д., Редько А.А. Эффективность инфракрасных водяных отопительных панелей // MOTROL. Commission of Motorization and Energetics in Agriculture. 2014. № 6. P. 69–76.

  6. Математическое моделирование сопряженного теплопереноса в системе с радиационным источником нагрева / Н.И. Куриленко, В.И. Максимов, Г.Я. Мамонтов, Т.А. Нагорнова // Науч.-техн. ведомости СПбГПУ. Наука и образование. 2012. № 2. С. 48–53.

  7. Sattari S., Farhanieh B. A parametric study on radiant floor heating system performance // Renewable Energy. 2006. V. 31. Is. 10. P. 1617–1626. https://doi.org/10.1016/j.renene.2005.09.009

  8. Olesen B.W. Radiant floor heating in theory and practice // ASHRAE J. 2002. V. 44. № 7. P. 19–24.

  9. Hu R., Niu J.L. A review of the application of radiant cooling & heating systems in Mainland China // Energy Build. 2012. V. 52. P. 11–19. http://www.doi.org/10.1016/j.enbuild.2012.05.030

  10. History of radiant heating & cooling systems. Part 1 / R. Bean, B.W. Olesen, K.W. Kim, W. Kwang // ASHRAE J. 2010. V. 52. № 1. P. 40–47.

  11. Rheea K.N., Kim K.W. A 50 year review of basic and applied research in radiant heating and cooling systems for the built environment // Build. Environ. 2015. V. 91. P. 166–190. https://doi.org/10.1016/j.buildenv.2015.03.040

  12. Бухмиров В.В., Крупенников С.А., Солнышкова Ю.С. Модификация зонального метода для решения задач радиационного теплообмена: основные положения // Вестник ИГЭУ. 2009. № 2. С. 61–63.

  13. Бухмиров В.В., Крупенников С.А., Солнышкова Ю.С. Алгоритм расчета систем лучистого отопления помещений // Вестник ИГЭУ. 2010. № 4. С. 23–25.

  14. Бухмиров В.В., Солнышкова Ю.С., Пророкова М.В. Экспериментальное исследование системы отопления с инфракрасными излучателями // Вестник ИГЭУ. 2011. № 3. С. 12–16.

  15. Бухмиров В.В., Солнышкова Ю.С. Зональные методы расчета радиационного и сложного теплообмена. Иваново: ИГЭУ, 2012.

  16. Р НП “АВОК” 4.1.6-2009. Рекомендации АВОК. Системы отопления с потолочными подвесными излучающими панелями. Утв. и введены в действие приказом президента НП “АВОК” от 28.09.2009. Дата введения 02.11.2009. М.: АВОК-Пресс, 2009.

  17. Сканави А.Н., Махов Л.М. Отопление. М.: АСВ, 2002.

  18. Богословский В.Н., Сканави А.Н. Отопление. М.: Стройиздат, 1991.

  19. Банхиди Л. Тепловой микроклимат помещений: расчет комфортных параметров по теплоощущениям человека: пер. с венг. В.М. Беляева; под ред. В.И. Прохорова и А.Л. Наумова. М.: Стройиздат, 1981.

  20. Fanger P.O. Thermal comfort. N.Y.: McGrow Hill Book Co., 1973.

  21. Богословский В.Н. Строительная теплофизика. 2‑е изд., перераб. и доп. М.: Высшая школа, 1982.

  22. Modest M.F. Radiative heat transfer. 3rd ed. San Diego: Academic Press, 2013.

  23. Lienhard IV J.H., Lienhard V J.H. A heat transfer textbook. 5th ed. Cambridge, Massachusetts, USA: Phlogiston Press, 2003.

  24. Табунщиков Ю.А., Бродач М.М. Математическое моделирование и оптимизация тепловой эффективности зданий. М.: АВОК-Пресс, 2002.

  25. Михеев М.А. Основы теплопередачи. М.: Госэнергоиздат, 1956.

  26. Мачкаши А., Банхиди Л. Лучистое отопление: пер. с венг. В.М. Беляева; под ред. В.Н. Богословского и Л.М. Махова. М.: Стройиздат, 1985.

Дополнительные материалы отсутствуют.