Вестник УКЦ АПИК: Мощность тепловой завесы с жидким теплоносителем в общем случае

0
337

Алексей Пухов
Алексей Пухов, технический директор компании «ТРОПИК»

В данной статье, являющейся изложением доклада, прочитанного в рамках выставки «МИР КЛИМАТА-2013», предлагается уточненный метод определения тепловой мощности воздушных завес с жидким теплоносителем. Этот эмпирический подход существенно повышает точность процедуры во всем рабочем диапазоне воздушных завес и является обобщением метода расчета, предложенного в статье «Методика оценочного расчета тепловых параметров водяных завес», опубликованной в журнале «МИР КЛИМАТА» № 76.

Предположим, что мощность, приходящаяся на конвекцию и излучение с поверхности корпуса работающей завесы, пренебрежимо мала по сравнению с мощностью, приходящейся на нагрев воздушной струи (W C + W R << W).

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

W=4,2(T г – T x )G               (1)

где Т г — температура поступающей воды, Т х — обратной, G — массовый расход воды, а множитель 4,2 кДж/кг °С — ее теплоемкость. Собственно, формула (1) выражает мощность не только воздушной завесы, но и любой теплообменной системы. Однако ниже будет показано, насколько может вводить в заблуждение механическое использование данного выражения, хотя само оно и не может быть подвергнуто сомнению. Выражение (1) к тому же не позволяет производить пересчеты мощности W при изменении температурных параметров даже в небольших окрестностях около некоторых условных начальных значений.

Рассмотрим воздушную завесу как черный ящик, входными параметрами которого являются G, Т г , Т х и Т 00 — воздушная температура окружения завесы), при этом Т г > Т х > Т 0. Единственным выходным параметром «черного ящика» является W.

Формулу (1) можно переписать в виде:

W=4,2(T г – T x )G   (при Т 0)                (2)

Мощность тепловой завесы с жидким теплоносителем в общем случае
Рис. 1

Выражение (2) означает, что мощность W производится данным «черным ящиком» при соответствующих входных параметрах Т г , Т х , G и Т 0 . Будем считать воздушный поток завесы постоянным. Рассмотрим теперь, как данный ящик будет отвечать на поочередное изменение всех входных параметров.

Для начала, зафиксировав Т г и Т 0 , попробуем уменьшить значение Т х , что приведет к увеличению разности (Т г – Т х ). Однако при фиксированных Т г и Т 0 уменьшить Т х можно, лишь уменьшив расход G. Опыт показывает, что, несмотря на то что разность (Т г – Т х ) возрастет, G уменьшится еще более значительно для любых воздушных завес при любых начальных условиях и мощность W снизится.

Теперь зафиксируем Т г , Т х и уменьшим окружающую температуру до значения Т 0 ‘< Т 0 . Что произойдет с формулой (2)? Она перестанет работать в представленном выше виде и сможет быть записана как:

W’=4,2(T г – T x )G’   (при Т 0 ‘)                (3)

где новые значения расхода G’ и мощности W’ будут больше своих начальных значений. И в этом случае очень мало практической пользы от выражений (1) или (2), так как нет никакой ясности в поведении функции G’ (Т 0 ‘).

Наконец, зафиксировав значения Т х и Т 0 , будем увеличивать Т г . Интуитивно мы будем ожидать увеличения мощности W в (2), однако опыт показывает, что она будет оставаться практически постоянной для всех воздушных завес, что явно противоречит здравому смыслу. Лишь при большом увеличении температуры прямой воды будет наблюдаться незначительное увеличение тепловой мощности W.

Приведенные опытные данные свидетельствуют, что наш «черный ящик» накладывает некоторые связи на свои входные параметры, и без понимания вида данных связей невозможно представить явную зависимость мощности от этих параметров. Формула (2) позволяет лишь снимать характеристики воздушных завес при некоторых условиях и предлагать их в виде таблиц со стандартными наборами температур.

Мощность тепловой завесы с жидким теплоносителем в общем случае
Рис. 2

Следует отметить еще одно интересное свойство этого «черного ящика». Оно является прямым следствием линейности уравнений теплопроводности (и вынужденной конвекции [2]) относительно температуры и заключается в существовании подобных температурных полей при фиксированном расходе. Поясним данные рассуждения. Пусть в (2) G постоянен. Изменим любые две температуры из Т г , Т х , Т 0 произвольным образом, но так, чтобы вода сохраняла свое агрегатное состояние. Пусть, например, Т г ‘ и Т 0 ‘ — это некоторые новые значения двух температур. Тогда при значении третьей температуры:

