@wwzhao 在OpenFOAM-2.1的时候，程序里写的就是“CrankNicholson” ，所以controlDict里也应该写“CrankNicholson”，虽然这个拼写确实是不正确的。但是从OpenFOAM-2.3开始，这个名字改成了正确的拼写：“CrankNicolson”，相差一个字母。所以这就导致了2.1和2.3的算例会在这里不兼容。

@搬运工不好当

在对拉普拉斯项离散的时候，如果使用非正交网格会引入误差，因此需要非正交修正。所以nNonOrthCorr的作用就是消除非正交网格引入的误差。

@搬运工不好当

压力泊松方程方程本身左边只有一项，这种结构一般压力给定场松弛因子，OpenOFAM里面速度预测方程速度给定方程松弛因子。系数的给定请参考楼上。

在压力修正方程中有关场松弛请参考：链接文本

PISO不需要给定松弛

你的arc定义错了，你觉得一个点能定义一条arc吗？

你得把错误内容贴出来，不然一句话看不出来。出错的原因太多了，而且不知道你修改了求解器，甚至case的哪部分...............

@搬运工不好当

多相求解器为了数值稳定性可能pimple更稳定。如果细看的话，piso在openfoam里面并不多。不是simple就是pimple。毕竟pimple可以实现piso。

下面是一个邮件里面的内容：

One comment I have is that the new solvers definitely requires PIMPLE to be used. Meaning, the number of outer correctors can hardly be 1 if you want a converged solution. In order to avoid iterating for too long, I usually fix required residuals. This will increase the computational time, since it takes between 5 and 20 iterations to converge.

虽然他建议使用pimple，不过在我的多相求解算例我用的piso。并且结果没有差别。
openfoam自带的算例好多都用的pimple。

使用pimple和piso我个人更感觉是一个数值稳定性的问题而不是迭代求解算法的问题。

transportpProperties没有声明。

这几天放假没在学校，才看到回复。谢谢分享 。我研究一下。万分感谢

是的，可以忽略。

并没有用过pyfoam

@xpqiu 谢谢！

在liggghts/src下面找到这个fix_heat_gran_conduction.cpp文件，里面讲了颗粒和壁面之间的静态导热的计算，看看代码应该可以知道的。

用surfaceCheck检查STL文件的结果呢？

diagonal是对角矩阵求解器。

fvScalarMatrix rhoEqn ( fvm::ddt(rho) + fvc::div(phi) == fvOptions(rho) );

从这个方程里面可以看出，只有时间项是隐式离散，fvc和fvOptions都进入了矩阵右边。因此离散之后是一个对角阵。直接除过去就可以。并不需要迭代。因此迭代数为0.

这个原项是怎么弄得？fvOptions，之前我记得都是MRFSimpleFoam，SRFSimpleFoam，现在都整合到fvOptions里面了。原理是什么啊？

@mark 这种“提示文件找不到”的都是路径搞错了，我看你已经编译出来了，仔细查看查看路径。

把fvsolutuon里的tolerance由e-6改为e-8对结果影响大不。

像大家请教一个问题，算低速流体（10m/s以下）且求解域很小（1cm1cm1cm量级）当用可压求解器求解好像CFL不管设多小，算足够长时间都会发散，除开用不可压求解器以外还有什么其他解决方式么？

@李东岳 想明白了，多谢东岳大神！

是中国人，还没改国籍。

互联网的某个角落 or 提问。

@搬运工不好当 每个边界都要赋值，格式大概是
nonuniform List<scalar>
xxx
(
aaa
bbb
ccc
.....
);
xxx 要等于边界面元的数目，括号里的数字的个数要等于xxx。

另一个用法是vtk可以存储部分数据，而不是整个的场。大型算丽非常有效。

另外有些paravirw的后处理需要vtk做源source.

@dyj19901127 hi.

我觉得一个好的习惯就是禁止复制构造函数。除非你每个类都写，但一般也用不到。henry应该是采取不写的做法吧。

这个在c++ primer plus 和effictive c加加里面提到过。书并不在手边你可以看看那个类里面存在指针的时候，复制构造函数怎么泄露的。我更觉得这是个和编程有关的问题和cfd关系不大？

@搬运工不好当 你先wclean一下试试:）

@搬运工不好当 shm生成的网格有的时候非常的tricky，不错如果你采用shm成功的生成了mesh okay的网格，openfoam大部分求解器都能算收敛。of里面单相求解器已经非常稳健了。

如果对于现有几何，一定要生成一个非常好的网格，我建议你

简化几何 换专业的网格生成程序 慢慢的，充分的调节shm里面的网格控制参数，主要是网格密度，迭代次数，甚至删除边界层，这个步骤需要很长时间。

话说回来，shm对于复杂几何生成的是混合网格。虽然单相求解器对混合网格非常稳健，但某些多相求解器还算不了。

@hongfu2233 不，湍流粘度并不是一个直接求解的量。并没有湍流粘度的方程，因此不需要初始值。

@liliang2009223 我都那个丸不是很了解。这俩个互溶么？液体

