CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Some problems with twoPhaseEulerFoam (https://www.cfd-online.com/Forums/openfoam-solving/58911-some-problems-twophaseeulerfoam.html)

 su_junwei May 1, 2008 09:59

Dear Group The bubble colum

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:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) 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"
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?

Thanks, su junwei

 sharonyue November 1, 2012 04:55

Quote:
 Originally Posted by su_junwei (Post 188931) 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
could anyone know how handle this question?~

 alberto November 2, 2012 02:12

It is quite simple: the turbulence model in twoPhaseEulerFoam is very simple, and not ideal to manage systems with phase separation.

You will need to implement a better model (see literature).

 All times are GMT -4. The time now is 07:43.