@mohui 我是基于那个方程进行修改添加一些，只不过它本身存在一些问题，当我设置1的时候温度正常

可以，so？

@李东岳 前辈，我一直用pyFoam来观察残差，我猜的这个bounded variable plot出现是因为计算出现了不稳定，如果这个bounded量计算没过界就不会出现。:confused:

另外观测压力等物理量的残差如下，这个就比较好理解了:cheeky:

0_1508035857600_WechatIMG15956.jpeg

@我也不知道吧 好的，我加你。

我认为这个force应该是用来维持湍流的并不是产生湍流，否则没有这一项的话湍流会被粘性耗散掉。如果初始的速度场没有用boxTurb产生随机场的话，而是均匀的速度为零的速度场的话，根本算不出湍流。

这个我还不确定，你可以不用boxturb，然后用dnsFoam算一下他那个tutorial看看有没有湍流出来。

想问一下这一段代码是如何定义force的？

目前只知道基于快速傅里叶变换，更深的代码还没看 :cheeky:

gv只能打开一个eps

@李东岳 除了添加力，还需要求解界面上表活剂输运方程，具体的介绍（James, Ashley J, Lowengrub, et al. A surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant[J]. Journal of Computational Physics, 2004, 201(2):685-722. ）。怎么求解输运方程这一块我不知道怎么处理。

http://www.cfd-china.com/topic/905/几个有意思的研究方向

@random_ran 这个词我一直在用，哈哈，逆压梯度

@李东岳 就是考虑粘性耗散项，粘性耗散热也是热源

类似这种求解不正确的问题很难判断 :confused:

这个看起来需要在你的局部区域做一个网格，然后mapField过去。

linux命令。同时将输出显示在屏幕上，你百度一下就明白了

@李东岳 这个装在一个cluster上，好久没维护了，不是很懂cluster不敢乱弄，怕万一弄崩了，没人会rebuild了。。。

不是point，有一个手动选取cell的按钮..

@wwzhao 好的，解决了

自己挖的坑自己填回去。。。。。。。。。。。。。

我今天发现是我自己定义的TELF函数出了问题，所以计算出来的结果有问题。
但是这种方法加入状态方程我实践过的确是可行的。
就是换了以后，连续性误差的地方会变的非常怪，查了continuityErr这个文档，发现里面计算连续性误差还是使用rho-thermo.rho所以会出现误差值十分大，
目前还在想着这么去改thermo下的东西，因为自己不能把黏性系数和普朗数设成固定值，
怎样把这些thermo文档下设置函数还在研究中。

@李东岳 多谢李老师！

@李东岳 谢谢大神的解答:happy:

需将2进制网格文件转换为ascii码格式

在controlDict下把format改为ascii，然后运行foamFormatConvert就可以转换了。