Уравнение теплопроводности в полярных координатах

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

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

Для расчета распределения температуры по поверхности пластины используется программа ПР-1. Она содержит цикл по времени с вложенными в него двумя циклами по i и по j, в которых перебираются все элементы пластины и вычисляются их температуры в последующие моменты времени. Результат решения задачи представлен на рис. 1, — на нем разными цветами изображены области с различными температурами. Изотермы (границы разноцветных областей) перпендикулярны теплоизолированным краям пластины и параллельны краям, температура которых поддерживается постоянной.

Рис. 1. Распределение температуры: двумерная среда.

Задача 2.

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

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

Для решения задачи используется программа ПР-2. В ней последовательно перебираются узлы двумерной сетки по строкам и по столбцам. Когда переменная naprav равна 0, то элементы перебираются по столбцам сверху вниз и снизу вверх. Когда переменная naprav равна 1, то элементы перебираются по строкам слева направо и справа налево. Нижний край пластины теплоизолирован, это задается циклом:

Все остальные края пластины поддерживаются при постоянной температуре. Результат работы программы представлен на рис. 2.

Рис. 2. Теплопроводность в неоднородной среде.

Задача 3.

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

Для решения задачи может быть использован алгоритм АЛ-1.

В программе ПР-3 используются два массива T[i] и TT[i], в которых сохраняются значения температуры элементов стержня в моменты t и t+1 соответственно. Расчет температуры осуществляется по выведенной выше формуле в цикле. Будем считать, что длина стержня l, его коэффициент температуропроводности k при x>0,2l равен 1,8, а при x

Рис. 3. Теплопроводность стержня.

Задача 4.

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

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

В конечных разностях получаем:

Для решения этой задачи используется программа ПР-4. Он содержит два вложенных цикла по i и j, в которых перебираются все узлы двумерной сетки и пересчитываются значения T[i,j] на следующем временном слое. При ее запуске на экране появляется цветное изображение, границы одноцветных областей соответствуют изотермам. Пример результата вычислений приведен на рис. 4.

Рис. 4. Решение уравнения теплопроводности в полярных координатах.

Задача 5.

Металлический шар нагрели до температуры T₀, а потом опустили в вязкую жидкость (конвекция отсутствует), имеющую температуру T₁. Шар некоторое время охлаждается, температура его поверхности уменьшается. После этого шар вынули из жидкости и поместили в газообразную среду с низкой теплопроводностью. За счет теплообмена с центральной частью шара температура его поверхности повышается, а температура в центре понижается. Через некоторое время из–за теплообмена с газообразной средой температура поверхности начинает медленно уменьшаться, шар охлаждается. Необходимо рассчитать температуру различных точек шара и среды в последовательные моменты времени.

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

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

Используется программа ПР–5. Чтобы учесть граничные условия сопряжения, введем функцию q(t), учитывающую теплопередачу от шара к среде. Будем считать, что на поверхности шара (точка с координатой r=R) имеется поглотитель тепла с отрицательной мощностью α(T_п-T_с), а в ближайшей точке среды имеется источник тепла с той же по модулю мощностью. Здесь T_п и T_с –– температуры поверхности шара и прилегающего к ней слоя среды, α –– коэффициент теплоотдачи.

На рис. 5.1 показано начальное распределение температуры. Пока шар находится внутри жидкости, он охлаждается (рис. 5.2), температура поверхности становится ниже, чем в центре. Когда шар помещают в воздух, то температура поверхности повышается за счет теплообмена с центральной частью шара; теплообмен с окружающей средой происходит существенно медленнее (рис. 5.3). Через некоторое время температура поверхности шара начинает уменьшаться за счет теплообмена с окружающей средой (рис. 5.4). Если жидкость, в которую опускают нагретый шар, имеет небольшую вязкость, то возникают конвективные потоки. При этом можно приближенно считать, что температура всех точек жидкости одинакова и равна некоторой средней температуре T_ср. Решение задачи упрощается.

Рис. 5. Результаты расчета поля температур при охлаждении шара.

Задача 6.

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

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

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

В силу симметрии поле температур не зависит от φ. Запишем уравнение теплопроводности в цилиндрической системе координат и конечных разностях:

