Ордена ленина




Скачать 303.14 Kb.
НазваниеОрдена ленина
Дата18.05.2013
Размер303.14 Kb.
ТипДокументы



РОИЙСКАЯ АКАДЕМИЯ НАУК

ОРДЕНА ЛЕНИНА

ИНСТИТУТ ПРИКЛАДНОЙ МАТЕМАТИКИ

имени М.В.КЕЛДЫША


А.Е.Коробочкин, В.И.Павлов, А.Х.Пергамент


ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ЛАЗЕРНЫХ НЕУСТОЙЧИВЫХ РЕЗОНАТОРОВ С ЗЕРКАЛАМИ ПРОИЗВОЛЬНОЙ ФОРМЫ В НЕОРТОГОНАЛЬНЫХ КООРДИНАТАХ


Москва

2002


А.Е.Коробочкин, В.И.Павлов, А.Х.Пергамент. ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ЛАЗЕРНЫХ НЕУСТОЙЧИВЫХ РЕЗОНАТОРОВ С ЗЕРКАЛАМИ ПРОИЗВОЛЬНОЙ ФОРМЫ В НЕОРТОГОНАЛЬНЫХ КООРДИНАТАХ.


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

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

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

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


А.Е.Коробочкин, В.И.Павлов, А.Х.Пергамент. The mathematical modeling of the laser beam propagation in the unstable resonator with the mirrors of the arbitrary form in nonorthogonal coordinat systems.


The mathematical modeling of the laser beam propagation in the unstable resonator with the mirrors of the arbitrary form is represented in the article. The method considered is based on using the non-orthogonal coordinate systems, which allows the adequate reduction of full wave boundary problem to the system of the evolution Shredinger’s equations.

The features of the method described are neglecting from the Siegman-Siclash transformations. As a result it doesn’t need to interpolate the field near the back mirror on the detailed grid, that is necessary generally for the solution of the equation describing the radiation propagation after the field reflection from the back mirror.

The spectral methods and iterations have been used for the numerical solution of the boundary problems. The telescope resonator with the cylindrical mirrors has demonstrated the features of the unstable resonators calculations when the active medium has been saturated.


Проблеме распространения лазерного излучения в открытых неустойчивых резонаторах, начиная со статьи Сигмена [1], посвящено большое количество теоретических и экспериментальных работ, основные результаты которых подытожены в известных монографиях [2,3]. Математические модели и основные численные методы изложены в работах [4,5].

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

В данной работе проведено четкое разделение особенностей двумерных и трехмерных задач, получаемых из соответствующих уравнений Гельмгольца [6,7]. Важной особенностью полученных эволюционных уравнений является наличие в них членов, характерных для уравнений переноса. Так, в частном, но важном случае телескопического резонатора в уравнении, описывающем распространение излучения от малого зеркала к большому, появляется член переносного типа.

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

В предлагаемой работе для численного решения краевых задач использованы спектральные методы [8,9] и итерационные процессы. При численном решении наличие достаточно резких изменений поля в резонаторе приводит к необходимости использования достаточно подробных сеток и проведения значительного числа итераций. Это, в свою очередь, приводит к необходимости точного учета кривизны зеркал.

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

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


Equation Section (Next) 1.Рассмотрим открытый оптический резонатор с зеркалами произвольной формы. Поверхности зеркал в декартовой системе координат описываются соответственно гладкими функциями

и , имеющими ограниченные вторые частные производные.

Во внутренней области резонатора электрическое поле удовлетворяет уравнению Гельмгольца:


, MERGEFORMAT

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

В произвольной системе координат уравнение (1) принимает вид:


. MERGEFORMAT

Здесь - метрический тензор, , . По повторяющимся индексам подразумевается суммирование.

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


.

или

. MERGEFORMAT

Для удобства перепишем (3) следующим образом:






. MERGEFORMAT

Здесь




,



,



,

,

,

. MERGEFORMAT


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



. MERGEFORMAT


Пусть электрическое поле излучения, падающего на отражающую поверхность , описывается в координатах . Тогда:


.

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

