Теплопроводность пластин
Все материалы способны проводить тепло, но по-разному! В данной задаче вам предстоит рассчитать теплофизические характеристики некоторых материалов. Закон теплопроводности полностью совпадает с законом Ома для участка цепи. Он был сформулирован французским физиком Жаном Батистом Жозефом Фурье в 1807 году.

Закон утверждает, что количество теплоты, протекающей через однородную пластину толщины \(h\) и площади \(S\) в единицу времени (т.е. мощность теплопередачи), определяется по формуле
\[q = \lambda S \frac{\Delta T}{h}\tag{1}\]
Коэффициент \(\lambda\) в этой формуле называется коэффициентом теплопроводности и является характеристикой вещества.
Замените в этой формуле температуру на электрический потенциал, поток теплоты на силу тока, коэффициент теплопроводности на удельную электрическую проводимость и... перед вами закон Ома.
Формула (1) справедлива, если коэффициент теплопроводности является постоянной величиной, если нет – разбивайте пластину на тонкие слои... и для каждого слоя пишите отдельные уравнения. Впрочем, удельное сопротивление также может изменяться, например, если проводник разогревается.
Во всех частях данной задачи мы будем рассматривать пластины одинаковых размеров \(S = 1{,}00\,\text{м}^2\), \(h = 10{,}0\,\text{см}\). Для воздушной пластины будем подразумевать, что воздух находится в коробке указанных размеров, масса воздуха внутри нее остается неизменной. Влиянием самой коробки на протекание рассматриваемых процессов будем пренебрегать. Также будем рассматривать одинаковые начальные условия: пластина имеет постоянную температуру \(T_0 = 0{,}0^\circ\text{C}\). Одна сторона пластины все время поддерживается при этой температуре, затем вторую сторону пластины приводят в тепловой контакт с очень массивным телом (термостатом), находящимся при температуре \(T_1 = 20{,}0^\circ\text{C}\). Координату холодной стороны полагаем равной нулю, координата второй - \(h\). По прошествии некоторого промежутка времени тепловой поток через пластину становится стационарным, то есть не зависящим от времени. Процесс установления равновесия достаточно сложный,
нас он интересоват не будет, основное внимание уделяется стационарному распределению температур и тепловым потокам в стационарном режиме. Часть 1. Металлические пластины.
В этой части рассматриваются две пластины указанных размеров, серебряная и свинцовая. Для металлов коэффициенты теплопроводности можно считать постоянными. Табличные характеристики выбранных металлов приведены в таблице.
Металл Плотность (\(\rho, \frac{\text{кг}}{\text{м}^3}\)) Удельная теплоемкость (\(C, \frac{\text{Дж}}{\text{кг} \cdot \text{К}}\)) Коэффициент теплопроводности (\(\lambda, \frac{\text{Вт}}{\text{м} \cdot \text{К}}\)) Серебро, \(Ag\) \(10{,}5 \cdot 10^3\) \(235\) \(429\) Свинец, \(Pb\) \(11{,}3 \cdot 10^3\) \(127\) \(35{,}1\) 1.1 Рассчитайте потоки теплоты через каждую пластину в установившемся режиме \(q_{Ag}\), \(q_{Pb}\).
1.2 Рассчитайте количество теплоты, которое пошло на нагревание каждой пластины \(Q_{Ag}\), \(Q_{Pb}\).

