applications/solvers/surfaceTracking/surfactantFoam/createVolFields.H

   1     // Create volume-to surface mapping object
   2     volSurfaceMapping vsm(aMesh);
   3
   4     volScalarField Cvf
   5     (
   6         IOobject
   7         (
   8             "Cvf",
   9             runTime.timeName(),
  10             mesh,
  11             IOobject::NO_READ,
  12             IOobject::AUTO_WRITE
  13         ),
  14         mesh,
  15         dimensionedScalar("0", dimless/dimLength, 0)
  16     );
  17
  18     vsm.mapToVolume(Cs, Cvf.boundaryField());
  19     Cvf.write();
  20
  21     volVectorField U
  22     (
  23         IOobject
  24         (
  25             "U",
  26             runTime.timeName(),
  27             mesh,
  28             IOobject::NO_READ,
  29             IOobject::AUTO_WRITE
  30         ),
  31         mesh,
  32         dimensionedVector("zero", dimVelocity, vector::zero)
  33     );
  34
  35     vsm.mapToVolume(Us, U.boundaryField());
  36     U.write();