Вестник Военного инновационного технополиса «ЭРА», 2023, T. 4, № 4, стр. 393-397

Управление показателем преломления линзовых элементов путем изменения температуры адаптивной оптики

Е. В. Медведев 1*, Д. С. Бурый 1, С. А. Гулин 1, В. Г. Страдинь 1, Р. В. Чумарин 1

1 Военный инновационный технополис “ЭРА”
Анапа, Россия

* E-mail: era_lab13@mil.ru

Поступила в редакцию 09.10.2023
После доработки 11.01.2024
Принята к публикации 20.01.2024

Полный текст (PDF)

Аннотация

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

ВВЕДЕНИЕ

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

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

В данной работе описана математическая модель коррекции волнового фронта ОС с помощью использования системы термонагревателей, которые размещаются на оптических (зеркальных или линзовых) элементах и управляются согласно заданному алгоритму. Следует отметить высокую точность тепловых методов коррекции. Управление перепадами температуры в диапазоне 0.03‒0.1°С позволяет создавать перемещения поверхности зеркальных элементов и, соответственно, фазовые коррекции ВФ с погрешностью до λ/40–λ/20 (λ – длина световой волны). Подразумевается, что при изначальных расчетах оптической системы будет сделана аберрационная коррекция предполагаемых величин искажения ВФ, а затем при помощи элемента с термоуправлением, который вводится в схему специально с целью высокоточного управления ВФ всей оптической системы, будет проводиться терморегулировка. Преимуществом теплового метода является возможность его применения как для управления формой зеркальных поверхностей, так и для изменения показателя преломления линзовых элементов. К недостаткам стоит отнести инерционность и долговременность, что ограничивает сферы применения адаптивной оптики на основе термоуправления [3].

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ

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

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

Представим ВФ оптического излучения W(x) в некоторый момент времени t в виде разложения по системе ортогональных функций Fm(x):

(1)
$W(x,t) = \sum\limits_{m = 1}^\infty {{{A}_{m}}(t){{F}_{{m\left( x \right)}}}} .$

Ортогональность функций Fm(x) означает, что

