МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

ИМЕНИ М. В. ЛОМОНОСОВА

ФАКУЛЬТЕТ ВЫЧИСЛИТЕЛЬНОЙ МАТЕМАТИКИ И КИБЕРНЕТИКИ

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

Жемухов Умар Хазреталиевич

РАВНОМЕРНЫЕ ЧИСЛЕННЫЕ МЕТОДЫ ДЛЯ

НЕКОТОРЫХ СИНГУЛЯРНО ВОЗМУЩЕННЫХ ЗАДАЧ

Специальность 01.01.07 вычислительная математика

АВТОРЕФЕРАТ

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

Москва 2013

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

Научный руководитель: доктор физико-математических наук, профессор кафедры вычислительных методов факультета ВМК МГУ имени М.В. Ломоносова Андреев Владимир Борисович.

Официальные оппоненты: доктор физико-математических наук, профессор кафедры математики Физического факультета МГУ имени М.В. Ломоносова Нефедов Николай Николаевич, кандидат физико-математических наук Коптева Наталья Викторовна, Department of Mathematics and Statistics University of Strathclyde, Glasgow, Scotland, UK.

Ведущая организация: Институт математики и механики имени Н. Н. Красовского Уральского отделения РАН.

Защита диссертации состоится 25 декабря 2013 г. в 15 ч. 30 мин. на заседании диссертационного совета Д 501.001.43 при Московском государственном университете имени М. В. Ломоносова по адресу: 119991, Москва, ГСП-1, Ленинские горы, 2-й учебный корпус, факультет ВМК, ауд. 685.

С диссертацией можно ознакомиться в научной библиотеке МГУ имени М. В. Ломоносова по адресу: 119992, Москва, Ломоносовский проспект, д. 27. С текстом автореферата можно ознакомиться на официальном сайте факультета ВМК МГУ http://cs.msu.ru/.

Автореферат разослан ноября 2013 г.

Ученый секретарь диссертационного совета, доктор физико-математических наук, профессор Е. В. Захаров

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

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

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

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

Это объясняется тем, что порядок вырожденного уравнения ниже, чем порядок исходного уравнения, из-за чего часть граничных условий оказывается лишней применительно к вырожденной задаче и эти неиспользованные условия приводят к появлению в окрестности границы пограничных слоев. Достаточно много примеров, иллюстрирующих сказанное, содержится в работах Э. Дулана, Дж. Миллера и У. Шилдерса1, Дж. Миллера, Е. О’Риордана и Г. И. Шишкина2.

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

Дулан Э., Миллер Дж., Шилдерс У. Равномерные численные методы решения задач с пограничным слоем. М.: Мир, 1983.

Miller J.J.H., O’Riordan E., Shishkin G. I. Fitted Numerical Methods For Singular Perturbation Problems: Error Estimates in the Maximum Norm for Linear Problems in One and Two Dimensions. World Scientic Co. Inc., Revised Edition, Roos H.-G., Stynes M., Tobiska L. Robust Numerical Methods for Singularly Perturbed Dierential Equations. Springer Series in Computational Mathematics. Vol. 24. 2nd ed. Springer-Verlag, Berlin, 2008.

Roos H.-G. Layer-adapted grids for singular perturbation problems. Z. Angew. Math. Mech. V. 78 (1998).

№ 5. P. 291–309.

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

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

Известно8, что гладкость искомого решения дифференциальной задачи оказывает существенное влияние на точность приближенного решения, найденного конечно-разностным методом. Между тем, если граница области, в которой рассматривается сингулярно возмущенная задача, содержит угловые точки и не предполагается выполнение в этих точках условий согласования9, то кроме пограничных слоев, решение или (и) его производные имеют и угловые особенности. А это часто приводит к существенному усложнению анализа погрешности численного решения или же к снижению порядка точности приближенного метода. Однако есть случаи10, когда рост производных в окрестностях угловых точек не приводит к (существенному) ухудшению погрешности приближенного решения, но обоснование этого факта требует дополнительных специальных исследований или модификаLinss T. Layer-Adapted Meshes for Reaction-Convection-Diusion Problems. Series: Lecture Notes in Mathematics, Vol. 1985. Springer, 2010.

