twoPhaseEulerFoam 的升力项怎么不在 动量方程里

[麦迪文]

不知我下面标的对不对，twoPhaseEulerFoam 里动量方程的升力项怎么不在动量方程中？是从哪里调用的？

         U1Eqn =
         (
             fvm::ddt(alpha1, rho1, U1) + fvm::div(alphaRhoPhi1, U1)   //对流扩散项
           - fvm::Sp(contErr1, U1)   //修正？？？
           + MRF.DDt(alpha1*rho1 + Vm, U1) //多坐标系
           + phase1.turbulence().divDevRhoReff(U1) // 粘度和湍流
          ==
           - Vm
            *(
                 fvm::ddt(U1)
               + fvm::div(phi1, U1)
               - fvm::Sp(fvc::div(phi1), U1)
               - DDtU2
             )  //虚拟质量力
           + fvOptions(alpha1, rho1, U1)
         );
         U1Eqn.relax();
         U1Eqn += fvm::Sp(Kd, U1); //曳力

[麦迪文]

不好意思，太马虎了，升力和壁面润滑力在peqn.h里，尴尬。。。。

     // Lift and wall-lubrication forces
     volVectorField F(fluid.F());

[Dingcy]

小师兄，不研究interFoam，要改行了？

[麦迪文]

@Dingcy interfoam告一段落了，接下来跟你们学习点先进的，整点夹杂物搁进去了

[Dingcy]

@麦迪文 可以可以，小师兄进度很快，我们老师给我博士课题分的要用这个求解器，我最近也刚刚开始做这个，找到小伙伴了

[麦迪文]

@Dingcy 哈哈哈哈，我就瞎看看，我主要做CFDEM方面