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
k
label k
Definition: LISASMDCalcMethod2.H:41
Foam::tmp
A class for managing temporary objects.
Definition: tmp.H:55
Foam
Namespace for OpenFOAM.
Definition: atmosphericBoundaryLayer.C:32
Foam::exp
dimensionedScalar exp(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:274
Foam::sqr
tmp< DimensionedField< typename outerProduct< Type, Type >::type, GeoMesh, Field >> sqr(const DimensionedField< Type, GeoMesh, PrimitiveField > &df)
Definition: DimensionedFieldFunctions.C:134
Foam::pow
tmp< DimensionedField< typename powProduct< Type, r >::type, GeoMesh, Field > > pow(const DimensionedField< Type, GeoMesh, PrimitiveField > &df, typename powProduct< Type, r >::type)
Definition: DimensionedFieldFunctions.C:61
Foam::Ek
tmp< scalarField > Ek(const scalar Ea, const scalar k0, const scalarField &k)
Definition: Ek.H:10