Дифференциальное уравнение конвективного теплообмена уравнение фурье кирхгофа

Теплоотдача (конвективный теплообмен)

4.2.1. Определение исходных понятий

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

Теплоотдача зависит и от характера движения жидкости – ламинарного (слоистого) или турбулентного (вихревого). Напомним, критерий перехода ламинарного течения в турбулентное оценивается безразмерным числом Рейнольдса Re = , где r — плотность жидкости, h — вязкость жидкости, v – средняя скорость потока жидкости, l – характерные размеры канала (например, диаметр трубы). При больших числах Рейнольдса наблюдается турбулентное течение. Экспериментально выявлено, что критическое значение числа Рейнольдса, когда в круглой трубе происходит переход от ламинарного течения в турбулентное, примерно2×10 3 . Характер течения жидкости (газа) в трубах разного сечения будет одинаков, если каждому сечению соответствует одно и то же значение числа Re.

При любом виде течения (турбулентном или ламинарном) скорость жидкости в тонком слое у стенки канала равна нулю. Этот слой обычно называют вязким подслоем. На рис. 4-6 приведен график изменения температуры в жидкости при турбулентном течении. Температура в основном изменяется по толщине вязкого слоя, а в турбулентном потоке изменение температуры незначительное.

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

Явление теплоотдачи определяется тепловыми и гидродинамическими процессами. Эти процессы описываются соответствующими дифференциальными уравнениями: уравнением теплоотдачи, уравнением теплопроводности, уравнением динамики течения жидкости (газа), уравнением непрерывности. Рассмотрим эти уравнения.

4.2.2. Уравнение теплоотдачи

Интенсивность теплоотдачи (конвективного теплообмена) с поверхности стенки площадью dS за время t описывается эмпирической формулой Ньютона-Рихмана

где: T_ст , T_ж — температура стенки и жидкости соответственно; dS — элементарная площадка стенки; a — коэффициент теплоотдачи. Коэффициент теплоотдачи имеет смысл тепла, передаваемого в жидкость за единицу времени единицей поверхности при разности температур в один кельвин (один градус 0 С).

Теплопередача осуществляется через вязкий подслой толщиной l, для которой справедлив закон Фурье

dQ= — l t dS ,

где — градиент температуры в вязком подслое.

Сравнивая законы Ньютона-Рихмана и Фурье, получим:

a = — . (23)

Уравнение (23) называется уравнение теплоотдачи. Уравнение теплоотдачи позволяет по известному температурному полю определить коэффициент теплоотдачи. Заметим, в технических расчетах процесса теплоотдачи именно коэффициент теплоотдачи и представляет наибольший интерес.

4.2.3. Уравнение теплопроводности Фурье-Кирхгофа

Выделим в потоке жидкости элементарный кубик с ребрами dx, dy, dz (рис. 4-7) и определим поток тепла через грани этого элементарного кубика.

Пренебрежем изменением давления p и будем считать, что параметры жидкости – коэффициент теплопроводности l, плотность r и удельная теплоемкость жидкости при постоянном давлении с_р – постоянные константы. Из первого начала термодинамики следует, что при p = const подведенная теплота равна приращению энтальпии в выделенный объем жидкости:

Здесь оператор d указывает, что теплота является одной из функций процесса переноса внутренней энергии (наряду с работой dA = pdV): элементарная порция подведенной теплоты dQ в общем случае не является полным дифференциалом. Однако так как энтальпия (H = U + pV) как функция состояния является полным дифференциалом, то припостоянном давлении подведенная элементарная порция тепла dQ формально становится полным дифференциалом, т.е. можно формально записать: dQ = dQ.

Поток тепла через грани кубика обусловлен градиентом температуры на противоположных гранях, находящихся в перпендикулярном направлении на расстоянии dx, dy и dz соответственно. Если в направлении x на левой грани (см. рис. 4-7) температура T, то на противоположной грани на расстоянии dx температура равна (T + dx), где – соответствующий градиент температуры. Аналогично по направлениям y и z.

За время dt в направлении оси x количество тепла, проходящее через левую грань, равно (закон Фурье)

Q_x1 = , где dS = dx dy – площадь грани,

а через правую грань Q_x2 = .

Имеем: dQ_x = Q_x1 — Q_x2 = .

Аналогично по осям y и z:

dQ_y = , dQ_z = .

За счет подвода тепла наблюдаем приращение энтальпии жидкости за время dt определится уравнением

dQ = dH = dQ_x + dQ_y + dQ_z = . (24)

С другой стороны подведенную теплоту dQ (изменение энтальпии dH) можно представить также уравнением

dQ =c_pr dx dy dz, (25)

где dm = r dx dy dz – масса рассматриваемого элементарного объема, c_p – теплоемкость жидкости при постоянном давлении, – изменение температуры рассматриваемого объема за время dt вследствие притока тепла.

