applications/solvers/surfaceTracking/bubbleInterTrackFoam/createFields.H

   1
   2 Info<< "\nReading field p\n" << endl << flush;
   3 volScalarField p
   4 (
   5     IOobject
   6     (
   7         "p",
   8         runTime.timeName(),
   9         mesh,
  10         IOobject::MUST_READ,
  11         IOobject::AUTO_WRITE
  12     ),
  13     mesh
  14 );
  15
  16
  17 Info<< "Reading field U\n" << endl << flush;
  18 volVectorField U
  19 (
  20     IOobject
  21     (
  22         "U",
  23         runTime.timeName(),
  24         mesh,
  25         IOobject::MUST_READ,
  26         IOobject::AUTO_WRITE
  27     ),
  28     mesh
  29 );
  30
  31
  32 #   include "createPhi.H"
  33
  34
  35 volScalarField rho
  36 (
  37     IOobject
  38     (
  39         "rho",
  40         runTime.timeName(),
  41         mesh,
  42         IOobject::NO_READ,
  43         IOobject::NO_WRITE
  44     ),
  45     mesh,
  46     dimensionedScalar("0", dimMass/dimVolume, 0)
  47 );
  48
  49
  50 freeSurface interface(mesh, rho, U, p, phi);
  51
  52
  53 volScalarField* fluidIndicatorPtr = NULL;
  54
  55 if(Pstream::parRun())
  56 {
  57     fluidIndicatorPtr = new volScalarField
  58     (
  59         IOobject
  60         (
  61             "fluidIndicator",
  62             runTime.timeName(),
  63             mesh,
  64             IOobject::MUST_READ,
  65             IOobject::AUTO_WRITE
  66         ),
  67         mesh
  68     );
  69 }
  70 else
  71 {
  72     fluidIndicatorPtr = new volScalarField
  73     (
  74         IOobject
  75         (
  76             "fluidIndicator",
  77             runTime.timeName(),
  78             mesh,
  79             IOobject::NO_READ,
  80             IOobject::NO_WRITE
  81         ),
  82         interface.fluidIndicator()
  83     );
  84 }
  85
  86 volScalarField& fluidIndicator = *fluidIndicatorPtr;
  87
  88
  89 rho = fluidIndicator
  90    *(
  91         interface.rhoFluidA()
  92       - interface.rhoFluidB()
  93     )
  94   + interface.rhoFluidB();
  95
  96
  97 volScalarField mu
  98 (
  99     IOobject
 100     (
 101         "mu",
 102         runTime.timeName(),
 103         mesh,
 104         IOobject::NO_READ,
 105         IOobject::NO_WRITE
 106     ),
 107     fluidIndicator
 108    *(
 109         interface.muFluidA()
 110       - interface.muFluidB()
 111     )
 112   + interface.muFluidB()
 113 );
 114
 115
 116 // Read interpolators if present
 117 if (interface.twoFluids())
 118 {
 119     interface.interpolatorAB();
 120 }
 121
 122
 123 #   include "setRefCell.H"
 124