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

Outlet boundary condition for wave flume with interFoam solver

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By Arnoldinho

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 17, 2010, 14:09
Unhappy Outlet boundary condition for wave flume with interFoam solver
  #1
Senior Member
 
Arne Stahlmann
Join Date: Nov 2009
Location: Hanover, Germany
Posts: 209
Rep Power: 18
Arnoldinho is on a distinguished road
Dear Foamers,

for the simulation of waves and currents in a flume, I'm now struggling with the problem of the outlet boundary condition. Ideally, waves and currents shall just leave the modelling domain without generating any kind of reflections which would otherwise interfere with the next incoming waves in the flume.

I already read a lot of posts in this forum concerning this topic, but its still not clear to me what the best solution could be for such an outlet boundary condition - and esp. how to implement it. Some posted their ideas (sometimes already months ago), e.g. in the direction of a 'sponge layer'. But I'm not sure if the ideas worked or not. At the moment I'm using numerical absorption by an extended mesh, which works quite well for waves in the flume, but not for currents.
So, if those who already fixed this issue for their simulations and made their own boundary conditions could give a short comment on this, I would be glad!

Another (maybe stupid) question: My understanding of the zeroGradient boundary condition was that it sets the regarded value at the boundary patch equal to the near-patch cell value. Is this right?
My problem is: For modelling constant flow in the flume (no waves) I set p, U and alpha as zeroGradient for the outlet.
For U and alpha at the inlet, groovyBC was used with constant and uniform values. What I got in the simulation is that the water level in the flume rose during runtime, so there was no or at least not enough 'outflow'. Could anyone explain this to me?
I could fix this problem by setting U and alpha for both inlet and outlet using groovyBC - but this only works for constant flow and not for other situations or waves. Which leads to the first topic of a non-reflecting outlet bc...

Hoping to see some of the topics a bit clearer soon,

Arne
ahparvin likes this.

Last edited by Arnoldinho; November 18, 2010 at 03:28.
Arnoldinho is offline   Reply With Quote

Old   November 24, 2010, 09:56
Default
  #2
Senior Member
 
Robert Sawko
Join Date: Mar 2009
Posts: 117
Rep Power: 22
AlmostSurelyRob will become famous soon enough
Hello Arnoldinho,

I am also struggling with similar issues in a two-phase flow, just a bit different geometry. Have you made any progress regarding the outlet boundaries? Any references perhaps?
AlmostSurelyRob is offline   Reply With Quote

Old   November 24, 2010, 10:24
Default
  #3
Senior Member
 
Arne Stahlmann
Join Date: Nov 2009
Location: Hanover, Germany
Posts: 209
Rep Power: 18
Arnoldinho is on a distinguished road
Hi Robert,

may I lead you to http://www.cfd-online.com/Forums/ope...tml#post284589?! Nevertheless, no improvement so far...
Arnoldinho is offline   Reply With Quote

Old   November 24, 2010, 10:39
Default
  #4
Senior Member
 
Robert Sawko
Join Date: Mar 2009
Posts: 117
Rep Power: 22
AlmostSurelyRob will become famous soon enough
Thanks for that! I was afraid that you will mention transmissive boundaries eventually which I believe are not coded in multiphase codes.

But please let also interest in a different post of mine:
http://www.cfd-online.com/Forums/mai...d-channel.html

I am doing 2D closed channel, so apart from increasing mesh size I was thinking of using two horizontal outlets. If set up properly they will become outlets for only one of the phases. But that workaround is not applicable in your case.
AlmostSurelyRob is offline   Reply With Quote

Old   November 25, 2010, 04:53
Default
  #5
Senior Member
 
Arne Stahlmann
Join Date: Nov 2009
Location: Hanover, Germany
Posts: 209
Rep Power: 18
Arnoldinho is on a distinguished road
Robert, regarding your 2D channel, is this what you mean ? -> http://www.cfd-online.com/Forums/ope...interfoam.html
Arnoldinho is offline   Reply With Quote

Old   October 16, 2011, 23:59
Default
  #6
Member
 
Albert Tong
Join Date: Dec 2010
Location: Perth, WA, Australia
Posts: 76
Blog Entries: 1
Rep Power: 16
tfuwa is on a distinguished road
Quote:
Originally Posted by Arnoldinho View Post
Dear Foamers,

for the simulation of waves and currents in a flume, I'm now struggling with the problem of the outlet boundary condition. Ideally, waves and currents shall just leave the modelling domain without generating any kind of reflections which would otherwise interfere with the next incoming waves in the flume.

