МОДЕЛИРОВАНИЕ ПРОЦЕССА ДЕЗОДОРАЦИИ ЖИРНЫХ КИСЛОТ
ХЛОПКОВОГО МАСЛА
Хамдамов Анвар Махмудович
канд. техн. наук, доцент, Наманганский инженерно-технологический институт, Республика Узбекистан, г. Наманган E-mail: [email protected]
Сарибаева Дилором Акрамжановна
старший преподаватель, Наманганский инженерно-технологический институт, Республика Узбекистан, г. Наманган E-mail: diloromsaribayeva@mail. ru
MODELING OF PROCESS DEZODORATION FATTY ACIDS COTTONSEED OIL
Anvar Khamdamov
Candidate of techniques sciences, Namangan Institute of engineering and technology, Uzbekistan, Namangan city
Dilorom Saribaeva
Senior lecturer,
Namangan Institute of engineering and technology, Uzbekistan, Namangan city
DOI: 10.32 743/UniTech.2020.80.11-2.94-98.
АННОТАЦИЯ
В статье рассматривается равновесие состояние между жидкой и паровой фазами, то есть коэффициенты 1,2,3,5 и 10 уровневых полиномов для диаграммы X-Y смеси жирных кислот C12 и С18 при давлении 10 кПа. Получено состояние межфазного равновесия, статистически построенных с использованием полиномов.
ABSTRACT
The article presents the equilibrium state between the liquid and vapor phases, that is, the coefficients 1,2,3,5 and 10 level polynomials for the X-Y diagram of a mixture of fatty acids C12 and C18 at a pressure of 10 kPa. The states of interfacial equilibrium are obtained, statistically constructed using polynomials.
Ключевые слова: летучие компоненты, свободные жирные кислоты, хлоповое масло, межфазное равновесие, T-X-Y и Y-X диаграммы, компьютерное моделирование, нейронная сеть, программа MATLAB, 10-уровневых полиномов для состояний равновесия смеси.
Keywords: Volatile components, free fatty acids, cottonseed oil, interfacial equilibrium, T-X-Y and Y-X diagrams, computer modeling, neural network, MATLAB program, 10-level polynomials for mixture equilibrium states.
Введение
Летучие компоненты хлопкового масла, включая свободные жирные кислоты, имеют давление пара при высоких температурах в вакууме. При одинаковых давлениях различные жирные кислоты имеют разную температуру кипения. Иными словами, при одинаковой температуре они имеют разную летучесть в одном и том же вакууме. Исходя из этого, с помощью метода ректификации можно отделить свободные жирные кислоты хлопового масла на фракции.
При производстве хлопкового масла отделяются множество свободных жирных кислот (из 1 тонны
неочищенного масла, кислотное число которого КОН 4,2 мг, отделяется 20-22 кг). Основная часть 95,2% уменьшается в процессе химической рафинации черного масла и 4,8 % в процессе дезодорации (перегонки). Свободные жирные кислоты в процессе ректификации отделяются на фракции, которые являются ценным сырьём в разных областях химической промышленности. В том числе, при производстве пластификаторов, детергентов, синтетических каучуков, прочных полистиролов и химических волокон требуются индивидуальные природные жирные кислоты [1].
Библиографическое описание: Хамдамов А.М., Сарибаева Д.А. Моделирование процесса дезодорации жирных кислот хлопкового масла // Universum: технические науки : электрон. научн. журн. 2020. 11(80). URL: https://7universum. com/ru/tech/archive/item/10981 (дата обращения: 25.11.2020).
№ 11(80)
7universum.com
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
ноябрь, 2020 г.
Состояние межфазного равновесия при разделении жирных кислот можно выразить следующим образом:
х* = / (Г, р, у) (1)
у* = / (Г, р, х) (2)
До настоящего времени состояние равновесия взяты из экспериментальных результатов, а диаграммы, представляющие состояние равновесия Т-Х^ и У-Х, строились с использованием традиционного метода, то есть графоаналитических расчетов процесса. При этом не учитывались изменения фазовых параметров (температуры и концентрации) в ходе процесса.
Для выполнения этого процесса путем современного компьютерного моделирования целесообразно использовать математические выражения или нейронные сети, полученные путем аналитических и экспериментальных расчётов состояния равновесия с высокой точностью [2].
Ниже мы рассмотрим использование математических выражений и нейронных сетей для аналитического определения состояния равновесия и его записи с высокой точностью.
Согласно закону Рауля, когда пар находится в равновесии с идеальной жидкой смесью, парциальное давление компонента в паре равно произведению концентрации этого компонента в жидкости на давление чистого пара:
Pi = xp
(3)
Фугутивность также влияет на природные смеси, которые отличаются от идеальных смесей.
Pi = xiPi7i
(4)
Известно, что существует определенная зависимость между давлением пара и температурой согласно законам термодинамики. Нернст выразил эту взаимосвязь с помощью следующей формулы [3]:
л
ln p =--+ 1,75ln T- — T + C'
RT
R
(5)
где р - давление паров данной жирной кислоты при температуре Т°К; Ло - молекулярная теплота испарения вещества при температуре абсолютного нуля; R - газовая постоянная 1,985 кал/град (1985*4,1868 дж/град); е - коэффициент температуры разности теп-лоемкостей вещества в жидком и паровом состояниях.
Уравнение (5) более удобно, если выразить её следующим образом:
lg P = -
Л
4,571T
+1,75 lg T -
4,571
T + C (6)
Используя значения Л0, е и С, можно рассчитать давление паров жирных кислот в широком диапазоне температур. По крайней мере, аналог позволяет рассчитать температуру кипения жирных кислот при заданном давлении. Для некоторых жирных кислот значения этих констант приведены в таблице 1[3].
Таблица 1.
Коэффициенты для расчета давления паров жирных кислот
Кислоты Число углерода Хо Е С
Каприновая 10 22920 0,0334 11,318
Лауриновая 12 23590 0,0349 11,408
Миристиновая 14 24960 0,0334 11,493
Пальмитиновая 16 26050 0,0326 11,568
Стеариновая 18 27280 0,0317 11,636
Для лауриновой (додекановой) жирной кислоты с 12 атомами углерода в молекуле значения этих констант составляют ^=23590, е=0,0349 и С=11,408. Для пальмитиновой (гексадекановой) жирной кислоты ^0=26050, £=0,0326 ва С=11,568.
Используя выше приведенные формулы и значения, можно рассчитать с достаточной точностью для моделирования значения параметров состояний жидкой и паровой фазовых равновесий компонентов смеси, поступающей в систему. Ниже приведена таблица для диаграммы Т^^ для различных давлений смеси лауриновой и пальмитиновой жирных кислот [4].
Если известны давление и концентрация легколетучих компонентов одной из двух фаз в смеси жирных кислот с помощью формулы (6), можно определить равновесную концентрацию соответствующей второй фазы. В таблице можно вычислить
значения состояний равновесия процессов разделения традиционным способом (используя формулу Нернста, которая представляет собой функциональную зависимость между парциальным давлением и температурой компонента в паровой фазе)
Соотношение давлений частиц дает относительную летучесть:
atj =-
(7)
Значение относительной летучести непостоянно, потому что оно изменяется с изменением температуры. До сих пор в расчетах было получено среднее значение, и расчет выполнялся без изменения этого значения:
—
а-.
= 4а1а2
(8)
Если значение и концентрация вещества в одной из фаз известны, равновесная концентрация в другой фазе может быть рассчитана следующим образом.
У =
аХ
1 -(а+1) X
(9)
Ниже приводится метод расчета равновесных концентраций жирных кислот на границе раздела фаз с числом атомов углерода в молекуле от Сю до С18 с использованием компьютера.
Следующая программа МАТЬАБ была разработана для расчета давления пара Р, соответствующего температуре Т, заданной уравнением 6.
Результаты можно наблюдать на графике, показанном на рисунке 1.
Используя структурированные программы МАТЬАБ, можно быстро, точно и легко выполнять
вычисления для соединений жирных кислот, получать графики и создавать динамические модели Simulink[5].
В некоторых случаях при моделировании процесса разделения необходимо проводить быстрые численные эксперименты с заданной точностью. В таких случаях взаимосвязь параметров равновесного состояния, полученных экспериментально или аналитически, целесообразно выражать с помощью п-уровневых уравнений, т.е. полиномов. Общий вид полинома степени п выглядит следующим образом:
у = ах" + а2х" 1+ ■■■ + а„_1х2 + апх + ап+1 (10)
Чтобы выразить состояние равновесия хлопкового масла в процессе разделения жирных кислот с помощью полиномов, ниже представлен метод расчета их коэффициентов с помощью компьютерной программы МАТЬАБ.
220
215
210
206 200 195
130
105
180
175
1 ) 4 1
-------г------* Г-----* 4 1 1 !
1 | |
! !
Ты
2 1 1 1 ^
1 1 1 ^^ 1 1 т
1
_1_
X
-1_
о
0.1 0.2 0 3 0.4 0.5 0 6 0.7 О.В 0.9
Рисунок 1. ^X^диаграмма смеси лауриновой и пальмитиновой жирных кислот при давлении 1,87 кПа (1-состав легкокипящего компонента в газовом фазе, 2-состав легкокипящего компонента в жидком фазе)
Коэффициенты 10-уровневых полиномов для состояний равновесия смеси лаурин-пальмитиновых
жирных кислот при давлении 1 кПа равны на следующие:
х* = / ( Р, У) ^1=351,6583252453378 а2=-165,1551092502495 аэ=3321,188566628895 Я4=-3723,533887560301 аз=2541,963170053172 а6=-1083,143380106641 ^7=283,9568381736603 а8=-43,12702533994424 ^9=3,501584063096183 аю=0,086220980919603 яц=0,000061497692188
I = / (Р, У ) а1=-2224,066396153840 а2= 10367,12782605067 аэ=-20726,68830503027 а4= 23121,89133556114 аз=-15718,72028978193 а6= 6671,138777123714 а7=-1745,373058345953 а8= 260,8025720254839 а9=-27,51090147217575 а10=-17,3829459094002 ап= 202,4286310317999
• 7universum.com
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ_ноябрь. 2020 г.
Использование полиномов для выражения состояний равновесия ускоряет вычисление компьютерной модели и позволяет проводить большее количество экспериментов за более короткое время.
Определено равновесие между жидкой и паровой фазами, то есть коэффициенты 1,2,3,5 и 10 уровневых полиномов для диаграммы X-Y смеси жирных кислот C12 и С18 при давлении 10 кПа. Они равны следующим показателям:
Таблица 1.
Показатели
Степень Коэфициенты Степень Коэфициенты
1 01=0.73407625469346 02=0.41718273155660 10 a1=-0.07437327333851 02=0.46561549878869 a3=-1.30813494173110 04=2.19074120624204 as=-2.45573147018139 Й6=1.96715181946477 Й7=-1. 17749114452463 08=0.54423929766038 Й9=-0.19790086625601 йю=0.05588384347727 a11=0.00000000301913
2 Й1=-1.39073 173195892 02=2.12480798665238 аэ=0.18771199578338
3 a1=2.32176077906796 Й2=-4.87337290056086 ЙЭ=3.51094560697153 04=0.07508338039080
5 a1=5.08552373008970 Й2=-16.27559179604100 Й3=20.71908424302194 Й4=-13.63462559516901 й5=5.10140732924744 a6=0.00896585166142
Рисунок 2. График состояния межфазного равновесия, статистически построенные
с использованием полиномов
Из этого графика можно увидеть, что, когда равновесие моделируется на кривой концентрации с использованием многочленов п-го уровня, то есть когда математическая функциональная связь определяется с помощью компьютера, точность фактического процесса увеличивается по мере увеличения
уровня уравнения. И по мере того, как обратные уравнения становятся проще, то есть с уменьшением уровня, уменьшается и ошибка. Когда в простейшем уравнении используется линейная функция первого порядка, ошибка становится очень большой.
Рисунок 3. График отличий полиномов от реальных процессов, представляющих состояния межфазного равновесия
Заключение
Графики ошибок вышеупомянутых функций в выражении фактической функции процесса приведены на графике выше. Видно, что погрешности результатов высокоуровневых функций близки к нулю. Из протестированных уравнений первой, второй, третьей, пятой и десятой степени ошибка уравнения 10-й степени была наименьшей, но точность уравнения пятой степени представляется достаточной для рассматриваемого нами процесса.
Таким образом, преимущества предложенного способа обеспечивают создание имитационной модели программы МАТЬАБ со скоростью и высокой точностью получения результатов фазового равновесия концентраций смеси жирных кислот, полученных экспериментальным путем.
Список литературы:
1. Артиков А., Хамдамов А., Хамидов Б., Маматкулов О., Маматкулов А. Математическая модел динамики распылительной перегонки инертным газом // Математические методы в технике и технологиях: Сб. тр. XVIII межд. науч. конф. - Казань, 2005. - С. 222-223.
2. Артиков А.А., Худайбердиев А.А., Хамдамов А.М. Математическое моделирование равновесного состояния экстракционного бензина и жирных кислоть // 20 международная научная конференция «Математические методы в технике и технологиях» (ММТТ-20), Ярославль, 28-31 мая, 2007: Сборник трудов. Т. 5. Секц. 11. -Ярославль: ЯГТУ, 2007. - с. 240-241.
3. Тютюнников Б.Н. Химия жиров. - М.: Пищевая промышленность, 1966. - 632 с.
4. Файнберг Е.Е. Ректификация природных жирных кислот и высших жирных спиртов. - М.: Пищевая промышленность, 1970. - 183 с.
5. The Mathworks Inc. Simulink - Simulation and Model-Based Design. http://www.mathworks.com/products/simulink