Вычислительная математика, численный анализ

Является периодическим научным изданием, отражающим тематику важнейших направлений теоретических исследований по математике и механике в МГУ имени М.В.Ломоносова. На его страницах печатаются оригинальные статьи, посвященные конкретным научным вопросам по всем основным направлениям теоретических и прикладных исследований.

Тематика журнала включает в себя: Процессы горения и взрыва. Математическое моделирование физико-химических процессов. Кластеры, кластерные системы и материалы. Межфазные слои и процессы взаимодействия в них. Квантово-химические расчеты. Нелинейные кинетические явления. Наноэлектронные приборы и устройства. Журнал включен в Реферативный журнал и Базы данных ВИНИТИ РАН.

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

Рассмотрены основные возможности математической системы MathCAD и применение ее для решения задач математического моделирования. Предназначено для бакалавров направления 09.03.02 «Информационные системы и технологии», изучающих дисциплину «Моделирование физических процессов», и магистров направления 09.04.02 «Информационные системы и технологии», изучающих дисциплину «Моделирование объектов, явлений, процессов и систем».

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

Является периодическим научным изданием, отражающим тематику важнейших направлений теоретических исследований по математике и механике в МГУ имени М.В.Ломоносова. На его страницах печатаются оригинальные статьи, посвященные конкретным научным вопросам по всем основным направлениям теоретических и прикладных исследований.

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

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

В курсе лекций рассмотрены вопросы математического моделирования связанных физико-механических задач и применения методов конечных элементов для их решения. Особое внимание уделено связанным задачам электроупругости, термоупругости, пороупругости и задачам взаимодействия деформируемых тел с акустическими средами. Акцентируются математические особенности рассматриваемых задач в их классических и обобщенных постановках. Теоретическая часть пособия подкреплена лабораторными работами по решению модельных задач пьезоэлектричества, термоупругости и задач о взаимодействии твердотельных тел с акустическими средами с использованием программного комплекса ANSYS версии 11.0 и выше. Пособие является переработанной версией в авторском переводе соответствующего англоязычного пособия: Nasedkin A. V., Nasedkina A. A. Finite element modeling of coupled problems: textbook / Rostov-on-Don: publishing house of Southern Federal University, 2015.

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

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

В кватернионной постановке рассмотрена задача оптимальной переориентации плоскости орбиты космического аппарата (КА). Управление (ускорение от вектора реактивной тяги, ортогональной плоскости орбиты) ограничено по модулю. Необходимо минимизировать затраты энергии на процесс переориентации плоскости орбиты КА. Рассмотрен актуальный частный случай задачи, когда орбита КА является круговой, а управление принимает постоянные значения на отдельных участках активного движения КА. Построен оригинальный генетический алгоритм нахождения траекторий оптимальных перелётов КА. При применении этого способа не требуется какая-либо информация о неизвестных начальных значениях сопряжённых переменных. Приведены примеры численного решения задачи для случая, когда отличие между начальной и конечной ориентациями орбиты КА составляет единицы градусов в угловой мере. При этом конечная ориентация плоскости орбиты КА соответствует ориентации плоскости орбиты спутников отечественной орбитальной группировки ГЛОНАСС.

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

Рассматривается задача дифракции стороннего электромагнитного поля на локально неоднородном теле, помещенном в свободном пространстве. Поставленная задача сводится к объемному сингулярному интегральному уравнению. Решение задачи производится параллельно численным методом Галеркина и численным методом коллокации. В связи с большой емкостью решение задачи численным методом Галеркина при различных параметрах было реализовано с использованием двух программных продуктов для суперкомпьютерных вычислительных комплексов: реализации MPI и программной системы x-com. Исследованы особенности выполнения задачи на суперкомпьютерном комплексе.

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

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

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

Рассматривается метод трехуровневого распараллеливания, применяемый в методике ТИМ-2D. На первом уровне осуществляется распараллеливание счета по математическим областям. На втором – распараллеливается счет внутри области по параобластям. На первых двух уровнях используется модель распределенной памяти и интерфейс передачи сообщений MPI. На третьем – распараллеливаются итерации счетных циклов в модели общей памяти с использованием OpenMP. Приводятся замеры эффективности распараллеливания на ряде тестовых расчетов

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

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

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

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

