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

unable to get parabolic velocity profile with pimplefoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 8, 2012, 01:47
Default unable to get parabolic velocity profile with pimplefoam
  #1
Member
 
HouKen
Join Date: Jul 2011
Posts: 67
Rep Power: 15
houkensjtu is on a distinguished road
Hi foamers!

In my last thread, I modified a pimplefoam tutorial case into a laminar, 2d channel flow.
The boundary conditions I applied could be summarized as:

inlet: fixedvalue velocity / zerogradient pressure
outlet: zerogradient velocity / fixed value pressure(=0 pa)
wall: fixedvalue velocity (0 0 0) / zerogradient pressure
frontandback: empty

From theory I knew for a fully developed channel flow, velocity profile on outlet should be parabolic with max velocity equals 1.5 times average velocity which should be the inlet velocity here, and this profile should have nothing to do with viscosity and density.

However, even I simulate very long time(3000 seconds), I can only get the following velocity profile:

444.jpg
555.jpg

Any ideas about this?
Also I want to know, how could I debug the pimplefoam using gdb, for example, view how the pressure field changing every time step. Is it necessary to modify the source code to add some print out code? (I have already compiled of with a debug flag)

thanks in advance
houkensjtu is offline   Reply With Quote

Old   October 8, 2012, 03:02
Default
  #2
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,267
Blog Entries: 1
Rep Power: 25
nimasam is on a distinguished road
1- whats the value of your inlet velocity? is it a laminar at all?
2-what the ration of length to width? as you know when you apply a uniform inlet,
there is an entrance length, after that flow would be developed
nimasam is offline   Reply With Quote

Old   October 8, 2012, 03:32
Default
  #3
Member
 
HouKen
Join Date: Jul 2011
Posts: 67
Rep Power: 15
houkensjtu is on a distinguished road
Quote:
Originally Posted by nimasam View Post
1- whats the value of your inlet velocity? is it a laminar at all?
2-what the ration of length to width? as you know when you apply a uniform inlet,
there is an entrance length, after that flow would be developed
Thanks for reply!
I also noticed may be the problem is entrance length.
If I raise the dynamic viscosity 100 times, to 1e-3, I could now view the "normal" parabolic velocity profile.
555.jpg
It seems that the lower viscosity, the longer entrance length will be.
However I read my textbook which says entrance length is usually about 116*D for laminar flow, where D is the channel diameter...seems no relation with viscosity...
houkensjtu is offline   Reply With Quote

Old   October 8, 2012, 04:10
Default
  #4
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 22
Bernhard is on a distinguished road
Quote:
Originally Posted by houkensjtu View Post
However I read my textbook which says entrance length is usually about 116*D for laminar flow, where D is the channel diameter...seems no relation with viscosity...
I suppose this length should be Reynolds depedent. What is your Reynolds number? Increasing viscosity by a factor 100, decreases Reynolds number by a factor 100. (Hence the first question about the laminarity)
Bernhard is offline   Reply With Quote

Old   October 8, 2012, 05:41
Default
  #5
Member
 
HouKen
Join Date: Jul 2011
Posts: 67
Rep Power: 15
houkensjtu is on a distinguished road
Quote:
Originally Posted by Bernhard View Post
I suppose this length should be Reynolds depedent. What is your Reynolds number? Increasing viscosity by a factor 100, decreases Reynolds number by a factor 100. (Hence the first question about the laminarity)
hmm...I will recheck the entrance length on text book.
Would you kindly comment on my second question, about how to debug and view data with gdb?
houkensjtu is offline   Reply With Quote

Reply

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
[UDF] Relative coordinates in logarithmic velocity profile cfdworker Fluent UDF and Scheme Programming 25 July 13, 2017 04:44
3D parabolic velocity profile for channel flow seetaratnam FLUENT 1 August 10, 2016 04:08
Half parabolic profile for the velocity inlet adele FLUENT 4 July 19, 2012 12:43
Logarithmic velocity profile cfdworker Fluent UDF and Scheme Programming 0 April 23, 2009 20:09
parabolic velocity profile? bssdyl FLUENT 4 March 22, 2006 12:32


All times are GMT -4. The time now is 17:37.