Механика жидкости и газа Вестник Нижегородского университета им. Н.И. Лобачевского, 2011, № 4 (3), с. 740-741
УДК 532.529
ВЗАИМОДЕЙСТВИЕ ФАЗ И ПРИСОЕДИНЕННАЯ МАССА ДИСПЕРСНЫХ ЧАСТИЦ В ПОТЕНЦИАЛЬНЫХ ПОТОКАХ ЖИДКОСТИ
© 2011 г. О.Б. Гуськов, Б.В. Бошенятов
Институт прикладной механики РАН, Москва
Поступила в редакцию 16.06.2011
С использованием метода самосогласованных полей для потенциальных течений несжимаемой жидкости рассмотрено движение монодисперсной двухфазной среды с учетом гидродинамического взаимодействия сферических дисперсных частиц. Для произвольного соотношения плотностей частиц и жидкости получены усредненные аналитические зависимости для относительной скорости и присоединенной массы частиц от их объемной концентрации.
Ключевые слова: гидродинамическое взаимодействие, дисперсная частица, присоединенная масса, самосогласованное поле, потенциальное течение.
При изучении движения гетерогенных сред возникает проблема описания внутрифазных и межфазных взаимодействий, вследствие чего для получения решаемых уравнений приходится использовать рациональные схематизации, упрощения и эмпирические зависимости. Так, например, инерциальные силы присоединенных масс часто определяют из потенциального обтекания одиночной сферы несущей (гомогенной) жидкостью. Ясно, что такое приближение справедливо лишь для достаточно малых концентраций дисперсных частиц. При более высоких концентрациях необходимо учитывать влияние дисперсных частиц друг на друга и их гидродинамическое взаимодействие.
В настоящем исследовании в рамках классических уравнений гидродинамики для потенциальных течений несжимаемой жидкости дано детальное описание мелкомасштабного движения двухфазной монодисперсной среды с учетом гидродинамического взаимодействия сферических дисперсных частиц. Затем путем усреднения по ансамблю находятся усредненные зависимости скорости и присоединенной массы частиц в макроскопическом масштабе. Задача решена при следующих упрощающих предположениях: внутренние движения частиц (деформация, вращение) отсутствуют; частицы распределены в пространстве статистически равномерно, кроме того, отсутствуют процессы дробления или слияния частиц, т.е. количество частиц неизменно. Поскольку действие тяготения (§) состоит просто в том, что к инерци-альным силам добавляется постоянная гидростатическая сила выталкивания, достаточно рассмат-
ривать движение с нулевой потенциальной энергией, т.е. при g = 0.
Пусть в начальный момент времени суспензия, включающая произвольное число N сферических частиц, имеющих радиус а, неподвижна и занимает полубесконечное пространство, ограниченное плоской бесконечной стенкой Ж; центры частиц расположены в точках х^) (г = 1, 2, ..., N у = 1, 2, 3), а их скорости и^. Определим среднюю скорость частиц вдали от поршня при его импульсном ускорении от нуля до скорости Ц;0). Задача решается в системе координат, связанной с подвижной стенкой Ж, ось Ох3 направлена перпендикулярно плоскости стенки. Безразмерный потенциал скоростей, удовлетворяющий граничным условиям на поверхностях частиц и стенки, ф = = ф*/Ш(0), где Ь — расстояние между центрами частиц; и(0) — скорость движения стенки в направлении Ох3 ищется в виде суммы возмущений ф(,), вносимых в заданный поток ф(0) всеми частицами, а также соответствующих им отражений от стенки ф^? [1]:
N
ф=ф(0)+Е (ф(г)+ф^).
г=1
Потенциалы возмущений от частиц ищутся во взаимосогласованном виде, что в результате приводит к определенному виду общего решения, а также к системе уравнений, связывающей заранее неизвестные тензорные коэффициенты в полученном представлении решения для потенциала ф. Предположение о малости параметра а = а/Ь позволяет упростить систему для тензорных ко -эффициентов и свести ее к рекуррентной системе
алгебраических уравнений, где каждое последующее приближение по малому параметру а для искомых функций определяется через предыдущие. Таким образом, зная нулевое приближение, можно определить все функции в любом приближении по а. Приведенные ниже аналитические уравнения получены с учетом членов разложения до степени а14 включительно.
При решении задач, имеющих большое количество частиц, целесообразно от дискретного (детального) описания движения перейти к осреднен-ному (макроскопическому) описанию дисперсных систем, используя процедуру усреднения по ансамблю [2]. Применение данной процедуры для скорости частиц u вдали от стенки Жс точностью до первой степени по объемной концентрации с дисперсной фазы приводит к следующему результату:
3и(0)
и = ,. ,.. [1 + к(Х)с],
где
(1 + 2А) 151 (1 -А)
к(А) =---------- -----^- +
72 (1+ 2А)
+
1415 (1 -А)2
5 (1 -А)3
5632 (1 + 2А)2 2112(1 + 2А)
3
личается лишь на 0.2% от соответствующего результата, полученного в работе [3].
Зная (1) и уравнение движения для частицы [3], находим коэффициент присоединенной массы дисперсной частицы:
т 1
ц(А) = —= - [1 + кт (А)с], рт 2
(2)
(1)
При X = р0/р = 1 плотность частиц р0 равна плотности несущей жидкости. В этом случае, как и следовало ожидать, и = и(0). При X ^ ^ скорость частиц и ^ 0 . Предел X ^ 0 соответствует случаю газовых пузырьков. В этом случае и = 3и(0)(1 — 1.848с), что совпадает с [1] и от-
где
к (А) = 111 _ 4245 (1 -А) + 40 (1 -А)2
т 48 11264 (1 + 2А) 11264 (1 + 2А)2 ’
т* - присоединенная масса частицы, т - объем частицы.
Аналитическая зависимость кт(А) в формуле (2) практически совпадает с данными численных расчетов [4], а при А = 0 - с результатом [3] для пузырьков газа.
Чтобы получить формулы, аналогичные (1) и (2) для более концентрированных суспензий, необходимо учесть следующие члены по степеням концентрации. Предварительный анализ показывает, что при определении коэффициента перед с2 принципиальных трудностей не возникает.
Список литературы
1. Струминский В.В, Гуськов О.Б., Корольков Г.А. // ДАН СССР. 1986. Т. 290, №4. С. 820-824.
2. Beenakker C.W.J., Mazur P. // Phys. Fluids. 1985. Vol. 28. P 767-769.
3. Wijngaarden L.V. // J. Fluid Mech. 1976. Vol. 77. Part 1. P. 27-44.
4. Felderhof B.U. // J. Fluid Mech. 1991. Vol. 225. P. 177-196.
INTERACTION OF PHASES AND VIRTUAL MASS OF DISPERSED PARTICLES IN POTENTIAL FLOWS OF FLUID
O.B. Guskov, B. V Boshenyatov
The motion of a monodisperse medium is considered on the bases of self-consistent field method, taking into account the hydrodynamic interaction of spherical dispersed particles. The averaged analytical expressions for the relative velocity and virtual mass of particles are derived as a function of volume concentration and density ratio of the particles and the fluid.
Keywords: hydrodynamic interaction, dispersed particle, virtual mass, self-consistent field, potential flow.