[求助]高手请留步

D调的华丽

案例描述：搅拌桶中的二次流现象的模拟。此处采用直搅拌浆而并非叶轮，流体在搅拌浆的作用下被甩向边壁，此时流体具有射流特性，在边壁上分上下两路运动，形成一对涡旋。在平面来看可看见搅拌浆的两侧各有一对涡旋。
问题：开始读入网格文件的时候并没提示体积为负，但是当运行一段时间之后突然冒出体积为负，不知道是什么原因；另外，不知道如何去定义角速度，因此此处采用的是二维网格，搅拌浆是绕着X轴旋转。
希望高手能帮我想想这个问题，感激不尽……

D调的华丽

这么多人，难道就没有一个高手，全部是菜鸟吗？？？
晕！