Здесь имеет смысл субстанциальной производной. Дело в том, что изменение температуры в текущей жидкости определяется двумя факторами: 1) изменением температуры во времени в рассматриваемом объеме жидкости; 2) изменением температуры вследствие перемещения этого объема от одной точки пространства к другой:

= + + + .

Здесь — локальное изменение температуры во времени (изменение температуры в данной точке пространства). Величина

+ + = + + ,

представляет собой конвективное изменение температуры, обусловленное перемещением элементов текущей жидкости (v_x, v_y, v_z – проекции скорости элементов жидкости). Таким образом, температурное поле в движущейся жидкости зависит и от распределения скоростей элементов жидкости.

Сравнивая (24) и (25) получим уравнение

c_pr dx dy dz = или

= . (26)

Уравнение (26) есть дифференциальное уравнение теплопроводности Фурье-Кирхгофа, где постоянная жидкости a = называется коэффициентом температуропроводности. Из уравнения Фурье-Кирхгофа следует, что коэффициент температуропроводности a характеризует скорость изменения температурного поля в среде. Уравнение (26) можно кратко записать посредством оператора Лапласа Ñ 2 = + + :

= . (26 * )

Запишем уравнение Фурье-Кирхгофа в развернутом виде

+ + + = . (26 ** )

Для твердых тел уравнение Фурье-Кирхгофа пример вид

= ,

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

4.2.4. Уравнение Навье-Стокса (дифференциальное уравнение

Движения жидкости)

Как было уже отмечено, из уравнения Фурье-Кирхгофа следует, что температурное поле зависит от распределения скоростей элементов жидкости (см., напр., уравнение 26 ** ). Динамика движения произвольного элемента жидкости массой dm описывается вторым законом Ньютона: dma = R, где R – результирующая сила, действующая на элемент.

Выделим в потоке жидкости элемент жидкости массой dm в форме кубика с ребрами dx, dy, dz (рис. 4-8). Жидкость течет по вертикали вниз. Ось 0x совпадает с левой стенкой. Скорость тонкого пристеночного слоя жидкост равна нулю. По мере удаления от стенки скорость v_x возрастает (на рисунке приведена эпюра скоростей жидкости).

На выделенный элемент жидкости действуют: 1) сила тяжести; 2) равнодействующая сил давления по вертикали; 3) равнодействующая сил вязкого трения на боковые поверхности выделенного элемента, и обусловленная градиентом скорости жидкости в направлении оси 0y.

Сила тяжести равна

где g_x – проекция ускорения свободного падения на ось 0х, r – плотность жидкости. Будем считать, что жидкость несжимаемая, т.е. r = const.

Равнодействующая сил давления. На верхней грани давление равно p,на нижней – (p + ). Площадь верхнего и нижнего граней dS = dy dz, следовательно, равнодействующая сил давления равна

p dy dz — (p + ) dy dz = — dy dz = — . (28)

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

F_t =h , (29)

где h – коэффициент вязкости (или просто – вязкость), единицей вязкости служит паскаль×секунда (Па×с).

Если на левой боковой грани вязкое трение равно F_t , то на правой грани – (F_t + ). Так как скорость жидкости растет слева направо, то сила трения на левой грани площадью dS = dx dz направлена вверх, а на правой той же площади – вниз. Равнодействующая сил трения равна:

(F_t + ) dx dz – F_t dx dz = dy dz = .

Воспользовавшись законом силы вязкого трения (29), получим при h = const:

= h dV . (30)

Если проекция скорости v_x изменяется по всем направлениям, т.е. проекция скорости v_x имеет градиент по всем направлениям, то равнодействующая вязкого трения в направлении оси 0х примет вид:

h + + dV или hÑ 2 v_x dV . (31)

Итак, проекция результирующей силы R, действующей на выделенный элемент массой dm и объемом dV, на ось 0х примет вид:

R_x = [r g_x — + h + + ] dV.

С другой стороны

dma_x = r dV = r + + + dV,

т.е. полное изменение скорости определяется: 1) локальной составляющей – изменением скорости во времени в данной точке пространства и 2) изменением в пространстве вследствие течения жидкости – конвективной составляющей + + . Итак:

r dV = r + + + dV.

Теперь можно записать второй закон динамики в развернутом виде:

r + r + + = r g_x — + h + + . (32)

Дифференциальное уравнение (32), описывающее динамику движения вязкой несжимаемой жидкости (r = const), называется уравнением Навье-Стокса. Уравнение Навье-Стокса для проекций на оси 0y и 0z имеют аналогичный вид:

r + r + + = r g_y — + h + + ; (32 * )

r + r + + = r g_z — + h + + . (32 ** )

Заметим, уравнение Навье-Стокса справедливо как для ламинарного, так и турбулентного потока жидкости.

4.2.5. Уравнение непрерывности (сплошности)