T x ‘=T 0 ‘ + k(T x -T 0 ), где k=(T г ‘ – T 0 ‘)/(T г -T 0)                (4)

очевидно, будет выполняться пропорциональность:

(T x ‘ -T 0 ‘)/(T x -T 0 )=(T г ‘ – T 0 ‘)/(T г – T 0 )=(T г ‘ – T x ‘)/(T г – T x )=k                (5)

и формулу (2) можно представить в виде:

W’=4,2k(T г – T x )G’ (при Т 0 ‘)                (6)

то есть мощность W’ увеличится в k раз по сравнению с мощностью W начального состояния (2). При граничных условиях, определенных (5), температурные градиенты в каждой точке теплообменника увеличатся в k раз, и равенство (6) не вызывает сомнений. Это свойство температурного подобия может быть использовано как при измерениях характеристик воздушных завес, так и в качестве необходимого проверочного критерия для любых эмпирических и теоретических соотношений. Заметим, что в процессах теплопередачи, где существенна роль свободной конвекции или излучения, например в радиаторах отопления, использовать подобие в виде (5) невозможно по причине нелинейности данных процессов.

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

Некоторое время назад специалистами компании ЗАО «Антарес ПРО» была предложена закономерность для мощности тепловой завесы в виде:

W’≈ 4,2k(T x – T 0 )C k                (7)

основанная на результатах измерений, где С к — некоторая определяющая тепловую завесу константа. Несмотря на свой приближенный характер, формула (7) в отличие от (2) производит пересчет тепловой мощности при изменении одного или нескольких входных параметров нашего «черного ящика». Действительно, за скобками выражения (7) остался расход теплоносителя G, который может служить свободным параметром при изменениях Т 0 , Т х и Т г , (само Т г явно в выражение (7) не входит, показывая, что зависимость W (Т г ) слабая). Заметим, что (7) не противоречит критериям температурного подобия (5) и (6). Данное условие не может быть выполнено наперед для «любого» заданного закона, например, гипотетические представления W ≈ 4,2 (Т х – Т 0 ) С к / (1+Т 0г ) или W ≈ 4,2С к / (Т 0 – Т х ) не удовлетворяют (5) и (6) и вследствие этого опровергаются априори до этапа экспериментального подтверждения.

Из совместного решения (7) и (2) можно получить искомую связь между входными параметрами «черного ящика» в виде:

G≈C k ((T x – T 0 )/(T г – T x ))                (8)

Данное выражение позволяет объяснить неочевидный на первый взгляд характер поведения мощности W в (2) при рассмотренных выше изменениях температур. Однако выражение (7), являясь произведением С к и некоторого температурного коэффициента, вводит еще и понятие относительного расхода G/С к . Один и тот же расход теплоносителя для некоторой воздушной завесы может быть большим, а для завесы со значительно большим С к (или, другими словами, с большей расчетной мощностью) — малым.

Что касается применения данного рассмотрения через С к для конкретных задач, то оно было использовано в упрощенной методике верификации воздушных завес на базе АПИК [3] (хотя само обозначение С к явно в методике и не использовалось). Напомним, что С к — это константа лишь условно, она некоторым образом меняется при изменении относительного расхода. Следовательно, чем меньше расход при измерениях отличается от расхода при расчетных параметрах, тем выше будет точность пересчета через С к .

Стремясь к значительно более высокой точности в новой редакции методики верификации [4], мы стали получать температуры, еще более близкие к расчетным, a незначительный пересчет производить по законам температурного подобия (5) и (6). И только минимальное температурное расхождение величиной обычно менее 1 °С, обусловленное невозможностью точной установки температурного подобия, пересчитывается согласно закономерности (7), что обеспечивает очень высокую точность итоговых расчетов.

Однако если для лабораторного определения параметров использование приблизительных закономерностей необходимо сводить к минимуму, то за пределами научных лабораторий (7) будет давать удовлетворительную точность в областях больших и средних расходов. Несколько худшую точность данная формула дает для малых расходов теплоносителя, и совершенно неприемлемой она становится для очень малых расходов, если С к для расчетов мощности был измерен в области больших расходов. Если же С к изначально измерить для малого расхода (смысл определений степеней расхода станет ясным непосредственно из приведенных ниже графиков), то эти данные могут давать погрешность для расчетов мощностей при больших расходах.

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

Предложим для этой цели уточненную эмпирическую формулу:

W=4,2k(0,12T г + 0,88T x – T 0 )C       (9)

Данная формула, так же как и (7), позволяет производить пересчеты мощности при изменении температурных условий (изменяется расход, который явно в (7) не входит). Совместное решение (9) и (2) дает для расхода зависимость:

