[su_junwei]

Dear Group

The bubble column test case in tutorials was carried out with twoPhaseEulerFoam solver.

With turbulence off( set in transportProperties file) , all went well. But if I turn turbulence on, the following messages appeared, and the execution aborted


DILUPBiCG: Solving for alpha, Initial residual = 0.560997, Final residual = 3.05437e-11, No Iterations 23
Dispersed phase volume fraction = 0.387282 Min(alpha) = -12.3385 Max(alpha) = 27.9207
DILUPBiCG: Solving for alpha, Initial residual = 0.985119, Final residual = 1.67429e-11, No Iterations 25
Dispersed phase volume fraction = 0.382111 Min(alpha) = -22.733 Max(alpha) = 29.9401
#0 Foam::error::printStack(Foam:stream&) in "/home/sujunwei/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/sujunwei/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xb7ef9420]
#3 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/home/sujunwei/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<foam::field,> const&, Foam::UPtrList<foam::lduinterfacefield> const&, Foam::Field<double> const&, unsigned char) const in "/home/sujunwei/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#5 Foam::GAMGSolver::Vcycle(Foam::PtrList<foam::lduma trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<foam::field<double> >&, Foam::PtrList<foam::field<double> >&, unsigned char) const in "/home/sujunwei/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/sujunwei/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#7 Foam::fvMatrix<double>::solve(Foam::Istream&) in "/home/sujunwei/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libfiniteVolume.so"
#8 main in "/home/sujunwei/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/twoPhaseE ulerFoam"
#9 __libc_start_main in "/lib/libc.so.6"
#10 Foam::regIOobject::readIfModified() in "/home/sujunwei/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/twoPhaseE ulerFoam"
float point exception



I also found that, the max phase fraction alpha>>1 see the message above.

I have read the thesis of Henrik Rusche.With the modified form of phase fraction equation, the problem is avoid. why it occured in my simulaiton.
Some things were omitted?

Any reply is appreciated!

Thanks, su junwei