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));