applications/solvers/heatTransfer/chtMultiRegionFoam/solid/createSolidFields.H

   1     // Initialise solid field pointer lists
   2     PtrList<volScalarField> rhos(solidRegions.size());
   3     PtrList<volScalarField> cps(solidRegions.size());
   4     PtrList<volScalarField> rhosCps(solidRegions.size());
   5     PtrList<volScalarField> Ks(solidRegions.size());
   6     PtrList<volScalarField> Ts(solidRegions.size());
   7
   8     // Populate solid field pointer lists
   9     forAll(solidRegions, i)
  10     {
  11         Info<< "*** Reading solid mesh thermophysical properties for region "
  12             << solidRegions[i].name() << nl << endl;
  13
  14         Info<< "    Adding to rhos\n" << endl;
  15         rhos.set
  16         (
  17             i,
  18             new volScalarField
  19             (
  20                 IOobject
  21                 (
  22                     "rho",
  23                     runTime.timeName(),
  24                     solidRegions[i],
  25                     IOobject::MUST_READ,
  26                     IOobject::AUTO_WRITE
  27                 ),
  28                 solidRegions[i]
  29             )
  30         );
  31
  32         Info<< "    Adding to cps\n" << endl;
  33         cps.set
  34         (
  35             i,
  36             new volScalarField
  37             (
  38                 IOobject
  39                 (
  40                     "cp",
  41                     runTime.timeName(),
  42                     solidRegions[i],
  43                     IOobject::MUST_READ,
  44                     IOobject::AUTO_WRITE
  45                 ),
  46                 solidRegions[i]
  47             )
  48         );
  49
  50         rhosCps.set
  51         (
  52             i,
  53             new volScalarField("rhosCps", rhos[i]*cps[i])
  54         );
  55
  56         Info<< "    Adding to Ks\n" << endl;
  57         Ks.set
  58         (
  59             i,
  60             new volScalarField
  61             (
  62                 IOobject
  63                 (
  64                     "K",
  65                     runTime.timeName(),
  66                     solidRegions[i],
  67                     IOobject::MUST_READ,
  68                     IOobject::AUTO_WRITE
  69                 ),
  70                 solidRegions[i]
  71             )
  72         );
  73
  74         Info<< "    Adding to Ts\n" << endl;
  75         Ts.set
  76         (
  77             i,
  78             new volScalarField
  79             (
  80                 IOobject
  81                 (
  82                     "T",
  83                     runTime.timeName(),
  84                     solidRegions[i],
  85                     IOobject::MUST_READ,
  86                     IOobject::AUTO_WRITE
  87                 ),
  88                 solidRegions[i]
  89             )
  90         );
  91     }