CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Some problems with twoPhaseEulerFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 1, 2008, 10:59
Default Dear Group The bubble colum
  #1
Senior Member
 
su_junwei's Avatar
 
su junwei
Join Date: Mar 2009
Location: Xi'an China
Posts: 151
Rep Power: 20
su_junwei is on a distinguished road
Send a message via MSN to 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
su_junwei is offline   Reply With Quote

Old   November 1, 2012, 04:55
Smile
  #2
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17
sharonyue is on a distinguished road
Quote:
Originally Posted by su_junwei View Post
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?~
sharonyue is offline   Reply With Quote

Old   November 2, 2012, 02:12
Default
  #3
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
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.
alberto is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
TwoPhaseEulerFoam convergence problems alberto OpenFOAM Running, Solving & CFD 74 June 7, 2016 06:26
TwoPhaseEulerFoam sara OpenFOAM Running, Solving & CFD 2 November 6, 2008 20:26
Bug in twoPhaseEulerFoam alberto OpenFOAM Bugs 2 May 20, 2008 22:25
TwoPhaseEulerFoam Bug alondono OpenFOAM Bugs 1 February 19, 2008 21:01
TwoPhaseEulerFoam newbee OpenFOAM 0 March 27, 2006 09:41


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