@东岳 抱歉，应该是我想当然了:mihu: ，如果有延迟的话，那中间的流体就被压缩到一起了，我用的是不可压流体求解器。

（1）应该可以在自己的用户文件夹中进行编译，etc/bashrc里面的安装路径最好改一下，免得和原来的冲突；
（2）icc也可以成功编译openfoam，之前用gcc编译openfoam3.0的时候，gcc编译openfoam的源码编译不成功（可能是gcc版本太低吧），就用icc替代了gcc；
（3）感觉主要注意一下环境变量的配置。

Ok 东岳老师，我问题解决了，谢谢您

PRESSURE GRADIENT TERM这个图中的$F=\frac{m}{V} \nabla p$，其中的$V$应该为$\rho$

const fvPatchList& patches = mesh.boundary(); forAll(patches, patch) { Info << patches[patch].name() << nl; const fvPatch& currPatch = patches[patch]; forAll(currPatch, face) { Info << abs_[0].boundaryField()[patch][face].Cf() << nl; } }

可是气体的扩散系数很小，靠扩散远远达不到10mg/m3的浓度. 之前试着用toposet设置了浓度为定值(10mg/m3)的面源，但扩散效果还是不明显。

你这个更像是设计问题。如果物理上靠扩散远远达不到10mg/m3的浓度，那应该靠对流的作用。具体效果如何，还取决于你的几何。只能建议更改几何看看在什么情况下效果最好。看起来像拓扑优化了。

不过是实验，不是模拟。

@东岳 现在并没有解决，我正打算找一个2维的算例跑一下试一试，谢谢老师！

@东岳 老师，DEM这块不是我编写的，工作主要是在CFD部分，据我了解DEM是计算单个颗粒的曳力，不用指定气泡直径大小吧；那我再看一下CFD部分曳力计算吧，谢谢您

大佬解决这个问题了吗？

如果有用，很高兴你们能引用这篇文章 Simulation of bubbly flows with special numerical treatments of the semi-conservative and fully conservative two-fluid model

@东岳 东岳老师你是不懂我（们？也许就我一个人这样干哈）这些底层的科研民工，比起完全自己创作我更愿意在一些现有的资源上进行修改

你的网格尺寸太大了，缩小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

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