CC BY
10
№ 12 (69)
UNIVERSUM:
технические науки
декабрь, 2019 г.
КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ РАВНОВЕСНОГО СОСТОЯНИЯ ПРОЦЕССА ТЕПЛОМАССООБМЕНА ПО ДАВЛЕНИЮ В ДВУХФАЗНОЙ СИСТЕМЕ
Машарипова Зулхумар Атабековна
PhD по техническим наукам, доцент кафедры "Информатика, автоматизация и управления "
Ташкентского химико-технологического института,
Узбекистан, г. Ташкент E-mail: zulxumor- 72@mail. ru
Сарболаев Фаррухбек Набиевич
ст. преп. кафедры «Информатика, автоматизация и управления» Ташкентского химико-технологического института,
Узбекистан, г. Ташкент E-mail: sarbo layev_f@mail. ru
Исламова Фарида Камилджановна
ст. преп. кафедры «Информатика, автоматизация и управления» Ташкентского химико-технологического института,
Узбекистан, г. Ташкент
COMPUTER SIMULATION OF THE EQUILIBRIUM STATE OF THE HEAT AND MASS TRANSFER PROCESS BY PRESSURE IN A TWO-PHASE SYSTEM
Zulkhumar Masharipova
PhD in Technical Sciences, associate professor "Informatics, Automation and Control" department of the Tashkent chemical-technological institute,
Uzbekistan, Tashkent
Farrukhbek Sarbolayev
senior lecturer of department "Informatics, automation and control" of the Tashkent chemical-technological institute,
Uzbekistan, Tashkent
Farida Islamova
senior lecturer of department "Informatics, automation and control" of the Tashkent chemical-technological institute,
Uzbekistan, Tashkent
АННОТАЦИЯ
Результаты расчета на компьютерной модели показывают, что в начальном этапе сушки продукта при атмосферном давлении равновесное давление системы имеют большое значение, около 200 кПа. По мере высушивания материала равновесное давление начинается уменьшаться, и динамический процесс завершается уравновешиванием равновесного давления к общему давлению системы (100 кПа). Соответственно уменьшается влага материала от 60 % до 41%. Таким образом, для расчета процесса сушки и других тепломассаобменных процессов можно за основу принимать равновесного давления.
ABSTRACT
Results based on a computer model shows that in the initial stage of drying of the product at atmospheric pressure, the equilibrium pressure of the system are important, about 200 kPa. As the drying of the equilibrium, pressure of the material begins to decrease, and the dynamic process is completed balancing the equilibrium pressure to the total pressure of the system (100 kPa). Accordingly, the moisture of the material is reduced from 60% to 41%. Thus, for the calculation of the drying process, and other processes heat mass transfer can take as the basis of the equilibrium pressure.
Ключевые слова: сушка, температура, равновесия, давления, летучий компонент, двухфазная система, тепломассообмен, движущая сила, математическая модель.
Keywords: drying, temperature, equilibrium, pressure, volatile component, two-phase system, heat and mass transfer, driving force, mathematical model.
Библиографическое описание: Машарипова З.А., Сарболаев Ф.Н., Исламова Ф.К. Компьютерное моделирование равновесного состояния процесса тепломассообмена по давлению в двухфазной системе // Universum: Технические науки : электрон. научн. журн. 2019. № 12(69). URL: http://7universum.com/ru/tech/archive/item/8509
№ 12 (69)
Д Ul\li
/ш. те:
UNIVERSUM:
технические науки
Из-за многих факторов анализ состояния равновесия процесса тепломассопереноса требует логического мышления и абстрактного подхода. Существуют различные методы расчета равновесия процесса тепломассопереноса в двухфазной системе. Они основаны на графоаналитических методах или методах последующего приближения с использованием табличных данных и номограмм [1]. С развитием информационных технологий разрабатываются методы автоматизированного расчета технологических процессов. Компьютерное отображение равновесия процесса тепломассопереноса в двухфазной системе отличается своей сложностью. Для удобного изложения мыслей мы проанализировали процессы дистилляции инертного газа, сушки материалов и т. Д. Используя метод многостадийного системного анализа [2], двухфазную систему анализируют путем
декабрь, 2019 г.
определения входных и выходных параметров. Элементами двухфазной системы являются жидкая и газовая фазы. Представляя их аналогично системам, имеющим входные и выходные параметры, определяется их взаимное влияние.
Одним из основных факторов в процессе тепло-и массообмена является давление в двухфазной системе. По мере увеличения давления в двухфазной системе сорбционные свойства будут возрастать. Это связано с парциальным давлением летучей составляющей и ее равновесным значением. Таким образом, можно приступать к определению парциального давления летучих компонентов. Общий вид дисплея компьютера (на основе использования пакета приложений МЛТЬЛБ) двухфазной системы показан на рис. 1.
Рисунок 1. Отображение двухфазной системы для изотермического случая
Для краткости значения входных параметров указываются одновременно на рисунке с обозначениями входных параметров. Анализ этих взаимных влияний включает в себя главным образом четыре фактора. В связи с этим движущая сила процесса [3] может быть охарактеризована сочетанием четырех факторов.
F(t,x,y,P) = 0 ,
(1)
это:
/ - температура в системе (в основном температура жидкой фазы);
х - концентрация летучего компонента в жидкой фазе;
у - концентрация летучего компонента в газовой фазе;
Р - давление в системе.
При анализе процесса в открытой системе на состояние равновесия влияют и другие факторы, в частности, перепад фазовой температуры и т.д.
В зависимости от величины равновесных значений может произойти десорбция от жидкой к газовой фазе или поглощение от газовой фазы летучей составляющей.
В этой статье мы ограничимся отображением процесса перехода системы в состояние равновесия с использованием фактора давления из четырех основных факторов, перечисленных выше. Используя комбинацию из четырех факторов, можно будет создать математическую и впоследствии компьютерную модель процесса перехода в равновесное состояние теп-ломассопереноса в двухфазной системе.
Мы представляем отображения последовательно по одному фактору давления в системе из четырех факторов.
Одним из факторов газовой фазы является концентрация летучего компонента, которая в большинстве случаев характеризуется (как известно) отношением числа молекул летучего компонента к общему количеству молекул газа фаза:
№ 12 (69)
А1
ЛОА
UNIVERSUM:
технические науки
декабрь, 2019 г.
y
N
N.
(2)
Робщ -—
У
(8)
и для тепломассообменного процесса в двухфазной системе, например, перегонка инертным газом или сушка материалов, концентрация летучего компонента в газовой фазе определяется из общеизвестно уравнения:
(
Gy + G M
\
Gy + Gy Gw - Gc
M
M.
(3)
где Оу- расход летучего компонента из жидкости;
О у- общий расход летучего компонента в газовой фазе;
Му- молекулярная масса летучего компонента;
Му- молекулярная масса воздуха;
Ооу- расход летучего компонента в первичном воздухе;
О о- общий расход высушиваемого материала
Здесь, расход испаренного летучего компонента записывается в виде:
G = g, 1 - ^
(4)
Тогда, общий расход летучего компонента в газовой фазе:
После ряда математических преобразований уравнение для фактического парциального давления летучей компоненты в газовой фазе получается в виде уравнения:
Gy + Gxy Gw - Goy
Mv
Mw
общ
Gy + Go Mw
(9)
Где : Оу- расход летучего компонента из жидкости;
О у- расход влаги в воздуха;
Му- молекулярная масса летучего компонента;
Му- молекулярная масса воздуха;
Ооу- расход летучего компонента в первичном воздухе;
Ру парциальное давление летучего компонента;
Робщ- общее давление
С другой стороны, для газовой фазы равновесная концентрация летучего компонента, участвующего в процессе массообмена, характеризуется состоянием жидкой фазы (ее температура, концентрация летучего компонента) и давления в системе.
Определение зависимости величины равновесного давления позволяет определить движущую силу в виде разницы между фактическим и равновесным давлением.
Лр= p*_p
(10)
G w — G + Gy
(5)
Равновесное парциальное давление летучего компонента определяется по показателям жидкой фазы. Его можно характеризовать концентрацией и температурой жидкой фазы.
P*—f(t,x)
(6)
Например, для системы аммофосной пульпы -газа из преобразования ее точки кипения [4] можно написать уравнение для равновесного давления летучего компонента в газовой фазе в виде:
Р*у = (0.571 -1.517 х -0.9635)26455/у (7) г ^ 50 - 100
Как известно c другой стороны, парциальное давление летучего компонента зависит от его молярной концентрации в газовой фазе и общего давления системы
Ру =у* Робщ (27)
Или, общее равновесное давления записывается в виде:
Переход в состояние равновесия можно охарактеризовать следующим уравнением массообмена:
dx dr
ßP ( p
p)
(11)
где: рр - коэффициент массообмена по давлению в системе
Если равновесное давление превышает фактическое давление, происходит процесс десорбции, в противном случае происходит процесс сорбции.
С учетом рассмотренных выше уравнений процесс перехода к системному равновесию по давлению может быть отображен следующей системой математических выражений:
Gy+ Gxy Gw - Goy
P *=P*y
ОБЩ J
t ^ t ж P ^ P
p ~ p ОБЩ
X ^ X ж dx
Mv
Mw
Gy+ Go Mw
d т
= ßv(P* - P)
№ 12 (69)
UNIVERSUM:
технические науки
декабрь, 2019 г.
Компьютерное отображение процесса перехода системы в равновесное состояние аналогично показанному на рис. 2.
Рисунок 2. Отображение компьютерного алгоритма расчета процесса по равновесному давлению в
системе
На компьютере проводится процесс изотермической десорбции при переходе в равновесное состояние под давлением в системе. Здесь для анализа компьютерный расчет выполняется по следующим значениям входных параметров: Gg =9 kg/c - расход газа, Go = 11 kg/c - расход пульпы, t=65 0 С - температура пульпы, х0=0.6 - концентрация летучего компонента в жидкой фазе, [ =0.5 еЛ -коэффициент мас-
сообмена по концентрации летучего компонента в газовой фазе, ^0=100 кПа - давление в системе. Компьютер автоматически рассчитывает следующие выходные параметры: у - концентрация летучего компо-
*
нента, у - равновесная концентрация летучего
компонента, Pp- давление летучего компонента в газовой фазе, Gw - расход летучего компонента, х - концентрация летучего компонента в жидкости, Gi - расход выходящей жидкости.
На рис. 3 показан характер перехода изобарической системы тепло-массообмена в равновесное состояние при равновесном давлении. Как видно из рисунка, изначально в неравновесной системе равновесное давление имеет большое значение (в данном случае около 200 кПа), с переходом в равновесное состояние оно становится близким к фактическому значению давления (в данном случае 100 кПа). Соответственно, уменьшается влажность материала и влажность газовой фазы. В общем, процесс становится статическим состоянием равновесия.
№ 12 (69)
UNIVERSUM:
технические науки
декабрь, 2019 г.
Рисунок 3. Характеристика динамики перехода системы по равновесному давлению. В правом верхнем углу показаны обозначения параметров ординат и соответствующие кривые
ЗАКЛЮЧЕНИЕ
Полученные результаты позволяют рекомендовать новый метод компьютерного расчета процесса тепломассообмена. Принимая за основной блок разработанного компьютерного отображения, можно формализовать модели и методы автоматизированного расчета различных процессов. Использование предложенных методов компьютерного отображения равновесия тепломассообменного процесса в двухфазной системе способствует анализу и автоматизированному расчету равновесия в двухфазной системе, принимая за основу равновесного давления.
ОБОЗНАЧЕНИЯ
Gy- расход летучего компонента из жидкости, kg/c;
Gy- расход воздуха, kg/c;
G0 -расход пульпы, kg/c;
G g-расход газа, kg/c;
Goy- расход летучего компонента в первичном воздухе, kg/c;
Gw - расход летучего компонента, kg/c;
Gi - расход выходящей жидкости, kg/c;
Mw- молекулярная масса летучего компонента,
кг;
Mv- молекулярная масса воздуха, кг;
N-молекулы летучего компонента;
Иобщ-обшие количества молекулы летучего компонента;
P - давление летучего компонента в газовой фазе, kPa;
P - действительное парциальное давление, kPa;
P*- равновесное парциальное давление, кПа;
Ру парциальное давление летучего компонента, кПа;
Робщ- общее давление, кПа;
Р*общ - равновесное общее давление
Р р - давление летучего компонента в газовой фазе, кПа;
t - температура в системе, о С;
^ действительная температура в жидкой фазе, о С; равновесная температура в жидкой фазе, о С;
tа.д- температурной депрессии при давлении Р = 100 кПа , о С;
t в- температура кипения воды, о С;
*
х - равновесная концентрация летучего компонента в жидкости, % ;
х - концентрация летучего компонента в жидкой фазе, % ;
у - концентрация летучего компонента в газовой фазе, % ;
у - действительная концентрация летучего компонента, % ;
*
у - равновесная концентрация летучего компонента, %;
[у - коэффициент массообмена по концентрации
летучего компонента газовой фазе, с"1;
3Р - коэффициент массообмена по давлению в системе, с"1;
3 - коэффициент массообмена по температуре в системе, с"1;
[х коэффициент массообмена по концентрации в жидкости, с"1.
^исок литературы:
1. Павлов К.Ф, Романков П.Г, Носков А.А., Примеры и задачи по курсу процессов и аппаратов химической технологии, Госхимиздат, 1987.- 576 с.
№ 12 (69)
UNIVERSUM:
технические науки
декабрь, 2019 г.
2. Артыков А. Multi-step method of computer model formalization with fuzzy sets application (Многоступенчатый метод формализации компьютерной модели с применением элементов нечетких множеств). WCIS-2004, world conference on intelligent systems for industrial automation, Tashkent-2004, TSTU
3. Юсуфбеков Н.Р, Артыков А., Юнусов И.И, Механизм перегонки с водяным паром и математическое моделирование процесса. Сб. Тепло и массоперенос. Том 4. Минск. 1972. С.325-327
4. Рейпназарова З.Д., Артиков А. Температура кипения аммофосной пульпы из фосфоритов Центральных Кызылкумов // Узб. Хим. ж.-Ташкент, 2007.-№4.-С. 34-37
