你的网格尺寸太大了，缩小1/10或者1/100试试

@chpjz0391 嗯 是的 。 其实是一样的 换汤不换药

有研究界面扩散模型的吗？大家可以相互交流交流哈:chigua:
qq:805759762

//2.4s label meshSize = this->thermo().rho()().size(); for(label i=0; i < this->nSpecie_; i++) { for(label icell = 0; icell < meshSize; icell++) { label izone = ZoneNumber[icell]; this->RR_[i][icell] = RRCCM[i][izone]*correctCH[icell]; } } //24s for(label i=0; i < this->nSpecie_; i++) { for(label icell = 0; icell < this->thermo().rho()().size(); icell++) { label izone = ZoneNumber[icell]; this->RR_[i][icell] = RRCCM[i][izone]*correctCH[icell]; } }

上边两种用法，计算耗时相差 10 倍！
改用 thermo.T() 之后计算时间又变成 2.4 s 了，所以原因应该是 rho 后边跟了两个括号，深层原因是啥呢？

//2.4s label meshSize = this->thermo().rho()().size(); for(label i=0; i < this->nSpecie_; i++) { for(label icell = 0; icell < this->thermo().T().size(); icell++) { label izone = ZoneNumber[icell]; this->RR_[i][icell] = RRCCM[i][izone]*correctCH[icell]; } }

此外，还可以使用 forAll(this->thermo().T(), icell)，效果和 for 循环几乎一样。

之前一直使用interFoam，求等温的液滴动力学问题。现在需要加入温度场求解温度用interTempFoam，手头只有of230版本的资料，跨版本编译of3和of1706编译通过，但是计算报错，，，

网格是非等间距的，of230不支持，不用blockMesh字典直接复制ployMesh后计算也报错，进退维谷这就很尴尬了:zoule:

有大神能分享下经验或者提供一下帮助吗？:chouchou:

昨晚看东岳大神说OF6发布了，，，可能会有类似的求非等温两相流求解器？

@aiweimo 非常感谢。

重新运行了算例, 提高了内存请求(200GB->251GB)。同样的结果。请求内存似乎不是问题所在，峰值 RSS 和 200G内存一样。

JobID Elapsed NCPUS AveRSS MaxRSS AveVMSize ReqMem ------------ ---------- ---------- ---------- ---------- ---------- ---------- 14389 06:00:15 192 251Gn 14389.batch 06:00:15 48 302K 194452190K 113468K 251Gn 14389.extern 06:00:15 192 115K 144K 107948K 251Gn

@东岳 谢谢东岳老师
之前看命令提示试过加括号了
但是没想到引号
非常感谢！

还没有专门看过这方面的研究

康师傅，你也在

@东岳 “......在计算域内部添加一个吹气孔进行吹气”
现在是要实现类似的吹气过程，将网格作为温度、质量（化学组分）源项，且网格赋予固定的速度以及压力，求解结果中反查这个网格的速度值和设定值是有偏差的（怀疑和组分浓度扩散相关，扩散系数存在导致网格有一定的附加“扩散速度”），结果的压力值误差更大（大于100%，猜测openfoam中应该用特殊的方式去限定这个压力约束）。

......BTW真正关心的问题来了，在fireFoam application源程序中没有找到质量运移的扩散系数D，也没找到施密特数Sc的设定位置（施密特数（Schmidt number, Sc）是一个无量纲的标量，定义为动黏滞系数和扩散系数的比值，用来描述同时有动量扩散及质量扩散的流体。施密特数的命名是为了纪念德国工程师 Ernst Heinrich Wilhelm Schmidt (1892-1975)。）

而且对于稍微复杂一点的结构，边界层经常出错。反正人家重心也不在这上面

这个问题我之前也琢磨过，

加vDotcAlpha是因为后面这个项原型应该为（1-alpha）*vDotcAlpha

其实在湍流模型的k方程里面也有类似的处理，我个人偏向于以下这个结论

尽量对各项做隐性处理，这样可以保证方程的稳定性，

在压力方程里面，mass transfer rate 这一项也是做隐性处理的哦

这些辐射系数都是和 wavenumber 相关, 在RT E中需要这些带量纲参数.
Thanks

@kingkong 您好，请问openfoam有自带的根据马赫数是否大于1选择的这种边界条件吗？

有人能讲解下interface composition的saturated和henry 这两个model吗。。到底怎么设置thermophysical 的文件？

Thank you. That is a great help.

看了一下，新版本的VOF算法被提出来了，更加的紧凑，尽可能的减少代码复用。在这种情况下，不同的算法具有不同的显性源项和隐性源项。在interFoam中，源项均为0，因此均初始化为zeroField。对于可压缩VOF，源项则初始化为其他，可参考compressibleInterFoam中的alphaSuSp.H。

@东岳
谢谢！用Surface Vectors filter就直接画出来了

@mr-h http://www.cfd-china.com/topic/1199 看这个，通过forAll操作边界