Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

07-Теплопроводность

.pdf
Скачиваний:
41
Добавлен:
18.04.2015
Размер:
458.09 Кб
Скачать

07-Теплопроводность 1

07-Теплопроводность

1.Уравнение Фурье. Условия однозначности.

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

Т а 2Т а х2Т2 у2Т2 z2Т2

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

o Геометрических условий, определяющих форму и размеры тела. o Физических параметров материала и среды.

oНачальных условий, описывающих распределения температуры в некоторый момент времени.

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

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

1.Граничные условия первого рода. Задается распределения температуры по всей поверхности тела и изменения этого распределения во времени Тпов=f(x,y,z,t).

2.Граничные условия второго рода. Задается распределения плотности теплового потока по всей поверхности тела и изменения этого распределения во времени qпов=f(x,y,z,t).

3.Граничные условия третьего рода. Задаются температура окружающей среды и закон теплообмена между средой и поверхностью тела. Этот способ задания граничных условий наиболее общий и часто встречающийся на практике. По существу задается некоторая связь между известной температурой среды и неизвестными температурами на поверхности тела и градиента перепада температуры в теле. Математической формулировкой граничных условий третьего рода в случае КТО на поверхности есть:

к Т f

Тпов

Т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х

 

 

 

 

 

 

 

 

 

 

 

 

Т

 

 

4

 

 

Т

 

 

 

4

 

 

 

 

г

 

 

пов

 

 

В случае лучистого теплообмена:

Cгм

 

 

 

 

 

 

 

 

Т

 

 

 

 

 

 

 

х

 

 

100

 

 

100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обычно в конвективной модели известны к , ,Т f

 

Тсреды , а неизвестны

07-Теплопроводность 2

Т Т Тпов Тцен , где Тпов задаем.

х

ГУ 4 рода.

Т

2

 

2T

 

2T

 

2T

 

 

а

Т а(

 

 

 

 

 

)

 

x2

y2

z2

 

 

 

 

 

ГУ 4 рода соответствуют теплообмену поверхности тела с окружающей средой ( теплообмен в пограничном слое при конвективном теплообмене). Среда отдаёт пограничному слою теплоту, которая конвекцией передаётся поверхности и дальше телу.

ГУ 4 рода.

Тпов ( ) Тсреды ( ) пов

 

 

Т

1

 

 

 

Т2

 

 

 

 

2

 

 

 

п

п

1

 

пов

 

 

пов

Уравнение Фурье совместно с условиями однозначности дают полную математическую формулировку конкретной задачи теплопроводности.

Поставленная задача может быть решена аналитическим, численным или экспериментальным путём.

В результате аналитического решения находится функция температуры от координат x, y,z,

Т=f(x,y,z)

которая должна удовлетворять самому дифференциальному уравнению теплопроводности (оно при подстановке в него этой функции должно обращаться в тождество), а также НУ и ГУ.

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

2.Решение уравнения Фурье при ГУ 3 рода.

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

07-Теплопроводность 3

Неограниченная пластина, толщиной 2δ, нагревается (или охлаждается) за счёт конвекции от окружающей среды с температурой Т f const

Вначальный момент времени (τ=0) темпе6ратура пластины по толщине одинакова Тн const .

Впроцессе нагрева температура в любой точке пластины будет увеличиваться, стремясь к температуреТ f .

Тн Т f

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

Математическая формулировка задачи.

Для определения функции Т(х,τ) имеем уравнение теплопроводности.

 

Т

2Т

 

 

 

 

 

 

а

 

 

 

 

 

(1)

 

 

х

2

 

 

 

 

 

 

 

 

 

 

НУ:Т(х, ) Тн

 

 

(2)

 

 

 

 

 

 

Т

 

ГУ: Т f

Т( , )

 

 

(3)

 

 

 

 

 

 

 

х х

 

Знак (+ ) перед левой частью ГУ условию на левой поверхности пластины, где х=-д, т. е. когда направление теплового потока совпадает с положительным направлением оси х. А знак (-) – условию на правой поверхности пластины при х=+д, когда тепловой поток направлен в противоположную сторону оси х.

Для решения этой задачи воспользуемся методом разделения переменных.

Сущность его заключается в том, что решение уравнения (1) в виде произведения 2 функций.

Т(х, ) Ф(х)*П( )

