applications/solvers/multiphase/settlingFoam/calcVdj.H

   1 if (VdjModel == "general")
   2 {
   3     Vdj = V0*
   4     (
   5         exp(-a*max(alpha - alphaMin, scalar(0)))
   6       - exp(-a1*max(alpha - alphaMin, scalar(0)))
   7     );
   8 }
   9 else if (VdjModel == "simple")
  10 {
  11     Vdj = V0*pow(10.0, -a*alpha);
  12 }
  13 else
  14 {
  15     FatalErrorIn(args.executable())
  16         << "Unknown VdjModel : " << VdjModel
  17         << abort(FatalError);
  18 }
  19
  20 Vdj.correctBoundaryConditions();