applications/solvers/multiphase/bubbleFoam/createRASTurbulence.H

   1     IOdictionary RASProperties
   2     (
   3         IOobject
   4         (
   5             "RASProperties",
   6             runTime.constant(),
   7             mesh,
   8             IOobject::MUST_READ,
   9             IOobject::NO_WRITE
  10         )
  11     );
  12
  13
  14     Switch turbulence
  15     (
  16         RASProperties.lookup("turbulence")
  17     );
  18
  19     dictionary kEpsilonCoeffs
  20     (
  21         RASProperties.subDict("kEpsilonCoeffs")
  22     );
  23
  24
  25     scalar Cmu
  26     (
  27         readScalar(kEpsilonCoeffs.lookup("Cmu"))
  28     );
  29
  30     scalar C1
  31     (
  32         readScalar(kEpsilonCoeffs.lookup("C1"))
  33     );
  34
  35     scalar C2
  36     (
  37         readScalar(kEpsilonCoeffs.lookup("C2"))
  38     );
  39
  40     scalar alphak
  41     (
  42         readScalar(kEpsilonCoeffs.lookup("alphak"))
  43     );
  44
  45     scalar alphaEps
  46     (
  47         readScalar(kEpsilonCoeffs.lookup("alphaEps"))
  48     );
  49
  50     dictionary wallFunctionCoeffs
  51     (
  52         RASProperties.subDict("wallFunctionCoeffs")
  53     );
  54
  55     scalar kappa
  56     (
  57         readScalar(wallFunctionCoeffs.lookup("kappa"))
  58     );
  59
  60     scalar E
  61     (
  62         readScalar(wallFunctionCoeffs.lookup("E"))
  63     );
  64
  65     nearWallDist y(mesh);
  66
  67
  68     Info<< "Reading field k\n" << endl;
  69     volScalarField k
  70     (
  71         IOobject
  72         (
  73             "k",
  74             runTime.timeName(),
  75             mesh,
  76             IOobject::MUST_READ,
  77             IOobject::AUTO_WRITE
  78         ),
  79         mesh
  80     );
  81
  82     Info<< "Reading field epsilon\n" << endl;
  83     volScalarField epsilon
  84     (
  85         IOobject
  86         (
  87             "epsilon",
  88             runTime.timeName(),
  89             mesh,
  90             IOobject::MUST_READ,
  91             IOobject::AUTO_WRITE
  92         ),
  93         mesh
  94     );
  95
  96
  97     Info<< "Calculating field nutb\n" << endl;
  98     volScalarField nutb
  99     (
 100         IOobject
 101         (
 102             "nutb",
 103             runTime.timeName(),
 104             mesh,
 105             IOobject::NO_READ,
 106             IOobject::AUTO_WRITE
 107         ),
 108         Cmu*sqr(k)/epsilon
 109     );
 110
 111     Info<< "Calculating field nuEffa\n" << endl;
 112     volScalarField nuEffa
 113     (
 114         IOobject
 115         (
 116             "nuEffa",
 117             runTime.timeName(),
 118             mesh,
 119             IOobject::NO_READ,
 120             IOobject::NO_WRITE
 121         ),
 122         sqr(Ct)*nutb + nua
 123     );
 124
 125     Info<< "Calculating field nuEffb\n" << endl;
 126     volScalarField nuEffb
 127     (
 128         IOobject
 129         (
 130             "nuEffb",
 131             runTime.timeName(),
 132             mesh,
 133             IOobject::NO_READ,
 134             IOobject::NO_WRITE
 135         ),
 136         nutb + nub
 137     );