Выделим в потоке жидкости мысленный неподвижный элементарный кубик с ребрами dx, dy, dz (рис. 4-9) и определим массу жидкости, протекающей через этот кубик за время dt.

В направлении оси 0x через левую грань с площадью dS = dydz втекает за время dt масса жидкости

а через правую грань той же площади вытекает масса

m_x2 = dy dz.

Имеем: dm_x = m_x2 — m_x1 = dx dy dz dt.

Аналогично по остальным осям:

dm_y = m_y2 — m_y1 = dx dy dz dt;

dm_z = m_z2 — m_z1 = dx dy dz dt.

Полное изменение массы жидкости в выделенном объеме dV = dx dy dz равно

dm = dm_x + dm_y + dm_z = + + dV dt.

Ясно, что изменение массы жидкости в выделенном объеме с течением времени может быть обусловлено только изменением плотности жидкости со временем. Следовательно, имеем:

+ + dV dt = — dV (33)

Знак минус в (33) указывает, что если dm > 0 (т.е. из объема dV вытекает большая масса жидкости, чем втекает), то 0, т.е. когда плотность жидкости в выделенном объеме увеличивается.

Уравнение (33) можно записать в виде:

+ + + = 0 (34)

Дифференциальное уравнение (34) и есть уравнение непрерывности (уравнение сплошности) в общем виде. Если плотность жидкости постоянна, (r = const) то уравнение непрерывности (34) примет вид

+ + = 0.

4.2.6. Условия однозначности математического описания

теплоотдачи (краевые условия). Подобие явлений

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

1) геометрические условия, которые характеризуют размеры и форму, например, диаметр и длина, гладкость или шероховатость стенок канала протекания жидкости и т.п.;

2) физические константы, характеризующие свойства среды, например, плотность, вязкость, коэффициент теплопроводности и т.д.;

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

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

Методология теории подобия основывается на трех утверждениях (трех теоремах подобия):

1) подобные явления имеют одинаковые критерии подобия (Бертран Ж.). В модельном эксперименте измеряют те величины, которые содержатся в критерии;

2) уравнения, описывающие процесс, могут быть представлены в виде зависимости между критериями подобия, определяющими процесс (Федерман А., Букингем Дж.). Результаты опыта представляют в виде критериальных уравнений;

3) подобны те явления, краевые условия которых подобны и для которых критерии подобия, составленных из краевых условий, равны (Кирпичев М.В., Гухман А.А.). Это утверждение устанавливает признак подобия процессов в модельном и рабочем объектах.

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

Число Эйлера Eu – это соотношение между перепадом силы давления в потоке F_p = Dpl 2 и поступательной силой инерции F_i = ma = rl 3 :

Eu = = = .

Число Эйлера характеризует соотношение между перепадом статического давления и динамическим давлением в потоке.

Число Рейнольдса Re определяет связь в потоке между поступательной силой инерции F_i = m и силой вязкого трения Ньютона F_h = h l 2 :

Re = = .

Число Рейнольдса характеризует гидродинамический режим вынужденного движения жидкости (газа).

Критерии теплового подобия получают из уравнений передачи тепла. Число Фурье Fo – отношение тепла, передаваемого теплопроводностью Q_l = l l 2 t и поступающим в тело массой m теплом Q = rl 3 cDT:

Fo = = = , где а – температуропроводность.

Число Фурье характеризует нестационарные процессы распространения тепла, определяемое скоростью изменения температурного поля, размерами тела, температуропроводностью.

Число Нуссельта Nu – отношение тепла передаваемого в процессе теплоотдачи Q_a = aDTl 2 t к теплу в процессе теплопроводности Q_l = l l 2 t:

Nu = .

Число Нуссельта характеризует интенсивность теплообмена на границе твердое тело (стенка) – жидкость.

Приложение:формулы и теоремы векторного исчисления

Векторные величины в тексте обозначены прямым полужирным шрифтом, скалярные – курсивом. Например, b – вектор, b – скаляр. Скалярное произведение векторов a и b записывается как(ab). Векторное произведение обозначено заключением сомножителей в квадратные скобки через запятую. Например, векторное произведение векторов a и b записывается как[a, b].

1.1. Скалярное произведение векторов abназывается произведение модулей векторов на косинус угла между этими векторами:

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

где i, j, k – единичные орты, ii = jj= kk = 1, ij= jk = ki =0.

1.2. Векторное произведение двух векторов[a, b]. В результате векторного произведения векторовaи b получается новый вектор c: [a, b] = c. Модуль результата векторного произведения – модуль вектора с— по определению равен c = ab sin (a^b). В декартовой системе координат векторное произведение имеет вид

[a, b] = .

1.3. Смешанное (скалярно-векторное) произведение трех векторов а[bc]. Результатом такого произведения является скалярная величина:

Смешанное произведение допускает циклическую (стековую) перестановку векторов:

