Численное моделирование самосогласованных структур в плазме и электронных потоках тема диссертации и автореферата по ВАК РФ 05.13.18, доктор физико-математических наук Юнаковский, Алексей Дмитриевич

• Юнаковский, Алексей Дмитриевич
• доктор физико-математических наук
• 2004

Оглавление диссертации доктор физико-математических наук Юнаковский, Алексей Дмитриевич

1 Сеточные методы решения одномерного нестационарного уравнения Шредингера.

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

1.2 Применение операторного компактного неявного метода для решения нестационарного уравнения Шредингера.

1.3 Схема повышенной точности.

2 Сеточные методы для двумерного нестационарного уравнения Шредингера

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

2.2 Применение метода Бунемана для решения нестационарного уравнения Шредингера.

2.2.1 Применение метода Бунемана в схеме расщепления при решении нестационарного уравнения Шредингера

2.2.2 Соединение операторного компактного неявного метода с методом Бунемана.

2.2.3 Случай цилиндрической симметрии.

2.3 Численное моделирование двумерной сильной леигмюровской турбулентности

3 Основные используемые функциональные методы

3.1 Метод осредненной операторной экспоненты для решения задачи Коши.

3.2 Проблемы атракторов траекторий непрерывного спуска для квадратичной формы полуограниченного оператора на единичной сфере

3.3 Применение преобразования всплесков при решении краевых задач.

3.3.1 Обозначения и определения конструкций всплесков

3.3.2 Постановка задачи.

3.3.3 Модельная задача.

3.3.4 Краевая задача с правой частью из W,L.

3.3.5 Краевая задача с правой частью из Vj.

4 Применение быстрого преобразования Фурье для решения нестационарного уравнения Шредингера

4/1 Приближенное решение НУШ во всём пространстве.

4.2 Приближённая операторная экспонента и схема расщепления по физическим процессам

4.3 Краевые задачи для нестационарного уравнения Шредингера

4.3.1 Постановка краевых задач.

4.3.2 Сохранение интегралов движения.

4.3.3 Приближенная операторная экспонента, оценка погрешности

4.4 Неустойчивость периодических решений одномерного иели

Ф иейного уравнения Шредингера.

4.5 Множественное дробление волновых структур в нелинейной среде.

5 Численное моделирование уравнений Захарова

5.1 Одномерные уравнения Захарова.

5.1.1 Сохранение интеграла энергии.

5.1.2 Тестовый пример.

5.2 Уравнения Захарова со сферической симметрией.

5.3 Двумерные векторные уравнения Захарова

5.3.1 Примеры и их анализ.

5.4 Моделирование уравнений Захарова с затуханием.'

5.4/1 Анализ способов введения затухания в уравнения Захарова 5.4.2 Уравнения Захарова с модельным затуханием Ландау

5.4.3 Построение приближенного решения уравнений Захарова с затуханием методом расщепления по физическим процессам.

5.4.4 Построение алгоритма с предиктором

5.4.5 Численный пример

5.4.6 Анализ результатов.

6 Моделирование пространственно-временных структур электронных потоков и ЭМ полей

6.1 Пространственные структуры электронных потоков в скрещенных полях.

6.1.1 Постановка математической задачи.

6.1.2 Ветвление решения.

6.1.3 Бриллюэновский поток электронов.

6.1.4 Виртуальный катод.17.

6.1.5 Решение в пролетном пространстве.

6.1.6 Случай магнитного поля, зависящего от радиальной координаты.

6.1.7 Ветвление решения в случае переменного магнитного поля.

6.1.8 Релятивистский случай.

6.2 Численное моделирование нестационарных процессов в лазерах на свободных электронах.

6.2.1 Уравнение взаимодействия электронного пучка с электромагнитной волной.

6.2.2 Вычислительная схема "FOAM" для интегрирования уравнения движения электронов.

6.2.3 Численное интегрирование уравнения для волны

6.2.4 Результаты тестирования и вычислительных экспериментов

6.3 Вычислительная схема "FROAM" для интегрирования уравнения движения электронов.

7 Обратные задачи для периодически гофрированных волно

Ш водов

7.1 Общая постановка задачи.

7.2 Азимутально-симметричный случай.

7.2.1 Сопротивление связи.

7.2.2 Нормальная составляющая электрического поля на поверхности гофрированного волновода.

7.2.3 Потери в стенках гофрированных волноводов

7.3 Применение метода гармонических возмущений к расчету периодически гофрированных волноводов.

7.3.1 Постановка задачи.

7.3.2 Круглый волновод.

7.3.3 Численные результаты.

7.3.4 Обратная задача.

8 Моделирование ускорительной секции суперколлайдера

Ъ 8.1 Постановка математической задачи рассеяния

8.2 Проблемы аналитического решения.

8.3 Задачи рассеяния в подводящих каналах.

8.3.1 Самосогласованная задача в каналах.

8.4 Оптимизация профиля каналов в в асимптотическом пределе плоской геометрии.

8.5 Основные идеи метода дискретных источников и особенности их применения для расчетов ускорительных секций суперкол-лайдеров

8.6 Функции Грина объемлющей области.

8.6.1 Случай граничных условий первого рода.

8.6.2 Случай граничных условий третьего рода.

8.6.3 Случай условий излучения Зоммерфельда

8.6.4 Случай плоской геометрии фк 8.7 Модельная задача рассеяния электромагнитных волн в канале коллайдера в случае плоской геометрии.

8.7.1 Метод дискретных источников в модельной задаче рассеяния

8.7.2 Процедура регуляризации.

8.7.3 Численные результаты и выводы

8.8 Обратная спектральная задача.

8.8.1 Тестовая задача

8.8.2 Численные Результаты.

8.9 Оптимизация структуры поля в параксиальной области

Введение диссертации (часть автореферата) на тему «Численное моделирование самосогласованных структур в плазме и электронных потоках»

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

Можно в самом общем случае сформулировать основные трббовннIIм. предъявляемые к алгоритмам при численном моделировании. К ним относятся: универсальность метода и его высокая разрешающая способность, что позволяет от идеализированных задач переходить к более реальным: высокая скорость при снижении точности и принципиальная возможность повышения точности до любой необходимой величины; алгоритмическая простота, удобство для программирования, доступность для программистов любой квалификации. Конечно, переход к ЭВМ нового поколения всякий раз существенным образом меняет наши представления о математических моделях, вычислительных методах, алгоритмах и программах. Однако общие требования простоты, универсальности, эффективности всей цепочки модель - алгоритм - программа оставались всегда неизменными. Менялись только пути их достижения. Упрощались алгоритмы и программы и совершенствовались модели, приближаясь к реальным физическим условиям.

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

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

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

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

Слово "структура" происходит от латинского устройство, построение. Понятие структуры широко распространено в химии, биологии, физике и т.п. Исследование процессов возникновения пространственных структ)ф и их разрушения — переход к хаосу — является фундаментальной задачей не только для физики плазмы и механики жидкости и газа, но и вообще для нелинейной физики. Изучению структуры газового разряда, структуры поля в волноводах, структуры волн, двигающихся без искажения - солитонов, -структуры пучка электронов и пространственного заряда, создаваемого ими в электронных приборах, структуры каверн при образовании и обрушении коллапсов плазменных волн, составляющих предмет исследования диссертационной работы, посвящено огромное количество работ. Разнообразие назначений и применений подобных структур и значительная математическая сложность описания происходящих в них физических процессов требуют постоянных теоретических и экспериментальных исследований. Если удается описать такую структуру аналитически, то полученное выражение является решением линейного или нелинейного уравнения или системы уравнений в частных производных. Чаще всего получившаяся функция обладает свойством бесконечной дифференцируемое™ или аналитичности. При создании эффективных алгоритмов поиска приближенного решения этих уравнений необходимо учитывать эту колоссальную априорную информацию. Традиционные методы конечных разностей и конечных элементов, применяемые при этом, почти не используют информацию о гладкости решения, т.е. это так называемые методы с насыщением.