При моделировании кинетики химических реакций, расчете электронных схем и электрических сетей и других важных приложений возникает необходимость решения задачи Коши для жестких систем обыкновенных дифференциальных уравнений. Для решения таких задач применяются L-устойчивые численные схемы. В таких методах при большой размерности системы дифференциальных уравнений основные вычислительные затраты приходятся на декомпозицию матрицы Якоби. Сокращения затрат достигают замораживанием матрицы Якоби, т. е. применением одной матрицы на нескольких шагах интегрирования. Дополнительного сокращения затрат добиваются за счет применения алгоритмов интегрирования на неоднородных схемах. В состав таких алгоритмов включаются явные и L-устойчивые методы. Эти алгоритмы сами распознают, является задача жесткой или нет. Эффективная численная схема выбирается на каждом шаге по критерию устойчивости. Здесь разработан неоднородный алгоритм интегрирования на основе L-устойчивого и явных двухстадийных методов. Построено неравенство для контроля устойчивости схемы Рунге-Кутта второго порядка точности. На основе стадий этого метода предложена численная формула первого порядка с расширенным до 8 интервалом устойчивости. На основе L-устойчивой (2, 2) -схемы и численных формул типа Рунге-Кутта первого и второго порядков точности разработан алгоритм переменной структуры, в котором эффективный метод выбирается на каждом шаге по критерию устойчивости. При расчетах по L-устойчивому методу допускается замораживание матрицы Якоби, которая может вычисляться как аналитически, так и численно. Алгоритм предназначен для решения как жестких, так и нежестких задач. Приведены результаты расчетов, подтверждающие эффективность построенного алгоритма.

Итерационный попеременно-треугольный кососимметричный метод (ПТКМ) используется для решения СЛАУ, полученной при аппроксимации центрально-разностной схемой первой краевой задачи конвекции–диффузии–реакции и использовании стандартного перебора узлов сеточной области. Для знакопеременного коэффициента реакции даны достаточные условия неотрицатльной определенности матрицы СЛАУ, полученной в результате такой аппроксимации. Это свойство гарантирует сходимость достаточно широкого класса итерационных методов, в частности ПТКМ. На тестовых задачах проверено соответствие теории вычислительному эксперименту и дано сравнение ПТКМ и SSOR.

Издается с 1978г. В данной серии публикуются оригинальные работы, обзоры и краткие сообщения по следующим научным направлениям: математическое моделирование физических процессов и свойств веществ, численные и аналитические методы решения прикладных задач математической физики и механики сплошной среды; вычислительная математика и применение математических методов и электронно-вычислительной техники в научных исследованиях; вопросы программирования; вопросы структуры алгоритмов и программ для современных ЭВМ; вопросы создания вычислительных комплексов и сетей ЭВМ. Главный редактор - д-р ф.-м. наук Р.М.Шагалиев

