Hi Researchers,

I'm doing a 3D simulation of two and three bubble rising using a modified interFoam solver and I need to bubbles centre position, velocity and surface area.

For a single bubble rising I used swak4Foam expressions for example for bubble centre position in Y as follows :

Code:

bubbleCentreY

{

type swakExpression;

valueType internalField;

verbose true;

variables (

"Vol= sum (alpha1 < 0.5 ? vol() : 0);"

"VolY= sum (alpha1 < 0.5 ? pos().y*vol() : 0);"

);

expression "VolY/Vol";

accumulations (

min

);

}

Is there a way to compute these for two or three bubble in this manner?

Thanks in advance,

Best Regards,

Arsalan.