jochemvandenbosch 
December 22, 2011 05:12 
2 Attachment(s)
Hi Stephen, my courant number is going like crazy indeed, but when I reduce my timestep with a factor 1000! to 0.0001 I get the same message after 54 iterations.
But when I look at those 54 iterations in ParaView I see no pressure distribution (see attachment) and also no speeds... This makes me believe there is something wrong in the boundary conditions....
Quote:
Time = 0.0054
Courant Number mean: 0.0418895 max: 355.16
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 4.47596e05, Final residual = 7.60243e08, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 6.83522e07, Final residual = 8.58481e08, No Iterations 125
DILUPBiCG: Solving for Uz, Initial residual = 9.32447e06, Final residual = 1.20745e08, No Iterations 1
DILUPBiCG: Solving for Alpha, Initial residual = 0.917158, Final residual = 6.78599e06, No Iterations 1001
Solid phase fraction = 22.4908 Min(Alpha) = 117512 Max(Alpha) = 92780.4
DICPCG: Solving for p_rgh, Initial residual = 0.996912, Final residual = 0.334048, No Iterations 1001
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.29575e+53, global = 2.09099e+53, cumulative = 6.60704e+65
DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 1.16875e16, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.51501e+150, global = 2.51501e+150, cumulative = 2.51501e+150
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam201/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam201/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam201/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#4
in "/opt/openfoam201/platforms/linuxGccDPOpt/bin/settlingFoam"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator*<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam201/platforms/linuxGccDPOpt/bin/settlingFoam"
#6
in "/opt/openfoam201/platforms/linuxGccDPOpt/bin/settlingFoam"
#7 __libc_start_main in "/lib/i386linuxgnu/libc.so.6"
#8
in "/opt/openfoam201/platforms/linuxGccDPOpt/bin/settlingFoam"
Floating point exception

