Обзор математических моделей, описывающих процесс транспорта примесей и одиночных частиц в потоке жидкости

ISSN 0868-5886 НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2015, том 25, № 4, c. 36-42

МАТЕМАТИЧЕСКИЕ МЕТОДЫ И МОДЕЛИРОВАНИЕ В ПРИБОРОСТРОЕНИИ

УДК 532.73-1+532.73-2+532.73-3

© В. Е. Курочкин, Б. П. Шарфарец, Е. Б. Шарфарец

ОБЗОР МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ, ОПИСЫВАЮЩИХ ПРОЦЕСС ТРАНСПОРТА ПРИМЕСЕЙ И ОДИНОЧНЫХ ЧАСТИЦ

В ПОТОКЕ ЖИДКОСТИ

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

Кл. сл.: поток жидкости, растворитель, растворенное вещество, дисперсные смеси, транспорт растворенного вещества, транспорт частиц

ВВЕДЕНИЕ

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

ДИСПЕРСНЫЕ СМЕСИ

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

Основой для описания течений в гетерогенных средах является гипотеза взаимопроникающих континуумов, предложенная Х.А. Рахматулиным [1]: многокомпонентная гетерогенная среда представляется совокупностью N сплошных сред,

каждая из которых описывается своей скоростью, плотностью, удельной внутренней энергией, давлением, температурой и т. д. Это положило начало использованию методов механики сплошных сред для описания различного рода смесей с привлечением идеи механики взаимопроникающих (многоскоростных) континуумов и определением взаимопроникающего движения составляющих смеси. Близкие результаты были получены К. Труссдел-лом в 1957 г. [2, 3].

Далее в работах Р.И. Нигматулина, в частности [4-6], эта идея последовательно развивается. Отметим также некоторые другие работы [7-12], посвященные этой тематике, библиография которой весьма обширна и достаточно полно приведена в упомянутых работах.

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

Согласно работам [4, с. 9] и [5, с. 15], наиболее изучены двухфазные дисперсные среды, одна из фаз которых, называемая дисперсной, — капли, пузырьки или твердые частицы. Различают следующие виды дисперсных смесей:

- суспензии — смеси жидкости с твердыми частицами;

- эмульсии — смеси жидкости с каплями другой жидкости;

- газовзвеси — смеси газа с твердыми частицами.

Окружающая несущая фаза (жидкость или газ) называется дисперсионной фазой.

В литературе наряду с динамикой многофазных сред рассматривается и динамика одиночных частиц в потоке жидкости. Изучение динамики одиночной частицы в потоке жидкости важно для

описания дисперснои среды в целом и представляет самостоятельный интерес [10]. Известно [12, с. 5, 6.], что траектория частицы в общем случае не совпадает с траекторией жидкой (лагранжевой) частицы из-за эффекта инерции. Если (массовая) концентрация С частиц достаточно мала (С < 102) и не происходит процессов перехода одной фазы в другую и химических реакций, то воздействие частиц на несущую жидкость осуществляется благодаря силе трения, возникающей из-за наличия проскальзывания, т. е. отличия скорости отдельно взятой частицы от локальной скорости жидкости.

Начнем с описания динамики одиночных частиц.

УРАВНЕНИЕ ДВИЖЕНИЯ ОДИНОЧНОЙ ЧАСТИЦЫ В НЕСЖИМАЕМОЙ ВЯЗКОЙ ЖИДКОСТИ

Теория

Изложенная ниже теория движения одиночной частицы в несжимаемой вязкой жидкости, а также подробный обзор публикаций приведены, например, в работе [12]. Изложим в общем виде строгую математическую постановку задачи о движении одиночной частицы в несжимаемой вязкой жидкости (попутно исправляя опечатки оригинала). Полагаем далее, что в отличие от [12] частица не обязательно является сферической.

Рассматривается частица с границей , плотностью рр , с координатой гр ^), движущаяся со скоростью

dr

f

dt

= V(t).

Пусть и(г, t) и р(г, t) — невозмущенные скорость течения и давление, которые имели бы место в однофазной жидкости, а U(г, t) и р(г, t) — соответствующие поля, возмущенные наличием частицы. Уравнение для скорости частицы имеет вид

dV г

mP 77 = mpg. + $ nj(- р8у + ° j)ds.

dt

Здесь 7,у = 1,2,3, х1 = х, х2 = у , х3 = г; тр —

масса частицы, легко рассчитываемая по ее плотности и объему; gi — 7 -я компонента ускорения

d

силы тяжести; — — производная вдоль траекто-dt

рии частицы; п — вектор внешней нормали к поверхности Бр. Тензор вязких напряжений с?у равен

С7 j = M

dU г + dU,

дх,. дх,.

где и — динамическая вязкость жидкости.

Далее используется гипотеза о несжимаемости жидкости. Поля скорости и давления U(г, t) и р (г, t) удовлетворяют уравнению Навье— Стокса, которое в координатной записи имеет вид

Pf

f dU i - dU Л + U k

dt

\

i = 1,2,3,

dx

k У

= -— + MU, + Pfgi

где р^ = const — плотность жидкости; по повторяющемуся индексу k = 1,2,3 осуществляется суммирование. Поле скорости жидкости удовлетворяет условию несжимаемости

V-U = 0 .

Граничные условия для поля скорости U на поверхности твердой частицы задаются условиями непроскальзывания (прилипания) жидкости

U

= V.

На расстояниях от поверхности частицы, многократно превышающих ее поперечный размер, скорость и давление равны

U

= U.

= р-

где U и р — невозмущенные поля течения при отсутствии частицы, удовлетворяющие системе уравнений

f

Pf

dUi

~дГ

+ U„

dU

\

дх.

k У

f + М&и + Pfgi •

i = 1,2,3,

V-U = 0.

Частица привносит локальное возмущение в течение жидкости благодаря условию прилипания на ее поверхности.

Принимаются следующие представления

и - и = V, р - р = ра.

Здесь поля (и', ра) представляют собой соответствующие возмущения исходных полей без частицы (и, р), вызванные наличием частицы. Тогда

поля и' и ра удовлетворяют следующей системе уравнений:

S=S

р

р

Р/

( ди* ~дГ

+и

еи;

дх,.

+и

еи,

дх.

+и

ди

*\

дх,

к у

1т+"А"' •

дх,.

СТ = ^

Гди дцЛ

V ^ дх,

ст ,, = ^

|ди* , диа, ^ -- +-=

дх,- дх,.

V 1 ,

§ «1(-+ (Г, ^ = П-д^ + ^Аи,

Л

дх

dV +

У

II

Л

дх

у

§ «1 (-р5, + 1 = °р (+ )|

Здесь гр — геометрический центр частицы. Из (3а) и (6) находим

др ди „ ди,.

Р =-~х + = Рf

р*, =-др_ + цАил { = дх

д*

-+и

дх,

к У

ди „ ди, * ди, ди, +и—-+и,* ——-+и,

* л

д*

дх

дх

дх

к у

Подстановка (9) в (2) с учетом последних равенств приводит к итоговому уравнению

V-Ц* = 0, Ц* = о, р*\ = 0.

Далее вычисляется поверхностная сила § «1 (-р8,, + СТ, )С.у в (2), действующая на частицу:

СК I ди, тт ди, Р?—Т7 = Р/1 — + и

dt

д*

дх,

к У

+

+ & ( Рр - Р/ ) +

§§ «1 (-Р$р + ст 1 = § «1 (-рду + СТу ^ +

+ § «1

где соответствующие тензоры вязких напряжений равны

После применения теоремы Гаусса— Остроградского к поверхностным интегралам в (8) получаем

для решения которого необходимо нахождение полей (и,р) и (и*, р*) из соответственно краевых задач (3а), (4а) и (6), (7).

Одним из случаев, когда решения граничной задачи для поля возмущения и * и р*

и выражение для силы F* (11) могут быть найдены аналитически, является случай, когда движение частицы происходит в стоксовом режиме, а это означает здесь, что наряду с условием малости диаметра частицы (* / L  1) должны выполняться следующие условия [13]:

- число Рейнольдса частицы в относительном

потоке мало

V, - и0

Яв = -°1

V

< 1,

- отношение е квадрата диаметра частицы * к квадрату характерного масштаба пограничного

vL

слоя в окрестностях частицы — мало

и

где О — объем частицы. Далее, используя предположение о том, что размер частицы много меньше характерного пространственного размера несущего течения L (для сферического включения диаметром * это условие имеет вид * / L  1), из (8а) получено

е =

ио *2 , —— << 1, IV

где ио и V, — характерные масштабы скорости жидкости и частицы соответственно. При этом

и01 и0*

числа Рейнольдса несущего течения - и -

V V

могут быть немалыми.

Решение задачи (6), (7) получено в работе [13]. Выражение для F * имеет вид

р * =-И£- ^ (г (*), *) +—А^

2 dt I 1 ^ ' ' 60

7* _ _/ а

2 с 18р

/

/

vQ (г (*), * )-

I С*' ^

9р

*

д.

С*

* \п Х

С*'

*

*

г=г

р

г=г

р

р

р

р

р

р

р

р

X

где

d2

Q

Слагаемые в правой части (15) соответствуют в порядке следования силам присоединенной массы, вязкого трения Стокса и силе Бассе (о силе Бассе см. в [14, с. 131, 132]).

Для анализа решения уравнения движения частицы (12) его удобно переписать в безразмерных переменных, например, нормированных на масштаб скорости и времени ио и Т = L / ио, в виде

„dV DU , ч , s d , ттч

s IF = s D+(<s -')g- 2*(V - U )-

- 9 (V - U)-x

i dp

dt'

Vt -11

где

s = Pp / Pf; g'=sgd2/(vUо);

DU=dU+( u-v) u

Dt dt v '

или в тензорной записи

DU dU

Dt

= —L+U dt k

dU. _i

dx.

(x1,x2,x3) = (x,y,z) .

поправки, пропорциональные

0 L

в формуле (16) (поправки Факсена), пренебрежимо малыми. Время релаксации частицы в стоксовом режиме определяется так:

d 2S

т =-

р 18v

смещение относительно линий тока несущей жидкости. Из (17) также следует, что характерное время, за которое происходит это смещение, оценивается как Т /11 - ^ .

Если произведение е5 ~ 1, то сила инерции частицы сравнима с силой вязкого трения и тяжести. В этом случае движение частицы определяется ее инерцией и силами вязкого трения и тяжести.

Рассмотрение частного случая

В работе [16, с. 116, 117] рассматривается одномерное движение сферической частицы радиусом а в следующем виде:

4 3 dV

Pp — = -6WaV + Fext •

3 ' p dt

Здесь V (t) — скорость частицы; -6n^aV — сила вязкого сопротивления Стокса при движении частицы со скоростью V; Fixt — некоторая внешняя сила. В предположении того, что Fext = const и V| = 0, в [16] получено следующее решение:

V (t ) = -Fext-6л^а

F„

1 - exp

9m

2Ppa

f

6лца

v ' p 1 - exp I —

//

Далее следуют некоторые результаты анализа уравнения (17), заимствованные из [15].

Диаметр частицы предполагается достаточно малым (см. условие (13)), что позволяет считать

й2

й2ли ~ ип%г

Конкретизируем уравнение (19). Примем допущение о том, что частица находится в стационарном плоском потоке жидкости, двигающемся с постоянной скоростью и . Тогда силу Fext можно принять равной соответствующей силе Стокса,

воздействующей течением на частицу ^ = 6Щаи,

или после подстановки этого выражения в (19) получаем

(

V (t) = U 1 - exp I —

t

\

Из анализа (17) следует, что при 5 = 0(1) сила вязкого трения (сила Стокса) преобладает над инерционными силами (связанными с градиентом давления и эффектом присоединенной массы). Это означает, что скорость частицы быстро за время релаксации хр << Т подстраивается под скорость

окружающей жидкости. В общем случае плотности частицы и несущей жидкости не равны рр Ф ру и инерция частицы обусловливает ее

В работе [16] показано, что фигурирующее

2p а

в решении (20) время релаксации т = -

9м

харак-

теризующее стабилизацию процесса при рр, примерно совпадающей с плотностью воды ру, и

а ~ 5 мкм равно х = 5 мкс . Тогда из (20а) следует, что через время Т ~ 2х = 10 мкс частица практически достигает стационарной скорости течения V = и.

x

t

Как видно из изложенного, уравнение (19) аппроксимирует общее уравнение (12) в стоксовом режиме при б = рр / pf = 0(1) и малом радиусе

частицы, когда сила вязкого трения (сила Стокса) преобладает над инерционными силами (первый член справа в (12)), а также, когда в (15) можно пренебречь эффектом присоединенной массы и силой Бассе, а в (16) поправкой Факсена. Кроме того, в (19) не учитывается поле силы тяжести, что, впрочем, легко устраняется.

Таким образом, при выполнении соответствующих условий, сложное уравнение для траектории частицы (12) вполне аппроксимируется упрощенным уравнением (19).

МЕХАНИКА МНОГОСКОРОСТНЫХ КОНТИНУУМОВ

Согласно [1-9 и др.], описание методами механики сплошной среды различного рода смесей, и гомогенных, и гетерогенных, связано с введением понятия многоскоростного континуума и определением взаимопроникающего движения составляющих смеси.

Согласно, например, [4], многоскоростной континуум представляет собой т континуумов, каждый из которых относится к своей составляющей (фазе или компоненте) смеси и заполняет один и тот же объем, занятый смесью. Для каждого из этих составляющих континуумов в каждой точке определяется обычным образом плотность (приведенная) pi (масса , -составляющей в единице объема среды), скорость V 1 (,' = 1,...,т ), а затем и другие параметры, относящиеся к своему континууму и своей составляющей смеси. Таким образом, в каждой точке объема, занятого смесью, будет определено т плотностей pi, т скоростей V { и т. д.

Вводятся также параметры, характеризующие смесь целиком, а именно плотность смеси р и среднемассовая скорость смеси V :

др т

-%: + V•{ptvt) = £1 , = 1,...,т .

от 1=1

Здесь JJi характеризует интенсивность перехода

массы из 1 -й в , -ю составляющую (или наоборот

из , -й в 1 -ю; тогда JJi < 0). С учетом

закон сохранения массы (23) записывается так [4, с. 15]

Р-

Если просуммировать (23) или (24), учитывая (22), то получается закон сохранения массы смеси в целом, имеющий обычный вид, как в односкоро-стном случае [4, с. 16]

дР + V•( pv ) = 0.

Закон сохранения импульса в дифференциальной форме для , -й фазы имеет вид [4, с. 16]

р |'дг+(V,-V) ^ 1=

т

V•стf + РА + £(Р -^V,).

1=1

Р=£ Р , РЛ? = £ Рг^'г .

Здесь — интенсивность обмена импульсом

между 1 -й и -й фазами; стк1 — тензор поверхностных сил (напряжения) , -й фазы; g — вектор ускорения силы тяжести.

т

Если принять ст = £ ст; , то после суммирова-

¿=1

ния (26) по т получаем уравнение импульсов среды в целом [4, с. 17]

дv

Р|¥ + (^ •V)^ " = = -V • стк + pg -£ V•( р'Я.'М,).

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

^ = V,-V, (£= 0).

Механика смесей строится на основе физических законов сохранения массы, импульса и энергии. Дифференциальный закон сохранения массы каждой составляющей имеет вид [4, с. 15]

Здесь ^^ — аффинор и операция вычисления дивергенции, определяется так:

V•(РIWIWг ^•(Р,^ ) ^ .

Уравнение импульсов среды в целом уже в отличие от уравнения неразрывности среды в целом зависит от относительного движения составляющих.

Закон сохранения энергии -й фазы сводится к выражению [4, с. 18]

f

dt

u. + — ' 2

2Л

= V( С - Ч, ) + Рг g •V +

/

+z

j=1

F - J

^ji Jji

f v]^

u + — 2

V

Здесь Ец — интенсивность обмена энергией между 7 -й и у -й фазами; ui — внутренняя энергия единицы массы 7 -й фазы; vi = ; ci — вектор плотности поверхностных сил (напряжений); qI — вектор плотности потока тепла в i -й фазе. После суммирования в (28) по всем фазам получается уравнение энергии смеси в целом [4, с. 18]:

Н Т7 т

Р^Е + ) =

Ш ,=1

= V • (c - q) + pg • v + g • £ (W,Pi),

где

v2 1 Ä ptWt2 F = u + — + — > г ' 2 p 1=1 2

c = Z С,

,=i

q = £ q,

w = W

Как видно, решение задачи определения взаимопроникающего движения составляющих смеси методами механики сплошной среды с введением понятия многоскоростного континуума представляет собой очень трудоемкую задачу решения системы связанных (в том числе нелинейных) уравнений в частных производных (23)-(25), (27) и (29). Ситуация существенно упрощается для так называемых бесконечно разбавленных растворов [17]. В этом случае система указанных уравнений сводится к уравнениям неразрывности (25) и На-вье—Стокса (27) для растворителя, уравнение сохранения энергии (29) трансформируется в уравнение теплопереноса, а закон сохранения массы примеси (23) трансформируется в уравнение мас-сопереноса. Обо всем этом подробно написано в работе [18].

Отметим, что теория разбавленных растворов обладает многими аспектами, которые лишь слегка видоизменяются в более общей теории концентрированных растворов [17, с. 276].

СПИСОК ЛИТЕРАТУРЫ

1. Рахматулин Х.А. Основы газодинамики взаимопроникающих движений сжимаемых сред // ПММ. 1956. Т. 20, № 2. С. 184-195.

2. Truesdell C. The Rational Mechanics of Materials (ed. by C. Truesdell). New York, London: Gordon&Breach Science Publisher, 1965.

3. Truesdell C., Toupin R.A. Handbuch der Physik. Bd III/1 (ed. S. Flügge). Berlin, Heidelberg: Springer, 1960. S. 226.

4. Нигматулин Р.И. Основы механики гетерогенных сред. М.: Наука, 1978. 336 с.

5. Нигматулин Р.И. Динамика многофазных сред. Ч. I. М.: Наука. 1987. 464 с.

6. Нигматулин Р.И. Динамика многофазных сред. Ч. II. М.: Наука, 1987. 360 с.

7. Пилюгин Н.Н., Тирский Г.А. Динамика ионизированного излучающего газа. М.: Изд-во МГУ, 1989. 311 с.

8. Manninen M., Taivassalo V. and Kallio S. On the mixture model for multiphase flow. Technical Research Center of Finland (VTT), 1996. 67 p.

9. Crowe C., Sommerfeld M. and Tsuji Y. Multiphase flows with droplets and particles. CRC Press, 1998. 471 p.

10. Шрайбер А.А., Гавин Л.Б., Наумов В.А., Яценко В.П. Турбулентные течения газовзвеси. М.: Наука. 1987. 239 с.

11. Soo S.L. Fluid dynamics of multiphase systems. Waltham, MA, Blaisdell Pub. Co. 1967. 524 p.

12. Дружинин О.А. Исследование динамики вихревых потоков и волн в дисперсных и стратифицированных средах. Дис.....докт. физ.-мат. наук. Нижний Новгород: ИПФ РАН, 2004. 300 с.

13. Maxey M.R., Riley J.J. Equation of motion for small rigid sphere in a non-uniform flow // Phys. Fluids. 1983. № 26. P. 883-889.

14. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т. 6. Гидродинамика. М.: Наука, 1988. 736 с.

15. Maxey M.R. On the advection of spherical particles in a non-uniform flow // Phil. Trans. R. Soc. of London A 333. 1990, no. 1631. P. 289-307. doi: 10.1063/1.864230.

16. Bruus H. Theoretical Microfluidics. Oxford University Press., 2008. 346 p.

17. Ньюмен Дж. Электрохимические системы. М.: Мир, 1977. 464 с.

18. Шарфарец Б.П. Обзор теории явлений переноса и поверхностных явлений применительно к решению некоторых задач научного приборостроения // Научное приборостроение. 2015. Т. 25. № 3. C. 45-64.

Институт аналитического приборостроения РАН, г. Санкт-Петербург)

ВЫВОДЫ

Таким образом, в работе рассмотрены поочередно математические модели транспорта одиночных частиц и модели транспорта дисперсных смесей в потоке жидкости.

Контакты: Шарфарец Борис Пинкутович, sharb@mail.ru

Материал поступил в редакцию: 7.07.2015

d

ISSN 0868-5886

NAUCHNOE PRIBOROSTROENIE, 2015, Vol. 25, No. 4, pp. 36-42

A REVIEW OF MATHEMATICAL MODELS DESCRIBING THE TRANSPORT OF SOLUTES AND SINGLE PARTICLES

IN FLUID FLOW

V. E. Kurochkin, B. P. Sharfarets, E. B. Sharfarets

Institute for Analytical Instrumentation of RAS, Saint-Petersburg, Russia

Discusses some mathematical models describing the transport of dispersed mixture and single particles in fluid flow. In relation to the movement of single particles under certain conditions, when the shift of the trajectory relative to the main fluid flow can be neglected. Noted the relative simplicity and usefulness of the of mass transfer theory of solutes in dilute solutions.

Keywords: fluid flow, solvent, solute, dispersed mixture, transport of solute, particle transport

REFERENCES

1. Rachmatulin Ch.A. [Fundamentals of gas dynamics of the interpenetrating movements of the squeezed environments]. Prikladnaya matematika i mechanika [Applied mathematics and mechanics], 1956, vol. 20, no. 2, pp. 184-195. (In Russ.).

2. Truesdell C. The Rational Mechanics of Materials (ed. by C. Truesdell), 1965, New York, London, Gordon&Breach Science Publisher.

3. Truesdell C., Toupin R.A. Handbuch der Physik. Bd 111/1, S. 226 (ed. S. Flügge), 1960, Berlin, Heidelberg, Springer.

4. Nigmatulin R.I. Osnovy mechaniki geterogennych sred [Bases of mechanics of heterogeneous environments], Moscow, Nauka Publ., 1978. 336 p. (In Russ.).

5. Nigmatulin R.I. Dinamika mnogofaznych sred [Dynamics of multiphase environments], Part I. Moscow, Nauka Publ., 1987. 464 p. (In Russ.).

6. Nigmatulin R.I. Dinamika mnogofaznych sred [Dynamics of multiphase environments], Part II. Moscow, Nauka Publ., 1987. 360 p. (In Russ.).

7. Pilyugin N.N., Tirskiy G.A. Dinamika ionizirovannogo iz-luchayuschego gaza [Dynamics of the ionized radiating gas], Moscow, MGU Pupl., 1989. 311 p. (In Russ.).

8. Manninen M., Taivassalo V. and Kallio S. On the mixture model for multiphase flow. Technical Research Center of Finland (VTT), 1996. 67 p.

9. Crowe C., Sommerfeld M. and Tsuji Y. Multiphase flows with droplets and particles. CRC Press, 1998. 471 p.

10. Shrayber A.A., Gavin L.B., Naumov V.A., Yazenko V.P. Turbulentnye techeniya gazovzvesi [Turbulent flows of a

Contacts: Sharfarets Boris Pinkusovich, sharb@mail.ru

gas-suspension]. Moscow, Nauka Publ., 1987. 239 p. (In Russ.).

11. Soo S.L. Fluid dynamics of multiphase systems. Waltham, MA, Blaisdell Pub. Co., 1967. 524 p.

12. Druzhinin O.A. Issledovanie dinamiki vichrevych potokov i voln v dispersnych i stratifizirovannych sredach. Dokt. Diss. [Research of dynamics of vortex streams and waves in the disperse and stratified environments. Dokt. Diss.]. Nizhniy Novgorod, IPF RAN, 2004. 300 p. (In Russ.).

13. Maxey M.R., Riley J.J. Equation of motion for small rigid sphere in a non-uniform flow. Phys. Fluids., 1983, no. 26, pp. 883-889.

14. Landau L.D., Lifshiz E.M. Teoreticheskaya fizika. T. 6. Gidrodinamika [Theoretical physics. Vol. 6. Hydrodynamics], Moscow, NaukaPubl., 1988. 736 p. (In Russ.).

15. Maxey M.R. On the advection of spherical particles in a non-uniform flow. Phil. Trans. R. Soc. of London, A 333, 1990, no. 1631, pp. 289-307. doi: 10.1063/1.864230.

16. Bruus H. Theoretical Microfluidics. Oxford University Press., 2008. 346 p.

17. N'yumen Dzh. Elektrochimicheskie sistemy [Electrochemical systems]. Moscow, Mir Publ., 1977. 464 p. (In Russ.).

18. Sharfarets B.P. [An overview of the theory of transport phenomena and surface phenomena in relation to the solution of some problems of analytical instrumentation]. Nauchnoe priborostroenie [Scientific Instrumentation], 2015, vol. 25, no. 3, pp. 45-64. doi: l0.l8358/np-25-3-i4564 .

Article received in edition: 7.07.2015