Термин "насыщение" введен К.И.Бабенко [20] . Явление насыщения метода приближения характеризуется тем, что пока гладкость приближаемой функции невелика, приближение в некотором выбранном классе по порядку равно наилучшему, т.е. улучшить порядок приближения нельзя. Этот класс является тем рубежом, при переходе через который дальнейшее увеличение гладкости не может повлиять на порядок приближения, и он остается таким же, как для функций этого класса. Его называют классом насыщения, а величину порядка приближения -порядком насыщения метода приближения.

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

Цитируя доклад Ингрид Добеши на Международном конгрессе математиков в Цюрихе, 1999г., можно сказать: "У математиков есть разные способы оценки новых теорий и конструкций. Одним важным критерием является эстетика - некоторые разработки просто "ощущаются" правильными, подходящими и красивыми. . Другим важным показателем ценности того или иного раздела математики является мера его полезности в приложениях: это почти единственный критерий, используемый нематематиками'.

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

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

Средством решения всех этих задач служат новые разработанные алгоритмы, обладающие указанным свойством "ненасыщаемости".

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

Л Переход к расчетам в образах Фурье и использование свойств этого математического аппарата делают его эффективным средством решения как нестационарных, рассмотренных в четвертой и пятой главах, систем уравнений в частных производных, так и стационарных, рассмотренных в седьмой главе. "Ненасыщаемость" алгоритма или, более точно, "слабая насыщаемость" алгоритма вычислений является следствием того факта, что при удвоении, как это принято в численном анализе, числа дискретных точек с N до 2N, точность нахождения, например, последнего из коэффициентов конечной суммы, пропорциональная отброшенному ряду, резко возрастает и для iV-ro коэффициента становится пропорциональной коэффициенту с номером 3N [208, 225]. Это свойство позволяет считать нестационарные задачи с большим числом шагов по времени. Для задач, рассмотренных в четвертой и пятой главах, где изучаются структура решения с особенностью, развивающейся за конечное время, число шагов возрастает за счет дробления при приближении к точке особенности. Подобные характерные особенности поведения решений этого уравнения требуют тщательного отбора классов их конечномерных приближений и способов описания эволюции этих прмбли-■■0 жений. Кроме того, выяснилось, что структуры у решений с особенностью особенно ярко проявляются именно в образах Фурье решения, где можно изучать возникающий "структурированный хаос".

При изучении структуры пространственного заряда, рассмотренного в шестой главе, была найдена "ненасыщаемая" аппроксимация, т.е. найдено подпространство, базис которого состоит из всплесков с компактным носителем - всплесков Добеши и Койфлетов [253, 254]. Процесс: ветвления решения в этом базисе сводится к переходу главной части соответствующего решения на более высокий уровень всплесков и обращению в нуль некоторых коэффициентов разложения на новом уровне. Краткое изложение кратномасштабного анализа всплесков приведено в разделе 3.3.

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

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

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

Принципы построения приближенного решения, развиваемые в диссертационной работе

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

