Методика расчёта диэлектрической проницаемости метанола и этанола в жидком состоянии

Канд. техн. наук, доц. МУЛ ЕВ Ю. В.

Белорусский ордена Трудового Красного Знамени политехнический институт

Докт. техн. наук СМИРНОВ С. Н.

Московский ордена Ленина и ордена Октябрьской Революции энергетический институт

 

Метанол и этанол являются добавками к моторному топливу. Контроль состояния их параметров — актуальный вопрос в технологии топ- ливопроизводства.

Безынерционный контроль параметров метанола и этанола может быть организован на основе измерения их диэлектрических проницаемостей.

Опубликованные данные по диэлектрической проницаемости метанола и этанола в жидком состоянии [1—3] получены для отдельных точек. Зависимости, описывающие эти результаты, отсутствуют.

Для теоретического описания диэлектрической проницаемости ε полярных веществ (метанол и этанол и др.) можно применить формулу Кирквуда [4]

 

  (ε-1)(2ε+1)= NAρ3M (α+μ2g3kT),     (1)

 

где NA — число Авогадро;

  ρ — плотность;

  М — молекулярная масса;

  α — поляризуемость смещения;

  k — постоянная Больцмана;

  T — температура.

  В (1) дипольный момент молекулы для расчетов принимается постоянным, а изменение его величины с температурой и давлением учитывается фактором корреляции g. Получить его значение при изменении плотности и температуры полярных веществ теоретическими расчетами невозможно [5].

  В современной молекулярно-статистической теории [6, 7] диэлектрическую проницаемость связывают с параметром у, определяющим вклад ориентационной поляризации,

 

  y=4/9πρNμ2MkT,    (2)

 

  Так, например, в [6] предложено

 

  ε-1ε+2=y-3y2+…  .    (3)

 

  Результаты расчета по (3) для метанола, как показано на рис. 1, меньше экспериментальных при малых у и значительно превышают εэкс

Рис. 1. Зависимость диэлектрической проницаемости метанола и этанола от параметра у:
  ——— - расчет ε метанола по (4); — — — — диэлектрическая проницаемость метанола по (3);
— · — —  диэлектрическая проницаемость метанола по (6); ▄ — экспериментальные точки для метанола; ▐ — для этанола

при y> 1,7. Большие расхождения с опытными данными связаны со сложностью нахождения потенциала взаимодействия, который у полярных диэлектриков анизотропен. Для метанола и этанола, являющихся к тому же и ассоциированными веществами, определение даже энергии водородной связи представляет большую проблему [5].

  Цель настоящей работы — создание методики, позволяющей рассчитывать диэлектрическую проницаемость метанола и этанола в жидком состоянии, на основе которой возможны интерполяция и экстраполяция с погрешностью, сопоставимой с точностью существующих экспериментальных данных.

  Поставленная цель достигается анализом экспериментальных зависимостей диэлектрических проницаемостей метанола и этанола от параметра у. Для их построения взяты результаты измерения диэлектрической проницаемости метанола из [1]. Их диапазон составил: по давлению 0,1÷250 МПа, по температуре 223,1÷298,1 К (всего около 50 экспериментальных значений). Нами использованы значения диэлектрической проницаемости метанола из [3], где приведены 14 точек для температур от 298 до 513 К- Для принятых экспериментальных значений ε метанола, температурный интервал составил 20 К, по давлению — 50 МПа.

Для этанола экспериментальные данные по диэлектрической проницаемости наиболее полно представлены в [2]. Температурный диапазон составил 163÷333 К и по плотности — 0,74÷0,9 г/см3 (всего около 20 экспериментальных точек). В [3] описаны восемь значений ε этанола для температур 303÷513 К. Для метанола значения дипольных моментов приняты равными 1,68, для этанола 1,74D).

Значения диэлектрической проницаемости метанола и этанола [1—3] заимствованы только для представленных табличных величин. Для междушагового пространства в последующих расчетах они исключались. Значения плотности представлены в описанных источниках.

Таблица 1

