Спектральный анализ подвижностей атмосферных частиц

ISSN 0868-5886

НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2014, том 24, № 4, c. 51-58

МАТЕМАТИЧЕСКИЕ МЕТОДЫ И МОДЕЛИРОВАНИЕ В ПРИБОРОСТРОЕНИИ

УДК 551.576 © А. Г. Варехов

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

В статье рассмотрен аналитический подход для оценки спектрального распределения подвижностей атмосферных частиц. За основу берется известная лабораторная практика анализа спектров с использованием вольтамперных характеристик аспирационных счетчиков. Поскольку при таких измерениях отсутствует статистический материал (выборка), делается целесообразное предположение относительно характера функции распределения. В качестве наиболее общей аналитической формы функции плотности вероятности выбрано гамма-распределение. Обоснованы и систематизированы представления о величине электрического заряда и подвижности частиц в диапазоне от 0.001 до 1.46 см2/(Вс). Предложены аналитические приемы, позволяющие определить параметры распределения на основе вольтамперной характеристики.

Кл. сл.: атмосферные частицы, заряд частиц, спектр подвижностей, гамма-распределение, оценка параметров

ВВЕДЕНИЕ

Инструментальный мониторинг атмосферных частиц на основе спектральных измерений размеров и подвижностей развивается в последнее время достаточно продуктивно. Интерес к изучению аэродисперсных систем быстро возрастает в связи с вредным воздействием частиц антропогенного аэрозоля на здоровье человека. Совместное изучение природных и антропогенных факторов, участвующих в образовании аэродисперсных систем, приводит в конечном итоге к расширению спектрального диапазона размеров частиц, который в большинстве случаев находится в пределах от 1 до 1000 нанометров.

Для нижней тропосферы над сушей крупных городов репрезентативным считается [1] так называемый континентальный тип одномодового (унимодального) распределения частиц по размерам, представленный на рис. 1. После резкого подъема до максимума около значения г = 0.01 мкм начинается медленный почти экспоненциальный спад плотности вероятности. Такое распределение соответствует хорошо известным [2] эмпирическим

формулам общего вида /

(а,Ь,р^ — эмпирические константы), справедливым для большинства конденсационных и дисперсионных аэрозолей.

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

ю5

ю3

§ УУ-

л

с Ii"

ю-51—

0.001 0.1 10 Радиус частицы, мкм

Рис. 1. Распределение по размерам частиц континентального аэрозоля

размеров, так и со стороны больших частиц, а модальный размер — передвигаться с течением времени в сторону больших частиц.

ПОСТАНОВКА ЗАДАЧИ

Большинство динамических свойств аэрозольной системы определяется подвижностью частиц, зависящей от их размеров и электрического заряда. Доля частиц, имеющих заряд q = ne (п — це-

лые положительные или отрицательные числа), в стационарном состоянии при симметричной биполярной ионизации описывается [2] больцма-новской функцией вида

Л\ -(пе)2/8п£0гквТ

п) = е • ' 0 в .

Задача изучения механизмов зарядки водных атмосферных частиц, поставленная еще А.И. Воейковым (1904), была последовательно решена Френкелем [3]. Исходное представление состояло в том, что плотная часть двойного электрического слоя на поверхности частицы состоит из диполь-ных гидроксилов, ориентированных отрицательными зарядами наружу частицы, а процесс заряжения обозначался как ориентированная адсорбция ионов на атмосферных ядрах. Оценка величины поверхностного потенциала частицы ф, располагалась в пределах 0.25 < ф, < 0.36 В, причем нижняя граница соответствовала скачку потенциала на двойном слое, а верхняя — значению, вычисленному исходя из соотношения для поверхностного потенциала ф,=па3рг / е0 (па, 3, рг — соответственно объемная концентрация диполей, толщина двойного слоя и проекция ди-польного момента на радиальное направление сферической частицы). Таким образом, поверхностный потенциал атмосферных водных капель был принят равным ф!, = 0.3 В. Это значение и сегодня [1] считается состоятельным. Радиус частицы, несущей один элементарный заряд, по Френкелю, был равен примерно 5 нм.

