В работе рассматриваются вопросы численного моделирования двухфазных течений с применением эйлеровой многоскоростной модели типа Баера-Нунциато. Представлено описание математической модели, подробно рассмотрен численный алгоритм решения задачи на основе разрывного метода Галеркина. Представлено описание разработанного программного комплекса, основное назначение которого - математическое моделирование двухфазных течений с прямым разрешением динамики границ раздела фаз. Особенностью предложенных алгоритмов является лимитирование простых и консервативных переменных, гарантирующее как отсутствие нефизичных осцилляций, так и допустимые значения физических полей с применением лимитера WENO-S. Основная цель работы - дать исчерпывающее описание предложенного комплекса алгоритмов для решения задач рассматриваемого класса и пригодного для решения задач в реалистичных постановках. На примере модельной задачи демонстрируется возможность его применения для решения задач на сетках больших размеров.