1-ая физическая лаборатория / 58 (приложение) - Учёт теплоёмкости образца
.pdfПРИЛОЖЕНИЕ К ЛАБОРАТОРНОЙ РАБОТЕ •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 всегда больше ¼. Поэтому ®2=®1 всегда пре- вышает 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) |
Время установления этого квазистационарного режима будет порядка 1=®2.
Ðèñ. 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ÑÐ).