CC BY
34
УДК 53 7
Т. В. Андрухова, Е.Ю. Потехин Динамика размеров и температуры аэрозольных частиц в мощных оптических полях при пониженном давлении
С развитием промышленного производства в крупных городах становится актуальной проблема загрязнения атмосферы промышленными аэрозолями. По данным экспериментов атмосферного контроля частицы образуют две группы - мелкие, менее 0,1 мкм и крупные, от 1 до 1000 мкм. Мелкие частицы образуются, главным образом, в результате химических процессов, а именно реакций горения. Их основными компонентами являются сернокислотные и сульфатные частицы, а также частицы, содержащие органические соединения. Частицы такого типа обладают ярко выраженными кислотными свойствами, являются источниками кислотных дождей. В отличие от них, крупные частицы образуются при механических разрушениях природных материалов.
Одним из источников антропогенного аэрозоля является сжигание топлива, которое влечет за собой выбросы в атмосферу частиц на основе углерода. Основными углеродосодержащими частицами в аэрозолях являются уголь и сажа, содержащие связанный углерод, золу (минеральный остаток), влагу, а также летучие вещества (водород, метан, смолы, которые выделяются при нагреве) . Одним из свойств атмосферных аэрозолей является их распространение на большие расстояния от источника образования. Аэрозольные дымки встречаются в настоящее время в любых областях земного шара, ранее считавшихся чистыми. Аэрозольные частицы имеют размеры, обеспечивающие наиболее эффективное рассеяние света. А это влияет на теплообмен между атмосферой и земной поверхностью, что может приводить как к охлаждению земной поверхности в результате экранирования солнечного света, так и нагреву верхних слоев атмосферы. В связи с вышеизложенным возникает потребность исследований процессов образования и разрушения аэрозолей, их поведения при воздействии мощного лазерного излучения (МЛИ) при различных давлениях, изучения возможности лазерного очищения атмосферы от аэрозолей.
Лазерные лучи можно использовать в различных слоях атмосферы, где давление меняется от одной атмосферы на уровне моря до вакуума на больших высотах. Интенсивное испарение твердой аэрозольной частицы под действием мощного лазерного излучения может сопровождаться образованием конденсационного аэрозоля. В допробойном интервале интенсивностей излучения явление переконденсации может оказаться решающим при анализе распространения МЛИ в аэрозольной среде.
Рассмотрим сферическую углеродную частицу радиуса а, помещенную в поле МЛИ и нагретую до температуры Т в условиях вакуума. Среда на большом расстоянии от частицы остается холодной, что приводит к большой неоднородности поля температур. Температура поверхности определяется плотностью потока энергии излучения, оптическими свойствами испаряющейся частицы и энергетическими потерями на испарение и переизлучение. Молекулы испаряющегося вещества могут находиться как в газообразной, так и в сконденсированной твердой или жидкой фазах. В работах [1, 2] рассмотрена газодинамическая модель испарения одиночной углеродной частицы в поле МЛИ в условиях вакуума с учетом переконденсации. Для решения газодинамических уравнений введена степень конденсации пара в виде с = гк / ( гк + г, ), где гк и гг - число молекул в единице объема, находящихся в конденсированной и газообразной фазах соответственно . Характерное время установления стационарных полей паров t = а / V , где V - скорость разлета пара. Скорость разлета испаренного углерода растет от значения скорости звука в газе до сверхзвуковых скоростей, так что при таких значениях мелкие частицы не успевают вырасти и далее разлетаются на бесконечность вместе с углеродным паром. При этом степень конденсации пара становится порядка 0,3, т. е. меньшая часть массы испаренного вещества превращается в частицы, а большая разлетается в виде газа. Температура двухфазной системы в окрестности частицы достаточно медленно убывает до расстояний г 20а, а затем резко стремится к температуре среды на бесконечности. Давление испаренного углерода у поверхности частицы вначале превышает атмосферное, но с увеличением расстояния от частицы быстро уменьшается и становится существенно ниже атмосферного. Степень конденсации испаренного вещества слабо зависит от температуры поверхности частицы или плотности потока энергии падающего излучения. Радиус вторичной частицы уменьшается со скоростью порядка СЦ /с^ = ш0^/гь, где аь - радиус вторич-
ной частицы, V - средняя тепловая скорость молекул пара, гь - плотность вещества вторичнои частицы, г - число атомов или молекул испаренного вещества в единице объема.
В области, где произошло образование зародышей, их величина а(г0) = а , т. е. приравнена к минимальному критическому размеру для того, чтобы он стал центром образования новой фазы. На всем интервале расстояний концентрация вторичных частиц убывает как 1 /г. При взаимодействии излучения с частицами единственным, как правило, известным параметром является плотность потока энергии. Плотность потока энергии от частицы складывается из энергопотерь за счет лучистого теплообмена и испарения с учетом переконденсации вещества. В [1, 2, 3] показано для сферически симметричного процесса испарения, что частица достигает максимальной температуры практически без изменения своего радиуса. Так при плотности потока энергии I = 109-1010Вт/м2 температура поверхности достигает значения 4000-5000 К практически независимо от размеров. В работах [1; 3-5] теоретически получены зависимости температуры поверхности и радиуса частицы от времени для различных размеров частиц. У крупных частиц с а^50 мкм температура поверхности после быстрого выхода на максимальное значение остается почти постоянной величиной до полного испарения частицы. Радиус частицы практически линейно убывает со временем. В случае испарения мелкой частицы с а0 10 мкм и менее наблюдается ее быстрое нагревание и быстрое охлаждение. Зависимости скорости испарения частицы от степени конденсации испарения вещества, от плотности потока энергии падающего излучения являются практически линейными функциями плотности потока энергии. Степень конденсации пара меняется в интервале исследуемых значений плотности потока энергии в 1,2 раза, а скорость испарения для I = (1-6)х109Вт/м2 - в 7,5 раза. В режиме интенсивного испарения скорость изменения радиуса Са/СЬ зависит от плотности потока энергии излучения в виде Са/СЬ = х1 , где х - размерный коэффициент, определяемый из уравнения теплового баланса.
В работе [6] экспериментально обнаружены высокодисперсные частицы в окрестности облучаемой лазером частицы Ш4С1, образованные в результате конденсации испаренного вещества. Для частиц хлорида аммония получено, что давление паров существенно больше давления насыщенного пара, поэтому уже вблизи частицы при испарении образуется существенное перенасыщение. Во всех опытах наблюдалась оболочка из частиц конденсационного типа, размер которой составлял (4-6) а. Центры конденсации могут представлять собой твердые частицы, которые всегда находятся в атмосфере, либо конденсационный аэрозоль (для вакуума), образующийся при значительном пересыщении пара. С целью наблюдения процесса испарения была сконструирована и
выполнена установка (рис. 1) :
1 Технологическая лазерная установка типа 7 "КИПР" с длиной волны, равной 10,6 мкм, полной
_____/ мощностью 900 Вт, работающая в непрерывном
режиме.
2 Вакуумная камера, имеющая окно из ЫаС1 -пластины, прозрачной для воздействующего мощного излучения, а также смотровые окна из стекла и отверстие для закрепления испаряемой частицы.
3. Поглотитель лазерного излучения, герметично закрепленный в нижней части вакуумной камеры.
4 Форвакуумный насос.
5. Вакуумметр термопарный типа ВТ-3 с подключенным к нему преобразователем ПМТ-4м.
4
5
2
3
Рис. 1. Принципиальная схема установки
-10-
мм рт. ст. с
Данная система позволяет измерять давление атмосферы в диапазоне (5^х10-1-погрешностью +20%.
6 Скоростная кинокамера для регистрации процесса испарения частицы.
7. Фокусирующая линза из ЫаС1 с Е = 20 см. Использование этой линзы позволило получить в области взаимодействия излучения и частицы диаметр пятна фокусировки 0,5 мм.
8. Электромагнитный затвор, время срабатывания которого составляет 25 мс.
Частица сажи или графита устанавливалась на кварцевую или нихромовую подложку. Влияние подложки на ход опыта было выявлено в предварительных экспериментах. Результаты опытов получены в виде кинограмм. Их обработка осуществлялась с помощью микрофотометра и микро-
скопа. Одними из объектов исследования являлись частицы сажи, графита, начальная форма которых близка к сферической. Использовалась сажа ПМ-100, в которой содержится: 99,1 - 99,3 % углерода, 0,3 - 0,4% водорода, 0,1% кислорода и химически чистый графит. Предварительно частицы отбирались и измерялись с помощью отсчетного микроскопа типа МПБ-2. Параметры некоторых опытов приведены в таблице.
Таблица.
№ Объект исследо- вания Размер, сі ,мкм Давление атмосферы, мм.рт.ст. Подложка Скорость, кадр/с Интенсивность, Вт/м2
1 сажа 425 3 х 10-1 нихром 1500 2,8 хг1 109
2 сажа 350 3 х 10-1 нихром 1500 3 х 109
3 сажа 300 750 нихром 1500 1,75 хГ109
4 сажа 350 2,51X10^ кварц 1200 3 х 108
5 сажа 1100 3 х10-1 нихром 750 2,5 Гх 109
6 сажа 1000 3 хґ 10-1 кварц 750 2,4 х 109
7 графит чистый 850 3 ґх 10-1 кварц 750 2,14 х 109
При воздействии мощного лазерного излучения с плотностью потока энергии излучения 109 Вт/м2 на частицу угля за время 1,3х10-3 с размер первичной частицы угля увеличивается в 1,4 раза, а при длительности воздействия 4, 0х10-2 с размер первичной частицы увеличивается 2-2, 5 раза. Процесс изменения размера частицы протекает неравномерно. Наблюдаются флуктуации до выхода размера первичной частицы на максимальное значение, соответствующее максимальной температуре 4000 - 5000 К. Это связано с нарушением строгой периодичности, присущей кристаллам. Частица претерпевает подобное изменение кристаллической структуры первоначально в поверхностном слое в момент начального воздействия на неё МЛИ. По мере прогревания частицы вглубь, сопровождающегося разлетом испаренного с поверхности вещества, происходит первоначально уменьшение размера частицы с дальнейшим его увеличением. Периодичностью данного процесса объясняются флуктуации размера частицы до ее полного прогревания. Тогда частицу можно рассматривать как переохлажденную жидкость с очень большим коэффициентом вязкости, у которой нарушения в кристаллической решетке произошли по всему её объёму. У крупных частиц температура поверхности после выхода на максимальное значение остается почти постоянной величиной до полного испарения частицы. Радиус частицы далее убывает практически линейно и за время 6, 67х10-2с уменьшается более чем в 1,3 раза. В дальнейшем происходит фрагментация угольной частицы. Каждый фрагмент является ядром образования конденсированного облака из вторичных частиц. Для сажистой частицы за время 2 10-2 с размер первичной частицы увеличивается в 2 и более раз. Так как сажа - это конгломират из мелких углеродных агрегатов, то
после нагрева частицы до температуры 4000 К происходит фрагментация сажистой частицы за счет кавитации внутри частицы и выхода испаренного вещества в результате разрушения полостей.
Динамика размеров испаряющейся частицы приведена на рисунке 2 (а, б), где сплошными линиями показаны результаты аппроксимации. Аппрокси-
Рис.2. Динамика размера сажистой частицы мация п°лученных ре-
а) а0 = 800 мкм, I = 2,8х108Вт/м2 зультатов описывается
б) а0 = 700 мкм, I = 2, 4х108Вт/м2 зависимостью вида
а('Ь) = A+Bt+Ct2+Dt3+Et4. В условиях I = 2,8 108Вт/м2, 2,4 ^108Вт/м2:
А = 800 мкм, 700 мкм;
В = -9 10-3м/с, 0,054м/с;
С = 0,74 м/с2, - 1,15 м/с2;
D = -9 м/с3, 10 м/с3;
Е = -28,8 м/с4, -30,7 м/с4, для рисунков 3 и 4 соответственно. Физический смысл параметров А,
В, С - начальный размер, скорость и ускорение испарения частипы соответственно. Измерения
температуры проводи-
лись с использованием пветового пирометра. К нему присоединялся запоминающий оспил-лограф в ждущем режиме с разверткой сигнала 2 с/см, на экране которого появлялась оспиллограмма пропес-са. При помощи градуировочных кривых, которые строились на установке "Проминь", снимались температурные характеристики пропесса. Результаты измерения температуры испаряющейся час-
Т,К
1500-
—I—
500
1000
2000
ІМС
а б
Рис. 3. Динамика температуры испаряющейся сажистой частиц
а) а0 = 800 мкм, I = 2,8х108 Вт/м2
б) а0 = 1000 мкм, I = 2,4х108Вт/м2
тицы приведены на рисунке 3 (а, б) , они являются типичными графиками, где сплошными линиями показаны результаты аппроксимации экспериментальных данных динамики температур, описываемых следующим выражением:
Т(У = А*+В^+С*^+Б*^+Е*^, где
А*= 1500 К, 1598 К;
В*= 2870 К/с, 4190 К/с;
С*= -544 К/с2, -11000 К/с2;
Б*= -1210 К/с3, 12 104 К/с3;
Е*= -386 К/с4, -4200 К/с4, для рисунков 5 и 6 соответственно. Физический смысл А*, В*, С* для уравнения Т('Ь) - начальная температура, скорость и ускорение ее изменения соответственно. Полученные экспериментально функции дают приблизительное решение уравнений [1, 2] :
ёа / Л = 10 /а0 {(го &0,Т8)у(а0,Т8)т0) /гь};
ЛТ / Л = І0 /(^Т ){3 1, I) / (4раС Гь)},
где а0 - начальный размер частицы, I - интенсивность излучения, п - первоначальная масса частицы Тз = Т/Тз, t = t / t0, а = а / а0 со следующими краевыми условиями: t = 1, а = 1, с точностью 3%.
Авторы выражают благодарность В.И. Букатому и Ю.Т. Лялину.
Литература
1. Лямкина Г.В., Шайдук А.М. Испарение и переконденсация тугоплавких частиц в мощном оптическом поле. Барнаул, 1990. Деп.в ВИНИТИ 11.11.90г. № 5653-ВУО.
2. Шайдук А.М. Моделирование процессов взаимодействия мощного лазерного излучения с дисперсными системами: Дисс. докт. физ.-мат.наук. Барнаул, 1998.
3. Лямкина Г.В., Шайдук А.М. Размеры и концентрация мелкодисперсной фракции аэрозоля при испарении тугоплавкой частицы в вакуум // Оптика атмосферы. 1991. Т 4. № 11.
4. Лямкина Г.В. Взаимодействие мощного лазерного излучения с тугоплавким аэрозолем в условиях вакуума с учетом переконденсации: Дисс. докт. физ.-мат.наук. Барнаул, 1992.
5. Букатый В.И., Суторихин И.А. Экспериментальные исследования воздействия излучения СО2-лазера на углеродные частицы // Физика горения и взрыва. 1982. № 2.
6. Свердлов М.Ю. Исследование переконденсации при воздействии мощного лазерного излучения на термически неустойчивые частицы / / Труды XI симпозиума по распространению лазерного излучения в атмосфере. Томск, 1991.
