FreeFOAM The Cross-Platform CFD Toolkit
kEpsilon.H
Go to the documentation of this file.
2 {
3  if (mesh.changing())
4  {
5  y.correct();
6  }
7 
8  dimensionedScalar k0("k0", k.dimensions(), SMALL);
9  dimensionedScalar epsilon0("epsilon0", epsilon.dimensions(), SMALL);
10 
12 
13  tmp<volTensorField> tgradU = fvc::grad(U);
14  volScalarField G = 2*mut*(tgradU() && dev(symm(tgradU())));
15  tgradU.clear();
16 
17  volScalarField Gcoef =
18  alphak*Cmu*k*(g & fvc::grad(rho))/(epsilon + epsilon0);
19 
20 # include "wallFunctions.H"
21 
22  // Dissipation equation
23  fvScalarMatrix epsEqn
24  (
25  fvm::ddt(rho, epsilon)
26  + fvm::div(phi, epsilon)
28  (
29  alphaEps*mut + mul, epsilon,
30  "laplacian(DepsilonEff,epsilon)"
31  )
32  ==
33  C1*G*epsilon/k
34  - fvm::SuSp(C1*(1.0 - C3)*Gcoef + (2.0/3.0*C1)*rho*divU, epsilon)
35  - fvm::Sp(C2*rho*epsilon/k, epsilon)
36  );
37 
38 # include "wallDissipation.H"
39 
40  epsEqn.relax();
41  epsEqn.solve();
42 
43  bound(epsilon, epsilon0);
44 
45 
46  // Turbulent kinetic energy equation
47 
48  fvScalarMatrix kEqn
49  (
50  fvm::ddt(rho, k)
51  + fvm::div(phi, k)
53  (
54  alphak*mut + mul, k,
55  "laplacian(DkEff,k)"
56  )
57  ==
58  G
59  - fvm::SuSp(Gcoef + 2.0/3.0*rho*divU, k)
60  - fvm::Sp(rho*epsilon/k, k)
61  );
62 
63  kEqn.relax();
64  kEqn.solve();
65 
66  bound(k, k0);
67 
68 
69  //- Re-calculate viscosity
70  mut = rho*Cmu*sqr(k)/epsilon;
71 
72 # include "wallViscosity.H"
73 }
74 
75 mu = mut + mul;
76 
77 // ************************ vim: set sw=4 sts=4 et: ************************ //