[scleakey]

OpenFOAM uses implicit methods so I thought I would see how far I could push the Courant number before it crashes - isn't the point of implicit methods to get large Courant numbers?

I'm running a multiphase case with interIsoFoam on six processors. Keeping everything else the same, I tried different values for maxCo...

1.0
2.0
2.3
2.5 crashed at t = 0.02 after high residual for p_rgh
2.7 crashed at t = 0.02 after high residual for p_rgh
3.0 crashed at t = 0.02 after high residual for p_rgh
5.0
10.0
20.0 crashed at t = 1.31, not sure exactly why

I don't understand how it ran for maxCo = 10 (with very sensible results) but crashed for maxCo = 2.5? I realise there are lots of reasons that it might crash but I thought that decreasing the Courant number would make it more stable, not less

Has anybody else experienced this?