Статьи Соросовского Образовательного журнала в текстовом формате
Рассмотрен термодиффузионный эффект: экспериментальные результаты, описание с помощью теории Онсагера, физический смысл тепла переноса, некоторые применения.
ТЕРМОДИФФУЗИЯБ. С. БОКШТЕЙН
Московский институт стали и сплавов
ВВЕДЕНИЕ
Термодиффузионный эффект был открыт Л. Дюфуром в 1873 году. Он открыл обратный, диффузионный термоэффект: при взаимной диффузии водорода и воздуха, а также воздуха и диоксида углерода через пористую перегородку возникала разность температур. Прямой термодиффузионный эффект (для газов) заключается в том, что при наличии постоянной разницы температур по обе стороны перегородки, разделяющей один и тот же газ, возникает, растет и становится постоянной разница давлений, даже если вначале ее не было. Таким образом, разница (градиент) температур приводит к разнице (градиенту) концентраций. Это и есть термодиффузия.
Позднее термодиффузионный эффект был обнаружен и измерен для конденсированных систем - жидких (эффект СорО) и твердых. Примером может служить образец сплава железа с углеродом, содержащий вначале 0,01% С по всей длине (рис. 1, а). Если поместить его в температурное поле, так что на левом конце образца поддерживается постоянная температура 710?С, а на правом - 670?С (рис. 1, б ), то концентрация углерода на левом, более горячем конце начнет расти, а на правом, более холодном - уменьшаться, пока не установится стационарное состояние, изображенное на рис. 1, в [1].
В указанном интервале температур образец однофазен и представляет собой твердый раствор углерода в железе, имеющий объемно-центрированную кубическую (ОЦК) решетку a-Fe. Это значит, что весь металлический образец состоит из повторяющихся кубиков (ячеек), в которых атомы железа расположены в вершинах и центре каждого куба, а атомы углерода - между ними, в так называемых междоузлиях (подробнее см. [2]). Длина ребра такого куба (период решетки) для чистого a-Fe составляет 0,286 нм. Углерод (0,01%) практически не меняет решетки.
Если нагреть образец выше 910?С, то решетка меняется: a-Fe превращается в g-Fe. Теперь атомы железа расположены в вершинах куба и центрах каждой грани. Такая решетка называется ГЦК (гранецентрированная кубическая). Ее период равен 0,365 нм. Если изучать термодиффузию углерода в g-Fe, то его концентрация в температурном поле не меняется, термодиффузионного эффекта нет.
Эти экспериментальные результаты непонятны и, по крайней мере на первый взгляд, противоречат законам переноса вещества или диффузии. Согласно первому закону Фика [3], диффузионный поток пропорционален градиенту концентрации, то есть вещество перемещается из области, где его больше, в область, где его меньше. Однако в описанных экспериментах давление газа было одинаково по обе стороны перегородки, равно как была одинакова концентрация углерода по длине железного образца. Тем не менее после "включения" температурного поля вещество (газ, углерод) начинало перемещаться, хотя диффузионный поток не должен, казалось бы, зависеть от разницы (градиента) температур. Также непонятно, почему небольшое (опять, казалось бы) изменение решетки железа так кардинально меняет результат: был поток и нет его. В чем же дело?
ПЕРЕКРЕСТНЫЕ ЭФФЕКТЫ
В первоначальной формулировке законов переноса потоки того или иного свойства пропорциональны градиенту этого свойства. Так, поток электрического заряда (электрический ток) пропорционален градиенту электрического потенциала (напряженности электрического поля) - закон Ома; поток тепла - градиенту температуры - закон Фурье; поток вещества - градиенту концентрации - закон Фика.
Однако обнаружено множество перекрестных эффектов, когда градиент одного свойства вызывает градиент и, следовательно, поток другого свойства. Наиболее известными среди них являются термоэлектрические эффекты: эффект Пельтье - появление разницы температур при протекании электрического тока через систему, состоящую из разнородных проводников, и термо-ЭДС - возникновение электродвижущей силы в системе из разнородных металлов, находящихся при разных температурах.
В 1968 году Ларсу Онсагеру была присуждена Нобелевская премия по химии за "формулировку соотношений взаимности... имеющих важное значение для термодинамики необратимых процессов". Термодинамика необратимых процессов справедлива при описании процессов в открытых системах (то есть в системах, обменивающихся теплом, работой и веществом с окружающей средой) при не слишком больших отклонениях от термодинамического равновесия [4] и при неизменных внешних условиях. В таких системах в конце концов устанавливается стационарное состояние - свойства системы перестают зависеть от времени [5].
В основе теории Онсагера лежат следующие положения.
1. Поток i-го свойства (ji ) линейно зависит от всех термодинамических сил, действующих в системе
где Xk - термодинамические силы, то есть градиенты свойств, а Lik - кинетические коэффициенты или коэффициенты Онсагера.
Кинетические коэффициенты могут быть диагональными Lii , отражающими влияние "своей" силы на "свой" поток (именно они содержались в законах Ома, Фурье и Фика), и перекрестными Lik (i ? ? k), отражающими влияние "чужой" силы на "свой" поток. Уравнения (1) иногда называют термодинамическими уравнениями движения, они являются обобщением опыта, который привел к законам Ома, Фурье и др.
2. Правило выбора термодинамических сил. Оно позволяет рассчитать термодинамические силы, действующие в системе. В случае векторных потоков эти силы всегда оказываются связанными с градиентами свойств: тепловая сила где WT - градиент температуры; диффузионная (для диффузии i-го компонента в разбавленном растворе) Xi = - RT W ln Ci , где WCi - градиент концентрации диффундирующего вещества, а R - газовая постоянная; электрическая Xe = - eWj, где e - заряд иона и Wj - градиент электрического потенциала (напряженность электрического поля E = - Wj) и т.д.
3. Соотношения взаимности или правило симметрии кинетических коэффициентов:
Lik = Lki .
Оно показывает, что k-я термодинамическая сила влияет на поток i-го свойства так же, как i-я сила на поток k-го свойства. В отличие от первых двух положений, являющихся постулатами, обобщением опыта, третье может быть выведено из принципов статистической механики, так называемого принципа микроскопической обратимости.
ОПИСАНИЕ ТЕРМОДИФФУЗИОННОГО ЭФФЕКТА
Теория Онсагера позволяет получить уравнение, описывающее термодиффузионный эффект, по крайней мере формально. Чтобы воспользоваться формулой (1), нужно ответить на вопрос, какие потоки возникают в нашей системе. Ответ представляется более или менее очевидным: поток вещества (углерода) j1 и поток тепла j2 . В соответствии с формулой (1)
j1 = L11X1 + L12X2 ,
j2 = L21X1 + L22X2 .
Подставляя в уравнения (3) и (4) выражения для термодинамических движущих сил и учитывая уравнение (2), то есть L12 = L21 , получим
В уравнении (6) мы учли, что L12 = L21 .
Из уравнения (5) следует, что поток вещества (углерода) состоит как бы из двух частей: j1 = j1C + j1T , где j1C = - L11RT W ln C - концентрационная часть, зависящая только от WC, а - тепловая часть, зависящая только от WT.
Поскольку мы берем образец (рис. 1, а), в котором углерод распределен равномерно по длине и то вначале возникает только тепловая часть потока углерода j1T . Мы не знаем знак L12 и не знаем заранее, в какую сторону направлен поток j1T . Поскольку мы знаем ответ, то предположим, что поток j1T направлен к более горячему концу.
Итак, под действием тепловой части потока j1T концентрация углерода на более горячем конце образца возрастает и становится больше, чем на более холодном. Следовательно, возникает градиент концентрации. Теперь WC ? 0 и соответственно j1C ? 0. Очевидно, что эта (концентрационная) часть потока направлена в обратную сторону, от горячего конца, где углерода больше, к холодному. Если | j1C | < < | j1T | , то концентрация углерода на горячем конце продолжает расти. Соответственно растет поток j1C до тех пор, пока не сравняется с j1Т . Вспомним, что потоки j1C и j1Т направлены в разные стороны. Поэтому если они сравнялись по величине, то суммарный поток углерода j1 = 0. Установилось стационарное состояние, концентрация углерода во всех точках перестает зависеть от времени.
Из уравнения (5) и условия j1 = 0 следует, что в стационарном состоянии
Из сравнения размерностей левой и правой частей уравнения (7) становится ясным, что отношение L12 / L11 имеет размерность энергии. Легко понять формальный смысл этой энергии, если разделить уравнение (6) на уравнение (5), приняв, что WT = 0. Тогда
то есть это отношение имеет смысл энергии, переносимой движущейся частицей (атомом углерода, молекулой газа) в направлении своего движения, но в изотермических условиях, когда WT = 0 и Т = const. Эту величину принято называть теплом переноса Q*. Таким образом, в стационарном состоянии
Из результатов эксперимента с термодиффузией в a-Fe ясно, что Q* < 0, а зная численные результаты, можно рассчитать по прямой на рис. 1, в, что Q* = - 96 кДж/моль.
Теплоты переноса разных примесей различаются и по величине, и по знаку даже для одного растворителя, тем более для разных. Экспериментальные значения колеблются от -121 кДж/моль (сурьма в серебре) до + 84 кДж/моль (кислород в b-цирконии).
ФИЗИЧЕСКИЙ СМЫСЛ ЧИСЛА ПЕРЕНОСА
Уравнение (9) позволяет предсказать, как перераспределяется вещество в температурном поле, если известны знак и величина тепла переноса. Но как их узнать? Для этого нужны физические модели. Рассмотрим несколько примеров.
В простой модели, описывающей термодиффузионный перенос атомов в твердом теле (углерода в железе), предполагается, что, переходя из одного узла (или междоузлия) в решетке в соседнее, атом затрачивает энергию в трех местах: уходя из узла (Н1), на середине пути (Н1/2) и занимая новое положение (Н2). Сравнение формулы (9) с результатом, полученным в этой модели [5], показывает, что Q* = Н1 - - Н2 . Очевидно, что при обычной изотермической диффузии [3] энергия активации процесса составила бы Е = Н1 + Н1/2 + Н2 . Для диффузии С в a-Fe E . . 100 кДж/моль. Сравнивая Q* = - 96 кДж/моль и Е . 100 кДж/моль, естественно предположить, что Н2 @ Н1 и Н2 @ Н1/2 , то есть основную энергию атом затрачивает в "месте приземления".
Напротив, для термодиффузии в g-Fe мы получили Q* © 0, а для диффузии С в g-Fe энергия активации Е . 135 кДж/моль. Естественно предположить, что либо Н1 © Н2 , либо Н1/2 @ Н1 и Н1/2 @ Н2 , то есть основную энергию при движении в более плотной ГЦК-решетке g-Fe-атом затрачивает посередине пути, расчищая себе дорогу и расталкивая соседей. О других попытках связать тепло переноса со свойствами системы можно прочитать в [5].
Для упомянутого выше переноса газа через пористую перегородку можно предположить, что молекулы растворяются в перегородке с одной стороны (при этом выделяется теплота растворения) и выделяются с другой (теплота поглощается). Таким образом, тепло переноса равно теплоте растворения газа в перегородке. Если в перегородке имеются узкие каналы (поры, капилляры), диаметр которых существенно меньше длины свободного пробега молекул, то при движении в перегородке молекулы не сталкиваются друг с другом, а только со стенками канала. В результате молекула как бы теряет одну поступательную степень свободы и связанную с ней энергию RT /2 (в расчете на один моль), входя в перегородку, и приобретает ее, выходя. Таким образом, Q* = - RT /2. Подставляя это значение в формулу (9) и заменяя в ней концентрацию С на давление пара р, получим после интегрирования:
Таким образом, отношение давлений газа в стационарном состоянии по обе стороны перегородки равно корню квадратному из отношения температур. Этот результат хорошо известен из опыта.
Рассмотрим перенос через перегородку, пропускающую пар и не пропускающую жидкость. Такая перегородка может быть сделана из пористого несмачиваемого жидкостью материала. Погрузим в жидкость трубку, на дне которой находится такая перегородка. Тогда жидкость, чтобы попасть в трубку, должна испаряться по одну сторону перегородки и конденсироваться по другую. Поэтому теплота переноса равна теплоте испарения жидкости. Если температура жидкости больше температуры пара в трубке, начнется термодиффузионный перенос жидкости через перегородку, то есть будет создан своеобразный тепловой насос. Расчет показывает, что в случае воды при разнице температур в 200?С процесс подъема прекратится, когда высота жидкости в трубке составит несколько километров. На такую высоту тепловой насос способен поднимать жидкость.
Необходимо еще раз отметить роль перекрестных эффектов. Термодиффузия - пример того, как разница температур приводит к переносу вещества. Есть много примеров практического применения термодиффузионного эффекта. Выше был описан тепловой насос. С помощью термодиффузии можно разделять вещества: различные вещества имеют разные теплоты переноса. Например, при разработке атомной бомбы было необходимо разделить изотопы урана - эта проблема была решена с помощью термодиффузии.
ЛИТЕРАТУРА
1. Шьюмон П. Диффузия в твердых телах. М.: Металлургия, 1966.
2. Шаскольская М.И. Кристаллография. М.: Высш. шк., 1982.
3. Бокштейн Б.С. Почему и как движутся атомы в металлах // Соросовский Образовательный Журнал. 1995. ╧ 1.
4. Маневич Л.И. Феноменология необратимости // Там же. 1997. ╧ 11.
5. Бокштейн Б.С. Диффузия в металлах. М.: Металлургия, 1978.
* * *
Борис Самуилович Бокштейн, доктор физико-математических наук, профессор, зав. кафедрой физической химии Московского государственного института стали и сплавов. Область научных интересов - физика и химия диффузии и дефектов в материалах. Автор около 200 статей и семи книг.