Электрическое поле излучения, отраженного от поверхности зеркала, описывается в других координатах, а именно имеет вид:

.

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

Тогда из соотношения (6) следует:


.

Окончательно:


. MERGEFORMAT

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

Выберем начало координат в точке фокуса зеркал конфокального неустойчивого резонатора. Примем за ось оптическую ось резонатора. Положительное направление оси - направление от малого зеркала к большому зеркалу.

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

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


.

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


.

Выберем такую систему криволинейных координат , в которой функция переходит в плоскость а функция в плоскость . Искомой системой криволинейных координат является:



где , .
  1. Поле в резонаторе представим в виде:





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


3.1Нечетные проходы.

Пусть источник отраженной от малого зеркала сферической волны находится в начале координат. Представим поле в виде: , где - фундаментальное решение оператора Гельмгольца [6,7]:


.

Здесь .

Компоненты тензора в координатах равны:


, , ,

, ,



.

Определитель матрицы метрического тензора : . Следовательно, согласно (5) имеем: , ,









,

, ,





.

Таким образом, в выбранной системе криволинейных координат уравнение (4) принимает вид:




. MERGEFORMAT

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


, при , .

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

Начальное условие получаем из граничных соотношений (7):




.

Здесь - коэффициент отражения от малого зеркала.

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


, ,

или в криволинейной системе координат :


, . MERGEFORMAT

Следовательно, учитывая условие (9), переходим от уравнения (8) к уравнению:


. MERGEFORMAT

3.2Четные проходы.

Поскольку резонатор телескопический, т.е. волна, отраженная от большого зеркала близка к плоской волне, естественно использовать решение оператора Гельмгольца в виде [6]:


.

Тогда, согласно (5) имеем: , ,








,

, ,

Таким образом, в выбранной системе криволинейных координат уравнение (4) принимает вид:


. MERGEFORMAT

Если выполняются условия (9), переходим в уравнении (11) к квазиоптическому приближению:

. MERGEFORMAT


Обсуждение краевых условий приведено выше. Начальное условие:




.

Здесь - коэффициент отражения от большого зеркала.


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

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

В системе декартовых координат , функция , определенная в области и описывающая поверхность малого зеркала имеет вид:


.

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


.

Выберем такую систему криволинейных координат , в которой функция переходит в плоскость а функция в плоскость . Искомой системой криволинейных координат является:




где , .

Аналогично методу, изложенному для трехмерного случая, поле в резонаторе представим в виде:

.

Здесь поле распространяется от малого зеркала к большому зеркалу (нечетный поход), а - от большого зеркала к малому (четный проход).


4.1Нечетные проходы.

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

Представим в виде , где - фундаментальное решение двумерного оператора Гельмгольца [6]:


.

Здесь - функция Ханкеля первого рода, .

Компоненты тензора в координатах равны:


, ,

,



Определитель матрицы метрического тензора : . Следовательно, согласно (5) имеем: ,




,



, MERGEFORMAT



. MERGEFORMAT

Таким образом, в выбранной системе криволинейных координат уравнение (4) принимает вид:


  1. . MERGEFORMAT

  2. Так как переменные и связаны соотношением , необходимо выразить функции и только в переменных .


Переходя к квазиоптическому приближению в уравнении (15) получаем:


. MERGEFORMAT

Выбор краевых условий обсужден выше. Начальное условие получаем из соотношений типа (7):


.

Здесь - коэффициент отражения от малого зеркала.


4.2Четные проходы.

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


.

Следовательно, согласно (5) имеем: ,



,

, .

Таким образом, в выбранной системе координат уравнение (4) имеет вид:


. MERGEFORMAT

Переход в уравнении (17) к квазиоптическому приближению приводит к уравнению:


. MERGEFORMAT

Выбор краевых условий обсужден выше. Начальное условие:


.

Здесь - коэффициент отражения от большого зеркала.


5.Рассмотрим математическую модель распространения лазерного излучения в генераторе на базе неустойчивого телескопического резонатора, заполненного активной нелинейной средой. Используем широко известную модель Ригрода [10] для описания усиливающей (активной) среды. Эффективный коэффициент усиления в этой модели представляется в виде:

