src/randomProcesses/turbulence/Ek.H

   1 #ifndef Ek_H
   2 #define Ek_H
   3
   4 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
   5
   6 namespace Foam
   7 {
   8
   9 inline tmp<scalarField> Ek
  10 (
  11     const scalar Ea,
  12     const scalar k0,
  13     const scalarField& k
  14 )
  15 {
  16     tmp<scalarField> tEk = Ea*pow(k/k0, 4.0)*exp(-2.0*sqr(k/k0));
  17
  18     /*
  19     scalarField& Ekf = tEk();
  20
  21     label i;
  22     forAll(Ekf, i)
  23     {
  24         if (k[i] < 2 || k[i] > 10)
  25         {
  26             Ekf[i] = 0.0;
  27         }
  28     }
  29     */
  30
  31     return tEk;
  32 }
  33
  34
  35 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  36
  37 } // End namespace Foam
  38
  39 #endif