applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/createFields.H

   1     Info<< "Reading thermophysical properties\n" << endl;
   2
   3     Info<< "Reading field T\n" << endl;
   4     volScalarField T
   5     (
   6         IOobject
   7         (
   8             "T",
   9             runTime.timeName(),
  10             mesh,
  11             IOobject::MUST_READ,
  12             IOobject::AUTO_WRITE
  13         ),
  14         mesh
  15     );
  16
  17     Info<< "Reading field p\n" << endl;
  18     volScalarField p
  19     (
  20         IOobject
  21         (
  22             "p",
  23             runTime.timeName(),
  24             mesh,
  25             IOobject::MUST_READ,
  26             IOobject::AUTO_WRITE
  27         ),
  28         mesh
  29     );
  30
  31     Info<< "Reading field U\n" << endl;
  32     volVectorField U
  33     (
  34         IOobject
  35         (
  36             "U",
  37             runTime.timeName(),
  38             mesh,
  39             IOobject::MUST_READ,
  40             IOobject::AUTO_WRITE
  41         ),
  42         mesh
  43     );
  44
  45 #   include "createPhi.H"
  46
  47 #   include "readTransportProperties.H"
  48
  49     Info<< "Creating turbulence model\n" << endl;
  50     autoPtr<incompressible::RASModel> turbulence
  51     (
  52         incompressible::RASModel::New(U, phi, laminarTransport)
  53     );
  54
  55     Info<< "Calculating field beta*(g.h)\n" << endl;
  56     surfaceScalarField betaghf("betagh", beta*(g & mesh.Cf()));
  57
  58     label pRefCell = 0;
  59     scalar pRefValue = 0.0;
  60     setRefCell
  61     (
  62         p,
  63         mesh.solutionDict().subDict("SIMPLE"),
  64         pRefCell,
  65         pRefValue
  66     );
  67
  68
  69     // Kinematic density for buoyancy force
  70     volScalarField rhok
  71     (
  72         IOobject
  73         (
  74             "rhok",
  75             runTime.timeName(),
  76             mesh
  77         ),
  78         1.0 - beta*(T - TRef)
  79     );