Моделирование многофазного потока

Добрый день! В кратце необходимо смоделировать течение газа через слой частиц для получения распределения скоростей и давлений в сечении. Особенность модели в том что слой выглядит в виде "пирога" с разными размерами твердых частиц. Сам слой частиц можно считать условно неподвижным. Проблема заключается в не физичных результатах. Дело в том что в нестационарной постановке задачи слои смешиваются, чего в реальном аппарате не происходит. Стандартное (ну по крайней мере описанное в туториале) решение задачи в Эйлеровой гранулярной постановке не учитывает модели вязкостного трения (Frictional Viscosity), которые для моделей с высокими концентрациями твердой фазы преобладающие. Но решение при включении моделей Schaeffer либо Johnson приводит к разваливанию решения на первых итерациях. Подобная ситуация наблюдается и при использовании опции Packed Bed, хотя я и выполняю требования по фиксации скоростей в каждой зоне. Какие есть ограничения по многофазной постановке в моем варианте. Может есть какие-то особенности инициализации или задания начальных значений?