CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

simpleFoam - Convergence problem - Simple rectangular prism domain

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

Like Tree1Likes
  • 1 Post By fredo490

LinkBack Thread Tools Display Modes
Old   April 26, 2013, 17:10
Default simpleFoam - Convergence problem - Simple rectangular prism domain
Join Date: Aug 2012
Posts: 74
Rep Power: 5
HakikiCanakkaleli is on a distinguished road

There is a convergence issue that I suspect I had done something wrong with the boundary conditions.

=== 1 ===
Brief description of the domain:

Domain is a rectangular prism like a simple channel through which fluid flows. All lateral surfaces are set to slip boundary condition. Fluid flows from inlet to outlet.

Brief description of the numerical setting:

simpleFoam is used with k-omega SST turbulence model.


=== 2 ===

Mesh is generated in GAMBIT. Exported with ".msh" extention. Converted with "fluentMeshToFoam".


=== 3 ===


In my opinion, the setting is OK.

=== 4 ===

Fluid flows at a constant speed of 1.73 m/s from inlet to outlet.
The "0" directory files are as follows:


=== 5 ===

Residuals up to 400-time-step
2 representative time-iteration detail from log file

The velocity field is constant through entire volume of the domain. However, the pressure field varies towards outlet. It means to me, considering dynamic pressure is constant due to constant speed, the static pressure changes, IMHO.

Plus, in my guess, there is something wrong with the streamwise velocity predictions as the number of iterations shoots up to 1000 at each time step:

smoothSolver:  Solving for Uz, Initial residual = 0.000153982, Final residual = 7.208e-05, No Iterations 1000
=== 6 ===
I have found some other forum pages which consider the same error message in a slightly different context. Therefore, I somehow couldn't adapt the given answers to my case.

I appreciate any help.

Many thanks in advance.
HakikiCanakkaleli is offline   Reply With Quote

Old   April 29, 2013, 04:00
Senior Member
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,147
Blog Entries: 1
Rep Power: 15
nimasam is on a distinguished road
which version of OpenFOAM do you use?
it seems there is a bug for simpleFoam at 2.0.1, so i suggest to update your version, if you had, then it would be solved!
Training Course on OpenFOAM at (
My Weblog (
nimasam is offline   Reply With Quote

Old   May 7, 2013, 13:08
Join Date: Aug 2012
Posts: 74
Rep Power: 5
HakikiCanakkaleli is on a distinguished road
Many thanks for your reply.

I have tried "fluent3DMeshToFoam" and it didn't work as well.

The version I am using is 2.1.0.

I suspect with the mesh itself. So that if anyone wants to have to look, I uploaded it to the post.

Rectangular prism GAMBIT mesh
HakikiCanakkaleli is offline   Reply With Quote

Old   May 8, 2013, 14:14
Join Date: Aug 2012
Posts: 74
Rep Power: 5
HakikiCanakkaleli is on a distinguished road

== 1 ==
I found the reason why I have been obtaining such a high level of residuals.

I have been using the script from "how to plot residuals". It seems the default script has been plotting the initial residuals rather than the final residuals. (In the script, -f13 needs to be used instead of -f9 to plot the final residuals.).

== 2 ==

This muddling, however, contributed to an interesting observation to me.

I have made 2 identical rectangular prism grids in GAMBIT and with blockMesh.

The grid has (height x width x length) = ( 10 x 10 x 100) and meshed with a node per 1m.

Case is simply steady-state and laminar flow where Re = 10 (as the characteristic length is 100).

This is the residual plot from GAMBIT grid case.

This is the residual plot from the blockMesh grid case.

I found it interesting and wonder why? Would anyone could comment on it? Otherwise, let's note it as an experience and leave it.

If anyone would glance at so that I uploaded both computed cases: GAMBIT&blockMeshGrid_Laminar_simpleFoam.

Many thanks.
HakikiCanakkaleli is offline   Reply With Quote

Old   May 9, 2013, 01:06
New Member
Join Date: May 2013
Posts: 9
Rep Power: 4
-Ali is on a distinguished road
Dear Hakiki

Did you checked a positive initial value for the pressure?
I mean although you are running an incompressible flow, setting p=0 may cause some problem.
I am working on partly similar problem but compressible flow using rhoCenteralFoam. I have wrote some more details here.
Do you have any idea for me?
-Ali is offline   Reply With Quote

Old   May 9, 2013, 04:59
Senior Member
Join Date: Jul 2010
Posts: 236
Rep Power: 8
fredo490 is on a distinguished road
Dear Ali, you cannot compare a compressible and an incompressible solver. The compressible solvers use the absolute pressure value while the incompressible solvers use the gauge value (or relative if you like).

Therefor, a compressible solver will crash if you set p=0 but not an incompressible solver as p=0 is actually equal to p=101325 (or whatever other pressure because the solver is independent of the ambient pressure).
-Ali likes this.
fredo490 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
CM+5 Convergence of the nonstationary problem ILYA87 STAR-CCM+ 0 May 22, 2011 04:35
Heat Transfer simulation: No convergence problem fiqs CFX 2 April 21, 2010 15:47
convergence problem limseokmin FLUENT 3 November 14, 2004 13:43
extremely simple problem... can you solve it properly? Mikhail Main CFD Forum 40 September 9, 1999 09:11
convergence problem with SIMPLER NURAY KAYAKOL Main CFD Forum 1 February 24, 1999 14:43

All times are GMT -4. The time now is 00:41.