По своему геометрическому смыслу смешанное произведение равно объему параллелепипеда, построенного на векторах а, bиc.

1.4. Двойное векторное произведение[а,[bc]]:

Здесь вначале выполняются скалярные произведения, стоящие в скобах (т.е. скалярные произведения acи ab), и только затем вектораb и c умножаются на соответствующие результаты скалярного произведения. Компоненты результата двойного векторного произведения имеют вид:

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

2.1. Векторный оператор набла Ñ (оператор Гамильтона) как символический вектор:

Ñ = + + .

2.2. Градиент скалярной функции (скалярного поля) j(x, y, z) – это вектор, направленный в сторону быстрейшего увеличения j и равный производной по этому направлению. В координатном представлении градиент имеет вид:

gradj = Ñj (в развернутом виде Ñj = + + ).

2.3. Дивергенция некоторого вектора Eв данной точке пространства– это поток вектора Eиз бесконечно малого объема dV, находящегося в этой точке пространства. Объем dV — источник или сток вектора E. Дивергенция вектора E – скалярная величина и определяется скалярным произведением оператора набла на вектор E:

div Eº ÑE. В развернутом виде ÑE = + + .

2.4. Ротор некоторого вектора E –это вектор, порождающий циркуляцию некоторого другого вектора по бесконечно малому контуру. Например, в уравнении Максвелла rot Е= вектор является ротором, который порождает вихревой вектор – напряженность вихревого электрического поля Е. Ротор вектора E – векторная величина и кратко записывается векторным произведением оператора набла на вектор E:

В развернутом виде [Ñ, E] = i + j + k).

2.5.Некоторые векторные тождества.

1. Ñ(Ñj) = Ñ 2 j = , где j — скалярная функция.

3. Ñ[Ñ, E] = 0, где E — векторная функция.

2.6.Теоремы векторного анализа.

1. Теорема Гаусса-Остроградского. Здесь объем V ограничен замкнутой поверхностью S. Вектор dS в данной точке замкнутой поверхности направлен по внешней нормали к поверхности.

а) = или в форме = ;

б) = или в форме =

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

а) = или в форме = ;

б) = или в форме = .

ЛИТЕРАТУРА

1. Савельев, И. В. Курс общей физики [Текст]: учеб. пособие. Кн.1., Кн.2., Кн.3, Кн.4 – М.: Наука, 2003.

2. Семенов, Н.А. Техническая электродинамика [Текст]: учеб. пособие для вузов. – М.: Связь, 1973. – 480 с.

3. Матвеев, А.Н. Электричество и магнетизм [Текст]: учеб. пособие для вузов. – М.: Высшая школа, 1983. – 463 с.

4. Алешин, Н.П. и др. Ультразвуковой контроль: учебное пособие/ под общ. ред. В.В. Клюева – М.: Изд. Дом «Спектр», 2011. – 224 с., ил.

5. Будадин, О.Н. и др. Тепловой контроль: учебное пособие/ под общ. ред. В.В. Клюева – М.: Изд. Дом «Спектр», 2011. – 176 с., ил.

Содержание

1. ОСНОВЫ ТЕОРИИ КОЛЕБАНИЙ И ВОЛН……………………………….7

1.1. Основные положения теории колебаний осцилляторов………………7

1.1.1. Уравнения колебаний гармонического осциллятора…………….7

1.1.2. Методы (способы) представления колебаний……………………10

1.1.3. Динамика гармонического осциллятора…………………………12

1.1.4. Динамика затухающих колебаний………………………………..15

1.1.5. Динамика вынужденных колебаний. Импеданс колебательной

1.1.6. Импеданс и фазовые соотношения между смещением и

вынуждающей силой. Резонанс смещения……………………….21

1.1.7. Импеданс и фазовые соотношения между скоростью и

вынуждающей силой. Резонанс скорости………………………..23

1.1.8. Устойчивость амплитуды вынужденных колебаний……………24

1.1.9. Добротность осциллятора…………………………………………26

1.1.10. Добротность и резонансная кривая поглощения

1.2. Исходные положения теории гармонических волн…………………….29

1.2.1. Волновое уравнение, описывающее бегущую волну в струне.

Фазовая скорость волны. …………………………………………29

1.2.2. Стоячие волны на струне………………………………………….33

1.2.3. Импеданс среды (на примере струны)…………………………….36

1.2.4. Отражение и прохождение волны на границе двух сред

с различными импедансами. Согласование импедансов………. 37

1.2.6. Затухание реальных волн в среде…………………………………41

2.1. Волновое уравнение плоской волны в газе и жидкости………………..42

2.2. Скорость акустической волны в твердом теле………………………….46

2.3. Поле плоского излучателя ультразвуковой волны……………………..47

3. ЭЛЕКТРОМАГНИТНОЕ ПОЛЕ……………………………………………..48

3.1. Электростатическое поле…………………………………………………48

