# Some problems with twoPhaseEulerFoam

 Register Blogs Members List Search Today's Posts Mark Forums Read

 May 1, 2008, 09:59 Dear Group The bubble colum #1 Senior Member     su junwei Join Date: Mar 2009 Location: Xi'an China Posts: 151 Rep Power: 20 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&, Foam::Field const&, Foam::Field const&, Foam::Field const&) const in "/home/sujunwei/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so" #4 Foam::GAMGSolver::scalingFactor(Foam::Field&, Foam::lduMatrix const&, Foam::Field&, Foam::FieldField const&, Foam::UPtrList const&, Foam::Field const&, unsigned char) const in "/home/sujunwei/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so" #5 Foam::GAMGSolver::Vcycle(Foam::PtrList const&, Foam::Field&, Foam::Field const&, Foam::Field&, Foam::Field&, Foam::Field&, Foam::PtrList >&, Foam::PtrList >&, unsigned char) const in "/home/sujunwei/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so" #6 Foam::GAMGSolver::solve(Foam::Field&, Foam::Field const&, unsigned char) const in "/home/sujunwei/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so" #7 Foam::fvMatrix::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

November 1, 2012, 03:55
#2
Senior Member

Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 844
Rep Power: 18
Quote:
 Originally Posted by 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
could anyone know how handle this question?~

 November 2, 2012, 01:12 #3 Senior Member   Alberto Passalacqua Join Date: Mar 2009 Location: Ames, Iowa, United States Posts: 1,912 Rep Power: 36 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). __________________ Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using.