想用双欧拉方法模拟水气两相流，请问twoPhaseEulerFoam中可以加入PBM模型来模拟气泡的聚合和破碎吗？该如何操作呢

@李东岳 算的固液两相流，李老师，我现在在做测试，就想搞清楚拉格朗日这面

@李东岳 李老师您好，还想请教一下，因为我的粒子直径相对于网格来说较大，所以用DPMFoam计算时发散了，取了几个时间步后的压力分布图发现有一个点压强突然变的很小，如下图：

我的粒子浓度相对来说应该是比较大的

但当我把kinematicCloudProperties里的alphamax从0.9变成0.7之后就可以正常计算了，请问这是什么原因呢？通过调整alphamax来进行计算可以吗？谢谢您

@李东岳 谢谢李老师！

@李东岳 但是我的模拟中粒子是通过kinematicCloudPositions来定义的，所以计算域内粒子数量是固定的。而且我的网格大小是6mm，粒子直径是3mm，对于这种情况是不是网格尺寸对粒子浓度影响较大啊？

各位大佬好，
本人最近在想是不是拉格朗日粒子的浓度是与网格大小相关的？因为在openfoam中浓度定义为：

theta[celli] += p.nParticle()*p.volume();
theta.primitiveFieldRef() /= mesh_.V();

当网格变小的时候，浓度不就自然而然变大了吗？如果是这样的话，该如何跟实验进行验证呢？

@zhe 在 MPPICInterFoam 中如何通过alphac得到alpha1和alpha2? 中说：

那么alpha1才应该是整个cell里的alpha1啊

没看懂这是什么意思。举个例子，比如alphac=1-theta=0.8, alpha1=0.5，那么phase1所占整个cell的体积分数就是0.8*0.5=0.4，phase2也是0.4，0.4+0.4+0.2=1，在MPPICInterFoam中粒子是考虑了体积分数的

大家好，我想请教一下调用momentum source函数的时候，返回值是哪一个？在KinematicCloudI.H中有这个函数的定义：

template<class CloudType>
 inline Foam::tmp<Foam::fvVectorMatrix>
 Foam::KinematicCloud<CloudType>::SU(volVectorField& U, bool incompressible)
 const
 {
     if (debug)
     {
         Pout<< "UTrans min/max = " << min(UTrans()).value() << ", "
             << max(UTrans()).value() << nl
             << "UCoeff min/max = " << min(UCoeff()).value() << ", "
             << max(UCoeff()).value() << endl;
     }
  
     dimensionSet dim(dimForce);
     if (incompressible)
     {
         dim.reset(dimForce/dimDensity);
     }
  
     if (solution_.coupled())
     {
         if (solution_.semiImplicit("U"))
         {
             volScalarField::Internal
                 Vdt(mesh_.V()*this->db().time().deltaT());
  
             if (incompressible)
             {
                 Vdt.dimensions() *= dimDensity;
             }
  
             return UTrans()/Vdt - fvm::Sp(UCoeff()/Vdt, U) + UCoeff()/Vdt*U;
         }
         else
         {
             tmp<fvVectorMatrix> tfvm(new fvVectorMatrix(U, dim));
             fvVectorMatrix& fvm = tfvm.ref();
  
             fvm.source() = -UTrans()/(this->db().time().deltaT());
  
             return tfvm;
         }
     }
  
     return tmp<fvVectorMatrix>::New(U, dim);
 }

在调用这个函数的时候，如果我用semiImplicit格式的话，返回的是

return UTrans()/Vdt - fvm::Sp(UCoeff()/Vdt, U) + UCoeff()/Vdt*U;

还是

return tmp<fvVectorMatrix>::New(U, dim);

呢？
我认为返回的应该是后者吧 tmp<fvVectorMatrix>::New(U, dim);
我想知道它是什么意思呢？

@bestucan 好的，非常感谢大佬的指点！