G=(C(0,12T г + 0,88T x – T 0 ))/(T г – T x)                (10)

Заметим, что данные выражения не противоречат и реализующемуся при постоянных расходах G критерию состояний температурного подобия (5) и (6).

Для определения С необходимо будет сделать всего одно измерение мощности при некоторых условиях. Заметим, что С ≠ С к . Если расход измеряется в кг/с, то:

C=(T г – T x)/(0,12T г + 0,88T x – T 0 )[кг/с]               (11)

Мощность тепловой завесы с жидким теплоносителем в общем случае
Рис. 3
Мощность тепловой завесы с жидким теплоносителем в общем случае
Рис. 4
Мощность тепловой завесы с жидким теплоносителем в общем случае
Рис. 5
Мощность тепловой завесы с жидким теплоносителем в общем случае
Рис. 6

Покажем теперь на графиках сравнительное поведение С к и С для воздушной завесы Тропик Т 109W [5] и трех воздушных завес Frico различных серий [6]. Для этого решим уравнения (7) и (9) относительно Ск и С соответственно для всех представленных в [5] и [6] данных и отобразим это в виде графиков (рис. 3–6) зависимостей данных величин от некоторого расходного параметра 0 < (Т х – Т 0 ) / (Т г – Т 0 ) < 1.

Из приведенных графиков следует, что С является константой уже во всех рассматриваемых областях, что и служит обоснованием формулы (9) для областей от больших до очень малых расходов. Физический смысл величины С — сколько килограммов теплоносителя за одну секунду данная воздушная завеса может «обработать» в области средних расходов, то есть при заметном падении величины Т х по сравнению с начальной температурой Т г . (Точнее, С = G при условиях Т г = 90 °С, Т х = 50,11 °С, Т 0 = 15 °С.) При оптимальном использовании воздушных завес расход теплоносителя воздушной завесы составляет величину более двух коэффициентов С для данной завесы.

Итак, получена закономерность для мощности и расхода воздушных завес во всем диапазоне изменения рабочих температур. Выражение (9) подтверждает слабую зависимость мощности воздушной завесы от температуры Т г . Покажем, как производить пересчет мощности завесы, если ее характеристики получены для стандартных условий: Т г = 95 °С, Т х = 70 °С, Т 0 = 15 °С. Как будет вести себя воздушная завеса при иных условиях работы? Когда можно и когда нельзя применять упрощенный подход через С к , предложенный в [1]? Для ответа на эти вопросы рассмотрим три случая возможного задания условий:

  1. Задается некоторый набор температур (Т г ‘, Т х ‘, Т 0 ‘), который однозначно определяет расход G’ и, следовательно, мощность W’. Рассчитаем сначала по (8) и (10) С к и С соответственно. Из (7), (8), (9) и (10) следует, что если мы хотим остаться в пределах погрешности 5 %, то для любой воздушной завесы пересчет мощности через (7) можно производить для случаев, когда Т х ‘ ближе к Т г ‘, чем к Т 0 ‘ или: Т х ‘ > (Т г ‘+Т 0 ‘) /2. Это же условие можно выразить через параметр абсолютного расхода, приблизительно как G’ > С к . Для расходов меньших, чем С к пересчет через С к даст заниженное выражение для мощности.
  2. Пусть теперь задан набор (Т г ‘, Т 0 ‘, G’), который определяет температуру обратной воды Т х ‘ и мощность W’. В этом случае область применимости расчетов через С к значительно расширяется: точности с погрешностью менее 5 % можно добиться и при значительно меньших расходах теплоносителя, и всю область применимости формулы (7) можно определить как G > 2/3Ск. Кажущееся противоречие с условием применимости для случая 1 легко разрешается, если вспомнить, что для всех воздушных завес при изменении величин Т х и G при фиксированных Т г и Т 0 , расход G изменяется сильнее, чем Т г – Т х . Заданный жестко расход G’ компенсирует возникающее расхождения слабо меняющегося параметра Т х ‘, что и расширяет границы применимости (7).
  3. Пусть теперь задан набор (Т х ‘, Т 0 ‘, G’), подразумевающий поиск температуры прямой воды, чтобы обратная была близка к (или не превосходила) Т х ‘. Учитывая условие, что Т г явно не входит в выражение мощности через С к , интуиция подсказывает, что в этом случае расхождения точной и приближенной зависимости обнаружатся при большем расходе, чем в первых двух случаях. Это действительно так: чтобы остаться в пределах 5 % погрешности, расчет через С к можно применять только при G > 1,15С к . В остальных случаях потребуется применение формулы (9).

