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

About flow division

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 25, 2009, 09:33
Default Hi, everyone, I have a geom
  #1
Member
 
Vivien
Join Date: Mar 2009
Posts: 52
Rep Power: 8
sunnysun is on a distinguished road
Hi, everyone,

I have a geometry which is roughly a bifurcation with cylinder tube(flow goes from cylinder to bifurcation, 1 inlet, 2 outlet). How can i set the flow division through boundary condisions? say I would like to have 30% of flow go in to one brance and 70% goes in to another? The inlet flow is pulsatile and I use timeVaringMappedFixedValue which read velocity from data file.

Thanks a lot!!

Vivien
sunnysun is offline   Reply With Quote

Old   February 25, 2009, 11:11
Default Seems to be a simple question.
  #2
Member
 
Vivien
Join Date: Mar 2009
Posts: 52
Rep Power: 8
sunnysun is on a distinguished road
Seems to be a simple question...but can somebody give me some ideas?
sunnysun is offline   Reply With Quote

Old   February 25, 2009, 11:11
Default Seems to be a simple question.
  #3
Member
 
Vivien
Join Date: Mar 2009
Posts: 52
Rep Power: 8
sunnysun is on a distinguished road
Seems to be a simple question...but can somebody give me some ideas?

Vivien
sunnysun is offline   Reply With Quote

Old   February 25, 2009, 12:28
Default Hi Vivien, For I know, ther
  #4
Member
 
Kati Laakkonen
Join Date: Mar 2009
Location: Espoo, Finland
Posts: 36
Rep Power: 8
kati is on a distinguished road
Hi Vivien,

For I know, there is no outlet BC in OpenFOAM which allows such predefined division of flow to different outlets (some codes do have such BC).

What you could do, is to set a flowRate BC on one of the outlets and define a negative flow rate that is about 30% of the inflow (there is also a timeVarying version of flowRate BC). Then you could use zero pressure on the other outlet for the rest of the outflow.

Regards,
Kati
kati is offline   Reply With Quote

Old   February 26, 2009, 04:59
Default Hi,Kati, Thanks a lot for y
  #5
Member
 
Vivien
Join Date: Mar 2009
Posts: 52
Rep Power: 8
sunnysun is on a distinguished road
Hi,Kati,

Thanks a lot for your answer.
What is still unclear for me is why should I use negative flow rate at onbe outlet?(say outlet1). Also, what should be the pressure BC for outlet1? Why zero pressure instead of zero gradient pressure BC at outlet2?

Thank you very much!

Vivien
sunnysun is offline   Reply With Quote

Old   March 2, 2009, 13:34
Default Vivien, Negative flow rate
  #6
Member
 
Kati Laakkonen
Join Date: Mar 2009
Location: Espoo, Finland
Posts: 36
Rep Power: 8
kati is on a distinguished road
Vivien,

Negative flow rate because positive would mean inflow.

Normal outlet BC is setting pressure and using zeroGradient for convected fields (actually inletOutlet instead of zeroGradient, more stable).

You could also fix flow rate at both boundaries and then use zeroGradient for pressure, but then you are fixing all flow rates on all boundaries, and I wouldn't do that because it leaves no freedom to the flow solution. It should in theory be okay with incompressible or steady state cases, but in compressible transient it is an error.

You could test different BC combinations and report the results here. I'd be interested to hear whether you had any problems with "fix all flow rates" approach.

Kati
kati 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
Division by zero vitke OpenFOAM Running, Solving & CFD 5 September 1, 2008 05:35
Rather obscure division by zero in triangleFuncsintersectAxesBundle gschaider OpenFOAM Bugs 2 August 13, 2008 07:08
Division by zero in Xoodles hannes OpenFOAM Bugs 3 August 4, 2008 11:04
ScalarField division maka OpenFOAM Pre-Processing 2 August 27, 2007 05:10


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