${{F}_{m}}{\text{|}}{{F}_{n}} = {{\left| {{{F}_{m}}} \right|}^{2}}{{{{\sigma }}}_{{mn}}},\quad {{{{\sigma }}}_{{mn}}} = \left\{ {\begin{array}{*{20}{c}} {1,~m = n;} \\ {0,~m \ne n.} \end{array}} \right.$

Введем обозначения для скалярного произведения функций:

(2)
${{F}_{m}}{\text{|}}{{F}_{n}} \geqslant \int\limits_D {F_{m}^{*}(x){{F}_{n}}(x)dD~} ,$
где $*$ – знак комплексного сопряжения, D – фазовое пространство системы (например, площадь апертуры S в рассматриваемом случае), x ϵ D – точка в пространстве D.

Пусть амплитуды Аm – моды искажений ВФ. Задача работы корректора ВФ состоит в обеспечении заданных значений аберраций $A_{m}^{0}$. Требуемый волновой фронт W0(x) записывается в виде

${{W}_{0}}(x) = \sum\limits_{m = 1}^\infty {A_{m}^{0}{{F}_{m}}(x)} .$

Для ошибки ΔW (отклонения ВФ от заданной формы) имеем

(3)
$\Delta W = W - {{W}_{0}} = \sum\limits_{m = 1}^\infty {[{{A}_{m}}(t) - A_{m}^{0}]{{F}_{m}}(x)} .$

Задача управления – обратить ΔW в значение, заданное ранее для коррекции величины искажения ВФ. Выполнить эту задачу в полном объеме не представляется возможным, так как число степеней свободы корректора ограничено. Пусть имеется возможность компенсации первых N аберраций [4], т.е. активный элемент (зеркало или линза) может вносить дополнительный фазовый сдвиг U в пределах апертуры:

(4)
$U(x,t) = \sum\limits_{m = 1}^N {{{B}_{m}}(t){{F}_{m}}(x)} .$

Тогда для ошибки $\Delta W$ можно записать:

(5)
$\begin{gathered} \Delta W = W - {{W}_{0}} - U = \\ \, = \sum\limits_{m = 1}^N {[{{A}_{m}}(t) - A_{m}^{0} - {{B}_{m}}(t)]{{F}_{m}}(x) + {{\varepsilon }}} , \\ \end{gathered} $
где ${{\varepsilon }}$ – остаточная неустранимая ошибка.

Цель работы корректора – обеспечение по каждой моде выработки управляющих воздействий ${{B}_{m}}\left( t \right)$ таких, чтобы $\Delta W$ принимало необходимые значения.

Необходимая величина сдвига ${{B}_{m}}$ должна обеспечиваться работой нагревателей, расположенных на активном элементе. Если, например, W представляет собой форму поверхности плоского круглого зеркала, то дополнительный фазовый сдвиг U будет представлять собой перемещение (деформацию) отражающей поверхности в поперечном к плоскости зеркала направлении.

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

МАТРИЦА СВЯЗИ МЕЖДУ ТЕМПЕРАТУРНЫМИ МОДАМИ И ОПТИЧЕСКИМИ АБЕРРАЦИЯМИ

Температурное поле T(x,t) в теле, имеющем объем V и ограниченном поверхностью L, при наличии теплового воздействия Q определяется уравнением теплопроводности:

(6)
$c{{\rho }}\frac{{dT}}{{dt}} = {{\chi \Delta }}T + Q,$
с начальным
(7)
$T(x,0) = f(x)$
и граничным условиями

(8)
${{\chi }}\frac{{dT(x,t)}}{{dn}} + {{\alpha [}}T(x,t) - {{\Theta (}}x,t{\text{)]}} = 0,\quad x \in L.$

Здесь введены $c$, $\rho $, ${{\chi }}$ – соответственно удельная теплоемкость, плотность и коэффициент теплопроводности тела, ${{\alpha }}$ – коэффициент теплообмена,$~{{\Theta }}$ – температура окружающей среды, n – внешняя нормаль к границе тела L.

Функцию Т можно разложить по системе ортонормированных (ортогональных и нормированных на единицу) функций [5]:

(9)
$T(x,t) = \sum\limits_{k = 1}^\infty {{{C}_{k}}(t){{{{\Phi }}}_{k}}(x)} ,~\quad {{{{\Phi }}}_{m}}{\text{|}}{{{{\Phi }}}_{n}} = {{{{\delta }}}_{{mn}}}.$

Реализуем синтез управления температурным полем из условия, что функции ${{\Phi }}$ удовлетворяют граничным условиям (8) и амплитуды Ck изменяются со временем независимо друг от друга. Если подставить выражение (9) в уравнение (6), то после умножения уравнения на ${{{{\Phi }}}_{m}}*$ и интегрирования получим

(10)
$\begin{gathered} c{{\rho }}\frac{{d{{C}_{m}}}}{{dt}} = {{\chi }}\sum\limits_{k = 1}^\infty {\left\langle {{{{{\Phi }}}_{m}}{\text{|}}\Delta {{{{\Phi }}}_{k}}} \right\rangle \Delta {{C}_{k}} + {{{{\Phi }}}_{m}}{\text{|}}Q} ,~ \\ m = 1,~2,~ \ldots ~~. \\ \end{gathered} $

В общем случае при недиагональности матрицы $\left\langle {{{{{\Phi }}}_{m}}{\text{|}}\Delta {{{{\Phi }}}_{k}}} \right\rangle $, Сk подчиняются бесконечной системе связанных друг с другом обыкновенных дифференциальных уравнений.

Потребуем, чтобы функции ${{\Phi }}$ удовлетворяли условиям

(11)
$\left\langle {{{{{\Phi }}}_{m}}{\text{|}}\Delta {{{{\Phi }}}_{k}}} \right\rangle = - {{{{\Lambda }}}_{m}}{{{{\delta }}}_{{mk}}},$
т.е. являлись бы собственными функциями оператора Лапласа с собственным значением $ - {{{{\Lambda }}}_{m}}$. Тогда из (10) получим уравнение для амплитуд Сm:

(12)
${{\dot {C}}_{m}} = - \frac{1}{{{{{{\tau }}}_{m}}}}{{C}_{m}} + {{P}_{m}},\quad m = 1,{\text{\;}}2{\text{\;}} \ldots \,{\kern 1pt} {\text{\;}}.$

Здесь введено обозначение для собственных времен тепловой задачи

(13)
${{{{\tau }}}_{m}} = \frac{{c{{\rho }}}}{{{{\chi }}{{{{\Lambda }}}_{m}}}}~$
и введены компоненты вектора модальных управляющих тепловых воздействий

(14)
${{P}_{m}} = \left\langle {{{{{\Phi }}}_{m}}{\text{|}}Q} \right\rangle = \int\limits_V {{{\Phi }}_{m}^{*}(x)Q(x,t)dV} .$

Размерность базиса для разложения температурного поля (число тепловых степеней свободы) NT должно быть не меньше N. Примем для простоты изложения, что NT = N. Тогда для T(x,t) имеем

(15)
$T(x,t) = \sum\limits_{k = 1}^N {{{C}_{k}}(t){{{{\Phi }}}_{k}}(x)} ,\quad \left\langle {{{{{\Phi }}}_{m}}{\text{|}}{{{{\Phi }}}_{n}}} \right\rangle = {{{{\delta }}}_{{mn}}},$
где функции ${{{{\Phi }}}_{k}}~$ удовлетворяют условиям (8), (9), (11). Амплитуды Сk, описываемые уравнением (12), назовем температурными модами.

Из ортонормированности функций ${{\Phi }}$ следует, что температурные моды представляют собой ортогональные проекции поля температуры на базисные векторы. Найдем в общем виде связь между оптическими модами Bm и температурными модами Ck. Прежде всего отметим, что между дополнительным фазовым сдвигом U, обусловленным работой корректора, и полем температур T(x,t) имеет место связь

(16)
$U(x,t) = \int\limits_V {\hat {G}(x,y)T(y,t)dy~} ,$
где линейный оператор ${\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{G} }$ – функция Грина для решения уравнений термоупругости.

Подставим в (16) разложение (4) для U и (15) для Т:

(17)
$\sum\limits_{m = 1}^N {{{B}_{m}}(t){{F}_{m}}(x)} = \sum\limits_{k = 1}^N {{{C}_{k}}(t)} \int\limits_V {\hat {G}(x,y){{{{\Phi }}}_{k}}(y,t)dy} ~{\kern 1pt} .$

Умножим обе части равенства (16) на $F_{m}^{*}$ (x) и проинтегрируем по х. В силу ортогональности функций F получим

$\begin{gathered} {{B}_{m}}\left( t \right){{\left| {{{F}_{m}}} \right|}^{2}} = \\ \, = \sum\limits_{k = 1}^N {{{C}_{k}}(t)} ~\int\limits_D {\int\limits_V {\hat {G}(x,y){{F}_{m}}(x){{{{\Phi }}}_{k}}(y)dxdy.~} } \\ \end{gathered} $

Если принять, что

(18)
$\left\langle {{{F}_{m}}{\text{|}}{{F}_{n}}} \right\rangle = {{{{\delta }}}_{{mn}}},$
то связь между заданным полем U и значения-ми В, обеспечивающими его, имеет следующий вид:
(19)
${{B}_{m}}(t) = \left\langle {{{F}_{m}}{\text{|}}U~} \right\rangle ,$
а связь между Сk и Bm может быть записана в виде

(20)
${{B}_{m}}(t) = \sum\limits_{k = 1}^N {{{O}_{{mk}}}{{C}_{k}}(t)} .$

Здесь введено обозначение для термооптической матрицы:

(21)
$\begin{gathered} {{O}_{{mk}}} = \left\langle {{{F}_{m}}{\text{|}}\hat {G}{\text{|}}{{{{\Phi }}}_{k}}} \right\rangle = \\ \, = \int\limits_D {\int\limits_V {\hat {G}(x,y){{F}_{m}}(x){{{{\Phi }}}_{k}}(y)dxdy} } ~. \\ \end{gathered} $
где интегрирование идет по объему тела V и апертуре D.

Таким образом, управляя температурными модами Ск, можно обеспечить требуемые значения оптических мод (фазовых сдвигов) Вm. И наоборот, необходимые значения температурных мод для получения требуемых величин Вm получаются с помощью обратной матрицы R = О1:

(22)
${{C}_{k}}(t) = \sum\limits_{m = 1}^N {{{R}_{{km}}}{{B}_{m}}(t)} .$

Для осуществления теплового управления разместим на активном элементе систему из Nq нагревателей, имеющих малые размеры по сравнению с размерами элемента. В математической модели объекта регулирования, каковым является активный элемент, нагреватели можно считать точечными, т.е. записать тепловое поле Q(x), создаваемое ими в виде

(23)
$Q(x) = \sum\limits_{a = 1}^{{{N}_{q}}} {{{q}_{a}}\delta (x - {{x}_{a}})} ,$
где δ(хха) – дельта-функция, ха – точка, в которой расположен нагреватель номер “a”, qa – поданная на него мощность.

Заметим, что величина Q, входящая в уравнение теплопроводности (6), имеет размерность плотности мощности (Вт/м2); с другой стороны, обобщенная функция δ(xxa), удовлетворяющая условию

(24)
$\int\limits_V {f(x){{\delta (}}x - {{x}_{0}}{\text{)}}dx = f({{x}_{0}})} ,$
также является размерной (м–3). Поэтому qa измеряется в ваттах и является непосредственно величиной мощности, выделяемой нагревателем с номером “а”.

Задача работы нагревателей – создать стационарное температурное поле, обеспечивающее требуемую форму ВФ [6]. При снятии теплового воздействия температурное поле изменится. Таким образом, речь идет о создании устойчивого положения в неравновесной системе за счет подвода энергии извне.

Подставляя выражение (23) в (14) с учетом (24) для компоненты вектора управляющих воздействий Рm, находим

${{P}_{m}} = \sum\limits_{a = 1}^{{{N}_{0}}} {{{{{\Phi }}}_{m}}({{x}_{a}}){{q}_{a}}} .~$

Если ввести матрицу Mma:

(25)
$\begin{gathered} {{M}_{{ma}}} = {{{{\Phi }}}_{m}}({{x}_{a}}),\quad m = 1,{\text{\;}}2,{\text{\;}} \ldots ,N,~ \\ a = 1,{\text{\;}}2,{\text{\;}} \ldots ,\quad {{N}_{q}}~,~ \\ \end{gathered} $
то для Pm получим

(26)
${{P}_{m}} = \sum\limits_{a = 1}^{{{N}_{q}}} {{{M}_{{ma}}}{{q}_{a}}} ,\quad ~m = 1,{\text{\;}}2,{\text{\;}} \ldots ,\quad N.$

Подставляя выражения для Pm в уравнение (12), а затем их решения в формулы (20), можно определить функцию отклика оптической системы на произвольное тепловое воздействие {qa}.

При синтезе управления требуется решить обратную задачу: на основании уравнения движения (13) определяются компоненты вектора Pm, обеспечивающие заданные значения амплитуд Cm, и на основании (26) рассчитываются необходимые значения мощностей на каждый нагреватель [7]:

(27)
${{q}_{a}} = \sum\limits_{m = 1}^N {{{S}_{{am}}}{{P}_{m}}} ,\quad a = 1,{\text{\;}}2,{\text{\;}} \ldots ,\;{{N}_{q}},~\quad S = {{M}^{{ - 1}}}~.$

Однако здесь возникает одна принципиальная проблема. Матрица S, обратная матрице М, определяется однозначно, только если Nq = N (квадратная матрица). Но в этом случае при произвольных величинах Рm получим для q как положительные, так и отрицательные значения. Отрицательные значения мощности соответствуют охлаждению, однако такая задача решается иным путем: берется число нагревателей больше числа управляемых мод и задается условием минимизации энергопотребления. Тогда задача решается однозначно, но возникают трудности, связанные с возрастанием температуры объекта, однако этот рост не будет значительным.

Остается найти алгоритм определения компонент вектора Pm на основании уравнений (12).

Обозначим значения Сk, необходимые для обеспечения требуемых величин амплитуд Bm, рассчитанных на основании (24), как $C_{k}^{0}$ и введем отклонения Ck(t) от требуемого значения:

(28)
${{\tilde {C}}_{k}}(t) = {{C}_{k}}(t) - C_{k}^{0}~.$

Величины ${{\tilde {C}}_{k}}\left( t \right)$ удовлетворяют уравнениям

(29)
$\frac{{d{{{\tilde {C}}}_{k}}}}{{dt}} = - \frac{1}{{{{{{\tau }}}_{k}}}}{{\tilde {C}}_{k}} + {{P}_{k}} - \frac{1}{{{{{{\tau }}}_{k}}}}C_{k}^{0}.~$

Задача воздействий Pk состоит в том, чтобы из любого начального состояния привести величины Ck к заданным значениям и поддерживать их в таком положении. Связь между Pk и Ck может быть записана в линейной форме посредством обратной передаточной функции $\hat {I}$:

(30)
${{P}_{k}}(t) = \int_0^t {{{I}_{k}}(t,t{\kern 1pt} '){{C}_{k}}(t{\kern 1pt} ')dt{\kern 1pt} '} ,$
где в качестве момента t = 0 выбрано начало регулирования.

Или в операторной форме

(31)
${{P}_{k}}(t) = \hat {I}_{{tt{\kern 1pt} '}}^{{(k)}}{{C}_{k}}(t{\kern 1pt} '),~$
где действие линейного оператора $\hat {I}$ определяется в (28). Полученные величины подставим в (27) и получим

(32)
$\begin{gathered} {{q}_{a}}(t) = \mathop \sum \limits_{k = 1}^n \sum\limits_{m = 1}^N {\{ {{S}_{{ak}}}\hat {I}_{{tt{\kern 1pt} '}}^{{(k)}}{{R}_{{km}}}\} ~{{B}_{m}}(t)} ~, \\ a = 1,{\text{\;}}2,{\text{\;}} \ldots ,\quad {{N}_{q}}. \\ \end{gathered} $

При помощи выражения (32) решается задача определения мощности, выводимой на каждый из нагревателей для обеспечения корректором заданного сдвига ВФ U. Время достижения tдост требуемых значений амплитуд Bm из произвольного начального состояния не превышает минимального из собственных тепловых времен объекта, учитываемых в управлении. Время τN определяется материалом и размерами, из которого изготовлен термоактивный элемент; для обычно используемых зеркал τN составляет от 1 до 30 мин. После достижения заданного значения амплитуд требуемый сдвиг ВФ может удерживаться сколь угодно долго.

ЗАКЛЮЧЕНИЕ

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

Список литературы

  1. Лукин В.П. // Интерэкспо ГЕО-Сибирь. 2019. № 8. С. 115.

  2. Шиховцев А.Ю., Лукин В.П., Ковадло П.Г. // Оптика атмосферы и океана. 2021. № 5 (388). С. 385.

  3. Русских И.В., Колобов Д.Ю., Киселев А.В. // Труды Международной Байкальской молодежной научной школы и XVI Конференции молодых ученых. 2019. С. 367.

  4. Гоголева Е.М., Фарафонтова Е.П. Прикладная оптика. Учебное пособие. Саратов, Екатеринбург: Профобразование. УрФУ, 2019. 187 с.

  5. Рукосуев А.Л., Белоусов В.Н., Мар Г.Н. и др. // Динамические процессы в геосферах. 2020. № 12. С. 156.

  6. Архипов С.А. // Контенант. 2019. Т. 18. № 4. С. 1.

  7. Заварзин В.И., Батшев В.И., Польщикова О.В. Компьютерные технологии и моделирование в оптотехнике. Учебное пособие. М.: МГТУ им. Н.Э. Баумана, 2017. 40 с.

  8. Дьяконов В.П. MATLAB. Полный самоучитель. ДМК-Пресс, 2017. 768 с.

Дополнительные материалы отсутствуют.