Материал: Д6768 Аналитическое описание процессаинестационарной теплопроводности

Тогда вместо формул (36)–(38), (39′) будем иметь

; (42)

. (44)

Частное решение дифференциального уравнения (41) можно написать, воспользовавшись результатом разобранного выше метода разделения переменных – решением уравнения (25):

Из условия симметрии следует А = 0, а именно:

вий охлаждения: распределение температуры должно быть симметричным относительно оси ординат (оси температур), следовательно, должно описываться четной функцией, т. е. в соотношении (45) должна быть только четная функция соs(kx).

Таким образом,

Удовлетворим соотношение (46) граничному условию (43):

26

откуда следует:

т. е. значений k существует бесчисленное множество. Следовательно, общее решение будет суммой всех частных ре-

шений:

n 1

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

будет представлена в виде ряда Фурье:

n 1

в котором коэффициенты определяются по соотношению

Решение исходной краевой задачи (36)–(39) получим, заменив по формуле (40) t(x, τ) на выражение

Тогда окончательное решение будет иметь следующий вид:

27

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

Замечание. Решение (52) является решением исходной краевой задачи (36)–(39) и одновременно решением задачи нахождения температурного поля в неограниченной пластине толщиной l R (0 < x < l), когда одна ее поверхность имеет тепловую изоляцию (при

противоположная поверхность x = l поддерживается при температуре tc. В начальный момент времени задано некоторое распределение температуры в виде функции f(x), причем в этом случае функция f(x) может быть любой, лишь бы она удовлетворяла условиям Дирихле.

Если начальное распределение температуры равномерное, т. е. f(x) = t(x, 0) = t0 = const, то после вычисления интеграла

Bn 2

2

решение (52) можно написать в безразмерном виде:

28

где θ = θ(X, Fо) – безразмерная (относительная избыточная) темпера-

число (критерий) Фурье (критерий гомохронности – безразмерное время – для процессов теплопроводности).

Некоторые выводы из анализа решения

1. Из анализа решения уравнения (53) следует, что относительная избыточная температура является функцией относительной координаты и числа Фурье: θ = θ(X, Fо).

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

Такой процесс теплообмена называется внутренним процес-

сом, а сама задача – внутренней задачей.

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

рис. 5.

Рис. 5. Распределение температуры в неограниченной пластине (внутренняя симметричная задача)

29

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

сходимости зависит, прежде всего, от величины Fo. Для достаточно больших Fo ряд сходится быстро, и для вычисления θ достаточно ограничиться одним–двумя членами ряда.

Пример

Пусть имеется пластина толщиной 2R = 0,2 м, t0 = 35 °C , tс = 5 °C. При этом коэффициент температуропроводности а = 1,4∙10–7 м2/с. Необходимо найти температуру в центре пластины (х = 0) через 10 часов после начала охлаждения.

Значение числа Fo достаточно велико, поэтому можно ограничиться только первым членом ряда:

t = 16,1 ° C.

2. На практике приходится решать обратную задачу: находить продолжительность охлаждения τ, т. е. время, через которое в определенной точке пластины будет достигнута необходимая температура. Данную задачу легко решить в том случае, если в уравнении (53) можно ограничиться первым членом, а остальные отбросить из-за их малой величины. Вычисления показывают, что это возможно при Fo > 0,4.

30