applications/solvers/multiphase/bubbleFoam/wallFunctions.H

   1 {
   2     labelList cellBoundaryFaceCount(epsilon.size(), 0);
   3
   4     scalar Cmu25 = ::pow(Cmu.value(), 0.25);
   5     scalar Cmu75 = ::pow(Cmu.value(), 0.75);
   6     scalar kappa_ = kappa.value();
   7
   8     const fvPatchList& patches = mesh.boundary();
   9
  10     //- Initialise the near-wall P field to zero
  11     forAll(patches, patchi)
  12     {
  13         const fvPatch& currPatch = patches[patchi];
  14
  15         if (isType<wallFvPatch>(currPatch))
  16         {
  17             forAll(currPatch, facei)
  18             {
  19                 label faceCelli = currPatch.faceCells()[facei];
  20
  21                 epsilon[faceCelli] = 0.0;
  22                 G[faceCelli] = 0.0;
  23             }
  24         }
  25     }
  26
  27     //- Accumulate the wall face contributions to epsilon and G
  28     //  Increment cellBoundaryFaceCount for each face for averaging
  29     forAll(patches, patchi)
  30     {
  31         const fvPatch& currPatch = patches[patchi];
  32
  33         if (isType<wallFvPatch>(currPatch))
  34         {
  35             const scalarField& nuw = nutb.boundaryField()[patchi];
  36
  37             scalarField magFaceGradU = mag(U.boundaryField()[patchi].snGrad());
  38
  39             forAll(currPatch, facei)
  40             {
  41                 label faceCelli = currPatch.faceCells()[facei];
  42
  43                 scalar yPlus =
  44                     Cmu25*y[patchi][facei]
  45                     *::sqrt(k[faceCelli])
  46                     /nub.value();
  47
  48
  49                 // For corner cells (with two boundary or more faces),
  50                 // epsilon and G in the near-wall cell are calculated
  51                 // as an average
  52
  53                 cellBoundaryFaceCount[faceCelli]++;
  54
  55                 epsilon[faceCelli] +=
  56                      Cmu75*::pow(k[faceCelli], 1.5)
  57                     /(kappa_*y[patchi][facei]);
  58
  59                 if (yPlus > 11.6)
  60                 {
  61                     G[faceCelli] +=
  62                         nuw[facei]*magFaceGradU[facei]
  63                         *Cmu25*::sqrt(k[faceCelli])
  64                         /(kappa_*y[patchi][facei]);
  65                 }
  66             }
  67         }
  68     }
  69
  70
  71     // perform the averaging
  72
  73     forAll(patches, patchi)
  74     {
  75         const fvPatch& curPatch = patches[patchi];
  76
  77         if (isType<wallFvPatch>(curPatch))
  78         {
  79             forAll(curPatch, facei)
  80             {
  81                 label faceCelli = curPatch.faceCells()[facei];
  82
  83                 epsilon[faceCelli] /= cellBoundaryFaceCount[faceCelli];
  84                 G[faceCelli] /= cellBoundaryFaceCount[faceCelli];
  85             }
  86         }
  87     }
  88 }