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

simpleFoam: Non-uniform mesh near inlet & outlet causes incorrect velocity profile?

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

LinkBack Thread Tools Search this Thread Display Modes
Old   January 28, 2015, 05:17
Default simpleFoam: Non-uniform mesh near inlet & outlet causes incorrect velocity profile?
New Member
Saumitra Joshi
Join Date: Dec 2012
Posts: 14
Rep Power: 13
Zaphod'sSecondHead is on a distinguished road
Hi all!

I'm running a very simple case of steady state, incompressible, laminar pipe flow using simpleFoam. I created the mesh using snappyHexMesh. The following three types of meshes were generated:

1) Fully uniform hex-mesh (uniform spacing between cells in flow direction)
2) Non-uniform hex mesh without layers (finer cells near inlet and outlet, no layers)
3) Non-uniform hex mesh with layers (finer cells near inlet and outlet, layers present) - image below.

Geometry: Straight cylindrical pipe
Diameter: 1 cm, Length: 1 m

Kinematic viscosity: 1e-6 m2/s

Boundary Conditions:
-> Inlet: Uniform Fixed Value ux = 0, uy = 0.1 m/s, uz = 0
-> Outlet: Zero Gradient
-> Wall: Uniform fixed value
-> Inlet, Wall: Zero Gradient
-> Outlet: Fixed Value = 0

I'm using the GAMG solver for pressure (Abs tolerance: 1e-8, Rel tolerance: 1e-3, other settings default), PBiCG solver for velocity (Abs tolerance: 1e-8, Rel tolerance: 1e-3, other settings default). Residual for pressure and velocity is 1e-6.

The Issue

For case 1, I get very nice results. (Parabolic profile, Ucent = 2*Uavg, correct pressure drop).

However for case 2 and 3, although the pressure drop is fine, the velocity profile along the center-line looks something like this:

Firstly, the velocity reaches only about 0.19 m/s instead of 0.2. Also, there is a weird kink towards the end. I checked if this is because of interpolation issues, but that is not the case.

I'm baffled. Is this because of the non-uniform mesh near inlet and outlet? If yes, why does it occur?
Zaphod'sSecondHead is offline   Reply With Quote


non uniform mesh problem, simplefoam, velocity field

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
LES supersonic free jet martyn88 OpenFOAM 22 April 17, 2015 06:00
Divergent temperature in chtMultiRegion(Simple)Foam akrasemann OpenFOAM Running, Solving & CFD 13 March 24, 2014 02:54
Inlet & Outlet Velocity BC issue naiter OpenFOAM 3 December 19, 2012 07:14
Need help with boundary conditions: open to atmosphere Wolle OpenFOAM 2 April 11, 2011 07:32
ATTENTION! Reliability problems in CFX 5.7 Joseph CFX 14 April 20, 2010 15:45

All times are GMT -4. The time now is 20:36.