Бахвалов H. C. К оптимизации методов решения краевых задач при наличии пограничного слоя // Ж. вычисл. матем. и матем. физ. 1969. Т. 9. № 4. С. 841-859.

Шишкин Г. И. Сеточные аппроксимации сингулярно возмущенных эллиптических и параболических уравнений. Екатеринбург: УрО РАН, 1992.

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

Волков Е. А. О дифференциальных свойствах решений краевых задач для уравнений Лапласа и Пуассона на прямоугольнике // Тр. МИАН СССР. 1965. Т. 77. С. 89–112.

Андреев В.Б. Равномерная сеточная аппроксимация негладких решений смешанной краевой задачи для сингулярно возмущенного уравнения реакции-диффузии в прямоугольнике // Ж. вычисл. матем.

и матем. физ. 2008. Т. 48. № 1. С. 90-114.

ции существующих методов. В этом направлении следует отметить работы В. Б. Андреева11 и Е. А. Волкова12.

Решения сингулярно возмущенных краевых задач для уравнений в частных производных в областях с угловыми точками могут иметь сложную структуру13, включающую регулярные, параболические (характеристические) и угловые пограничные слои вместе с угловыми особенностями. Для анализа погрешности аппроксимации дискретной задачи, как известно, важно иметь поточечные оценки производных искомого решения. Обычно, эти оценки в случае регулярных (без малого параметра) эллиптических и параболических уравнений получаются как следствие гельдеровых оценок14 вплоть до границы. Но в сингулярно возмущенном случае, чтобы в полной мере выявить характерные свойства решения и производных (в том числе и связанные с угловыми особенностями), целесообразно сначала строить декомпозицию искомого решения на регулярную, погранслойную и угловую составляющие, которая позволяет получить оценки решения и производных отдельно для каждой компоненты. Такое разложение впервые было построено в вышеупомянутой работе Н. С. Бахвалова для получения оценок производных решения сингулярно возмущенной эллиптической задачи. Позже были предложены и другие варианты декомпозиции15, 16 решения сингулярно возмущенных задач.

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

Целью диссертационной работы является исследование равномерной по малому параметру сходимости в равномерной метрике разностных схем на сгущающихся сетках, аппроксимирующих сингулярно возмущенные краевые задачи для двумерного эллиптического уравнения конвекциидиффузии и уравнения теплопроводности при наличии у производных исАндреев В. Б. Равномерная сеточная аппроксимация негладких решений сингулярно возмущенного уравнения конвекции-диффузии в прямоугольнике // Дифференц. уравн. 2009. Т. 45. № 7. С. 954–964.

Волков Е. А. О дифференциальных свойствах решений уравнений Лапласа и Пуассона на параллелепипеде и эффективных оценках погрешности метода сеток// Тр. МИАН. 1969. Т. 105. С. 46–65.

Kellogg R. B., Stynes M. Corner singularities and boundary layers in a simple convection-diusion problem // J. Di. Equ. 2005. V. 213. № 1. P. 81-120.

Ладыженская О. А., Уральцева Н. Н. Линейные и квазилинейные уравнения эллиптического типа.

Изд. 2, М.: Наука, 1973.

O’Riordan E., Shishkin G. I. A technique to provide parameter-uniform convergence for a singularly perturbed convection-diusion equation // J. Comput. Appl. Math. 2007. V. 206. P. 136–145.

O’Riordan E., Shishkin G.I Parameter uniform numerical methods for singularly perturbed elliptic problems with parabolic boundary layers // Appl. Numer. Math. 2008. V. 58. P. 1761-1772.

комых решений особенностей в угловых точках области из-за несогласованности входных данных.

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

Научная новизна. В диссертации исследована равномерная по малому параметру сходимость численных методов для ряда сингулярно возмущенных задач при существенно более слабых предположениях о гладкости искомых решений по сравнению с известными результатами. Именно, для сеточного решения сингулярно возмущенной эллиптической задачи, входные данные которой в угловых точках прямоугольника удовлетворяют лишь условиям согласования нулевого порядка, получена равномерная по оценка сходимости O(N 3/2 ln2 N ) в сеточной норме Lh для всех (0, 1].

