@bestucan 看了官方的，主要就是怕需要一些依赖的包。。安装一下试试，看看新功能爽一下

何止有啊，简直一大堆，搜搜 OpenQBMM。:chigua:

@胖小帅 速度是不是要设置成slip?你试试

@东岳 李老师。我是要在壁面附近一层网格区域加源项，每个网格需要加的质量计算的公式，需要乘以壁面这个网格的面积，再除以网格的体积，请问我需要怎么写呀

@bestucan 感谢前辈的回复，目前问题已经解决:xinxin:

在编译OpenFOAM之前，需要先编译ThirdParty，但OpenFOAM WiKi中的Ubuntu18.04版安装教程里面好像并没有“编译ThirdParty这一步”，希望以后大家注意吧

@东岳 :xinxin: 谢李老师！

@液固两相流 我是在计算域内初始化一定量的颗粒，输出计算过程中每个颗粒的位置和速度信息，这个可以实现吗？

@东岳 谢谢李老师！:ok:

箭头是对指针的成员的引用，找到这个类的定义：turbulenceModel.H

epsilon()函数声明在205行，是个虚函数，所以在 turbulenceModel.C 里没有定义，而是在使用其作为模板创建新类的时候再各自定义

在上面那个窗口里搜 epsilon() 函数：
Screenshot from 2020-09-15 16-37-07.png

看你的情况是用的哪种定义，其实里面 epsilon() 函数定义大差不差，
你用的 LES，咋找上层流的了，你找的那个确实是层流的对于 epsilon() 的一个实现，在我发的那个图片里也有

脉动速度都用上了，那应该不是层流了，速度分成脉动速度和平均速度这就是用了湍流模型了

最近看了以下jasak的论文，找到了相关描述
屏幕截图 2020-09-14 160918.jpg

@random_ran 这位大佬，我也尝试了您说的这个方法，是不是这样产生的网格都是大小一致的，没办法局部加密网格？

感谢分享: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的。

@闻久STU 谢谢哦，您这个方法我也记在小本本上了，以后肯定能用上:xiezuoye:

简单验证了一下，可以用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

@cccrrryyy 受教受教，授我以渔，万分感谢~ 帖子链接里开篇就阐明了我的疑惑，也难怪k-epsilon可以在工程领域运用的如此广泛~ 我就先对照着openfoam tutorials和CFDEM案例里的k-epsilon取值作参考，计算后再做调整，再次感谢指点~

@东岳 非常感谢老师:xinxin: ！ :xiexie:
就是这样迂回的话好像有些难，我还要再思考思考，再尝试尝试 :144:

@东岳 好滴

Already solved.

@bestucan Thanks for your attention.

谢谢老师

@东岳 这个当时手一抖就点上了，发现没法删掉，真是尴尬:mihu:

@莫晓柒 编译之后得到的库和求解器程序都在 platforms 下面，所以只能重新安装了。

@东岳 我不算运营人员，只是参与了号主的一些活动，有时候帮帮忙。。。

@upc_ngh
1.是的，$\mathbf{U}$代表的是KinematicParcel中的$\mathbf{U_c}$，也就是连续相的在颗粒所处位置的速度插值，具体的插值方式OF里面有cell，cellpoint那些，$\mathbf{U}_\rm{d}$表示的就是颗粒的速度。
2.是的，系数乘以速度差。

新知识点

@沐风 因为是百度云，所以没有下载:xinlei: 不过我有几个小建议：
计算域是否够大，圆柱和上下前后边界的距离是否够远，去除边界的影响；
snGrad建议用corrected；
我印象中上下边界应该用symmetry，用noslip是有什么原因或者依据吗。

@xpqiu 感谢

@chien 是的谢谢🙏，放弃编译了😹

@小考拉 maxU会突然增大只是发散的迹象，不是发散的原因，所以这种问题很难debug :135:

我也是每次都重新设置，很麻烦，前来学习

min(alpha.water) 都出现负值了。

时间步自己调整的特别小是为了让库朗数降低，结果时间步那么小，库朗数还那么大，可能网格有问题

尝试解释一下:mihu:

就是0.3图里的涡在0.4不见了，然后1里又出来了？

0.3里刚过柱体，柱体后面压力小，然后左右流体往中间挤，挤不动了，要么往上要么往下，往下是墙，只能往上冒，就是0.3里的涡

0.4里不见了，因为这个柱体下面有个空隙，流体在过空隙的时候很挤，过了空隙就释放压力，下边是墙，是能往上释放，然后空隙里出来的流体就把0.3里的涡往上托，托到液面托没了就成0.4了

然后0.6，为了释放压力往上托涡，但又因为惯性，释放过头了，涡把流体都往上带，下面就出现低压区了，然后左右的流体都往中间流，在这里碰头，可以从0.6看见他们碰头后很纠结，不知道往上走还是往下走

然后是1，往上的那部分变成了0.3的涡尾，往下的那部分，又卷出来个涡。

要是水足够深，应该能同时看见0.3和1的涡共存。

再形象点，就是0.3里的涡是柱体后的低压区造成。1里的涡就像柱体下的狭缝放的“屁”，一个屁把0.3里的涡挤走了……:yes: