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

pressure driven flow in interFoam

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

Like Tree1Likes
  • 1 Post By deniggo

Reply
 
LinkBack Thread Tools Display Modes
Old   May 9, 2011, 13:24
Default pressure driven flow in interFoam
  #1
Member
 
Nico T
Join Date: Aug 2010
Location: Leipzig, Germany
Posts: 39
Rep Power: 7
deniggo is on a distinguished road
Hi,
I'm working on open channel flow using interFoam. A case with fixedValue-BC for U at the inlet and p_rgh fixedValue-BC at outlet is working properly (zeroGradient BC at outlet for U and inlet for p_rgh).

But is it possible to implement a pressure drop by specifying inlet AND outlet pressure fixedValue-BC (channel slope), without a fixedValue-BC for U? Thus, U should be developed corresponding to the pressure gradient.

I tried zeroGradient and pressureInletOutletVelocity for U-inlet and outlet, but the runs failed. Which BC for U would be appropriate for thise case?

Thanks for your help,

Nico
joris.hey likes this.
deniggo is offline   Reply With Quote

Old   May 9, 2011, 15:25
Default
  #2
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 418
Rep Power: 14
santiagomarquezd will become famous soon enough
What about totalPressure for p_rgh and zeroGradient for U and alpha1 at inlet?

Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Post-doctoral Fellow
Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL
T.E.: 54-342-4511594 Ext. 1005
Güemes 3450 - (3000) Santa Fe
Santa Fe - Argentina
http://www.cimec.org.ar
santiagomarquezd is offline   Reply With Quote

Old   May 11, 2011, 08:40
Default
  #3
Member
 
Nico T
Join Date: Aug 2010
Location: Leipzig, Germany
Posts: 39
Rep Power: 7
deniggo is on a distinguished road
Hi, thanks for your answer,
I tried also totalPressure for inlet and outlet of p_rgh with no success. totalPressure is also set for the top of the channel = atmosphere.
For alpha1, I forgot to mention, I use groovyBC at inlet for time varying waterlevel.
Maybe that's the problem?

Regards.
deniggo is offline   Reply With Quote

Old   June 24, 2011, 11:11
Default
  #4
New Member
 
Benjamin Mandt
Join Date: Sep 2009
Location: Germany
Posts: 8
Rep Power: 7
bmandt is on a distinguished road
Hi,

I have the same problem trying to calc a flooding inside a ship, only with pressure given at the inlet side.

In a personal discussion at another forum, I received this answer (the discussion was in german, I try to translate it):

"
Hydrostatic pressure boundary condition calculated as
pRefValue + rho*g.(x - pRefPoint)
where rho is provided and assumed uniform.
"
for the inlet pressure: "uniformDensityHydrostaticPressure"
the velocity at the inlet: "pressureInletUniformVelocity"
water at the inlet: alpha/gamma = fixedValue; value uniform 1;
at the outlet "zeroGradient"for using Piso


unfortunately the author mentioned piso, not interfoam, but it may be, and I hope this, this information helps you. I was not able to set up an interfoam case with this information, the reason may be, that I am a very beginner at openfoam and cfd...


I hope you will give feedback.


Best regards


Benjamin


bmandt is offline   Reply With Quote

Old   June 25, 2011, 08:13
Default
  #5
New Member
 
HuyHoang
Join Date: Mar 2011
Posts: 7
Rep Power: 6
bongbongxanh is on a distinguished road
i suggest convert the pressure drop at the boundary condition into a pressure source, i.e. a constant gradient of pressure in the desired flow direction and add it to momentum equation in the solver. the velocity at the boundary can be set at zeroGradient or inletOutlet, the velocity itselft will then be driven by the pressure source and develop over the length of the channel.
Hope this idea can help.
bongbongxanh is offline   Reply With Quote

Old   June 25, 2011, 15:42
Default
  #6
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,123
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
look into channelFoam, you can download channelInterFoam here, but i cant promise its totally true

http://www.4shared.com/file/75DpcP1g...InterFoam.html
nimasam is offline   Reply With Quote

Old   October 11, 2011, 06:45
Default
  #7
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 275
Rep Power: 7
Andrea_85 is on a distinguished road
Hi Nico,
Have you found a solution? Have you tried channelFoam? I need to specify a pressure gradient between the inlet and the outlet of my domain and i would like that the flow was driven by this gradient.

Best

andrea
Andrea_85 is offline   Reply With Quote

Old   October 11, 2011, 16:59
Default
  #8
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,123
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
use a fan boundary condition between outlet and inlet so the pressure gradient* channel length will be the value of pressure jump
nimasam is offline   Reply With Quote

Old   October 12, 2011, 03:52
Default
  #9
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 275
Rep Power: 7
Andrea_85 is on a distinguished road
Hi nima and thanks for reply.
I never used fan BC, can you be a bit more specific about how to use them for my case? or can you direct me to some reference or tutorial?

thanks again

andrea
Andrea_85 is offline   Reply With Quote

Old   October 13, 2011, 04:51
Default
  #10
Member
 
Nico T
Join Date: Aug 2010
Location: Leipzig, Germany
Posts: 39
Rep Power: 7
deniggo is on a distinguished road
Hi nimasam,
thanks for your help.
I tried to use fan BC. But I use groovyBC for defining the water level at the inlet. Therefore, cyclic fan BC does not work, because "inout" (cyclic patch of inlet & outlet) can not be defined by groovyBC?!

Nico
deniggo is offline   Reply With Quote

Old   October 13, 2011, 05:01
Default
  #11
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,123
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
use the funkySetFields to initiate the level of the water at time zero!
nimasam is offline   Reply With Quote

Old   October 13, 2011, 06:00
Default
  #12
Member
 
Nico T
Join Date: Aug 2010
Location: Leipzig, Germany
Posts: 39
Rep Power: 7
deniggo is on a distinguished road
ok! defining water level at time 0 by using setFields works. water level stays constant through the simulation.
thank you very much!
deniggo is offline   Reply With Quote

Old   October 18, 2011, 07:28
Default
  #13
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 275
Rep Power: 7
Andrea_85 is on a distinguished road
Hi again,
I understand how to define a "inout" cyclic patch between inlet and outlet. Now my question is which is the correct boundary condition for alpha1 in "inout"?
i have

0/U
inout
{
type cyclic;
value uniform (0 0 0);
}

0/p_rgh
inout
{
type fan;
patchType cyclic;
f List<scalar> 1(1000);
value uniform 0;
}

0/alpha1
inout ??

Can i use the inletOutlet patch? I would like to have fixedValue at the inlet (uniform 1) and something like zeroGradient or inletOutlet at the outlet.

best

andrea
Andrea_85 is offline   Reply With Quote

Old   October 18, 2011, 08:07
Default
  #14
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,123
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
no you can't , when you use the class cyclic, you can only use cyclic boundary condition or its sub derived class like fan
but the main question is that how ur alpha treats? is it cyclic or not? if it is not cyclic why you are going to use cyclic BC for ur simulation?
nimasam is offline   Reply With Quote

Old   October 18, 2011, 08:37
Default
  #15
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 275
Rep Power: 7
Andrea_85 is on a distinguished road
Hi nima,
honestly i'm trying different types of boundary condition to see what is the best choice to reproduce an experiment. In the experiment the flow is driven by a pressure difference between inlet and outlet and I'm loking for something similar in my simulation (at least i would like to not specify the pressure at outlet, because i'm interesting in what happens at the breakthrough and i think fixed pressure BC affects the results). Probably you are right and my alpha1 is not cyclic because i want only phase1 which enters at the inlet and not a mix of the two (which is what flows out at the outlet). the problem remains!
if you have other ideas, of course are welcome.

best

andrea
Andrea_85 is offline   Reply With Quote

Old   September 13, 2013, 04:57
Default
  #16
Senior Member
 
Mojtaba.a's Avatar
 
Mojtaba Amiraslanpour
Join Date: Jun 2011
Location: Zanjan, Iran
Posts: 233
Rep Power: 7
Mojtaba.a is on a distinguished road
Send a message via Yahoo to Mojtaba.a
This maybe an old thread, but did you found any solution for this?
__________________
Complex Heat & Flow Simulation Research Group
If you can't explain it simply, you don't understand it well enough. "Richard Feynman"
Mojtaba.a is offline   Reply With Quote

Old   September 13, 2013, 08:51
Default
  #17
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,123
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
Quote:
Originally Posted by Mojtaba.a View Post
This maybe an old thread, but did you found any solution for this?
solution for what?
__________________
Training Course on OpenFOAM at (http://www.isme.ir/)
My Weblog (http://openfoam.blogfa.com/)
nimasam is offline   Reply With Quote

Old   September 13, 2013, 12:36
Default
  #18
Senior Member
 
Mojtaba.a's Avatar
 
Mojtaba Amiraslanpour
Join Date: Jun 2011
Location: Zanjan, Iran
Posts: 233
Rep Power: 7
Mojtaba.a is on a distinguished road
Send a message via Yahoo to Mojtaba.a
Quote:
Originally Posted by nimasam View Post
solution for what?
To solve for pressure driven flow using interFoam.
I mean the boundary conditions.

I don't know why not typical boundaries converge for this case.

P value must be greater that a specific value (i.e 20000 Pa) to be able to drive the flow. But the one I want it to be is about 1000 Pa. With this lower value I can not converge the solution.
__________________
Complex Heat & Flow Simulation Research Group
If you can't explain it simply, you don't understand it well enough. "Richard Feynman"

Last edited by Mojtaba.a; September 13, 2013 at 13:58.
Mojtaba.a is offline   Reply With Quote

Old   September 13, 2013, 13:21
Default
  #19
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,123
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
i have no idea, you may want to setup a simple test case and show this, which may help other OpenFOAM user help you to find out how to solve this problem
__________________
Training Course on OpenFOAM at (http://www.isme.ir/)
My Weblog (http://openfoam.blogfa.com/)
nimasam 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
Compressible flow, no data at the outlet mireis FLUENT 1 July 28, 2010 05:22
Simulating a high pressure flow through a valve Kromagnsss FLUENT 8 July 2, 2010 05:20
pressure BC in buoynacy driven flow Sasha FLUENT 3 October 11, 2006 10:08
fluid flow fundas ram Main CFD Forum 5 June 17, 2000 21:31
Hydrostatic pressure in 2-phase flow modeling (CFX4.2) HB &DS CFX 0 January 9, 2000 14:19


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