Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

THE THEORY OFFLUID HEAT CV

Pavlov A.M. 1
1 S. Amanzholov East Kazakhstan State University
This paper proposes a formula for calculating the heat capacity Cv monatomic liquids. Formula is derived using the formula of the potential energy W molecules in the liquid, taken from [7]. Decomposing the W in a series in powers of y = r/a, managed to get a quasi-rigid molecular «springs» and calculate the dependence of the oscillation frequency on the relationship of its own volume of the molecule to the volume per particle. Besides the terms making up the potential energy of the oscillating particles in the expansion has the energy W(0), which the particle is in equilibrium. This energy, according to our calculations, makes a significant contribution to the heat capacity of the liquid. It depends on the coordination number f, and the heat capacity depends on f and df/dT. Calculation of these values was made on the assumption that the cluster is a liquid drop; and the surface tension acts on the surface of the drop. Moreover, the surface energy is associated with W(0).
heat capacity
interaction energy
fluctuations
frequency
1. Belaja M.L., Levadnyj V.G. Molekuljarnaja struktura vody // Znanie. 1987. № 3, ser. «Fizika».
2. Vasserman A.A., Kazavchinskiy Ya.Z., Rabinovich V.A. Teplofizicheskie svojstva zhidkogo vozduha i ego komponentov. M. Izd. komiteta standartov, mer i izmeritel’nyh priborov pri Sov. min. SSSR, 1968. 239 р.
3. Kikoin I.K., Kikoin A.K. Molekuljarnaja fizika. M.: Gos. izd. fiz. mat. lit., 1963. 500 р.
4. Krokston K. Fizika zhidkogo sostojanija. Statisticheskoe vvedenie. M.: Mir, 1978. 400 р.
5. Sakenova R.E., Pavlov A.M. O teorii teplojomkosti zhidkostej. // Sb. materialov I Rossijsko-Kazahstanskoj molodjozhnoj nauchno-tehnicheskoj konferencii. Barnaul 2013. рр. 51–56.
6. Frenkel’ Ja.I. Kineticheskaja teorija zhidkosti. L.:Nauka, 1975. 592 р.
7. Cjao Sjuje-Sen’. Fizicheskaja mehanika. M.: Mir, 1965. 544 р.

Обычно в литературе, например в [6], утверждается, что теплоёмкость жидкости близка к 3R, как у твёрдых тел. Однако таблицы теплоёмкости Cp, приведенные в [2], показывают, что эта теплоёмкость далека от 3R (в связи с малой расширяемостью жидкости Cv ≈ Cp). Так, теплоёмкость жидкого аргона при давлении 20 бар и температуре 85°К равна 4,7R, а при Т = 120°К – 6,32R, а при давлении 40 бар и Т = 145°К теплоёмкость Ср достигает значения 11,45R. Оказалось, что теплоёмкость жидкости зависит как от Т, так и от Р. При повышении температуры при том же давлении Р теплоёмкость увеличивается, а при повышении давления при Т = const она уменьшается.

На мой взгляд, изменение теплоёмкости связано с изменением числа ближайших соседей, т.е. координационного числа. В статье [5] была получена формула теплоёмкости Сv с учетом изменения координационного числа. Изменение этого числа при изменении температуры, т.е. df/dT, подсчитывалось с использованием идеи Я.И. Френкеля о разбиении жидкости на Z областей, каждая из которых содержит f частиц. Подсчитав увеличение энтропии и свободной энергии, там была получена следующая формула для f:

pavlov01.wmf, (1)

где US – энергия единицы поверхности; n – число частиц в единице объёма. Однако полученная из (1) формула для df/dT давала завышенные результаты. В данной статье сделана попытка найти df/dT с учетом другой связи US и f, а также зависимости US от σ и Т. Зависимость σ(Т), приведенная в [3] не дала желаемого результата. Поэтому использовалась формула σ(Т), взятая в [4].

Вывод формулы для Сv

В [7] приводится следующая формула для потенциальной энергии молекулы жидкости, окруженной f соседними молекулами:

pavlov02.wmf (2)

где rm – равновесное расстояние в потенциале Леннарда-Джонса; «a» – равновесное расстояние между молекулами в жидкости, ε – глубина потенциальной ямы в потенциале Леннарда-Джонса. Эта формула получена в результате усреднения потенциала взаимодействия частицы, находящейся в центре сферы, на которой находится f молекул, по всем направлениям колебания выбранной частицы. Если круглые скобки разложить в ряд по степеням pavlov03.wmf, то, учитывая малость y, получим

pavlov04.wmf (3)

где pavlov05.wmf (4)

Выражение в квадратной скобке можно представить в ином виде:

pavlov06.wmf (5)

Как известно, pavlov07.wmf – есть коэффициент жёсткости молекулярных квазипружин:

pavlov08.wmf (6)

с помощью которого определяется частота колебаний pavlov09.wmf

Хотя наличие слагаемого с r4 указывает на то, что колебания молекул нелинейные с частотой ω ≠ ω0, тем не менее будем считать, что все молекулы колеблются с одинаковой частотой ω = ω0, т.е. воспользуемся моделью А. Эйнштейна.

В этом случае свободная энергия, приходящаяся на одну частицу, имеет вид

pavlov10.wmf (7)

Тогда, учитывая, что U = F + TS и pavlov11.wmf, для внутренней энергии, приходящейся на одну частицу, получаем:

pavlov12.wmf (8)

Дифференцируя (8) по Т с учётом того, что W(0) зависит от координационного числа f, являющегося функцией температуры, получим:

pavlov13.wmf (9)

В W(0) и С можно сделать замену: pavlov14.wmf – собственный объём молекулы, a3 = V1 – объём, приходящийся на одну частицу. В этом случае

