WWW.DISS.SELUK.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА
(Авторефераты, диссертации, методички, учебные программы, монографии)

 

На правах рукописи

ПАНИН АЛЕКСАНДР АНАТОЛЬЕВИЧ

ОЦЕНКИ ТОЧНОСТИ ПРИБЛИЖЁННЫХ РЕШЕНИЙ

И ИХ ПРИМЕНЕНИЕ В ЗАДАЧАХ

МАТЕМАТИЧЕСКОЙ ТЕОРИИ ВОЛНОВОДОВ

Специальность 01.01.03

МАТЕМАТИЧЕСКАЯ ФИЗИКА

АВТОРЕФЕРАТ

диссертации на соискание учёной степени кандидата физико-математических наук

Москва — 2009

Работа выполнена на кафедре математики физического факультета МГУ имени М.В.Ломоносова.

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

Официальные оппоненты: доктор физико-математических наук, профессор М. Л. Гольдман, доктор физико-математических наук, профессор А. С. Ильинский

Ведущая организация: Институт Прикладной Математики Российской Академии Наук

Защита диссертации состоится « » 2009 г. в часов на заседании Диссертационного совета Д 501.002.10 при Московском государственном университете имени М.В.Ломоносова по адресу:

119991, Москва, ГСП-1, Ленинские горы, МГУ, д. 1, стр. 2, физический факультет, ауд. №.

С диссертацией можно ознакомиться в библиотеке физического факультета МГУ.

Автореферат разослан « » 2009 г.

Учёный секретарь Диссертационного Совета Д 501.002. доктор физико-математических наук Ю. В. Грац

Общая характеристика работы

Актуальность. С ростом возможностей вычислительной техники, позволивших решить весьма сложные краевые задачи математической физики, все более ощущается отсутствие простых алгоритмов, позволяющих получить оценки точности приближённого решения. Например, в настоящее время существуют ориентированные на широкую аудиторию коммерческие комплексы программ, такие как PDE Toolbox в MatLab’е и FEMLab, реализующие алгоритм метода конечных элементов для всех основных линейных краевых задач в ограниченных областях размерности от 1 до 3, однако в этих программах не встроено вообще никакого механизма оценки точности полученного приближённого решения и пользователю не остается ничего другого, как сгустить сетку и визуально оценить скорость сходимости, хотя достоверность такой оценки никак не обоснована. Сходимость проекционных методов для линейных задач математической физики была строго доказана в 1960–70-х годах, однако в теоретических работах по численным методам обычно ограничиваются доказательством самого факта сходимости и нахождению её порядка.

Большинство известных алгоритмов оценки погрешности применимо только к краевым задачам для уравнения Пуассона или, более общо, к задачам с положительно определёнными операторами, но, к сожалению, ни один из них пока не реализован в виде общедоступных комплексов программ. Вероятно, наиболее простым в реализации будет алгоритм, основанный на апостериорных оценках, полученных недавно С. И. Репиным (2000 г. и далее) и не требующих вычисления каких-либо общих констант, кроме оценки сверху константы в неравенстве Пуанкаре—Фридрихса. Однако имеются принципиальные препятствия для перенесения этого метода на задачи с незнакоопределёнными операторами и, в частности, на задачи для уравнения Гельмгольца. Известные оценки для отклонения приближённого решения для всех основных уравнений математической физики были суммированы в работах М. Накао (2001—2007 гг.); в настоящей диссертации на их основе получены улучшенные оценки, построен и реализован алгоритм для оценки точности вычисления собственных значений оператора Лапласа и построен значительно более сложный алгоритм оценки точности решения для уравнения Гельмгольца, где коэффициент k 2 может быть переменным. Предложенные алгоритмы пригодны не только для метода конечных элементов, но и для любого проекционного метода, для которого можно ввести аналог шага сетки.

Ещё одной из актуальных проблем является построение специальных проекционных методов. Так, в работах И. Бабушки и М. Меленка предложен так называемый обобщённый метод конечных элементов (generalized FEM), суть которого в выборе особых базисных функций (функций формы), а в работах С. Саутера показано, как такой специальный выбор позволяет улучшить сходимость. С другой стороны, интересно и практически ценно так называемое явление суперсходимости (superconvergence phenomenon), которое заключается в том, что в некоторых точках расчётной области (положение которых можно определить) приближённое решение сходится к точному быстрее, чем гарантируется общей оценкой в интегральной или равномерной норме.

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