Теперь соединим обе пластины. Сначала с холодной стороны находится свинцовая пластина (а серебряная с теплой), потом наоборот.
Для каждой составной пластины:
1.3 Рассчитайте потоки теплоты в установившемся режиме \(q_1\), \(q_2\).
1.4 Выразите значения \(q_1\), \(q_2\) через значения потоков \(q_{Ag}\), \(q_{Pb}\), найденных в п.1.1.
1.5 Рассчитайте количества теплоты \(Q_1, Q_2\), которые пошли на нагревание пластин до установившегося режима (с одной стороны температура \(T_0\), с другой \(T_1\)).
1.6 Постройте графики зависимости температуры внутри пластин от координаты (на одном бланке.
Часть 2. Воздушная пластина.
Теплоемкость воздуха зависит от температуры, кроме того, внутри коробки может происходить и перераспределение плотностей воздуха, поэтому, на первый взгляд, задача становится не решаемой! Но попробуйте! Конвективные явления не учитывать!
Коэффициент теплопроводности воздуха при температуре \(T_0 = 0{,}0^\circ\text{C}\) равен \(\lambda_0 = 0{,}024 \frac{\text{Вт}}{\text{м} \cdot \text{К}}\).
Молярная теплоемкость воздуха при постоянном объеме равна \(C_V = \frac{5}{2} R\). Давление воздуха в коробке при начальной температуре \(T_0 = 0{,}0^\circ\text{C}\) равно \(P_0 = 1{,}0 \cdot 10^5 \text{ Па}\).
Математическая подсказка: при малых изменениях \(\Delta x\) справедливо равенство
\[x \Delta x = \frac{1}{2} \Delta(x^2)\]
Согласно упрощенной модели теплопроводности газов, их коэффициент теплопроводности должен определяться по формуле
\[\lambda = \frac{1}{3} \langle v \rangle \langle l \rangle \rho c_V\tag{2}\]
где \(\langle v \rangle\) — средняя скорость теплового движения молекул, \(\langle l \rangle\) — средняя длина свободного пробега молекулы, \(\rho\) — плотность газа, \(c_V\) — удельная теплоемкость газа в изохорном процессе.
2.1 Длина свободного пробега молекул зависит от эффективного диаметра молекулы \(d\) и концентрации молекул \(n\). Установите вид зависимости длины свободного пробега от указанных параметров.
2.2 Коэффициент теплопроводности пропорционален абсолютной температуре в некоторой степени \(\lambda = \beta T^z\). Определите показатель степени \(z\) в этой формуле.
Независимо от полученного вами результата далее считайте, что теплопроводность воздуха прямо пропорциональна абсолютной температуре, т.е. \(z = 1\).
Итак, одна сторона воздушного слоя имеет температуру \(t_0 = 0{,}0^\circ\text{C}\), вторая \(t_1 = 20{,}0^\circ\text{C}\). Поток теплоты стал стационарным, газ в равновесии.
2.3 Рассчитайте поток теплоты через воздушный слой \(q_a\), в предположении, что коэффициент теплопроводности воздуха постоянен (такой как при \(t_0 = 0{,}0^\circ\text{C}\)).
2.4 Нарисуйте схематический график зависимости температуры воздуха от координаты в установившемся режиме. (Здесь вам нет необходимости получать точную формулу, и даже вид зависимости).
2.5 Рассчитайте поток теплоты через воздушный слой, учитывая зависимость коэффициента теплопроводности от температуры.
2.6 Рассчитайте количество теплоты, которое пошло на нагревание воздуха от начальной температуры до установившегося распределения температур, при условии, что коэффициент теплопроводности зависит от температуры.

▶ПодсказкаСкрыть: подсказка
ПодсказкаЗакон Фурье аналогичен закону Ома: тепловые сопротивления пластин при последовательном соединении складываются, температура на стыке находится как «падение напряжения»; для воздуха теплопроводность растёт с температурой — разбейте слой на тонкие слои.
▶ОтветСкрыть: ответ
Ответ1.1) \(q_{Ag}=8{,}58\cdot10^{4}\) Вт, \(q_{Pb}=7{,}02\cdot10^{3}\) Вт; 1.2) \(Q_{Ag}=2{,}47\cdot10^{6}\) Дж, \(Q_{Pb}=1{,}43\cdot10^{6}\) Дж; 1.3–1.4) \(q_1=q_2=\left(\dfrac1{q_{Ag}}+\dfrac1{q_{Pb}}\right)^{-1}=6{,}49\cdot10^{3}\) Вт; 1.5–1.6) температура стыка \(18{,}5\,^\circ\)C (Pb на холодной стороне) и \(1{,}5\,^\circ\)C (Ag на холодной стороне), \(Q_1=6{,}07\cdot10^{6}\) Дж, \(Q_2=1{,}73\cdot10^{6}\) Дж; 2.1) \(\langle l\rangle=\dfrac1{\pi d^2n}\); 2.2) \(z=\dfrac12\); 2.3) \(q_a\approx4{,}80\) Вт; 2.5) \(q_a=\dfrac{\lambda_0S}{2T_0l}(T_1^2-T_0^2)\approx4{,}98\) Вт; 2.6) \(Q=\dfrac54P_0V\dfrac{\Delta T}{T_0}\approx915\) Дж.
▶РешениеСкрыть: решение
Решение1.1 Для расчетов потоков теплоты надо просто подставить численные значения в заданную формулу и провести расчет Для серебра
\[q_{Ag} = \lambda S \frac{\Delta T}{h} = 8{,}58 \cdot 10^{4} \text{ Вт}\tag{1}\]
Для свинца
\[q_{Pb} = \lambda S \frac{\Delta T}{h} = 7{,}02 \cdot 10^{3} \text{ Вт}\tag{2}\]
1.2 Так как коэффициент теплопроводности постоянен в установившемся режиме температуру изменяется по линейному законы. Поэтому при расчетах можно считать, что среднее значение изменения температуры равно \(\delta T = \frac{1}{2} \Delta T\). После этого опять считаем Для серебра
\[Q_{Ag} = \frac{1}{2} c m \Delta t = \frac{1}{2} c \rho S h \Delta t = 2{,}47 \cdot 10^{6} \text{ Дж}\tag{3}\]
Для свинца
\[Q_{Pb} = \frac{1}{2} c \rho S h \Delta t = 1{,}43 \cdot 10^{6} \text{ Дж}\tag{4}\]
1.3 Для решения этого пункта воспользуемся упорно навязываемой аналогией с законами постоянного тока. Так для каждой пластины можно ввести и рассчитать тепловое сопротивление
\[R = \frac{1}{\lambda} \frac{h}{S}\tag{5}\]
Так как пластины соединены последовательно (теплота проходит через одну, потом через другую), то общее сопротивление равно сумме сопротивлений. Поэтому поток теплоты через сдвоенные пластины будет равен
\[q_1 = q_2 = \frac{\Delta T}{\frac{h}{\lambda_{Ag}S} + \frac{h}{\lambda_{Pb}S}}\tag{6}\]
Этой формуле можно придать интересный вид
\[q_1 = q_2 = \left( \frac{1}{q_{Ag}} + \frac{1}{q_{Pb}} \right)^{-1}\tag{7}\]
Расчеты по этой формуле дают значение
\[q_1 = q_2 = 6{,}49 \cdot 10^3 \, \text{Вт}\tag{8}\]
Как видно поток определяется потоком через плохо проводящий слой.
Для расчета количества теплоты, потраченное на нагревание необходимо сначала найти распределение температуры, т.е. решить следующие пункты задачи.
1-4 – 1.6 Распределение температуры в каждом слое, по-прежнему, линейное. Поэтому достаточно найти температуру стыка. Опять воспользуемся аналогией с протеканием электрического тока. Обозначим свойства вещества примыкающего к холодной стороне индексом «х», а к горячей стороне индексом «г». Тогда можно записать (здесь \(T_z\) - температура стыка):
\[T_z - T_0 = q R_x = \frac{\Delta T}{\frac{h}{\lambda_{Ag}S} + \frac{h}{\lambda_{Pb}S}} \cdot \frac{h}{\lambda_x S} = \frac{T_1 - T_0}{\frac{1}{\lambda_x} + \frac{1}{\lambda_z}} \frac{1}{\lambda_x}\tag{9}\]
Откуда следует, что
\[T_z = T_0 + \frac{T_1 - T_0}{\frac{1}{\lambda_x} + \frac{1}{\lambda_z}} \frac{1}{\lambda_x}\tag{10}\]
Проведем расчеты. В первом случае (\(Pb \to Ag\))
\[T_{z1} = T_0 + \frac{T_1 - T_0}{\frac{1}{\lambda_{Ag}} + \frac{1}{\lambda_{Pb}}} \frac{1}{\lambda_{Pb}} = 18{,}5^\circ C.\tag{11}\]

