ПОИСК Статьи Рисунки Таблицы Особенности алгоритмов расчета из "Моделирование промышленных процессов полимеризации" Применительно к реакторам идеального смешения при использовании типовых моделей статические характеристики определяют по табл. где приведены формулы вычисления основных характеристик X в случае микросмешения. Уравнения решают последовательно для каждого к-го реактора (от к=1 до k=N) по следующей схеме. Находят концентрации мономера из уравнения п-й степени (п=2, или 3, или 5, в зависимости от типа обрыва), затем находят концентрации активных цепей единичной длины и моментов активных цепей в порядке увеличения их номера (нулевой, первый, второй...). Далее находят моменты неактивных -цепей в порядке еличения их номера. Размерность вектора основных координат X применительно к типовым модулям ограничена восемью, но при необходимости число координат может быть увеличено. [c.135] Этот вопрос ранее специально не обсуждался, поскольку он в большей степени относится к моделированию производства в целом как сложной химико-технологической схемы. Ограничимся здесь кратким рассмотрением. Условимся различать следующие модули межреакторных связей и аппаратов соединительный, смесительный, разделительный. [c.139] Такова роль соединительных трубопроводов, промежуточных емкостей. Их не учитывают при анализе стационарных режимов, однако если нас интересуют капитальные затраты, их учет при моделировании обязателен. При анализе нестационарных режимов наличие соединительных блоков также не может не учитываться. [c.139] Равенство (1У.4) записано с учетом аддитивности основных параметров моделей X в таком виде, который логически следует из самого способа их введения (моменты не нормированы и не центрированы). [c.139] Следует подчеркнуть, что такая проверка выполнения ограничений с разбивкой расчета на два этапа (расчет статических характеристик материальных потоков и проверка физической реализуемости теплового режима) значительно облегчает и ускоряет расчет по сравнению с одновременным расчетом по уравнениям материального и теплового баланса (о способах решения которого будет сказано несколько позже, после рассмотрения аналогичных характеристик установившегося режима для полностью сегрегированных систем). [c.140] Как следует из уравнений (1У.2), вычисление характеристик каскадно-реакторных схем для условий макросмешения нельзя свести к решению простых алгебраических уравнений с готовыми аналитическими решениями. Искомые характеристики получают последовательным численным интегрированием, причем процесс решения разбивают на несколько этапов, охарактеризованных ниже. [c.140] В качестве примера приведем схему метода Рунге — Кутта (наиболее часто используемого для этих целей) при вычислении по одной из формул, погрешность которой пропорциональна Л (/I — шаг интегрирования) на каждом шаге. Ограничимся здесь рассмотрением системы второго порядка. [c.140] Аналогичный вид имеют формулы для систем более высокого порядка, какими чаще всего и являются наши модели. [c.141] Анализ совместного поведения функций / ( , 0) и X t) показал, что пределы интегрирования уравнения (1У.2) с достаточной точностью могут быть установлены от О до 40, что существенно экономит затраты машинного времени при вычислении. [c.141] Сходимость метода Ньютона сильно зависит от хорошего начального приближения, поэтому иногда при расчете прибегают к разнообразным вариантам итеративных расчетов с запоминанием результатов предыдущих итераций, среди которых наиболее известен метод Вольфа [20]. [c.142] Вернуться к основной статье