CC BY
32
УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА Том 150, кн. 4 Физико-математические пауки 2008
УДК 532.529.6
РАСЧЕТ ХАРАКТЕРИСТИК ТОНКОСТЕННОГО ПРОБООТБОРНИКА ПРИ БОЛЬШИХ ЗНАЧЕНИЯХ ОТНОШЕНИЯ СКОРОСТЕЙ ВЕТРА И АСПИРАЦИИ
А.К. Гильфшпов
Аннотация
Рассмотрена задача об аспирации аэрозоля в тонкостенную трубку из высокоскоростного потока. Поле течения несущей среды найдено па основе метода граничных элементов и метода конечных объемов для моделей потенциального и вязкого течений несжимаемой жидкости соответственно. Исследовано поведение коэффициента аспирации при высоких значениях отношения скоростей ветра и аспирации. Построены кривые зависимости коэффициента аспирации от числа Стокса и от отношения скоростей. Дано сравнение рассчитанных кривых с экспериментальными результатами и приближенными формулами. Ключевые слова: пробоотбор аэрозоля, коэффициент аспирации.
Введение
Для определения концентрации аэрозольных частиц в воздушной среде широко применяются пробоотборники устройства для отбора проб аэрозоля. Отбор аэрозольных проб часто осуществляется из движущегося воздуха (отбор в атмосфере в условиях сильных ветровых потоков) в иеизокинетических условиях, когда концентрация частиц внутри пробоотборника может заметно отличаться от концентрации частиц в изучаемом аэрозоле. Для коррекции искажений в теории пробоотбора аэрозолей вводится коэффициент аспирации А, представляющий собой отношение средней концентрации в пробоотборнике к счетной концентрации в невозмущенной среде. Обзор теоретических и экспериментальных работ, посвященных определению коэффициента аспирации для различных типов пробоотборников и условий пробоотбора. содержится в монографии [1].
Наиболее изученными в настоящее время являются пробоотборники в виде тонкостенной трубки, осуществляющей пробоотбор из движущегося газа. Обзор экспериментальных исследований аспирации в трубку приводится в работах [1. 2]. В работе [3] коэффициент аспирации для тонкостенной трубки рассчитывался в приближении осесимметричного потенциального течения несущей среды на основе численного решения уравнения для функции тока методом конечных разностей. Потенциальная модель поля течения несущей среды была использована в [4] для приближенного расчета коэффициента аспирации в трубку на основе решения линеаризованных уравнений движения частиц, пригодных или для больших или для малых чисел Стокса. Решение задачи аспирации аэрозоля в трубку в приближении вязкого течения было получено в [5 8].
В настоящей работе описывается математическая модель и приводятся результаты исследования поведения коэффициента аспирации тонкостенной трубки для случая, когда отношение Еа скорости ветра к скорости аспирации может достигать больших значений. Такая ситуация реализуется при очень высоких скоростях
Рис. 1. Схема течения аэрозоля при аспирации в трубку при Ка > 1
ветра (пробоотбор в атмосфере) или низких значениях скорости аспирации. Поведение коэффициента А в области больших значений Еа недостаточно изучено. В работе [9] экспериментально исследовался пробоотбор в трубку при Ка > 1, вплоть до значений Еа = 50, и было показано, что зависимость коэффициента аспирации от Ка имеет нелинейный характер при больших Ка. Коэффициент аспирации в трубку в рамках модели осесимметричного потенциального течения при высоких значениях отношения Еа скоростей ветра и аспирации рассчитывался в [10]. В настоящей статье для описания поля течения несущей среды используются модели потенциального и вязкого течений несжимаемой жидкости.
1. Постановка задачи
Рассматривается течение аэрозоля при аспирации в тонкостенную трубку кругового сечения радиуса Я1 (рис. 1). Вдали от пробоотборника несущая среда движется равномерно со скоростью и0. Направление вектора скорости ветра и0 совпадает с осью трубки и направлением скорости аспирации Па. Введем безразмерный параметр Ка = и0/иа. Осевая симметрия позволяет рассматривать задачу в цилиндрической системе координат (Х,Я) в меридиональной плоскости. Осевая координата X направлена вдоль оси пробоотборной трубки. Пусть
и Кс ординаты разделительной линии тока и предельной траектории в меридиональной плоскости вдали от пробоотборника. Разделительная линия тока определяет области течения воздуха, попадающего и не попадающего в пробоотборник. Предельной будем называть траекторию, разделяющую дисперсную фазу на потоки аспирируемых частиц и частиц, проходящих мимо пробоотборника. Используя условие баланса массы для несущей фазы, величину Я3 можно выразить как = Я,^ЩЩ = Щ/у/Щ.
В предположение об однородности параметров невозмущенного потока задача вычисления коэффициента аспирации сводится к нахождению в меридиональной плоскости предельной траектории, ордината Ес которой позволяет вычислить коэффициент аспирации А по формуле
А = Я2/Я2 = Яс2Яа/Яг2- (1)
При малых концентрациях частиц задача моделирования течения аэрозоля сводится к двум отдельным: определению поля скоростей несущей среды и расчету траекторий в найденном поле скоростей. В пренебрежении всеми силами, кроме аэродинамического сопротивления, представим систему уравнений, описывающих движение аэрозольной частицы в цилиндрических координатах, в безразмерном
ВИД6
дюх их — ух диг иг — уг дх дг
(2)
Л Яі ’ 'II Яі ’ ей А ’ ей
где ух, уг и их, иг - компоненты скоростей частицы и несугцей среды, х, г -координаты частицы, St = С2 ри0/\8рС - число Стокса, й р - диаметр и плотность частицы, ^ - коэффициент динамической вязкости среды, Б = 2Д4. Радиус трубки Кг и скорость тесу щей среды П0 выбираются в качестве масштаба для безразмерных координат и скоростей.
При скоростях, характерных для аспирации аэрозоля в типичных условиях, сжимаемостью воздуха можно пренебречь. Пусть движение несущей среды представляет собой осесимметричное потенциальное течение несжимаемой жидкости. Введем функцию тока фр(х,г), удовлетворяющую уравнению
Тогда безразмерные составляющие их, иг скорости газа в меридиональной плоскости (х, г) определяются по формулам
где ф3 (х,г) - функция тока течения, создаваемого аспирацией газа из неподвижной среды вдали от трубки при единичной скорости внутри трубки на значительном удалении от входного отверстия. Функция ф8(х,г) удовлетворяет уравнению (3) и граничным условиям
Течение газа, определяемое функцией ф3(х,г), формируется распределенными кольцевыми вихрями интенсивности ^(0 на поверхности трубки (£ — абсцисса положения кольца). В этом случае функция ф8(х,г) может быть выражена в ви-
где (З(.гуг) = (2/к — к)К(к)— 2Е(к)/к, к = 2^г/[х1 + (г + I)2], К (к) и Е(к) эл-
липтические интегралы первого и второго рода, которые нами находятся численно. Из условия (6) получается интегральное уравнение для определения ^(0
Уравнение (8) решается численно методом граничных элементов. После нахождения 7(£) поле скоростей (их,иг) определяется дифференцированием функции (4).
д2фр д2фр 1 дфр _
0 г,----1---7Г~П--------------о— — и
дх2 дг2 г дг
(3)
1 дфр 1 дфр
^X о 1 О •
г дг г дх
Функция фр(х, г) ищется в виде суммы двух функций
1. и ЦУр
г дг
(4)
ф8(х,0) = 0, —то < х < +то.
ф8(х, 1) = 0 < х < +оо.
(5)
(6)
де [И]
(7)
о
(8)
о
Рис. 2. Разделительные линии тока при различных Ra : 1 — модель потенциального течения. 2 модель вязкого течения
В приближении осесимметричного ламинарного течения несжимаемого газа описывается двумерными уравнениями Навье Стокса. В качестве расчетной области выбирается прямоугольная область в верхней части меридиональной плоскости. На левой и верхней границах расчетной области задается скорость невозмущенного внешнего потока. На правой границе выше трубки давление полагается равным нулю. В выходном сечении трубки задано параболическое распределение скорости по высоте, соответствующее развитому вязкому течению в канале круглого сечения со средней скоростью равной скорости аспирации. На внешней стенке н поверхности трубки принимаются условия прилипания, а на нижней границе расчетной области условия симметрии.
2. Результаты расчетов
Поле скоростей несущей среды в приближении вязкого течения несжимаемой жидкости рассчитывалось на основе численного решения уравнений Навье Стокса методом конечных объемов с помощью программы FLUENT. При численном интегрировании уравнений движения частиц (2) значения составляющих скорости газа находились интерполированием в поле (ux,ur), найденном из решения уравнений Навье Стокса. Вдали от трубки скорость частиц полагалась равной скорости газа.
В результате решения уравнений (2) находилась предельная траектория, разделяющая поток частиц на аспирируемые и проходящие мимо пробоотборника. Значение rc = Rc/Rt начальной радиальной координаты предельной траектории позволяет найти коэффициент аспирации по формуле A = r"^Ra.
Проведены исследования задачи аспирации в трубку с Rt = 0.005 м при скорости ветра U0 = 1 м/с. Скорость аспирации изменялась в диапазоне U0 £ [0.02,0.5] м/с. Соответствующее число Рейнольдса газового потока оказывается меньше 600, что обеспечивает ламинарный режим течения. Размеры расчетной области выбираются достаточно большие, чтобы дальнейшее увеличение области не влияло на решение. Сеточное разбиение расчетной области обеспечивает сходимость численного решения с точностью до величины 10-9.
На рис. 2 показаны разделительные линии тока газа, рассчитанные в рамках
Ra
чие разделительных линий тока незначительно и становится заметным при больших Ra. При Ra > 1 точка N разделения аспирируемого и уходящегося мимо пробоотборника потоков располагается с внутренней стороны трубки. С увеличе-
а б
Рис. 3. Линии тока (а) и траектории частиц (6) при Ка = 50 и St = 0.5 (штриховые линии разделительные лилии тока несущей среды, штрпх-пупктнрпая кривая предельная траектория)
Рис. 4. Зависимость Л(Ка) при различных числах Стокса: 1 — модель потенциального течения. 2 модель вязкого течения
нием Еа точка N удаляется от края внутрь трубки. При этом внутри трубки образуется возвратное течение, которое оказывает заметное влияние на процесс аспирации частиц при больших Еа.
Для учета влияния возвратного течения в качестве предельной выбирается траектория, попадающая в точку разделения газового потока N (рис. 3). Отметим также, что при Ка > 1 па внешней стенке трубки образуется вихревое течение.
Результаты параметрических расчетов коэффициента аспирации показаны на рис. 4-6. На рис. 4 приведены зависимости Л(Еа) для водяных капель с диаметрами d = 13, 34, 89.5 мкм (St = 0.051, 0.351, 2.43). Выбранные значения диаметров капель соответствуют эксперименту [9], результаты которого обозначены символами (о, □, Д). Расчетные кривые Л(Да), полученные в рамках двух моделей, близки друг к другу, то есть учет вязкости не приводит к заметному уточнению значений коэффициента аспирации. Наблюдается также удовлетворительное согласие расчетных и экспериментальных данных при St = 0.051 и 2.43. Для промежу-= 0.351
экспериментальных значений при увеличении Ка.
St
Рис. 5. Зависимость А(8^ при различных Ка: 1 — модель потенциального течения, 2 модель вязкого течения
Табл. 1
А — 1 + fi(^Ra — 1) Параметр (3 Работа
Аг 1 l + (2 + 0.62/fia)St. Беляев Левин (1974) [12]
А2 1 1 1 + St./(1 -I- 0.418/St.) Zliang. Liu (1989) [13]
Аз . 1 Paik. Vincent (2002) [9[
1 + (2 + 0.62/Ra - O.dRO-^St
Зависимости А (Я!;) при некоторых Еа даны на рис. 5. В целом расчетные и экспериментальные данные лучше согласуются при меньших Ка.
Представляется интересным сравнить результаты расчетов зависимостей А (Я!;) в области больших Ка с известными приближенными формулами для коэффициента аспирации тонкостенных пробоотборников. К настоящему времени предложены несколько формул для зависимости коэффициента аспирации от числа Стокса и параметра Еа , рекомендуемых для применения в области Еа > 1. Все формулы имеют общий вид
А = 1 + в(Яа — 1)- (9)
В табл. 1 приведены выражения для параметра в, предлагаемые в различных работах.
Наиболее известная формула Беляева Левина [12] хорошо описывает коэффициент аспирации в диапазоне 0.18 < Ка < 6. В [9] было показано, что формула Беляева Левина дает завышенные значения коэффициента аспирации при больших значениях Еа. На рис. 6 дается сравнение приближенных зависимостей коэффи-А А( )
модели течения, при Ка = 50. Наиболее близкими к результатам численной модели оказываются значения А, полученные то формулам для А^ и Аз. Значения коэффициента аспирации по формуле Беляева Левина оказываются завышенными. При меньших Ка значения коэффициента аспирации, найденные по формулам и из численной модели, сближаются. Таким образом, формулы из работ [9, 13] удовлетворительно описывают коэффициент аспирации в области больших значений Еа .
St
Рис. 6. Сравнения зависимостей A (St) при Ra = 50, полученных по приближенным формулам из табл. 1
Заключение
Развита математическая модель и проведены численные исследования коэффициента аспирации аэрозоля в тонкостенную трубку из движущегося газа при высоких значениях отношения скоростей ветра и аспирации. Получено удовлетворительное согласие расчетных данных с экспериментальными данными и приближенными формулами.
Автор выражает благодарность за оказанную помощь научному руководителю профессору Зарипову Шамилю Хузеевичу и профессору Маклакову Дмитрию Владимировичу.
Работа выполнена при финансовой поддержке РФФИ (проект X- 07-07-00183).
Summary
А.К. Gilfanov. The Calculation of Thin-Walled Sampler Performance in the Case of High Ratio of Wind and Aspiration Velocities.
The problem of aerosol sampling from a high-speed flow into a tliin-walled tube is considered. The velocity field of air flow is found using Boundary Element and Finite Volume methods within the potential and viscous flow of incompressible fluid models. The aspiration efficiency in the case of high ratio of wind and sampling velocities is studied. Curves of dependence of the aspiration efficiency on the Stokes number and velocities ratio are plotted. The comparison with experimental results and approximate formulas is given.
Key words: aerosol sampling, aspiration efficiency.
Литература
1. Vincent J. Aerosol sampling: science, standards, instrumentation and applications. John
Wiley&Sons, 2007. 616 p.
2. Grinshpun S., Lipatov G.N., Sutugin A.G. Sampling errors in cylindrical nozzles // Aerosol Sci. Techn. 1990. V. 12. P. 716 740.
3. Vituls V. Theoretical limits of errors due to anisokinetic sampling of particulate matter // J. Air Pollution Control Association. 1966. V. 16, No 2. P. 79 84.
4. Волощук B.M. Введение в гидродинамику грубодисперсных аэрозолей. Л.: Гидро-
метеоиздат, 1971. 208 с.
5. Rader D.J., Marple V.A. A study of the effects of anisokinetics sampling // Aerosol Sci. Techn. 1988. V. 8. P. 283 299.
6. Liu Zhang Z.Q., Kuehn Т.Н. A numerical study of inertial errors in anisokinet.ic
sampling // J. Aerosol Sci. 1989. V. 20, No 3. P. 367 380.
7. Medvedev A.A. Numerical study of aerosol particle sampling from a low-speed flow // Atmospheric and Oceanic Optics. 2002. V. 15, No 8. P. 663 666.
8. Kramer М., Afshine A. Sampling characteristics of inlets operated at low U/Uo ratios: new insights from computational fluids dynamics (CFX) modelling // J. Aerosol Sci. 2004. V. 35. P. 683 694.
9. Paik, S., Vincent, J.H. Aspiration efficiency for thin-walled nozzles facing the wind and for high velocity ratios // J. Aerosol Sci. 2002. V. 33. P. 705 720.
10. Dunnett S.J. A numerical study of the aspiration efficiency of a thin-walled sampler facing the wind for high velocity ratios // J. Aerosol Sci. 2005. V. 36. P. Ill 121.
11. Batchelor G.K., F.R.S. An introduction to fluid dynamics. Cambridge: Cambridge
Univ. Press, 2000. 631 p.
12. Беляев С.П., Левин, Л.М. Корректирование ошибок отбора проб аэрозолей // Изв. АН СССР. Физика атмосферы и океана. 1974. Т. 10, .V 5. С. 512 518.
13. Zhang Z.Q., Liu B.Y.H. On the empirical fitting equations for aspiration coefficients for thin-walled sampling probes // J. Aerosol Sci. 1989. V. 20. P. 713 720.
Поступила в редакцию 06.10.08
Гильфанов Артур Камилевич аспирант НИИ математики и механики им. Н.Г. Чеботарева Казанского государственного университета.
E-mail: artur.gilfanoveksu.ru
