- •Модуль 1 Основы теории разностных схем
- •Комплексная цель модуля
- •Введение
- •1.3. Сетки и сеточные функции. Вводные понятия
- •1.4. Построение разностных схем
- •1.4.1. Интегро-интерполяционный метод
- •1.4.2. Консервативные и не консервативные разностные схемы
- •1.4.3. Исследование аппроксимации дифференциального уравнения
- •1.4.4 . Исследование погрешности аппроксимации граничного условия
- •1.4.5. Понятие о погрешности аппроксимации, устойчивости и сходимости разностной схемы на примере краевой задачи для оду второго порядка
- •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. Тест рубежного контроля
- •Бланк правильных ответов
Построение энергетической эквивалентности оператора второй разностной производной с переменными коэффициентами
До сих пор мы рассматривали простейшие операторы, не зависящие от пространственных переменных. На практике это редко встречающаяся ситуация. Поэтому получим оценки для оператора разностной производной второго порядка с переменными коэффициентами.
;
.
В двумерном случае рассматривается задача Дирихле.
;
,
В этих задачах мы приходим к оператору второго порядка с переменными коэффициентами.
(1) .
Будем использовать:
(2) .
Получим по возможности более точные оценки вида:
(3) , где константы и не могут быть улучшены.
Опр. 1
Постоянные в неравенстве (3) называются постоянными энергетической эквивалентности или энергетическими постоянными.
Воспользуемся первой разностной формулой Грина для финитных сеточных функций.
(4) .
В соответствии с (4), скалярное произведение (2) запишется в виде:
(*)
. В соотношении (1), если иметь в виду аппроксимацию ОДУ II-го порядка .
Следовательно, , тогда:
(*)
Воспользуемся еще раз формулой (4) в которой , прочитав ее справа на лево и умножив на (-1).
.
Вспомним оценку сверху для , т.е.
(5) , тогда
.
Итак, получили оценку:
(6) .
Будем считать, что - точная константа в (5).
___________________________________________________________
?? Можно ли улучшить оценку (6), уменьшив .
‼ Получить самостоятельно оценку снизу (7).
___________________________________________________________
(7) , где .
Приведем без вывода соответствующие оценки – укажем постоянные энергетической эквивалентности для оператора (двумерный случай):
(8) (8) - дискретный аналог (9):
(9) ; .
(10) , где
(11) ,
(12)
(13) , где
(14)
___________________________________________________________
‼ Получить самостоятельно оценки для: а) и ; б) уточнить формулы (12) и (13).
___________________________________________________________
Корректность операторно-разностных уравнений
1.7.1. Общие сведения
Построенные разностные аналоги краевой задачи для ОДУ II-го порядка, а также задачу Дирихле для уравнения Пуассона, могут быть представлены в виде операторного уравнения:
(1) .
Здесь - известная функция правой части из гильбертова пространства сеточных функций, - искомая сеточная функция, которая может считаться финитной при соответствующем переопределении граничных условий и правой части.
,
,
,
.
Пусть нормы и - нормы, введенные на пространстве функций, являющихся решением (1), норма для , а норма введена на пространстве функций, являющихся правыми частями в (1), т.е. для .
В простейшем случае эти нормы могут совпадать.
Можно считать, что решение – это множество финитных функций.
Опр. 1
Операторное уравнение (1) называется корректным, если:
решение задачи (1) существует и единственно для всех - семейство функций;
для любой пары существует , не зависящая от , такая, что выполняется неравенство:
(2) .
Гильбертово пространство финитных функций обозначим .
С пунктом 2) определения корректности, мы встречались ранее, называя его условием устойчивости по правой части. Новым здесь является то, что в левых и правых частях 2) могут фигурировать разные нормы, а также граничные условия, если они есть, учтены в правой части. Неравенство (2) означает, что малым изменениям входных данных соответствует малое изменение возмущений. Т.е. малое по норме изменение функции приводит к малому изменению нормы .
В формулировке, как и ранее, подразумевается линейность и ограниченность оператора .
Докажем последнее утверждение, т.е., на ряду с задачей (1) рассмотрим задачу (3):
(3) .
Вычтем из (3) (1), получим
, операторы задач (1), (3), (4) совпадают.
(4) .
.
Поскольку , то, следовательно, .