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

Floating point exception (core dumped), running a new solver

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 13, 2015, 20:24
Default Floating point exception (core dumped), running a new solver
  #1
New Member
 
anonymous
Join Date: Jan 2012
Location: Canada
Posts: 24
Rep Power: 6
Mahyar Javidi is on a distinguished road
Hi everybody
I am running a new solver that I developed based on interFoam by adding a new field (volume fraction) to it and using a new viscosity model. It workes up to a certain time then I get this error: Floating point exception
I tried to change the solvers in fvSolution and schems in fvSchemes but no thing really changed.
I just copy/paste the last part that I get from solving a case by this solver, which may help to get a clue.



MULES: Solving for alpha1
Phase-1 volume fraction = 0.516683 Min(alpha1) = -6.43313e-18 Max(alpha1) = 1.1454
MULES: Solving for alpha1
Phase-1 volume fraction = 0.516631 Min(alpha1) = -1.54198e-21 Max(alpha1) = 1.1454
DICPCG: Solving for p_rgh, Initial residual = 0.00254906, Final residual = 8.62831e-05, No Iterations 3
time step continuity errors : sum local = 8.62297e-05, global = -2.64316e-08, cumulative = -2.53353e-05
DICPCG: Solving for p_rgh, Initial residual = 0.00031195, Final residual = 1.17906e-05, No Iterations 28
time step continuity errors : sum local = 1.61631e-05, global = 1.28852e-06, cumulative = -2.40467e-05
DICPCG: Solving for p_rgh, Initial residual = 0.000107948, Final residual = 9.00722e-09, No Iterations 82
time step continuity errors : sum local = 1.41569e-08, global = 1.41857e-09, cumulative = -2.40453e-05
GAMG: Solving for Qhi, Initial residual = 0.00241765, Final residual = 8.19049e-05, No Iterations 1
ExecutionTime = 87.7 s ClockTime = 89 s

Courant Number mean: 0.0534621 max: 0.428816
Interface Courant Number mean: 0.00291851 max: 0.357529
deltaT = 0.000111111
Time = 0.00222222

#0 Foam::error:rintStack(Foam::Ostream&) in "/home/mahyar/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/home/mahyar/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 in "/lib/x86_64-linux-gnu/libm.so.6"
#4 Foam:ow(Foam::Field<double>&, Foam::UList<double> const&, double const&) in "/home/mahyar/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:ow<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::dimensioned<double> const&) in "/home/mahyar/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleTransportModels.so"
#6 Foam::viscosityModels:articleVolumFraction::calc Nu() const in "/home/mahyar/OpenFOAM/mahyar-2.1.x/platforms/linux64GccDPOpt/lib/libuserparticleVolumFraction.so"
#7 Foam::viscosityModels:articleVolumFraction::corr ect() in "/home/mahyar/OpenFOAM/mahyar-2.1.x/platforms/linux64GccDPOpt/lib/libuserparticleVolumFraction.so"
#8 Foam::twoPhaseMixture::calcNu() in "/home/mahyar/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleTransportModels.so"
#9
in "/home/mahyar/OpenFOAM/mahyar-2.1.x/platforms/linux64GccDPOpt/bin/interFoamParticle"
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11
in "/home/mahyar/OpenFOAM/mahyar-2.1.x/platforms/linux64GccDPOpt/bin/interFoamParticle"
Floating point exception (core dumped)


I would appreciate any help
Mahyar
Mahyar Javidi is offline   Reply With Quote

Old   January 14, 2015, 08:55
Default
  #2
Senior Member
 
anonymous
Join Date: Aug 2014
Posts: 198
Rep Power: 4
ssss is on a distinguished road
So as you can see in the log, there is an error in the calcNu() function, when you call pow(). Try searching in your code if there can be a divition by zero or if there is a sqrt which can be imaginary
ssss is offline   Reply With Quote

Old   January 14, 2015, 11:34
Default
  #3
New Member
 
anonymous
Join Date: Jan 2012
Location: Canada
Posts: 24
Rep Power: 6
Mahyar Javidi is on a distinguished road
Thanks for the reply. I will check the code and let you know what was the issue.

Regards
Mahyar
Mahyar Javidi is offline   Reply With Quote

Old   January 15, 2015, 22:19
Default
  #4
New Member
 
anonymous
Join Date: Jan 2012
Location: Canada
Posts: 24
Rep Power: 6
Mahyar Javidi is on a distinguished road
Hi

The problem was division by the zero. Solver is running perfectly now. Thanks for the help.

Mahyar
Mahyar Javidi is offline   Reply With Quote

Old   January 16, 2015, 03:52
Default
  #5
Senior Member
 
anonymous
Join Date: Aug 2014
Posts: 198
Rep Power: 4
ssss is on a distinguished road
You are welcome
ssss is offline   Reply With Quote

Old   August 4, 2015, 04:01
Default
  #6
New Member
 
Kittipos
Join Date: Feb 2015
Posts: 7
Rep Power: 3
Kittipos is on a distinguished road
Hi

I used rhoPorousSimpleFoam to running simple case but it has an error

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam221/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #2 Uninterpreted: #3 Foam::heRhoThermo<Foam::rhoThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() in "/opt/openfoam221/platforms/linuxGccDPOpt/lib/libfluidThermophysicalModels.so" #4 Foam::heRhoThermo<Foam::rhoThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() in "/opt/openfoam221/platforms/linuxGccDPOpt/lib/libfluidThermophysicalModels.so" #5 in "/opt/openfoam221/platforms/linuxGccDPOpt/bin/rhoPorousSimpleFoam" #6 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6" #7 in "/opt/openfoam221/platforms/linuxGccDPOpt/bin/rhoPorousSimpleFoam" Floating point exception (core dumped)

Could you explain me pls
Kittipos 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
Floating Point Exception Error nyox FLUENT 8 January 22, 2016 08:53
solving a conduction problem in FLUENT using UDF Avin2407 Fluent UDF and Scheme Programming 1 March 13, 2015 03:02
floating point exception Error AmirBaqa1987 OpenFOAM 1 October 10, 2013 01:29
interDyMFoam- Floating point exception Error marhamat OpenFOAM Running, Solving & CFD 2 December 26, 2012 14:40
Gmsh and samplesurface touf Open Source Meshers: Gmsh, Netgen, CGNS, ... 2 December 10, 2007 03:27


All times are GMT -4. The time now is 14:35.