applications/utilities/preProcessing/engineSwirl/createFields.H

   1 Info<< "Reading combustion properties\n" << endl;
   2
   3 IOdictionary engineGeometry
   4 (
   5     IOobject
   6     (
   7         "engineGeometry",
   8         runTime.constant(),
   9         mesh,
  10         IOobject::MUST_READ,
  11         IOobject::NO_WRITE
  12     )
  13 );
  14
  15 vector swirlAxis
  16 (
  17     engineGeometry.lookup("swirlAxis")
  18 );
  19
  20 vector swirlCenter
  21 (
  22     engineGeometry.lookup("swirlCenter")
  23 );
  24
  25 dimensionedScalar swirlRPMRatio
  26 (
  27     engineGeometry.lookup("swirlRPMRatio")
  28 );
  29
  30 dimensionedScalar swirlProfile
  31 (
  32     engineGeometry.lookup("swirlProfile")
  33 );
  34
  35 dimensionedScalar bore
  36 (
  37     engineGeometry.lookup("bore")
  38 );
  39
  40 dimensionedScalar rpm
  41 (
  42     engineGeometry.lookup("rpm")
  43 );
  44
  45
  46 Info<< "Reading field U\n" << endl;
  47 volVectorField U
  48 (
  49     IOobject
  50     (
  51         "U",
  52         runTime.timeName(),
  53         mesh,
  54         IOobject::MUST_READ,
  55         IOobject::AUTO_WRITE
  56     ),
  57     mesh
  58 );
  59
  60 vector zT = swirlAxis;
  61 vector yT = vector(0, zT.z(), -zT.y());
  62 vector xT = vector(zT.y()*zT.y() + zT.z()*zT.z(), -zT.x()*zT.y(), -zT.x()*zT.z());
  63
  64 // if swirl is around (1, 0, 0) we have to find another transformation
  65 if (mag(yT) < SMALL)
  66 {
  67     yT = vector(zT.y(), -zT.x(), 0);
  68     xT = vector(-zT.x()*zT.z(), -zT.y()*zT.z(), zT.x()*zT.x() + zT.y()*zT.y());
  69 }
  70
  71 //swirlAxis doesn't have to be of unit length.
  72 xT /= mag(xT);
  73 yT /= mag(yT);
  74 zT /= mag(zT);