DPMFoam计算粒子分布



  • 在对DPMFoam自带的例子Goldschmidt修改后出现的几个问题希望大佬们看一下:chouchou:
    1、修改了turbulenceProperties.air,将其内部的层流模型改成了湍流LES模型,粒子的速度为什么没有变化,
    只有流体的速度发生了变化,DPM不是四向耦合的吗;
    2、修改了流体性质为水的性质,计算出来粒子速度非常大,有些达到了50m/s,这个量级是对的吗?



  • 上面的问题有点想明白了,现在又遇到了新问题,因为这是液固系统必须要考虑虚拟质量力(这应该也是为什么粒子速度非常大的原因),但我在添加virtualMass之后,出现了下面的报错,进源文件里面看也没有用到U的地方吧,烦劳了解的老师同学解答一下,不胜感激:xinxin:
    /home/lx/图片/2019-10-26 15-40-37屏幕截图.png



  • 怎么上传不了截图,好吧,直接上代码吧

      request for volVectorField U from objectRegistry region0 failed
        available objects of type volVectorField are 1(U.air)
    
        From function const Type& Foam::objectRegistry::lookupObject(const Foam::word&) const [with Type = Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>]
        in file /home/lx/OpenFOAM/OpenFOAM-5.x/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 193.
    


  • 你在设置虚拟质量立的时候需要指定关键词

    U     U.air;
    

    否则默认读取的是U,它提示你

    request for volVectorField U from objectRegistry region0 failed
        available objects of type volVectorField are 1(U.air)
    


  • @东岳 谢谢东岳老师,我刚刚也发现了,好像求虚拟质量力时候需要先求解pressureGradient这个力,求解这个力的时候需要用到U再次感谢东岳老师:xinxin: :xinxin:



  • 你好,我想知道你这个虚拟质量是怎么加的,方便透露一下吗?最好能截图展示一下,万分感谢!



  • @ZYj

      interpolationSchemes//添加离散格式
       {
           rho.air         cell;
           U.air           cellPoint;
           mu.air          cell;
          DUcDt             cellPoint;
       }
    
      pressureGradient//添加虚拟质量力
            {
                U       U.air;
           } 
           virtualMass
           {
           Cvm   0.5;
           }
    

Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2