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

Размещено на http://www.allbest.ru/

УДК 517.968

Томский политехнический университет

О ПОВЫШЕНИИ ДОСТОВЕРНОСТИ РЕШЕНИЯ ИНТЕГРАЛЬНОГО УРАВНЕНИЯ ПРИ БЕЗОТБОРНОЙ ДИАГНОСТИКЕ НАНОЧАСТИЦ

Мышкин В.Ф., д.ф.-м.н., профессор

Власов В.А., - д.ф.-м.н., профессор

Хан В. А., - д.т.н.

Бурдовицын А. Н., - аспирант

Аннотация

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

Ключевые слова: ПОВЫШЕНИЕ ДОСТОВЕРНОСТЬ РЕШЕНИЕ ИНТЕГРАЛЬНОЕ УРАВНЕНИЕ ДИАГНОСТИКА НАНОЧАСТИЦА

Гранулометрический состав аэрозолей оказывает существенное влияние на перенос оптического излучения в атмосфере, определяет физико-химические процессы в гетерогенных (газодисперсных) системах.

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

Основные вопросы взаимодействия оптического излучения с одиночной дисперсной частицей и соотношения для расчета интенсивности рассеянного в разных направлениях света на малых частицах впервые рассмотрены в работе [1]. Теоретической основой интегральных оптических методов диагностики дисперсности аэрозолей является решение интегрального уравнения Фредгольма первого рода. Задача сводится к численному решению с помощью вычислительных машин обратной задачи рассеяния зондирующего оптического излучения.

К обратным методам оптики аэрозолей относятся: метод спектральной прозрачности (МСП), метод малых углов (ММУ), метод обращения полной индикатрисы рассеяния (МПИ) и лидарные дистанционные методы диагностики гранулометрического состава аэрозолей.

Например, задача обработки экспериментальных индикатрис рассеяния сводится к решению интегрального уравнения Фредгольма первого рода:

наночастица лазерный спектр дисперсный

,(1)

где (,) - ядро интегрального уравнения - эффективное сечение рассеяния дисперсных частиц; f () - функция распределения дисперсных частиц по размерам; I () - индикатриса рассеяния. При неточно известных значениях величин ядра интегрального уравнения и его правой части интегральные уравнения могут иметь множество решений, т.е. задача некорректна [2]. В операторном виде уравнение (1) имеет вид:

Ay = fy Y, f F,(2)

где Y, F - некоторые метрические пространства, А - непрерывный оператор, отражающий метрические пространства Y на F.

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

,(3)

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

Вектор регуляризованного решения S, минимизирующий функционал Тихонова, удовлетворяет системе линейных алгебраических уравнений:

.(4)

Здесь А^* - матрица, транспонированная к матрице А.

Сглаживающая матрица Н имеет вид:

,(5)

где р_о, р₁ - масштабные множители, соответственно равные (r₂ - r₁) и (r₂ - r₁)³, r₁, r₂ - границы интервала радиусов дисперсных частиц.

Важным и достаточно сложным вопросом при построении регуляризованного решения является выбор величины параметра регуляризации . При малых значениях решаемое уравнение близко к исходному. С другой стороны, при очень малом , в силу малой информативности оптических измерений и наличия ошибок, полученное решение будет сглаженным (или вовсе не иметь физического смысла). Параметр регуляризации в уравнении (3) зависит от уровня ошибок экспериментальных данных. Наиболее простыми являются: выбор параметра регуляризации по невязке, поиск квазиоптимального значения .

На основе известных алгоритмов решения обратной задачи нами был сформирован программный комплекс для обработки экспериментальных спектров коэффициентов пропускания различных аэрозолей [3].

Тестирование программы обработки экспериментальных данных проводили в два этапа. На первом этапе выполнялись расчеты по схеме: дисперсность спектр пропускания дисперсность при различных уровнях вносимых в спектр пропускания случайных и систематических ошибок. Установлено, что при суммарной экспериментальной ошибке, не превышающей 10 %, расчетная погрешность, вносимая программой обработки спектров, не превышает 25 %. На втором этапе сопоставлялись данные о дисперсности в счетном объеме, получаемые из анализа процессов седиментации различных фракций ультрадисперсного порошка (УДП) меди в циклогексане и обращением интегрального уравнения с регистрируемыми в те же моменты времени коэффициентами пропускания этой суспензии.