Во втором случае (\(Ag \to Pb\))
\[T_{z2} = T_0 + \frac{T_1 - T_0}{\frac{1}{\lambda_{Ag}} + \frac{1}{\lambda_{Pb}}} \frac{1}{\lambda_{Ag}} = 1{,}5^\circ C.\tag{12}\]
Во втором случае можно было просто вычислить \(T_{z2} = \Delta T - T_{z1}\). Графики зависимостей температуры от координаты показаны на рисунке. Видно, что основное «падение
температуры» происходит на свинцовой пластине. Также по графику видно, что основное различие в количестве теплоты возникает из-за того, что серебряная пластина в первом случае нагревается больше.
Для расчетов теплот учтем, что пластина, примыкающая к холодной стороне, в среднем нагрелась на \(\frac{1}{2}T_Z\); а вторая - на величину \(\frac{1}{2}(T_z + T_1)\).
Поэтому количество теплоты, которой пойдет на нагревания будет равно:
В первом случае (\(Pb \to Ag\)):
\[Q_1 = \frac{1}{2}c_{Pb}\rho_{Pb}ShT_Z + \frac{1}{2}c_{Ag}\rho_{Ag}Sh(T_Z + T_1) = 6{,}07 \cdot 10^6 \text{ Дж}\tag{13}\]
Во втором случае (\(Ag \to Pb\))
\[Q_1 = \frac{1}{2}c_{Ag}\rho_{Ag}ShT_Z + \frac{1}{2}c_{Pb}\rho_{Pb}Sh(T_Z + T_1) = 1{,}73 \cdot 10^6 \text{ Дж}.\tag{14}\]
Часть 2. Воздушная пластина.
В данной части задачи все температуры следует рассчитывать в шкале Кельвина.
2.1 Для определения вида зависимости длины свободного пробега от указанных параметров, заметим, что движущаяся молекула заметает на своем пути площадь равную \(\pi d^2\). На длине свободно пробега в цилиндре с указанной площадью поперечного сечения и высотой равной длине свободного пробега в среднем должна находиться одна молекула. Поэтому можно записать
\[\pi d^2 \langle l \rangle n = 1 \quad \Rightarrow \quad \langle l \rangle = \frac{1}{\pi d^2 n}.\tag{15}\]
2.2 Подставляя найденную зависимость в формулу для теплопроводности, получим выражение
\[\lambda = \frac{1}{3} \langle v \rangle \langle l \rangle \rho c_V = \frac{1}{3} \sqrt{\frac{8kT}{\pi m}} \cdot \frac{1}{\pi d^2 n} \cdot nm \cdot C_V M.\tag{16}\]
Из которого следует, что коэффициент теплопроводности воздуха пропорционален корню квадратному из абсолютной температуры. Т.е. показатель степени \(z = \frac{1}{2}\). Важно также отметить, что теплопроводность не зависит от плотности (концентрации).
Экспериментальные данные показывают, что при температурах близких к комнатным показатель степени ближе к единице. Поэтому выбранная модель также имеет право на существование.
2.3 Используя формулу для потока теплоты, получаем, что поток теплоты через воздушную прослойку (в приближении постоянства теплопроводности) равен
\[\tilde{q}_a = \lambda S \frac{\Delta T}{h} = 4{,}80 \text{ Вт}\tag{17}\]

