applications/solvers/multiphase/bubbleFoam/wallViscosity.H

   1 {
   2     scalar Cmu25 = ::pow(Cmu, 0.25);
   3
   4     const fvPatchList& patches = mesh.boundary();
   5
   6     forAll(patches, patchi)
   7     {
   8         const fvPatch& currPatch = patches[patchi];
   9
  10         if (isType<wallFvPatch>(currPatch))
  11         {
  12             scalarField& nutw = nutb.boundaryField()[patchi];
  13
  14             forAll(currPatch, facei)
  15             {
  16                 label faceCelli = currPatch.faceCells()[facei];
  17
  18                 // calculate yPlus
  19                 scalar yPlus =
  20                     Cmu25*y[patchi][facei]
  21                    *::sqrt(k[faceCelli])
  22                    /nub.value();
  23
  24                 if (yPlus > 11.6)
  25                 {
  26                     nutw[facei] =
  27                         yPlus*nub.value()*kappa
  28                        /::log(E*yPlus)
  29                       - nub.value();
  30                 }
  31                 else
  32                 {
  33                     nutw[facei] = 0.0;
  34                 }
  35             }
  36         }
  37     }
  38 }