[sepponen]

Hello,

I would like to calculate the mass of a phase as function of time using twoPhaseEulerFoam. Could this be done without modifying the solver, maybe with a function object? Is there a function object which integrates a variable over the domain?

Thank you for your attention,
Simo

[tomf]

Hi,

If the density of the phase is constant you can follow the volume fraction of the phase that you are interested in in the log of the solver. If you want to plot it you can use some bash to extract the values from the log. I typically use something like this:

grep fra twoPhaseEulerFoam.log | awk '{print $5}' > volumeFraction

The grep command gives you something like this for every occurrence of "fra":

alpha.particles volume fraction = 0.0005668175397  Min(alpha.particles) = 0  Max(alpha.particles) = 0.05349701286

The awk command gives the 5th column, which is the overal volume fraction of the particles phase. You can than plot the file volumeFraction.

awk also allows multiplication to have the entire mass using:

awk '{print rho*$5*domainvolume}'

where rho is your density and domainvolume is the entire volume of the domain.

If your density is not constant you maybe can use the volRegion functionObject to obtain a weighted volume average or a volume integral.

Use the following to find some more information:

find $FOAM_SRC -iname *volRegion*