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

pimpleDyMFoam error message

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

Like Tree2Likes
  • 1 Post By louisgag
  • 1 Post By Artur

Reply
 
LinkBack Thread Tools Display Modes
Old   September 27, 2013, 02:07
Default pimpleDyMFoam error message
  #1
Member
 
laurentb's Avatar
 
Laurent B
Join Date: Jun 2009
Location: Lille, FRANCE
Posts: 70
Rep Power: 8
laurentb is on a distinguished road
Hi all,

I want to simulate a mixer (rushton turbine) and I m using pimpleDyMFoam solver but after some iterations i've got an error message. Can you help me understand why ?

Thanks
------------
Courant Number mean: 0.000809178 max: 0.483624
deltaT = 8.08409e-05
Time = 0.0321407

solidBodyMotionFunctions::rampedAxisRotationMotion ::transformation(): Time = 0.0321407 transformation: ((0 0 0) (0.99999 (0 0 0.00450739)))
AMI: Creating addressing and weights between 25130 source faces and 25130 target faces
AMI: Patch source weights min/max/average = 1.0565e-06, 1.00177, 0.834999
AMI: Patch target weights min/max/average = 0, 1.00183, 0.83511
PIMPLE: iteration 1
[2] #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[2] #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[2] #2 in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[2] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/pimpleDyMFoam"
[2] #5
[2] in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/pimpleDyMFoam"
[2] #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #7
[2] in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/pimpleDyMFoam"
[laurent-Precision-WorkStation-T7400:17917] *** Process received signal ***
[laurent-Precision-WorkStation-T7400:17917] Signal: Floating point exception (8)
[laurent-Precision-WorkStation-T7400:17917] Signal code: (-6)
[laurent-Precision-WorkStation-T7400:17917] Failing at address: 0x3e8000045fd
[laurent-Precision-WorkStation-T7400:17917] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7f92ee10b4a0]
[laurent-Precision-WorkStation-T7400:17917] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7f92ee10b425]
[laurent-Precision-WorkStation-T7400:17917] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7f92ee10b4a0]
[laurent-Precision-WorkStation-T7400:17917] [ 3] /opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKdRK NS_5UListIdEE+0x24) [0x7f92ef33a034]
[laurent-Precision-WorkStation-T7400:17917] [ 4] pimpleDyMFoam(_ZN4FoamdvINS_12fvPatchFieldENS_7vol MeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKNS_ 11dimensionedIdEERKS8_+0x2a0) [0x4480a0]
[laurent-Precision-WorkStation-T7400:17917] [ 5] pimpleDyMFoam() [0x41aa93]
[laurent-Precision-WorkStation-T7400:17917] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7f92ee0f676d]
[laurent-Precision-WorkStation-T7400:17917] [ 7] pimpleDyMFoam() [0x41eebd]
[laurent-Precision-WorkStation-T7400:17917] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 2 with PID 17917 on node laurent-Precision-WorkStation-T7400 exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------
laurentb is offline   Reply With Quote

Old   October 18, 2013, 05:48
Default
  #2
Senior Member
 
Artur's Avatar
 
Artur
Join Date: May 2013
Location: Southampton, UK
Posts: 226
Rep Power: 6
Artur is on a distinguished road
Take a look at the fragments highlighted in red:

Quote:
Originally Posted by laurentb View Post
Hi all,

AMI: Creating addressing and weights between 25130 source faces and 25130 target faces
AMI: Patch source weights min/max/average = 1.0565e-06, 1.00177, 0.834999
AMI: Patch target weights min/max/average = 0, 1.00183, 0.83511

mpirun noticed that process rank 2 with PID 17917 on node laurent-Precision-WorkStation-T7400 exited on signal 8 (Floating point exception).
Your sliding mesh explodes due to incorrect weighting (there is at least one face with zero weighting which is unphysical). There have been many posts on this forum giving guidance on how to get rid of that. In short, you should aim at having all your AMI weightings as close to 1.0 as possible.

