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
Foam::sqr
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Definition: dimensionedSymmTensor.C:49
k
label k
Boltzmann constant.
Definition: LISASMDCalcMethod2.H:41
Foam::Ek
tmp< scalarField > Ek(const scalar Ea, const scalar k0, const scalarField &k)
Definition: Ek.H:10
Foam::exp
dimensionedScalar exp(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:274
Foam::pow
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Definition: dimensionedScalar.C:73
Foam::tmp
A class for managing temporary objects.
Definition: PtrList.H:53
Ea
scalar Ea(const scalar p, const scalar T) const
Definition: HtoEthermo.H:20
Foam
Namespace for OpenFOAM.
Definition: atmBoundaryLayer.H:213