CC BY
14
УДК 621.319.7.001.6
Стохастическая модель индукционной зарядки аэрозольных частиц
Семенов В.К., д-р техн. наук, Сорокин А.Ф., канд. техн. наук
Предлагается стохастическая математическая модель контактной зарядки аэрозольных частиц в аппаратах электронно-ионной технологии, основанная на уравнении Фоккера-Планка, позволяющая рассчитать не только средние заряды частиц, но их флуктуации.
Ключевые слова: индукционная зарядка, аэрозольная частица, электрод, электрическое поле.
Aerosol particle inductive charging operation probabilistic model
The probabilistic mathematical model of aerosol particles contact charging operation is suggested to be used in electron-ion technology operation. Based on formula of Foker-Planka, the model makes possible to calculate not only average particle charges, but also their fluctuations.
Keywords: induction charging operation, aerosol particle, electrode, electrical field.
Индукционная зарядка частиц используется во многих электротехнологических процессах -электросепарация, электропрядение, электроокраска, электроворсование и пр. Для удобства решения задачи допускаем, что частицы имеют форму полусферы или полуэллипсоида. При этом площадь контакта равна площади большого круга, а проводимость в месте контакта принимаем равной проводимости самой частицы. Однако оба эти предположения не соответствуют действительности: частицы имеют самую разнообразную форму, а переходное сопротивление в месте контакта частицы с электродом столь велико, что процесс зарядки частицы является квазистатическим, т.е. в каждый момент времени заряд, полученный частицей, распределяется по ней в соответствии с законами электростатики. Корректное решение задачи по расчету предельного заряда частицы приведено в [1]. Поскольку частицы могут иметь самую разнообразную форму, то довольно хорошей аппроксимацией их форм является эллипсоид вращения, изменением соотношения между полуосями которого можно изменять форму частицы в довольно широких пределах. В [1] показано, что заряд частицы, оторвавшейся от электрода, удовлетворяет следующему соотношению:
12 + 15k 2^ 4 + 15k 2сг
-a Eo < Я ^ , a2Eo,
(1)
1 + 32к2 " ° " 8 + 2к
где а - большая полуось частицы; Е0 - напряженность внешнего поля; к - коэффициент формы, зависящий от соотношения полуосей эллипсоида а и Ь.
Для вытянутого вдоль поля эллипсоида имеем [2]
1 e—, e Л-f*f.
3|n^ - 2e V faJ
1 - e
(2)
Для сплюснутого эллипсоида вращения с осью, перпендикулярной к вектору напряженности электрического поля, имеем
k=
2e3
Зл/Г+ё21^(1 + e2 )arctge - e
e =
(3)
Для определения зависимости заряда частицы от времени запишем уравнение кинетики зарядки
dNe
dt
■ = EaS,
(4)
где Е - напряженность электрического поля в месте контакта частицы с электродом; N - заряд частицы, измеряемый числом элементарных зарядов е; 5 - площадь контакта частицы с электродом; ст - электропроводность переходного сопротивления контакта.
После элементарных вычислений получим
dN
dt
где
N
1 -
Nm =
N
Nm
4neo (4 -
e(8 + 2k)
(5)
15k) 2 /—
a Eo - предельный заряд
4пєп (4 - к)
частицы; т = ——--------- - характерное время
(8 + 2к )оБ
зарядки частицы.
Интегрируя последнее уравнение, находим закон накопления заряда частицы во времени:
N = N„
1 - exp(—
(6)
Для электрофизических технологий представляет интерес не только величина зарядов, полученных частицами, но и флуктуации в процессе зарядки. В таком случае на процесс зарядки следует смотреть как на стохастический и исходить из вероятностных представлений. Поскольку величина заряда и время его приобретения в одном элементарном акте соответственно малы, по сравнению со средними значениями, определяющими макроскопический механизм процесса, стохастический процесс будет медленным. Для него вероятность
данного состояния определяется только вероятностью предшествующего и не зависит от предыстории процесса. Процессы, удовлетворяющие названному условию, называются марковскими. Как известно, непрерывный марковский процесс подчиняется уравнению Фокера-Планка для функции распределения /) [3]:
зг N, і) ді
0]+д_[S(N)f(N, I)], (7)
где A(N), B(N) - кинетические коэффициенты, соответственно представляющие собой среднее и среднеквадратичное изменение заряда частицы за единицу времени:
A(N) = |ю(^q)qdq, B(N,q) = 1 |ю(^q)q2dq. (8)
Коэффициенты A(N) и B(N) определяются вероятностью перехода т^, q), которая представляет собой средний поток зарядов на частицу, когда ее заряд фиксирован и равен N.
Для практики важно знать не сами функции распределения частиц по зарядам, а средние заряды и их флуктуации, которые можно найти методом моментов. Найдем уравнение для среднего заряда частицы <^. С этой целью умножим левую и правую части уравнения Фоккера-Планка на N и проинтегрируем по всевозможным значениям заряда. После несложных преобразований получим
б < N > бі
=< ) >.
(9)
Полученное уравнение имеет очевидный физический смысл: скорость изменения среднего заряда частицы определяется средним потоком элементарных зарядов. Однако чтобы им воспользоваться, нужно перейти от среднего потока к потоку от среднего заряда . Для этого разложим A(N) в ряд Тейлора вблизи <М>:
А^) * А(< N >) + А'(< N >)^- < N >) +
+1 А"(< N >)^-< N >)2.
Здесь штрихом обозначена производная по N. Усредняя по N, получим
< А(^) >= А(< N >) +1 А"(< N >)А. (10)
Так как зависимость A(N) линейна, то второе слагаемое, содержащее дисперсию А, исчезает:
d < N > .. .. . Nm (. < N >
--------= А(< N >) = т I 1--------
dt х I Nm
(11)
Именно это уравнение, совпадающее с (5), определяет кинетику зарядки частиц на детерминированном уровне описания.
Выведем теперь уравнение для дисперсии распределения:
А =<N-< N >)2 =<N2 >-< N >2. (12)
Из уравнения (11) имеем
б < N >2
М
= 2 < N > Л(< N >).
(13)
Умножая все члены уравнения (7) на N и интегрируя по N с учетом (12) и (13), получим
dА = 2(< ЫА{Ы > - < N >< А(М) >) + 2В(< N >).
Л
Раскладывая NA(N) и усредняя по Л,
имеем
< NA(N) >*< N > А(< N >) + 2 [А^)] N=<N> • А. (14)
Подставляя это разложение в предыдущую формулу, получим уравнение для дисперсии функции распределения:
— = 2А'(< N >) • А + 2В(< N >). (15)
dt
Так как поток заряда A(N) явно от времени не зависит, то вместо t введем новую переменную <М>, определяемую по уравнению (11):
dА Л ^ г, л2, «, 1 ^ В(< N >
--------= А-|!п А (< N >1 + 2—-------------------. (16)
d < N > d < N >^ -I А(< N >)
Полученное уравнение является линейным и интегрируется в квадратурах:
<ы >
А = 2 Л
В(< N >)
б < N >.
0 Л3(< N >) Подставляя сюда
. Nm (. <N >
Л(< N >) = —т 1-----------
т [ Nm
л, ч 1 ^ Л < N >
в(< N >)=—т і-------------
2 т I Nm
после интегрирования получим
(17)
(18)
А =<N>11 -
< N >
Nn
(19)
Максимума дисперсия достигает при <М> / Nm = 1/2, а дисперсия в максимуме -А = 0,25Nm (см. рисунок).
0,2
0,4 0,6 0,8 <Л>/Лт
Зависимость дисперсии распределения от заряда частицы
Список литературы
1. Волков В.Н. К теории зарядки частиц на электродах // Новые методы исследования в теоретической электротехнике и инженерной электрофизике: Межвуз. сб. науч. тр. / Иван. энерг. ин-т им В.И. Ленина. - Иваново, 1974. - Вып. 3. - С. 67-71.
2. Ландау Л.Д., Лифшиц Е.М. Электродинамика сплош- 3. Лифшиц Е.М., Питаевский Л.П. Физическая ки-
ных сред. - М.: Наука, 1982. нетика. - М.: Наука, 1979.
Семенов Владимир Константинович,
Ивановский государственный энергетический университет, доктор технических наук, профессор кафедры атомных электростанций, телефон (4932) 38-57-78, e-mail: npp@aes.ispu.ru
Сорокин Александр Федорович,
Ивановский государственный энергетический университет, кандидат технических наук, доцент, профессор кафедры электрических систем, телефон (4932) 41-60-10 e-mail: deaneef@eef.ispu.ru
