CC BY
35
- ПРОЦЕССЫ, ОБОРУДОВАНИЕ И АППАРАТЫ
ПИЩЕВЫХ ПРОИЗВОДСТВ -
УДК 544.015; 517.958; 544.275.5; 544.275
З.Н. Есина, А.М. Мирошников, М.Р. Корчуганова
ПРОГНОЗИРОВАНИЕ ФАЗОВЫХ РАВНОВЕСИЙ В СИСТЕМАХ ГЛИКОЛЕЙ И ЭФИРОВ ГЛИКОЛЕЙ
Приводятся аналитические выражения для энтальпии испарения и энтальпии плавления гликолей и эфиров гликолей, полученные с помощью теории подобия. Показано, что прогнозирование этих параметров позволяет рассчитать фазовые диаграммы «жидкость - твердое» и «жидкость - пар», а также эвтектические и азео-тропные параметры в системах «эфир гликоля - вода» и «эфир гликоля - алкан».
_______Энтальпия испарения, энтальпия плавления, термодинамическое подобие, эвтектика, азеотроп.________
Введение
Энтальпия испарения и плавления для некоторых членов гомологических рядов неизвестна или измерена недостаточно точно. В работе [1] отмечается, что расчет параметров эвтектических и азеотропных смесей, используемых в технологических процессах кристаллизации и ректификации, осложняется отсутствием точных данных об энтальпии фазового перехода чистых компонентов. Особенно следует отметить эту проблему при расчете термодинамических параметров растворов насыщенных жирных кислот, а также эфиров гликолей, имеющих важное значение при разработке составов медицинских и косметических препаратов, тормозных жидкостей, антифризов, теплоаккумуляторов.
Температура кипения Т кип и температура плавления Тпл, давление Р, плотность р, энтальпия испарения Нисп и теплота плавления Нпл имеют определяющее значение для расчета термодинамических свойств чистых компонентов и растворов. Прогнозирование свойств соединения предполагает их определение на основе структурной формулы. Теория, связывающая структуру молекулы известного состава с макроскопическими параметрами, в настоящее время не завершена, поэтому необходимо обобщить эмпирические данные о свойствах соединений. Теоретической основой обобщения является метод термодинамического подобия, представляющий собой часть общей теории подобия [2].
В качестве основного математического метода изучения фазового равновесия в реальных растворах взят метод термодинамических потенциалов.
Объект и методы исследования
В теорию подобия в зависимости от природы изучаемого объекта вводится система единиц измерения, и затем делается переход к безразмерным, так называемым приведенным величинам. Для изменения масштаба, т.е. перехода от одной системы к другой, необходимо совершить линейное преобразование соответственных величин:
х1 = сх х1,
где сх - постоянная, называемая множителем подо-
х1 '
бия; х и х' - соответственные величины.
Закономерности, выявленные для одной системы, справедливы для группы подобных систем, что позволяет получить информацию о свойствах неизвестных систем.
Метод термодинамического подобия кратко можно охарактеризовать как способ выделения определенного числа параметров, однозначно определяющих термофизические свойства веществ. В качестве параметров приведения применяются энтальпия плавления Н пл и испарения Н исп, относительная молекулярная масса М, температура плавления Т пл и
кип
кипения Т .
Теплоту фазового перехода жидкость - пар можно найти по зависимости давления от температуры на соответствующей кривой двухфазного равновесия, которая выражается уравнением Клаузиуса -Клапейрона:
1п Р = -Нисп / КТ + свтг. (1)
Уравнение (1) получено при условии, что энтальпия испарения не зависит от температуры, пар является идеальным газом и объем жидкости мал по сравнению с объемом газа. При отсутствии данных о давлении пара энтальпия испарения может быть вычислена по эмпирическим уравнениям зависимости от температуры кипения или от критических параметров. Группа веществ, имеющих близкое химическое строение, характеризуется параметрами подобия, которые варьируются при переходе к другой группе. Соотношения подобия связаны с понятием критического состояния. В [3] приводится эмпирическая формула
4 [(-Тр )(кр-Тр ) + гкр
Т - Т
кр
Т - Т
кр тр
где Нтр - энтальпия испарения в тройной точке; Ткр -критическая температура; Ттр - температура тройной точки; гкр - определяющее число подобия.
В работе [3] утверждается, что широкий класс жидкостей, за исключением квантовых жидкостей, имеет универсальное поведение на кривой испарения.
В данной работе молярная энтальпия плавления чистого компонента для гомологических рядов: алка-нов, спиртов, кислот, гликолей и эфиров гликолей -представлена как функция числа углеродных атомов N, относительной молекулярной массы М, темпера-
кип
туры кипения Т и универсальной газовой постоянной К:
Н„сп = а^)М кт — + ), (2)
Ж
где а(ЩМ / NR - определяющий параметр подобия, безразмерный коэффициент; В( N) = р( N )(-1)N +у( N) -функция, введенная для более точного моделирования энтальпии испарения членов гомологического ряда с четным и нечетным числом атомов углерода; Р^), у(Ю имеют размерность Дж/моль.
В табл. 1 представлен коэффициент уравнения
(2): а(Ы) как функции от числа атомов углерода в молекуле, для гликолей и эфиров гликолей Р( N) = у( N) = 0.
Таблица 1
Параметры уравнения для расчета энтальпии испарения членов гомологических рядов гликолей и эфиров гликолей
Гомологический ряд а( N), Дж/моль К
Гликоли N (N + 1,1) 1,7(N -1)
Эфиры гликолей N (N +1,1) 1,15(N -1)
На рис. 1 приведены экспериментальные данные и результаты прогнозирования энтальпии испарения для гомологического ряда этиловых эфиров гликолей.
иосюо ■■ хооооо ■
л
О 90000 2
80000
С
" 70000
X
60000 50000
4 6 8
N
Рис. 1. Энтальпия испарения гомологического ряда
этиловых эфиров гликолей: ряд 1 - нисп = -(N+11^MT ,
мод 1,15(N -1)
Л гг исп ряд 2 - Нэксп
Для подобных членов различных гомологических рядов в данной работе предлагается отношение:
НГ - Вп (N) ая (N)МяТякип
НГ - Вп (N) а, (N)МЯТГ ’
где В1 (N) = р,. (Ж) (-1)+ у,. (N, , = я, п.
Для связи энтальпии испарения различных членов одного гомологического ряда с числом углеродных атомов N1, N2 можно использовать соотношение:
Н™ -В(N2) а(N2)М(N^
Н™ - В(N1) а(N. )М(N. )Т™ ,
где а(N2)/ а(N1) - отношение определяющих параметров подобия.
Энтальпия плавления чистого компонента может быть представлена для молекулы с числом углеродных атомов N как функция молярной массы М, температуры плавления Т пл и универсальной газовой постоянной К:
Н1Ш = «(^Ж КТ 1Ш + А( N), (3)
ЫК
где A(Ы) = Ь(И)(-1)* + с^). Здесь а(И)М / Ж -определяющий параметр подобия, безразмерный коэффициент, Ь^), сN имеют размерность Дж/моль.
В табл. 2 приведены коэффициенты уравнения
(3): а(Щ, b(Ы), c(Ы) как функции от числа атомов углерода в молекуле. Это представление не является единственно возможным и зависит от точности экспериментальных данных, используемых для аппроксимации энтальпии.
Таблица 2
Параметры уравнения для расчета энтальпии плавления членов гомологических рядов гликолей и эфиров гликолей
На рис. 2 приведены экспериментальные данные и результаты прогнозирования энтальпии плавления для гомологического ряда гликолей.
Рассматривается подобие энтальпии плавления, взятой с точностью до постоянной A(N), соответствующих членов гомологических рядов гликолей и эфиров гликолей:
НПГ - Ап (N) = «я (N)МЯТЯШ НЯ - Ап (N) а,п (N)МпТТ’
Гомологический ряд а(М), Дж/моль К ът Дж/моль е(Ю, Дж/моль
Гликоли и эфиры гликолей (1 < N < 10) N2 +1 26,5 0 -320 N+10755
где А1 (N = (N (-1)Ы + с1 (N, , - номер ряда. Определяющим параметром подобия является безразмерный комплекс: a(Ы)M/.
О
s
*
d
Ряді
Ряд2
Рис. 2. Энтальпия плавления гликолей: ряд 1 -
H” = N +1MT - 320N +10755 = ряд 2 - Hl
мод 26,5N
пл
эксп
Для подобия энтальпии плавления различных членов одного гомологического ряда с числом углеродных атомов N1, N справедливо соотношение:
Н!? - А(N2) а(N2)М(N2 )TлПл
Н” - А(N2) а(N1)М(N1)TЛПI
Выражение (4) позволяет рассчитать энтальпию плавления по известной энтальпии испарения
H- b(N)(-1)N - c(N) a(N) Tn
Hисп -P(N) -y(N) a(N) Tк
(4)
В табл. 3 приведено отношение параметров подобия a(N)/a(N) для такого сравнения.
Таблица З
Отношение параметров подобия для членов гомологических рядов гликолей и эфиров гликолей
Гомологический a( N )
ряд a( N)
Гликоли (N2 +1)(N -1)
(1 < N < 10) 15,59 N (N +1,1)
Эфиры гликолей ( N2 +1)( N +1)
(1 < N < 10) 2З,04 N ( N +1,1)
Результаты и их обсуждение
PCEAS (Phase Chart Eutectic and Azeotropic Systems) [4] является моделью уравнения состояния и основана на минимизации избыточной энергии Гиббса по параметру сольватации X, характеризующему отношение числа молекул A к числу молекул B в молекулярном соединении.
Результатом работы является расчет фазовых равновесий «жидкость - твердое» и «жидкость -пар» в бинарных системах, содержащих эфиры гли-колей с применением программы PCEAS [5]. Энтальпия плавления и испарения для эфиров гликолей и алканов, необходимые для расчета, взяты по моделям, предложенным в данной работе. Фазовые равновесия c учетом ассоциации и сольватации молекул изучались при нормальном атмосферном давлении.
На рис. З-4 приведены кривые ликвидуса и соли-дуса в системах метилцеллозольв - вода и метилцел-лозольв - циклогексан.
х1, мае, д.
Рис. 3. Фазовая диаграмма «жидкость - твердое» системы «метилцеллозольв - вода»
Рис. 4. Фазовая диаграмма «жидкость - твердое» системы «метилцеллозольв - циклогексан»
В табл. 4 приводится результат расчета эвтектических и азеотропных параметров в системах на основе эфиров гликолей.
N
Таблица 4
Расчетные и экспериментальные эвтектические и азеотропные параметры в системах, содержащих эфиры гликолей
Раствор Хэвт, масс., расч. t , oC, эвт расч. Чаз , масс., расч. ,C ч. o сч , а аз р t q , масс., ^аз ’ ’ эксп. [б] *аз ^ эксп. [б]
Метилцеллозольв - вода 0,959 -102,87 0,21 99,77 0,212 99,75
Метилцеллозольв - циклогексан 0,452 -111,бб 0,09 77,5б 0,080 79,80
Метилцеллозольв - гептан 0,404 -102,40 0,18 91,88 0,230 92,50
Метилцеллозольв - октан 0,809 -88,73 0,329 105,59 0,480 110,00
Метилкарбитол - вода 0,948 -91,09 не азеотр. не азеотр. не азеотр. не азеотр.
Метилкарбитол - ундекан 0,959 -85,б5 0,340 173,59 0,400 (об.) 178,70
Этилцеллозольв - вода 0,909 -80,57 0,297 99,74 0,300 99,30
Этилцеллозольв - гептан 0,291 -97,74 0,147 93,б9 0,140 94,50
Этилцеллозольв - октан 0,801 -78,54 0,294 110,11 0,380 11б,00
Этилцеллозольв - нонан 0,б0б -79,91 0,4б2 120,85 0,50 128,00
Результаты расчета, приведенные в табл. 4, согласуются с экспериментально полученными параметрами критических точек [6].
Выводы и рекомендации
Прогнозирование термодинамических параметров гликолей и эфиров гликолей на основе теории подобия позволяет рассчитать фазовые диаграммы бинарных систем, для которых отсутствуют экспериментальные данные об энтальпии плавления и энтальпии испарения чистых компонентов. Метод может быть рекомендован при выборе компонентов и состава с целью получения требуемых свойств раствора: температуры и энтальпии плавления, температуры и энтальпии испарения, что позволит значительно сократить время и средства на проведение экспериментов.
Список литературы
1. Kenisarin, M. Solar energy storage using phase change materials / M. Kenisarin, K. Mahkamov // Renewable and Sustainable Energy Reviews. - 2007. - № 11. - P. 1913-1965.
2. Филиппов, Л.П. Прогнозирование теплофизических свойств жидкостей и газов / Л.П. Филиппов. - М.: Энерго-атомиздат, 1988. - 168 с.
3. Meyra, A.G. Universal behavior of the enthalpy of vaporization: an empirical equation / A.G. Meyra, V.A. Kuz, G.J. Zarragoicoechea // Fluid Phase Equilibria. - 2004. - № 218. - P. 205-207.
4. Есина, З.Н. Математическое моделирование фазового перехода жидкость - твердое / З.Н. Есина, М.Р. Корчугано-ва, В.В. Мурашкин // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. - 2011. - № 3 (16). - С. 13-23.
5. А. с. 2012618394 Российская Федерация. Phase Chart Eutectic and Azeotropic System (PCEAS) / З.Н. Есина, В.В. Мурашкин, М.Р. Корчуганова. - № 2012616324; заявл. 25.07.2012; рег. 17.09.2012.
6. Огородников, С.К. Азеотропные смеси: Справочник / С.К. Огородников, Т.М. Лестева, В.Б. Коган; под ред. проф. В.Б. Когана. - Л.: Химия, 1971. - 849 с.
ФГБОУ ВПО «Кемеровский технологический институт
пищевой промышленности», 650056, Россия, г. Кемерово, б-р Строителей, 47.
Тел/факс: (3842) 73-40-40, е-шаП: office@kemtipp.ru
На рис. 5 приведена зависимость температуры кипения от массовой доли первого компонента в системе «мелилцеллозольв - вода».
390
386
382 380 378 376 374 37? J
0.1 0.2 0.3 0.4 0.5 0.6 0.7 08 0.9 1
х1, мае. д.
Рис. 5. Фазовая диаграмма «жидкость - пар» системы «метилцеллозольв - вода»
SUMMARY
Z.N. Esina, A.M. Miroshnikov, M.R Korchuganova PREDICTION PHASE EQUILIBRIA IN THE GLYCOLS AND GLYCOL ETHERS
The paper presents analytical expressions for the enthalpy of vaporization and the enthalpy of fusion of glycols and glycol ethers obtained by using the theory of similarity. It is shown that the prediction of these parameters allows us to calculate the phase diagrams of liquid-solid and liquid-vapor, and eutectic and azeotrope parameters in systems of glycol ether-water and glycol ether-alkane.
Enthalpy of vaporization, enthalpy of melting, thermodynamic similarity, eutectic, azeotrope.
Kemerovo Institute of Food Science and Technology 47, Boulevard Stroiteley, Kemerovo, 650056, Russia.
Phone/fax: +7 (3842) 73-40-40, е-mail: office@kemtipp.ru
Дата поступления: 08.07.2013
- -
