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

multiphaseEulerFoam: exception related to nu

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 21, 2014, 22:28
Default multiphaseEulerFoam: exception related to nu
  #1
New Member
 
Join Date: Jun 2011
Posts: 22
Rep Power: 14
timbojones is on a distinguished road
Hello Foamers,

I'll post more details on the specifics of my error later today, but wanted to post the general question if someone might have a quick fix.

I'm simulating a drum with one inlet and two outlets. The simulation runs fine with all quantities solving with reasonable values. However, after about 1300 time steps, the code will stop owing to aa what appears to be a division by zero error involving nu. My LES model has been set as laminar. Anyone have a similar experience and knows how to fix my mistake?

Sorry for the sloppy post. I'm not near my computer, but wanted to pose the general question. I'll post the exact error stream in a short while.

Thanks,
Timbojones
timbojones is offline   Reply With Quote

Old   September 22, 2014, 07:16
Default Printstack
  #2
New Member
 
Join Date: Jun 2011
Posts: 22
Rep Power: 14
timbojones is on a distinguished road
As promised, here is the printstack associated with the floating-point exception I discussed in my previous post.

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2
at sigaction.c:0
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField>(Foam::FieldField< Foam::fvPatchField, double>&, Foam::FieldField<Foam::fvPatchField, double> const&, Foam::FieldField<Foam::fvPatchField, double> const&) in "/opt/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libmultiphaseSystem.so"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libmultiphaseSystem.so"
#6 Foam::multiphaseSystem::nu() const in "/opt/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libmultiphaseSystem.so"
#7 Foam::incompressible::LESModels::laminar::nuSgs() const in "/opt/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libincompressibleLESModels.so"
#8 Foam::incompressible::LESModel::nut() const in "/opt/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libincompressibleLESModels.so"
#9
in "/opt/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/multiphaseEulerFoam"
#10 __libc_start_main in "/lib64/libc.so.6"
#11
in "/opt/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/multiphaseEulerFoam"
[1]+ Floating point exceptionmultiphaseEulerFoam > log.mpef

I set the LES model to be laminar for this case so I'm a bit confused by the references in lines #7 and #8 which seem to refer to nut() (turbulent viscosity) and nuSgs (viscosity related to LES). However, it does seem to be nu() that starts the whole mess. I've set values for nu for both phases in transportProperties as well as the dummy entry at the bottom of transportProperties. I'm not sure what's going wrong here and would be very grateful for some help.

Thanks and Regards,
Timbojones
timbojones is offline   Reply With Quote

Old   September 22, 2014, 08:20
Default
  #3
Senior Member
 
Gerhard Holzinger
Join Date: Feb 2012
Location: Austria
Posts: 339
Rep Power: 28
GerhardHolzinger will become famous soon enoughGerhardHolzinger will become famous soon enough
The Lines 3 and 4 of the stacktrace tell you that division cause the floating point exception, most probably division by zero.

The method Foam::multiphaseSystem::nu() computes the kinematic viscosity nu. In the return statement, there is a division ( return tmu/rho(); ).

Did your solution diverge otherwise, which caused the mixture density to become zero?

What does the solver output say? The stacktrace alone is too little information to determine the cause of the solver crash.
GerhardHolzinger is offline   Reply With Quote

Reply

Tags
exception, multiphaseeulerfoam

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
Floating Point Exception Error nyox FLUENT 11 November 30, 2018 12:31
Adding Boussinesq Approximation to multiphaseEulerFoam? dschmidt OpenFOAM Programming & Development 1 January 16, 2017 16:40
MPI Error - simpleFoam - Floating Point Exception scott OpenFOAM Running, Solving & CFD 3 April 13, 2012 16:34
Cannot Open .sim (Floating Point Exception) trex930 STAR-CCM+ 1 July 30, 2010 06:51
Received exception: Disk full or disk related problem thekay ANSYS Meshing & Geometry 4 June 20, 2010 10:22


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