Development and research of algorithms for estimating optical element parameters operating for transmission of infrared radiation

Cover Page

Cite item

Full Text

Abstract

In this article we consider the problem of estimating the absorption coefficient of the optical element having the shape of a circular disk and the irradiated laser beam, by measuring temperature at the periphery, and the heat exchange coefficient with the environment. To do this, we considered the quality criterion function, which has a trading structure. The solution of the direct thermal conductivity problem in the form of an infinite Fourier-Bessel series and the inverse problem is obtained using the coordinate descent method, the fastest descent method, the modified Golden section method, and the Gauss method with two step corrections. Also, a comparative analysis was made on the criteria of accuracy and execution time of the described methods based on the results of accurate and noisy measurements.

Full Text

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

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

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

Описание физического эксперимента

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

Оптический элемент помещен в вакуум-камеру, поскольку в воздухе очень сильно действие коэффициента теплообмена. Он имеет большую величину порядка 102 ~ 5×103 и является нестабильным. Поэтому есть необходимость уменьшить влияние коэффициента теплообмена. В вакуум-камере ситуация стабильная, влияние изменения окружающей среды минимально и погрешность меньше. Вакуум-камера поддерживает стабильный вакуум и параметр теплообмена мало меняется (значение  становится близкой к константе). В дальнейшем он предполагается постоянной величиной.

Лазерный луч направляют через окно вакуум-рна оптический элемент. Оптический элемент частично поглощает лазерное излучение, вследствие чего нагревается. В силу теплопроводности происходит распространение тепла от зоны засветки до периферии. На периферии ( r=R) внутри оправы установлены датчики температуры, обладающие высокими теплоизолирующими свойствами. С помощью датчиков температуры проводят измерения теплового поля в процессе эксперимента.

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

Оптический элемент представляет собой диск с параллельными гранями кругового сечения радиуса R.  из однородного материала толщиной l. Также он характеризуется коэффициентом поглощения β.

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

Математическая постановка прямой и обратной задачи

Математическая модель теплового процесса, протекающего в диске, описывается следующей системой [5]:

 

сut=kTΔru2αlu+Irβ,       u|t=0=0,                                               urr=R=0,                                               0<rR, 0tT;0rR;                     0tT;                     

 

где u=ur,t – разница между температурой оптического элемента и температурой окружающей среды,

α – коэффициент теплообмена,

kT – коэффициент теплопроводности материала диска,

l – толщина оптического элемента,

c=cyρ – объемная теплоёмкость,

Ir – интенсивность лазерного излучения.

 

Рис. 1. Вид установки

 

Интенсивность лазерного излучения имеет гауссовское распределение по радиусу диска:

Ir=Pπa2era2,                                         

где  a – эффективный радиус лазерного пучка,

P – мощность излучения.

Решение задачи запишется в виде:

ur,t=n=0+βPJ0μnrReμn2a24R2πR2J02μnμn2R2+2αkTlkT1eμn2R2+2αkTlkTct.

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

Суть решаемой обратной задачи теплопроводности состоит в вычислении значений параметров уравнения теплопроводности (коэффициента поглощения β и коэффициента теплоотвода  α с граней оптического элемента во внешнюю среду) по результатам измерений температуры на периферии оптического элемента при r=R.

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

В качестве критерия оптимальности оценок параметров использовался следующий среднеквадратический критерий качества [6]:

Gα,β=1Kk=1KuR, tk, α,βϑk2minα,β,

где uR, tk, α,β – решение задачи теплопроводности в моменты времени tk при некоторых значения α,β.

ϑk – результаты измерения температуры периферии элемента в моменты времени  tk, k=1,K¯.

Предполагается, что после включения лазера в течение времени T производится регистрация температуры с помощью тепловых датчиков, размещенных на поверхности r=R оптической пластины. При вычислении значения критерия качества Gα,β однократно решается прямая задача теплопроводности при значениях параметров α и β уравнения, выбираемых в соответствии с алгоритмом оптимизации. 

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

Исследование функции критерия качества

Для расчета α и β уравнения теплопроводности проведем исследование функции цели Gα,β. Здесь при проведении вычислительных экспериментов по исследованию критерия оптимальности в качестве «измерений» ϑk  брались результаты расчета теплового поля при значениях параметров α=0,005 и β=0,008 с отсутствием шума (рис. 2).

Исследование функции цели Gα,β проводилось методом сечений плоскостями α=const, для конкретного α из набора значений вычисляется функция Gα,β  для различных β. Данное исследование показало, что функция цели имеет форму криволинейного оврага [7]. Далее была осуществлена минимизация функции Gα,β.

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

 

