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

1-ая физическая лаборатория / 58 (приложение) - Учёт теплоёмкости образца

.pdf
Скачиваний:
6
Добавлен:
25.04.2014
Размер:
169.11 Кб
Скачать

ПРИЛОЖЕНИЕ К ЛАБОРАТОРНОЙ РАБОТЕ •58

Учет теплоемкости образца.

1. Точное решение уравнения теплопроводности.

Формулы (4)-(6) основного текста справедливы, когда теплоемкостью образца C0 можно пре- небречь по сравнению степлоемкостью приемника. Если считать C0 = 0 (тепло не поглощается образцом), то количество тепла, втекающее за время dt в каждый элементарный объем dV , равно

количеству тепла, вытекающего из него. Поэтому дивергенция вектора плотности теплового потока ~q должна в пределах образца равняться нулю.

Åñëè æå C0 =6 0, то необходимо написать уравнение теплового баланса для элементарного объема. Рассмотрим элементарный объем dV = dxdydz с центром в точке x; y; z (рис.1). Через левую

грань в него втекает за врем dt количество тепла, равное qx(x ¡ dx2 ; y; z) dydzdt, через правую - вытекает - qx(x + dx2 ; y; z) dydzdt. Разность этих количеств равна ¡@q@xx dxdydzdt. Рассматривая ана-

логичным образом другие пары граней, находим, что объем dV за время dt получает количество

 

³

@qy

 

´

 

 

 

 

 

 

 

 

 

тепла ¡

@qx

@qz

dV dt = div~q dV dt. Это тепло идет на нагревание вещества в данном объ-

@x

+ @y

+ @z

еме, температура

 

которого в результате повышается на

 

 

div~q dV dt

 

 

T

dT = ¡ C0dV , ãäå C0

- теплоемкость

 

 

 

 

 

 

 

 

 

единицы объема вещества. Учитывая, что qx = æ @T@x ; ¢ ¢ ¢ , получим

 

 

 

 

 

 

C0 @x

= æ µ @x2

+ @y2 +

@z2

:

 

 

(1)

 

 

 

 

 

 

@T

@2T

 

@2T

@2T

 

 

 

 

Это уравнение обычно называют уравнением теплопроводности (хотя оно включает также и уравнение теплового баланса).

В нашем случае T зависит только от одной координаты z, и уравнение (1) упрощается:

C0

@T

= æ

@2T

:

(2)

@t

@z2

 

 

 

 

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

(z = 0) : T (0; t) = Tí = const

 

(z = d) : T (d; t) = Tõ;

Cõ @Tõ

= ¡æ

@T (z; t)

jz=d

 

 

 

 

 

S @t

 

@z

(3)

(4)

Здесь Cх - теплоемкость всего приемника, площадь верхней (примыкающей к образцу) поверхности которого равна S, Tн - температура нагревателя и Tх - температура приемника. Кроме того нужно задать начальное условие распределение температуры в начальный момент t = 0. Ôèçè-

чески ясно, что задание начальных и граничных условий однозначно определяет функцию T (z; t)1;

соответствующие теоремы о существовании и единственности доказываются в курсе матеметики. Нетрудно убедиться, что уравнение (2) имеет следующие частные решени

T1(z; t) = C ; T2(z; t) = D z ;

 

T3(z; t) = A e¡®t Sin¯z ; T4(z; t) = B e¡®t Cos¯z ;

(5)

ãäå A, B, C è D произвольные постоянные, а константы ® è ¯ связаны уравнением

 

C0® = æ ¯2 :

(6)

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

В силу линейности уравнения (2) любая линейная комбинация функций (5) также является реше-

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

(3) и (4) и начальному условию

t = 0 : T (z; 0) = T0(z) ;

(7)

 

ãäå T0(z) - произвольная функция (начальное распределение температуры может быть любым).

Ðèñ. 1:

Условие (3) будет выполняться при любых t, åñëè C = Tí, B = 0. Ïðè t ! 1 решение должно иметь вид T (z; t) = Tõ = Tí, поэтому D = 0. Условие (4) дает

Cõ

