- •Модуль 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.5. Принцип максимума и следствия из него
1.5.1. Аппроксимация задачи Дирихле для уравнения Пуассона в прямоугольнике
В области требуется найти решение уравнения:
, (1)
если заданы граничные условия:
. (2)
Будем предполагать необходимую гладкость функции правой части и граничных условий.
Рассмотрим частный случай, когда
. (3)
Для задачи (1), (2) при условии (3) справедлив принцип максимума:
. (4)
Максимальное по модулю решение этой задачи достигается на границе.
. (5)
(4) является частным случаем (5), которое является неравенством определяющим, устойчивость задачи (1-3) по граничным условиям.
При построении дискретных аналогов задачи (1), (2), будем стремиться к наследованию свойства (4).
Для аппроксимации задачи (1), (2), область покроем равномерной сеткой с шагами по оси Ох1 и по оси Ох2 .
.
Внутренностью сетки будем называть множество точек:
, с
граничными узлами:
.
Будем использовать пятиточечный шаблон, включающий в себя точки:
.
Так как уравнение с постоянными коэффициентами, то будем использовать непосредственную аппроксимацию.
, .
.
Подставим это все в (1), получим.
(6) .
(7) .
В системе n уравнений и n (n = N1N2 ) неизвестных, следовательно задача (6), (7) имеет и притом единственное решение (будет показано далее).
Так же будет показано, что решение данной задачи непрерывным образом зависит от граничных условий и правой части (устойчивость по правой части и граничным условиям).
1.5.2. Каноническая форма сеточных уравнений
Запишем (6) предыдущего параграфа в виде:
(1) .
Систему (1) можно записать несколько иначе, если граничные условия (7) учесть в правой части.
Если обратиться к шаблону, задачи, то выделяется две совокупности узлов, т.е. выделяется узел и остальные узлы, которые будем называть окрестностью узла .
Уравнение (1) вместе с граничными условиями является частным случаем уравнения вида:
(2) ;
– центр шаблона,
– окрестность центра шаблона,
Выделим исходное семейство сеточных уравнений.
Коэффициенты, участвующие в уравнении (2) должны удовлетворяют следующим условиям:
(5)
Для простоты, будем предполагать, что Р из внутренних узлов сетки.
Сетка = является связной относительно используемого шаблона, если:
Опр. Для любых узлов , найдется конечное множество узлов каждый из которых принадлежит сетке , таких, что:
.
Первоначально предполагается, что хотя бы один из узлов имеет не пустую окрестность( не граничный узел).
1.5.3. Каноническая форма общего вида для сеточного уравнения
Разностная схема, аппроксимирующая задачу Дирихле для уравнения Пуассона в прямоугольнике:
(1) ;
(2) .
В предыдущем параграфе была записана в канонической форме:
Напомним, что сеточные уравнения имеют следующий канонический вид:
(3) ;
(4) ..
Сравнивая с (3), получаем:
,
, .
Очевидно, что
(6) .
Разностное условие (4) можно трактовать, как частный случай записи (3) в которой , если , в этом случае .
Возможна и другая трактовка граничного условия (4).
- произвольные положительные числа,
, .
Кроме того, сеточное уравнение (3) в канонической форме может быть записано в виде:
(3′) .
Формально граничные узлы сетки можно определить двумя способами (в случае условия Дирихле):
те узлы сетки, значения в которых задано, образуют множество .
.