,

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

Тогда процесс распространения лазерного излучения описывается системой нелинейных уравнений для медленно меняющихся комплексных амплитуд и : - нечетный проход излучения, - четный проход излучения.


5.1В трехмерном случае с учетом эффективного коэффициента усиления активной среды имеем систему уравнений:




,



.

Условие отражения от малого зеркала:




.

Условие отражения от большого зеркала:




.

Здесь - коэффициент поглощения, - коэффициент усиления слабого поля, - интенсивность насыщения.


5.2В двумерном случае имеем систему уравнений:




, MERGEFORMAT

. MERGEFORMAT

Условие отражения от малого зеркала:


.

Условие отражения от большого зеркала:


.


6.Рассмотрим подробнее двумерный случай. Коэффициенты и в уравнении (15) зависят согласно (13), (14) как от функции Ханкеля, так и от ее производных, что создает определенные неудобства при численном исследовании задачи распространения лазерного излучения в резонаторе. Воспользуемся приближенной формулой [11]:


, MERGEFORMAT

где :

.

Ограничимся тремя членами в разложении (21):


. MERGEFORMAT

Подставляя (22) в выражения для коэффициентов и , получаем:




,







.

Представим уравнения (19) и (20) в виде, удобном для численного решения. Поскольку справедливы соотношения:


, ,

после тождественных преобразований уравнения (19), (20) перепишем следующим образом:






,



.

Введем безразмерные переменные: , .

В переменных уравнение для нечетного прохода:






.

Уравнение для четного прохода:



.

По переменной применяем преобразование Фурье. Для нечетного прохода получаем:






.

Для четного прохода имеем:






Здесь нижним индексом обозначены соответствующие Фурье компоненты функций, например, , , и т.д.

По переменной применяем схему Кранка-Николсона [12] с итерациями. Для нечетного прохода:








. MERGEFORMAT

Для четного прохода:






. MERGEFORMAT

Уравнения (23)-(24), условия отражения на зеркалах резонатора:


,

.

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


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




где , .

Компоненты тензора в координатах равны:

, ,

,



Определитель матрицы метрического тензора : . Следовательно, имеем: ,

,

, .

Таким образом, в выбранной системе координат уравнение (4) имеет вид:

. MERGEFORMAT

Переходим в уравнении (25) к квазиоптическому приближению:


. MERGEFORMAT

Выбор краевых условий обсужден выше. Начальное условие с учетом апертуры малого зеркала:




Поскольку поле распространяется от поверхности до плоскости вне активной среды, согласно (23) имеем:




. MERGEFORMAT


Решая уравнение (26) получаем установившееся значение поля на произвольной плоскости в декартовой системе координат.

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

Приведены результаты расчета поля в неустойчивом телескопическом резонаторе, имеющем следующие геометрические характеристики: радиус малого зеркала равен 1120мм. Радиус большого зеркала равен 3360мм. Длина волны затравочного излучения в генераторе – 10.7мкм. Расчет наиболее добротных мод и соответствующих им собственных значений проведен для диапазона изменения чисел Френеля от 0.5 до 5.0.

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

На Рис.1.б изображено поле основной моды в ближней зоне в декартовой системе координат .

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

