Результаты аппроксимации зависимостей диэлектрических проницаемостей от параметра у для метанола и этанола

  В результате расчетов зависимостей диэлектрических проницаемостей метанола и этанола от параметра у на основе массива из [1—3] построены зависимости (рис. 1). Очевидно, что нанесенные экспериментальные точки образуют в большей части линейные участки. Предельное значение параметра у, полученное в расчете, для рассмотренных веществ равно 2,02. Для этанола нижнее значение параметра у — 0,8. Для меньших его значений экспериментальные данные отсутствуют. Для метанола при y<0,6 наблюдается отклонение точек от линейной зависимости. Это свидетельствует о существовании при таких значениях у паровой фазы.

Таблица 2

Диэлектрическая проницаемость метанола εм и этанола εэ при плоскостях 0,6 и 0,65 г/см3

Продолжение таблицы 2

  Если для этанола все экспериментальные точки расположены на прямой, то для метанола выделяются три точки (рис. 1), которые явно отклонены от линейного вида, не имеют связи ни с температурой, ни с плотностью и их можно отнести к промахам.

  Полученные зависимости аппроксимированы прямой линией. Для рассмотренных веществ ее параметры определены методом наименьших квадратов [8]. Расчеты проведены для у с учетом и без учета составляющей поляризации смещения. В результате получены уравнения, показанные в табл. 1.

  Из табл. 1 видно, что лучшее совпадение с экспериментом (наименьшие среднеарифметические и среднеквадратичные отклонения) для значений у метанола и этанола имеют уравнения (5) и (7), т. е. учет составляющей поляризации смещения приводит к уменьшению отклонения экспериментальных точек от линейности.

Практическое использование полученных уравнений (5) и (7) проиллюстрировано интерполяцией при плотностях 0,6 и 0,65 г/см3 (табл. 2) и экстраполяцией в неисследованные области параметров (табл.3).

 

ВЫВОДЫ

1.  Современная молекулярно-статистическая теория диэлектриков не позволяет описать зависимость ε=f(y) метанола и этанола в широком диапазоне параметров.

2.  Предложенная методика дает возможность проводить расчет диэлектрической проницаемости метанола и этанола закритических плотностей в пределах у до 2,02 с погрешностью, сопоставимой с точностью существующих экспериментальных данных.

 

Таблица 3

Диэлектрическая проницаемость метанола и этанола при плоскостях 0,9 и 1,0 г/см3

ЛИТЕРАТУРА

1.  Ledwig R., Wurflinger А. Dielectric Measurements at High Pressures and Low Temperatures Dielectric and pVT Data of Metanol//Zeitschrift for Physikalische Chemie, 1982,—Bd 132,—S. 21—27.

2.  Timmermans J. The phys-chemist constants of binary Systems in concentrated solvente J. Phys. Chem.—1959.— V. 2.—• P. 1272.

3.  Dannhauser W., В a h e L. W. Dielectric Constant of Hydrogen Bonded liquids. 3. Superbead alcohols//J. Chem. Phys.—1964.— V. 40.— N 10.— P. 3058—3066.

4.  Kirkwood J. G. On the theory of dielectric polarization // J. Chem. Phys.— 1936,—V. 4,—N 3,—P. 592—601.

5.  Смирнов С. H. Расчетно-экспериментальное определение диэлектрической проницаемости сжатых газов, газовых смесей и таблицы для автоматизации контроля теплофизических свойств: Дис. ... докт. техн. наук.— М.: МЭИ, 1986.—600 с.

6.  Р a t е у G. N. An integral equation theory for the dence dipolar hard-sphere fluid // Mol. Phys.—1977,—V. 34,—N 2,— P. 427—440.

7.  Adams D. J. Computer simulation of highly polar liguids: the hard spheres plus point dipole potential // Mol. Phys.'—1980.— V. 40,-—N 5.— P. 1261—1271.

8.  Гусак А. А. Элементы методов вычислений.— Мн.: Изд-во БГУ, 1982.—166 с.

 

Представлена НТС кафедры ВТ и СУ  Поступила 6.05.1991