[maybee]

hi, am I right that the following code of file "zonePhaseVolumes.H" calculates the phase volumes for each phase in each cell?

{
    2     const scalarField& V = mesh.V();
    3 
    4     forAll(mesh.cellZones(), czi)
    5     {
    6         const labelList& cellLabels = mesh.cellZones()[czi];
    7 
    8         forAllConstIter(PtrDictionary<phaseModel>, fluid.phases(), iter)
    9         {
   10             const volScalarField& alpha = iter();
   11             scalar phaseVolume = 0;
   12 
   13             forAll(cellLabels, cli)
   14             {
   15                 label celli = cellLabels[cli];
   16                 phaseVolume += alpha[celli]*V[celli];
   17             }
   18 
   19             reduce(phaseVolume, sumOp<scalar>());
   20 
   21             Info<< alpha.name()
   22                 << " phase volume in zone " << mesh.cellZones()[czi].name()
   23                 << " = " << phaseVolume*1e6 << " ml " << endl;
   24         }
   25     }
   26 }

with reduce(...) see:

http://foam.sourceforge.net/docs/cpp...7e34b70aa75f91

I think sometimes a brief description of the content and especially the results of included files in solvers would be great. With results I mean what is calculated and where is it stored in the end...

greetings
maybee