Ранее такая оценка была получена для гладкого случая в предположении выполнения в угловых точках условий согласования до 2-го порядка и при условии, что CN 1. В случае параболического уравнения при тех же предположениях о согласованности входных данных, что и для эллиптической задачи получена поточечная оценка погрешности сеточного решения O( +N 2 ln2 N ) ln(j+1), отличающаяся лишь логарифмическим множителем от известной оценки, полученной при условии достаточной гладкости искомого решения.

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

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

Апробация результатов. Результаты, полученные в диссертации, докладывались и обсуждались на 4 конференциях: XVIII Международная научная конференция студентов, аспирантов и молодых учёных Ломоносов 2011 (Москва), Международная научная конференция Fifth Conference on Numerical Analysis and Applications (Болгария, 2012 г.), Научная конференция Ломоносовские чтения (Москва, 2012 г.), Научная конференция Тихоновские чтения (Москва, 2013 г.).

Публикации. Основные результаты диссертации опубликованы в 7 работах, список которых приводится в конце автореферата. Из них три статьи [1–3] в изданиях, рекомендованных ВАК.

Структура и объем диссертации. Диссертация состоит из введения, трех глав и библиографии, включающей 89 наименований. Общий объем диссертации составляет 75 страниц.

Во введении дается обзор литературы по теме диссертации, обосновывается актуальность темы исследования. В нём сформулирована цель диссертации, описана структура диссертации и перечислены основные результаты.

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

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

В области = (0, 1)2 с границей = D N рассматривается смешанная краевая задача для сингулярно возмущенного уравнения конвекции-диффузии:

где a = const > 0, q = const > 0, n единичный вектор внешней нормали к N, а (0, 1] малый параметр. Правая часть и граничные функции предполагаются достаточно гладкими.

Граница области состоит из частей D = 1 3, N = 2 4, где k стороны квадрата, пронумерованные против хода часовой стрелки, Naughton A., Stynes M., Regularity and derivative bounds for a convection-diusion problem with Neumann boundary conditions on characteristic boundaries // Z. Anal. Anwend. 2010. V. 29. № 2. P. 163-181.

начиная с 1 = { (x, y) | x = 0}. Положим, что g1 (y) = g(0, y), g2 (y) = g(1, y) и 1 (x) = (x, 0), 2 (x) = (x, 1).

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

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

В разделе 1.2 ставится разностная задача на сетке = x y, являющейся тензорным произведением двух кусочно-равномерных одномерных сеток Шишкина на отрезке [0, 1]. Сетка в направлении x сгущается вблизи правой границы и состоит из двух частей [0, 1 x ] и [1 x, 1] с шагами H1 = 2(1 x )N 1 и h1 = 2x N 1 соответственно, а в направлении y сетка сгущается в окрестностях обоих концов отрезка и состоит и трех частей соответственно, где x = min{1/2; c ln N }, y = min{1/4; c ln N }, N число узлов сетки в каждом координатном направлении.

Разностная схема18, используемая для аппроксимации задачи (1)–(4), имеет повышенный порядок аппроксимации и является неоднородной в том смысле, что ее вид не одинаков в различных частях сеточной области и зависит от величины параметра. Именно, при < aH1 /2 внутри регулярного погранслоя конвективный член аппроксимируется симметричной первой разностью 2-го порядка точности, а за его пределами применяется аппроксимация односторонней разностью в полуцелых узлах. Если же aH1 /2, то схема однородная и первая производная в уравнении аппроксимируется на всей сетке симметричной первой разностью. В обоих случаях для аппроксимации оператора Лапласа используется классический пятиточечный Clavero C., Gracia J.L., Lisbona F., Shishkin G.I. A robust method of improved order for convectiondiusion problems in a domain with characteristic boundaries // ZAMM. Z. Angew. Math. Mech. 2002. V.

82. № 9. P. 631-647.

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

Основная теорема первой главы:

Теорема 1. Пусть u(x, y) решение исходной задачи (1)-(4), а uh реij шение соответствующей разностной задачи (см. в диссертации (1.10)– (1.13)) на кусочно-равномерной сетке Шишкина. Тогда при (0, 1] справедлива оценка сходимости Доказательство теоремы 1 получается путем объединения результатов теорем 2 и 3 из последующих двух разделов.

В разделе 1.3 исследуется проблема сходимости сеточного решения задачи (1)–(4) при малых значениях параметра ( < aH1 /2) и основное внимание концентрируется на изучении погрешности сеточных функций характеристического и углового слоев, содержащих также угловые особенности.

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

соответствующей разностной задачи (см. в диссертации (1.10)–(1.13)) при < aH1 /2. Тогда при N > N0, где N0 целое положительное число и не зависящее от, верны оценки где В разделе 1.4 рассматривается случай aH1 /2 и используется разностная схема с аппроксимацией конвективного члена симметричной первой разностью. Для обоснования сходимости сеточного решения в этом случае применяется тот же подход, что и в предыдущем разделе, но при получении оценки погрешности приближенного решения гладкой составляющей Andreev V. B. Pointwise approximation of corner singularities for singularly perturbed elliptic problems with characteristic layers // Internat. J. of Num. Analysis and Modeling. 2010. V. 7. № 3. P. 416-428.

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

соответствующей разностной задачи (см. в диссертации (1.10)–(1.13)) при aH1 /2. Тогда в справедлива равномерная по параметру оценка где В разделе 1.5 приводятся результаты численного эксперимента.

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

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

В области QT = {(x, t)| 0 < x < l, 0 < t T } рассматривается первая начально-краевая задача где a постоянная, f (x, t), (x), µ1 (t), µ2 (t) заданные достаточно гладкие функции.

Известно, что при классическом подходе для построения приближенного метода решения задачи (8)–(10) с точностью O( +h2 ) требуется непрерывность производных 2 u/t2, 4 u/x4 вплоть до границы. А это, при условии достаточной гладкости входных функций, подразумевает выполнение, например в угловой точке (0, 0), условий согласования:

Требование выполнения условий согласования не всегда является естественным и может быть слишком обременительным.

В точке (0, 0) мы отказываемся от условия (12), а в (l, 0), для простоты изложения, предполагаем, что аналогичное условие выполнено.

Для искомого решения справедливо представление Функции из (13) являются решениями следующих начально-краевых задач:

Регулярная составляющая решения U (x, t) удовлетворяет в угловых точках условиям согласования до 1-го порядка и является достаточно гладкой.

Решение задачи (15) выписывается в явном виде:

где erfc(z) = 1 erf(z) дополнительная функция ошибок.

Чтобы изучить характер угловой особенности, из этого представления находятся производные функции V (x, t) :

Бижанова Г. И. Решение в пространствах Гельдера краевых задач для параболических уравнений при рассогласовании начальных и граничных данных // Современная математика. Фундаментальные направления. 2010. Т. 36. С. 12–23.

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

Для численного решения задачи (8)–(10) используется неявная четырехточечная разностная схема на равномерной сетке h = h, где 0,..., M }. Множество внутренних узлов обозначается через h = h = {(xm, tj ) : m = 1,..., N 1, j = 1,..., M }.

В соответствие (8)–(10) ставится разностная задача где yt,m = (ym ym )/, yxx,m = (ym+1 2ym + ym1 )/h2, fm = f (xm, tj ).

Основной результат второй главы:

Теорема 4. Пусть ym решение разностной задачи (17), а u(x, t) решение дифференциальной задачи (8)–(10), (11). Тогда справедлива оценка погрешности Доказательство данной теоремы приводится в следующих разделах.

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

Решение разностной задачи (18) представляется в виде Справедливы оценки Из (19), используя последнее неравенство, получим априорную оценку.

Теорема 5. Пусть Fm сеточная функция, убывающая на бесконечности. Тогда для решения задачи (18) справедлива поточечная оценка В разделе 2.3 строится дискретный аналог декомпозиции (13) и для соответствующих разностных схем, аппроксимирующих задачи (14) и (15), доказываются оценки сходимости.