pavlov15.wmf (10)

pavlov16.wmf (11)

Объём V1 находится из таблиц, например [2] по удельному объёму:

pavlov17.wmf (12)

где NA – число Авогадро.

Зависимость координационного числа от температуры

В статье [5] для получения такой зависимости, как выше было сказано, использовалась идея Я.И. Френкеля, которая дала завышенное значение теплоёмкости. Здесь предлагается использовать зависимость коэффициента поверхностного натяжения от температуры. Будем исходить из того, что жидкость кластеризована и кластер представляет собой капельку жидкости. Если энергию одной частицы, находящейся на поверхности, обозначить через W, то энергия всей поверхности будет

U = WNпов,

где Nпов – число поверхностных молекул. Энергия единицы поверхности равна

pavlov18.wmf

где Rc – радиус кластера. Пусть эффективный радиус молекулы (радиус Вигнера ‒ Зейтца) равен R0. Тогда

Rc = R0n 2/3,

где n – число частиц в кластере.

Следовательно

pavlov19.wmf

но 4n 2/3 = Nпов. Поэтому pavlov20.wmf

Энергия молекулы на поверхности и энергия молекулы, находящейся внутри кластера, связаны некоторым образом. Будем считать, что энергия поверхностной молекулы составляет часть α от энергии внутренней частицы. Тогда

pavlov21.wmf

Дифференцируя это равенство по Т, нетрудно получить:

pavlov22.wmf (13)

Согласно [4]

pavlov23.wmf (14)

а

pavlov24.wmf (15)

где Т0 – критическая температура; σ – коэффициент поверхностного натяжения; σ0 и μ – константы, причем μ = 1,24.

Используя равенства (13)‒(15), можно определить df/dT. Так как

pavlov25.wmf

то pavlov26.wmf

или

pavlov27.wmf (16)

Поскольку US пропорционально f, то это же выражение будет определять и df/dT:

pavlov28.wmf

или, подставив μ = 1,24, можно записать

pavlov29.wmf (17)

После интегрирования получается следующее выражение для f:

pavlov30.wmf (18)

Постоянную С можно определить, если известно f при какой-либо температуре. Согласно [1] координационное число у аргона равно 10,3. Хотя температура там не указана, будем считать Т = 85°К. Тогда получается С = 11,18 и при Т = 120°К f = 9,06, а при 140°К f = 7,15. При 85°К (17) оказалось равным 0,0023, т.е. величиной, существенно меньшей единицы. Это позволяет считать, что df/dT практически не сказывается на частоте колебаний.

Вычисление теплоёмкости CV с помощью (9) и (17)

Чтобы воспользоваться формулой (9), необходимо определить частоту колебаний, которая зависит от отношения pavlov31.wmf и от координационного числа. Используя (12), найдем частоту колебаний:

pavlov32.wmf

Обозначим диаметр молекулы буквой D. Тогда V* = D3, pavlov33.wmf и pavlov34.wmf. При 85°К и Р = 20 бар получаем: ω0 = 3,847∙1012; pavlov35.wmf; pavlov36.wmf; f = 10,3. Тогда первое слагаемое в (9) равно 1,42, а второе слагаемое ‒ 2,975. Следовательно, теплоёмкость CV = 4,39R.

Произведем эти же вычисления при Р = 20 бар и Т = 120°К. Первое слагаемое в (9) равно 2,64, а второе – 3,044. Таким образом, CV = 5,68R.

Табличное значение CP при этих температурах равно 4,7R и 6,3R соответственно. Таким образом, при 85°К погрешность вычислений составляет 6,5 %, а при 120°К – 9,8 %. Это в том случае, если CV ≈ CP. На самом деле CV несколько меньше, чем СР и, следовательно, погрешность будет меньше. При Р = 40 бар и Т = 140°К наши вычисления дают CV = 8,35R, а табличное значение CV = 9R, т.е. погрешность равна 7,2 %.

Погрешность формулы, возможно, определяется неточностью вычисления координационного числа. Так, если при 85°К считать число соседних молекул равным 12, то постоянная С в (18) будет иметь значение 12,9 и координационное число при Т = 120°К и Т = 140°К увеличивается до 10,47 и 8,25 соответственно. Погрешность формулы (9) при этом сокращается. Поэтому определение координационного числа в жидкости и его зависимости от температуры остаётся актуальной задачей. Кроме этого, уменьшить погрешность можно, подбирая число μ.

Выводы

  1. Вклад колебательного движения молекул жидкости в теплоёмкость составляет 3R, если частота колебаний не зависит от температуры.
  2. Существенный вклад в теплоёмкость жидкостей даёт потенциальная энергия взаимодействия молекул. При низких температурах он составляет более 1,5R, а при высоких температурах – от 3R до 6R и более.
  3. Разработанный способ подсчёта зависимости координационного числа от температуры даёт удовлетворительные результаты. Погрешность подсчёта теплоёмкости с помощью рассчитанных значений f и df/dT не превышает 10 %.

Рецензенты:

Плотников С.В., д.ф.-м.н., профессор, заведующий кафедрой «Техническая физика», Восточно-Казахстанский государственный технический университет им. Д. Серикбаева, г. Усть-Каменогорск;

Алонцева Д.Л., д.ф.-м.н., профессор кафедры «Программирование и автоматизация технологических процессов», Восточно-Казахстанский государственный технический университет им. Д. Серикбаева, г. Усть-Каменогорск;

Бичурин М.И., д.ф.-м.н., профессор, зав. кафедрой «Проектирование и технология радиоаппаратуры», Новгородский государственный университет, г. Великий Новгород.

Работа поступила в редакцию 30.04.2014.