Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,074

МОДЕЛИРОВАНИЕ ПРОЦЕССА ПОЛУЧЕНИЯ ДВУХФАЗНОЙ СИСТЕМЫ ЖИДКОСТЬ – ТВЕРДОЕ В АППАРАТЕ РОТОРНОГО ТИПА

Ширина Н.Ю. 1 Сидоров В.Н. 1 Мурашов А.А. 1
1 ГОУ ВПО «Ярославский государственный технический университет»
Для получения двухфазной системы жидкость – твердое предполагается использование аппарата роторного типа. В процессе диспергирования происходят два совмещенных процесса: перемешивание системы жидкость – твердое и дезагрегация твердых частиц. Для моделирования процесса смешения предлагается использовать диффузионную модель. Расчет производится исходя из предположения, что процесс перемещения твердых частиц в основном, происходит в направлении, перпендикулярном движению жидкости. Поток твердых частиц пропорционален площади циркуляционного контура и градиенту концентрации. В соответствии с принятой моделью получено уравнение диффузии твердых частиц в циркуляционном потоке жидкости. Коэффициент диффузии зависит от размера твердых частиц, числа оборотов ротора и свойств жидкости. Граничные условия для уравнения диффузии строятся из тех предположений, что на стенке аппарата поток концентрации твердых частиц отсутствует. Начальное условие строится из тех предположений, что в начале процесса твердые частицы подаются сверху на движущийся слой жидкости и, следовательно, полностью заполняют все циркуляционные контуры, граничащие со стенкой аппарата. Для решения уравнения диффузии используется разностная схема с последующим использованием сеточного метода. При моделировании процесса дезагрегации используется экспоненциальная зависимость среднего размера твердых частиц от времени. Параметры модели определяются методами статистического анализа на основе экспериментальных данных в программе Statistica 6.0. Качество процесса перемешивания жидкости и твердых частиц оценивается коэффициентом неоднородности распределения твердых частиц по объему жидкости. Модель использована для расчета получения двухфазной системы вода – кальцит с процентным содержанием кальцита 40 %. Представлены результаты регрессионного анализа процесса дезагрегации частиц кальцита; сопоставление теоретических и экспериментальных данных по качеству смеси.
двухфазная система
жидкость
твердые частицы
диспергирование
перемешивание
дезагрегация
аппарат
роторный тип
диффузионная модель
циркуляционный контур
концентрация
средний размер
разностная схема
коэффициент неоднородности
1. Кафаров В.В., Меньшиков В.В. Применение рециркуляции для повышения эффективности химико-технологических процессов // ТОХТ. – 1995. – Т.24, № 6. – С. 565–574.
2. Кухтин И.В. Усовершенствование лакокрасочного оборудования для производства ЛКМ // Лакокрасочные материалы и их применение. – 2009. – № 1–20. – С. 70.
3. Соловьев В.П. Современное диспергирующее оборудование для производства лакокрасочных материалов // Лакокрасочные материалы и их применение. – 1996. – № 10. – С. 37–38.
4. Степин С.Н. Организация производства и оборудование для получения пигментированных лакокрасочных материалов: учеб. пособие. – Казань: Казан. гос. технол. ун-т, 2005. – 104 с.
5. Таршис М.Ю. К расчету смесителя сыпучих материалов со сложным движением эластичной рабочей камеры // Изв. ВУЗов, «Химия и химическая технология». – 2008. – т. 51, – Вып. 8. – С. 75–77.

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

Процесс получения двухфазной системы жидкость – твердые частицы в аппарате роторного типа осуществляется в следующей последовательности. Вначале в рабочую емкость вводят жидкий компонент. Затем при медленном вращении вводят порошкообразный компонент. Затем скорость вращения вала резко увеличивают до образования воронки (эффект Dougnut). Расчетная схема движения представлена на рис. 1.

pic_36.tif

Рис. 1. Расчетная схема движения в аппарате роторного типа

Для удобства описания процесса смешения перейдем в систему координат, вращающуюся вместе с валом. В выбранной системе координат при вращении турбинной мешалки выделяются два циркуляционных контура: верхний и нижний (рис. 1). Верхний циркуляционный контур можно аппроксимировать в виде треугольника. Соответственно, линии тока в данном контуре будут представлять собой также треугольники.

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

Поток массы твердых частиц через поверхность циркуляционного контура определяется выражением [2]:

shirina01.wmf (1)

где c – концентрация твердой фазы; ⟨D⟩ – среднестатистический коэффициент макродиффузии; x – координата, определяющая положение циркуляционного контура; S(x) – площадь поверхности циркуляционного контура, которая определяется уравнением

shirina02.wmf (2)

где shirina03.wmf shirina04.wmf shirina05.wmf (3)

l – максимальное значение координаты x, которое будет определяться центром циркуляции.

Изменение количества твердой фазы в единицу времени в циркуляционном слое толщиной dx будет определяться разностью потоков массы:

shirina06.wmf (4)

Выполняя предельный переход, получим следующее уравнение:

shirina07.wmf (5)

Граничное условие составляется из тех соображений, что на границе циркуляционных контуров у твердой стенки перенос твердой фазы во внешнюю среду отсутствует, и записывается в виде

shirina08.wmf (6)

Начальное условие строится из предположения, что в начале процесса твердые частицы подаются сверху на движущийся слой жидкости и, следовательно, полностью заполняют все циркуляционные контуры, граничащие со стенкой аппарата, и записывается в виде

