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



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

 


Владельцы патента RU 2562153:

Федеральное государственное унитарное предприятие "Центральный аэрогидродинамический институт имени профессора Н.Е. Жуковского" (ФГУП "ЦАГИ") (RU)

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

 

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

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

1. Наиболее известные фирмы Lavision (www.lavision.com), Dantec (www.dantecdynamics.com), а также отечественная фирма «Полис» (http://www.itp.nsc.ru) используя самые современные достижения в науке и технике (зарубежные исследования), выпускают приборы, определяющие характеристики и структуры многофазных течений. Известны способы измерения концентрации частиц (весовые нефелометры, нефелометры, основанные на поглощении излучения и др.) и капель в потоке, когда количественные измерения проводят только в отдельных точках. Трехмерную структуру течения восстанавливают путем томографического многократного фотографирования перемещающейся световой плоскости лазерного листа (напр., www.lavision.com), интенсивность излучения которой меняется во всей этой плоскости. Это существенно усложняет обработку изображения и может привести к ошибкам измерения. При этом, для построения лазерного листа используют цилиндрическую линзу или вращающиеся зеркала. Это неизбежно приводит к расхождению плоскости лазерного листа и, как следствие, к снижению точности измерений.

2. Известен способ PDPA - Particle Doppler Phase Analyzer, заключающийся в измерении скорости, размеров и концентрации дисперсной фазы в области пересечения двух лазерных лучей путем анализа изменения интерференционной картины (см., напр., С.М. Sipperley and W.D. Bachalo Triple Interval Phase Doppler Interferometry: Improved Dense Sprays Measurements and Enhanced Phase Discrimination ILASS Americas, 25th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 2013). Недостатком этого способа является пространственная ограниченность измерений - возможность получать значения параметров потока лишь в одной точке пространства.

3. Известен весовой способ определения запыленности и водности среды, заключающийся в выделении частиц из известного объема аэрозольной среды с последующим их взвешиванием. Устройство состоит из (Волков Э.П., Гаврилов Е.И., Дужин Ф.П. Газоотводящие трубы ТЭС и АЭС // М.: Энергоатомиздат. - 1987. - 280 с., Волков Э.П. Контроль загазованности атмосферы выбросами ТЭС // М.: Энергоатомиздат. 1986. - 255 с.) заборной трубки (при осаждении пыли вне газорасхода), устройства для осаждения частиц, устройства для измерения расхода отбираемых газов и средства для отсоса газов. Этот способ обладают рядом недостатков, к которым можно отнести внесение возмущения в аэрозольный поток, невозможность непрерывного автоматического контроля распределения концентрации частиц в аэрозольном потоке и длительность отбора проб частиц и их анализа.

4. Известен способ построения лазерного листа с помощью вращающегося зеркала (см., напр., www.laserpribor.ru). Метод вращающегося зеркала неизбежно приводит к расхождению лучей излучения и трудностям обработки полученных изображений.

5. Информация о числовой концентрации частиц и капель может быть также получена из анализов данных измерения способом масляной ловушки (использованного, напр., в Миллер А.Б., Потапов Ю.Ф., Токарев О.Д., Яшин А.Е. Экспериментальные исследования образования барьерного льда на обычных и наномодифицированных поверхностях авиационных материалов // Материалы XXIV научно-технической конференции по аэродинамике в пос. им. В. Володарского. - М.: Изд-во ЦАГИ. - 2013. - С. 179), заключающийся во внесении в поток импактора-пластины, покрытой маслом, в которое попадают частицы или капли, застревая в нем. Данный способ предназначен для измерения распределения частиц и капель в потоке по их размерам, а погрешность определения концентрации частиц и капель может быть существенна.

6. Информация о числовой концентрации частиц в аэрозольном потоке может быть также получена путем анализа Шлирен изображений, полученных при зондировании дисперсного потока (см., напр., William D. Bachalo, Chad Sipperley, and Gregory Payne Aircraft Icing Research: Challenges in Cloud Simulation and Characterization ILASS Americas, 23-nd Annual Conference on Liquid Atomization and Spray Systems, Ventura, CA., May 2011). Данный способ позволяет лишь качественно оценивать распределение числовой концентрации дисперсной фазы, т.к. на пути света оказывается все пространство заполненное частицами, что может привести к существенным ошибкам измерений.

Прототипом принят способ определения концентрации дисперсной фазы в аэрозольном потоке методом лазерного листа (Василевский Э.Б., Безменов В.Я., Боровой В.Я., Горелов В.А., Жилин Ю.В., Казанский Р.А., Мошаров В.Е., Чирихин А.В., Яковлева Л.В. Экспериментальное исследование течения, теплообмена и электрооптических явлений при обтекании тел сверхзвуковым аэродисперсным потоком. // «ЦАГИ - основные этапы научной деятельности 1993-2003» - М.: Физматлит, - 2003. - С. 452-457). С помощью цилиндрической линзы пучок лазерного излучения разворачивают в расходящуюся световую плоскость. В соответствии с расхождением этой плоскости находят пространственное распределение интенсивности зондирующего излучения в пространстве вдали от источника. По интенсивности рассеянного частицами излучения на фотоснимке лазерного листа определяют распределение концентрации дисперсной фазы вблизи обтекаемого тела в относительных величинах. Устройство для реализации способа содержит источник лазерного излучения, цилиндрическую линзу, фотоприемник и цифровую аппаратуру обработки изображения лазерного листа.

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

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

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

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

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

Фиг. 1. Схема оптического зондирования обтекания тела аэрозольным потоком.

Фиг. 2. Схема фотосъемки лазерного листа при оптическом зондировании обтекания тела аэрозольным потоком (ось x потока перпендикулярна плоскости фигуры).

Фиг. 3. Оптическая схема устройства построения лазерного листа постоянной ширины (вид сбоку).

Фиг. 4. Оптическая схема устройства построения лазерного листа постоянной ширины (вид спереди).

Фиг. 5. Зависимость максимального угла β расхождения светового пучка после его прохождения через цилиндрическую линзу от радиуса r этой линзы при различных значениях показателя преломления оптического стекла nλ.

Фиг. 6. Зависимость расстояния L от показателей преломления оптических стекол линз системы mλ и nλ малой цилиндрической и большой собирающей линз соответственно.

Фиг. 7. Зависимость радиуса кривизны большой собирающей линзы от угла γ при различных значениях показателя преломления материалов линз оптической системы mλ=nλ.

Фиг 8. Визуализация методом лазерного листа газопылевого и суспензионного потоков.

Фиг. 9. Схема использования информации при интегрировании вдоль пути распространения излучения. Ось x потока перпендикулярна плоскости фигуры.

Фиг. 10. Верификация способа обработки изображения лазерного листа.

Фиг. 11. Схема количественного определения пространственного распределения концентрации частиц в аэрозольном потоке.

Фиг. 12. Индикатрисы рассеивания, рассчитанные по теории Ми. , a - радиус частицы, λ=0.63 мкм - длина волны зондирующего излучения.

Рассмотрим использование лазерного листа при зондировании осесимметричного обтекания затупленного тела сверхзвуковым аэрозольным потоком (фиг. 1) с пространственным распределением в единице объема (в дальнейшем, концентрации) n(x, y, z) одинаковых шаровых частиц из известного материала и радиуса а: пример геометрических и физических параметров указаны в табл. 1 и 2 соответственно.

На фиг. 1 изображены 1 - плоскость лазерного листа, 2 - направление излучения 3 - направление аэрозольного потока, 4 -обтекаемое тело, 5 - область, охватываемая фотоприемником, 6 - объектив фотоприемника (светочувствительная матрица). Плоскость объектива фотоприемника параллельна плоскости лазерного листа, плоскость лазерного листа проходит через ось симметрии аэрозольного потока x и принадлежит плоскости z=0. Каждая частица обладает своим коэффициентом рассеивания Qsca, коэффициентом поглощения Qa и коэффициентом ослабления (экстинкции) Qext=Qsca+Qa падающего излучения интенсивностью I0 с длиной волны λ, которые вычисляются с помощью известной теории Ми (или формуле Релея в случае малых частиц) в зависимости от электромагнитных свойств материала частиц и их размеров. В пренебрежении многократным рассеиванием энергии излучения получим формулу, связывающую излучение в каждой ячейке фотоснимка световой плоскости (лазерного листа), и концентрацию частиц n в зондируемом объеме.

Интенсивность излучения I A i ( x , y ) в каждой точке A(x, y) плоскости z=0 лазерного листа (фиг. 1) пропорциональна концентрации частиц n(x, y, 0) в малой окрестности этой точки, толщине ножа h, а также коэффициенту рассеивания излучения πa 2Qsca(a) при заданной длине волны. Здесь индекс i означает incident (падающий).

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

Каждая точка плоскости лазерного листа взаимно-однозначно отображается на фотоснимке. Из каждой светящейся (в результате рассеивания) ячейки световой плоскости листа лучистая энергия, затухая, попадает, пройдя через линзу объектива (фиг. 2), на соответствующую точку S(xM, yM) матрицы фотоприемника, а выражение для интенсивности сигнала Ie(xxM(х), yM(y)) в произвольной точке S(xM, yM) на фотоснимке будет следующим:

где s - путь, пройденный светом от точки плоскости лазерного листа A(x, y), до точки S(xM(x), y0(y)), плоскости светочувствительной матрицы: фотоснимка. Малая площадка A1A2 и отрезок AL образуют телесный угол Ω, внутри которого происходят потери лучистой энергии (фиг.2), попавшей в фотоприемник из элемента A1A2 плоскости листа.

На фиг. 2 изображены 1 - лазерный лист, 2 - направление излучения, 4 - обтекаемое тело, 6 - светочувствительная матрица фотоприемника: изображение объекта, 7 - линза фотоприемника. Весь путь от точки A до точки S определяется выражением AS = ( x x M ) 2 + ( y y M ) 2 + Z 2 , где Z - расстояние от плоскости лазерного листа до фотоприемника. Концентрация частиц при интегрировании вдоль пути s определяется следующим выражением n(s)=n(ξ (s), η(s), ζ(s)). Входящие в него величины находятся из геометрии фиг. 1: ξ(s)=s sinεcosφ, η(s)=ssinεsinφ, ζ(s)=scosε. Из геометрии (Фиг. 1) также находим входящие в них углы: t g ε = x 2 + y 2 / Z , tgψ=y/x.

На Фиг. 3 и Фиг. 4 изображены схемы (фиг. 3 - вид сбоку и фиг. 4 - вид спереди) устройства для построения лазерного листа. Устройство для построения лазерного листа 1 постоянной ширины состоит из цилиндрической линзы (постоянного радиуса кривизны r) 8, большой собирающей линзы с переменным радиусом кривизны 9, диафрагмы 10, щели 11. Пучок лазерного излучения 12, пройдя сквозь описанную выше оптическую систему, преобразуется в плоскопараллельный поток лучей 13 (имеющий элипсообразную форму в своем сечении), который после обрезания диафрагмой преобразуется в тонкий лазерный лист постоянной ширины.

Запишем известный закон Снеллиуса для преломления лучей света, прошедших сквозь линзы (см. фиг. 3):

Здесь nλ и mλ - показатели преломления материалов, из которых сделаны первая и вторая линзы соответственно. Максимальные значения углов α и β определяются диаметром входного лазерного луча d: и радиусом кривизны r первой (цилиндрической) линзы:

На фиг. 5 показаны зависимости значений максимального угла β при различных значениях показателя преломления оптического стекла цилиндрической линзы.

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

O2D+FD-O1F=O2D+FD-[O1B+FB]=O1O2=L, где L - расстояние между центрами кривизны цилиндрической и большой собирающей линз,

или

Здесь углы α и β зависят от угла γ:

Величину L находим исходя из диаметра входного пучка d и ширины лазерного листа D:

Или

Таким образом, L зависит только от показателей преломления линз nλ и mλ.

На фиг. 6 показана зависимость значений длины L от показателей преломления оптических стекол, из которых изготовлены линзы.

А радиус кривизны выражается следующей формулой:

Здесь nλ - показатель преломления материала, из которого изготовлены линзы, R(γ) - радиус кривизны линзы; L - расстояние между центрами кривизны линз.

На фиг. 7 приведена зависимость радиуса R(γ) кривизны большой собирающей линзы при различных значениях показателя преломления материала, из которого изготовлены линзы: nλ=1.1-2.0.

Поток параллельных лучей, полученный при прохождении света через большую собирающую линзу 9, имеет вид эллипса 13 в плоском сечении (фиг. 3, 4). Для создания плоского лазерного листа 1 (фиг. 3) этот эллипс нужно «обрезать» пропустив его через диафрагму 10 (фиг. 3, 4) со щелью 11 (фиг. 3, 4). Меняя ширину щели 11, подбирают оптимальную толщину плоскости лазерного листа для эффективного оптического зондирования многофазных потоков.

На фиг. 8 показана визуализация методом лазерного листа газопылевого и суспензионного потоков.

Опишем метод обработки изображения. Положим Ie - поле интенсивности излучения на фотоснимке, который получен при фотографировании плоскости лазерного листа, зондирующего обтекание тела аэрозольным потоком (фиг. 1). Разобьем фотоснимок - расчетную область на Nx и Ny ячеек вдоль осей x и y соответственно. Пусть i, j - номера ячеек вдоль координатных осей x и y соответственно (i ∈ [1, Nx+1], j ∈ [1, Ny+1]), a Ie(i, j) - значение интенсивности излучения в ячейке с номерами i, j. Требуется найти пространственное распределение концентрации частиц вблизи затупленного тела при помощи описанного выше фотоснимка: решить обратную задачу. Существуют методы решения таких задач (Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач // М.: Наука. Гл. ред. физ.-мат. лит. - 1986. - 288 с.; Алифанов О.М. Обратные задачи теплообмена // М.: Машиностроение. - 1988. - 280 с.; Бек Дж., Блэкуэлл Б., Сент-Клэр Ч. Некорректные обратные задачи // М.: Мир. - 1989. - 312 с.), основанные, главным образом, на использовании дополнительной априорной информации и последовательном (итерационном) приближении к решению.

Основные принципы разработанного способа обработки изображения состоят в том, что

1. компьютерную обработку изображения производят вдоль пути распространения энергии зондирующего излучения на фотоснимке;

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

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

На фиг. 9 изображена схема использования информации при обработке изображения лазерного листа, где плоскость лазерного листа 1, направление падающего излучения 2, направление сканирования при обработке изображения 3, направление излучения, попадающего в фотоприемник 4, линии равных значений концентрации дисперсной фазы 5, линза фотоприемника 6, светочувствительная матрица: изображение объекта 7. Опишем численный метод решения обратной задачи, подставив (1) в (2) и записав (2) в следующем виде:

Здесь β=πa 2Qext - ослабление (экстинкция) излучения, коэффициент α связывает интенсивность излучения на снимке Ie с концентрацией дисперсной фазы n в потоке при оптическом зондировании аэрозольного потока.

В верхней области лазерного листа энергия света не ослабла в результате рассеивания вдоль оси y (фиг. 9), т.к.

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

или

Это значение меньше того значения концентрации частиц, которое требуется получить и которое впоследствии будет найдено методом последовательных приближений: методом итераций. Множитель 1 2 нужен для того, чтобы найденное в первом приближении значение концентрации частиц не превосходило искомого значения концентрации частиц, которое требуется найти. Отладка алгоритма показала, что в противном случае решение не будет сходиться с ростом числа приближений. Во втором приближении (итерации) концентрация частиц вычисляется, используя значение концентрации частиц в данной ячейке в первом приближении n ( i ,2 ) 1 . При этом ослабление излучения между первой и второй ячейкой равно произведению β на шаг разбиения расчетной области по вертикали Δy и на среднее значение концентрации частиц между этими ячейками. Концентрация частиц во втором приближении равна

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

Устремляя число итераций p→∞, получаем значение концентрации дисперсной фазы n в ячейке (i, 2) с точностью до численной аппроксимации. Опыт предшественников (Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач // М.: Наука. Гл. ред. физ.-мат. лит. - 1986. - 288 с.; Алифанов О.М. Обратные задачи теплообмена // М.: Машиностроение. - 1988. - 280 с.) свидетельствует об эффективности итерационных методов решения обратных задач. Продолжая действовать аналогичным образом, находят значение концентрации дисперсной фазы в третьей ячейке по вертикали, заметив при этом, что нумерация ячеек отсчитывается от места попадания излучения в расчетную область, т.е. сверху вниз:

в первом приближении,

- во втором и последующих приближениях.

Проводя сканирование по направлению распространения излучения s аналогичным образом, уловив закономерность, находят значения концентрации дисперсной фазы в остальных js=Ny +1-j ячейках:

в первом приближении,

- во втором и последующих приближениях.

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

В верхней области лазерного листа интенсивность излучения не ослабла в результате рассеяния как вдоль оси y, так и вдоль оси z (см. фиг. 9), т.к. оптические пути (вдоль y и z), пройденные излучением, равны нулю

Интеграл вдоль оси z равен нулю, т.к. при расстоянии от оси x, большем значения H, дисперсная фаза отсутствует (см. фиг. 9). Поэтому концентрация частиц в верхней ячейке расчетной области равна или В следующей ячейке, ниже на малую величину Δy концентрация частиц в первом приближении (первой итерации) равна

Где Δ12 - первая и единственная (см. фиг. 4) ячейка вдоль оси z, которая соответствует второй ячейке сканирования

Здесь

Полученное значение концентрации частиц в этой ячейке меньше того значения концентрации частиц, которое требуется получить и которое впоследствии находят методом последовательных приближений: методом итераций. Множитель 1 2 нужен для того, чтобы найденное в первом приближении значение концентрации частиц не превосходило значение концентрации частиц, которое требуется найти. Отладка алгоритма показала, что в противном случае решение не будет сходиться с ростом числа приближений. Во втором приближении (итерации) концентрацию частиц вычисляют, при этом используют значение концентрации частиц в этой ячейке в первом приближении n ( i ,2 ) 1 . При этом ослабление излучения по вертикали равно произведению β на шаг разбиения расчетной области по вертикали Δy и на среднее значение концентрации частиц между этими ячейками, а ослабление на пути z равно произведению β на шаг разбиения расчетной области вдоль z Δ12 и на среднее значение концентрации частиц между этими ячейками. Концентрация частиц во втором приближении равна

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

Устремляя число итераций p→∞ получают значение концентрации дисперсной фазы n в ячейке (i, 2) с точностью до численной аппроксимации. Продолжая действовать аналогичным образом, находят значение концентрации дисперсной фазы в третьей ячейке по вертикали:

- в первом приближении,

- во втором и последующих приближениях.

Здесь Δ j s ,k = ( tanα j s k 1 tanα j s k ) y j s - величина ячейки ВС (фиг. 9), y j s = H ( j s 1 ) Δy ,

Таким образом, величина произвольной ячейки ВС равна ( H ( j 1 ) Δ y ) Δ j s , k . Проводят сканирование по направлению распространения излучения s и используют информацию о симметрии относительно оси x, аналогичным образом, уловив закономерность, находят значения концентрации дисперсной фазы в остальных js=Ny+1-j ячейках:

в первом приближении,

во втором и последующих приближениях.

Отладка алгоритма показала, что итерационный процесс сходится (максимальная ошибка менее 0.0001%) не более чем за 10 итераций в диапазоне параметров, соответствующем выполнению закона Бугера. В противном случае, итерационный процесс может привести к бесконечно большой ошибке. При значительном ослаблении излучения где характерная длина пробега излучения - толщина сжатого слоя вблизи обтекаемой сферы) метод может дать значительную ошибку, т.к. известный экспоненциальный закон ослабления Бугера-Ламберта-Бера не выполняется.

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

Покажем иллюстративно работу алгоритма в численном эксперименте.

На фиг. 10 показана проверка способа обработки изображения: 1 - распределение концентрации частиц ( n ( x = 1.25 R , y 2 + z 2 ) ) (нормированных на концентрацию частиц n0 у точки торможения), 2 - зарегистрированный фотоприемником оптический сигнал IeМ, yM), 3 - распределение концентрации частиц, полученное в результате обработки «изображения» 2, 4 - ослабление зондирующего излучения в плоскости лазерного листа, 5 - ослабленное зондирующее излучение на пути от листа до фотоприемника. В результате численного эксперимента, в котором по полю интенсивности излучения на «фотоснимке» (полученному в результате решения прямой задачи (фиг. 10, кривая 2): расчет интенсивности излучения 1е на фотоснимке по заданному распределению концентрации частиц n (фиг. 10, кривая 1), которое качественно описывает поведение дисперсной примеси в экспериментах. (Василевский Э.Б., Безменов В.Я., Боровой В.Я., Горелов В.А., Жилин Ю.В., Казанский Р.А., Мошаров В.Е., Чирихин А.В., Яковлева Л.В. Экспериментальное исследование течения, теплообмена и электрооптических явлений при обтекании тел сверхзвуковым аэродисперсным потоком. // «ЦАГИ - основные этапы научной деятельности 1993-2003» - М.: Физматлит, - 2003. - С. 452-457), в котором методом лазерного листа определялось поле концентрации частиц (фиг. 10, кривая 3) с помощью описанного выше метода и сравнивалось с заданным распределением (фиг. 10, кривая 1). Проведенный расчет (фиг. 10) и отладка описанного выше алгоритма свидетельствуют о ничтожности относительной погрешности определения распределения концентрации дисперсной фазы в аэродисперсном потоке по интенсивности рассеянного ими зондирующего излучения; тем не менее, результат существенно зависит от мелкости разбиения расчетной области: числа ячеек. При сканировании в направлении распространения излучения s необходимым условием сходимости решения обратной задачи является подбор шага разбиения расчетной области таким образом, что при монотонном возрастании пройденного светом пути s(i, j+1)>s(i, j) в плоскости лазерного листа.

Калибровку измерительной системы осуществляют путем сопоставления яркости рассеянного излучения на фотоснимке лазерного листа и измеренной концентрации дисперсной фазы весовым нефелометром. На фиг. 11 представлена схема количественного определения пространственного распределения концентрации частиц в аэрозольном потоке. 2 - направление распространения зондирующего излучения, 3 - направление аэрозольного потока, 4 - обтекаемое тело, 17 - сопло, 18 - контрольные объемы, в которых измеряют распределение частиц по размерам в аэрозольном потоке, 19 - построитель лазерного листа постоянной толщины, 20 - область потока, предназначенная для калибровки измерительного комплекса, 21 - весовой нефелометр с датчиком скорости потока. В верхней части оптически зондируемой области 20 (там, где лучистая энергия меньше всего ослабла в результате ее поглощения и рассеивания частицами) устанавливается весовой нефелометр 21 с датчиком расхода потока. Частицы ловятся в нефелометр и взвешиваются. Зная площадь входного сечения нефелометра, расход потока и массу попадающих в него частиц в единицу времени, измеряют концентрацию частиц в абсолютных единицах. Этой концентрации соответствует яркость области 20 на фотоснимке лазерного листа, зондирующего обтекание тела аэрозольным потоком. Коэффициент ослабления зависит только от материала частицы, ее радиуса а и длины волны излучения λ: β=πa 2Qext(a, λ).

Распределение частиц по размерам находим с помощью ореольного нефелометра (Лагунов А.С., Байвель Л.П., Гусев Б.А. Определение дисперсного состава аэрозолей методом рассеяния света под малыми углами // Журнал прикладной спектроскопии. - Том XI. - Вып. 1. - 1969. - С. 98-103.; Рудаков В.П., Стасенко А.Л., Флаксман Я.Ш. Определение массового спектра частиц в газодисперсном потоке методом малоуглового рассеяния // Ученые записки ЦАГИ. - Т. XXIV. - №2. - 1993. - С. 114-122. Горчаков Г.И., Токарев О.Д. Оценка микроструктуры аэрозоля по данным измерений ореольных индикатрис рассеяния света // Сб. Атмосферная радиация и фотометрия. - Томск. - 1988. - С. 22-25; Любовцева Ю.С. Статистические характеристики светового ореола в тумане // Изв. АН СССР. - Физика атмосферы и океана. - 1971, т. 7. - №4. - С. 465-467; Голиков В.И. Прибор для измерения микроструктуры облаков и туманов методом малых углов // Труды Главной геофизической обсерватории им. А.И. Воейкова. Л.: 1964. - Вып. 152. - С. 142-159), основанного на известном методе малоуглового рассеивания (по индикатрисе рассеивания), определяется распределение частиц и капель по их размерам. Или с помощью известного метода PDPA (particle Doppler Phase Analyser), основанного на известных физических эффектах интерференции и Допплера ((см., напр., С.М. Sipperley and W.D. Bachalo Triple Interval Phase Doppler Interferometry: Improved Dense Sprays Measurements and Enhanced Phase Discrimination ILASS Americas, 25th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 2013).

Коэффициенты рассеивания, ослабления и поглощения, соответственно, находят из известных формул Густава Ми (Mie G. Beiträge zur Optik Medien speziell kolloidaler Metallösungen. - Ann. Phys. 25. - P. 377-445. - 1908) в зависимости от физических свойств материала частиц аэрозоля и дифракционного параметра x=2πa/λ, зависящего от радиуса частицы а и длины волны λ зондирующего излучения.

Здесь i = 1 .

Зависимость индикатрис рассеянного частицами излучения от их размеров при заданном значении длины волны изображена на фиг. 12. Рассчитав светорассеивающие характеристик отдельных частиц с помощью изложенных выше формул Ми для коэффициентов рассеивания в уравнении (3) остается только одно неизвестное - нормировочный коэффициент α:

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

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

Линзы могут быть изготовлены на заводах оптического стекла. Большая линза переменного радиуса кривизны может быть изготовлена из параболической собирающей линзы или собирающей линзы постоянного радиуса кривизны путем шлифовки на токарном станке. Линзы и диафрагма на соответствующих держателях могут быть установлены на оптическую скамью. После преобразования пучка лазерного излучения в плоскопараллельную световую плоскость последнюю направляют на исследуемую область. С помощью цифрового аппарата регистрируют изображение рассеянного в плоскости лазерного листа излучения. Обработка изображения цифрового фотоснимка лазерного листа с использованием описанного выше алгоритма может быть реализована с помощью ЭВМ.

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

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

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



 

Похожие патенты:

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

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

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

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

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

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

Изобретение относится к области техники, а именно автоматизации измерений при анализе взвешенных наночастиц в газах. Для этого используют устройство для определения спектра размеров взвешенных наночастиц в газах, содержащее размещенные по ходу анализируемого потока газа входное сопло с каналами подачи; диффузионные батареи сетчатого типа для пропускания аэрозольных частиц определенного размера; укрупняющее устройство конденсаторного роста; счетный объем; вакуумный насос; температурные датчики, нагреватель, охладитель и микроконтроллер для управления процессами нагревания и охлаждения в укрупняющем устройстве конденсаторного роста; оптическую систему, включающую импульсный источник излучения, осветитель и объективы для фокусировки оптического излучения в области счетного объема потока частиц и формирования изображений на матрице ПЗС; аналогово-цифровой преобразователь и ЭВМ для управления микроконтроллером термостатирования, ваккумным насосом и обработки шести изображений укрупненных частиц для анализа спектра их размеров.

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

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

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

Предложенный способ позволяет измерять распределение по фракциям и концентрации твердых и жидких частиц аэрозоля в интервале размеров частиц: от 0,8 мкм до 2 мкм, от 2 мкм до 5 мкм, от 5 мкм до 10 мкм и более 10 мкм при помощи полупроводниковых кондуктометрических сенсоров по изменению проводимости. При применении заявляемого способа измерения концентрации частиц аэрозоля в реальном времени и устройства для его осуществления, состоящего из забора пробы из окружающей среды, распределения по дисперсности исследуемого потока аэрозоля по соответствующим камерам и измерения концентрации аэрозоля в каждой камере, измерения проводятся полупроводниковыми кондуктометрическими сенсорами по изменению их проводимости после предварительного прогрева, для чего проводится предварительная калибровка сенсоров по их проводимости на атмосферном воздухе и на воздухе, содержащем аэрозоль с известной концентрацией и дисперсностью, которая принимается постоянной вплоть до окончания измерений, причем калибровка осуществляется для каждого сенсора, далее с учетом калибровки значения проводимостей по заданному алгоритму пересчитываются в концентрации. Техническим результатом является снижение стоимости измерения и потребления энергии, а в части устройства - уменьшение его массогабаритных характеристик. 2 н. и 2 з.п. ф-лы, 1 ил.

Изобретение относится к исследованию физико-механических свойств сталей и сварных соединений и может применяться в различных отраслях промышленности. Сущность: по окончании процесса воздействия на образец ударным изгибом предварительно готовят микрошлиф образца. С помощью оптического микроскопа формируют изображение микроструктуры участка разрушения, по которому определяют величину пакетов бейнитной структуры, и находят параметрические характеристики вязкого и хрупкого разрушения по следующим зависимостям: dя=0,1 dб dфтс=0,6-0,8 dб, для трубных сталей класса прочности до К60, dфтс=0,4-0,6 dб, для трубных сталей класса прочности свыше К60 включительно, где: dя - средний размер ямки в участке вязкого разрушения (мкм), dфтc - средний размер фасетки транскристаллитного скола в участке хрупкого разрушения (мкм), dб - средний размер пакета бейнита (мкм). Технический результат: сокращение времени исследования и упрощение процесса оценки искомых характеристик. 18 ил.

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

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

Изобретение относится к области оптических методов измерения физико-химических характеристик аэрозольных сред и может быть использовано при разработке лидарных комплексов для дистанционного контроля дисперсного состава аэрозольных облаков стойких токсичных химикатов (ТХ) при возникновении запроектных аварий в местах хранения и уничтожения химического оружия (УХО) и на других химически опасных объектах. В способе проводится зондирование полидисперсного аэрозольного облака ТХ многочастотным лазерным излучением ультрафиолетового, видимого и инфракрасного диапазонов спектра и регистрируются интенсивности сигналов обратного упругого аэрозольного рассеяния. В процессе хранения ТХ осуществляется контроль их оптических констант (коэффициента преломления и показателя поглощения). По результатам спектральных измерений создается база данных характеристик аэрозольного рассеяния ТХ на основе многопараметрических рядов, включающих относительные характеристики обратного аэрозольного рассеяния с использованием инструментально измеренных значений мнимой и действительной частей комплексного показателя преломления ТХ, а также медианного диаметра и дисперсии распределения логарифмически нормального закона распределения аэрозоля ТХ по дисперсному составу. При этом контроль дисперсного состава аэрозолей ТХ осуществляют в рамках теории распознавания образов по минимальному значению меры близости сигналов аэрозольного рассеяния, полученных в эксперименте с помощью дистанционного средства, и данных многопараметрических рядов в составе базы данных средства локации. Изобретение обеспечивает дистанционный контроль размеров тонкодисперсных аэрозолей стойких ТХ с логарифмически нормальным законом распределения частиц по дисперсному составу для оценки масштабов и последствий аварийных выбросов ТХ на объектах УХО. 3 табл.

Изобретение относится к технике измерений, может использоваться в автомобильной, сельскохозяйственной, авиационной, нефтеперерабатывающей и других отраслях промышленности, где необходимо проводить оперативный анализ качества моторного масла. Устройство анализа загрязненности моторного масла двигателя внутреннего сгорания дисперсными частицами включает кювету с чистым маслом, измерительный канал анализа угарных частиц, расположенный на высоте минимального уровня масла в картере, и измерительный канал анализа металлических частиц, расположенный внизу масляного поддона картера двигателя. Также устройство включает лазер в качестве источника зондирующего излучения, три смотровых окна, три светоделителя (полупрозрачных зеркала), световую ловушку, три объектива, датчик температуры и три ультразвуковых излучателя каналов анализа угарных, металлических частиц. Кроме того, устройство также включает эталонный канал, усилитель, четыре аналого-цифровых преобразователя, цифроаналоговый преобразователь, генератор ультразвуковых колебаний, коммутатор, электронно-вычислительную машину. Также устройство дополнительно содержит три ПЗС-матрицы и три DSP-процессора. Техническим результатом является повышение точности измерения параметров угарных и металлических частиц, а также повышение информативности данных для оценки концентрации, размера и формы взвешенных металлических и угарных дисперсных частиц, находящихся в масле, в частности дает возможность контролировать качество работы двигателя, оставшийся ресурс работы масла до его замены. 1 ил.

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

Изобретение относится к области исследования частиц с помощью ИК спектроскопии, в частности к методам экспресс-анализа полимерных композитов. В способе определения ориентации анизометричных частиц наполнителя в объеме полимерной матрицы при выполнении условия |nМ-nН|>0, где nМ и nН - показатели преломления матрицы и наполнителя соответственно, производится регистрация ИК спектров пропускания при облучении композитов под разными углами источником ИК излучения. Далее, преобразованием полученных спектров в кривые распределения частиц наполнителя по размерам и идентификацией экстремумов кривых распределения частиц наполнителя по размерам определяют угол ориентации частиц в объеме полимерной матрицы. Техническим результатом является разработка ИК спектроскопического экспресс-способа определения ориентации анизометрических частиц наполнителя в объеме полимерной матрицы. 2 табл., 4 ил.
Наверх