Нестационарная теплопроводность

Содержание

Слайд 2

Условия однозначности процесса задаются в виде: – физических параметров – коэффициента

Условия однозначности процесса задаются в виде:
– физических параметров – коэффициента

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

– граничных условий в виде:

– значение соответствующей координаты.

Иногда начальные и граничные условия объединяются под названием «краевые условия».

Слайд 3

Охлаждение неограниченной пластины Рассмотрим неограниченную пластину (толщина пластины значительно меньше ее

Охлаждение неограниченной пластины
Рассмотрим неограниченную пластину (толщина пластины значительно меньше ее длины

и ширины) толщиной

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

Таким образом, задача может рассматриваться как одномерная.

К охлаждению плоской неограниченной пластины

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

(граничные условия второго рода). Начальное распределение температуры в пластине запишем в виде:

Введем обозначение

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

Слайд 4

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

При начальных условиях

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

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

при

при

Для решения уравнения теплопроводности используем метод разделения переменных Фурье, т.е. ищем решение в виде произведения двух функций, одна из которых является функцией только времени, а вторая – только координаты:

Тогда:

Или

Или

Слайд 5

Постоянную выбираем из граничных условий, а знак «-» выбираем из физических

Постоянную

выбираем из граничных условий, а знак «-» выбираем из

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

Или

Или

При

Принимая

Или

Или

И окончательно:

Тогда:

Полученное решение принадлежит дифференциальному уравнению теплопроводности при любых

Слайд 6

Учитывая граничные условия, имеем: При Или Это возможно только при ,

Учитывая граничные условия, имеем:
При

Или

Это возможно только при

, т.е. при


Таким образом, частное решение

должно быть отброшено, как не удовлетворяющее заданным граничным условиям. Если обозначить

, тогда решение будет иметь вид:

Далее рассмотрим граничное условие при

Тогда

Умножая числитель и знаменатель на

, получаем:

Принимая во внимание, что

(число Био), и обозначая

Слайд 7

Получаем трансцендентное уравнение: При каждом значении существует бесконечное множество решений. Для

Получаем трансцендентное уравнение:

При каждом значении

существует бесконечное множество решений.

Для поиска решения используется графоаналитический метод (рис.). Обозначим

Пересечение графиков этих функций дает значение корней характеристического уравнения

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

соответствует свой набор значений

Если

возрастает, то прямая

стремится к оси абсцисс; если

то прямая

совпадает с осью абсцисс. При

прямая

совпадает с осью ординат. Таким образом, любому значению

Слайд 8

будет соответствовать свое распределение температуры: . . . . . .

будет соответствовать свое распределение температуры:

. . . . . . .

. . . . . . . . . . .

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

Слайд 9

Откуда определим начальное распределение температуры, учитывая начальные условия Уравнение есть не

Откуда определим начальное распределение температуры, учитывая начальные условия

Уравнение есть не

что иное, как разложение четной функции в ряд Фурье, откуда:

Умножим (1) на

(1)

и проинтегрируем по

от

до

или

Слайд 10

Тогда выражение для примет вид: После преобразований: Уравнение определяет температурное поле

Тогда выражение для

примет вид:

После преобразований:

Уравнение определяет температурное поле

для любого момента времени при любом начальном распределении температуры. Если при

температура

, то:

Слайд 11

Отсюда Распределение температуры может быть записано: Переходя к безразмерным координатам (критерий

Отсюда

Распределение температуры может быть записано:

Переходя к безразмерным координатам

(критерий или

число Фурье – безразмерное время), получим:

Анализируя, видим, что

представляет собой ряд возрастающих чисел. Чем больше

Слайд 12

тем быстрее убывают члены ряда. При ряд становится быстросходящимся и с

тем быстрее убывают члены ряда. При

ряд становится быстросходящимся и с

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

Или, введя обозначение

Учитывая, что

то величина

является функцией только

и может быть определена для различных значений

.

Рассмотрим характер изменения температуры при заданных значениях

(рис.). Если

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

(рис.а), то на поверхности пластины

Слайд 13

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

Зависимость распределения температуры от

При малых

(рис.б) температура пластины незначительно отличается

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

Количество теплоты, передаваемой пластиной за время от

до

равняется изменению ее внутренней энергии за период полного охлаждения до

-масса, удельная теплоемкость, плотность, площадь поперечного сечения пластины.