applications/solvers/discreteMethods/dsmc/dsmcFoam/createFields.H

   1
   2     Info<< nl << "Reading field boundaryT" << endl;
   3     volScalarField boundaryT
   4     (
   5         IOobject
   6         (
   7             "boundaryT",
   8             runTime.timeName(),
   9             mesh,
  10             IOobject::MUST_READ,
  11             IOobject::AUTO_WRITE
  12         ),
  13         mesh
  14     );
  15
  16     Info<< nl << "Reading field boundaryU" << endl;
  17     volVectorField boundaryU
  18     (
  19         IOobject
  20         (
  21             "boundaryU",
  22             runTime.timeName(),
  23             mesh,
  24             IOobject::MUST_READ,
  25             IOobject::AUTO_WRITE
  26         ),
  27         mesh
  28     );
  29
  30     Info<< nl << "Reading field rhoN (number density)" << endl;
  31     volScalarField rhoN
  32     (
  33         IOobject
  34         (
  35             "rhoN",
  36             runTime.timeName(),
  37             mesh,
  38             IOobject::MUST_READ,
  39             IOobject::AUTO_WRITE
  40         ),
  41         mesh
  42     );
  43
  44     Info<< nl << "Reading field rhoM (mass density)" << endl;
  45     volScalarField rhoM
  46     (
  47         IOobject
  48         (
  49             "rhoM",
  50             runTime.timeName(),
  51             mesh,
  52             IOobject::MUST_READ,
  53             IOobject::AUTO_WRITE
  54         ),
  55         mesh
  56     );
  57
  58     Info<< nl << "Reading field rhoNdsmc (dsmc particle density)" << endl;
  59     volScalarField dsmcRhoN
  60     (
  61         IOobject
  62         (
  63             "dsmcRhoN",
  64             runTime.timeName(),
  65             mesh,
  66             IOobject::MUST_READ,
  67             IOobject::AUTO_WRITE
  68         ),
  69         mesh
  70     );
  71
  72     Info<< nl << "Reading field momentum (momentum density)" << endl;
  73     volVectorField momentum
  74     (
  75         IOobject
  76         (
  77             "momentum",
  78             runTime.timeName(),
  79             mesh,
  80             IOobject::MUST_READ,
  81             IOobject::AUTO_WRITE
  82         ),
  83         mesh
  84     );
  85
  86     Info<< nl << "Reading field linearKE (linear kinetic energy density)"
  87         << endl;
  88
  89     volScalarField linearKE
  90     (
  91         IOobject
  92         (
  93             "linearKE",
  94             runTime.timeName(),
  95             mesh,
  96             IOobject::MUST_READ,
  97             IOobject::AUTO_WRITE
  98         ),
  99         mesh
 100     );
 101
 102     Info<< nl << "Reading field internalE (internal energy density)" << endl;
 103     volScalarField internalE
 104     (
 105         IOobject
 106         (
 107             "internalE",
 108             runTime.timeName(),
 109             mesh,
 110             IOobject::MUST_READ,
 111             IOobject::AUTO_WRITE
 112         ),
 113         mesh
 114     );
 115
 116     Info<< nl << "Reading field iDof (internal degree of freedom density)"
 117         << endl;
 118
 119     volScalarField iDof
 120     (
 121         IOobject
 122         (
 123             "iDof",
 124             runTime.timeName(),
 125             mesh,
 126             IOobject::MUST_READ,
 127             IOobject::AUTO_WRITE
 128         ),
 129         mesh
 130     );
 131
 132     Info<< nl << "Reading field q (surface heat transfer)" << endl;
 133     volScalarField q
 134     (
 135         IOobject
 136         (
 137             "q",
 138             runTime.timeName(),
 139             mesh,
 140             IOobject::MUST_READ,
 141             IOobject::AUTO_WRITE
 142         ),
 143         mesh
 144     );
 145
 146     Info<< nl << "Reading field fD (surface force density)" << endl;
 147     volVectorField fD
 148     (
 149         IOobject
 150         (
 151             "fD",
 152             runTime.timeName(),
 153             mesh,
 154             IOobject::MUST_READ,
 155             IOobject::AUTO_WRITE
 156         ),
 157         mesh
 158     );
 159
 160     Info<< nl << "Constructing dsmcCloud " << endl;
 161
 162     dsmcCloud dsmc("dsmc", boundaryT, boundaryU);