Прохождение жидкости со случайной скоростью через среду, заполненную частицами-гранулами

НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ, Н. Э. БАУМАНА

НАУКА и ОБРАЗОВАНИЕ

Эл № ФС77 - 48211. ISSN 1994-0408

Прохождение жидкости со случайной скоростью через среду, заполненную частицами-гранулами # 08, август 2014

Б01: 10.7463/0814.0724713 Морозов А. Н.1, Скрипкин А. В.1а

УДК 532.542

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

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

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

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

Предлагаемая работа посвящена исследованию потока жидкости с флуктуирующей составляющей скорости через цилиндрическую трубу заполненную частицами-гранулами. Найдены статистические характеристики скорости жидкости, если первоначально они соответствовали характеристикам белого шума. Заметим, что вопросы распространения жидких и газообразных веществ в средах, содержащих крупные частицы, а также в пористых средах, активно обсуждаются в научной литературе последние несколько десятков лет. В частности, в работе [7] исследуется течение жидкости в среде с порами, которые рассматриваются как ловушки; статья [8] рассматривает частицы-гранулы как соответствующую броуновскую решетку. При этом происходит взаимодействие частиц движущейся жидкости и таких броуновских частиц. Заметим также, что в последнее время вместо частиц-гранул активно изучаются другие типы объектов, находящихся в жидкой среде, например, волокна [9].

Рассмотрим прямую трубку постоянного кругового сечения 5, заполненную жидкостью, в общем случае с непостоянной плотностью, и направим ось Ох вдоль оси трубки (рис. 1).

Рис 1. К постановке задачи

Будем считать, что при х > 0 в ней содержатся частицы-гранулы, создающие раствор суспензии, причем в общем случае концентрация частиц-гранул будет зависеть от расстояния вдоль трубки: N = ^х) (зависимостью от времени пренебрегаем). Величиной, способной охарактеризовать степень пористости, может являться отношение

Ф) = ^, (1)

где (х) - площадь сечения, не занятая гранулами, 5 - общая площадь сечения трубки. В случае гомогенной среды величину в(х) можно считать постоянной.

В работе [10], исходя из уравнения непрерывности, показано, что в указанных условиях диффузионное уравнение для концентрации частиц движущейся в трубке жидкости имеет вид

^х)*Ы= о , (2)

д дх2

где В - коэффициент диффузии (который, вообще говоря, тоже зависит от расстояния вдоль трубки, но мы этим фактом пренебрегаем).

Начальные и граничные условия для уравнения (2) имеют вид:

«ь(0, (3)

и(х,0) = 0. (4)

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

Дополним выражения (2) - (4) уравнением для потока жидкости через сечение х = 0, при этом учтем тот факт, что течение жидкости всегда сопровождается случайной составляющей, вызванной, например, изменениями давления или температуры, а также плотности жидкости. Поэтому для указанного потока числа частиц q (?) запишем соотношение

дп(х, /)

q(t ) = -D x=0

+ (t).

Здесь ^ (^) - указанный случайный поток, который мы будем считать белым шумом с ин

тенсивностью v.

Если распределение гранул в среде не зависит от расстояния вдоль оси трубки ( s(x) = s = const), то поставленная задача (2) - (5) может быть полностью решена. Решение дифференциального уравнения (2) с граничным и начальным условием (3) и (4) имеет вид [11]

n( x, t) =

nD f (t-T)3

n0(x)dT, x > 0.

Продифференцировав последнее выражение по координате x, получим уравнение

для частной производном плотности жидкости

dn( x, t)

откуда для потока (5) находим

) = -,IE f.1 ^t + 5, (t),

ns^ л/t-t dx

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

Заметим, что формулы (6) и (7) соответствуют обычному диффузионному процессу в полупространстве с коэффициентом диффузии О/в, что связано с наличием частиц-гранул в среде.

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

