请教各位老师，有没有研究涉及OpenFOAM与PFC3D的耦合的。老板让用PFC3D研究DEM-CFD coupling，感觉很少有研究这个的，目前看的资料也就是github上一个itasca公司的老师编写的一系列代码和一篇今年6月份发表的文献，资料有点少，也不知道这是否是一个正确的方向，还望各位老师指点指点，谢谢各位老师！

@mohui 是的，interface处是有间断的，但是one-field形式的控制方程体现不出来间断
压力分布不对，看一下理论解就明白了～

不是少部分，是一点点负值，比如说-0.00006。如果你的方程有一些要对浓度做开方，那一点点负值都不行。

@李东岳 CFDEMcoupling-PUBLIC-3.0.x用这个，耦合的是OpenFoam中的PisoFoam。createPatch这个是ICEM中的，还是OpenFoam中的？

我从OpenFOAMl里面导出了算例网格和数据，用Fluent进行后处理时，找不到浓度、含率参数，这是什么原因？？？？！！！。其中C和alpha.water对应Fluent参数序号如下图。0_1510822319990_QQ图片20171116165029.jpg

好的，我试试

0_1510557038187_1510484187969-捕获-resized.jpg

想用二相欧拉模型做一个固液二相流问题，下层是固体（体积分数0.6），上层是水，入口处水给一个固定流量，入口的边界条件该怎么给呢？示意图：
0_1510474822975_1.PNG

@李东岳 感谢东岳老师，非常详细的步骤。我会好好看看这块

@李东岳 您好！感谢您的提议。首先rhoCentralFoam应该是可以算湍流的，我是基于这个前提，想知道除了在constant文件下turbulenceModel选择RAS或LES以后，如何改动其system文件夹下的fvsolution和fvscheme字典以及在0文件夹下新加相应的变量字典来计算湍流。不知道我的这种想法是不是正确的。

东岳大神说了开发的问题。我来说说我的看法？清楚物理模型的物理意义貌似是最基本的也是最容易出错的部分。openfoam方面可以参考The Finite Volume Method in Computational Fluid Dynamics。还有就是最笨的办法：读代码

@asininno 我没有详细对比过具体有多大影响。关于非正交修正，建议你仔细看一下公式，《The Finite Volume Method in Computational Fluid Dynamics An Advanced Introduction with OpenFOAM and Matlab》这本书有很详细的说明。说白了就是采用和正交网格上一样的方法去做离散，对于非正交网格会有误差，那么再加上一些措施进行修正。
但是还是建议做好网格，控制网格的正交性和弯曲度；毕竟修正只是一种补救措施。

@dzw05
哦了，有空整理整理。

好吧，你说的很对，看来就是：
\begin{equation}
\frac{\nabla p}{\rho}=\frac{\mathrm{D} \mathbf{U}}{\mathrm{D} t}
\end{equation}

of中多相流这块不可压好像都没有能量方程，在模仿着可压的能量方程添加时发现，剪切力做功是用湍流粘性为形参的库函数表示的，现在我不调用库函数，怎么表示这一项呢？谢谢

@wwzhao
谢谢赵老师的回答，我尝试下写入不同的值看是否能正常写入文档

有关Eo数的定义：http://47.93.231.95/doku.php?id=cfd百科#eoetvoes_number

谁那见过三维Eo数的计算公式和文献出处？OpenFOAM植入的这个公式不知道从哪来的。

@李东岳 额，从controlDict里的时间来看，多相计算的物理时间要远大于单相，但还是没能得到合理结果。我再算下去看看。

@dzw05 好的，真的非常感谢:happy: :happy:

@李东岳 嗯，调分辨率就可以了。东岳老师我最近在用二相欧拉模型求解一个水沙运动问题，但是算到一定时间就会蹦掉，您有空帮我看下我的边界条件设置是不是有问题吗？