- •Моделирование физических процессов Учебное пособие
- •Содержание
- •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
5.5.2. Распространение тепла в пластине радиатора.
Рассмотрим
пластину радиатора при температуре
окружающего воздуха, равной
(задача аналогична рассмотренной в п.
3.2). Можно считать, что температура в
поперечном сечении, перпендикулярном
оси пластины (оси x),
постоянна. Температура Ts
у левого края пластины (при х = 0) известна.
Уравнение теплопроводности можно записать в виде
,
где
а – коэффициент температуропроводности,
;
S
– площадь поперечного сечения; Р –
периметр поперечного сечения;
– коэффициент конвективной теплопередачи;
– коэффициент теплопроводности.
Граничные
условия при x
= 0 Т = Ts,
а при x
= L
перенос тепла обусловлен только
конвекцией и, следовательно,
.
Если
ввести безразмерные переменные
и
,
то дифференциальное уравнение и граничные
для него условия запишутся следующим
образом:
,
,
,
где
,
а
.
Для
решения этого уравнения используем
метод прогонки и введем уравнения
и
.
Тогда уравнения для прямой прогонки запишутся в виде
и
,
а
граничные условия –
и
.
Граничные значения при x = L получим в виде
и
.
Проинтегрировав
от
до
,
найдём
и
.
Результаты для нескольких значений
и
приведены в табл. 5.2, в которую также
включены значения
,
найденные по точной формуле.
Таблица 5.2
Сравнение метода прогонки с точным решением
|
|
|
Метод прогонки |
Точное решение |
|
|
|
|
|
||
|
2,0 |
0,5 |
0,3494 |
0,1747 |
0,3494 |
|
1,5 |
0,2364 |
0,3546 |
0,2364 |
|
|
3,0 |
0,1592 |
0,4775 |
0,1592 |
|
|
1,0 |
0,5 |
0,4693 |
0,2347 |
0,4693 |
|
1,5 |
0,3025 |
0,4537 |
0,3025 |
|
|
3,0 |
0,1973 |
0,5919 |
0,1973 |
|
Известно, что точное решение этой граничной задачи даётся формулой
.
Откуда
при
получается
.
На
этом заканчивается этап прямой прогонки.
Чтобы получить распределение температуры,
нужно использовать обратную прогонку,
приняв за начальные условия
и
из таблицы.
5.6. Решение краевых задач методом пристрелки
5.6.1. Метод пристрелки.
Итерационные методы решения граничных задач (как линейных, так и нелинейных) могут быть разделены на две большие группы, а именно методы пристрелки и методы конечных разностей.
В метоле пристрелки выбирается недостающее (незаданное) начальное условие в начальной точке интервала, а потом решается задача Коши – численно интегрируется дифференциальное уравнение до конечной точки интервала. Затем, сравнивая вычисленное значение зависимой переменной в конечной точке с заданным, проверяют правильность выбора недостающего начального условия. Если эти значения не совпадают, то выбирают другое значение недостающего начального условия, и процесс повторяется до тех пор, пока разность между вычисленным и заданным условиями в конечной точке не станет достаточно малой.
Наиболее часто используются три таких метода: метод Ньютона, метод параллельной пристрелки и метод квазилинеаризации. Все эти методы по существу основаны на одном и том же принципе – ньютоновском методе решения нелинейных алгебраических уравнений – и поэтому в равной мере обладают двумя важными свойствами: монотонной сходимостью и квадратичной сходимостью. Для большинства задач методы Ньютона и квазилинеаризации одинаково эффективны. Но из-за того что в методе Ньютона подгоняется лишь недостающее начальное условие, в то время как в методе квазилинеаризации систематически подгоняются значения функции во всех точках интервала, можно ожидать, что последний обладает лучшей сходимостью. Метод параллельной пристрелки обычно применяется в тех задачах, где решение весьма существенно зависит от выбора производной в начальной точке.
