mpi problem with twoPhaseEulerFoam
 

Hello,

I think there is bug in OpenFOAM when using twoPhaseEulerFoam and try to run it parallel.

I was trying to run "bed2" tutorial on 2 CPU's but it crashed. It works fine with single CPU. The description of the problem can be found here.

It looks like the problem is coming from alpha equation. Below there is a printout from the log file. First the alpha equation is solved and max alpha=0.8, after correction of alpha equation, the min and max values of alpha are: Min(alpha) = -0.684932 Max(alpha) = 1.10975, which generates further problems. This is not the case when running it on single CPU.

Courant Number mean: 0.0280557 max: 0.36132
Max Ur Courant Number = 0.461666
Calculating averages

Time = 0.9655

DILUPBiCG: Solving for alpha, Initial residual = 0.00161798, Final residual = 3.43713e-11, No Iterations 4
Dispersed phase volume fraction = 0.275303 Min(alpha) = -2.12575e-18 Max(alpha) = 0.802811
DILUPBiCG: Solving for alpha, Initial residual = 0.000113185, Final residual = 5.23969e-11, No Iterations 3
Dispersed phase volume fraction = 0.275303 Min(alpha) = -2.11835e-18 Max(alpha) = 0.802806
GAMG: Solving for p, Initial residual = 0.000399759, Final residual = 8.44465e-09, No Iterations 12
time step continuity errors : sum local = 3.62591e-09, global = -7.10158e-10, cumulative = 2.70152e-07
DILUPBiCG: Solving for alpha, Initial residual = 8.22487e-05, Final residual = 6.99567e-11, No Iterations 3
Dispersed phase volume fraction = 0.275303 Min(alpha) = -2.11584e-18 Max(alpha) = 0.802824
DILUPBiCG: Solving for alpha, Initial residual = 8.09589e-06, Final residual = 0.00755801, No Iterations 1001
Dispersed phase volume fraction = 0.275245 Min(alpha) = -0.684932 Max(alpha) = 1.10975
GAMG: Solving for p, Initial residual = 0.00143021, Final residual = 7.14171e-09, No Iterations 19
time step continuity errors : sum local = 3.07254e-09, global = 5.55504e-10, cumulative = 2.70707e-07
DILUPBiCG: Solving for epsilon, Initial residual = 0.00673291, Final residual = 2.54617e-06, No Iterations 5
DILUPBiCG: Solving for k, Initial residual = 0.00751057, Final residual = 3.37159e-06, No Iterations 5


Regards,
Peter