В частности, погрешность численного решения zm = ym u(xm, tj ) представляется в виде суммы Первое слагаемое этой суммы предствляет собой решение разностной заk j дачи, аналогичной (18), с правой частью m = LV (xm, tj ) LVm, где L разностный оператор из (17). А второе слагаемое есть решение разностной схемы где погрешность аппроксимации m = LU (xm, tj ) LUm.

Такое разложение позволяет провести анализ погрешности приближенного решения отдельно для каждой компоненты.

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

Теорема 6. Пусть V (x, t) решение дифференциальной задачи (15) на отрезке [0, l], а Vm решение соответствующей разностной схемы (см.

в диссертации (2.20)). Тогда справедлива оценка Оценка для j гладкой компоненты погрешности численного решения получается при помощи стандартных рассуждений с применением принципа максимума.

В разделе 2.4 приводятся результаты численных расчетов.

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

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

В области QT рассматривается начально-краевая задача (8)–(10) с малым параметром 2 ( (0, 1]) перед старшей производной, вместо коэффициента a2. Правая часть и граничные функции предполагаются достаточно гладкими, последние удовлетворяют только условиям согласования нулевого порядка:

Решение обсуждаемой задачи при 0 имеет вдоль границ x = 0 и x = l параболические пограничные слои, где решение быстро меняется, а старшие производные имеют в угловых точках особенности порядка c t1.

Для численного решения указанной задачи используется четырехточечная неявная схема на равномерной по времени и кусочно-равномерной по пространству сетке: h = h, где = {tj = j, = T /M, j = 0,..., M } равномерная сетка по временной переменной, а h – кусочно равномерная сетка Шишкина с шагами h = 4N 1 на отрезках [0, ], [l, l] и H = 2(l 2)N 1 на [, l ], где = min{l/4, 2 ln N }. Множество внутренних узлов обознчается через h = h = h QT.

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

где U (x, t) – регулярная составляющая, V (x, t) – сингулярная часть, состоящая из суммы двух погранслойных функций V1 (x, t) и V2 (x, t), которые, в свою очередь, включают в себя угловые особенности в точках (0, 0) и (l, 0) соответственно, а u(x, t) – остаточный член, являющийся достаточно гладким, без угловых особенностей и пограничных слоев.

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

Теорема 7. Для функций U (x, t), V (x, t) и u(x, t) из вышеприведенного разложения, определенных в QT, при 1, 0 n + 2m 4, n, m Z, n, m 0 справедливы оценки При доказательстве данной теоремы пользуемся принципом максимума, методом расширения области, а также мотодом функции Грина.

Для остаточного члена справедливы следующие оценки:

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

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

В частности, для погрешности регулярной составляющей сеточного решения получена оценка Аналогичная оценка справедлива и для сеточной аппроксимации остаточного члена.

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

Приведем эти оценки для дискретной функции V1,i, аппроксимирующей погранслойную функцию V1 (xi, tj ), содержащую также угловую особенность в точке (0, 0).

Лемма 1. Пусть V1,i решение разностной задачи (см. в диссертации (33) при k = 1), а V1 (x, t) есть решение соответствующей дифференциальной задачи. Тогда вне области пограничного слоя при 2, N N справедлива оценка погрешности При доказательстве этой леммы используется оценка для V1 (xi, tj ) из теоремы 7 и специальная сеточная барьерная функция, учитывающая поведение погранслойной функции.

Лемма 2. Пусть V1,i решение разностной задачи (см. в диссертации (33) при k = 1), а V1 (x, t) есть решение соответствующей дифференциальной задачи. Тогда в области пограничного слоя справедлива оценка погрешности Доказательство леммы 2 получается путем разбиения погрешности аппроксимации схемы на две части, одна из которых характерезует поведение погрешности вблизи угловой точки, а другая вне конечной окрестности этой точки.

Таким образом, объединяя оценки (22)–(24), получим основной результат третьей главы:

