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

Unknown error when changing turbulence model

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 13, 2013, 03:57
Default Unknown error when changing turbulence model
  #1
Senior Member
 
Bobi
Join Date: Oct 2012
Posts: 286
Rep Power: 5
babakflame is on a distinguished road
Dear buddies
I changed my turbulence model from kepsilon to realizableKE
and this error happens during running:

Code:
Starting time loop

Time = 1

DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 1.76289432e-08, No Iterations 6
DILUPBiCG:  Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for Uz, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.000958120731, No Iterations 40
time step continuity errors : sum local = 0.0906109149, global = -0.00473111621, cumulative = -0.00473111621
DILUPBiCG:  Solving for H, Initial residual = 0.999999999, Final residual = 1.40752135e-08, No Iterations 9
Updating look-up table extractions...
Updating mass fraction extractions...
DILUPBiCG:  Solving for csi, Initial residual = 1, Final residual = 3.12448826e-08, No Iterations 10
DILUPBiCG:  Solving for csiv2, Initial residual = 1, Final residual = 5.93143858e-09, No Iterations 4
#0  Foam::error::printStack(Foam::Ostream&) in "/home/babak/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/home/babak/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/babak/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4  void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/babak/OpenFOAM/babak-2.1.x/platforms/linux64GccDPOpt/bin/turbulentFlameletRhoSimplecFoam"
#5   at realizableKE.C:0
#6  Foam::compressible::RASModels::realizableKE::correct() in "/home/babak/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#7  
 in "/home/babak/OpenFOAM/babak-2.1.x/platforms/linux64GccDPOpt/bin/turbulentFlameletRhoSimplecFoam"
#8  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9  
 in "/home/babak/OpenFOAM/babak-2.1.x/platforms/linux64GccDPOpt/bin/turbulentFlameletRhoSimplecFoam"
Floating point exception
would somebody plz hint me what is this error? and how can be fixed?


Regards
Bobi
babakflame is offline   Reply With Quote

Old   October 13, 2013, 17:15
Default
  #2
Senior Member
 
Lieven
Join Date: Dec 2011
Location: Mol, Belgium
Posts: 295
Rep Power: 13
Lieven will become famous soon enough
Dear Bobi,

This looks like you get a division by 0 in the evaluation of the fields in the turbulence model which (always) throws a Floating Point Exception.

I guess you set the initial value for k and/or epsilon to 0.0. If this is indeed the case, change this to something else. That will most likely solve the issue.

Cheers,

L
Lieven is offline   Reply With Quote

Old   October 13, 2013, 23:41
Default
  #3
Senior Member
 
Bobi
Join Date: Oct 2012
Posts: 286
Rep Power: 5
babakflame is on a distinguished road
Greetings Lieven

My case works with k-epsilon model. I tried to work with variants solvers and found that GAMG solver and Gauss Seidel smoother make this error when trying to move to for instance kOmegaSST model. I employed PCG abd PBiCG solvers then the problem removed.

Thanks for your hint.

Regards
Bobi
babakflame 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
Fluent :- turbulence Model prince_pahariaa FLUENT 6 November 25, 2014 13:22
Centrifugal Pump and Turbulence Model Michiel CFX 12 January 25, 2010 04:20
turbulence model equation Andy Chen FLOW-3D 4 January 1, 2010 22:45
Discussion: Reason of Turbulence!! Wen Long Main CFD Forum 3 May 15, 2009 09:52
Fan heater model: what turbulence source to use? andy20 Main CFD Forum 0 March 2, 2008 13:46


All times are GMT -4. The time now is 16:22.