3.1.1. Объекты электростатики…………………………………………..48

3.1.2. Закон Кулона. Напряженность электростатического поля и

интегральная теорема Гаусса……………………………………..49

Дифференциальное уравнение конвективного теплообмена

Для вывода дифференциального уравнения конвективного теплообмена в установившемся потоке жидкости (газа) выделяют элементарный параллелепипед с гранями dx, dy, dz (рис. 7.2) и составляют для него тепловой баланс, приняв физические параметры λ, с_р и ρ постоянными. Скорости движения жидкости в направлении осей x, y и z, соответственно – w_x, w_y и w_z. Температура жидкости t изменяется вдоль граней параллелепипеда. Теплоперенос в жидкости осуществляется путем конвекции и теплопроводности. Все подведенное к параллелепипеду тепло затрачивается только на изменение его энтальпии.

Количество тепла, которое вводится жидкостью путем конвекции по направлению оси х через грань dydz за время dτ,

.

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

,

.

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

.

Аналогично в направлении осей y и z

;

.

Полная разность между количествами введенного и выведенного путем конвекции тепла в объеме параллелепипеда за время dτ

.

На основании уравнения неразрывности потока (1.34)

.

Тогда конвективная составляющая теплового потока

,

где – объем параллелепипеда.

Количество тепла, которое вводится в параллелепипед жидкостью путем теплопроводности за время dτ, в соответствии с уравнением (7.21)

.

Суммарное количество тепла, подведенное конвекцией и теплопроводностью

Это тепло dQ, согласно первому началу термодинамики, равно изменению энтальпии жидкости в объеме параллелепипеда

.

.

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

,

, (7.66)

где – коэффициент температуропроводности.

Полученное уравнение (7.66) является дифференциальным уравнением конвективного теплообмена, или уравнением Фурье-Кирхгофа. Оно выражает в наиболее общем виде распределение температур в движущемся потоке жидкости (газа).

Для твердых тел w_x = w_y = w_z = 0 и уравнение (7.66) превращается в дифференциальное уравнение теплопроводности (7.23).

При установившемся процессе теплообмена , тогда уравнение конвективного теплообмена для этого случая

. (7.67)

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

Тепловое подобие

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

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

.

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

.

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

. (7.68)

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

Аналогом критерия Нуссельта при нестационарном теплообмене между твердым телом и жидкостью или газом является критерий Био:

, (7.69)

где – коэффициент теплопроводности твердого тела.

Условия подобия в ядре потока определяются из уравнения Фурье-Кирхгофа (7.66). Для одномерного движения потока жидкости его можно представить в виде

. (7.70)

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

;

.

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

. (7.71)

Комплекс является мерой соотношения между теплом, переносимым конвекцией и теплопроводностью при конвективном теплообмене, и носит название критерия Пекле:

. (7.72)

Необходимым условием теплового подобия являются предшествующие ему гидродинамическое и геометрическое подобия. Гидродинамическое подобие определяется критериями гомохронности Но, Рейнольдса Re и Фруда Fr, а геометрическое – постоянством отношений основных геометрических размеров поверхности (стенки) к некоторому характерному размеру L. В качестве характерного линейного размера для трубных поверхностей, которые часто используются в качестве теплопередающих, принимают диаметр трубы L = d, иногда – длину трубы L = l, радиус кривизны изогнутой трубы L = R и т.д. В результате критериальное уравнение конвективного теплообмена выражается функцией вида

, (7.73)

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

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

. (7.74)

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

,

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

По физическому смыслу критерий Pr определяет соотношение полей скоростей и температур в потоке. Для капельных жидкостей он зависит от температуры (с повышением температуры увеличивается), а его численные значения лежат в пределах 3 ÷ 300. Для газов критерий Прандтля не зависит ни от температуры, ни от давления, являясь постоянной величиной для газов одинаковой атомности (для одноатомных газов Pr = 0,67; для двухатомных Pr = 0,72; для трехатомных Pr = 0,8; для четырехатомных и более Pr = 1).

С введением критерия Прандтля уравнение (7.74) принимает вид:

. (7.75)

При установившихся тепловых процессах и одинаковой геометрической форме потоков из уравнения (7.75) исключаются критерии Fo и Ho, содержащие время, и геометрические симплексы, тогда

. (7.76)

При рассмотрении конкретных задач теплообмена уравнение (7.76) может быть видоизменено.

Так, при вынужденном движении, когда влияние сил тяжести на гидродинамику потока пренебрежительно мало, критерием Фруда можно пренебречь и критерий Нуссельта будет являться функцией лишь двух критериев:

. (7.77)

При свободном движении жидкости (в условиях естественной конвекции) из критериального уравнения исключают критерий Рейнольдса:

. (7.78)

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

,

где ρ и ρ₀ – плотности холодной и нагретой жидкости (газа).

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

