@dzw05 嗯，谢谢你，这种方法的确比较简单，但我更想用上面那种方法来实现，因为之后会加非稳态风场等等，后续会比较方便一点:xinlei:

@dzw05 谢谢！！

静候大作！

用旧版本的paraview可以，我现在用新版的也有这个问题。paraFoam确实慢。
没有解决办法，只能心中自有云图咯 :mihu:

耦合求解和分离求解的结果有什么不同呢？

俩种求解策略。高超那面比较用户耦合求解。不可压缩低速这面SIMPLE PISO近似垄断。这不好说哪个结果好哪个结果不好。你看SCI，发表出来的，结果都很好。算法上的区别，耦合求解构建一个大矩阵，压力速度放在一起求解，一锅粥出来了。分离求解分开构建矩阵，比如压力方程、速度方程。

@东岳 谢谢李老师！

@linhan-ge 改过来也没有发生变化，据说难以保证气泡的形状精准。

附上温度限制：

temperatureLimit { type limitTemperature; active true; limitTemperatureCoeffs { selectionMode all; min <Tmin>; max <Tmax>; } }

div改成upwind有所改善，但不能根本上解决问题。
如果温度使用uniformFixedValue将压力从小到大随时间增大，为何实际算的时候也不起作用？
初始压力极小如绝对压力为1，即使温度也很低，则能正常计算，压力和温度同事增大或减小才可以计算，这是什么原因？

谢谢各位指教

@东岳 李老师好，目前发散问题已解决，但是由于epsilon会偏离正常值范围，所以会导致流场混乱，且层流算出来的流场和湍流k-e模型算出来的基本相似，尝试用低雷诺数模型SA计算，则无法观察到流场中的涡。

对于场，IOobject类型里面如果是AUTO_WRITE就是输出。对于颗粒，目前不确定是不是IOobject ..

@东岳 谢谢东岳老师，本来我的solver是基于OpenFOAM4.0版本，根据您的建议，我换用了OpenFOAM-extend4.0来编写solver，成功通过编译。
如您所说，错误确实很有可能是因为粘弹性的库不匹配OF的库函数。
谢谢东岳老师的提示。

Zalesak disk可以参考isoAdvector里面的算例配置：
https://github.com/isoAdvector/isoAdvector/tree/master/OpenFOAM/run/prescribedU/notchedDiscInSolidBodyRotation

@winsway_zero这个问题解决了吗？

在某些湍流模型中，需要计算指向壁面的面矢量，OpenFOAM提供几种计算方法，但这些方法还需要进行详细的测试。

比如对于下面的2D管道，下图是采用meshWave方法计算的指向壁面的面矢量，可见是正确的。

0_1535517133002_correct.jpg

下图是采用Poisson方法计算的指向壁面的面矢量，管道中心的指向壁面的面矢量是错误的。

0_1535517141915_noCorrect.jpg

需要上床，不能粘贴~

各位老师，你们好。我正在学习openfoam，我在重复tutorials里面的第二个例子stress of a plate with a hole时，细化网格后，用singleGraph把（0 0.5 0.25）到（0 2.0 0.25）这条线上的正应力提取出来时，发现只提取了0.5这个点的值。我在paraFoam用plot over line看到0.5到2.0这条线上的值都是存在的，为什么用singleGraph提取会出现这种问题呢

想过将物体的几何模型分成立柱和底下的浮筒，这样可以设置立柱和浮筒的单独的边界条件，这样就可以监测力，但是物体是整体运动的，所以要用createPatch将立柱和浮筒组合成平台的边界条件，这样立柱和浮筒的边界条件又没有了

0_1533986941754_2018-08-11 19-28-30屏幕截图.png 貌似有答案了。。。。

但求解结果文件中并没有发现创建好的通量场

你强制写一下看看？field.write();