Работа посвящена математическому и численному моделированию течения среды в ступени центробежного насоса. Поставленная задача представляет собой сложную физико-математическую проблему, т.к. изучается течение в трехмерном канале, где имеются повороты потока на большие углы и вращающиеся элементы. В работе сформулирована начально-краевая задача для уравнений в частных производных (уравнений Навье-Стокса, осредненных по Рейнольдсу), описывающих 3-х мерное нестационарное вязкое турбулентное течение среды в ступени насоса. В качестве модели турбулентности использовалась модель. В качестве численного метода использован конечно-разностный метод второго порядка точности. Проведенные расчеты позволили получить детальное распределение параметров в ступени и оценить возможные критические ситуации в работе устройства. Результаты численного моделирования проанализированы, получены их интегральные характеристики и сопоставлены с имеющимися экспериментальными данными и упрощенными аналитическими моделями. Сравнение с экспериментом показало хорошее качественное и количественное соответствие результатов математического моделирования и эксперимента.