УДК 535.361 535.43.07
ЭКСПЕРИМЕНТАЛЬНАЯ УСТАНОВКА ДЛЯ ОПРЕДЕЛЕНИЯ ДИСПЕРСНЫХ ПАРАМЕТРОВ АЭРОЗОЛЬНЫХ СУБМИКРОННЫХ
СРЕД
Э.А. Мецлер, С.С. Титов, Е.В. Муравлев, А.А. Павленко, В.А. Архипов
Рассмотрена новая экспериментальная установка для определения концентрации, среднего размера частиц и функции распределения частиц по размерам аэрозольной среды в динамическом режиме. Представлена методика определения дисперсных характеристик аэрозольных сред. Приведены результаты экспериментального исследования параметров аэрозольной среды.
Ключевые слова: аэрозоль, обратная задача оптики аэрозолей, размер частиц, функция распределения частиц по размерам, спектральная оптическая плотность
ВВЕДЕНИЕ
Работа направлена на развитие и модификацию турбидиметрического метода (метода спектральной прозрачности [1]) путем перехода от использования непрерывного спектра к набору дискретных длин волн зондирующего излучения. При замене источника излучения с непрерывным спектром на набор монохроматических источников излучения (лазеров) необходима замена приемника излучения (спектрометра) на набор фотодиодов с системой селекции. К преимуществам разработанного метода относятся:
- возможность диагностики дисперсных сред с увеличенной (до нескольких метров) длиной оптического пути за счет использования коллимированного излучения с большой пространственной плотностью энергии;
- возможность исследования оптически плотных сред с обеспечением выполнимости закона Бугера-Ламберта-Бера;
- повышение точности измерения интенсивности зондирующего излучения за счет увеличения отношения сигнал/шум при использовании фотодиодов в качестве приемников излучения;
- компактность и мобильность измерительной установки.
ОПРЕДЕЛЕНИЕ ДИСПЕРСНЫХ ПАРАМЕТРОВ АЭРОЗОЛЬНОЙ СРЕДЫ
Восстановление функции
распределения частиц по размерам при дистанционной диагностике дисперсных сред по экспериментальным данным об
ослаблении оптического излучения этой средой осуществляется путем решения обратной задачи оптики аэрозолей. В качестве основного уравнения для решения данной обратной задачи оптики аэрозолей используется уравнение Фредгольма первого рода [2]:
/(!) = 10 (!) ехр[-г(А)], (1)
nCJ
\Q (ZD, m (Я)! D2f (D ) dD
где т(Я) =
где I(!) - поток излучения, прошедший через дисперсную среду на длине волны ! ; 10 (!) - поток излучения, падающий на дисперсную среду на длине волны ! ; т(!) -спектральная оптическая плотность на длине волны ! ; С„ - счетная концентрация частиц; I - длина оптического пути; < - фактор эффективности ослабления; о - диаметр частицы; ш(Х) = (п - гк) - комплексный
показатель преломления материала частиц ( п - показатель преломления; к - показатель поглощения); / (о) - функция счетного
распределения частиц по размерам.
Искомая функция распределения частиц по размерам априорно задается в виде логарифмически нормального распределения частиц по размерам (ЛНР):
(1п о - 1п о)
f ( D ) =
Djlä 1 ист
exp
2 (lna)2
(2)
где 1ист - логарифм среднеквадра-тического отклонения диаметров частиц;
4
о
lnD - логарифм математического ожидание диаметров частиц.
Для решения обратной задачи оптики аэрозоля используется вариационный подход к построению регуляризирующего алгоритма А.Н. Тихонова [3]. При этом на рассматриваемое множество функций распределения частиц по размерам {f} (
f=f (D)) накладывается ограничение,
используя априорную информацию о среднем объемно-поверхностном диаметре ( D32). Задача сводится к поиску минимума
функционала Тихонова minM%[f] на
множестве G , где G - замкнутое множество, % - параметр регуляризации. В соответствии
с алгоритмом решения некорректно поставленных задач осуществляется минимизация функционала Тихонова [3]:
M %[f ]= J
CJ Jß [f,m(X)) f (D) D2dD'2
d X-
(3)
+U
0
dD
f 2 (D)-( f'(D ))2 _
ts (X) - экспериментально определенная
спектральная оптическая плотность исследуемой дисперсной среды с
погрешностью S ; Xn, Xx - минимальная и максимальная длина волны зондирующего излучения.
Функция распределения частиц по размерам f (D), которая соответствует минимуму функционала Тихонова minMf], является решением обратной
f^G L J
задачи.
Однако в предложенной модификации турбидиметрического метода определения дисперсных характеристик при переходе на ограниченное число длин волн зондирующего излучения возникает неопределенность при восстановлении функции распределения частиц по размерам, связанная с привлечением меньшего количества измерительных данных об исследуемой среде. Уменьшение этой неопределенности достигается промежуточным этапом, который заключается в определении осредненного размера исследуемой дисперсной среды D32. Уменьшение неопределенности реализуется вследствие привлечения дополнительной информации, которая характеризует особенности искомой функции
распределения частиц по размерам. Также осредненный размер () функционально связан с параметрами распределения частиц по размерам. Это позволяет конкретизировать решение поставленной задачи, то есть ограничить область решений обратной задачи оптики аэрозолей. Средний объемно-поверхностный диаметр частиц для ЛНР по размерам определяется выражением [2]
т: 5.
ln d32 = ln D + - (!пст)
(5)
Средний объемно-поверхностный
диаметр частиц, определяемый выражением
D32 =
JP'f (D)dD
г о
\D2f (D)dD
(6)
Этот размер широко используется для описания дисперсных систем.
Определение среднего объемно-поверхностного (Б32) осуществляется с использованием модифицированного метода спектральной прозрачности (ММСП) [5], который основан на измерении спектральной оптической плотности аэрозольной среды для двух фиксированных длин волн зондирующего лазерного излучения. Основным уравнением метода определения является выражение для спектральной оптической плотности слоя равномерно распределенных полидисперсных частиц:
3ст1 ё (а, рз2)
2РВ32
где Сш - массовая концентрация частиц; ё (А, Д2) - усредненный фактор эффективности ослабления излучения; р -плотность материала частиц.
Усредненный фактор эффективности ослабления ё является важнейшей характеристикой, определяющей оптические свойства совокупности полидисперсных частиц, определяется выражением
- i°ß( ß X D32 ) = ^
tTD/
7/l
m
(X)) f (D )
D2dD
(9)
) D2dD
! / (с)'
Как следует из уравнения (1), отношение для двух соответствующих длин волн А и является однозначной функцией В32:
ё (А2, Оз2 )
X2
Т, ß (4, D32)
= F (D32).
(10)
2
X
Т
л
А
2
Путем математического моделирования установлено, что функция Р (Б32) зависит от
ширины функции распределения частиц по размерам (а), это может приводить к возникновению методической погрешности в определении Б32 . На рисунке 1 приведены
зависимости Р(Б32) для ЛНР с а = 1,2 и а
= 1,6 для длин волн 3,39 мкм и 0,405 мкм
(водный аэрозоль). Оптические константы для частиц водного аэрозоля взяты из [5].
(М
со
С) 1,6
1,4
1,2
0,8
0,6
0,4
0,2
1 J
Тд2 V / /
1 Д032 2
л
—« L
0,5
1,5
2,5
3,5
4.5
D
5,5 32
Рисунок 1 - Зависимости Р (Щ2) для
ЛНР с различной шириной распределения 1 - ЛНР для а = 1,2; 2 - ЛНР для а = 1,6; 3 - среднее
Таким образом, используя информацию об ослаблении оптического излучения аэрозольной среды решается обратная задача оптики аэрозоля по восстановлению функции распределения частиц по размерам. Область решения обратной задачи ограничивается диапазонам Д032. Далее осуществляется минимизация функционала Тихонова. Выбирается / ( Б)
соответствующая положению минимума. Затем по известным параметрам / (Б) рассчитывается Б32 и Сш .
ЭКСПЕРИМЕНТАЛЬНАЯ УСТАНОВКА
Структурная схема экспериментальной установки (рисунок 2) состоит из семи источников лазерного излучения (1): гелий-неоновый лазер Я = 3,39 мкм и Я = 0,632 мкм, полупроводниковые лазеры Я = 0,905 мкм,
Х- 0,850 мкм и Я = 0,780 мкм, диодные лазеры Я- 0,532 мкм и Я- 0,405 мкм. Сведение лазерного излучения один световой пучок в диапазоне длин волн от 0,405 мкм до 0,905 мкм происходит с использованием блока дихроических зеркал. Затем объединение светового пучка с лазерным излучением Я- 3,39 мкм
происходит с использованием
плоскопараллельной пластины из германия (Ge).
Для гелий-неонового лазера Я- 3,39 мкм осуществляется контроль мощности, так как у этого источника во времени происходит изменение мощности излучения. Канал контроля мощности гелий-неонового лазера Я - 3,39 мкм состоит из светоделительной плоскопараллельной пластины из кремния (Si) и фотоприемника. Лазерное излучение, пройдя диафрагму (3), попадает в измерительный объем (4), где ослабляется аэрозольной средой. Далее ослабленное лазерное излучение попадает в фотоприемный блок, где происходит разделение лазерного пучка на исходные составляющие с использованием
дифракционной решетки (5).
Лазерное излучение, ослабленное облаком аэрозоля, регистрируется фотоприемниками отдельно для каждой длины волны. Для регистрации лазерного излучения в диапазоне длин волн от Я -0,405 мкм до Я-0,905 мкм используются кремниевые фотоприемники ФДУК-11. Для регистрации лазерного излучения с Я-3,39 мкм используются фотоприемные устройства AMP48 (компании АИБИ) на основе гетероструктуры
GaSb/GalnAsSb/GaAlAsSb с однокаскадным термоэлектрическим модулем. AMP48 предназначены для регистрации оптического излучения в диапазоне от 0,8 мкм до 4,9 мкм с частотой от 150 Гц до 350 кГц. Механический обтюратор (2) производит модуляцию оптического излучения с частотой 480 Гц и скважностью 2. Регистрация сигналов с фотоприемников осуществляется аналого-цифровым преобразователем (АЦП) L-502 (фирмы L - Card), имеющим 16 дифференциальных входов с частотой дискретизации до 400 кГц.
На рисунке 3 приведен внешний вид экспериментальной установки.
Рисунок 2 - Структурная схема измерительной установки 1 - лазерные источники излучения; 2 - механический модулятор света; 3 - диафрагмы; 4 - измерительный объем; 5 - дифракционная решетка; Si - плоскопараллельная пластина из кремния; Ge - плоскопараллельная пластина из германия
Рисунок 3 - Внешний вид экспериментальной установки
ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ
Проведена серия экспериментов по определению динамики аэрозольной
среды, при этом первоначально для каждого момента времени был рассчитана функция распределения частиц по размерам. Далее по известному значению Д32 была рассчитана массовая концентрации
частиц Ст. Определение дисперсных характеристик аэрозольной среды осуществлялось при заполнении
измерительного объема облаком аэрозоля, создаваемым распылителем «Муссон-2». В качестве модельной жидкости для создания аэрозоля использовалась дистиллированная вода. Генерация аэрозоля осуществлялась в течении 5 с. Измерительный объем имеет размеры 8х10х16 см. Окна измерительного
объема изготовлены из лейкосапфира (АЮэ), который имеет спектральную прозрачность в диапазоне длин волн зондирующего излучения (0,3 мкм - 6 мкм). На рисунке 4 приведена типичная динамика Б32 аэрозоля в измерительном объеме.
<М 3,5 ГО
О 3
1 2 3 4 5
I
1 1
\
1
60 С
Рисунок 4 - Динамика Б32 водного аэрозоля в замкнутом объеме
На рисунке 5 приведено изменение массовой концентрации частиц С .
Рисунок 5 - Динамика концентрации водного аэрозоля
Эволюция аэрозоля определяется совокупным воздействием трех основных факторов: коагуляция, испарение и гравитационное осаждение. Совокупность этих основных факторов определяет динамику концентрации и размеров частиц аэрозольной среды [6].
Внешний вид зависимостей,
приведенных на рисунке 4 и 5, позволяет сделать следующие выводы. Разброс значений Б32 до 10 с (рисунок 4) объясняется нестационарным образованием облака аэрозоля. Далее наблюдается
тенденция плавного снижения размеров ВЪ2, соответствующая экспоненциальному
уменьшению концентрации частиц. Для объяснения полученной зависимости БЪ2{г) была проведена оценка влияния гравитационного осаждения частиц. В стационарном режиме скорость осаждения частицы и определяется уравнением [7]:
и = (Р^ з
18ц (11)
где ра - плотность воздуха; g - ускорение свободного падения; ц - коэффициент динамической вязкости воздуха.
На рисунке 6 приведена зависимость скорости гравитационного осаждения от диаметра частиц.
Рисунок 6 - Зависимость скорости гравитационного осаждения от диаметра частиц
Время т похождения частицей длины пути ь при гравитационном осаждения определяется уравнением:
т=Ь.
и
(12)
На рисунке 7 приведена зависимость времени прохождения пути ь =12 см (так как измерение дисперсности производилось на высоте 4 см от дна измерительного объема) от диаметра частиц.
Рисунок 7 - Зависимость времени прохождения пути 12 см при осаждении частиц от их диаметра
Анализ приведенных результатов показал, что за время измерения мелкодисперсная фракция в <4 мкм практически не осаждается. Уменьшение размеров частиц при г >50 а можно объяснить процессом испарения частиц.
Для моментов времени (1-5) выделенных на рисунке 4 приведены функции распределения частиц по размерам (рисунок 6).
Рисунок 6 - Функции распределения частиц
по размерам 1 - 2,3 с; 2 - 8,3 с; 3 - 24,3 с; 4 - 55,3 с; 5 -59,3 с
Внешний вид зависимостей,
приведенных на рисунке 6, демонстрирует деформацию функции распределения частиц по размерам в сторону меньших диаметров в течении времени эксперимента.
ВЫВОДЫ
• Разработана экспериментальная лазерная установка для определения
концентрации, среднего размера частиц и функции распределения частиц по размерам аэрозольной среды в динамическом режиме.
• Разработан математический аппарат решения обратной задачи оптики аэрозолей, позволяющий определять функцию распределения частиц по размерам с использованием априорной информации об экспериментально измеренном значении среднего объемно-поверхностного диаметра частиц исследуемой среды.
• Проведена серия экспериментальных исследований по определению динамики среднего объемно-поверхностного диаметра частиц аэрозольной среды, также для каждого момента времени определена функция распределения частиц по размерам и массовая концентрация частиц. Показана возможность исследования эволюции различных аэрозольных сред с использованием разработанной установки.
• Результаты измерений позволяют судить о различных процессах эволюции аэрозоля, таких как коагуляция, испарение, гравитационное осаждение частиц.
Работа выполнена при финансовой поддержке РФФИ (грант мол_нр 16-3250111).
СПИСОК ЛИТЕРАТУРЫ
1. Титов С.С., Павленко А.А., Кудряшова О.Б., Архипов В.А., Бондарчук С.С., Турбидиметрический метод определения параметров субмикронных аэрозольных сред // Автометрия. - 2012. - Т. 48, № 3. - С. 68-74.
2. Архипов В.А., Бондарчук С.С. Оптические методы диагностики гетерогенной плазмы продуктов сгорания: учеб. пособие. - Томск: Изд-во Том. ун-та, 2012. - 265 с.
3. Ягола, А. Г. Обратные задачи и методы их решения. Приложения к геофизике [Электронный ресурс] / А. Г. Ягола, Ван Янфей, И. Э. Степанова и др. - 2-е изд. (эл.). -М. БИНОМ. Лаборатория знаний, 2014. -216 с.
4. Архипов В.А., Ахмадеев И.Р., Бондарчук С.С., Ворожцов Б.И., Павленко А.А., Потапов М.Г. Модифицированный метод спектральной прозрачности измерения дисперсности аэрозолей // Оптика атмосферы и океана. - 2007. - Т. 20, №1. - С. 48-52.
5. Refractive index database: электронная база данных. URL: http://refractiveindex.info (дата обращения: 02.09.2016).
6. Кудряшова О.Б., Ворожцов Б.И., Антонникова А.А. Физико-математическая модель динамики функции распределения частиц по размерам с учетом процессов коагуляции, испарения и осаждения // Вестник Томского государственного университета. Математика и механика. -2012. - №1(17). - С. 81-90.
7. Архипов В.А, Усанина А.С. Движение аэрозольных частиц в потоке: Учеб. пособие. - Томск: Изд-во Том. ун-та, 2013. - 92 с.
Мецлер Эдуард Андреевич, младший научный сотрудник лаборатории физики преобразования энергии высокоэнергетических материалов Федерального государственного бюджетного учреждения науки Института проблем химико-энергетических технологий Сибирского отделения Российской академии наук (ИПХЭТ СО РАН), E-mail: mezlered@mail.ru Ул. Социалистическая, 1, г. Бийск, Россия. Тел. 8 (3854) 30-58-47.
Титов Сергей Сергеевич, кандидат технических наук, ученый секретарь Федерального государственного бюджетного учреждения науки Института проблем химико-энергетических технологий
Сибирского отделения Российской академии наук (ИПХЭТ СО РАН), E-mail: titov@ipcet.ru Ул. Социалистическая, 1, г. Бийск, Россия. Тел. 8 (3854) 30-17-25.
Муравлев Евгений Викторович,
старший научный сотрудник лаборатории физики преобразования энергии
высокоэнергетических материалов
Федерального государственного бюджетного учреждения науки Института проблем химико-энергетических технологий
Сибирского отделения Российской академии наук (ИПХЭТ СО РАН), E-mail: admin@ipcet.ru Ул. Социалистическая, 1, г. Бийск, Россия. Тел. 8 (3854) 30-58-47.
Павленко Анатолий Александрович, доктор физико-математических наук, заведующий лабораторией физики преобразования энергии высокоэнергетических материалов Федерального государственного бюджетного учреждения науки Института проблем химико-энергетических технологий Сибирского отделения Российской академии наук (ИПХЭТ СО РАН), E-mail: admin@ipcet.ru Ул. Социалистическая, 1, г. Бийск, Россия. Тел. 8 (3854) 30-14-43.
Архипов Владимир Афанасьевич, доктор физико-математических наук, главный научный сотрудник лаборатории физики преобразования энергии высокоэнергетических материалов Федерального государственного бюджетного учреждения науки Института проблем химико-энергетических технологий Сибирского отделения Российской академии наук (ИПХЭТ СО РАН), E-mail: leva@niipmm.tsu.ru Ул. Социалистическая, 1, г. Бийск, Россия.