Численное исследование методов поиска многомерных солитонов Текст научной статьи по специальности «Математика»

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

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

Библиографические ссылки

1. Рубан А. И. Глобальная оптимизация методом усреднения координат / Краснояр. гос. техн. ун-т. Красноярск, 2004.

2. Кузнецов А. В. Генераторы последовательностей равномерно распределенных точек // Вестн. Краснояр. гос. техн. ун-та. Мат. методы и моделирование. 2005. № 37. С. 148-155.

References

1. Ruban A. I. Globalnaya optimizatsiya metodom us-redneniya koordinat (Global optimization by averaging coordinates). Krasnoyarsk, KGTU, 2004. 302 p.

2. Kuznetsov A. V. Vestnik KGTU, 2010, № 37, рр. 148-155.

© Рубан А. И., Кузнецов А. В., 2013

УДК 519.634

ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ МЕТОДОВ ПОИСКА МНОГОМЕРНЫХ СОЛИТОНОВ

Н. П. Савенкова, В. С. Лапонин

Московский государственный университет имени М. В. Ломоносова Россия, 119991, ГСП-1 Москва, Ленинские горы. E-mail: [email protected]

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

Ключевые слова: солитонное решение, численный метод, сходимость.

NUMERICAL RESEARCH OF METHODS FOR MULTIDIMENTIONAL SOLITONS SEARCH

N. P. Savenkova, V. S. Laponin

Moscow State University named after M. V. Lomonosov GSP-1 Lenin Hills, Moscow, 119991, Russia. E-mail: [email protected]

Over the last years, along with the development of optical computing systems, the mathematical modeling of laser femtosecond pulses in nonlinear media has been actively developing. Among the proposed methods for finding solitons the most widespread methods are the method of the inverse problem, spectral methods and other methods. The authors present two iterative methods, namely M1 and M2, for finding soliton solutions in the propagation of optical radiation in a medium with a cubic nonlinearity in the axially symmetric case.

Keywords: soliton solution, numerical method, convergence.

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

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

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

частотах одновременно существуют и распространяются вместе оптические волны вдоль нелинейной среды. Эволюция этих солитонов описывается системами связанных уравнений Шредингера. Интерес к этим солитонам в литературе постоянно сохраняется в связи с многочисленными потенциальными приложениями их в задачах передачи информации оптическими методами. Экспериментальные исследования [1-3], выполненные в различных лабораториях, показали практическую реализацию солитонов данного типа.

Ниже предлагаются два итерационных метода М1 и М2 для поиска солитонных решений, подробно описанные в [5] и [6] соответственно, которые хорошо себя зарекомендовали при нахождении солитонных решений одномерных уравнений таких, как уравнение Кортевега-де Фриза, уравнение 8ш-Гордона и нелинейного уравнения Шредингера. В настоящей работе методы М1 и М2 применяются к системе нелинейных уравнений Шредингера, описанной в [7], и приводится сравнение результатов, полученных с помощью предложенных методов, с результатами, полученными в работе [7].

Математическая постановка задачи. Распространение оптического излучения в среде с кубичной нелинейностью в аксиально-симметричном случае описывается следующей системой нелинейных уравнений Шредингера [7]:

дА

д2 А

— + ЮД ГА + іБ1^- + 1уА" Бе—Ь +

дг

ді 2

+ іаА(| А |2 +2 | Б |2) = 0,

дБ .Ю. 0 .п д2Б -----+і—Д ГБ +1&2 —2

дг 2 г 2 ді2

+ і чА2Єкі +

(1)

+ 2іаБ(2| А |2 +1 Б |2) = 0,

0 < г < Я, 0 < г < 1,

0 < і < Т, к = к2 - 2к1,

Д г =-

На входе в нелинейную среду задается распределение импульса:

А(Л, г, г = 0) = А0 (Л, г),

Б(Л, г, г = 0) = Б0(Л, г).

Граничные условия имеют вид

А \і =0,Т = 0,

дА

г — |г—>0 = 0, А \г=я = 0,

дг

Б Л=0,т = 0, гд 1г ^0 = 0, Б 1г=я = °. дг

Обобщение итерационного метода M1 к решению многомерной задачи. В системе (1) сделаем замену переменных С = Бв~гЬ и получим нелинейную систему уравнений

дА — д2 А . *

----+/БДГА +/Б1 —~—+ /уА С +

г 1 дл2

дг

+ іаА(| А |2 +2 | С |2) = 0,

дС Ю

~ + ікС +Д гС + іЮ дг 2 г

д2С . ,2 2 —г + чЛ +

(2)

+ 2/а С (2 | А |2 +1 С |2) = 0.

Далее произведем замену переменных

|^1 = г,

1^2 = Л - У2, и получим следующую систему уравнений

-V— + ЮД Е А + іЮ + іуА*С +

щ>2 51 1 д^2

(3)

