MPPICFoam/denseParticleFoam 对于重力的考虑

[sxz0823]

（1）如果在求解时加入g（0 0 -9.81）则流场outlet处会出现不合适的速度同时有压差出现：


（2）改变g（0 0 0），流场速度正常，没有出现其他变化，
以下附上p，U.air文件
p

inlet
{
    type            zeroGradient;
}

outlet
{
    type            fixedValue;
    value           uniform 0;
}

walls
{
    type            symmetry;
    value           uniform 0;
}

top
{
    type            symmetry;
    value           uniform 0;
}

bottom
{
    type            zeroGradient;
}

U.air文件

inlet
{
	type  		fixedValue;
	value 		uniform (11 0 0);
}

outlet
{
    type            zeroGradient;
     
}

walls
{
    type            symmetry;
}

top
{
    type            symmetry;
}

bottom
{
    type        fixedValue;
    value       uniform (0 0 0);
}

但是在我没有添加kinematicCloudProperties以及相关颗粒离散格式、求解方式设置时这个边界条件是正常的，如果我在 kinematicCloudProperties 文件中打开了 gravity，在constant文件中g为（0 0 0），那么例子所受重力还存在吗（这个我还没有测试）。同时我在查阅以往的记录里面发现了相似的问题（链接如下，不过时间有点遥远），有没有相关大佬遇到过这方面问题，可以指点一下吗，谢谢？
https://www.cfd-china.com/topic/394/如果用dpmfoam求解稀相流会怎么样-误差大么

[李东岳]

你把网格处理少一点，然后再openfoam10里面复现一下错误，有错误你发给我（算例清理干净）我算一下

[sxz0823]

@李东岳 2bdedb7f-a996-40c4-84db-584bbf26f12f-TP.zip

[sxz0823]

@李东岳 感谢李老师

[sxz0823]

@李东岳 注射时间应该是一直注射的，不是6s

[李东岳]

那你重新发一个

[sxz0823]

@李东岳 TP.zip

[sxz0823]

@李东岳 麻烦李老师了

[李东岳]

openfoam10的denseParticleFoam没跑起来

[sxz0823]

@李东岳 李老师你好，of10用的比较少，但是我切换了一下of10，使用的MPPICFoam跑动了，还是和上面一样。denseParticleFoam没有使用过

[sxz0823]

@李东岳 TP_of10.zip 东岳老师，这是of10能运行的文件
https://www.cfd-china.com/topic/394/如果用dpmfoam求解稀相流会怎么样-误差大么/47
和这个一样的情况，在不改变g情况下应该怎么修改呢，因为我的粒子需要考虑重力

[李东岳]

工作日给你看一下 :-)

[sxz0823]

@李东岳 非常感谢东岳老师

[sxz0823]

@李东岳 因为像是我的场属于大气流动，有明显的压力梯度，应该把p换成p_rgh，应该是把DPMFoam求解器里面p改为p_rgh，我是这样考虑的，不过还要您修改一下看看，我这个想法的实现过程还不太明白

[sxz0823]

@李东岳 东岳老师，您好，我还有个考虑，我的是气固流动颗粒在70um，可不可以将连续项里面的重力加速度项去掉 phicForces里面的 rAUcf*(g & mesh.Sf())

这样会影响到我粒子的重力计算吗，此时我的g仍然给（0 0 -9.81）

[李东岳]

不能去掉。去掉就不对了。

[sxz0823]

@李东岳 那东岳老师 ，您在这篇帖子中提到的办法是指什么呢，1方法换成p_rgh，“直接将将pEqn.H改为p_rghEq.H，是无密度速度方程+prgh方程，虽然流场预测正确，但这种植入在密度不为1的时候耦合项不正确导致速度过低。（引用了您在这篇帖子的结论）”；所以我想着有没有办法去掉动量方程中的重力影响，但是保留颗粒相的重力（在不可压缩的前提下）
https://www.cfd-china.com/topic/394/如果用dpmfoam求解稀相流会怎么样-误差大么/59 ，

[李东岳]

这个已经过去太久了，我已经忘了怎么回事了.. 等我处理个bug之后再给你看看吧

[sxz0823]

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

[sxz0823]

@李东岳 东岳老师，我把rAUcf*(g & mesh.Sf())注释掉以后，除了压力场为0，似乎没有造成其它影响，但是肯定不符合正常要求了，那还是要在MPPIC里面添加p_rghEqn.H吗？东岳老师可以简单指导一下吗