Для расчета нестационарного поля температур на следующем временном слое организуют два цикла по i и j, в которых рассчитываются температура в узлах сетки (процедуры R1 и R2) сначала в одном, а потом в другом направлениях. На экране рисуется половина осевого сечения цилиндрического сосуда с помощью точек, цвет которых зависит от температуры. Та часть вещества, которая находится в жидком состоянии (T больше T_K) изображается красным цветом. По мере охлаждения эта область постепенно уменьшается и исчезает.

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

For r:=1 to N do For j:=1 to M do begin if (j 15) then V:=V+3.1415*(r*r-(r-1)*(r-1))*dz; end;

Рис. 6. Образование углубления при отвердевании парафина.

Здесь учтено, что объем кольца с внешним радиусом r и внутренним r-Δr и высотой Δz равен 3,14*(r 2 -(r-Δr) 2 )Δz. Зная значение V в предыдущий момент времени t-1, можно найти объем вещества, затвердевшего в течение последнего шага по времени: ΔV=V t-1 -V t . Обозначим высоты верхней и нижней точек расплавленной области через h_a и h_b (рис. 6.1), тогда эффективная площадь горизонтального сечения S расплавленной области S=V/(h_a-h_b). Так как объем вещества за счет кристаллизации уменьшился на αΔV, то высота h_b уменьшится на Δh=αΔV/S, где α=0,02–0,05. Программа может нарисовать горизонтальную прямую AG, где G –– точка на высоте h_a, в которой жидкая фаза переходит в твердую.

Используется программа ПР–6; результаты расчета формы поверхности представлены на рис. 6.2. Видно, что в центральной части сосуда образуется углубление. Этот результат устойчиво получается и при небольших вариациях параметров системы. Можно усложнить модель и учесть, что коэффициент теплопроводности вещества в жидком состоянии существенно больше, чем в твердом. Например, в узлах сетки, для которых T больше T_K, коэффициент теплопроводности k=0,8, а в остальных узлах k=0,2. При этом получаются аналогичные результаты (рис. 6.2 и 7). Модель имеет несущественный недостаток: при решении уравнения теплопроводности мы пренебрегаем искривлением верхней границы расчетной области, считая ее горизонтальной плоскостью. Учет этого фактора приведет к изменению поля температур; изотермы останутся перпендикулярны верхней искривленной границе, решение задачи в центральной области будет более точным.

Рис. 7. Поле температур при остывании вещества в цилиндрическом сосуде.

Задача 7.

Рассмотрим одномерную двухфазную задачу Стефана. Имеется одномерная теплопроводящая среда, находящаяся в жидком и твердом состояниях. В начальный момент ее температура меньше температуры плавления T_пл. В точках с координатами x меньше 2 имеются источники тепла, под воздействием которых происходит плавление на границе раздела фаз, протекающее при температуре T_пл. Необходимо найти распределение температуры в последовательные моменты времени и определить расположение границы фаз.

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

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

где λ –– удельная теплота плавления, ρ –– плотность среды, x=ξ(t) –– уравнение границы разделяющей жидкую и твердую фазы, k₁, k₂ –– коэффициенты температуропроводности жидкой и твердой фазы соответственно, T₁ и T₂ –– температуры вблизи границы раздела x=ξ(t) со стороны жидкой и твердой фаз.

Рассмотрим решение задачи Стефана для одномерной среды длиной L, когда температура левого конца все время выше T_пл, правый конец теплоизолирован, а граница раздела фаз смещается вправо (рис. 8). На каждом шаге определяется число узлов j, для которых T_i больше T_пл (i не превосходит j). Температура в j–том узле незначительно отличается от температуры плавления T_пл, поэтому переменной T_j присваивается значение T_пл. Количество теплоты, приходящее в j-тый узел, на шаге t

Здесь k_j-1 и k_j+1 –– коэффициенты температуропроводности слева и справа от границы раздела двух фаз.

Рис. 8. Одномерная двухфазная задача Стефана.

