@Jacob-Z 在 CFDEM安装 中说： @fubianhanshu CFDEM®coupling 3.8.0=LIGGGHTS® 3.8.0+OF 5.x 好的，谢谢！我也尝试下，以后一起交流

@CFDngu 我不做波浪所以不清楚一般大家都怎么做，所以我觉得看一下湍流粘度比值会有帮助。RANS计算的话，对一般的流动而言这个比值应该较大，如果比值很小那说明其实没有太多湍流的作用。也许真的像你们说的，其实用层流就没问题。 网格我是想问你有没有做过无关性分析之类的，另外这种场里面突然有一个小点突然算出来很大的量以至于发散，有可能是网格正好碰到了什么奇怪的数值点。当然究其根本肯定是程序哪里有什么问题，但有时候换一套网格有可能（极小概率的哈）能避开这个问题。

@东岳 这是我44核心的log文件的最后几个时间步：之前的迭代次数也差不多： Time = 9.8 Courant Number mean: 0.733368 max: 1.27353 GAMG: Solving for p, Initial residual = 0.0998056, Final residual = 0.00082233, No Iterations 6 time step continuity errors : sum local = 4.80691e-07, global = 1.99762e-20, cumulative = 1.59798e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000123775 GAMG: Solving for p, Initial residual = 0.0809999, Final residual = 0.000584259, No Iterations 6 time step continuity errors : sum local = 3.45394e-07, global = 5.63939e-20, cumulative = 1.60362e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000120435 GAMG: Solving for p, Initial residual = 0.0189678, Final residual = 9.4483e-07, No Iterations 39 time step continuity errors : sum local = 5.51633e-10, global = 3.45318e-20, cumulative = 1.60707e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000119638 ExecutionTime = 170.51 s ClockTime = 171 s Time = 9.9 Courant Number mean: 0.73339 max: 1.27886 GAMG: Solving for p, Initial residual = 0.100098, Final residual = 0.000824461, No Iterations 6 time step continuity errors : sum local = 4.82558e-07, global = 8.66441e-20, cumulative = 1.61574e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000123823 GAMG: Solving for p, Initial residual = 0.0809731, Final residual = 0.000577791, No Iterations 6 time step continuity errors : sum local = 3.42103e-07, global = 4.07868e-20, cumulative = 1.61981e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000120471 GAMG: Solving for p, Initial residual = 0.0189619, Final residual = 9.31189e-07, No Iterations 39 time step continuity errors : sum local = 5.44568e-10, global = 7.2382e-20, cumulative = 1.62705e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000119674 ExecutionTime = 172.22 s ClockTime = 172 s Time = 10 Courant Number mean: 0.733413 max: 1.27933 GAMG: Solving for p, Initial residual = 0.0999699, Final residual = 0.000818608, No Iterations 6 time step continuity errors : sum local = 4.79761e-07, global = 8.13808e-20, cumulative = 1.63519e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.00012387 GAMG: Solving for p, Initial residual = 0.081007, Final residual = 0.000576888, No Iterations 6 time step continuity errors : sum local = 3.41793e-07, global = 5.45899e-20, cumulative = 1.64065e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000120508 GAMG: Solving for p, Initial residual = 0.0189675, Final residual = 9.02577e-07, No Iterations 39 time step continuity errors : sum local = 5.28191e-10, global = 1.77208e-20, cumulative = 1.64242e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.00011971 ExecutionTime = 173.99 s ClockTime = 174 s End 这是我32核心的： Time = 9.8 Courant Number mean: 0.733369 max: 1.27363 GAMG: Solving for p, Initial residual = 0.0998053, Final residual = 0.000728134, No Iterations 6 time step continuity errors : sum local = 4.26347e-07, global = 1.16708e-20, cumulative = 5.86486e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000123776 GAMG: Solving for p, Initial residual = 0.0810209, Final residual = 0.00054928, No Iterations 6 time step continuity errors : sum local = 3.24835e-07, global = 4.55333e-20, cumulative = 5.9104e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000120434 GAMG: Solving for p, Initial residual = 0.0189497, Final residual = 8.88034e-07, No Iterations 38 time step continuity errors : sum local = 5.18253e-10, global = 2.24608e-20, cumulative = 5.93286e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000119638 ExecutionTime = 167.35 s ClockTime = 167 s Time = 9.9 Courant Number mean: 0.73339 max: 1.27894 GAMG: Solving for p, Initial residual = 0.100093, Final residual = 0.00073128, No Iterations 6 time step continuity errors : sum local = 4.28756e-07, global = 1.40838e-20, cumulative = 5.94694e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000123823 GAMG: Solving for p, Initial residual = 0.0809929, Final residual = 0.000542798, No Iterations 6 time step continuity errors : sum local = 3.21522e-07, global = -2.66223e-20, cumulative = 5.92032e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000120471 GAMG: Solving for p, Initial residual = 0.0189464, Final residual = 8.9418e-07, No Iterations 38 time step continuity errors : sum local = 5.22767e-10, global = 2.26167e-20, cumulative = 5.94294e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000119674 ExecutionTime = 169.06 s ClockTime = 169 s Time = 10 Courant Number mean: 0.733413 max: 1.27934 GAMG: Solving for p, Initial residual = 0.0999826, Final residual = 0.000732193, No Iterations 6 time step continuity errors : sum local = 4.29802e-07, global = 1.57268e-20, cumulative = 5.95866e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.00012387 GAMG: Solving for p, Initial residual = 0.0810305, Final residual = 0.00054485, No Iterations 6 time step continuity errors : sum local = 3.22957e-07, global = 2.11478e-20, cumulative = 5.97981e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000120508 GAMG: Solving for p, Initial residual = 0.0189583, Final residual = 8.97181e-07, No Iterations 38 time step continuity errors : sum local = 5.24979e-10, global = 2.66036e-20, cumulative = 6.00642e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.00011971 ExecutionTime = 170.82 s ClockTime = 171 s End

