Предложен неявный численный метод, основанный на совместном интегрировании уравнений газовой фазы и фракций частиц, который позволяет проводить расчеты сверхзвукового двухфазного течения при больших весовых долях конденсата и малых размерах частиц. Задача рассмотрена в рамках многожидкостной модели сплошной среды с лагранжевым представлением процессов взаимодействия частиц друг с другом.