ПОИСК Статьи Рисунки Таблицы Расчет температуры кипения из "Машинный расчет парожидкостного равновесия многокомпонентных смесей" Для уяснения методики расчета обратимся к блок-схеме, представленной на рис. V-1. Вначале вводятся исходные данные (подпрограмма INPUT) и задается начальное приближение по температуре для первой итерации. Несмотря на неточность этого приближения, оно совершенно не сказывается на конечных результатах расчета. Коэффициенты фугитивности компонентов паровой фазы полагаются равными единице (для первой итерации) вводятся заданные параметры — давление и состав жидкой фазы. [c.59] Вновь полученная величина суммы концентраций SUMY сравнивается с предыдущим значением. Если сумма изменяется, то опять вызывается подпрограмма PHIMIX для расчета коэффициентов фугитивности, и так до тех пор, пока при двух последующих итерациях значения сумм концентраций не будут постоянными (с заданной точностью), но не обязательно равными единице. [c.61] Как только будет получено постоянное значение величины SUMY, последняя сравнивается с единицей. Если сумма концентраций равна единице в пределах допустимой ошибки, то полученная температура отвечает условиям равновесия в противном случае вновь корректируется температура и выполняется еще одна итерация по температуре. Последующие уточнения температуры производятся по способу Ньютона ., В программе предусмотрен соответствующий контроль сходимости решения. [c.61] Детальное описание программы дано в главе VII. Примеры расчета приведены в Приложении. [c.62] Вернуться к основной статье