@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也可以。

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

例如下面的两相界面处的jump condition，在表面张力系数σ不为常值函数时，沿界面方向的表面张力作用不能被忽略
[μ∇u · n] = ∇s σ
其中∇s =(I - n⊗n)∇，表示表面梯度算子
在interFoam中表面张力系数为常值函数，用CSF模型计算法向jump condition，想知道在多数CFD计算中，这种jump condition都会怎么处理?

好的 谢谢 我去看一下