DPMFoam计算粒子分布

upc_ngh

在对DPMFoam自带的例子Goldschmidt修改后出现的几个问题希望大佬们看一下：
１、修改了turbulenceProperties.air，将其内部的层流模型改成了湍流LES模型，粒子的速度为什么没有变化，
只有流体的速度发生了变化，ＤＰＭ不是四向耦合的吗；
２、修改了流体性质为水的性质，计算出来粒子速度非常大，有些达到了５０m/s，这个量级是对的吗？

upc_ngh

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

upc_ngh

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

　　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)

upc_ngh

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

ZYj

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

upc_ngh

@ZYj

  interpolationSchemes//添加离散格式
   {
       rho.air         cell;
       U.air           cellPoint;
       mu.air          cell;
      DUcDt             cellPoint;
   }

  pressureGradient//添加虚拟质量力
        {
            U       U.air;
       } 
       virtualMass
       {
       Cvm   0.5;
       }