. (7.79)

Таким образом, критерий Gr представляет собой определяющий критерий теплового подобия при естественной конвекции, когда движение жидкости целиком обусловлено самим процессом теплообмена, а уравнение конвективного теплообмена для этого случая

. (7.80)

Для газов при Pr ≈ 1 уравнение (7.79) превращается в функциональную зависимость критерия Нуссельта от критерия Грасгофа.

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

критерий Стантона , (7.81)

который отражает соотношение количества тепла, передаваемого конвекцией, и тепла, переносимого движущимся потоком;

критерий Релея ; (7.82)

критерий Грэтца , (7.83)

(G – массовый расход жидкости) характеризует конвективный теплоперенос при ламинарном режиме.

.

В соответствии с полученными критериальными уравнениями обрабатываются опытные данные по переносу тепла конвекцией и представляются в виде степенных уравнений. Например, уравнение (7.77) – в виде ; уравнение (7.79) – в виде и т.д.

При решении практических задач по найденному из соответствующего критериального уравнения значению критерия Нуссельта определяют коэффициент теплоотдачи:

.

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

ЛЕКЦИЯ 10 ПОТЕНЦИАЛ ПЕРЕНОСА. УРАВНЕНИЕ ФУРЬЕ-КИРХГОФА ЭЛЕМЕНТЫ ТЕОРИИ ПОДОБИЯ В ТЕПЛООБМЕНЕ. Потенциал переноса

Фёдор Коновалов 3 лет назад Просмотров:

1 ЛЕКЦИЯ 10 ПОТЕНЦИАЛ ПЕРЕНОСА. УРАВНЕНИЕ ФУРЬЕ-КИРХГОФА ЭЛЕМЕНТЫ ТЕОРИИ ПОДОБИЯ В ТЕПЛООБМЕНЕ Потенциал переноса Рассмотрим газовую или жидкую сплошную среду. Примем: все точки среды находятся в неравновесном состоянии. Это приводит к возникновению полей концентраций, температур, давлений, а наличие градиентов этих параметров вызывает перенос массы и энергии. Выделим элемент объема движущейся жидкости в неоднородном поле некоторого потенциала переноса. Под потенциалом переноса понимают удельную массу или энергию (отнесѐнную к единице объѐма). (x, y, z) — скалярная величина. Из курса математики известно, что скалярная функция называется потенциалом векторной функции, если между ними существует связь в форме: q grd Под градиентом скалярной функции (10.1) grd подразумевают векторную функцию: grd i j k (10.) x y z где i, j, k базисные векторы или орты. В дальнейшем будем понимать связь между q и как пропорциональность q grd (10.) Таким образом, поток переносимой субстанции (массы или энергии) является векторной величиной q. В случае переноса массы под потенциалом переноса обычно понимают концентрацию компонента в смеси: mi (10.3) i где mi масса i-го компонента в объѐме, кг i; i концентрация i-го компонента в смеси, кг i/м 3. 1

2 При переносе энергии в качестве потенциала переноса рассматривают энтальпию единицы объема среды: c P cp (10.4) Здесь: ср изобарная теплоѐмкость среды, Дж (кг К) ; температура, К; — плотность, кг/м 3 ; объем, м 3 c Дж c Дж P P 3 м В рассматриваемой среде могут существовать, так называемые, объемные (непрерывно распределѐнные по объѐму) источники или стоки массы и энергии. В химической технологии под ними подразумеваются химические превращения. Известно, что процессы тепло- и массообмена осуществляются двумя основными механизмами: молекулярным и конвективным. Молекулярный перенос (диффузия, теплопроводность) возникает в результате стремления системы к термодинамическому равновесию, а конвективный вызывается наличием поля скоростей в жидком или газовом объѐме. Следует отметить, что в случае переноса энергии в форме теплоты существует ещѐ и радиантный перенос (тепловое излучение), вклад которого учитывают при достаточно высоких температурах. Молекулярный перенос является определяющим в неподвижных средах, хотя он вызывает естественную конвекцию и практически всегда ею сопровождается. Процессы молекулярного переноса массы и энергии описываются соответствующими феноменологическими уравнениями, являющимися, как правило, линейными градиентными законами. закону Фика: Молекулярный перенос массы (молекулярная диффузия) подчиняется первому qmm -D grd i (10.5) где коэффициент молекулярной диффузии, м с; q Mm — плотность массового потока, кг i /м с. Молекулярный перенос энергии в форме теплоты описывается законом Фурье: q M — grd (10.6)

