Ek.H
Go to the documentation of this file.
1#ifndef Ek_H
2#define Ek_H
3
4// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
5
6namespace Foam
7{
8
9inline 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
Ea
scalar Ea(const scalar p, const scalar T) const
Definition: HtoEthermo.H:32
k
label k
Definition: LISASMDCalcMethod2.H:41
Foam::tmp
A class for managing temporary objects.
Definition: tmp.H:65
Foam
Namespace for OpenFOAM.
Definition: atmBoundaryLayer.C:34
Foam::exp
dimensionedScalar exp(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:261
Foam::sqr
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Definition: dimensionedSymmTensor.C:51
Foam::pow
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Definition: dimensionedScalar.C:75
Foam::Ek
tmp< scalarField > Ek(const scalar Ea, const scalar k0, const scalarField &k)
Definition: Ek.H:10