applications/solvers/multiphase/twoPhaseEulerFoam/liftDragCoeffs.H

   1     volVectorField Ur(Ua - Ub);
   2     volScalarField magUr(mag(Ur));
   3
   4     volScalarField Ka(draga->K(magUr));
   5     volScalarField K(Ka);
   6
   7     if (dragPhase == "b")
   8     {
   9         volScalarField Kb(dragb->K(magUr));
  10         K = Kb;
  11     }
  12     else if (dragPhase == "blended")
  13     {
  14         volScalarField Kb(dragb->K(magUr));
  15         K = (beta*Ka + alpha*Kb);
  16     }
  17
  18     volVectorField liftCoeff(Cl*(beta*rhob + alpha*rhoa)*(Ur ^ fvc::curl(U)));