shirina09.wmf (7)

где значение координаты x1 определяется из равенства

shirina10.wmf (8)

Нетрудно убедиться, что уравнение (5) полностью соответствует условию сохранения массы

shirina11.wmf (9)

где M – масса твердой фазы; ρт – истинная плотность твердой фазы.

Среднестатистический коэффициент макродиффузии ⟨D⟩ существенно зависит от размера частиц твердой фазы. Чем меньше размеры частиц твердой фазы, тем выше их подвижность в хаотическом движении. В связи с этим величина ⟨D⟩ может быть представлена в виде

shirina12.wmf (10)

где D0 – коэффициент макродиффузии при размере частиц размером d0; shirina13.wmf – средний размер частиц твердой фазы, k – некоторая константа, определяемая в рамках модели. Экспериментальные исследования показывают, что величина коэффициента макродиффузии D0 определяется числом оборотов ротора и конструктивными особенностями диспергатора.

Предполагаем, что процесс смешения практически не влияет на процесс диспергирования твердой фазы, так как за счет высокой турбулентности распределение частиц по объему жидкой фазы практически равномерно.

В работе [2] показано, что процесс диспергирования твердой фазы может быть описан уравнением

shirina14.wmf (11)

где kg – коэффициент скорости диспергирования; n – константа, определяемая природой материала твердой фазы; d0 – начальный средний размер частиц твердой фазы. Параметры уравнения (10) определяются регрессионными методами на основе обработки экспериментальных данных.

Насосная производительность в верхнем контуре определяется известной формулой [1]:

shirina15.wmf (12)

где n – число оборотов ротора; d – диаметр мешалки. В соответствии с формулой (9) время движения в верхнем циркуляционном контуре будет равно

shirina16.wmf (13)

где Vв – объем верхнего циркуляционного контура.

Для решения дифференциального уравнения в частных производных (5) относительно c(x, t) будем использовать численный анализ для значений x ∈ [0, l] и t ∈ [0, T] методом сеток, сделав следующие разбиения x и t с шагами, которые определяются следующими формулами:

hx = l/m; ht = τв, (14)

где m – соответственно число точек для переменных x и t, которое определяется качеством смеси, о котором будет сказано позднее; определяется формулой (13).

В соответствии с численным анализом используем приближенные выражения для производных:

shirina17.wmf (i = 0, 1, 2, ..., n; j = 0, 1, 2, ..., m); (15)

shirina18.wmf (i = 0, 1, 2, ..., n; j = 0, 1, 2, ..., m); (16)

shirina19.wmf (i = 0, 1, 2, ..., n; j = 0, 1, 2, ..., m). (17)

Подставляя формулы (15)–(17) в уравнение (5), получим:

shirina20.wmf (18)

Решение модели осуществляется численно стандартным сеточным методом. Время перемешивания T (число N) определяется качеством смеси, которое оценивается по величине коэффициента неоднородности:

shirina21.wmf (19)

где k – количество проб; ci – концентрация твердой фазы в i-м циркуляционном контуре; c0 – среднее значение концентрации твердой фазы по объему аппарата.

На рис. 2 представлено сопоставление теоретических и расчетных данных по качеству смеси при различных значениях времени процесса смешения.

pic_37.wmf

Рис. 2. Зависимость качества смешения от времени процесса перемешивания

Расчет производился для получения двухфазной системы вода ‒ кальцит с процентным содержанием кальцита 40 %. В процессе диспергирования отношение shirina22.wmf составило 0,76. Процесс производился при следующих параметрах аппарата: D = 1,5 м, H = 2,0 м, d = 0,4 м, h = 0,15 м. При числе оборотов ротора n = 1500 об/мин коэффициент диффузии составил ⟨D⟩ = 0,57⋅10–3, k = 1,27.

Параметры модели дезагрегации твердых частиц определяются методами статистического анализа на основе экспериментальных данных в программе Statistica 6.0. В результате получено уравнение

shirina23.wmf (20)

где τ – время диспергирования (с).

Значение коэффициента неоднородности рассчитывалось по массе твердой фазы. Пробоотборником отбиралось 30 проб из различных точек объема смеси, положение которых определялось расстоянием от стенки аппарата. Все точки находились на одной высоте над мешалкой со значением равным величине 0,3 м. Масса кальцита в пробе определялась путем полного выпаривания воды. Дисперсность кальцита определялась микроскопическим методом. Расхождение теоретических и экспериментальных данных не превышает 12 %.

Рецензенты:

Бачурин В.И., д.ф.-м.н., профессор кафедры «Высшая и прикладная математика» Ярославского филиала МИИТ, г. Ярославль;

Смирнов А.Д., д.ф.-м.н., профессор Ярославского государственного университета им. П.Г. Демидова, г. Ярославль.

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


Библиографическая ссылка

Ширина Н.Ю., Сидоров В.Н., Мурашов А.А. МОДЕЛИРОВАНИЕ ПРОЦЕССА ПОЛУЧЕНИЯ ДВУХФАЗНОЙ СИСТЕМЫ ЖИДКОСТЬ – ТВЕРДОЕ В АППАРАТЕ РОТОРНОГО ТИПА // Фундаментальные исследования. – 2014. – № 8-6. – С. 1341-1344;
URL: http://www.fundamental-research.ru/ru/article/view?id=34765 (дата обращения: 14.12.2019).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074