Л ИТЕРАТУРА



  1. Siegman A.E. Proc. IEEE, 1965, v. 53, p. 277.

  2. Ананьев Ю.А. Оптические резонаторы и проблемы расходимости лазерного излучения. – М.: Наука, 1979.

  3. Ананьев Ю.А. Оптические резонаторы и лазерные пучки. – М.: Наука, 1990.

  4. Карамзин Ю.Н., Конев Ю.Б. Численные методы исследования работы неустойчивых телескопических резонаторов с учетом дифракции и эффекта насыщения в активной среде. – Препринт ИПМ им.М.В.Келдыша АН СССР, 1974, № 57.

  5. Елкин Н.Н., Напартович А.П. Прикладная оптика лазеров. – М.:ЦНИИатоминформ, 1989.

  6. Владимиров В.С. Уравнения математической физики. - М.: Наука, 1971.

  7. Тихонов А.Н., Самарский А.А. Уравнения Математической физики. – М.: Наука, 1972.

  8. Карамзин Ю.Н., Цветкова И.Л. Спектральный метод решения нелинейных квазиоптических задач. – Препринт ИПМ им.М.В.Келдыша АН СССР, 1979, № 115.

  9. Карамзин Ю.Н., Цветкова И.Л. О сходимости спектрального метода решения одной задачи нелинейной оптики. – ЖВМ и МФ, 1982, т.22, с. 235-240.

  10. Rigrod W. Gain saturation and output power of optical mazers. – J. Appl. Phys., 1963, v. 34, p. 2602.

  11. Ватсон Г.Н. Теория бесселевых функций, часть первая. – М.: Из-во Иностранной литературы, 1949.

  12. Самарский А.А. Теория разностных схем. – М.: Наука, 1977.

  13. Фаддеев Д.К., Фаддеева В.Н. Вычислительные методы линейной алгебры. – М.-Л.: ГИФМЛ, 1963.

  14. Уилкинсон Дж. Х. Алгебраическая проблема собственных значений. – М.: Наука, ГИФМЛ., 1970.

  15. Siegman A.E., Miller H.Y. Unstable Optical Loss Calculations Using the Prony Method. – Applied Optics, 1970, v. 9, № 12, p. 2729.














Похожие:

Ордена ленина iconМолекулярные аспекты функционирования рибосомных рнк
Московский ордена ленина, ордена октябрьской революции и ордена трудового красного знамени государственный университет имени М. В....
Ордена ленина iconОтчете, подлежит раскрытию в соответствии с законодательством Российской Федерации о ценных бумагах
Открытое акционерное общество "Ярославский ордена Ленина и ордена Октябрьской Революции шинный завод"
Ордена ленина iconМетодические указания к учебно-исследовательской лабораторной работе №13 по дисциплине «охрана труда»
Московский ордена ленина ордена трудового красного знамени институт инженеров железнодорожного транспорта
Ордена ленина iconОтчет открытое акционерное общество «Ярославский ордена Ленина и ордена Октябрьской Революции шинный завод»
Открытое акционерное общество «Ярославский ордена Ленина и ордена Октябрьской Революции шинный завод»
Ордена ленина iconОрдена Ленина Институт прикладной математики им. М. В. Келдыша удк 517. 958 М. П. Галанин, А. П. Лотоцкий, В. Ф. Левашов
Расчет электродинамического ускорения плоских пластин в лабораторном магнитокумулятивном генераторе
Ордена ленина iconМосковский ордена ленина государственный университет
Л 92001. Подписано к печати 24/1 1948 г. Зак. №1100. М. Г. У. №503—47. Печ л. 66V2 +8 вклеек = 703/4 п л. Уч ав л. 6072. Тираж 3000....
Ордена ленина iconРоссийская Академия Наук ордена ленина институт прикладной математики им. М. В. Келдыша А. В. Ахтёров, А. А. Кирильченко
Задача обнаружения подвижных объектов при информационном мониторинге динамической среды распределённой мобильной системой
Ордена ленина iconОрдена ленина главмосстрой при мосгорисполкоме мосоргстрой технологическая карта
Технологическая карта на устройство внутриквартальной фекальной канализации из полиэтиленовых труб
Ордена ленина icon,43 + 549,3 отьшдные здщю-никшене руды йоко-дошренского расслоешюго їшут0на (состав и условия образования) Специальность
Академия наук СССР ордена Ленина Сибирское отделение Бурятский филиал Геологический институт
Ордена ленина iconОрдена Ленина Институт прикладной математики им. М. В. Келдыша Галанин М. П., Щеглов И. А
Делоне и методы исчерпывания. Приведены варианты алгоритмов для каждого из указанных методов. Обсуждены особенности построения сеток...
Разместите кнопку на своём сайте:
Библиотека


База данных защищена авторским правом ©lib.znate.ru 2014
обратиться к администрации
Библиотека
Главная страница