applications/solvers/multiphase/bubbleFoam/liftDragCoeffs.H

   1     volVectorField Ur = Ua - Ub;
   2     volScalarField magUr = mag(Ur);
   3
   4     volScalarField CdaMagUr =
   5         (24.0*nub/da)*(scalar(1) + 0.15*pow(da*magUr/nub, 0.687));
   6
   7     volScalarField CdbMagUr =
   8         (24.0*nua/db)*(scalar(1) + 0.15*pow(db*magUr/nua, 0.687));
   9
  10     volScalarField dragCoef
  11     (
  12         "Cd",
  13         0.75*(beta*rhob*CdaMagUr/da + alpha*rhoa*CdbMagUr/db)
  14     );
  15
  16     volVectorField liftCoeff = Cl*(beta*rhob + alpha*rhoa)*(Ur ^ fvc::curl(U));