附上温度限制：

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();

多谢，已更新，

应该把继续把去掉。刚看了下后面还没改，等我回来之后改一下。明天坐飞机去喀纳斯

@昝磊 在solid2中，你通过读取DT文件注册了volScalarField类型的DTsolid2，所以objectRegistry中并没有DT。你可以试试将DT文件改为DTsolid2。

@东岳 据说速度没有原生的ubuntu快，但是我没有亲自测试，window界面不习惯，不过能算也行了

starFrom:latestTime

有个问题想请教一下大家。我使用interDyMFoam得到数据后，想获得沿某一条直线上的相分布。但是导入paraview之后使用plot over line得到的是一根绝对坐标系下静止直线上的数据，事实上，我想得到的是在运动的物体上固定的一条线的数据。不知道有没有方法可以让这根线跟着物体运动？或者有没有方法可以在计算的过程中直接输出沿这条线上的数据？
感谢大家！

@cyw 你好，请问你是在哪个版本中添加的传热模型呢？
网上只找到了of230的方法。
我看过模拟加热气泡的文献，晚几天把名字发给你可以参考下。
话说，方便细聊么？我最近也在弄这方面的东西
vx：zhigu4610

@vamoszhi 您的声波扰动方程添加成功了吗？我也准备将声学扰动方程加入OF，请问是否方便交流下