MPPICInterFoam中静置情况下水流沸腾

zhe

如题，在图中底下的凹槽里填满里粒子，红色为水域，蓝色为空气。在静置的状态下，水流会像开锅一样沸腾。有遇到过类似问题的同学吗？希望能指出问题所在，谢谢。

星星星星晴

@zhe 啥也不设置自己沸腾？压力？

zhe

@星星星星晴对，什么设置都没有。观察了流速和压强状况，在粒子堆与纯水的交界处最强烈。那个地方的不同就是由于粒子的存在，alphac会出现梯度。但是在动量方程里又找不到这一项东西。有没有想法，有几个猜想也是好的。先感谢！

李东岳

@zhe 是不是粒子到了packing limit了

Tens

颗粒和水的密度比大吗，是不是有可能是颗粒在浮力作用下运动引起的波动呢

zhe

@星星星星晴 @李东岳 @Tens 找到原因了。就像之前观测发现的现象，在粒子群和纯水的界限处流速和压强极强。问题是出在压力泊松方程上。MPPICInterFoam应该按着interFoam的连续方程，考虑alphac的影响。但目前MPPICInterFoam是按照单项流的DPMFoam的方程。ddt(alphac)这一项是造成水沸腾的原因。把现有的

    surfaceScalarField Dp("Dp", alphacf*rAUcf);

    fvScalarMatrix p_rghEqn
    (
        fvm::laplacian(Dp, p_rgh)
     ==
        fvc::ddt(alphac) + fvc::div(alphacf*phiHbyA)
    );

改成：

    surfaceScalarField Dp("Dp", alphacf*rAUcf);

    fvScalarMatrix p_rghEqn
    (
        fvm::laplacian(alphacf*Dp, p_rgh) == fvc::div(alphacf*phiHbyA)
    );

这样就不会出现“无风自动”的现象了。是不是对别的地方会造成影响，还需要再研究。（其中出现的alphac*alphac是因为我的粒子浓度太高，不能忽略压力项。李东岳老师在其他帖子里也讨论过这个问题。根据自己的粒子粒径或占网格比来定。

李东岳

@zhe