Научная статья на тему 'ЭФФЕКТИВНОСТЬ УЛАВЛИВАНИЯ ИНЕРЦИОННЫХ АЭРОЗОЛЬНЫХ ЧАСТИЦ В ПЕРИОДИЧЕСКОЙ ЯЧЕЙКЕ С ПОРИСТЫМ ЦИЛИНДРОМ'

ЭФФЕКТИВНОСТЬ УЛАВЛИВАНИЯ ИНЕРЦИОННЫХ АЭРОЗОЛЬНЫХ ЧАСТИЦ В ПЕРИОДИЧЕСКОЙ ЯЧЕЙКЕ С ПОРИСТЫМ ЦИЛИНДРОМ Текст научной статьи по специальности «Физика»

CC BY
19
3
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ВЗВЕШЕННЫЕ ЧАСТИЦЫ / ПОРИСТЫЙ ЦИЛИНДР / МОДЕЛЬ СТОКСА - БРИНКМАНА / ЭФФЕКТИВНОСТЬ УЛАВЛИВАНИЯ ЧАСТИЦ

Аннотация научной статьи по физике, автор научной работы — Хазиев Айрат Ринатович, Зарипов Шамиль Хузеевич, Марданов Ренат Фаритович, Пилюгин Александр Геннадьевич

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

EFFICIENCY OF CAPTURE OF INERTIA AEROSOL PARTICLES IN A PERIODIC CELL WITH A POROUS CYLINDER

In this study, aerosol flow in circular and rectangular periodic cells of a regular row of porous cylinders is simulated. The hydrodynamic field of the carrier medium flow outside and in the domain of the porous cylinder is described in the approximation of the Stokes-Brinkman model for an incompressible gas. For the model of a circular periodic Kuwabara cell, an analytical solution is obtained; in the case of a rectangular cell, the boundary value problem for the flow stream function is solved numerically using the boundary element method. In the determined velocity fields of the carrier medium, the Lagrangian equations of motion of suspended particles are integrated, and the dependence of the particle capture efficiency as a result of the inertial impact and the interception mechanism on the Stokes number at different cell density and the Darcy number of a porous cylinder is calculated. The results of the calculations are found to be in a good agreement with those from other works carried out using the CFD package. Parametric studies of the efficiency of capture of aerosol particles as a function of the Stokes number for various values of the Darcy number of a porous cylinder, the periodic cell density, and the interception parameter are performed. It is shown that when a porous cylinder is flown around, the particle capture efficiency tends to a non-zero value that increases depending on how high the Darcy number and the density of the periodic cell are. In the region of small Stokes numbers, the capture of particles is provided by both the interception mechanism and the flow through the porous cylinder. A formula is derived for estimating the efficiency of particle trapping under the effect of particle interception on a porous cylinder.

Текст научной работы на тему «ЭФФЕКТИВНОСТЬ УЛАВЛИВАНИЯ ИНЕРЦИОННЫХ АЭРОЗОЛЬНЫХ ЧАСТИЦ В ПЕРИОДИЧЕСКОЙ ЯЧЕЙКЕ С ПОРИСТЫМ ЦИЛИНДРОМ»

УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО УНИВЕРСИТЕТА. СЕРИЯ ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ

ISSN 2541-7746 (Print) ISSN 2500-2198 (Online)

2022, Т. 164, кн. 2-3 С. 181-193

ОРИГИНАЛЬНАЯ СТАТЬЯ

УДК 532.529.6 <!ок 10.26907/2541-7746.2022.2-3.181-193

ЭФФЕКТИВНОСТЬ УЛАВЛИВАНИЯ ИНЕРЦИОННЫХ АЭРОЗОЛЬНЫХ ЧАСТИЦ В ПЕРИОДИЧЕСКОЙ ЯЧЕЙКЕ С ПОРИСТЫМ

ЦИЛИНДРОМ

А.Р. Хазиев, Ш.Х. Зарипов, Р.Ф. Марданов, А.Г. Пилюгин

Казанский (Приволжский) федеральный университет, г. Казань, 420008, Россия

Аннотация

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

Ключевые слова: взвешенные частицы, пористый цилиндр, модель Стокса-Бринк-мана, эффективность улавливания частиц

Введение

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

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

Рис. 1. Круговая периодическая ячейка с пористым цилиндром

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

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

1. Аналитическое решение для круговой ячейки Кувабара

Рассмотрим течение в периодической круговой ячейке с пористым цилиндром (рис. 1). Задача формулируется в безразмерном виде. В качестве масштаба длины и скорости выбраны радиус пористого цилиндра Яс и скорость невозмущенного потока и. Расчетная область П = Пе У П® состоит из однородной области Пе внешнего течения радиуса г = Н и пористой среды П® внутри цилиндра радиуса г = 1. Пусть фе, ф® - функция тока течения вне и внутри цилиндра. Здесь и далее индексы " е" и " г" соответствуют внешней и внутренней областям течения.

Течение во внешней области Пе описывается в приближении течения Стокса. В этом случае функция тока фе(х,у) внешнего течения будет удовлетворять би-гармоническому уравнению [10]

Д2фе = 0 (1)

и условиям на внешней границе круговой ячейки Ге и линии раздела Г® между однородной и пористой областями. Для круговой ячейки граничные условия на внешней границе ЛВ запишем в виде [11]

фе = у, = 0, (2)

где = — Дфе - завихренность. На участках ЛВ и СВ выполняется условие симметрии

фе = 0, =0. (3)

В пористой области П* функция тока ф* (х,у) внутреннего течения в цилиндре удовлетворяет уравнению Бринкмана

Д2ф* - Б2Дф* = 0, (4)

где Б = Ис/л/к - безразмерный параметр, к - проницаемость пористого цилиндра. На участке ВОС задано условие симметрии

ф* = 0, и* = 0. (5)

На линии BC раздела двух сред (границе пористого цилиндра) выполняются условия [12]

uex = uX, uey = uy, pe = pi, = Jí, (6)

где ux, uy - компоненты вектора скорости течения, p - давление. Решение краевой задачи (1)-(6) позволяет определить функции фе,1(х,у], ue,i(x,y), uexi (х,у),

uyi(x,y).

Решения для функций тока внешнего и внутреннего течений представляются в виде

фe(r,e) = ^ — + Bir + Cir ln r + Dir3^ sin в, (7)

ф'1(т,в) = (B2r + C2I1(Sr)) sin в, (8)

где r, в - полярные координаты, Ii - модифицированная функция Бесселя первого рода первого порядка. Подставив (7) и (8) в граничные условия (5) и (6), получим систему линейных алгебраических уравнений для нахождения неизвестных Ai, Bi, Ci, Di, B2 , C2 :

aAi + Bi - 1Па Ci + 1 Di = 1,

2 a

aCi + 4Di = 0,

Ai + Bi + Di - B2 - Ii(S)C2 =0,

-Ai + Bi + Ci + 3Di - B2 + (Ii(S) - SIo(S))C2 = 0,

2Ci - 8Di - S2B2 = 0,

2Ci +8Di - S2Ii(S)C2 = 0,

где Io - модифицированная функция Бесселя первого рода нулевого порядка. Компоненты скорости течения выражаются через найденные фe(r, в) и фi (r, в) по формулам:

1 dфe (A

uer М) = - = ( ^ + Bi + Ci ln r + Dir2 ) cos в, r dti \r2 '

4 M) = i ^ = (B + C2 ^Jcoe в, r dti \ r 1

u«M) = - = ( ^T - Bi - Ci(ln r + 1) - 3Dir2 ) sin в,

rlf \ f2 '

d^e = í at

dr r2

?(r, в) = - di = - (B2 + C2S (Io(Sr) - ^J ) sin в.

Рис. 2. Прямоугольная периодическая ячейка с пористым цилиндром

2. Прямоугольная периодическая ячейка

Модель круговой периодической ячейки Кувабара позволяет получить аналитическое решение задачи Стокса - Бринкмана, однако она дает приближенное поле скоростей течения. Для обеспечения большей точности описания поля течения в ряду или упаковке пористых цилиндров рассмотрим задачу Стокса - Бринкмана для прямоугольной периодической ячейки (рис. 2). В области Qe внешнего течения функция тока фе(х,у) удовлетворяет бигармоническому уравнению (1). На границах AF и DE задаются граничные условия

Г(-Ьл,у) = у, фе(Н1,у) = у, ше(-Н1,у)=0, ше(Ниу)=0. (9)

Условия симметрии на границе EF задаются в виде

фе = h, ле = 0 (10)

и на оси симметрии ячейки (границы AB и CD) в виде

фе = 0, ле =0. (11)

В области пористого цилиндра Q® функция тока ф®(х,у) внутреннего течения удовлетворяет уравнению Бринкмана (4). На границе BOC выполняется условие симметрии

фi = 0, Л =0. (12)

Условия на границе Г® однородной и пористой областей имеют вид (6).

Приведенная краевая задача (1), (4), (6), (9)-(12) решается с помощью метода граничных элементов [12].

3. Уравнения движения аэрозольных частиц

В приближении стоксова аэродинамического сопротивления уравнения движения взвешенной частицы записываются в виде [13]

¿V и — V дГ _ . ,

, (13)

3Ь т дЬ

где и, V - скорость несущей среды и частицы соответственно, г - радиус-вектор частицы, т = ррдр/\8/ время релаксации сферической частицы, / - вязкость среды, рр и др - плотность и диаметр частицы. Для рассматриваемой двумерной задачи уравнения (13) представляются в безразмерном виде

dvx ux — vx dvy uy — vy dx d'у

dt St ' dt St ' dt xx dt

(14)

v

y

где их, иу и ух, уу - декартовы составляющие скорости несущей среды и частицы соответственно, St = Пот/Кс - число Стокса. Решение задачи Коши для системы обыкновенных дифференциальных уравнений (14) в найденном поле скоростей несущей среды с соответствующими начальными условиями позволяет рассчитать траекторию отдельной частицы. Значения составляющих скорости газа в текущей точке траектории частицы находятся на основе билинейной интерполяции во множестве значений, полученных из численного решения уравнений движения несущей среды. В области решения вводится равномерная прямоугольная сетка {Х®,у^}, в узлах которой рассчитывается множество значений декартовых компонент скорости течения их. ., иу. .. Затем текущие значения компонент скорости их(х,у), иу(х, у) вычисляются по формулам

их = 1 [пх^ (1 — 0(1 — п) + ихг+1„ (1 + 0(1 — п) +

(1+о(1+п)+их«,3+1 (1 — оа + п)],

иу =1 [иуг, з (1 — оа — п)+иуг+1,. (1+оа — п)+

+ пУi+1, 3 + 1

(1+оа+п)+иуг 3+1 (1 — оа + п)],

£ = 2——х--1, п = 2——---1.

х®+1 — х® у®+1 — у®

4. Результаты расчетов

Плотность периодической ячейки определяется как отношение площади области внешнего течения ко всей площади ячейки. В случае прямоугольной ячейки под областью внешнего течения понимается квадратная область размера 2Н х 2Н. Для круговой и квадратной ячеек плотность ячейки а находится по формулам а = Н-2 и а = п/АН? соответственно. Для пористого цилиндра вводится безразмерный параметр - число Дарси Ба = к/К2с, который связан с параметром й формулой Ба = 1/Бс.

Задача Коши для системы уравнений (14) движения частиц решается численно методом Рунге-Кутта для полей течений, определенных в разд. 1 и 2. В случае прямоугольной ячейки расчет введется в области с Н1 = 5Н. В качестве начальных условий вдали от цилиндра задаются положение частиц и начальная скорость частиц равная скорости газа. Проведены расчеты эффективности Е улавливания частиц в результате инерционной импакции в зависимости от числа Стокса. Для этого определялась предельная траектория, разделяющая поток частиц, попадающий на цилиндр и проходящий мимо. Безразмерная ордината предельной траектории представляет собой эффективность осаждения частиц на цилиндре. Для расчета эффективности осаждения в результате как инерционной им-пакции, так и механизма зацепления определяется предельная траектория, которая попадает на цилиндр радиуса 1 + 5 (5 = др/2Кс - параметр зацепления).

С целью тестирования развитого метода проведены сравнения поля течения и эффективности осаждения частиц в прямоугольной ячейке для сплошного цилиндра с результатами работы [14]. На рис. 3 показаны распределения их(0,у) скорости течения несущей среды. Наблюдается хорошее согласие рассчитанных распределений из двух работ. На рис. 4, а приведены кривые эффективности улавливания частиц за счет инерционной импакции и механизма захвата для сплошного

- /1=2 -— /1=3.33 ....... /1=10 * Muller et a L, 2014 [14]

f / / / / i

* *

0.1 0.2 0.3 0.4 0.5

У/2 h

Рис. 3. Распределения скорости течения их(0,у) для различных расстояний между цилиндрами

цилиндра в зависимости от числа Стокса для 6 = 0.1 и различных Н. Кривые для различных Н могут быть сведены в одну при введении модифицированного числа Стокса = 0.5St ', /' = (4а/п)1/3 [14] (рис. 4, б). Сравнение кривых, полученных по развитой численной модели течения несущей среды и на основе численного решения в СЕБ пакете в [14], показывает хорошее согласие.

На рис. 5 представлены зависимости для различных чисел Дарси, рас-

считанные по описанным аналитической и численной моделям и на основе расчета в среде пакета А^УВ/ЕЬиЕМТ в [9]. Видно, что все модели дают близкие зависимости, что говорит о применимости аналитической модели круговой ячейки в приближении Стокса -Бринкмана для описания поля течения в периодической ячейке с пористым цилиндром. В случае сплошного цилиндра, начиная с некоторого числа Стокса, эффективность захвата частиц становится равной нулю, так как все частицы с малыми числами Стокса огибают цилиндр вместе с несущей средой. Для пористого цилиндра часть газа начинает протекать через пористое тело, увлекая за собой взвешенные частицы. Это приводит к улавливанию частиц при малых числах Стокса. Эффективность улавливания частиц при St ^ 0

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

В работе [15] получена формула для расчета эффективности Е$ улавливания частиц за счет зацепления для сплошного цилиндра в ячейке Кувабара

¿2

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

E + Д)2/3(1-а) ,

k(a) = -0.5lnа - 0.75 + а - 0.25а2.

- /1=10 -— h=5 ....... /1=2 • Muller et al.. 2014 Г141 . *" *

/

—*-*—♦ i» * ,*i

Рис. 4. Зависимости Е(Я^ (а) и Е(Б^) (б) для различных Н

Из аналитического решения (7) выразим эффективность улавливания безынерционных частиц ^ = 0) за счет зацепления на пористом цилиндре в виде

ESp(a, S, Da) = ф(1 + S, п/2) - ф(1, п/2) =

1 + S

AIo(S) + BIi(S) CIo(S) + DIi(S)

-16Ii(S)

CIo(S)+ DIi(S)'

A = 4S(a - 1)7,

B = 7[S2(2 + a7) - 8(1 + a)] - 4S2(1 + S)2 ln(1 + S) - 16, C = -4S(a - 1)2,

D = S2(2ln a + 3 - 4a + a2) - 8(1 + 2a - a2),

(16)

S = 1/%/Dä, y = S(2 + S).

1

Ю"1 10° ю1

Рис. 5. Зависимость Е(Я^ для пористого цилиндра при а = 0.05 и различных Ба

- СГ = 0.01

----а = 0.05

....... сг = 0.1 * Круговая ячейка .. ■ г.*-»-- >~ . *........ »■ ..............—

>'' г/

* / у

у' У

10° ю1

Рис. 6. Зависимость для пористого цилиндра при Ба =0.1 для различных а

Сравнение зависимостей (а, 5) и Е^р(а,5, Ба) при а = 0.05 показано на рис. 7. В случае пористого цилиндра эффективность улавливания частиц за счет зацепления больше, чем эффективность улавливания сплошного цилиндра. Разница растет с ростом числа Дарси. Зависимости эффективности Е(Б1т) улавливания аэрозольных частиц пористым цилиндром, учитывающей как инерционную импакцию, так и эффект зацепления, для 5 = 0.1; 0.5 и Ба = 0; 0.1 при а = 0.05 показаны на рис. 8. Для сравнения приведены соответствующие зависимости для сплошного цилиндра из [14], которые практически совпадают с полученными в настоящей работе для Ба = 0 .В области малых чисел Стокса улавливание частиц обеспечивается как механизмом зацепления, так и течением через пористый цилиндр. Для больших чисел Стокса величина Е(Я^) стремится к величине 1 + 5 как для сплошного, так и пористого цилиндра.

б

Рис. 7. Зависимости Е$(а.,5) и Е$р(а.,5, Юа) при а = 0.05

---- 0а=0.1 - Ра=0 -

* МиИег^а! , 2014 [14] < 6 = 0.5

/ / в = 0.1

/ //У'

/ / /

/7

Ю"1 10° ю1 ю2

эй

Рис. 8. Зависимость Е(Я^) для 5 = 0.1; 0.5 и Юа = 0; 0.1 при а = 0.05

Заключение

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

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

Благодарности. Работа выполнена при финансовой поддержке Российского научного фонда в рамках научного проекта № 22-21-00176.

Литература

1. Кирш В.А. Аэрозольные фильтры из пористых волокон // Коллоидный журн. -1996. - Т. 58, № 6. - С. 786-790.

2. Стечкина И.Б. Сопротивление пористых цилиндров в потоке вязкой жидкости при малых числах Рейнольдса // Изв. АН СССР. Механика жидкости и газа. - 1979. -№ 6. - C. 122-124.

3. Kirsh V.A. Stokes flow in periodic systems of parallel cylinders with porous permeable shells // Colloid J. - 2006. - V. 68, No 2. - P. 173-181. - doi: 10.1134/S1061933X06020086.

4. Deo S., Filippov A., Tiwari A., Vasin S., Starov V. Hydrodynamic permeability of aggregates of porous particles with an impermeable core // Adv. Colloid Interface Sci. -2011. - V. 164, No 1-2. - P. 21-37. - doi: 10.1016/j.cis.2010.08.004.

5. Yu P., Zeng Y., Lee T.S., Bai H.X., Low H.T. Wake structure for flow past and through a porous square cylinder // Int. J. Heat Fluid Flow. - 2010. - V. 31, No 2. - P. 141-153. -doi: 10.1016/j.ijheatfluidflow.2009.12.009.

6. Kirsh V.A. Stokes flow in model fibrous filters // Sep. Purif. Technol. - 2007. - V. 58, No 2. - P. 288-294. - doi: 10.1016/j.seppur.2007.07.002.

7. Bhattacharyya S., Dhinakazan S., Khalili A. Fluid motion around and through a porous cylinder // Chem. Eng. Sci. - 2006. - V. 61, No 13. - P. 4451-4461. - doi: 10.1016/j.ces.2006.02.012.

8. Kirsch V.A. Inertial deposition of aerosol particles in a model filter with dust-loaded fibres // Trans. Filtr. Soc. - 2002. - V. 2, No 4. - P. 109-113.

9. Zaripov S.K., Solov'eva O.V., Solov'ev S.A. Inertial deposition of aerosol particles in a periodic row of porous cylinders // Aerosol Sci. Technol. - 2015. - V. 49, No 6. -P. 400-408. - doi: 10.1080/02786826.2015.1036834.

10. Хаппель Дж., Бреннер Г. Гидродинамика при малых числах Рейнольдса. - М.: Мир, 1976. - 630 с.

11. Kuwabara S. The forces experienced by randomly distributed parallel circular cylinders or spheres in a viscous flow at small Reynolds numbers //J. Phys. Soc. Jpn. - 1959. -V. 14, No 4. - P. 527-532. - doi: 10.1143/JPSJ.14.527.

12. Mardanov R.F., Zaripov S.K., Sharafutdinov V.F., Dunnett S.J. A Stokes-Brinkman model of the fluid flow in a periodic cell with a porous body using the boundary element method // Eng. Anal. Boundary Elem. - 2018. - V. 88. - P. 54-63. - doi: 10.1016/j.enganabound.2017.12.006.

13. Райст П. Аэрозоли. Введение в теорию. - М.: Мир, 1987. - 278 с.

14. Miiller T.K., Meyer J., Kasper G. Low Reynolds number drag and particle collision efficiency of a cylindrical fiber within a parallel array //J. Aerosol Sci. - 2014. - V. 77. -P. 50-66. - doi: 10.1016/j.jaerosci.2014.07.007.

15. Lee K.W., Gieseke J.A. Note on the approximation of interceptional collection efficiencies // J. Aerosol Sci. - 1980. - V. 11, No 4. - P. 335-341. - doi: 10.1016/0021-8502(80)90041-5.

Поступила в редакцию 07.02.2022

Хазиев Айрат Ринатович, аспирант кафедры аэрогидромеханики Казанский (Приволжский) федеральный университет

ул. Кремлевская, д. 18, г. Казань, 420008, Россия E-mail: aira-6991@mail.ru

Зарипов Шамиль Хузеевич, доктор физико-математических наук, профессор, заведующий кафедрой моделирования экосистем

Казанский (Приволжский) федеральный университет

ул. Кремлевская, д. 18, г. Казань, 420008, Россия E-mail: shamil.zaripov@kpfu.ru

Марданов Ренат Фаритович, кандидат физико-математических наук, доцент кафедры аэрогидромеханики

Казанский (Приволжский) федеральный университет

ул. Кремлевская, д. 18, г. Казань, 420008, Россия E-mail: Renat.Mardanov@kpfu.ru

Пилюгин Александр Геннадьевич, старший преподаватель кафедры моделирования экосистем

Казанский (Приволжский) федеральный университет

ул. Кремлевская, д. 18, г. Казань, 420008, Россия E-mail: pag@kpfu.ru

ISSN 2541-7746 (Print) ISSN 2500-2198 (Online) UCHENYE ZAPISKI KAZANSKOGO UNIVERSITETA. SERIYA FIZIKO-MATEMATICHESKIE NAUKI (Proceedings of Kazan University. Physics and Mathematics Series)

2022, vol. 164, no. 2-3, pp. 181-193

ORIGINAL ARTICLE

doi: 10.26907/2541-7746.2022.2-3.181-193

Efficiency of Capture of Inertia Aerosol Particles in a Periodic Cell with a Porous Cylinder

A.R. Khaziev*, Sh.Kh. Zanpov**, R.F. Mardanov***, A.G. Pilyugin****

Kazan Federal University, Kazan, 420008 Russia E-mail: *aira-6991@mail.ru, **shamil.zaripov@kpfu.ru, ***Renat.Mardanov@kpfu.ru,

**** pag@kpfu.ru

Received February 7, 2022 Abstract

In this study, aerosol flow in circular and rectangular periodic cells of a regular row of porous cylinders is simulated. The hydrodynamic field of the carrier medium flow outside and in the domain of the porous cylinder is described in the approximation of the Stokes-Brinkman

model for an incompressible gas. For the model of a circular periodic Kuwabara cell, an analytical solution is obtained; in the case of a rectangular cell, the boundary value problem for the flow stream function is solved numerically using the boundary element method. In the determined velocity fields of the carrier medium, the Lagrangian equations of motion of suspended particles are integrated, and the dependence of the particle capture efficiency as a result of the inertial impact and the interception mechanism on the Stokes number at different cell density and the Darcy number of a porous cylinder is calculated. The results of the calculations are found to be in a good agreement with those from other works carried out using the CFD package. Parametric studies of the efficiency of capture of aerosol particles as a function of the Stokes number for various values of the Darcy number of a porous cylinder, the periodic cell density, and the interception parameter are performed. It is shown that when a porous cylinder is flown around, the particle capture efficiency tends to a non-zero value that increases depending on how high the Darcy number and the density of the periodic cell are. In the region of small Stokes numbers, the capture of particles is provided by both the interception mechanism and the flow through the porous cylinder. A formula is derived for estimating the efficiency of particle trapping under the effect of particle interception on a porous cylinder.

Keywords: suspended particles, porous cylinder, Stokes-Brinkman model, particle trapping efficiency

Acknowledgments. This study was supported by the Russian Science Foundation (project no. 22-21-00176).

Figure Captions

Circular periodic cell with a porous cylinder. Rectangular periodic cell with a porous cylinder.

Distributions of flow velocity ux(0,y) for different distances between the cylinders. Dependencies E(St) (a) and E(Stm) (b) for different h. Dependence E(St) for a porous cylinder at a = 0.05 and different Da. Dependence E(St) for a porous cylinder at Da = 0.1 for different a. Dependencies Eg(a,S) and Egp(a,5, Da) for a = 0.05. Dependence E(Stm) for 5 = 0.1; 0.5 and Da = 0; 0.1 at a = 0.05.

References

1. Kirsh V.A. Aerosol filters made of porous fibers. Colloid J., 1996, vol. 58, no. 6, pp. 737740.

2. Stechkina I.B. Drag of porous cylinders in a viscous fluid at low Reynolds numbers. Fluid Dyn., 1979, vol. 14, no. 6, pp. 912-915. doi: 10.1007/BF01051997.

3. Kirsh V.A. Stokes flow in periodic systems of parallel cylinders with porous permeable shells. Colloid J., 2006, vol. 68, no. 2, pp. 173-181. doi: 10.1134/S1061933X06020086.

4. Deo S., Filippov A., Tiwari A., Vasin S., Starov V. Hydrodynamic permeability of aggregates of porous particles with an impermeable core. Adv. Colloid Interface Sci., 2011, vol. 164, nos. 1-2, pp. 21-37. doi: 10.1016/j.cis.2010.08.004.

5. Yu P., Zeng Y., Lee T.S., Bai H.X., Low H.T. Wake structure for flow past and through a porous square cylinder. Int. J. Heat Fluid Flow, 2010, vol. 31, no. 2, pp. 141-153. doi: 10.1016/j.ijheatfluidflow.2009.12.009.

6. Kirsh V.A. Stokes flow in model fibrous filters. Sep. Purif. Technol., 2007, vol. 58, no. 2, pp. 288-294. doi: 10.1016/j.seppur.2007.07.002.

7. Bhattacharyya S., Dhinakazan S., Khalili A. Fluid motion around and through a porous cylinder. Chem. Eng. Sci., 2006, vol. 61, no. 13, pp. 4451-4461. doi: 10.1016/j.ces.2006.02.012.

Fig. 1.

Fig. 2.

Fig. 3.

Fig. 4.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Fig. 5.

Fig. 6.

Fig. 7.

Fig. 8.

8. Kirsch V.A. Inertial deposition of aerosol particles in a model filter with dust-loaded fibres. Trans. Filtr. Sac., 2002, vol. 2, no. 4, pp. 109-113.

9. Zaripov S.K., Solov'eva O.V., Solov'ev S.A. Inertial deposition of aerosol particles in a periodic row of porous cylinders. Aerasal Sci. Technal., 2015, vol. 49, no. 6, pp. 400-408. doi: 10.1080/02786826.2015.1036834.

10. Happel J., Brenner H. Gidradinamika pri malykh chislakh Reinal'dsa [Low Reynolds Number Hydrodynamics]. Moscow, Mir, 1976. 630 p. (In Russian)

11. Kuwabara S. The forces experienced by randomly distributed parallel circular cylinders or spheres in a viscous flow at small Reynolds numbers. J. Phys. Sac. Jpn., 1959, vol. 14, no. 4, pp. 527-532. doi: 10.1143/JPSJ.14.527.

12. Mardanov R.F., Zaripov S.K., Sharafutdinov V.F., Dunnett S.J. A Stokes-Brinkman model of the fluid flow in a periodic cell with a porous body using the boundary element method. Eng. Anal. Baundary Elem., 2018, vol. 88, pp. 54-63. doi: 10.1016/j.enganabound.2017.12.006.

13. Reist P.C. Aerazali. Vvedenie v teariyu [Introduction to Aerosol Science]. Moscow, Mir, 1987. 278 p. (In Russian)

14. Muller T.K., Meyer J., Kasper G. Low Reynolds number drag and particle collision efficiency of a cylindrical fiber within a parallel array. J. Aerasal Sci., 2014, vol. 77, pp. 50-66. doi: 10.1016/j.jaerosci.2014.07.007.

15. Lee K.W., Gieseke J.A. Note on the approximation of interceptional collection efficiencies. J. Aerasal Sci., 1980, vol. 11, no. 4, pp. 335-341. doi: 10.1016/0021-8502(80)90041-5.

Для цитирования: Хазиев А.Р., Зарипов Ш.Х., Марданов Р.Ф., Пилюгин А.Г. Эффективность улавливания инерционных аэрозольных частиц в периодической ячейке \ с пористым цилиндром // Учен. зап. Казан. ун-та. Сер. Физ.-матем. науки. - 2022. -Т. 164, кн. 2-3. - С. 181-193. - doi: 10.26907/2541-7746.2022.2-3.181-193.

For citation: Khaziev A.R., Zaripov Sh.Kh., Mardanov R.F., Pilyugin A.G. Efficiency / of capture of inertia aerosol particles in a periodic cell with a porous cylinder. Uchenye \ Zapiski Kazanskogo Universiteta. Seriya Fiziko-Matematicheskie Nauki, 2022, vol. 164, no. 2-3, pp. 181-193. doi: 10.26907/2541-7746.2022.2-3.181-193. (In Russian)

i Надоели баннеры? Вы всегда можете отключить рекламу.