Моделирование ротационного вискозиметра

Добрый день!

Помогите, пожалуйста, решить задачу для университета. 

Суть задачи состоить в моделировании работы ротационного вискозиметра. Принцип его работы заключается в измерении динамической вязкости жидкости через измерение момента сопротивления вращению цилиндра в этой жидкости. Я моделирую обратную задачу - через известную динамическую вязкость найти момент сопротивления вращению. 

Задача следующая: есть внешний цилинд диаметром 1000 мм и высотой 1000 мм, заполненный жидкостью с известной динамической вязкостью 40 Пас. В него погружается цилиндр диаметром 400 мм на глубину 360 мм. Внутренний цилиндр вращается с постоянной угловой скоростью 300 об/мин. Необходимо определить момент сопротивления и распределение скоростей. 

Я решаю задачу так: моделирую цилиндр диаметром 1000 мм и высотой 360 мм с цилиндрическим отверстием 400 мм.

Строю сетку и назначаю слудующие граничные условия:

1. Домен fluid -> без вращения -> ламинарная модель течения жидкости

1. Inlet - Mass Flow Rate = 0.00001 kg/s (тк осевого движения жидкости нет)

2. Outlet - Average Static Presure = 0 Pa

3. Вращение врутренней стенки - Wall -> Boundary Details -> Option = No Slip Wall -> Wall velocity -> Rotating wall = 0 rev min^-1

4. Нужная стенка статическая - 

При решении сходимость обеспечивается, кроме P-Mass Imbalance, которая +-100.

По итогу распределение скоростей получается верным, а расчетный момент в 2 раза больше теоретического (ansys = 460 Hm ,расчет по формулам 227 Hm).

Не понимаю в чем может быть ошибка. Подскажите, пожалуйста, для меня это очень важно)