Для тестирования программы обработки спектров был разработан экспериментальный стенд (рисунок 1), состоящий из источника широкого спектра (лампа накаливания), объектива "Гелиос-44", прямоугольной кварцевой кюветы с анализируемой суспензией, фотоприемника (ФЭУ-28), установленных на оптической оси монохроматора. Рабочая и контрольная кюветы размерами 1501010 мм располагались на каретке, имеющей возможность перемещения перпендикулярно оптической оси системы. Зондирование осуществлялось лучом с диаметром 5 мм на отметке 100 мм ниже верхней свободной поверхности суспензии. Питание лампы накаливания, используемой в качестве источника света, осуществлялось стабилизированным напряжением.

Размещено на http://www.allbest.ru/

Измерение уровня напряжения ФЭУ проводилось с помощью цифрового вольтметра и осциллографа. При подготовке суспензии ультрадисперсный порошок меди насыпался в кювету с жидкостью. Суспензия размешивалась и разбавлялась так, чтобы оптическая плотность в кювете не превышала значения 2,0 во всем спектральном интервале измерения. Далее кювета с суспензией облучалась ультразвуком для разбивания возможных конгломератов в течение 30 мин. Кювета помещалась на стенд, предназначенный для измерения коэффициентов пропускания. Определялись значения интенсивности прошедшего суспензию излучения на дискретных длинах волн в видимой области, по которым строился нулевой спектр коэффициентов пропускания.

После обработки суспензии в центрифуге в течение 15 с кювета устанавливалась в колориметрическую установку. Измерялся первый спектр. Далее кювета с суспензией помещалась в центрифугу и процедура повторялась. Для каждого отсчета измерялись уровни напряжений ФЭУ на выбранных длинах волн с пустой кюветой и с кюветой с суспензией, по которым определялся коэффициент пропускания суспензии во всем рабочем диапазоне спектра. Весь тестовый эксперимент с суспензией продолжался непрерывно в течение 5 часов. Экспериментальные данные, снятые в выборочные моменты времени, приведены ниже. Рассчитанная по нулевому спектру поглощения (методом МСП) дисперсность порошка меди принималась за исходный гранулометрический состав при расчете по седиментационной формуле.

Влияние коагуляции частиц с размерами 0,01-0,5 мкм в циклогексане оценивалось по отношению интенсивностей рассеяния в направлении углов 3 и 175° на длине волны 0,63 мкм после размешивания ультразвуком и по истечении времени не менее времени регистрации коэффициентов пропускания. Для этого суспензия наливалась в стеклянный цилиндр высотой 140 мм, а зондирование осуществлялось на расстоянии 120 мм от свободной поверхности. Установлено, что при первом измерении отношение интенсивностей рассеянного вперед и назад излучения составило 1,60, а при втором - 1,62. Поэтому влиянием коагуляции за время измерения можно было пренебречь.

Экспериментально определены коэффициенты поглощения суспензией излучения в видимой области спектра (рисунок 2). Номера гистограмм соответствуют порядковому номеру измерения от начала эксперимента. За время измерения наблюдается изменение спектра поглощения. Оптическая плотность в измеряемом спектральном интервале находится в диапазоне 0,1-1,9 (коэффициент пропускания 0,9-0,16). Уменьшение оптической плотности на каждой длине волны происходит не монотонно, а имеет участки возрастания и убывания. Значения коэффициентов пропускания, определяемых с погрешностью менее 2,5 %, увеличиваются с увеличением длины волны.

Размещено на http://www.allbest.ru/

На рисунке 3 приведены функции распределения дисперсных частиц по размерам в различные моменты времени, полученные обработкой на ЭВМ экспериментальных спектров поглощения суспензии. Номера гистограмм приведены в соответствии с рисунком 2.

В начальный момент времени функция распределения дисперсных частиц по размерам имеет распределение с модовым диаметром 0,08 мкм. Скорость седиментации частиц разных размеров оценивалась соотношением [4]:

,(6)

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

Размещено на http://www.allbest.ru/

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

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

Размещено на http://www.allbest.ru/

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

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

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

Размещено на http://www.allbest.ru/

С целью поиска критерия для составления наиболее эффективной выборки были проведены измерения коэффициентов пропускания модельного аэрозоля. На рисунке 5 приведен спектр пропускания взвеси ультрадисперсного порошка оксида алюминия в воде. Коэффициент пропускания растет с длиной волны. При этом спектр имеет участки с локальными экстремумами.

Для определения дисперсности УДП из всего набора (рисунок 5) было составлено 15 выборок по коэффициентам пропускания на десяти длинах волн. Эти длины волн выбирались в областях некоторых локальных экстремумов графика спектра пропускания. При выполнении расчетов для каждого значения проводилось несколько итераций (iter). Расчеты останавливались, когда изменение невязки между расчетным и экспериментальным коэффициентами пропускания проходило через экстремум.

При проведении расчетов проводили оптимизацию также по значениям комплексного показателя преломления и диапазона размеров дисперсных частиц. Три результата обработки (из множества случайных выборок) по десяти значениям длин волн из спектра пропускания приведены в таблице 1. Использованы следующие обозначения: R - радиус дисперсной частицы, m - комплексный показатель преломления материала частицы, - начальное значение параметра регуляризации, iter - число итераций решения интегрального уравнения с заданным значением , U - невязка между расчетными и экспериментальными (k_пр) коэффициентами пропускания на выбранных длинах волн, Z_rach - расчетная гистограмма распределения дисперсных частиц по размерам.

Таблица 1 - Данные, полученные обращением интегрального уравнения

При обработке различных выборок из экспериментальных значений коэффициентов пропускания появляются отрицательные значения числа частиц (таблица 1). Для сопоставительного анализа гистограмм составлялись выборки, при обработке которых отрицательное значение не превышало 3 % от максимальной величины. В таблице 1 жирным шрифтом выделены параметры гистограммы в максимуме распределения по размерам.

Модальные радиусы получаемых распределений по размерам имели следующие значения: 0,200; 0,241; 0,255; 0,274; 0,283; 0,320; 0,324; 0,327; 0,330; 0,335; 0,341; 0,366; 0,444; 0,491; 0,514. Эти данные для наглядного представления приведены на числовой оси (рисунок 6).

Из рисунка 6 видно, что модальные радиусы шести результатов расчетов объединяются в компактно расположенную группу. Это позволяет утверждать (с высокой достоверностью), что максимум кривой распределения по размерам находится в области 0,32-0,34 мкм. Расчеты по остальным выборкам были остановлены в связи с попаданием результатов расчета в один из многочисленных локальных минимумов функционала Тихонова.

Авторы выражают искреннюю признательность д.ф.-м.н., профессору ТПУ Корепанову В. И. за оказанную помощь при проведении экспериментальных исследований.

Выводы

Результаты показывают, что путем обращения интегрального уравнения с использованием экспериментального спектра пропускания аэрозолей, регистрируемого на большом числе длин волн зондирования, можно получать достоверные данные о гранулометрическом составе дисперсной фазы УДП, даже при достаточно высокой экспериментальной погрешности измерения интенсивностей зондирующих потоков.

Литература

1. Ван де Хюлст Г. Рассеяние света малыми частицами.- М.: Изд. Иностр. лит., 1961. - 320с.

2. Верлань А.Ф., Сизиков В.С. Интегральные уравнения: методы, алгоритмы, программы. Спр. пос. Киев. “Наукова Думка”. 1986.- 544с.

3. Мышкин В.Ф., Тихомиров И.А., Цимбал В.Н. и др. Лазерная диагностика гранулометрического состава дисперсной фазы плазмы горения пиротехнических составов // Заводская лаборатория. 1999.- №3.- Т.65.- С.24-27.

4. Диагностика металлических порошков/ В. Я. Буланов, Л. И. Квартер, Т. В. Долгаль, Т. А. Угольникова. В. Б. Акименко.- М.: Наука, 1983.- 278 с.

Размещено на Allbest.ru