2012
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА Математика и механика
№ 1(17)
МЕХАНИКА
УДК 531.351
М.А. Бубенчиков
ДВИЖЕНИЕ ЧАСТИЦ КСЕНОНА В ЦИКЛОННОЙ КАМЕРЕ
В работе описана вычислительная технология расчета среднего (без учета броуновских флуктуаций) перемещения атомов тяжелого инертного газа ксенона в воздушной среде. Проведены исследования движения таких частиц в циклонной камере. Найдены режимы, позволяющие улавливать указанные частицы при достаточно больших расходных характеристиках аэродинамической центрифуги.
Ключевые слова: ксенон, наночастица, циклонная камера, аэродинамика, движение частиц, режимы улавливания.
Ксенон - простое вещество - инертный одноатомный газ без цвета, вкуса и запаха. Ксенон - весьма редий элемент. При нормальных условиях 1000 м3 воздуха содержат около 87 см3 ксенона. Диаметр атома ксенона ё = 0,44 нм, а атомный вес М = 131 у.е. Ксенон используют для наполнения ламп накаливания, а также мощных газоразрядных и импульсных источников света (высокая атомная масса газа в колбах ламп препятствует испарению вольфрама с поверхности нити накаливания). Ксенон, как в чистом виде, так и с небольшой добавкой паров цезия, является высокоэффективным рабочим телом для электрореактивных (главным образом, ионных и плазменных) двигателей космических аппаратов. В 1999 году ксенон был разрешен к медицинскому применению в качестве средства для общего ингаляционного наркоза. В наши дни ксенон проходит апробацию в лечении зависимых состояний. Ксенон получают как побочный продукт производства жидкого кислорода на металлургических предприятиях.
Цель настоящей работы - оценить возможность получения ксенона в лабораторных условиях методом центрифугирования.. Для того чтобы решить задачу, необходимо знать сопротивление таких частиц в рассматриваемой среде и уметь рассчитывать аэродинамику циклонной камеры при высокой скорости вращения камеры: п ~ 6000-12000 об/мин. Вопросы сопротивления компактных наночастиц решены в работах [1, 2], а вычислительная технология, используемая для расчета аэродинамики при высоких частотных характеристиках циклонной камеры, представлена в [3]. Поэтому перейдем непосредственно к движению наночастиц (атомов ксенона) в циклонном аппарате.
Физическая область течения
Так как течение является осесимметричным и независимых переменных для расчета аэродинамики только две - г и г, то можно рассматривать любое осевое сечение цилиндрической камеры, а область интегрирования в этом случае будет плоской фигурой, в нашем примере прямоугольником.
На рис. 1 показана область изменения независимых параметров, являющаяся продольным сечением камеры. Здесь сплошные линии представляют непроницаемые границы, пунктирные - проницаемые и условную ось симметрии.
Л
Численное решение задачи динамики частиц
Наночастица участвует в сложном движении. Во-первых, она переносится потоком газа и, во-вторых, перемещается относительно потока несущей среды под действием сил инерции. Поэтому по теореме сложения скоростей для проекций абсолютной скорости частицы можем записать
ёг ё*
г1к = W{{*^г{*)) + ™{*).
(1)
(2)
(3)
Здесь и, У, V - проекции скорости газа на оси цилиндрических координат; и, V, w
- проекции скоростей относительного движения частиц на те же оси.
При этом относительное движение частицы определяется следующими уравнениями:
=а7{{* ^г {*))-2ви {/);
V2 ({*),г{*))
— = аг ( {*^ г {*))+— - 2Pv {*);
ш г
dw
= ^ {{*),г{*))-2Pw{*),
(4)
(5)
(6)
где
а,7 =-У-Уи , аг = -У-УУ , а,„ =-У ,
V = (и,у,ж), у = Г—+] —+
дг дг г дф
Здесь а,, аг, а^ - проекции переносного ускорения; V2/г - центробежное ускорение; 2ри, 2pv, 2pw - проекции ускорения (замедления), обусловленного действием силы сопротивления среды на частицу.
При * = 0 частицы стартуют с позиции во входном сечении (г е [0, Ьвх ], г = г1) с нулевой относительной скоростью по отношению к несущей среде:
Численное интегрирование системы уравнений (1) - (6) с начальными условиями (7), (8) будем проводить с использованием схемы Рунге-Кутта четвертого порядка точности и неявного этапа для промежуточных значений и, V, w, найденных по уравнениям (4) - (6). Из-за осесимметричного характера течения несущей среды в приводимой ниже численной модели движения частиц будут отсутствовать окружные перемещения, поскольку они не будут использованы в дальнейшем физическом анализе. Согласно подходу Рунге-Кутта, опирающемуся на идею пересчета, для того чтобы вычислить значения искомых величин на новом слое по
времени (в момент *к+1), необходимо предварительно найти правые части выпи-
Г к к+1 1
санных уравнений в четырех точках интервала * е I * , * I или, что то же самое,
в четырех позициях на фрагменте траектории, отвечающем указанному интервалу времени:
Начальные условия
г(0) = z0, г(0) = ^ ф(0) = ф0 ; и (0) = 0, V(0) = 0, w (0) = 0.
(7)
(8)
Схема вычислений
= гк, г1 = гк, и1 = ик, V = Vй, и = и(],г1), V = V(гх,г1); (9)
и2 = ----------------------------------------- , V-, = ------------------------------------------------------------------- ;
2 (1 + рИ) 2 (1 + рИ)
г3 = Г + 2 (и 2 + и2 ) г3 = гк + 2 ( + ^ ), и3 = и (^ г3 ), 13 = 1 (^ г3 );
(11)
(1+ РИ)
(12)
(13)
(1+ РИ)
¿4 = гк + И (и3 + и3), г4 = гк + И V + v3), и4 = и (4, г4), І4 = V (4, г4); (14)
(1 + 2РИ)
(1 + 2РИ)
Здесь ю^, г )=— ’ - угловая скорость квазитвердого вращения сплошной
г
среды, И - шаг по времени.
Тогда значения искомых величин на новом слое по времени найдутся по формулам
/+1 = 2к + И ( + м) + 2 (и 2 + и2) + 2 (и3 + и3 ) + (и4 + и4 ); (16)
6
гк+1 = гк + И((V + VI) + 2 (2 + У2) + 2(Уз + Уз) + (V. + У4 ) ; (17)
6
и = — (1 + 2^2 + 2^з + м.); (18)
6
Ук+1 =1 ( + 2у2 + 2у3 + у4 ). (19)
6
Соотношения (16) - (19) представляют собой формулы явного определения цилиндрических координат и скорости частицы по технологии Рунге-Кутта. При этом промежуточные значения относительных скоростей частицы найдены с использованием неявного представления силы сопротивления в разностных уравнениях (4) - (6). Полунеявный этап расчетов выражен соотношениями (11), (13), (15).
Проводя последовательно расчеты по формулам (9) - (19), находим численное решение задачи о движении частицы в закрученном потоке, полученное с использованием построенной явно-неявной схемы. Нетрудно видеть, что для систем без трения (в = 0) технология (9) - (19) превращается в классический вариант системы Рунге-Кутта четвертого порядка точности.
Интерполяция
Частицы, двигаясь по своим траекториям, разве только случайно могут пройти через узел разностной сетки, используемой для решения уравнений аэродинамики. Поэтому необходима интерполяция значений скорости и переносного ускорения в каждую из промежуточных позиций на траектории, необходимых для расчета искомой величины на новом слое по времени по технологии Рунге-Кутта.
Рис. 2. Интерполяция во внутреннюю точку по четырем значениям прямоугольной разностной ячейки
Если точка с координатами &, гх) попадает в указанную на рисунке ячейку, то интерполяция проводится по следующим линейным зависимостям:
и г О' )) = и ( (), г (}))+ ) (и ( (^г ())-и ((), г (})))^
2 (/ + 1)-2 (/)
и (( г и +1)) = и ((), г (} +1))+ (2+1)г(г ) (и И* +1), г ( +1))-и и(' X г (} + Ш
г(/ +1)-г(/) и(^г1) = иги)) + (Г1+1)(]\ ■)(и(иг(^+1))-игО'))). г и +1)-г и)
Аналогичным образом проводится интерполяция для У, Ш (либо ю), аг, а2 в точки
(22, г2), (23, г3), (24, г4).
Результаты расчетов
Во всех представленных ниже вариантах движутся с одной угловой скоростью либо две границы Г2 и Г4, либо они же совместно с границами Г1 и Г3 (торцевыми стенками камеры). Во всех вариантах скорость радиальной подачи воздуха Увх = 1 м/с. Предварительно был рассчитан вариант прямоточного течения (п = 0), для которого траектории частиц ксенона в точности повторили линии тока несущей среды. Ниже представлены различные варианты закрученных течений, отличающиеся степенью и характером закрутки. Как видим из рис. 3 - 5, существуют принципиально различные режимы движения наночастиц в циклонной камере. Первый режим (рис. 3) без вращения торцевых стенок камеры, имеющий умеренные частотные характеристики (п = 6000 об/мин) и следующие размеры: г1 = 0,4 м, г2 = 0,405 м, Ьвх = Ьвых = 0,1Ь, Ь = 0,1 м (длина камеры), характеризуется выносом всех поступающих в камеру частиц. Второй режим (рис. 4) характеризуется следующими геометрическими параметрами: г1 = 0,3 м, г2 = 0,35 м,
Ьвх = Ьвых = 0,1Ь, Ь = 0,4 м и в большей степени определяется вращением торцевых стенок камеры. В этом случае радиально входящие в камеру частицы сразу же устремляются в область пониженного давления, находящуюся у левой торцевой стенки. Попадая в конечном счете на эту стенку, они продолжают испытывать влияние действующей на них центробежной силы и медленно перемещаются в угловую зону камеры, ограниченную левой стенкой и внешней цилиндрической поверхностью аппарата. Третий режим (рис. 5, 6) определяется следующими параметрами камеры: г1 = 0,4 м, г2 = 0,5 м, Ьвх = Ьвых = 0,1Ь, Ь = 0,4 м. Существование этого режима в большей степени обязано появлению в потоке внутренних сингулярных зон высокой радиальной скорости несущей среды (на рис. 6 это зоны сгущения линий тока), которые ограничивают движение частиц в аксиальном направлении. Из-за специфического влияния этих зон на характер движения наночастиц их целесообразно назвать «шторками». Правая шторка является более протяженной и полностью перекрывает частицам путь к выходному сечению. Поэтому относительно более тяжелые атомы ксенона (в сравнении с молекулами Ы2 и 02) вынуждены циркулировать в ограниченном объеме камеры, выполняющем роль пылесборника. Аналогичный эффект наблюдается в циклонных пылесосах, не имеющих мешков для сбора мусора.
Рис. 3. Траектории частиц ксенона в режиме без вращения торцевых стенок камеры: п = 6000 об/мин
Рис. 4. Траектории частиц ксенона в режиме с вращением торцевых стенок камеры: п = 6000 об/мин
д Л
0,49-
0,48 - \ «ШШВВнН
0,47- 0,46-
0,45 - ті ВіиАІ
0,44- °'44 імИ іі 1
0,43- ИР І II о ¡¡І
0,42- 0,42 ~ ЖшЩ 1
0,41 -
0 4Л- пал--'''' — ~ ^ ^
Рис. 5. Траектории частиц ксенона в режиме без вращения торцевых стенок камеры: п = 12000 об/мин
Рис. 6. Линии тока несущей среды. Тот же режим, что и на рис. 5
Заключение
В работе описана вычислительная технология расчета динамики наночастиц, испытывающих значительное удельное сопротивление (приходящееся на единицу массы частицы). Найдены режим быстрого выхода частиц на вращающуюся торцевую стенку аппарата и режим внутренней циркуляции наночастиц в рабочей зоне устройства (с последующим выходом на внешнюю стенку вращающейся камеры). Оба режима могут быть использованы для улавливания наночастиц в центробежном аппарате. Эффект улавливания в двух найденных режимах локализации ксенона достигается за счет совместного действия механизма конвективного
переноса инерционных частиц и их диффузионного перемещения под действием центробежной силы.
ЛИТЕРАТУРА
1. Бубенчиков М.А. Механическое сопротивление компактных наночастиц в воздушной среде // Изв. вузов. Физика. 2011. № 1. С. 92-96.
2. Потекаев А.И., Бубенчиков М.А. Седиментация наночастиц в поле центробежных сил // Изв. вузов. Физика. 2011. № 2. С. 37-42.
3. БубенчиковМ.А., Иванова И.А. Расчет аэродинамики циклонной камеры // Вестник Томского государственного университета. Математика и механика. 2011. № 1(13). С. 46-51.
Статья поступила 01.03.2011г.
Bubenchikov M.A. MOTION OF XENON PARTICLES IN A CYCLONE CHAMBER. In this work, the computing calculation technology of an average (without regard to Brown fluctuations) motion of heavy inert xenon atoms in an air environment is described. Researches of the motion of such particles in a cyclonic chamber are carried out. Modes that make it possible to catch the particles at enough big account characteristics of an aerodynamic centrifuge are found.
Keywords: xenon, nanoparticle, cyclonic chamber, aerodynamics, motion of particles, catching modes.
BUBENCHIKOV Mikhail Alekseevich (Tomsk State University)
E-mail: michael121@mail.ru