Рис. 2. Графики сечений функции цели плоскостями

 

Рис. 3. График линии дна оврага функции

 

Рис. 4. Оптиматизация методом покоординатного спуска

 

Метод покоординатного спуска

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

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

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

Рассмотрим работу покоординатного метода для исходной задачи из начальной точки α=0,01 и β=0,01 (рис. 4).

Метод наискорейшего спуска

Формализация алгоритма наискорей-шего спуска выглядит следующим образом:

Шаг 1. Задаются начальное приближение для α, β и точность ε.

Шаг 2. Находим значение градиента по формулам: Gα=Gu×uα  и , где Gβ=Gu×uβ

 

uβ=n=0+PJ0μnrReμn2a24R2πR2J02μnμn2R2+2αkTlkT1eμn2R2+2αkTlkTct;

uα=n=0+2βPJ0μnrRteμn2a24R2πR2J02μnμn2R2+2αkTlkTcleμn2R2+2αkTlkTct

2βPJ0μnrReμn2a24R2πR2J02μnμn2R2+2αkTl2kT2l1eμn2R2+2αkTlkTct.                   

Шаг 3. Пересчитываем α, β по формулам:

αk+1=αkλkGα,

βk+1=βkλkGβ, где λk увеличивается в p раз, пока значение функции цели уменьшается.

Шаг 4.

Если выполняется αk+1αkαk+1+α*ε и βk+1βkβk+1+β*ε, где в качестве α* и β* берется малое значение 107, то возвращаем текущие значения αk+1, βk+1.

Иначе возврат к шагу 3.

Продемонстрируем работу метода наискорейшего спуска (рис. 5).

 

Модифицированный метод золотого сечения

Алгоритм поиска оптимальных значений параметров уравнения теплопроводности для нашей задачи является адаптацией одномерного метода золотого сечения к задаче двумерной оптимизации. Оптимальные значения оцениваемых коэффициентов ищутся как по одной, так и по другой переменной, однако здесь не идёт речь о классической покоординатной оптимизации. Оптимизация осуществлена вдоль линии оврага. Сначала по правилу золотого сечения выбираются точки на оси a без вычисления значений критерия качества. Затем для каждой точки на оси a по одномерному методу золотого сечения осуществляется спуск на дно оврага функции Gα,β. При этом оптимизация проводится по параметру с использованием опять же одномерного метода золотого сечения. В полученных таким образом точках на дне оврага вычисляются значения критерия Gα,β. В результате сравнения этих значений производится сужение интервалов неопределённости вдоль линии дна оврага и, следовательно, одновременно происходит сужение интервалов неопределенности как по параметру a, так и по параметру β [7]. Продемонстрируем работу метода золотого сечения (рис. 6).

 

Рис. 5. Оптиматизация методом наискорейшего спуска

 

Рис. 6. Оптиматизация модифицированным методом золотого сечения

 

Рис. 7. Оптиматизация методом Гаусса

 

Рис. 8. Оптиматизация методом Гаусса с корректировкой шага, основанной на аппроксимации параболой

 

Рис. 9. Оптиматизация методом Гаусса с корректировкой шага, основанной на методе наискорейшего спуска

 

Метод Гаусса

Метод Гаусса относится к семейству градиентных методов, как и метод наискорейшего спуска. Итерационный метод Гаусса будет иметь следующий вид (рис. 7):

 

Однако неудачно выбранное начальное приближение может значительно увеличить число итераций или даже вызвать зацикливание. Можно заметить, что метод Гаусса стреляет с большим шагом, из-за чего иногда случаются «перелеты» и растет время работы программы. Этого можно избежать с помощью добавления корректировки для шага «выстрела».

Первая корректировка основана на аппроксимации параболой. Её суть заключается в том, чтобы уменьшать или увеличивать шаг в 2 раза и держать в памяти компьютера три оптимальные точки функции Gx¯. Из графика критерия качества заметно, что функция цели похожа на параболу, экстремум которой будет являться оптимальным значением. Для аппроксимации воспользуемся уравнением Gx¯=ax¯2+bx¯+c. Подствляя полученные программой точки, решаем систему из трех уравнений и находим коэффициенты a, b, c, после чего находим экстремум полученной параболы и сравниваем значение функции цели в точке минимума Gx¯,  с предложенной методом Гаусса, и выбираем наименьшую. Продемонстрируем работу метода Гаусса с применением данной корректировки шага из точки α=0,01, b=0,02 (рис. 8).

В качестве эксперимента была также реализована вторая корректировка шага, которая применялась ранее в методе наискорейшего спуска. Продемонстрируем работу метода с такой корректировкой шага из точки α=0,02, b=0,01 (рис. 9).

Обе корректировки не дали большого выигрыша во времени, но позволили методу Гаусса быстрее добраться до линии оврага.

Сравнительный анализ методов

В результат проведения вычислительных экспериментов на «точных» измерениях все методы позволили достичь оптимума с высокой точностью, хотя за разное время. Стоит отметить, что различие между полученными значениями α, β и αtrueβtrue  не превосходили 0,0005% при использовании модифицированного метода золотого сечения и 0,0001% при использовании метода Гаусса. Поскольку при проведении натурных экспериментов всегда присутствуют ошибки измерений, важно было проверить работу реализованных методов оптимизации на зашумленных измерениях. Для имитации зашумленных измерений использовалась аддитивная добавка Гауссовского шума с нулевым средним и среднеквадратическим отклонением, равным σ2. В данной работе было рассмотрено σ в диапазоне 0,1;0,3, что соответствует параметрам современных динамических датчиков температуры. Формула критерия оптимальности оценок параметров с учетом шума v:

Gx,ξ=ξuxTξuxminx,           

где x=αβ,ξ= utruex+vux=uR,t1,α,βuR,t2,α,βuR,tM,α,β.

                                                        

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

  

Таблица 1

Сравнение средней погрешности разных методов оптимизации на зашумленных измерениях

Метод

Погрешность по α  (εα ,%)

Погрешность по β , ( εβ,%)

σ=0,1

σ=0,2σ=0,3

σ=0,1

σ=0,2σ=0,3

Метод золотого сечения

   0,73

  2,04

  3,25

   0,46

  0,99

  1,87

Метод Гаусса

   0,59

  1,74

  2,91

   0,38

  0,85

  1,55

Таблица 2

Сравнение среднего времени работы разных методов оптимизации на зашумленных измерениях

Используемый метод

    σ=0,1

     σ=0,2

      σ=0,3

Метод покоординатного спуска

450 секунд

500 секунд

600 секунд

Метод наискорейшего спуска

250 секунд

300 секунд

400 секунд

Метод золотого сечения

115 секунд

115 секунд

115 секунд

Метод Гаусса

10 секунд

15 секунд

18 секунд

Метод Гаусса с корректировкой шага

3 секунды

5 секунд

8 секунд

 

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

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

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

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

×

About the authors

Alyona Konstantinovna Shlyakhova

Samara University

Author for correspondence.
Email: gastingz@gmail.com

graduate student of the Faculty of Informatics

Russian Federation, 443086, Russia, Samara, Moskovskoye Shosse, 34

Alexander Alexandrovich Degtyarev

Samara University

Email: aadegt@gmail.com

assistant professor of the Department of Technical Cybernetics of the Samara University

Russian Federation, 443086, Russia, Samara, Moskovskoye Shosse, 34

Supplementary files

Supplementary Files
Action
1. JATS XML

Copyright (c) 2020 Proceedings of young scientists and specialists of the Samara University

Creative Commons License
This work is licensed under a Creative Commons Attribution-ShareAlike 4.0 International License.

Proceedings of young scientists and specialists of the Samara University

ISSN 2782-2982 (Online)

Publisher and founder of the online media, journal: Samara National Research University, 34, Moskovskoye shosse, Samara, 443086, Russian Federation.

The online media is registered by the Federal Service for Supervision of Communications, Information Technology and Mass Communications, registration number EL No. FS 77-86495 dated December 29, 2023

Extract from the register of registered media

Regulation of the online media

Editor-in-chief: Andrey B. Prokof'yev, Doctor of Science (Engineering), associate professor,
head of the Department of Aircraft Engine Theory

2 issues a year

0+. Free price. 

Editorial address: building 22a, room 513, Soviet of Young Scientists and Specialists, 1, Academician Pavlov Street, Samara, 443011, Russian Federation.

Address for correspondence: room 513, building 22a, 34, Moskovskoye shosse, Samara, 443086, Russian Federation.

Tel.: (846) 334-54-43

e-mail: smuissu@ssau.ru

Domain name: VMUIS.RU (Domain ownership certificate), Internet email address: https://vmuis.ru/smus.

The previous certificate is a printed media, the journal “Bulletin of Young Scientists and Specialists of Samara University”, registered by the Office of the Federal Service for Supervision of Communications, Information Technologies and Mass Communications in the Samara Region, registration number series PI No. TU63-00921 dated December 27, 2017.

© Samara University

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies