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

МОДЕЛЬ НЕСТАЦИОНАРНОГО ПОТОКА СУСПЕНЗИИ В ПОРИСТОЙ СРЕДЕ

Жеглова Ю.Г. 1 Сафина Г.Л. 1
1 Национальный исследовательский Московский государственный строительный университет
Изучение фильтрации суспензии в пористой среде является актуальной проблемой при проектировании и строительстве туннелей, подземных и гидравлических сооружений, систем очистки сточных вод и жидких промышленных отходов. Рассматривается математическая модель переноса и осаждения твердых частиц при фильтрации нестационарного потока монодисперсной суспензии в пористой среде с механико-геометрическим механизмом захвата частиц. Построено точное решение задачи на входе пористой среды. Проводится расчет полного осадка на входе пористой среды. Рассмотрены различные распределения концентрации взвешенных частиц на входе пористой среды для линейного блокирующего коэффициента фильтрации. Показано, что предельное значение полного осадка не зависит от формы распределения концентрации взвешенных частиц на входе фильтра. Величина предельного суммарного осадка определяется объемом взвешенных частиц, впрыскиваемых в пористую среду. Этот максимум достигается за конечное время для конечного распределения времени.
фильтрация
суспензия
пористая среда
математическая модель
численное решение
1. C. Tien. Principles of filtration. – Elsevier, Oxford, 2012.
2. D.J. Logan. Transport Modelling in Hydrogeochemical Systems. – Springer, N.Y., 2001.
3. Bedrikovetsky P. Upscaling of Stochastic Micro Model for Suspension Transport in Porous Media // Transport in Porous Media. – 2008. – vol. 75. – Р. 335–369.
4. You Z., Bedrikovetsky P., Kuzmina L. Exact Solution for Long-Term Size Exclusion Suspension-Colloidal Transport in Porous Media // Abstract and Applied Analysis. – 2013. – vol. 2013, iss. Mathematical and Computational Analyses of Flow and Transport Phenomena, 9 p.
5. You Z., Osipov Y., Bedrikovetsky P., Kuzmina L. Asymptotic model for deep bed filtration // Chemical Engineering Journal. – 2014. – vol. 258. – Р. 374–385.
6. Osipov Yu., Kuzmina L. Calculation of the filtration of polydisperse suspension with a small rate // Matec Web of Conferences. – 2017. – vol. 117, 00131. – 6 p.
7. Kuzmina L.I., Osipov Yu.V. Asymptotic solution for deep bed filtration with small deposit // Procedia Engineering. – 2015. – vol. 111. – Р. 491–494.
8. Kuzmina L., Osipov Yu. Filtration model of the unsteady suspension flow in a porous medium // Matec Web of Conferences. – 2017. – vol. 117, 00097. – 6 p.
9. Kuzmina L.I., Osipov Yu.V., Galaguz Yu.P. A model of two-velocity particles transport in a porous medium // International Journal of Non-linear Mechanics. – 2017. – vol. 93. – Р. 1–6.
10. Galaguz Y.P., Safina G.L. Modeling of Particle Filtration in a Porous Medium with Changing Flow Direction // Procedia Engineering. – 2016. – vol. 153. – Р. 157–161.
11. Galaguz Y.P., Safina G.L. Modeling of fine migration in a porous medium // MATEC Web of Conferences. – 2016. – vol. 86, 03003. – 6 p.
12. Galaguz Y., Safina G. Calculation of the filtration in a heterogeneous porous medium // Matec Web of Conferences. – 2017. – vol. 117, 00052. – 6 p.
13. Galaguz Y., Safina G. Modeling of filtration of 2-types particles suspension in a porous medium // Matec Web of Conferences. – 2017. – vol. 117, 00053. – 6 p.

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

В работе рассматривается математическая модель движения частиц в фильтре, основанная на механико-геометрическом взаимодействии частиц с пористой средой [3]. Предполагается, что частицы свободно проходят через поры большого размера и застревают на входе пор, если диаметр частицы больше диаметра поры. Осажденная частица не может быть выбита из поры потоком жидкости или другой частицей, и навсегда остается в поре. С ростом осадка количество свободных малых пор сокращается и скорость образования осадка замедляется. Если все поры малого размера заблокированы частицами, концентрация осадка S достигает максимально возможного предельного значения SM. Коэффициент пропорциональности Λ(S) между скоростью роста осадка и концентрацией взвешенных частиц, обращающийся в ноль при блокировании всех малых пор S = SM, называется блокирующим коэффициентом фильтрации.

Для различных моделей фильтрации стационарного потока суспензии в пористой среде получены точные и асимптотические решения [4–9]. При отсутствии аналитических решений используются численные расчеты [10–13].

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

Постановка задачи