Рассмотрим конкретные примеры расчетов завес. Пусть некоторая воздушная завеса продемонстрировала для данных температурных условий мощность W = 27,5кВт при расходе G = 0,262кг/с [7]. Можно рассчитать характеризующие ее константы С = 0,113кг/с и С к = 0,119кг/с.

  1. Пусть для рассмотренной завесы необходимо найти мощность для условий Т г ‘ = 65 °С, Т х ‘ = 35 °С, Т 0 ‘ = 25 °С. Так как 35 ближе к 25, чем к 65, пересчет мощности необходимо произвести, используя С. При помощи (10) найдем G’ = 0,0512 кг/с и для мощности получим W’ = 6,45кВт. (Расчет через С к дал бы заметно заниженные расход и мощность, соответственно 0,0397 кг/с и 5 кВт.)
  2. Пусть заданы температуры прямой воды и помещения и расход теплоносителя: Т г ‘ = 65 °С, Т 0 ‘ = 25 °С, G’ = 0,2 кг/с. Так как расход значительно больше, чем 2/3С к = 0,0793 кг/с, то для расчета мощности применима упрощенная закономерность через С к . Решая уравнение (8) относительно температуры холодной воды, получим: Т х ‘ ≈ 50,1 °С и для итоговой мощности: W’ ≈ 12,5кВт. (В этом случае пересчитывать через С нет необходимости, при большей трудоемкости мы получим значения: Тх’ = 49,9 °С и W’ = 12,7 кВт, то есть почти те же значения, что и выше, если учесть начальные погрешности измерений, проведенных для расчета С и С к .)
  3. Пусть теперь заданы температуры помещения, обратной воды и расход воды Т 0 ‘ = 25 °С, Т х ‘ = 40 °С, G’ = 0,05 кг/c. В этом случае расход значительно меньше, чем 1,15С к = 0,137 кг/с. То есть, пересчет произведем, используя (9). Сначала решим уравнение (10) относительно температуры горячей воды: Т г ‘ = 86,5 °С. Для мощности получим W’ = 9,8 кВт. Пересчет через С к привел бы к сильно заниженным выражениям для температуры и мощности: Т г ‘ ≈ 75,7 °С и W’ ≈ 7,5 кВт.

Примечание: Все утверждения данной статьи относятся к стандартным воздушным завесам с жидким теплоносителем. Измерения и данные, используемые в статье, относятся к максимальной скорости воздушного потока, когда коэффициенты рекуперации относительно низкие, и приведены для воздушных завес шириной около 1 м. Расходы теплоносителя в терминах (Т х – Т 0 ) / (Т г – Т 0 ) варьировались от достаточно больших (0,7) до малых (0,25). Выводы данной статьи нельзя экстраполировать на другие типы теплообменных систем «вода — воздух».

Поиск и представление формулы для мощности воздушной завесы в виде: W = 4,2 (αТ г + (1 – α) Т х – Т 0 ) С были обусловлены тем, что данный вид удовлетворяет критериям температурного подобия при любых α. Обработка большего объема опытных данных сможет привести к некоторым уточнениям α для различных типов воздушных завес. В этом случае для данных завес в формулах (9), (10) и (11) вместо коэффициентов 0,12 и 0,88 необходимо будет использовать коэффициенты α и 1 – α соответственно. Как показано выше, коэффициенты С в широком диапазоне расходов теплоносителя с хорошей точностью являются для воздушных завес константами. Однако, при ничтожно малых расходах теплоносителя G произведения Сα также станут малы. Это обусловлено физическим смыслом вышеуказанной формулы для мощности воздушной завесы W при данном предельном переходе.

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

Литература

  1. Методика оценочного расчета тепловых параметров водяных завес. Лысцев С. А., Азин А. В. Журнал «МИР КЛИМАТА» № 76.
  2. Основы теплопередачи. Михеев М. А. Издательство «Энергия», 1977 г., с. 55.
  3. Методика верификации воздушных завес с водяным нагревом, результаты испытаний. Журнал «МИР КЛИМАТА» № 71.
  4. Обновленная методика верификации воздушных завес с водяным нагревом. Журнал «МИР КЛИМАТА» № 77, с. 172.
  5. Протоколы измерений на стенде АПИК воздушной завесы Тропик Т109W от 18.01.13.
  6. Данные из каталога Frico для воздушных завес серий AD210W, AD310W и AD410W (см., например, www.frico-sales.ru).
  7. Паспорт верификации воздушной завесы Тропик Х525W. Журнал «МИР КЛИМАТА» № 77, с. 177.