- •Моделирование физических процессов Учебное пособие
- •Содержание
- •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. Метод пристрелки.
Итерационные методы решения граничных задач (как линейных, так и нелинейных) могут быть разделены на две большие группы, а именно методы пристрелки и методы конечных разностей.
В метоле пристрелки выбирается недостающее (незаданное) начальное условие в начальной точке интервала, а потом решается задача Коши – численно интегрируется дифференциальное уравнение до конечной точки интервала. Затем, сравнивая вычисленное значение зависимой переменной в конечной точке с заданным, проверяют правильность выбора недостающего начального условия. Если эти значения не совпадают, то выбирают другое значение недостающего начального условия, и процесс повторяется до тех пор, пока разность между вычисленным и заданным условиями в конечной точке не станет достаточно малой.
Наиболее часто используются три таких метода: метод Ньютона, метод параллельной пристрелки и метод квазилинеаризации. Все эти методы по существу основаны на одном и том же принципе – ньютоновском методе решения нелинейных алгебраических уравнений – и поэтому в равной мере обладают двумя важными свойствами: монотонной сходимостью и квадратичной сходимостью. Для большинства задач методы Ньютона и квазилинеаризации одинаково эффективны. Но из-за того что в методе Ньютона подгоняется лишь недостающее начальное условие, в то время как в методе квазилинеаризации систематически подгоняются значения функции во всех точках интервала, можно ожидать, что последний обладает лучшей сходимостью. Метод параллельной пристрелки обычно применяется в тех задачах, где решение весьма существенно зависит от выбора производной в начальной точке.