7. Нестационарный теплообмен

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

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

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

На рис. 7.24 показан характер изменения температур на поверхности t_пови в глубине телаt_гл, внесенного в среду с более высокой температуройt_ж.

Рисунок 7.24 – Характер изменения температуры тела во времени

Любой процесс нагревания или охлаждения тела можно разделить на три режима. Первый из них охватывает начало процесса и характеризуется постепенным изменением температуры в глубь тела (участок τ₁). При этом скорость изменения температуры в отдельных точках различна, а температурное поле зависит от начального распределения температур в теле. Этот режим носит названиенеупорядоченного. В дальнейшем влияние начального распределения температур в теле исчезает, и относительная скорость изменения температуры в каждой точке тела становится постоянной величиной, зависящей от формы и размеров тела, его теплофизических свойств, условий теплообмена на поверхности тела. Такой упорядоченный режим нагрева названрегулярным(участок τ₂). По истечении длительного времени (теоретически по истечении бесконечно большого времени) наступает третий,стационарныйрежим, характерной особенностью которого является постоянство распределения температур во времени, в каждой точке тела. Если температура во всех точках тела одинакова и равна температуре окружающей среды, то наступает режимтеплового равновесия.

Решение задачи нестационарной теплопроводности заключается в определении зависимости изменения температуры во времени для любой точки тела tи количества подведенной или отведенной теплоты. Это может быть осуществлено аналитическим путем в результате решения дифференциального уравнения теплопроводности (7.23) приQ= 0 совместно с условиями однозначности.

Для технических целей в большинстве случаев ограничиваются рассмотрением течения процесса лишь в одном каком-либо направлении – направлении Х, например. В этом случае

.

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

1) распределением температуры на поверхности тела t_повв любой момент времени τ – граничное условие 1-го рода;

2) распределением плотности теплового потока по поверхности тела q_повво времени – граничное условие 2-го рода;

3) распределением температуры окружающей среды t_жи коэффициентом теплоотдачиот поверхности тела к этой среде:– граничные условия 3-го рода.

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

Для плоской неограниченной стенкитолщиной 2δ, омываемой жидкостью с температуройt_ж,

. (7.199)

Краевые условия при этом:

;;

;,

где t₀– температура в начальный момент времени, равная температуре поверхности стенки.

Для сплошного цилиндранеограниченной длины

. (7.200)

Краевые условия:

;;

;,

(R– радиус сечения цилиндра).

Для сплошного шара с радиусом R

. (7.201)

Краевые условия:

;;

;.

Решение уравнений (7.199) – (7.201) ввиду их сложности приводятся в специальной литературе.

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

или. (7.202)

Величина mв уравнении (7.202) является положительным числом, сохраняющим одно и то же значение для любой точки тела. Это число характеризует собой скорость охлаждения тела и называетсятемпомохлаждения; оно полностью определяется формой и размерами тела, значением тепловых параметров (a,,), интенсивностью теплообмена с окружающей средой и не зависит от начальных условий. Для ее определения достаточно измерить температуры телаt₁иt₂в произвольной точке тела в два определенных момента времени τ₁и τ₂:

.

В соответствии с теорией регулярного режима, разработанной Г.М. Кондратьевым, темп охлаждения m

, (7.203)

