Приборы и техника эксперимента, 2021, № 3, стр. 82-87

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

А. Р. Гасанов a*, Р. А. Гасанов a, Р. А. Ахмедов a, М. В. Садыхов a

a Азербайджанская национальная академия авиации
АZ1045 Баку, просп. Мардакан, 30, Азербайджан

* E-mail: afig.gasanov.51@mail.ru

Поступила в редакцию 06.10.2020
После доработки 30.10.2020
Принята к публикации 02.11.2020

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

Аннотация

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

ВВЕДЕНИЕ

Акустооптические линии задержки (а.о.л.з.) применяются для решения широкого круга радиотехнических задач, прежде всего, для обработки сигналов во временной области [16]. Различают а.о.л.з. с прямым детектированием и гетеродинного типа [7]. При этом результаты исследований а.о.л.з. с прямым детектированием в большинстве случаев могут быть экстраполированы на а.о.л.з. гетеродинного типа.

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

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

ОБЩИЕ СВЕДЕНИЯ

В а.о.л.з. с прямым детектированием (рис. 1) обрабатываемый сигнал uin(t) модулирует в амплитудном модуляторе (АМ) сигнал генератора высокой частоты (ГВЧ) по амплитуде (в зависимости от характера решаемой задачи могут быть применены амплитудная модуляция, балансная амплитудная модуляция или амплитудная манипуляция). Модулированное высокочастотное колебание с несущей частотой Ω передается на прикрепленный к торцу фотоупругой среды (ФУС) электроакустический преобразователь (ЭАП) длиной L и шириной H. В зависимости от рабочей частоты толщина ЭАП может составлять от единиц микрометров до единиц миллиметров.

Рис. 1.

Блок-схема а.о.л.з. с прямым детектированием. АМ – амплитудный модулятор; ГВЧ – генератор высокой частоты; ФУС – фотоупругая среда; ЭАП – электроакустический преобразователь; ФД – фотодетектор; Ф – фильтр.

Ячейку, состоящую из ФУС и прикрепленного к ее торцу ЭАП, называют акустооптическим модулятором (АОМ). АМ и ГВЧ обеспечивают перенос низкочастотного спектра обрабатываемого сигнала в область рабочих частот АОМ. Очевидно, что при обработке радиочастотного сигнала, который имеет несущую частоту, равную центральной частоте АОМ, необходимость в применении АМ и ГВЧ отпадает.

ЭАП возбуждает в ФУС акустическую волну длиной Λ, мощностью Pa(t), которая распространяется со скоростью $v$, в ~105 раз меньшей скорости распространения электромагнитной волны. Временная зависимость акустической мощности Pa(t) обусловлена входным сигналом uin(t).

Рисунок 1 поясняет режим дифракции Брэгга, в котором лазерное излучение интенсивностью I0, с длиной волны λ и частотой ν = ω/2π попадает в апертуру ФУС под углом θB = arcsin(0.5λ/Λ). В результате взаимодействия акустической и оптической волн в ФУС часть светового пучка с интенсивностью I1 отклоняется в первый дифракционный порядок и через отверстие в диафрагме падает на светочувствительную поверхность фотодетектора (ФД). Выделенный с помощью фильтра (Ф) выходной сигнал uout(t) является копией входного сигнала uin(t) и отстает от него на время

(1)
$\tau = {{x}_{0}}{\text{/}}{v},$
где x0 – расстояние от ЭАП до области акустооптического взаимодействия.

ДИАПАЗОН ПЛАВНО УПРАВЛЯЕМОЙ ЗАДЕРЖКИ СИГНАЛОВ

Существенная задержка τ (десятки микросекунд) сигнала обусловлена относительно низкой скоростью распространения акустической волны в ФУС, качество которой оценивают коэффициентом добротности:

$M = {{n}^{6}}{{p}^{2}}{\text{/}}(\rho {{{v}}^{3}}),$
где n, ρ и p – соответственно коэффициент преломления, плотность и фотоупругая постоянная среды взаимодействия [8].

