关于fluent重力失效的问题

yekai339

我在做的是蒸发器内的单相水流场分析，计算后发现流场压降和基本估计的不一致：流场中存在一段竖直圆形截面通道，长约6m，液相水从上往下流动。估算该段的摩擦压降和形阻压降之和不超过10kpa，6m的高度能产生约合58kpa，所以估算的进出口总压差应该为-48kpa。但是计算结果显示，进出口压降为1kpa左右。我初步估计是重力没有起作用，但是我有勾选重力选项，并在相应坐标上设定了重力加速度。

为了找出问题所在，我建立了一个简单竖直圆形截面管道进行研究。

简单竖直管道模型描述：管道长5m，管径0.5m，流体常物性，密度677.9034 kg/m3，粘度80 upa*s。估算得第一层网格高度0.8-4.8mm，建立了第一层高度为2mm的边界层网格，全局非结构化网格，质量0.3以上。X方向重力加速度9.81m/s2（方向是对的），使用standard k-e模型，standard-wall function。进口速度边界，速度0.176839m/s，湍流强度2.7052%，水力直径1m。出口压力边界，压力为0，回流湍流强度2.7052%，水力直径1m。SIMPLE算法，均采用二阶离散格式。

对简化模型进行了计算，发现进出口总压差还是和估计的严重不符。Ps，整个计算最终没有错误和警告提醒，连续性和速度残差都是十的负六次方，k和e残差十的负四次方，进出口质量守恒，进口总压稳定。查看一篇国外的CFD-Online的帖子，（帖子地址：http://www.cfd-online.com/Forums ... lation-results.html）里面的楼主遇到了和我一样的问题，他最后的解决办法是在operating condition里设置operating density为0，而且他认为这是fluent中的bug（我觉得这对fluent来说简直难以置信）。

我也进行了同样尝试，发现进出口的压差确实和估计的比较一致了。但是计算结果显示进出口的总压都是负数，尽管SIMPLE只考虑压力的相对值，但是我不知道这样做是否正确。不知道论坛里有没有哪位前辈有遇到过类似的问题，还望能多多指教。

ps.在论坛里深水区发过这个帖子了，但是好像看到的人比较少，就在这再发一遍了。