® A e¡®t Sin¯d = +æ ¯ A e¡®t Cos¯d :

(8)

S

 

 

Это второе уравнение, связывающее ® è ¯. Сокращая на A e¡®t и подставляя ® из (6), имеем

tg¯d =

S C0

:

(9)

Cõ¯

 

 

 

Это трансцендентное уравнение удобно решать графически, построив функции y1 = tg¯d è y2 =

d S C0

1

в зависимости от x = ¯d. Решение найдетс как точка пересечения графиков y1(x) è y2(x). Èç

Cõ

¯d

рис.2 видно, что имеетс бесконечное множество решений1. При малых S C0

d наименьшее решение

 

 

Cõ

 

¯1d будет лежать в области малых x, и его можно найти, используя разложение тангенса в ряд

Мак-Лорена

 

 

 

(¯1d)3

 

 

 

 

S C0

 

 

 

1

 

 

 

 

 

 

 

¯1 d +

+ ¢ ¢ ¢ =

d

 

:

 

 

(10)

 

 

3

 

Cõ

 

¯1d

 

 

Ограничиваясь первым членом, находим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

 

 

 

 

 

 

 

 

 

¯1 d ' r

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S 0

d ; ®1

'

 

 

 

S :

 

 

 

(11)

 

 

Cõ

d Cõ

 

 

 

Остальные значения ¯d будут очень близки к ¼, 2¼, ¢ ¢ ¢ , òî åñòü

 

 

 

 

¯

n '

(n ¡ 1) ¼

;

®

 

 

 

(n ¡ 1)2¼2

 

æ

 

(n = 2; 3;

¢ ¢ ¢

) :

(12)

 

 

 

 

C0

 

d

n '

 

 

 

d2

 

 

 

 

 

 

 

 

Таким образом, мы получили решение нашей задачи в форме

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T (z; t) = Tí +

 

 

 

Ame¡®mt Sin¯mz :

 

 

(13)

m=1

1Подумайте, почему можно не рассматривать область ¯ < 0.

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

X1

T0(z) ¡ Tí = Am Sin¯mz :

(14)

m=1

 

Модно показать, что любая физически осуществимая функция T0(z) ¡ Tí, обращающаяся в нуль при z = 0, может быть прадставлена в виде такого ряда. Иначе говоря, функция (13) удовлетворяет

произвольным начальным условиям, то есть является общим решением задачи о распространении тепла в системе, описываемой уравнениями (2)-(4).

2. Переходный процесс.

Из рис.2 видно, что ¯1d всегда меньше ¼=2, à ¯2d всегда больше ¼. Поэтому ®21 всегда пре- вышает 4, а при выполнении приближения (11)

®2

¼2

 

Cõ

 

À 10 :

®1

S 0 d

Поэтому все члены с m > 1 в (13) убывают намного быстрее первого члена, и через достаточно длительное врем после начала опыта мы получим

T (z; t) = Tí + A1 e¡®1t Sin¯1z ; T0 = Tí + A1 e¡®1t Sin¯1t :

(15)

Время установления этого квазистационарного режима будет порядка 12.

Ðèñ. 2:

3. Измерение теплопроводности.

Приближение (11) соответствует замене синусоидального распределения температур в образце

линейным. Улучшенное приближение можно получить, переписывая (10) в форме

¯1 d

µ1 +

3

= Cõ

d ;

2

2

¯12d2

S C0

 

и подставляя (11) в поправочный член в скобках. Тогда

d S

+

 

 

 

 

' r

 

Cõ

µ ¡

3Cõ

3

 

 

 

 

C0 S d

1

C02 d2 S2

; ®1 =

 

 

æ

 

 

=

Отсюда æ = ®1 d Cõ ÝÔÔ

 

 

 

2

¡

Cõ

 

C0 d

¢

¯10d

 

 

 

 

 

 

 

S , где эффективна теплоемкость приемника

1

Cõ ÝÔÔ = Cõ + 3C0 d S :

æ S

:

(16)

d Cõ ÝÔÔ

Это приближение будет хорошим, если выполнено условие

C0 d S

¿ 1 ;

