applications/solvers/compressible/rhoCentralFoam/createFields.H

   1 Info<< "Reading thermophysical properties\n" << endl;
   2
   3 autoPtr<basicPsiThermo> pThermo
   4 (
   5     basicPsiThermo::New(mesh)
   6 );
   7 basicPsiThermo& thermo = pThermo();
   8
   9 volScalarField& p = thermo.p();
  10 volScalarField& e = thermo.e();
  11 const volScalarField& T = thermo.T();
  12 const volScalarField& psi = thermo.psi();
  13 const volScalarField& mu = thermo.mu();
  14
  15 bool inviscid(true);
  16 if (max(mu.internalField()) > 0.0)
  17 {
  18     inviscid = false;
  19 }
  20
  21 Info<< "Reading field U\n" << endl;
  22 volVectorField U
  23 (
  24     IOobject
  25     (
  26         "U",
  27         runTime.timeName(),
  28         mesh,
  29         IOobject::MUST_READ,
  30         IOobject::AUTO_WRITE
  31     ),
  32     mesh
  33 );
  34
  35 #include "rhoBoundaryTypes.H"
  36 volScalarField rho
  37 (
  38     IOobject
  39     (
  40         "rho",
  41         runTime.timeName(),
  42         mesh,
  43         IOobject::NO_READ,
  44         IOobject::AUTO_WRITE
  45     ),
  46     thermo.rho(),
  47     rhoBoundaryTypes
  48 );
  49
  50 volVectorField rhoU
  51 (
  52     IOobject
  53     (
  54         "rhoU",
  55         runTime.timeName(),
  56         mesh,
  57         IOobject::NO_READ,
  58         IOobject::NO_WRITE
  59     ),
  60     rho*U
  61 );
  62
  63 volScalarField rhoE
  64 (
  65     IOobject
  66     (
  67         "rhoE",
  68         runTime.timeName(),
  69         mesh,
  70         IOobject::NO_READ,
  71         IOobject::NO_WRITE
  72     ),
  73     rho*(e + 0.5*magSqr(U))
  74 );
  75
  76 surfaceScalarField pos
  77 (
  78     IOobject
  79     (
  80         "pos",
  81         runTime.timeName(),
  82         mesh
  83     ),
  84     mesh,
  85     dimensionedScalar("pos", dimless, 1.0)
  86 );
  87
  88 surfaceScalarField neg
  89 (
  90     IOobject
  91     (
  92         "neg",
  93         runTime.timeName(),
  94         mesh
  95     ),
  96     mesh,
  97     dimensionedScalar("neg", dimless, -1.0)
  98 );
  99
 100
 101 surfaceScalarField phi("phi", mesh.Sf() & fvc::interpolate(rhoU));
 102
 103 Info<< "Creating turbulence model\n" << endl;
 104 autoPtr<compressible::turbulenceModel> turbulence
 105 (
 106     compressible::turbulenceModel::New
 107     (
 108         rho,
 109         U,
 110         phi,
 111         thermo
 112     )
 113 );