- •Модуль 1 Основы теории разностных схем
- •Комплексная цель модуля
- •Введение
- •1.3. Сетки и сеточные функции. Вводные понятия
- •1.4. Построение разностных схем
- •1.4.1. Интегро-интерполяционный метод
- •1.4.2. Консервативные и не консервативные разностные схемы
- •1.4.3. Исследование аппроксимации дифференциального уравнения
- •1.4.5. Понятие о погрешности аппроксимации, устойчивости и сходимости разностной схемы на примере краевой задачи для оду второго порядка
- •1.4.4 . Исследование погрешности аппроксимации граничного условия
- •1.4.6. Разностные операторы и их свойства
- •1. Разностные тождества
- •1.1. Формулы разностного дифференцирования произведения
- •1.2. Первая разностная формула Грина
- •1.3. Вторая разностная формула Грина
- •2. Оператор разностной производной второго порядка и его собственные значения
- •3. Обобщение оценок для скалярных произведений, содержащих оператор второй разностной производной с переменными коэффициентами
- •1.4.7. Построение разностных схем для уравнения теплопроводности
- •1.4.8. Связь аппроксимации, устойчивости, сходимости для линейных разностных задач
- •1.4.9. Исследование устойчивости разностных задач методом гармоник
- •1.4.10. Неявная схема для уравнения теплопроводности
- •1.4.11. Схема с весами
- •1.4.12. Трехслойные разностные схемы
- •Исследование устойчивости методом гармоник однородного уравнения с нулевыми граничными условиями. Необходимое условие устойчивости
- •1.4.13. Примеры трехслойных разностные схем для уравнения теплопроводности
- •1.5. Принцип максимума и следствия из него
- •1.5.1. Аппроксимация задачи Дирихле для уравнения Пуассона в прямоугольнике
- •1.5.2. Каноническая форма сеточных уравнений
- •1.5.3. Каноническая форма общего вида для сеточного уравнения
- •Разностная схема, аппроксимирующая задачу Дирихле для уравнения Пуассона в прямоугольнике:
- •Принцип максимума
- •1.5.4. Следствия из принципа максимума
- •1.5.5. Теорема сравнения . Мажоранта
- •1.5.6. Следствие из теоремы сравнения
- •1.5.7. Оценка решения неоднородного уравнения
- •1.5.8. Исследование устойчивости разностных схем на основе принципа максимума и следствий из него
- •1.5.9. Примеры решения задач на принцип максимума и следствий из него
- •1.6. Энергетический метод исследования устойчивости
- •1.6.1. Некоторые сведения из теории операторов
- •Знакоопределенные операторы
- •1.6.2. Свойства оператора симметричной разностной производной второго порядка в гильбертовом пространстве
- •Свойства разностных эллиптических операторов
- •Построение энергетической эквивалентности оператора второй разностной производной с переменными коэффициентами
- •Корректность операторно-разностных уравнений
- •1.7.1. Общие сведения
- •1.7.2. Краткие сведения об обратимости оператора в линейных нормированных пространствах
- •1.6.3. Свойства операторов разностных производных первого порядка
- •Каноническая форма двухслойных схем и ее применение к исследованию устойчивости
- •1.7.4. Необходимое и достаточное условие устойчивости в на
- •Энергетическое пространство
- •1.7.5. Исследование устойчивости в случае несамосопряженного оператора
- •Экономичные схемы для нестационарных задач математической физики
- •Схемы переменных направлений
- •Факторизованные схемы
- •Понятие о монотонных схемах
- •Понятие асимптотической устойчивости.
- •Понятие о трехслойных разностных схемах.
- •Введение в метод конечных элементов
- •Кусочно-линейные восполнения сеточных функций.
- •Понятие о методе конечных элементов.
- •10. Общее описание метода.
- •20. Пример построения и исследования конечно-элементой схемы.
- •Свойства приближенного решения мнк.
- •10. Вводные сведения.
- •20. Вариационные свойства приближенного решения.
- •30. Формулы вложения.
- •Сходимость приближенного решения к точному решению.
- •1.9. Заключение по модулю
- •Вопросы для самоконтроля
- •1.10. Проектное задание
- •Контрольная работа №1.
- •Контрольная работа №2.
- •1.11. Тест рубежного контроля
- •Бланк правильных ответов
Принцип максимума
Напомним определение связности сетки.
Опр.
Сетка
называется связной, если
найдется множество узлов:
,
что выполняется:
(7)
.
Некоторые из этих точек, попадающих в окрестность, могут быть и граничными.
Теорема 1. Принцип максимума.
Пусть для сеточной
функции
не являющейся тождественной константой
для любых
выполняется:
,
(считаем, что
сеточный оператор с коэффициентами,
удовлетворяющими условию (5), а сетка
-связная), тогда сеточная функция
не может достигать наибольшего
положительного значения во внутренних
узлах сетки.
,
тогда сеточная функция
не может достигать наименьшего
отрицательного значения во внутренних
узлах сетки.
Доказательство.
Предположим
противное: существует такая точка
.
Тогда найдется хотя бы одна точка
.
Так как
,
среди всех
зафиксируем
.
В силу связанности сетки существует конечное множество узлов:
.
Рассмотрим уравнение
(3′) в точке
.
(8)
.
В силу (5, 6):
в силу (5):
.
Если использовать
«сильное» определение, то точка
- точка в которой сеточная функция
достигает наибольшего значения, то
,
то неравенство (8) не выполняется, т.к.
,
а
.
Более интересным
является применение другого «слабого»
определения наибольшего значения:
- точка в которой сеточная функция
достигает своего наибольшего значения,
то
при условии, что найдется хотя бы одна
точка
.
Будем ориентироваться на это более содержательное определение.
Тогда
,
,
следовательно
.
Выберем точку
в качестве центра шаблона.
Имеем:
.
Так как
,
то
.
Следовательно,
.
Аналогично проделав
для остальных узлов неравенства (3′)
тоже самое, приходим к цепочке равенств:
.
Теперь рассмотрим левую часть неравенства (8).
.
Полученное неравенство противоречит предположению .
1.5.4. Следствия из принципа максимума
Следствие 1.
Пусть выполнены условия предыдущей теоремы, кроме одного может быть тождественной константой на сетке.
Сеточная функция
,
заданная
неотрицательна на границе
,
т.е.
и имеет место
.
Тогда функция .
Доказательство.
Пусть
;
Рассмотрим первый пункт.
Выберем в качестве
точки
так называемый приграничный узел, т.е.
.
Рассмотрим уравнение (3′) в точке .
(9)
1)
,
тогда
,
т.к.
.
2)
,
тогда
,
,
.
Отсюда вытекает,
что левая часть (9) отрицательна, что
противоречит условию, что
.
Рассмотрим второй пункт.
Пусть существует
точка
,
тогда в силу предыдущей теоремы 1
,
следует, что нет такой точки
.
Если точка
не единственная точка сетки, в которой
значение функции отрицательно, то среди
всех этих точек
следует выбрать такую точку, что
,
которая удовлетворяет неравенству.
Если
не одна, то выбирается любая.
Следствие 2.
Однородное
уравнение
с
однородными граничными условиями
,
,
имеет только тривиальное решение.
Доказательство.
Из теоремы
1)
,
если
;
2)
, если
Из 1) и 2) следует, что и , следовательно .
Следствие 3.
Исходное сеточное
уравнение
при условии, что выполняются ограничения
(5), (6), а также сетка связная,
имеет единственное решение.
Это следствие вытекает из известной теоремы алгебры, т.к. соответствующая однородная система имеет тривиальное решение в силу следствия 2, следовательно неоднородная система имеет единственное решение.
Опр.
Сеточное уравнение вида:
(1)
будем называть
уравнением, а оператор L
– оператором из исходного семейства,
если (1) определено во внутренних
узлах сетки
,
а
- связное множество относительно шаблона.
Следствие 4. Сеточный принцип максимума.
Однородное
уравнение
с
неоднородными граничными условиями
,
,
,
где
- заданная функция из исходного семейства
обладает свойством:
,
т.е. максимальное
по модулю значение функции
являющейся решением задачи (2), (3)
достигается на границе
сетки
.
Доказательство.
Результат прямо следует из теоремы 1, т.к. для задачи (2), (3) справедливы одновременно случай 1) и 2) теоремы 1.
1) , следовательно наибольшее положительное значение достигается на границе сетки .
2) , следовательно наименьшее отрицательное значение достигается на границе сетки .
Также для следствия
1 можно сформулировать его аналог для
.
Следствие 1′.
Пусть дано однородное уравнение из исходного семейства вида:
(4) , с граничным условием Дирихле:
(5)
,
.
Тогда решение
задачи (4), (5) не положительно на
,
т.е.
.
