@李东岳 在 有关侵入边界法。 中说：

Hi，

OpenFOAM-extend里面有IBM方法，Extend和官方的OpenFOAM侧重不同，几个我了解到的不同主要是IBM，粘弹流体，动网格表面跟踪。

还有可压缩流动高分辨率格式。https://github.com/Unofficial-Extend-Project-Mirror/Unofficial-Extend-Project-Mirror.github.io/wiki

http://www.tfd.chalmers.se/~hani/kurser/OS_CFD/HrvojeJasak/DiscretisationBestPractice.pdf

看了这个，有点懂了，问题应该在于离散格式和梯度限制器啥的。
Finite Volume Discretisation in OpenFOAM Best Practice Guidelines

Hi,

并没有用过这个，不过cfd-online有关waves2Foam的讨论已经多达2000贴。从常识来看，附加湍流波高减小是可以理解的。

This questions is not clear. What is 液力变矩器?

@搬运工不好当 说：

求速度方程之后为什么没有求压力方程，而

Hi，

我觉得能量方程放在压力方程之后和之前结果是相同的。就像湍流求解，标量传输放在压力方程之前和压力方程之后也是相同的。

你可以自己测试一下并反馈？

@李东岳 谢谢东岳大神。祝论坛越来越好。

@李东岳 东岳兄，这个是一个基于OpenFOAM和LAMMPS的求解器，有一篇文章叫：Sun_2016_SediFoam A general-purpose, open-source CFD-DEM solver for particle-laden，可以一看。这个求解器其实功能是和CFDEM差不多，也就是求解密相流动的。其不属于CFDEM一部分。

Hi，

这和你的求解方程有关，例如：

fvm::ddt(rho) + fvc::div(U)
和
fvm::laplacian(nu, p) = fvc::div(phi)

是对称矩阵。更进一步的，第一个是对角阵。

fvm::ddt(U) + fvm::div(U)

则是非对称矩阵。

of用的是自己的容器

@xiaohsu2006 是根据温度和压力算出来的，　constant 参数后面的值是根据什么给的呢？

“\”好像是Linux的语法。 我的理解，是换行的意思，表示后面还有代码需要执行。你看最后一行就没有""了

musgs= rho*nusgs alphasgs = musgs/prt

对于mvCoeff 和　mvCoeff 两个变量，当他是可压缩的时候，成为密度的函数，该怎样修改呢？？

很感谢！！！

omegaWallFunction 这个边界条件继承自 fixedValueFvPatchField类，其构造函数如下
0_1460115657665_upload-47e7e05b-ec74-4419-8ff9-d55730af0dfa

这个构造函数里，需要注意两点，第一是传给基类fixedValueFvPatchField的参数，也就是 fixedValueFvPatchField<scalar>(p, iF, dict)，这一句里，基类会通过参数iF 去读取 类似 value uniform 2; 这样的语句，来讲当前边界的值初始化为 2。不过，看下面最后一句，this->operator==(patchInterField());，这里将当前边界的值重新赋值为临近边界的第一层网格中心的值。

所以，综合起来，value uniform 2; 这一句不起实际作用。

@搬运工不好当

这种情况很常见，如果使用第三方网格生成程序比如ICEM做的网格需要使用periodic严格定义才可以。如果你确定你的面试匹配的话，可以加大tolerance来忽略这个提醒，也不会有错误。这就是网格的原因。

这里的errors是误差的意思，而不是错误。可参考cfd-online上的这个帖子。

谢谢回答

fvm::xxx(u)，其中的u是需要求出来的，fvc::xxx(u)，其中的u为当前时间步的值，其返回一个场。所以，需要你要求某个场u，就用fvm。