Поверхностный потенциал сферической частицы, а также ее заряд и радиус связаны кулонов-ским соотношением Фs = q / 4кее0г, откуда следует равенство q = 4кее0ф/ . Подставляя это равенство (при £ = 1) в выражение для функции / (п) , получаем аналогичное выражение в зависимости от радиуса частицы:

/

= ¿>~2п£оФ»г/квТ

г ) = е

При этом должно выполняться условие нормировки

2 м -2ж£оф1 .

2п£оф2 М е -ВТ" ¿г = 1,

КТ

где Х = 2ж£0ф1 / квТ — нормирующий множитель. Отыскание среднего

М

Г = ге~1г d г

приводит при ф^ = 0.3 В к противоречивому результату, т. е. к значениям среднего радиуса Г = 1/ Х = 0.827 -10~9м = 1 нм и дробного заряда q = 0.21 е . Физический смысл имеет минимальный заряд частицы q = 1е и соответственно радиус частицы гт1П = 4 -10~9м (4 нм). Таким образом, для интегрирования от начального значения г = гт1п необходимо перенормировать функцию / (г), т. е. определить дополнительный нормирующий множитель ам из соотношения

аыА | е кг d г

откуда следует ам = еАгт1п . Вычисления дают значения аы = 126.06 и средний радиус г = = 104.25 -10~9м = 0.1мкм . При этом средний заряд частиц составит д = 22е. Ограничение области интегрирования сверху, например значением г = 1мкм, не изменяет сколько-нибудь заметно результат.

Отметим другие приближения для оценки заряда частицы, определяющие предельный заряд: рэлеевское приближение поверхностного натяжения [4]; приближение начала спонтанной автоэлектронной эмиссии в электрическом поле у поверхности заряженной частицы [4]; ударную зарядку во внешнем электрическом поле [5]; диффузионную зарядку в среде с заданной концентрацией газовых ионов [5]. Два последних механизма зарядки характерны для технологии электрофильтрации, тогда как два первых описывают фундаментальные атмосферные процессы. Наконец, ориентированная адсорбция ионов на атмосферных ядрах дает значения заряда, обоснованные выше.

Для всех пяти механизмов зарядки величина предельного заряда в функции радиуса частиц представлена в табл. 1.

Как и следовало ожидать, диффузионная зарядка и ориентированная адсорбция ионов дают наиболее близкие закономерности по величине зарядов частиц. Отметим, что диффузионная зарядка описывается кинетически, т. е. содержит множитель, зависящий от времени. Напротив, представление об ориентированной адсорбции ионов основано только на свойствах сорбирующей поверхности и потому может рассматриваться как более предпочтительное.

Вычисления механической и электрической подвижностей производились в соответствии с общим выражением к = пеСс тг / 6пгц, в котором

1

г

Табл. 1. Величина предельного заряда частицы в единицах е = 1.610-19 Кл для различных механизмов зарядки

Механизм зарядки Радиус частицы, мкм

0.001 0.01 0.1 1.0

Рэлеевское приближение поверхностного натяжения 126 4103 1.3105 4106

Приближение начала спонтанной автоэлектронной эмиссии 1 70 7103 7105

Ударная зарядка в земном электрическом поле Е = 130 Вм1 0 0 0 0

Диффузионная зарядка 0 1 12 125

Ориентированная адсорбция ионов на атмосферных ядрах по Френкелю 1 1-2 22 220

Табл. 2. Результаты расчетов электрической подвижности

Радиус частицы, мкм Поправки Каннингэма

0.0005 156.275 (Ст) 0.9161015 1 1.4662

0.001 78.137 (Ст) 0.229 1015 1 0.3664

0.002 39.068 (Ст) 0.572-1014 1 0.0916

0.004 19.534 (Ст) 0.143 1014 1 0.0229

0.006 13.023 (Ст) 0.636 1013 1 0.0102

0.01 11.787 (Сс) 0.345 1013 1-2 0.0055

0.1 1.842 (Сс) 0.540 1011 22 0.0019

1.0 1.076 (Сс) 0.3151010 220 0.0011

5.0 1.015 (Сс) 0.595 109 1100 0.00105

зависящие от радиуса частиц поправки Каннингэ-ма Сс или молекулярно-кинетическая Ст вычислялись по формулам [2]

С,

1+ — г

1.142 + 0.558ехр

\ Хт /

и

