compressilbeinterfoam质量守恒问题?



  • 请问compressibleinterfoam每一相的质量不守恒是什么原因?比如一个膨胀气泡,一个周期后气相(alpha2)的质量增加了,比如图b,一个气泡膨胀收缩后,质量增加了很多,这是求解器本身的问题吗,是哪里导致了质量的不守恒呢。
    0_1495070776057_TIM图片20170518092538.png


  • 管理员

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



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



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


  • 版主

    @litong189456compressilbeinterfoam质量守恒问题? 中说:

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

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



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



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



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



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

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



  • @李东岳
    好的,谢谢东岳老师。


登录后回复
 

与 CFD中文网 的连接断开,我们正在尝试重连,请耐心等待