@mark dgdt分析已经上传

对于第一类边界条件，OF里面壁面法向导数只有一介精度吧？

Foam::tmp<Foam::Field<Type> > Foam::fvPatchField<Type>::snGrad() const { return (*this - patchInternalField())*patch_.deltaCoeffs(); }

请问有没有二阶精度的sngrad?

请教一下，hooke/stiffness model怎么使用，有成功的吗？能否分享下。

@搬运工不好当 说：

在某文献里看到用OpenFoam划分高速列车流体网格，六面体

目前OpenFOAM处理这种非常简单，snappyHexMesh是一种自动网格划分的生成器。特别擅长画复杂几何。不过从个人经验来看：

你不需要过分追求细节； OpenFOAM不是纯六面体网格；

因此要是发文章我更倾向与使用ICEM

能不能在水槽内部形成一个固体边界？不紧挨着边界处

If you dig deeper, 你会发现ThirdParty里面不是OpenFOAM基金会开发的。例如：paraFoam是paraview的插件，paraview属于Kitware，scotch属于法国的一个...看这里：链接文本。所以称作第三方包。

你是指的取断面为边界入口吗？这样的话可以采用timeVaryingMappedFixedValue；

@搬运工不好当 说：

得到一个特定时间步结果

不得不说如果调整了时间格式为何还要期望获得一样的数据结果？不能期望每一个格式获得同样的数据。

Hi,

这里我并没有发现错误？

Hi,

你的twoPhaseProperties没有定义，如果你是在interFoam中调用mu，在interFoam的createFields.H中已经定义了：

Info<< "Reading transportProperties\n" << endl; incompressibleTwoPhaseMixture twoPhaseProperties(U, phi);

所以可以这样写：

volScalarField& muc(twoPhaseProperties.mu());

试一下。我没有测试，有问题再联系。

@mark Hi mark,

错误提示，缺少twoPhaseMixtureThermo.H文件。可以看看这个文件有没有在你的Make/options中包含

Hi,

我曾经调用过Fortran库来计算矩阵特征值，不过后来我采用Alberto的库而不用Fortran了。在OpenFOAM里面调用Fortran混编需要安装gfortran，可以在终端下通过：

sudo apt-get install gfortran

来安装。安装后下载你所需要的Fortran源代码（目前我还不会写fortran，我用的是现成的包），里面通过makeFile来制定规则，在终端通过make编译就可以。编译出来的是.so文件。OpenFOAM可以直接调用。

也可以在./Allwmake里面把make写为脚本。

使用入口边界条件，turbulentInlet。
或者延长入口使用，maped，类似于前面加一个自循环的pipe效果。
或者使用，timeVaryingMapped，提供入口出的数据。

很早之前用过ssh client之类的东西，可以利用图形界面进行登录，文件拷贝之类的东西，不过运行程序还是需要命令行。

不过还是使用ssh命令行靠谱啊，时间长了效率会比图形界面高的。

这个倒是一个问题，因为最近做到了流化床的燃烧，那么里面肯定是可压缩的，因为燃烧放热，密度变化等，气体为可压缩的。CFDEM在2.9.0好像是可以支持压缩的，但是感觉：流速比较低，处于低马赫数，气体压缩效应并不明显，所以我看绝大多数的冷态CFD-DEM计算都用不可压缩的。我们这里有做高马赫数（激波）等，这时候应该要考虑可压缩性了。

Hi，

你可以试试scotch方法，从他的介绍来说是自动的处理网格疏密问题。可以这样测试一下：分解之后用paraview看一下各个块怎么样。我一直用scotch分解，不过还没测试过。这几天有空我测试下。

@mark Hi,

等我有空把文档写完，那个暂时一句俩句说不清楚。

因为流体力学中常用的是前五个量纲，因此可以默认忽略最后两个量纲，详见openFoamUserManual_PFM第19页。

Hi,

目前我知道的有在windows下面有使用GPU加速计算的，例如Symscape，还有这个：RapidCFD。

前几天看到在CUDA上求解矩阵的算法，感觉很有意思：CUDA。我的理解GPU计算就是在CUDA上用CUDA支持的语言编程来求解矩阵？

暂时知道这么多

Hi，

采用MULES方法求解的相场不应该存在越界问题，不过守恒倒是有一些问题。能否详述一下？

OpenFOAM里面相场越界这个问题已经是一个很重要的问题了。从OpenFOAM-1.4开始OpenCFD公司（目前已注销）就着手采用MULES来处理标量场的守恒问题，并逐步完善比如添加LTS功能等。

@搬运工不好当

Hi,

PMV目前在海尔科学研究院有一套标准，国际上也有一套ISO9000x的标准。在这个标准中，除了公式的值都需要自己设定。然后在PMV框架下迭代求解。可以google一下PMV ISO标准。

Hi,

请注意epsilon的不能赋值为0，因为在kEpsilon方程中的湍流粘度定义为：

\begin{equation}
\nu_t=C \frac{k^2}{\varepsilon}
\end{equation}

明显的如果epsilon为0，湍流粘度不可以计算。