Полученные в диссертации оценки точности вычисления собственных значений позволяют для широкого класса волноводов со сложной геометрией решить вопрос о существовании ловушечных мод прямым расчетом. Теоретические исследования волноведущих систем ведутся довольно давно. Их начало положено классическими работами А. Н. Тихонова и А. А. Самарского. Дальнейшее активное развитие этой области математической физики в нашей стране связано с именами Г. В. Кисунько, П. Е. Краснушкина, Е. И. Моисеева, А. Г. Свешникова, Р. В. Хохлова, В. П. Шестопалова и ряда других учёных. Из зарубежных специалистов можно назвать Ф. Реллиха, Д. Джонса, П. Вернера, П. Экснера.



Введём некоторые термины. Полубесконечной трубой или просто трубой будем называть множество вида T = R+, где — ограниченная односвязная область в R1 или R2. Волноведущей системой V назовём связную область в R2 или R3, вне некоторого шара (круга) представляющую собой объединение конечного числа непересекающихся труб.

Задача о возбуждении такой системы гармоническим током f (x)eit в скалярном приближении имеет вид u + k 2 q(x)u = f, где функции f и q 1 финитны, а аргумент x обозначает вектор всех пространственных координат. Также важно исследовать наличие ловушечных мод, то есть таких функций u 0, что Значения k 2, при которых такое ненулевое решение u существует, называются собственными значениями, а u — собственными функциями спектральной задачи (1). Эта тема восходит к работам Реллиха, который указал на возможность существования собственных значений у спектральных задач в неограниченных областях. Джонсом введено понятие непрерывного спектра и получен ряд результатов относительно его непустоты, а также относительно наличия собственных значений. Оказалось, что непрерывный спектр волновода образуют частоты, для которых k 2 [1 ; +). Здесь 1 — наименьшее собственное значение задачи Дирихле для оператора Лапласа на сечениях полубесконечных труб: Квадратные корни из этих собственных значений, то есть числа название частот отсечки, потому что при переходе частоты возбуждения регулярного и локально нерегулярного волновода через каждую из них добавляется новая мода распространяющихся в волноводе волн. Кроме того, при гармоническом возбуждении регулярного волновода на частотах, равных частотам отсечки, за исключением случая, когда возбуждение ортогонально соответствующей собственной функции сечения, в волноводе не устанавливается режим гармонических колебаний, а происходит рост амплитуды поля пропорционально t. Если же частота возбуждения отлична от частоты отсечки, то устанавливается гармоническое поле, то есть распространяющиеся волны. Поскольку нормальным режимом работы волновода как передатчика энергии (информации) является не резонансный режим, а распространение волн, то именно последний случай представляет наибольший практический интерес. Чтобы гарантировать отсутствие резонансного режима в регулярном волноводе при данной частоте возбуждения, достаточно доказать её отличие от всех частот отсечки. Таким образом, если для первых N собственных значений сечения {i }N найдены интервалы [i ; i ] (то есть показано, что i [i ; i ]), причём для данного k верно k 2 < N и ни при каком i = 1,..., N k 2 [i ; i ], то можно гарантировать (и в этом состоит одна из целей диссертации), что k не совпадает ни с одной из резонансных частот и, следовательно, режима резонанса в регулярном волноводе на данной частоте не будет.

Резонансное множество регулярного волновода исчерпывается частотами отсечки, а в нерегулярных (деформированных или со вставкой, то есть при q 1) волноводах могут существовать упомянутые выше ловушечные моды. С конца 80-х — начала 90-х годов начинается серия работ, посвящённых наличию ловушечных мод у изогнутых волноводов постоянного сечения, само появление которых оказалось некоторой неожиданностью. Так, в статье Экснера и Шебы это показано для плоского (двумерного) волновода с достаточно гладкой границей, а в диссертации Крейцирика — для двумерной полосы, которая изогнута в пространстве.

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

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

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

Научная новизна. В диссертации предложены:

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

• Метод прямого численного исследования ловушечных мод волноведущих систем.

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

• Метод вычисления оценки погрешности приближённого решения, полученного проекционным методом, для уравнения Гельмгольца (в L2 - и H 1 -нормах).

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

Практическая ценность. Подход, предлагаемый в работе, даёт возможность численно установить существование ловушечной моды в волноведущей системе в случаях, когда это непосредственно не следует из известных теоретических результатов. В частности, такой системой может быть волновод с резонатором сложной формы, для собственного значения которого нельзя получить хорошую оценку сверху, вписав в него стандартную область с известным спектром: параллелепипед, шар и т. п. Другой пример — «изломанный» волновод с негладкой границей (поверхностью). Этот подход также позволяет найти частоты, на которых гарантируется работа волновода в режиме распространения гармонических колебаний, то есть не происходит резонанса. Полученные в работе двусторонние оценки собственных значений области (в том числе и задачи «с весом») представляют и самостоятельный интерес, а подпрограмма, написанная для вычисления параметров дискретизации задачи методом конечных элементов в пакете PDE Toolbox, может найти широкое применение для оценки погрешности приближённых решений.

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

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