Согласно вышеприведенному описанию работы а.о.л.з., диапазон плавно управляемой задержки сигналов полностью определяется типом ФУС и ее длиной l. Во всех случаях должно выполняться условие

(2)
$l > {{x}_{0}} = {v}\tau .$

Из анализа формул (1) и (2) легко заключить, что для расширения диапазона плавно управляемой задержки сигналов необходимо использовать фотоупругую ячейку большой длины l и с малой скоростью распространения v упругой волны в материале ячейки. Ячейку большой длины (10–20 см) можно изготовить из стеклообразных фотоупругих материалов. С точки зрения малой скорости распространения упругой волны парателлурит TeO2 считается наилучшим материалом для ФУС. Скорость медленной сдвиговой волны в TeO2 равна ≈616 м/с.

ЧАСТОТА СРЕЗА

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

(3)
${{f}_{s}} = \ln (9){\text{/}}(2\pi {{\tau }_{{\text{н}}}}) = 0.35{\text{/}}{{\tau }_{{\text{н}}}},$
где τн – время нарастания переходной характеристики а.о.л.з.

Как правило, в а.о.л.з. используется слабое акустооптическое взаимодействие. Поэтому можно считать, что интенсивности падающего и отклоненного света связаны следующим примерным равенством [9]:

(4)
${{I}_{1}}(t) \approx {{I}_{0}}\eta (t),$
где η(t) – дифракционная эффективность, которая определяется по формуле

(5)
$\eta (t) = [{{\pi }^{2}}M{{P}_{a}}(t)L]{\text{/}}(2{{\lambda }^{2}}H).$

На основе формул (4) и (5) можно написать следующее выражение для мощности отклоненного в первый порядок светового пучка с равномерным распределением интенсивности в пределах площади его поперечного сечения S1:

(6)
${{P}_{1}}(t) = {{S}_{1}}{{I}_{1}}(t) = {{S}_{1}}{{I}_{0}}\eta (t) = {{S}_{1}}\eta (t){{P}_{0}}{\text{/}}{{S}_{0}},$
где S0 и P0 = S0I0 – соответственно площадь поперечного сечения и мощность падающего в апертуру АОМ светового пучка.

При амплитудной модуляции колебаний ГВЧ напряженность электрического поля падающего на фоточувствительную поверхность ФД оптического сигнала может быть представлена в виде

(7)
$e(t) = E[1 + ms(t)]\cos (\omega --\Omega )t,$
где m – индекс амплитудной модуляции; s(t) = = uin(t)/uinmax – нормированный входной сигнал – модулирующий процесс, который может меняться в пределах $ \pm \,1$; uinmax – максимальное значение входного сигнала.

Отметим, что, согласно формуле (7), частота отклоненного света меньше на величину Ω по сравнению с частотой падающего света ω, что обусловлено эффектом Доплера.

Для определенности примем, что в качестве ФД используется фотоэлектронный умножитель. Ввиду того что ФД реагирует на интенсивность света, пропорциональную |e(t)|2, ток фотокатода iф(t) будет определяться выражением:

(8)
$\begin{gathered} {{i}_{{\text{ф}}}}(t)\sim {{E}^{2}}{{[1{\text{ }} + ms(t)]}^{2}}{{\cos }^{2}}(\omega --\Omega )t = \\ = 0.5{{E}^{2}}[1{\text{ }} + {\text{ }}2ms(t) + {{m}^{2}}{{s}^{2}}(t)]\left\{ {1 + \cos [2(\omega --\Omega )t]} \right\}. \\ \end{gathered} $

Коэффициент пропорциональности в (8) определяется как среднее значение тока фотокатода:

