各位前辈好，关于外流场计算中的动压、总压等有些疑问，想和大家请教一下：
对于不可压缩的流体，总压=静压+动压。OF的多相流计算中，初始条件中需要给一个参考压力p_rgh，然后经计算得到总压p。
那么，以一个初始位置在水位以上的结构物（即没有静水压力）为例，当结构物受到波浪冲击时，结构物所受到的流体压力可以认为全部是动压吗？或者说，OF计算出的在自由液面处的压力都是动压？

在计算多相流动网格案例时，计算一段时间后库朗数突然增大并导致计算发散。起初以为是网格变形过大导致的发散，但查看日志文件和结果发现，计算并非是结构物受到波浪冲击后达到最大变形时终止的，而是在逐渐回到原位置时终止的，此时相分数alpha.water超过了1并在急剧增大。这时结构物附近有一定的越流（如下图所示，截图截得不太整齐请大家见谅），且计算停止时（44.2s）的气液交界面显然是有很大问题的。因此我认为应该是相方程计算不收敛进而导致的整个计算终止。

日志文件中相分数在44.076s开始大于1：
最后计算崩溃：

因此，想请问大家是否在多相流、动网格的计算中是否遇到过类似的问题，有什么办法可以解决这种相分数越界的问题吗？

想了一下，虽然网格有变形但整个过程中网格拓扑不变，是否有办法根据该单元的拓扑来提取这个单元的压力时程呢？

以结构物受波浪冲击导致结构物产生变形为例，想要得到距离结构物最近的一层流体网格上某一单元的压力变化情况。在以往的非动网格案例中，我是在system目录下设置一个sampleDict文件并在其中指定提取方式、提取点坐标，但是动网格案例中位置会变，这样的做法显然不行。请问各位是否有比较好用的提取方法？

@guohuiqun
目前，我们是通过渗透试验测出孔隙介质的压差、流速，拟合出I=AU+BU^2中的A和B，然后反算出来a和b的。至于c，OlaFlow作者在其论文中提到该参数的影响很小，故一般是取经验值0.34。

@bestucan

如果是计算时要实时计算边界条件，这个可能符合你的要求：
http://openfoamwiki.net/index.php/Contrib/groovyBC#Usage_of_the_boundary_condition

感谢回复，需要的确实是在计算中实时计算的边界条件，我先了解一下

目前在尝试用质量流率\体积流率边界flowRateInletVelocity来实现，但是发现官方给出的说明中好像只能施加一个常数的质量\体积流率边界，然后根据此流率边界得出此边界上的速度，无法根据边界上的流速得到体积流率（和我要实现的恰巧相反。。）
请问OF中有无根据边界上的流速来增加流体体积的方法呢？

各位好，想请教一个关于多相流的流固耦合问题。具体要模拟的情况就是一个孔隙介质挡板，左侧有水右侧无水，左侧的水进入挡板并从右侧渗出。

目前已在固体求解器中算出了固体边界上的变形和渗流速度并将二者加到了OpenFOAM流体域的相应边界上，请问怎样才能通过该渗流速度在边界上相应的位置补充上因渗流而增加的水呢？

再和各位前辈请教一个相关的问题，在用snappyHexMesh处理网格量比较大的案例时（网格数量200w左右，模型和上图类似，有一些拐角）：

1. Explicit feature edge refinement设置成level 1就会出现FOAM FATAL ERROR: attempt to read beyond EOF；
2. refinementSurfaces设置到level ( 1 1 )还能正常生成网格，提高到( 2 2 )就报错face 17 area does not match neighbour by 0.0484736%。
   请问这两种情况是否就是网格太多导致的？有没有什么解决办法呢？

@李东岳
东岳前辈好，我也没找到那个帖子中的Use VTK Polyhedron(通过paraview和paraFoam打开都没有)，但是有这个Decompose polyhedra，默认是勾选上的。如果将它取消勾选，显示的网格会有一些不同，但是扭曲的部分还是存在，如下面的对比。


另外，即使用了那个帖子中的第2、3点，显示的结果也是没什么区别。