Если количество теплоты dq t , поступающее в j-тый слой, меньше заданного значения теплоты плавления q_пл=λρΔx, которую необходимо затратить, чтобы расплавить слой толщиной h=Δx единичной площади, то фронт плавления не смещается. ЭВМ вычисляет температуру в 2, 3, . (j-1)-том, а затем в (j+1), (j+2), . (N-1)-ом узлах. Температура T_j остается равной T_пл; значение dq t накапливается в переменной dq. Если dq t превышает q_пл, то j–тый слой плавится, и граница раздела жидкой и твердой фаз смещается вправо на 1 узел. Величина dq t уменьшается на затраченную энергию q_пл. Переменной T_j+1 присваивается значение T_пл, и все повторяется снова.

Рассмотренный алгоритм реализован в программе ПР–7. В нашем случае теплопроводность жидкой фазы в раза больше, чем у твердой фазы. На рис. 9.1 показаны получающиеся распределения температуры в последовательные моменты времени, если температура левого конца одномерной среды непрерывно растет. На рис. 9.2 приведен результат, соответствующий ситуации, при которой температура левого конца стержня повышается до некоторой величины, а затем остается постоянной. В обоих случаях время плавления зависит от количества теплоты q_пл, необходимого для плавления слоя толщиной h=Δx единичной площади, взятого при температуре плавления.

Рис. 9. Моделирование плавления стержня.

Рассмотренные выше задачи в общем–то известны; методика их решения изложена в научной и учебной литературе. Некоторые из них обсуждались профессором В.А.Сараниным на конференции “Учебный физический эксперимент: Проблемы и решения”, проводимой в Глазовском педагогическом институте. Решение других задач можно найти в Интернете.

Тексты программ находятся в zip-архиве, файл gl-7.pas.

Уравнение теплопроводности

В данном разделе приводится вывод уравнений теплопроводности в виде дифференциальных уравнений в декартовых системах координат. Наиболее удобна дифференциальная форма уравнения теплопроводности. Тепловое уравнение изотропного материала. Рассмотрим бесконечно малое пространство с размерами Dx, boo, bg, iso -. Ферментация в 3-мерной системе координат x, y, R. 2-4.It также учитывается нестационарность, то есть изменение температуры в момент времени t.

