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

Outlet boundary condition for wave flume with interFoam solver

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

Reply
 
LinkBack Thread Tools 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: 8
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

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: 13
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: 8
Arnoldinho is on a distinguished road
Hi Robert,

may I lead you to transmissive BC / numerical beach?! 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: 13
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:
Pressure outlet in two-phase flow in horizontal 2D channel

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: 8
Arnoldinho is on a distinguished road
Robert, regarding your 2D channel, is this what you mean ? -> Different boundary conditions between OF-1.6 and OF-1.7.1 for interFoam
Arnoldinho is offline   Reply With Quote

Old   October 16, 2011, 22:59
Default
  #6
Member
 
Albert Tong
Join Date: Dec 2010
Location: Perth, WA, Australia
Posts: 76
Blog Entries: 1
Rep Power: 6
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, 02:48
Default
  #7
Senior Member
 
Arne Stahlmann
Join Date: Nov 2009
Location: Hanover, Germany
Posts: 209
Rep Power: 8
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, 05:35
Default no ourflow with interFoam
  #8
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 185
Rep Power: 7
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, 06:25
Default
  #9
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 185
Rep Power: 7
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, 365 views)
File Type: jpg impactOnOutflowT0,15_withXCoordsmallerZero.jpg (14.6 KB, 43 views)

Last edited by vonboett; June 14, 2012 at 09:09.
vonboett 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
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 08: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 03:27


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