Апробация работы. Основные результаты диссертации докладывались на Международной конференции «Дифференциальные уравнения и смежные вопросы», посвящённой памяти И. Г. Петровского (в 2007 г.), Международной конференции «Дифференциальные уравнения, теория функций и приложения», посвящённой 100-летию со дня рождения акад. И. Н. Векуа (в 2007 г.), и Международной конференции студентов, аспирантов и молодых учёных «Ломоносов» в 2006 и 2008 годах, на научном семинаре кафедры математики физического факультета МГУ, на научном семинаре «Численные методы электродинамики» физического факультета МГУ под рук. проф. А. Г. Свешникова и проф. А. С. Ильинского (в 2008 г.).

Публикации. По теме диссертации опубликовано 8 печатных работ.

Структура диссертации. Диссертация состоит из введения, четырёх глав, списка литературы (69 наименований) и 18 рисунков; изложена на 90 страницах.

Содержание работы В начале первой главы вводятся основные понятия, используемые в работе. Пусть — произвольная ограниченная область в евклидовом пространстве Rn, граница которой далее считается достаточно гладкой с тем, чтобы было возможно применять формулу Гаусса—Остроградского. Зададимся произвольным конечным набором функций {1,..., N }, носитель которых лежит в, и образуем натянутую на них линейную оболочку SN. Даже если эти функции бесконечно гладкие, произвольную функцию v из L2 () или H0 () нельзя приблизить элементом из SN относительно норм этих пространств, поскольку и в том, и в другом пространстве имеются элементы, ортогональные ко всем функциям {1,..., N }. Поэтому имеет смысл говорить лишь о близости SN и некоторого линейного подпространства X() пространства H0 (). Очевидно, что для дальнейшего нужно, чтобы этому пространству принадлежали обобщенные решения основных краевых задач математической физики, рассматриваемых в, а именно задачи Дирихле и задачи на собственные значения для оператора Лапласа. Поэтому вполне естественно взять в качестве X() линейное пространство, образованное функциями из H0 (), являющимися решениями однородных краевых задач Дирихле для уравнения Пуассона, правые части которых принадлежат L2 (). Тогда близость SN и X() будем характеризовать наименьшей константой C(N ), удовлетворяющей неравенству где. обозначает норму в L2 ().

Константа C(N ) для произвольного проекционного метода играет ту же роль, которую в проекционно-разностных методах играет шаг сетки. Например, если — прямоугольник и в качестве {1,..., N } используются кусочно билинейные конечные элеh менты на квадратной сетке шага h, то C(N ) =. В диссертации приведены и другие примеры конкретных пространств конечных элементов, используемых на практике и удовлетворяющих условию (2). В частности, показано, что ему также удовлетворяют пространства, связанные с разложением по собственным функциям, которые могут быть использованы для решения краевой задачи в том случае, если эти собственные функции известны. Для каждого из рассмотренных конечномерных подпространств указана оценка сверху для константы C(N ).

Для удобства обозначим через PN проектор, ставящий каждой функции v H0 () функцию, реализующую точную нижнюю грань в левой части формулы (2). Его существование легко обосновать, например с помощью теоремы Рисса. Тогда (2) принимает вид Нетрудно показать, что из требования (3) следует оценка верная для всех u H0 (). В то же время, можно показать, что выполняются неравенства для любой функции v из X() и для любой u H0 (), где Cp — константа в неравенстве Пуанкаре—Фридрихса верном при всяком w H0 (). Существенно, что если для некоторого пространства SN удаётся доказать (2), то отсюда автоматически следует (4) и, следовательно, для дискретизации в любом таком подпространстве выполняются полученные в диссертации оценки.

В этой же главе поставлена и решена проблема вычисления двусторонней оценки использующая лишь свойство (4) пробного подпространства и не требующая предварительного вычисления каких-либо других величин. Полученная оценка имеет вид Здесь m — оценка сверху m-го собственного значения сверху проекционным методом в пространстве SN, 1 = sup, C(N ) — константа из (2). Основная идея метода состоит в использовании для пробного функционального подпространства оценки (4) в сочетании с минимаксной характеризацией собственных значений, предложенной Р. Курантом:

где в обеих формулах {1,..., m1 } — наборы из (m 1) функций, принадлежащих H0 ().