(9)
${{\bar {i}}_{{\text{ф}}}} = {{P}_{1}}e\eta {\kern 1pt} {\text{'/}}(h\nu ),$
где P1/(hν) – число фотонов, которые падают на фоточувствительную поверхность ФД за 1 с; h = = 6.63 ⋅ 10–34 Дж · с – постоянная Планка; е = 1.6 × × 10–19 Кл – заряд протона; η' – квантовый выход фотокатода (среднее число электронов, испускаемых фотокатодом при падении на него одного фотона).

Подставляя (9) в (8) получаем:

(10)
$\begin{gathered} {{i}_{{\text{ф}}}}(t) = ({{P}_{1}}{\text{/}}(h\nu ))e\eta {\text{'}} \times \\ \times \,\,[1{\text{ }} + {\text{ }}2ms(t){\text{ }} + {{m}^{2}}{{s}^{2}}(t)]\{ 1 + \cos [2(\omega --\Omega )t]\} . \\ \end{gathered} $

Модулирующий процесс из (10) выделяется фильтром Ф, на выходе которого формируется ток

(11)
${{i}_{{{\text{out}}}}}(t) = 2Gme\eta {\text{'}}{{P}_{1}}s(t){\text{/}}(h\nu ),~$
где G – коэффициент усиления фотоэлектронного умножителя.

Следует выделить две области формирования тока на выходе Ф: область стационарного режима и область установления стационарного режима. В стационарном режиме мощность P0 падающего в апертуру АОМ светового пучка и отношение поперечных сечений падающего и отклоненного световых пучков S1/S0 постоянны, т.е. P0 = const, S1/S0 = const. Поэтому временная зависимость iout(t) обусловлена только входным сигналом uin(t), т.е. модулирующим процессом s(t). При отсутствии модуляции, т.е. при s(t) = 0, определяемый соотношением (11) выходной ток Ф будет равен нулю: iout(t) = 0.

В контексте решаемой задачи процесс установления стационарного режима рассматривается при отсутствии модулирующего процесса, т.е. при s(t) = 0. При этом вид функциональной зависимости S1(t) определяется временем и характером вхождения упругой волны в световой пучок. В этих условиях процесс установления стационарного режима нормируется выражением

(12)
$g(t) = {{S}_{1}}(t){\text{/}}{{S}_{0}}.$

Отметим, что в данном случае S0 – площадь поперечного сечения отклоненного светового пучка в стационарном режиме, т.е. наибольшее значение S1(t).

С учетом того, что АОМ является широкополосным устройством (полоса пропускания составляет 40–60% от центральной частоты), функциональное выражение (12) может рассматриваться как переходная характеристика а.о.л.з. c прямым детектированием.

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

(13)
${{S}_{1}}(x)\, = \,\int\limits_x^{x + {v}{{\tau }_{i}}} {H{\kern 1pt} '(x)dx} \,\,{\text{при}}\,\,{{x}_{0}} \leqslant x \leqslant {{x}_{0}} + d + v{{\tau }_{i}},$
где
(14)
$\begin{gathered} H{\text{'}}(x) = 2{{[d(x--{{x}_{0}})--{{(x--{{x}_{0}})}^{2}}]}^{{0.5}}} \\ {\text{при}}\quad {{x}_{0}} \leqslant x \leqslant {{x}_{0}} + d. \\ \end{gathered} $
Здесь H'(x) – длина линии взаимодействия переднего фронта упругого волнового пакета с лазерным пучком в плоскости x0z (см. рис. 2), x – текущая координата, x0 – расстояние от ЭАП до точки акустооптического взаимодействия.

Рис. 2.

Проекция акустооптического взаимодействия в плоскости, перпендикулярной направлению распространения лазерного пучка.

В области стационарного режима, т.е. при выполнении условия x > x0 + d, площадь акустооптического взаимодействия будет иметь максимальное значение и определяться как:

(15)
${{S}_{1}}(x){{{\text{|}}}_{{x = {{x}_{0}} + d}}} = {{S}_{0}} = \int\limits_{{{x}_{0}}}^{{{x}_{0}} + d} {H{\kern 1pt} {\text{'}}(x)\,dx} = 0.25\pi {{d}^{2}}.$

В выбранной системе отсчета текущая координата x связана с текущим временем t равенством $x = vt$. Поэтому соотношения (13) и (14) переносятся на временную плоскость в следующей форме:

${{S}_{1}}(t) = v\int\limits_\tau ^t {H{\kern 1pt} {\text{'}}(t)\,dt} \quad {\text{при}}\quad \tau \leqslant t \leqslant \tau + d{\text{/}}v;$
$\begin{gathered} H{\kern 1pt} {\text{'}}(t) = 2v\sqrt {\frac{d}{v}t - {{t}^{2}}} \left[ {\sigma (t - \tau ) - \sigma \left( {t - \tau - \frac{d}{v}} \right)} \right] \\ {\text{при}}\quad \tau \leqslant t \leqslant \tau + d{\text{/}}v. \\ \end{gathered} $

На основе последних выражений конечное уравнение для площади поперечного сечения отклоненного светового пучка приобретает следующий вид:

(16)
$\begin{gathered} {{S}_{1}}(t) = 2{{v}^{2}}\int\limits_\tau ^t {\sqrt {\frac{d}{v}(t - \tau ) - {{{(t - \tau )}}^{2}}} \,dt} \\ {\text{при}}\quad \tau \leqslant t \leqslant \tau + d{\text{/}}v. \\ \end{gathered} $

Подставляя (15) и (16) в (12) получаем следующее выражение для переходной характеристики а.о.л.з. с прямым детектированием:

(17)
$g(t) = \frac{{4v}}{{\pi {{d}^{2}}}}\int\limits_\tau ^t {H{\kern 1pt} {\text{'}}(\xi )d\xi } \quad {\text{при}}\quad \tau \leqslant t \leqslant \tau + {{\tau }_{0}}$
или
(18)
$\begin{gathered} g(t) = \frac{8}{{\pi \tau _{0}^{2}}}\int\limits_\tau ^t {\sqrt {{{\tau }_{0}}(\xi - \tau ) - {{{(\xi - \tau )}}^{2}}} d\xi } \\ {\text{при}}\quad \tau \leqslant t \leqslant \tau + {{\tau }_{0}}, \\ \end{gathered} $
где τ0 = d/v время установления стационарного режима.

Время нарастания τн переходной характеристики определяется как время, в течение которого функция g(t) изменяется от 0.1 до 0.9 от своего максимального значения. В соответствии с этим определением на основе выражения (18) составим следующую систему уравнений для вычисления τн:

(19)
$\left\{ \begin{gathered} \frac{8}{{\pi \tau _{0}^{2}}}\int\limits_\tau ^{{{t}_{1}}} {\sqrt {{{\tau }_{0}}(\xi - \tau ) - {{{(\xi - \tau )}}^{2}}} d\xi } = 0.1; \hfill \\ \frac{8}{{\pi \tau _{0}^{2}}}\int\limits_\tau ^{{{t}_{2}}} {\sqrt {{{\tau }_{0}}(\xi - \tau ) - {{{(\xi - \tau )}}^{2}}} d\xi } = 0.9; \hfill \\ {{\tau }_{{\text{н}}}} = {{t}_{2}} - {{t}_{1}}. \hfill \\ \end{gathered} \right.$

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ

Вычисленный по формуле (18) расчетный график переходной характеристики а.о.л.з. с параметрами $v$ = 3600 м/с, d = 2.2 мм, τ = 0.5 мкс изображен на рис. 3.

Рис. 3.

Расчетный график переходной характеристики а.о.л.з. с параметрами $v$ = 3600 м/с, d = 2.2 мм, τ = = 0.5 мкс.

Время нарастания τн переходной характеристики вычисляется по формулам (19), которые легко моделируются в среде Mathcad. Расчетное значение составило τн = 0.419 мкс и однозначно совпадает со значением, определяемым по графику на рис. 3.

Вычисленная по формуле (3) частота среза амплитудно-частотной характеристики а.о.л.з. составила  fs = 835 кГц.

Численный анализ зависимостей времени нарастания от скорости распространения упругой волны в ФУС при различных диаметрах светового пучка и от диаметра светового пучка при различных скоростях распространения упругой волны в ФУС проведен по уравнениям (19). Результаты численного анализа представлены в виде графиков на рис. 4 и на рис. 5 соответственно.

Рис. 4.

Зависимости времени нарастания переходной характеристики от скорости распространения упругой волны в ФУС при двух значениях диаметра светового пучка.

Рис. 5.

Зависимости времени нарастания переходной характеристики от диаметра светового пучка при двух значениях скорости распространения упругой волны в ФУС.

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

ЭКСПЕРИМЕНТАЛЬНАЯ АПРОБАЦИЯ

Вышеприведенные интерпретации физических процессов в а.о.л.з. с прямым детектированием, установленные на их основе положения и результаты численного анализа апробированы экспериментально на макете установки, схема которой приведена на рис. 6. Прямоугольный импульс с заданными параметрами от генератора Г5-54 модулирует колебание высокочастотного генератора Г4-107 (работает в режиме внешней импульсной модуляции) и синхронизирует осциллограф MSO4052. Частота колебания генератора Г4-107 выбирается равной центральной частоте АОМ, что в наших экспериментах составляет 80 МГц.

Рис. 6.

Схема экспериментальной установки. АОМ – акустооптический модулятор.

Осциллограммы напряжений на входе и выходе а.о.л.з. с прямым детектированием с параметрами $v$ = 3600 м/с, d = 2.2 мм, τ = 0.7 мкс приведены на рис. 7. Длительность входного импульса (определяется по осциллограмме на уровне 0.5 от максимального значения) равна τi ≈ 1.07 мкс. Время нарастания, т.е. время, в течение которого напряжение на выходе изменяется от 0.1 до 0.9 от своего максимального значения, равно примерно 419 нс, что полностью совпадает с временем нарастания переходной характеристики на рис. 3. Осциллограмма 2 на рис. 7б в целом хорошо коррелирует с данными расчетных графиков на рис. 4 и 5. Небольшие отклонения обусловлены в основном неидеальностью входного воздействия и инерционностью нагрузочной цепи.

Рис. 7.

Осциллограммы импульсов на входе (1) и на выходе (2) а.о.л.з. с прямым детектированием с параметрами $v$ = 3600 м/с, d = 2.2 мм, τ = 0.7 мкс. Масштаб: по вертикали – 2 В/деление (1), 500 мВ/деление (2); по горизонтали – 200 нс/деление.

ВЫВОДЫ

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

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

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

  1. Shakin O.V., Nefedov V.G., Churkin P.A. // Prog. of Conf. “Wave Electronics and its Application in Information and Telecommunication Systems”. (St. Petersburg, Russia) St. Petersburg: State University of Aerospace Instrumentation, 2018. NOV 26–30. P. 340.

  2. Yushkov K.B., Molchanov V.Ya., Ovchinnikov A.V., Chefonov O.V. // Physical Review. 2017. V. 96. Iss. 4. № 043866. https://doi.org/10.1103/PhysRevA.96.043866

  3. Schubert O., Eisele M., Crozatier V., Forget N., Kaplan D., Huber R. // Optics Letters. 2013. V. 38. P. 2907. https://doi.org/10.1364/OL.38.002907

  4. Chandezon J., Rampnoux J.-M., Dilhaire S., Audoin B., Guillet Y. // Optics Express. 2015. V. 23. P. 27011-27019. https://doi.org/10.1364/OE.23.027011

  5. Okon-Fafara M., Kawalec A.M., Witczak A. // Proc. of SPIE. 2019. V. 11055. https://doi.org/10.1117/12.2525032

  6. Diewald A.R., Steins M., Müller S. // Advances in Radio Science. 2018. V. 16. P. 203. https://doi.org/10.5194/ars-16-203-2018

  7. Гасанов А.Р., Гасанов Р.А. // ПТЭ. 2017. № 5. С. 112. https://doi.org/10.7868/S0032816217050081

  8. Christofer C.D. Lasers and Electro-optics. Cambridge University Press, 2014. https://doi.org/10.1017/CBO9781139016629

  9. Lee J.N., van der Lugt A. // Proc. IEEE. 1989. V. 77. № 10. P. 158.

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