|
楼主 |
发表于 2016-6-22 10:07:22
|
显示全部楼层
太感谢您了,您的解答对我很有启发。这也和我之前的看法比较一致。动量方程里面的压力梯度项里的压力指的是静压,help里面说到勾选重力后除了引入boussinesq假设外都会使用操作密度ρ0,那重力项就变成了(ρ-ρ0)g。而在fluent中又会进行压力的重新定义,P'=P-ρ0g 。当流体稳定时,P'=0的。 也就是说fluent中即使勾选了重力,但是由于fluent中对压力的重新定义,使得求解得到的静压和没有重力时是一样的。而如果需要确定真实的静压,只需要在计算得到的静压的基础上加上重位压力就行了。
fluent这么处理让人有点摸不到头脑,但是我觉得也并不是没有原因的:加入重力后,如果不对压力进行重新定义,那么使用SIMPLE算法时,会发现在求p'时存在的常数项b里面会加入重力的影响。而对于单相水来说,由于其密度较大,可能其对b的贡献会比其他项大很多,以此产生了数值误差?所以fluent为了保证数值模拟的精确度,在计算静压时采取了那样的处理。
我尝试过简单管道勾选重力但是将operating density改为0的情况,计算结果虽然进出口的压力差和估算的差不多,但是流体从进口进入流域后存在不合理的速度上升,计算已经收敛,但是出现了不合理的流场,不就应该是发生了数值假扩散?
不知道您对我的分析有没有什么看法?
|
|