[mnikku]

Hi,
you haven't set any residualControl to PIMPLE as I suggested earlier:

PIMPLE: no residual control data found. Calculations will employ 3 corrector loops

From my experience, this leads to high initial residuals and funky results. If you check your average particle volume fraction (alpha.particles volume fraction = x), you see that it increases constantly, ie. your system is generating mass!

Again, to reach a converged solution, add residualControl to PIMPLE, there are several examples available for example here

[foamiste]

Thank you for your fast reply .

Please find attached my fvsolution, I've modified it, is it ok or i should add something else?

[foamiste]

Yes it seems that my system is generating mass!!! so what can I do to solve this problem?

[mnikku]

You need to make your simulation converge. Try it with your new fvSolution settings and you should see a difference.

[foamiste]

Yes it converges now!! thank you very much

Please do you have any idea about how to calculate torque and forces in multiphase, because with simpleFoam I use

functions
(
forces
{
type forces;
functionObjectLibs ("libforces.so");
outputControl timeStep;
outputInterval 1;
patches (wall_impeller_1 wall_shaft-fluid_tank wall_shaft-fluid_mrf_1);
pname p;
Uname U;
rhoName rhoInf;
log true;
rhoInf 998;
CofR (0 0 0);
}
);

But when I use it with multiphase case, I get this error:

--> FOAM FATAL ERROR:
No valid model for viscous stress calculation

From function forces::devRhoReff()
in file forces/forces.C at line 258.

FOAM exiting