(17)

Cõ

(точнее, если квадрат этой величины мал по сравнению с относительной погрешностью измерения æ). В общем случае мы можем получить из (6) и (9)

 

 

 

æ =

®1 C0

=

®1 Cõ d tg¯1 d

:

 

 

(18)

 

 

 

 

 

 

 

 

 

 

 

 

 

¯12

 

S

 

¯1 d

 

 

 

 

таким образом, точное значение эффективной теплоемкости равно Cõ =

tg¯1 d

. Измеряя

®1 =

1

dTõ

æ, åñëè Cõ ÝÔÔ известно.

 

 

 

 

 

¯1 d

 

 

Tí¡Tõ

dt мы найдем

 

 

 

 

 

 

 

 

4. Измерение эффективной теплоемкости приемника.

При измерении теплоемкости приемника включается печка нагрева приемника. С учетом мощности W, выделяемой этой печкой, граничное условие (4) заменяется на

(z = d) : T (d; t) = Tõ ;

Cõ @Tõ

=

¡

æ

@T

+

W

:

(19)

 

 

 

 

S @t

@z

S

 

 

 

 

 

 

Остальные условия остаются неизменными. Таким образом, в (5) по-прежнему C = Tí, B = 0. Ïðè

t ! 1 должно установитьс стационарное распределение температур, @T@t = 0, и тогда (19) дает D = æWS . Поэтому решение следует искать в виде

 

 

 

 

1

 

 

 

 

W

X

 

T (z; t) = Tí +

æ S

z + Am e¡®mt Sin¯mz :

(20)

 

 

 

 

 

m=1

 

Подстановка в (19) дает

 

 

 

 

 

 

1

 

Cõ

®m Sin¯md + æ ¯m Cos¯md= 0 :

(21)

m=1 Am e¡®mt

µ¡ S

X

 

 

 

 

 

 

Это уравнение выполняется при всех t, если выражение в скобках равно нулю при любом m. Таким образом мы получаем прежний набор ®m è ¯m, а функция (20) при подходящем выборе Am

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

 

 

 

 

 

 

 

 

Считая, что с момента начала эксперимента прошло достаточное время, отбросим все члены

ñ m > 1. Далее, пусть в некоторый момент t0 температуры с обеих сторон образца сравнялись,

Tí = Tõ. Тогда

 

 

 

 

 

 

 

W

 

 

 

 

 

 

 

 

 

 

A1 e¡®t0 Sin¯1d = ¡

d :

 

 

(22)

 

 

 

 

 

æ

 

 

В этот же момент, согласно (19),

W ¡

 

µ æ

 

 

 

 

 

 

 

 

 

S dt

 

 

1

 

 

 

1

1

 

 

Cõ

 

dTõ

=

 

æ

 

W

+ A

 

e¡®1t0 ¯

 

Cos¯ d

:

(23)

 

 

 

 

 

 

 

Подставляя A1 из (22), находим

1

Cõ

tg¯1d dTõ

= W = Cõ ÝÔÔ

1

 

dTõ

 

(24)

 

 

¢

 

 

 

 

 

;

S

¯1 d

dt

S

 

dt

ãäå Cõ ÝÔÔ имеет в точности то же значение, что и в предыдущем эксперименте. Таким образом,

мы можем не заботиться о поправках на теплоемкость образца. Достаточно определить Cõ ÝÔÔ =

Cõ tg¯1d

¯1 d из (24) и подставить в (18).

5. Задание для самостоятельной работы.

При решении нашей задачи мы считали, что приемник изготовлен из материала с бесконечной теплопроводностью. Попробуйте рассмотреть самостоятельно случай, когда это приближение недопустимо. Для этого необходимо решать два уравнени теплопроводности в материале образца и материале приемника. Условие (4) заменится условием равенства температур и тепловых потоков с обеих сторон от плоскости z = d. На второй границе приемника можно принять условие теплоизо-

ляции (q = 0) или условие теплоотдачи в среду с постоянной температурой TÑÐ (q » T ¡ TÑÐ).

Соседние файлы в папке 1-ая физическая лаборатория