В данной работе для моделирования процессов размыва связного грунта и распространения поверхностных волн используется нестационарная неоднородная система уравнений Навье—Стокса с переменными вязкостью и плотностью. Значение плотности определяется при помощи уравнения конвекции-диффузии. Для решения полученной системы уравнений используется алгоритм, состоящий из схемы расщепления по физическим факторам для системы уравнений Навье—Стокса и метода предиктора–корректора для уравнения переноса. Система решается на разнесенной сетке методом сеток. Представлены результаты двух- и трехмерных расчетов.