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: 12 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, 04:55
#2
Senior Member

Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 754
Rep Power: 10
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, 02:12 #3 Senior Member   Alberto Passalacqua Join Date: Mar 2009 Location: Ames, Iowa, United States Posts: 1,910 Rep Power: 28 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. OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using.

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post alberto OpenFOAM Running, Solving & CFD 74 June 7, 2016 05:26 sara OpenFOAM Running, Solving & CFD 2 November 6, 2008 20:26 alberto OpenFOAM Bugs 2 May 20, 2008 21:25 alondono OpenFOAM Bugs 1 February 19, 2008 21:01 newbee OpenFOAM 0 March 27, 2006 08:41

All times are GMT -4. The time now is 19:10.