CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Bugs

OpenFOAM-2.1 instability problems?

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

LinkBack Thread Tools Display Modes
Old   May 25, 2012, 07:35
Default OpenFOAM-2.1 instability problems?
New Member
Ali Q Raeini
Join Date: Feb 2010
Posts: 21
Rep Power: 8
aliqasemi is on a distinguished road
I had some problems with using OpenFOAM-2.1, and I was wondering if anybody else had similar problems, or if they have worked on finding the solution.

I have noticed when I run my two phase simulations using OpenFOAM-2.1 they are more prone to instability, specially when run in parallel. When using the OpenFOAM-1.6-extend, however, my simulation are more stable,

The instabilities occur after some time from the beginning of the simulation and appear in terms of some non-physical velocities slowly increasing and finally destroying the simulation results. When I run my code in serial if I just use smaller time-steps then the simulations are less likely to diverge. I haven't seen these instabilities when using OF-1.6-ext. So my first question is what has made the OF-1.6-ext more stable. It should be some improvements in discretization algorithms, but it would be great if we know more exactly what are these improvements.

There is another problem which is more an issue for me. It happens when using OF-2.1 in parallel, where I have seen instabilities mainly occurring in the processor boundaries. This time the problem can be due to a looser coupling of the linear equation solvers (I use GAMG most of the times), or maybe improper handling of the boundary conditions. I don't have this problem when using OF-1.6-ext, with identical algorithm, so this is not a problem with my code. My simulations are two-phase flow at low capillary numbers, and these instabilities happen when capillary forces are active on the processor boundaries.

At this time, I use OF-2.1 for preprocessing and post processing, while OF-1.6-ext for running my simulations. So it is not that bad, but it would be better if I could figure out what is the problem.

Is there anybody else having similar problems?
aliqasemi is offline   Reply With Quote

Old   May 25, 2012, 07:55
Senior Member
Join Date: Jun 2009
Location: France, grenoble
Posts: 263
Rep Power: 10
olivierG is on a distinguished road

If i remember correctely, Co calculation has changed with the 2.0/2.1 version, so this may be a starting point to investigate : try with lower time step and compare with 1.6-ext .

olivierG is offline   Reply With Quote

Old   May 25, 2012, 16:42
Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,679
Blog Entries: 39
Rep Power: 103
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
Greetings to both of you!

I knew I had seen reports about this sometime ago in bug tracker and here they are:
Best regards,
wyldckat is offline   Reply With Quote

Old   May 26, 2012, 15:07
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,910
Rep Power: 27
alberto will become famous soon enoughalberto will become famous soon enough
What solver(s) is showing this problem? InterFoam?
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats.
OpenQBMM - An open-source implementation of quadrature-based moment methods.

To obtain more accurate answers, please specify the version of OpenFOAM you are using.
alberto is offline   Reply With Quote

Old   May 26, 2012, 16:38
New Member
Ali Q Raeini
Join Date: Feb 2010
Posts: 21
Rep Power: 8
aliqasemi is on a distinguished road
Dear Alberto, Bruno and Olivier

Thanks for your help. The modification on Courant number calculation can explain my problem, or part of it at least. I didn't know about it and were using the same maxCo for both OF1.6 and the OF2.1. Nevertheless, I have put an additional constraint on time step size myself, so I need to check to see whether fixing this solves the problem.

The code is something based on the interFoam code, but heavily modified.

I will test this out further and post back.

Thanks again,
aliqasemi is offline   Reply With Quote

Old   September 24, 2012, 13:23
Default update
New Member
Ali Q Raeini
Join Date: Feb 2010
Posts: 21
Rep Power: 8
aliqasemi is on a distinguished road
Just to confirm that this problem is NOT due to the time-step size. I have to mention that this happens for parallel cases; the difference is not that significant for the single-CPU simulations. I suspect it is due to the linear eqn solver (the computed pressure is different from the beginning of the simulation), or I may have done something with the boundary conditions which is not consistent with the new version. Anyway, I rather stick to the old OF for now; I might spend more time on this if I have got the time.

Last edited by aliqasemi; September 24, 2012 at 14:09.
aliqasemi is offline   Reply With Quote

Old   May 17, 2016, 10:46
New Member
Join Date: Feb 2016
Posts: 1
Rep Power: 0
mara61 is on a distinguished road
Dear all,

I am experiencing similar problem with my solver (laminar melting with buoyant forces), it works in serial and crashes at some point in parallel. I tried modifying linear solvers, but without much success. Have you solved the problem? Do you have any suggestions?

mara61 is offline   Reply With Quote


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
Can openfoam solve problems in cylindrical and spherical coordinate systems? GGerber OpenFOAM Running, Solving & CFD 6 February 5, 2015 11:37
New OpenFOAM Forum Structure jola OpenFOAM 2 October 19, 2011 06:55
64bitrhel5 OF installation instructions mirko OpenFOAM Installation 2 August 12, 2008 18:07
OpenFOAM installation on SUSE92 problems derath OpenFOAM Installation 1 April 3, 2006 06:44
Problems of numerical instability because of high gradients Xiangyang Ye Main CFD Forum 4 September 28, 1998 03:48

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