File: zonePhaseVolumes.H

package info (click to toggle)
openfoam 4.1%2Bdfsg1-1
  • links: PTS, VCS
  • area: main
  • in suites: stretch
  • size: 163,028 kB
  • ctags: 58,990
  • sloc: cpp: 830,760; sh: 10,227; ansic: 8,215; xml: 745; lex: 437; awk: 194; sed: 91; makefile: 77; python: 18
file content (26 lines) | stat: -rw-r--r-- 718 bytes parent folder | download | duplicates (3) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
 
{
    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;
        }
    }
}