颗粒曳力的孔隙率alphac是怎么读取的(插值计算还是网格中心值)？

hotairballon1997

ErgunWenYuDragForce.C中，计算颗粒曳力的代码如下：

     scalar alphac(alphac_[p.cell()]);
  
     if (alphac < 0.8)
     {
         return forceSuSp
         (
             Zero,
             (mass/p.rho())
            *(150.0*(1.0 - alphac)/alphac + 1.75*Re)*muc/(alphac*sqr(p.d()))
         );
     }
     else
     {
         return forceSuSp
         (
             Zero,
             (mass/p.rho())
            *0.75*CdRe(alphac*Re)*muc*pow(alphac, -2.65)/(alphac*sqr(p.d()))
         );
     }

这里第一行的alphac(alphac_[p.cell()])，是插值计算网格的alphac值(cellPoint方式)，还是只读取网格中心值(cell方式)？因为没有查到关于这个[p.cell()]的解释，所以想请教一下论坛里的前辈

李东岳

@hotairballon1997 p.cell()是拉格朗日粒子所在单元的标号，那就是你说的只读取网格中心值(cell方式)

hotairballon1997

@李东岳谢谢李老师解答。

那是不是可以理解为，在比如下图中（打开颗粒对流体的影响），颗粒的中心点在网格1，但是颗粒的实际大小有如图显示那么大。现在就只有网格1受到颗粒的影响，而2,3,4号网格都完全感受不到颗粒的存在？

不知道这对总体结果会造成多大的影响？

李东岳

应该影响不大。这种颗粒应该不多，不会所有的颗粒都恰好那样。所以对alpha的影响也很小。除非颗粒直径非常大。不过直径非常大的话，那就会引起收敛性的问题。