3 где коэффициент теплопроводности, Вт/ м K; q M — плотность теплового потока, Дж /с м = Вт/ м В более общей форме закон Фурье можно переписать следующим образом: q M — cp grd cp — grd cp c grdc P P (10.7) Здесь коэффициент температуропроводности, м /с. c P Следует обратить внимание, что коэффициенты диффузии D и температуропроводности а имеют одинаковую размерность (м /с) и называются молекулярными коэффициентами переноса. Таким образом, молекулярный перенос массы и энергии описываются одинаковыми по форме законами, и они могут быть обобщены следующим выражением: q M -k grd (10.8) При конвективном переносе масса и энергия транспортируются макроскопическим путѐм, движущейся со скоростью v средой. Плотность конвективного потока массы и энергии на каждом участке поверхности А можно выразить следующим образом: v A q K v (10.9) A где А участок поверхности, ориентированный перпендикулярно вектору скорости v. Размерность q плотность потока массы или энергии, соответственно. q K v c P (10.10) Таким образом, в случае молекулярного и конвективного переноса общая плотность потока массы или энергии складывается из двух векторных величин: q q M q K (10.11) 3

4 Рис К выводу балансовых уравнений переноса В газовой или жидкой среде, находящейся в движении, выделим произвольный объѐм, ограниченный поверхностью А (рис. 10.1). На поверхности А выделим элемент поверхности dа и представим его в векторной форме, умножив на единичный вектор, нормальный к этому элементу и направленный из объѐма,. Составим балансовое уравнение по типу: Накопление внутри объёма = Вход Выход + Образование Примем, что в произвольном объеме нет источников субстанции или стоков, т.е. образование равно нулю. Плотность потока субстанции через элементарную площадку будет Знак в этом произведении делает входящие потоки положительными, а выходящие отрицательными. Результирующий поток массы или энергии (Вход минус Выход) будет получен суммированием всех потоков через замкнутую поверхность A: A q da (10.1) Таким образом, физически этот интеграл представляет разницу между входящими и выходящими потоками субстанции через всю поверхность А. 4

5 Если в объѐме происходит накопление субстанции, то это вызовет изменение d потенциала переноса во времени, которое для элементарного объѐма d можно dt d представить как d, а для всего объема как интеграл: dt M d d dt (10.13) Приравняв выражения (10.1) и (10.13), получим: q da A d d dt (10.14) Согласно теореме Остроградского-Гаусса, дающей преобразование интеграла, взятого по объѐму, ограниченному поверхностью A, в интеграл, взятый по этой поверхности, будем иметь: q da A div q d С учѐтом (10.15) соотношение (10.14) примет вид: (10.15) div q d t (10.16) Интеграл, взятый по произвольному объѐму, может быть равен нулю только в случае равенства нулю подынтегральной функции: div q 0 t (10.17) Полученное выражение (10.17) и есть основное дифференциальное уравнение переноса субстанции массы или энергии, как будет показано далее. В случае изотропных сплошных сред с его помощью можно получать поля температур или концентраций в однофазной среде. Искомой величиной является плотность потока субстанции, которая определяет удельный поток массы или энергии. 5

6 Дифференциальное уравнение конвективного теплообмена (Уравнение Фурье-Кирхгофа) Дифференциальное уравнение конвективного теплообмена является частным случаем уравнения переноса энергии в форме теплоты в однофазной сплошной изотропной среде. При этом теплоѐмкость с P, теплопроводность и плотность среды считаются постоянными; отсутствует также перенос энергии в форме теплового излучения и объѐмное источники (стоки) теплоты. Как было отмечено выше, потенциалом переноса теплоты является энтальпия единицы объѐма среды (ур. 10.4): Где q M вид: q c P (10.4) Тогда с учѐтом выражений (10.7) и (10.9), будем иметь: q q q grd c c (10.18) M K v P плотность потока теплоты, представляющая векторную сумму молекулярной и конвективной q K компонент. Основное уравнение переноса субстанции (10.17) в этом случае примет следующий cp t t div div При = const получим: c c v 0 P grd P P (10.19) grd (10.0) v 0 div grd (10.1) — оператор Лапласа. Дифференциальная операция скалярную функцию сопоставляет скалярную функцию и учѐтом div grd (10.) x y z Примем также, что гидродинамически среда является стационарной, тогда div v уравнение неразрывности имеет вид: 0 (10.3) с 6

7 Дивергенцию от v, как произведения векторной и скалярной величины, можно представить в виде: div v divv v grd (10.4) С учѐтом (10.1), (10.3) и (10.4) выражение (10.0) примет вид: v grd t 7 (10.5) Полученное выражение (10.5) называется дифференциальным уравнением конвективного теплообмена или уравнением Фурье-Кирхгофа. Оно является частным случаем дифференциального баланса энергии в форме теплоты в движущейся среде, где имеет место перенос энергии теплопроводностью. Полная форма уравнения конвективного теплообмена в скалярном виде будет: t v x vy vz (10.6) x y z x y z Левая часть этого соотношения представляет собой субстанциональную производную: D Dt vx vy vz t x y z Решением дифференциального уравнения конвективного теплообмена в общем виде является функция температур в движущейся среде. В неподвижной среде t, которая представляет собой нестационарное поле и выражение (10.5) принимает вид: (10.7) Уравнение (10.7) описывает нестационарное температурное поле в неподвижной среде, может применяться также для твѐрдых тел и нестационарной теплопроводности. имеет вид: называется уравнением Уравнение для установившегося процесса в неподвижной среде или в твердом теле 0 (10.8) Отметим, что согласно уравнению (10.7), локальное изменение температуры пропорционально коэффициенту температуропроводности, который, таким образом, характеризует теплоинерционные свойства среды. При прочих равных условиях быстрее нагреется или охладится то тело, которое имеет бόльший коэффициент температуропроводности. t

