Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,074

МОДЕЛИРОВАНИЕ ОСЛАБЛЕНИЯ ОПТИЧЕСКОГО ИЗЛУЧЕНИЯ В ГАЗОДИСПЕРСНОЙ СИСТЕМЕ «УГОЛЬНАЯ ПЫЛЬ ‒ ВОЗДУХ»

Лисаков С.А. 1 Павлов А.Н. 1 Сыпин Е.В. 1 Леонов Г.В. 1
1 Бийский технологический институт (филиал) Алтайский государственный технический университет имени И.И. Ползунова
Проведено моделирование ослабления оптического излучения в газодисперсной системе «угольная пыль ‒ воздух». Определены основные аналитические зависимости, описывающие ослабление излучения в газодисперсной системе. Выполнена проверка адекватности компьютерной модели. Проведено исследование ослабления оптического излучения на базе компьютерного моделирования. Моделирование выполнялось на основе теории Ми и алгоритма BHMIE. При проверке адекватности модели установлено, что результаты расчетов на базе теории Ми и экспериментальные данные, описывающие поглощательную способность, хорошо согласуются друг с другом. Врезультате исследования ослабления оптического излучения в газодисперсной системе «угольная пыль ‒ воздух» установлено, что ослабление оптического излучения может быть охарактеризовано через поглощающую способность среды. Промежуточная среда в виде газодисперсной системы «угольная пыль ‒ воздух» существенно ограничивает применение оптико-электронных датчиков, что необходимо учитывать при их использовании для контроля возникновения очагов возгорания при различных технологических процессах. Результаты моделирования будут использоваться для повышения помехоустойчивости оптико-электронных датчиков.
оптическое излучение
газодисперсная система
ослабление
1.Блох А.Г. Теплообмен в топках паровых котлов. – Л.: Энергоатом-издат. Ленингр. отд-ние, 1984. – 240 с., ил.
2.Блох А.Г. Тепловое излучение в котельных установках. – Л.: Энергия, 1967. – 328 с.
3.Блох А.Г. Основы теплообмена излучением. – Л. Госэнергоиздат, 1962. – 332 с.
4.Борен Г.Ф. Поглощение и рассеяние света малыми частиц-ми: пер. с англ./ Г.Ф. Борен, Д.Р. Хафмен. – М.: Мир, 1986. – 664 с.: ил
5.Дейрменджан Д. Рассеяние электромагнитного излучения сферическими полидисперсными частицами: пер. с англ. – М.: Мир, 1971. – 164 с.
6.Патент РФ №2459269: МПК7 G 08 B 17/12. Пирометрический датчик координат очага возгорания с цилиндрическими линзами/ Е.В. Сыпин и др.; заявитель и патентообладатель ГОУВПО АлтГТУ им. И.И. Ползунова// Бюл. – 2012.
7.Патент РФ №2536355: МПК: G08B17/00. Пирометрический датчик координат очага возгорания с полевой диафрагмой/ Е.В. Сыпин и др.; заявитель и патентообладатель ГОУВПО АлтГТУ им. И.И. Ползунова// Бюл. – 2014.
8.Повернов Е.С. Оптико-электронный датчик координат очага возгорания на базе одноэлементного позиционно-чувствительного детектора/ Е.С. Повернов, А.Н. Павлов, Е.В. Сыпин// Южно-Сибирский научный вестник. – 2012. – №2. – С. 157–160. – Режим доступа: http://s-sibsb.ru/images/articles/2012/2/39_157-160.pdf.
9.Повернов Е.С. Лабораторный образец быстродействующего многопорогового прибора контроля аварийных и предаварийных ситуаций во взрывоопасной атмосфере// Ползуновский вестник. – 2013. – №2. – С. 172–178.
10.Тупикина Н.Ю. Оптическая система трёхканального оптико-электронного прибора двух спектральных отношений с раздельными объективами каналов/ Н.Ю. Тупикина, Е.В Сыпин// Южно-Сибирский научный вестник. – 2013. – №1 (3). – С. 25–30. – Режим доступа: http://s-sibsb.ru/images/articles/2013/1/6_25-30.pdf.
11.Laven, P. MiePlot. A computer program for scattering of light from a sphere using Mie theory & the Debye series [электронный ресурс]. – 2013. – Режим доступа: http://www.philiplaven.com.
12.Lisakov S.A. Determination of Optical Radiation Attenuation in Dispersed System «Coal Dust-Air»// 16th International Conference of Young Specialists on Micro/Nanotechnologies and Electron Devices EDM 2015: Conference proceedings. – Novosibirsk: NSTU, 2015. – P. 353–358.
13.Mikhanoshina J.L. The study of influence of gas-dispersed medium like coal dust-air on the optical radiation attenuation/ Mikhanoshina, A.N. Pavlov, E.V. Sypin// International conference and seminar on micro/nanotechnologies and electron devices EDM 2012: Conference proceedings. – Novosibirsk: NSTU, 2012. – P. 226–229.
14.Pavlov A.N. Optoelectronic system for determination of ignition center three-dimensional coordinates at initial stage// International conference and seminar on micro/nanotechnologies and electron devices EDM’2010: Conference proceedings. – Novosibirsk: NSTU, 2010. – P. 417–419.
15.Yakusheva O.Y. Experimental research backscattering in the disperse system/ O.Y. Yakusheva, A.N. Pavlov, E.V. Sypin// 14th International Conference of Young Specialists on Micro/Nanotechnologies and Electron Devices EDM 2013: Conference proceedings. – Novosibirsk: NSTU, 2013. – P. 246–248.

