CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM

buoyantBoussinesqSimpleFoam divergence problems

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 5, 2010, 04:39
Default buoyantBoussinesqSimpleFoam divergence problems
  #1
New Member
 
Andrea Beretti
Join Date: Nov 2009
Posts: 13
Rep Power: 7
beretti_a is on a distinguished road
Hello foamers,

I'm studying the mixing between two streams at different temperatures of an incompressible fluid.
For the analysis I've used the solver buoyantBoussinesqSimpleFoam.

The calculation has worked fine for about 2600 iterations and the solution is reliable. Also comparing the results obtained with Fluent I get a good agreement.
All of a sudden the solution has started to diverge shutting off the calculation. This is the printout I've obtained.

#0 Foam::error:rintStack(Foam::Ostream&) in "/home/andrea/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/andrea/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib/libc.so.6"
#3 Foam:ICPreconditioner::calcReciprocalD(Foam::Fie ld<double>&, Foam::lduMatrix const&) in "/home/andrea/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 Foam:ICPreconditioner:ICPreconditioner(Foam::l duMatrix::solver const&, Foam::dictionary const&) in "/home/andrea/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#5 Foam::lduMatrix:reconditioner::addsymMatrixConst ructorToTable<Foam:ICPreconditioner>::New(Foam:: lduMatrix::solver const&, Foam::dictionary const&) in "/home/andrea/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#6 Foam::lduMatrix:reconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) in "/home/andrea/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#7 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/andrea/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#8 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/andrea/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.so"
#9 main in "/home/andrea/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/buoyantBoussinesqSimpleFoam"
#10 __libc_start_main in "/lib/libc.so.6"
#11 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
Floating point exception


Does anybody have an ideas?

Thanks in advance,

Andrea
beretti_a is offline   Reply With Quote

Old   February 8, 2010, 09:25
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 7,518
Blog Entries: 33
Rep Power: 74
wyldckat is a jewel in the roughwyldckat is a jewel in the roughwyldckat is a jewel in the roughwyldckat is a jewel in the rough
Greetings Andrea,

It's good to know you've got OpenFOAM up and running

OK, by what I can estimate, you should check the output residuals in the final iterations. They will probably indicate issues with the propagation of really small or large values, leading to the busted simulation.

I remember I've come against such an issue some time ago, but it was because I hadn't checked the units of the geometry... I was using an STL in millimetres instead of meters... so the solver gave in to the massive small lunar object on the surface of the earth

Additionally, you could also have triggered a bug. Try following the instructions on the wiki here, to try to isolate the problem!

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   February 8, 2010, 12:13
Default
  #3
New Member
 
Andrea Beretti
Join Date: Nov 2009
Posts: 13
Rep Power: 7
beretti_a is on a distinguished road
Greatings Bruno,

It's good to hearing from you.
Yeah, few weeks ago I asked to my IT dept for a brand new machine and using Debian as OS I haven't had any further problems. Sincerely I still don't know what were the troubles with SLES 10. Probably I was missing some critical packages as you suggested.

Coming back to my simulation, I can exclude any mistake about the units. Please consider that the calculation reached convergence (checking both the residuals and the velocity in certain points of the domain) and gave me a physical solution comparable to the one I've gotten using Fluent.
As a sudden the solution starts to diverge, giving me higher residuals that kept the solution to be busted.

Anyway in these days I'm going to perform a debugging as you suggest. I'm not so expert and I guess I'll need a while to get through the problem. Nevertheless I'll keep you informed as soon as I found something.

Best regards,
Andrea
beretti_a is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Divergence problem for species transport model MY FLUENT 3 January 11, 2014 05:46
Divergence Problems Ben FLUENT 2 January 5, 2007 00:18
divergence problems Kaatje FLUENT 3 January 23, 2006 14:50
problems simulation ideal gas, divergence in AMG S Ralf Schmidt FLUENT 11 October 1, 2005 14:21
LES NITA giving divergence problems anindya FLUENT 0 June 20, 2005 08:23


All times are GMT -4. The time now is 09:48.