Have a look at these to start with:

Problem using AMI

pimpleDyMfoam simulation keeps blowing up

Without seeing your case I'd recommend increasing the number of feature snap iterations and solver iterations in your snap controlls. Sometimes coarsening the mesh and increasing the AMI tolerance also helps.

Best wishes,

A
Artur is offline   Reply With Quote

Old   October 22, 2013, 08:12
Default
  #3
Member
 
laurentb's Avatar
 
Laurent B
Join Date: Jun 2009
Location: Lille, FRANCE
Posts: 70
Rep Power: 8
laurentb is on a distinguished road
If you want to have a look to my case (200 Mo) you can download it at :

https://www.dropbox.com/s/allf2zv7v1...L131003.tar.gz

I'm working on an hexaedral mesh in order to avoid zero wheighting. Does someone can explain exactly what are "patch source and target weights" ?
laurentb is offline   Reply With Quote

Old   May 13, 2015, 03:45
Default
  #4
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 208
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
For those reading on, an quick way to get around this weighting problem is sometimes to use a constant timestep (if you're using a variable one), and ensure it adds up correctly to the period of the mesh movement cycle and/or verify using the moveDynamicMesh command that you will never hit that null weighting zone...
-Louis
Artur likes this.
louisgag is offline   Reply With Quote

Old   May 13, 2015, 04:16
Default
  #5
Senior Member
 
Artur's Avatar
 
Artur
Join Date: May 2013
Location: Southampton, UK
Posts: 226
Rep Power: 6
Artur is on a distinguished road
Hi,

Thanks for sharing the tip!

Since I last posted in this thread I've also found that sometimes it's a very particular angular position that causes a crash, even if the overall AMI mesh is good. I've found that a hack that helps to deal with that sometimes is to save your solution JUST before the crash point (1 time step), let it crash, and then restart with a bigger time step (just temporarily) in order to "jump over" the problematic conformation of the interpolation patches. It's not a very neat way of doing it and it's quite tedious but helps if you need to get the results and don't have the time to re-mesh and re-run.

All the best,

A
louisgag likes this.
Artur is offline   Reply With Quote

Old   May 13, 2015, 04:19
Default
  #6
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 208
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Hi Artur,
I've noted the same thing, but how do you save just before the crash when you don't know beforehand where it will crash ?
Regards,
-Louis
louisgag is offline   Reply With Quote

Old   May 13, 2015, 04:31
Default
  #7
Senior Member
 
Artur's Avatar
 
Artur
Join Date: May 2013
Location: Southampton, UK
Posts: 226
Rep Power: 6
Artur is on a distinguished road
In general, for big simulations I tend to save my solution quite often anyway but use the purgeWrite parameter in controlDict to only keep the last 2 or 3 saves. Then, if a problem occurs, I restart and save every time step (it's terrible, I know, but as I said this is a hack). Then when it crashes I have the data for just before it did and hence can restart with a big time step to perform the "jump".

In the past it would usually happen to me more-or-less once or twice per revolution so it's bearable to deal with it this way. If it happened more often for a particular mesh I would probably just try to re-do it.

A
Artur is offline   Reply With Quote

Old   May 13, 2015, 05:48
Default
  #8
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 208
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Great, thanks for the tip about purgeWrite.
-Louis
louisgag 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
pimpleDyMFoam issue giovanidiniz OpenFOAM Running, Solving & CFD 1 July 5, 2013 07:25
Problem with transient simulation of pump inner flow by pimpleDyMFoam on 1.6-ext renyun0511 OpenFOAM Running, Solving & CFD 0 April 10, 2013 03:26
pimpleDyMFoam issue for pitching foil PLD OpenFOAM Running, Solving & CFD 2 November 23, 2012 03:59
pimpleDyMFoam samiam1000 OpenFOAM 2 September 19, 2012 10:11
Error with pimpleDyMFoam samiam1000 OpenFOAM 2 June 11, 2012 06:21


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