Для обеспечения пожарной безопасности взрывоопасных технологических процессов, связанных с образованием газодисперсных систем (в частности, при добыче угля в шахтах), используются автоматические системы противопожарной защиты и взрывоподавления. Такие системы, как правило, строятся на основе быстродействующих оптико-электронных датчиков (ОЭД) [8–10, 14].

Излучение от очага возгорания приходит на ОЭД ослабленным за счет влияния промежуточной среды, в которой оно распространяется. Промежуточная среда оказывает существенное влияние на работу ОЭД, приводя к значительным погрешностям или делая их работу невозможной. Поэтому учет влияния запыленной среды на работу оптико-электронных датчиков является актуальной задачей [13, 15] и одной из основных при выполнении проекта «Разработка научных основ построения системы предотвращения и локализации взрывов на потенциально опасных промышленных объектах» (грант РФФИ №15-08-06719 а). Врамках указанного проекта в работе [12] выполнена постановка задачи моделирования ослабления оптического излучения в газодисперсной системе «угольная пыль ? воздух» и сформулированы исходные данные для проведения моделирования, приведенные в табл.1.

В результате моделирования ослабления оптического излучения в газодисперсной системе «угольная пыль ? воздух» необходимо:

–вычислить коэффициенты рассеяния и поглощения для полидисперсной системы в зависимости от длины волны падающего излучения;

Таблица 1

Исходные данные для проведения моделирования

Параметр

Значение

Диапазон длин волн излучения l, мкм

0,7–5

Диаметр частиц пыли x, мкм

1–100

Плотность распределения fm(x) массы пыли по диаметрам частиц x

Распределение Розина ? Раммлера

fm(x)=baxa–1exp(–bxa)

с коэффициентами a=2 и b=0,005

Комплексный показатель преломления m=n –iχ частиц пыли каменного угля марки Г

n(λ)=1,6+0,08λ0,5

χ(λ)=0,33 – 0,025λ0,5 – 0,003(λ – 1)2

Концентрация частиц угольной пыли, г/м3

0,06–30

–получить индикатрисы рассеяния оптического излучения для полидисперсной системы частиц;

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

В связи с вышеизложенным была сформулирована цель: выполнить моделирование ослабления оптического излучения в газодисперсной системе «угольная пыль ? воздух».

Для достижения поставленной цели необходимо решить следующие задачи:

–определить основные аналитические зависимости, описывающие ослабление излучения в газодисперсной системе «угольная пыль ? воздух», и на их основе построить компьютерную модель;

–выполнить проверку адекватности компьютерной модели;

–провести исследование ослабления оптического излучения в газодисперсной системе «угольная пыль – воздух» на базе компьютерного моделирования.

Математическая модель ослабления оптического излучения в газодисперсной системе «угольная пыль ? воздух»

В работе [12] сделан вывод о том, что для расчета коэффициента ослабления в газодисперсной системе «угольная пыль ? воздух» может использоваться теория Ми.

К основным параметрам теории рассеяния Ми, характеризующим перенос оптического излучения в промежуточной среде, относятся: спектральные факторы эффективности: поглощения – Ka(λ), рассеяния – Ks(λ) и ослабления K(λ)=Ka(λ)+Ks(λ); а также спектральные коэффициенты: поглощения α(λ), рассеяния β(λ) и ослабления k(λ)=α(λ)+β(λ).

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

