DPMFoam求解中，重力居然会加速液体相？

upc_ngh

@zhangxc0223 在 DPMFoam求解中，重力居然会加速液体相？中说：

manualInjection
好的，明白了，多谢多谢

zhangxc0223

@东岳根据那个帖子的后续讨论，想提两个问题。

按照我的算例，重力是顺气流方向，如果我只想把重力作用在颗粒上，不作用在流体上，能不能直接删掉phiForces里面的 “+ rAUcf*(g & mesh.Sf())”项？直接删掉的话会有什么不合理的地方吗？源代码如下：

surfaceScalarField phicForces
(
   fvc::flux(rAUc*cloudVolSUSu/rhoc) + rAUcf*(g & mesh.Sf())
);

if (pimple.momentumPredictor())
{
    solve
    (
        UcEqn
     ==
        fvc::reconstruct
        (
            phicForces/rAUcf - fvc::snGrad(p)*mesh.magSf()
        )
    );

    fvOptions.correct(Uc);
}

我试着用这样修改后的求解器进行一些简单的算例（比如一个颗粒流体中，颗粒流体初速度均为0，让其在重力作用下加速）结果看起来没有问题。但是不知道会不会有其他的隐性问题，想请问一下各位。

还有想问一下，这个DPMFoam，考虑重力作用在流体上的初衷是什么？

李东岳

标准的拉格朗日离子方程需要通过连续相的压力梯度计算浮力。如果你假定连续相的不可压缩，这一项需要在拉格朗日粒子中同时重新植入一下。对于粒子项密度更大的时候，或许可以忽略

批注 2019-11-16 144145.png

你看下拉格朗日那面浮力怎么计算的。手头没有openfoam有点忘了

zhangxc0223

@东岳有个问题还是没明白，在竖直向下的颗粒-空气流中（重力方向向下，只考虑重力和拽力），速度关系应该是：固体颗粒的速度>颗粒流中气体的速度>纯空气流动中的气体速度。但是在DPMFoam求解出的结果中，无论总是得到：固体颗粒的速度<颗粒流中气体的速度。不知道这个是不是一个bug？还是说DPMFoam的求解初衷不适合这种流动？

zhangxc0223

@zhangxc0223 在 DPMFoam求解中，重力居然会加速液体相？中说：

@东岳有个问题还是没明白，在竖直向下的颗粒-空气流中（重力方向向下，只考虑重力和拽力），速度关系应该是：固体颗粒的速度>颗粒流中气体的速度>纯空气流动中的气体速度。但是在DPMFoam求解出的结果中，无论总是得到：固体颗粒的速度<颗粒流中气体的速度。不知道这个是不是一个bug？还是说DPMFoam的求解初衷不适合这种流动？

抱歉句子没说完整：
无论有没有这个phiforce的重力项，总是得到：固体颗粒的速度<颗粒流中气体的速度。不知道这个是不是一个bug？还是说DPMFoam的求解初衷不适合这种流动？

李东岳

@zhangxc0223 在 DPMFoam求解中，重力居然会加速液体相？中说：

固体颗粒的速度<颗粒流中气体的速度

为什么？DPMFoam在求气泡的时候气泡颗粒速度肯定大于液体速度呀

zhangxc0223

@东岳确实求出的结果就是这样，我也纳闷。然后我试了很多组模拟，发现在总体速度稍低（Ub≈3）的颗粒流中，颗粒速度总是会大于气体速度，然而在总体速度稍高（Ub≈12）的颗粒流中，颗粒速度只有在近壁区域会大于气体速度，其他区域都是小于气体速度。如下两图所示，其中纵坐标是颗粒或者气体的绝对速度，横坐标是距离中心线的距离，r/R=1表示管壁。我的模拟都是重颗粒（颗粒空气密度比为1000）

我一开始以为是网格的问题，然后根据速度的不同，按相同比例细化之后，发现仍然是这个趋势。不知道是什么原因导致的？会是DPMFoam中对颗粒速度处理的问题嘛？

我查找过相关文献，但是确实在稍高速度的条件下的研究很少，所以也无法确定是模拟的问题，还是颗粒流的真实情况就是这样？但是我的理解是颗粒收重力的作用，任何情况下都会比气体快的，这个理解没错吧？

李东岳

这是竖直管道颗粒流？边界条件都是什么样的？稀疏的？St数有多少。你这个截面位置在什么位置

zhangxc0223

@东岳对，竖直管道，重力和就是流向方向，出入口是周期边界，然后用添加了一个Ubar控制总体速度，初始速度场是从充分发展湍流的管流映射过来的。颗粒-气体体积分数是1o^-4量级，St在5左右。这个截面也是充分发展之后的结果。

李东岳

@zhangxc0223 我正在用矩方法算这个耶.. 针对你的问题，感觉跟这个很像，https://www.jianguoyun.com/p/DSS2YxIQ9s3ZBhjX_6oC 看图3，二者速度无差异

CFD中文网

DPMFoam求解中，重力居然会加速液体相？