东岳大神说了开发的问题。我来说说我的看法？清楚物理模型的物理意义貌似是最基本的也是最容易出错的部分。openfoam方面可以参考The Finite Volume Method in Computational Fluid Dynamics。还有就是最笨的办法：读代码

@asininno 我没有详细对比过具体有多大影响。关于非正交修正，建议你仔细看一下公式，《The Finite Volume Method in Computational Fluid Dynamics An Advanced Introduction with OpenFOAM and Matlab》这本书有很详细的说明。说白了就是采用和正交网格上一样的方法去做离散，对于非正交网格会有误差，那么再加上一些措施进行修正。
但是还是建议做好网格，控制网格的正交性和弯曲度；毕竟修正只是一种补救措施。

@dzw05
哦了，有空整理整理。

好吧，你说的很对，看来就是：
\begin{equation}
\frac{\nabla p}{\rho}=\frac{\mathrm{D} \mathbf{U}}{\mathrm{D} t}
\end{equation}

of中多相流这块不可压好像都没有能量方程，在模仿着可压的能量方程添加时发现，剪切力做功是用湍流粘性为形参的库函数表示的，现在我不调用库函数，怎么表示这一项呢？谢谢

@wwzhao
谢谢赵老师的回答，我尝试下写入不同的值看是否能正常写入文档

有关Eo数的定义：http://47.93.231.95/doku.php?id=cfd百科#eoetvoes_number

谁那见过三维Eo数的计算公式和文献出处？OpenFOAM植入的这个公式不知道从哪来的。

@李东岳 额，从controlDict里的时间来看，多相计算的物理时间要远大于单相，但还是没能得到合理结果。我再算下去看看。

@dzw05 好的，真的非常感谢:happy: :happy:

@李东岳 嗯，调分辨率就可以了。东岳老师我最近在用二相欧拉模型求解一个水沙运动问题，但是算到一定时间就会蹦掉，您有空帮我看下我的边界条件设置是不是有问题吗？

https://github.com/OpenFOAM/OpenFOAM-2.1.x/blob/master/src/finiteVolume/finiteVolume/fvc/fvcMeshPhi.C

看这个代码

注入质量

如果你指的是往计算域里面添加液体的话，设置进口边界条件不就好了？

这图应该指的是在网格分辨率足够的情况吧，如果分辨率不够，在0cell和1cell之间不会只夹着一个0，1之间的cell

计算完毕后，规定任意一条直线及点数，如何插值求出直线上任意一点的值？

特别是计算区域分为几个小的区域，每个区域内有不同的物性参数的情况如何计算？

@李东岳
的确奇怪，我现在都不知道为什么纯粹进行值之间的交换会出现问题，又没有显示命名的重复，
不过我把发生错误的变量名字改了以后，错误就不在出现，虽然解决了，但还是对于错误的原因存疑。

你好，请问你这个问题解决了吗？

你在MRFProperties里面或者fvoptions里面需要定义旋转的区域，paraview是显示的那个效果那个变不了。你可以把结果当做稳态时候的结果。或者搅拌器叶片转到某一个时刻，流场就是这个样子的。

@wwzhao 谢谢BZ！我试试~

好像以前见过一个类似的帖子，下面那个value 是没用的，主要用于paraview读取，你随便写100也可以。

这个问题通常是由于内存不够引起的崩溃，你机器的内存多大？算例的网格多大？