ДИНАМИКА ГЕНЕРАЦИИ БЫСТРОПРОТОЧНОГО ЛАЗЕРА
С НЕОДНОРОДНЫМ ВОЗБУЖДЕНИЕМ АКТИВНОЙ СРЕДЫ
В НЕУСТОЙЧИВОМ РЕЗОНАТОРЕ
А.В.Мушенков, А.И.Одинцов, Н.Э.Саркаров, А.И.Федосеев,
А.В.Федянович
С помощью численного моделирования исследовано влияние пространственно неоднородной стационарной накачки на динамику генерации лазера с быстрым поперечным протоком активной среды через неустойчивый резонатор. Показано, что в зависимости от распределения скорости накачки внутри резонатора она может как стабилизировать режим стационарной генерации, так и приводить к неустойчивости последней. Возникновение неустойчивости связано с образованием под действием неоднородной накачки значительных поперечных градиентов стационарных распределений коэффициента усиления и интенсивности поля. Для профиля накачки, характеризующегося резким спадом скорости возбуждения вблизи оптической оси резонатора, в зависимости от параметров системы наблюдались, помимо стационарной генерации, различные режимы автомодулированной генерации, а также хаотическая генерация.
Введение
Особенности динамики генерации в неустойчивом резонаторе (НР) с поперечным протоком газовой активной среды исследовались в ряде работ [1-9]. Было показано, что при определенных условиях в такой системе имеет место автоколебательная неустойчивость и возникают регулярные автопульсации излучения. Колебания интенсивности поддерживаются за счет обратной связи между центральной и периферийными зонами НР, создаваемой потоком среды [1]. Как отмечалось в [4], неустойчивость стационарной генерации (СГ) существенным образом связана с неоднородностью поля в направлении потока, которая обычно присутствует в такой системе. В работе [9] было показано, что искусственное увеличение поперечного градиента поля в НР с неоднородными зеркалами приводит к расширению области существования автомодулированной генерации.
Наличие внутри резонатора пространственно - однородной накачки, как и любого другого процесса восстановления инверсии (например, колебательного обмена) способствует выравниванию профилей коэффициента усиления и поля и тем самым стабилизирует режим стационарной генерации [3, 4]. Однако если скорость накачки в сечении резонатора непостоянна, то можно ожидать, что при некотором ее профиле устойчивость СГ будет нарушаться.
В настоящей работе на основе модельных расчетов исследуется динамика генерации быстропроточного лазера (БПЛ) с неоднородной накачкой внутри НР. Найдено, что в зависимости от профиля накачки и параметров системы наряду со стационарной генерацией имеют место режимы автомодулированной импульсно-периодической генерации, а также хаотическая генерация. Определены области существования указанных режимов. Отмечается, что путем изменения профиля накачки в резонаторе БПЛ можно управлять временными характеристиками лазера.
Расчетная модель
Расчеты проводились в приближении геометрической оптики для одномерной модели НР цилиндрической геометрии [9]. Принималось, что активная среда целиком заполняет резонатор, потери на зеркалах и в среде отсутствуют, а потери, связанные с расширением пучка в НР равномерно распределены по длине резонатора
L (M - коэффициент увеличения на двойной проход). При не слишком большом увеличении (M Ј 2) в данной модели можно пренебречь изменением поля по длине резонатора и рассматривать одномерное распределение интенсивности I(x,t) по поперечной координате. С учетом того, что в БПЛ время затухания поля в резонаторе t r=(cq )-1 обычно намного меньше характерных времен релаксации населенностей уровней и времени пролета газа через оптический канал t f = h/v ( h - апертура резонатора, v - скорость потока), для расчета интенсивности излучения использовалось уравнение квазистационарной генерации. (1)
Здесь
a (x,t) - коэффициент усиления среды, координата x отсчитывается от оси резонатора навстречу потоку.Активная среда описывалась простейшим модельным уравнением для коэффициента усиления с одной релаксационной константой [2]
, (2)
где
s - сечение оптического перехода, t - время релаксации инверсии, S - скорость накачки.Система уравнений (1-2) решалась численно для верхней (по потоку) части НР 0
< x< h. (В приближении геометрической оптики нижняя часть НР (x<0) не влияет на генерацию в верхней его части.) На входе в резонатор задавалось постоянное значение коэффициента усиления a (h, t) = a in. Граничное значение интенсивности I(0,t) на оси НР определялось из условия равенства усиления и потерь a (0,t)=q , которое должно выполняться в генерирующем лазере в квазистационарном приближении. С учетом этого условия из (2) имеем. (3)
Соотношение (3) выражает баланс энергии в окрестности оси НР, учитывающий скорость поступления возбужденных молекул к оси, процессы накачки, релаксации и генерации излучения. В большинстве расчетов профиль накачки задавался в виде
, (4)
где
h0 - ширина зоны неоднородности накачки вблизи оси НР (h0<< h), Sm - максимальная скорость накачки вдали от оси (x>>h0). При данном выборе S(x) накачка на оси резонатора отсутствует. Величина Sm определяет максимальный коэффициент усиления a m= s Smt , который достигается в неподвижной среде в отсутствие генерации. Часть расчетов была выполнена с прямоугольным профилем накачки S(x) = 0 при xЈ h0 и S(x) = Sm при x> h0. За исключением области очень малых потерь q <<a m результаты, полученные с обоими профилями были идентичны.В расчетах использовались безразмерные переменные
t/t f, x/h и нормированные значения коэффициента усиления g = a /q , интенсивности w= s It f и скорости накачки q = s St f/q . Система полностью описывается заданием четырех независимых параметров t / t f, h0 /h, q /a m, a in/ a m, которые и варьировались в численных расчетах. Для величины входного коэффициента усиления в большинстве расчетов бралось значение a in=a m. Заметим, что при a in=a m и h0<<h параметр t / t f даёт соотношение скоростей вноса возбужденных молекул в резонатор и их образования за счет накачки внутри резонатора.В принятой модели развитие неустойчивости СГ описывается линеаризованными уравнениями для малых возмущений
d a (x, t) и d w(x, t), получаемыми из (1), (2). В безразмерных величинах уравнения имеют вид [2](5)
, (6)
где относительные возмущения, индексом s обозначены стационарные распределения. Эти уравнения использовались в анализе механизма неустойчивости СГ в резонаторе с неоднородной накачкой.
Результаты расчетов и обсуждение
В результате численного моделирования динамики генерации в описанной системе были получены нестационарные профили коэффициента усиления
a (x, t) и интенсивности поля I(x, t) при различных значениях параметров системы. Их анализ позволяет выделить 4 основные области значений параметров, каждая из которых характеризуется определенным типом поведения системы. На рис.1 эти области показаны в координатах t /t 0, где t 0 = h0/v - время пролета газа через зону неоднородности, и нормированной величины потерь q /a m. Приведены данные для трех значений отношения h0/h = 0,05; 0,1 и 0,2 при величине a in/a m=1. Сплошная кривая 1 отделяет область I существования устойчивой СГ для значения h0/h = 0,2, кривые 2 и 3 соответственно для значений h0/h = 0,1 и 0,05. Со стороны меньших значений t /t 0 ниже границ устойчивости СГ в области II наблюдается режим регулярных автопульсаций с относительно высокой частотой повторения импульсов (с периодом Т< < t f). Справа к границам устойчивости СГ примыкают зоны нерегулярной (хаотической) генерации III (заштрихованы). При дальнейшем увеличении t /t 0 хаотическая генерация переходит в регулярные автопульсации с периодом близким к времени пролета t f (область IV). Кривая 4 показывает пороговое значение потерь резонатора q , выше которого генерация невозможна.Результаты расчетов характеристик генерации в различных ее режимах представлены на рис. 2-5. Размер зоны неоднородности везде взят равным
h0=0,1h, для величины нормированных потерь принято значение q /a m=0,26. На рис. 2 показаны типичные временные зависимости интенсивности генерации на выходе из резонатора I(h, t). В области II период автопульсаций равен Т=0,17t f (рис. 2а), что заметно превосходит время пролета среды через зону неоднородности t 0=0,1t f. В области IV (рис.2в) генерируются короткие импульсы с периодом следования Т=1,02 t f. (На рисунке 2в их длительность определяется шириной сглаженного края зеркал НР). При переходе от данного режима к хаотической генерации (рис.2б) в узкой переходной зоне между областями III и IV наблюдалось удвоение периода, являющееся одной из промежуточных стадий развития хаоса [12]. Оно выражалось в неглубокой модуляции интенсивности автопульсаций с периодом 2Т.Форма границы устойчивости СГ и наличие двух областей автомодулированной генерации с разными частотами пульсаций указывают на наличие двух различающихся механизмов автомодуляции.
Первый из них реализуется в области IV (рис.1) при
t > > t f, когда восстановление инверсии обусловлено сменой газа в резонаторе. Масштаб пространственной неоднородности стационарных распределений коэффициента усиления a s(x) и интенсивности поля Is(x), с наличием которой связана неустойчивость СГ, в этом случае определяется апертурой резонатора h, а период автопульсаций близок к времени пролета t f.Второй механизм автомодуляции имеет место при
t <<t f в области II (рис.1). Здесь восстановление инверсии в приосевой зоне, в которой накачка отсутствует, происходит за счет процессов накачки в соседних зонах и переноса возбужденных молекул к оси резонатора. Условия на входе потока в резонатор и полное время пролета t f в этом случае не влияют на динамику генерации в "управляющей" приосевой зоне. Период следования импульсов определяется как временем релаксации t , так и временем пролета через зону неоднородности t 0 , причем всегда T>t 0 (Заметим, что в рассматриваемой системе величины t 0 и t не могут сильно отличаться). Отношение t /t 0 однозначно определяет границу устойчивости СГ по величине потерь q , что видно из подобия кривых 2 и 3 (рис.1) в этой области. (Если эти кривые представить в зависимости от t /t f, то они будут близки друг к другу в области t >t f). Что касается кривой 1, то для нее указанное подобие нарушается, так как условие t 0<<t f выполняется хуже.Чтобы убедиться в том, что неоднородность накачки играет определяющую роль в механизме неустойчивости СГ и возникновении режима автопульсаций в области
t <t f, были сделаны расчеты с пространственно-однородной накачкой. Было найдено, что в области t <t f в этом случае всегда наблюдается устойчивая СГ, в то время как в области IV (рис.1) сохраняется автомодуляция излучения.Производился также расчет для НР с однородной накачкой и полностью однородными распределениями коэффициента усиления и поля СГ. В рамках нашей расчетной модели такая система реализуется при использовании "согласованного" значения коэффициента усиления на входе, равного потерям резонатора
a in= q .* Аналитическое исследование устойчивости такой системы показывает, что возмущения в ней всегда затухают. В численных расчетах устойчивая СГ наблюдалась при всех значениях параметров системы. При a in<q для всех значений скорости накачки, при которых достигается превышение порога, также имела место стационарная генерация.Неустойчивость СГ в области II (pис.1) обусловлена видом формирующихся под действием накачки стационарных pаспpеделений коэффициента усиления и поля в приосевой зоне pезонатоpа, которые быстро возрастают от оси резонатора. Такая геометрия системы способствует усилению относительных возмущений
(x, t) и (x, t) в направлении навстречу потоку. Это качественно видно из уравнений (5), (6), в которых члены, ответственные за взаимодействие возмущений, пропорциональны ws и gs, соответственно. Так как сама накачка в приосевой зоне почти отсутствует, то релаксация относительных возмущений инверсии в потоке, описываемая первым членом правой части (5), оказывается сильно замедленной. Благодаря этому перенос к оси HР возникающих на некотором расстоянии от нее более сильных возмущений сpеды пpоисходит достаточно эффективно и способен пpиводить к pаскачке автоколебаний.Таким обpазом, наблюдающаяся неустойчивость СГ в области
t <t f связана с выбоpом специального пpофиля накачки, хаpактеpизующегося pезким спадом в центральной зоне НР. Пpи данном пpофиле влияние накачки на динамику генеpации в большей степени связано с фоpмиpованием наpастающих пpофилей gs(x) и ws(x) и в меньшей степени с pелаксацией возмущений в потоке. Пpи ином пpофиле накачки, в частности, в случае одноpодной накачки главную pоль игpает влияние накачки на затухание возмущений в потоке, что ведет к стабилизации СГ [4, 5].Различное поведение системы в зависимости от пpофиля накачки пpоявляется и в pазличном влиянии на динамику генеpации величины потеpь pезонатоpа
q , опpеделяемых коэффициентом увеличения М. В pяде pабот отмечалось, что уменьшение M, как и увеличение скоpости накачки, способствует стабилизации СГ [4, 5]. Это согласуется с хаpактеpом зависимости pелаксационного члена в (5) от указанных величин. В пpотивоположность этому в HР с пpофилем накачки типа (4) автопульсации пpи t <t f возникают только пpи достаточно малых потеpях (pис.1). В соответствии со сказанным выше, это объясняется возpастанием гpадиентов gs и ws пpи уменьшении q . Заметим, что пpи фиксиpованном отношении a in/a m величины S и q входят в число паpаметpов системы только в комбинации S/q .Пpиведенные качественные сообpажения подтвеpждаются численными pасчетами пpостpанственных pаспpеделений коэффициента усиления и поля в HР с неодноpодной накачкой.
Особый интерес представляет хаотический режим генерации (рис.2б). В отличие от работ, где подобные режимы получались при внешней модуляции параметров системы [10, 11], в данном случае какая-либо модуляция параметров отсутствует. Хаотическая генерация имеет место в промежуточной области соотношения характерных времен системы
t ~ t f (рис.1), где восстановление инверсии и формирование профилей a (x, t) и I(x,t) обусловлено как сменой газа, так и процессами накачки и релаксации. Можно предположить, что усложнение динамического поведения системы в этой области связано с взаимодействием указанных двух механизмов с различающимися характерными временами. На рис.3 приведен фазовый портрет хаотической генерации в координатах I, ¶ I/¶ t. Природа хаоса и сценарий его развития в данной работе специально не исследовались. Наблюдавшаяся на границе областей III и IV (рис.1) модуляция огибающей импульсов позволяет предположить, что развитие хаоса идет через последовательность удвоения периода автопульсаций [12].В серии расчетов моделировалась динамика переходного процесса развития неустойчивости СГ и формирования различных нестационарных режимов генерации. Вначале для выбранных значений параметров системы моделировалась СГ, которая искусственно стабилизировалась путем фиксирования самосогласованного значения интенсивности
I0 на оси резонатора, при котором выполняется условие равенства на оси усиления и потерь. После установления стационарного режима ограничение на величину I0 снималось и далее она рассчитывалась по формуле (3). Результаты расчетов приведены на рис.4. В линейной стадии малые возмущения носят гармонический характер. По мере увеличения их амплитуды появляются нелинейные искажения, при этом период пульсаций несколько возрастает. Дальнейшее поведение системы зависит от выбранных значений параметров.В области II (рис.1) быстро устанавливается режим регулярных автопульсаций (рис.4а). В области III вначале наблюдаются колебания с пеpиодом Т
» 0,2t f , слабо пpомодулиpованные с пеpиодом близким к t f (pис.4б), причем значения обоих периодов не строго кратны друг другу. В нелинейной стадии pегуляpность колебаний наpушается и постепенно pазвивается хаотическая генеpация. (Hа pис.4б показано начало этого пpоцесса. В дальнейшем генеpация пpинимает вид, изобpаженный на pис.2б.) Интересную особенность обнаруживает переходный процесс в области IV. Здесь наблюдается промежуточная стадия нерегулярной генерации, внешне имеющей вид хаоса, которая затем переходит в регулярные автопульсации с периодом Т » t f (рис.4в).На рис.4г более подробно показан начальный этап развития нерегулярностей в этой области. Видно, что на фоне малого гармонического возмущения с периодом Т
= 0,84t f постепенно возникают мелкомасштабные искажения с более высокими инкрементами нарастания. Образующиеся выбросы и провалы повторяются с периодом равным t f. Это говорит о том, что механизм усиления мелкомасштабной структуры связан с наличием резкого градиента (спада) поля на краю апертуры зеркала. По мере усложнения указанной структуры и ее нелинейного искажения генерация приобретает нерегулярный характер. Последующая регуляризация, по-видимому, является результатом нелинейной синхронизации автопульсаций.На рис.5 приведены значения инкрементов развития неустойчивостей СГ в областях I-IV (в расчете на время пролета
t f), которые определялись на линейной стадии малых возмущений. Видно, что вдали от границ устойчивости СГ развитие или затухание возмущений происходит достаточно быстро.Представляет интерес сравнение энергетической эффективности различных режимов генерации. В принятой модели квантовая эффективность использования энергии возбуждения
h для верхней (по потоку) половины НР может быть найдена по соотношению, (7)
где угловые скобки означают усреднение по времени.
Выражение (7) учитывает потери энергии на релаксацию в объеме и на вынос возбужденных молекул. Согласно полученным предварительным результатам переход от искусственно стабилизированной СГ в области неустойчивости последней к режимам автопульсаций или к хаотической генерации при неизменных параметрах системы сопровождается некоторым возрастанием величины
h . Это означает, что осуществление того или иного режима работы системы при заданном наборе ее параметров происходит в соответствии с критерием "энергетической выгодности". При малых размерах зоны неоднородности накачки h0<<h энергетическая эффективность автомодулированной генерации в рассматриваемой системе может быть достаточно высока (h і 0,6) и сравнима с той, которая достигается в НР с неоднородными зеркалами [9].
Заключение
Выполненные модельные расчёты характеристик генерации быстропроточного лазера с неоднородной накачкой внутри НР показали, что, в зависимости от пространственного профиля накачки и параметров системы могут реализовываться различные динамические режимы генерации лазера.
Накачка оказывает двоякое влияние на устойчивость СГ. С одной стороны, ускоряя процесс генерации, она ведет к затуханию возмущений инверсии в среде и стабилизации СГ. С другой стороны, накачка косвенно влияет на устойчивость СГ, изменяя вид стационарных распределений коэффициента усиления и интенсивности поля. Однородная накачка, выравнивая указанные распределения, повышает устойчивость СГ. В противоположность этому накачка с профилем, быстро спадающим к оси НР, за счет увеличения пространственных градиентов поля и среды стимулирует развитие неустойчивости СГ.
В практически важном случае, когда время пролёта намного превосходит время релаксации среды,
t f>>t , динамическое поведение системы полностью определяется видом профиля накачки в “управляющей” приосевой зоне НР.Сложный характер взаимодействия поля с движущейся активной средой в пространственно-неоднородной системе при наличии двух конкурирующих процессов восстановления инверсии (внутренняя накачка и внос возбужденной среды) с различными характерными временами приводит к возникновению в промежуточной области
t ~ t f хаотического режима генерации. Природа хаоса и его характеристики нуждаются в дальнейшем исследовании.В заключение отметим, что влияние профиля накачки на динамику генерации проточного лазера с НР может дать удобный способ управления временными характеристиками лазерного излучения. По сравнению с ранее предложенными способами получения автомодуляционных режимов генерации [8, 9] данный способ обладает определенными преимуществами благодаря относительной простоте технической реализации изменения профиля накачки в лазере с секционированными электродами.
Данная работа выполнена при поддержке Российского фонда фундаментальных исследований (код проекта 96-02-18609).
Московский государственный университет им. М. В. Ломоносова,
физический факультет, 119899, Москва, Воробьевы горы, Россия.
A.I. Fedosejev, A.V. Fedjanovich, A.V. Mushenkov, A.I. Odintsov, N.E. Sarkarov. Operation dynamics of fast flow laser in the case of inhomogeneous excitation within an unstable resonator.
The temporal characteristics of the fast flow CW gas laser has been numerically investigated under the condition of the spatial inhomogeneity of the excitation within an unstable resonator. The nonuniform distribution of the pumping rate across the resonator was found to cause the instabilities of the output power. The instabilities arise due to the transverse gain and power gradients stimulated by the inhomogeneous pumping. Various regimes of the laser action were found in the case of the spatial pumping profile having the abrupt fall near the optical axis. They are the continuous, chaotic and periodical selfpulsing generation depending on the parameters of the laser system.
Рисунки
Рис.1 Области существования различных режимов генерации (пояснения в тексте).
Рис.2 Временные зависимости интенсивности генерации в областях II (а), III (б) и IV (в). Значения t /t f равны: 0,13 (а), 2 (б), 4 (в).
Рис.3. Фазовый портрет хаотической генерации. Параметры те же, что на рис.2 б.
Рис.4. Развитие неустойчивостей и формирование режимов генерации в областях II (а), III (б) и IV (в, г). Значения
t /t f равно 0,11 (а), 2 (б), 4 (в, г).Рис.5. Инкременты развития неустойчивостей СГ.