[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

[Chethan]

Hi,
You are right. You can access the zone phase volumes. In fact it is automatically displayed onto your log during run time. The only pre- requisite is that, you need to have defined cell zones in your mesh.

You can define it in your blockMeshDict when you define your block. You can verify the existence of your zones in the poly mesh directory after running blockMesh.

Another hint for tracking cell zone based values would be through the use of libsofunctionobjects in controlDict, again a cell zone must be defined while creating the mesh.


Hope this helps.

Regards,
Chethan