Факторы эффективности Ka(λ), Ks(λ) и K(λ) являются безразмерными величинами, а коэффициенты α(λ), β(λ) и k(λ) имеют размерность м–1.

На основании электромагнитной теории Г.Максвелла Ми были теоретически установлены формулы, связывающие факторы эффективности рассеяния и ослабления с диаметром частицы х, длиной волны падающего излучения λ и комплексным показателем преломления частицы m(λ). Дифракционные свойства частиц учитываются путем использования в расчетных формулах безразмерного параметра (дифракции) ρ=πx/λ.

Факторы эффективности рассеяния и ослабления излучения определяются в соответствии с формулами [5]:

lisakov01.wmf

lisakov02.wmf

где n – положительные целые числа.

Коэффициенты Ми an, bn рассчитываются в соответствии с рекуррентными выражениями, представленными в форме, удобной для вычисления на ЭВМ (алгоритм BHMIE [4]):

lisakov03.wmf

lisakov04.wmf

где An(ρ·m(λ)) и ξn(ρ) являются комплексными рекуррентными соотношениями и рассчитываются по выражениям

lisakov05.wmflisakov06.wmflisakov07.wmf

lisakov08.wmf lisakov09.wmf.

Содержание в дисперсной системе частиц различных размеров характеризуется числовой концентрацией частиц N0 и функцией их распределения по размерам f(x).

Коэффициенты поглощения α(λ), рассеяния β(λ) и ослабления k(λ) для полидисперсной системы частиц определяются зависимостями вида

lisakov10.wmf

lisakov11.wmf

lisakov12.wmf

где γ – плотность угольной пыли (1,8?106г/м3); μ – массовая концентрация частиц, г/м3.

Массовая концентрация связана с числовой концентрацией следующей зависимостью:

lisakov13.wmf

Из приведенных в табл.1 исходных данных видно, что параметр дифракции ρ>1, поэтому рассеяние в системе «угольная пыль – воздух» происходит на средних и больших частицах. Рассеянная такими частицами энергия концентрируется в узком направленном вперед пучке и практически не вычитается из основного потока [1, 2, 3]. При этом ослабление при достаточно больших концентрациях или толщинах слоя происходит главным образом за счет многократного поглощения на частицах, расположенных на пути пучка.

Полная поглощательная способность среды характеризуется коэффициентом ослабления лучей в поглощающей среде. Таким образом, поглощательная способность a(λ) зависит как от коэффициента поглощения частиц α(λ), так и от их рассеивающей способности, которая характеризуется коэффициентом рассеяния β(λ) и индикатрисой рассеяния γ(ϑ) [1].

Рассеянное частицей излучение, как правило, распределяется по различным направлениям ϑ неравномерно, т.е. является анизотропным. Анизотропию рассеяния характеризуют индикатрисой рассеяния γ(ϑ, λ):

lisakov14.wmf

в которой величина i1 характеризует интенсивность рассеяния излучения, поляризованного в направлении, перпендикулярном плоскости рассеяния, а величина i2 характеризует интенсивность рассеяния излучения, поляризованного в плоскости рассеяния [1].

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

Поглощательная способность слоя рассчитывается по формуле [1]:

lisakov15.wmf (1)

lisakov16.wmf

lisakov17.wmf,

где Sc(λ)=Ks(λ)/(Ka(λ)+Ks(λ)) – критерий Шустера; l – толщина поглощающего слоя, м; ηнз – параметр, характеризующий долю энергии, рассеянной частицей назад в полусферическом телесном угле [2].

Параметр ηнз рассчитывается согласно формуле

lisakov18.wmf

Направление ϑ=180° соответствует направлению падающего излучения.

Зависимость (1) для поглощательной способности слоя учитывает, как соотносятся фактор эффективности рассеяния с фактором эффективности поглощения через критерий Шустера Sc (для больших частиц Sc≈0,5)[1].

Для расчета ослабления оптического излучения в газодисперсной системе «угольная пыль ? воздух» по приведенным аналитическим зависимостям использовался программный пакет MiePlot v4.4, поскольку он позволяет рассчитывать коэффициенты поглощения и рассеяния для полидисперсных систем частиц [11]. ВMiePlot реализован классический алгоритм BHMIE, используемый для расчета характеристик рассеяния на однородном шаре, опубликованный в [4].