Актуальность и цель. В статье К. И. Бабенко «О некоторых задачах теории приближений и численного анализа»1 среди ряда важных проблем вычислительной математики были сформулированы две проблемы: 1) вычисление поперечников Колмогорова и Бабенко на классе Q(Ω,M)(класс Q(Ω,M)со-

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

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

Гиперсингулярные интегралы в настоящее время находят все большие области применения – аэродинамика, теория упругости, электродинамика и геофизика. При этом их вычисление в аналитическом виде возможно лишь в весьма частных случаях. Поэтому приближенные методы вычисления гиперсингулярных интегралов являются актуальной задачей вычислительной математики. Этой задаче посвящено много работ. В частности, И. В. Бойковым и Ю. Ф. Захаровой опубликованы циклы работ по построению оптимальных методов вычисления гиперсингулярных интегралов. В 1975 г. в докладах АН СССР (т. 221, № 1) опубликована статья К. И. Бабенко, в которой он сообщил об открытии им принципиально новых – ненасыщаемых численных методах. Отличительной особенностью последних является способность автоматически подстраиваться под классы корректности решений рассматриваемых задач. Анализ известных квадратурных и кубатурных формул вычисления гиперсингулярных интегралов показал, что они являются насыщаемыми. Поэтому является актуальной задача построения ненасыщаемых алгоритмов вычисления гиперсингулярных и полигиперсингулярных интегралов. Этой задаче посвящена данная работа.

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

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

Издание рассматривает методики, предназначенные для решения задач механики сплошной среды на неструктурированных многоугольных и многогранных сетках. Это методики ДМК, ТИМ и МЕДУЗА, применяемые в Институте теоретической и математической физики РФЯЦ-ВНИИЭФ.

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

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

Издается с 1978г. В данной серии публикуются оригинальные работы, обзоры и краткие сообщения по следующим научным направлениям: математическое моделирование физических процессов и свойств веществ, численные и аналитические методы решения прикладных задач математической физики и механики сплошной среды; вычислительная математика и применение математических методов и электронно-вычислительной техники в научных исследованиях; вопросы программирования; вопросы структуры алгоритмов и программ для современных ЭВМ; вопросы создания вычислительных комплексов и сетей ЭВМ. Главный редактор - д-р ф.-м. наук Р.М.Шагалиев

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

Пособие содержит комплекс работ по моделированию объектов классической механики. Существенное внимание уделено фундаментальным положениям механики.

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

СибЖВМ - единственный общероссийский журнал по вычислительной математике, издающийся за Уралом с привлечением авторов и рецензентов со всего СНГ.Основные направления журнала:- вычислительная математика;- математическое моделирование;- прикладная информатика;- автоматизация научных и прикладных исследований.Статьи публикуются на русском и английском языках, в зависимости от языка оригинала.

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

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

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

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

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

Отражены результаты анализа известных литературных источников по использованию вероятностных методик расчета термомеханики твэлов в различных режимах и выбора методики расчета методом Монте-Карло. Разработан программный код (комплекс) CaPpaPI (Complex of Programs for Probability Investigations) вероятностных расчетов, реализующий алгоритм вероятностных расчетов поведения твэла на основе метода Монте-Карло с использованием кодов START-3 и RAPTA-5.2 и теплогидравлических условий, рассчитываемых кодом RELAP/SCDAPSIM. Описана архитектура разработанного кода. Представлено описание реализации многовариантного расчета решения задачи и графический интерфейс. Приведены результаты работы кода. Код учитывает неопределенности входных параметров, погрешности эмпирических корреляций и погрешности расчетных моделей. Код позволяет получать реалистические оценки основных критериальных параметров твэлов в авариях типа LOCA, например, максимальной температуры оболочки, максимальной эквивалентной степени окисления оболочки и др. с оценками их неопределенностей в виде границ доверительного интервала с заданным уровнем доверия.

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

СибЖВМ - единственный общероссийский журнал по вычислительной математике, издающийся за Уралом с привлечением авторов и рецензентов со всего СНГ.Основные направления журнала:- вычислительная математика;- математическое моделирование;- прикладная информатика;- автоматизация научных и прикладных исследований.Статьи публикуются на русском и английском языках, в зависимости от языка оригинала.

Пособие посвящено одному из важнейших этапов конструирования плазменных установок различного назначения – математическому моделированию комплекса процессов, протекающих как в самой высокотемпературной среде – плазме, так и в элементах конструкции, обеспечивающих работоспособность технических устройств. Рассмотрены методы аналитического и численного решения систем уравнений различных типов, приведены решения ряда практически важных задач, которые встречаются студентам при выполнении домашних заданий по основным дисциплинам специальностей 140403 и 140505, курсовом и дипломном проектировании. Пособие основано на материалах лекций, семинарских и лабораторных занятий по методам математического моделирования процессов в плазменных установках, проводимых автором в течение ряда лет в МГТУ им. Н.Э. Баумана

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

Цель работы: изучение математической модели распространения поверхностных электромагнитных ТЕ-волн в плоском неоднородном диэлектрическом волноводе, заполненном средой с нелинейностью, выраженной законом Керра. Материал и методы исследования: проблема сводится к исследованию нелинейного интегрального уравнения с ядром в виде функции Грина. Существование решений интегрального уравнения доказано с помощью метода сжимающих отображений. Для численного решения задачи предложены два метода: итерационный алгоритм (доказана его сходимость), а также метод, основанный на решении вспомогательной задачи Коши (метод пристрелки). Результаты: доказано существование корней дисперсионного уравнения - постоянных распространения волновода. Получены условия, когда могут распространятся k волны, указаны области локализации соответствующих постоянных распространения. Выводы: полученные результаты свидетельствуют о наличии волноводного режима распространения электромагнитных волн в нелинейной среде.