Рассмотрим основные уравнения, описывающие динамику взвешенных и осажденных частиц суспензии в фильтре длиной 1. Обозначим через C(x, t) концентрацию взвешенных частиц, пусть S(x, t) – концентрация осажденных частиц. Концентрации частиц удовлетворяют уравнению непрерывности, которое с учетом неподвижности осадка в простейшем случае имеет вид

geg01.wmf. (1)

Рост осадка пропорционален концентрации взвешенных частиц, из которых он образуется, а также зависит от текущей величины осадка S(x, t): чем больше пор малого диаметра закрыто частицами, тем меньше скорость роста осадка

geg02.wmf. (2)

Здесь Λ(S) > 0 при S ≥ 0.

Система уравнений (1), (2) рассматривается в области geg03.wmf. Краевые условия для системы (1), (2) ставятся на входе фильтра x = 0 и в начальный момент времени t = 0:

geg04.wmf; (3)

geg05.wmf; geg06.wmf. (4)

Условия (4) означают, что в начальный момент пористая среда не содержит взвешенных и осажденных частиц. Суспензия переменной концентрации начинает поступать на вход фильтра в момент t = 0 и движется в пористой среде со скоростью v = 1. Характеристика t = x – фронт концентраций взвешенных и осажденных частиц является подвижной границей суспензии. В области пустой пористой среды geg07.wmf решение нулевое; в области суспензии geg08.wmf решение положительно.

Решение на входе фильтра

Уравнение (2) на входе фильтра x = 0 имеет вид

geg09.wmf. (5)

Делим обе части уравнения (5) на Λ(S) и интегрируем по переменной t:

geg10.wmf. (6)

Используя условие (4), преобразуем интеграл в левой части (6):

geg11.wmf. (7)

Формула (7) задает концентрацию осажденных частиц на входе фильтра.

Рассмотрим модель с линейным блокирующим коэффициентом фильтрации:

geg12.wmf. (8)

В этом случае согласно (7) концентрация осадка на входе фильтра

geg13.wmf. (9)

Согласно [3, 4] для коэффициента фильтрации (8) концентрация осажденных частиц в области geg14.wmf равна

geg15.wmf. (10)

Расчет полного осадка

Объем осажденных частиц (полный осадок) в пористой среде в момент времени t равен

geg16.wmf. (11)

Если суспензия впрыскивается в пористую среду в течение конечного промежутка времени t1 ≤ t ≤ t2, то Sf(t) = 0 при t ≤ t1. Поскольку частицы проходят от входа до выхода пористой среды за время t0 > 1, то Sf(t) возрастает при t1 < t < t2 + 1 и Sf(t) = Sf(t2 + 1) = const при t > t2 + 1. Если впрыск суспензии происходит при всех значениях времени t > 0, то функция Sf(t) возрастает при t > 0.

Численный расчет выполнен для линейного коэффициента фильтрации (8) при geg17.wmf для трех распределений pi(t) концентрации взвешенных частиц на входе фильтра, заданных на рис. 1.

gekl1.wmf

Рис. 1. Концентрации взвешенных частиц на входе фильтра

gekl2a.tif gekl2b.tif

gekl2c.tif

Рис. 2. Полный осадок в пористой среде для трех граничных распределений

На рис. 2 изображены графики полного осадка geg18.wmf.

Общее количество взвешенных частиц, впрыскиваемых в пористую среду, является одинаковым для всех трех рассматриваемых распределений pi(t), поскольку geg19.wmf (рис. 1). Для финитных распределений концентрации взвешенных частиц p1(t), p2(t) на входе фильтра полный осадок достигает максимального значения при t = 2; для распределения p3(t) процесс фильтрации длится неограниченно, полный осадок возрастает и стремится к предельному максимальному значению при t → ∞ (рис. 2).

Заключение

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

Форма распределения концентрации взвешенных частиц на входе фильтра pi(t) определяет вид графика geg20.wmf. Поскольку geg21.wmf при 0 < t < 1, то на этом интервале график geg22.wmf выше графика geg23.wmf.

Согласно рис. 2 предельный полный осадок не зависит от формы распределения концентрации взвешенных частиц на входе фильтра. Величина предельного полного осадка определяется объемом взвешенных частиц, впрыскиваемых в пористую среду.

Найденное точное решение позволит упростить нахождение параметров процесса фильтрации в лабораторных и полевых условиях.


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

Жеглова Ю.Г., Сафина Г.Л. МОДЕЛЬ НЕСТАЦИОНАРНОГО ПОТОКА СУСПЕНЗИИ В ПОРИСТОЙ СРЕДЕ // Фундаментальные исследования. – 2017. – № 11-1. – С. 63-66;
URL: https://fundamental-research.ru/ru/article/view?id=41899 (дата обращения: 02.12.2021).

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

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