Проверка адекватности компьютерной модели

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

В источниках [1, 2, 3] приведены экспериментальные данные о поглощающей способности газовой среды, содержащей зольную пыль с концентрацией от 10 до 220г/м3 и средней удельной поверхностью пыли от 0,1 до 0,4м2/г, в диапазоне длин волн от 1 до 5мкм. На основании этих данных приведена зависимость с использованием закона Бугера для расчета поглощающей способности частиц золы, содержащейся в газовом потоке:

lisakov19.wmf

где k(λ) – коэффициент ослабления излучения; F – средняя удельная поверхность пыли, м2/г; μ – концентрация пыли, г/м3; l – толщина поглощающего слоя, м.

Для расчета величины k(λ) дается следующая эмпирическая формула [1]:

lisakov20.wmf (2)

где ρ – параметр дифракции; A′ – коэффициент, характеризующий влияние на коэффициент k(λ) электрооптических свойств вещества и формы пылевых частиц.

Эмпирическая зависимость (2) характерна для случая, когда ослабление происходит главным образом за счет поглощения K(λ)≈Ka(λ) [1].

В табл.2 приведены значения коэффициента A′ в зависимости от рода пыли [1, 3].

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

Данные для проверки адекватности соответствуют параметрам дисперсной системы, для которой была получена зависимость (2) (табл.3).

Таблица 2

Значения коэффициента A′

Род пыли

A′

Зола печорского угля

0,25–0,20

Зола воркутинского и кузнецкого угля

0,23

Зола подмосковного бурого и донецкого газового углей; зола эстонского сланца, зола ирша-бородинского бурого угля

0,15–0,2

Зола экибастузского угля

0,12

Зола антрацитового штыба

0,08

Угольная пыль различных видов топлива

0,1–0,15

Таблица 3

Исходные данные для проверки адекватности

Параметр

Значение

Диапазон длин волн излучения l, мкм

1–5

Диаметр частиц пыли x, мкм

1–50

Плотность распределения f(x) частиц золы ирша-бородинского бурого угля по диаметрам частиц x

lisakov21.wmf с модальным диаметром частиц xm=6,25мкм

Комплексный показатель преломления m=n – iχ частиц пыли золы ирша-бородинского бурого угля

l, мкм

n(λ)

χ(λ)

0,7

1,674

0,038

1,2

1,691

0,1

1,7

1,65

0,194

2,2

1,891

0,302

2,7

1,993

0,27

3,2

1,881

0,119

3,7

1,889

0,169

4,2

1,846

0,132

4,7

1,839

0,215

5

1,913

0,225

Концентрация частиц угольной пыли, г/м3

10

Данные, приведенные в табл.3, близки к данным, которые будут использованы при моделировании (табл.1), и могут быть использованы при проверке адекватности компьютерной модели.

На рис.1 показаны зависимости поглощательных способностей от длины волны излучения, полученные на основе теоретических и эмпирических зависимостей для концентрации 10г/м3, средней удельной поверхности пыли 0,13м2/г и толщины поглощающего слоя 1м.

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

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

Можно сделать вывод о том, что разработанная компьютерная модель на базе теории Ми применима для расчета ослабления оптического излучения в газодисперсной системе «угольная пыль –воздух».

pic_36.tif

Рис.1. Зависимости поглощательных способностей от длины волны излучения, полученные на основе теоретических и эмпирических зависимостей

Исследование ослабления оптического излучения в газодисперсной системе «угольная пыль –воздух» на базе компьютерного моделирования

Исходные данные, использованные для проведения исследования ослабления оптического излучения в газодисперсной системе «угольная пыль – воздух», приведены в табл.1. Поскольку дисперсный состав пыли представлен функцией плотности распределения fm(x) массы пыли по диаметрам частиц x, а в формулах для расчета коэффициентов поглощения α(λ), рассеяния β(λ) и ослабления k(λ) используется функция плотности распределения частиц по размерам f(x), необходимо выполнить пересчет с помощью выражения [2]:

lisakov22.wmf

Коэффициент поглощения α(λ) рассчитывается согласно формуле [2]:

lisakov23.wmf

Коэффициенты рассеяния β(λ) и ослабления k(λ) рассчитываются аналогично.

С учетом приведенных зависимостей проведено компьютерное моделирование ослабления оптического излучения в газодисперсной системе «угольная пыль – воздух».

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

pic_37.tif

