zonePhaseVolumes.H

Go to the documentation of this file.

{
    const scalarField& V = mesh.V();
 
    forAll(mesh.cellZones(), czi)
    {
        const labelList& cellLabels = mesh.cellZones()[czi];
 
        forAllConstIter(PtrDictionary<phaseModel>, fluid.phases(), iter)
        {
            const volScalarField& alpha = iter();
            scalar phaseVolume = 0;
 
            forAll(cellLabels, cli)
            {
                label celli = cellLabels[cli];
                phaseVolume += alpha[celli]*V[celli];
            }
 
            reduce(phaseVolume, sumOp<scalar>());
 
            Info<< alpha.name()
                << " phase volume in zone " << mesh.cellZones()[czi].name()
                << " = " << phaseVolume*1e6 << " ml " << endl;
        }
    }
}