![](/user_photo/2706_HbeT2.jpg)
- •Федеральное агентство по образованию
- •2404000000-35 Удк 681.142:519.6
- •Оглавление
- •1. Основы работы в Mathcad 10
- •2. Роль численных методов 36
- •3. Методы аппроксимации и интерполирования 42
- •4. Лабораторная работа № 1. Интерполирование степенными многочленами 44
- •5. Лабораторная работа № 2. Параметрическая идентификация математических моделей методами аппроксимации 65
- •6. Решение систем линейных алгебраических уравнений 93
- •7. Лабораторная работа № 3. Решение систем линейных уравнений приближенными методами 109
- •8. Лабораторная работа № 4. Решение нелинейных уравнений приближенными методами 142
- •9. Лабораторная работа № 5. Решение систем нелинейных уравнений приближенными методами 174
- •10. Лабораторная работа № 6. Решение обыкновенных дифференциальных уравнений первого порядка численными методами 197
- •11. Лабораторная работа № 7. Численное интегрирование 212
- •12. Лабораторная работа № 8. Моделирование реактора идеального вытеснения для многостадийной химической реакции с линейной кинетикой 226
- •13. Лабораторная работа № 9. Расчет моделей процессов диффузии и теплопроводности с помощью явной разностной схемы 249
- •Предисловие
- •1. Основы работы в Mathcad
- •1.1. Панели инструментов
- •1.2. Ввод и вывод данных
- •1.3. Осуществление несложных вычислений
- •1.4. Построение и настройка графиков
- •1.5. Программирование в Mathcad
- •1.5.1. Программирование без программирования
- •1.5.2. Язык программирования Mathcad
- •1.5.3. Создание программы (Add Line)
- •1.5.4. Редактирование программы
- •1.5.5. Локальное присваивание ()
- •1.5.6. Условные операторы (if, otherwise)
- •1.5.7. Операторы цикла (for, while, break, continue)
- •1.5.8. Возврат значения (return)
- •1.5.9. Перехват ошибок (on error)
- •1.5.10. Примеры программирования
- •2. Роль численных методов
- •2.1. Этапы решения задачи на компьютере
- •2.2. Математические модели
- •2.3. Численные методы
- •3. Методы аппроксимации и интерполирования
- •4. Лабораторная работа № 1. Интерполирование степенными многочленами
- •4.1. Постановка задачи
- •4.2. Порядок выполнения работы
- •4.3. Краткие теоретические сведения
- •4.3.1. Метод неопределенных коэффициентов
- •4.3.2. Интерполяционный многочлен Лагранжа
- •4.3.3. Интерполяционные формулы Ньютона для равностоящих узлов
- •4.4. Примеры выполнения
- •4.4.1. Интерполирование степенными многочленами с использованием метода неопределенных коэффициентов
- •4.4.2. Интерполирование степенными многочленами с использованием второй интерполяционной формулы Ньютона
- •4.5. Требования к отчету
- •4.6. Контрольные вопросы и задания
- •4.7. Задания
- •5. Лабораторная работа № 2. Параметрическая идентификация математических моделей методами аппроксимации
- •5.1. Постановка задачи
- •5.2. Порядок выполнения работы
- •5.3. Краткие теоретические сведения
- •5.3.1. Метод выбранных точек
- •5.3.2. Метод средних
- •5.3.3. Метод наименьших квадратов
- •5.4. Примеры выполнения
- •5.4.1. Аппроксимация с использованием метода выбранных точек
- •5.4.2. Аппроксимация с использованием метода средних
- •5.4.3. Аппроксимация с использованием метода наименьших квадратов
- •5.4.4. Сравнительный анализ методов аппроксимации
- •5.5. Требования к отчету
- •5.6. Контрольные вопросы и задания
- •5.7. Задания
- •6. Решение систем линейных алгебраических уравнений
- •6.1. Общие положения
- •6.2. Точные методы решения систем линейных уравнений
- •6.2.1. Метод Крамера
- •6.2.2. Метод Гаусса
- •6.2.3. Метод обращения матриц
- •7. Лабораторная работа № 3. Решение систем линейных уравнений приближенными методами
- •7.1. Постановка задачи
- •7.2. Порядок выполнения работы
- •7.3. Краткие теоретические сведения
- •7.3.1. Математическое описание реактора идеального смешения непрерывного действия
- •7.3.2. Математическое описание кинетических закономерностей химических превращений
- •7.3.4. Приближенные методы решения систем линейных уравнений
- •7.3.4.1. Метод простых итераций
- •7.3.4.2. Метод Зейделя
- •7.4. Примеры выполнения
- •7.4.1. Пример выполнения задания точным методом
- •7.4.2. Пример выполнения задания методом итераций и методом Зейделя
- •7.5. Требования к отчету
- •7.6. Контрольные вопросы и задания
- •7.7. Задания
- •8. Лабораторная работа № 4. Решение нелинейных уравнений приближенными методами
- •8.1. Постановка задачи
- •8.2. Порядок выполнения работы
- •8.3. Краткие теоретические сведения
- •8.3.1. Этапы решения нелинейного уравнения
- •8.3.4. Метод деления отрезка пополам (вилки, дихотомии)
- •8.3.5. Метод Ньютона (метод касательных)
- •8.3.6. Метод простых итераций
- •8.4. Пример выполнения задания методом итераций
- •8.5. Требования к отчету
- •8.6. Контрольные вопросы и задания
- •8.7. Задания
- •9. Лабораторная работа № 5. Решение систем нелинейных уравнений приближенными методами
- •9.1. Постановка задачи
- •9.2. Порядок выполнения работы
- •9.3. Краткие теоретические сведения
- •9.3.1. Метод Ньютона
- •9.3.2. Метод итераций
- •9.4. Примеры выполнения
- •9.4.1. Метод Ньютона
- •9.4.2. Метод итераций
- •9.5. Требования к отчету
- •9.6. Контрольные вопросы и задания
- •9.7. Задания
- •10. Лабораторная работа № 6. Решение обыкновенных дифференциальных уравнений первого порядка численными методами
- •10. 1. Постановка задачи
- •10.2. Порядок выполнения работы
- •10.3. Краткие теоретические сведения
- •10.3.1. Метод Эйлера
- •10.3.2. Модифицированный метод Эйлера
- •10.3.3. Метод Эйлера-Коши
- •10.3.4. Метод Рунге-Кутта 4-го порядка
- •10.4. Примеры выполнения
- •10.4.1. Реализация метода Эйлера в математическом редактореMathcad
- •10.4.2. Решение обыкновенного дифференциального уравнения 1-го порядка с помощью функции rkfixed
- •10.5. Требования к отчету
- •10.6. Контрольные вопросы и задания
- •10.7. Задания
- •11. Лабораторная работа № 7. Численное интегрирование
- •11. 1. Постановка задачи
- •11.2. Порядок выполнения работы
- •11.3. Краткие теоретические сведения
- •11.3. 1. Метод прямоугольников
- •11.3.2. Метод трапеций
- •11.4. Пример выполнения
- •11.5. Требования к отчету
- •11.6. Контрольные задания
- •1 Таблица 101.7. Задания
- •12. Лабораторная работа № 8. Моделирование реактора идеального вытеснения для многостадийной химической реакции с линейной кинетикой
- •12. 1. Постановка задачи
- •12.2. Порядок выполнения работы
- •12.3. Краткие теоретические сведения
- •12.3.1. Математическая модель реактора идеального вытеснения
- •12.3.2. Численное решение систем дифференциальных уравнений
- •12.4. Пример выполнения
- •12.5. Проверка расчета с помощью функции rkfixed
- •12.6. Требования к отчету
- •12.7. Контрольные вопросы и задания
- •1 Таблица 112.8. Задания
- •13. Лабораторная работа № 9. Расчет моделей процессов диффузии и теплопроводности с помощью явной разностной схемы
- •13.3.2. Решение уравнений в частных производных
- •13.3.3.Метод сетки
- •13.3.4. Явная разностная схема
- •13.3.5. Условия устойчивости явной разностной схемы
- •13.4. Пример выполнения
- •13.5. Требования к отчету
- •13.6. Контрольные вопросы и задания
- •13.7. Задания
- •Библиографический список
- •Использование
7. Лабораторная работа № 3. Решение систем линейных уравнений приближенными методами
Цель работы:
Изучение приближенных методов решения систем линейных уравнений, а также получение навыков реализации и исследование эффективности работы этих методов на ЭВМ.
Получение практических навыков разработки математической модели технологического объекта на примере моделирования реактора идеального смешения, использование результатов моделирования для исследования процессов, протекающих в реакторах указанного типа непрерывного действия.
7.1. Постановка задачи
Для химического реактора идеального смешения непрерывного действия, функционирующего в изотермическом стационарном режиме, разработать математическую модель.
Используя методы простых итераций и Зейделя, а также один из точных методов (согласно варианту), рассчитать концентрации реагентов на выходе реактора и исследовать зависимость последних от времени пребывания веществ в зоне реакции.
Осуществить сравнительный анализ эффективности использованных методов.
Вид химической реакции, протекающей в реакторе, константы скоростей отдельных стадий, концентрации исходных реагентов и точные методы решения систем линейных уравнений получить согласно варианту (пункт 7.7, табл. 7). Время пребывания изменяется в диапазоне 1<<10 (с шагом 1), степень точности получения решения численными методами =10-3.
7.2. Порядок выполнения работы
Изучить методические указания и ответить на контрольные вопросы.
Получить у преподавателя номер варианта.
Разработать математическую формулировку реактора идеального смешения непрерывного действия, работающего в изотермическом режиме (составить уравнения материального баланса).
Осуществить расчет концентраций реагентов на выходе реактора как функций времени пребывания. Для этого осуществить решение системы линейных уравнений.
Осуществить решение системы точным методом. Выбрать метод в соответствии с вариантом.
Преобразовать систему в нормальный вид. Осуществить проверку сходимости.
Осуществить решение системы, используя методы простых итераций и Зейделя.
Составить схемы алгоритмов и написать программы в среде Mathcadили Си++ (в помощь см. пункт 1.4).
Отладить программы и получить результаты расчетов.
Результаты расчетов представить графически.
Провести анализ полученных результатов (на предмет соответствия хода полученных кривых виду химической реакции, а также сравнения точности и скорости сходимости использованных при выполнении задания методов).
Оформить отчет.
7.3. Краткие теоретические сведения
7.3.1. Математическое описание реактора идеального смешения непрерывного действия
Реактор смешения представляет собой чаще всего цилиндрическую емкость со сферическим или коническим днищем, снабженную мешалкой (рис. 42).
Рис. 42. Внешний вид реактора
Внутри может происходить химическая реакция. Часть аппарата покрыта рубашкой, в которую поступает теплоноситель с целью отвода теплоты от содержимого или, наоборот, его нагрева.
Под идеальностью смешения понимается отсутствие внутри реактора градиента концентраций и температур, т. е. в любой точке в объеме реактора в конкретный момент времени концентрация вещества или его температура одинаковы.
Реактор идеального смешения может работать в нестационарном и стационарном режимах.
Математическое описание реактора смешения можно получить, исходя из уравнений модели идеального смешения, если учесть скорость образования продуктов в реакционной зоне. Если положить, что в процессе химического превращения число молей реагентов не изменяется, то можно записать:
(35)
где V –объем
реактора смешения, м3;,
– вектор концентраций продуктов на
входе и выходе из реактора соответственно,
доли;t– время, ч;– скорость подачи продукта в зону
реакции, м3/ч;W – скорость
образования вещества (продукта), ч-1.
При наличии теплового эффекта реакции и теплообмена с внешним носителем необходимо использовать соотношение, определяющее изменение температуры в зоне реакции. Уравнение теплового баланса модели идеального смешения представляется в следующем виде:
(36)
где ср– удельная теплоемкость реагирующей смеси, Дж/(м3С);Т0– температура реагентов на входе,C;T– температура реагентов на выходе, C;Тх– температура теплоносителя,C;Q– тепловой эффект реакции, Дж/(м3ч);kт – коэффициент теплопередачи, Дж/(м2чС);F– площадь поверхности теплообмена, м2.
Уравнения (35) и (36) представляют собой математическое описание реактора идеального смешения в случае нестационарных режимов. При совместном решении этих уравнений можно получить графики зависимости концентрации реагентов и температуры в реакторе от времени в нестационарных условиях.
Стационарные условия работы реактора смешения достигаются, если концентрации веществ на входе реактора постоянны во времени. Тогда концентрации веществ внутри и на выходе аппарата также постоянны во времени Процесс является изотермическим, если температура внутри реактора T=const. Тогда можно записать:
. (37)
В этом случае уравнения (35) и (36) преобразуются к виду:
, (38)
, (39)
где
– среднее время пребывания реагентов
в аппарате, ч.
Задачей, возникающей при расчете реакторов смешения, работающих в стационарном режиме, является определение концентраций веществ на выходе реактора для заданных концентраций на входе и времени пребывания. В случае линейной кинетики (для реакций первого и псевдопервого порядка) и при изотермическом режиме работы реактора (T=const) соотношения (38) представляют собой систему линейных уравнений. В результате решения полученной системы уравнений будут рассчитаны концентрации реагентов на выходе.