感谢分享:xiexie:

谢谢李老师啦，我先看一下这几个求解器，了解一下思路:xinxin:

@东岳 ok，明白了，谢谢东岳老师

@Samuel-Tu
記得好像是用巨集 #calc

@benqing 会不会是因为固体域没有时间步，所以网格不发生变化？:xiezuoye:

今天在编译liggghts和cfdem的时候也出现如下报错：
185e1dd1-371e-43a4-9c79-1be84f630d3e-image.png

结果证明需要如下依赖包：

sudo apt-get install libvtk6-dev

另外好像OF8没法完成和liggghts的耦合，我用的是5.x的。

简单验证了一下，可以用shm，没问题:xinxin:

@walten 在 如何CFD自学成才 中说：

@闻久STU
乐观点， 六年:high: 这样才学得深入

不能提:135: 长期抗战是肯定的了

抱歉问的问题可能太泛了，具体举个例子来说。比如说在直径D=0.01m的空气管流中，速度Uf=15m/s，雷诺数Re=10000。颗粒的斯托克斯数St=10，直径是大概dp=50×10^(-6)m，密度为ρp=1000kg/m^3，颗粒-空气的质量比是Φ=0.1，双向耦合，只考虑颗粒曳力。

在这种情况下，我查看OpenFOAM计算得到的kinamaticCloude:UTrans文件中，动量转移都是在10^(-14)量级（单位是kg×m/s），请问这合理吗，为什么只有这么小。。？

据我所知，这个UTrans文件表征的就是颗粒→流体相的动量转移，也就是颗粒相通过源项Sp，施加在流体相的曳力作用。不知道这个理解对不对？

最近多次使用foamDictinary感觉真是太用了，自己在文件里面加的keys也能识别出来并修改

各位大佬问一下有关拉格朗日颗粒喷射的问题，我如何往喷射模型中（比如manualInjection模型）动态的传递喷射点坐标，坐标每一个时间步都会发生变化，怎么实现内部传递呢。

@zhe 在 请问，MPPIC这个方法的particle distribution function的运作 中说：

固相是由particle distribution function表示的，并且固相的相分数也是通过它进行求解的。那么这个方程在OpenFOAM里有体现吗？

没有，你说的这个方程是一个非常基本的控制方程，其实就是玻尔兹曼方程，很多模型都是从这个模型上展开对他进行封闭的。比如

双流体模型，就是粒子分布函数传输方程取2次矩的模型。

矩方法就是取高阶矩的模型。

MPPIC在这个方程上对碰撞进行BGK封闭。可以更简单的处理碰撞。看这个文献的方程2

A model for collisional exchangeing as /liquid/solidfluidizedbeds

《无痛苦N-S方程笔记》4,8节也有讨论。这个方程没人去解。自由度太多了。

可以彼此融会贯通 左右逢源一下 :144:

@马乔 有没有什么办法输出这个误差呢，一般大家是怎么判断pisoFoam或者pimpleFoam是否收敛呢?
我发现误差降到１e-5以下就很困难了。。。所以我是设置的监测点，发现关注的变量变化不大了就认为收敛。。

你这是什么版本啊，为嘛那一通运算没有量纲检查。forAll这个只会遍历内场。

看论文发现有作者说fluent稳态计算残差降到1e-7则认为收敛停止计算，但是我用of的simpleFoam的时候发现残差根本降不到那么低，连降到１e-６都已经是极限了。我都不敢在论文里写残差判断标准了，一般你们怎么说收敛标准呢

@bestucan 好滴，我试一下

@xpqiu 好滴

@Samuel-Tu 用了，确实快，但是有时候分网不成功，而且网格质量很难保证

本来想说试试加深水槽。 静水流深嘛，一看是浅水方程……
如果边界条件设定为带反射的，那就应该不停波动。如果是自由边界，不停波动不正常

在cfd online 里逛了一下，有些收货。
似乎non-orthogonal corrector 是在组装矩阵Ax=b的时候用的，它会把扩散项中的正交部分放入系数矩阵Ａ中，而把非正交部分放入方程右边的源项b中，使得系数矩阵更稀疏且主对角线控制，使得计算更稳定？？cfd online