Рис.2. Коэффициенты рассеяния и поглощения для полидисперсной системы в зависимости от длины волны падающего излучения при различных концентрациях угольной пыли

Из данных, представленных на рис.2, видно, что коэффициент рассеяния β(λ) излучения газодисперсной системой по порядку величины соизмерим с коэффициентом поглощения α(λ) и немного превышает α(λ) во всем рассматриваемом диапазоне длин волн.

Индикатрисы рассеяния оптического излучения полидисперсной системой частиц для длин волн 1, 2, 3, 4 и 5мкм представлены на рис.3. Функция γ(ϑ, λ) является безразмерной (интенсивности i1 и i2 в относительных единицах), а ее значения определяются в зависимости от фактора эффективности ослабления Ks(m(λ), ρ) и параметра дифракции ρ.

С уменьшением длины волны падающего излучения (увеличивается ρ) происходит существенная деформация индикатрис рассеяния. Рассеяние вперед значительно возрастает по ходу распространения падающего излучения.

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

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

Зависимость поглощательной способности газодисперсной системы «угольная пыль – воздух» от длины волны при различных концентрациях угольной пыли (толщина поглощающего слоя 1м) представлена на рис.5,а.

pic_38.tif

Рис.3. Индикатрисы рассеяния оптического излучения полидисперсной системой частиц для длин волн 1, 2, 3, 4 и 5 мкм

pic_39.tif

Рис.4. Зависимость параметра ηнз от длины волны падающего излучения

pic_40.tif pic_41.tif

а б

Рис. 5. Зависимость поглощательной способности газодисперсной системы«угольная пыль –воздух»:а – от длины волны при различных концентрациях угольной пыли; б – от произведения толщины поглощающего слоя и концентрации угольной пыли

Из данных представленных на рис.5,а, видно, что поглощательная способность газодисперсной системы a(λ) слабо зависит от длины волны падающего излучения. Например, при концентрации пыли 5г/м3 в диапазоне длин волн от 0,7 до 4мкм происходит увеличение поглощательной способности с ростом длины волны – a(0,7)=0,612 и a(4)=0,646. При этом в диапазоне от 4 до 5мкм поглощательная способность уменьшается с ростом длины волны – a(4)=0,64. Для различных концентраций характер зависимости поглощательной способности от длины волны сохраняется.

Поглощательная способность рассматриваемой газодисперсной системы в значительной мере зависит от концентрации угольной пыли и толщины поглощающего слоя (рис.5,б). Например, при толщине поглощающего слоя 10м и концентрации 1,4г/м3 поглощательная способность превышает 0,9.

Можно сделать вывод о том, что промежуточная среда в виде газодисперсной системы «угольная пыль – воздух» существенно ограничивает применение оптико-электронных датчиков, что необходимо учитывать при их использовании для контроля возникновения очагов возгорания.

Выводы и заключение

В результате выполнения работ проведено моделирование ослабления оптического излучения в газодисперсной системе «угольная пыль – воздух».

Решены следующие задачи:

–определены основные аналитические зависимости, описывающие ослабление излучения в газодисперсной системе «угольная пыль – воздух»;

–выполнена проверка адекватности компьютерной модели;

–проведено исследование ослабления оптического излучения в газодисперсной системе «угольная пыль – воздух» на базе компьютерного моделирования.

Моделирование ослабления оптического излучения в газодисперсной системе «угольная пыль – воздух» выполнялось на основе теории Ми и алгоритма BHMIE.

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

В результате исследования ослабления оптического излучения в газодисперсной системе «угольная пыль – воздух» установлено, что ослабление оптического излучения может быть охарактеризовано через поглощающую способность среды.

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

Дальнейшая работа направлена на проведение экспериментального исследования ослабления оптического излучения в газодисперсной системе «угольная пыль ? воздух».

Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта №15-08-06719 а.


Библиографическая ссылка

Лисаков С.А., Павлов А.Н., Сыпин Е.В., Леонов Г.В. МОДЕЛИРОВАНИЕ ОСЛАБЛЕНИЯ ОПТИЧЕСКОГО ИЗЛУЧЕНИЯ В ГАЗОДИСПЕРСНОЙ СИСТЕМЕ «УГОЛЬНАЯ ПЫЛЬ ‒ ВОЗДУХ» // Фундаментальные исследования. – 2015. – № 12-2. – С. 288-296;
URL: http://www.fundamental-research.ru/ru/article/view?id=39407 (дата обращения: 12.11.2019).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074