@东岳 是想借助于充分发展湍流作为入口条件来做模拟，或者有什么别的途径可以实现？

请问一下，有哪位大佬用过datToFoam这个命令吗，*.dat文件需要写成什么格式呢？

https://packages.vmware.com/tools/esx/latest/windows/

上面那个公式引用是公式（8）

你这个直接在ICEM里面就应该能看出问题吧，在ICEM里面先看看网格质量之类的

当年，，人人网还叫校内网的时候。。。。我也特别喜欢这个图，我还用这个图做过封面。。做过头像。。

@东岳 我重新编译试过了，还是这个问题呢

就是从定义来求解，k还好说，epsilon和omega牵涉到一个特征长度的问题，这个取起来就不容易了，只能多调试。

运用pimpleFoam，求解一个平板的运动问题，想让整个网格一起运动，程序顺利运行了，但是结果明显是错误。猜测是边界条件设置的问题。想请教一下各位大佬，这种情况下，边界条件如何设置？或者是否有其他方面的问题。

奇怪了，虽然用起来没什么问题，但是打开的时候还是　会有warning Invalid $PV_PLUGIN_PATH and -plugin-path= not defined No supplementary ParaView/OpenFOAM reader modules Using builtin reader: paraFoam -vtk Created temporary 'platePimple.foam'

@xpqiu 还真是，我用paraview5.6.3可以完美显示出来了

链接文本，研究了一下用法，大家可以看看

@sibo 并行限制2核，网格限制是52万

李老师太客气了。免费看了这么久，作为新手其他的也没办法出力，能找几个小bug，也算是贡献微薄之力。 159页5.4.4中的“参见图5.12”应为”图5.11“。 我上贴最后一条写错了。不过新版UserGuide里5.4.5小节后增加了一小节网格细化功能的介绍，所以李老师更新时这些小bug应该会注意到。 161页5.4.6，“图5.14有一个例子”，应为“图5.13”。5.4.7，“参见图5.15”应为“参见图5.14”。 163页5.4.7，“边界层划分通过snappyHexMesh”应为“snappyHexMeshDict”。表中“minMedianAxisAngle 中轴点角度 130”应为“90”。 164页5.4.8， 表遗漏了minTetQuality关键字。（是因为V5版本后引入新功能取代该关键字的缘故吗？）

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

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

@东岳 是的李老师，p2是一个体标量场，0是整型，那这个地方应该怎么处理呢，求指点

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

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

在利用per phase求解湍流方程时，可以用以下语句获取各相指针。 Domain *dp = Get_Domain(2); /primary phase domain pointer/ Domain *ds = Get_Domain(3) ;/secondary phase domain pointer/ Thread *tp = Lookup_Thread(dp, 7); Thread *ts = Lookup_Thread(ds, 7); 2,3分别是各相domain_id，7为区域ID，可以从用户图形界面获得

@bestucan 感谢前辈的回复，目前问题已经解决 在编译OpenFOAM之前，需要先编译ThirdParty，但OpenFOAM WiKi中的Ubuntu18.04版安装教程里面好像并没有“编译ThirdParty这一步”，希望以后大家注意吧

