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

How to specify BCs to solve a Poiseuille flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 11, 2008, 15:01
Default Hello, I was trying to solv
  #1
kar
Senior Member
 
Kārlis Repsons
Join Date: Mar 2009
Location: Latvia
Posts: 111
Rep Power: 8
kar is on a distinguished road
Hello,

I was trying to solve a parabolic flow in 2D square-like domain.

* In blockMeshDict and boundary file I specified just "patch" for inlet and outlet,
* In /0/U I set zeroGradient on in/out, (0 0 0) on walls + empty
* In /0/p I set fixedValue on in/out, zeroGradient on walls + empty.

But.... solution obviously diverge. Aren't those BCs sufficent to describe a problem?
I think, I don't understand how to set BCs in OpenFOAM, for example - why should I use physicalType? Maybe you can advice something useful at least for that poiseuille flow?
kar is offline   Reply With Quote

Old   February 11, 2008, 15:55
Default Oh, I see a parabolic velocity
  #2
kar
Senior Member
 
Kārlis Repsons
Join Date: Mar 2009
Location: Latvia
Posts: 111
Rep Power: 8
kar is on a distinguished road
Oh, I see a parabolic velocity distribution now and it feels so good! The problem was - I forgot to take into account a fact, that for icoFoam pressure must be divided by density. So I tried to solve for 1000 times greater pressure!

Maybe I'll see somewhere else why physicalType should be specified...
kar is offline   Reply With Quote

Old   February 11, 2008, 17:14
Default I checked solution more carefu
  #3
kar
Senior Member
 
Kārlis Repsons
Join Date: Mar 2009
Location: Latvia
Posts: 111
Rep Power: 8
kar is on a distinguished road
I checked solution more carefully and .. trouble doesn't stop - I have almost exactly 2x greater U values comparing with analytical!

So in SI:
nu = 1e-4,
geometry: dx=0.01, dy=0.1, dz=0.001,
flow is in y direction,
delta p = 0.016 (divided by ro)

According to wikipedia, Umax = 0.01. But solution gives 0.02!! How come?
kar is offline   Reply With Quote

Old   February 11, 2008, 18:12
Default Hi Karlis The Poiseuille-fl
  #4
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,594
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Hi Karlis

The Poiseuille-flow I could find on wiki (please add link in the future), was for a circular pipe flow. Your setup is for a 2D case, i.e. Poiseuille flow between two flat plates. The ratio between u_max in these two cases is exactly a factor of 2.

/ 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   February 12, 2008, 04:07
Default Thank you, Niels! I was think
  #5
kar
Senior Member
 
Kārlis Repsons
Join Date: Mar 2009
Location: Latvia
Posts: 111
Rep Power: 8
kar is on a distinguished road
Thank you, Niels!
I was thinking in a quick way - take a plane of cylinder and calculate for it.. Actually different problem is solved in such way.
kar is offline   Reply With Quote

Reply

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
Incompressible laminar poiseuille flow samarth OpenFOAM Running, Solving & CFD 20 July 13, 2010 03:38
Basic Poiseuille flow simulation question Ashish Main CFD Forum 0 October 2, 2007 13:05
Poiseuille flow Kirikou Main CFD Forum 0 November 23, 2006 14:32
Simple Poiseuille flow Maggie CFX 0 May 27, 2004 10:58
Poiseuille flow problem Rosie FLUENT 1 December 6, 2002 17:52


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