Существует много разных методов и численных схем нахождения приближенного решения нестационарных систем уравнений в частных производных, начиная от классического метода прямых (см., например, [211]), когда дифференциальный оператор по пространственным переменным заменяется конечно-разностным, а получившаяся при этом конечномерная система обыкновенных дифференциальных уравнений решается методом типа Рунге-Кутта (сюда же относится и схема предиктор - корректор), - и кончая комбинированным или гибридным методом, когда на конечной стадии образования особенности моделирование ведётся методом многих частиц [257|.

В случае двух пространственных переменных часто используется метод переменных направлений. Применение этого метода к нестационарному уравнению Шредингера описано в книге [186]. Однако у этого метода есть явно выделенное "первое" пространственное направление, что мешает на больших временных интервалах сохранять заложенную в начальных данных симметрию, а также контролировать результаты, когда структура получаемого решения теряет изначальную симметрию и становится слоистой. Пример такого поведения решения рассмотрен в разделе 2.1

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

Использованный в работе метод "расщепления по физическим процессам" [161] состоит в сведении исходной эволюционной задачи, описывающей сложный физический процесс, к решению последовательности задач, описывающих процессы более простой физической природы. Приближённо это удаётся сделать на основе аддитивности рассматриваемых процессов в малом. В четвертой и пятой главах показано, что схема "расщепления по физическим процессам" строится путем последовательных операторных замен неизвестной функции на малом интервале. Как спектральные методы, развиваемые в четвертой и пятой главах, так и сеточные, применявшиеся в первой и второй главах, могут быть использованы для реализации отдельных операторов схемы расщепления. Этот путь позволяет наиболее просто сделать обобщения на случай двух пространственных переменных. При проведении замен на первом этапе мы стремимся осуществить переход от уравнения или системы уравнений с неограниченными по пространственным координатам операторами к уравнениям с ограниченными операторами. Такой переход существенно расширяет возможности последующего численного моделирования.

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

На заданном малом отрезке (to,to + At) мы стремились построить такое приближенное решение задачи, которое кроме высокого порядка аппроксимации обладает ещё и тем свойством, что в начальной и конечной точках рассматриваемого интервала производные по времени приближенного решения совпадают с производными точного решения, определяемыми по исходному уравнению. Естественно, что в конечной точке t0 + At мы попадём на другую, близкую траекторию, но с производной, соответствующей Е h 2 t0 +А/

Рис. 0.1 Характер построенного приближенного решения именно этой траектории ( см. рис. 0.1. ). При этом ошибка на шаге, т.е. разность между точным и приближенным решением обладает в начальной и конечной точках шага счета нулевой производной по времени. В сочетании с высоким порядком аппроксимации это свойство, как и при сплайн - аппроксимации, существенно уменьшает накопление ошибки. При этом если у решения уравнения Е(т) , to < т < to + At есть инвариант 1(Е(т)) = const , то автоматически выполняются равенства для приближенного решения Еар{т) в начальной и конечной точках рассматриваемого интервала.

При построении неявного алгоритма мы стремились в случае нелинейного уравнения получить в точке tQ + At наиболее простую алгебраическую (трансцендентную) систему уравнений относительно E(to + At) и отыскать ее решение. Наиболее эффективными являются алгоритмы, в которых удается найти аналитическое решение каждого уравнения схемы расщепления.

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

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

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

3. Одной из целей диссертационной работы являлся учет характерных особенностей поведения решений НУШ уравнения для тщательного отбора классов их конечномерных приближений и способов моделирования эволюции этих приближений. Изучение получаемых в результате эволюции структур является еще одной целью, для достижения которой потребовалсь разработка целого ряда "слабо насыщаемых" алгоритмов, дающих возможность реализации большого числа временных шагов.

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

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

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

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

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

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

Экономичность вычислительного процесса в широком смысле слова может быть понята только в рамках всей технологической цепочки вычислительной математики. Главный недостаток ЭВМ состоит в том, что каждая задача, для которой она может быть применена, должна быть приведена, часто довольно утомительным способом, к последовательности арифметических задач. С учётом больших затрат на программирование, окончательная оценка качества алгоритма должна определяться не только количеством времени расчёта на ЭВМ, но и затратами человеческого времени на программирование, отладку, тестирование и усовершенствование программы. Кроме того, работа программиста, отлаживающего большую программу, напоминает работу авиадиспетчера: и тот и другой должны обладать способностью видеть "картинку" целиком, интуитивно предвидя возможные ослол<нения. Особенно важно наличие такого "диспетчера" при работе коллектива программистов над программным комплексом.

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

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

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

2. Разработаны новые схемы счета для двумерных по пространству систем уравнений, включающих в себя нестационарное уравнение Шредингера, обладающие свойством "слабой насыщаемости" за счет комбинации операторного компактного неявного метода с методом Бунемана, что позволило повысить порядок аппроксимации, сохранив минимальное число операций (NlnN) на каждом временном шаге.

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

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

Дородницын А.А. Информатика: предмет и задачи. - Природа, 1985. No 2, с. 26 - 29.

5. Разработаны новые операторные методы нахождения приближенных решений нелинейного уравнения Шредингера, в том числе и описывающего самовоздействие широкого класса волн, поверхности волновых векторов у которых имеют седловую точку (например, гравитационные волны на глубокой воде [110] , плазменные колебания в замагниченной плазме [152] ).

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

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

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

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

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

Диссертация состоит из введения, восьми глав и заключения.

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Юнаковский, Алексей Дмитриевич

Выводы и перспективы гл. 8

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

Во-вторых, фазовые отклонения численного решения от асимптотического при dmin < А/8 составляют проценты.

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

3. Разработан и реализован алгоритм метода дискретных источников решения двумерного эллиптического уравнения типа уравнения Гельмгольца в цилиндрической геометрии для установления возможности создания конструкции системы колец с одной проточкой в канале. В качестве источников предложено выбирать функцию Грина задачи рассеяния минимальной объемлющей области. Данный алгоритм использован для расчетов параметров секции суперколлайдера в проведенном в ИПФ РАН эксперименте.

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

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

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

6. Предложено использовать дополнительный демпфирующий желоб не W только для предотвращения пробоя, но и как дополнительный резонатор на используемой основной временной частоте. Расположенная на достаточно малом радиусе, верхняя граница этого желоба сможет минимизировать полную микроволновую энергию в получившейся системе из двух резонаторов, не уменьшая амплитуды требуемой гармоники. В качестве конечного результата вычислительного эксперимента была сконструирована система связанных резонаторов различной конфигурации на одну и ту же частоту. По результатам расчетов в отделе электроники больших мощностей Института прикладной физики РАН изготовлена секция суперколлайдера и проведен "холодный" эксперимент по исследованию электродинамических характеристик полученной системы. Эксперимент подтвердил хорошие резонансные свойства исследуемой связанной системы резонаторов.

7. Проведено исследование возможности появления аномалии Вуда во внешней по отношению к системе колец области и вносимые ею изменения в

• структуру полей. При ее появлении во внешней части колец только у магнитного поля образуется стоячая волна большой амплитуды. У электрического поля этой составляющей нет.

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

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

• При этом в случае, когда зазоры узкого и широкого каналов dmin и dmax стремятся к нулю, имеет место сходимость численного решения задачи рассеяния к асимптотическому (см. Рис. 8.10).

9 Заключение

Основные научные результаты работы заключаются в следующем:

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

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

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

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

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

• Разработана методика построения схем повышенной точности, предложен экономичный аналога метода Бунемана прямого решения системы разностных уравнений. Метод требует (N\ogN)2 действий. Здесь N - число узлов по каждой из пространственных координат. Экономичность предложенного метода делает его перспективным для работы в непрямоугольных областях, где возможно применение альтернирующего метода Шварца.

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

• Исследована структура стационарного СВЧ разряда в пересекающихся волновых пучках с помощью разработанных методов . Проведено сравнение численных результатов с экспериментальными данными. Показано их хорошее качественное согласие.

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

• Разработанные сеточные методы использовались в качестве контроля для всех уравнений и систем глав 3-5 на начальной стадии рассматриваемых процессов.

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

• Доказана возможность использования усредненного по времени оператора правой части при решении задачи Коши для получения приближенного решения с точностью О (At3). Этот теоретический результат используется в главах 4 и 5 при построении приближенной операторной экспоненты.

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

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

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

• Проведен анализ возможности использования БПФ для различных типов граничных условий и исследованы возможности получения алгоритмов без насыщения.

• Проанализированы численные аналоги двух первых интегралов движения нестационарного уравнения Шредингера.

• Построено приближенное решение нестационарного уравнения Шредингера путем проведения последовательных операторных замен неизвестной функции. Замена проводится таким образом, чтобы в начальной точке производная по времени приближенного решения совпадала с производной точного решения, а в конечной точке интервала с соответствующей производной для той близкой к решению траектории, на которую выходит приближенное решение. Поэтому разность между точным и приближенным решениями в начальной и конечной точках временного шага обладает нулевой производной по времени. Сделанные оценки погрешности не только имеют порядок по времени 0(Д£3), но и являются лучшими из известных по пространственным переменным за счет аналитического вычисления с помощью БПФ операторных экспонент. Предложенный подход применен для нахождения приближенного решения систем уравнений, включающих нестационарное уравнение Шредингера.

• Исследована неустойчивость периодических решений уравнения Шредингера с кубичной нелинейностью, так называемых кноидальных волн, построена математическая схема возникновения "волн-убийц", аналогичная схеме развития так называемых аномалий Вуда. В линеаризованной задаче образуется ячейка Жордана, соответствующая нулевому собственному значению оператора, из присоединенных векторов которой и возникают аномальные решения.

• Продемонстрирована возможность успешного применения метода, основанного на схеме расщепления с использованием дискретного преобразования Фурье при нахождении решения нелинейного уравнения Шредингера гиперболически — параболического типа. Обнаружен новый тип самовоздействия волновых полей в нелинейных средах — множественное дробление волновых структур в нелинейной среде.

5. Предложен и реализован метод расчета уравнений Захарова, основанный на аналитическом вычислении операторных экспонент с помощью БПФ. Метод достаточно хорошо сбалансирован по временной и пространственным переменным, что невозможно сделать при использовании традиционных сеточных методов. В случае одной пространственной переменной метод протестирован на известном примере системы сталкивающихся солитонов.

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

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

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

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

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

Изучены процессы быстрой хаотизации вихревой части поля на стадии разрушения (схлопывания) каверны. Четко прослежены образование и рост кольцевых структур в образах Фурье у вихревой части поля и проанализирован механизм их образования. Это говорит о том, что для адекватного описания вихревой части поля необходимо вести его моделирование в образах Фурье. Более точно - при распараллеливании процесса вычислений требуется использовать локализованные тригонометрические базисы типа базисов Габора где функции wk(t) - оконные функции, хорошо локализованные в пространстве [94]. Применение этих базисов стало развиваться с появлением вейвлет - преобразований.

Перспективным является использование Быстрого Преобразования Ханкеля [228] .

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

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

6. Найдена в работах [217, 31, 221, 231] связь между дифференциалом дуги огибающей семейства характеристик уравнения Власова и так называемым "приведенным пространственным заряда" в уравнении Пуассона, что позволяет свести задачу описания стационарного пространственного заряда, создаваемого потоком электронов в скрещенных полях к граничному интегральному уравнению относительно приведенного пространственного заряда.

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

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

Изложенный подход к задаче является новым. Для установления существования и исследования устойчивости решения системы квазилинейных интегральных уравнений с дополнительными нелокальными граничными условиями адекватным задаче кажется аппарат теории всплесков (wavelets). Компактность носителя образующей всплесков и возможность построения такой системы всплесков, для которой дополнительные нелокальные граничные условия являлись бы естественной связью всплесков разных уровней, позволяют надеяться на построение конечномерной приближенной модели задачи с виртуальным катодом при небольшом числе моделирующих элементов (базисных всплесков). Постоянное число элементов, связанных дополнительным нелокальным условием, облегчает проведение предельных переходов. Это высказывание хорошо согласуется с утверждением (см. [196]), что плотность пространственного заряда, как и всплесков Добеши, обладает сплошным широким спектром.

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

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

• Рассмотрена задача о возмущении граничного интегрального уравнения для азимутально симметричных мод азимутально симметричных волноводов с плавной периодической гофрировкой путем гармонических возмущений его профиля. Найдено уравнение для возмущения плотности азимутального магнитного тока на граничной поверхности. Для исходного круглого волновода найдена явная связь между гармониками возмущенной плотности и возмущенного профиля волновода, а также выражение для возмущенного волнового числа.

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

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

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

Во-вторых, фазовые отклонения численного решения от асимптотического при c/min < А/8 составляют проценты.

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

Разработан и реализован алгоритм метода дискретных источников решения двумерного эллиптического уравнения типа уравнения Гельмгольца в цилиндрической геометрии для установления возможности создания конструкции системы колец с одной проточкой в канале. В качестве источников предложено выбирать функцию Грина задачи рассеяния минимальной объемлющей области. Данный алгоритм использован для расчетов параметров секции суперколлайдера в проведенном в ИПФ РАН эксперименте.

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

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

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

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

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

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

Список литературы диссертационного исследования доктор физико-математических наук Юнаковский, Алексей Дмитриевич, 2004 год

1. Агошков В.И., Дубовский П.Б., Шутяев В.П. Методы решения задач математической физики// М.: Физматлит. -2002. -320с.

2. Айдагулов Г.Р. Метод подвижной сетки для решения нестационарного уравнения Шредингера // Вычислительные методы и программирование. -2004. -Т.5. № 1. С. 22-34.

3. Алгазин С.Д. Численные алгоритмы без насыщения в классических задачах математической физики// М.: Научный мир. -2002. -156с.

4. Алексидзе М.А. Фундаментальные функции в приближенных решениях граничных задач// М.: Наука. -1991. 352с.

5. Альшина Е.А., Калиткин Н.Н., Панченко С.Д. Численное решение краевых задач в неограниченной области // Мат. моделирование. -2002. -Т.14. № 11. С. 10-22.

6. Анфиногентов В.Г., Храмов А.Е. Неавтономные колебания электронного потока с виртуальным катодом в плоском диодном промежутке // Изв.вузов "ПНД". -1997. -Т.5. № 6. С.61-75.

7. Анфиногентов В.Г., Храмов А.Е. Исследование колебаний в электронном потоке с виртуальным катодом в виркаторе и виртоде // Изв. вузов "ПНД". -1999. -Т.7. № 2,3. С. 33-55.

8. Арансон И,С., Горшков К.А., Рабинович М.И., Юнаковский А.Д. Стохастизация и расплывание волновых пакетов в осциллирующем потенциале.

9. Препринт № 256. Горький, ИПФ АН СССР. 1989. 20с.

10. Арнольд В.И. О матрицах, зависящих от параметров// Успехи мат. наук. -1971. Т. 26. С. 101-114.

11. Арсеньев А.А. Об особенностях продолжения и резонансных свойствах решения задачи рассеяния для уравнения Гельмгольца // ЖВМиМФ. -1972. -Т.12. № 1. С.112-138.

12. Арсеньев А.А. Глолбальное существование слабого решения системы уравнений Власова // ЖВМиМФ. -1975. -Т.15. № 1. С. 131-143.

13. Арсеньев А.А. О существовании резонансных полюсов и резонансов при рассеянии в случае краевых условий II и III рода // ЖВМиМФ. -1976. -Т.16. № 3. С.718-724.

14. Арсеньев А.А. Об асимптотике энергии: переданной почти периодическим источником колебаний открытому резонатору за большое время // Мат. сборник. -2004. -Т. 195. № 3. С. 3-14.

15. Архипов А.В., Богданов Л.Ю., Воскресенский С.В., Левчук С.А., Лукша О.И., Сомин-ский Г.Г. Исследование колебаний объемного заряда и формирования пространственных структур в электронном потоке с магнитным удержанием. Часть II // Известия вузов

16. Звездочкой помечены работы по теме диссертации.

17. ПНД" -1995. -Т. 3 №5. С. 35-54.

18. Ахманов С.А., Сухоруков А.П., Хохлов Р.В. Самофокусировка и дифракция света в нелинейной среде // УФН. -1967. -Т.93. № 1. С. 19- 70.

19. Бабенко К.И. Теоретические основы и конструирование численных алгоритмов задач математической физики// М.: Наука. -1979. -744с.

20. Бахвалов Н.С. Жидков Е.П. Каданцева Е.П. Рубин С.Б. Сердюкова С.Н. Численное решение модельной самосогласованной электродинамической задачи // ЖВМиМФ. 1979. -Т.19. № 6. С. 1228-1236.

21. Беллюстин С.И. ЖЭТФ. -1939. -Т.9. № 7. С.840.

22. Белов В.Е., Родыгин Л.В., Родыгина Л.С., Юнаковский А.Д. Явления на минимуме потенциала // Электроника сверхвысоких частот /тезисы докладов VI межвузовской конференции по электронике СВЧ. МРТИ, Минск. -1969. С.124. '

23. Белов В.Е., Родыгин Л.В., Юнаковский А.Д. К теории статического режима магнетрона с учетом начальных скоростей // Аннотации и тезисы докладов секции электроники Всесоюзной научн. сессии НТОРЭ и С им. А. С. Попова,. Москва. -1968. С.14.

24. Белов В.Е. Юнаковский А.Д. Задача дифракции акустических волн на упругом цилиндре в слоистом волноводе // Спектральш i еволюцшш задачь Тез! доповщей. Випуск 2 С1мферополь. -1993. С. 156-157.

25. Бесараб М.А. Кравченко В.Ф. Численно-аналитические методы решения задач дифракции и рассеяния электромагнитных волн на объектах сложной формы // Успехи современной радиоэлектроники. -2004. . № 4. С.4-71.

26. Богомолов Я.Л., Демехов А.Г., Трахтенгерц В.Ю., Шер Э.М., Юнаковский А.Д. О динамике квазилинейной релаксации быстрых частиц в плазме при наличии источников и стока частиц и волн// Препринт № 211. Горький, ИПФ АН СССР. 1988. 30с.

27. Богомолов Я.Л., Сергеев А.С., Юнаковский А.Д. Построение схем численного интегрирования для обыкновенных дифференциальных уравнений специального вида// Препринт №125. Горький, ИПФ АН СССР. 1985. 32с.

28. Богомолов Я.Л., Юнаковский А.Д. Метод построения явных разностных схем для обыкновенных дифференциальных уравнений второго порядка // X Всесоюзная школа "Теории прикладных проблем вычислительной математики и математической физики". Рига 1985г. С. 183.

29. Братман В.Л., Гинзбург Н.С., Петелин М.И. Лазеры на свободных электронах: перспективы продвижения классических электронных генераторов в коротковолновые диапазоны // Известия АН СССР, серия физическая. 1980. -Т.44. С. 1593-1602.

30. Бугаев С,П., Ильин В.П., Кошелев В.И. и др. Формирование сильноточных релятивистских электронных пучков для мощных генераторов и усилителей СВЧ.// Горький, ИПФ АН СССР. 1979. С. 5.

31. Бургейн Ж. Недавний прогрес в квазипериодических решеточных операторах Шредингера и Гамильтоновых дифференциальных уравнениях с частными производными // Успехи мат. наук. 2004. -Т.59. № 2. С.37-52.

32. Быков Н.М., Губанов В.П., Гунин А.В. и др. Релятивистские импульсно-периодические СВЧ-генераторы сантиметрового диапазона длин волн.// Релятивистская высокочастотная электроника, -вып.5/ Горький, ИПФ АН СССР. 1988. С.115.

33. Быковский В.А. О вычислительной сложности Дискретного Преобразования Фурье // ДАН. -2003. -Т.389. № 3. С. 293-298.

34. Вайнштейн Л,А. Электро-магнитные волны// М.: Радио и связь. 1988. 440с.

35. Веденяпин В.В. Граничные задачи для стационарного уравнения Власова // ДАН СССР. 1986. -Т.290. С. 777-780.

36. Вихарев А.Л., Гильденбург В.Б., Иванов О.А., Степанов А.И. СВЧ разряд в пересекающихся пучках электромагнитных волн // Физика плазмы. 1984. -Т.19. № 1. С. 165-168.

37. Вишик М.И., Люстерник Л.А. Решение некоторых задач о возмущении в случае матриц и самосопряженных и несамосопряженных дифференциальных уравнений // Успехи мат. наук. 1960. -Т. 15. № 1. С. 3-80.

38. Владимиров B.C. Уравнения математической физики// -М.: Наука,. 1976. -528с.

39. Власенко В.А., Лаппа Ю.М. Метод анализа погрешностей вычисления в алгоритмах Быстрого Преобразования Фурье // Радиотехника и электроника. 1986. -Т.31. № 7. С. 1348-1351.

40. Власов С.Н., Петрищев В.А., Таланов В.И. Усреднённое описание волновых пучков в линейных и нелинейных средах (метод моментов) // Изв. вузов Радиофизика. 1971. -Тю14. 9. С. 1353 -1364.

41. Власов С.Н., Таланов В,И. Распределённый волновой коллапс в модели нелинейного уравнения Шредингера// В кн. Нелинейные волны. Динамика и эволюция. /Отв.ред. акад. А.В.Гапонов-Грехов, д.ф.-м.н. М.И.Рабинович.// М.: Наука,. 1989. С. 218 - 227.

42. Габдулхаев Ю.Г. Прямые методы решения сингулярных интегральных уравнений первого рода. Казань: изд-во Казанского ун-та. 1994. 288с.

43. Гадылыпин P.P. Об амплитуде колебаний для резонатора Гельмгольца // ДАН СССР. 1990. -Т.310. № 5. С.1094-1097.

44. Гадылыпин P.P. О квазистационарном режиме резонатора Гельмгольца // Прикладная матем. и механика. 1993. -Т.57. № 5. С.52-59.

45. Гадылыпин P.P. Расщепление полюсов резонатора Гельмгольца // Изв. РАН Сер. матем. 1993. 57. № 5. С.44-74.

46. Гадылыпин P.P. Системы резонаторов // Изв. РАН Сер. матем. 2000. -Т.64. № 3. С.51-96.

47. Гадылыпин P.P. О рассеянии плоской волны на ловушке в критическом случае // Прикладная математика и механика. 2002. -Т.66. № 1. С.52-61.

48. Гилмор Р. Прикладная теория катастроф, кн. 1,2/1 М.: Мир,. 1984. -352с., 288с.

49. Гильденбург В.Б. Неравновесный высокочастотный разряд в полях электромагнитных волн// В кн. / Нелинейные волны. Распространение и взаимодействие. М.: Наука. 1981. С. 87-96.

50. Гильденбург В.Б., Жарова Н.А., Позднякова В.И. Нелинейная динамика волновых пучков и пакетов в ионизируемой среде // Изв. вузов. Радиофизика. 2003. -Т.46. № 5-6. С. 360-376.

51. Гинзбург Н.С., Кузнецов С.П., Федосеева Т.Н. Теория переходных прцессов в релятивистских JIOB // Изв.вузов. Радиофизика. 1978. -Т.21. № 7. С. 1037-1045.

52. Годунов С.К. Уравнения математической физики// М.: Наука. 1971. -365с.

53. Годунов С.К. Современные аспекты линейной алгебры// Новосибирск: Научная книга. 1997. -390с.

54. Головнин Б.А. Параллельные вычислительные системы// М.: Наука. 1980. -520с.

55. Горбачук В.И., Горбачук M.JI. Граничные задачи для дифференциально операторных уравнений// - Киев: Наук, думка. 1989. - 284с.

56. Горшков К.А., Островский Л.А., Папко В.В. Журн. Эксперим. и теорет. физики. 1976. -Т.71. С. 585.

57. Гохберг И.Ц., Крейн М.Г. Введение в теорию несамосопряженных операторов// М.: Наука. 1965. -448с.

58. Градштейн И.С., Рыжик И.М. Таблице интегралов, сумм, рядов и произведений// М.: ГИФМЛ. 1963. -1100с.

59. Григорьев Ю.Н., Вшивков В.А. Численные методы "частицы-в-ячейках"// Новосибирск: Наука. Сибирская издательская фирма РАН. 2000. -184с.

60. Громов Е.Н., Тютин В.В. Короткие солитоны огибающей, изотропные и анизотропные среды // Изв. вузов. Радиофизика. 2003. -Т.46. № 5-6. С. 429-442.

61. Гюнтер Н.М. Теория потенциала и ее применение к основным задачам математической физики// М.: ГИТТЛ. 1953. -416с.

62. Давыдов А.С. Теория твердого тела// М.: Наука. 1976.

63. Дегеярёв Л.М., Захаров В.Е., Рудаков Л.И. Два примера коллапса ленгмюровских волн // ЖЭТФ. 1975. -Т.66. № 1. С.115 126.

64. Дегтярёв Л.М., Копа-Овдиенко А.А. Скалярная модель ленгмюровского коллапса. // Физика плазмы. 1984. -Т. 10. № 1. С. 9 -20.

65. Дегтярёв JI.M., Крылов В.И. Метод численного решения задач динамики волновых полей с особенностями // ЖВМиМФ. 1977. -Т.17. № 6. С. 1523 1530.

66. Деммель Дж. Вычислительная линейная алгебра. Теория и приложения// -М.: Мир. 2001. -430с.

67. Джеффрис В., Джонсон Г.В. Операторное исчисление Фейнмана для семейств некомму-тирующих операторов: тензорные произведения, упорядоченные носители и выпутывание экспотенциального множителя // Мат. заметки. 2001. -Т.70. № 6. С. 815-838.

68. Ингрид Добеши. Всплески и другие методы локализации в фазовом пространстве// В кн. Международный конгресс математиков в Цюрихе 1994■ М.: Мир. 1999. С. 84-108 .

69. Додд Р., Эйлвен Дж., Гиббон Дж., Моррис X. Солитоны и нелинейные волновые уравнения// М.:Мир,. 1988. -696с.

70. Дубинов А.Е. Особенности динамики электронов в виркаторе с магнитной пробкой // Радиотехника и электроника. 2000. -Т.45. № 7. С. 875-879.

71. Дубинов А.Е., Ефимова И.А. Колебания сжатого состояния электронных пучков в виркаторе на встречных потоках // Электромагнитные волны и электронные системы. 2003. № 11 -12. С. 55-57.

72. Дубинов А.Е., Сайков С.К., Селемир В.Д. К динамике метастабильных электронов в виртуальном катоде // Изв.вузов.Радиофизика. 2002. -Т.45. № 4. С. 349-353.

73. Жаблон К. Симон Ж.К. Применение ЭВМ для численного моделирования в физике// -М.: Наука,. 1983. -238с.

74. Жарова Н.А., Литвак А.Г., Миронов В.А. Самовоздействие сверхкороткого лазерного импульса в среде с нормальной дисперсией групповой скорости // Изв. вузов. Радиофизика. 2003. -Т.46. № 5-6. С. 331-341.

75. Жиженкова Е.Ф. Жислин Г.М. О существовании минимума некоторых квадратичных функционалов в неограниченной области // Труды Московского матеметического общества. 1960. -Т.9. С. 121 129.

76. Захаров В.Е. Устойчивость периодических волн конечной амплитуды на поверхности глубокой жидкости // ПМТФ-. 1968. -Т.6. № 3. С. 86 -94.

77. Захаров В.Е. Коллапс плазменных волн // ЖЭТФ. 1972. -Т.62. № 5, С. 1745 1759.

78. Захаров В-Е., Кузнецов Е.А. Квазиклассическая теория трёхмерного волнового коллапса. // ЖЭТФ. 1986. -Т.91. № 4(10). С. 1310 1324.

79. Захаров В.Е., Рубенчик A.M. Журн. эксперим. и meopem. физики. 1973. -Т.65. С. 997.

80. Земанян А.Г. Интегральные преобразования обобщенных функций// М.: Наука, Гл. ред. физ.мат. лит. 1974. -400с.

81. Ильинский А.С., Кравцов В.В., Свешников А.Г. Математические модели электродинамики/ / М.: Высшая школа. 1991. - 224с.

82. Инфельд Э. Письма в ЖЭТФ. 1980. -Т.32. С. 97.

83. Кадомцев Б.Б. Колллективные явления в плазме М.:Наука. 1988. -304с.

84. Калинин Ю.А., Кузнецов Н.Н., Украинская Т.Н. Исследование широкополостных шумо-подобных колебаний в интенсивных пучках заряженных частиц в режиме образования виртуального катода // Изв. вузов "ПНД'Ш2. -Т.10. № 5. С.32-35.

85. Камоцкий И.В., Назаров С.А. Аномалии Вуда и поверхностные волны в задаче рассеяния на периодической границе. 1. // Мат. сборник. 1999. -Т.190. № 1. 109-138.

86. Камоцкий И.В., Назаров С.А. Экспонентно затухающие решения задачи дифракции на жесткой периодической границе // Мат.заметки. 2003. -Т.73. № 1. С.138-140.

87. Канторович J1.B. Функциональный анализ и прикладная математика // УМН. 1948. -Т.З. № 6. С.89-185.

88. Кашкин В.Б., Носков М.В., Осипов Н.Н. Вариант Дискретного Преобразования Фурье на параллелепипедальных сетках // ЖВМиМФ. 2001. -Т.41. № 3. С. 355-359.

89. Ким Ы., Ковалев Н.Ф., Фильченков С.Е. Коаксиальный волновод с прорезями на внутреннем проводнике Радиотехника и электроника 2003. -Т.48 №4 С. 1-9.

90. Кингсеп А.С. Сильная ленгмюровская турбулентность и турбулентный нагрев плазмы// В кн. Итоги науки и техники, серия "физика плазмы." /Под ред. Шафранова В.Д. вып. 4 М.: ВИНИТИ. 1983. С. 48-113.

91. Кингсеп А.С. Введение в нелинейную физику плазмы// М.: МЗ Пресс. 2004. -264с.

92. Кириленко А.А., Сенкевич C.JI. Сравнение эффективности четырех методов решения волноводных задач // Радиотехника и электроника. 1984. -Т.29. № 6. С. 1089-1097.

93. Кирпичев С.Б., Поляков П.А. Постановка начальной задачи для системы релятивистских заряженных частиц // Электромагнитные волны и электронные системы. 2004. -Т.9. № 6. С. 8-13.

94. Ковалев Н.Ф., Петелин М.И., Райзер М.Д. и др. Приборы типа О, основанные на индуцированном черенковском и переходном излучениях релятивистских электронов.// В кн. Релятивистская высокочастотная электроника Горький ИПФ АН СССР. 1979. С. 76.

95. Колоколов А.А., Суков А.Н. О векторной теории самофокусировки. // ПМТФ. 1977. № 6. С. 3 5.

96. Копосова Е.В., Петелин М.И. Изв. вузов. Радиофизика. 1989. Т.32. №2 С. 1178.

97. Костюченко А.Г., Саргсян И.С. Распределение собственных значений// М.: Наука. 1979. -400с.

98. Кудряшов В.Р., Михайловский А.Б. ЖЭТФ. 1986. -Т.90. С. 1656.

99. Кузнецов Е.А. Коллапс электромагнитных волн в плазме // ЖЭТФ. 1974. -Т.66. № 6. С. 2037 2047.

100. Кузнецов Е.А. Интегральные критерии волновых коллапсов // Изв. вузов. Радиофизика. 2003. -Т.46. № 5-6. С. 342-359.

101. Купрадзе В.Д. О приближенном методе решения задач математической физики // Успехи мат. наук. 1967. -Т.22. № 2. С.58-109.

102. Кураев А.А. Мощные приборы СВЧ. Методы анализа и оптимизации// М.: Радио и связь. 1986. 208с.

103. Курант Р. Уравнения с частными производными// М.: Мир. 1964. -822с.

104. Кучмент П.А. Теория Флоке для дифференциальных уравнений в частных производных // УМН. 1982. -Т.37. № 4. С. 3 52.

105. Кюркчан А.Г., Минаев С.А. Решение задач дифракции с использованием техники вейвлетов // Радиотехника и электроника. 2003. -Т.48. № 5. С.552-557.

106. Кюркчан А.Г., Минаев С.А., Соловейчик A.JI. Модификация метода дискретных источников на основе априорной информации об особенностях дифракционного поля // радиотехника и электроника. 2001. -Т.46. № 6. С.666-672.

107. Кюркчан А.Г., Суков А.И., Клеев А.Г. Особенности волновых полей и численные методы решения краевых задач для уравнения Гельмгольца. // Успехи современной радиоэлектроники. Зарубежная электроника. 2000. № 5. стр . 14-33.

108. Ландау Л.Д., Лифшиц Е.М. Электродинамика сплошных сред// М.: Наука. 1982. -620с.

109. Ланцош К. Практические методы прикладного анализа// М,: ГИФМЛ. 1961. -524с.

110. Леонтович М.А Изв. АН СССР, Сер. физич. 1944. -Т.8. С. 16.

111. Леус В.А. О временной сложности параллельной реализации численного преобразования Фурье // Мат. моделирование. 1995. -Т.7. № 5. С. 99-110.

112. Сергеев А.С., Литвак А.Г. Сб. Высокочастотный нагрев плазмы// Горький, ИПФ АН СССР. 1983.

113. Литвак А.Г., Таланов В.И. Изв. высш. уч. зав., сер. Радиофизика. 1967. -Т.9. № 10. С. 539.

114. Лифанов И.К. Метод сингулярных интегральных уравнений и численный эксперимент// -М.: ТОО "Янус". 1995. -520с.

115. Малкин В.М. ЖЭТФ,. 1984. -Т.87. С. 433.

116. Вычислительные процессы и системы. Под ред. Марчука Г.И.// № 2. М.: Наука. 1985.

117. Марчук Г.И. Методы расщепления// -М.: Наука,. 1988. -264с.

118. Марчук Г.И. Методы вычислительной математики// -М.: Наука. 1989. -608с.

119. Марчук Г.И., Шайдуров В.В. Повышение точности решений разностных схем// М.: Наука. 1979. -320с.

120. Месяц Г.А. Импульсная энергетика и электроника// М.: Наука. 2004. -704с.

121. Миллер Р. Введение в физику сильноточных пучков заряженных частиц// М.: Мир. 1984. -432с.

122. Молотков И.А. Аналитические методы в теории нелинейных волн// М.: Физматлит. 2003. -208с.

123. Морозов А.И., Соловьев Л.С. Движение заряженных частиц в электромагнитных полях// В кн. Вопросы теории плазмы, вып. 2 /Под ред. Леонтовича М.А. М.: Атомиздат. 1963. с. 177.

124. Назайкинский В., Стернин В., Шаталов В. Методы некоммутативного анализа// М.: Техносфера. 2002. -234с.

125. Новиков И.Я., Стечкин С.Б. Основы теории всплесков, // УМН. 1998. -Т.53. № 6. С. 55-128.

126. Овсянников Д.А., Егоров Н.В. Математическое моделирование систем формирования электронных и ионных пучков// Издательство Санкт-Петербургского университета. 1998. -276с.

127. Окуно Ю., Икуно X. Метод Ясууры для решения граничных задач уравнения Гельмгольца в однородных средах // Радиотехника и электроника. 1992. -Т.37. № 10. С. 1744-1763.

128. Павленко В.Л., Петвиашвили В.И. Физика плазмы. 1982. -Т.8. С. 206.

129. Парангишвили И.В., Виленкин С.Я., Медведев И.Д. Параллельные вычислительные системы с общим управлением// М.: Энергоиздат. 1983. 312с.

130. Пелиновский Е.Н., Слюняев А.В., Талипова Т.Г., Хариф К. Нелинейное параболическое уравнение и экстремальные волны на морской поверхности // Изв,вузов.Радиофизика. 2003. -Т.46. № 7. С.499-512.

131. Рид М., Саймон Б. Методы современной математической физики.: Т.2. Гармонический анализ. Самосопряженность.// М.:Мир. 1987. -400с.

132. Рид М., Саймон Б. Методы современной математической физики.: Т.2. Анализ операторов.// М.:Мир. 1982. -430с.

133. Рихтмайер Р., Мортон К. Разностные методы решения краевых задач// М.: Мир. 1972. -420с.

134. Рошаль А.С. Моделирование заряженных пучков// М.: Атомиздат. 1979. -234с.

135. Рудаков Л.И. Торможение электронных пучков в плазме с высоким уровнем лэнгмюров-ской турбулентности // ДАН. 1973. -Т.207. № 4. 821 824.

136. Савостьянов Д.В., Тыртышников Е.Е. О случае алгебраической эквивалентности метода коллокаций и метода Галеркина // ЖВМиМФ. 2004. -Т.44. № 4. С.686-693.

137. Самко С.Г., Килбас А.А., Маричев О.И. Интегралы и их производные дробного порядка// Минск: Наука и техника. 1987. 688с.

138. Самарский А.А. Теория разностных схем// М.: Наука. 1989. -616с.

139. Самарский А.А., Гулин А.Н. Устойчивость разностных схем// -М.: Наука. 1973. -416с.

140. Сидоров Н.А., Синицын А.В. Стационарная система Власова Максвелла в ограниченных областях// В кн. Нелинейный анализ и нелинейные дифференциальные уравнения. -М.: ФИЗМАТЛИТ. 2003. С. 50-88.

141. Сигалов А.Г. Уравнение непрерывного спуска для квадратичной формы полуограниченного оператора на единичной сфере. // ЖВМиМФ. 1971. -Т.П. № 4. С. 809 828.

142. Скубачевский А.Л. Эллиптические задачи с нелокальными условиями вблизи границы// Мат. сб. 1986. -Т.181 № 2. С.179-202.

143. Скубачевский А.Л. Модельные нелокальные задачи для эллиптических уравнений в двугранных углах // Дифференциальные уравнения. 1990. -Т.26. № 1. С.120-131.

144. Сыровой В.А. Геометризованная теория пространственных электростатических электронных пучков // Радиотехника и электроника. 2002. -Т.47. № 3. С. 372-382.

145. Таланов В.И. О самофокусировке волновых пучков в нелинейных средах // Письма ЖЭТФ-. 1965,. -Т.2. № 5. С. 218 -222.

146. Тарасов Р.П. Счетная устойчивость граничных уравнений I рода в задачах дифракции вблизи резонансов. // ЖВМиМФ. 1999. -Т.39. № 6. С. 943-969.

147. Теория магнетрона (по Бриллюэну). Сборник переводов -М.: Сов. радио. 1946.

148. Тихонов А.Н., Самарский А.А. Уравнения математической физики// М.: ГИТТЛ. 1953. 68с.

149. Трубецков Д.И., Храмов А.Е. Лекции по высокочастотной электронике для физиков т. /.// М.: ФИЗМАТЛИТ. 2003. -496с.

150. Уилкинсон Дж. X. Алгебраическая проблема собственных значений// М.: Наука. 1970. -564с.

151. Финк Л.М. Сигналы, помехи, ошибки.// М.: Радио и связь. 1984. -256с.205. под ред. проф. Я.Н.Фельда. Справочник по волноводам, пер. с англ.// -М.: Советское радио. 1952. -432с.

152. СМБ Функциональный анализ. Под ред. С.Г.Крейна// -М: Наука. 1964. -424с.

153. Фурсиков А.В. Оптимальное управление распределенными системами. Теория и прило

154. Новосибирск: Научн. кн. 1999 . -352с.

155. Хемминг Р.И. Численные методы// М.: Наука. 1968. -400с.

156. Хокни Р., Джессхоуп К. Параллельные ЭВМ// М.: Радио и связь. 1986. -392с.

157. Хокни Р., Иствуд Дж. Численное моделирование методом частиц// М.: Мир. 1987. -640с.

158. Холодниок М., Клич А., Кубичек М., Марек М. Методы анализа нелинейных динамических моделей// М.:Мир. 1991. 365с.

159. Шварц JI. Комплексные многообразия. Эллиптические уравнения.// М.: Мир. 1964. -212с.

160. Юнаковский А.Д. Математическое моделирование взаимодействия электромагнитного излучения с плазмой: Дис. . канд. физ.-мат. наук: 05.13.16 /ИПФ АН СССР. Горький, 1988,- 214с.

161. Юнаковский А.Д. Математическое моделирование физических процессов /вычислительные аспекты/ Методические указания для студентов специальностей 0645, 2001. Часть I Нижегородский политехнич. ин-т. Нижн. Новгород,. -1990,

162. Яненко Н.Н. Метод дробных шагов решения многомерных задач математической физики// Наука. Сибирское отделение. Новосибирск. -1967. -196с.

163. Abramyan L.A., Yunakovsky A.D. Application of Petviashvili's Method for Determination of Three-dimentional Langmuiere's Multipoles // Spectral and evolutional problems. Vol.4, Simferopol. 1996. Crimea, Ukraine. 1996. C.153

164. Ahner J.F. A scattering trinity: the reproducing kernel, null field equations and modified Green's Function. // Quart. J.Mech. and Appl. Math. 1986. -V.39. № 1. p.151-162.

165. Akhiezer N.I., Glazman I.M. Theory of Linear Operators in Hilbert Space// Frederick Ungar Publishing Co., New York. Vol.1 1961., Vol.11 1963.

166. Al-Abawi H., Hopf F.A., Moore G.T., Scully M.O. Lasers: Laser Letqargy and Coherence Brightening // Optics Communications. 1979. -V.30. № 2. p. 231-235.

167. Batt J., Berestycki H., Deyond P. Parthame В/. Sjme Families of Solutions of the Vlasov -Poisson System // Arch. Rational Mech. Anal. 1977. -V.104. p.79-103.

168. Batt J. Fabian K. Stationary Solutios of the Relativistic Vlasov- Maxwell System of Plasma Physics // Chin. Ann. Math, (ser В.) 1993. -V.14. p.253-278.

169. Beale J.T. Scattering Frequencies of Resonators // Comm. Pure and Appl. Math. 1973. -V.26. № 4. p.549-563.

170. Luc Berge. Dynamical stability analysis of strong/weak wave collapses // Journal of Mathematical Physics. 1994. -V.35. № 11. p.5765 5781.

171. Beylkin G., Coifman R., Rokhlin V. Fast wavelet transforms and numerical algorithms. // Comm. Pure Appl. Math. 1991. -V.44. p.141-183.

172. Beylkin G. Keiser J.M. On the Adaptive Numerical Solution of Nonlinear Partial Gifferential Equations in Wavelet Bases // Journ. of Сотр. Phys. 1997. -V.133. p.233-259.

173. Chiao R.Y., Garmire E., Townes C.H. Self-Trapping of Optical Beams // Phys.Rev.Lett. 1964. -V.13. p.469 472.

174. Ciment Melvyn, Leventhal Stephen H., Weinberg Bernard C. The Operator Compact Implicit method for parabolic equations. // Journ. of Сотр. Phys. 1978. -V.28. № 1. p.135 166.

175. Daubechies I. Orthonormal bases of compactly supported wavelets, // Comm. Pure Appl. Math. 1988. -V.41. p.909-996.

176. Daubechies I. The wavelet transform, time-frequency lokalisation and signal analysis, // IEEE Trans. Inform Theory. 1990. -V.36. p.961-1005.

177. Deacon D.A.G,. Elias L.R., Madey J.M.J., Ramian G.J., Schwettmah H.A., Smith T.I. First Operation of Free Electron Lasers // Phys. Rev. Lett. 1977. -V.38. p. 892-894.

178. DiPerna R.J., Lions P.L. Global Weak Solutions of Vlasov Maxwell System // Comm. Pure and Appl. Math. 1989. V. 42. №6. p. 729-757.

179. Dyachenko A.I., Pushkarev A.N., Shvets V.E., Zakharov V.E. Computer simulation of Langmuir collapse // Physica D. 1991. -V.52. p.78 102.

180. Fadeev D.K., Fadeeva V.N. Computational Methods of Linear Algebra. Freeman, San Francisco, California. 1963.

181. Feit M.D., Fleck Y.A., Steiger. Solution of the Schrodinger Equation by a Spectral Method. // J. Сотр. Phys. 1982. -V.47. № 4. p.412 433.

182. Feynman R. An operator calculus having applications in quantum electrodynamics // Phys. Rev. 1951. -V.84. p.108-128.

183. Hardin R.H., Tappert F.D. Application of the split-step Fourier method tj the numerical solution of nonlinear and variable coefficient wave equation // SI AM Rev. Chronicle. 1973. -V.15. p.423.

184. Infeld E, Rowlands G. Z. Phys. B. 1980. -V.37. p.277.

185. Infeld E., Ziemkiewicz J. Acta Phys. Polonica. 1981. -V.59. p.255.

186. Keller J.B. Singularities and Rayleigh's hypothesis for diffraction gratings // J. Opt. Soc. Am. A. 2000. -V.17. № 3. p. 456-457.

187. Kosmatov N.E., Shvets V.E., Zakharov V.E. Computer simulation of wave collapses in the nonlinear Schrodinger equation // Physica D. 1991. -V.52. № 1. p.16 35.

188. Kovalev N.F., Nechaev V.E., Petelin M.I., Zaitsev N.I. A Scenario for Output Pulse Shorteningm Microwave Generators Driven by Relativistic Electron Beams // IEEE Trans, on Plasma Science. 1998. -V.26. № 3. pp. 246-251.

189. Lonngren K.E., Pecseli H.L., Juui.Rasmussen J., Thomsen K. Transient effects of nonlinear wave propagation in magnetized plasmas. // Phys. Scr. 1982, T2: 2: Int. Conf. Plasma Phys., Goteborg, June 9 15. 1982. pp.541 - 545.

190. Lonngren K.E., Pecseli H.L., Juui.Rasmussen J. On the transient effects of nonlinear wave propagation in magnetized plasmas. // Ibid, pp.546 547.

191. Mallat S.G. Multifrequency channel decomposition of images and wavelets models. // IEEE Trans. Acoust. Speech Signal processes. 1989. -V.37. № 12. pp. 2091-2110.

192. Mallat S.G., Hwang W.L. Singularity detection and processing with wavelets // IEEE Trans, of Information Theory. 1991. -V.31. № 3. .

193. Marcuvitz N. Waveguide. Handbook./j Dover Publications Inc. 180 Varick Street New York, N.Y. 10115. -951. 428p.

194. Mesyats G.A., Proskurovsky D.I. Pulsed Electrical Discharge in Vacuum// Springer, Heidelberg. 1989.

195. Рецензии. Meyer Y. Ondelettes et operateurs // Алгебра и анализ. 1991. -Т.З. № 2. С. 253-264.

196. Meyer Y. Ondelettes sur rintervalle. // Rev.Mat. Iberoamericana. 1992, -V.7. pp.115-133.

197. Meyer Yves. Wavelets Algorithms & Applications// SIAM Philadelphia. 1993. 134p.

198. Myra J.R., Liu C.S. Phys. Fluids. 1980. -V.23. p.2258.

199. Newton P.K. Wave Interaction in the Singular Zakharov Sistem // J. Math. Phys. 1991. -V.32. № 2. .

200. Payne G.L., Nicholson D.R., Downie R.M. Numerical solution of the Zakharov Equations. // J. Сотр. Phys. 1983. -V.50. № 1. pp.482 498.

201. Petelin M.I. High power sources of coherent microwave radiation (gyrotrons, free electron lasers) and their applications // Annales des Journees Maxwell, Bordeaux. 1995. pp.277-286.

202. Petelin M.I. Quasi-Optical Collider Concept.// American Institute of Physics. 2002. pp. 459-468.

203. Petelin M.I., Suvorov E.V., Kovalev N.F., Filchenkov E.C., Smirnov A.I. Quasi optical Diffraction Gril for Exitation of Lower - hybrid Waves in Tokamaks// Plasma Phys. Control. Fusion 1996/, -V.38 pp. 593-610.

204. Petit R., Cadilhac M. Sur la difraction d'une onde plane par un reseau infiniment conducteur // C.R Acad.Sci, Paris, Ser. B. 1966. -V.262. pp. 468-471.

205. Richtmyer Robert D. Principles of Advanced Mathematical Physics. Vol. 1// Springer-Verlag New York Inc. 1978.

206. Ruth R. D. et. al. The Next Linear Collider Test Accelerator // Proc. of the IEEE Particle Accelerator Conference, Washington DC, May 1993,. Washington DC. 1993. pp.543-545.

207. Sack Ch., Schamel H. SUNION An Algorithm for One - Dimentional Laser Plasma Interaction. // J. Сотр. Phys. 1984. -V.53. № 3. pp.395 - 428.

208. Sami G., Tantawi et. al. The Next Linear Collider Test Accelerator's RF pulse Compression and Transmission Systems // Advanced Accelerator Concepts workshop, Lake Tahoe, California, October 12-18, 1996. 398. AIP proceedings. 1996. pp.229-317.

209. Schauder J.Zur. Theorie stetiger Abbildungen in Funktionalraumen // Math. Zeitschr. 1927. -V.26. pp.47-65.

210. Shwinger J. Discontinuites in waveguardes// Gordon Breach Science publischers, N -Y., London, Paris. 1969.

211. Simon B. Schrodinger operators in the twentieth century // J. Math. Phys. 2000. -V.41. № 6. pp.3523-3555.

212. Skubachevskii A.L. Nonlokal Elliptic Problems and Multidimentional Diffusion Prjcesses // Russian Journal of Mathematical Physics. 1995. -V.3. pp.327-360.

213. Steinberg J. Numerical Treatment of Singularities in Elliptic Problems. // J. Inst. Math. Applies. 1980. -V.26. pp.349-368.

214. Swartz В.К. The Construction of Finite Difference Analjgs of Some Finite Elements. Mathematical Aspects of Finite Elements in Partial Differential Equations, ed. by G.De Boor. Academic Press, New York. 1974. pp.412 433.

215. Tarakanov V.P. User's manual for code Karat. Bercley Research Associate Inc., Springfield, VA. 1992.

216. Weitzner Harold Longitudinal plasma oscillations // Proceedings of a Symposium in Appl. Math, of AMS, New York City, April 13-15 1965. ed.by Harold Grad. pp. 127-151.

217. Weitzner Harold. Initial Value-Boundary Value Problem for Bernstein Modes // Comm. Pure and Appl. Math. 1969. XXII. № 4. pp. 415-455.

218. Wilson P.V. Application of High-Power Microwave Sources to TeV Linear Colliders. 229-317. Artech House, Boston к London. 1994.

219. Wilson P.V. В кн: Proceedings ofthe ITP Conferenceon Future High EnergyColliders. University of California, Santa Barbara. 1996.

220. Yunakovsky A.D. Mathematical Simulation of Pitch-angular Diffusion in Raduiative Fields // Спектральные и эволюционные задачи вып. 3 Симферополь -. 1993. Крым, Украина. 1993. С. 104-106.

221. Yunakovsky A.D. Approximate Evaluation of Operator Exponential by Operator Compact Implicit Method // Spectral and evolutional problems. Vol.5, Simferopol. 1996. Simferopol, Crimea, Ukraine. 1996. C. 202-203.

222. Yunakovsky A.D. Simulation of the electromagnetic field collapse dynamics. // Ibid. C. 143-145.

223. Yunakovsky A.D. About attractors of trajectories of continuons descent for a square-law form of semibounded operator on a unit sphere // Spectral and Evolutional Problems.- vol.8, Symferopol. 1998. pp.166-170.

224. Yunakovsky A.D. Attractors of trajectories of a continuous descent for a square-law form of semibounded operator on a unit sphere // Nonlinear boundary value problems. 1998. № 8. pp.243-252.

225. Yunakovsky A.D. Wavelets and boundary value problems, j j Nonlinear boundary value problems. Book of abstracts. International conference Dedicated to J.P.Schauder, Lviv, August 23-29 1999. p.218.