(4)

Ф зависит от координаты х. П – от координаты ф.

 

 

 

 

 

 

 

 

 

 

 

 

07-Теплопроводность

4

Подставим функции в уравнение (1) и получим:

 

 

Т

Ф

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

2Т

 

П

d2Ф

 

 

 

х2

2

 

 

 

 

 

 

Ф

П

аП

d2Ф

 

(5)

 

2

 

 

 

 

 

 

 

 

 

Левая часть выражения (5) является функцией времени ф, а правая только функцией координаты х.

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

Обозначим постоянную величину через (2).

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

1 * ак2

Пd

d2Ф к2 *ф 0 2

(6)

(7)

Решение уравнения (6) имеет вид:

П( ) С1 *ехр( ак2 )

 

 

 

 

(8)

А решение уравнения (7) имеет вид:

(9)

Ф(х) С2 *cos(kx) C3 *sin(kx)

 

 

В (8) и (9) С1, С2, С3 – произвольные постоянные.

 

Подставим (8) и (9) в (4) и получим, что

 

Т(х, ) С *ехр( ак2 )* С

2

cos(kx) C

3

sin(kx)

 

1

 

 

 

Обозначим С12

 

 

 

 

 

С13=D

 

 

 

 

и получим

 

 

 

 

 

Т(х, ) ехр( ак2 )* Сcos(kx) Dsin(kx)

(10)

В виде (10) и будем искать решение уравнения (11), для чего нужно определить значение констант С, D, к из НУ и ГУ.

D определяется просто:

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

f(x)=f(-x)

РИСУНОК!!!!!!!!!!!!!!

f(-x)=-f(-x)

sin(kx) является нечётной функцией, поэтому D=0 для любой

07-Теплопроводность 5

симметричной задачи. И тогда решение уравнения (1) нужно искать в виде:

Т(х, ) С*ехр( ак2 )*сos(kx)

Для решения задачи введём избыточную температуру (x, )

(х, ) Т f Т(х, )

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

 

 

а

2

 

 

 

 

 

 

х2

 

 

 

 

 

 

 

 

НУ: (х,0) Т f

Тн

н

 

 

 

 

 

 

 

ГУ: ( , )

 

 

 

 

 

 

 

 

 

х х

Знаки в левой части ГУ изменились на обратные, т. к.

Тх х

(12)

(12а)

(12б)

Так как уравнение (12) для избыточной температуры имеет тот же вид, что и (1) для абсолютной температуры, то решение (12) будем искать в виде

(11), т. е.

(х, ) С*ехр( ак2 )*сos(kx)

(13)

Константы С и к определим из краевых условий задачи (12а) и (12б).

В начале используем ГУ (12б), например, на левой поверхности пластины, где х=-д.

Подставим решение (13) в уравнение (12б) и получим

*С*ехр( ак2 )*cos(k ) *C*k *exp( ak2 )*sin(k )

сos(k ) ctg(k ) k sin(k )

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

ctg(k ) k

Из этого уравнения найдём к. Для удобства будем искать не саму к, а величину кд.

кд=м

м – характеристическое число, а комплекс

Bi - число Био (критерий Био).

Для отыскания одной из констант к, уравнение (11), а именно число м, в которое оно входит, мы получим характеристическое уравнение.

ctg

 

(14)

 

Bi

Так как решение (16) должно удовлетворять не только граничным условиям (12б), но и начальным условиям (12а), то Сn определим из начальных условий (12а). Для этого положим в решение (16), что критерий

07-Теплопроводность

6

Фурье равен нулю, а это значит время равно нулю.

n Сn cos( n X) (17)

n 1

Данное выражение представляет собой разложение величины n в ряд с заданными значениями n , которое определяется уравнением (14).

Используя уравнение (17), можно показать, что сos( n X) и сos( m X) ( m - любое характерное число, выбранное из бесконечного множества значенийт )являются ортогональными функциями, что означает если

cos( n X) cos( m X)dx 0

при любых значениях n не равных значениям m, то этот интеграл равен нулю. И только при n=m является величиной конечной.

Используем это обстоятельство для определения значения Сn. Для чего умножим обе части выражения (17) на сos( n X) и проинтегрируем по х в пределах толщины пластины, и поменяем порядок суммирование и интегрирование и получим выражения

