Гомогенные равновесная и неравновесная (релаксационная) модели используются для численного моделирования вскипающих потоков в соплах. Расчеты выполнены с помощью разработанного авторами CFD-кода ANES. Тестирование реализованной математической модели и модифицированных для рассматриваемого класса задач алгоритмов кода ANES проведено с привлечением известных экспериментальных данных. Представлены результаты тестовых расчетов, а также данные, полученные применительно к режимным параметрам и конструкции сопла и расширителя системы парогенератор-сепаратор теплоутилизационного комплекса ЗАО НПВП "Турбокон". Показана работоспособность алгоритма SIMPLE на задачах трансзвукового и сверхзвукового истечения вскипающей жидкости. Продемонстрировано, что учет неравновесности с помощью релаксационной модели позволяет получить лучшее согласие с экспериментальными данными по распределению доли паровой фазы вдоль оси сопла, при этом различие 1D- и 2D- моделей для рассматриваемого класса течений невелико.