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

Pressure reference in cyclic interFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 9, 2011, 04:58
Default Pressure reference in cyclic interFoam
  #1
Senior Member
 
Robert Sawko
Join Date: Mar 2009
Posts: 107
Rep Power: 12
AlmostSurelyRob will become famous soon enough
Dear Foamers,

I am seeking some advice on interFoam, pressure reference cells and cyclic boundaries conditions. I am trying to modify interFoam to run cyclic boundary conditions in a similar manner to channelFoam and then apply this model to simulate stratified flow in a pipeline.

I set up pipeline geometry, I put cyclic patches, I write velocity and pressure corrections to account for specified flow rates for each of the two phases.

When I was preparing my case I removed all pressure boundaries so eventually the only patches I have are wall and cyclic. interFoam would protest when I run it so I had to add

Code:
pRefCell 0
pRefValue 0
to system/fvSolution. My simulations seem to converge but after several hundred iterations they begin to oscillate.

So I modified the code to output cell residuals for each quantity. The most exciting action happens at the interface between phases (expected!) and at pressure reference cell (surprising?). For some reason velocity residuals are always very high at pressure reference cell. Please have a look at the attachment.

It would be great if you had the time to lend me some assistance on this matter. I can send/attach some case files if think it necessary.
Attached Images
File Type: jpeg pRefCell.jpeg (55.0 KB, 52 views)
AlmostSurelyRob is offline   Reply With Quote

Old   February 11, 2011, 09:42
Default >1 values of phase fraction!
  #2
Senior Member
 
Robert Sawko
Join Date: Mar 2009
Posts: 107
Rep Power: 12
AlmostSurelyRob will become famous soon enough
Something curious, which seems to be related to this problem, are increased values of phase fraction. In the region around the pressure reference cell I get values of alpha1 greater than 1. What is the significance of it?
AlmostSurelyRob is offline   Reply With Quote

Old   February 13, 2011, 14:48
Default
  #3
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,048
Blog Entries: 1
Rep Power: 13
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
hi friend
i had the same problem in cyclic boundary in interFoam, i dont know its for pressure refrence or not but it cant be solved with finner mesh or finner delta time step or alpha correction or ...... alpha residual is up rising !!!!!!!!
nimasam is offline   Reply With Quote

Old   February 14, 2011, 06:03
Default
  #4
Senior Member
 
Robert Sawko
Join Date: Mar 2009
Posts: 107
Rep Power: 12
AlmostSurelyRob will become famous soon enough
Thanks for your interest. I don't know what type of multiphase case you're doing. I am trying a stratified flow in a simple horizontal pipeline. It's quite possible that I have chosen inputs that do not give steady solutions so I am trying to find some experimental/theoretical cases for which I am sure that the interface is smooth.


I was wrong about the pressure cell. The curious behaviour of high alpha values occurs somewhere else. You can see my geometry, pressure reference cell and high alpha values in the attached picture. Also the Maximum value of alpha is increasing up to 1.14 and then decreasing again.
Attached Images
File Type: jpg screen.jpg (86.1 KB, 39 views)
Attached Files
File Type: pdf alpha.pdf (29.9 KB, 13 views)
AlmostSurelyRob is offline   Reply With Quote

Old   February 14, 2011, 06:56
Default interFoam - channelFoam - Periodic Flows
  #5
Member
 
Javier Basurco
Join Date: Jan 2010
Location: Rio de Janeiro, Brazil
Posts: 33
Rep Power: 6
Ingenierias2003 is on a distinguished road
Dear Robert,

I am working too in two-phase flow regarding periodic flow. I am trying to modify the interFoam solver in according to periodic flow (channelFoam).
I was not sure if you got to do it. I trying to do that my be could you help with any doubts.

My kind regards.

Respectfully
Javier Basurco
Ingenierias2003 is offline   Reply With Quote

Old   February 14, 2011, 08:10
Default
  #6
Senior Member
 
Robert Sawko
Join Date: Mar 2009
Posts: 107
Rep Power: 12
AlmostSurelyRob will become famous soon enough
Quote:
Originally Posted by Ingenierias2003 View Post
I am working too in two-phase flow regarding periodic flow. I am trying to modify the interFoam solver in according to periodic flow (channelFoam).
Well, in this case we are in the same boat or maybe pipeline, channel... or other conduit which is relevant to your application. ;-)

For now, I can recommend some literature and other posts on this forum.

Single phase:
1. J.Y.Murthy, S.Mathur "Periodic flow and heat transfer using unstructured meshes", Int. Journal for Numerical Method in Fluid, Vol 25,659-677(1997)

2. S.V.Patankar and D.B.Spalding "A calculation procedure for ehat, mass and momentum transfer in three dimensional parabolic flows" Int. J. of Heat Mass Transfer, Vol 15, pp. 1787-1806, 1972

There is also a splendid book which explains the basics of SIMPLE method:
"Numerical Heat Transfer and Fluid Flow" by S.V. Patankar, 1980

Two phase
Level set method:
3. G.Son "A numerical method for incompressible two-phase flow with open or periodic boundaries.", Numerical Heat Trasnfer Part P, 2001


There is also this post which you might find helpful
Cyclic boundaries importing Gambit msh file and interFoam

We could have a chat about it to exchange some ideas, but this will be a bit off topic. Tell me what you think.
AlmostSurelyRob is offline   Reply With Quote

Old   February 14, 2011, 08:52
Default
  #7
Member
 
Javier Basurco
Join Date: Jan 2010
Location: Rio de Janeiro, Brazil
Posts: 33
Rep Power: 6
Ingenierias2003 is on a distinguished road
Dear Robert,

Thanks for your quick reply. Definitively, I would like to exchange ideas with you. My e-mail is javierbas82@hotmail.com or javier.basurco@gmail.com. If you could add me to msn, I think is the best way to talk about that or maybe do you have a skype account.

I a look a forward to hearing with you

My deep regards

Respectfully
Javier Basurco
Ingenierias2003 is offline   Reply With Quote

Old   February 16, 2011, 09:58
Default
  #8
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,048
Blog Entries: 1
Rep Power: 13
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
i use a cyclic boundary condition, you can find my case here! i simulate a bubble train flow!
convergence in interFoam
nimasam is offline   Reply With Quote

Reply

Tags
cyclic boundaries, interfoam, pressure reference cell

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
specify a reference pressure point in FD method? Rui Main CFD Forum 10 August 30, 2010 08:50
Compiling new Solver with wmake lin123 OpenFOAM 3 April 13, 2010 14:18
Error with Wmake skabilan OpenFOAM Installation 3 July 28, 2009 00:35
Calculation of the Governing Equations Mihail CFX 5 July 25, 2008 17:29
Does star cd takes reference pressure? monica CD-adapco 1 April 19, 2007 11:26


All times are GMT -4. The time now is 01:09.