Particle dynamics with a viscous liquid through a porous medium cell with a pressure gradient

Full Text

Введение Проблема моделирования течения вязкой жидкости в пористой среде является одной из актуальных в механике жидкости. Причем это обусловлено как математической сложностью самой задачи, так и чрезвычайно широкой областью возможных приложений результатов моделирования: извлечение нефти и газа из пласта, водоотведение и водоотбор, создание бытовых и промышленных фильтров для очистки жидкостей и т. п. В мировой практике используются различные подходы для моделирования движения таких сред. Наиболее распространенный подход основан на уравнении Дарси и его модификациях [1-4]. Полученная система уравнений решается аналитически или численно при заданных значениях коэффициентов, характеризующих параметры среды, и которые определяются экспериментально. Однако для реальных пористых сред экспериментальное нахождение этих параметров - достаточно сложная задача, так как такие системы уникальны по своим свойствам, и, кроме того, достаточно трудно воспроизвести в лаборатории условия, характеризующие такую систему в естественном состоянии. Развитие компьютерных технологий дало возможность использовать и другие подходы к моделированию фильтрации жидкости в пористой среде. Это численное моделирование, позволяющее рассматривать течения в пористой среде на масштабе пор [5-7]. В этом случае задача сводится к прямому численному решению уравнений гидродинамики в некоторой пористой среде, геометрия которой воспроизводит структуру реального пористого материала. При этом используются известные методы численного решения уравнений гидродинамики и методы компьютерной томографии образца материала для моделирования его структуры. Проблема данного подхода в том, что для его реализации требуются значительные вычислительные ресурсы. Так, например, для цифрового представления реального образца пористого материала в виде кубика размерами в 600 мкм по каждой стороне требуется 109 пикселей [8]. Помимо этого, реализация самих численных методов интегрирования уравнений гидродинамики с граничными условиями на сложных геометрических поверхностях также является сложной задачей. Сложность задачи возрастает, когда в жидкости присутствуют подвижные диспергированные частицы и расчет ведется на реальном масштабе пор. Даже в последнее время, когда возможности вычислительной техники выросли, применение такого метода для большого числа обтекаемых поверхностей является трудной задачей, требующей большого объема вычислений. Поэтому, как указывается в работе [8], получение уравнений фильтрационных течений вязкой жидкости с использованием средних величин, характеризующих поток в пористой среде, по-прежнему является актуальной задачей. Фактически это означает, что необходимо исследовать процессы массо- или теплопереноса на масштабе пористой ячейки, а затем получить средние значения параметров на масштабе, много большим, чем размер ячейки, но много меньшим, чем размер пласта. Таким образом, в настоящее время не существует достаточно точного и при этом несложного метода, позволяющего моделировать фильтрацию жидкости с частицами с учетом их взаимодействия между собой и пористой структурой. Трудности моделирования возрастают при увеличении числа частиц, взаимодействующих между собой, что сказывается на практической реализации вычислительных схем известных методов. Поэтому получение новых аналитических и численных результатов в этой области по-прежнему остается актуальной задачей. В работе [12] на основе результатов, ранее полученных авторами, предлагается метод, основанный на представлении динамики частиц в пористой среде, как динамики системы частиц при учете как внутренних сил, так и сил гидродинамического взаимодействия между ними. Пористая среда моделируется большим числом неподвижных полидисперсных частиц, имитирующих структуру пористой среды. Движение частиц, моделирующих включения, происходит в результате течения жидкости через модельную структуру. В расчет берется взаимодействие модельных частиц между собой и с частицами структуры. Постановка задачи и метод решения Рассматривается модель фиктивного грунта, состоящего из N1 твердых сферических частиц, образующих случайную структуру в объеме V, и помещенных в несжимаемую жидкость плотности ρ и вязкости η. Частицы имеют эффективный радиус a. Предполагается, что на частицы наложены связи, препятствующие их поступательному и вращательному движению. Как возможный вариант постановки задачи - наличие вдали или вблизи от рассматриваемой структуры частиц в жидкости плоской стенки в качестве модели непроницаемой поверхности, вдоль которой происходит течение со скоростью U. Внутри структуры помещаются N2 частиц меньшего радиуса b, чем эффективный радиус частиц a, и которые могут поступательно и вращательно перемещаться в жидкости. Эти частицы моделируют наличие второй фазы в жидкости. Для расчета течения жидкости внутри рассматриваемой структуры и динамики системы свободных частиц записываются уравнения движения жидкости и каждой частицы с учетом всех сил и моментов, связанных с гидродинамическим взаимодействием всех частиц и сил реакции связей. Рассматривается случай малых чисел Рейнольдса. Система уравнений для вязкой жидкости имеет вид: (1) На поверхности неподвижной частицы с номером должны выполняться следующие граничные условия: (2) На поверхности подвижной частицы с номером должны выполняться следующие граничные условия: (3) Далеко от частиц имеет место затухание возмущений: (4) Здесь - скорость жидкости, и - соответственно, поступательная и угловая скорости частицы с номером , - невозмущенное давление в жидкости, - вектор положения произвольной точки среды относительно центра частицы с номером . Уравнения динамики для подвижных частиц записываются в виде (квазистационарное приближение): (5) здесь - сумма сил, - сумма моментов сил, действующих на подвижную частицу с номером со стороны жидкости. Для неподвижных частиц имеем соотношения: (6) здесь - сумма внешних сил, - реакции связей, удерживающих частицы в заданном положении, - сумма моментов внешних сил, - сумма моментов реакций связей, наложенных на неподвижную частицу с номером Метод решения системы уравнений (1), (5), (6) с граничными условиями (2), (3), (4), а также процедура получения гидродинамических сил и моментов, действующих на частицы со стороны жидкости с учетом гидродинамического взаимодействия всех частиц, представлены в работах [9-11]. Скорости подвижных частиц находятся из условия (5), а реакции связей - из (6). Ввиду большого количества уравнений 6(N1+N2)+4 был применен специальный программный комплекс, аналогичный тому, что использовался в работах [10-12]. Полученная система уравнений решается численным способом, и по найденным числовым значениям находятся скорости свободных частиц, строятся траектории их движения частиц, а также определяются усилия, возникающие в связях, удерживающих частицы в структуре. Моделирование динамики частиц при фильтрационном течении Моделирование динамики частиц при течении вязкой жидкости через пористую среду основывается на методе подобия, разработанном в работах [10]. Вместо реальных размеров частиц a и таких параметров жидкости, как её вязкость η и скорость течения, при компьютерном моделировании задаются модельные значения размеров частиц , вязкости жидкости и скорости внешнего потока . Это связано с тем, что значения реальных параметров достаточно малы, что затрудняет нахождение численного решения системы уравнений и визуализацию результатов. Поскольку плотность жидкости явным образом не входит в уравнения движения жидкости и частиц, то ее значение считается одинаковым для модельной и реальной жидкостей. Условие подобия реального и модельного течений означает, что число Рейнольдса должно быть одинаковым в обоих случаях. Откуда получаем следующее равенство: (7) Из соотношения (7) получаем: (8) Поскольку программа для моделирования течения жидкости через упаковку частиц использует промежуток времени, то необходимо найти связь между такими параметрами, как характерное время в реальном процессе фильтрации t=a/v и при моделировании . Используя соотношение (8), получим такую зависимость между указанными параметрами: . Зная время движения частиц в пористой структуре при моделировании, можно найти время движения реальной частицы при фильтрации. Для определения модельной скорости течения жидкости внутри структуры использовался следующий подход. Из решения системы уравнений для модельных частиц, образующих пористую структуру, получаем значения составляющих реакций связей по оси k для частицы с номером i, удерживающих частицы на месте. Причем эти значения учитывают гидродинамическое взаимодействие всех частиц, как неподвижных N, так и подвижных n. Для каждой неподвижной частицы записываем соотношение: (9) Здесь - составляющая скорости потока жидкости для частицы с номером i, учитывающая гидродинамическое взаимодействие всех частиц. Другими словами, - скорость самосогласованного течения жидкости, образованная основным потоком и возмущениями от всех остальных частиц. Поскольку для каждой частицы из структуры можно получить такого рода течения, то для всей структуры в целом можно считать, что имеется течение со скоростью потока жидкости внутри структуры, определяемое как средняя скорость течений для всех частиц, которая находится из соотношения: (10) С учетом равенства (9) из (10) получим выражение для средней скорости потока жидкости внутри структуры: (11) Аналогичные выражения можно записать и для величин реального фильтрационного потока. Учитывая равенства (8), (9), можно получить следующее соотношение для величин и в модельном и реальном течениях: (12) Выражение (11) можно представить и в другом виде. Вводя такие величины, характеризующие структуру, как пористость , удельную поверхность частиц , составляющую напряжение на поверхности частиц , получим соотношение: . Здесь учтены следующие соотношения: . Необходимо отметить, что параметр пористости для модельной и реальной структуры одинаковый, так как не зависит от масштаба, а только от геометрии структуры. Возможны два варианта моделирования: зная реальные значения параметров, при моделировании можно задавать значения двух параметров, а значение третьего находится из условия подобия или, задавая значения всех трех модельных параметров, можно найти один из трех реальных параметров, зная значения двух из них. Результаты моделирования Для моделирования пористой структуры и динамики частиц внутри ее была использована программа, аналогичная той, что использовалась в [10-12]. В заданном объеме пространства, ограниченном координатами , , , случайным образом располагаются частицы. Из них N частиц неподвижны и имеют радиус а одна, моделирующая подвижную частицу, располагается внутри структуры и имеет радиус . X3 X1 X2 Рисунок 1 - Структура пористой среды, образованной неподвижными частицами Рассматривались два модельных случая N и подвижной частицы: N=449, k = 0,3; N=477, k = 0,2. На рисунке (рис. 1) показана характерная структура, образуемая неподвижными частицами. Задавались следующие модельные значения: скорость потока вдоль оси OX2, вязкость жидкости , начальное положение подвижной частицы и промежуток времени Δ t= 0,1. Скорость течения выбрана таким образом, чтобы на границах структуры , имитирующих плоские поверхности, жидкость покоилась. Заданному потоку жидкости соответствует градиент давления вдоль оси OX2, равный . Затем запускалась программа. Результаты расчетов динамики частиц представлены ниже на рисунках. На рис. 2 представлена динамика частицы для случая N=449, k = 0,3. Начальное положение подвижной частицы задавалось координатами (0; -7,5; 0). Рисунок 2 - Траектория движения подвижной частицы при k = 0,3. Общее время счета Рисунок 3 - Сечение структуры с движущейся частицей при k = 0,3 Траектория движения подвижной частицы в случае N=477, k = 0,2 представлена на рисунке (рис. 4). Рисунок 4 - Траектория движения подвижной частицы при k = 0,2. Общее время счета Начальное положение подвижной частицы задавалось координатами (-0,6; -7,3; 0,6). Рисунок 4 - Сечение структуры с движущейся частицей при k = 0,2 В рассмотренных случаях имеется составляющая перемещения подвижной частицы внутри структуры в направлении, противоположном течению жидкости без структуры. Такое поведение частицы кардинальным образом отличается от того, которое имеет место при однородном течении вязкой жидкости через аналогичные структуры. Как установлено в работе [12], в этом случае перемещение подвижной частицы происходит вдоль и перпендикулярно направлению течения жидкости (рис. 5). Рисунок 5 - Траектория движения подвижной частицы при k = 0,3 (однородный поток). Общее время счета Поскольку в одних случаях имеется градиент давления, а в других нет, а скорость однородного поток бралась равной скорости параболического при , то можно предполагать, что такое существенное отличие объясняется только наличием градиента давления. Заключение На основе предложенной ранее модели фильтрации двухфазной жидкости, основанной на представлении пористой структуры, как системы неподвижных частиц, а диспергированной фазы - как системы подвижных частиц, находящихся внутри пористой структуры, рассмотрен случай фильтрации с заданным градиентом давления. Учитывается гидродинамическое взаимодействие всех частиц, как подвижных, так и неподвижных. Компьютерное моделирование рассмотренных течений проводилось с помощью разработанного программного комплекса. Проведен численный расчет течения жидкости в пористой среде на примере 2 модельных структур, образованных, соответственно, из 450 и 478 частиц эффективного размера. Размеры дисперсной частицы, помещенной в вязкую жидкость, составляли 0,3 и 0,2 от размера эффективных частиц. Результаты численного моделирования показывают, что динамика подвижных частиц в потоке с градиентом давления количественно и качественно отличается от динамики частицы в однородном потоке. В частности обнаружено, что перемещение подвижной частицы в градиентном потоке имеет составляющую, противоположную направлению течения жидкости вне структуры. Это означает, что использование градиента давления для вытеснения двухфазной жидкости из пористой среды может приводить не к увеличению, а к уменьшению массопереноса диспергированной фазы. Необходимо более точно учитывать влияния микротечений на перенос каждой из фаз, что позволит более качественно предсказывать массоперенос каждой фазы и определять способы его увеличения.

About the authors

Sergei I. Martynov

Yugra State University

Author for correspondence.
Email: martynovsi@mail.ru

Doctor of Philosophy in Physics and Mathematics, Professor of the Department of Higher Mathematics of the Institute of Technical Systems and Information Technologies 

Russian Federation, 16, Chehova street, Khanty-Mansiysk, 628012

Maria S. Deryabina

Yugra State University

Email: deryabinams@mail.ru

Postgraduate student of the Department of Higher Mathematics of the Institute of Technical Systems and Information Technologies

Russian Federation, 16, Chehova street, Khanty-Mansiysk, 628012

References

Supplementary files