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

Steady pipe flow mean velocity higher than inlet velocity

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 6, 2008, 03:59
Default Hi, I am trying to simulate
  #1
New Member
 
Anita K.
Join Date: Mar 2009
Posts: 25
Rep Power: 17
anita is on a distinguished road
Hi,

I am trying to simulate a turbulant pipe flow with simpleFoam k-epsilon (with settings from pitzDaily tutorial).

My pipe has a diameter of 43.1 mm. At the inlet is a flow of 10 m/s.

In 2D:
If I take the velocity profile after a pipe length of 20 * diameter, there is a flat "parabolic" profile with mean velocity of 9.7 and a max velocity of 11.1 m/s. This is a factor 0.88.
My book says the factor should be around 0.8.

In 3D:
There is a mean velocity of 10.8 m/s and a max velocity of 12.5 m/s.

My problem is that the mean velocity is higher than the velocity at the inlet.
Furthermore the profile is developed already near the inlet and doesn't change along the pipe. In 2D I can see how the flow develops slowly.

Is this possible ?

Anita
anita is offline   Reply With Quote

Old   March 6, 2008, 04:40
Default Hi Anita, Could you post so
  #2
Member
 
lillberg's Avatar
 
Eric Lillberg
Join Date: Mar 2009
Location: Stockholm
Posts: 80
Rep Power: 17
lillberg is on a distinguished road
Send a message via Skype™ to lillberg
Hi Anita,

Could you post some pictures describing your problem. Also, what boundary conditions are you using on the inlet and outlet planes for U,P,k and epsilon.

Regards

/Eric
lillberg is offline   Reply With Quote

Old   March 6, 2008, 06:29
Default Hi Eric, boundary conditions
  #3
New Member
 
Anita K.
Join Date: Mar 2009
Posts: 25
Rep Power: 17
anita is on a distinguished road
Hi Eric,
boundary conditions from 3D case

Inlet
U fixedValue, uniform(0 10 0)
k fixedValue, uniform 0.375
p zeroGradient
epsilon fixedValue, uniform 8.7549

Outlet
U zeroGradient
k zeroGradient
p fixedValue, uniform 0
epsilon zeroGradient

profile development in 2D:

Cut through pipe:

Comparing of 2D and 3D results:

anita is offline   Reply With Quote

Old   March 6, 2008, 08:34
Default Hi Anita, The 3D case, what
  #4
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,900
Rep Power: 37
ngj will become famous soon enoughngj will become famous soon enough
Hi Anita,

The 3D case, what is the cross section like? Any way, the difference in that maximum velocity is dependent on the flow geometry.
What on the other hand is worrying me is that your mean velocity is not the same as your input mean velocity. Could you please try to evaluate it along the axis and see if this deviation is varying along the axis or not? This might tell you if you have done something wrong when evaluating U_mean.

Best regards,
Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is offline   Reply With Quote

Old   March 6, 2008, 11:11
Default Hi Niels, thanks for your h
  #5
New Member
 
Anita K.
Join Date: Mar 2009
Posts: 25
Rep Power: 17
anita is on a distinguished road
Hi Niels,

thanks for your help. The cross section is a circle.

When I evaluated U_mean again, I found out I was my fault. I calculated U_mean over a line on the cross section (like in 2D) instead over the area of cross section.

Now I get a U_mean of 9.82 m/s.

My second fault was that I defined inlet velocity in the wrong direction. So the flow was coming in through the outlet and going out though inlet.

I will calculate again with the right inlet conditions and evaluate U_mean along axis.

Anita
anita is offline   Reply With Quote

Old   March 7, 2008, 08:58
Default Hi, after rerunning the 3D
  #6
New Member
 
Anita K.
Join Date: Mar 2009
Posts: 25
Rep Power: 17
anita is on a distinguished road
Hi,

after rerunning the 3D case U_mean is between 9.94 and 9.96 m/s along the axis. 10 m/s is the inlet velocity. I think now this simulation should be ok.

Regards
Anita
anita is offline   Reply With Quote

Old   April 12, 2010, 07:29
Default
  #7
Member
 
Javed
Join Date: Mar 2010
Location: Mumbai,India
Posts: 32
Rep Power: 16
Javed is on a distinguished road
Hi Anita...
R u using axisymmetric geometry?

Thanks
Javed
Javed is offline   Reply With Quote

Old   September 25, 2012, 06:35
Thumbs down
  #8
New Member
 
abhi
Join Date: Sep 2012
Posts: 11
Rep Power: 13
ashgun is on a distinguished road
can u tell me how u calculated mean velocity?
I calculate the velocity value by mass weighted average value over cross section and that value is very less. close to 0.7 times Umax
ashgun 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
velocity profile as an velocity inlet condition KEES Main CFD Forum 0 April 15, 2008 12:26
reversed flow at velocity inlet / mass flow inlet ib FLUENT 1 March 26, 2007 14:11
SimpleFoam pipe flow velocity blowup gzink OpenFOAM Running, Solving & CFD 0 June 9, 2006 17:52
velocity or mass flow inlet? Silvia FLUENT 1 August 12, 2005 14:55
Can reversed flow occur at velocity inlet Bcs. kanok FLUENT 1 September 30, 2003 12:54


All times are GMT -4. The time now is 04:27.