В процессе подготовки работы написана подпрограмма на MatLab’е, вычисляющая оценку сверху для C(N ) в случае аппроксимации задачи в выпуклой многоугольной области в пространстве кусочно линейных конечных элементов. Именно такое пространство используется в PDE Toolbox. На вход программы поступают данные о сетке, которые могут быть непосредственно экспортированы из PDE Toolbox’а в рабочее пространство MatLab’а. Результат использовался для вычисления оценок собственных значений конкретных областей, что было существенно использовано ниже в главе 2.

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

В главе 2 полученные оценки были применены к задачам математической теории волноводов. Как уже было сказано выше, собственные значения {i } поперечного сечения волновода определяют его основные свойства, а именно частоты отсечки i, на которых происходит добавление новой моды распространяющихся волн и которыми исчерпывается резонансное множество регулярного волновода, и непрерывный спектр k 2 [1 ; +), на частотах ниже которого излучение не происходит. Вычислив или оценив эти собственные значения, можно найти диапазоны частот одно-, двух-, трёхмодового режима и т. д. как интервалы, лежащие между частотами отсечки. Более того, к изолированным ловушечным модам (не вложенным в непрерывный спектр) применм уже упоминавшийся принцип минимакса. Применяя формулы (6)—(7) к задаче (1), можем записать:

где 1,..., m1 H0 (V ). Поэтому если первые M (M = 0, 1, 2,... ) значений µ, доставляемые формулой (8), оказались меньше 1, то можно утверждать, что в волноводе существуют ловушечные моды, квадраты частот которых равны µ. Для численного моделирования важно, что приближённые собственные значения, вычисленные произвольным конформным проекционным методом согласно формуле (9), с необходимостью µj < 1. Отсюда следует, что нахождение M таких приближённых собственных значений {j }M, что µj < µ j=1 1, гарантирует существование в волноведущей системе не менее M ловушечных мод, лежащих ниже непрерывного спектра. Поскольку, в свою очередь, 1 может быть оценено снизу способом главы 1, во многих конкретных случаях оказывается возможным численно доказать существование таких ловушечных мод. С целью иллюстрации этого подхода в диссертации сначала рассмотрено несколько плоских волноводов, причём для того, чтобы триангуляция не нарушала конформность дискретных подпространств, вогнутый участок границы был заменён ломаной из касательных с достаточно большим числом звеньев. Далее метод был применён к конкретным трёхмерным волноведущим системам сложной геометрии и для них установлено существование ловушечных мод. В частности, для резонатора с присоединённой к нему полубесконечной трубой наличие ловушечной моды теоретически обосновано для достаточно больших резонаторов. Поэтому для конкретных систем вопрос остаётся открытым, если собственное значение задачи Дирихле для оператора Лапласа в полости резонатора не оказалось меньше 1. В работе рассмотрено два таких случая.

В обоих сечением полубесконечной трубы является правильный треугольник, двусторонние оценки собственных значений которого получены в главе 1. В частности, для первого собственного значения имеем 50,9760 1 52,6382. В первом случае резонатор имеет форму прямоугольного параллелепипеда размерами 1 1 0,5, для которого 59,2176 < 1 = 6 2 < 59,2177. Как видно, этой оценки «не хватает». В то же время, взяв ром случае резонатор имеет форму куба со стороной 1, поэтому существование первой ловушечной моды тривиально: 1 = 3 2 29,6088. Но 2 = 6 2, поэтому для второй ловушечной моды требуется взять собственное значение большей подобласти, откуда имеем оценку µ2 50,8783 < 1 = 50,9760. Приближённо вычисленные ловушечные моды для всех примеров изображены на рисунках, некоторые из которых приведены в данном автореферате.

Ловушечная мода. Пример 6, симметричная мода. µ = 106,262, 1 > 106,810.

Ловушечная мода. Пример 6, антисимметричная мода. µ2 = 106,718, 1 > 106,810.

Ловушечная мода. Пример 9. Резонатор. µ = 47,5330, 1 > 50,9760.

Ловушечная мода. Пример 10. Резонатор, вторая мода. µ2 = 50,8783, 1 > 50,9760.

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

В главе 3 для краевой задачи в ограниченной области ности приближённого решения проекционным методом, пробное пространство которого удовлетворяет условию (2). В диссертации сформулированы явные алгоритмы вычисления оценки погрешности.

Для изложения результатов потребуется ввести некоторые обозначения. Пусть C(N ) и Cp — константы, введённые выше; пусть, далее, оператор A ставит в соответствие произвольной функции u H0 () обобщённое решение Au H0 () краевой задачи Дирихле для оператора Лапласа Пусть матрицы G и D определяются так:

