MPPICInterFoam中静置情况下水流沸腾
-
如题,在图中底下的凹槽里填满里粒子,红色为水域,蓝色为空气。在静置的状态下,水流会像开锅一样沸腾。有遇到过类似问题的同学吗?希望能指出问题所在,谢谢。
-
@zhe 啥也不设置自己沸腾?压力?
-
@星星星星晴 对,什么设置都没有。观察了流速和压强状况,在粒子堆与纯水的交界处最强烈。那个地方的不同就是由于粒子的存在,alphac会出现梯度。但是在动量方程里又找不到这一项东西。有没有想法,有几个猜想也是好的。先感谢!
-
@zhe 是不是粒子到了packing limit了
-
颗粒和水的密度比大吗,是不是有可能是颗粒在浮力作用下运动引起的波动呢
-
@星星星星晴 @李东岳 @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是因为我的粒子浓度太高,不能忽略压力项。李东岳老师在其他帖子里也讨论过这个问题。根据自己的粒子粒径或占网格比来定。
-