Кто-нибудь проводил расчет газожидкостного равновесия (ГЖР) для газоконденсатной смеси по уравнению состояния и компонентному составу? Там надо вычислять константы равновесия в итерационном цикле, при этом надо на каждой итерации решать итерационным способом уравнение Рашфорда-Райса или Рэкфорда-Райса (Rachford and Rice) для нахождения молярной доли газовой фазы. В книгах пишут что нужно применять итерационный метод Ньютона или Ньютона-Рафсона. А на деле у меня этот метод иногда расходится. Применял также метод деления пополам (бисекции) - то же самое. Просто на каком-нибудь шаге итерационного цикла получаются такие значения констант равновесия что уравнение Рашфорда-Райса не имеет корней и асимптотически приближается к нулю и имеет точку разрыва. Посоветуйте пожалуйста как правильно построить алгоритм.
Кто-нибудь проводил расчет газожидкостного равновесия (ГЖР) для газоконденсатной смеси по уравнению состояния и компонентному составу?
- Давно уже написан софт (не один) по этому поводу. Странно, что кто-то еще занимается такими расчетами. Если не секрет для кого "программу" пишите?
Для себя, а какие программы есть, ссылки приведите пож-та если можно.