Теорема 8. Пусть yi решение исходной разностной задачи, а u(x, t) решения сингулярно возмущенного варианта дифференциальной задачи (8)–(10), (21). Тогда справедлива оценка погрешности В разделе 3.4 проводится анализ результатов численного эксперимента.

Основные результаты, выносимые на защиту 1. Для сингулярно возмущенного эллиптического уравнения, вырождающегося при = 0 в уравнение первого порядка по одной из переменных, со смешанными краевыми условиями на сторонах прямоугольника и негладким в окрестностях угловых точек решением доказано, что решение неоднородной монотонной разностной схемы на кусочноравномерной сетке Шишкина сходится в норме Lh равномерно относительно малого параметра со скоростью O(N ln N ) для всех (0, 1], где N количество узлов сетки в каждом направлении.

2. Доказано, что решение четырехточечной неявной разностной схемы, аппроксимирующей на отрезке начально-краевую задачу для уравнения теплопроводности при наличии угловой особенности у производных решения, сходится на равномерной сетке к точному решению со скоростью O( + h2 ) ln(j + 1), где h шаг сетки по пространственной переменной, а tj = j узлы сетки по временной переменной.

3. Для решения четырехточечной неявной разностной схемы, аппроксимирующей на равномерной по времени и кусочно-равномерной по пространственной переменной сетке Шишкина начально-краевую задачу для сингулярно возмущенного уравнения теплопроводности на отрезке при условии, что входные данные удовлетворяют в угловых точках только условиям согласования нулевого порядка, получена равномерная по малому параметру поточечная оценка погрешности O( +N 2 ln2 N ) ln(j +1), где N количество узлов пространственной сетки, а tj = j узлы сетки в направлении временной переменной.

Благодарности Автор выражает искреннюю признательность своему научному руководителю, профессору Андрееву Владимиру Борисовичу за постановку задач и постоянное внимание к работе.

Публикации автора по теме диссертации 1. Жемухов У. Х. Равномерная сеточная аппроксимация негладких решений с улучшенной сходимостью для сингулярно возмущенного уравнения конвекции-диффузии с характеристическими слоями // Ж. вычисл.

матем. и матем. физ. 2012. Т. 52. № 9. С. 1633-1654.

2. Жемухов У. Х. О сходимости численного решения начально-краевой задачи для уравнения теплопроводности при наличии угловой особенности у производных решения // Вестн. Моск. ун-та., Серия 15. Вычисл.

матем. и киберн., 2013. № 4. С. 9-18.

3. Zhemukhov U. Kh. Uniform Grid Approximation of Nonsmooth Solutions of a Singularly Perturbed Convection - Diusion Equation with Characteristic Layers // Numerical Analysis and Its Applications ( Series: Lecture Notes in Computer Science), V. 8236, P. 562-570, Springer Berlin Heidelberg, 2013.

4. Жемухов У. Х. Равномерная по параметру оценка погрешности неявной четырехточечной разностной схемы для сингулярно возмущенного уравнения теплопроводности с угловыми особенностями // Препринт, ISBN 978-5-317-04584-5, М.: МАКС Пресс 2013, 20 с.

5. Жемухов У. Х. Равномерная сеточная аппроксимация негладких решений с улучшенной сходимостью для сингулярно возмущенного эллиптического уравнения с характеристическими слоями // Сборник тезисов XVIII Международной конференции студентов, аспирантов и молодых ученных Ломоносов–2011, секция Вычислительная математика и кибернетика. М., 2011, C. 123-124.

6. Zhemukhov U. Kh. Uniform grid approximation of nonsmooth solutions to singularly perturbed convection-diusion equation with characteristic layers // The Abstracts of NAA’12: Fifth Conference on Numerical Analysis and Applications, Lozenetz, Bulgaria, Organized by University of Rousse, June 15-20, 2012, P. 59.

7. Жемухов У. Х. Об оценке погрешности численного решения сингулярно возмущенного уравнения теплопроводности с угловыми особенностями // Сборник тезисов ежегодной научной конференции Тихоновские чтения, секция Математическое моделирование и вычислительные методы. М.: МАКС Пресс 2013, С. 42-43.