частиц жидкости определяется равенством^ («)=|~(«)2, где ) - преобразование Лапласа по времени функции q (?), р - параметр преобразования, принимая во внимание также, что спектральная плотность мощности белого шума равна его интенсивности, т.е.

|~00®)2 =а, для спектральной плотности мощности флуктуаций потока жидкости получим

Gq(«) = Оа®. (9)

Таким образом, в этом случае флуктуации потока будут относиться к синему шуму, являющемуся производной от фликкер-шума.

Заметим, что аналогичные выражениям (6) и (7) соотношения были получены в задачах теплопроводности и диффузии вблизи плоской поверхности в другой работе [12]

Спектральная плотность мощности флуктуаций изменения концентрации жидкости (связанная с изменением плотности) определяется с помощью метода, используемого выше для получения функции (9), и у нулевого сечения имеет вид

оп(«)=—,--.. (10)

« «+л--ъ о / в£

«I V в^ )

Эта формула также аналогична полученным в работе [11] выражениям, относящимся к другим задачам. Эта аналогия может быть обоснована тем, что если касательную силу,

действующую со стороны жидкости на единицу площади поверхности, ¥ = Миз [12]

заменить на соотношение Q = £—, где Q - общее количество частиц жидкости, прошедшее за время Л через сечение трубки, то формулы переходят одна в другую. Из формулы (10) легко видеть, что при малых площадях сечения трубки спектральная плотность Сп («) принимает характер фликкер-шума. Отметим также, что спектральная плотность (10) сильно зависит от отношения О/ £, ввиду чего при плавном увеличении площади поперечного сечения трубки мощность флуктуаций количества проходящих частиц жидкости может проявлять характер некоторых пульсаций (это, в частности, видно из рис. 2, на котором изображены графики функции (10) при разных 5).

Рис. 2. Графики спектральной плотности флуктуаций числа проходящих по трубке частиц при S = 1,5 мм

(1), S = 1,6 мм 2 (2), S = 1,75мм2 (3)

Заметим, наконец, что если функция s(x) крайне медленно зависит от координаты х: s(x) = a+bx, |b| << 1, то все полученные выше формулы верны в первом приближении и в

этом случае. В них следует заменить величину £ на записанные выражения для s(x). Разумеется, все это справедливо при не очень больших х.

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

Список литературы

SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU

SCIENCE and EDUCATION

EL № FS77 - 48211. ISSN 1994-0408

Passing of Liquid from the Random Velocity Through Medium Filled by Particles-Granules # 08, August 2014

DOI: 10.7463/0814.0724713 A.N. Morozov1, A.V. Skripkin1,a

askripkin@bmstu.ru

The article discusses in passing the particles generally inhomogeneous density fluid through the tube in which the particles are granules, the dimensions of which are much exceed size liquid particles. It is assumed that the liquid is exposed to external interference, which leads to the presence of the velocity fluctuations of liquid particles. The aim of this work is to statistical characteristics of liquid density in these conditions, including the power spectral densities of the density fluctuations, and to identify the statistical characteristics considered depending on the cross-sectional area of the tube. For this problem, we obtain an equation similar to diffusion. His solution for additional initial and boundary conditions led to a stochastic integral equation for the density of the liquid. Using the method of description of stochastic systems described by linear integral transforms it possible to obtain the necessary characteristics of the process. It turned out that the low-frequency fluctuations in the density of the liquid have the character of flicker noise (pink noise). Moreover, due to the strong dependence of the behavior of the fluid and the diffu-sivity ratio of the cross sectional area of the tube, the power spectral density of the density fluctuations are not dependent on the cross-sectional area continuously, but undergoes a kind of "oscillation". It is also shown that if the distribution of particle beads in an medium is not homogeneous and is weakly dependent on the distance along the tube axis, it is still possible to use the first approximation results for the homogeneous distribution of the particles of the granules. To do this, the formulas obtained granules distribution parameter should be replaced with the corresponding function. The results obtained can be used in problems of fluid flow in pipes and suspensions, as well as issues of hydraulic engineering.

References