- •Модуль 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.4. Построение разностных схем
1.4.1. Интегро-интерполяционный метод
Различные физические процессы (распространение тепла или колебаний в среде, процессы газодинамики и аэродинамики ) характеризуются некоторыми интегральными законами сохранения (тепла, массы, количества движения, энергии). При выводе дифференциальных уравнений математической физики обычно исходят из некоторого интегрального соотношения (уравнения баланса), выражающего закон сохранения для малого объема. Дифференциальное уравнение получается из уравнения баланса при стягивании объема к нулю в предположении существования непрерывных производных, входящих в уравнение.
Одно из основных свойств разностной схемы является ее консервативность. Разностная схема называется консервативной, если для нее выполняются разностные аналоги физических законов сохранения. Для получения консервативных разностных схеме естественно исходить из уравнений баланса, записанных для элементарных объемов (ячеек) сеточной области. Входящие в эти уравнения баланса интегралы и производные следует заменить приближенными разностными выражениями. В результате получаем однородную разностную схему. Такой метод получения консервативных схем будем называть интегро-интерполяционным методом (методом баланса).
Итак, метод позволяет строить разностные схемы, удовлетворяющие основным законам сохранения, если этим законам удовлетворяют соответственные непрерывные модели.
Метод рассмотрим на примере краевой задачи для ОДУ второго порядка вида:
(1)
(2)
(3)
В уравнении (1) и условии (2) искомой является функция . Функции и постоянные , , задаются. Относительно этих функции будем предполагать, что
, и ≥ 0 .
Кроме того, будем предполагать, что все функции имеют необходимую по ходу изложения гладкость. При указанных условиях решение задачи (1), (2), (3) существует и единственно.
Уравнение (1) можно рассматривать как уравнение установившегося распределения температуры u(x) в стержне длины l , на одном конце которого поддерживается заданная температура, а на другом происходит теплообмен с окружающей средой.
Для построения решения разностной схемы будем использовать равномерную сетку: ;
;
;
.
Проинтегрируем обе части уравнения (1) по отрезку
.
. (4)
Будем предполагать, что температура. Поток тепла внутри тела
.
Уравнение (2) является несколько измененной формой теплообмена по закону Ньютона.
Введем функцию теплового потока:
. (5)
Тогда уравнение (4) может быть переписано в виде:
. (4׳)
Используем приближенные замены (6-8) :
. (6)
. (7)
, (8)
получим из (4׳) приближенное равенство, которое заменим на точное, не меняя для простоты прежние обозначения.
. (9)
Выразим значение функции , используя значения функции (х) в узлах сетки. Из уравнения (5) получаем, что
.
Интегрируя соотношение на отрезке , имеем:
.
.
Здесь
. (10)
Переходя к точным равенствам, получаем:
, (11)
аналогично
. (12)
Подставим (11) и (12) в равенство (9) получаем:
.
Разделим обе части полученного равенства на :
. (13)
Разностная аппроксимация уравнения (1) построена.
По сравнению с непосредственной аппроксимацией возникли особенности в вычислении коэффициентов разностных уравнений:
.
.
.
Аппроксимация граничного условия (3) очевидна:
. (14)
С граничным условием (2) необходимо аппроксимировать. Для получения его дискретного аналога воспользуемся интегро – итерполяционным методом.
Проинтегрируем обе части уравнения (1) по отрезку , тогда
.
Из условия (2) : .
Из условия (12) : . Подставляя значения в предыдущее уравнение, получаем
; (15)
Где ; .
Формула (15) аппроксимирует граничное условие (2). Таким образом, разностная схема (13), (14), (15) построена.
Полученное разностное граничное условие (15) можно записать:
. (16).
Т.о. краевое условие (2) имеет разностным аналогом уравнение (16).
В дальнейшем решение разностной задачи в отличии от решения дифференциальной задачи будем обозначать буквой y, так что yi = y( xi ),
Объединяя уравнений (13), (16) получаем дискретную модель для задачи (1) – (3) :
(17)
; yN = µ2 .
При исследовании дискретной модели (17) возникают следующие вопросы:
1) о существовании и единственности решения системы линейных алгебраических уравнений (17); 2) каким методом надо решать эту систему; 3) вопрос об аппроксимации дифференциальной задачи (1) – (3) разностной схемой (17); 4) сходиться решение разностной задачи к решению дифференциальной задачи при ? На первые два вопроса можно ответить, преобразуя систему (17) к системе с ленточной треугольной матрицей:
Где , , и .
Граничные условия можно записать в виде:
здесь
; и
,
Из условия ai >0 , β≥0 , di ≥0 следует, что т . е. выполнены условия устойчивости прогонки. Систему (17) можно решать методом скалярной прогонки. На последние два вопроса ответы даются в следующих параграфах.