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

simpleFoam: simple 1-D channel flow, yet very strange convergence behavior

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

LinkBack Thread Tools Display Modes
Old   June 15, 2013, 08:44
Exclamation simpleFoam: simple 1-D channel flow, yet very strange convergence behavior
New Member
Join Date: Jun 2013
Posts: 19
Rep Power: 5
kishpishar is on a distinguished road
Dear FOAMers,

I am new to OpenFOAM - so to start with, tried simpleFoam on one of the simplest of the flow situations - namely laminar flow between two parallel plates. The domain is simple rectangular - 0.5 m long and 0.05 m wide. I started with a 10 x 5 grid (10 cells lengthwise and 5 cells breadthwise) and noticed a very weird convergence behavior. The issue is that the case converges for some inlet velocities, but diverges for some others - even for this coarse grid. More specifically:

Converged perfectly for inlet velocities (V_in) upto 0.05 m/s. Beyond that:

Diverged for V_in upto 0.1 m/s (I tried 0.06, 0.07, 0.08, 0.09)

Diverged for V_in = 0.1 m/s, but CONVERGED for V_in = 0.1085 m/s

Diverged for 0.15 m/s, but converged for 0.2 m/s, again Diverged for 0.25 m/s

Diverged for 0.5 m/s, but Converged for 0.49 m/s

Converged for 1.0 m/s, but Diverged for 0.99 m/s

I see similar pattern for a 30 x 8 grid, only the above numbers are different.

Boundary conditions:

Inlet: fixedValue on U, zeroGradient on p

Outlet: zeroGradient on U, fixedValue (uniform 0) on p

Can someone help me understand what is going on here? Many thanks indeed.

- Kumar

P.S: I've attached my blockMesh dict and fvSchemes/fvSolution.

Sorry: I forgot to mention that the fluid used was water with kinematic viscosity 1.0 x 10^-6
Attached Files
File Type: txt blockMeshDict.txt (1.6 KB, 5 views)
File Type: txt fvSchemes.txt (1.9 KB, 7 views)
File Type: txt fvSolution.txt (1.6 KB, 7 views)

Last edited by kishpishar; June 16, 2013 at 03:57. Reason: To add one piece of info
kishpishar is offline   Reply With Quote

Old   June 19, 2013, 04:16
Senior Member
Tushar Chourushi
Join Date: Jul 2009
Location: IIT-Indore, India
Posts: 319
Blog Entries: 1
Rep Power: 10
Tushar@cfd is on a distinguished road

You should have provided with the Boundary condition files. Anyways, in your case a proper initialization of the solver is not done.

By initialization I mean to say like:

For velocity = 0.1 m/s,
BC: Inlet: fixedValue on U, zeroGradient on p.
And for initialization: internalField uniform (0.1 0 0);

The problem which you have mentioned is something unrealistic in CFD unless some parameters are missing.
Tushar@cfd is offline   Reply With Quote

Old   June 20, 2013, 13:55
New Member
Join Date: Jun 2013
Posts: 19
Rep Power: 5
kishpishar is on a distinguished road
Hi Tushar,

Many thanks for your time and help. Sorry I didn't include the 0/U file - I used U(0 0 0) to initialize; my p and U attached for your reference.

You are absolutely right - the case converges without any glitch for the initial U field (0.1 0 0), in a small no. of iterations. Isn't it because the initial guess is so close to the converged solution? Shouldn't the code be able to seek the right solution even for a bad guess, say (0 0 0). Field initialization based on the inlet value may not be ideal always (what if the pipe has a 90 degree bend or worse, if the flow recirculates in the domain?)

I was particularly confused why the instabilities develop so quickly for certain velocities and not for others? for comparison, see the cases of 0.1 (blows up) and 0.1085 m/s (nicely converges)

Thanks and regards
Attached Files
File Type: txt pandU.txt (2.4 KB, 4 views)
kishpishar 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
Behavior of k and epsilon in a channel flow mn1363 Main CFD Forum 1 January 11, 2013 12:49
Yacht in Open Channel Flow andreimour FLUENT 1 October 14, 2010 23:54
Complex problem: Open Channel Flow lett FLUENT 17 April 9, 2009 10:56
Strange Solution for a simple pipe flow!! shekharc Main CFD Forum 4 May 9, 2005 09:21
compressible channel flow.. R.D.Prabhu Main CFD Forum 0 July 17, 1998 17:23

All times are GMT -4. The time now is 23:56.