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

multiphaseInterFoam floating point exception

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 14, 2019, 07:24
Unhappy multiphaseInterFoam floating point exception
  #1
New Member
 
Jad
Join Date: Jun 2013
Posts: 14
Rep Power: 10
JTDN is on a distinguished road
Hello foamers,



So I'm trying to model a simple separating tank (oil and water) using multiphaseInterFoam (unmodded solver from OF 18.12). As i'm relatively new at OF I have roughly indexed my case on the damBreak4phase tutorial (notably for transportProperties, sigmas, etc)

As soon as I run, I get a floating point exception (after one or two iterations) with the following error:



#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in /lib/x86_64-linux-gnu/libc.so.6
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#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&) at ??:?
#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&) at ??:?
#6 Foam::multiphaseMixture::solve() at ??:?
#7 ? at ??:?
#8 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#9 ? at ??:?
Floating point exception (core dumped)



I understand there's a divide by zero somewhere and it comes from "multiphaseMixture" but I can't get a handle on it. The pre-processing part up to and including setFields seems to work fine.



Troubleshooting that I've tried:

1) Maybe the issue came from way i defined the two-phase inlet (50% water, 50% oil) so I tried a one phase inlet - same issue

2) Tried modding the alpha.water file so it takes into account surface tension when dealing with the walls - same issue

3) Tried reducing inlet speed as i noticed Co number (not interface Co, the "normal" Co) was rising quickly. This mod did lower Co but the case still crashed.



I think the geometry and volume mesh are OK since I ran them first on interFoam (with water only) just to check the basic BC, the Co numbers I can expect, etc.


Any help would be much appreciated. I've attached my log file as well as the alpha.water, alpha.oil and alpha.air and finally fvSolutions.

Thanks a lot!
Jad
Attached Files
File Type: txt alpha.oil.txt (1.3 KB, 1 views)
File Type: txt alpha.air.txt (1.6 KB, 0 views)
File Type: txt alpha.water.txt (1.3 KB, 2 views)
File Type: txt log.testRun4.txt (7.1 KB, 1 views)
File Type: txt transportProperties.txt (1.5 KB, 1 views)
JTDN is offline   Reply With Quote

Old   May 3, 2021, 13:01
Default
  #2
New Member
 
Mahmoud
Join Date: Nov 2020
Posts: 4
Rep Power: 2
Mahmoud Abbaszadeh is on a distinguished road
Quote:
Originally Posted by JTDN View Post
Hello foamers,



So I'm trying to model a simple separating tank (oil and water) using multiphaseInterFoam (unmodded solver from OF 18.12). As i'm relatively new at OF I have roughly indexed my case on the damBreak4phase tutorial (notably for transportProperties, sigmas, etc)

As soon as I run, I get a floating point exception (after one or two iterations) with the following error:



#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in /lib/x86_64-linux-gnu/libc.so.6
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#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&) at ??:?
#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&) at ??:?
#6 Foam::multiphaseMixture::solve() at ??:?
#7 ? at ??:?
#8 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#9 ? at ??:?
Floating point exception (core dumped)



I understand there's a divide by zero somewhere and it comes from "multiphaseMixture" but I can't get a handle on it. The pre-processing part up to and including setFields seems to work fine.



Troubleshooting that I've tried:

1) Maybe the issue came from way i defined the two-phase inlet (50% water, 50% oil) so I tried a one phase inlet - same issue

2) Tried modding the alpha.water file so it takes into account surface tension when dealing with the walls - same issue

3) Tried reducing inlet speed as i noticed Co number (not interface Co, the "normal" Co) was rising quickly. This mod did lower Co but the case still crashed.



I think the geometry and volume mesh are OK since I ran them first on interFoam (with water only) just to check the basic BC, the Co numbers I can expect, etc.


Any help would be much appreciated. I've attached my log file as well as the alpha.water, alpha.oil and alpha.air and finally fvSolutions.

Thanks a lot!
Jad

I have got the same issue. Have you figured out how to resolve the error?
Mahmoud Abbaszadeh is offline   Reply With Quote

Reply

Tags
floating point exception, multiphase flow, multiphaseinterfoam

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
icoFoam floating point exception (8) leizhao512 OpenFOAM Running, Solving & CFD 7 November 1, 2018 11:43
A floating point exception has occurred: floating point exception [Overflow]. starlight STAR-CCM+ 4 May 4, 2016 09:08
A floating point exception - SEM Model yansheng STAR-CCM+ 1 April 4, 2016 04:57
Floating point exception from twoPhaseEulerFoam openfoammaofnepo OpenFOAM Running, Solving & CFD 1 March 19, 2016 13:56
Floating point exception (core dumped) for GAMG solver yuhou1989 OpenFOAM Running, Solving & CFD 2 March 24, 2015 19:28


All times are GMT -4. The time now is 18:28.