这是一个非常常见的现象就是求解的量越界了。CFD里面，尤其是有限体积法，需要满足若干限制。主要的有：守恒、有界。守恒就是开始1kg的水，后来还是1kg。有界就是比如温度不能出现负值。你这个算例，没有温度下降以及上升的进出口以及源项，计算后低于300k就是越界了。你用的什么数值格式？你全部用Upwind试试。

你是正版Fluent用户么，联系ANSYS咨询没有，这种有界性boundedness的CFD问题很偏理论，JCP上一堆文章在处理这个问题，不知道ANSYS Fluent那面怎么解决的。

@sunwinter 平台上除了fluent,还有CFX,star-ccm+等，另外还开发了openfoam线上仿真，提供界面

液体射流喷射雾化后，在高温高压的环境中燃烧

不是很熟悉你这个物理过程。普遍性原则是你们是否需要某一内部坐标点的数量密度函数，如果需要就需要重构。比如蒸发过程就需要当内部坐标为0时候的权重，因此需要NDF重构。

感觉跟充气机翼挺像的，，，

谢谢您的回答，fluent里不是只有压力和密度求解器吗，问题有些浅显，不要见怪。还有怎么设置流体域为真空。

@李东岳 哦，明白，谢谢李老师

我的模型是计算一个泵的外特性以及泵的内部流场，这个泵由一个叶轮加一个导叶组成，叶轮叶片数为3，导叶叶片数为4，现在为了加快计算速度，采用了周期性边界，选取了一片叶轮加一片导叶来计算，所以这个流量的大小应该怎么设置？

@灬Xin 16.1,我将同样的case发给两个朋友，分别是15.0,16.0，都没问题。然而,诡异的是，当我个人电脑和服务器(linux)换成15.0,上述情况仍在。

@深海矿工 非常感谢您的回答，我在Fluent中已设定每个时间步都保存整支叶片zone的压力数据，后面再处理，缺点是数据会占较大的空间

是不是udf里面没有保存运动的数据

@下里巴人 流线图知道怎么画了，也能画出整体的矢量图，可是对于同一水平线的矢量图还不知道怎么画，能讲下么，谢谢

@下里巴人 最近稍微试了一下均匀流下粗糙网格的3组试验，一组时间步长给定为旋转2度的时间，一组为0.2度，一组Fluent自调整。因为时间原因仅做了一小段时间内的对比（未达到稳定周期性变化），关注的两个物理量（一个轴向力，一个转矩）在三组试验中趋势几乎一致，数值误差小于3%.另外对一种非均匀流下的精细网格进行了同样的试验（同样时间较短），趋势几乎一致，但不如均匀流下趋势符合得那么好，误差在9%以内。
受时间限制，对比不严谨，仅作为前期考察

参考的文献是网上的“Parallel EDEM CFD Coupling for ANSYS FLUENT - version 1.0 - User Guide”，以及是用的网上下载的耦合包，所使用的EDEM版本为EDEM2.6，以及FLUENT中的耦合接口出现了。

@李东岳 除了phase的图像，其他压力图，密度图等都没有问题。0_1504141697200_ea5fe613-9544-43ae-acce-f76b5a6f726e-image.png 这是密度图。我在网上有人提到会不会是由于相与相之间的界限不明显，而产生这个问题。老师，请问您觉得这种说法可靠么？谢谢

比如调成2的话，这个时候的时间步长是多少？

你的时间步的2倍。

根据这幅图的说法，隐式耦合计算时，非稳态的courant我也可以设置的比1大吗？

隐式理论上可以使用任意大的时间步，因此没有库狼数限制。所以可以。

Gosman无处不在

应该可以写udf输出吧

@李东岳 非常感谢。
又读了几篇类似的文献，真是用时均处理的仿真模型找不到，高雷诺数下的这种流动都是在用实验测试了，但是我需要计算一个比较准确的结果，合作单位的同志希望能找到统一的规律，形成经验公式类似的结论，然后用作产品开发:crying: 任重道远啊，有个师兄推荐用DES试试，目前在测试，希望别掉进坑了。

周期性边界条件主要是这俩个边界的通量一正一负，速度相同。

对称性边界主要是此边界通量为0，边界法向速度为0，只存在切向方向，对称性边界对于标量等于零法向梯度边界。

我用openfoam也出现这个问题，有种说法是不可压假设导致的，我还没试过