Проверка расчета с помощью функции rkfixed
Проведем расчет профиля концентраций в реакторе для схемы реакции, приведенной в примере:
Реактор описывается системой (152):
Преобразуем (152) к виду:
Решим систему (157) с помощью встроенной в Mathcad функции rkfixed, которая использует метод Рунге-Кутта. 1. Зададим исходные данные:
2. Рассчитаем количество точек на интервале решения:
3. Зададим начальные условия:
4. Введём правую часть системы дифференциальных уравнений:
Расчеты в функции rkfixed производятся с использованием двумерных массивов. Чтобы функция работала правильно, нужно ввести замену в системе (157): XA заменить на x0, XP – на x1, XD – на x2. В данном случае элементы массивов x0, x1 набираются с использованием кнопки Matrix®Subscript (Матрицы®Нижний индекс) или горячих клавиш Shift+] (или Shift+ъ при русской раскладке клавиатуры). 5. При вызове функции возврат рассчитанных значений производится в массив Z:
6. Зададим ранжированную переменную i, рассчитаем сумму концентраций, переобозначим массивы, выведем результаты. Построим графики, отражающие изменение концентраций веществ по длине реактора (рис. 100).
Рис. 100. Профиль концентраций в реакторе идеального вытеснения
|
.
(157)