где – средний перепад температур по поверхности;– средний перепад температур по объему тела;FиV– поверхность и объем тела;с_р– удельная теплоемкость тела;– безразмерный коэффициент пропорциональности, называемыйкоэффициентом неравномерности распределения температуры в теле и зависящий для данного тела от условий теплообмена между его поверхностью и окружающей средой, т.е. от критерия Био (Bi = (l/); l – характерный линейный размер тел.

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

. (7.204)

Коэффициент Сзависит только от геометрической формы и размеров тела:

– для плоской неограниченной стенки толщиной 2δ ;

– для шара радиусом R;

– для цилиндра радиуса Rи длинойl;

– для параллелепипеда со сторонами b₁,b₂иb₃

.

По экспериментально найденным температурам t₁иt₂в произвольной точкетела определяют коэффициент его температуропроводности а, затем, рассчитавпо уравнению (7.203) величинуm, находят коэффициент теплоотдачи и скорость нагревания или охлаждения тела.

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

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

критерий Био ;

критерий Фурье

;

симплекс геометрического подобия

.

Безразмерная температура в любой точке тела (t₀– температура, принятая за масштаб температур) может быть выражена обобщенной зависимостью

. (7.205)

Эта функция постоянна для всех подобных процессов нестационарной теплопроводности.

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

и,

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

В качестве примера на рис. 7.25 приведены зависимости идля плоской пограничной стенки.

Расчет по таким графикам производится следующим образом. По известным значениям λ,с_р,ρнаходят величину коэффициента температуропроводностии рассчитывают критерий Фурье. Затем рассчитывают значение критерия Био, принимая за определяющий размер толщину стенкиδ. При этом величина коэффициента теплоотдачиαдолжна быть задана или рассчитана (при расчете α задаются температурой поверхности тела, которую потом проверяют, т.е.αнаходят методом последовательных приближений).

а

б

Рисунок 7.25 – Зависимости θ_пови θ_срот Вi и Fо для плоской неограниченной стенки:а–;б–

По пересечению перпендикуляра, проведенного из точки на оси абсцисс (рис. 7.25), соответствующей величине Вi и кривой, отвечающей значению Fо, определяют температуры θ_пови θ_ср. По этим температурам находят температуры на поверхностиt_пови в среднем сечении стенкиt_ср.

Результаты интегрирования уравнения теплопроводности можно представить также в виде зависимости

, (7.206)

где Q_τ– количество тепла, переданное за время τ;Q– полное количество тепла, переданное телу до наступления равновесия (до полного нагревания или охлаждения).

Расчет производят по соотношениям:

. (7.207)

В уравнении (7.207) V– объем тела;t– средняя по объему тела температура, достигнутая за время нагрева (эта температура задается или принимается).

На рис. 7.26 представлена зависимость от критериев Био и Фурье для плоской неограниченной стенки.

Рисунок 7.26 – Зависимость от Bi и Fо для плоской неограниченной стенки

Для определения времени нагрева стенки до температуры t необходимо из точки на оси абсцисс (рис. 7.26), соответствующей данному значению Bi, восстановить перпендикуляр до пересечения с ординатой, отвечающей значению отношения . Точка пересечения определит величину, по которой определяется искомое время нагрева τ.

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

Рисунок 7.27 – К расчету нестационарного теплообмена между двумя жидкостями

Если расход жидкости, протекающей через змеевикG₁, ее удельная теплоемкостьс₁, начальная температура (на входе в змеевик) постояннаt_1н= соnst, а конечная (на выходе)t_1кизменяется во времени, при этомt_1к<t_1н, то количество тепла, отданное этой жидкостью за времяdτ

. (7.208)

Это же тепло перейдет через стенку змеевика к другой жидкости:

, т.е., (7.209)

где Δt_ср– средняя разность температур жидкостей, переменная во времени.

Температура жидкости в аппарате при теплообмене изменяется от t_2ндоt_2к(t_2к>t_2н), ее значение в любой произвольный момент времениt₂. Тогда разность температур жидкостей у входа горячего теплоносителя в змеевикt_1н–t₂, а у выходаt_1к–t₂.

Следовательно, Δt_срв данный момент

. (7.210)

Подставив полученное значениеΔt_срв уравнение (7.209), получим

,

откуда . (7.211)

Если принять средние значения коэффициента теплопередачи Ки удельной теплоемкостиспостоянными, то.

Из уравнения (7.211) можно найти переменную конечную температуру жидкости в змеевике:

. (7.212)

При подстановке t_1киз уравнения (7.212) в уравнение (7.208) количество переданного теплаdQвыразится равенством

. (7.213)

Это же количество тепла, перешедшее в жидкость в аппарате, можно записать в виде

, (7.214)

где G₂ис₂– количество и удельная теплоемкость жидкости в аппарате, соответственно.

Приравнивая правые части последних уравнений, получим

. (7.215)

Разделив переменные и проинтегрировав уравнение(7.215)

,

найдем время нагревания жидкости в аппарате

(7.216)

либо

.

Последнее выражение может быть использовано также для определения температуры, до которой нагревается жидкость в аппарате t_2кза время τ_нагр.

Начальную и конечную температуры жидкости в змеевике можно рассчитать, пользуясь уравнением (7.212).

Так, при t₂=t_2ни разности температур жидкостейt_1н–t_2н

. (7.217)

При t₂=t_2ки разности температурt_1к–t_2к

. (7.218)

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

,

откуда . (7.219)

Приведенный расчет справедлив также и для случая охлаждения жидкости в аппарате. При охлаждении начальная, конечная и средняя температуры охлаждающей жидкости рассчитываются по уравнениям (7.217) – (7.219), соответственно, в которых знак сложения или вычитания изменяется на противоположный.

Рисунок 7.28 – К расчету нестационарного теплообмена между конденсирующимся паром и перемешиваемой жидкостью

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

Начальная температура жидкости в аппарате t_жн, конечная –t_жк. Температура греющего пара –t_п. За промежуток времениdτ количество тепла, переданное паром жидкости через поверхность нагреваF,составит

,

где K– коэффициент теплопередачи;t_ж– температура жидкости в данный момент времени.

Это тепло расходуется на подогрев жидкости

(G_ж,с_ж–количество и удельная теплоемкость жидкости соответственно).

Следовательно,

.

В результате разделения переменных и интегрирования последнего выражения, получим

;, (7.220)

откуда время нагревания жидкости в аппарате от t_жндоt_жк

. (7.221)

Уравнение (7.220) может быть использовано также для определения температуры нагреваемой жидкости t_жв любой момент времени τ:

. (7.222)

Расход греющего пара определяется по тепловому балансу

, (7.223)

где G_п– расход пара;r– удельная теплота парообразования.