- •Швыдкий Владимир Серафимович дифференциальные уравнения и численные методы
- •Екатеринбург
- •Введение
- •Методы исследования
- •Физический смысл дифференциальных уравнений
- •Глава 1. Дифференциальные уравнения в частных производных
- •Основные положения
- •1.1.1. Характеристики дучп
- •Интерпретация на физической основе
- •1.2. Гиперболические дифференциальные уравнения в частных производных
- •1.3. Параболические дифференциальные уравнения в частных производных
- •Интерпретация с помощью характеристик
- •Интерпретация на физической основе
- •1.4. Эллиптические дифференциальные уравнения в частных производных
- •Глава 2. Математическое описание физических процессов
- •2.1. Общий вид законов сохранения
- •2.2. Характеристическая макроскопическая скорость и диффузионный поток
- •2.3. Баланс массы
- •2.4. Уравнение сохранения количества движения
- •2.5. Уравнение баланса энергии
- •2.6. Обобщённое дифференциальное уравнение
- •2.7. Выбор координат
- •Односторонние и двухсторонние координаты
- •Глава 3. Методы дискретизации
- •Сущность численных методов
- •Р ис. 3.2. Изменение относительной температуры по длине трубы
- •Методы получения дискретных аналогов
- •Метод формальной замены производной конечно-разностными отношениями. Этот метод основан на разложении в ряд Тейлора достаточно гладких функций. Рассмотрим узловые точки, показанные на рис. 3.3.
- •Запишем, например, следующее разложение по формуле Тейлора
- •Анализ методов аппроксимации. Вернёмся к модельной задаче, сформулированной выше, и выявим связь точности численного решения со способом аппроксимации.
- •Основные правила построения дискретных аналогов
- •Заключение
- •Глава 4. Численное моделирование процессов теплопроводности
- •4.1. Стационарная одномерная теплопроводность
- •Нестационарная одномерная теплопроводность
- •Двух- и трёхмерные задачи
- •4.4. Некоторые геометрические соображения
- •4.5. Заключение
- •Глава 5. Конвекция и диффузия
- •5.1. Рассматриваемая задача
- •5.2. Установившиеся одномерные конвекция и диффузия
- •5 .3. Дискретный аналог для двухмерных задач
- •5.4. Дискретный аналог для трёхмерных задач
- •5.5. Односторонняя пространственная координата
- •5.6. Схемная искусственная диффузия
- •5.7. Заключение
- •Глава 6. Расчёт поля течения
- •6.1. Обоснование необходимости в специальной методике
- •6.2. Трудности расчёта поля давления
- •6.3. Переменные "завихренность – функция тока"
- •6.4. Шахматная сетка
- •6.5. Поправки скорости и давления
- •6.6. Уравнение для поправки давления
- •6.7. Модифицированный алгоритм simpler
- •6.8. Заключение
- •Оглавление
- •Глава 1. Дифференциальные уравнения в частных производных...................................12
- •Глава 2. Математическое описание физических процессов............................................28
- •Глава 3. Методы дискретизации.........................................................................................49
- •Глава 4. Численное моделирование процессов теплопроводности.................................72
- •Глава 5. Конвекция и диффузия..........................................................................................109
- •Глава 6. Расчёт поля течения................................................................................................134
- •Библиографический список
1.4. Эллиптические дифференциальные уравнения в частных производных
В теплофизике (и не только в этой предметной области) эллиптические ДУЧП ассоциируются со стационарными проблемами. Простейший пример эллиптического ДУЧП даёт уравнение Лапласа
(1.32)
которое определяет стационарное распределение температур в длинной заготовке прямоугольного поперечного сечения, потенциальное течение несжимаемой жидкости, фильтрацию грунтовых вод в почве и многое другое. Кроме того, слагаемые левой части уравнения (1.32) входят составной частью в уравнения других типов, например, движения газов.
Пусть (x, y) характеризует значение температуры в сечении квадратной заготовки (в данном случае у – пространственная координата). Поскольку уравнение (1.32) имеет второй порядок по обеим координатам, то для его решения необходимо задать четыре граничных условия. Зададим на нижней грани (см. рис. 1.8) распределение температуры, соответствующее зависимости (х, 0) = sin(x), а на верхней грани (x, 1) = sin(x)exp( ). На левой и правой гранях зададим нулевые температуры: (0, у) = (1, у) = 0. При этих граничных условиях уравнение (1.32) допускает точное решение
(x, y) = sin(x)exp( y) (1.33)
в области 0 x 1, 0 y 1.
Для
эллиптического ДУЧП второго порядка,
представленного в форме уравнения
(1.2), существует важный принцип
максимума,
сущность которого состоит в утверждении:
как
максимальное, так и минимальное значения
достигаются
на границе Г, если
не считать того тривиального случая,
когда величина
постоянна. Этот принцип весьма важен
для проверки правильности решения.
Как уже отмечалось ранее, эллиптические ДУЧП не имеют вещественных характеристик, что, естественно, должно наложить отпечаток на физические свойства решения.
Интерпретация на физической основе
Наиболее важная особенность, касающаяся эллиптических ДУЧП, состоит в том, что возмущение, внесённое в точке Р (см. рис. 1.8), оказывает влияние на все другие точки области исследования, хотя вдалеке от точки Р это влияние будет малым. Тем самым косвенно утверждается, что при попытке отыскания численных решений эллиптических задач необходимо рассматривать глобальную область, т.е. всю область исследования сразу. В противоположность этому параболические и гиперболические ДУЧП можно решать прогрессивно-маршевым методом, отталкиваясь от начальных условий. Разрывы в граничных условиях для эллиптических ДУЧП сглаживаются внутри области.
Отметим ещё одно обстоятельство. При рассмотрении параболических и гиперболических ДУЧП путём применения разложения в ряд (например, Тейлора) всегда можно получить некое локальное решение в непосредственной близости от границы или от некоторой точки, относительно которой у исследователя имеется подробная информация о закономерностях изменения искомого параметра. Попытка проделать то же самое, имея дело с эллиптическими ДУЧП, приводит, как правило, к получению бесконечно возрастающего решения, т.е. локального решения эллиптических задач не существует в принципе; имеются лишь глобальные решения.
