Теории Юнга и Лапласа
В 1804 г. Томас Юнг [7] обосновал теорию капиллярных явлений на принципе поверхностного натяжения. Он также наблюдал постоянство угла смачивания жидкостью поверхности твердого тела (краевого угла) и нашел количественное соотношение, связывающее краевой угол с коэффициентами поверхностного натяжения соответствующих межфазных границ. В равновесии контактная линия не должна двигаться по поверхности твердого тела, а значит, говорил
(1)
где sSV, sSL, sLV - коэффициенты поверхностного натяжения межфазных границ твердое тело - газ (пар), твердое тело - жидкость, жидкость - газ соответственно, q - краевой угол. Это соотношение теперь известно как формула Юнга. Эта работа все же не оказала такого влияния на развитие науки в этом направлении, какое оказала вышедшая несколькими месяцами позже статья Лапласа (Pierre Simon Laplace). Это, по-видимому, связано с тем, что Юнг избегал использования математических обозначений, а пытался описывать все словесно, отчего его работа кажется запутанной и неясной. Тем не менее он считается сегодня одним из основателей количественной теории капиллярности.
Явления когезии и адгезии , конденсация пара в жидкость, смачивание твердых тел жидкостями и многие другие простые свойства вещества - все указывало на наличие сил притяжения, во много раз более сильных, чем гравитация, но действующих только на очень малых расстояниях между молекулами. Как говорил Лаплас, единственное вытекающее из наблюдаемых явлений условие, налагаемое на эти силы, состоит в том, что они «неощутимы на ощутимых расстояниях».
Силы отталкивания создавали больше хлопот. Их наличие нельзя было отрицать - они должны уравновешивать силы притяжения и препятствовать полному разрушению вещества, но их природа была совершенно неясной. Вопрос осложнялся двумя следующими ошибочными мнениями. Во-первых, часто считалось, что действующей силой отталкивания является тепло (как правило, мнение сторонников теории теплорода), поскольку (такова была аргументация) жидкость при нагревании сначала расширяется и затем кипит, так что молекулы разъединяются на гораздо большие расстояния, чем в твердом теле. Второе ошибочное мнение возникло из уводящего назад к Ньютону представления, согласно которому наблюдаемое давление газа происходит вследствие статического отталкивания между молекулами, а не из-за их столкновений со стенками сосуда, как тщетно доказывал Даниель Бернулли.
На этом фоне было естественно, что первые попытки объяснить капиллярность или вообще сцепление жидкостей основывались на статических аспектах вещества. Механика была хорошо понимаемой теоретической ветвью науки; термодинамика и кинетическая теория были еще в будущем. В механическом рассмотрении ключевым было предположение о больших, но короткодействующих силах притяжения. Покоящиеся жидкости (в капиллярной ли трубке или вне ее) находятся, очевидно, в равновесии, а потому эти силы притяжения должны уравновешиваться силами отталкивания. Поскольку о них можно было сказать еще меньше, чем о силах притяжения, их часто обходили молчанием, и, говоря словами Рэлея, «силам притяжения предоставлялось исполнять немыслимый трюк уравновешивания самих себя». Лаплас первым удовлетворительно разрешил эту проблему [8], полагая, что силы отталкивания (тепловые, как он допускал) можно заменить внутренним давлением, которое действует повсеместно в несжимаемой жидкости. (Это предположение приводит временами к неопределенности в работах XIX в. в отношении того, что строго понимается под «давлением в жидкости».) Приведем расчет внутреннего давления по Лапласу. (Этот вывод ближе к выводам Максвелла [2] и Рэлея [10]. Вывод приводится по [9] .)
Оно должно уравновешивать силы сцепления в жидкости, и Лаплас отождествлял это с силой на единицу площади, которая оказывает сопротивление разделению бесконечного жидкого тела на два далеко разъединяемых полубесконечных тела, ограниченных плоскими поверхностями. Приведенный ниже вывод ближе к выводам Максвелла и Рэлея, чем к оригинальной форме Лапласа, но существенного различия в аргументации нет.
Рассмотрим два полубесконечных тела жидкости со строго плоскими поверхностями, разделенные прослойкой (толщины l) пара с пренебрежимо малой плотностью (рис. 1), и в каждом из них выделим элемент объема. Первый находится в верхнем теле на высоте r над плоской поверхностью нижнего тела; его объем равен dxdydz. Второй находится в нижнем теле и имеет объем , где начало полярных координат совпадает с положением первого элементарного объема. Пусть f(s) - сила, действующая между двумя молекулами, разделенными расстоянием s, а d - радиус ее действия. Поскольку это всегда сила притяжения, имеем
Если r - плотность числа молекул в обоих телах, то вертикальная составляющая силы взаимодействия двух элементов объема равна
(2)
Полная сила притяжения, приходящаяся на единицу площади (положительная величина), есть
(3)
Пусть u(s) - потенциал межмолекулярной силы:
(4)
(5)
Рис. 1.
Интегрируя по частям еще раз, получаем
(6)
Внутреннее давление Лапласа K есть сила притяжения на единицу площади между двумя плоскими поверхностями при их контакте, т.е. F(0):
(7)
где - элемент объема, который можно записать как . Поскольку u(r) по предположению всюду отрицательно или равно нулю, то K положительно. Лаплас полагал, что K велико по сравнению с атмосферным давлением, но первую реалистическую численную оценку предстояло сделать Юнгу.
Приведенный выше вывод основан на неявном допущении, что молекулы распределены равномерно с плотностью r, т.е. жидкость не обладает различимой структурой в шкале размеров, соизмеримых с радиусом действия сил d. Без этого предположения нельзя было бы написать выражения (2) и (3) в такой простой форме, а надо было бы выяснить, как присутствие молекулы в первом элементе объема влияет на вероятность наличия молекулы во втором.
Натяжение на единицу длины вдоль произвольной линии на поверхности жидкости должно быть равным (в соответствующей системе единиц) работе, затраченной на создание единицы площади свободной поверхности. Это следует из опыта по растяжению пленки жидкости (рис. 2).
Рис. 2.
На проволочной рамке держится жидкая пленка, прикрепленная правым краем к свободно перемещаемой проволочке. Сила F, необходимая для уравновешивания натяжения в двусторонней пленке, пропорциональна длине L. Пусть F = 2sL. Смещение проволочки на расстояние x требует работы Fsdx = sdA, где dA - увеличение площади. Таким образом, натяжение на единицу длины на отдельной поверхности, или поверхностное натяжение s, численно равно поверхностной энергии на единицу площади.
Величина этой работы может быть сразу получена из выражения (6) для F(l). Если взять два полубесконечных тела в контакте и развести их на расстояние, превышающее радиус действия межмолекулярных сил, работа на единицу площади будет определяться как
(8)
При разделении образуются две свободные поверхности, и потому затраченную работу можно приравнять удвоенной поверхностной энергии на единицу площади, которая равна поверхностному натяжению:
(9)
Таким образом, K есть интеграл от межмолекулярного потенциала, или его нулевой момент, а H - его первый момент. В то время как K недоступно прямому эксперименту, H может быть найдено, если мы сможем измерить поверхностное натяжение.
Пусть - плотность когезионной энергии в некоторой точке жидкости или газа, т.е. отношение dU/dV где dU - внутренняя энергия малого объема V жидкости или газа, содержащего эту точку. Для молекулярной модели принимаем
(10)
где r - расстояние от рассматриваемой точки. Рэлей отождествлял лапласовское K с разностью этого потенциала 2 между точкой на плоской поверхности жидкости (значение 2S) и точкой внутри (значение 2I). На поверхности интегрирование в (10) ограничено полусферой радиуса d, а во внутренней области проводится по всей сфере. Следовательно, S есть половина I, или
(11)
Рассмотрим теперь каплю радиуса R. Расчет fI не изменяется, но при получении fS интегрирование теперь проводится по более ограниченному объему из-за кривизны поверхности. Если - угол между вектором и фиксированным радиусом , то
(12)
Тогда внутреннее давление в капле есть
(13)
где H определяется уравнением (9). Если бы мы взяли не сферическую каплю, а порцию жидкости с поверхностью, определяемой двумя главными радиусами кривизны R1 и R2 , то получили бы внутренне давление в виде
(14)
По теореме Эйлера сумма равна сумме обратных радиусов кривизны поверхности вдоль любых двух ортогональных касательных.
Так как K и H положительны и R положительно для выпуклой поверхности, то из (13) следует, что внутреннее давление в капле выше, чем в жидкости с плоской поверхностью. Наоборот, внутреннее давление в жидкости, ограниченной вогнутой сферической поверхностью ниже, чем в жидкости с плоской поверхностью, поскольку R в этом случае отрицательно.
Эти результаты составляют основу теории капиллярности Лапласа. Уравнение для разности давлений (давление жидкости внутри сферической капли радиуса R) и (давление газа снаружи) теперь называют уравнением Лапласа:
(15)
Достаточно трех идей - натяжения у поверхности, внутреннего давления и краевого угла, а также выражений (1) и (15), чтобы решить все задачи обычной равновесной капиллярности методами классической статики. Таким образом, после работ Лапласа и Юнга основы количественной теории капиллярности были заложены.
Результаты Юнга были получены позже Гауссом вариационным методом. Но все эти работы (Юнга, Лапласа и Гаусса) обладали одним общим недостатком, изъяном, если можно так выразиться. Об этом недостатке будет рассказано позже.
При расчете давления внутри искривленной жидкой поверхности был введен потенциал Рэлея 2 (10); попутно было отмечено, что I является плотностью когезионной энергии. Впервые это полезное понятие в 1869 г. ввел Дюпре, который определил его как работу дробления куска вещества на составляющие его молекулы (la travail de dйsagrйgation totale - работа полной дезагрегации).
Рис. 3
Направленная внутрь сила, действующая на молекулу на глубине r < d, противоположна по знаку направленной наружу силе, которая бы возникла со стороны молекул в заштрихованном объеме, если бы он был заполнен равномерно с плотностью .
Он приводит [12] вывод, проделанный его коллегой Ф. Ж. Д. Массье следующим образом. Сила, действующая на молекулу у поверхности по направлению к объему жидкости, противоположна по знаку силе, возникающей от заштрихованного объема на рис. 3, поскольку внутри жидкости сила притяжения от шарового объема радиуса равна нулю из симметрии. Таким образом, сила, направленная внутрь, есть
(16)
Эта сила положительна, так как f(0 < s < d) < 0 и F(d) = 0 из-за нечетности функции f(s). Никакая сила не действует на молекулу, если только она не находится в пределах расстояния d по ту или иную сторону от поверхности. Следовательно, работа удаления одной молекулы из жидкости равна
(17)
поскольку u(r) - четная функция. Эта работа равна минус удвоенной энергии на молекулу, необходимой для дезинтеграции жидкости (удвоенной, чтобы не считать молекулы дважды: один раз при их удалении, другой раз - как часть среды):
(18)
Это простое и понятное выражение для внутренней энергии U жидкости, содержащей N молекул. Отсюда следует, что плотность когезионной энергии дается выражением (10), или
(19)
что совпадает с (11), если убрать индекс I. Сам Дюпре получил тот же результат окольным путем. Он рассчитывал dU/dV через работу против межмолекулярных сил при однородном расширении куба жидкости. Это дало ему
(20)
Поскольку K имеет форму ((7) и (11)), где постоянная a дается выражением
(21)
то интегрирование (20) снова приводит к (19).
Рэлей критиковал вывод Дюпре [10]. Он считал, что рассмотрение работы однородного расширения от состояния баланса когезионных и отталкивающих межмолекулярных сил при учете только когезионных сил было необоснованным; прежде чем предпринять подобный шаг, следовало бы располагать лучшим знанием вида сил отталкивания.
Мы видим, что в этом выводе, как и в выводах Юнга, Лапласа и Гаусса, существенным образом используется предположение о скачкообразном изменении плотности числа молекул вещества на границе раздела фаз. В то же время, чтобы проведенные рассуждения описывали реальные явления в веществе, необходимо предполагать, что радиус действия межмолекулярных сил в веществе много больше характерного расстояния между частицами. Но при этом предположении граница раздела двух фаз не может быть резкой - должен возникнуть непрерывный переходный профиль плотности, иначе говоря, переходная зона.
Были предприняты попытки обобщить эти выводы на непрерывный переходный профиль. В частности, Пуассон, пытаясь пойти по такому пути, пришел к ошибочному выводу, что при наличии переходного профиля поверхностное натяжение должно вообще исчезнуть. Позже Максвелл показал ошибочность такого вывода.
Однако, само предположение о том, что радиус действия межмолекулярных сил в веществе много больше характерного расстояния между частицами не соответствует экспериментальным данным. В действительности, эти расстояния одного порядка. Поэтому механистическое рассмотрение в духе Лапласа является, говоря современным языком, теорией среднего поля. Таковой же является не описанная здесь теория Ван-дер-Ваальса, давшая знаменитое уравнение состояния реальных газов. Во всех этих случаях точный расчет требует учета корелляций между плотностями количества частиц в различных точках. Это делает задачу очень сложной.
К 1819 г. он был занят детальным обсуждением межмолекулярных сил отталкивания, которые, хотя и приписывались еще теплоте или теплороду, обладали существенным свойством уменьшаться с расстоянием быстрее, чем силы притяжения.
Юнг упоминал наличие градиента плотности в конечном по толщине слое, но отбросил этот эффект, посчитав его несущественным.