applications/solvers/lagrangian/icoUncoupledKinematicParcelFoam/createFields.H

   1     Info<< "\nReading transportProperties\n" << endl;
   2
   3     IOdictionary transportProperties
   4     (
   5         IOobject
   6         (
   7             "transportProperties",
   8             runTime.constant(),
   9             mesh,
  10             IOobject::MUST_READ_IF_MODIFIED,
  11             IOobject::NO_WRITE
  12         )
  13     );
  14
  15     dimensionedScalar rhoInfValue
  16     (
  17         transportProperties.lookup("rhoInf")
  18     );
  19
  20     volScalarField rhoInf
  21     (
  22         IOobject
  23         (
  24             "rho",
  25             runTime.timeName(),
  26             mesh,
  27             IOobject::NO_READ,
  28             IOobject::AUTO_WRITE
  29         ),
  30         mesh,
  31         rhoInfValue
  32     );
  33
  34     Info<< "Reading field U\n" << endl;
  35     volVectorField U
  36     (
  37         IOobject
  38         (
  39             "U",
  40             runTime.timeName(),
  41             mesh,
  42             IOobject::MUST_READ,
  43             IOobject::AUTO_WRITE
  44         ),
  45         mesh
  46     );
  47
  48     #include "createPhi.H"
  49
  50     Info<< "Creating turbulence model\n" << endl;
  51
  52     singlePhaseTransportModel laminarTransport(U, phi);
  53
  54     const volScalarField nu(laminarTransport.nu());
  55
  56     autoPtr<incompressible::turbulenceModel> turbulence
  57     (
  58         incompressible::turbulenceModel::New(U, phi, laminarTransport)
  59     );
  60
  61     volScalarField mu
  62     (
  63         IOobject
  64         (
  65             "mu",
  66             runTime.timeName(),
  67             mesh,
  68             IOobject::NO_READ,
  69             IOobject::AUTO_WRITE
  70         ),
  71         nu*rhoInfValue
  72     );
  73
  74     word kinematicCloudName("kinematicCloud");
  75     args.optionReadIfPresent("cloudName", kinematicCloudName);
  76
  77     Info<< "Constructing kinematicCloud " << kinematicCloudName << endl;
  78     basicKinematicCollidingCloud kinematicCloud
  79     (
  80         kinematicCloudName,
  81         rhoInf,
  82         U,
  83         mu,
  84         g
  85     );
  86
  87     IOobject Hheader
  88     (
  89         "H",
  90         runTime.timeName(),
  91         mesh,
  92         IOobject::MUST_READ,
  93         IOobject::AUTO_WRITE
  94     );
  95
  96     autoPtr<volVectorField> HPtr;
  97
  98     if (Hheader.headerOk())
  99     {
 100         Info<< "\nReading field H\n" << endl;
 101
 102         HPtr.reset(new volVectorField (Hheader, mesh));
 103     }
 104
 105     IOobject HdotGradHheader
 106     (
 107         "HdotGradH",
 108         runTime.timeName(),
 109         mesh,
 110         IOobject::MUST_READ,
 111         IOobject::AUTO_WRITE
 112     );
 113
 114     autoPtr<volVectorField> HdotGradHPtr;
 115
 116     if (HdotGradHheader.headerOk())
 117     {
 118         Info<< "Reading field HdotGradH" << endl;
 119
 120         HdotGradHPtr.reset(new volVectorField(HdotGradHheader, mesh));
 121     }
 122
 123     #include "createNonInertialFrameFields.H"