applications/solvers/lagrangian/reactingParcelFoam/createFields.H

   1     Info<< "Reading thermophysical properties\n" << endl;
   2
   3     autoPtr<psiChemistryModel> pChemistry
   4     (
   5         psiChemistryModel::New(mesh)
   6     );
   7     psiChemistryModel& chemistry = pChemistry();
   8
   9     hsCombustionThermo& thermo = chemistry.thermo();
  10
  11     SLGThermo slgThermo(mesh, thermo);
  12
  13     basicMultiComponentMixture& composition = thermo.composition();
  14     PtrList<volScalarField>& Y = composition.Y();
  15
  16     const word inertSpecie(thermo.lookup("inertSpecie"));
  17
  18     if (!composition.contains(inertSpecie))
  19     {
  20         FatalErrorIn(args.executable())
  21             << "Specified inert specie '" << inertSpecie << "' not found in "
  22             << "species list. Available species:" << composition.species()
  23             << exit(FatalError);
  24     }
  25
  26     volScalarField& p = thermo.p();
  27     volScalarField& hs = thermo.hs();
  28     const volScalarField& T = thermo.T();
  29     const volScalarField& psi = thermo.psi();
  30
  31     volScalarField rho
  32     (
  33         IOobject
  34         (
  35             "rho",
  36             runTime.timeName(),
  37             mesh,
  38             IOobject::NO_READ,
  39             IOobject::AUTO_WRITE
  40         ),
  41         thermo.rho()
  42     );
  43
  44     Info<< "\nReading field U\n" << endl;
  45     volVectorField U
  46     (
  47         IOobject
  48         (
  49             "U",
  50             runTime.timeName(),
  51             mesh,
  52             IOobject::MUST_READ,
  53             IOobject::AUTO_WRITE
  54         ),
  55         mesh
  56     );
  57
  58     #include "compressibleCreatePhi.H"
  59
  60     DimensionedField<scalar, volMesh> kappa
  61     (
  62         IOobject
  63         (
  64             "kappa",
  65             runTime.timeName(),
  66             mesh,
  67             IOobject::NO_READ,
  68             IOobject::AUTO_WRITE
  69         ),
  70         mesh,
  71         dimensionedScalar("zero", dimless, 0.0)
  72     );
  73
  74     Info<< "Creating turbulence model\n" << endl;
  75     autoPtr<compressible::turbulenceModel> turbulence
  76     (
  77         compressible::turbulenceModel::New
  78         (
  79             rho,
  80             U,
  81             phi,
  82             thermo
  83         )
  84     );
  85
  86     Info<< "Creating field DpDt\n" << endl;
  87     volScalarField DpDt
  88     (
  89         "DpDt",
  90         fvc::DDt(surfaceScalarField("phiU", phi/fvc::interpolate(rho)), p)
  91     );
  92
  93     multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields;
  94
  95     forAll(Y, i)
  96     {
  97         fields.add(Y[i]);
  98     }
  99     fields.add(hs);
 100
 101     DimensionedField<scalar, volMesh> chemistrySh
 102     (
 103         IOobject
 104         (
 105             "chemistry::Sh",
 106             runTime.timeName(),
 107             mesh,
 108             IOobject::NO_READ,
 109             IOobject::NO_WRITE
 110         ),
 111         mesh,
 112         dimensionedScalar("chemistry::Sh", dimEnergy/dimTime/dimVolume, 0.0)
 113     );