I already read a lot of posts in this forum concerning this topic, but its still not clear to me what the best solution could be for such an outlet boundary condition - and esp. how to implement it. Some posted their ideas (sometimes already months ago), e.g. in the direction of a 'sponge layer'. But I'm not sure if the ideas worked or not. At the moment I'm using numerical absorption by an extended mesh, which works quite well for waves in the flume, but not for currents.
So, if those who already fixed this issue for their simulations and made their own boundary conditions could give a short comment on this, I would be glad!

Another (maybe stupid) question: My understanding of the zeroGradient boundary condition was that it sets the regarded value at the boundary patch equal to the near-patch cell value. Is this right?
My problem is: For modelling constant flow in the flume (no waves) I set p, U and alpha as zeroGradient for the outlet.
For U and alpha at the inlet, groovyBC was used with constant and uniform values. What I got in the simulation is that the water level in the flume rose during runtime, so there was no or at least not enough 'outflow'. Could anyone explain this to me?
I could fix this problem by setting U and alpha for both inlet and outlet using groovyBC - but this only works for constant flow and not for other situations or waves. Which leads to the first topic of a non-reflecting outlet bc...

Hoping to see some of the topics a bit clearer soon,

Arne
Hi Arne,

Have you solved the problem? If yes, can you please tell me how?

Cheers,
Albert
tfuwa is offline   Reply With Quote

Old   October 17, 2011, 03:48
Default
  #7
Senior Member
 
Arne Stahlmann
Join Date: Nov 2009
Location: Hanover, Germany
Posts: 209
Rep Power: 18
Arnoldinho is on a distinguished road
Hi Albert,

you could have a look at this thesis http://www.google.de/url?sa=t&source...VJ_amw&cad=rja, which has, with some modifications, also worked for most of my cases.

Arne
Arnoldinho is offline   Reply With Quote

Old   May 10, 2012, 06:35
Default no ourflow with interFoam
  #8
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 240
Rep Power: 17
vonboett is on a distinguished road
...this might be some error on my side, but I did dozends of simulations both with setting type zeroGradient; for the outlet boundary condition in alpha1 as well as using type inletOutlet with inletValue uniform 0 for the atmosphere boundary, and both worked fine using OF 1.7.1 . But now with 2.1.x. I face the same problem as reported in different threads that the phase alpha1 gets reflected at the outflow, no matter if dense or coarse grid, zeroGradient or inletOutlet....
vonboett is offline   Reply With Quote

Old   May 23, 2012, 07:25
Default
  #9
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 240
Rep Power: 17
vonboett is on a distinguished road
Ok turning old searching the cause of alpha1 being reflected at the outflow, I finally got it. Maybe this is a bug dependent on ubuntu version, but it is quite relevant. The difference between the two pictures below showing an outflow of a channel is only that I moved the grid from positive x quadrant to negative x quadrant. When the whole grid lies at a position that the x-coordinates are smaller than 0 the outflow works! This strange behavior vanishes when using zeroGradient for p_rgh instead of fixedValue 0.
Attached Images
File Type: jpg impactOnOutflowT0,15.jpg (13.2 KB, 459 views)
File Type: jpg impactOnOutflowT0,15_withXCoordsmallerZero.jpg (14.6 KB, 123 views)

Last edited by vonboett; June 14, 2012 at 10:09.
vonboett is offline   Reply With Quote

Old   July 10, 2018, 06:15
Default
  #10
Member
 
hua1015's Avatar
 
Hua
Join Date: May 2012
Posts: 31
Rep Power: 14
hua1015 is on a distinguished road
Hi Arne,

I cannot open this link including the thesis you recommended. Would you mind tell me the name of this thesis?

Cheers,
Hua
hua1015 is offline   Reply With Quote

Reply


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
Selective outlet boundary condition for InterFoam caw OpenFOAM Running, Solving & CFD 5 February 7, 2012 14:48
Outlet boundary condition =/ AdidaKK CFX 4 December 4, 2009 09:12
Outlet Boundary Condition ? Maria FLUENT 3 August 20, 2008 09:30
Outlet boundary condition for pd in InterFoam gopala OpenFOAM Running, Solving & CFD 0 March 19, 2008 10:26
massflow outlet boundary condition ?? Charles Renard Main CFD Forum 6 August 20, 2001 04:27


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