applications/solvers/multiphase/settlingFoam/correctViscosity.H

   1 {
   2     mul = muc +
   3         plasticViscosity
   4         (
   5             plasticViscosityCoeff,
   6             plasticViscosityExponent,
   7             alpha
   8         );
   9
  10     if (BinghamPlastic)
  11     {
  12         volScalarField tauy = yieldStress
  13         (
  14             yieldStressCoeff,
  15             yieldStressExponent,
  16             yieldStressOffset,
  17             alpha
  18         );
  19
  20         mul =
  21             tauy/
  22             (
  23                 mag(fvc::grad(U))
  24               + 1.0e-4*
  25                 (
  26                     tauy
  27                   + dimensionedScalar
  28                     (
  29                         "deltaTauy",
  30                         tauy.dimensions(),
  31                         1.0e-15
  32                     )
  33                 )/mul
  34             )
  35           + mul;
  36     }
  37 }