如果用DPMFoam求解稀相流会怎么样?误差大么


  • 管理员

    将noSlip改为

    type  fixedValue;
    value uniform (0 0 0);
    


  • @东岳 谢谢老师,流场确实收敛了。我注意到您做的更改中关掉了重力g(设为(0 0 0)),出于什么考虑呢?测试的考虑重力后就不合理了.


  • 管理员

    @alvin如果用DPMFoam求解稀相流会怎么样?误差大么 中说:

    出于什么考虑呢?测试的考虑重力后就不合理了.

    我需要深入研究一下。主要起源于体积力重力的处理,从期刊看到的公式来看,重力项的处理不太一样,比如下面这俩个:

    0_1523417739053_捕获.PNG

    0_1523417745705_捕获2.PNG

    我怀疑还是粒子压力梯度的处理问题,http://www.cfd-china.com/topic/1488 在弄清楚之后,可以从代码上进而在结果上解释原因。

    目前的猜测,DPMFoam连续性方程里面的phiForces引起的重力驱动流动。后续我会更新。



  • 谢谢您耐心的讲解,受教了,可是如果关掉重力,DPmFoam求解器是不是就不能考虑颗粒受重力的情况了,感觉这个求解器的处理不是很理解


  • 管理员

    问题是phiForces导致,比如重力向下,那么phiForces将为负值,其产生一个向下的速度。目前我自己的求解器也遇到了这个问题。还需要进一步研究。还不能说是DPMFoam的bug,因为我用这个模拟过气泡流,结果不错。

    我要看一下密度为常数的情况下,速度和压力的平衡处理。


  • 管理员

    一个解决办法是将浮力项和重力项进行下面的转换(参考其他求解器):$\nabla p - \rho \mathbf{g}=\nabla p_{\mathrm{rgh}}+\mathbf{g} \mathbf{h} \nabla \rho$,我植入看看。



  • 真是太谢谢您了,不过排除了案例本身设置的问题外,从物理直观上不好理解同样求解不可压缩流场,simpleFoam不考虑重力,DPMFoam流体动量方程考虑重力就不合理。


  • 管理员

    0_1523837840346_2018-04-16 08-15-57屏幕截图.png

    pEqn.H改为p_rghEq.H之后,初步看起来是合理的,上图左侧为一个渐进的,下面越来越大的水压,右侧速度为0(1e-10)。这只一种单相附加重力的静态模拟。不过还需要进一步测试。

    http://dyfluid.com/buoyantPimpleFoam.html 我在这里添加了点内容,跟这个有关



  • @东岳 参考您提供的网址 http://dyfluid.com/buoyantPimpleFoam.html 中的部分解析,下面这个方程替换若应用于DPMFoam求解器中不可压缩流场的求解,动量方程中等效于不考虑重力

    0_1523842729945_捕获.PNG

    您提到“OpenFOAM中并没有附加重力的单相流求解器”,显然“buoyantPimpleFoam是OpenFOAM的传热求解器之一,其用于求解瞬态的、由于温度变化导致的密度变化、浮力驱动流动。”,受重力(浮力)驱动。
    物理上讲,不管流体可压与否,重力做功会对竖直方向上的流动产生影响。最好还是能够了解到DPMFoam求解器中流场代码植入之所以是现在这个样子,它的设计思想及应用范围,目前测试来看,在求解不可压缩单相流场时,至少它不是通用的。


  • 管理员

    谢谢。

    OpenFOAM中并没有附加重力的单相流求解器

    已更改为

    OpenFOAM中并没有附加重力的恒密度单相流求解器。例如在恒密度的icoFoam和simpleFoam中,是压力差引起的流体流动。

    最好还是能够了解到DPMFoam求解器中流场代码植入之所以是现在这个样子,它的设计思想及应用范围,目前测试来看,在求解不可压缩单相流场时,至少它不是通用的。

    DPMFoam主要用于气固流动,连续相密度远小于固相,空气的密度无法引致这么大的压力差。这不同于气液(连续相密度大于离散相)。


登录后回复
 

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