Файл: Смольский Б.М. Нестационарный теплообмен.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 22.07.2024

Просмотров: 127

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Авторами работы [24] обнаружено, что с ростом тол­ щины стенки образца коэффициент нестационарного теплообмена повышается. Возможные причины такого результата состоят в следующем.

При разработке метода последовательных интерва­ лов использовалось условие Fo^0,5, что для длинных образцов приводит к необходимости при обработке тем­ пературных кривых разбивать процесс на интервалы вре­ мени по 20—30 сек и больше. Расчет дает всего одну точку для такого интервала (среднее значение для всего интервала). Такое усреднение может привести к сущест­ венным погрешностям, особенно на первых интервалах, где изменение теплового потока относительно велико.

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

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

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

128


ляет собой в этом случае некоторую эффективную вели­ чину, которая, как оказалось (рис. 32) [24], уменьшает­ ся практически в течение всего процесса нагрева и не стремится к постоянному квазистацибнарному значению.

Впроцессе нагрева (охлаждения) с уменьшением скорости изменения температуры должен наступить ин­ тервал времени, начиная с которого нагрев будет квазистационарным. В работах [34, 36—38] квазистационарный режим при охлаждении образцов в воздухе наступал либо через 5—15 сек,'либо процесс был квазистационарным с самого начала, но так или иначе он неизбежно возникал.

Несмотря на несоответствие некоторых результатов работы [24] результатам других работ, основные зако­ номерности, установленные в этой , работе, важны для понцмания механизма нестационарного теплообмена.

Висследованиях [105, 117] получено подтверждение найденной в работе [24] зависимости нестационарного коэффициента теплообмена при нагреве вырезки из плоскопараллельной пластины потоком воды от времени, размера и свойств материала образца. Для небольших

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

Установленные закономерности можно объяснить, используя рассмотренную схему процесса.

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

В процессе нагрева шара или вырезки из пластины в термостате с водой нестационарное распределение тем­ пературы в воде у поверхности шара в каждый момент времени будет иметь больший наклон (профиль будет более наполненным), чем квазистационарное распреде­ ление для тех же значений температуры стенки и жидко­ сти вдали от поверхности. Тепловой поток и коэффици­ ент теплообмена будут в'этом случае больше, чем в слу­ чае стационарного распределения температуры, что следует соответственно из законов Фурье и Ньютона.

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

9. Зак. 1284

129


также могут быть объяснены исходя из механизма не­ стационарного теплообмена. Если сравнивать зависимо­ сти температуры поверхности от времени двух шаров из одинаковых материалов, но с разной толщиной стенки, то для одного и того же момента времени скорости из­ менения температуры во времени будут для них различ­ ными, что вызывает отставание распределения темпера­ туры в жидкости различной степени. Это в свою очередь приведет к различным значениям теплового потока и ко­ эффициента теплообмена в рассмотренных случаях. При­ чем при увеличении б и ср замедляется скорость роста температуры поверхности в начальный период времени (см. рис. 40, а) и уменьшается qu/qСт или Nuh/Nuct- Этот качественный вывод согласуется с-полученными количе­ ственными результатами (см. рис. 47). Аналогичное рас­ суждение приводит к выводу, согласующемуся с опыт­ ными данными, который состоит в том, что с ростом вре­ мени величина Nuh/NuCt должна уменьшаться, стремясь к единице. Таким образом, в процессе нагрева тел дол­ жен-наступать момент, соответствующий определенной скорости изменения температуры поверхности, после которого коэффициент теплообмена остается постоян­ ным, равным его квазистационариому значению (рис. 42). Промежуток времени до этого момента зависит от ср, б (3.21). Очевидно, что коэффициент теплообмена в случае квазистационарпого процесса нагрева практиче­ ски равен коэффициенту теплообмена в соответствую­ щих стационарных условиях.

