CFD中文网

@队长别开枪 您好，最近在用您的PLIC工具，发现了一个问题，目前我还没有找到原因，不知您是否遇见过。
当我计算的时候，发现在纯水或者纯气体区域会出现另一相，二维和轴对称都出现了这个问题，用您的输出界面的writeplicfaces功能也发现，在本应该是纯水的区域产生了气体。黑色的是writeplicfaces输出的界面，粉色的是体积分数1e-8到9.9999999的单元，本来应该只在界面上。

Cellcells()不能并行，东岳老师是怎么做的？@东岳

@东岳
东岳前辈，我看您的一篇文章《Simulation of bubbly flows with special numerical treatments of the semi-conservative and fully conservative two-fluid model》，其中关于可压缩相方程boundness这里有一点没看懂，可以请您详细解释一下吗？

@队长别开枪 谢谢

@队长别开枪 您是做带有几何重构的VOF算法的，我想请问一下，现有的算法都是在不可压缩流体中使用，对于可压缩流体，这些几何重构算法（比如PLIC）还能适用吗，为什么现在没有相应的文章呢，是没有人做，还是说在算法上有难以解决的问题呢，这个问题困扰我好久了。

大家好，关于文章《Simulation of bubbly flows with special numerical treatments of the semi-conservative and fully conservative two-fluid model》中关于相方程有界性有一些疑问想要请教。

，

请问为什么速度散度是0，相分数就是0到1之间呢?

公式30的等号右侧在每一个相分数是0的时候，并不等于0啊。

请问有详细的讲这方面的书或者文献吗？

@yhdthu 看了一下他的文章，发现就是基于不可压缩流体来推导的方程，您觉得能够推广到可压缩流体吗？

我还得更深入的阅读代码了，太感谢了。

@李东岳
好的，谢谢东岳老师。

@程迪
您写的博客实在太详细了，关于openfoam的piso算法，

1.这个是由piso.correct()控制的，而piso.correct()是由变量label corrPISO_控制的。
2.重复的时候，UEqn所引用的U是有更新的。所以用UEqn构造的UEqn.A(), UEqn.H()也可能会有所更新。

之前有帖子曾经讨论过这个问题，http://www.cfd-china.com/topic/589/关于piso循环的二次修正的一个疑问/15 。

@李东岳 在 关于piso循环的二次修正的一个疑问 中说：

piso循环的二次修正是建立在第一次修正的基础上的。第二次修正的连续性方程是涉及到u**。上面的代码里，UEqn并没有做更新。如果忽略函数fvc::ddtCorr(U, phi)，pEqn也是没有变动的。
问题1： 这是不是就意味着对同一个pEqn重复一次迭代运算？

虽然构建的速度矩阵Ueqn没有更新，但是HbyA有变化，因此求解的pEqn不是同一个pEqn

问题2：我对问题1的理解是，在icoFoam的算法里，速度u通过连续性方程被隐性的耦合到了pEqn里面。第二次的修正是没有必要的。这与原始的PISO算法有区别。请问这样的理解对吗？

不是很清楚你说的哪个二次修正，如果你说的是楼上指的非正交修正，一般情况下不需要进行，除非在那个potentialFoam中。因为这个求解器没有进行速度-压力耦合迭代。

为什么Ueqn没有更新，但是HbyA有变化，HbyA不是等于UEqn.H()/UEqn.A()吗？ Ueqn没有更新，那么UEqn.H()和UEqn.A()应该没变啊，为什么HbyA反而有变化。

@李东岳
感谢东岳老师。

• 主要原因是教材都非常老了，经典教材第一版写的时候都是90年代。同理，教材都是采用的交错网格。虽然同位网格有介绍。如果2018年有人写一本新的CFD巨著，很有可能用的就是压力方程而不是压力修正方程。

您的意思是压力方程和压力修正方程都一样，还是说压力方程相对于压力修正方程更有优势？

• 但为了遵守SIMPLE算法的基本思想，压力修正方程的介绍还是有必要。但是你看PISO直接用的就是压力泊松。但相对于PISO，SIMPLE更容易理解因此教材中主要介绍原始的SIMPLE压力修正。

我看computational methods for fluid dynamics等教材, PISO也是用压力修正方程推导的。

是不是对于可压缩的压力基计算，压力修正方程更适合一些呢？
对于压力方程：

只有是rho=a*p这种形式的状态方程才可以直接把对流项里面的rho换成p。
但是对于压力修正方程：
对于任何状态方程都是可以用的。
谢谢您！

请问一个压力基的的问题，压力方程和压力修正方程的区别到底是什么，为什么几乎所有的教材在讲述simple和Piso算法的时候都是用的压力修正方程，但是openfoam里面却全是压力方程？

@mohui
考虑过，但是刚开始压力发生这种令人不解的变化，之后压力又能够恢复正常，气泡脉动在流场产生的压力又在合理的范围内。但是我对程序进行定点调试，就发现经过我图3的代码之后，压力就突然变了。

@李东岳
气泡内部初始压力是高于外界的，所以气泡会膨胀。但是发现布置的probe测点压力不太对的时候，就让时间步特别小1e-9，结果两个时间步之后，压力就不对了。

请教一个问题，compressibleinterfoam计算气泡膨胀的时候压力场问题。我的时间步选择是1e-9，按理说经过有限几个（4-5个）时间步，由于气泡的膨胀产生的压力扰动还根本没有传播到远处，一定范围外的压力应该是初始时候的压力，即压力随着水深增加。 1e-9时刻，压力分布还是合理的，如图1。但是2e-9的时候，流场的压力产生了一个很大的变化，如图2。
图1

图2

经过我对程序的调试，发现在2e-9解压力方程之后，压力发生的突变，解压力方程代码如图3。没经过这段代码之前，p_rgh合理分布（图4），经过图3的代码之后，p_rgh发生了很大的变化（图5），
图3

图4

图5

请问这是什么问题？

@lxjcfd
请问您那里有第三篇文章的原文吗，在网上怎么找也找不到。[3] Weller, H.G., 2002. Derivation, modelling and solution of the conditionally averaged two-phase flow equations. Technical Report TR/HGW/02. Nabla Ltd

@李东岳
好的，谢谢东岳老师。

@李东岳
请问您说的Bug的意思是说这是算法本身的问题吗？

@random_ran
如果仅仅想要加密网格的话，计算量太大。我想请教这个质量不守恒是什么原因导致的，是算法本身的问题吗？ 有人说是使用了状态方程（而不是解能量方程），也有人说是数值扩散(Numerical Diffusion)。因为我也刚刚接触发现这个问题，还不是很懂。

@李东岳
谢谢东岳老师，我算的是轴对称的例子。