@李东岳 这一项和interpolationSchemes改为linear就可以运行，用MUSCL和venLeer就不可以

@聆听whatever

Hi,

可以看看这个：链接文本。欢迎补充交流。:baby:

我的并行出问题了，不并行倒是可以，我想具体学一学怎么用Debug模式来调试，能否举个例子教一下我呢？谢谢各位了~

@lhzhu 在 如何拷贝/创建 volVectorField 在每一个时间步文件夹中 中说：

只要0目录里面有一个场文件，decomposePar 都会在processorX里面生成子区域的场文件。是不是你写的0/HAG_BC 的file header 里面的object 不是 HAG_BC ?

decomposePar 后 在“processorX” 中的 "0" 目录里有这个HAG_BC, 问题是下一时间步的时候没法创建这个HAG_BC啊，
所以会有那个error.
不知道怎么解决

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

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; 这一句不起实际作用。