n cos( m X)dx Сn соs( m X) cos( n X)dx

n 1

В правой части последнего выражения все интегралы равны нулю кроме одного, соответствующему случаю n=m. Поэтому это выражение получает вид:

 

 

 

 

 

 

 

n cos( n X )dx Сn соs 2 ( n X )dx

 

n 1

 

 

 

 

 

Отсюда определим Сn:

 

 

 

 

 

 

 

 

2

 

 

 

cos( n X)dx 2 cos( n X)dx

sin n

 

 

 

 

0

 

 

n

 

 

 

 

 

 

 

 

 

 

n

sin n cos n

cos2 ( n X)dx 2 cos2 ( n X)dx

 

n

 

0

 

 

 

 

 

Cn 2 n

 

sin n

(18)

 

 

n

sin n cos n

 

 

Решение (16) принимает вид:

 

 

2sin n

 

 

 

 

 

 

(19)

n

 

 

cos( n X)e n2 Fo

 

sin n cos n

n 1 n

 

 

 

 

 

 

Следовательно, задача решена, так как из (19) можно найти значение избыточной теплоты, а значит и абсолютные температуры в любой точки пластины в любой момент времени. В решении (19) все независимые переменные фигурируют в безразмерном виде. Роль координаты играет безразмерная величина X x , а роль времени – число Фурье. Поэтому

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

07-Теплопроводность 7

безразмерном виде. Для этого нужно разделить обе части уравнения (19) на начальную избыточную температуру n , которая известна из начальных условий, тогда получим решение задачи в безразмерном виде:

 

 

 

Tf

T

 

 

 

sin n

 

cos( n X) e

n2Fo

(20)

 

n

T

f

T

 

n

sin

n

cos

 

 

 

 

 

 

n

 

 

 

 

n

 

 

- относительная безразмерная избыточная температура. Общий вид этого решения расписан так:

f (Bi;Fo;X) (21)

В (21) Fo и Х безразмерные независимые переменные, а Bi – параметр задачи. Термин «безразмерное время», относящееся к ..., для пластин различной толщины и с различными коэффициентами теплопроводности, но при одинаковых значениях критерия Био, температурные кривые построены в координатах f (X) будут совпадать при равенстве чисел Фурье, то есть не

......... те же моменты времени, а сходственные, которые определяются из равенства чисел Фурье. Конкретный вид решения (20) получим для неограниченной плоской пластины, но принципиально в таком же виде, то есть в виде уравнения (21) получаем решения и для тел любой другой геометрической формы. Изменяется только смысл вкладываемый в понятие «характерный размер тела». В случае цилиндра и шара роль характерного размера играет радиус.

Рассмотрим физический смысл критерия Био. Для чего представим его в

виде Bi

 

 

, где

интенсивность внешнего теплообмена, а

 

 

 

 

 

 

 

 

 

 

 

характеризует интенсивность внутреннего теплообмена между поверхностями тела и его внутренними слоями. Критерий Био можно

представить в виде Bi

 

 

, где - внутреннее тепловое сопротивление, а

 

1

 

 

 

 

 

1

его внешнее тепловое сопротивление. Поэтому можно сделать вывод о

 

 

 

 

 

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

Для выявления влияния Био не процесс нагрева (охлаждения) тела, сделаем анализ решения (20) при очень малых значениях, Био стремится к нулю, и при очень больших значениях, Био стремится к бесконечности.

При Био, стремящемуся к нулю.

Такие тела называют «тонкими» в тепловом отношении телами. Угловой

коэффициент 1Bi

07-Теплопроводность 8

в пороговой части выражения (14) стремится к

бесконечности. И, следовательно, прямая на графике совпадает с осью координат. А корни характеристического уравнения пишут значения 0,

,2 ....(п 1) .

Bi при Вi 0

Тогда из (18) видно, что все коэффициенты Сn, кроме первого, обращается в ноль. А для первого получается неопределенность типа 0/0, которую можно раскрыть.

Обозначим коэффициент, стоящий в (20) уравнении под знаком суммы через:

Dn

 

2sin n

n sin n cos n

 

 

Dn пропорциональна коэффициенту μn, как это следует из выражения (18). Раскроем с помощью предельного перехода неопределенность для D1 при

1 0

lim D1

lim

2sin 1

lim

2 1

 

 

 

0

0 1 sin 1 cos 1

0 1 1 1

При малых значениях 1 имеем: sin 1

1 , tg 1 1

Следовательно, решение (20) для этого случая примет вид:

cos( 1X)e n2Fo

(22)

 

 

 

 

 

 

 

В решении (22) нужно определить конкретный вид связи между 1 и Bi.

Уравнение (14) принимает вид:

1

 

1

 

Следовательно, 1

 

, а решение

 

Bi

 

Bi

(22) примет вид

 

1

 

 

 

 

 

 

 

 

 

 

 

cos(

 

 

 

 

 

 

 

 

 

 

 

BiX)e BiFo

 

 

 

 

 

 

 

 

lim cos(

 

 

 

 

 

 

 

 

 

BiX) 1, так как 0 X 1

 

 

 

 

 

Bi 0

 

 

 

 

 

 

 

 

 

 

 

И окончательно получаем e BiFo

(23)

 

 

Следовательно,

для случая Bi 0 получим равномерное распределение

Bi
Bi

07-Теплопроводность 9

температуры по объему тела на протяжении всего периода нагрева или охлаждения тела, так как в выражении (23) не зависит от Х. Получаем, что температура во всех точках тела растет или падает во времени по экспоненциальному закону.

«Тонкими» в тепловом отношении телами, которые нагреваются или охлаждаются по всему объему, называются тела, у которых Bi 0.25. При нагреве «тонких» тел лимитирующим звеном в процессе нагрева является внешний теплообмен.

Лекция 14

3. Решение уравнения Фурье для твердого тела при граничных условиях первого рода

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

По формуле закона Ньютона Tf

Tпов

 

q

 

разность температур стремится к

 

 

 

 

 

бесконечности. При Bi (практически

при Bi 100) поверхность тела

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

 

07-Теплопроводность

10

Математическая формулировка задачи.

д

а

д2

 

 

 

 

 

 

д дх2

 

Начальные условия: (х,0) н

 

Граничные условия: ( , ) 0

 

Т f Т

н Т f Тн

 

Решение этой задачи получим как частный случай решения (20) с задачей граничного условия третьего рода при Bi .

Характеристическое уравнения (14), в этом случае, принимает вид ctg 0, следовательно, прямая линия на графике совпадает с осью абсцисс и значение

корня уравнения (14)

получают

 

;

3

;...(2п 1)

 

(24). Поэтому sin n 1 при

 

 

 

 

 

 

 

2

2

2

 

 

 

 

нечетных n, а sin n 1

при четных n, то есть sin n ( 1)n 1 . Так как cos( n ) 0, то

 

2( 1)

n 1

 

 

 

 

 

 

 

2( 1)

n 1

Cn

 

и решение этой задачи принимает вид:

 

cos( n X)e n2Fo (25)

n

 

n

 

 

 

 

 

 

 

 

 

 

n 1

 

 

Так как в решении (25) n определяется последовательностью (24), и, следовательно, n f (Bi), то есть при граничных условиях первого рода общий вид решения (25) имеет вид: f (Fo, X)

Вэтом случае Bi не является параметром задачи и лимитирующим звеном

впроцессе нагрева или охлаждения, является внутренним теплообменом.

4. Решение уравнения Фурье при граничных условиях второго рода. Математическая формулировка задачи.

 

дТ

а

д2Т

 

 

 

 

 

д

дх2

 

 

 

 

 

 

 

 

Начальные условия: Т(х,0) Тн

 

 

 

 

 

дТ

 

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

 

 

qпов

 

 

 

 

 

 

дх х

 

Особенность этой задачи: при граничных условиях первого и третьего родов температурное поле с течением времени стремится к равновесию, то есть температура в любой точке тела стремится к температуре среды при граничных условиях третьего рода или к температуре поверхности при граничных условиях первого рода. Для решения этой задачи введем вместо температуры новую переменную плотность теплового потока q(х, ), связанную с температурой формулой закона Фурье:

q(х, ) дТ(х, )

дх

Продифференцируем по х обе части уравнения дТ а д2Т .

д дх2

д

дТ

 

д

 

2

Т

 

а

 

д

 

 

 

 

 

 

 

 

2

 

 

 

 

дх

 

дх

д

 

дх

 

 

Соседние файлы в предмете Физика