8 Граничное условие Уравнения Фурье-Кирхгофа на практике используется совместно с граничным условием, т.е. условием на границе среды у неподвижной твердой стенки. Вблизи твердой стенки теплота передается только теплопроводностью внутри пограничного слоя. Следовательно, по закону Фурье: Q da grd da n (10.9) В тоже время, количество теплоты, передаваемой из ядра потока к твердой стенке, можно выразить законом Ньютона (уравнение теплоотдачи): Q ст da (10.30) Если перенос тепла стационарный, это один и тот же тепловой поток: Q grd da ст da grd ст (10.31) Это и будет граничным условием, дополняющим уравнение Фурье-Кирхгофа. Элементы теории подобия в теплообмене Рассмотрим гидродинамически одномерный поток жидкости. Запишем уравнение Фурье-Кирхгофа: t vz z (10.3) Получим приближенное решение этого уравнения методами теории подобия. Для этого зададим константы подобия, выражающие отношения величин, входящих в уравнение Фурье-Кирхгофа:, t. v. Умножим каждый из элементов дифференциального уравнения (10.3) на соответствующую константу подобия, причем последняя как постоянная величина, выносится за знак дифференциала. v v z (10.33) t t z Для сохранения тождественности полученного и исходного уравнений необходимо выполнение следующего условия: 8

9 t v (10.34) Разделим поочередно дроби на правую дробь. Отношения будут равны единице, т.к. они все являются индикаторами подобия, а у подобных явлений индикаторы равны единице.после деления заменим константы подобия их значениями: = 1/, t = t1/t, а = а1/ а, = 1/. Первый комплекс при делении / 1 t Fo 9 t (10.35) Критерий Фурье Fo является мерой отношения между количеством теплоты, вызывающей изменение температуры в данной точке движущейся среды и количеством теплоты, передаваемой данной движущейся среде теплопроводностью. Равенство критериев Fo в сходственных точках подобных систем — необходимое условие подобия нестационарных процессов. Критерий Fo — аналог критерия Ho в гидродинамике. Второй комплекс при делении v / 1 v Pe (10.36) Критерий Пекле Pe характеризует отношение количества теплоты, передаваемой конвекцией, к количеству теплоты, передаваемой теплопроводностью. Критерий Pe можно представить так: v Pe Re Pr где: Pr (10.37) Критерий Прандтля Pr — мера соотношения между толщиной гидродинамического пограничного слоя δгидр и толщиной теплового пограничного слоя δтепл. Критерий Прандтля Pr составлен только из физических параметров. В газах при Pr = 1 поля температур и скоростей подобны, толщины теплового и гидродинамического слоев соизмеримы по масштабу δтепл δгидр. При Pr = 0,7 1, толщины теплового и гидродинамического слоев практически равны по величине. Здесь Pr мало зависит от температуры и давления.

10 В жидкостях Pr = Поэтому в капельных жидкостях толщина теплового слоя δтепл меньше толщины гидродинамического δгидр слоя, δтепл 11 3 g Ar 1 (10.40) 1 Зависимость плотности от температуры можно записать как: 1 1 ( 1 ) (10.41) где β — коэффициент объемного расширения (1/град). Отсюда 1 ( 1 ) (10.4) Критерий Грасгофа Gr: 1 g Gr 3 (10.43) Gr — критерий Грасгофа равен отношению подъемной силы, определяемой разностью плотностей в разных точках потока и 1, к силе внутреннего трения в неизотермической движущейся среде. В большинстве случаев Δ в критерии Грасгофа определяют как положительную разницу температур между ядром потока и стенкой. Поэтому, в условиях свободной конвекции стационарного процесса переноса теплоты, критериальное уравнение может быть записано так: Nu = f( Gr, Pr, Г1, Г. ) (10.44) Основным видом критериальных зависимостей, применяемых в инженерных расчетах, является степенная функция вида: Nu = ARe m Pr n Gr p Г q. (10.45) Коэффициенты: A, m, n, p, q получают, проводя экспериментальные исследования в определенных границах изменения параметров для определенной группы подобных явлений. Полученная зависимость применяется в инженерных расчетах для получения коэффициентов теплоотдачи α в указанном диапазоне величин. 11