где А( г, г, і), Б( г, г, і) - комплексные амплитуды гармоник; г - нормированная на длину пучка основной волны продольная координата; г - безразмерная координата; і - безразмерное время в сопровождающей импульс основной волны системе координат; 1 д ( д

-I г— | - оператор Лапласа по координате г ;

+ іаА(| А |2 +2 | С |2) = 0,

д 2С

дС Ю К п

—V------+ ікС +1—Д? С + іЮл —~—+

д^2 2 2 %2

(4)

г дг ^ дг

1 дк,

Б, =-----------— - коэффициенты, характеризующие дис-

' 2 дю, _ персию второго порядка , = 1,2 ; Б - коэффициент, характеризующий дифракцию (в выбранной нормировке Б = 1); к,, ю, - волновое число и частота волны соответственно , = 1,2; у - коэффициент нелинейной связи взаимодействующих волн; а, - коэффициенты самовоздействия волн , = 1,2; 2 - безразмерная длина нелинейной среды; я - поперечный размер нелинейной среды; т - безразмерное время, в течение которого анализируется рассматриваемый процесс.

+ іуА2 + 2іаС(2 | А |2 + | С |2) = 0.

В системе (4) проинтегрируем первое и второе уравнение по |2 и получим систему

уА = /Б [ Д£ А^2 + /Б1— +

д^2

+ /у| А* Сё|2 + /а| А(| А |2 +21 С12)ё|2,

уС = / [кСё§2 + / — \Д,Сё^2 + /Б2 — +

•' 2 : 1 д|2

+ /у| А2ё§2 + 2/а|С(21А |2 +1 С |2) ё§2.

Задачу (5) будем решать разностным методом в области О = {-Ь < |1 < Ь, -Ь < |2 < Ь}. Величина Ь

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

дС

(5)

ченного по столбцам. Далее для системы (5) запишем итерационный процесс:

(

АП+1 = АП + т

—£ д^ АпкН2 + Б к=1

Ап - Ап

Аз А,-1

+/у£ (АП )*С1Н2 + /а^А (| А |2 +21 СП ^ ^

к=1 к=1

/'~»п+1 ,оп . _

С, = С.^- +т

^/к £^7 + £ Д^1 С^П^2 +

V т=1 2 к=1

СП - Сп

+ /Б~,

3 з-1 I V,

'уЕ (Ак”+1)2 к2

к=1

+ 2/а£Скп(2 | Ап+1 |2 +1 Сп |2)^2

по следующей формуле:

уп+1 =1

Г (ап+1, АП) -1

(Ап, Ап)

л

Итерационный процесс (6) прекращает свою работу, когда будет выполнено условие

уП+1 - уП

< 6.

(6)

к=1

где

Д Аз = 72 ё“(^ з+0,5 Аз +1 -к

- (^,+0,5 +^,-0,5 ) А, +^,-0,5 Ау-1),

А0,С0 - задаются изначально, уп+1 вычисляется

(7)

Результаты расчетов. В качестве начальных приближений А0, С0 для действительной и комплексной части брались функции вида «домик» (рис. 1) и функция «параллелепипед» (рис. 2) или функция «двойной домик» (рис. 3) для получения солитонного решения второй моды.

Расчеты производились при следующих значениях управляющих параметров: Б1 = 0,08, Б2 = 0,14, Б = 0,1, а = 5, у = 20 . Шаг по пространству 71 = к2 = 0,02 , параметр т подбирается экспериментально, в данном

случае т = 10-5. Результат работы метода изображен на рис. 4, а распределение фазы волны представлено на рис. 5. Время численного расчета составило 40 минут, число итераций п = 2 521, невязка у = 0,04. Также на рис. 6 приведено солитонное решение второй моды, для которого шаг по пространству был взят 7 = к2 = 0,01, время расчета составило около 5 часов, число итераций п = 6 725, невязка у = 0,02. Стоит отметить, что для получения солитонного решения второй моды можно использовать в качестве начального приближения функцию «двойной параллелепипед».

Рис. 1.

Рис. 2

Рис. 3

Рис. 4

п

у

2

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

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

Метод М1 был реализован на многопроцессорных компьютерах с использованием программного средства МР1. Решение исходной задачи методом М1 при к = 0,01 на различных машинах дало следующие результаты:

ПК (1 процессор (1 ядро), общая память) - 140 мин.

ПК (2 процессора (4 ядра), общая память) - 45 мин.

Кластер (16 процессоров (32 ядра), разделенная память) - 8 мин.

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

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

§1 = г,

§2 = Л - У^

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

укАЩ+1 = рп (Ащ, Сщ) АП+1,

Лп+1 _ ПП(~ЛП+17^п\7^п+1

УкСк = V- (Ак , Ск )Ск , где Рп и Qn - матрицы размерности N2 х N2, полученные из (6) АЩ - к-й собственный вектор длины N2; уП - соответствующее собственное значение.

Алгоритм М2 для данной системы будет выглядеть следующим образом

1) по АЩ и СП строим матрицу Рп;

2) находим собственное значение уЩ1 и собственный вектор АЩ+1 матрицы Рп;

3) по АЩ+1 и СП строим матрицу Оп;

4) находим собственное значение уЩ1 и собственный вектор СЩ+1 матрицы Оп ;

5) сравниваем результаты на данной и предыдущей итерациях.

А0, С0 - задаются изначально. Метод прекращает свою работу, когда будет выполнено условие

| уЩ+1 - УкП|<6.

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

Б1 = 0,08, Б2 = 0,14, Б = 0,1, а = 5, у = 20. Шаг по пространству к1 = к2 = 0,1, параметр т подбирается экспериментально, в данном случае т = 10-5. Результат работы метода изображен на рис. 7. Метод завершил работу за 15 минут, произведя 27 итераций, невязка составила 0.13, распределение фазы волны соответствует рис. 5. Отметим также, что если в ходе итерационного метода М2 рассматривать второе собственное значение и соответствующий ему собственный вектор, то можно получить солитонное решение второй моды, аналогично для к-й моды.

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

Рис. 7

преобразований и могут быть применимы по стандартной схеме, описанной в [5; 6], что значительно упрощает решение задачи. Особенностью методов M1 и M2 является сходимость к солитонному решению для любых начальных приближений, схожих с описанными выше, что позволяет обойти проблему выбора начального приближения, которую не удается так легко избежать в работе [7].

Для решений задачи (1) на сетке, предложенной выше с шагом h = 0,01 на ПК, методу M1 требуется около двух часов, а методу, предложенному в [7], требуется более суток. Метод выигрывает в 10 раз на ПК, а если использовать параллельные вычислительные системы, то можно добиться выигрыша в несколько сотен раз. Если учесть, что данную задачу необходимо решать для большого количества значений управляющих параметров а, у, чтобы найти область существования солитонных решений, то метод M1 с использованием параллельных вычислительных систем позволит найти область за несколько дней, в зависимости от области и шага.

В заключение отметим, что предлагаемые методы M1 и M2 можно успешно использовать при математическом моделировании оптических компьютеров. Возникающие при определенных значениях управляющих параметров солитонные решения, сохраняющие форму, соответствует сохранению формы сигнала [4]. Отсутствие солитонного решения (нарушение сходимости методов) соответствует потере сигнала, т. е. сбою передачи данных в оптическом компьютере.

Библиографические ссылки

1. Brull L., Lange H. Stationary, oscillatory and solitary wave type solution of singular nonlinear Schrodinger equations // Math. Meth. in Appl. Sci. 1986. Vol. 8, № 4. P. 559-575.

2. Rozanov N. N., Fedorov S. V., Shatsev A. N. Incoherent weak coupling of laser solitons // Optics & Spectroscopy. 2007. Vol. 102, № 1. P. 83-85.

3. Кернер Б. С., Осипов В. В. Автосолитоны. М. : Наука, 1991.

4. Кившарь Ю. С., Агравал Г. П. Оптические соли-тоны. М. : Физматлит, 2005. С. 622-634.

5. Лапонин В. С., Савенкова Н. П., Ильютко В. П. Численный метод поиска солитонных решений // Прикладная математика и информатика : сб. фак. ВМК. 2011. № 38. С. 69-80.

6. Дорохова Т. В., Савенкова Н. П., Трофимов В. А. Численное моделирование солитонных решений в задаче распространения фемтосекундного импульса в среде с кубической нелинейностью // Прикл. математика и информатика : сб. фак. ВМК. 1999. № 2. С. 63-68.

7. Матусевич О. В., Трофимов В. А. Численный метод нахождения 3Б-солитонов нелинейного уравнения Шредингера в аксиально-симметричном случае // Журн. вычисл. математики и мат. физики. 2009. Т. 49, № 11. С. 1-12.

Referens

1. Brull L., Lange H. Stationary, oscillatory and solitary wave type solution of singular nonlinear Schrodinger equations. Math. Meth. in Appl. Sci. 1986, vol. 8, № 4, pp. 559-575.

2. Rozanov N. N., Fedorov S. V., Shatsev A. N. Incoherent weak coupling of laser solitons // Optics & Spectroscopy. 2007, vol. 102, № 1, pp. 83-85.

3. Kerner B. S., Osipov V. V. Avtosolitony (Autosolitons). Moscow, Nauka, 1991.

4. Kivshar Yu. S., Agraval G. P. Opticheskiye solitony (Optical solitons). Moscow, Fizmatlit, 2005, pp. 622-634.

5. Laponin V. S., Savenkova N. P., Ilyutko V. P. Prikladnaya matematika i informatika. 2011, № 38, pp. 69-80.

6. Dorokhova T. V., Savenkova N. P., Trofimov V. A.

Prikladnaya matematika i informatika. 1999, № 2, pp. 63-68.

7. Matusevich O. V., Trofimov V. A. Zhurnal vychis-litel'noy matematiki i matematicheskoy fiziki. 2009, vol. 49, № 11, pp. 1-12.

© Савенкова Н. П., Лапонин В. С., 2013