compressilbeinterfoam质量守恒问题？

[litong189456]

请问compressibleinterfoam每一相的质量不守恒是什么原因？比如一个膨胀气泡，一个周期后气相（alpha2）的质量增加了，比如图b，一个气泡膨胀收缩后，质量增加了很多，这是求解器本身的问题吗，是哪里导致了质量的不守恒呢。

[赵一铭]

这些线都是什么？坐标坐标轴是什么啊？你用compressibleInterFOam模拟的单个气泡由于压力膨胀压缩？》

[litong189456]

@赵一铭
是的，抱歉没写清楚，横轴是时间，纵轴是当前时间气体的质量与初始气体质量的比值，发现气体的质量是逐渐增大的。红绿蓝分别表示网格从粗到密，发现网格越粗，质量的不守恒越严重。

[李东岳]

我在两年前也遇到过这个问题并反馈给基金会并没有解决，直观来看你的问题应该属于bug。或许你可以汇报给基金会。你的算例网格多少？2D么？

[random_ran]

@litong189456 在 compressilbeinterfoam质量守恒问题？ 中说：

格越粗，质量的不守恒越严重

如果是这样的结果，不正好说明了CFD计算结果很依赖网格的精度吗？随着你加密网格，质量不守恒的现象逐步消失，说明你的计算结果向准确值逼近。这正是网格敏感性检查要达到的目的？

[litong189456]

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

[litong189456]

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

[litong189456]

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

[李东岳]

类似的bug我在2014年遇到过，就是质量不守恒，并且我的验证表示和网格多少有关、是否并行有关，体现就是本来计算域里面有五杯水，算着算着变成四杯水了，后来反映给基金会也没处理掉，henry当时说是随机误差导致。再后来就没有结果了。后来在计算中我强制进行守恒了。类似的操作就是没个网格单元乘以一个乘数保证守恒。不过这个不守恒并不是很严重，我验证的不超过10%之内，一半就是2-3%左右。不知道你这个是否可以乎略。

不过你那个采用网格比较密的时候，可以认为近似守恒了。 @random_ran 说的也有道理，有可能网格比较糙，把质量数值扩散到周围去了 :cheeky:

[litong189456]

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