Потребуются также константы: операторная норма PN APN |SN H0, равная наибольшему собственному значению задачи а также M = (min— соб. знач.G=D ||)1. Отметим, что в последних двух случаях возникли обобщённые задачи на собственные значения, в которые вошла матрица Грама D базиса {i }N пробного пространства SN в смысле скалярного произведения (., поскольку базис не является ортонормированным.

Пусть ещё Тогда можно получить следующие оценки погрешности приближённого решения (где f — правая часть, u — точное решение, v —приближённое):

где предполагается (это проверяется в процессе расчёта), что < 1.

Более того, константа показывает, во сколько раз приближённое решение v более удалено от точного решения u, чем наилучшее приближение последнего в подпространстве SN, поскольку верна оценка нения градиента. Прежде подпространство H 1 (). Далее каждая компонента v проектировалась на SN. Полученная векторv. Затем, после вычисления L2 -норм функция была обозначена вычислялась оценка погрешности которая в разных случаях может оказаться как более, так и менее точной, чем первая из оценок (12).

В главе 4 описан метод построения специального проекционного метода для краевой задачи для ОДУ 2-го порядка где f L2 (0; 1), c L (0; 1), а b(x) «кусочно-W », что означает, что отрезок [0; 1] можно разбить на конечное число отрезков [xk ; xk+1 ], на каждом из которых b(x) принадлежит W (xk ; xk+1 ). Непрерывности b на всём [0; 1] не требуется. Ограничений на знаки коэффициентов тоже нет.

К задаче (13) применяется метод Галёркина—Петрова, т. е. решение ищется в виде линейной комбинации пробных функций {n }M, а условие ортогональности невязки поверочным функциям {m }M записывается в виде Основная идея состоит в использовании в качестве поверочных функций проекционного метода функций, удовлетворяющих однородному сопряжённому уравнению.

Точнее, пусть дано множество точек 0 = x0 < x1 <... < xN 1 < xN = 1, в которых мы хотим знать решение без погрешности. В это множество необходимо включить точки разрыва коэффициента b(x). Предположим, что на каждом интервале (xn1 ; xn ), i = 1,..., N, однородная краевая задача Дирихле для оператора L имеет только тривиальное решение, т. е. число 0 не является собственным значением этих операторов ни на одном из этих отрезков. (Это равносильно однозначной разрешимости задачи для оператора L. Если это условие не выполняется, можно добиться его выполнения, разбив отрезки, где оно нарушается, новыми точками.) Тогда однозначно Поверочные функции выбирались следующим образом, обобщающим обычные для традиционного метода конечных элементов кусочно аффинные функции. Пусть для Доказана следующая теорема:

1. дифференциальная задача (13) в обобщённой постановке однозначно разрешима;

2. задачи (14) для всех отрезков [xn1 ; xn ], n = 1,..., N, имеют только тривиальное 3. пробные функции принадлежат H0 [0; 1] (отметим, что поверочные принадлежат этому пространству по построению);

удовлетворяющие условиям m (xn ) = mn, m, n = 1,..., N 1 (отсюда, в частности, следует линейная независимость как системы указанных линейных комбинаций, так и исходной системы функций {n (x)});

5. СЛАУ метода Галёркина с выбранными функциями {n (x)}N 1 и {m (x)}N 1, построенными по формуле (15), однозначно разрешима.

Тогда верны равенства то есть приближённое решение совпадает с точным в выбранных точках.

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

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

Для сравнения приведены фрагменты графиков решений при одном и том же числе узлов N = 1000. Задача решалась на отрезке [0; 1], для наглядности на графиках представлен отрезок [0,1; 0,2]. График с кружочками представляет приближённое решение, без них — точное; кружочками отмечены значения приближённого решения в узлах сетки. На втором графике (решение предлагаемым методом) обе кривые практически совпадают.

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

Заключение Приведём основные результаты работы.

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

• Построены оценки частот ловушечных мод волноведущих систем.

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

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

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

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

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

Основные публикации 1. Боголюбов А. Н., Малых М. Д., Панин А. А. Временная асимптотика поля, возбуждаемого в волноводе гармоническим током. // Журнал вычисл. мат. и мат.

физики. 2005. Т. 45. № 12. С. 2219—2231.

2. Боголюбов А. Н., Малых М. Д., Панин А. А. Принцип предельной амплитуды для волновода. // Вестник Моск. ун-та. Сер. 3. Физ. Астрон. 2006. № 5. С. 9—13.

3. Панин А. А. О проблеме суперсходимости алгоритмов метода конечных элементов. // Журнал вычисл. мат. и мат. физики. 2008. Т. 48. № 12. С. 2180—2185.

4. Боголюбов А. Н., Панин А. А. Об оценке погрешности приближённого решения эллиптических уравнений с некоэрцитивной билинейной формой. // Вычислительные методы и программирование. 2009. Т. 10. С. 34—48.

5. Панин А. А. Временная асимптотика поля, возбуждаемого в волноводе гармоническим током. // Сборник тезисов международной конференции студентов, аспирантов и молодых учёных «Ломоносов—2006», секция «Физика». М.: Физический факультет МГУ, 2006. Т. 1, С. 133—135.

6. Боголюбов А. Н., Малых М. Д., Панин А. А., Свешников А. Г. О спектральных свойствах краевой задачи для оператора Лапласа. // Международная конференция «Дифференциальные уравнения и смежные вопросы», посвящённая памяти И. Г. Петровского. Сб. тезисов. Москва, 2007. С. 44.

7. Боголюбов А. Н., Малых М. Д., Панин А. А., Свешников А. Г. О спектре краевой задачи Дирихле для оператора Лапласа. // Международная конференция «Дифференциальные уравнения, теория функций и приложения», посвящённая 100-летию со дня рождения акад. И. Н. Векуа. Сб. тезисов. Новосибирск, 2007.

С. 98—99.

8. Панин А. А. Об эффективности одного класса апостериорных оценок. // Сборник тезисов международной конференции студентов, аспирантов и молодых учёных «Ломоносов—2008», секция «Физика». М.: Физический факультет МГУ, 2008.

С. 76—78.





Похожие работы:

«ГОЛУБ Константин Юрьевич РОССИЙСКАЯ ПРАВОВАЯ ПОЛИТИКА В СФЕРЕ ВНЕШНЕЭКОНОМИЧЕСКОЙ ДЕЯТЕЛЬНОСТИ 12.00.01 – теория и история права и государства; история учений о праве и государстве АВТОРЕФЕРАТ диссертации на соискание учной степени кандидата юридических наук Саратов – 2012 2 Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Саратовская государственная юридическая академия Научный руководитель – доктор...»

«Карабан Вадим Михайлович ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ТЕПЛОФИЗИЧЕСКИХ ПРОЦЕССОВ ТЕРМОСТАТИРОВАНИЯ 01.04.14 - Теплофизика и теоретическая теплотехника Автореферат диссертации на соискание ученой степени кандидата физико-математических наук Томск – 2009 2 Работа выполнена на кафедрах прикладной информатики и математики, конструирования и производства радиоаппаратуры Томского государственного университета систем управления и радиоэлектроники. Научный руководитель : доктор...»

«Степанов Вилен Степанович МЕТОДИКА ПРОЕКТИРОВАНИЯ ПРИВОДА НА ОСНОВЕ ВОЛНОВОЙ ПЕРЕДАЧИ С ТЕЛАМИ КАЧЕНИЯ Специальность: 05.02.02 Машиноведение, системы приводов и детали машин Автореферат диссертации на соискание ученой степени кандидата технических наук Москва 2009 г. Работа выполнена на кафедре Системы приводов авиационнокосмической техники Московского авиационного института (государственного технического университета) Научный руководитель : д.т.н., профессор Самсонович Семен...»

«ВИШНИЦКИЙ АЛЕКСАНДР ФЕДОРОВИЧ Исследование и разработка технологии создания микромодулей бесконтактной идентификации для электронных документов Специальность 05.27.06 – Технология и оборудование для производства полупроводников, материалов и приборов электронной техники. АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Москва – 2009 Работа выполнена на кафедре микроэлектроники Московского государственного института электронной техники...»

«Лисеенко Наталья Владимировна СИНТЕЗ И СВОЙСТВА КЕРАМИЧЕСКИХ ПИГМЕНТОВ НА ОСНОВЕ СИСТЕМЫ CaO–RO(R2O3)–SiO2 C ИСПОЛЬЗОВАНИЕМ ПРИРОДНОГО И ТЕХНОГЕННОГО МИНЕРАЛЬНОГО СЫРЬЯ Специальность 05.17.11 – технология силикатных и тугоплавких неметаллических материалов АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Томск 2011 2 Работа выполнена на кафедре технологии силикатов и наноматериалов ФГБОУ ВПО Национального исследовательского Томского...»

«Милованова Любовь Анатольевна Семантико-грамматические свойства и отношения предлога за1, оформляющего винительный падеж, и предлога за2, оформляющего творительный падеж, в современном русском языке Специальность 10.02.01 – Русский язык АВТОРЕФЕРАТ диссертация на соискание ученой степени кандидата филологических наук Челябинск - 2009 Работа выполнен на каф на федре русс ского язы и мето ыка одики пре еподавани русско ия ого яз зыка ГОУ ВПО Челябинск госуда У кий арственны...»

«БУЛАВИНА Мария Александровна ПРАВОВОЕ ПОЛОЖЕНИЕ КАТОЛИЧЕСКОЙ ЦЕРКВИ В РОССИИ В XVIII ВЕКЕ 12.00.01 – теория и история права и государства; история учений о праве и государстве Автореферат диссертации на соискание ученой степени кандидата юридических наук Москва, 2008 2 Диссертация выполнена на кафедре истории государства и права Ставропольского государственного университета и рекомендована к защите на кафедре государственного строительства и права Федерального...»

«Матвеев Евгений Леонидович ОПТИМИЗАЦИЯ КВАНТИЛЬНОГО КРИТЕРИЯ ПРИ ВЫПУКЛОЙ ЦЕЛЕВОЙ ФУНКЦИИ С ПОМОЩЬЮ СТОХАСТИЧЕСКОГО КВАЗИГРАДИЕНТНОГО АЛГОРИТМА Специальность 05.13.01 Системный анализ, управление и обработка информации (авиационная и ракетно-космическая техника) АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук Москва, 2010 Работа выполнена на кафедре Теории вероятностей Московского авиационного института (государственного технического...»

«Завалишина Ольга Валентиновна ПЕДАГОГИЧЕСКАЯ ПОДДЕРЖКА ПОДРОСТКОВ, СКЛОННЫХ К ИНТЕРНЕТ-ЗАВИСИМОСТИ 13.00.01 – общая педагогика, история педагогики и образования АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата педагогических наук Курск 2012 Работа выполнена на кафедре коррекционной психологии и педагогики ФГБОУ ВПО Курский государственный университет доктор педагогических наук, профессор Научный руководитель кафедры коррекционной психологии и педагогики ФГБОУ ВПО...»

«БАЧУРИН Александр Борисович ГИДРОАВТОМАТИКА РЕГУЛИРУЕМОЙ ДВИГАТЕЛЬНОЙ УСТАНОВКИ (РАЗРАБОТКА И ИССЛЕДОВАНИЕ) Специальность: 05.04.13 – Гидравлические машины и гидропневмоагрегаты АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук УФА 2014 Работа выполнена в ФГБОУ ВПО Уфимский государственный авиационный технический университет на кафедре прикладной гидромеханики Научный руководитель : доктор технических наук, профессор Целищев Владимир Александрович...»

«Шаповалов Михаил Сергеевич Палестинский аспект ближневосточной политики Великобритании в 1914-1931 гг. 07.00.03 - всеобщая история АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата исторических наук Томск - 2008 Работа выполнена на кафедре всеобщей истории ГОУ ВПО Омский государственный университет имени Ф.М. Достоевского Научный руководитель : доктор исторических наук, профессор Фоменко Светлана Владимировна Официальные оппоненты : доктор исторических наук,...»

«Тихоненко Алексей Митрофанович ФУНКЦИОНАЛЬНАЯ УСТОЙЧИВОСТЬ И ЭЛЕКТРОМАГНИТНАЯ ЗАЩИЩЕННОСТЬ АВТОМАТИЗИРОВАННЫХ ИДЕНТИФИКАЦИОННЫХ СИСТЕМ ДЛЯ МОНИТОРИНГА И УПРАВЛЕНИЯ ДВИЖЕНИЕМ СУДОВ НА ВНУТРЕННИХ ВОДНЫХ ПУТЯХ НА ОСНОВЕ ПОСЛЕДОВАТЕЛЬНОПАРАЛЛЕЛЬНЫХ СЛОЖНЫХ СИГНАЛОВ Специальность 05.13.06 Автоматизация и управление технологическими процессами и производствами Автореферат диссертации на соискание ученой степени кандидата технических наук Санкт-Петербург 2013 Диссертация выполнена...»

«Горин Антон Анатольевич РАЗВИТИЕ ЛИЧНОЙ И РОДОВОЙ ИДЕНТИФИКАЦИИ СЛУЖИЛОГО СОСЛОВИЯ МОСКОВСКОГО ГОСУДАРСТВА XIV–XVII ВЕКОВ: ГЕРАЛЬДИЧЕСКИЙ АСПЕКТ Специальность 07.00.02 – Отечественная история АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата исторических наук Казань – 2012 Работа выполнена на кафедре отечественной истории Казанского (Приволжского) федерального университета. Научный руководитель : доктор исторических наук, доцент Ибнеева Гузель Вазыховна...»

«Крайнова Екатерина Александровна РАЗРАБОТКА ТЕХНОЛОГИИ СУЛЬФОКАТИОНИТОВ МЕТОДОМ СЕРНОКИСЛОТНОГО ОБУГЛИВАНИЯ РАЗЛИЧНЫХ ВИДОВ ЦЕЛЛЮЛОЗОСОДЕРЖАЩИХ ТВЕРДЫХ ОТХОДОВ 05.17.01 – Технология неорганических веществ АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Москва 2008 Работа выполнена в Российском химико-технологическом университете им. Д.И. Менделеева и Объединенном институте ядерных исследований (ОИЯИ), Дубна Московской обл. Научный руководитель :...»

«ПОЛИТОВ Михаил Сергеевич ЭКСПЕРИМЕНТАЛЬНО-АНАЛИТИЧЕСКИЙ МЕТОД ОЦЕНКИ И ПРОГНОЗИРОВАНИЯ УРОВНЯ ЗАЩИЩЁННОСТИ ИНФОРМАЦИОННЫХ СИСТЕМ НА ОСНОВЕ МОДЕЛИ ВРЕМЕННЫХ РЯДОВ Специальность 05.13.19 – Методы и системы защиты информации, информационная безопасность АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Уфа - 2010 Работа выполнена в ГОУ ВПО Челябинский государственный университет на кафедре вычислительной механики и информационных технологий Научный...»

«Соломыков Василий Сергеевич МЕТОД КОРРЕКЦИИ КОНТУРА ШРИФТОВОГО ЗНАКА ДЛЯ ВЫВОДНЫХ УСТРОЙСТВ С НИЗКИМ РАЗРЕШЕНИЕМ Специальность 05.13.06. – Автоматизация и управление технологическими процессами и производствами (полиграфические средства информации и информационные системы) АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Москва – 2013 2 Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального...»

«Ахмедшин Рамиль Линарович КРИМИНАЛИСТИЧЕСКАЯ ХАРАКТЕРИСТИКА ЛИЧНОСТИ ПРЕСТУПНИКА Специальность 12.00.09. – уголовный процесс, криминалистика и судебная экспертиза, оперативно-розыскная деятельность Автореферат диссертации на соискание ученой степени доктора юридических наук Томск - 2006 Работа выполнена в государственном образовательном учреждении высшего профессионального образования Томский государственный университет на кафедре криминалистики. Научный консультант :...»

«Савченко Андрей Иванович УПРАВЛЕНИЕ ИННОВАЦИОННЫМ РАЗВИТИЕМ ПРОМЫШЛЕННОГО КОМПЛЕКСА РЕГИОНА Специальность 08.00.05 – Экономика и управление народным хозяйством (экономика, организация и управление предприятиями, отраслями, комплексами – промышленность; управление инновациями) АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата экономических наук Ижевск - 2010 1 Диссертация выполнена в ГОУ ВПО Удмуртский государственный университет Научный руководитель : доктор...»

«ДАВЫДОВ Андрей Александрович ЗРЕЛИЩНОСТЬ КАК КУЛЬТУРНАЯ ОСНОВА АНТИЧНОЙ КАРТИНЫ МИРА Специальность 24.00.01 – Теория и история культуры АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата культурологии Иваново 2014 Работа выполнена на кафедре культурологии и литературы ФГБОУ ВПО Ивановский государственный университет, Шуйский филиал Научный руководитель :доктор философских наук, профессор, ФАТЕНКОВ АЛЕКСЕЙ НИКОЛАЕВИЧ Официальные оппоненты : КОЧЕРОВ СЕРГЕЙ...»

«Ерофеева Наталья Валерьевна ИССЛЕДОВАНИЕ СЕГРЕГАЦИИ ГРУЗА НА ЛЕНТОЧНОМ КОНВЕЙЕРЕ ПОД ВОЗДЕЙСТВИЕМ УДАРНЫХ ИМПУЛЬСОВ Специальность 05.05.06 – Горные машины Автореферат диссертации на соискание ученой степени кандидата технических наук Кемерово 2011 2 Работа выполнена в Государственном образовательном учреждении высшего профессионального образования Кузбасский государственный технический университет Научный руководитель : доктор технических наук, доцент Захаров Александр...»








 
2014 www.av.disus.ru - «Бесплатная электронная библиотека - Авторефераты, Диссертации, Монографии, Программы»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.