Слайд 1Теория тепловых процессов
Лекция №1
«Основные положения учения о теплопроводности»
Слайд 2Теплопередача или теплообмен — учение о самопроизвольных необратимых процессах распространения теплоты
в пространстве.
Слайд 3Методы изучения физических явлений
феноменологический
статистический
Плюсы метода:
-позволяет сразу установить общие связи между параметрами,
характеризирующими процесс;
-использовать экспериментальные данные, точность которых предопределяет и точность самого метода.
Минусы:
-сложный;
-дорогостоящий.
Плюсы метода:
-позволяет получить феноменологические соотношения на основании заданных свойств микроскопической структуры среды без дополнительного эксперимента.
Минусы:
-сложность;
-требуются знания ряда параметров, для которых требуются знания специальных разделов физики.
Слайд 4Температурное поле
Теплопроводность обусловлена движением микрочастиц вещества.
Математическое выражение температурного поля:
Слайд 5Уравнение для одномерного поля:
Простой вид уравнения одномерного стационарного температурного поля:
Слайд 6Температурный градиент
Изотермической поверхностью называется геометрическое место точек в температурном поле, имеющих
одинаковую температуру.
Изотермы
Возрастание температуры в направлении нормали к изотермической поверхности характеризуется градиентом температуры.
Градиент температуры:
Слайд 7
Проекции вектора grad t на оси Ox, Oy, Oz:
Слайд 8Тепловой поток.
Закон Фурье.
Для передачи теплоты теплопроводностью необходимо неравенство нулю температурного
градиента в различных точках тела.
Гипотеза Фурье:
Плотностью теплового потока:
Слайд 9
Изотермы и линии
теплового тока
dF-элемент изотермической поверхности.
Слайд 11Коэффициент теплопроводности
Порядок значений коэффициентов теплопроводности различных веществ
Коэффициент теплопроводности,
Вт/(м*К):
Слайд 12Коэффициент
теплопроводности газов
Слайд 13Коэффициент теплопроводности жидкости
Коэффициент теплопроводности капельных жидкостей лежит примерно в пределах от
0,07 до 0,7 Вт/(м*К).
При повышении давления коэффициенты теплопроводности жидкостей возрастают.
Слайд 14Коэффициент теплопроводности твердых тел
Слайд 16Дифференциальное уравнение теплопроводности
Слайд 17Если ограничиться двумя первыми членами ряда, то:
Количество теплоты, отведенное через
противоположную грань элементарного параллелепипеда в направлении оси Ох, запишется как:
Разница количеств теплоты:
или
Количество теплоты dQ:
Т. о.
Функция является непрерывной в рассматриваемом интервале dx и может быть разложена в ряд Тейлора:
Слайд 18Если система тел не содержит внутренних источников тепла (q=0), тогда выражение
принимает форму уравнения Фурье:
Дифференциальное уравнение теплопроводности превращается в уравнение- Пуассона:
Для стационарной теплопроводности и отсутствия внутренних источников теплоты выражение принимает вид уравнения Лапласа:
Слайд 19Лекция №2
«Теплопроводность при стационарном режиме»
Слайд 20Передача теплоты
через плоскую стенку ( qv=0 )
или
если qv=0, то:
Слайд 21Граничные условия первого рода.
Если ось Ох направить, как показано на
рис. 2-1, то температура в направлении осей Оу и Оz будет оставаться постоянной:
Граничные условия задаются в виде:
Слайд 22После второго интегрирования:
Постоянные С1 и С2 определяются из:
Слайд 23Закон распределения температуры:
- текущий температурный напор
- полный температурный напор
- безразмерный температурный
напор
- безразмерная координата
Получим:
Слайд 24Согласно:
Получаем:
Количество теплоты:
В случае когда:
Слайд 27Граничные условия третьего рода (теплопередача).
Плотность теплового потока от горячей жидкости
к стенке определяется уравнением
При стационарном тепловом режиме тот же тепловой поток пройдет путем теплопроводности через твердую стенку:
Тот же тепловой поток передается от второй поверхности стенки к холодной жидкости за счет теплоотдачи:
Таким образом получим:
Слайд 28Сложив равенства получим:
Отсюда плотность теплового потока:
Обозначим:
Тогда:
Полное термическое сопротивление однослойной стенки :
Слайд 29Если стенка состоит из слоев, то полное термическое сопротивление теплопередачи через
такую стенку будет равно:
Отсюда
Плотность теплового потока через многослойную стенку, состоящую из n слоев, будет равна:
Тепловой поток Q, Вт, через поверхность F твердой стенки:
Слайд 30Температура при граничных условиях третьего рода:
Слайд 31Графический метод определения температур на поверхностях слоев неоднородной стенки
Общее термическое сопротивление
теплопередачи через стенку равно:
Из подобия треугольников C1B1B2 и C1C2E1 следует:
Т. к.
тогда
Слайд 32Граничные условия второго и третьего рода
В силу стационарности теплового режима можно
записать следующие уравнения:
на внешней правой поверхности
на внешней левой поверхности
на поверхности между слоями m — 1 и m
Слайд 33Передача теплоты через цилиндрическую стенку
1) Граничные условия первого рода
Цилиндрическая система
координат:
Т. к. поле одномерное, то:
Изотермические поверхности являются цилиндрическим, сл-но:
Граничные условия:
Слайд 34Введем новую переменную:
Тогда:
Интегрируем и получаем:
Граничные условия:
Слайд 35Решение уравнения:
или
Согласно закону Фурье:
Т. к.
то
Слайд 36Тепловой поток, отнесенный к единице длины трубы, измеряется в Вт/м и
называется линейной плотностью теплового потока.
Среднеинтегральное значение:
Уравнение закона Фурье:
Слайд 372) Граничные условия третьего рода (теплопередача)
Температурный напор:
Тогда:
Слайд 38Величина ki называется линейным коэффициентом теплопередачи, он измеряется в Вт/ (м
• К).
Линейное термическое сопротивление теплопередачи:
Плотность теплового потока:
где
Слайд 39На практике часто встречаются цилиндры, толщина стенок которых мала по сравнению
с диаметром. В этом случае при расчетах можно пользоваться упрощенными формулами Для получения таких формуй поступим следующим образом.
Если
то
Если стенка тонкая:
Для многих технических расчетов ошибка, не превышающая 4%, вполне допустима Обычно в инженерных расчетах при d2/d1≤1,8 используются формулой.
Ошибку можно уменьшить, если в качестве расчетной поверхности брать поверхность, со стороны которой α меньше:
Слайд 40В случае теплопередачи через многослойную цилиндрическую стенку система равенств должна быть
заменена системой, учитывающей сопротивление теплопроводности всех слоев:
Слайд 41В случае задания граничных условий первого рода можно рассматривать как предельный
случай граничных условий третьего рода:
Температура между граничными слоями:
Слайд 42Пути интенсификации теплопередачи
а) Интенсификация теплопередачи путем увеличения коэффициентов теплоотдачи
Из уравнения теплопередачи
Для
плоской стенки:
Слайд 43Из уравнения видно, что коэффициент теплопередачи не может быть больше α.
Проследим это на числовых примерах.
Слайд 44б) Интенсификация теплопередачи за счет оребрения стенок
При передаче теплоты через цилиндрическую
стенку термические сопротивления определяются не только значениями коэффициентов теплоотдачи, но и размерами самих поверхностей. При передаче тепла через шаровую стенку влияние диаметров оказывается еще сильнее. Отсюда следует, что если α мало, то термическое сопротивление теплоотдачи можно уменьшить путем увеличения соответствующей поверхности. Такой же результат можно получить и для плоской стенки, если одну из поверхностей увеличить путем оребрения. Последнее обстоятельство и положено в основу интенсификации теплопередачи за счет оребрения. При этом термические сопротивления станут пропорциональными величинам:
Строгое аналитическое решение задачи о распространении тепла в ребре связано со значительными трудностями. В основу решения поэтому кладут некоторые допущения, которые позволяют сравнительно простым путем получить нужный результат. Ниже рассмотрим метод решения задач о теплопроводности в ребрах простейших геометрических форм.
Слайд 45Теплопроводность круглого ребра постоянной толщины
Рассмотрим расчет теплопроводности круглого ребра постоянной толщины
Избыточная температура ребра будет:
Составим уравнение баланса энергии для кольцевого элемента ребра
получаем дифференциальное уравнение вида:
Уравнение представляет собой уравнение Бесселя, имеющее общее решение вида
Свойства функции:
Слайд 46Если температурой с торца круглого ребра пренебречь:
для температуры на конце ребра
для
количества теплоты
Количество теплоты, которое будет отдаваться поверхностью круглого ребра постоянной толщины
Слайд 47Лекция № 3
Основные положения учения о конвективном теплообмене.
Слайд 48Основные понятия и определения.
Конвекция
Свободная конвекция.
движение в рассматриваемом объеме жидкости возникает за
счет неоднородности в нем массовых сил.
Вынужденная конвекция.
Вынужденное движение рассматриваемого объема жидкости происходит под действием внешних поверхностных сил, приложенных на его границах за счет предварительно сообщенной кинетической энергии .
Слайд 49Физические свойства жидкости.
Процесс теплообмена зависит от физических свойств жидкостей:
коэффициент теплопроводности Л,
удельная теплоемкость ср,
плотность р,
коэффициент температуропроводности а
коэффициент вязкости μ .
Слайд 50Реальные жидкости обладают вязкостью; между частицами или слоями, движущимися с различными
скоростями, всегда возникает сила внутреннего трения, противодействующая движению.
Согласно закону Ньютона эта касательная сила s, Па
Коэффициент μ называется динамическим коэффициентом вязкости или просто коэффициентом вязкости; его единица измерения Н-с/м2. При
Численно
Слайд 51отношение вязкости μ к плотности р, называемое кинематическим коэффициентом вязкости и
обозначаемое буквой v, м2/с:
У капельных жидкостей вязкость почти не зависит от давления, но значительно уменьшается при повышении температуры.
Слайд 52
На теплоотдачу оказывает влияние сжимаемость жидкостей. Изотермической сжимаемостью или коэффициентом сжатия
тела при
t=const называют величину
Для капельных жидкостей изотермическая сжимаемость чрезвычайно мала. Для воды
Слайд 53Помимо изотермической сжимаемости для конвективного теплообмена большое значение имеет тепловое расширение
жидкости, характеризующееся температурным коэффициентом объемного расширения, определяемым уравнением (р=const)
Для некоторых жидкостей, например для воды при t<4°С, коэффициент р может иметь отрицательное значение.
Для идеального газа температурный коэффициент объемного расширения есть величина, обратная абсолютной температуре газа,
Слайд 54Дифференциальные уравнения конвективного теплообмена.
Из уравнения:
следует, что плотность теплового потока в
любой точке жидкости для каждого момента времени однозначно определяется, если известны поля температур, удельной энтальпии и скорости.
Слайд 55Для реальной жидкости i=i(T, р), и согласно понятию о полном дифференциале
Отсюда
Слайд 56в предположении о несжимаемости жидкости. (р=const) с достаточной степенью точности можно
принять (di/dp)T—0 т. е. пользоваться соотношением, справедливым для термодинамически идеального газа
Слайд 57Вывод уравнения энергии
Где
Согласно уравнению проекции плотности теплового потока q на
координатные оси Ох, Оу и Oz равны:
Слайд 58Подставляя значений qx, qy и qz в уравнение (1-25), можно получить:
Для
несжимаемых жидкостей p=const
Слайд 59Тогда
или, если
если t=t(т, х, у, z)t то на основании понятия
о полной производной имеем:
где
Слайд 60Здесь dt/dx характеризует изменение температуры во времени в какой-либо точке жидкости,
т. е. является локальным изменением t;
Применяя обозначение
Уравнение энергии можно записать
если
Уравнение энергии переходит в уравнение теплопроводности
Слайд 61При стационарных процессах конвективного теплообмена dt/dт=0.
упрощенный вывод этого уравнения для
случая одномерного течения несжимаемой жидкости.
Выделим в потоке вязкой жидкости элементарный объем с размерами ребер dx, dy и dz
Вывод уравнения движения основан на втором закона Ньютона: сила равна массе, умноженной на ускорение
Слайд 62Силы, действующие на рассматриваемый элемент жидкости, можно разделить на массовые (или
объемные) и поверхностные
проекция на ось Ох равна произведению проекции ускорения свободного падения gx на массу элемента:
Если на верхней грани элемента давление жидкости равно р, то на площадку dydz действует сила р dy dz.
На нижней грани давление с точностью до второго члена разложения в ряд Тейлора равно
и на эту грань действует сила
Слайд 63Равнодействующая этих сил равна алгебраической сумме:
Подставляя
получаем
Согласно второму закону механики эта равнодействующая
равна произведению массы элемента на его ускорение dwx/dt и учитывает силы инерции (б):
Суммируя dfi, df2 и dfs, получаем проекцию на ось Ох равнодействующей всех сил, приложенных к объему (а):
Слайд 64Приравнивая правые части уравнений (а) и (б) и производя сокращения, окончательно
имеем уравнение движения вдоль оси Ох:
трехмерного движения несжимаемой жидкости с постоянными физическими параметрами скоростное поле описывается тремя уравнениями движения, каждое соответственно в проекциях сил на оси Ох, Оу и Oz: для оси Ох
Слайд 65ДЛЯ оси ОY
Для оси OZ
На основании понятия о полной производной имеем:
Аналогично
и для других осей:
Слайд 66Лекция № 4
Теплоотдача при вынужденном поперечном омывании труб и пучков труб.
Слайд 67Основные положения
Омывание трубы поперечным неограниченным потоком жидкости характеризуется рядом особенностей. Плавное,
безотрывное обтекание цилиндра
Слайд 68Теплоотдача при поперечном омывании пучков труб
Слайд 69Если учитывается только сила тяжести, то
Свободное движение возникает за счет неоднородного
распределения в рассматриваемой жидкости массовых (объемных) сил.
Такими силами являются:
Сила тяжести;
Центробежная сила;
Сила за счет наведения в жидкости электромагнитного поля высокой напряженности.
Слайд 70Теплоотдача при свободном движении жидкости в большом объеме.
Теплоотдача при свободном
ламинарном движении вдоль вертикальной пластины
Вертикальная пластина с неизменной температурой поверхности равной tc , находится в жидкости или газе. Температура жидкости вдали от пластины постоянна и равна t0.
tc > t0
Слайд 71Будем полагать, что температура в движущемся слое жидкости изменяется по уравнению:
Слайд 72Для стационарного течения я с учетом ранее принятых допущений уравнение движения
упрощается, В результате будем иметь:
При линейной зависимости плотности от температуры, где
Отсюда
Слайд 73Подставляя значение υ
в уравнение и учитывая последнее соотношение для плотности,
уравнение движения можно написать следующим образом:
Интегрирование уравнения движения дает:
Слайд 74Примем следующие граничные условия для скорости: wx=0 как при у=0, так
и при у= δ. Отметим, чго, строго говоря, при у= δ(υ =0) скорость может быть не равна нулю. Из принятых условий следует:
Подставив значение C1 и C2 в уравнение и преобразовав получим уравнение распределения скоростей в движущемся слое жидкости:
Максимум скорости соответствует значению:
Слайд 75Среднеинтегральная скорость равна:
среднюю температуру жидкости в слое определим приближенно как среднеинтегральную
по сечению слоя
Слайд 76Подставляя значение в уравнение, получаем:
Расход жидкости
через поперечное сечение слоя равен:
Слайд 77Можно считать, что в среднем жидкость нагревается до температуры
. На этот нагрев затрачивается теплота
Из уравнения (е) следует, что
Приравнивая правые части уравнений , получаем дифференциальное уравнение, описывающее изменение δ по высоте стенки:
Слайд 78Постоянную интегрирования с найдем из условия, что при х=0 δ =0.
Отсюда с=0.
Из уравнения следует, что
Согласно уравнению α=2λ/δ
Подставляя сюда значение δ, получаем:
Слайд 79Приведем уравнение к безразмерному виду, для чего левую и правую части
уравнения умножим на х и разделим на λ.
После некоторых преобразований получим:
Из уравнения следует, что .При этом
Слайд 80Для расчета местных коэффициентов при свободном ламинарном течении вдоль вертикальных стенок
можно использовать формулу:
при условии, что
Тогда расчетная формула для средних коэффициентов теплоотдачи будет:
Слайд 81Теплоотдача при свободном турбулентном движении вдоль вертикальной пластины
развитое турбулентное течение
наступает при числах
.
Для местных коэффициентов теплоотдачи при развитом турбулентном течении в предложена формула:
Линейный размер входит в числа Nu и Gr :
Слайд 82Теплоотдача при переходном режиме свободного движения вдоль вертикальной пластины
переходный режим имеет
место примерно при
Переходный режим отличается неустойчивостью процесса течения и теплоотдачи и, как следствие, большим разбросом опытных точек.
Слайд 83Теплоотдача при свободном движении около горизонтальной трубы
Слайд 84Теплоотдача при очень малых значениях комплекса GrPr
Слайд 85Лекция № 5
Основные законы излучения.
Слайд 86Описание процесса. виды лучистых потоков, вектор излучения.
Тепловое излучение представляет собой
процесс распространения внутренней энергии излучающего тела путем электромагнитных волн.
Слайд 87Виды лучистых потока: вектор излучения.
Поверхностное (полусферическое) излучения.
Тело излучает энергию в виде
непрерывного(сплошного)
Или прерывистого спектра по длинам волн.
Интегральный поток, испускаемый с единицы поверхности, носит название поверхностной плотности потока интегрального излучения
Лучистый поток со всей поверхности выразится интегралом:
Слайд 88Отношение плотности лучистого потока, испускаемого в бесконечном малом интервале длин волн,
к величине этого интервала длин волн называется спектральной плотностью потока излучения:
Количество энергии, испускаемое в определенном направлении /, определяемым углом -ψ с нормалью к поверхности п единицей элементарной площадки в единицу времени в пределах элементарного телесного угла dω, называется угловой плотностью излучения.
Слайд 89Из этих соотношений следует, что
Интенсивностью излучения называется количество лучистой энергии, испускаемое
в направлении угла - ψ в единицу времени элементарной площадкой в пределах единичного элементарного телесного угла, отнесенное к проекции этой площадки на плоскость, ортогональную к направлению излучения
Слайд 90Потоки интегрального и монохроматического излучения связаны следующими зависимостями:
Плотность потока поглощающей лучистой
энергии
Часть падающей энергии, которую поверхность данного тела отражает обратно окружающим его телам, носит название потока отраженного излучения. Плотность потока отраженного излучения E отр, Вт/м2, равна:
Слайд 91Часть падающей энергии излучения, проходящая сквозь тело, называется плотностью потока пропускаемого
излучения Е проп, Вт/м2:
Суммарная величина плотностей потоков собственного и отраженного излучения, испускаемого поверхностью данного тела, называется плотностью эффективного излучения:
Слайд 92В первом способе (Нуссельта) поток результирующего излучения определяется из теплового баланса
относительно поверхности а-а расположенной внутри тела вблизи его поверхности с учетом, что Е проп = 0
Второй способ (О. Е. Власова) состоит в определении фрез из баланса относительно воображаемой поверхности б— б, находящейся вне тела, но вблизи его поверхности
В общем случае плотность потока результирующего излучения определяется разностью .встречных потоков излучения, падающих на условную поверхность в — в (Ю. А. Суринов)
Слайд 93если выразить Е эф из способа Власова, а Е пад из
способа Нуссельта
и подставить в способ Суринова, то получим:
Слайд 94Объемное излучение
Для среды, которая заполняет некоторый объем системы и может быть
излучающей, поглощающей и рассеивающей, характерными являются объемные плотности потоков излучения
Интегральной и монохроматической объемными плотностями потоков собственного излучения называются лучистые потоки, испускаемые единицей объема среды в единицу времени по всем различным направлениям в пределах пространственного угла ω=4π:
Слайд 95Лучистые потоки, отнесенные ко всему объему, выразятся зависимостями
Потоки монохроматического и интегрального
излучений связаны соотношением
Плотность поглощенного объемного излучения
и плотность рассеянного объемного излучения
Слайд 96плотностью эффективного объемного излучения называется суммарная величина плотностей потоков собственного и
рассеянного излучений
плотность потока результирующего объемного излучения выражается зависимостью
Последние две зависимости могут быть использованы для получения уравнений, связывающих плотности потоков результирующего и эффективного объемных излучений, аналогичных
Слайд 97Вектор излучения.
Вектор излучения определяет направление наиболее интенсивного переноса лучистой энергии в
рассматриваемой точке поля излучения.
, получим:
Составляющие плотности потока результирующего излучения в направлении осей координат OX OY OZ являются компонентами вектора излучения
Слайд 99Законы теплового излучения.
Закон Планка.
Слайд 100Закон Релея-Джинса
Переходит в соотношение, выражающее закон Релея-Джинса
Слайд 101Закон смещения Вина
Решение этого уравнения приводит к соотношению
Откуда
Слайд 102Закон Планка в безразмерной форме.
Закон смещения Вина позволяет привести закон Планка
к безразмерному виду: