CC BY
62
___________УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО УНИВЕРСИТЕТА
Том 156, кн. 3 Естественные науки
2014
УДК 504.75.05+532.529.6
РАСЧЕТ ВДЫХАЕМОЙ ФРАКЦИИ ДИСПЕРСНЫХ ЗАГРЯЗНЕНИЙ С УЧЕТОМ ЗАЩИТНОЙ МАСКИ
И.Т. Мухаметзанов
Аннотация
Представлена математическая модель и результаты расчетов вдыхаемой фракции дисперсных загрязнений при дыхании человека с защитной маской. Для идеализированной сферической модели головы человека воздушное течение через защитную маску описывается течением через пористый слой. Проведены параметрические исследования потоков частиц, попадающих на маску и проходящих в щель между маской и лицом, в зависимости от проницаемости маски.
Ключевые слова: защитная маска, вдыхаемая фракция, пористый слой, математическая модель.
Введение
Для защиты органов дыхания от дисперсных воздушных загрязнений, особенно инфекционной природы, широко используются различные лицевые маски и респираторы. Ношение защитных масок является обязательным для людей, работающих в области здравоохранения. Однако эффективность респиратора или маски зависит не только от эффективности материала фильтра, но и от плотности прилегания к лицу. При неплотном облегании образуются щели между маской и лицом, через которые могут проникать пылевые частицы и попадать в дыхательные пути через ротовое или носовое отверстия. Эксперименты, проведенные в работах [1, 2], показали, что даже небольшие щели между респиратором и поверхностью лица могут существенно снизить эффективность респиратора. Характеристики респиратора с отверстием в фильтрующем слое для манекена головы человека исследовались в [1, 3-7]. В [8] была развита математическая модель течения воздуха при дыхании человека через маску для реальной геометрии головы человека. Европейский FFP2 и американский N95 классы защиты широко используемых респираторов предполагают эффективность фильтрации запыленного воздуха не менее 0.95 в области размеров частиц порядка 0.1 мкм и боковую протечку воздуха менее 8% [9].
Детальные исследования эффективности используемых защитных масок при наличии щелей остаются актуальными для оптимизации параметров масок и возможной их модернизации. Учет реальной геометрии манекенов головы человека осложняет проведение необходимого числа расчетов при изменении ключевых параметров задачи (проницаемость материла маски, размеры частиц и щели и т. п.) из-за большого компьютерного времени единичного расчета. Создание адекватной приближенной модели могло бы способствовать решению задачи.
133
134
И.Т. МУХАМЕТЗАНОВ
Ранее в [10] было показано, что вдыхаемая фракция пылевых загрязнений хорошо описывается коэффициентом аспирации идеализированного сферического пробоотборника.
Целью настоящей работы является создание математической модели и проведение параметрических исследований вдыхаемой фракции дисперсных загрязнений при дыхании человека через защитную маску. Для описания геометрии головы человека предлагается двумерная идеализированная модель, представляющая собой сферу с аспирирующим отверстием, перед которым расположен пористый слой - защитная маска. В приближении малых концентраций дисперсной фазы двухфазное течение описывается уравнениями движения несущей среды и аэрозольных частиц. Проведены численные исследования вдыхаемой фракции при различных значениях параметров, характеризующих геометрию и материал фильтрующего слоя маски.
Математическая постановка задачи
Рассматривается стационарное двумерное осесимметричное течение газа с взвешенными частицами при обтекании сферы радиуса Rs с аспирирующим отверстием и фильтрующим слоем перед ним (рис. 1). Вдали от сферы воздух движется со скоростью U0. Дыхание моделируется аспирацией через круговое отверстие радиуса ra со скоростью Ua. В отсутствие пористого слоя рассматриваемая постановка полностью соответствует задаче аспирации аэрозоля в сферический пробоотборник [11, 12].
Рис. 1. Г еометрия сферического пробоотборника с пористым слоем и предельные траектории
В предположении малости концентраций частиц влиянием дисперсной фазы на газовое течение пренебрегается. Движение несущей среды в однородной и пористой областях в приближении ламинарного вязкого течения газа описывается уравнениями [13]
V-U = 0, (1)
=-vp+^au - ь HU, s к (2)
где U - вектор скорости газа, р - коэффициент динамической вязкости, р -плотность воздуха, P - давление, е - пористость цилиндра, k - проницаемость
РАСЧЕТ ВДЫХАЕМОЙ ФРАКЦИИ ДИСПЕРСНЫХ ЗАГРЯЗНЕНИИ... 135
пористой среды. Коэффициент b задается равным нулю вне пористого слоя и единице внутри него. Скорость газа U, осредненная по объему порового пространства, связана со скоростью фильтрации Uf соотношением U = sU^ .
В однородной области (е = 1, b = 0) уравнения (1), (2) переходят в уравнения Навье - Стокса несжимаемой жидкости. Внутри пористого тела (1), (2) представляют собой расширенные уравнения Дарси - Бринкмана, учитывающие дополнительное сопротивление в пористой среде.
Уравнения движения газовой среды течения (1), (2) дополняются условиями на границах расчетной области, заключенной между двумя полуокружностями: границей сферы радиуса Rs и внешней границей радиуса Rb. Вдали от сферической поверхности на левой половине граничной окружности задается скорость ветра U0, на правой границе избыточное давление принимается равным нулю. На сфере задаются условия прилипания U = 0. В поперечном сечении входного отверстия пробоотборника задается постоянная скорость аспирации Ua. Область, занятая маской, определяется как пористая среда с проницаемостью k.
В общем случае на движение частиц могут влиять различные силы: аэродинамическое сопротивление, сила тяжести, диффузия, электростатические силы и т. п. Диффузия для частиц со значительной инерцией мала по сравнению с конвективной составляющей их движения. Поэтому, как правило, в задачах аспирации не учитывают диффузию частиц [14]. Уравнения движения взвешенной частицы в предположении стоксового закона сопротивления газовой среды без учета силы тяжести запишем в виде [15, 16]
dV
dt
fi(U - V),
drp _
-T = V, dt
(3)
где V - вектор скорости частицы, p - вектор-координата положения частицы, в = 3n^d/ mp , dp - диаметр частицы, mp - масса частицы.
Для количественной оценки вклада двух путей проникновения частиц вводится величина P, которая представляет собой отношение потока Qp частиц, проходящих через щель и попадающих в ротовое отверстие, к потоку Qp частиц,
попадающих на пористый слой. Для определения P рассчитываются предельные траектории, разделяющие указанные потоки частиц (рис. 1). Пусть rm0 - радиус трубки траекторий вдыхаемых частиц, Гр - радиус трубки траекторий частиц, попадающих на пористый слой. Введенные величины rm0 и Гр позволяют рассчитать потоки Qp и Qp по формулам
Qf = coUof Q,p = rp-rno)
/0-1
(4)
где C0 - концентрация частиц в невозмущенном потоке. Тогда величина P определяется из формулы
P =
QL
QP
( Y
V rf0 J
-1.
(5)
136
И.Т. МУХАМЕТЗАНОВ
Для оценки эффективности работы защитных масок используется понятие коэффициента защиты Pf [17], представляющего собой отношение концентрации С0 частиц в невозмущенной среде к концентрации частиц C, попадающих в дыхательные пути после прохождения защитной маски: р = C0 / Ci. В то же
время в теории пробоотбора аэрозольных частиц вводится понятие коэффициента аспирации А как отношение концентрации С, частиц, попадающих в измерительное устройство, к концентрации частиц С0 в невозмущенной среде: A = C / C . Таким образом, коэффициент защиты может быть выражен как обратная величина коэффициента аспирации А
р = 1/A. (6)
Отметим, что коэффициент аспирации соответствует понятию вдыхаемой фракции, определяющей долю вдыхаемых человеком пылевых частиц [18]. В общем случае величина A зависит от характеристик самой частицы (размера, плотности, формы), свойств газового потока (отношения скоростей ветра и аспирации Ra = U0 /Ua), геометрии ротового или носового отверстия (размера, формы),
ориентации человека относительно направления ветра и направления силы тяжести. В случае учета защитной маски к параметрам, влияющим на вдыхаемую фракцию, добавляются геометрические размеры маски, свойства фильтрующего материала, форма и размер щели между маской и лицом человека.
В предположении осевой симметрии коэффициент аспирации сферического пробоотборника может быть вычислен по формуле [12]:
C U жг2 (г ^
^_C0U0'l'm0 _R rm0
C U жг2 a г C0Ua'U а V Га У
(7)
Определим коэффициент аспирации для сферического пробоотборника с учетом пористого слоя как отношение потоков частиц через фильтрующий слой и через щель к потоку безинерционных частиц через аспирирующее отверстие:
Qf^f + Qf4
C0 Ua^a
C0U0ЖГ/0^ + C0U0Ж (Гт0 - =R rf0lf + (Гт0 -
C0 Ua*ra
-,(8)
г
a
где nf и ni - коэффициенты пропуска частиц через пористый слой и через щель соответственно. В предположении идеального фильтрующего слоя (nf = 0) и полного пропуска частиц через щель (nl = 1) коэффициент защиты из (8) может быть записан в виде
Pf =
Ra (Гт20 - Г/о)
2
r
a
(9)
Результаты расчетов
Для экспериментальных исследований вдыхаемой фракции с учетом и без учета защитной маски использовались манекены головы человека с размерами, близкими к реальным размерам головы. Так, средняя высота и ширина головы манекена выбирались в диапазоне: 0.15-0.16 и 0.14-0.15 м [19], 0.216 и 0.142 м [20],
РАСЧЕТ ВДЫХАЕМОЙ ФРАКЦИИ ДИСПЕРСНЫХ ЗАГРЯЗНЕНИИ... 137
Рис. 2. Сеточное разбиение в окрестности пористого слоя
0.225 и 0.134 м [21]. Остановимся на геометрических параметрах из [21]. Поставим в соответствие эллипсу с высотой и шириной головы манекена из [21] круг радиуса 0.09 м такой же площади 0.0236 м2. Площади ротового отверстия Sa = 0.00016 м2 соответствует круговое отверстие радиусом 0.007 м. Таким образом, для численных исследований выбраны следующие значения геометрических параметров сферического пробоотборника: Rs = 0.09 м, ra = 0.007 м.
Защитный слой представляет собой сферический сегмент радиуса b = 0.057 м и высотой a = 0.05 м (рис. 1). Площадь внешней поверхности фильтрующего слоя составляет Sf = 0.0179 м2. Толщина слоя равномерна по всей площади 0.003 м. Приведенные размеры пористого слоя приближенно соответствуют параметрам широко используемых респираторов указанных ранее европейского FFP2 и американского N95 классов защиты. Размер щели e по краю защитной маски регулируется перемещением сегмента вдоль оси симметрии х. Размеры расчетной области внешнего течения подбирались достаточными, чтобы обтекаемое тело -сфера - не влияло на течение воздуха вдали от него: Rb = 20 м.
Различают три основных режима дыхания человека с различным расходом q вдыхаемого воздуха: спокойное, умеренное и интенсивное (при больших нагрузках). Для расхода воздуха, соответствующего спокойному дыханию q = 8 л/мин, скорость аспирации через принятое круговое отверстие будет равна Ua = 1.73 м/с. Скорость ветра в соответствии с характерными скоростями воздуха внутри помещений выбрана равной U0 = 0.1 м/с.
Система уравнений течения несущей среды (1), (2) решается с помощью метода конечных объемов в среде пакета ANSYS/Fluent [22]. Расчетная область вне и внутри пористой зоны разбивается на неструктурированные четырехугольные элементы с обеспечением более подробного разбиения на границе пористой зоны (рис. 2). Задача Коши для (3) в найденном поле скоростей газа решается численно на основе метода Рунге - Кутта.
138
И.Т. МУХАМЕТЗАНОВ
St
Рис. 3. Зависимости A (St) для сферического пробоотборника: Ra = 0.38
Рис. 4. Траектории частиц диаметром 1 (а) и 30 (б) мкм
В целях тестирования проводились расчеты коэффициента аспирации в сферический пробоотборник без защитной маски. Наблюдается хорошее согласие полученной зависимости коэффициента аспирации A от числа Стокса St = /3~lU0 /ra при Ra = 0.38 с расчетными данными из [11] (рис. 3).
Пример траекторий частиц диаметром dp = 1 и 30 мкм при проницаемости маски к = 8.9310-11 м2 дан на рис. 4. При дыхании через защитную маску вдыхаемый воздух проходит через фильтрующий слой. Часть вдыхаемого воздуха устремляется в щель между маской и лицом, что приводит к попаданию частиц в ротовое отверстие мимо фильтрующего слоя. Рис. 4 демонстрирует различное поведение частиц двух диаметров после попадания в щель. Часть частиц с dp = 30 мкм после прохождения сечения щели будет оседать на сферической поверхности, то есть не все частицы, проходящие в щель, будут попадать в ротовое отверстие.
Указанное обстоятельство учитывалось при вычислении величины Р. Различались отношения потоков частиц, проходящих через щель и через ротовое отверстие к потоку частиц, попадающих на пористый слой: Р1 и Р2 соответственно. Зависимость Р1 и Р2 от диаметра частиц для трех проницаемостей (к = 5• 10-11, 8.93• 10-11, 5• 10-10 м2) при размере щели e = 1 мм представлена на рис. 5. Выбранные значения проницаемостей в целом находятся в диапазоне значений, соответствующих проницаемостях, применяемых в материалах защитных масок.
РАСЧЕТ ВДЫХАЕМОЙ ФРАКЦИИ ДИСПЕРСНЫХ ЗАГРЯЗНЕНИИ... 139
Рис. 5. Зависимость Р1 (1, 3, 5) и Р2 (2, 4, 6) от диаметра частиц для e = 1 мм: к = 510 11 м2 (1, 2), к = 8.93 • 10-11 м2 (3, 4), к = 5-10-10 м2 (5, 6)
IQ-
10
1
2
3
I /
Ч
100
d, мкм
Рис. 6. Зависимость Pf (dp)
В зависимости от значения проницаемости материала маски доля частиц, проходящих через щель, может быть значительной. Так, для к = 8.93-1011 м2 Р\ и Р2 достигают значений, равных приблизительно 8, в области малых чисел Стокса. Близкие значения Р наблюдались в экспериментах [2]. Различие между Pi и Р2 наблюдается для частиц размером более 25 мкм. Предельный размер вдыхаемых через щель частиц составляет 25-35 мкм, в то время как в щель могут попадать частицы размером до 80-100 мкм. По мере уменьшения проницаемости маски доля частиц, проходящих через щель, возрастает.
Наряду с величиной Р рассчитывался коэффициент защиты маски Pf по формуле (8). Полученные зависимости Pf (dp) представлены на рис. 6. Видно, что значение коэффициента защиты для частиц с аэродинамическим диаметром менее 25 мкм мало и находится в диапазоне от 1 до 2. Как было отмечено во введении, принятые для классов FFP2/N95 степени фильтрации запыленного воздуха должны быть не менее 0.95. Таким образом, допустимыми значениями степени коэффициента защиты респираторов считаются величины, большие 20. Для рассматриваемой задачи приемлемый уровень защиты обеспечивается для
140
И.Т. МУХАМЕТЗАНОВ
частиц размером выше ~30 мкм. Примечательно, что для выбранных значений проницаемостей и коэффициентов пропуска фильтра и щели (nf = 0, pl = 1) максимальную защиту обеспечивает наиболее проницаемая маска. Однако для такой маски предположение пу = 0 может оказаться недостоверным. Часть частиц, попадающих на поверхность фильтровального слоя, может пройти сквозь пористый слой, не задерживаясь. Доля частиц, которые не улавливаются фильтром, будет возрастать с увеличением проницаемости. В общем случае следует в формуле (8) учитывать реальные значения коэффициента пропуска пу и вычислять коэффициент защиты с учетом новой формулы для коэффициента аспирации. Кроме того, необходимо отметить, что для всех трех проницаемостей величина Ру намного меньше 20, следовательно, эффективность защитной маски определяется в первую очередь размером щели.
Заключение
Предложена математическая модель и проведены численные исследования вдыхаемой фракции дисперсных загрязнений при дыхании человека с защитной маской. Модель включает в себя уравнения ламинарного течения несущей среды и уравнения движения взвешенных частиц с учетом аэродинамического сопротивления. Для идеализированной сферической модели головы человека воздушное течение через защитную маску описывается течением через пористый слой. Численное решение уравнений модели достигается в среде CFD пакета FLUENT. Проведены параметрические исследования потоков частиц, попадающих на маску и проходящих в щель между маской и лицом, и коэффициента защиты маски в зависимости от ее проницаемости. Показано, что в области значений проницаемостей, соответствующих материалам реальных респираторов, низкие значения эффективности защиты определяются в первую очередь размером щели между маской и лицом. Созданная модель и результаты расчетных исследований могут способствовать проектированию защитных масок с улучшенными характеристиками.
Работа выполнена при поддержке Российского фонда фундаментальных исследований (проекты № 12-01-00333, 14-01-31118).
Литература
1. Chen C.C., Willeke K. Characteristics of face seal leakage in filtering facepieces // Am. Ind. Hyg. Assoc. J. - 1992. - V. 53, No 9. - P. 533-539.
2. Grinshpun S.A., Haruta H., Eninger R.M., Reponen T., McKay R.T., Lee S.A. Performance of an N95 filtering facepiece particulate respirator and a surgical mask during human breathing: two pathways for particle penetration // J. Occup. Environ. Hyg. -2009. - V. 6, No 10. - P. 593-603. - doi: 10.1080/15459620903120086.
3. Chen C.C., Ruuskanen J., Pilacinski W., Willeke K. Filter and leak penetration characteristics of a dust and mist filtering facepiece // Am. Ind. Hyg. Assoc. J. - 1990. - V. 51, No 12. - P. 632-639.
4. Hinds W.C., Kraske G. Performance of Dust Respirators with Facial Seal Leaks: I. Experimental // Am. Ind. Hyg. Assoc. J. - 1987. - V. 48, No 10. - P. 836-841.
5. Hinds W.C., Bellin P. Performance of Dust Respirators with Facial Seal Leaks: II. Predictive Model // Am. Ind. Hyg. Assoc. J. - 1987. - V. 48, No 10. - P. 842-847.
РАСЧЕТ ВДЫХАЕМОЙ ФРАКЦИИ ДИСПЕРСНЫХ ЗАГРЯЗНЕНИИ... 141
6. Rengasamy S., Eimer B.C. Nanoparticle penetration through filter media and leakage through face seal interface of N95 filtering facepiece respirators // Ann. Occup. Hyg. -
2012. - V. 56, No 5. - P. 568-580. - doi: 10.1093/annhyg/mer122.
7. Rengasamy S., Eimer B.C. Total inward leakage of nanoparticles through filtering facepiece respirators // Ann. Occup. Hyg. - 2011. - V. 55, No 3. - P. 253-263. - doi: 10.1093/annhyg/meq096.
8. Lei Z., Yang J., Zhuang Z., Roberge R. Simulation and evaluation of respirator faceseal leaks using computational fluid dynamics and infrared imaging // Ann. Occup. Hyg. -
2013. - V. 57, No 4. - P. 493-506. - doi: 10.1093/annhyg/mes085.
9. Бобрик А.В., Хорошев П.В. Современные маски и респираторы в системе инфекционного контроля и обеспечения безопасности персонала в ЛПУ. - М.: Акварель, 2010. - 20 с.
10. GilfanovA.K., Mukhametzanov I.T., Zaripov S.K. Mathematical study of inhalable fraction in low velocity conditions // Digital Abstract book and Conference Handbook: European Aerosol Conf. - Granada, 2012. - C-WG09S1P33.
11. Dunnett S.J., Vincent J.H. A mathematical study of aerosol sampling by an idealised blunt sampler oriented at an angle to the wind: the role of gravity // J. Aerosol Sci. -2000. - V. 31, No 10. - P. 1187-1203. - doi: 10.1016/S0021-8502(00)00024-0.
12. Galeev R.S., Zaripov S.K. A theoretical study of aerosol sampling by an idealized spherical sampler in calm air // J. Aerosol Sci. - 2003. - V. 34, No 9. - P. 1135-1150.
13. Bhattacharyya S., Dhinakaran 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.
14. Wen X., Ingham D.B. Aspiration efficiency of a thin-walled cylindrical aerosol sampler at yaw orientations with respect to the wind // J. Aerosol Sci. - 2000. - V. 31, No 11. -P. 1355-1365. - doi: 10.1016/S0021-8502(00)00036-7.
15. Райст П. Аэрозоли. Введение в теорию. - М.: Мир, 1987. - 278 с.
16. Фукс Н.А. Механика аэрозолей. - М.: Изд-во АН СССР, 1955. - 353 с.
17. Liu B.Y., Lee J., Mullinsb H., Danisch S. Respirator leak detection by ultrafine aerosols: A predictive model and experimental study // Aerosol Sci. Technol. - 1993. - V. 19, No 1. - P. 15-26.
18. Мухаметзанов И.Т., Зарипов Ш.Х., Гильфанов А.К. Теоретическое исследование вдыхаемой фракции дисперсных воздушных загрязнений // Учен. зап. Казан. ун-та. Сер. Естеств. науки. - 2013. - Т. 155, кн. 1. - С. 50-60.
19. Yang J., Dai J., Zhuang Z. Simulating the interaction between a respirator and a headform using LS-DYNA // Computer-Aided Design Appl. - 2009. - V. 6, No 4. - P. 539-551.
20. Anthony T.R. Contribution of facial feature dimensions and velocity parameters on particle inhalability // Ann. Occup. Hyg. - 2010. - V. 54, No 6. - P. 710-725. - doi: 10.1093/annhyg/meq040.
21. KennedyN.J., Hinds W.C. Inhalability of large solid particles // J. Aerosol Sci. - 2002. -V. 33, No 2. - P. 237-255. - doi: 10.1016/S0021-8502(01)00168-9.
22. ANSYS Fluent. - URL: http://www.ansys.com.
Поступила в редакцию 18.07.14
Мухаметзанов Ильнар Тальгатович - лаборант кафедры моделирования экологических систем, Казанский (Приволжский) федеральный университет, г. Казань, Россия. E-mail: ilnar.mukhametzanov@mail.ru
142
И.Т. МУХАМЕТЗАНОВ
•к к к
CALCULATION OF INHALABLE FRACTION OF DISPERSE POLLUTION WITH REGARD TO PROTECTIVE MASK
I. T. Mukhametzanov Abstract
A mathematical model and calculation results for a fraction of disperse pollution inhaled while
breathing through a protective mask are presented. For an idealized spherical model of the human head,
the air flow moving through the protective mask is described by a flow moving through a porous layer.
Parametric studies of the particle fluxes penetrating through the porous layer into the leakage between
the mask and the face are conducted depending on the permeability of the mask material.
Keywords: protective mask, inhalable fraction, porous layer, mathematical model.
References
1. Chen C.C., Willeke K. Characteristics of face seal leakage in filtering facepieces. Am. Ind. Hyg. Assoc. J, 1992, vol. 53, no. 9, pp. 533-539.
2. Grinshpun S.A., Haruta H., Eninger R.M., Reponen T., McKay R.T., Lee S.A. Performance of an N95 filtering facepiece particulate respirator and a surgical mask during human breathing: two pathways for particle penetration. J. Occup. Environ. Hyg., 2009, vol. 6, no. 10, pp. 593-603. doi: 10.1080/15459620903120086.
3. Chen C.C., Ruuskanen J., Pilacinski W., Willeke K. Filter and leak penetration characteristics of a dust and mist filtering facepiece. Am. Ind. Hyg. Assoc. J., 1990, vol. 51, no. 12, pp. 632-639.
4. Hinds W.C., Kraske G. Performance of Dust Respirators with Facial Seal Leaks: I. Experimental. Am. Ind. Hyg. Assoc. J., 1987, vol. 48, no. 10, pp. 836-841.
5. Hinds W.C., Bellin P. Performance of Dust Respirators with Facial Seal Leaks: II. Predictive Model. Am. Ind. Hyg. Assoc. J., 1987, vol. 48, no. 10, pp. 842-847.
6. Rengasamy S., Eimer B.C. Nanoparticle penetration through filter media and leakage through face seal interface of N95 filtering facepiece respirators. Ann. Occup. Hyg., 2012, vol. 56, no. 5, pp. 568-580. doi: 10.1093/annhyg/mer122.
7. Rengasamy S., Eimer B.C. Total inward leakage of nanoparticles through filtering facepiece respirators. Ann. Occup. Hyg., 2011, vol. 55, no. 3, pp. 253-263. doi: 10.1093/annhyg/meq096.
8. Lei Z., Yang J., Zhuang Z., Roberge R. Simulation and evaluation of respirator faceseal leaks using computational fluid dynamics and infrared imaging. Ann. Occup. Hyg., 2013, vol. 57, no. 4, pp. 493-506. doi: 10.1093/annhyg/mes085.
9. Bobrik A.V., Khoroshev P.V. Modern Masks and Respirators in the System of Infection Control and Personnel Safety of Health Facilities. Moscow, Akvarel, 2010. 20 p. (In Russian)
10. Gilfanov A.K., Mukhametzanov I.T., Zaripov S.K. Mathematical study of inhalable fraction in low velocity conditions. Digital Abstract book and Conference Handbook: European Aerosol Conf., Granada, 2012, C-WG09S1P33.
11. Dunnett S.J., Vincent J.H. A mathematical study of aerosol sampling by an idealised blunt sampler oriented at an angle to the wind: the role of gravity. J. Aerosol Sci., 2000, vol. 31, no. 10, pp. 11871203. doi: 10.1016/S0021-8502(00)00024-0.
12. Galeev R.S., Zaripov S.K. A theoretical study of aerosol sampling by an idealized spherical sampler in calm air. J. Aerosol Sci., 2003, vol. 34, no. 9, pp. 1135-1150.
13. Bhattacharyya S., Dhinakaran 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.
14. Wen X., Ingham D.B. Aspiration efficiency of a thin-walled cylindrical aerosol sampler at yaw orientations with respect to the wind. J. Aerosol Sci., 2000, vol. 31, no. 11, pp. 1355-1365. doi: 10.1016/S0021-8502(00)00036-7.
15. Raist P. Aerosols. Introduction into Theory. Moscow, Mir, 1987. 278 p. (In Russian)
16. Fuks N.A. The Mechanics of Aerosol. Moscow, Izd. AN SSSR, 1955. 353 p. (In Russian)
РАСЧЕТ ВДЫХАЕМОЙ ФРАКЦИИ ДИСПЕРСНЫХ ЗАГРЯЗНЕНИИ... 143
17. Liu B.Y., Lee J., Mullinsb H., Danisch S. Respirator leak detection by ultrafine aerosols: A predictive model and experimental study. AerosolSci. Technol., 1993, vol. 19, no. 1, pp. 15-26.
18. Mukhametzanov I.T., Zaripov Sh.Kh., Gilfanov A.K. Theoretical study of inhalable fraction of dust air pollution. Uchenye Zapiski Kazanskogo Universiteta. Seriya Estestvennye Nauki, 2013, vol. 155, no. 1, pp. 50-60. (In Russian)
19. Yang J., Dai J., Zhuang Z. Simulating the interaction between a respirator and a headform using LS-DYNA. Computer-Aided Design Appl., 2009, vol. 6, no. 4, pp. 539-551.
20. Anthony T.R. Contribution of facial feature dimensions and velocity parameters on particle inhala-bility. Ann. Occup. Hyg, 2010, vol. 54, no. 6, pp. 710-725. doi: 10.1093/annhyg/meq040.
21. Kennedy N.J., Hinds W.C. Inhalability of large solid particles. J. Aerosol Sci., 2002, vol. 33, no. 2, pp. 237-255. doi: 10.1016/S0021-8502(01)00168-9.
22. ANSYS Fluent. Available at: http://www.ansys.com.
Received July 18, 2014
Mukhametzanov Ilnar Talgatovich - Laboratory Assistant, Department of Modeling of Ecological Systems, Institute of Ecology and Geography, Kazan Federal University, Kazan, Russia.
E-mail: ilnar.mukhametzanov@mail. ru