При решении задачи нестационарного теплообмена между потоком жидкости и участком стенки трубы авто­ ры работы [3] получили выражение для удельного теп­ лового потока на поверхности:

q ^ А ехр

(4.2)

Сравнение формул (3.12) и (4.2) показывает, что ха­ рактер зависимости теплового потока от основных пара­ метров в этих выражениях близок. Таким образом, полу­ чено экспериментальное подтверждение вида зависимо­ сти нестационарного теплового потока от т, ср и б, установленного путем теоретического решения задачи.

Из найденных зависимостей <?и/(7ст= /:(т) и NuH/NuCT= = /(/C TJ (рис. 48 и 50) следует, что с ростом параметра

ІЗО

/<т разница в интенсивностях теплообмена в нестацио­ нарных и стационарных условиях увеличивается, начи­ ная с некоторого значения параметра. Затем рост отно­ шения N u h / N u c t с увеличением /Ст замедляется. Уста­ новленные закономерности согласуются с результатами работы [17], полученными для случая теплообмена при течении газа в трубе.

При обработке экспериментальных данных оказалось возможным выразить обобщенную зависимость интен­ сивности теплообмена двумя способами: в зависимости от ср и б, а также в зависимости от /Ст (3.16) и (3.19). В указанных формулах левые части, как отмечалось вы­ ше, равны. Из этого следует, что параметр Ят однознач­ но определяется величинами б, ср и т. В частности, для граничных значений Ят между ними имеется зависи­ мость, определяемая формулой (3.22).

Следует подчеркнуть, что параметры б и ср оказыва­ ют влияние на нестационарный коэффициент теплооб­ мена, определяя скорость изменения температуры по­ верхности. При таком подходе к этому вопросу противо­ речия между основными выводами работ [24] и [17] становятся скорее формальными.

Можно также отметить, что из формул (3.16) и (3.19) более удобной для использования является зависимость N u h / N u c t и л и <7нА7ст о т 6, ср и X (3.16), так как значения температуры поверхности и ее производной, входящие в Кт (3.17), обычно заранее не известны.

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

В работах [36—38], как уже отмечалось выше, не бы­ ло обнаружено отклонение коэффициента теплообмена от его квазистационариого значения. Коэффициент а в этих условиях не зависел от скорости охлаждения, раз­ мера и материала шара и являлся постоянным для дан­ ного значения критерия Re.

Отсутствие в указанных опытах нестационарных эф­ фектов можно объяснить, по нашему мнению, небольши­ ми скоростями изменения температуры во времени.

На начальной стадии, процесса они непревышали 0,5—— 2,0°С/сеіс, что особенно в случае теплообмена с воздухом

5* 131


давало возможность пограничному слою перестраивать­ ся в соответствии с температурой поверхности. Измене­ ние размера и значения параметра ср шаров при соблю­ дении условия Ві<0,1 мало влияло на данную величину. Этот качественный вывод подтверждается результатами расчета коэффициента теплообмена по данным работы [38] с помощью разработанного нами метода (2.40). ^

Используя выражение (2.40), приведенное в гл. II, был подсчитан тепловой поток на поверхности для одно­ го из случаев, приведенных в работе, а именно для свин­ цового шара диаметром d = 29,9 мм (Я,—34,9 вт/м- °С,

а = 0,084 м2/час, ^0=ЮО°С, £ИІ = 21,9°С). Функция, описы­ вающая избыточную температуру на поверхности этого шара, имеет вид

£п= 78[1 — ехр (—0,011т)].

Подсчитанная по этой формуле температура на поверх­ ности дает отклонение от экспериментальных данных не более чем на 1,4%. Средний по поверхности коэффици­ ент теплообмена свинцового шара при скорости набегаю­ щего потока Шоо= 11,7 м/сек оказался равным а =

=81,2 вт]мг- °С и не зависел от времени, что согласуется

сданными, полученными в работе [38]. Это является косвенной проверкой разработанного нами метода опре­ деления условий на границе теплообмена между шаром

иокружающей средой.

Вработах [34, 46] при нагреве металлических шаров в потоке воздуха в отличие от работ [36, 38] было обна­ ружено отклонение коэффициента теплообмена от его квазистационарного значения.

Эти на первый взгляд противоречивые результаты могут быть объяснены следующим образом.

Тела, использованные в работах [34, 46], были тон­

костенными, и начальный нагрев их составлял 180 °С. Это позволяло получить в первые секунды охлаждения более высокие скорости изменения температуры во вре­ мени, чем в случае массивных шаров, использованных в работе [38], при начальном нагреве их до 80 °С. Кроме того, в исследованиях [34," 46] обращалось особое вни­ мание на регистрацию температуры в первые 2—4 сек процесса. Зарегистрированные скорости изменения тем­

пературы достигали 5—6°С/сек [34,

36]. В работах [37,

38] они не превышали 0,5—2 °С/сек.

При изменениях

132-


скорости в воздухе 5—6 °С/сек возможно появление не­ стационарных эффектов [16], ио полученные в работах [34, 36] отношения Nu„/NuCT, достигающие 2—2,5, види­ мо, завышены.

Установленный на основании опытов [38] квазистационарный характер теплообмена можно объяснить сле­ дующим образом. Перестройка распределения темпера­ туры в жидкости, как отмечалось выше-, в случае нестационарного теплообмена отстает от перестройки распределения температуры в образце. Это отставание тем значительнее, чем больше различаются между собой коэффициенты температуропроводности этих двух сред-. Поэтому рассматриваемый эффект труднее обнаружить в случае теплообмена металлических тел с воздухом, ко­ эффициенты температуропроводностей которых отлича­ ются всего на один порядок. Заметный эффект в этом случае будет при достаточной скорости изменения тем­ пературы поверхности во времени, что следует из рас­ смотренного выше механизма процесса. В работе [17] было обнаружено отклонение коэффициента теплообме­ на потока воздуха с металлическим телом от квазистационарного значения, достигающее 60% при скорости изменения температуры поверхности 200°С/се/с и 5% при скорости изменения температуры поверхности 5— 10 °С/сек.

Эффективное значение критерия Re и скорости потока жидкости, при которых проводились опыты в работе [117], были определены с помощью зависимости Крамерса для одиночных шаров, омываемых водой в стацио­ нарных условиях [114]:

NuCT= 2 ,0 + І.ЗРг0,15 + 0,66Pr°'31 Re0,5 .

Из опытов по стационарному теплообмену следует, что N u c t = - 2 4 7 . Число Рг в данном случае равно 3,02. Таким образом, значение критерия Re равно ~ ІО5 и скорость потока Шоо = 0,3 м/сек.

Значения этих парамеров были также оценены с по­ мощью зависимости для поперечного обтекания одиночт иых цилиндров потоком'жидкости [91]:

■NuCT= 0,25Re0,6 Рг«38 (Ргж/Ргп)0'25.

В этом случае значение Re оказалось также равным

~ ІО5.

133

По данным работы [17], отношение N uh/N uCt не за­ висит от числа Re. Это позволяет считать, что зависимо­ сти (3.16), (3.19) справедливы не только для числа Re=105, но и по крайней мере в некотором диапазоне чи­ сел Re.

Условия, в которых проведены исследования, и ре­ зультаты основных имеющихся в литературе работ при­ ведены в табл. 8. Приведенные данные подтверждают предложенное объяснение механизма процесса нестацио­ нарного теплообмена.

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

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

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

где

N = f (х±,

. . ., хп).

С учетом (4.3) выражение (2.39) приводится к виду (4.4). При дифференцировании выражения (2.39) неза­ висимыми переменными считаем величины A ii2, k u2, ti,2, /'i,2, X . Влияние последних двух переменных оценено прн-

134