Ek.H
Go to the documentation of this file.
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
dimensionedSymmTensor sqr(const dimensionedVector &dv)
label k
Boltzmann constant.
tmp< scalarField > Ek(const scalar Ea, const scalar k0, const scalarField &k)
Definition: Ek.H:10
dimensionedScalar exp(const dimensionedScalar &ds)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
A class for managing temporary objects.
Definition: PtrList.H:54
Namespace for OpenFOAM.