- •Моделирование физических процессов Учебное пособие
- •Содержание
- •3. 8. Математическая модель массообмена ……………………………. 54
- •5.1.1. Формулы конечного дифференцирования …………………… 74
- •5.6.2. Метод Ньютона ………………………………………………. 95
- •Введение
- •VII этап
- •2. Структурные математические модели
- •Основы структурного моделирования
- •2.2. Типовые элементы
- •2.2.1. Способы математического описания звеньев и систем.
- •2.2.2. Типовые динамические звенья.
- •2.3. Схемы соединения типовых элементов
- •2.4. Структурные модели теплообменного оборудования
- •2.4.1. Составление математических моделей тепловых объектов.
- •3. Физические математические модели
- •3.1 Использование математических моделей
- •Под углом 57 к горизонту со скоростью 20 м/с
- •Изменение температуры при остывании тела
- •3.2 Математическая модель теплопроводности
- •. Математическая модель гидродинамики
- •3.4 Математическая модель теплообменника
- •3.5 Математическая модель парогенератора
- •С естественной циркуляцией
- •Математическая модель горения
- •Математическая модель образования токсичных
- •3. 8. Математическая модель массообмена
- •Дифференциальные уравнения диффузионного пограничного слоя
- •Моделирование теплоотдачи диффузией.
- •3.9. Математическая модель парового котла
- •3.9.1. Математическая модель динамических процессов.
- •3.9.2. Регулирование давления пара и тепловой нагрузки
- •3.10. Пример построения математической модели объекта
- •4. Планирование эксперимента
- •Полный трехфакторный эксперимент
- •Численные методы
- •5.1. Интерполяция
- •5.1.1. Формулы конечного дифференцирования
- •5.1.2. Сплайны
- •Способы получения краевых условий:
- •5.2. Решение систем алгебраических линейных уравнений
- •Метод прогонки. Метод прогонки используется для решения систем специального вида
- •5.3. Решение систем алгебраических нелинейных уравнений
- •5.3.1. Метод итераций
- •5.3.2. Метод Ньютона
- •5.4. Решение дифференциальных уравнений
- •5.4.1. Постановка задачи. Разностные схемы.
- •5.4.2. Разностные схемы. Метод Рунге – Кутта.
- •Пример решения дифференциального уравнения методом Рунге-Кутта.
- •Результаты решения уравнения методом Рунге-Кутта
- •5.5. Решение краевых задач методом прогонки
- •5.5.2. Распространение тепла в пластине радиатора.
- •Сравнение метода прогонки с точным решением
- •5.6. Решение краевых задач методом пристрелки
- •5.6.1. Метод пристрелки.
- •5.6.2. Метод Ньютона. Рассмотрим граничную задачу, определяемую дифференциальным уравнением второго порядка
- •5.7. Решение краевых задач методом конечных разностей
- •5.7.1. Метод конечных разностей.
- •5.7.2. Линейные дифференциальные уравнения.
- •Моделирующие программы
- •Смеси газов в реакторе
- •Перемещение свободной поверхности жидкости
- •Скоростей течения испаряющегося метана
- •Рекомендованный библиографический список
- •Пасько Петр Иванович «моделирование физических процессов»
- •346360, Г. Волгодонск, ул. Ленина 73/94
Метод прогонки. Метод прогонки используется для решения систем специального вида
содержащих n+1 уравнение с неизвестными z0, z1, ..., z n .
Построим матрицу системы:
Основная система называется трехточечной разностной схемой, а первое и последнее уравнения – краевыми (граничными) условиями.
Если выполняются следующие условия:
и ,
то находим неизвестные:
,
,
где ; ; i = n – 1, n – 2, ..., 2, 1.
Произвольные системы линейных уравнений легче всего решать методами исключения с применением обращения матриц. Итерационные методы эффективны лишь в том случае, если матрицы содержат много нулевых элементов.
5.3. Решение систем алгебраических нелинейных уравнений
Решение нелинейных уравнений. Нахождение корней уравнения. Рассмотрим нелинейное уравнение , где – непрерывная функция на отрезке [а, b]. Прежде чем вычислять корни с необходимой точностью, нужно или определить их границы, или найти приближенные значения корней. Одним из самых простых методов является метод половинного деления, который применим, если , т.е. функция на концах интервала имеет разные знаки. Разделим отрезок [а, b] пополам и возьмем ту половину, на концах которой имеет разные знаки. В свою очередь разделим эту половину пополам и поступим далее аналогично. В результате получаем последовательность вложенных друг в друга отрезков. Процесс продолжают, пока не будет выполнено условие , где – точность, с которой нужно найти корень. Середина отрезка [аn, bn] дает нужное приближение к корню. Сходимость к корню в этом методе обязательна и от требуется лишь непрерывность.
5.3.1. Метод итераций
Рассмотрим уравнение вида на отрезке [a; b] при условии, что данная функция дифференцируема на [a; b]. Тогда на этом отрезке функция имеет единственный корень, который может быть найден как предел последовательности приближений:
,
где x0 – начальное приближение; .
Погрешность метода простых итераций:
,
где все xk [x0; x0+r]; – константа (0 < < 1). Сходимость имеет место со скоростью геометрической прогрессии со знаменателем меньше 1.
Геометрический смысл метода итераций заключается в отыскании точки пересечения графиков функций y = x и y = (x). Процесс сходимости показан на рис. 5.2,а если функция убывает (-1 < < 0), и на рис. 5.2,б если функция возрастает (0 < <1).
а) б)
Рис. 5.2. Метод итераций
Условие 0 < < 1 является достаточным, но не необходимым: если > 1, то итерации могут разойтись, а уравнение может иметь решение. (Например функция (x) = bx, где b > 1. Итерации разойдутся, если в качестве начального приближения взять любое x0 0.)
Метод итераций можно применять и для решения нелинейных систем уравнений. Рассмотрим пространство Rn с введенной на этом пространстве нормой x. Замкнутым шаром S(x(0), r) с центром в точке x(0) и радиусом r называется множество всех n-мерных векторов x таких, что расстояние между x и Y не превосходит r: .
Рассмотрим нелинейное уравнение , где – вектор-функция, а х – n-мерный вектор, т.е. рассмотрим систему:
Определим матрицу .
Если функция непрерывно дифференцируема в шаре S и в шаре S выполняются условия: и , то для решения нелинейных систем уравнений можно применять метод итераций.