applications/solvers/lagrangian/uncoupledKinematicParcelFoam/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 rho
  10     (
  11         IOobject
  12         (
  13             "rho",
  14             runTime.timeName(),
  15             mesh,
  16             IOobject::NO_READ,
  17             IOobject::AUTO_WRITE
  18         ),
  19         thermo.rho()
  20     );
  21
  22     Info<< "\nReading field U\n" << endl;
  23     volVectorField U
  24     (
  25         IOobject
  26         (
  27             "U",
  28             runTime.timeName(),
  29             mesh,
  30             IOobject::MUST_READ,
  31             IOobject::AUTO_WRITE
  32         ),
  33         mesh
  34     );
  35
  36     #include "compressibleCreatePhi.H"
  37
  38     Info<< "Creating turbulence model\n" << endl;
  39     autoPtr<compressible::turbulenceModel> turbulence
  40     (
  41         compressible::turbulenceModel::New
  42         (
  43             rho,
  44             U,
  45             phi,
  46             thermo
  47         )
  48     );
  49
  50     const word kinematicCloudName
  51     (
  52         args.optionLookupOrDefault<word>("cloudName", "kinematicCloud")
  53     );
  54
  55     Info<< "Constructing kinematicCloud " << kinematicCloudName << endl;
  56     basicKinematicCloud kinematicCloud
  57     (
  58         kinematicCloudName,
  59         rho,
  60         U,
  61         thermo.mu(),
  62         g
  63     );