applications/solvers/multiphase/bubbleFoam/liftDragCoeffs.H

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