我现在做出模拟结果但是后处理做不出湍动能，大家帮帮孩子吧，该怎么弄呀。

OpenFOAM中有pbm群体平衡模型的求解器吗？？？求指点

@纵浪大化 好的，常常交流联系

@东岳 谢李老师！

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

@东岳 不不不，您才是大佬，我是仰慕小菜鸡

@dzw05 谢谢

@东岳 谢谢李老师！

箭头是对指针的成员的引用，找到这个类的定义：turbulenceModel.H epsilon()函数声明在205行，是个虚函数，所以在 turbulenceModel.C 里没有定义，而是在使用其作为模板创建新类的时候再各自定义 在上面那个窗口里搜 epsilon() 函数： 看你的情况是用的哪种定义，其实里面 epsilon() 函数定义大差不差， 你用的 LES，咋找上层流的了，你找的那个确实是层流的对于 epsilon() 的一个实现，在我发的那个图片里也有 脉动速度都用上了，那应该不是层流了，速度分成脉动速度和平均速度这就是用了湍流模型了

@leilei 我用RSM做过旋流场气液流模拟，收敛性还可以，只是刚开始微调一点就行，但是计算速度比很多其他RANS慢很多。

我使用的buoyantBoussinesqSimpleFoam，边界条件设置如下 出现的问题如下

最近看了以下jasak的论文，找到了相关描述

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

文献有点大，上传不了，麻烦东岳老师下载吧（1）http://dx.doi.org/10.1016/j.ijsrc.2017.09.006（2） 10.1002/2013JF002911文献里都是用DEM方法模拟的，绝大多数文献里MPPIC都是用于气固或气液，想着试试用它做液固，跟别人的对比看看如何，增加些模型什么的，目前刚开始接触欧拉拉格朗日法，没什么思路

@雾雨魔理沙 不好意思哈，回复的比较晚...我的qq是891931269

@Calf-Z-DNS 你好有关于对流传热的问题想要咨询下你，就是对流出口边界条件的设置问题，是只对U使用的，还是全部都用呢？

@刘雄国 在 并行效率疑问 中说： 大牛说 下一代cfd 需要 大规模计算机网格 + Algorithm + Coding 整体考虑 我觉得很有道理 我说错了, 下一代cfd 需要 大规模计算机网格 + CFD Algorithm + CFD Coding + physical网格 整体考虑

李老师，您好，我是中国石油大学博1研究生，硕士阶段一直在用fluent做数值模拟。后来由于课题的需要，接触到了openfoam，才发现了openfoam的魅力。刚开始由于自己的计算流体力学底子薄弱，遇到了很多问题。我参考了李老师的很多资料，cfd中文网对我帮助很大，让我得到了提升，特别感谢。目前主要是以空泡动力的研究为主，希望自己能获得李老师的支持，谢谢李老师！ 李东岳：APPROVED，价格锁定45（京东），第二本京东当当均无货，到货后补

@东岳 那老师对这个问题有什么看法呢？standard k omega 在 Re=3000 时并不能表现出网格无关性，但相同的条件，在 Re=10000 下可以表现出网格无关性。

感谢分享

@东岳 哈哈，明白，谢谢李老师~

你需要系统看一下流体力学的相关知识，可能你没有这方面的背景所以不太理解。对物理空间中的一块区域而言，其所包含的物理量有几种变化方式，要么它自己随时间在变化（变化率，比如密度因为温度升高降低了），要么因为有流动带它进来或者出去（对流，比如能量），要么它自己在向外扩散或者外部在向这块区域扩散（扩散，比如这一块是高浓度区域而周围是低浓度），要么这块区域中这个物理量自己在生成或者消灭（源项，比如有化学反应）。所有的这些变化，都需要满足物理学的基础即守恒定律，其数学表达就是所谓的“广义”NS方程。这里“广义”是指这个方程描述了物理量的一般过程。 phi等于1，这个方程代表质量守恒；phi等于速度，这个方程代表动量守恒；phi等于内能，这个方程代表能量守恒；phi等于某个物质的浓度，那么这个方程就代表该物质的质量守恒，等等。

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

有幸参加过2020年1月的OKS课程，解开了我对求解器内部代码的困惑，虽然我还是菜鸡，但是我现在用商软，选择一个按钮，我会想到在OpenFoam里面都是对应哪个方程，方程的哪一项等等。早学早开窍啊 。

@西湖冷月 nice兄弟，d=====(￣▽￣*)b