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

how to get parabolic velocity distribution in a channel?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 29, 2011, 20:42
Default how to get parabolic velocity distribution in a channel?
  #1
New Member
 
Roland
Join Date: Jul 2011
Posts: 5
Rep Power: 6
F1r3 is on a distinguished road
Hi,

I am new to openfoam and I want to simulate a channel flow with simpleFoam. The simulation runs but i have no parabolic velocity distribution when I look at my results. It seems that everywhere is the same velocity. I want to use the kOmegaSST model.

I'm searching in this forum but I didn't find a similar problem.
Could anybody help me please?

Here are my case folder and pictures of the velocity and pressure of my channel.

greetz Roland
Attached Images
File Type: jpg velocity.jpg (55.1 KB, 105 views)
File Type: jpg pressure.jpg (82.3 KB, 68 views)
Attached Files
File Type: gz Kanal.tar.gz (2.6 KB, 27 views)
F1r3 is offline   Reply With Quote

Old   July 30, 2011, 04:38
Default
  #2
Senior Member
 
Amir's Avatar
 
Amir
Join Date: May 2009
Location: Shiraz, Iran
Posts: 739
Blog Entries: 1
Rep Power: 14
Amir is on a distinguished road
Dear Roland,

you are using turbulent model and in turbulent flow the velocity profile is not parabolic. (it's not laminar)
Indeed, your velocity magnitude is very high, consequently, you achieved small boundary layer. Also, to reach such profiles, the flow has to be fully developed.

Amir

Last edited by Amir; July 30, 2011 at 07:26.
Amir is offline   Reply With Quote

Old   July 30, 2011, 08:10
Default
  #3
New Member
 
Roland
Join Date: Jul 2011
Posts: 5
Rep Power: 6
F1r3 is on a distinguished road
Hi Amir,

thanks for your reply.
So the kOmegaSST model is only for turbulent flow. When I change it to laminar I would get a parabolic velocity? Our could you tell me what to change in my case to reach a parabolic velocity? I've tried a simulation with a laminar model and set in the RASProperties turbulence to off but I don't reach a parabolic velocity. Furthermore I made my channel longer so that the flow can fully develop.
To your information I've never worked with CFD before...
Later on I want to simulate the flow through a hydraulic valve. There I think I have to simulate with the kOmegaSST model...

Roland
F1r3 is offline   Reply With Quote

Old   July 30, 2011, 08:56
Default
  #4
Senior Member
 
Amir's Avatar
 
Amir
Join Date: May 2009
Location: Shiraz, Iran
Posts: 739
Blog Entries: 1
Rep Power: 14
Amir is on a distinguished road
To reach parabolic profile, use laminar model and turn off other turbulent parameters; also, it doesn't need to extend geometry to achieve fully developed region, just set boundary condition of velocity to zeroGradient at both inlet and outlet of your pipe.


Amir
Amir is offline   Reply With Quote

Old   July 30, 2011, 09:42
Default
  #5
New Member
 
Roland
Join Date: Jul 2011
Posts: 5
Rep Power: 6
F1r3 is on a distinguished road
Thanks for your quick reply.

Now I've reached a parabolic velocity. I've changed my viscosity to a higher value (8.95e-04) only for testing. With the higher viscosity I see the parabolic profil. I think with a lower viscosity I have to simulate with a fine mesh cause the change in velocity is only at the wall and in a very thin region.

Now when I set the velocity for the beginning to (0 0 0) I get bad results. What do I have to change to get better results (faster convergence) cause I can't calculate the velocity for my valve later on...

A colleague told me that he reaches good results (good convergence) within 20 iterations with ANSYS CFX. Is it possible to reach that with OpenFOAM too or do have to run a lot more iterations?

Roland
Attached Images
File Type: jpg pressure.jpg (36.6 KB, 72 views)
File Type: jpg velocity.jpg (40.7 KB, 68 views)
F1r3 is offline   Reply With Quote

Old   July 30, 2011, 11:34
Default
  #6
Senior Member
 
Amir's Avatar
 
Amir
Join Date: May 2009
Location: Shiraz, Iran
Posts: 739
Blog Entries: 1
Rep Power: 14
Amir is on a distinguished road
Dear Roland,
Obviously, the results which you attached are incorrect and the most probable reason is that it's not converged yet.
Many factors can affect rate of convergency such selecting proper solver,schemes,good initialization, good solution control and so on. For your case, I think that it's more appropriate to initialize with average velocity.

PS: I have used both commercial software (FLUENT) and OpenFOAM; it can be proved that if the case is managed precisely, openFOAM would be much more efficient in rate of convergency.


Regards,


Amir
Amir is offline   Reply With Quote

Old   July 30, 2011, 12:52
Default
  #7
New Member
 
Roland
Join Date: Jul 2011
Posts: 5
Rep Power: 6
F1r3 is on a distinguished road
Hi Amir,

thanks a lot for your help. So I've to try with diffrent solvers and schemes for simulating my valve. Do you have any suggestions which solvers and schemes are good for hydraulic flow?

greetz Roland
F1r3 is offline   Reply With Quote

Old   July 30, 2011, 15:18
Default
  #8
Senior Member
 
Amir's Avatar
 
Amir
Join Date: May 2009
Location: Shiraz, Iran
Posts: 739
Blog Entries: 1
Rep Power: 14
Amir is on a distinguished road
Quote:
Originally Posted by F1r3 View Post
Hi Amir,

thanks a lot for your help. So I've to try with diffrent solvers and schemes for simulating my valve. Do you have any suggestions which solvers and schemes are good for hydraulic flow?

greetz Roland
Hi,

I believe that it depends on your specific case (flow and geometry); there are many solvers which you can use (refer to manual 3.5) and also different schemes, since this is not my field I prefer to let this discussion open for more experienced users.


Bests,


Amir
Amir is offline   Reply With Quote

Old   July 30, 2011, 19:26
Default
  #9
New Member
 
Roland
Join Date: Jul 2011
Posts: 5
Rep Power: 6
F1r3 is on a distinguished road
Dear Amir,

thanks for your help!
I will try with my case and if I don't find a solution I think it's best to start a new thread for a new topic.

greetz Roland
F1r3 is offline   Reply With Quote

Old   August 5, 2011, 06:19
Default
  #10
Senior Member
 
Mohammad
Join Date: Feb 2010
Location: Shiraz, Iran
Posts: 108
Rep Power: 7
m2montazari is on a distinguished road
hi Roland,
getting parabolic velocity profile is very simple! first you should specify what is the case you want to run? if you want to model a pipe with a fixed velocity at inlet and have a parabolic velocity at outlet, pipe length should be long enough to have fully develope flow at outlet. in this case you should specify a small velocity value (to have reynolds number less than 2300 to have laminar flow) at inelt and zerogradient at outlet; for pressure just use zerogradient at inlet and fixedvalue,value=0 at outlet.
but if you want to specify pressure gradient in pipe and have parabolic velocity profile, as amir said, just set zerogradient for velocity in both inlet and outlet, and fixedvalue for pressure(pressure of inlet should be a bit more than outlet).
you can use simplefoam solver for steady state and icofoam for unsteady(laminar) and piso foam for unsteady(turbulent and laminar). shemes should be linear or upwind and you should achieve converged solution in a minute!

I hope it helps,
Mohammad
m2montazari is offline   Reply With Quote

Reply

Tags
channel, parabolic velocity, pipe

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
Parabolic Velocity profile in a channel mohsen cheraghi OpenFOAM Running, Solving & CFD 3 April 13, 2011 04:21
Velocity profiles for channel and rectangular duct Jeff Main CFD Forum 0 November 21, 2005 20:46
problem in 3d parabolic velocity profile Lokesh FLUENT 8 August 11, 2005 05:36
PARABOLIC DISTRIBUTION OF VELOCITY pradeep sivakumar FLUENT 4 June 25, 2001 07:38
Exponential or parabolic velocity in inlet Kim, Hongju CFX 1 December 18, 2000 18:12


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