compressibleinterfoam的压力突变问题。

[litong189456]

请教一个问题，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

请问这是什么问题？

[李东岳]

compressibleinterfoam计算气泡膨胀

什么原因引起的膨胀？

[litong189456]

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

[mohui]

@litong189456 考虑过这个求解器本身存在一些问题吗？

[litong189456]

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