2.4 Так как теплопроводность воздуха возрастает с ростом температуры, то график температурного профиля будет выгнут вверх. Действительно, в тех местах, где теплопроводность больше, градиент температуры меньше.
2.5 Так как теплопроводность воздуха зависит от температуры, то она изменяется и при переходе от одной точки к другой. Поэтому закон теплопроводности Фурье следует записывать только для тонких слоев (т.е. в дифференциальной форме). Выразим зависимость теплопроводности от температуры посредством формулы
\[\lambda = \lambda_0 \frac{T}{T_0}.\tag{18}\]
И запишем уравнение Фурье для тонкого слоя воздуха \(\Delta x\):
\[q = \lambda_0 \frac{T}{T_0} \frac{\Delta T}{\Delta x} S.\tag{19}\]
Воспользуемся очевидной подсказкой и запишем
\[q \Delta x = \frac{\lambda_0 S}{2 T_0} \Delta(T^2)\tag{20}\]
Теперь легко просуммировать по всем слоям:
\[q l = \frac{\lambda_0 S}{2 T_0} (T_1^2 - T_0^2)\tag{21}\]
Откуда получаем окончательно для потока теплоты
\[q = \frac{\lambda_0 S}{2 T_0 l} (T_1^2 - T_0^2).\tag{22}\]
Расчет дает следующее численное значение: \(q_a = 4{,}98\) Вт, что несколько больше, чем получено ранее в приближенном расчете.
2.6 Так как объем воздуха не изменяется, то все полученная теплота пойдет на увеличение внутренней энергии воздуха. В установившемся режиме теплопередачи температура будет изменяться от точки к точке, постоянным будет оставаться давление! Поэтому выразим внутреннюю энергию воздуха через его давление:
\[U = \frac{5}{2} R \nu T = \frac{5}{2} P V.\tag{23}\]
Следовательно, изменение внутренней энергии воздуха (и равное ему количество полученной теплоты) выражается через изменение давления
\[Q = \Delta U = \frac{5}{2} V \Delta P.\tag{24}\]
Таким образом, нам необходимо рассчитать установившееся давление. Для этого воспользуемся условием постоянства числа молекул воздуха. Для этого разобьем воздушную пластину на тонкие слои \(\Delta x_i\). Концентрацию частиц в каждом слое можно выразить из уравнения состояние газа
\[n_i = \frac{P}{k T_i}.\tag{25}\]
Тогда общее число молекул можно получить, если просуммировать числа частиц во всех слоях
\[N = \sum_i \frac{P S}{k T_i} \Delta x_i.\tag{26}\]
Это же число можно выразить через начальные условия
\[N = \frac{P_0 S l}{k T_0}.\tag{27}\]
Таким образом, для определения давления получаем уравнение
\[P \sum _ { i } \frac { \Delta x _ { i } } { T _ { i } } = \frac { P _ { 0 } l } { T _ { 0 } } .\tag{28}\]
Прямое вычисление суммы требует расчета распределения температуры и последующего интегрирования. Можно поступить проще: выразим \(\Delta x _ { i }\) из уравнения теплопроводности (19):
\[q = \lambda _ { 0 } \frac { T } { T _ { 0 } } \frac { \Delta T } { \Delta x } S \quad \Rightarrow \quad \Delta x _ { i } = \lambda _ { 0 } \frac { T _ { i } } { T _ { 0 } } \frac { \Delta T _ { i } } { q } S\tag{29}\]
И после подстановки получим
\[\sum _ { i } \frac { \Delta x _ { i } } { T _ { i } } = \sum _ { i } \frac { 1 } { T _ { i } } \lambda _ { 0 } \frac { T _ { i } } { T _ { 0 } } \frac { \Delta T _ { i } } { q } S = \frac { \lambda _ { 0 } S } { q T _ { 0 } } \sum _ { i } \Delta T _ { i } = \frac { \lambda _ { 0 } S } { q T _ { 0 } } ( T _ { 1 } - T _ { 0 } ) .\tag{30}\]
Подставим найденное выражение для потока теплоты (22)
\[\sum _ { i } \frac { \Delta x _ { i } } { T _ { i } } = \frac { \lambda _ { 0 } S } { q T _ { 0 } } ( T _ { 1 } - T _ { 0 } ) = \frac { \lambda _ { 0 } S } { \displaystyle \frac { \lambda _ { 0 } S } { 2 T _ { 0 } l } ( T _ { 1 } ^ { 2 } - T _ { 0 } ^ { 2 } ) / T _ { 0 } } ( T _ { 1 } - T _ { 0 } ) = = \frac { 2 l } { ( T _ { 1 } + T _ { 0 } ) } .\tag{31}\]
Теперь из уравнения (28) находим
\[P \sum _ { i } \frac { \Delta x _ { i } } { T _ { i } } = \frac { P _ { 0 } } { T _ { 0 } } \quad \Rightarrow \quad P = P _ { 0 } \frac { T _ { 1 } + T _ { 0 } } { 2 T _ { 0 } }\tag{32}\]
Наконец, подставляем в формула для количества теплоты
\[Q = \frac { 5 } { 2 } V \Delta P = \frac { 5 } { 2 } P _ { 0 } V \left( \frac { T _ { 1 } + T _ { 0 } } { 2 T _ { 0 } } - 1 \right) = \frac { 5 } { 4 } P _ { 0 } V \frac { \Delta T } { T _ { 0 } }\tag{33}\]
Численное значение поглощенной теплоты равно \(Q = 915\) Дж.
Прочитали решение — спросите себя, как можно было до этого додуматься. А через день-два решите задачу заново, не подглядывая.
Похожие задачи
- Транзистор на пластине Всесоюзная 1985/86, заключительный этап, 9 кл.
- Остывание воды после кипятильника Белорусская 1989/90, заключительный этап, 9 кл.
- Пластина на солнце Белорусская 1990/91, заключительный этап, 9 кл.
- Два стержня Белорусская 1995/96, заключительный этап, 11 кл.
- Похолодание «Туймаада» 1995/96, заключительный этап, 11 кл.
- Остывание чайника МОШ 1997/98, заключительный этап, 8 кл.