сотрудник
Исследования проводились с целью математического моделирования процесса фильтрования тонкодисперсных неньютоновских суспензий в намывных фильтрах периодического действия. Для описания реологического состояния среды использована модель степенной жидкости. В начале цикла фильтрования намывается слой вспомогательного материала, который используется в качестве фильтрующей перегородки. Тонкодисперсная суспензия в процессе фильтрования образует второй слой, по мере нарастания которого возрастает сопротивление, оказываемое потоку жидкой фазы. Фильтровальный цикл завершается после достижения критических значений давления, либо производительности в зависимости от режима работы аппарата. Уравнения фильтрации записаны в цилиндрической системе координат отдельно для каждого слоя. На границе раздела слоев учтены скачки давления и скорости, возникающие из-за отличия порозности. Для расчета толщины образуемого слоя осадка построена дифференциальные уравнение. Показано, что построенная математическая модель может быть использована для расчета фильтрования суспензий с линейной моделью реологического состояния, если выполнить предельный переход n→1. Также показано, что после предельного перехода R→∞, результаты работы позволяют моделировать процесс фильтрования неньютоновских сред в фильтрах пластинчатого типа. Численные расчеты проведены для режима фильтрования с постоянной скоростью. Исследованы закономерности влияния эквивалентной вязкости, степени нелинейности среды, а также исходной концентрации дисперсной фазы на процесс фильтрования. Показано, что при увеличении этих параметров предельное значение давления в аппарате достигает раньше, что приведет к сокращению времени фильтровального цикла.
фильтр, степенная модель жидкость, фильтрование суспензии, слой осадка
Введение. Фильтрование жидкостей с тонкодисперсными примесями встречается во многих отраслях народного хозяйства [1, 2, 3]. В качестве примера можно назвать проблему восстановления технических масел, подготовку технической воды промышленности, очистку масел растительного происхождения и сгущения фруктовых соков в пищевой промышленности. Одним из методов фильтрования тонкодисперсных суспензий является использование вспомогательных фильтрующих материалов. Вспомогательный фильтрующий материал может быть добавлен в исходную суспензию с целью повышения проницаемости слоя образуемого осадка. Тогда рассматривается процесс фильтрования трехфазной гетерогенной среды. Вспомогательный фильтрующий материал может быть использован в качестве фильтрующей перегородки. В этом случае он намывается отдельно в виде нижнего слоя в начале цикла фильтрования. Тонкодисперсная суспензия фильтруется через нижний слой, образуя слой тонкодисперсных частиц. При этом возникает задача расчета фильтрации среды через двойной слой осадка, когда толщины верхнего слоя является переменной величиной. В настоящее время изучение процесса фильтрования неньютоновских сред с тонкодисперсными частицами на многослойных проницаемых подложках остается актуальной задачей [4, 5, 6].
Цель исследований - математическое моделирование и численный расчет процесса фильтрования тонкодисперсных неньютоновских суспензий, подчиняющихся степенной модели реологического состояния среды, в намывных фильтрах периодического действия.
Условия, материалы и методы. Многие аномально вязкие среды хорошо описываются моделью степенной жидкости. При этом обобщенное уравнение фильтрации имеет вид [7, 8, 9]
(1)
В данном уравнении параметр n указывает на степень нелинейности среды. При значении n=1 получаем классическое уравнение фильтрации Дарси. Через коэффициенты Ki – обозначены проницаемости слоев пористого осадка. Величина Ψ называется эквивалентной вязкостью. Она зависит не только от реологических параметров жидкости, но и от пористости и проницаемости соответствующих слоев осадка. Для вычисления эквивалентной вязкости известны разные формулы, которые справедливы при определенных диапазонах изменения скорости и конкретных моделей пористой среды [9].
Существуют разные конструкции намывных фильтров периодического действия – пластинчатые, трубчатые [10, 11, 12]. Возможны разные режимы их работы: фильтрование при постоянной разности давлений, фильтрование при постоянной скорости, а также фильтрование при переменных разности давлений и скорости. В этих аппаратах предварительно намывается нижний фильтрующий слой с достаточно высокой проницаемостью, через которого фильтруется тонкодисперсная суспензия. По мере увеличения толщины второго слоя возрастает сопротивление, оказываемое
потоку жидкой фазы. В режиме фильтрования при постоянной разности давлений уменьшается производительность, в режиме фильтрования при постоянной скорости увеличивается разность давлений. После достижения их предельных значений цикл работы аппарата прекращается. Далее выполняется промывка фильтрующего элемента, и далее цикл повторяется заново.
Результаты и обсуждение. Пусть на поверхности трубчатого фильтра образован слой вспомогательного материала толщины δ = const. В процессе фильтрования тонкодисперсной суспензии происходит накопление второго слоя осадка h(t). Уравнения фильтрации (1) для соответствующих слоев в цилиндрической системе координат имеют вид
(2)
(3)
Запишем уравнения сохранения массы
(4)
(5)
Для решения уравнений фильтрации и сохранения массы необходимо задавать граничные условия. В рабочей зоне фильтра создается давление P, которое определяется с учетом режима работы аппарата. В математической модели это давление задается на поверхности верхнего слоя осадка
(6)
На границе раздела двух слоев задаются условия сшивания для давления и скорости. Поскольку пористости разных слоев осадка различаются, на границе раздела произойдет скачкообразное изменение давления и скорости фильтрации. Эти условия могут быть записаны в виде [13]
(7)
где εδ, εh- пористости нижнего и верхнего слоев осадка.
На поверхности фильтрующего трубчатого элемента задается выходное давление
(8)
Решение системы уравнений (2), (3), (4) и (5) при заданных граничных условиях (6), (7), (8) не представляют труда. После интегрирования уравнения (4) имеем
(9)
После подстановки зависимости (9) можно проинтегрировать уравнения фильтрации (2)
Постоянного интегрирования С2 определим из условия (8). Тогда последнее соотношение примет вид
Отсюда найдем
(10)
Аналогично проинтегрируем уравнения (3) и (5). Если интеграл уравнения (5) представим в виде
то решения уравнения (3) ,с учетом граничного условия (6), запишется так
(11)
Условие сшивания скоростей εhVh=εδVδ на границе r=R+δ дает зависимость εhC3=εδC1. При этом постоянная интегрирования C3 выражается через константу C1 в виде:
Тогда зависимость (11) можно представить следующим образом
(12)
Полученные соотношения для давлений (10) и (12) должны выполнять граничное условие (7), с помощью которого определяется постоянное интегрирования С1. После определения С1, в результате несложных преобразований, для расчета скоростей фильтрации окончательно получим
(13)
(14)
Уравнения (13), (14) позволяют вычислять скорость фильтрации сплошной фазы суспензии в аппаратах трубчатого типа. Полученные уравнения могут быть преобразованы для расчета фильтровальных аппаратов пластинчатого типа. Для этого необходимо выполнить предельный переход R→∞. Предварительно проведем необходимые преобразования.
Несложно показать, что
Выполним предельный переход R→∞, используя замечательный предел
после замены переменных 1/R=x или 1/(R+δ)=x. С учетом того, что при условии R→∞ имеет место предел r/R→1, получим
Тогда, уравнение (13) для фильтровальных аппаратов пластинчатого типа примет вид
(15)
а уравнение(14) остается в силе.
Если фильтруемая среда подчиняется линейному реологическому закону состояния, то уравнения (13) также необходимо преобразовывать, используя предельный переход n→1. С помощью замечательного предела
можно показать, что имеет место равенства
Тогда уравнение (13) преобразуется к виду
(16)
Построим уравнение для расчета толщины слоя осадка. При заданной концентрации дисперсных частиц άT, их приток на элементарную площадку ∆S=(R+δ+h)∆l составляет
Продифференцируем данное соотношение с учетом условия (1-άT)V=εhVh на границе осадок-суспензия
(17)
С другой стороны за время ∆t объем слоя осадка на единичной площадке ∆S изменится на величину ∆q=(1—εh)∆S∆h
Разделим обе части данного равенства на ∆t и выполним предельный переход ∆t→0:
(18)
Приравнивая правых частей соотношений (17) и (18), для расчета изменения толщины верхнего слоя осадка получим дифференциальное уравнение
(19)
Построенное уравнение решается с начальным условием h=0 при t=0.
Численные расчеты по построенным математическим моделям были проведены в безразмерных переменных
где V* - характерная скорость, которая вычисляется через число Рейнольдса. В частности, после перехода к безразмерным переменным, уравнение (13) примет следующий вид (знаки тильда ~ опущены)
,
где
При этом вид уравнения (19) не изменится.
Некоторые результаты расчетов, для режима фильтрования с постоянной скоростью, приведены на рисунках 1-3. В режиме фильтрования с постоянной скоростью, по мере увеличения толщины осадка, происходит возрастание давления в аппарате.
На процесс фильтрации сильно влияют реологические параметры суспензии Ψi и n. Эквивалентная вязкость среды обратно пропорциональна числу Рейнольдса Re. На рис. 1 показано влияние параметра Re на интенсивность возрастания давления. Как видим, более интенсивный рост давления наблюдается при малых значениях числа Re. Следовательно, при увеличении эквивалентной вязкости среды, предельное значение давления в аппарате достигает раньше, что приведет к сокращению времени фильтровального цикла. Влияние степени нелинейности среды n показано на рис. 2. При увеличении этого параметра давление в аппарате растет быстрее, следовательно, сокращается время достижения критического давления.
При фильтровании суспензий с повышенной концентрации нарастание слоя осадка происходит более интенсивно. В результате этого быстрее растет давление в аппарате и сокращается длительность фильтровального цикла (рис.3).
Выводы. Построена математическая модель фильтрования тонкодисперсных суспензий, подчиняющихся степенной модели реологического состояния среды, в намывных фильтрах периодического типа. Для режима фильтрования с постоянной скоростью приведены результаты численных расчетов. Установлены закономерности влияния эквивалентной вязкости, степени нелинейности среды, а также исходной концентрации дисперсной фазы, на длительность фильтровального цикла.
.
1. Ахмадиев Ф.Г., Фарахов М.И., Бекбулатов И.Г., Исянов Ч.Х. Математическое моделирование процесса фильтрования двухфазных суспензий в трубчатых фильтрах в неизотермических условиях // Теоретические основы химической технологии. - 2016. - Т. 50. - № 1 - С. 44.
2. Ибятов, Р. И. О моделировании случайных процессов в агропромышленном комплексе / Р. И. Ибятов, Б. Г. Зиганшин // Вестник Казанского государственного аграрного университета. - 2022. - Т. 17. - № 1(65). - С. 50-55.
3. Багайсков Ю.С. Разработка состава композиционного материала трубчатых фильтров для очистки технических жидкостей // Технология машиностроения. - 2013. - № 6. - С. 41-43.
4. Булышев Е.М., Худобин Л.В. Высокопроизводительные намывные фильтры // Водоочистка. - 2019. - № 6. - С. 51-62.
5. Ибятов Р.И., Холпанов Л.П., Ахмадиев Ф.Г. Течение многофазной среды по проницаемой поверхности с образованием осадка // Инженерно-физический журнал. - 2005. - № 2. - С. 65-72.
6. Ибятов Р.И., Холпанов Л.П., Ахмадиев Ф.Г., Бекбулатов И.Г. Математическое моделирование течения многофазной гетерогенной феды по проницаемой трубе // Теоретические основы химической технологии. - 2005. - № 5. - С. 533-541.
7. Ю. Ф. Лачуга, Р. И. Ибятов, Б. Г. Зиганшин. Моделирование траектории движения зерна по рабочим органам пневмомеханического шелушителя // Российская сельскохозяйственная наука. - 2020. - № 4. - С. 73-76.
8. Жуков В.Г., Чесноков В.М., Лукин Н.Д. Получение кумулятивной функции удельной поверхности пористых сред по заданной гистограмме частотного распределения размеров пор // Достижения науки и техники АПК. - 2019. - Т. 33. - № 11. - С. 82-87.
9. Бернадинер М.Г., Ентов В.М. Гидродинамическая теория фильтрации аномальных жидкостей. М.: Наука, 1975. - 200 с.
10. Балашов В.А., Тябин Н.В. Фильтрация неньютоновских жидкостей, подчиняющихся степенному реологическому закону // Теоретические основы химической технологии. - 1989. - № 6. - С. 844-846.
11. Севинс Дж. Неньютоновское течение в пористой среде // Механика: Сборник переводов. М.: Мир. 1974. Вып. 2. - С. 59-115.
12. Гольдштик М.А. Процессы переноса в зернистом слое. Новосибирск, 1984. - 164 с.
13. Ю. Ф. Лачуга, Р. И. Ибятов, Ю. Х. Шогенов // Метод расчета траектории движения зерна в пневмомеханическом шелушителе // Российская сельскохозяйственная наука. - 2021. - № 6. - С. 64-67.