Ст = (A + Q

г

где Хт — длина свободного пробега газовых молекул при нормальных условиях, Хт = 66.5 нм; сумма A + Q= 1.175 соответствует зеркальному отражению при столкновениях молекул и частиц; кроме того, динамическая вязкость п =

= 18.1-10-6 Па с. Результаты расчетов электрической подвижности приведены в табл. 2. и графически представлены на рис. 2.

Кнутсон и Уитби [6] описали в 1975 г. ставший основой для этих разработок электростатический

классификатор частиц по подвижности, который был предложен Чапменом [2] еще в 1937 г. Однако для преодоления технических трудностей уже в ранних вариантах таких устройств [2] требова-

0.0001 0,001 0,01 0.1 1

Радиус частицы, мкм

Рис. 2. Электрическая подвижность (см2/(В• с)) частиц в функции их радиуса (мкм) по данным табл. 2

лась предварительная симметричная зарядка частиц. Кроме того, необходимо было удалять из исследуемого воздуха легкие газовые ионы с подвижностью выше 1 см2 • В-1 • с-1, искажающие результаты измерений для относительно малоподвижных крупных частиц. Предварительная зарядка частиц при помощи ^-активного ионизирующего источника или мягкого (9.5 кэВ) рентгеновского излучения используется и в настоящее время, в частности в крупных центрах по исследованию атмосферных аэрозолей [7, 8]. В этих исследованиях электростатическая сепарация частиц используется для выделения достаточно узких фракций частиц для последующего анализа другими, чаще всего оптическими методами. Предпринимались одиночные попытки обойтись без предварительной зарядки [9].

Ясно, что эти экспериментальные приемы приводили к искажению зарядовой структуры исходного аэрозоля и к обеднению спектра. Совокупность этих причин, как отмечал уже на ранних этапах исследований Фукс [2], вела к тому, что во многих работах по физике высокодисперсных аэрозолей были допущены существенные ошибки.

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

МЕТОДИКА И РЕЗУЛЬТАТЫ

Рассмотрим цилиндрический аспирационный конденсатор, между электродами которого в осевом направлении протекает ламинарный поток анализируемого воздуха. Нетрудно показать, что для осаждения частицы (иона), имеющей начальную координату г>г1 и подвижность к (г), должно удовлетворяться неравенство

к ( г )>

1п ( г2/ г1 ) 2^и

( г2 - г12),

электроду. Распределение частиц по подвижно-стям характеризуется некоторой нормированной функцией / (к) плотности вероятности.

При определенном значении напряжения и для начальной радиальной координаты г' осаждаются ионы с подвижностью к > к (г') . Разобьем интервал г2 - г1 на некоторое не очень большое число одинаковых участков Аг, в пределах каждого из которых будем считать координату г постоянной. Рассмотрим для определенности деление межэлектродного расстояния на 10 равных участков. Число участков принципиального значения не имеет. Очевидно, из первого участка, наиболее удаленного от осадительного электрода, осаждаются ионы с подвижностью к > к (г'), из

второго — с подвижностью к > к (г '-Аг), из третьего — с подвижностью к > к(г' - 2Аг) и т. д. При заданном напряжении и осаждаются все ионы с подвижностью к > к (г'), 9/10 ионов с подвижностью к > к (г' - Аг) , 8/10 ионов с подвижностью к > к (г' - 2Аг) и т. д. Ток аспирационного конденсатора 1 (и) определяется суммарным зарядом частиц различного радиуса и заряда. Доля осаждающихся из каждого участка частиц некоторого радиуса и заряда определяется функцией плотности вероятности распределения ф(к), получающейся из функции / (к) после умножения ее значений на коэффициенты 1;9/10; 8/10 и т. д. Функция плотности вероятности ф (к) в отличие от общей функции / (к) может быть названа частной. При увеличении напряжения и измеряемый ток определяется возрастающей функцией ф (к) . При и ^ итах этот ток определяется функцией / (к) . При неизменности / (к) этот ток есть наибольший и равный току насыщения 1н = 1 (итах) .

Приведенные рассуждения позволяют получить следующее соотношение:

где г1, г2 — радиусы соответственно внутреннего (осадительного) и внешнего электродов конденсатора; ^ — время пребывания частицы в конденсаторе; и — приложенное напряжение. Для определенности предполагается, что внутренний электрод имеет положительный потенциал относительно нулевого потенциала внешнего электрода, а отрицательно заряженные частицы движутся по криволинейным траекториям к внутреннему

1н - 1 (и) = 1н

22

| / (к) d к/ (к) d к

22 г2 - г2

так что

22 г2 - г2

/ (к) ггЛ = ф (к).

г2 - г1

Кроме того, выполняется условие нормировки

функции /(к) :

со

| / ( к ) d к = 1.

При г = г2 и соответственно при к > к (г2) получаем равенство ф (к) = / (к) . Из выражения (2) имеем

1н - 1 (и) = 1н

к

| / ( к ) dk+ | / ( к ) d к-

к

к

2 2

* г — г я

| / (к ) -2Г-Г dk - | / (к ) dк

22 г - г

к

к

^^[к(г2)]- | /(к)^ак

22 г - г

22 г - г

- г2 к

перепишем (3) в виде

1н - 1 (и ) =1н ^ [к ( г2 )]-

0 <

< 1.

вольтамперная характеристика конденсатора, однако она не дает непосредственно информации о спектре подвижностей. Поскольку статистический материал (выборка) отсутствует, целесообразно сделать некоторые предположения относительно функции плотности вероятности.

В качестве наиболее общей аналитически формы распределения может быть выбрано двухпара-метрическое гамма-распределение

/

ва

Г

ка-1е-вк

где F [к (г2)] — соответствующее значение функции распределения F(к). Поскольку на основании (1) выполняется равенство

, к(г2) 1

[к

Наконец, удобно ввести безразмерный ток ас-пирационного конденсатора

где а > 0; / > 0; к > 0, причем форма кривой определяется параметром а, а параметр / носит характер масштабного множителя; Г(а) — гамма-функция Эйлера.

В зависимости от величины а функция / (к) видоизменяется от экспоненциального распределения (а = 1) до почти симметричного (нормального) распределения с ростом а. Отметим, что гамма-распределение хорошо соответствует эмпирическим формулам вида /(г) = агрехр(-Ьг*),

отмеченным выше. По этой причине эмпирические формулы этого вида иногда обозначают как обобщенное гамма-распределение [10].

На рис. 3 показано распределение подвижно-стей частиц, приблизительно соответствующее распределению частиц по размерам (рис. 1). Математическое ожидание и дисперсия подвижности даются в следующей простой форме:

к=а, в (к)=4. в* К} в2

Функция распределения F (к) дается выражением

Тогда уравнение, связывающее вольтамперную характеристику конденсатора с функцией плотности вероятности распределения подвижностей / (к ), запишется в виде

т1=1 [к

Нетрудно показать, что соотношение (5) удовлетворяет в интервале напряжений 0 < и < итах

обоим очевидным условиям 1 (0) = 0 и 1 (итах ) =

= 1н .

Единственным источником информации при измерениях аспирационными приборами является

F

1

:Г вк (а),

Г(а) вк

где Грк (а) — неполная гамма-функция, определяемая интегралом

Грк (а) = | ^-^ .

Следовательно, для вольтамперной характеристики с учетом равенства (7) получаем из уравнения (5) соотношение

'

1-

Г вк

к Г рк

Г

к

г

2

н

н

г

k

1

k (Г2 ) k (Г2 ) в

3

1 а

ктах в

из выражения (8), вводя аналогичные обозначения

1 - = ¿2; k ( Г2 ) = ¿0

получаем оценку для асимптотическои части характеристики:

3

Г Ар (а ) (а / в) Г А0 (а +1) Г(а) к0 Г(а +1) .

Таким образом, значения параметров распределения а ив могут быть получены при совместном решении уравнений системы

0,01 0.1 1 Подвижность, см2В-1с-1

Рис. 3. Функция плотности вероятности распределения подвижностей (гамма-распределение).

f (к) = 864к2ехр (-вк); а = 3; в = 12; к = 0.25;

о = \[В = 0.144; коэффициент вариации V = \/4а ; коэффициент асимметрии 2/-\/а= 1.15; эксцесс 6 / а = 2; максимум (мода) при значении подвижности к = (а -1) / в = 0.167

Используем теперь начальный (и ^ 0) и асимптотический (и ^ итах ) участки вольтамперной характеристики. При и ^ 0, т. е. при к (г2 ктах (для наибольшего значения подвижности; в частности, для распределения рис. 3 значения к = 2), из выражения (8) получаем

3

, Г вктх

1--+-•-,

Г

Г вк0

Г

к

Г

Первое из уравнений системы (12) приводится к виду (10). Вводя далее в (1) обозначение

( Г22 - Г12) Ь ^

2к 1

где Хс — коэффициент пропорциональности, зависящий только от параметров аспирационного конденсатора, получаем из (10) простое соотношение

3

кк

к

ти

тах с

Введем обозначение г (0) / 1н=31 — произвольное, но малое значение безразмерного тока; при этом к ( г2 ) = ктах — наибольшее значение подвижности. Поэтому перепишем (9) в виде

в котором напряжение и1 соответствует принятому малому уровню безразмерного тока. Таким образом, независимо от характера функции плотности вероятности f (к) математическое ожидание

подвижности определяется значением производной

d

dU I ^ У и=0

то есть

к=Х

Малая величина 31 отсчитывается непосредственно по вольтамперной характеристике; соответствующее значение подвижности определяется на основании (1) при напряжении и1 , соответствующем уровню 31 .

При и ^ итах (и 1н ) , т. е. при к (Г2 к0

dU

ни)

V У

Преобразуем теперь второе уравнение в (12). Для этого запишем его в форме

3

1

Г

1

Г вк0

н

3

2

а

н

d

=0

ßko

ßk0

1 Н*0 1 н^о

--—— f za-e^zdz--— f zV"zdz .

Г(a) J ßkor(а) 0

Произведем оценку величины /Зк0 . Для относительно крупных частиц (малых подвижностей) величина 2 в интегралах выражения (15) мало отличается от единицы для всех 2 из области интегрирования, и потому из (15) после вычислений приближенно получаем неразрешимое в элементарных функциях соотношение

S2 =

( ßko )а

Г(а + 2) .

i (U) 1 ßk

r = 1" F (k) + ^ f f (k)dk_

Тн ßk 0

_ 1 .Га!^) + Га

Г( а ) Г( а + 1)'

степенным рядом. Для так называемой регуляри-зированной формы неполной гамма-функции используем ряд [11]

(x,a)! _ ax+l ( x, да)! ( x + 1)!

f

1 - а 1!

x + 1 + а x + 2 2!

x + 1 x + 3

а

зГ

x + 1 x + 4

- + .

Таким образом, вообще говоря, необходимые параметры распределения а, в получаются при совместном решении уравнений (10) и (16). Трудность заключается в том, что асимптотический участок характеристики может быть сильно растянутым по причине слишком медленного роста ас-пирационного тока для крупных (малоподвижных) частиц. Для таких частиц (г > > 1 мкм; к< 0.001 см2/(В• с), см. табл. 2) и при условии а> 1 величина д2 становится чрезмерно малой, как это следует из анализа соотношения (16), и ее отсчет непосредственно по вольтампер-ной характеристике становится затруднительным.

Для определения параметров распределения достоверной точечной оценкой будем считать математическое ожидание подвижности, определяемое отношением к = а / в. Это значение может быть использовано для анализа распределения на основе экспериментальной вольтамперной характеристики. Для этого необходимо сначала определить значение приложенного напряжения и соответствующего осаждению на центральном электроде конденсатора всех частиц с подвижностью к(г2) = к . Далее необходимо в соответствии с (5) и при учете, как это следует из (7) и (8), записи подвижности в безразмерной форме (вк) получить безразмерный аспирационный ток в виде следующего соотношения:

Правая часть (17) является функцией только параметра распределения а , однако также неразрешима аналитически. Для оценки параметра можно воспользоваться представлением неполной гамма-функции достаточно быстро сходящимся

V ' ' - . J

где x = а; а + 1 для расчетов соответственно второго и третьего слагаемых в правой части соотношения (17). При значении параметра а = 2 и при учете первых десяти слагаемых разложения расчет дает значение безразмерного тока, равное i (U) / 1н = 1 - 0.595 + 0.233 = 0.638. Аналогично

для значения а = 3 получаем i (U) / IH = 1 -- 0.998 + 0.182 = 0.184. Эта вычислительная процедура может быть организована программным путем.

ЗАКЛЮЧЕНИЕ И ВЫВОДЫ

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

СПИСОК ЛИТЕРАТУРЫ

1. Ивлев Л.С., Довгалюк Ю.А. Физика атмосферных аэрозольных систем. СПб.: НИИХ СПбГУ, 1999. 194 с.

2. Фукс Н.А. Механика аэрозолей. М.: Изд. АН СССР, 1955. 114 с.

3. Френкель Я.И. Теория явлений атмосферного электричества. М.-Л.: ГИТТЛ, 1949. 155 с.

4. Райст П. Аэрозоли. Введение в теорию. М.: Изд. Мир, 1987. 280 с. (Reist P.C. Introduction to aerosol science. McMillan Publ. Company, NY, 1984).

5. Kirsch A.A., Zagnit'ko A.V. Diffusion charging of submicrometer aerosol particles by unipolar ions // J. Colloid Interface Science. 1981. Vol. 80, No. 1. P. 111-117.

6. Knutson E.O., Whitby K.T. Aerosol classification by electric mobility: Apparatus, theory, and applications // J. Aerosol Sci. 1975. No. 6. P. 443-451.

7. URL: (www.cas.manchester.ac.uk/restools/ instruments/aerosol/differential/schem_DMPS/ index.html.

8. URL: (www.iac.ethz.ch/edu/courses/master/ lab_fiel_work/atmospheric_physics_lab_work/ DMPS_en.pdf.

x

9. Laakso L., Gagne S., Petaja T. et al. Detecting charging state of ultra-fine particles: instrumental development and ambient measurements // Atmos. Chem. Phys. 2007. Vol. 7. P. 1333-1345.

10. Архипов В.А., Усанина А.С. Движение аэрозольных частиц в потоке. Томск: Изд. дом Томского гос. университета, 2013.

11. Янке Е., Эмде Ф. Таблицы функций с формулами и кривыми. М.: Гос. изд. физ.-мат. лит., 1959. 110 с.

Государственный университет аэрокосмического приборостроения, г. Санкт-Петербург

Контакты: Варехов Алексей Григорьевич, varekhov@mail. ru

Материал поступил в редакцию: 20.06.2014

SPECTRAL ANALYSIS MOBILITY OF ATMOSPHERIC PARTICLES

A. G. Varekhov

State University of Aerospace Instrumentation, Saint-Petersburg, RF

There is analytical approach for estimate of spectral distribution of mobility of atmospheric particles considered in the article. The known laboratory practice of the analysis of ranges with use of volt-ampere characteristics of aspiration counters is taken as a basis. Because at such measurements there is no statistical material (sampling), the expedient assumption concerning nature of function of distribution becomes. As the most general analytical form of function of density of probability gamma distribution is chosen. Values of electric charge and mobility of particles in the range from 0.001 to 1.46 cm2V-1s-1 are proved and systematized. The analytical receptions allowing to determine distribution parameters on the basis of the volt-ampere characteristic are proved.

Keywords: atmospheric particles, charge of particles, spectrum of mobility, gamma distribution, estimate of parameters

REFERENСES

1. Reist P.C. Introduction to aerosol science. McMillan Publ. Company, NY, 1984.

2. Kirsch A.A., Zagnit'ko A.V. Diffusion charging of submicrometer aerosol particles by unipolar ions. J.

Colloid Interface Science, 1981, vol. 80, no. 1, pp. 111-117.

3. Knutson E.O., Whitby K.T. Aerosol classification by electric mobility: Apparatus, theory, and applications. J. Aerosol Sci., 1975, no. 6, pp. 443-451.

4. URL: (www.cas.manchester.ac.uk/restools/

instruments/aerosol/differential/schem_DMPS/ index.html.

5. URL: (www.iac.ethz.ch/edu/courses/master/ lab_fiel_work/atmospheric_physics_lab_work/ DMPS_en.pdf.

6. Laakso L., Gagne S., Petaja T. et al. Detecting charging state of ultra-fine particles: instrumental development and ambient measurements. Atmos. Chem. Phys, 2007, vol. 7, pp. 1333-1345.

Contacts: Varekhov Aleksey Grigor'evich, varekhov@mail. ru

Article arrived in edition: 20.06.2014