Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
математические методы моделирования физических...doc
Скачиваний:
21
Добавлен:
13.11.2019
Размер:
802.82 Кб
Скачать

Если на границе происходит теплообмен, то

λдТ/дх + α(Т – Т0 . (10)

Данное условие носит название граничного условия 3-го рода (закон Ньютона-Рихмана).

1.4 Стационарное распределение температуры в среде

Стационарными называются такие задачи, в которых функция распределения температуры не зависит от времени и является функцией только координаты. Частная производная дТ/дt=0. В этом случае уравнение теплопроводности (6) принимает вид:

d2T(x)/dx2=0, (11)

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

Требуется найти решение уравнения (11), если на границах среды поддерживаются постоянные температуры Т1 и Т2 (рис.1). Граничные условия записываются так:

Т(0) = Т1, Т(1) = Т2. (12)

Общее решения дифференциального уравнения (11) имеет вид:

Т(х) = С1х + С2,

где С1 и С2 – произвольные постоянные. Для их определения необходимо воспользоваться граничными условиями (12). Частное решение удовлетворяющее этим условиям, выглядит так:

Т(х) = (Т2 - Т1)/1+ Т1.

Температура вдоль оси х меняется по линейному закону (рис. 1).

Рассмотрим случай, когда внешняя поверхность стенки толщиной 1 контактирует с массивным твердым телом, вследствие чего на внешней поверхности поддерживается постоянная температура. На внутренней поверхности плоской стенки задана плотность теплового потока Q (рис. 3).

T2 твердое Рис. 3

тело

Q

0 1 х

внутренняя внешняя

поверхность поверхность

Уравнение теплопроводности имеет вид (11), а граничные условия записываются так:

, Т(1) = Т2.

Общее решение, как и в предыдущем случае имеет вид:

Т(х) = С1х + С2.

Производная dT/dx = С1 и из первого условия получаем:

1 = Q, С1 = -Q/.

Второе условие, с учетом найденного значения С1, дает:

С2 = Т2+ (Q/)1.

Подставляя найденные выражения для С1 и С2 в общее решение, получаем:

Т(х) = (Q/)(1 – x) + T2.

Температурный перепад на стенке

Т(0) – Т(1) = (Q/)1.

Рассмотрим теперь стенку толщиной 2d, в которой равномерно распределены источники теплоты, объемной плотностью q. На внутренней и внешней плоскостях имеет место теплообмен с окружающей средой с коэффициентом теплоотдачи , температура среды в дали от поверхности постоянна и равна Т0.

В стационарном случае и уравнение, описывающее распределение температуры в стенке принимает вид:

.

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

,

.

Постоянные интегрирования находятся из граничных условий. В силу постоянства точности источников теплоты в стенке и одинаковости условий охлаждения на внутренней и внешней поверхностях распределение температуры симметрично относительно ее середины. Удобно в этом случае принять за x = 0 середину стенки. Тогда граничные условия имеют вид:

;

.

После подстановки в граничные условия полученного общего решения дифференциального уравнения

.

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

.

В силу симметрии максимальная температура стенки достигается в середине стенки, она равна:

.

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

.

Тепловой поток на охлаждающей поверхности соответственно равен:

,

Откуда температура стенки:

,

А перепад температуры между середины стенки и ее поверхностью равен:

.

Если коэффициент теплоотдачи положить  = , то полученное решение будет представлять собой температурное в стенке при граничных условиях I рода, так как Т(d) = Т0 при  = . В этом случае:

.

При этом температура в середине стенки (x = ):

,

А перепад температуры остается равным:

.