По закону теплопроводности Фурье. тепло, поступающее по оси Х в основной объем, можно описать следующим образом: dQn =- Величину теплового тока вне получить, только первые 2 части ДК разобрать Рисунок 2-4.To вывод уравнения теплопроводности. Около объем вдоль оси x, который возможен с рядом Тейлора, сохраняется в качестве достаточного приближения： приращение теплового потока, вызванного теплопроводностью в направлении x, равно: =(2-7） Как и в формулах (от 2 до 7), 2 формулы в направлениях y и 2 могут быть записаны таким же образом. Сумма приращений теплового потока — это количество тепла, которое должно накопиться в объеме. [а-(Л£ -)+£(Я5 -)+е-( «>)] ’* .

Разновидностью теплотехники является теплоэнергетика. Другим из ответвлений общей теплотехники — строительная теплотехника. Людмила Фирмаль

Если количество теплоты Q ’(x, y, z, t) в единицу времени единичного объема происходит, то накопление теплоты в основном объеме происходит следующим образом: Ох (2-9） Тепло, остающееся в элементе объема по проводимости[уравнение (2-8)] и тепло, выделяемое самим объемом [уравнение (2-9)], увеличивают тепловую энергию элемента объема. За счет такого увеличения тепловой энергии изменяется теплоемкость объемного элемента, которую можно описать следующим образом: cp8×8z / 8з -, (2-10） Где c-удельная теплоемкость. p-это плотность. X-это время. Энергетический баланс объемного элемента может быть изменен путем уравнивания изменения количества тепла объемного элемента с тепловым потоком, вызванным теплом, генерируемым теплопроводностью и самим элементом.£ ) H (x£)+£(g£)+ s’ — 2- » >Вот что следует отметить А = ч(х, гг З. ы. т)\ с = р(х, у, Z, T)и р = р (х, г, З. ы. т).

Таким образом, (2-11) справедливо для изотропных гетерогенных сред. Если можно опустить термин, обозначающий тепло в объеме(нет источника в теле), (2-11), то в проекции на 3 координатные оси можно записать как: — B + 4-10b + b(2-12) эта формула более распространена и служит в разрезе анизотропного материала. Формула(2-11) может быть упрощена при применении к изотропным однородным материалам и когда теплопроводность считается постоянной. Подобный этому Х / ФА сложно? Он имеет размерность линейного измеренного значения в квадрате, деленного на время, и называется коэффициентом температуропроводности А. характеризует свойства материала.

Тепловое уравнение цилиндрической системы координат (с использованием hi Рис. 2-5.Цилиндрическая система координат. Рисунок 2-6.Сферическая система координат. Система координат (2-13) может быть описана в более удобной форме для цилиндрической системы. Таким образом, рисунок 2-5 х = rcos6; г / = г грех б; з = з \ ДТ udChₜ1 | 1(Д1/\, м ’/ₙ1 + Р «ПФП + ДЗ) ПК (21⁴） Уравнение теплопроводности для сферической системы координат. Аналогичное преобразование сферической системы(рис. 2-6) приводит к следующей формуле: х = г sinφcos п; г = р sinφsin п; г = gsozf; ■1 Д /(&Г1d \ РТ). д / Я 1d2/ 1, м ’ ДТ [гр гду **. ГСИН у у ’ duJ ’g2sin2u ’ ПК’ (2-15） Уравнение теплопроводности анизотропного material. In в предыдущем разделе получено уравнение теплопроводности для изотропных сред.

Некоторые технологически продвинутые слоистые материалы имеют теплопроводность 6Т и существенно изменяются в зависимости от направления теплового потока, проходящего через корпус. Материалы этой категории включают: кристаллический материал, древесину, составные плиты и металл используемые для анкеров в трансформаторах, и plywood. To примените уравнения теплопроводности к этим анизотропным материалам, они должны быть соответствующим образом изменены. Обычная форма этого уравнения очень сложна[L. 5], поэтому она не рассматривается в этой книге. Но、 Рисунок 2-7.Тепловой поток в анизотропной среде. Для 2-мерных измерений учитываются основные понятия. В случае 2-мерных измерений теплопроводность распределяется так, что максимальные и минимальные значения распределяются так, что они называются»приоритетной»осью или»главным валом».

Величина теплопроводности в другом направлении имеет промежуточное значение. Распределение теплопроводности может быть представлено эллипсом с осями, соответствующими максимальному и минимальному значениям теплопроводности. представьте себе объект в системе координат (x, y) (рис.2-7). он образует угол 0 с главной осью теплопроводности материала. Система координат ( £ , m]) совпадает с главной осью теплопроводности. Тепло течет по телу в направлении координат|и Т）： −1. dt. Поток в сторону хны выглядит так: (2-16） Градиенты температуры могут быть преобразованы в градиенты X и y с помощью следующих соотношений: ДТ dt ДХ ИДТ сделать df «dhdG «’ duod ^’ ДТ _ _ _ _ д / ДХ. д / д Дри ДХ ДТ] ’у д（\ И в соответствии с геометрией фигуры Р / = грешить£= — qcosp; х = $ КГУ = —qsinp.

Если подставить эти значения в Формулу (2-16) и преобразовать ее, то получится тепловой поток. Ях = — потому что₽+ грехах ’ п)^ — (iₜ— потому что грех п п; Общая форма теплового потока [уравнение (2-12)] уравнение теплопроводности [уравнение (2-17)] может быть использовано для описания двумерного уравнения теплопроводности анизотропного материала в следующем виде: 57 = \ COS в ’п + грехах’ п） +(sinsinр+ x, cosp)g +(я£-X₄)sin2pd.(2

18） Для изотропных сред L.= LT. и p =0.In при этих условиях уравнение(2-18) сводится к 2-мерному уравнению(2-13). Если пластины из анизотропного материала зажаты между изотермическими поверхностями испытуемой теплопроводной системы, и Образец готовят таким образом, что его главная ось образует изотермическую поверхность и угол р, после чего определяют измеренную теплопроводность (в зависимости от того, была ли она выполнена вдоль направления x или y). Xia3 как COS Ильин; = грех₽ — ч cos3.

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

Если геометрическая ось анизотропного тела совпадает с главной осью теплопроводности, то уравнение (2-18） Упрощенный: д / ч-л ч мск Для древесины с различной теплопроводностью по окружности 0, вдоль волокна z, вдоль волокна g идея формулы(2-19) может быть применена следующим образом: Рисунок 2-8.Вектор теплового потока слоистого материала. Уравнение(2-14) помещает ось z вдоль линии центра дерева, игнорируя термин, обозначающий источник тепла. Поэтому в цилиндрических координатах (2-20） Пример 2-1.Пластина из слоистого материала (пластина) используется для экспериментов по теплопроводности. Слой под углом с гладкой поверхности образца(рис. 2-8).

Поверхность A является изотермической поверхностью для поддержания постоянной, но разной temperature. It необходимо рассчитать угол, который складывается за счет тепла N потоков, перпендикулярных изотермической плоскости. ЮЖД рис. 2-8, С ^ С ^ А ДТ / ДТС тгт〜\ ДТ / дл■ Но основываясь на предыдущих расчетах dt. В. С.?^⁺сентенция ДТ. Сделать это. ДТ, ДТ. «Чо?»ДТ⁺ы? В. Отсюда ATⱼ[соз£(ДТ [ды) — грех п(dtjdx)] у- ^ [что? (dtfdx) — грех п (dtjdy）] * (ля） Однако в системе, которую мы рассмотрим здесь, интерфейс остается изотермическим. То есть координатная ось y находится в направлении перпендикуляра n, а x — в изотермическом plane. So, d / / dx = 0, Формула (a)принимает вид: тгр = ^ — Ctgf = ^ ТГА. Таким образом, y a, А Вектор теплового потока не перпендикулярен изотермической плоскости, как в случае, когда материал является изотермическим. Если X ^ * = = 2X ^ среда-это дерево и 0 = 45°、 1 tgT = » 2 -> Задачи 2-1.

Предложен способ измерения теплопроводности жидкого металла при высоких температурах temperatures. It схематически показана основная часть устройства и определена погрешность предлагаемых измерений. 2-2.Используя цилиндрическую систему координат, разложим уравнение (2-11) по модели уравнения (2-14). 2-3.Используя сферическую систему координат и малые объемные элементы этой системы координат, мы расширяем уравнение (2-15) в соответствии с моделью уравнения (2-11). 2-4. Создайте тепловое уравнение для трехмерной анизотропной среды.

Образовательный сайт для студентов и школьников

Копирование материалов сайта возможно только с указанием активной ссылки «www.lfirmal.com» в качестве источника.

Дифференциальное уравнение теплопроводности

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

(1.12)

где с_р, Дж/(кг×К) – изобарная теплоемкость; r, кг/м 3 – плотность; l, Вт/(м×К) – коэффициент теплопроводности; w_х, w_y, w_z – проекции вектора скорости движения жидкости; q_v , Вт/м 3 – объемная плотность внутреннего тепловыделения жидкости.

Уравнение (1.12) записано для случая l=const.

Дифференциальное уравнение температурного поля для твердых тел называется дифференциальным уравнением теплопроводности и может быть получено из (1.12) при условии w_х= w_y= w_z=0, с_р= с_v=с:

где — коэффициент температуропроводности, характеризует скорость изменения температуры в теле. Значения а = f (t) для различных тел приводятся в справочниках.

Дифференциальное уравнение теплопроводности

(1.13)

описывает нестационарное температурное поле твердых тел с внутренним тепловыделением (с внутренними источниками тепла). Такими источниками тепла могут быть: джоулева теплота, выделяемая при прохождении электрического тока по проводникам; теплота, выделяемая ТВЭЛами ядерных реакторов и т.д.

Дифференциальное уравнение теплопроводности (1.13), записанное в декартовых координатах, можно представить в цилиндрических (r, z, φ) и сферических (r, φ, ψ).

В частности, в цилиндрических координатах (r –радиус; φ – полярный угол; z — аппликата) дифференциальное уравнение теплопроводности имеет вид

(1.14)

Условия однозначности

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

· геометрические условия, характеризующие форму и размеры тела;

· физические условия, характеризующие свойства участвующих в теплообмене тел;

· граничные условия, характеризующие условия протекания процесса на границе тела;

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

При решении задач теплопроводности различают:

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

· граничные условия второго рода, когда задается плотность теплового потока на поверхности тела:

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

В соответствии с законом Ньютона-Рихмана тепловой поток, передаваемый с 1м 2 поверхности в среду с температурой t_ж,

В то же время этот тепловой поток подводится к 1м 2 поверхности из глубинных слоев тела теплопроводностью

Тогда уравнение теплового баланса для поверхности тела запишется в виде

(1.15)

Уравнение (1.15) является математической формулировкой граничных условий третьего рода.

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

Контрольные вопросы и задания

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

2. Почему в правой части уравнения (1.3) стоит минус?

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

4. Как определяется тепловой поток (Q, Вт) при конвективной теплоотдаче, теплопроводности, тепловом излучении?

5. Запишите дифференциальное уравнение теплопроводности в декартовых координатах, описывающее трехмерное стационарное температурное поле без внутренних источников теплоты.

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

2. ТЕПЛОПРОВОДНОСТЬ И ТЕПЛОПЕРЕДАЧА
ПРИ СТАЦИОНАРНОМ РЕЖИМЕ

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

Дано:плоская однородная стенка толщиной δ (рис. 2.1) с постоянным коэффициентом теплопроводности λ и постоянными температурами t₁ и t₂ на поверхностях.

Определить:уравнение температурного поля t=f (x) и плотность теплового потока q, Вт/м 2 .

Температурное поле стенки описывается дифференциальным уравнением теплопроводности (1.3) при следующих условиях:

· т. к. режим стационарный;

· т.к. отсутствуют внутренние источники теплоты;

· т.к. температуры t₁ и t₂ на поверхностях стенки постоянны.

Температура стенки является функцией только одной координаты х и уравнение (1.13) принимает вид

(2.1)

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

Граничные условия первого рода:

при х=0 t= t₁ ,	(2.2)
при х= δ t= t₂.	(2.3)

Выражения (2.1), (2.2), (2.3) являются математической постановкой задачи, решение которой позволит получить искомое уравнение температурного поля t= f (x).

Интегрирование уравнения (2.1) дает

При повторном интегрировании получим решение дифференциального уравнения в виде

t=с₁х+с₂.

(2.4)

Из уравнения (2.4) при условии (2.2) получим

а при условии (2.3)

Подстановка констант интегрирования с₁ и с₂ в уравнение (2.4) дает уравнение температурного поля

(2.5)

по которому можно рассчитать температуру по толщине стенки на любой координате 0 2 ; t₂, t₃.

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

	(2.8)
	(2.9)
	(2.10)

(2.11)

Сложив левые и правые части уравнений (2.11), получим расчетную формулу для плотности теплового потока, передаваемого через трехслойную стенку

(2.12)

Температуры на границах слоев t₂ и t₃ можно рассчитать по уравнениям (2.8) – (2.10) после того, как найдена плотность теплового потока (q) по (2.12).

Общий вид уравнения (2.12) для многослойной плоской стенки, состоящей из п однородных слоев с постоянными температурами на наружных поверхностях и , имеет вид

(2.13)

Средний коэффициент теплопроводности многослойной стенки называют эффективным (λ_эф). Он равен коэффициенту теплопроводности однородной стенки, толщина и термическое сопротивление которой равны толщине и термическому сопротивлению многослойной стенки

(2.14)

2.2. Теплопроводность цилиндрической стенки
при граничных условиях первого рода

Дано:Однородная цилиндрическая стенка (стенка трубы) с внутренним радиусом r₁, наружным – r₂, длиной , с постоянным коэффициентом теплопроводности λ, с постоянными температурами на поверхностях t₁ и t₂.
(рис. 2.3).

Определить: уравнение температурного поля
t = f (r), тепловой поток, передаваемый через стенку
Q, Вт.

Дифференциальное уравнение теплопроводности в цилиндрических координатах (1.14) для условий данной задачи:

(2.15)

Граничные условия первого рода:

при r=r₁ t=t₁ ,	(2.16)
при r=r₂ t=t₂ .	(2.17)

Порядок решения системы уравнений (2.15) – (2.17) тот же, что и в случае плоской стенки: находится общий интеграл дифференциального уравнения второго порядка (2.15), который содержит две константы интегрирования
с₁ и с₂ . Последние определяются с помощью граничных условий (2.16) и (2.17) и после подстановки их значений в решение дифференциального уравнения (общий интеграл) получаем уравнение температурного поля цилиндрической стенки t = f (r) в виде

(2.18)

где r₁ r r₂ – текущий радиус.

Нетрудно убедиться, что при подстановке в (2.18) r= r₁ получим t=t₁ , при r=r₂ получим t=t₂. Распределение температуры по толщине цилиндрической стенки, в соответствии с (2.18) подчиняется логарифмическому закону (рис. 2.3).

Для определения теплового потока воспользуемся законом Фурье:

(2.19)

Если взять производную от правой части уравнения (2.18) и подставить в (2.19), получим расчетную формулу для теплового потока цилиндрической стенки

(2.20)

В технических расчетах часто тепловой поток вычисляется для 1 м длины трубы:

и называется линейной плотностью теплового потока.

Запишем уравнение (2.20) в виде

где – термическое сопротивление теплопроводности цилиндрической стенки.

Для трехслойной цилиндрической стенки (трубы, покрытой двумя слоями тепловой изоляции) с известными постоянными температурами поверхностей (t₁ и t₄), с известными геометрическими размерами (r₁ , r₂, r₃, r₄ , ) и коэффициентами теплопроводности слоев (λ₁, λ₂, λ₃) (рис. 2.4) можно записать следующие уравнения для теплового потока Q:

(2.21)

Совместное решение системы уравнений (2.21) дает расчетную формулу для теплового потока, передаваемого через трехслойную стенку при заданных температурах на поверхностях,

(2.22)

Температуры на границах слоев (t₂, t₃) можно рассчитать по уравнениям (2.21).

Для многослойной цилиндрической стенки, состоящей из п слоев, формулу (2.22) можно записать в общем виде

(2.23)

Эффективный коэффициент теплопроводности для многослойной цилиндрической стенки, как и для многослойной плоской стенки, определяется из равенства суммы термических сопротивлений многослойной стенки термическому сопротивлению однородной стенки той же толщины, что и многослойная. Так, для двухслойной тепловой изоляции трубы
(рис. 2.4) эффективный коэффициент теплопроводности (λ_эф) определ ится из равенства

2.3. Теплопроводность плоской и цилиндрической стенок
при граничных условиях третьего рода (теплопередача)

Граничные условия третьего рода состоят в задании температуры жидкости (t_ж) и коэффициента теплоотдачи ( ) между поверхностью стенки и жидкостью.

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

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

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

Рассмотрим теплопередачу через плоскую и цилиндрическую стенки при условии только конвективного теплообмена на поверхностях. Теплопередача с радиационно-конвективным теплообменом (сложным теплообменом) на поверхностях будет рассмотрена позже.

Плоская стенка(рис. 2.5)

Дано:

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

– от горячей жидкости к стенке

(2.24)

(2.25)

– от стенки к холодной жидкости

(2.26)

Записав уравнения (2.24) – (2.26) в виде

(2.27)

и сложив почленно правые и левые части уравнений (2.27), получим формулу для расчета теплопередачи (q, Вт/м 2 ) через плоскую стенку в виде

(2.28)

Величины называются термическими сопротивлениями теплоотдачи. Они прямо пропорциональны перепадам температур .

Температуры на поверхностях стенки t₁ и t₂ можно рассчитать по уравнениям (2.24) – (2.26) после того, как определена плотность теплового потока (q) по уравнению (2.28).

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

(2.29)

где — коэффициент теплопередачи плоской стенки,характеризует интенсивность процесса теплопередачи.

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

(2.30)

Цилиндрическая стенка(рис. 2.6)

Дано:

Для цилиндрической стенки, по аналогии с плоской стенкой, можно записать следующую систему уравнений:

	(2.31)
	(2.32)
	(2.33)

где — площади внутренней и наружной поверхностей трубы.

Записав уравнения (2.31) – (2.33) относительно разностей температур, а затем сложив правые и левые части уравнений, получим формулу для расчета теплопередачи (Q, Вт) через цилиндрическую стенку в виде

(2.34)

Температуры на поверхностях стенки t₁ и t₂ рассчитываются по уравнениям (2.31) – (2.33).

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

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

Для металлических труб с можно пренебречь кривизной стенки и теплопередачу рассчитать по формулам для плоской стенки:

источники:

http://lfirmal.com/uravnenie-teploprovodnosti/

http://megaobuchalka.ru/8/46137.html